ABAQUS仿真模拟源文件 三维岩石试样压裂仿真 试样尺寸d50mm,h100mm 试样参数岩石 工作条件一端固定另一端15N 注意是ABAQUS源文件此处应有应力分布示意图自己跑完仿真记得截图哈最近在折腾ABAQUS做岩石压裂模拟发现这玩意儿比想象中带劲儿。就拿这个直径50mm、高100mm的岩石圆柱体来说想用15N载荷搞出裂缝效果实操时遇到不少坑。直接上干货分享下建模关键步骤和对应代码片段。模型骨架搭建在ABAQUS/CAE里画圆柱体谁都会但inp文件里的几何定义才是真本事。看这段实体生成代码*Part, nameSample *Node 1, 0., 0., 0. 2, 25., 0., 0. ...节点坐标省略 *Element, typeC3D8R 1, 1, 2, 3, 4, 101, 102, 103, 104 ...单元连接省略 *End Part这里用C3D8R八节点线性减缩积分单元既能保证计算精度又不至于让电脑冒烟。注意节点编号规律——前四位是底面节点后四位对应顶面这种排列在后续加载时找节点集会方便很多。让石头硬起来岩石材料参数设置是门玄学实验室数据没到手的可以先用经典参数顶着*Material, nameRock *Elastic 5e9, 0.25 # 弹性模量5GPa泊松比0.25 *DruckerPrager 30., 30. # 内摩擦角30°膨胀角30° *DruckerPrager Hardening 0., 1e6 # 初始屈服应力1MPa这里Drucker-Prager模型比Mohr-Coulomb更适合三维模拟。重点说下膨胀角——这参数控制材料屈服后的体积变化设为和内摩擦角相同值能较好模拟岩石的剪胀现象。锁死底座玩加载ABAQUS仿真模拟源文件 三维岩石试样压裂仿真 试样尺寸d50mm,h100mm 试样参数岩石 工作条件一端固定另一端15N 注意是ABAQUS源文件边界条件处理有个骚操作不用直接固定底面而是创建参考点施加约束*Boundary BaseSet, 1, 6 # 约束底面所有自由度 *Cload TopRef, 2, -15 # 在顶部参考点Y方向施加-15N载荷这里用参考点加载比直接选面更稳定。注意力的方向是Y轴负向所以数值是负的。有兄弟问为啥不用位移加载——对于脆性材料力加载更容易触发失稳破坏这个后面再说。网格的刀工划网格时在中间区域加密用这个骚操作*Nset, nsetMidZone, generate 501, 600, 1 # 中间100个节点 *Elset, elsetMidElements, generate 201, 300, 1 # 中间100个单元就像切香肠时中间多切几刀裂缝更容易从这里开始。单元尺寸建议控制在5mm左右太细了计算时间指数爆炸太粗了裂缝路径会跑偏。让电脑开始燃烧分析步设置暗藏杀机*Step, nlgeomYES *Static 0.1, 1., 1e-5, 1. *Output, field, variablePRESELECT *Output, history, frequency10把最大增量步设为1秒初始增量0.1秒。重点在nlgeomYES打开几何非线性——别小看这个开关岩石在破坏前会有微小变形这个设置能让接触算法更准确。跑完仿真打开ODB文件用这个Python脚本快速提取最大主应力from odbAccess import openOdb odb openOdb(fracture.odb) lastFrame odb.steps[Step-1].frames[-1] stress lastFrame.fieldOutputs[S].values maxPrincipal max([s.maxPrincipal for s in stress]) print(f最大主应力{maxPrincipal/1e6} MPa)这个脚本直接读取最后一个分析步的数据适合快速判断是否达到抗拉强度。如果输出值超过岩石抗拉强度比如10MPa说明裂缝要开始发育了。翻车预警遇到过不收敛的情况试试这两招在材料定义里加*Density设置密度2500kg/m³动力分析更稳定把加载方式改成位移控制用Boundary配合Amplitude渐进加载最后提醒这个inp文件需要配合合适的材料参数才能出效果。建议先做单轴压缩试验校准参数再玩压裂模拟。源文件已打包解压后注意文件路径别带中文ABAQUS这货对中文路径过敏。伸手党直接私信拿文件
岩石压裂应力云图](https://example.com/fracture_simulation.png
ABAQUS仿真模拟源文件 三维岩石试样压裂仿真 试样尺寸d50mm,h100mm 试样参数岩石 工作条件一端固定另一端15N 注意是ABAQUS源文件此处应有应力分布示意图自己跑完仿真记得截图哈最近在折腾ABAQUS做岩石压裂模拟发现这玩意儿比想象中带劲儿。就拿这个直径50mm、高100mm的岩石圆柱体来说想用15N载荷搞出裂缝效果实操时遇到不少坑。直接上干货分享下建模关键步骤和对应代码片段。模型骨架搭建在ABAQUS/CAE里画圆柱体谁都会但inp文件里的几何定义才是真本事。看这段实体生成代码*Part, nameSample *Node 1, 0., 0., 0. 2, 25., 0., 0. ...节点坐标省略 *Element, typeC3D8R 1, 1, 2, 3, 4, 101, 102, 103, 104 ...单元连接省略 *End Part这里用C3D8R八节点线性减缩积分单元既能保证计算精度又不至于让电脑冒烟。注意节点编号规律——前四位是底面节点后四位对应顶面这种排列在后续加载时找节点集会方便很多。让石头硬起来岩石材料参数设置是门玄学实验室数据没到手的可以先用经典参数顶着*Material, nameRock *Elastic 5e9, 0.25 # 弹性模量5GPa泊松比0.25 *DruckerPrager 30., 30. # 内摩擦角30°膨胀角30° *DruckerPrager Hardening 0., 1e6 # 初始屈服应力1MPa这里Drucker-Prager模型比Mohr-Coulomb更适合三维模拟。重点说下膨胀角——这参数控制材料屈服后的体积变化设为和内摩擦角相同值能较好模拟岩石的剪胀现象。锁死底座玩加载ABAQUS仿真模拟源文件 三维岩石试样压裂仿真 试样尺寸d50mm,h100mm 试样参数岩石 工作条件一端固定另一端15N 注意是ABAQUS源文件边界条件处理有个骚操作不用直接固定底面而是创建参考点施加约束*Boundary BaseSet, 1, 6 # 约束底面所有自由度 *Cload TopRef, 2, -15 # 在顶部参考点Y方向施加-15N载荷这里用参考点加载比直接选面更稳定。注意力的方向是Y轴负向所以数值是负的。有兄弟问为啥不用位移加载——对于脆性材料力加载更容易触发失稳破坏这个后面再说。网格的刀工划网格时在中间区域加密用这个骚操作*Nset, nsetMidZone, generate 501, 600, 1 # 中间100个节点 *Elset, elsetMidElements, generate 201, 300, 1 # 中间100个单元就像切香肠时中间多切几刀裂缝更容易从这里开始。单元尺寸建议控制在5mm左右太细了计算时间指数爆炸太粗了裂缝路径会跑偏。让电脑开始燃烧分析步设置暗藏杀机*Step, nlgeomYES *Static 0.1, 1., 1e-5, 1. *Output, field, variablePRESELECT *Output, history, frequency10把最大增量步设为1秒初始增量0.1秒。重点在nlgeomYES打开几何非线性——别小看这个开关岩石在破坏前会有微小变形这个设置能让接触算法更准确。跑完仿真打开ODB文件用这个Python脚本快速提取最大主应力from odbAccess import openOdb odb openOdb(fracture.odb) lastFrame odb.steps[Step-1].frames[-1] stress lastFrame.fieldOutputs[S].values maxPrincipal max([s.maxPrincipal for s in stress]) print(f最大主应力{maxPrincipal/1e6} MPa)这个脚本直接读取最后一个分析步的数据适合快速判断是否达到抗拉强度。如果输出值超过岩石抗拉强度比如10MPa说明裂缝要开始发育了。翻车预警遇到过不收敛的情况试试这两招在材料定义里加*Density设置密度2500kg/m³动力分析更稳定把加载方式改成位移控制用Boundary配合Amplitude渐进加载最后提醒这个inp文件需要配合合适的材料参数才能出效果。建议先做单轴压缩试验校准参数再玩压裂模拟。源文件已打包解压后注意文件路径别带中文ABAQUS这货对中文路径过敏。伸手党直接私信拿文件
