有限元教程5.pptVIP

五、等参数单元 * §5.1 什么是等参数单元 背景： 解决边界条件复杂的问题，要获得较高的精度 什么是等参数单元？ 实际单元 局部坐标系－母单元 建立一一对应关系 用局部坐标系中的插值函数在（ξ,η）值和已知结点（实际单元结点）的坐标（x1,y1）… （x4,y4）线性组合得到p（x,y） （1）坐标变换 以上过程叫做坐标变换 （2）位移插值 在母单元中，按照： 位移模式选择→求解待定系数→形函数 如果： →等参数单元 证明:能否通过坐标变换将实单元与母单元之间建立一一对应关系？ 线？ 整体坐标系中3－4边 同理： ξ－η中每一条直线都唯一对应着实单元中一直线→通过坐标变换将实单元与母单元之间建立一一对应关系，在局部坐标系中研究问题 §5.2 实际单元的位移插值函数 形函数矩阵[N] 形函数 §5.3几何矩阵 根据隐函数求导 §5.4单刚 根据隐函数求导得到[B]中元素 单刚。全部元素为（ζη）函数，很难得到解析积分结果→借助高斯积分 数值积分是利用被积函数在若干点上的函数值来近似计算定积分的方法。 例 梯形求积公式 从几何上看，求定积分就是求曲边梯形的面积。最简单的办法，可以用梯形面积近似代替曲边梯形面积。从而有 O a b x y y=f(x) y=p1(x) f(a) f(b) A B §5.5高斯积分 数值积分的普遍形式 近似计算定积分准确吗？→代数精度 能使上式精确成立的多项式的最高幂次 权系数 求积结点 高斯积分结点个数可取2、3、4、5…，相应的代数精度为3、5、7、9…。 高斯积分点及权系数 n=2时高斯积分点 二维高斯数值积分 n=3时高斯积分点(9) 高斯积分求解四结点等参元单刚过程： 输入实际单元结点坐标 计算[J]中系数 A、B、β1～ β4 计算[D] 确定高斯积分阶数n For I=1:n for J=1:n [J]= for k=1:4 [B]k end [k]=[k]+[B]T×[D] × [B] ×[J] ×hi ×hj end end 练习六： 图式四结点等参数单元，材料弹性模量E=200000Pa，u=0.2，板厚为单位1，平面应力问题。请计算结点位移 提示：可以用有限元软件计算结果验证 *