MultiFracs输入数据格式及输出结果数据格式
输入“/”快速插入内容
MultiFracs输入数据格式及输出结果数据格式
MultiFracs输入数据格式及输出结果数据格式
用户2475用户24752025年9月25日修改
---------------------------------参数设置-------------------------------------
*Material=elemSetName
第一行（总共60个数据,实体单元用到的参数）
IeType,密度(初始时的天然密度),E(弹性模量),u(泊松比),pn(法向接触罚参数),ps(切向接触罚参数),接触内摩擦角,-7,kf11(渗透率),kf22,kf33,kf12,kf23,kf31,-14,Biot系数,孔隙率,应力渗透率耦合参数,突变系数,-19,kT11(热传导系数),kT22,kT33,kT12,kT23,kT31,-26,Cp(比热容),tec(热膨胀系数),htc(接触换热系数),-30,kM11(湿度传导系数),kM22,kM33,kM12,kM23,kM31,-37,msf(湿度收缩系数),-39,ft(实体单元抗拉强度),fc(实体单元抗压强度),-42,x(裂缝中冰水混合物损失质量占总冰水混合物质量占比),-44,-45,-46,-47,-48,-49,-50,-51,-52,-53,-54,-55,-56,-57,-58,-59,
第二行（总共40个数据,节理单元用的参数） 连接单元
IbkFlag,c(粘聚力,必须大于0),phi(内摩擦角),ft(抗拉强度,>0),fs(取1e15),pfn(节理单元法向罚参数，100E),pfs(切向罚参数，100E),GfI(I型断裂能),GfII(II型断裂能),-9,a0(初始张开度),amin(张开度最小值),amax(张开度最大值),-13,hTfr(水岩换热系数),-15,hTnbj(未断裂节理单元换热系数),dhTnbj(节理单元断裂时hnbj的折减系数),-18,hFnbj(未断裂节理单元换流系数),hFpof(孔隙-裂隙渗流在裂隙处的换流系数),-21,hMnbj(未断裂节理单元换湿系数),dhMnbj(断裂节理单元换湿系数的浙减系数),-24,该组节理处两侧附近三角形单元的接触内摩擦角,imat1,imat2,eflag,-29,am0(力学虚拟张开量>=0),-31,dd(3D时为倾向，2D为自x轴到层理面的夹角angle，逆时针为正,x轴水平,y轴竖直),dip(3D时为倾角,2D填写-33即可),c_min(垂直于层理粘聚力),phi_min(内摩擦角),ft_min(抗拉强度),GIf_min,GIIf_min,-39,kTc(延裂缝的热扩散系数),-41,kmc(延裂缝的湿度扩散系数),-43,-44,-45,-46,-47,-48,-49,-50,-51,-52,-53,-54,-55,-56,-57,-58,-59,
第三行（总共40个数据，实体单元塑性本构用到的参数）
CModelType(本构模型类别),K(体积模量),G(剪切模量),c(粘聚力),phi(内摩擦角),dilation(剪胀角),ft(抗拉强度),-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
(kf-渗透率张量,kT-热传导系数张量,kM-湿度传导系数张量)
(IeType=1,平面应力单元,IeType=2,平面应变单元,IeType=3,平面刚体单元,
IeType=11 4节点可变形四面体单元,IeType=12 4节点刚体四面体单元,IeType=13 10节点可变形四面体单元, Ietype=14 10节点刚体四面体单元)
(IbkFlag=1表示节理单元完全断裂,当接触开启时IbkFlag后的力学参数不起作用,当接触未开启时IbkFlag后的除pfn,pfs,内摩擦角以外的力学参数不起作用；
IbkFlag=0节理单元未断裂,IbkFlag后的力学参数均起作用
IbkFlag=-1节理单元未断裂但是允许裂隙流体流过，用于裂隙渗流等效表征孔隙渗流模型 或者 能够透水但仍然有一定强度裂缝
)
(
CModelType=0全量线弹性本构
CModelType=1增量线弹性本构
CModelType=2横观各向同性本构 二维情况下
2,E1(平行于层理面的弹性模量),E3(垂直于层理面的弹性模量),v12(平行于层理面的泊松比),v13(垂直作用于层理面的泊松比),G13(垂直于层理面的平面内的剪切模量),angle(自x轴到层理面的夹角，逆时针为正,x轴水平,y轴竖直),-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
CModelType=2横观各向同性本构 三维情况下
2,E1(平行于层理面的弹性模量),E3(垂直于层理面的弹性模量),v12(平行于层理面的泊松比),v13(力垂直作用于层理面的泊松比),G13(垂直于层理面的平面内的剪切模量),dd(层理面的倾向),dip(层理面的倾角),-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
CModelType=3莫尔库伦本构
3,K(体积模量),G(剪切模量),c(粘聚力),phi(内摩擦角),dilation(剪胀角),ft(抗拉强度),-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
CModelType=4 Drucker-Prager本构
4,K(体积模量),G(剪切模量),kphi(参数),qphi(参数),qpsi(参数),ft(抗拉强度),-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
CModelType=5 vonmises本构
5,K(体积模量),G(剪切模量),strength(屈服强度),H(塑性模量),-5,-6,-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
CModelType=101 经典粘弹性本构-马克斯威尔体
101,K(体积模量),G(剪切模量),eta(粘性系数),-4,-5,-6,-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
CModelType=102 粘弹性本构-power model
102,K(体积模量),G(剪切模量),A1(power常数),A2(power常数),n1(power指数),n2(power指数),rs1(参考应力1),rs2(参考应力2),-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,-25,-26,-27,-28,-29,-30,-31,-32,-33,-34,-35,-36,-37,-38,-39,
)
(GfI,GfII两者的初始取值可为2,10)
(对于第二行-31后的 dd(>=0),c_min(垂直于层理的粘聚力),phi_min(内摩擦角),ft_min(抗拉强度),GIf_min,GIIf_min，当
当dd>=0，后面的参数才起作用，此时-31前面的那组c,phi,ft,GIf,GIIf表示平行于层理方向的参数取值
当dd<0，后面的参数不起作用，表示没有层理
）
*JointBetweenMat=imat1,imat2
IbkFlag,c(粘聚力,必须大于0),phi(内摩擦角),ft(抗拉强度,必须大于0),fs(取1e15不用修改),pfn(节理单元法向罚参数),pfs(节理单元切向罚参数),GIf(I型断裂能),GfII(II型断裂能),-9,a0(初始张开度),amin(张开度最小值),amax(张开度最大值),-13,hTfr(水岩换热系数),-15,hTnbj(未断裂节理单元换热系数),dhTnbj(节理单元断裂时换热系数hnbj的折减系数),-18,hFnbj(未断裂节理单元换流系数),hFpof(孔隙-裂隙渗流在裂隙处的换流系数),-21,hMnbj(未断裂节理单元换湿系数),dhMnbj(断裂节理单元换湿系数的折减系数),-24,该组节理处两侧附近三角形单元的接触内摩擦角,imat1,imat2,eflag,-29,am0(力学虚拟张开量>=0),-31,dd(3D时为倾向，2D为自x轴到层理面的夹角angle，逆时针为正,x轴水平,y轴竖直),dip(3D时为倾角,2D填写-33即可),c_min(垂直于层理粘聚力),phi_min(粘聚力),ft_min(抗拉强度),GIf_min,GIIf_min,-39,kTc(延裂缝的热扩散系数),-41,kmc(延裂缝的湿度扩散系数),-43,-44,-45,-46,-47,-48,-49,-50,-51,-52,-53,-54,-55,-56,-57,-58,-59,
*Joint=lineElemSetName