有限元分析基础.pptVIP

*第五章等参元或其中(5-1)(5-2)实际单元与基本单元的对应关系可写为用同样的形状函数来插值单元内任意一点(x,y)的位移*第五章等参元(5-3)即：为此单元的结点位移列向量，为形状函数矩阵。这里采用了同样的形状函数式(5-2)，用同样的结点插值函数表示出单元的几何坐标x、y与u、v，这种单元称为等参单元。*第五章等参元5.2平面八结点等参元类似地可以推广到具有更多结点的单元，如图5-2所示(a)直角坐标系与实际单元(b)自然坐标系与基本单元图5-2八结点等参单元该基本单元的位移函数可取为(5-4)*第五章等参元其中在顶角结点与边中点上的形函数分别为(5-6)(5-5)*第五章等参元(5-7)5.3单元刚度矩阵首先给出单元内的应变列向量，对平面问题，应有按坐标变换关系式(5-1)，有*第五章等参元01写成矩阵表达式为：02(5-8)03由式(5-8)可解出04其中05称为坐标变换的雅可比（Jacabian）矩阵，其中06(5-9)*第五章等参元合写成矩阵形式有(5-10)将式(5-3)代入式(5-7)中，则有为应变转换矩阵，按结点分块表示，有*第五章等参元将式(5-9)代入式(5-11)，即可得出此单元的应变转换矩阵，进而求出。单击此处添加小标题单击此处添加小标题同上，单元内的应力可表示为单击此处添加小标题单元刚度矩阵由虚功原理求得，即单击此处添加小标题(5-11)单击此处添加小标题而i=1,2,3,4*(5-12)上述积分在自然坐标系内进行，得刚度矩阵(5-13)一般参数单元的计算都采用数值积分求式(5-13)的近似值，同时，为了减少计算点的数目和便于编写程序，多采用高斯数值积分方法。二维积分法的高斯求积公式为(5-14)第五章等参元*第五章等参元式中，为对应的坐标位置、值，、为权重系数，L、M为沿、方向的积分点数目。5.4ANSYS平面结构计算示例一个长方形面板，如图5-3所示，其高AB=1m，宽BC=1.5m，板厚b=0.04，孔半径R＝0.2m，长方形面板的弹性模量E=210GPa，泊松比μ=0.3，约束条件：在长方形底边约束全部自由度。载荷：BC边施加垂直向下均布载荷F/m。5.4.1问题描述图5-3长方形板结构*5.4.2ANSYS求解操作过程打开Ansys软件，在Ansys环境下做如下操作。（1）选择单元类型运行PreprocessorElementTypeAdd/Edit/Delete，弹出ElementType对话框，如图5-4所示。单击Add，弹出LibraryofElementType窗口，如图5-5所示，选择PLANE82。图5-4单元类型对话框图5-5单元类型库对话框第五章等参元*第五章等参元厚度的。运行PreprocessorRealConstantsAdd/Edit/Delete。弹出如图5-7所示对话框，点击Add，弹出如图5-8所示对话框，点击OK，弹出如图5-9所示对话框。在THK选项栏中设置板厚度为0.04m。设置完毕单击OK按钮完在ElementTypes对话框中单击Options对话框，弹出如图5-6所示对话框，设置K3选项栏为Planestrsw/thk，设置K5选项栏为Nodalstress，设置K6选项栏为Noextraoutput。表示单元是应用于平面应力问题，且单元是有图5-6PLANE82单元类型选项对话框成设置。*（2）设置材料属性运行PreprocessorMaterialPropsMaterialModels，弹出如图5-10所示对话框，双击Isotropic，弹出图5-11所示对话框，在EX选项栏中设置数值2.1e11，在PRXY选项栏中设置数值0.3。设置完毕单击OK按钮。图5-7实常数对话框图5-8选择要设置实常数的单元类型图5-9PLANE82实常数对话框第五章等参元*第五章等参元图5-10选择材料属性对话框图5-11设置材料属性对话框1建立模