ANSYSMechanical的強大非线性分析能力.doc

ANSYS/Mechanical的强大非线性分析能力 ? / 《ANSYS Mechanical ─ A Powerful Nonlinear Simulation Tool》，Grama R. Bhashyam，2002年9月） 1 前言 随着有限元算法理论、计算机硬件和软件技术、实际工业需求等的进步，现代CAE技术的应用已逐步由以线性模拟为主向以非线性模拟为主快速发展。一个好的非线性CAE程序必须在这样一些技术领域有完美的解决方案：(1)有限元单元技术；(2)材料本构模型；(3)接触算法和装配分析；(4)针对大型复杂非线性问题的有效解算算法；(5)良好的应用程序组织结构（易用性和可靠性）；等等。本文概要讲述ANSYS/Mechanical程序的部分主要的和独特的非线性技术，以让大家有一个概略性的了解。 2 ANSYS单元技术：实现有效的分析仿真的基础 早在1994年的5.3版本中，ANSYS就已经具备了一个功能强大、适应面宽的单元库。为了适应复杂多变的材料本构和宽范围应用的需求，ANSYS公司随后集中技术力量开发了一组全新的单元（180系列单元），这些单元具有如下显著特点： ????????? 功能丰富 ????????? 融合了最新的基础理论和最先进的算法 ????????? 柔性的单元构架 传统的全积分等参实体单元应用范围很有限，在线性或非线性分析中，可能会发生严重的“锁定”等问题。作为一个通用的分析工具，ANSYS的有限元单元面对的是一个范围非常宽广的应用范畴。因而，180系列实体单元（182、183、185、186和187）综合采用了如下的一些全新的单元技术来保证其宽范围的应用需求： ????????? 选择性降阶积分：避免低阶单元在分析近似不可压材料时发生体积锁定； ????????? 增强应变列式：避免低阶单元的过刚（弯曲状态）、剪切和体积锁定（近似不可压材料）； ????????? 一致性降阶积分：更有效地避免在分析近似不可压材料时发生体积锁定，避免低阶单元在弯曲状态的体积锁定，采用沙漏控制避免降阶积分的零能模式； ????????? 混合u-P列式：除位移外，将静水压力作为附加自由度，有基于罚函数和基于拉格郎日乘子等两种u-P列式，适用于近似不可压超弹材料、近似不可压弹塑性材料和完全不可压超弹材料等。 180系列梁单元（188和189）具有如下独特特性： ????????? 梁截面的真实描述：任意形状的截面（可用ANSYS二维建模功能建立）均可作为梁截面，程序自动计算所有截面特性参数； ????????? 截面计算的网格数可随意设定； ????????? 梁截面可由任意多种不同材料构成； ????????? 支持端面约束扭转和传统的无约束扭转； ????????? 全三维真实形状可视结果； ????????? 支持所有的弹塑性、超弹性、粘弹性等材料本构模型； 在最近的版本中，还将从所有几何特性上真正支持锥型梁（变截面梁），而不是象以往那样通过变化截面特性参数来近似模拟。 同样，在180系列壳单元（181）中，也融合了一系列独特特性： ????????? 横向剪切基于Bathe-Dvorkin假设，可用一致性降阶积分或全积分（采用非协调模式增强膜特性）； ????????? 支持多种弹塑性、超弹性和粘弹性材料本构模型； ????????? 支持层状复合材料结构（带层间剪切应力恢复）； ????????? 提供专用“壳截面定义”工具来定义壳单元的几何和材料特性：用函数方式定义壳的厚度变化，直接定义各复合材料层的厚度、各向异性方向、积分点数目，等等； 单元技术是有限元分析程序的基础和重中之重，如上提到的和大量未提到的新技术的运用使得ANSYS全新的180系列有限元单元具有超强的性能和功能，诸如：满足宽广范围的应用需求（例如，对于材料非线性，现在在所有新单元中都可以采用任意的材料本构关系而不象以前那样特定的材料本构只能由特定的单元来模拟，同时，也可更加方便地采用具有独特应用的新材料模型）；具有良好的结构适应性（例如：对于必须采用退化单元的复杂结构模型，ANSYS程序自动选择合适的单元形函数和单元积分规则；不同阶次的单元自动满足相容性条件等）；等等。 2 ANSYS的材料非线性 在工程设计和分析中，必须很好地理解和准确地描述材料行为特性，但要做到这一点是非常困难的，Lemaitre和Chaboche曾经用如下很形象的例子描述了材料特性的复杂性： “在室温下的一块钢材，可以认为是： 常规结构分析中的线弹性材料； 振动阻尼问题中的粘弹性材料； 刚体； 在非常大的载荷作用下，是理想塑性材料； 要精确计算永久变形时，它是具有硬化特性的弹塑性材料； 对于应力松弛问题，又有弹性粘塑性特性； 当计算成型极限时，它是延