ANSYS ls-dyna在节理岩石爆破裂纹损伤数值模拟中的应用：从建模到后处理的全方位解析与...

YrqnxehxDo

279人浏览 · 2026-03-29 15:30:00

YrqnxehxDo · 2026-03-29 15:30:00 发布

ANSYS/ls-dyna包含不同倾斜角度节理岩石爆破裂纹损伤数值模拟 1.CAD-ANSYS模型信息化建立，从原理角度进行建模，模型建立简单化，自由度高。 2.网格优化处理，网格设计技巧，实现最优裂纹效果。 3.节理创建、材料参数、边界条件等定义，快速完成关键字的定义。 4.不同角度节理修改、实用ls-prepost前处理技巧。 5.后处理云图数据操作、出图技巧，输出各类云图、裂纹演化图。课程对该案例模型从建模到后处理全过程进行了讲解，过程中还演示了许多与案例相关的ls-prepost实用小技巧。附件中资料齐全，包含爆破常用资料、软件操作指南、材料参数、软件安装等，适合入门及想深入了解软件的同学学习。

在岩土爆破模拟的坑里摸爬滚打两年多，总算搞明白了怎么用LS-DYNA玩转带节理的岩石模型。今天咱们直接上干货，手把手教你怎么用不同倾斜角度的节理模型模拟爆破裂纹损伤——别被专业名词吓到，其实操作比吃火锅还简单！

建模就像搭积木

打开CAD画个20x20x20cm的岩体模型，记住一定要用布尔运算切割出节理面。这里有个骚操作：把节理面参数写成EXCEL表格，用APDL脚本批量导入ANSYS。这样下次要改倾斜角度？直接改表格数据重新跑脚本，模型自动更新，比手动调参快十倍！（别问我怎么知道的，都是熬夜改模型的血泪经验）

!APDL批量生成节理面示例
*do,i,1,5    !生成5组不同角度节理
  plane_angle=30*i-15  !角度递增
  csys,1      !激活局部坐标系
  local,11,cyl,,,plane_angle  !创建倾斜坐标系
  vsbw,all,,sepo       !切割岩体
*enddo

网格就该这么玩

爆破裂纹最怕网格瞎划分，分享个独门口诀："四面体粗、六面体细、节理处加密"。在LS-PrePost里用快捷键Ctrl+Shift+M调出网格质量检查，把雅可比系数控制在0.7以上。偷偷告诉你，在节理面附近用sweep网格划分法，配合2mm的局部加密，裂纹扩展效果直接拉满！

材料参数别踩雷

JH-2模型参数设置是重灾区，这里给出爆破老司机都在用的黄金配方：

*MAT_JOHNSON_HOLMQUIST_CONCRETE
$#    mid        ro        gi        ...（其他参数按实际填写）
  2.65e-9   0.67e11   !密度2.65g/cm³，剪切模量67GPa
*EOS_GRUNEISEN        !状态方程必选这个
$#   eosid        c         s    gammao    ...
      2    4500e3    1.35      0.5

特别注意损伤参数D1~D4别超过0.04，否则模型直接碎成渣（别问，问就是试过）。

前处理神操作

遇到要改30°→45°节理角度？不用重建模型！LS-PrePost里选中节理面，F12调出旋转面板，直接输入角度值。更骚的是用Python脚本批量修改关键字：

with open('joint.k','r+') as f:
    lines = f.readlines()
    for i in range(len(lines)):
        if 'INCLINATION_ANGLE' in lines[i]:
            lines[i] = lines[i].replace('30','45')  #30度改45度
    f.seek(0)
    f.writelines(lines)

后处理装逼指南

爆破完别只会看应力云图，试试这些操作：