7 平面分析-等参单元-1幻灯片.pptVIP

式（Δ）是形状函数，与位移函数一样，在i节点，Ni=1,在其它节点，Ni=0，该形状函数与位移函数一样求解。 为简便起见，将（Δ）式写成如下通用公式 7.4等参单元-任意直边四边形单元（4节点） * * 平面问题有限元分析-等参单元 位移函数： 7.4等参单元-任意直边四边形单元（4节点） * * 平面问题有限元分析-等参单元 坐标变换函数与位移函数采用相同的形状函数 7.4等参单元-任意直边四边形单元（4节点） 对比 节点位移个数与节点坐标个数一样，在选取单元的位移函数时，可取单元的自然坐标作为自变量，取图形（坐标）变换式的形状函数为位移函数的形状函数。 * * 平面问题有限元分析-等参单元 像上述这样，坐标变换与位移函数采用相同的节点，且取相同的插值函数（形状函数）的变换叫做等参变换。这种单元叫等参元。 在等参元的坐标变换中，局部坐标（自然坐标）中的正方形或立方体称为母单元，而整体坐标内曲边形称子单元。 7.4等参单元-任意直边四边形单元（4节点） * * 平面问题有限元分析-等参单元 MATHCAD例子 每个班由负责人压缩成一个压缩包后发到 caoguohua@cumt.edu.cn 每个人提交的文件名为： 班级-姓名-学号（序号） 比如： 机自08（1）班-某某-***（1） 截止时间： 2010年12月5日 作业： 采用MATHCAD求解 如图所示，该模型中有两个四边形单元，弹性模量为 ＝210GPa，厚度为 ＝0.025m，泊松比 ＝0.3， ＝1kN，求单元所受应力。 * * 平面问题有限元分析-等参单元 * * 平面问题有限元分析-等参单元 平面问题有限元分析-等参单元 平面问题有限元分析-等参单元 * 平面问题有限元分析-等参单元 * 7 平面问题有限元分析等参单元 曹国华 7.1四节点矩形单元位移 函数 7.2四节点矩形单元应变与应力矩阵 7.3四节点矩形单元刚度矩阵 7.4等参单元 虽然三角形单元具有很好的“适应性”，几乎任何复杂边界的弹性体总可以划分为三角形，并且三角形单元计算公式简单，但精度较低。 7.1四节点矩形单元位移函数 三角形单元间虽然能够保证位移连续，但应力的精度较差，不能很好的反映弹性体内应力的准确分布规律。为了提高计算精度，准确反映弹性体内的应力状态，可以采用一些较精密的单元类型。 本节将介绍常用的矩形单元，它采用了比常应变三角形单元次数更高的位移模式，因而可以更好地反映弹性体中的位移状态和应力状态。另外，对一些边界比较规则且呈直线的平面结构的分析，采用矩形单元较合适。这时单元总数可以减少，相应的原始数据准备工作和单元特征计算工作均可节省。 * * 平面问题有限元分析-等参单元 如图所示的矩形单元，不失一般性，令矩形单元的长、宽分别为2a、2b。矩形单元有4个节点，共8个自由度，即共有8个节点位移，采用类似三角形单元的分析方法，同样可以完成对矩形单元的力学特性分析。 7.1四节点矩形单元位移函数 * * 平面问题有限元分析-等参单元 这里引入一个局部坐标系?、?，这样可以推出比较简洁的结果。如图所示，取矩形单元的形心o为局部坐标系的原点，?和?轴分别与整体坐标轴x和y平行，两坐标系存在有以下的坐标变换关系 式中： 、 ——矩形形心处坐标。 矩形形心处坐标以及矩形长、宽可由下式计算 6.1四节点矩形单元位移函数 * * 平面问题有限元分析-等参单元 7.1四节点矩形单元位移函数 3(1,1) 2(1,-1) 1(-1,-1) 4(-1,1) 在局部坐标系中，节点i的坐标是 ，其值分别为±1。如节点1在局部坐标系下的坐标为（－1，－1）。 * * 平面问题有限元分析-等参单元 由于矩形有4个节点，共8个自由度，可以选择有8个待定参数的位移模式，如下 　　该函数称为双线性函数。将节点的局部坐标值代入上式，可列出四个节点处的位移分量，即两组四元联立方程，由此可求得位移模式中的8个未知参数?1，?2，…，?8 7.1四节点矩形单元位移函数 * * 平面问题有限元分析-等参单元 * * 平面问题有限元分析-等参单元 3(1,1) 2(1,-1) 1(-1,-1) 4(-1,1) 7.1四节点矩形单元位移函数 求出α1, α2, α3, α4；α 5, α 6 , α7 , α8 7.1四节点矩形单元位移函数 * * 平面问题有限元分析-等参单元 式中： ——矩形单元的形函数，i＝1，2，3，4； ——形函数矩阵； ——单元节点位移列阵， ，i＝1，2，3，4。 7.1