变分与弹塑性力学PPT6.pptVIP

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