ls-dyna 常见错误

bateren

作为一名LS-DYNA初学者，经常遇到各种报错，但又解决无门。这里整理一下我遇到过的报错，希望能为其它初学者提供帮助。如果有用的话，还希望点个赞。

报错查看方法：使用txt打开message文件或ls.run.out文件都可以查看。

1、*** Error 70245 (OTH+245) Memory allocation failed while allocating
问题：内存分配失败。

解决办法：内存不够。可以增大内存，减小网格数量，在DYNA左上角的内存分配处不使用“auto”，手动调小内存。



2、 *** Error 10144 (KEY+144) CHECKING VOID PARTS part set for initial void parts. # 1 is not found. This is #1 in the input order)
问题：找不到parts，我这里是流固耦合设置不对。

解决方法：我遇到这个问题时，是流固耦合的part设置的不对，本来应该选取两个component，却只选择了一个。解决方法是建立一个新的part set，然后在流固耦合中重新选择。



3、*** Error 41438 (SOL+1438) Write error on file d3plot , which=773 where= 1520128 wrote 0 of 53528. Handle =1
问题：写入文件失败。

解决方法：硬盘空间不足。加大或预留更多硬盘空间即可。



4、要求使用双精度版本，具体报错代码记不得了。
解决方法：使用双精度版本。


在SOLVER处选择双精度版本（一般在文件名中带有‘d’，我这里是lsdyna_dp.exe），并在Preset中选择双精度版本（图中选择的是单精度），再进行计算。



5. *** Error 70022 (OTH+22)
Please file a Service Request (SR) via the ANSYS Customer
Portal at https://support.ansys.com
Customers using an academic license may seek assistance
at https://forum.ansys.com
E r r o r t e r m i n a t i o n
问题：许可证有问题，我遇到的是端口冲突了。

解决方法：先打开这个：


打开之后页面大概是这样，但是license没有running:


具体解决方法可以观看以下视频，来自up主 口算N-S方程：

【ANSYS 2022R2许可证管理器失效解决方案】 ANSYS 2022R2许可证管理器失效解决方案_哔哩哔哩_bilibili



6、*** Error 40509 (SOL+509) negative volume in solid element # 1 cycle 136560 time 1.0000E+00
问题：负体积。如材料发生大变形时，可能导致单元的体积成为负值。还有一种情况，是模型本身就出了问题，例如我有一次就遇到单元的ID没有更新导致一开始算就负体积了。调试的时候一定要回到负体积发生的前一步去看发生了什么问题。

解决方法：这个问题的解决过程个人感觉比较玄学，下面的方法是在网上找到的解决方法，可以根据自己的实际模型都试一试。其中，在括号前带有^号的是我自己尝试过，应该有用的方法。

^(1)在 *CONTROL_TIMESTEP里 面设置ERODE选项为1，而且 在*CONTROL_TERMINATION里设置DTMIN项为任何非零的值，在这种情况下，出现负体积的单元会被删掉而且计算继续进行(大多数情况)。有时即使ERODE和DTMIN按上面说的设置了，负体积可能还是会发生从而导致因错误终止。

(2)简单的把材料应力－应变曲线在大应变时硬化, 这种方法会非常有效。

(3)有时候修改初始网格来适应特定的变形场将阻止负体积的形成。此外，负体积通常只对非常严重的变形情况是个问题，而且特别是仅发生在像泡沫这样的软的材料上面。

^(4)减小时间步缩放系数(timestep scale factor), 缺省的0.9可能不足以防止数值不稳定。

(5)避免用全积分的体单元(单元类型2和3)，它们在包含大变形和扭曲的仿真中往往不是很稳定。全积分单元在大变形的时候鲁棒性不如单点积分单元，因为单元的一个积分点可能出现负的Jacobian而整个单元还维持正的体积。在计算中用全积分单元因计算出现负的Jacobian而终止会比单元积分单元来得快。

^(6)用缺省的单元方程(单点积分体单元)和类型4或者5的沙漏(Hourglass)控制(刚化响应)。对泡沫材料首先的沙漏方程是：如果低速冲击type 6，系数1.0; 高速冲击type 2或者3。

(7)对泡沫用四面体(tetrahedral)单元来建模，使用类型10体单元。

(8)增加DAMP参数(foam model 57)到最大的推荐值0.5。

(9)对包含泡沫的接触，用*CONTACT选项卡B来关掉shooting node logic。

(10)使用*CONTACT_INTERIOR卡用*SET_PART来定义需要用*CONTACT_INTERIOR来处理的parts，在*SET_PART Card 1的第5项DA4来定义*CONTACT_INTERIOR类型。缺省类型是1，推荐用于单一的压缩。在版本970里，ELFORM = 1的体单元可以设置type=2，这样可以处理压缩和减切混合的模式。

(11)如果用mat_126，尝试ELFORM=0;

(12)尝试用EFG方程(*SECTION_SOLID_EFG)。因为这个方程非常费时，所以只用在变形严重的地方，而且只用于六面体单元。

此方法(12)即为传说中的“无网格伽辽金法”，其表面上看起来的效果和自适应网格比较像，本质上已经修改了模型的求解方式。

注：“缺省”的意思即为“默认”。在看LS-Dyna的教程时才第一次见到这个词，当时特别迷惑，不知道什么意思。后来用多了才知道其实就是“默认”的意思。不知道为啥要这样翻译，也不知道为啥DYNA的教程里这么喜欢用这个词。

7.初始穿透问题

在Tools中找到PenetrationCheck


在左边出现的页面中选择要检查的Component

点OK，点Check，就可以自动检查穿透问题。


上面的几个选项可以自动或手动调整修复穿透问题。

如何查找单元？
例如在 *** Error 40509 (SOL+509)

negative volume in solid element # 49221 cycle 1 time 5.2982E-05

有的时候，我们会想要查看报错中的单元具体在哪里。



通过Tool面板里的find，在单元类型（此处为nodes）处选择by id即可。

LSPP中的操作为EleTol->Find->Nd/El/Part ID


后处理查看单元距离？
LSPP中的操作为EleTol->Measur


单元质量检查
hm中为Mesh-＞Check-＞Elements


编辑于 2025-09-29 09:19・黑龙江