第三章等参单元ppt课件.ppt

《第三章等参单元ppt课件.ppt》由会员分享，可在线阅读，更多相关《第三章等参单元ppt课件.ppt（59页珍藏版）》请在三一办公上搜索。

1、船舶结构有限元分析Finite Element Analysis of Ship Structure,上海海事大学商船学院江国和,第三章 等参单元,3.1等参单元的引入,采用等参变换的单元称之为等参单元，即单元几何形状的变换和单元内的场函数，采用相同数目的节点参数，及相同的插值函数进行变换。基本思想：现在具有规则形状的单元（区域）上构造位移插值函数，然后把这个具有规则形状的坐标变换映射为物理平面上的一个形状比较复杂的单元，因而，等参单元也被称为映射单元。,3.2四节点矩形双线性单元3.21位移插值函数,对于某些有规则边界的区域，可以把所讨论的区域划分成若干个矩形单元，设一个中心位于坐标原点，变

2、长为2的正方形，如图 取正方形的四个角点作为单元节点，单元的位移模式为,式中所示的位移插值函数的特点：固定x，它是y的线性函数；固定y，它是x的线性函数。因而，我们称它为双线性插值函数，把采用这种位移插值函数的四节点矩形单元称为矩形双线性单元。 矩形双线性单元有着常应变三角形单元所不具有的反映单元弯曲变形的能力，从而计算精度有望提高。 由于上式满足常应变准则即完备性条件，又满足位移连续性条件即协调性条件。因而，矩形双线性单元是完备的协调单元。,3.2.2形函数的确定,建立以节点位移为基本未知量的有限元方程，把u和v表示成形函数与节点位移乘积之和的形式 形函数确定的条件： （1）每个形函数都仍然

3、是坐标的双线性函数； （2）Ni在结点i处等于1，在其他节点处等于0。,根据条件2，可得根据矩形单元各个节点的坐标，可将形函数Ni统一表示为 （i=1,2,3,4）对任意双线性函数f(x , y),都有 为f(x , y)在节点i的值，分别取f(x , y)=1， f(x , y)=x和f(x , y)=y，得：,形函数是定义在给定单元上的，所以上式对单元的每一点（x , y）都成立。单元内一点的坐标通过单元节点坐标 的表达式，即坐标插值函数。 位移插值函数和坐标插值函数具有完全相同的构造，都是以节点处的函数值作为参数，并具有相同的形函数。因而矩形双线性单元也是一种等参数单元。,3.2.2矩形

4、双线性单元的完备性和协调性,完备性：所采用的位移插值函数的确能反映任意给定的刚体位移和常应变，简称为满足常应变准则或完备性条件上式是任何等参数单元满足常应变准则的充分必要条件,协调性：只要相邻单元在其公共边界上具有共同节点，相邻单元间公共边界上任一点的位移就必定相同，所以位移连续条件（即协调性条件）得到满足。在x=1边界上，N1=0，N4=0，有同理，在x=1的边界有 。可知，矩形双线性单元的边界在手里变形后，仍保持为直线边界。此协调性论证适用于任何等参单元。,3.3四节点四边形等参单元,前面介绍的四节点矩形双线性单元只适用于规则区域。本届介绍可用于不规则区域的四节点任意四边形单元。3.3.1

5、四节点任意四边形单元对任意四边形单元，不能再直接采用双线性插值函数。通过适当的坐标变换，把物理平面（xy平面）四边形单元变换成平面上以坐标原点为中心的边长为2的正方形单元，xy平面上的节点1，2，3，4分别对应正方形单元的节点1，2，3，4。如图所示,对于a所示的四边形单元，将各边的等分点用之线连接，并规定它与局部坐标系下单元的相应对边、等分点的连线对应，这样，就建立了图a所示四边形单元和图b所示正方形单元个点间的一一对应关系。,3.3.2自然坐标系下的位移插值函数,自然坐标系下的位移插值函数：式中形函数 （i=1,2,3,4）有可知，插值函数在自然坐标系下满足常应变准则和协调性条件,3.3.

6、3坐标变换,整体坐标和自然坐标之间的坐标变换关系式：由形函数在节点i为1，在其他节点为0的性质可知，上式把正方形的四个顶点变换成整体坐标下的节点（Xi，Yi）。其次，考察母单元某一条边，如12边=-1,上式是一条以为参数的直线方程，因而变换式把母单元的12边变换成整体坐标下连接节点的直边。对其他边也有同样结论。位移插值公式和坐标变换公式具有完全相同的构造，因此这种四节点任意四边形单元也是等参单元，简称为四节点等参数单元，在整体坐标下也必然满足常应变准则。位移插值函数在自然坐标系下的协调性自动保证了坐标变换式的合理性及位移插值函数在整体坐标系下的协调性今后，对每个单元，我们把立足点放在自然坐标上

7、，一切计算都转换到自然坐标系下进行,位移插值函数可用矩阵记号记为：3.3.4应变矩阵B建立对平面应力问题,上式B是一个三行8列矩阵 （i=1，2，3，4）利用复合函数求导法则，得写成矩阵形式，,式中坐标变换矩阵（雅可比矩阵）由矩阵求逆法则得式中,3.3.5应力应变关系,S=DB为应力矩阵3.3.6单元刚度矩阵及等效节点载荷1.面元变换公示自然坐标系下的无限小面积元，经过坐标变换式变换为整体坐标系下无限小平行四边形，如图,变换后微元面积记为dA，则2.单元刚度矩阵把积分区域变换到自然坐标系平面上的相应区域,A单元所占区域t单元厚度,3.等效节点载荷,讨论作用在单元边界上的表面里引起的等效节点载荷

8、,在12边上，=-1，N3=0,N4=0,3.3.7等参变换的条件,为了确保坐标变换式能确定整体坐标与自然坐标间一一对应关系，要求：所以，为保证等参数变换得以进行，所谓的任意四边形所有内角都必须小于180当内角接近180，丨J丨接近与零，因为丨J丨出现在求逆矩阵的分母，将导致计算误差迅速增大，因而计算实践不应使四边形单元过于歪斜。,3.4等参单元的收敛性,为了保证协调性，相邻单元在这些公共边（或面）上应有完全相同的节点，同时每一单元沿这些边(或面)的坐标和未知函数应采用相同的插值函数加以确定。,关于单元的安全性，对于 型单元，要求插值函数中包含完全的线性项，本章讨论的所有单元在自然坐标中是满足

9、此要求的。研究经过等参变换后，笛卡尔坐标下有：将上式代入，就得到单元内的函数表达式可以看出，如果插值函数满足 ，说明单元满足完备性要求,某些情况下，坐标插值多项式次数低于位移插值多项式，计算更为便捷，这样的单元称为次参数单元，如图为典型的二维次参数单元。坐标插值函数的阶数高于位移插值函数的阶数的单元，称为超参数单元，可以更准确的描述壳体的复杂几何形状,3.5数值积分方法,3.5.1插值求积法在求积区间-1，1上取n个分点i（i=1,2n）,n 个求积节点的分布可等距也可不等距，记根据节点上的函数值Fi，构造一个n-1次插值多项式Pn(),使得 ，称为F()的插值多项式或拉格朗日插值多项式。写成

10、如下形式：根据Ni()的性质可知,记则有插值多项式可写为于是有 若记 就得到插值求积公式,式中的系数Hi(i=1,2,3n)称为求积系数或加权系数。它的值与被积函数无关，只依赖于求积节点的个数和位置。当F()本身是次数不超过n-1的多项式时， 即具有n个求积节点的插值求积公式至少具有n-1次的代数精确度。3.5.2高斯求积法对求积节点的位置予以优选，使插值求积公式对任何次数不超过2n-1次的多项式函数都精确成立。这时求积公式的最高代数精确度为2n-1次，称为高斯求积公式,设F()为次数不超过2n-1次的任意多项式，因此有只需要选取i(i=1,2,，n)，使以下n个等式成立，则插值求积公式即为高

11、斯求积公式。积分点的个数n，称为数值积分的阶,N次勒让德多项式在区间-1，1上的n个根就是高斯求积公式中的n个求积节点。利用此性质可确定节点i和系数Hi。因此高斯求积公式又叫高斯勒让德求积法。 n=1情形，得一阶高斯求积公式式中1=0，H1=2，此为中点矩形插值公式，对任意一次多项式都精确成立 n=2情形，二阶高斯求积公式式中 ,它对任意三次多项式函数都精确成立。,下表列出n=1至n=4时高斯求积公式中的求积节点坐标和相应的求积系数,3.5.3二维和三维高斯求积公式因为故二维高斯求积公式为：同理，三维高斯求积公式为：对有限元分析而言，由于几何各向同性要求，个方向的求积节点个数取为相同，故有上述

12、公式,3.6数值积分阶次的选择,选择积分阶次的原则：1.保证积分的精度缩减积分：高斯积分阶数低于被积函数所有项次精确积分所需要阶数的积分方案。采用缩减积分可以取得较完全精确积分更好的精度。理由如下： .精确积分常常是由插值函数中非完全项的最高次方要求，而决定有限元精度的是完全多项式的次方。 .有限元计算模型具有较实际结构偏大的整体刚度，缩减积分方案将使模型刚度有所降低，因此提高计算精度,3.2坐标变换与平面四结点等参元,图3-1(a)为一个任意四边形单元，称为实际单元。在实际单元内以对边的中点连线建立起一个局部坐标系，通过坐标转换把实际单元“映射”为如图3-1(b)所示的一个正方形，此坐标系称

13、为单元的自然坐标系或等参数坐标系，正方形称为基本单元，基本单元内任一点P( ， )与实际单元内的一点P(x，y)唯一对应。,图3-1(a),(a) 直角坐标系与实际单元 (b) 自然坐标系与基本单元,图3-1(b),实际单元与基本单元的对应关系可写为,(3-1),(3-2),用同样的形状函数来插值单元内任意一点(x, y)的位移,(3-3),为此单元的结点位移列向量 为形状函数矩阵。这里采用了同样的形状函数式(3-2)，用同样的结点插值函数表示出单元的几何坐标x、y与u、v，这种单元称为等参单元。,类似地可以推广到具有更多结点的单元，如图3-2所示,(a) 直角坐标系与实际单元 (b)自然坐标

14、系与基本单元,图3-2 八结点等参单元,该基本单元的位移函数可取为,(3-4),其中在顶角结点与边中点上的形函数分别为,(3-5),(3-6),3.3 单元刚度矩阵,首先给出单元内的应变列向量，对平面问题应有,(3-7),按坐标变换关系式(3-1)，有,写成矩阵表达式为：,(3-8),由式(3-8)可解出,其中,(3-9),J称为坐标变换的雅可比（Jacabian）矩阵，其中,合写成矩阵形式有,(3-10),将式(3-3)代入式(3-7)中，则有,为应变转换矩阵，按结点分块表示，有,而 i=1,2,3,4,(5-11),将式(3-9)代入式(3-11)，即可得出此单元的应变转换矩阵 ,进而求出

15、 。,同上，单元内的应力可表示为,(3-12),上述积分在自然坐标系内进行，得,刚度矩阵,(3-13),单元刚度矩阵由虚功原理求得，即,式中， 为对应的坐标位置 、 值， 、 为权重系数，L、M为沿 、 方向的积分点数目。,一般参数单元的计算都采用数值积分求式(3-13)的近似值，同时，为了减少计算点的数目和便于编写程序，多采用高斯数值积分方法。二维积分法的高斯求积公式为,(3-14),3.4 ANSYS平面结构计算示例,3.4.1 问题描述,一个长方形面板，如图5-3所示，其高AB=1m，宽BC=1.5m，板厚b=0.04，孔半径R0.2m，长方形面板的弹性模量E=210GPa，泊松比=0.

16、3，约束条件：在长方形底边约束全部自由度。载荷：BC边施加垂直向下均布载荷F10000000N/m。,图3-3 长方形板结构,3.4.2 ANSYS求解操作过程,打开Ansys软件，在Ansys环境下做如下操作。 （1）选择单元类型运行PreprocessorElement TypeAdd/Edit/Delete，弹出Element Type对话框，如图3-4所示。单击Add，弹出Library of Element Type窗口，如图3-5所示，选择PLANE82。,图3-4 单元类型对话框,图3-5单元类型库对话框,在Element Types对话框中单击Options对话框，弹出如图3-

17、6所示对话框，设置K3选项栏为Plane strsw/thk，设置K5选项栏为Nodal stress，设置K6选项栏为No extra output。表示单元是应用于平面应力问题，且单元是有厚度的。,图3-6 PLANE82 单元类型选项对话框,运行PreprocessorReal Constants Add/Edit/Delete。弹出如图3-7所示对话框，点击Add，弹出如图3-8所示对话框，点击OK，弹出如图3-9所示对话框。在THK选项栏中设置板厚度为0.04m。设置完毕单击OK按钮完成设置。,图3-10选择材料属性对话框 图3-11 设置材料属性对话框,（2）设置材料属性运行Pre

18、processorMaterial Props Material Models，弹出如图3-10所示对话框，双击Isotropic，弹出图3-11所示对话框，在EX选项栏中设置数值2.1e11，在PRXY选项栏中设置数值0.3。设置完毕单击OK按钮。,（3）建立模型 选择Preprocessor Modeling Create Area RectangleBy 2 Corners；弹出如图3-12所示对话框，设置参数，WP X选项栏中填写0，WP Y选项栏中填写0，Width选项栏中填写1.5，Height选项栏中填写1，单击OK，设置完毕。继续运行Preprocessor Modeling

19、Create Area CircleSolid Circle；得到如图3-13所示对话框，在WP X选项栏中填写0.75，WP Y选项栏中填写0.5，在Radius选项栏中填写0.2，设置完毕点击OK按钮。得到如图3-14所示图形。,图5-14 长方形板模型,（4）划分网格运行MeshingSize Cntrls ManualSizeAreasAll Areas，弹出如图3-15所示对话框，在SIZE选项栏中填写0.05，点击OK按钮；运行MeshAreasFree划分网格，网格划分如图3-16所示。,图5-15 设置网格尺寸对话框,图5-16划分网格后板的有限元模型,图5-17 对线施加全约

20、束,（5）施加约束 选择菜单SolutionDefine LoadsApplyStructureDisplacementOn Lines，选择长方形底边，弹出图5-17所示对话框，选择All DOF，单击OK。,（6）施加载荷 选择菜单SolutionDefine LoadsApply StructurePressureOn Lines，弹出如图3-18所示对话框。拾取长方形上边，单击OK按钮。弹出如图3-19所示对话框。在VALUE选项栏中填写10000000。设置完毕点击OK完成设置。,图3-18 选取要施加载荷的边,图3-19 施加载荷对话框,（7）求解 选择SolutionSolveC

21、urrent LS，弹出如图3-20所示对话框。单击OK按钮，开始计算，计算结束会弹出计算完毕对话框，单击Close。关闭对话框计算完毕。,图3-20 求解对话框,（8）后处理 运行General PostprocPlot Results Contour PlotNodal Solu。弹出如图3-21所示对话框，运行DOF SolutionDisplacement vector sum和Stressvon Mises stress，分别显示长方形面板的位移云图和应力云图。,图3-21 云图显示对话框,结果显示如图3-22和图3-23所示。,图3-22 位移变形云图,图3-23 等效应力云图,从图3-22长方形面板的位移云图可知，最大位移发生在圆孔的上部，最大位移为0.75710-4m。 从图3-23长方形面板的应力云图可知，最大应力发生在圆孔的两侧，最大应力为34.6MPa。,