ANSYS/LS-DYNA 隧道光面爆破数值模拟（CAD+LS-DYNA）课程说明：模型建立、...

ANSYS/LS-DYNA隧道光面爆破数值模拟（CAD+LS-DYNA）的课程说明 1.通过CAD完成光面爆破模型的建模导入ANSYS划分过渡网格，大量减少网格数量和网格划分时间。 2.炸药part分区后如何设置延期时间，ls-prepost实用前处理操作技巧。 3.后处理输出应力云图、损伤轮廓、时程曲线等。

隧道光面爆破数值模拟实战：从建模到后处理

搞隧道爆破模拟的兄弟应该都懂，LS-DYNA这玩意儿虽然功能强，但网格划分和参数设置要是没整明白，分分钟能把人逼疯。今天咱就唠唠怎么用CAD+ANSYS+LS-PrePost一条龙搞定光面爆破模拟，重点解决网格效率、炸药延期设置和后处理这三个坑点。

一、CAD建模：别让网格拖后腿

光面爆破模型最大的特点就是掏槽眼和周边眼的密集布孔，直接用ANSYS画六面体网格？怕不是得等到天荒地老。这里有个骚操作：先用CAD把隧道轮廓和炮孔位置画利索，导出为.step格式，再丢进ANSYS里用扫掠网格（Sweep）搭配局部细化。

举个例子，导入模型后，用ANSYS的APDL命令控制过渡区网格密度：

! 定义炮孔周边为局部细化区域 ESIZE, 0.02 ! 全局网格尺寸20mm LESIZE, 25, 0.01 ! 指定边25的网格尺寸10mm VMESH, ALL ! 执行体网格划分

这波操作的核心是“该密的地方密，该糙的地方糙”。比如炮孔周围用10mm网格捕捉应力集中，而远离区域的网格直接放大到20mm，整体网格量能减少30%以上，计算时间直接砍半。

二、炸药分区与延期起爆：K文件里的小心机

光面爆破的核心是不同炮孔的起爆时间差，但LS-DYNA里可没有“右键设置延期”这种好事。得手动在K文件里搞事情。

ANSYS/LS-DYNA隧道光面爆破数值模拟（CAD+LS-DYNA）的课程说明 1.通过CAD完成光面爆破模型的建模导入ANSYS划分过渡网格，大量减少网格数量和网格划分时间。 2.炸药part分区后如何设置延期时间，ls-prepost实用前处理操作技巧。 3.后处理输出应力云图、损伤轮廓、时程曲线等。

假设炸药按区域分了3组（Part 101~103），每组间隔25ms起爆，对应的关键字长这样：

*INITIAL_DETONATION $ 单元组ID 起爆时间(s) 101 0.0 102 0.025 103 0.05

但这里有个坑：炸药单元必须提前用*SET_PART分组！在LS-PrePost里按住Shift多选炮孔单元，用Main Menu > Group > Create直接打包成不同Part，再在K文件里引用。

如果想更骚一点，可以用ALE多物质算法避免网格畸变：

*SECTION_SOLID_ALE 2 ! 算法类型 *MAT_HIGH_EXPLOSIVE_BURN 3 1.6 ! 炸药材料参数

ALE适合大变形，但算得慢；如果是小规模模型，用拉格朗日+JWL状态方程也行。

三、后处理：云图不够，曲线来凑

算完了别急着关软件，LS-PrePost的Ctrl键全家桶才是精华。

1. 应力云图：按Ctrl+B调出截面工具，切个隧道纵剖面，再Ctrl+F调出应力分布。想让图更专业？在Fringe Level里勾选Average，避免锯齿状云图。

1. 损伤轮廓提取：在命令行直接敲：
   bash
splot dmg all # 显示所有单元损伤因子

   阈值设0.99（完全损伤），用Cut Plane工具切出爆破漏斗轮廓，导出为.dxf给CAD比对设计误差。

1. 时程曲线：比如监测拱顶位移：
   k
DATABASENODOUT
0.0001 ! 输出时间间隔
DATABASEHISTORY_NODE
1 ! 输出节点历史
1148 ! 节点编号

   算完后进History > Time History，选节点1148的Y向位移，直接导出Excel做振动分析。

最后说句大实话：

光面爆破模拟最怕“理想很丰满，现实很骨感”。CAD导进来的模型经常有破面，ANSYS里记得用Tool > Check Model修几何；延期时间设错了？直接K文件里全局替换时间参数重算，别头铁反复点界面。数值模拟嘛，有时候糙点反而更快出活——毕竟甲方爸爸只关心你的图和实际爆破效果像不像，对吧？