第6章平面问题高阶单元(简)技巧.pptVIP

第6章 平面问题高阶单元 ; 要想提高计算精度，必须细分网格，增加单元数和点数，因而加大输入数据的工作量。;h;选择位移模式时，第2章提到要考虑解的收敛性，即要考虑到位移模式的完备性和协调性。实际操作中，一般应考虑位移模式的对称性。这是因为，有限元位移模式的选择实际是以帕斯卡（Pascal）三角形基础上的(如图6-2所示)，由低价至高阶，顺序选取，组成多项式。多项式中的项数等于单元节点自由度数。如三节点三角形单元，位移模式取完全一次式，共3项。 六节点三角形单元，位移模 式取完全二次式共6项。如 果某一阶次不能全取，则 应按对称性原则适当选取。; 1 x y x2 xy y2 x3 x2y xy2 y3 x4 x3y x2y2 xy3 y4 ; 在上式表示的位移模式中，a1, a2, a3, a5, a6, a7， a8反映了单元的刚体位移和常 应变。在单元的边界(x=±a 或y =±a)上（或），位移是 按线性分布的。因此，相邻 单元在公共边上的位移是连 续的。这样，位移模式满足 了解答收敛性的充分条件。 ;式中形函数为 ：;则式（6-3）可简写为;（6-7）; 由此可见，[B]是?、 ?的函数，即是x、y的函数。因此单元中的应变不再是常数。 ;其 中 ;在四节点矩形单元中，[k]是一个8×8的矩阵。将[k]写成分块形式： ;（6-17） ; 对本问题给定的位移函数，若体积力是重力的情形（设重度为?），单元等价节点载荷列阵为：; 在三角形单元i, j, m的各边中点增设一个节点，使每个单元具有6个节点， 得到图6-4所示的六节点 三角形单元。 这种单元具有12个 自由度，可以采用完全 二次多项式的位移模式： ;（6-20）;2、面积坐标 ;其中A为三角形ijm的面积，Ai, Aj, Am分别为三角形的Pjm, Pmi, Pijd的面积。 这三个比值Li, Lj, Lm称为P点的面积坐标。 ;（6-23）; 将式（6-24）的3个式子分别乘以xi, xj, xm，然后相加，并利用关系式（6-23a），有;（6-26）;在三角形单元上进行积分时，有;?; 再考虑三角形的边中点，例如i?点，面积划分如图6-7所示。显然有： ;形函数确定后，单元中任意一点的位移可以表示为：; 式中应变矩阵[B]为：;单元中的应力仍可表示为： ;（5）等价节点力向量;单元表面力引起的结点力计算公式仍为： ;如在ij边上受到x方向的三角形分布面力，其集度在i点为ps，在j点为0。对应的节点力向量为（图3-9）