岩石爆破数值模拟全流程指南：LSDYNA显式分析从建模到调参实践

模拟岩石爆破说穿了就是用算法还原爆炸那一下岩石爆破这个事搞过现场的人都知道一炮下去炮孔周围的岩石从完整岩体到破碎抛出全过程只有几毫秒到几十毫秒。人眼还没反应过来应力波已经走完了一圈。而这短短一瞬决定了爆堆形态、块度分布、飞石距离和周边围岩损伤范围。想靠经验公式吃遍所有工况地质条件一变就露馅。于是计算机模拟成了绕不开的选项而在爆破数值仿真这个领域LSDYNA是绝大多数人绕不开的求解器。LSDYNA做岩石爆破模拟核心优势就一个显式中心差分算法天生就是干这个的。爆炸冲击加载下的岩石属于超高应变率问题动载持续微秒级、压力达到GPa量级隐式算法根本收敛不了而显式算法不需要迭代求解大型刚度矩阵每个时间步只要满足稳定条件就往前走正好匹配爆炸荷载这种短、猛、烈的载荷特征。这篇文章根据我自己的爆破建模经验把从建模准备、材料参数标定到边界条件设置、常见报错排查的完整链路捋一遍。适合正在做岩石爆破数值模拟的研究生、从事爆破设计的工程师以及想从零搭建爆破仿真流程的入门者。内容不绕弯子都是实操中验证过的东西你拿过去就能用。1. 为什么爆破这类高速冲击问题必须走显式这条路1.1 隐式算法在爆炸荷载面前的无力感很多刚接触仿真的人会有个疑问ANSYS经典界面里用隐式算法做静力分析也做得挺好啊为什么换个求解器要这么大费周章核心区别在时间尺度。静力问题看的是最终状态加载过程有无穷长的时间来平衡。隐式算法每个增量步里要反复迭代求解平衡方程确认结构内力跟外力对上了才肯前进。可爆炸这个问题里爆轰波扫过炸药只需要几微秒压力直接从0冲到几万个大气压之后又快速衰减。这种载荷压根不给结构平衡的机会结构的响应本身就是非平衡的——应力波还在传播岩石还没来得及变形到位。隐式算法在这种场景下的表现基本就是迭代发散、计算不收敛、单元畸变要么算不动要么算出来的结果谁都不敢信。这么说吧隐式算法像背着书包慢走的人每一步都得踩实了才迈下一步显式算法则是踩着滑板冲下坡只要脚底下能保证足够的小步幅满足稳定条件就能一路狂奔。爆破问题正好是下坡冲刺的活LSDYNA这种显式代码才是对路的工具。1.2 LSDYNA求解方式的三个关键支撑LSDYNA能扛住爆破模拟背后有几个关键技术支撑理解了这些建模的时候才知道哪些参数是动不得的。第一个是集中质量矩阵配合显式中心差分。整个模型的质量都被集中到节点上每个时间步的运动方程变成了对角矩阵不需要组装大型刚度矩阵更不需要迭代求解。程序每推进一步就是简单地做完一次加速度、速度、位移的递推算得快且内存开销小。第二个是单点积分加沙漏控制。爆炸问题网格量轻松破百万如果用全积分单元光是单元级别的积分点计算就能把算力拖垮。LSDYNA默认用单点积分但单点积分会引入一种诡异的变形模式——沙漏模式网格像搓麻将一样扭曲却又不产生任何应变能。所以必须靠沙漏控制算法把这个模态压住这个后面在控制卡片环节细说。第三个是ALE多物质组描述。爆炸和岩石接触本质上是炸药爆生气体Euler描述和固体岩石Lagrange描述之间的相互作用。炸药在爆轰过程中变形极其剧烈如果让岩石和炸药都用Lagrange网格炸药单元很快会严重畸变导致计算终止。LSDYNA的ALEArbitrary Lagrangian-Eulerian方法让炸药的网格可以独立于物质运动辅以流固耦合把爆炸压力传导给身后的岩石结构这才是爆破模拟能实现的基石。因此当你拿到一个爆破模拟任务时先别急着画网格先想清楚你要模拟的是耦合装药还是空气不耦合装药炸药近区要不要精细研究这决定了你用Lagrange接触方案还是ALE流固耦合方案。多数工程级岩石爆破模拟我都建议走ALE路线。2. 动工前的三件大事单位制、坐标系与几何取舍2.1 量纲统一cm-g-us这套单位制怎么来的这是新手最容易翻车的地方没有之一。LSDYNA的求解器本身没有内置量纲系统它只认数字。你输入一个密度值它不知道这个数字是千克每立方米还是克每立方厘米。如果你把长度用米、时间用秒、质量用吨混着来算出来的压力就会差几个数量级。爆破模拟里最常用的量纲体系是cm-g-us厘米-克-微秒。为什么是这套因为岩石密度大概在2.4~2.7水的密度是1.0炸药密度在1.0~1.6这些数字在cm-g-us体系下非常直观。弹模用GPa即10^10 g/(cm·s²)压力用Mbar也能直接对应上常用量纲表。物理量cm-g-us单位制对应国际单位制长度cmm质量gkg时间μss密度g/cm³kg/m³应力/压力MbarGPa速度cm/μsm/s换算关系记住一句话1 cm/μs 10000 m/s1 Mbar 100 GPa。JWL状态方程参数表里动不动就是A 3.7 Mbar这种你输入时候得想清楚这个数是哪套单位制来的。很多参数表直接给你GPa加mm/ms体系照抄之前一定要先完成换算。2.2 几何简化该精细的地方精细该糙的地方糙爆破模拟的几何模型并不追求长得像而是追求力学行为对。实际工程有几个炮孔就画几个孔没必要。除非你研究孔间协同作用否则单孔或双孔模型才是性价比之王。装药结构是决定几何简化的第一要素。耦合装药炸药直接填满炮孔在几何上就是岩石上一个圆柱孔、孔内一根炸药柱炸药和岩石之间可以共用节点网格或设置接触。径向空气不耦合装药炸药直径小于孔径则需要额外建一层空气环壳这层空气承担着爆生气体的传压作用网格要细一些不然爆生气体膨胀过程的压力衰减会被严重抹平。几何上还需要考虑你的计算域边界离爆源多远理论上应力波碰到边界会产生反射反射波叠加上后续的爆生气体膨胀压力会让结果失真。常见的做法有两种要么把计算域取足够大让边界反射波在关注时刻之前到不了目标区域要么在边界上施加无反射边界条件。我个人做法是分析时间在1ms以内的近区模型计算域取爆源到边界的最小距离至少4~5倍最大关注距离同时再叠加无反射边界双保险。坐标系的约定也很重要。岩石爆破模型通常是圆柱坐标或直角坐标我建议模型沿炮孔轴向对齐全局坐标系的Z轴这样后处理提取不同截面的压力云图、振速监测点布置都会方便得多。别小看这个习惯等到你需要在10万个节点里挑出某个深度截面的数据时坐标对齐能帮你省半天时间。2.3 建模顺序与零件命名习惯LSDYNA的K文件本质是一堆关键字卡片每个PART对应一个零件。前处理阶段养成规范的命名习惯后面调试K文件的时候能救命。我的命名格式一般像这样ROCK_1岩石主体EXPLOSIVE炸药柱AIR_COUPLING耦合空气层AIR_OUTER外部空气域每个PART对应独立的材料号和单元段SECTION材料号用PID的尾号同步编排这样在LS-PrePost里一眼就能看出哪个PART用的是哪张材料卡。接触定义也遵循这个规律接触名称里写明接触类型和主从面比如BLAST-ROCK_2D对应炸药与岩石的流固耦合。3. 岩石、炸药、空气三大件建模与材料卡定义3.1 岩石材料模型HJC与MAT_272的取舍岩石在爆炸荷载下的行为不是你拿一张线弹性卡就能糊弄过去的。爆炸近区岩石会经历压碎、裂隙扩展、损伤累积远区则是应力波衰减后的弹性响应。这意味着一个合适的材料模型需要同时具备以下几个特征能反映高应变率下的强度增强应变率效应能描述累积损伤引起的刚度退化能处理拉伸和压缩不对称岩石抗拉远弱于抗压破坏后残余强度的合理表达在LSDYNA里常用的两个高手一个是MAT_111HJC模型Holmquist-Johnson-Cook另一个是MAT_272RHT模型。RHT模型可以说是HJC的全面升级版它引入了更细致的失效面描述和残余强度面对岩石爆破模拟的匹配度更高。MAT_003塑性运动硬化虽然也有人用但它本质上是金属材料模型缺少压力相关的屈服行为和损伤演化算出来的爆腔形态跟实测差距很大我不建议在这种场景用。材料模型应变率效应损伤累积残余强度适用性MAT_003简单无无不推荐用于爆破MAT_111 HJC有有压缩为主可接受参数易获取MAT_272 RHT完整完整完整推荐爆破近区更准用RHT模型时有一组核心参数需要填初始密度、剪切模量、抗压强度、抗拉强度、孔隙度以及一大堆失效面参数。这些参数怎么来最好的来源是同批岩石的霍普金森杆试验SHPB数据做不同应变率下的单轴压缩测试拟合出失效面参数。但多数工程场景下你没有这个条件那就用经验参数加上文献数据校核模型标定过程中重点调整的是失效面常数B和压缩应变率指数这两个参数直接决定爆炸应力波衰减的速率。如果用了MAT_272还需要指定它的侵蚀应变。这个参数的作用是单元应变达到该值时单元被从模型中删除防止网格过度畸变导致计算步长崩溃。对于岩石爆破侵蚀主应变我一般设0.5~1.0之间压碎区设得偏低裂隙区可以稍高但注意别设太低否则岩石还没碎呢单元就没了损伤云图会糊成一片。3.2 炸药与JWL状态方程压力的来源炸药在LSDYNA里最经典的定义方式是MAT_HIGH_EXPLOSIVE_BURNMAT_008加EOS_JWL。为什么要有状态方程因为爆轰过程就是化学能瞬间转变成高温高压气体膨胀的过程要模拟出这个压力-体积-内能之间的函数关系就得靠JWL状态方程来翻译。JWL方程的标准形式是P A(1 - ω/R₁V)e^(-R₁V) B(1 - ω/R₂V)e^(-R₂V) ωE₀/V其中P是压力V是爆轰产物的相对体积当前体积比初始体积E₀是初始内能密度A、B、R₁、R₂、ω都是炸药特性参数。这些参数的来源是炸药圆筒试验测出爆速、爆压和膨胀做功数据后反求拟合出来的。爆速D、CJ压力P_CJ与状态方程参数之间不是孤立的。实际参数表里A通常是CJ压力的3~7倍B是较小的数值在GPa量级R₁在4~5之间R₂在1~1.5附近。R₁控制爆轰波前沿附近高压区的衰减速度R₂控制中等压力区的行为。参数表里A、B都是绝对值很大的数字输入错一位小数计算结果就会差得离谱。参数乳化炸药典型值TNT典型值密度 (g/cm³)1.151.63爆速 D (cm/μs)0.450.693A (Mbar)2.23.7377B (Mbar)0.0180.0373R₁4.54.15R₂0.90.9ω0.350.35E₀ (Mbar)0.0650.06还要设置一个点火点Initiation Point坐标。这里有个容易误解的地方MAT_HIGH_EXPLOSIVE_BURN并不能模拟炸药的逐步爆轰传播它只是给定一个爆炸起始位置和时间程序按爆速计算爆轰波到达每个炸药单元的时间之后按JWL方程控制压力释放。所以点火点的位置直接影响爆轰波扫过药柱的时序工程上雷管起爆点在孔底的就在药柱底部设点火点多点起爆则需要设置多个点火点序列。3.3 空气材料的定义要点空气在爆破模拟里不是没用的东西它是爆炸能量传递的介质。用MAT_NULLMAT_009配线性多项式状态方程是最常见的做法。MAT_NULL的意思是这个材料不提供剪切刚度只有体积响应正好对应气体的力学行为。线性多项式EOS的压力表达式是P C₀ C₁μ C₂μ² C₃μ³ (C₄ C₅μ C₆μ²)E其中μ ρ/ρ₀ - 1。对于理想气体C₀C₁C₂C₃C₆0C₄C₅γ-1γ为绝热指数空气取1.4即P (γ-1)E/μ。初始内能密度E₀在标准大气压下约2.5×10⁻⁶ Mbar量级。很多人会忽略这个值但如果你做的是空气冲击波场分析这个E₀大小决定了初始大气背景压力飞石轨道和冲击波到达时间的计算都会受影响。4. JWL状态方程与材料参数那些不可告人的调参内幕4.1 为什么同一款炸药不同文献里的参数不一样接触爆破模拟久了你会发现一个让人头疼的现象同一个牌号的乳化炸药在A论文和B论文里的JWL参数居然不一样有的甚至就差得离谱。原因在于JWL参数本质上是拟合产物。不同研究者做圆筒试验时测压探头的位置、炸药装填工艺、约束条件存在差异拟合出来的参数自然不同。再加上不同厂家炸药成分、密度差异同一牌号不同批次也不一样。这就要求你建立自己的参数库每次做实验或查阅文献时记录的不只是参数值还要记录对应炸药密度、试验方法、适用场景。如果你手上只有爆速和装药密度没有圆筒试验数据怎么办还有些估算途径。可以用CJ理论从爆速、密度算出CJ压力P_CJ≈ρ₀D²/4然后按经验比例关系推导A、B的初始值再通过单发模型试验比如水泥砂浆试块爆破的爆腔半径来反标定R₁、R₂。这个过程比较费功夫但对工程级计算的精度提升是肉眼可见的。4.2 初始内能怎么算一个手把手的例子不少教程里给的E₀都是现成的数但如果你换一种炸药怎么自己估算初始内能以乳化炸药为例。JWL方程中的E₀是指炸药的初始爆轰产物的单位体积内能密度工程上常用以下粗略思路估算总化学能的一部分在CJ面释放并转化为内能。取爆热Q单位质量炸药释放的化学能一般为3.5~5.5 kJ/g初始内能密度E₀ ≈ ρ₀ × Q 换算到cm-g-us单位制下。举个例子乳化炸药密度ρ₀ 1.15 g/cm³取比爆热q 4.5 kJ/g 4.5×10¹⁰ erg/g而1 erg 1 g·cm²/s²cm-g-us制中压力单位是Mbar 10¹² dyne/cm² 10¹² g/(cm·s²)能量密度单位是press单位。换算过来E₀ 1.15 × 4.5×10¹⁰ erg/cm³ 5.175×10¹⁰ erg/cm³对应Mbar量级的数值是0.05175 Mbar。这个值和经验参数表里的0.06 Mbar量级一致说明估算思路合理。注意爆轰波后的压力场不只是E₀决定的A、B、R₁、R₂共同控制整个压力衰减曲线E₀的误差影响的主要是爆生气体的膨胀后期行为。4.3 一个容易坑人的单位陷阱前面强调过cm-g-us体系但在JWL参数这里还有个大坑很多文献直接把A、B列成GPa单位比如A 373.77 GPa。如果你在LSDYNA里直接输入373.77而其他几何尺寸用的是cm那这个参数就大了10000倍因为1 GPa 0.01 Mbar。算出来的压力必然会爆表轻则单元直接负体积重则解完全发散。我自己的对照组做法是在K文件写好之后专门花十分钟把材料卡里每个参数的单位制在表里过一遍。密度对不对应该在1~3之间、弹模对不对应该几百到几千kbar级别、JWL的A/B值对不对应该几Mbar量级。这个习惯救了我很多次。5. 网格、边界条件与控制卡片决定算得准不准的三道坎5.1 网格尺寸不是越细越好爆破模拟里网格尺寸的选取是一门艺术。太粗了压力波传播被数值耗散严重抹平爆腔尺寸偏小太细了显式算法的时间步长被CFL条件锁死——时间步长跟最小单元尺寸成正比单元每小一半计算时间翻好几倍。怎么把握这个度经验法则是在一个应力波波长内至少保证6到8个单元。对于岩石爆破载荷频率主要在1 kHz到100 kHz之间在典型花岗岩中应力波速约4500 m/s到5500 m/s那么100 kHz对应波长约5 cm理论上6~8个单元对应的单元尺寸是6~8 mm。但爆炸近区压力梯度极陡近区的单元尺寸建议控制在1~2 mm远区再过渡到5 mm到10 mm。具体操作中我通常把炮孔周围1倍孔径范围内做局部加密往外逐步放大。LSDYNA前处理里用*PART和网格自动划分配合岩石用六面体网格ALE空气区域用六面体或结构化网格避免四面体带来的数值刚度偏大问题。如果实在画不出全六面体模型至少保证爆炸近区是六面体主导远区有一些楔形/四面体单元影响没那么大。网格尺寸还有一个需要跟侵蚀应变配合的点。如果岩石单元侵蚀应变设得太大单元变形拉长后时间步长会急剧缩小导致计算龟速爬行。这也是为什么侵蚀应变不能设太高的原因防止你被自己设的参数拖死。5.2 边界条件无反射边界不是万能的模型边界处理直接影响应力波反射。如果不处理边界应力波到达模型边缘会像鞭梢一样甩回去形成拉伸波可能人为制造出虚假的拉伸损伤区。这在近区爆破模拟中经常被误认为反射拉伸致裂。LSDYNA里对应的解决工具是*BOUNDARY_NON_REFLECTING。这个关键字在边界面上施加黏性阻尼力近似模拟波穿过边界进入无限域的效果。用的时候注意几点第一无反射边界对掠射角大的波存在反射残留。也就是说波近乎平行边界射过来时吸收效果会打折扣。所以即使加了无反射边界计算域边界距离爆源也不能太近。第二无反射边界对体波有效但对表面波的吸收效果很差。如果模型表面是自由面比如露天台阶爆破自由面本身是反射面那是真实的物理反射不能加无反射边界。只有截断的人工边界才需要吸收处理。第三如果模型是对称结构建议用对称边界而不是把整个模型都建出来。比如单孔耦合装药可以用1/4模型对称面上设对称约束计算量直接省四分之三。对称面不需要无反射边界因为它本身就是反射对称。5.3 控制卡片时间步、沙漏、输出频率LSDYNA的K文件里有一组控制卡片看起来不起眼实际决定了计算能不能收敛、结果能不能用。核心的几个如下时间步长缩放因子*CONTROL_TIMESTEP里的TSSFAC。显式中心差分的稳定时间步由最小单元尺寸和材料波速决定TSSFAC就是在这个稳定上限上乘一个安全系数一般取0.67默认值但很多人不知道这个值的含义。ALE部分由于是流固耦合时间步长相对更严格有时候要降到0.5左右宁可慢一点求稳。沙漏控制*CONTROL_HOURGLASS。单点积分导致的沙漏能必须在总内能的10%以内。工程上判断标准是看GLSTAT文件里的沙漏能和总能量比值超过10%说明沙漏太严重结果不可信。控制方式一般是把沙漏控制系数调到0.05~0.1并选IHQ4基于刚度的沙漏控制它对爆破这类冲击问题效果最好。注意沙漏系数不是越大越好过大的沙漏控制会引入额外刚度相当于在模型里加了隐形弹簧。输出控制*DATABASE_BINARY_D3PLOT。默认输出的时间间隔可能太粗导致你拿到的应力云图漏掉了峰值。我做爆破模拟时先跑一次粗间隔的摸底计算确认压力峰值出现的时刻再在那个时刻附近设置加密输出间隔确保峰值时刻被捕捉到。终止时间*CONTROL_TERMINATION。爆破模拟的分析时间取决于研究目标只看近区爆腔和应力波传播的话几百微秒足够看爆生气体准静态膨胀阶段和裂隙最终形态需要几毫秒。我习惯分两步走第一步算到应力波离开关注区通常200~500 μs第二步把第一步的d3dump重启动继续算到气体膨胀稳定2~5 ms。这样避免前期时间步太小导致整段计算时间过长。6. 算崩了怎么排查负体积、穿透、压力飙升的实例链路6.1 现象一负体积Negative Volume导致计算终止这是爆破模拟里出现频率第一的报错。表现为计算进行到某一步突然终止message文件里提示某个单元负体积。负体积的实质是单元被压到翻转体积变成负值这在物理上不可能属于网格畸变的极端情况。爆破模拟里最常见的原因是炸药单元在JWL压力下体积膨胀至极或者ALE网格运动导致网格翻转。我的排查链路是固定的第一步打开d3plot看负体积发生前最后几步的变形图确定出问题的位置在炸药区还是岩石区。如果在炸药区大概率是炸药单元压力过高且ALE网格没有及时重分如果在岩石区大概率是岩石单元侵蚀设置不当允许单元变形过大。第二步如果是炸药区负体积检查JWL参数。一个典型错误A输入成了GPa而不是Mbar导致压力大了10000倍。把参数重新校对一遍80%的问题出在这里。第三步参数没错那就检查网格质量。炸药近区单元长宽比不能太大扭曲度过大的单元在高压下更容易形成负体积。给炸药区做网格加密或重划分。第四步如果前三步都查不出来考虑改用ALE方法而不是纯Lagrange接触。ALE网格本身允许物质流过网格边界能从根本上避免炸药大变形下的网格畸变。6.2 现象二岩石损伤范围大面积偏大或偏小算出来的压碎区比实测大了一倍或者裂隙压根没扩展出来这是最让人头疼的结果类问题。这种问题一般不是单个参数引起的而是材料参数组合失配。排查链路先看压力云图和损伤云图。如果压力峰值的空间分布基本合理但损伤过大优先怀疑侵蚀应变设得过大——单元在失效前承受了超出实际的变形人为放大了损伤区。把侵蚀应变从1.0调到0.5压碎区通常会明显收敛。如果压力峰值分布本身就不对——比如爆腔附近的压力比CJ压力小了一个量级——那要回头检查JWL状态方程里E₀的单位是否换算正确以及点火点是否设置成功。来个很常见的错忘记设置点火点或者点火点坐标写错位置炸药单元永远不会起爆程序只能按初始压力慢慢释放压力峰值会低得离谱。如果损伤裂隙形态呈接近圆球形而实测是沿径向优势裂隙发展的非对称形态那多半没有给岩石模型设置初始各向异性或节理弱面。RHT模型本身能通过设置不同的拉伸和压缩强度来间接体现这种差异但如果岩石有显著的层理面建议考虑加入节理单元或使用能显式描述层理面的材料模型。6.3 现象三能量曲线异常——内能飙升或总能量不守恒打开GLSTAT文件看能量曲线如果总能量随时间明显上升或者沙漏能在总能量里占比越来越高说明求解过程存在问题。总能量上升往往伴随异常的质量缩放。我的处理办法是第一步在控制卡片里看*CONTROL_TIMESTEP的DT2MS设置。如果设了一个非零的DT2MS值启用质量缩放程序会人为增大某些单元的密度来放大时间步长这会导致总质量增加内能曲线异常。爆破问题宁可不做质量缩放让计算慢一点也别为速度牺牲能量守恒。第二步看沙漏能。如果沙漏能占比超过20%说明网格变形模式已经严重失控。此时除了加强沙漏控制系数还要检查网格是否存在长条薄片单元。爆破近区的单元长宽比控制在3以内是一个比较稳妥的约束。6.4 排查思路的通用化先定位层级再动手总结一下排查经验爆破模拟出问题不要一上来就改参数按参数层-网格层-接触层三层结构逐层定位。参数层的检查顺序依次是单位制、密度、弹模、JWL参数、初始内能、点火点。网格层检查顺序是单元类型、最小尺寸、长宽比、加密区域位置。接触层检查顺序是主从面定义、接触类型、接触刚度系数、侵蚀设置是否与材料匹配。这个三层排查法我用了很多年基本能覆盖九成以上的爆破模拟问题。具体到每个参数怎么改、改成多少需要配合计算结果和实际经验来微调但至少能让你的踩坑过程从盲人摸象变成按图索骥。7. 算完之后后处理里必须看的三个量爆破模拟的后处理不光是看看彩色的云图爽一下需要做定量验证。我在每个爆破模型算完后的标准检查动作是看三个量第一个是压力时程曲线。在距离爆源不同距离比如5倍、10倍、20倍孔径布置监测点输出压力随时间变化的历史曲线。峰值压力随距离的衰减规律应该符合近区陡降、远区趋缓的形状。跟经验公式或实测值对比误差在30%以内说明参数基本靠谱超过这个范围就要回到参数标定环节。第二个是岩石损伤云图的体积统计。在LS-PrePost里设定损伤变量的阈值通常以损伤变量D0.2视为损伤统计损伤区体积换算成等效损伤半径。这个值可以和现场爆破后的声波测试孔里得到的松动圈半径做对比。这是我做参数反标定最常用的一组数据。第三个是质点振速PPV。爆破振动的安全评估主要看质点峰值振速后处理中提取地表或保护目标的振动速度换算成PPV后与安全规程里的控制标准对比。这个量也常用来做远区模型的校核。最后一个提醒爆破模拟毕竟是数值试验它的价值在于做方案对比和机理分析而不是抠绝对精度。同一套建模流程和参数算几个不同装药结构的方案看它们的相对差距这个相对结论的可靠性比单个方案的绝对数值高得多。这也是我在实际项目里最推荐的工作方式——仿真数据的意义在于告诉你A方案比B方案好多少而不是告诉你爆堆飞出去多远的精确数值。我个人的体会是LSDYNA岩石爆破模拟的入门门槛不低但真正的难处从来不是软件操作而是对爆炸力学、材料本构和数值方法的综合理解。把单位制、材料参数、边界条件这三件事做扎实了大部分问题都能在前处理阶段避免。如果建模的时候贪快省事后处理阶段迟早要把这些省下来的时间加倍还回去。