应用有限元分析工程实例3.pptxVIP

第 3 讲4 实体单元及其特点和应用; 4节点单元： 在单元内，位移随坐标线性变化，应变和应力为常数。 8节点单元（ 9节点单元）： 在单元内，位移为坐标的二次函数，应变和应力随坐标线性变化。可以模拟复杂边界，精度高。 5-7节点单元： 用于单元网格过度。 退化单元（3-6节点单元） 易于网格剖分，单元网格过度。 （2）按适用的力学问题，二维实体等参单元可分为： 平面应力单元、平面应变单元、平面轴对称单元。 平面应力单元又分为单位厚度和制定厚度两种。 ;（3）二维等参元的位移函数 * 坐标转换 两个单元的4个节点和4条边一一对应。 采用坐标变换式，把实际单元映射为 一个边长为2的正方形单元。 为形函数矩阵 用实际单元的节点坐标 ，通过插值方法来表示实际单元内任意一点的坐标（x,y）。 因为形函数是用局部坐标定义的，式(4.1)建立了正方形单元内任意点 与实际单元内任意点（x,y）之间的意义对应关系。 （4.2）中形函数构造规律：直线方程; * 单元位移函数 为形函数矩阵 为单元节点位移列阵 * 等参元定义（按物理场和几何描述） 等参元：坐标变换和单元位移函数采用相同的节点数，且采用相同 的形函数 超参元：坐标变换采用相同的节点数多于单元位移函数的节点数 亚参元：坐标变换采用相同的节点数少于单元位移函数的节点数 ;（4）单元应变、单元应力 * 单元应变 * 单元应力 （5）程序实现 平面应力单元、平面应变单元、平面轴对称单元可以采用同一类单元计算公式，并且在有限元程序中，这三种问题一般采用的是同一个程序路径和同一组子程序。 ; （6）单元刚度矩阵 分别表示单元的体积、面积、厚度。按复合函数链式求导规则 雅可比矩阵 式（4.12）两边左乘 ，得 被积函数是局部坐标 的函数，积分域是边长为2的正方形，积分限是-1到1，即积分域规则，积分限简单。 这里刚度矩阵是 8x8 ???矩阵，其计算一般采用高斯数值积分。;（7） 二维等参单元的等效节点载荷 计算原理和方法与三角形节相同，区别主要在于积分是在局部坐标系下进行的。 * 体积力: 体积力等效的节点载荷公式： A是单元面积，t 是单元厚度。 体积力列矢量： 因为 ，所以，积分域由实际单元的面积A，变为局部坐标系下的边长为2的规则正方形，并且被积分函数是局部坐标 的函数。; * 体积力: 体积力等效的节点载荷公式： A是单元面积，t是单元厚度。 体积力列矢量： 被积分函数是局部坐标 的函数。 ;* 表面力 等效节点载荷计算公式 表示载荷作用的边界 载荷作用边的法矢量（假作用在边12上）： 表面力表达式 微弧长： 为8x1的列阵，但是，只有与载荷 所在边的节点对应的行才有等效节点载荷。 ;例4.1 计算图所示单元的表面力的等效节点载荷（假定单元厚度为 ） 表面力的等效节点载荷可按式（4.25）计算 因为积分只沿 的边进行 在该边上 ， ,所以 平均分配到节点1、2的负y方向上。 在其他边上亦可类似方法处理。 沿边界变化的面力 ; 计算公式 251边对应于 的边 ， 积分只沿 的边进行， 只有 不等于零 ; ;平面单元算例 实际问题常简化为平面问题，对于一个平面问题，可以采用不同的平面单元进行有限元分析，可以采用不同节点数的三角形单元、四边形单元，也可以将三角形单元和四边形单元组合使用。 在力学模型确定之后，选择什么样的单元和如何划分网格是必须考虑的问题。 下边通过几个算例，以便进一步了解不同单元的特点。 例4.3 一悬矩形截面臂梁，长400mm， 高40mm，厚10mm， 左端固定，右端部作用1000N的负 方向 的集中力，求其根部的最大正应力。 ;例4.4 分析带有圆孔的平板在单向拉伸下的应力分布（本题与例2.4相同） 理论