应用有限元分析工程实例学习教案.pptx

《应用有限元分析工程实例学习教案.pptx》由会员分享，可在线阅读，更多相关《应用有限元分析工程实例学习教案.pptx（54页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1应用有限元分析应用有限元分析(fnx)工程实例工程实例第一页，共54页。（3 3）二维等参元的位移函数）二维等参元的位移函数*坐标转换坐标转换两个单元的两个单元的4 4个节点和个节点和4 4条边一一对应。条边一一对应。采用坐标变换式，把实际单元映射为采用坐标变换式，把实际单元映射为一个边长为一个边长为2 2的正方形单元。的正方形单元。为形函数矩阵为形函数矩阵(j zhn)(j zhn)用实际单元的节点坐标用实际单元的节点坐标 ，通过插值方法来表示实际单元内任意一点的坐标（，通过插值方法来表示实际单元内任意一点的坐标（x,yx,y）。）。因为形函数是用局部坐标定义的，式因为形函数是用局部

2、坐标定义的，式(4.1)(4.1)建立了正方形单元内任建立了正方形单元内任意点意点与实际单元内任意点（与实际单元内任意点（x,yx,y）之间的意义对应关系。）之间的意义对应关系。（4.24.2）中形函数构造规律：直线方程）中形函数构造规律：直线方程规则(guz)单元-局部系实际(shj)单元-总体系第1页/共54页第二页，共54页。*单元位移函数单元位移函数 为形函数矩阵为形函数矩阵 为单元节点位移列阵为单元节点位移列阵 *等参元定义等参元定义(dngy)(dngy)（按物理场和几何描述）（按物理场和几何描述）等参元：坐标变换和单元位移函数采用相同的节点数，且采用相同等参元：坐标变换和单元位移

3、函数采用相同的节点数，且采用相同 的形函数的形函数 超参元：坐标变换采用相同的节点数多于单元位移函数的节点数超参元：坐标变换采用相同的节点数多于单元位移函数的节点数 亚参元：坐标变换采用相同的节点数少于单元位移函数的节点数亚参元：坐标变换采用相同的节点数少于单元位移函数的节点数 第2页/共54页第三页，共54页。（4 4）单元应变）单元应变(yngbin)(yngbin)、单元应力、单元应力 *单元应变单元应变(yngbin)(yngbin)*单元应力单元应力 （5 5）程序实现）程序实现 平面应力单元、平面应变平面应力单元、平面应变(yngbin)(yngbin)单元、平面轴对称单元、平面轴

4、对称单元可以采用同一类单元计算公式，并且在有限元程序中，单元可以采用同一类单元计算公式，并且在有限元程序中，这三种问题一般采用的是同一个程序路径和同一组子程序。这三种问题一般采用的是同一个程序路径和同一组子程序。第3页/共54页第四页，共54页。（6 6）单元刚度）单元刚度(nn d)d)矩阵矩阵 分别表示单元的体积、面积、厚度。按复合函数链式求分别表示单元的体积、面积、厚度。按复合函数链式求导规则导规则雅可比矩阵雅可比矩阵式（式（4.124.12）两边左乘）两边左乘 ，得，得 被积函数是局部坐标被积函数是局部坐标 的函数，积分域是边长为的函数，积分域是边长为2 2的正方形，的正方形，积分限是

5、积分限是-1-1到到1 1，即积分域规则，积分限简单。，即积分域规则，积分限简单。这里刚度这里刚度(nn d)d)矩阵是矩阵是 8x8 8x8 的矩阵，其计算一般采用高的矩阵，其计算一般采用高斯数值积分。斯数值积分。高斯(o s)数值积分第4页/共54页第五页，共54页。（7 7）二维等参单元的等效节点载荷二维等参单元的等效节点载荷 计算原理和方法与三角形节相同，区别主要在于积计算原理和方法与三角形节相同，区别主要在于积分是在局部坐标系下进行的。分是在局部坐标系下进行的。*体积力体积力:体积力等效的节点载荷公式：体积力等效的节点载荷公式：A A是单元面积是单元面积(min j)(min j)，

6、t t 是单元厚度。是单元厚度。体积力列矢量：体积力列矢量：因为因为 ，所以，积，所以，积分域由实际单元的面积分域由实际单元的面积(min j)A(min j)A，变为局部坐标系下的，变为局部坐标系下的边长为边长为2 2的规则正方形，并且被积分函数是局部坐标的规则正方形，并且被积分函数是局部坐标 的函数。的函数。第5页/共54页第六页，共54页。*体积力体积力:体积力等效的节点载荷公式体积力等效的节点载荷公式(gngsh)(gngsh)：A A是单元面积，是单元面积，t t是单元厚度。是单元厚度。体积力列矢量：体积力列矢量：被积分函数是局部坐标被积分函数是局部坐标 的函数。的函数。（7 7）二

7、维等参单元的等效节点二维等参单元的等效节点(ji din)(ji din)载荷载荷 计算原理和方法与三角形节相同，区别主要在于积分计算原理和方法与三角形节相同，区别主要在于积分是在局部坐标系下进行的。是在局部坐标系下进行的。第6页/共54页第七页，共54页。*表面力表面力表面力表面力等效节点载荷计算公式等效节点载荷计算公式等效节点载荷计算公式等效节点载荷计算公式 表示载荷作用表示载荷作用表示载荷作用表示载荷作用(zuyng)(zuyng)的边界的边界的边界的边界载荷作用载荷作用载荷作用载荷作用(zuyng)(zuyng)边的法矢量（假作用边的法矢量（假作用边的法矢量（假作用边的法矢量（假作用(

8、zuyng)(zuyng)在边在边在边在边1212上）：上）：上）：上）：表面力表达式表面力表达式表面力表达式表面力表达式微弧长：微弧长：微弧长：微弧长：为为为为8x18x1的列阵，但是，只有与载荷的列阵，但是，只有与载荷的列阵，但是，只有与载荷的列阵，但是，只有与载荷所在边的节点对应的行才有等效节点载荷。所在边的节点对应的行才有等效节点载荷。所在边的节点对应的行才有等效节点载荷。所在边的节点对应的行才有等效节点载荷。第7页/共54页第八页，共54页。例例例例4.1 4.1 计算图所示单元的表面力的等效节点载荷（假定单元厚度为计算图所示单元的表面力的等效节点载荷（假定单元厚度为计算图所示单元的

9、表面力的等效节点载荷（假定单元厚度为计算图所示单元的表面力的等效节点载荷（假定单元厚度为 ）表面力的等效节点载荷可按式（表面力的等效节点载荷可按式（表面力的等效节点载荷可按式（表面力的等效节点载荷可按式（4.254.25）计算）计算）计算）计算因为因为因为因为(yn wi)(yn wi)积分只沿积分只沿积分只沿积分只沿 的边进行的边进行的边进行的边进行在该边上在该边上在该边上在该边上 ，,所以所以所以所以平均分配到节点平均分配到节点平均分配到节点平均分配到节点1 1、2 2的负的负的负的负y y方向上。方向上。方向上。方向上。在其他边上亦可类似方法处理。在其他边上亦可类似方法处理。在其他边上亦

10、可类似方法处理。在其他边上亦可类似方法处理。沿边界变化的面力沿边界变化的面力沿边界变化的面力沿边界变化的面力 第8页/共54页第九页，共54页。计算公式计算公式251251边对应于边对应于 的边的边 ，积分只沿积分只沿 的边进行的边进行(jnxng)(jnxng)，只有只有 不等于零不等于零*单元等效节点载荷（单元等效节点载荷（5-85-8节点等参单元）节点等参单元）计算方法与计算方法与4 4节点等参单元类似。节点等参单元类似。举例说明如何计算表面举例说明如何计算表面(biomin)(biomin)力等效节点载荷力等效节点载荷例例4.2 4.2 如图所示的单元，如图所示的单元，q=q=常数。常

11、数。（假定单元厚度为（假定单元厚度为 t t）。）。计算结果表明：载荷分配到节点(ji din)上的比例 为1:4:1第9页/共54页第十页，共54页。均布载荷向节点等效的两种分配方法对比按公式(gngsh)计算时，各节点的分配的载荷比例，凭直观想象是不易理解的。数值计算表明，如果单元网格较密，两种方法对计算结果影响不大。第10页/共54页第十一页，共54页。平面单元算例平面单元算例平面单元算例平面单元算例 实际问题常简化为平面问题，对于一个平面问题，可以采用不同实际问题常简化为平面问题，对于一个平面问题，可以采用不同实际问题常简化为平面问题，对于一个平面问题，可以采用不同实际问题常简化为平面

12、问题，对于一个平面问题，可以采用不同的平面单元进行有限元分析，可以采用不同节点数的三角形单元、四的平面单元进行有限元分析，可以采用不同节点数的三角形单元、四的平面单元进行有限元分析，可以采用不同节点数的三角形单元、四的平面单元进行有限元分析，可以采用不同节点数的三角形单元、四边形单元，也可以将三角形单元和四边形单元组合使用。边形单元，也可以将三角形单元和四边形单元组合使用。边形单元，也可以将三角形单元和四边形单元组合使用。边形单元，也可以将三角形单元和四边形单元组合使用。在力学在力学在力学在力学(l xu)(l xu)(l xu)(l xu)模型确定之后，选择什么样的单元和如何划分网格模型确定

13、之后，选择什么样的单元和如何划分网格模型确定之后，选择什么样的单元和如何划分网格模型确定之后，选择什么样的单元和如何划分网格是必须考虑的问题。是必须考虑的问题。是必须考虑的问题。是必须考虑的问题。下边通过几个算例，以便进一步了解不同单元的特点。下边通过几个算例，以便进一步了解不同单元的特点。下边通过几个算例，以便进一步了解不同单元的特点。下边通过几个算例，以便进一步了解不同单元的特点。例例例例4.3 4.3 4.3 4.3 一悬矩形截面臂梁，长一悬矩形截面臂梁，长一悬矩形截面臂梁，长一悬矩形截面臂梁，长400mm400mm400mm400mm，高高高高40mm40mm40mm40mm，厚，厚，

14、厚，厚10mm10mm10mm10mm，左端固定，右端部作用左端固定，右端部作用左端固定，右端部作用左端固定，右端部作用1000N1000N1000N1000N的负的负的负的负 方向方向方向方向的集中力，求其根部的最大正应力。的集中力，求其根部的最大正应力。的集中力，求其根部的最大正应力。的集中力，求其根部的最大正应力。第11页/共54页第十二页，共54页。例例例例4.4 4.4 分析带有圆孔的平板在单向拉伸下的应力分布（本题分析带有圆孔的平板在单向拉伸下的应力分布（本题分析带有圆孔的平板在单向拉伸下的应力分布（本题分析带有圆孔的平板在单向拉伸下的应力分布（本题(bnt)(bnt)与例与例与例

15、与例2.42.4相同）相同）相同）相同）理论解理论解理论解理论解 FEM FEM8节点(ji din)单元24个单元(dnyun)93个节点第12页/共54页第十三页，共54页。例例例例4.5 4.5 受均布内压长圆筒的应力分析（内压受均布内压长圆筒的应力分析（内压受均布内压长圆筒的应力分析（内压受均布内压长圆筒的应力分析（内压100MPa100MPa）文件文件文件文件(wnjin)(wnjin)：W4-5a.txt(4W4-5a.txt(4节点轴对称单元节点轴对称单元节点轴对称单元节点轴对称单元)文件文件文件文件(wnjin)(wnjin)：W4-5.txt(8W4-5.txt(8节点轴对称

16、单元节点轴对称单元节点轴对称单元节点轴对称单元)第13页/共54页第十四页，共54页。4.2 4.2 4.2 4.2 三维实体等参单元三维实体等参单元三维实体等参单元三维实体等参单元（1 1 1 1）按单元的节点数，三）按单元的节点数，三）按单元的节点数，三）按单元的节点数，三 维实体等参元可分为：维实体等参元可分为：维实体等参元可分为：维实体等参元可分为：8 8 8 8节点单元节点单元节点单元节点单元 20 20 20 20节点单元节点单元节点单元节点单元 8-20 8-20 8-20 8-20节点单元节点单元节点单元节点单元 退化单元退化单元退化单元退化单元 6 6 6 6节点单元节点单元

17、节点单元节点单元 14 14 14 14节点单元节点单元节点单元节点单元 4 4 4 4节点单元节点单元节点单元节点单元 10 10 10 10节点单元节点单元节点单元节点单元 采用单元的原则：采用单元的原则：采用单元的原则：采用单元的原则：应力变化、几何形状应力变化、几何形状应力变化、几何形状应力变化、几何形状(xngzhun)(xngzhun)(xngzhun)(xngzhun)、疏密过度、疏密过度、疏密过度、疏密过度 特点：统一格式、同一程序路径特点：统一格式、同一程序路径特点：统一格式、同一程序路径特点：统一格式、同一程序路径第14页/共54页第十五页，共54页。（2 2）三维等参元的

18、位移函数）三维等参元的位移函数(hnsh)(hnsh)点、棱、面一一对应点、棱、面一一对应坐标变换公式坐标变换公式第15页/共54页第十六页，共54页。积分域是边长为积分域是边长为2 2的立方体，积分限是的立方体，积分限是-1-1到到1 1，即积分域规则，即积分域规则，积分限简单。积分限简单。这里刚度矩阵这里刚度矩阵 是是24x24 24x24 的矩阵，其计算一般采用高斯的矩阵，其计算一般采用高斯数值积分。数值积分。是在三个方向是在三个方向(fngxing)(fngxing)的高斯积分权系数的高斯积分权系数 是在三个方向是在三个方向(fngxing)(fngxing)的高斯积分点数的高斯积分点

19、数（3）三维等参元的刚度(n d)矩阵总体坐标系和局部坐标系之间的转换关系单元刚度(n d)矩阵可以写为被积函数是局部坐标 的函数。第16页/共54页第十七页，共54页。（4 4）三维等参单元的等效节点载荷）三维等参单元的等效节点载荷 单元等效节点载荷的计算原理和方法与单元等效节点载荷的计算原理和方法与2.2.42.2.4节相同节相同(xin tn)(xin tn)，计算公式也具有相同，计算公式也具有相同(xin tn)(xin tn)的形式，区别主要在于等参单元的等的形式，区别主要在于等参单元的等效节点载荷计算中的有关积分是在局部坐标系下进行的。效节点载荷计算中的有关积分是在局部坐标系下进行

20、的。1.1.体积力体积力2.2.表面力表面力（1 1）确定表面力表达式）确定表面力表达式 第17页/共54页第十八页，共54页。（2）确定(qudng)有表面力作用面的法线方向（3）确定(qudng)表面力列阵 表面力垂直于单元表面，假定压力为正，拉力为负，则将式(4.57)、(4.58)代入式（4.55）,得表面力等效节点载荷计算公式：第18页/共54页第十九页，共54页。（5）算例及应用 有限元求解过程中应该注意的几个问题：（1）网格的边长比不宜太大，形状不能过度扭曲。平面网格越接近正方形，三维网格越接近正立方体，单元性能越好。（2）要保证单元边界上的位移协调。既使使用协调单元，如单元编号

21、不恰当，也可能产生单元边界上的位移不协调情况，如对图4-19中的平面单元，如果单元1的节点号是1、3、8、6、2、5、7、4 单元2的节点号是6、8、10、9 两个单元的交界处位移不协调。单元2的节点号应该是6、8、10、9、7（5节点元），则协调。（3）如果出现单元的雅可比行列式的值为负的情况，计算就无法(wf)进行，单元过 分扭曲、翘曲、单元内角接近于零度或180度、节点编号反向等都会造成 这种情况。（4）高斯积分点的应力精度高于节点的应力精度。（5）通过节点重合形成的退化单元，不能计算节点上的应力。第19页/共54页第二十页，共54页。例例4.6 图示圆筒，底部固定，图示圆筒，底部固定，

22、高度高度0.82m，筒的外半径，筒的外半径(bnjng)0.25m，筒壁厚度筒壁厚度 0.005m。筒内装满水，。筒内装满水，求在静水压力下，筒壁中径处的求在静水压力下，筒壁中径处的径向位移和外壁环向应力。径向位移和外壁环向应力。第20页/共54页第二十一页，共54页。（6 6）高斯积分阶数的选取。从应用的角度讲，高斯积分阶数的选取主要应考两个方面，）高斯积分阶数的选取。从应用的角度讲，高斯积分阶数的选取主要应考两个方面，一是计算精度的需要，二而是计算量和结果输出量。决定单元精度的是完全多项式的方一是计算精度的需要，二而是计算量和结果输出量。决定单元精度的是完全多项式的方次，对于含有非完全多项

23、式的单元，并非积分阶数越高精度也越高。次，对于含有非完全多项式的单元，并非积分阶数越高精度也越高。例例4.7 4.7 四十齿齿轮强度有限元计算四十齿齿轮强度有限元计算 图示为齿轮有限元网格，材料为图示为齿轮有限元网格，材料为40CrNiMoA40CrNiMoA，齿轮轴孔表面固定。载荷为施加在齿轮上的啮合力，在其中一个齿的齿顶分别施加切向齿轮轴孔表面固定。载荷为施加在齿轮上的啮合力，在其中一个齿的齿顶分别施加切向力力74.2574.25吨，径向力吨，径向力(指向圆心指向圆心(yunxn)34.62(yunxn)34.62吨。对其进行强度计算与校核。吨。对其进行强度计算与校核。第21页/共54页第

24、二十二页，共54页。5 杆系结构单元杆系结构单元杆只承受轴力杆只承受轴力一维杆一维杆 一维杆单元一维杆单元(dnyun)(dnyun)平面杆系平面杆系 平面杆单元平面杆单元(dnyun)(dnyun)空间杆系空间杆系 空间杆单元空间杆单元(dnyun)(dnyun)梁能够承受轴力、弯矩、扭矩梁能够承受轴力、弯矩、扭矩 与坐标轴平行的梁与坐标轴平行的梁 梁单元梁单元(dnyun)(dnyun)平面梁平面梁 平面梁单元平面梁单元(dnyun)(dnyun)空间梁空间梁 空间梁单元空间梁单元(dnyun)(dnyun)在有限元法中有：在有限元法中有：考虑轴力的梁单元考虑轴力的梁单元(dnyun)(d

25、nyun)公式公式不考虑轴力的梁单元不考虑轴力的梁单元(dnyun)(dnyun)公式公式第22页/共54页第二十三页，共54页。5.1 5.1 杆单元杆单元5.1.1 5.1.1 一维杆单元一维杆单元 单元位移函数单元位移函数杆单元的横截面积为杆单元的横截面积为A A，长度为，长度为 L L，弹性摸量为，弹性摸量为 E E。单元的节点位移单元的节点位移:单元位移函数单元位移函数单元节点位移条件单元节点位移条件(tiojin):(tiojin):由式（由式（6.36.3）解出）解出 ，代入（，代入（6.26.2）得）得 单元应变单元应变应变矩阵：应变矩阵：单元应力单元应力 单元刚度矩阵单元刚度

26、矩阵应力(yngl)矩阵：单元刚度(n d)方程：单元刚度矩阵：单元节点力：第23页/共54页第二十四页，共54页。刚度矩阵的物理意义刚度矩阵的物理意义杆的轴向刚度：杆的轴向刚度：单元单元(dnyun)(dnyun)刚度方程：刚度方程：当当 时，时，当当 时，时，第24页/共54页第二十五页，共54页。5.1.2 5.1.2 平面和空间平面和空间平面和空间平面和空间(kngjin)(kngjin)杆单元杆单元杆单元杆单元平面杆单元平面杆单元平面杆单元平面杆单元平面杆单元：可以处在平面内的任意平面杆单元：可以处在平面内的任意平面杆单元：可以处在平面内的任意平面杆单元：可以处在平面内的任意方位，与

27、坐标轴平行，或不平行。方位，与坐标轴平行，或不平行。方位，与坐标轴平行，或不平行。方位，与坐标轴平行，或不平行。取两个直角坐标系：取两个直角坐标系：取两个直角坐标系：取两个直角坐标系：夹角为夹角为夹角为夹角为 。（1 1）局部系下节点力与节点位移关系）局部系下节点力与节点位移关系）局部系下节点力与节点位移关系）局部系下节点力与节点位移关系令令令令局部坐标系下单元刚度方程局部坐标系下单元刚度方程局部坐标系下单元刚度方程局部坐标系下单元刚度方程（单元节点力与节点位移关系）（单元节点力与节点位移关系）（单元节点力与节点位移关系）（单元节点力与节点位移关系）第25页/共54页第二十六页，共54页。（2

28、 2）两种坐标系下节点力与节点位移关系）两种坐标系下节点力与节点位移关系(gun x)(gun x)局部坐标系下，节点力与节点位移关系局部坐标系下，节点力与节点位移关系(gun x)(gun x)：两种坐标系下，节点位移关系两种坐标系下，节点位移关系(gun x)(gun x)：两种坐标系下，节点力关系两种坐标系下，节点力关系(gun x)(gun x)：将式（将式（5.6)5.6)和和(5.7)(5.7)代入式（代入式（5.5)5.5)，得，得第26页/共54页第二十七页，共54页。用 左乘式（5.8）两端(lin dun)，得总体坐标系下、任意方位平面杆单元的节点力与节点位移关系：其中：5

29、.1.3 空间杆单元第27页/共54页第二十八页，共54页。例例 5.1 图示为一平面桁架，假定弹性模量图示为一平面桁架，假定弹性模量(tn xn m lin)为为 E，各杆的横截面积为各杆的横截面积为 A，为为45O，在点，在点1的负的负Y 方向作用力为方向作用力为F。求点。求点1的竖向位移。的竖向位移。(W5-1.txt)3个单元的个单元的刚度矩阵：刚度矩阵：因节点因节点2、3、4的位移等于零，只需组装与节点的位移等于零，只需组装与节点1有关的方程：有关的方程：所求结果与材料力学理论解相同。所求结果与材料力学理论解相同。令令 ,则则 ，有限元解：有限元解：节点(ji din)1节点(ji

30、din)1节点1节点2节点3节点4第28页/共54页第二十九页，共54页。5.2 5.2 梁单元梁单元梁单元梁单元与坐标轴平行的平面梁单元与坐标轴平行的平面梁单元与坐标轴平行的平面梁单元与坐标轴平行的平面梁单元能承受轴向力能承受轴向力能承受轴向力能承受轴向力NN、弯矩、弯矩、弯矩、弯矩MM和横向剪力和横向剪力和横向剪力和横向剪力QQ的梁单元的梁单元的梁单元的梁单元1.1.单元位移函数单元位移函数单元位移函数单元位移函数每个节点每个节点每个节点每个节点(ji di(ji di n)n)有有有有3 3个广义位移个广义位移个广义位移个广义位移(自由度自由度自由度自由度)：沿沿沿沿 x x轴和轴和轴和

31、轴和 y y轴方向的位移轴方向的位移轴方向的位移轴方向的位移 、绕、绕、绕、绕z z 轴的转角轴的转角轴的转角轴的转角:与与与与x x向位移有关的节点向位移有关的节点向位移有关的节点向位移有关的节点(ji di(ji di n)n)位移有位移有位移有位移有2 2个个个个 ,单元的单元的单元的单元的x x向位移取向位移取向位移取向位移取坐标的线性函数。坐标的线性函数。坐标的线性函数。坐标的线性函数。与与与与 y y向位移有关的节点向位移有关的节点向位移有关的节点向位移有关的节点(ji di(ji di n)n)位移有位移有位移有位移有4 4个：个：个：个：,单元的单元的单元的单元的y y 向位移

32、取坐标的向位移取坐标的向位移取坐标的向位移取坐标的3 3次函数，即次函数，即次函数，即次函数，即常数常数常数常数 由单元节点由单元节点由单元节点由单元节点(ji di(ji di n)n)位移条件确定：位移条件确定：位移条件确定：位移条件确定：由（由（由（由（5.115.11）解出）解出）解出）解出 由（由（由（由（5.125.12）解出）解出）解出）解出 代入式（代入式（代入式（代入式（5.105.10）得到：）得到：）得到：）得到：第29页/共54页第三十页，共54页。将式位移u,v两式合写在一起，以矩阵形式(xngsh)表示，则有第30页/共54页第三十一页，共54页。2.单元单元(dn

33、yun)应变和单元应变和单元(dnyun)应力应力单元单元(dnyun)应变应变 应变可以分为两部分应变可以分为两部分 ：拉压应变拉压应变 和弯曲应变和弯曲应变 。略去剪切变形的影响，则单元略去剪切变形的影响，则单元(dnyun)应变为应变为令令式（式（5.14）可写为）可写为单元单元(dnyun)应力应力第31页/共54页第三十二页，共54页。3.轴线与坐标轴平行的梁单元刚度矩阵轴线与坐标轴平行的梁单元刚度矩阵单元刚度方程：单元刚度方程：令单元刚度矩阵令单元刚度矩阵单元刚度矩阵显式表达式：单元刚度矩阵显式表达式：式中，式中，A是单元的横截面面积，是单元的横截面面积，I 为横截面惯性矩，为横截

34、面惯性矩，E 为材料为材料(cilio)的弹性模量。的弹性模量。第第1行、第行、第1列对应于列对应于ui,.。第32页/共54页第三十三页，共54页。其中，其中，为剪切影响系数，为剪切影响系数，A A为梁的横截面面积，为梁的横截面面积，E E和和GG是弹性摸量和剪切摸量，是弹性摸量和剪切摸量，为截面剪切因子。矩形和圆形为截面剪切因子。矩形和圆形截面，文献截面，文献33建议建议 分别取分别取5/65/6和和10/910/9。对于高为对于高为h h 宽为宽为b b 的矩形截面，剪切影响系数的矩形截面，剪切影响系数 ，所以，当梁的高度所以，当梁的高度 相对于跨度相对于跨度 很小时，可以很小时，可以(

35、ky(ky)忽略剪切忽略剪切变形的影响。变形的影响。4.考虑剪切变形时的单元刚度矩阵考虑剪切变形时的单元刚度矩阵当梁的高度当梁的高度(god)相对于跨度不是很小时，需要考虑剪切变形的影响，为此，需要对相对于跨度不是很小时，需要考虑剪切变形的影响，为此，需要对式（式（5.21）的刚度矩阵做如下修正）的刚度矩阵做如下修正第33页/共54页第三十四页，共54页。例例5.2 按考虑和不考虑剪切变形影响两种情况，计算图示矩形截面悬臂梁在端部横向力按考虑和不考虑剪切变形影响两种情况，计算图示矩形截面悬臂梁在端部横向力 F作作用下的端点挠度。截面高为用下的端点挠度。截面高为 h，宽为，宽为b。忽略轴向变形。

36、忽略轴向变形。用一个单元，只有两个节点。用一个单元，只有两个节点。忽略轴向变形，每个节点只有忽略轴向变形，每个节点只有2个自由度，个自由度，即即y 方向的移动和绕方向的移动和绕z 轴的转动。轴的转动。节点节点1 固定，固定，总体平衡方程只有与节点总体平衡方程只有与节点2对应的对应的2个方程。个方程。按不考虑剪切变形影响，按不考虑剪切变形影响，根据式（根据式（6.63）的单元刚度矩阵，可得总体平衡方程：）的单元刚度矩阵，可得总体平衡方程：使用一个单元，就可以得到与材料力学的理论使用一个单元，就可以得到与材料力学的理论(lln)解相同解相同的位移结果，这是因为这种梁单元的挠度位移函数是的位移结果，

37、这是因为这种梁单元的挠度位移函数是三次的，悬臂梁的挠度函数也是三次的。三次的，悬臂梁的挠度函数也是三次的。考虑剪切变形影响考虑剪切变形影响总体平衡方程：总体平衡方程：矩形截面：矩形截面：当当 ，括号中第二项等于，括号中第二项等于0.54%。第34页/共54页第三十五页，共54页。平面梁单元和三维梁单元平面梁单元和三维梁单元平面梁单元：可以处在平面内的任意平面梁单元：可以处在平面内的任意方位，与坐标轴平行，或不平行。方位，与坐标轴平行，或不平行。取两个直角坐标系：取两个直角坐标系：夹角为夹角为 。（1）节点力与节点位移）节点力与节点位移局部系局部系总体系总体系(tx)（2）节点力与节点位移转换关

38、系）节点力与节点位移转换关系第35页/共54页第三十六页，共54页。令局部系下的节点位移用总体系下的表示同理，局部系下的节点力用总体系下的表示由式（5.25）知局部系下的刚度方程可写为为将式（5.26）和（5.27）代入（5.28），得到用 左乘式（5.29）两端，得总体坐标系下、任意方位(fngwi)平面梁单元的节点力与节点位移关系：第36页/共54页第三十七页，共54页。令局部系下的节点(ji din)位移用总体系下的表示同理，局部系下的节点(ji din)力用总体系下的表示由式（5.25）知局部系下的刚度方程可写为:将式（5.26）和（5.27）代入（5.28），得到用 左乘式（5.29

39、）两端，得令 则总体坐标系下、任意方位平面梁单元的节点(ji din)力与节点(ji din)位移关系：第37页/共54页第三十八页，共54页。（1 1）局部坐标系下的三维梁单元）局部坐标系下的三维梁单元三维梁单元能够承受：三维梁单元能够承受：轴力：轴力：剪力：剪力：弯矩：弯矩：扭矩：扭矩：每个节点每个节点6 6个自由度（广义）：个自由度（广义）：三个方向的位移：三个方向的位移：绕三个轴的转角：绕三个轴的转角：刚度刚度(nn d)d)方程方程式（式（5.205.20）给出了局部系下节点力）给出了局部系下节点力 与节点位移与节点位移 的关系的关系同理，可以给出局部系下节点力同理，可以给出局部系下

40、节点力 与节点位移与节点位移 的关系的关系局部系下空间梁单元刚度局部系下空间梁单元刚度(nn d)d)方程（节点力与节点位移关系）：方程（节点力与节点位移关系）：节点力：节点力：节点位移：节点位移：刚度刚度(nn d)d)矩阵：矩阵：第38页/共54页第三十九页，共54页。（2）总体坐标系下的三维梁单元）总体坐标系下的三维梁单元总体坐标系下的梁单元刚度矩阵总体坐标系下的梁单元刚度矩阵 可通过如下可通过如下(rxi)转换得到转换得到 T是三维转换矩阵，可以根据单元的总体节点坐标求出，具体计算方法可见文献是三维转换矩阵，可以根据单元的总体节点坐标求出，具体计算方法可见文献1。第39页/共54页第四

41、十页，共54页。5.3 5.3 板壳单元板壳单元板壳单元板壳单元(dnyun)(dnyun)*对于板壳结构对于板壳结构,原则上说，可以采用二维或三原则上说，可以采用二维或三维实维实 体单元对其进行结构分析，但是，可能体单元对其进行结构分析，但是，可能会出现：会出现：（1 1）单元的三个方向的尺寸相差很大，从而）单元的三个方向的尺寸相差很大，从而导致形成的有限元总体方程为病态或奇异，最后导致形成的有限元总体方程为病态或奇异，最后使解丧失精度，或无法求解。使解丧失精度，或无法求解。（2 2）为避免上述问题）为避免上述问题(wnt)(wnt)，若使单元在，若使单元在各个方向尺寸相近，将导致单元总数和

42、节点总数各个方向尺寸相近，将导致单元总数和节点总数过分庞大，也会使计算无法进行。过分庞大，也会使计算无法进行。*与板理论和壳理论对应，在有限元法中，有板与板理论和壳理论对应，在有限元法中，有板单元和壳单元。单元和壳单元。*对板壳结构应该采用板壳单元。对板壳结构应该采用板壳单元。板与壳的区别：板与壳的区别：板只考虑弯曲变形，壳考虑弯曲和中曲面内的变板只考虑弯曲变形，壳考虑弯曲和中曲面内的变形。形。本章只介绍比较简单的板单元和壳单元。本章只介绍比较简单的板单元和壳单元。第40页/共54页第四十一页，共54页。5.3.1 5.3.1 薄板的弯曲变形薄板的弯曲变形薄板的弯曲变形薄板的弯曲变形工程中平板

43、构件的长度或宽度要小得多，如果厚度与板面内的其他尺工程中平板构件的长度或宽度要小得多，如果厚度与板面内的其他尺工程中平板构件的长度或宽度要小得多，如果厚度与板面内的其他尺工程中平板构件的长度或宽度要小得多，如果厚度与板面内的其他尺寸之比小于寸之比小于寸之比小于寸之比小于1/51/5，可以认为是薄板。，可以认为是薄板。，可以认为是薄板。，可以认为是薄板。位于板厚度中点的平面称为板的中面，位于板厚度中点的平面称为板的中面，位于板厚度中点的平面称为板的中面，位于板厚度中点的平面称为板的中面，薄板所受的载荷可以分为两类：薄板所受的载荷可以分为两类：薄板所受的载荷可以分为两类：薄板所受的载荷可以分为两类

44、：作用在中面内的纵向载荷，认为沿厚度均匀分布，可以按平面应力进作用在中面内的纵向载荷，认为沿厚度均匀分布，可以按平面应力进作用在中面内的纵向载荷，认为沿厚度均匀分布，可以按平面应力进作用在中面内的纵向载荷，认为沿厚度均匀分布，可以按平面应力进行计算，行计算，行计算，行计算，垂直于中面的所谓横向载荷。横向载荷使板弯曲，可以按薄板弯曲问垂直于中面的所谓横向载荷。横向载荷使板弯曲，可以按薄板弯曲问垂直于中面的所谓横向载荷。横向载荷使板弯曲，可以按薄板弯曲问垂直于中面的所谓横向载荷。横向载荷使板弯曲，可以按薄板弯曲问题进行计算。题进行计算。题进行计算。题进行计算。这里的薄板弯曲是指弹性薄板小挠度弯曲。

45、这里的薄板弯曲是指弹性薄板小挠度弯曲。这里的薄板弯曲是指弹性薄板小挠度弯曲。这里的薄板弯曲是指弹性薄板小挠度弯曲。在忽略横向剪切变形的前提下，在忽略横向剪切变形的前提下，在忽略横向剪切变形的前提下，在忽略横向剪切变形的前提下，克希霍夫克希霍夫克希霍夫克希霍夫(Kirchhoff)(Kirchhoff)薄板弯曲理论假设：薄板弯曲理论假设：薄板弯曲理论假设：薄板弯曲理论假设：*直法线直法线直法线直法线(f(f xin)xin)假设，即变形前垂直于中面的法线假设，即变形前垂直于中面的法线假设，即变形前垂直于中面的法线假设，即变形前垂直于中面的法线(f(f xin)xin)，变形后仍垂直于中面。变形后

46、仍垂直于中面。变形后仍垂直于中面。变形后仍垂直于中面。*板的中面只发生弯曲变形，没有面内的板的中面只发生弯曲变形，没有面内的板的中面只发生弯曲变形，没有面内的板的中面只发生弯曲变形，没有面内的 伸缩变形，即中面无变形。伸缩变形，即中面无变形。伸缩变形，即中面无变形。伸缩变形，即中面无变形。*忽略沿厚度方向的挤压变形，忽略沿厚度方向的挤压变形，忽略沿厚度方向的挤压变形，忽略沿厚度方向的挤压变形，即即即即 ：第41页/共54页第四十二页，共54页。由上述(shngsh)假设，可以写出薄板弯曲问题的位移、应变、应力和内力矩：（1）位移 ，z方向挠度只是 x,y的函数，板内各点的位移为(2)应变由于

47、w仅是坐标x,y 的函数，虑 三个所以只需要考应变分量根据几何方程和式，应变可以表示为反映了弯曲程度的合成为曲率：（3）应力第42页/共54页第四十三页，共54页。（4）内力（矩）内力（矩）单位宽度板上的内力单位宽度板上的内力弯矩弯矩弯矩弯矩扭矩扭矩矩阵矩阵(j zhn)表示的内力表示的内力薄板弯曲弹性矩阵薄板弯曲弹性矩阵(j zhn)（试与梁的弯曲刚度比较）（试与梁的弯曲刚度比较）应力(yngl)与内力的关系最大应力在板上下表面(biomin)（与梁的弯曲应力比较，排在一起的梁与板比较）第43页/共54页第四十四页，共54页。5.3.2 三节点三节点 三角形弯曲三角形弯曲(wnq)板单元板单

48、元三角形板单元易于模拟复杂边界。三角形板单元易于模拟复杂边界。三角形板单元种类较多。三角形板单元种类较多。节点位移（广义自由度）节点位移（广义自由度）Z向挠度：向挠度：；绕绕x轴转角：轴转角：；绕绕y轴转角：轴转角：3个节点，个节点，9个位移分量：个位移分量：如果单元位移函数取如果单元位移函数取 x、y的多项式，则多项式最多只能取的多项式，则多项式最多只能取9项，而项，而 x、y的三次的三次完全多项式有完全多项式有10项，即项，即*必须减少一项。必须减少一项。前前6项代表刚体位移和常应变，必须保留。项代表刚体位移和常应变，必须保留。*后后4项去掉一项则破坏对称性，有人尝试将项去掉一项则破坏对称

49、性，有人尝试将 合并，但当单元是两个边合并，但当单元是两个边 分别平行于分别平行于x、y 轴的等腰三角形时，确定未知参数轴的等腰三角形时，确定未知参数 的系数矩阵奇异，故不可的系数矩阵奇异，故不可行。行。*采用面积坐标可以构成不同的单元位移函数，以克服上述困难。采用面积坐标可以构成不同的单元位移函数，以克服上述困难。第44页/共54页第四十五页，共54页。三节点(ji din)三角形弯曲板单元刚度方程第45页/共54页第四十六页，共54页。*这种位移模式在相邻单元间的挠度和沿边界切线方向的斜率是连续的，但是 其 法线方向的斜率是不连续的，因此，这种单元是一种完备的非协调单元。*采用不同形状的单

50、元、采用不同的节点数、采用不同的单元位移函数等方法，可以构造(guzo)出不同的板单元。例如，对于本节的三节点三角形单元，可以构造(guzo)出在相邻单元边界上完全协调的基于离散Kirchhoff理论的板单元2,3。此外，还有其它形式的板单元3,6。*板单元有协调元和非协调元，只要非协调单元能通过分片试验，则单元是可以 使用的，对于许多工程问题，用非协调元得到的解的精度还是能满足要求的，*有时候还能给出比协调元更好一些的结果，这是因为利用最小位能原理求得的 近似解一般使结构显得过于刚硬，而非协调元未精确满足最小势能原理的要求，使结构趋于柔软，正好抵消了上述过于刚硬带来的误差。*对于薄板弯曲问题