ANSYS非线性基础培训手册Basic5a

1、Training Manual几何非线性基础几何非线性基础Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual大应变理论大应变理论 该附录中所包含的知识对于成功地使用该附录中所包含的知识对于成功地使用 ANSYS 中的几何非线性不是必需的中的几何非线性不是必需的, 因此因此, 通常不在课程通常不在课程中讲解中讲解. 它作为附加的背景知识提供给那些希望更深入地它作为附加的背景知识提供给那些希望更深入地理解

2、理解 ANSYS 大位移特征的用户大位移特征的用户.Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual将非线性应变定义推广到一般的三维情况将非线性应变定义推广到一般的三维情况 在二维和三维中在二维和三维中, 当一个元件经历大应变变形当一个元件经历大应变变形时时, 不仅长度元素改变不仅长度元素改变, 厚度厚度、面积和体积也改面积和体积也改变变.AA0PBasic Structural Nonlinear

3、ities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual运动和变形运动和变形 物体在外载荷的作用下会移动和变形物体在外载荷的作用下会移动和变形. 如果考察该物体上某一点的运动如果考察该物体上某一点的运动, 它的初始位置它的初始位置为为 X , 最终位置为最终位置为 x , 则位移量为则位移量为 u , 图中图中 XY X u x XxuBasic Structural Nonlinearities 11.0Basic Structural Nonlinear

4、ities 11.0Basic Structural Nonlinearities 11.0Training Manual变形梯度变形梯度 变形梯度是物体变形程度的度量变形梯度是物体变形程度的度量, 定义为定义为 变形梯度变形梯度 F 包含如下信息包含如下信息: 体积改变体积改变 转动转动 由于应变而引起的形状改变由于应变而引起的形状改变 XxFBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual 变形

5、梯度变形梯度 注意由定义注意由定义, 变形梯度变形梯度 F 消除了平动消除了平动; 是应变定是应变定义的必要条件义的必要条件. 当定义应变时当定义应变时, 还应排除转动部分还应排除转动部分 (因为它对应变因为它对应变没有贡献没有贡献), 并提取形状改变部分并提取形状改变部分. 这可以通过应这可以通过应用用极分解定理极分解定理 实现实现.Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual极分解极分解 变

6、形梯度变形梯度 F 可以用极分解定理分解为转动部分和变形部可以用极分解定理分解为转动部分和变形部分分: F = R U R U = 转动矩阵转动矩阵, 包含物质点象刚体一样转动的量和方向包含物质点象刚体一样转动的量和方向的信息的信息.= 伸长矩阵伸长矩阵, 包含在物质点处物体的应变信息包含在物质点处物体的应变信息.Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual应变定义中应变定义中 U 的构造的构造

7、 知道了伸长矩阵知道了伸长矩阵 U , 对一维对数应变对一维对数应变 e el 和一维和一维 Green-Lagrange应变应变 e eG 进行推广进行推广, 可以构造出三可以构造出三维一般应变维一般应变. Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training ManualHencky应变应变 对数对数 (Hencky) 应变按下式计算应变按下式计算:式中式中 e eH 是按矩阵形式表示的应变张量是按矩阵形式表示的

8、应变张量. 在此情况下在此情况下, e eH 是一维对数或真实应变是一维对数或真实应变 e el 的三的三维等效维等效. UHlneBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training ManualGreen-Lagrange应变应变 在三维中在三维中, 可以按下式所示直接由伸长矩阵可以按下式所示直接由伸长矩阵 U 计算计算 Green-Lagrange 应变应变: 该度量从应变场估算中直接忽略转动矩阵该度量从应变场

9、估算中直接忽略转动矩阵 R . e eG 可以用位移场的梯度项重写为下式可以用位移场的梯度项重写为下式: 前两项是线性小应变项前两项是线性小应变项, 最后一项是对应变度量最后一项是对应变度量的非线性贡献的非线性贡献. IUUTG21e XuXuXuXuTTGeBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual将非线性应力定义推广到一般三维情况将非线性应力定义推广到一般三维情况 与一维中的情况一样与一维

10、中的情况一样, 在二维和三维中有共轭应在二维和三维中有共轭应力度量力度量, 它可以对每一种非线性应变定义它可以对每一种非线性应变定义.Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training ManualCauchy应力应力 Cauchy 或真实应力张量或真实应力张量 t t (此处写为矩阵形式此处写为矩阵形式) 给出变形构件中单位变形面积的当前力给出变形构件中单位变形面积的当前力. 如果令如果令 那么那么, 在三维中在

11、三维中, Cauchy应力张量应力张量 t t 把把 dP 与与 dA 联系起来联系起来 Cauchy 应力是容易解释的物理量应力是容易解释的物理量.= 定义变形体中单元面积分量的矢量定义变形体中单元面积分量的矢量= 作用的相应单元力作用的相应单元力 dA dP dAdPtBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual 第二第二 Piola-Kirchhoff 应力应力 令令 代表由变换形式代表由

12、变换形式 推导出的推导出的力分量力分量 第二第二 Piola-Kirchhoff 应力张量应力张量 S 把把 和和 联联系起来系起来 S 是对称应力张量是对称应力张量, 常常用于有限应变弹性公式中常常用于有限应变弹性公式中, 是是 Green-Lagrange 应变应变 e eG 的共轭应力张量的共轭应力张量. S 是是个非物理的应力张量个非物理的应力张量 (一个伪应力张量一个伪应力张量). 不能直接对不能直接对S 进行物理解释进行物理解释. Pd dPFPd1= 定义未变形体中单元面积的矢量定义未变形体中单元面积的矢量, 这里这里 0dAdAdAndeformatio0 Pd0dA 0dAS

13、PdBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manualt t 和和 S 的关系的关系 物理的物理的 Cauchy t t 应力可以通过下式直接与非物应力可以通过下式直接与非物理的第二理的第二 Piola-Kirchhof f伪应力伪应力 S 联系起来联系起来 : TFSFFdet1Basic Structural Nonlinearities 11.0Basic Structural Nonlinea

14、rities 11.0Basic Structural Nonlinearities 11.0Training Manual完全一致非线性切向刚度矩阵完全一致非线性切向刚度矩阵 众所周知众所周知, 在大多数非线性结构问题中在大多数非线性结构问题中, 应用一应用一致非线性切向刚度矩阵可以迅速提高基于致非线性切向刚度矩阵可以迅速提高基于 Newton-Raphson 求解过程的收敛速度求解过程的收敛速度. 一致或完全切向刚度矩阵通常在迭代求解过程中一致或完全切向刚度矩阵通常在迭代求解过程中产生二次方的收敛速度产生二次方的收敛速度.Basic Structural Nonlinearities 11

15、.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual何谓一致非线性刚度矩阵何谓一致非线性刚度矩阵? 一致非线性刚度矩阵一致非线性刚度矩阵 通过对离散化的有限通过对离散化的有限元方程求导得到元方程求导得到, 它是单元内力矢量它是单元内力矢量 和单元和单元施加的载荷矢量施加的载荷矢量 的函数的函数. nleK inteF aeFBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Ba

16、sic Structural Nonlinearities 11.0Training Manual离散化的非线性静态有限元方程离散化的非线性静态有限元方程 求解的离散化的非线性静态有限元方程可以在单元层次上求解的离散化的非线性静态有限元方程可以在单元层次上描述为描述为:式中式中 01eNeaeinteTnFFT= 单元总数单元总数= 单元坐标系中的单元内力矢量单元坐标系中的单元内力矢量= 转换矩阵将转换矩阵将 变换到全局坐标系变换到全局坐标系= 全局坐标系中全局坐标系中, 在单元层次上施加的载荷矢量在单元层次上施加的载荷矢量eNinteFaeFnTinteFBasic Structural N

17、onlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual单元内力矢量单元内力矢量 单元内力矢量单元内力矢量 由下式由下式给出给出inteF式中式中按照上面给出的内力定义按照上面给出的内力定义, 离散化的非线性有限元方程离散化的非线性有限元方程 (力平衡力平衡) 可以重写为可以重写为: 01eeNeaeVevTnFdVBT eTinteveVFBdV= 单元应变单元应变 - 节点位移矩阵节点位移矩阵= 单元应力矢量单元应力矢量= 单元体积单元体

18、积vB eVBasic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Basic Structural Nonlinearities 11.0Training Manual推导增量非线性刚度矩阵推导增量非线性刚度矩阵 一致非线性刚度矩阵一致非线性刚度矩阵 通过对离散化的有限通过对离散化的有限元方程求导获得元方程求导获得,如下所示如下所示:nleK aeueeinceaeinteTnnleKKKKFFTuK式中式中 eNlVvNkninceincedVjullkBkiTjiKKeu ,11 eNlVvNkneedVljulkBkiTjiKKeu 11, eNlVvNknueuedVllkBjukiTjiKKeu,11 uFKaeaeComponentsStressofNumberNFreedomofDegreesElementofNumberNuuBasic Structural Nonlinearities 1