内容发布更新时间 : 2024/6/7 12:23:23星期一 下面是文章的全部内容请认真阅读。

第10章 非线性动力有限元法 ....................................................................................................... 1

10.1 几何非线性问题的有限元法 .......................................................................................... 2

10.1.1 几何非线性问题的牛顿迭代法 ........................................................................... 2 10.1.2 典型单元的切线刚度矩阵 ................................................................................. 4 10.2 材料非线性问题的有限元法 .......................................................................................... 8

10.2.1 弹/粘塑性问题的基本表达式 .............................................................................. 8 10.2.2 粘塑性应变增量和应力增量 ............................................................................... 9 10.2.3 弹/粘塑性平衡方程 ............................................................................................ 10 10.3 材料非线性问题的动力有限元法 ................................................................................ 11 10.4 应用举例 ........................................................................................................................ 14

10.4.1 粘弹粘塑性动力有限元分析举例 ................................................................... 14 习题......................................................................................................................................... 15

第10章 非线性动力有限元法

当机械结构受到较大的外载荷，或受到持续时间较短的冲击载荷作用时，结构会产生过大的变形, 以至于必须考虑结构几何大变形对结构整体刚度及固有频率的影响，即所谓的几何非线性影响。另外, 对于多数非线性动力学问题，还需要考虑材料非线性、接触非线性等方面的影响。

非线性动力学分析求解的基本方程有如下形式

???I?P?0 (4.141) Mu??Ku为粘性效应项，考虑阻尼、粘塑、粘弹等效应。P为外部激励。 式中，I?Cu对于考虑各种非线性效应的动力学问题求解，需要对动力学方程进行直接时间积分。即非线性动力有限元分析具有如下特点：（1）问题分析过程需要考虑时间积分效应，不必做模态分析，不必提取固有频率；（2）采用直接积分方法求解非线性动力学方程，需要对时间作积分计算，因此计算量远远大于线性模态动力学方法；（3）非线性动力学分析中可以施加不同类型的载荷，包括结点力、非零位移、单元载荷；（4）在每个时间步上，进行质量、阻尼、及刚度的集成，采用完整矩阵，不涉及质量矩阵的近似；（5）可以同时考虑几何、材料和接触等多种非线性效应。

非线性动力有限元分析程序常采用隐式Hilber-Hughes-Taylor法进行时间积分运算。这种方法适于模拟非线性结构的动态问题，对于冲击、地震等激发的结构动态响应以及一些由于塑性或粘性阻尼造成的能量耗散，隐式算法特别有效。隐式积分方法需要对刚度矩阵求逆计算，并通过多次迭代求解增量步平衡方程。隐式Hilber-Hughes-Taylor时间积分算法为无条件稳定，对时间步长没有特别的限制。

采用子空间法也可以对动力学平衡方程作时间积分运算。子空间法是提取模态分析得到的各阶特征模态，并采用与线性模态动力学分析方法相近的分析方式进行求解。对于带有微小非线性效应的问题，如材料小范围进行入屈服、结点转角不大的情况，子空间法效率比进接积分法要高。

此外，非线性动力有限元分析还可以采用显式动态算法，如中心差分法。显式时间积分算法为有条件稳定，其临界稳定时间步长限制了时间步长的大小，与有限元模型最小单元尺寸、材料应力波速等有关。显式时间积分法适于模拟高速冲击、接触等问题。

上述方法的选择需要综合考虑计算量、分析问题的规模、单元限制等多方面因素，需要丰富的有限元模拟的理论、经验和实践知识。以下以几何非线性问题和材料非线性问题为例介绍非线性有限元法，其中粘弹粘塑性非线性材料问题的分析是典型的非线性动力有限元的求解思想。

10.1 几何非线性问题的有限元法

几何非线性问题一般是指物体经历大的刚体位移和转动，但固连于物体坐标系中的应变分量仍假设为小量, 即大位移小应变情况。

10.1.1 几何非线性问题的牛顿迭代法

由数值分析技术可知，求解非线性方程组的数值方法的常规方法是Newton-Raphson法，即牛顿迭代法，这是一种近似线性化迭代求解方法。对于非线性方程?(x)?0，具有一阶导数，在xn点作一阶泰勒级数展开，它在xn点的线性近似为

?(x)??(xn)?(d?)n(x?xn) (4.142) dx因此，非线性方程?(x)?0在xn附近似为线性方程：

?(xn)?(当(d?)n(x?xn)?0 (4.143) dxd?)n? 0时，由上式求得n步的修正项 dx?Xn?1???(xn)/(Newton-Raphson方法的迭代公式为

d?)n (4.144) dxXn?1?xn??Xn?1 (4.145)

在几何非线性有限元法中，结构的刚度矩阵与其几何位置有关，平衡方程由变形后的位形描述，因此，结构的刚度矩阵是几何变形的函数。设变形为δ, 结构的平衡方程式

K(δ)δ?R?0 (4.146)

为一个非线性方程组。记非线性方程

ψ(δ)?Kδ?R?0 (4.147)

用Newton-Raphson方法求ψ(δ)?0的根时，迭代公式分别为

δn?1?δn?δn?1 (4.148)

其中,

?δn?1满足下式

KT?δn?1?R?K(δn)δn (4.149)

式中, KTn称为切线刚度矩阵，表达式为

dψ(δ))n (4.150) dδ 在每一个迭代步中，通过求解切线刚度矩阵KTn，进而用?δn?1进行迭代求解，称为

KTn?(Newton-Raphson方法，又称切线刚度法。

牛顿法的收敛性是好的。但是某些非线性问题中，使用牛顿法迭代时，若KTn出现奇异或病态，则对KT的求逆出现困难。关于这一点也可以采用其它修正办法，如引入阻尼因子。

对于已经建立的有限元方程，设ψ表示内为和外力矢量的总和，有

δ*Tψ??ε*TσdV?δ*TR?0 (4.151)

V式中, R为载荷列阵；δ为虚位移；ε为虚应变

用应变的增量形式dε?Bdδ代入上式，消去δ项，可以得到非线性问题的一般平衡方程式为

***ψ(δ)??BTσdV?R?0 (4.152)

V该式不论位移或应变的大小与否均成立。

在有限变形中，应变和位移之间的关系是非线性的，即B矩阵是δ的非线性函数。但是，近似地可将进行如下分解：

B?B0?BL (4.153)

式中, B0为线性应变分析的部分； BL为由非线性变形引起的，与δ有关。

假定应力应变关系为线弹性，于是有

σ?D(ε?ε0)?σ0 (4.154)

式中 [D]为材料的弹性矩阵； {?0}为初应变列阵；{?0}为初应力列阵

对于式(4.152)的非线性平衡方程式，可用Newton-Raphson方法进行迭代求解。对该式微分，有

dψ??dBTσdV??BTdσdV (4.155)

VV不考虑初应变和初应力的影响，得

dσ?Ddε?DBdδ

并且

dB?dBL