搜索

曹金凤老师:从零开始学Abaqus到接触非线性分析高阶应用实用小技巧

  • 分类:行业动态
  • 作者:
  • 来源:
  • 发布时间:2022-05-19
  • 访问量:0

【概要描述】Abaqus 是国际上先进的大型通用非线性有限元分析软件,它在技术、品质以及可靠性等方面具有非常卓越的声誉,拥有世界上最大的非线性力学用户群。对于工程中各种线性和非线性问题,Abaqus 都能够提供完美的解决方案,它每年都会更新版本,注入新“血液”、增加新功能,引领着全世界非线性有限元技术的发展。

曹金凤老师:从零开始学Abaqus到接触非线性分析高阶应用实用小技巧

【概要描述】Abaqus 是国际上先进的大型通用非线性有限元分析软件,它在技术、品质以及可靠性等方面具有非常卓越的声誉,拥有世界上最大的非线性力学用户群。对于工程中各种线性和非线性问题,Abaqus 都能够提供完美的解决方案,它每年都会更新版本,注入新“血液”、增加新功能,引领着全世界非线性有限元技术的发展。

  • 分类:行业动态
  • 作者:
  • 来源:
  • 发布时间:2022-05-19
  • 访问量:0
详情

Abaqus的非线性分析包括几何非线性、边界非线性和材料非线性等,其中接触非线性分析最为复杂,原因是分析过程中接触状态的变化都会引起约束条件的变化,使得搜索满足所有条件的收敛解的难度大大增加。

 

  • 数值奇异(numerical singularity)

 

缺乏边界条件约束是 Abaqus / Standard 静力分析中最容易犯的错误之一,这时往往会在 MSG 文件中出现数值奇异(numerical singularity)或零主元(zero pirot)的警告信息,此处将重点介绍数值奇异的相关知识,应该高度重视这方面的问题。

 

1 定义了位移载荷的模型

 

1所示模型是一个二维平面应变问题,包含两个弹性体:圆筒和平板。在圆筒中心的圆孔内壁上定义了固支边界条件,在平板顶部中央的 A 点上给定了位移 U2 = -2,希望使平板向正下方移动,和圆筒发生接触。

 

提交分析后,计算可以完成,但在分析结果中看到平板发生了异常的位移,偏移到了一侧(如图2所示),这是什么原因引起的?

 

2后处理时看到平板发生了异常的位移

 

  • 错误原因解析

 

要养成好的建模习惯,在提交分析作业后一定要查看 DAT 文件中是否出现了错误信息(error),MSG 文件中是否出现了警告信息(warning)。对于 Abaqus / Explicit 分析,还应查看 STA 文件中的信息。

 

打开本实例的 MSG 文件 Job-Singularity-Displacment-Wrong.msg,可以看到下列警告信息:

 

 

其含义为:平板(PLATE-1)在自由度1和2上都出现了数值奇异。这说明平板在这两个方向上都缺乏边界条件的约束。

 

对于三维模型,每个部件都有3个平动自由度和3个转动自由度;对于二维模型,每个部件都有2个平动自由度和1个转动自由度。在建立静力分析模型时,必须在模型每个实体的所有平动和转动自由度上定义足够的边界条件,以避免它们出现不确定的刚体位移,否则会在 MSG 文件中看到上述警告信息,这时分析往往无法收敛,即使能够收敛,分析结果也往往是错误的。

 

本实例中两个部件实例的边界条件为:

 

1)圆筒:在多个节点上定义了固支边界条件,不会出现刚体位移;

2)平板:在 x 方向上没有定义任何边界条件,因此在 x 方向上的刚体位移是不确定的;在 y 方向上,只在一个节点(A点)上给定了位移 U2,这时整个平板仍然可以绕 A 点做刚体转动,即除了 A 点之外,平板上的其他各个节点的 U2 都是不确定的。

 

虽然整个模型并没有施加让平板发生转动或 x 方向平动的载荷,直观感觉上该模型的边界条件设定是正确的,但这样的模型不符合有限元分析的要求。这种“因为没有受力,所以不会移动”的因果关系,只是我们根据生活经验在头脑中进行逻辑分析时的思路,而 Abaqus / Standard 的求解过程恰恰与此相反,其过程是:迭代尝试各种可能的位移状态,检验它们是否能够满足静力平衡方程。

 

在本实例中,无论平板发生多大的转动或 x 方向的平动,都可以满足静力平衡方程,即符合静力平衡条件的位移解有无限多个,因此会出现“数值奇异”。有限元是一种数值计算方法,计算过程中的微小数值误差会导致平板在缺乏约束的自由度上发生刚体运动,因此会看到如图2所示的异常结果。

 

  • 解决方案

 

『解决方法1

 

本实例中的模型左右对称,因此圆筒和平板都应该只取一半建模,在整个对称面上定义对称边界条件(即 U1 = 0),这样平板就不会再发生转动或在 x 方向上产生平动。

 

 Abaqus/CAE 里定义对称边界条件的方法是:将边界条件的类型设为 Symmetry / Antisymmetry / Encastre,然后根据模型的具体情况,选择 XSYMM、YSYMM 或 ZSYMM。

 

建立有限元模型之前,应该考虑的第一件事是:该模型是否具备对称性,是否可以只取1/2、1/4甚至1/8进行建模。这样做有多方面的重要意义:

 

1)在对称面上定义对称边界条件,有助于避免刚体位移约束不足问题;

2)可以大大减小模型的自由度和计算规模,缩短计算时间;

3)接触面上的节点减少一半,接触分析更容易收敛;

4)施加了对称边界条件,整个模型的支撑状态变得更加稳固,可能出现的位移状态大大减少,Abaqus / Standard 不必反复尝试不满足对称性的位移解,就更容易找到正确的位移解,会使复杂的非线性分析更容易收敛。

 

有限元分析的一个重要原则是:尽量把所有不会发生位移的节点都固定住,不要让求解器再去通过迭代计算来确定这些节点的位移。

 

需要注意的是:一个模型是否具有对称性,不仅取决于它的几何形状,还要看材料、载荷、边界条件、接触等是否都对称,即变形后的模型是否对称。本实例中,如果平板左侧受到 x 方向的载荷作用,则不允许利用对称性建模。

 

『解决方法2』

 

如果模型不具有对称性,则根据具体情况添加适当的边界条件,以消除不确定的刚体位移。本实例中,可以在平板中央对称线上定义边界条件 U1 = 0。需要注意的是:不能只定义 A 点的 U1 = 0,因为这样整个平板仍然可以绕 A 点做刚体转动。

在上述正确模型的 DAT 和 MSG 文件中,仍然会看到很多警告信息,但都是分析过程中的正常提示信息,并不意味着模型存在错误。本书第20.1.7节“MSG 文件中的正常提示信息”和第20.2.12节“DAT 文件中的正常提示信息”将对这些信息的含义做出详细解释。

 

  • 从零开始学Abaqus到接触非线性分析高阶应用

 

对于动力分析,不需要在所有自由度上定义足够的边界条件,因为动力分析会考虑惯性力,可以避免产生无限大的瞬时运动。如果在动力分析时看到“数值奇异”的警告信息,往往是由于模型中存在其他问题,例如“过度塑性”等。

 

对于这样实用技巧还有很多,想了解更多相关内容欢迎点击在线咨询了解更多。

 

南京复创是达索知名的合作伙伴,专注于为中国制造业数字化和信息化建设提供咨询服务和解决方案,主要涉及汽车、轨道交通、机械、电子、金融、教育、通讯、医疗等各大领域。

关键词:

扫二维码用手机看

Recommended news