[工学]北航有限元第4讲 等参元和高斯积分.ppt

第4讲 等参单元和数值积分 实际问题常常需要使用一些几何形状不太规整的单元来逼近原问题。直接研究这些不规整单元的表达式比较困难（在整体坐标系下构造位移插值函数，则计算形状函数矩阵、单元刚度矩阵及等效节点载荷列阵时十分冗繁）。事实上，形状不规整的单元和形状规整的单元（矩形单元、正六面体单元）可以建立一种映射关系，使得物理坐标系中的整体坐标和自然坐标系中的局部坐标一一对应。 等参单元的提出为有限元法成为现代工程实际领域最有效的数值分析方法迈出了极为重要的一步。 4.1 等参单元 简单杆系问题分析的新途径 等参单元定义的给出 平面问题四边形等参单元计算公式 三维问题六面体等参单元计算公式 采用等参单元的优点 简单杆系问题分析之新途径 关于坐标系 当点的运动轨迹已知时，通常采用自然法确定点的运动规律、速度、加速度。 在自然坐标系中表示质点速度，是非常简单的，因为无论质点处在什么位置上速度都只有切向分量，而没有法向分量。 等参单元定义的给出 等参单元：用同样的节点和相同的形状函数通过插值的方式表示出单元的几何坐标与位移的单元，称为等参单元。 如果坐标变换节点数多于位移插值的节点数，称为超参变换。反之，如果坐标变换节点数少于位移插值的节点数，则称为亚参变换。 等参单元的插值函数用自然坐标给出。 平面问题四边形等参单元的推导 四边形等参单元形状要求 三维问题六面体等参单元的计算公式 采用等参单元的优点 借助于等参元可以对于一般的任意几何形状的工程问题方便地进行有限元离散。 等参元的插值函数是用自然坐标给出的，等参元的一切计算都是在自然坐标系中规格化的母单元内进行，相关运算大大简化。 不管各个积分形式的矩阵的被积函数如何复杂，都可以采用标准化的数值积分方法计算，从而使工程问题的有限元分析纳入了统一的通用化程序。 4.2 数值积分 数值积分及其基本思想 Newton-cotes积分公式 Gauss-Legendre积分公式 等参元中积分阶次的选择 关于数值积分 计算刚度矩阵及等效节点载荷列阵的元素时，往往涉及到复杂函数的定积分，在有限元分析中广泛采用数值积分方法。 数值积分方法是一种近似的方法。 数值积分的基本思想 求积公式—插值法 Newton-cotes求积公式 如果n个结点 等距分布，则前面的插值型求积公式称为Newton-cotes求积公式。 Newton-cotes求积公式具有n-1次代数精度 几个常用求积公式 梯形公式，n=1 Simpson公式，n=2 Gauss-Legendre求积公式 n个插值结点非等距分布 结点和积分权系数可以查表 高斯积分方法预先定义了积分点和相应的加权系数，求出被积分的函数在指定积分点上的数值，加权后求和，就得到了该函数的积分。 高斯积分方法具有最高的计算精度。采用n个积分点的高斯积分可以达到2n-1阶的精度，也就是说，如果被积分的函数是2n-1次多项式，用n个积分点的高斯积分可以得到精确的积分结果。 等参元高斯求积公式的一般形式 等参元中积分阶次的选择 积分阶次的选择直接影响计算的精度和计算工作量。 积分阶次的选择必须保证积分的精度。（完全精确积分） 很多情况下，实际选取的高斯积分点数低于精确积分的要求，往往可以取得较完全精确积分更好的精度。（减缩积分） 一个函数的定积分可以通过n个结点的函数值的加权组合来表示 至少具有n-1次代数精度 北京航空航天大学 北京航空航天大学 金朝海 jch666@ 途经1：在整体直角坐标系下进行单元分析（参看第3讲内容） 途径2：建立局部自然坐标系进行单元分析 直角坐标系( x , y , z) 极坐标(r,?) ，2维 球坐标系(r,θ, ?) 柱坐标系 (?, ?, z) 自然坐标系 自然坐标系: 选轨迹上任一点O为原点 用轨迹长度S 描写质点位置 O m S 质点沿切线前进方向的单位矢量为 切向单位矢量(tangential unit vector) 质点与切向正交且指向轨迹曲线凹侧的 单位矢量为法向单位矢量(normal unit vector) 新途径：建立局部自然坐标系进行单元分析 坐标插值函数： 局部自然坐标和整体直角坐标可以建立一种映射关系 节点条件： xi xj 单元内坐标由节点坐标插值表示 局部坐标到物理坐标的变换 单元位移函数： 节点条件： 观察： 单元的几何坐标与位移用同样的节点和相同的形状函数通过插值的方式表示。 形状函数是用自然坐标给出的，表达式很简单 单元应变能： 单元刚度矩阵 单元外力功 等效节点力 对于本例 自然坐标系下的分析结果与整体直角坐标系下的分析结果完全相同。 忽略单元间作用力 整体直角坐标 单元局部自然坐标 （一般四边形） （规格化