ls-dyna *** Error 30073 (INI+73)

Question

请教各位。ls-dyna中求解k文件，发生这个错误是怎么回事呢。

*** Error 30073 (INI+73)
     negative or zero determinant - solid element 1399485
       node 1521129 ( -3.6452E+00, -2.4552E+00, -2.3000E+01)
       node 1521125 ( -3.6483E+00, -2.4643E+00, -2.3000E+01)
       node 1521126 ( -3.6342E+00, -2.4653E+00, -2.3000E+01)
       node 1521130 ( -3.6350E+00, -2.4572E+00, -2.3000E+01)
       node 1626750 ( -3.6483E+00, -2.4643E+00, -2.2500E+01)
       node 1626746 ( -3.6452E+00, -2.4552E+00, -2.2500E+01)
       node 1626745 ( -3.6350E+00, -2.4572E+00, -2.2500E+01)
       node 1626749 ( -3.6342E+00, -2.4653E+00, -2.2500E+01)

Answer 1

让阿豪来帮你解答，本回答参考chatgpt3.5编写提供，如果还有疑问可以评论或留言

这个错误提示意味着在求解k文件时，ls-dyna遇到了一个实体元素的determinant（行列式）值为负数或0的情况。这很可能是由于实体元素弯曲或畸变造成的，进而导致解方程的求解失效。这个错误提示的具体信息包括了造成错误的节点编号和坐标，在进行问题分析时可以从这些信息入手进行排查。 解决这个问题的一种方法是检查模型中的实体元素是否存在问题。可以在ls-dyna中使用反演求解器（INVERT）来定位这些节点，可以通过修改节点位置或元素连接选项来解决这个问题。另一个方法是将实体元素转换为壳元素，但这可能会导致模型的实体细节部分丢失。最好的方法是检查几何图形并手动调整其中的问题部分。 案例： 在一个计算机模拟中，ls-dyna在计算实体元素时出现了30173错误。在分析问题之前，我们首先需要了解以下具体信息：这个错误是由于实体元素的determinant（行列式）值为负数或0造成的，并且涉及到了某些节点的坐标。 我们使用ls-dyna的反演求解器进行了进一步的分析，并最终确定了造成问题的节点和元素。该元素存在严重的畸变，超出了ls-dyna的计算极限。我们手动调整了模型几何图形，以解决这个问题。 代码： 下面是一个使用ls-dyna反演求解器来定位问题节点的简单示例代码： *DEFINE_CURVE_TITLE 1 Fix Top of Column *LILBULK_TITLE $ CURVE ID X-AXIS --------------------------- Y-AXIS ---------------------------$ 1 0.0000000 -0.4000000 1.00000000 1 2.0000000 0.6000000 1.00000000 *MAT_ELASTIC .....（省略部分内容） *SOLID_TITLE 1399485 7 11216 11217 11222 11223 11224 11229 0.000E+000.000E+000.000E+00 1521129 -3.6452 -2.4552 -23.0000 0.0 0.000E+000.000E+000.000E+00 1521125 -3.6483 -2.4643 -23.0000 0.0 0.000E+000.000E+000.000E+00 1521126 -3.6342 -2.4653 -23.0000 0.0 0.000E+000.000E+000.000E+00 1521130 -3.6350 -2.4572 -23.0000 0.0 0.000E+000.000E+000.000E+00 1626750 -3.6483 -2.4643 -22.5000 0.0 0.000E+000.000E+000.000E+00 1626746 -3.6452 -2.4552 -22.5000 0.0 0.000E+000.000E+000.000E+00 1626745 -3.6350 -2.4572 -22.5000 0.0 0.000E+000.000E+000.000E+00 1626749 -3.6342 -2.4653 -22.5000 0.0 0.000E+000.000E+000.000E+00 *END *CONTROL_TERMINATION_TITLE *STOP_OPTION 0 1.1898E+001 4.2570E+003 0 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.117E+000.000E+00 0.100E+000.000E+000.000E+000.000E+00 0.000E+00 以便解决这个问题，我们需要使用ls-dyna的反演求解器，第一步就是确定错误的节点的编号和坐标。通过查看错误提示信息，我们可以看到问题涉及到了节点编号为：1521129、1521125、1521126、1521130、1626750、1626746、1626745、1626749，它们的坐标分别为（-3.6452，-2.4552，-23.0000）、（-3.6483，-2.4643，-23.0000）、（-3.6342，-2.4653，-23.0000）、（-3.6350，-2.4572，-23.0000）、（-3.6483，-2.4643，-22.5000）、（-3.6452，-2.4552，-22.5000）、（-3.6350，-2.4572，-22.5000）、（-3.6342，-2.4653，-22.5000）。 接下来，我们可以使用INVERT关键字来定位问题节点和元素。以下是相应的代码段： *CONTROL_SOLUTION_TITLE *INVERT $ Invert Title 1000 0 0 $ PLOT TITLE FACTOR NODES/ELM

    49          1           1

*END 我们可以将这个代码段放在输入文件的开头位置，从而使ls-dyna在进行求解之前就判断了问题节点是否存在。 这个示例演示了如何使用ls-dyna的反演求解器来定位问题节点。如果我们找到了问题元素，那么我们就可以对其进行调整或替换，以解决30173错误。