8 平面分析-等参单元-2幻灯片.pptVIP

* 平面问题有限元分析-等参单元 * * 平面问题有限元分析-等参单元 * * 平面问题有限元分析-等参单元 * 积分点数的选取与被积函数f(?)有关，当f(?)是m次多项式时，取 当f(?)不是多项式时，需要通过一些试算来判断选取适当的n值。积分点数n(2～10)不能取得过大，否则将随着积分点的增多计算量急剧增加 高斯积分法 高斯求积公式的积分点坐标与加权系数表 n ξi Wi 2 ±0.57735 1 3 0 ±0.77460 8/9 5/9 4 高斯积分法 MATHCAD举例 平面问题有限元分析-等参单元 平面问题有限元分析-等参单元 * 平面问题有限元分析-等参单元 * 8 平面问题有限元分析等参单元 曹国华 8.1等参单元刚度（4节点） 8.2等参单元等效节点力（4节点） 8.3矩形单元（8节点） 8.4等参单元（8节点） 8.5高斯积分法 坐标变换函数与位移函数采用相同的形状函数 对比 * * 平面问题有限元分析-等参单元 等参单元刚度（4节点） （i=1，2，…，4） 等参单元刚度（4节点） 根据形函数，几何矩阵为 * * 平面问题有限元分析-等参单元 * 平面问题有限元分析-等参单元 * 已有 若 求z对(?,?)偏导数 应力等分析需要整体坐标下的形函数，怎样求解？ 等参单元刚度（4节点） * 平面问题有限元分析-等参单元 * 等参单元刚度（4节点） 求不出来 在上式中， 是ξ、η的函数，因此必须用坐标变换式来转换导数关系，根据复合函数的求导法则，有 （i=1，2，…，4） 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 式中的J称为 雅可比（Jacobi）矩阵 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 求出 关键 雅可比（Jacobi）矩阵的逆矩阵为 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 可求 将单元应变代入平面问题的物理方程式，就得到平面4节点等参单元的应力列阵，为 （i=1，2，…，4） 等参元刚度矩阵仍可用下式表示 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 面积微元的变换 计算单元刚度矩阵时，需将对整体坐标的积分变成对局部坐标（自然坐标）的积分，其面积微元也由直角坐标中的矩阵dxdy变成自然坐标中的曲边四边形 以下求直角坐标中的曲边面积元dA，用局部坐标中的面积微元 表示的表达式。 0 x y 由图 (a)得 图 (a) 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 由图(b)，在oxy坐标系内任一微矢量 x y o 图(b) 由于x,y是 的函数，所以 故 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 x y o 图(b) 由于 是沿 坐标变化的，所以它只是 的函数，与 无关。可得： 同理 得 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 雅可比行列式 Jacobi 因为 曲边面积元dA: 则等参元刚度矩阵为 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 由以上分析可知，在划分单元时，只需确定单元节点的整体坐标值，而不必画出等参元的具体形状，因为在计算中实际使用的只有单元四个节点在整体坐标系下的位移值。 等参元变换的条件为 ，因此在有限元网格划分时，要特别注意这一点。 由此式可计算出单元的刚度矩阵，此处的积分一般得不到显式形式，因此要采用数值积分法（如Gauss积分法），该积分对任何形状的四边形单元具有相同的积分表达式和相同的积分限。上式积分运算是在母单元（自然坐标）内进行的，但节点力与节点位移的方向均沿整体坐标x，y方向。所以[k]e是整体坐标内的刚度矩阵。 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 单元内某点受到集中载荷P=[Px Py]T，移置到单元节点上的等效节点力为： 等参单元等效节点力（4节点） （1）集中力引起的单元节点载荷 其分量形式： * * 平面问题有限元分析-等参单元 （2）体积力引起的单元节点载荷 单元的单位体积力为P=[X Y]T，移置到单元节点上的等效节点力为： * * 平面问题有限元分析-等参单元 高斯积分 （3）表面力引起的单元节点载荷 0 x y * * 平面问题有限元分析-等参单元 单元的某