做LS-DYNA水射流破岩模拟的兄弟应该都有这种体会写K文件不算难真正磨人的是调试。尤其是水射流这种动静耦合的奇葩工况水是流体岩石是固体俩介质在一个模型里互相撕扯稍微一个参数没配对算着算着就给你抛个error——不是负体积就是时间步塌缩要么就是水粒子穿透岩石飞得满天都是。最近我刚好在搞固定/移动水射流破岩的数值模拟把几个关键代码段拆开揉碎了重新过了一遍踩了不少坑也翻出了不少平时不太会注意的魔鬼细节。这篇就写给做同类模拟的同仁做个参考从模型搭建、材料参数到动静耦合的代码细节再到调试排错一条龙说清楚。1. K文件调试到底在调什么1.1 K文件不是代码但调试逻辑和调bug一样先说个基础认知。K文件是LS-DYNA的关键字输入文件从*KEYWORD开头后面跟着一堆关键字卡片描述节点、单元、材料、接触、边界、载荷、输出控制这些信息。它本质上是模型的“图纸”不是可执行程序所以我们常说的K文件调试和写代码编译跑通的逻辑还不完全一样。代码调试是找程序的bugK文件调试则是让模型在三个层面上同时成立语法层关键字有没有写错卡片参数有没有放对列单元编号有没有重复这些属于“能不能跑起来”的问题。物理层材料参数是否真实合理边界条件是否符合试验工况单位制是否统一这些属于“算得对不对”的问题。数值层时间步是否稳定沙漏能是否超标质量缩放是否引入了过多虚假质量这些属于“算得稳不稳”“结果可信不可信”的问题。很多新手容易卡在第一层但真正的麻烦往往在第二层和第三层。K文件调试更像是做体检——你得能分辨出哪个参数异常了、哪个环节在拖后腿而不是看到一个error就慌。1.2 水射流破岩的动静耦合为什么是硬骨头水射流破岩这个工况难就难在它把两种性质完全不同的介质强行放在了一起。水是流体可变形、可流动、能飞溅岩石是固体有强度、有刚度、能破坏。两者之间还有高速冲击、应力波传播、裂纹萌生扩展、岩屑剥离这一连串物理过程。用LS-DYNA做动静耦合本质上是在求解一个瞬态流固耦合问题。水射流以几百米每秒的速度撞向岩石表面冲击区域的应力瞬间就能超过岩石的抗压强度岩石在局部发生压碎破坏紧接着水进入裂纹尖端形成水楔效应加剧裂纹扩展破裂的岩屑又被高速水流冲走。这个过程时间尺度在毫秒级甚至微秒级空间尺度从微米级裂纹到厘米级的破碎坑都有涉及。从数值模拟角度看难点集中在三个方面一是流体介质的极度变形需要用SPH光滑粒子流体动力学或ALE任意拉格朗日-欧拉方法描述二是岩石材料在高应变率、高压下的本构关系和失效准则是出了名的难定参数三是流体和固体之间的耦合界面怎么处理才能既让水“漏”不进岩石里又让岩石的破坏能真实传递到流体上。这些坑逐个拆开说都够写一篇。2. 模型体系与材料参数选型2.1 单位制第一道魔鬼关卡调K文件踩得最频繁的坑就是单位制。LS-DYNA不像某些软件那样强制你选单位制它只看数值所有计算都是无量纲的。但你输入的时候必须保持一致否则算出来的结果轻则差几个数量级重则直接发散。水射流破岩这种工况我强烈建议用mm-ms-kg体系。在这个体系下长度单位是毫米时间单位是毫秒质量单位是千克对应的导出单位如下物理量单位数值关系长度mm1 m 1000 mm时间ms1 s 1000 ms质量kg1 t 1000 kg速度mm/ms1 mm/ms 1 m/s加速度mm/ms²1 mm/ms² 1000 m/s²密度kg/mm³1 kg/mm³ 1e6 kg/m³应力/压力GPa1 N/mm² 1 MPa 1e-3 GPa能量N·mm1 J 1000 N·mm注意看速度那一行在mm-ms-kg体系下1 mm/ms恰好等于1 m/s这是个很方便的巧合很多老手在设置射流初速时都直接写数不用换算。但一旦混用了单位比如长度用了mm、密度却填了kg/m³那后面怎么调都白搭因为岩石强度、水密度、声速这些参数全错位了。我自己的习惯是在K文件头部的注释区$开头写清楚单位体系把密度、弹性模量、强度这几个容易被搞混的量换算好然后建一个参数清单表贴在K文件里调试时对着表查。另外*CONTROL_UNITS这个卡片在部分LS-DYNA版本支持单位换算但它治标不治本手动保证参数自洽永远是第一位的。2.2 水的材料MAT_NULL EOS_GRUNEISEN水射流里的“水”怎么建模低速问题可以直接用*MAT_NULL加一个声速截断但高速射流冲击岩石时水的可压缩性、冲击波传播都不可忽略必须用状态方程EOS描述压力-密度关系。常用组合是*MAT_NULL*EOS_GRUNEISEN。Null材料本身只有密度和压力截止压力截止设成负值表示允许拉伸。Gruneisen状态方程则给出水的压力响应*MAT_NULL $ MID RO PC MU TEROD CEROD YM PR 1 1.0E-06 -1.0E-06 0.0 0.0 0.0 0.0 0.0 *EOS_GRUNEISEN $ EOSID C S1 S2 S3 GAMMA0 A E0 1 1647.0 1.921 -0.096 0.0 0.35 0.0 2.9E-04在mm-ms-kg体系下水的密度RO 1.0e-6 kg/mm3也就是1000 kg/m³。Gruneisen参数里C是声速截距水大约是1647 mm/ms这就是理论上水中的声速约1647 m/sS11.921是一阶冲击波参数S2-0.096是二阶项GAMMA00.35是Gruneisen系数。这些参数参考的是标准水的冲击Hugoniot数据初次使用可以直接按这个填基本没问题。初动能E0的初始内能这个参数经常被忽略。E0是单位初始体积的内能对可压缩流动影响很大我在多个模型里试过E0设为0和设为2.9e-4在远场压力上会差几个百分点靠近冲击区差异更明显。要跟试验对标时别省这一步。2.3 岩石材料从弹性到HJC的取舍岩石材料是另一个大坑。如果只用MAT_ELASTIC那就完全无法模拟破坏用MAT_JOHNSON_COOK的话又缺少对静水压力敏感性的描述。水射流破岩这个问题里岩石处在三向受压、高应变率、大变形条件下我常用的是HJC模型*MAT_JOHNSON_HOLMQUIST_CONCRETE它对混凝土和岩石类材料的冲击压碎、损伤累积、压溃-软化行为描述得比较到位。HJC的关键参数如下*MAT_JOHNSON_HOLMQUIST_CONCRETE $ MID RO G A B C N FC 2 2.6E-06 14.8 0.79 1.60 0.007 0.61 4.0E-05 $ T EPS0 EFMIN SFMAX PC UC PL UL 4.0E-06 1.0E-06 0.0100 7.0000 1.6E-05 8.0E-06 1.0E-03 1.1E-03 $ D1 D2 K1 K2 K3 FS 0.040 1.0000 85.0 -171.0 208.0 0.0RO是密度2.6e-6 kg/mm³对应2600 kg/m³G是剪切模量单位是GPaFC是单轴抗压强度4.0e-5 GPa就是40 MPa这算中等强度岩石T是抗拉强度通常取FC的十分之一左右D1、D2是损伤常数控制材料在塑性应变累积下的劣化速度K1、K2、K3是压力-体积应变多项式系数。HJC有个特点它的归一化抗压强度公式里含了C和N两个参数分别控制应变率效应和压力硬化指数。水射流冲击的应变率可以达到10³/s以上C0.007配合应变率项能明显提高动态强度这符合岩石类材料的率相关特性。关于参数标定岩石种类千差万别HJC参数最好用SHPB分离式霍普金森压杆、三轴压缩等试验数据标定没有试验条件的话可以参考同类岩石的文献值但一定要做敏感性分析。我调过几次之后发现FC、T、SFMAX这三个参数对破碎坑形态影响最大先用它们做基准拟合后面再动损伤相关参数。2.4 网格尺度与射流速度的匹配关系水射流模拟里水的粒子大小和岩石的单元尺寸不是随便画的。一般的经验是射流直径方向至少要有10个SPH粒子岩石在冲击区域的最小单元尺寸要和粒子尺寸接近避免出现“粗岩石单元吃不消细粒子冲击”的数值穿透。举个例子如果射流直径是2 mm粒子间距取0.2 mm左右岩石冲击区单元尺寸也控制在0.2 mm级别远离冲击区可以用渐变网格放大到1-2 mm既保证精度又控制计算量。射流速度方面普通水射流在100-300 m/s超高压水射流水刀可以到600-900 m/s。速度越高冲击压力越大越需要更小的时间步长。用Courant条件简单估算一下水声速约1500 m/s粒子间距0.2 mm时间步长约0.13微秒也就是0.00013 ms。几毫秒的计算时长就要算上万个循环网格一细化计算量直接指数级上升。这也是为什么水射流模拟要在精度和效率之间反复权衡。3. 动静耦合关键代码段逐段拆解3.1 控制卡片时间步、终止、SPH三板斧K文件调试最先要看的就是头部控制卡片这些卡片决定整个求解过程的“气质”。*CONTROL_TERMINATION $ ENDTIM ENDCYC DTMIN ENDNEG ENDMAS NOSOLVR 5.0E-01 *CONTROL_TIMESTEP $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MS1ST 0.0 0.600 0 0.0 -1.0E-06 0.0 0.0 0.0 *CONTROL_SPH $ NCBS BOXID DT START VMAX CONT IDIM MEMORY 100 0 0.1 0.0 2.0E-03 0.0 2 0*CONTROL_TIMESTEP里的TSSFAC是时间步缩放系数SPH和瞬态冲击问题不推荐用默认的0.9我一般取0.6甚至0.5稳定性优先。DT2MS设成负值表示启用质量缩放负号代表“允许质量缩放”后面的1.0e-6就是目标时间步。这里要小心质量缩放是双刃剑后面单独说。*CONTROL_SPH里NCBS100表示每步最多允许100个新的粒子邻近搜索IDIM2是三维模拟的关键选项1是二维2是三维MEMORY0表示自动分配SPH粒子相关的内存。VMAX是SPH粒子的最大速度限制设太大会导致粒子飞出去无法收敛设太小又会人为截断物理速度2e-3 mm/ms就是2 m/s在mm-ms单位下的表示高压水射流600 m/s远大于这个值所以这里一般按实际最高速度稍微放大一点。3.2 SPH方法把水离散成粒子SPH是无网格方法把水离散成一个个粒子每个粒子携带质量、速度、压力等信息通过核函数插值计算粒子间的相互作用。它比ALE方法天然更适合模拟水的飞溅、破碎和自由表面大变形的场景。SPH单元的截面定义*SECTION_SPH $ SECID CSLH HMIN HMAX SPHINI 1 1.20 0.500 2.000 0CSLH1.2是粒子光滑长度的比例系数它决定核函数的支持域半径。过小会让粒子间相互作用太弱容易出现空洞过大会让结果过于平滑丢失局部冲击细节。HMIN和HMAX是光滑长度的下限和上限粒子被压缩时光滑长度自适应减小、被拉伸时自适应增大但要限制在0.5-2.0倍初始值范围内防止粒子拉伸时密度计算失真。SPH粒子的初始速度施加常用的方式是*INITIAL_VELOCITY_GENERATION*INITIAL_VELOCITY_GENERATION $ STYP EID VX VY VZ IVATN ICID 2 0 0.0 0.0 -600.0 1STYP2表示按Part施加EID0配合后面的Part列表使用。VX/VY/VZ是速度分量射流方向是Z轴负方向的话就写0, 0, -600单位mm/ms即600 m/s。这个卡片对固定射流和移动射流的区别在于固定射流可以在一开始把所有水粒子都设一个初速移动射流则要根据喷头位置和运动速度给不同时刻进入计算域的粒子设置不同的速度方向。这里有个非常容易踩的坑水射流出流孔的位置如果在模型内那些处于射流管壁内的SPH粒子如果也给了初速它们会在计算开始时直接撞向管壁引发压力震荡和粒子穿透。正确的做法是把初始速度只赋给射流入口到岩石表面这段“自由射流区”的粒子或者在射流管壁处设置让粒子“滑过”的边界条件。3.3 流固耦合核心卡片CONSTRAINED_LAGRANGE_IN_SOLID水射流破岩里如果水用SPH粒子、岩石用Lagrange网格二者之间的耦合其实是粒子-节点接触但如果水用ALE网格、岩石用Lagrange网格就需要*CONSTRAINED_LAGRANGE_IN_SOLID简称CLIS把两者绑在一起。它是LS-DYNA流固耦合的核心卡片也是动静耦合调试里最折磨人的一张卡片。*CONSTRAINED_LAGRANGE_IN_SOLID $ SLAVE MASTER SSTYP MSTYP NQUAD CTYPE DTECT MCTC 3 1 0 0 2 4 1 1 $ PFAC FRIC IGNORE NORMDIR ITYPE IGAP DELINK 0.100 0.0 1 2 1 0 1.0E-10SLAVE是Lagrange部分的Part ID这里是岩石MASTER是ALE/SPH部分的Part ID这里是水。CTYPE4是罚函数耦合通过主从面之间的相对穿透计算耦合力CTYPE5则在罚函数基础上叠加了法向压力更适合水对结构冲击的工况。DTECT1表示自动检测耦合方向MCTC1表示在主从面之间按最近点搜索。PFAC是罚因子默认0.1。罚因子太小水会渗进岩石里出现非物理的“湿透”现象罚因子太大又会在界面处引入高频振荡导致时间步降低甚至发散。IGNORE1很重要它允许移除初始穿透。模型里水和岩石如果有初始重叠部分网格没完全贴合不加这个选项会直接报穿透错误加了之后程序会把这些初始穿透忽略掉但在物理上被忽略的穿透区域在计算中相当于已经“贴合”了所以在建立模型时还是尽量让网格对齐不要全靠IGNORE兜底。DELINK1.0e-10是解除耦合的阈值。岩石单元被侵蚀删除后被删除单元上的流固耦合载荷自动转移到剩余单元上这个阈值控制单元删除后耦合关系解除的灵敏度设太小会造成数值噪声设太大又会让新暴露的面来不及建立耦合。我在调试CLIS时最常遇到的症状是水粒子整体穿透岩石或者水在岩石表面“弹开”而不是侵入裂纹。前者通常是PFAC太小或DTECT选择错误后者通常是耦合方向设反了。遇到这类问题建议先缩小模型做一个单射流单元的二维试算把耦合参数调稳了再上三维全模型。3.4 岩石破坏侵蚀接触与失效准则的配合岩石在射流冲击下要发生破坏、删除单元、形成破碎坑这需要组合使用侵蚀接触和材料失效。*CONTACT_ERODING_SURFACE_TO_SURFACE $ SSID MSID SSTYP MSTYP SBOXID MBOXID SPR MPR 2 1 3 3 0 0 0 0 $ FS FD DC VC VDC PENCHT BT DT 0.100 0.100 0.0 0.0 0.0 0.0 0.0 1.0E20*CONTACT_ERODING_SURFACE_TO_SURFACE让接触面会随着单元删除而自动更新这是模拟冲击破碎的基本配置。面面接触的FS和FD分别是静/动摩擦系数在这里对结果影响不大保持0.1左右即可。真正决定破碎形态的是失效准则。HJC材料本身有FS参数但它只能给出基于压力或损伤的失效现实中的岩石破坏是拉压混合的所以我经常再加一个*MAT_ADD_EROSION做补充*MAT_ADD_EROSION $ MID EXCL MXPRES MNPRES EFFEPS VOLEPS RATIO SIGP1 2 $ SIGVM MXEPS EPSSH SIGTH IMPULSE FAILTM 0.0 0.0 0.1000 0.0040 0.0 0.0EPSSH0.1是最大剪切应变失效SIGTH0.004是最大主应力失效GPa单位4 MPa。失效准则设得多不一定是好事——失效太容易岩石会像豆腐一样被冲掉太苛刻破碎坑又出不来。我的调试经验是先把SIGTH按抗拉强度的1.5-2倍设EPSSH从0.2开始往下压看破碎坑形态和试验照片的吻合度再迭代。另外MXPRES最大压力和MNPRES最小压力只设一个阈值防止冲击区压力异常时单元不删导致计算崩溃。3.5 移动射流实现方式与边界处理固定射流相对简单水射流位置固定持续冲击岩石同一点。移动射流则是射流喷头沿着一个方向移动模拟实际水刀切割或钻头的作业过程。移动射流的实现有两种常见思路方式一给射流入口加移动边界。把水粒子的生成区域或ALE入口边界设为一个移动的刚体Part通过*BOUNDARY_PRESCRIBED_MOTION_RIGID给它一个恒定速度让喷头以指定速度沿轨迹移动水粒子从移动的入口不断注入。方式二射流不动岩石移动。利用相对运动原理固定水射流、让岩石试件以逆方向速度移动。这个方式收敛性更好因为水和岩石的相对速度不变数值实现也简单。方式二在实际调试中更稳我一般推荐先用它做参数研究。关键代码是给岩石底部节点集合加一个恒定速度*BOUNDARY_PRESCRIBED_MOTION_SET $ NSID DOF VAD LCID SF VID DEATH BIRTH 10 2 2 0 500.0 0 1.0E28 0DOF2表示Y方向移动VAD2是加载速度时间导数SF500.0就是移动速度单位mm/ms即0.5 m/s。实际切割速度通常在这个量级100-500 mm/s即0.1-0.5 m/s。移动射流的边界处理有个关键细节岩石网格范围要足够长移动试件时不能让射流跑到岩石边界外面去。我在初版模型里吃过亏射流还没“切完”就跑出了岩石右边界结果边界附近单元变形畸变结果完全不可用。后来把岩石长度方向加长了同时把两端边界设成非反射边界*BOUNDARY_NON_REFLECTING才算把端部效应降下来。4. 魔鬼细节调试实录与排查思路4.1 负体积和粒子飞溅最容易让人崩溃的错误SPH粒子在高速冲击下如果粒子间距被压得过小或者粒子飞出主计算域经常会报negative volume错误。这个错误在Lagrange网格中表现为单元畸变在SPH里则表现为粒子密度异常、光滑长度塌缩。我遇到最多的情景是水粒子以600 m/s速度撞上岩石后在坑内高速回流、粒子密集堆叠SPH核函数计算出的密度出现非物理峰值然后程序直接崩溃。排查步骤供参考先看报错位置是不是在射流冲击区如果是说明粒子初始速度或粒子尺寸分布有问题把粒子间距放大10%-20%再试。检查*CONTROL_SPH的VMAX设置如果设得太小粒子飞溅速度被截断会导致堆积。适当提高*CONTROL_TIMESTEP中TSSFAC的稳定裕度从0.9降到0.6很多时候负体积只是时间步过大导致的数值失稳。给水粒子之间加*CONTACT_SPH自接触能缓解粒子过度穿透的现象但要小心自接触会增加不少计算量。4.2 时间步塌缩一杯“鸡尾酒”引发的宕机LS-DYNA显式计算中时间步由最小单元尺寸和材料声速决定。水射流模型里SPH粒子和岩石单元尺寸都很小时间步通常在微秒量级但只要局部有个畸变单元或速度过高的粒子时间步就会雪崩式缩小出现一根红得发黑的曲线——老玩家叫它“cocktail”鸡尾酒因为时间步曲线长得像一杯分层鸡尾酒。应对“cocktail”的常规手段用*CONTROL_TIMESTEP的DT2MS设定时间步下限负值时启用质量缩放。但注意质量缩放的副作用是把虚假质量加到临界网格上可能改变冲击波传播特性所以我一般把质量缩放限制在模型总质量的2%以内可以通过*DATABASE_GLSTAT输出的total mass变化检查。在*CONTROL_TERMINATION里配合DTMIN设置最小时间步低于这个值程序直接终止避免无意义的死循环。优化网格岩石冲击区细网格、非冲击区粗网格中间用平滑过渡别放突变。网格尺寸突变本身就是时间步问题的引信。4.3 沙漏模式算完了才发现结果全废沙漏hourglass是低阶单元在单点积分下产生的零能模式。岩石在冲击下会经历大变形如果沙漏能占比超过总内能的10%破碎坑形态、裂纹走向基本都不可信。这个坑很隐蔽因为我见过不少人在曲线图上看应力云图觉得“看起来挺对”实际沙漏能早就超标了。排查沙漏要盯住*DATABASE_GLSTAT输出的沙漏能hourglass energy和*DATABASE_MATSUM里各Part的沙漏能。如果超标可以在岩石Part上设置合适的沙漏控制方式*CONTROL_HOURGLASS $ IHQ QH QB QS Q1 Q2 QB QW 4 0.050 0.0 0.0 0.0 0.0 0.0 0.0IHQ4是Flanagan-Belytschko刚度形式对低速冲击问题效果不错水射流这种高速冲击问题我更喜欢用IHQ6基于压力平均的粘性-刚度混合形式对抑制冲击区的沙漏更有效。QH是沙漏系数默认0.05如果沙漏能比例太高可以调到0.1但调太高会人为增加结构刚度导致破碎坑偏小。4.4 质量缩放的两面性加速和失真的权衡水射流模拟的计算成本很大时间步小、计算时长长。很多人为了省时间把DT2MS设到1e-5甚至更大速度快了但结果可能悄悄变差。质量缩放的原理是把时间步低于目标值的那些单元或粒子的密度人为调大从而降低声速、稳住时间步。但密度一调惯性就变了冲击问题中惯性力占据主导地位密度失真带来的误差会被直接放大。我实测过一组对比同一个水射流模型DT2MS分别设为5e-7几乎无质量缩放和2e-6质量缩放明显加快计算后者破碎坑深度差了约18%且岩石内部的应力波到达时间明显偏移。所以调试阶段可以适当放宽但正式出货的模型一定要把质量缩放控制在极小范围最好同时对比一两个关键结果和原始模型确认误差在工程接受范围内。4.5 固定射流和移动射流在调试上的差异固定射流调试相对简单因为载荷和边界都是稳态的出现问题容易定位。移动射流则是在固定射流的基础上叠了“时变边界”多了一层复杂度。移动射流调试排错我摸索出一个固定流程先冻结移动把移动速度设为0当成固定射流跑通确认水和岩石的参数组合没问题。再给速度从低移动速度比如0.1 m/s开始逐步提高到目标值观察破碎坑形态如何演变。重点关注射流驻留时间移动射流中射流在某一位置的驻留时间变短冲击坑可能从“深坑”变成“浅槽”这和试验趋势要一致。检查射流轨迹是否偏航喷头移动速度过快时射流本身的惯性会让冲击点滞后于喷头位置滞后量可以通过对比射流粒子轨迹和喷头位置来判断。移动射流调试最常见的问题不是崩溃而是结果不符合物理直觉。如果你发现破碎槽深浅不一、有周期性变化别急着改参数先看是不是喷头移动速度和求解步长之间有共振效应——把时间步缩小一半再算一次如果槽深分布变化很大说明是数值效应不是物理现象。5. 常见问题速查表与实践建议K文件调试的问题千奇百怪但归纳起来就那几类。我整理了一个速查表按“现象-可能原因-排查方法”列出来碰到问题可以先对着查现象可能原因排查/解决办法计算开始就报错单位制混乱、关键字参数位置错检查单位制核对关键字卡片格式和字段位置负体积/单元畸变时间步过大、网格突变调小TSSFAC优化网格过渡检查粒子间距水粒子穿透岩石罚因子太小、耦合方向错误调大CLIS的PFAC检查DTECT和NORMDIR设置破碎坑出不来岩石失效准则过严降低SIGTH/EPSSH检查HJC强度参数沙漏能超标积分方式不合适更换IHQ类型调高QH检查是否有单点积分单元时间步雪崩局部畸变单元、质量缩放失控查METSUM找罪魁祸首优化局部网格计算时间过长时间步太小、网格过细合理质量缩放非冲击区网格放大移动射流轨迹偏航喷头速度过快、射流惯性大降低移动速度步进检查滞后量边界处破碎异常边界反射、网格截断加非反射边界延长岩石范围几条实践建议都是真金白银换来的教训调试K文件不要在一开始的模型上加满所有细节。我强烈建议先做一个2D平面模型或者缩小版3D模型把射流直径、材料参数、耦合方式都验证了一遍确认关键物理过程能复现之后再逐步往全尺寸模型加东西。这能省下大量反复试错的机时。每次改参数都要留档。K文件调试经常是改一个参数、跑一遍、看曲线、再改很容易陷入“盲调”的泥潭。我习惯在每次提交前把当前K文件按“v01、v02...”命名同时在文件头的注释里记录本次改了什么、为什么改。等发现结果异常时回溯起来非常管用。最后再分享一个我最近才真正重视起来的细节输出控制卡片的设置。调试阶段我建议把*DATABASE_BINARY_D3PLOT的DT设得小一点比如每10个或20个循环输出一帧把*DATABASE_GLSTAT、*DATABASE_MATSUM、*DATABASE_SPHOUT都打开。很多问题比如沙漏能异常、质量缩放影响、时间步塌缩发生的位置在动画里根本看不出来但看能量曲线和时间步曲线一眼就明白了。数据密度够排查才能快。