第七章等参元解说.ppt

第七章 等参元 子单元与母单元的变换 高斯积分 等参元实际运用 教学程序学习阶段 子单元向母单元的变换 上式求逆： 回代后得： 利用上式可以把形函数Ni（?，?）对x, y求导问题化为对?，?求导的问题。 子单元向母单元的变换 为了计算单元刚度矩阵及等效节点载荷，还要把总体坐标下的微 元面积 dA 换为局部坐标上去。 子单元向母单元的变换 等参变换受什么限制？ 对雅克比矩阵进行展开观察 由此可看出，只要（xi, yi）i=1,2,3,4给出，J就确定，且仅是?，? 的线性函数 子单元向母单元的变换 为简化，令 于是 为了求得雅可比变换列阵的逆矩阵J-1及单元刚度矩阵积分式中的微元面积dA,要求变换行列式在整个单元上不等于零,而且不能为负值。 子单元向母单元的变换 具体到单元网格上： 上述限制虽然是由四节点斜四边形单元分析得出来,但其结论可推广到一般等参元. 为了能确保进行等参变换,在总体坐标下所划分的斜四边形单元必须是凸四边形,而不能有一内角大于或等于180度的四边形.如下图: 积分运算的转化 当将一个子单元转换为母单元后，原来对子单元的积分就变成在母单元上进行： 由此可见被积函数的表达式必将非常复杂，进行积分运算需要较大的工作量. 由多元微积分的变量替换规则可知 积分运算的转化 采用高斯积分（Gauss积分）可以解决这一问题 ： 高斯积分的精度？ 在单元内选某些点做为积分点，把这些点的坐标值代入，算出其被积函数值，在乘以加权系数，然后求其总和就得到近似的积分值。 基本步骤 由此必然会带来一定的误差。一般来讲，数值积分的精确度与所采用的求积分公式及所取的积分点数有关，对同一种积分方法，积分点数越多则产生的误差越小，但所需的计算工作量也就越大。两者（精度与计算时间）需做统一安排。 积分运算的转化 高斯积分的精度可以保证 数值积分的精确度与选取积分点的数目和位置以及HK中的插值表达式密切相关，若对积分?K及HK中的插值表达式相关，求积公式为 积分点 坐标 积分点数 加权系数 在选择适当的?K 和HK之后，对任何次数不超过2n-1次的多项式函数f（?）均能精确成立，这就是高斯积分公式， ?K 称为Gauss点 。 高斯积分公式的积分点及加权系数可以根据勒让德多项式的n个不同的实根求得，通常应用时，根据积分点数的n可查表求得。 积分运算的转化 一维高斯积分的高斯点及权函数值表 积分运算的转化 二维高斯积分 首先做? 向的积分而保持?不变 再做向?的积分 积分运算的转化 二维高斯积分点 在实际计算中，为了保证计算的精度而又不过分增加计算的工作量，通常高斯积分点的数目n可根据等参元中的节点个数选取，若有需要，在不同坐标方向上也可选用不同的积分点数 当二维情况积分点数为3时，9个GAUSS点(3?3) 的位置如： * 为什么会有等参元 来源一： 在有限元的网格划分中常用的一些单元，像三角形、矩形、六面体单元等，都是形状很规则的单元。对于形状规则的连续体，用这些单元来离散可以获得比较好的结果。 但是，对于一些几何形状比较复杂的连续体，再用这些单元离散就比较困难了。于是出现了坐标变换的方法来解决这个问题。通过一一对应的坐标变换，把规则的单元转变成形状不规则的单元，就可以用它们来离散几何形状复杂的连续体。 为什么会有等参元 来源二： 有些求解问题的域的几何形状比较规则，那么采用原来的坐标进行积分运算就不是很复杂，但是有些几何形状比较畸形，会使运算处理很麻烦，于是就有这样的想法，能不能去把这个畸形的形状转换到一个比较规则，比较普遍的通用的形体上，在这个形体上去研究它的性质，却不改变原来的问题。 等参元基本概念 例如： 将原来的任意四边形单元（子单元）变换为一个正方形单元（母单元）。子单元来自结构，它代表了真实结构的几何特征，然后通过坐标变换转换到母单元上去，进行一系列有限元运算。 等参元基本概念 在有限元分析中，恰当的选择位移函数是整个方法中最重要的 部分。位移函数的描述通常可用两种方式，其中一种是采用含有若干待定系数?i（也称为广义坐标）的简单多项式，其二是可以采用形函数（插值函数）直接描述，形函数通常是用插值多项式表示， 对于二维位移场来说，可以写成： 式中Ni（x,y）由具体单元和节点确定。 等参元基本概念 对于高阶插值单元来说，所选择的位移函数除了应满足完备和协调条件外，除待定系数的数目还应和单元的节点自由度都一致外（以便由节点位移值唯一的确定全部系数），在此多项式中还应满足几何各向同性条件，即位移函数对于坐标的任何线性变换都能保持不变。 采用高阶单元能够得到较高的精度，得阶数过高会导致计算复杂、计算量大为增加，故一般选用不高于三次的多项式函数。 多项式是位移函数