大变形问题的有限单元法讲义.ppt

2000.4 哈尔滨建筑大学 王焕定教授制作 第五章 大变形问题的有限单元法 1. 弹性大变形问题的有限元法 2. 弹性分支点稳定问题有限元分析 3. 物质描述大变形增量问题的T.L 、U.L法 1. 弹性大变形问题的有限元法 2. 弹性分支点稳定问题有限元分析 3.1 物质描述大变形增量问题的T.L法 3.2 物质描述大变形增量问题的U.L法 验证δALθ=ALδθ 推证式（A） 式（5. 96）的推证 t时刻应力引起的等效结点力矩阵 t+Δt时刻荷载引起的等效结点力矩阵 或 将其按集成规则集装后可得 再引入如下记号 将 和 的表达式代入，可得 几何或非线性应变增量刚度矩阵 或 利用这些关系，非线性平衡方程可写为 为求解上述方程，尚需解决如下两方面问题 首先假设 然后将ΔS和ΔE的关系线性化。根据本构关系 则有 为使其线性化，设( t时刻的材料性质矩阵) 在做了上述两方面处理后，可得 由此出发，用非线性方程的相关解法，即可解决大变形非线性（材料非线性）问题。 将其代回非线性平衡方程，可得 讲义上给出了T.L法的求解步骤，可供大家编程序参考。 因t+Δt的位移是用t时刻位形为基准度量的，因此 在[t,t+Δt]间隔内，以t时刻位形为参考位形，其增量位移为 象T.L法一样，设单元和位移的描述为 但需指出的是，式中形函数N是t时刻单元自然坐标的函数。在计算 等导数时，要先作坐标变换(Xi应换为xi)。 To 26 类似地，用矩阵来表示则有 在时刻t和t+Δt的格林应变是以t时刻位形定义的，因而它们分别为 式中算子矩阵象T.L法一样，但应将Xi换为xi。 再次强调，式中算子符号象T.L法一样，但应将Xi换为xi。 基于上述说明，象T.L法一样可导得 关于应力的处理也和T.L法一样，对t 时位形定义的t和t+Δt时刻的克希荷夫应力分别为 象T.L法一样由虚位移原理虚功方程可导得 象T.L法一样推导，引入如下符号定义 式中 其中 分别为t时刻位形定义的单元体积、应力表面、体力和表面力。 几何或非线性应变增量刚度矩阵 荷载引起的等效结点力矩阵 则可得 t时刻应力引起的等效结点力矩阵 t+ Δt时刻的非线性平衡方程 象T.L法一样，为求解上述方程也需解决线性化问题。首先讨论ΔS的计算。因为 因此 可改写为 由第四章已知 Chap4 64 式中 Chap4 44 Chap4 26 可得 由如下两式消去 并利用 ,且注意到Vij对称、Ωij反对称 GoTo 49 又由于 、 和 ，因此 GoTo 50 上式最后一项将使本构张量不对称，对金属类不可压缩介质，这一项可略去，也即 在有限的克希荷夫应力和格林应变增量之间仍认为 这就是U.L法的本构关系线性化。 与T.L法一样，除本构关系线性化外，还需解决几何方面的线性化。因为 为对上述三项积分作几何方面的线性化，可设 这一关系可有两种矩阵表达方式：其一是 其另一方式是 引入讲义上式(5. 102)算子矩阵 ，则有 基于上述讨论，积分 时，作如下处理 式中 为由 组成的6×6矩阵， 矩阵为 结 束 引入如下两个矩阵 则t+ Δt时刻的非线性平衡方程 可改写成 需要指出的是，要提高精度可减小Δt或采用适当的修正技术。 * 弹性大变形问题，需要考虑变形的非线性项和变形对平衡的影响。 若以初始自然平衡状态作初始位形，则物质描述的格林应变为 式中 线性部分 非线性部分 为便于计算机编程，将张量转换为矩阵： 格林应变矩阵和张量的分量间有如下关系 对应的克希荷夫应力矩阵和张量分量间关系为 引入两个算子矩阵 式中 再引入位移梯度向量的记号 3阶单位矩阵 在上述符号基础上，格林应变由位移表为 则单元格林应变为 其中线性应变矩阵B和有限元(I)一样 设单元位移场和有限元(I)一样为 非线性部分“应变矩阵”为 式中G为如下9×3m的矩阵 式中AL为如下6×9的矩阵 由(AL)可见，格林应变-位移关系是非线性的。 非线性部分因为 所以 为用虚位移原理建立单元特性方程，还得建立应变增量和位移增量间的关系。对线性部分 综上所述，格林应变增量为 验证 如果记 ，则