FE-Ch02.5矩形单元.pptVIP

结点位移： (i=1,2,3,4) 单元结点位移: 位移模式： 求解 ?： 将?代入位移模式， 可以得到单元位移通过单元结点位移的表达式。 其中 为单元插值函数，令 则， ?，?是矩形单元的自然坐标 自然坐标的原点取在形心上。 又： 且插值函数满足 因此， 又： 显然， 由于ξ0和η0中都包含变量 ξ 和η ， 和 均不是常数矩阵。 高次三角形单元可以借助 pascal三角形(杨辉三角) 系统的求解。 1）二次单元（平面问题：6个结点，12个自由度) 利用6个结点坐标 xi ,yi (i=1,2…6) 求解广义坐标。 2）三次单元（三角形) 10个结点*2=20自由度 xi ,yi (i=1,2…10) 即广义坐标的个数为20 计算量大、繁琐。 1.面积坐标（三角形单元） 定义：三角形内任一点P的面积坐标可表示为： 令： 即： 2.面积坐标的性质 1)与j-m 边平行线上的三角形内点有相同的值Li 2)角点坐标为 3)形心的坐标为 4)每边的方程为： 5) 三个坐标不是独立的，只有两个独立的坐标。 3.面积坐标与直角坐标的转换关系 从而 ? 由于以上两种坐标变换是线性变换： 所以， 面积坐标表示的多项式 ?? 直角坐标中的同阶多项式。 4.面积坐标的微积分计算 1)导数计算（复合函数求导） 同理： 2)积分计算 面积积分： 已知 如，均质等厚单元的自重问题 线积分（三角形的某一条边i-j边上，边长为l） 注：此处不能有 ，?在i-j边上， 例： n=3, 即单元结点数等于3，一次三角形单元，k=2，对于结点1： 例：一次三角形单元，对于角结点1，23边方程为：L1=0 * 设边长分别为2a和2b。 取坐标原点在单元形心上， 轴分别平行两对对边。 单元共有8个自由度，广义坐 标取8。 即，广义坐标： 位移多项式： 位移模式是具有完全一次式，非完全二次式。 因为它在x,y方向呈线形变化，所以称之为双线形位移模式。 §2.5 矩形单元与高精度三角形单元 一. 矩形单元 单元等效结点载荷列阵 二.高精度三角形单元 三.面积坐标为自然坐标的插值函数 j i m y x P ?Pij ?Pim ?Pjm 其中 5.用面积坐标给出单元插值函数（线性单元与三角形高次单元） 线性单元： 插值函数构造点i的单元插值函数一般方法 研究一个每边有k个节点（等间距）的高次单元,单元中结点的总数为: P 为 (0,1,…k-1) 由图可类推出一个数列： 插值函数的次数等于k-1次.对二次单元, 每边有k-1=2个节点和结点总数n=6. 单元插值函数为 同理： 例 对于三角形高次单元, n=6，每边有3个结点，k=3 构造点i的单元插值函数Ni的几何解释: 单元插值函数 类似地 1、角节点--选择若干直线,这些直线应该经过除点i以外的其他点并且平行（或经过）点i的对边. 2、边节点--选择非此边节点所在的其他两直线. 3、将这些直线方程的左端项相乘并除以点I的面积坐标(到相应直线的距离),即可构造点i的单元插值函数Ni。