《第4讲等参元和高斯积分.ppt》由会员分享,可在线阅读,更多相关《第4讲等参元和高斯积分.ppt(30页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第4讲等参元和高斯积分 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望北京航空航天大学l实际问题常常需要使用一些几何形状不太规整的单元来逼近原问题。直接研究这些不规整单元的表达式比较困难(在整体坐标系下构造位移插值函数,则计算形状函数矩阵、单元刚度矩阵及等效节点载荷列阵时十分冗繁)。事实上,形状不规整的单元和形状规整的单元(矩形单元、正六面体单元)可以建立一种映射关系,使得物理坐标系中的整体坐标和自然坐标系中的局部坐标一一对应。l等参单元等参单元的提出为有限元法成
2、为现代工程实际领域最有效的数值分析方法迈出了极为重要的一步。北京航空航天大学4.1 等参单元l等参单元定义的给出等参单元定义的给出l平面问题四边形等参单元计算公式平面问题四边形等参单元计算公式l三维问题六面体等参单元计算公式三维问题六面体等参单元计算公式l采用等参单元的优点采用等参单元的优点北京航空航天大学等参单元定义的给出l等参单元:用同样的节点和相同的形状函数通过插值的方式表示出单元的几何坐标与位移的单元,称为等参单元。l等参单元的插值函数用自然坐标给出。l如果坐标变换节点数多于位移插值的节点数,称为超参变换。反之,如果坐标变换节点数少于位移插值的节点数,则称为亚参变换。北京航空航天大学平
3、面问题四边形等参单元的推导整体直角坐标整体直角坐标单元局部自然坐标(一般四边形)(规格化的矩形)映射坐标映射北京航空航天大学映射节点条件:构造插值函数北京航空航天大学北京航空航天大学北京航空航天大学节点条件:位移函数北京航空航天大学同理可得:北京航空航天大学北京航空航天大学单元的几何坐标与位移用同样的节点和相同的形状函数通过插值的方式表示。形状函数用自然坐标给出。北京航空航天大学?北京航空航天大学偏导数变换雅可比矩阵:北京航空航天大学北京航空航天大学四边形等参单元形状要求不能有重节点不能出现内角大于180o的情况内角最好介于30o-150o之间(有限变形的情况)避免出现北京航空航天大学三维问题
4、六面体等参单元的计算公式北京航空航天大学北京航空航天大学北京航空航天大学采用等参单元的优点l借助于等参元可以对于一般的任意几何形状的工程问题方便地进行有限元离散。l等参元的插值函数是用自然坐标给出的,等参元的一切计算都是在自然坐标系中规格化的母单元内进行,相关运算大大简化。l不管各个积分形式的矩阵的被积函数如何复杂,都可以采用标准化的数值积分方法计算,从而使工程问题的有限元分析纳入了统一的通用化程序。北京航空航天大学4.2 数值积分l数值积分及其基本思想lNewton-cotes积分公式lGauss-Legendre积分公式l等参元中积分阶次的选择北京航空航天大学关于数值积分l计算刚度矩阵及等
5、效节点载荷列阵的元素时,往往涉及到复杂函数的定积分,在有限元分析中广泛采用数值积分方法。l数值积分方法是一种近似的方法。一个函数的定积分可以通过n个结点的函数值的加权组合来表示北京航空航天大学数值积分的基本思想北京航空航天大学求积公式插值法至少具有n-1次代数精度北京航空航天大学Newton-cotes求积公式l如果n个结点 等距分布,则前面的插值型求积公式称为Newton-cotes求积公式。lNewton-cotes求积公式具有n-1次代数精度l几个常用求积公式梯形公式,n=1Simpson公式,n=2北京航空航天大学Gauss-Legendre求积公式ln个插值结点非等距分布l结点和积分
6、权系数可以查表北京航空航天大学l高斯积分方法预先定义了积分点和相应的加权系数,求出被积分的函数在指定积分点上的数值,加权后求和,就得到了该函数的积分。l高斯积分方法具有最高的计算精度。采用n个积分点的高斯积分可以达到2n-1阶的精度,也就是说,如果被积分的函数是2n-1次多项式,用n个积分点的高斯积分可以得到精确的积分结果。北京航空航天大学等参元高斯求积公式的一般形式北京航空航天大学等参元中积分阶次的选择l积分阶次的选择直接影响计算的精度和计算工作量。l积分阶次的选择必须保证积分的精度。(完全精确积分完全精确积分)l很多情况下,实际选取的高斯积分点数低于精确积分的要求,往往可以取得较完全精确积分更好的精度。(减缩积分减缩积分)北京航空航天大学线性单元完全精确积分 二次单元减缩积分