4.Lagrangian网格研讨.ppt

* TL和UL之间的转换。 * 材料时间导数：如速度是位移的材料时间导数，即当材料坐标固定，对时间求偏导数。 * 初始条件回到材料坐标上。 * 材料时间导数：如速度是位移的材料时间导数，即当材料坐标固定，对时间求偏导数。 * * 杂交元：在单元内部独立假设应力场，在节点处独立假设位移场。 * * 局部减缩积分也是选择性减缩积分 N与x之间不可逆，前面是x与ξ之间可逆 * 共旋列式方法就是从单元总变形位移中扣除刚体平动与刚体转动的部分，在共旋坐标系下基于小变形假设利用现有性能优良的线性小变形单元来计算单元内力和刚度矩阵，再通过变换得到总体坐标系下的单元内力和切线模量。即在保证客观性前提下处理小应变大转动问题。 * * 5 编制程序 数值积分 节点力、质量矩阵和其它单元矩阵的积分不是解析计算的，而是应用数值解答，称为数值积分。最广泛应用的是Gauss积分 式中nQ个积分点的权重wQ 和坐标值ξQ有表可查；见附录3,p545。 指定方程在母单元域上进行积分，其积分区间为[－1，1]。 一个二维单元的Gauss积分为 在非线性分析中，采用积分点数的规则一般基于在线性分析中的相同规则；对于一个规则的单元，积分点数目的选择是能恰好积分内部节点力。单元的规则形式，是指仅通过母单元的拉伸而不是剪切能得到的形式，例如一个矩形二维等参单元。 对于一个4节点四边形单元，如何选择内部节点力的积分点数目？ 完全积分 由于速度是双线性的，单元中的变形率D 和B 矩阵是线性的。如果应力与变形率线性相关，那么它在单元内线性变化。内部节点力的被积函数是近似为二次的，因为它是B矩阵和应力的乘积。在Gauss积分中，对于一个二次函数的精确求解在每一方向上需要两个积分点，所以对于线性材料，需要2×2个积分点得到内部节点力的精确解。 对于线性本构方程的积分，几乎得到内部节点力精确解的积分公式，称为完全积分。 完全积分 4节点四边形单元 母单元节点坐标 形状函数和节点变量 单元与节点速度 单元与节点坐标 B矩阵 内部节点力（a是厚度） 完全积分 5 编制程序 局部减缩积分（选择性减缩积分） 对于完全不可压缩或接近不可压缩的材料，运动必须是等体积的 0 J = 1 式中K是体积模量，? 是剪切模量。在任意的等体积运动中，单元的整个体积将保持常数，即在整个单元中的运动必须是等体积的， 否则，当K是一个非常大的数时(材料接近于不可压缩)，任何非零体积应变将吸收几乎全部的能量。 内部节点力的完全积分可能引起单元的自锁，即出现很小的位移而不收敛或收敛得非常慢。考虑一种线性材料，如果将线性弹性应变能分解为静水和偏量部分，可以写为 5 编制程序 局部减缩积分 为了克服这个困难，最容易的方法是使用局部减缩积分。在局部减缩积分中，压力为不完全积分，而应力矩阵的其余部分为完全积分。为此，将应力张量分解为偏斜部分和静水部分 体积自锁源于单元没有能力准确地表示一个等体积运动。为了消除自锁，必须设计应变场，这样在假设的应变场中整个单元的膨胀为零：为了避免自锁，对于任意保持单元体积的速度场，整个单元的应变场必须是等体积的。对于等体积运动的四边形单元，为了克服沙漏模式，在整个单元中膨胀必须为零。 5 编制程序 局部减缩积分 在局部减缩积分中，压力为不完全积分，以此确定与外力平衡的应力场。如杂交单元（奠基人为卞学鐄先生）。 5 编制程序 局部减缩积分 注意到膨胀部分和偏斜部分是彼此正交的，因此内部虚功率为 将变形率表示为形函数形式，膨胀和偏斜部分的被积函数分别为 局部减缩积分包含在偏斜功率上的完全积分和在膨胀功率上的减缩积分。对于一个4节点四边形单元的局部减缩积分为 关于内力的局部减缩积分表达式为 5 编制程序 局部减缩积分 只有四边形和六面体单元才能采用减缩积分；而所有楔形、四面体和三角形实体单元只能采用完全积分。 关于减缩积分，针对静水部分的单个积分点是在母单元的质心（质量均匀在形心）。偏斜部分是通过完全积分方式积分的。这种方法类似于应用不可压缩材料线性分析的算法。关于其它单元的局部减缩积分算法，可以采取对于线性有限元类似的修正方法建立起来。 单元力和矩阵转换 问题：应用三角坐标建立三角形单元（也称为面积坐标和重心坐标） 解：一个具有线性位移场的3节点单元；单元厚度a。在母单元中，节点以逆时针方向编号。在初始构形中，节点编号必须是逆时针的，否则初始构形和母单元域之间的映射行列式将成为负值。 例4.1：三角形3节点单元 具有线性位移三角形的形状函数是三角坐标： 将空间坐标表示为三角坐标 将空间坐标表示为三角坐标的形式： 这里附加上三角形单元坐标之和等于1的条件。上式的逆矩阵