[工学]第三章 平面问题有限单元法3 有限单元法与程序设计 教学课件.ppt

第三章 平面问题有限单元法 3. 等参数单元 9) 数值积分： b) 一维Gauss积分： 积分点的选取：n个积分点不是等间距布置，具体的?i由下列条件确定： 其中： 第三章 平面问题有限单元法 3. 等参数单元 9) 数值积分： b) 一维Gauss积分： 有： 近似被积函数（2n-1次多项式）： 其中： y (米) 1.50 1.00 0.50 0.00 -0.50 -1.00 -1.50 CST (Error) -249 (23) -191 (-11.5) -108 (-18.8) -28 (-28) 52 （-37.2） 134 （-45.5） 210 （-62） T6 (Error) -272.7 (-0.7) -180.5 (-1.0) -89.2 (0.0) -0.6 (-0.6) 89.1 （-0.1） 179.6 （0.1） 271.2 （-0.8） Exact -272.0 -179.5 -89.2 0.00 89.2 179.5 272.0 在 x = 0 截面处的 ?x (t/m2) 10 t/m 3m 9m x y T6 单元网格(182DOF) x y o 2a 2b 2. 四结点矩形单元 (x0 ,y0) o ? ? 第三章 平面问题有限单元法 为便于表达引入局部坐标系（?，?） 八、几种常用的平面单元 且有坐标变换式： 可看出：?,?是坐标的无量纲数，在(-1，1)之间取值。 该局部坐标系又称为正则坐标系或自然坐标系。 2. 四结点矩形单元 设位移模式取为： x y o 2a 2b 1 2 3 4 (x0 ,y0) o ? ? 第三章 平面问题有限单元法 代入结点坐标解出待定参数，可以得到用结点位移表示的单元位移模式： 其中： 也可用试凑法直接写出形函数，如对结点1： 2. 四结点矩形单元 第三章 平面问题有限单元法 而： 将1点自然坐标值及N1值带入，得： 所以： 其中： 2. 四结点矩形单元 第三章 平面问题有限单元法 形函数矩阵： 用矩阵表示： 2. 四结点矩形单元 第三章 平面问题有限单元法 应变矩阵： 其中： 非常应变 2. 四结点矩形单元 第三章 平面问题有限单元法 应力矩阵： 其中： 1. 四结点矩形单元 第三章 平面问题有限单元法 单元刚度矩阵： 可得显式： y x L L 例：方板纯弯曲问题 CST 0.3116 0.6045 0.4126 0.8088 0.4683 0.9121 0.4899 0.9574 R4 0.4461 0.9363 0.4797 0.9770 0.4931 0.9894 0.4978 0.9947 Mesh(1/4Plate) ( 1?1 ) ( 2?2 ) ( 4?4 ) ( 8?8 ) Element Euc/?0L ?xc/?0 Euc/?0L ?xc/?0 Euc/?0L ?xc/?0 Euc/?0L ?xc/?0 Exact Euc/?0L=0.5000 ?xc/?0=1.0000 纯弯曲方板在C点的位移和应力 A B D C y x L L A B D C 3. 等参数单元 第三章 平面问题有限单元法 八、几种常用的平面单元 1) 问题的提出： 构造曲线/曲面边界单元的必要性 如何构造曲线/曲面边界单元： a) 在总体坐标系下直接进行曲边单元分析 b) 将已知的自然坐标下的正方形单元和三角形单元通过坐标变换，转换成曲边单元 3. 等参数单元 第三章 平面问题有限单元法 2)等参单元的概念： 如何描述单元几何形状？ 利用自然坐标下的已知单元构造曲边单元需解决的两个问题： 如何描述单元内任一点的物理量？ 3. 等参数单元 第三章 平面问题有限单元法 2)等参单元的概念： 如何描述单元几何形状？ 将已知单元（母单元）经坐标变换成曲边单元（子单元），建立两个单元上点的一一对应关系。 3. 等参数单元 第三章 平面问题有限单元法 2)等参单元的概念： 由于母单元形函数具有在本点值为1，在其它结点值为0的性质，可利用它建立两种单元坐标的对应关系（且为一一对应） 3. 等参数单元 第三章 平面问题有限单元法 2)等参单元的概念： 将母元中的坐标线映射到子元上，则得到一个斜角坐标系或曲线坐标系?? ，称其为局部坐标，原xy坐标系称为整体坐标。 当母元确定后（即结点数和形函数确定），坐标变换式只与子元结点坐标有关，而子元坐标与网格划分有关，故一个母元与一族子元相对应。 试证明这种坐标变换形成的