MARC 出错调解

更新时间:2023-11-02 03:07:01 阅读量: 综合文库 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

MARC 提供的求解器包括:直接 profile 求解器;直接稀疏求解器;稀疏迭代求解 器;单元迭代求解器等等。

本文采用 Iterative Sparse-稀疏迭代求解器:支持 MARC 的绝大多数分析功能的

求解。支持实对称、正定系数矩阵的求解,内存占用很少。支持 In-Core 和 Out-of core。 性迭代求解。

RUN 提交作业并开始运行。该模型有 20000 左右单元,计算时间视计算机硬件不同

而有差别。如果最后输出数字(EXIT NUMBER)为 3004,即表示求解成功。如果存在边 界条件给的不正确;输出的是别的数字,则代表不成功。最常见的输出数字包括:13, 1009,3002,3009,3015 等。如果有网格重划分,还会出现 33、5059、5063 等错误。 需要参照 EXIT MESSAGE 提示和 OUTPUT FILE 去发现问题和寻求解决问题的方法。然后

去反复修改模型,直至提交输出数字正确为止。 常见的输出数字可以做如下处理:

(1)输出 13 号错误,Marc 输出信息如下:

Data errors have been detected during data input. Refer to output for

location of error. Likely causes are misspelled keywords, mistyped lines or invalid input options。

如果模型中不存在网格重划分,代表几何模型有错误,这种情况下应该检查建模过 程中可能出现的失误,一般有可能是以下原因:

① 存在畸形的单元,应该重新划分,尽量将单元形状划分规则。

② 存在冗余节点或单元没有排列好,应该用 sweep 命令清除掉多余的节点;用 check 命令检查出来,然后 flip element;最后用 renumber 命令重新进行编 号。

③ 有可能是材料问题或接触刚体的问题。 ④ 约束加的不合适。

⑤ 约束或接触体的定义名称给的不合适也会出现退出号 13 的情况。

在网格重划中出现该项错误,比较常见的是因为使用了 herrmann 单元。该种单元 主要用于橡胶和金属等体积不可压缩材料,在模拟该类材料大变形过程中,要注意求解 类型的选定。如是橡胶类材料可选用 elasticity1 和 elasticity2,但更多的是选用 elasticity2,因为该模式求解效率更高。而对于金属类材料的大塑性变形问题,只能 选用 plasticity5。特别对于三维网格重划问题,目前 mentat2003 只提供了 157 号单 元,该单元就是 herrmann 单元,其求解类型只能用 elasticity2 和 plasticity5,否 则就会出现错误。

以上各种求解类型分别对应于如下情况: ①

asticity:jobs-analysis opions-rubber elasticity procedure-large

strain-total lagrange。 ②

elasticity2:jobs-analysis opions-rubber elasticity procedure-large strain-update lagrange。 ③

plasticity5:jobs-analysis opions-plasticity procedure-large strain multiplicate。

plasticity3:jobs-analysis opions-plasticity procedure-large strain additive。 36

(2)33 号错误。Marc 输出信息如下:

There are more boundary nodes than the upper bound declared. Increase the upper bound to the number of boundary nodes on the CONTACT option。

这个错误提示一般出现在网格重划分时,是因为网格重新细化以后,模型中的网格 数目和节点数目急剧增大,Job 中设置的网格重划参数小于实际中网格重划后生成的参 数,改正的方法是将 Job 中设置的网格重划参数中相应的单元数和节点数增大。即只要 将该选项中的三项(即 max element;max node ;max contanct)数值加大,在设置 这几个数值之前应该大概估计一下将会生成的单元数目和节点数目,有目的地去修改。 (3)输出 3009,Marc 输出信息如下:

The time step size has become too small to continue the analysis due to too many time step cut-backs in the current increment.

It is not possible to obtain convergence with the specified allowed maximum number of cut-backs.

以上提示代表计算结果不收敛,有可能是以下原因: ①

收敛准则选的不合适或收敛容差给的不合适,对待不同的问题应该有不同

的收敛准则,有些用位移收敛准则好,有些用力收敛准则好,有些用应变能收敛准则好; 当然收敛的容差数值也应该给的合适。 ②

有可能是步长给的不合适,应该选用自适应载荷步。使用 adaptive steps auto step 并把 parameters 中的#cut backs allow 值设置大些。 ③

减少 loadcases 总时间。 ④

增大时间步长。

(4)输出 1009,Marc 输出信息如下:

Error encountered in stress recovery. The output reveals which element has a particular problem. The error is usually due to excessive deformation in the element. This can occur during the iterative process, so that it is not always possible to visualize the excessive deformation.

Check the material behavior and the magnitude of the incremental loads (5)输出 5063,Marc 输出信息如下:

The element size for remeshing is found to be zero. You need to specify the element size or number of elements needed for remeshing. Check input.

以上提示说明网格重划中没有设置单元边的平均长度值,这个长度值是网格重划分 模型目标的平均长度,是网格能够得以重新划分成功必不可少的条件。 (6)5059 错误,Marc 输出信息如下:

The mesher failed to create a mesh. Check error messages in the output file (jidname.out) or in the log file.

网格重划分失败,有可能是网格重划分的准则或网格重划分参数设置的不好,应该

去调节,主要包括网格重划分的单元长度和总单元、总节点数目。 元建模及其仿真实现。

本文来源:https://www.bwwdw.com/article/q7v2.html

Top