考虑一个质量等于体积*密度=2*2*2*0.07=0.56的4节点壳单元。

由此得出，每个节点的质量为0.56/4=0.14 。

此外，还有两个通过*element_mass定义的集中质量：

● 在节点2处，集中质量=0.10

● 在节点3处，集中质量=0.20

这意味着：

● 节点1:0.14（节点结构质量）

● 节点2:0.14+0.10=0.24（节点结构+集中质量）

● 节点3:0.14+0.20=0.34（节点结构+集中质量）

● 节点4:0.14（节点结构质量）

此外，节点1和3通过*constrained_nodal_rigid_body定义在节点刚体中。

//在D3HSP文件中，可以找到以下消息

总质量为 0.56 + 0.10 + 0.20 = 0.86 (结构+集中质量 = 总质量)

NRB质量=2*0.14+0.20=0.48（节点结构+集中质量）

这是未包含在节点刚体中的壳体材料的质量（总质量-节点刚体的质量）

// “质量总结”部分

这是节点2处的集中质量。节点3处的集总质量包括在节点刚体中，因此不在结构的可变形部分中。

第一个数字是壳体材料质量（结构质量）；第二个数字是包括在节点刚体中的质量的部分。集中质量不包含在里面。

代号为NaN（非一个数字）的速度报告表明，由于多种原因，分析变得不稳定，追踪造成NaNs不稳定的根本原因非常困难。尽管如此，仍可以通过激活ISNAN（*CONTROL_SOLUTION），在消息文件中报告输出具有力或力矩阵列的节点（超出范围…）。

频繁地将状态数据转储到d3plot数据库会导致计算不稳定。它可以通过从D3DUMP或RUNRSF文件重新启动，并通过*CONTROL_BINARY_D3PLOT，或通过改变输出间隔来实现。频繁的绘图状态输出允许用户查看模型从最初到不稳定性的演变。

造成这个现象的根本原因是接触破裂，或由非物理负载、不真实的材料参数输入或严重的沙漏引起的严重单元变形。非物理阻尼参数也会导致不稳定发展，时间步长过大也是如此。对于预计会发生较大物理变形的零件，最好坚持默认的单元类型公式，因为这些公式往往是最稳健的。

如果将可变形点焊（梁类型9）与*CONTACT_SPOTWELD_TORSION一起使用，建议对壳体零件调用*DAMPING_PART_STIFFNESS（在阻尼系数为0.1时使用）。

如果涉及非常薄的外壳，自动接触可能会中断。在这些情况下，可以增加薄壳的接触厚度（请参见*PART_CONTACT）。

如果您怀疑复合材料会产生问题，可尝试用MAT1（*MAT_ELASTIC）代替MAT54/55，看看是否仍会出现不稳定性。

使用*PART_INERTIA关键字卡时，必须定义惯性矩（Ixx、Iyy等）。此外，代码必须知道参考坐标轴。有时，惯性值在局部坐标系和全局坐标系的坐标系中提供给分析人员也是不同的。在这种情况下，IRCS＝1允许用户直接输入这些惯性值，而无需首先将它们转换为全局系统。LS-DYNA将自动进行转换。无论在输入中使用全局参考系还是局部参考系，该参考系的原点始终位于刚体的质心处。

动态松弛不适用于一般的准静态分析。当预加载只产生小的弹性应变时，施加预加载是可以的，或者将系统初始化为规定的几何结构，但这并不是对其它方面都有好处。

您可以用常规的显式模拟进行准静态分析，根据需要调用时间或质量缩放，在合理的时间范围内得出结果，但这种方法可能很棘手。您必须关注系统中的动能，因为您想最大限度地减少惯性效应。基本上，动能相对于内能应该保持较小。（按时间缩放，我的意思是比准静态实验更快地施加载荷，以减少模拟时间。）有关质量缩放的更多信息，请参阅文件“mass_scaling”。

或者，您可以尝试使用LS-DYNA v.960进行隐式静态分析。请参阅命令*CONTROL_IMPLICIT_...，以及950用户手册中的附录M。在ls-dyna/example目录中，我们的用户ftp站点上有一些隐式分析的示例。

另请参阅

● 质量缩放

● 长运行时间

注意

初始化为指定的几何图形

1. 编写第一次运行最终状态的节点位移文件。要以必要的格式获取此数据，请按如下方式使用LS-TAURUS：

● ls-taurus g=d3plot<执行1000<进入最终状态deform<如上所述写入文件＜结束 ls-taurus

● 请注意，d3plot不包含节点旋转自由度，因此旋转自由度被写成零，这可能是壳和梁初始化的一个实际问题。

● LS-PrePost可以选择使用Output>Nodal Displacements写入位移，但输出是i8,3e16，而不是所需的i8,3e15，因此建议使用LS-TAURUS。LS-TAURUS不适用于Windows电脑。它是免费的，可用于Unix和Linux工作站。

● 如果您执行常规的动态松弛运行以达到初始化状态，DR阶段结束时将自动写入指定位移和旋转自由度的文件。

2. 在第二次运行中，快速初始化到步骤 1 中写入的指定几何形状。您需要在 *CONTROL_DYNAMIC_RELAXATION 中设置 IDRFLG=2，并在执行行中包含 m = filename，其中 filename 是步骤 1 中创建的文件。在瞬态运行开始之前，LS-DYNA 将自动运行 100 个时间步的预分析，其中节点根据 filename 中的数据发生位移。

更多内容分享，欢迎继续关注我们！

How to | LS-DYNA使用技巧（二）

How to | LS-DYNA使用技巧（三）

How to | LS-DYNA使用技巧（四）

LS-DYNA电池行业应用介绍（三）

LS-DYNA电池行业应用介绍（二）

LS-DYNA电池行业应用介绍（一）