等参数单元等参元学习教案

《等参数单元等参元学习教案》由会员分享，可在线阅读，更多相关《等参数单元等参元学习教案（65页珍藏版）》请在装配图网上搜索。

1、会计学1等参数等参数(cnsh)单元等参元单元等参元第一页，共65页。 对于矩形单元，由于它采用了双线性位移模式，使得单元内的应力和应变不是常量而是按线性变化，它比常应变三角形单元能较好地反映出结构(jigu)的实际应力分布状态，但是它很难适应曲线边界和非正交的直线边界；同时在划分单元时，改变单元的大小也很困难，即不便于在不同部位采用大小不同的单元，因为已把每个单元的边长之半 作为常量而引入单元刚度矩阵中（见式(2.48)）。因此，矩形平面单元未能在实际中得到广泛的应用。 为此，我们希望找到一种单元，一方面它具有较高次的位移模式，能更好地反映结构(jigu)的复杂应力分布状态，即或是单元网格划

2、分的比较疏些，也可以得到比较好的计算精度；另一方面，它又能很好地适应曲线边界和非正交的直线边界。等参元就具备了上述两条优点，因而得到广泛应用。ba,第1页/共65页第二页，共65页。niiiniiivNvuNu11,niiiniiiyNyxNx11, 前面已谈到：无论是三角形单元还是矩形单元，其单元内位移用形函数(hnsh)表示为实际上不难证明：单元内任一点(y din)的坐标同样有上述关系，即 (3-1) (3-2)第2页/共65页第三页，共65页。 可见，常应变三角形单元和矩形单元内任一点的位移函数插值公式与该点的位置坐标变换式，都具有完全相同的形式。它们都是用同样个数的相应结点值（结点位

3、移值或坐标值）作为参数，并且用完全相同的形函数作为这些(zhxi)结点值前面的系数项。当参数取为结点位移时就得到位移函数插值公式；当参数取为结点坐标时，就得到位置坐标插值公式（或位置坐标变换式）。 常应变三角形单元和矩形单元的这种位移函数插值公式与位置坐标变换式之间的对应协调关系，就是等参元的基本特征。所以，等参元的基本概念可简单概括成：一个单元的位移函数插值结点数与其位置坐标变换结点数相等，其位移函数插值公式与位置坐标变换式都用相同的形函数与结点参数进行插值者，称为等参元。 显然，常应变三角形单元和矩形单元就是两种最简单的等参元。但是，本章所要研究的等参元，并不是这种单元，而是4结点任意四边

4、形等参元和8结点曲边四边形单元。第3页/共65页第四页，共65页。xoy4 , 3 , 2 , 1,iyxii1o1,1111,1331,144 由前述知，具有双线性位移模式的矩形单元只适用于正交的、规则形状的结构。对于(duy)非正交的、不规则形状，可以用任意四边形单元代替矩形单元进行有限元分割。 在直角坐标系（又称整体坐标系） 中，任取一任意四边形单元1,2,3,4，四边形的四个角点取为结点，各结点的直角坐标值为 。对于(duy)这种任意四边形等参元，可令其实际形状所构成的单元为子单元，把子单元的各边中点连线做一个局部坐标系（或称自然坐标系） ，且令单元各结点的局部坐标系分别是： ； ；

5、； 。这样，就把子单元影射到局部坐标系上，而成为正方形单元，称此正方形单元为母单元。整体坐标系适用于所有单元，即适用于整个求解区，而局部坐标系只适用于每一个单元。,1212一位移插值函数式及坐标(zubio)变换式第4页/共65页第五页，共65页。1oxoy 在子单元上再作各对边的等分线，这些等分线影射(yngsh)到母单元上，也必然是母单元各对应边上的等分线。这样，母单元与子单元之间的相应点存在着一一对应的关系。这种对应关系说明，在母单元 平面上平行于 或 的直线，在 平面内的子单元上仍然是相对应的直线。第5页/共65页第六页，共65页。niiiiiivNvuNu141),(,),(4141

6、41),(1,),(,),(iiiiiiiiNNN4141),(,),(iiiiiiyNyxNx 因此，我们就可以把矩形单元的位移函数插值式(3-1)，单元内任一点的坐标变换式(3-2)，以及局部坐标变换式（类似坐标变换式），用在任意(rny)四边形等参元上，并重新写成(3-3)(3-4) (3-5)第6页/共65页第七页，共65页。 式(3-3)和(3-4)是任意四边形在局部坐标系下的位移插值函数和单元内任一点局部坐标插值公式，而式(3-5)是每个单元的局部坐标系与结构的整体坐标系之间的坐标变换式。由这些公式看出，任意四边形单元符合等参元条件，它当然是等参元。 由于任意四边形单元的位移插值函

7、数(3-3)，在局部坐标系下满足形容条件，因此坐标变换式(3-5)也就满足相容条件，从而使得式(3-3)在整体坐标下满足相容条件。也就是说，在两相邻任意四边形单元公共边上的位移是连续的，坐标变换后仍然是连续的，两相邻单元公共边上的公共点在坐标变换后仍为公共点，决不会出现重叠和开裂(ki li)现象。第7页/共65页第八页，共65页。xoy BiNxy ATedxdytBDBkiNxydxdy),(iN Bdxdy 利用任意四边形等参元分析平面(pngmin)问题时，有了该单元的位移插值函数式(3-3)和坐标变换式(3-5)，就可以应用第二章已导出的一系列公式去求解。但是，这一系列公式都是在整体

8、坐标 下导出的，其中，应变矩阵 的每个元素都是各结点形函数 对整体坐标 进行重积分，而任意四边形等参元的形函数 又是针对局部坐标的，因此需要对 和 进行坐标变换。这样(zhyng)，就引出了坐标变换矩阵和变换行列式。 和 的偏导数；单元(dnyun)刚度矩阵 的每个元素又是各结点形函数 对整体坐标 和 的偏导数的乘积，再对二坐标变换矩阵及变换行列式第8页/共65页第九页，共65页。),(, ),(yxvvyxuu),(),(, ),(),(141vvNvuuNuniiiiii),(),(, ),(),(141yyNyxxNxniiiiii 设任意四边形在整体坐标下的位移(wiy)插值函数式为而

9、该单元在局部(jb)坐标系下的位移插值函数式(3-3)可以写成： 这两种形式的位移插值函数式通过坐标变换式(3-5)联系(linx)起来。为了方便，把式(3-5)写成： (3-6) (3-7) (3-8)第9页/共65页第十页，共65页。yyvxxvvyyvxxvvyyuxxuuyyuxxuu,yyxxyyxx, 根据(gnj)复合函数的求导法则，(3-6)，(3-7)，(3-8)三式之间有如下关系由式(3-9)可抽象(chuxing)出 (3-9)第10页/共65页第十一页，共65页。yxyxyxJJyxyx),(把上面(shng min)二式写成矩阵形式，得令 (3-11) (3-10)第

10、11页/共65页第十二页，共65页。JyxJ1Jyx 称 为坐标变换矩阵或雅克比矩阵，它是局部(jb)坐标的函数。因此式(3-10)变成故有 (3-13) (3-12)第12页/共65页第十三页，共65页。1JxxyyJJ11J ATedxdytBDBkdxdy 式中 称为(chn wi)坐标变换矩阵或雅克比逆阵，它也是局部坐标的函数。式中的 是坐标(zubio)变换行列式。 另外，为了把 化成对局部坐标的重积分，还需把微分(wi fn)面积 做相应的变换。 (3-14)第13页/共65页第十四页，共65页。P,ba ,ayx ,yxaa ,dyadxayx, 设任意(rny)四边形等参元1,

11、2,3,4内任一点 沿局部坐标 方向的微分矢量为 ，由于在 方向上只有 变化，而 不变，故微分矢量 在整体坐标系的 轴上的投影 分别为第14页/共65页第十五页，共65页。byx ,dybdxbyx,ba ,abFdydxdydxbbaaFyxyxab 同理，由于在 方向上只有(zhyu) 变化而 不变，故 在轴上的投影(tuyng)：两个(lin )微分矢量所构成的微小平行四边形面积 ：第15页/共65页第十六页，共65页。abFdydxFabddJddyxyxdydxdydxdydxyxyxJ而 又可以(ky)看成是在整体坐标中的微分面积，故有式中 (3-15) （3-15）第16页/共6

12、5页第十七页，共65页。 ek B, JJJ41414141),(),(),(),(iiiiiiiiiiiiyNxNyNxNJ 为了进一步阐明和计算任意四边形等参元的单元(dnyun)刚度矩阵及应力矩阵 （或应变矩阵 ），需要(xyo)把前述的 及 展成具体形式的表达式，为此，将(3-5)式代入(3-12)式，得第17页/共65页第十八页，共65页。)1)(1 (41),(iiiN,J41414141)()()()(41iiiiiiiiiiiiiiiiiiiiiiiiyyxxyyxxJ 再将对然后(rnhu)代入 中得 (3-16)分别(fnbi)求偏导数，第18页/共65页第十九页，共65页

13、。ixJiyi,4 , 3 , 2 , 1, 4413412411414141,ayaxayaxByAxiiiiiiiiiiiiiiiiiiii 式(3-16)表明，只有给定整体坐标下的单元四个结点的坐标值 和 就完全由单元的局部坐标来决定了，而且它的每个元素(yun s)都是 和 的线性函数，令常数项分别为第19页/共65页第二十页，共65页。BaAaBaAaJ432141)()()(16134213241BaAaAaBaaaaaJ则(3-16)式可写成由式(3-16)可以(ky)得出雅克比行列式 (3-16) (3-17)第20页/共65页第二十一页，共65页。,JAaAaBaBaBaAa

14、AaBaaaaaAaAaBaBaJJ13243421324113241)()()()()(4)()(411J, 它也是 的线性函数。由(3-16)和(3-17)式可以(ky)直接写成的逆阵：式中的每个元素(yun s)变成 (3-18)的较复杂(fz)的函数。第21页/共65页第二十二页，共65页。1J0J0180 由式(3-18)看出，为了(wi le)确保 的存在，必须要求变换行列式 ，这个条件的实质是，要求任意四边形等参元在整体坐标下的形状必须是凸的四边形，而不能有一个内角等于或大于 。否则，在单元上将得不到整体坐标与局部坐标之间的一一对应的变换关系，而使计算方法失效。因此，所谓任意四边

15、形等参元，其任意性还是(hi shi)有一定限度的，要求四边形的任意两对边不能通过适当的延伸而在单元内出现交点。通常，在实际有限元计算中，为了尽量使其形状接近于正方形比较好，但可以大小不一样。第22页/共65页第二十三页，共65页。JJ1J ek Bdydx及 表达式中的有关 及 都换成局部坐标的函数表达式。此时，任意四边形等参元的一切计算都可以(ky)立足在局部坐标系下进行了。 根据上述已求得的 ， 及 等函数(hnsh)表达式，就可以将首先(shuxin)，由式(3-13)引出：第23页/共65页第二十四页，共65页。iiiiNNJyNxN1iiiiNNAaAaBaBaGNNAaAaBaB

16、aBaAaAaBaaaaa1324132434213241)()()()()()()(4第24页/共65页第二十五页，共65页。iiiiiiNAaGNAaGyNNBaGNBaGxN)()()()(1324)()()(434213241BaAaAaBaaaaaG所以(suy)式中 (3-19)第25页/共65页第二十六页，共65页。 B B, B ATedxdytBDBk ek, 1111ddJtBDBkTe ek 然后(rnhu)，将式(3-19)代入 中，就把 的各元素化成 的函数；再将式(3-17)代入式(3-15)，并将式(3-15) 及 代入，就把的重积分(jfn)，其被积函数(hns

17、h) 应该指出， 中的每个元素都含有对 和 的重积分，尽管其积分区域变得十分简单，而其被积函数都比较复杂，需要采用数值积分（通常是采用高斯求积法），由于任意四边形等参元的应力 是 和 的函数，因此在求解单元应力时，必须指明是求哪一点的应力，而且各单元之间的应力是不连续的。的每个元素化成对局部坐标的函数式。都是 和 的复杂函数，对于各单元的应力 也可以化成是 和第26页/共65页第二十七页，共65页。 四节点任意四边形等参元尽管比矩形单元好，比三角形单元的精度高，但是它对结构的曲线边界仍然要以许多小直线段去逐渐逼近，计算精度仍不够理想，为了(wi le)进一步提高计算精度，可在四节点任意四边形等

18、参元的基础上，增加结点个数，选用高幂次结构模式的等参元。一般常用的是八结点曲边四边形等参元。一位移插值函数(hnsh)及坐标变换第27页/共65页第二十八页，共65页。xoy8 , 2 , 1,iyxii1o 左图是8结点平面等参元在整体坐标 下的实际形状，它除了(ch le)四个角点1，2，3，4之外，又在每边中点选一个结点5或6，7，8，各结点的整体坐标值为 。类似于4结点任意四边形等参元，这种子单元映射到局部坐标系 上，就变成边长为2的8结点正方形母单元。第28页/共65页第二十九页，共65页。2162152141321211109282726524321aaaaaaaavaaaaaaa

19、au821,aaai)8 , 2 , 1(ii 现在，我们首先(shuxin)考虑在局部坐标系下的8结点正方形母单元。由于它的8个结点共有16个位移分量，故必须选择局部坐标的双二次多项式，做为它的位移模式。式中 是由单元(dnyun)8个结点的局部坐标值 和来决定(judng)的待定常数。 (3-20)第29页/共65页第三十页，共65页。iivu ,ii,1621,aaa8181),(,),(iiiiiivNvuNu),(iN)8 , 3 , 2 , 1(ii 可用类似于以前用过的方法，将8结点(ji din)的局部坐标值代入式(3-20)，得到以 为未知数，以 及 为已知数的16个联立方程

20、组，求解这个方程组即得 的表达式。然后再回代到(3-20)中，经整理 式中 是第 (3-21)个结点(ji din)的形函数。1621,aaa第30页/共65页第三十一页，共65页。iNiNiN 应该指出，用上式求解联立方程组的方法来导出式(3-21)，比较麻烦，特别(tbi)是当待定系数比较多时，欲导出形函数 的显式是非常烦琐的。为了简化求导 的过程，我们可以利用 的特点来决定各结点的形函数显式。iNiNiNvu ,iNvu ,iN由第二章知，各种( zhn)单元的形函数都有两个(lin )特点时， 或 是 的二次函数；当固定 时， 或 又是的二项函数。 1) 是形如单元位移模式的同幂次多项

21、式，对于8结点曲边形等参元，它的 一定是局部坐标 和 的双二次多项式，固定 第31页/共65页第三十二页，共65页。iNi)(ijj0),(jjiN4 , 3 , 2 , 1),(iN8 , 3 , 275,43,328 , 3 , 201,01,01) 1)(1)(1(8 , 3 , 2) 1)(1)(1() 1)(1)(1(),(1kN2) 在 结点(ji din)的值为1，在其余结点(ji din) ， 根据形函数的这两个特点，先求8结点(ji din)曲边四边形等参元四个角点 的形函数。对于角点1，根据第2)个特点， 在结点(ji din) 的值全为零，而直线 通过 等七个结点(ji

22、din)，故这三条直线的方程分别为 同时，由 组成的函数，在 结点上的值恰好都等于零；另外，根据(gnj)第1)个特点，函数也恰好是双二次函数。因此，可设第32页/共65页第三十三页，共65页。k111) 1,1(1N41k) 1)(1)(1 (41),(1N) 1)(1)(1 (41),() 1)(1)(1 (41),() 1)(1)(1 (41),(432NNN式中 是待定常数。将1结点的局部(jb)坐标值 和 代入上式，并考虑 这一特点，可求得同理可得 (3-22)，于是(ysh)第33页/共65页第三十四页，共65页。8 , 7 , 6 , 5),(5N8 , 7 , 6 , 4 ,

23、3 , 2 , 143,32,2101,01,01) 1)(1(28 , 7 , 6 , 4 , 3 , 2 , 1) 1)(1(2iN) 1)(1(),(215kN1k101)0,1(5N211k)1)(1 (21),(25N 再来求8结点曲边四边形单元各边中点 的形函数。对于(duy)5结点，根据第2)个特点， 在 结点的值全应等于零，而直线 通过上述七点，这三条直线的方程分别为 于是，函数 在结点 上的值恰好都等于零；另外，函数 也符合 的第1)个特点(tdin)。因此，可设 式中 也是待定常数(chngsh)。将5结点的局部坐标值 和 代入上式，并考虑 这一特点，可求得 ，于是第34页

24、/共65页第三十五页，共65页。)1)(1 (21),()1)(1 (21),()1)(1 (21),(282726NNNii00,8 , 2 , 1)1)(1)(1 (21)1)(1)(1 (21) 1)(1)(1 (41),(22022202220000iNiiiiiii 同理可得 (3-23)，则可将(3-22)和(3-23)式合并(hbng)成一个通式 (3-24)如令第35页/共65页第三十六页，共65页。iNuvuvuv 从式(3-24)看出， 是双二次函数 ，从而使单元任一点(y din)的位移插值函数 和 式(3-21)也是双二次函数：单元每一条边上的和 是 或 的二次函数，它

25、完全由边上的3结点的函数值唯一(wi y)决定，而且在相邻两单元的公共边上，其三个结点有相同的函数值。因此，这种单元的位移插值函数 和 ，以及形函数能完全满足变形连续性条件和相容条件，结构变形后，各单元之间和每个单元都不能出现开裂和重叠现象。 由于 是双二次函数，它对结构线性函数都是精确成立的，故可用类似于4结点任意(rny)四边形等参元曾用过的方法，直接写出单元内任一点的局部坐标的线性插值式。iN8181),(,),(iiiiiiNN (3-25)第36页/共65页第三十七页，共65页。81),(1iiN8181),(,),(iiiiiiyNyxNxiN且有 以及该点的整体(zhngt)坐标

26、与局部坐标之间的变换式由于 的相容(xin rn)性，式(3-27)也满足相容(xin rn)条件。(3-26)(3-27)第37页/共65页第三十八页，共65页。J1JJxyx,y81818181,iiiiiiiiiiiiyNyyNyxNxxNx 8结点曲边四边形等参元的坐标变换矩阵 和其逆阵 以及变换行列式 ，仍可以(ky)采用(3-12)、(3-14)和(3-15)，但这三式中及项，需做如下(rxi)改变 (3-28)所共同(gngtng)包含的第38页/共65页第三十九页，共65页。iNiN8 , 2 , 1)1 ()1 (21)1 ()1 ()2()1 (41)1 ()1 (21)1

27、 ()1 ()2()1 (4122220220202222022020iNNiiiiiiiiiiiiiiiiiiii而式(3-28)中的和，可将式(3-24)分别(fnbi)对,求得： (3-29)的偏微分而第39页/共65页第四十页，共65页。 最后(zuhu)，由式(3-13)可直接引出8 , 2 , 1,1iNNJyNxNiiii (3-30)第40页/共65页第四十一页，共65页。J1JJxNiyNi ek 将式(3-29)代入式(3-28)，再将式(3-30)分别代入式(3-12)，(3-14)及(3-15)中，即可求得8结点(ji din)曲边四边形等参元的坐标变换矩阵和其逆阵 ，

28、以及坐标变换行列式 的具体(jt)表达式。这些表达式都是局部坐标 和 的函数。最后，再把式(3-29)代入式(3-30)中，就把 和 转化成为局部坐标的函数，这对于立足于局部坐标去计算单元应力 和单元刚度(n d)阵 是非常方便的。 应该指出，为了保证8结点曲边四边形等参元的坐标变换能顺利进行，对整体坐标下8个结点位置 的配置必须做一定的限制。既不能使单元太偏斜，又要求任意两条对边经过适当延伸也不能在单元内出现交点。通常，在划分单元网格时，尽量把每个单元配置成接近正方形。)8 , 2 , 1,(iyxii第41页/共65页第四十二页，共65页。 将式(3-21)代入平面(pngmin)问题的几

29、何方程中，便得出8结点曲边四边形单元的应变分量计算式 81821iiieeBBBBBxvyuyvxu (3-31)二单元(dnyun)分析第42页/共65页第四十三页，共65页。iB B8 , 2 , 100ixNyNyNxNBiiiii e 8 , 2 , 1821ivuTiiTiTTTTe式中 单元应变(yngbin)矩阵 的第i个子矩阵单元(dnyun)结点位移列阵第43页/共65页第四十四页，共65页。 eTeexyyxSSSSBD821 S iS8结点(ji din)曲边四边形单元的应力表达式式中 和应力矩阵及其子矩阵，对于平面(pngmin)应变问题(3-32) 8 , 2 , 1

30、212112ixNyNyNxNyNxNEBDSiiiiiiii第44页/共65页第四十五页，共65页。 应用虚功方程(fngchng)，仍可以导出这种单元的刚度矩阵。 1111ddJtBDBdxdytBDBkTATe (3-33) ijk 8, 2 , 1,1111 jiddJtBDBkjTiij把式(3-33)写成分块矩阵，可分成(fn chn)88个子矩阵，每个子矩阵都是22阶矩阵 ： (3-34)第45页/共65页第四十六页，共65页。J1JJJ ijk 将式(3-29)代入式(3-28)中，再将式(3-28)代入式(3-12)、(3-14)和(3-15)即可求得8结点(ji din)曲

31、边四边形等参元的 , 及 ；然后将式(3-29)代入式(3-30)中，再将式(3-30)及已求得 的代入式(3-34)中，经过局部坐标的积分，可得到 ，其中被积函数 xNxNyNyNyNxNxNyNxNyNyNxNyNyNxNxNEBDBjijijijijijijijijTi2121212112 (3-35)第46页/共65页第四十七页，共65页。 ek ek K eF F FK e e e或 的具体数值；最后(zuhu)把各单元的 组集结构整体刚度矩阵，把各单元的结点力列阵(li zhn) 组集成整个结构的结点力列阵(li zhn) ，并组成结构(jigu)刚度方程 ，再考虑结构(jigu)

32、约束条件，即可求解出离散结构(jigu)上各结点的位移分量列阵 和各单元的结点位移分量列阵 。 求得了各单元的 后，再把 和已计算过的式(3-30)一起代入式(3-32)中，且要给出各单元需要求应力的局部坐标值，就可以求得各单元内需要求应力那些点的各应力分量值。第47页/共65页第四十八页，共65页。 应该指出，无论是8结点曲边四边形等参元或者是4 结点任意四边形等参元，其单元内每一点的应力 都是不相同的，且 是点的局部坐标复杂函数，因而(yn r)这种单元的精度比较高。但是，在相邻单元公共边上的应力函数仍然是不连续的（其位移是连续的）。因此，通常是求单元各高斯积分点处的应力。（见3.4） 8

33、结点(ji din)曲边四边形等参元基本和4结点(ji din)任意四边形等参元的等效结点(ji din)力计算方法相同。现以8结点(ji din)四边形等参元为例，讨论如何把单元上的载荷化成等效结点(ji din)力。三等效结点(ji din)力计算第48页/共65页第四十九页，共65页。 TyxPPP P 8 , 2 , 100iPPNNPNFFFyxiicieiyixeiciNiNiNciN P 1. 集中力：设单元(dnyun)上任意点c受有集中载荷，则 被移置到单元各有关受载结点上的等效结点力，可按第2章2-8节中讲过的方法(fngf)直接写成 (3-36)式中 代表 在集中力作用点

34、处c的取值，把c点的局部坐标值代入形函数(hnsh) ，再去计算 。实际计算时，应尽量把集中力作用点取为结点，而把 直接加到该结点上。第49页/共65页第五十页，共65页。 TyxQQQ Q 8 , 2 , 11111 iddJtQQNFFFyxieiyixei 2. 体积力：设单元(dnyun)上的体力为 ， 移置到单元(dnyun)各结点上的等效结点力。按式(2-38)可写成式中t是单元(dnyun)厚度。 (3-37)第50页/共65页第五十一页，共65页。 Tyxqqq dstqqNFFFyxieiyixeids3. 面力：设单元(dnyun)的某边界上受有面力上有关(yugun)结点

35、的等效结点力按式(2-39)可写成式中 是单元(dnyun)作用有面力的边界域； 是在 边界域的一个微弧长；i为受面力边界上的结点号码。，这条边界(3-38)第51页/共65页第五十二页，共65页。xqyqnqtqnqtqnqtq 应指出，式(3-38)中所给出的面力分量 和 ，实用时不太方便。在实际结构上往往是给出沿单元曲线边界的法线和切线方向(fngxing)的面力 和 ，故需对式(3-38)做适当修改。 规定：法向面力 以沿边界(binji)曲线的外法线方向为负，切向面力 以沿单元边界(binji)逆时针方向前进者为正。图中指出的 和都是正的。第52页/共65页第五十三页，共65页。48

36、3nqtqnq090dsdxdsdy)(cos;)(sindsdxqdsdyqqqqdsdyqdsdxqqqqntntyntntx)2(sin)(sin)2(cos)(cos 设图中8结点平面(pngmin)等参元的 边界上受有面力 及 ，且与x轴的夹角为 。则 与x轴的夹角就是 。由图知而第53页/共65页第五十四页，共65页。xqyq 8 , 4 , 3idxqdyqdyqdxqNtdstdsdxqdsdyqdsdyqdsdxqNFFFntntintntieiyixei),(, ),(),(81yyxxNxiiidydydydxdxdx,把 和 代入(3-38)式得由(3-21)知，故x或

37、y对和的重积分(jfn)为(3-39)第54页/共65页第五十五页，共65页。4831dydydxdx,dxdy 1111dxqyqyqxqNtFFFntntieiyixei 对于(duy)图中等参元的 边界，其局部坐标 ， 是变化的，因此x或y对局部坐标的全微分应为 将 和 代入式(3-39)得 (3-40)xNi,y1由于(yuy)，及都是的复杂(fz)函数（因的积分也要用数值积分法（常用高斯求积法来求解）。），故式(3-40)第55页/共65页第五十六页，共65页。 TT0110 dydxtDBHHHTieiyixei0iB DddJTtNNEdydxTtNNEHHHyixiyixiei

38、yixei 1111,11 4. 温度荷载(hzi)考虑温度变化产生的初应变则任意结点(ji din)上的等效结点(ji din)力是将 和 代入上式，可以(ky)写成第56页/共65页第五十七页，共65页。 TeeTEDBDBDBDBD01118210此时，计算应力(yngl)的（3-32）式改写为第57页/共65页第五十八页，共65页。 111111111111),(,),(,)(dddfddfdfff 在前二节的刚度矩阵和等效结点(ji din)力的计算公式中，都需要做如下形式的积分运算 显然，被积函数 一般是很复杂的，往往不能得出(d ch)它的显式。因此，在有限单元法的计算中都用数值

39、积分。我们在单元内选出某些点（称为积分点），算出被积函数 在这些积分点处的函数值，然后用一些加权系数乘上这些函数值，再求出总和作为近似的积分值。高斯求积法就是数值积分法中具有较高精度的方法。第58页/共65页第五十九页，共65页。niiifHdf111)()(iiH一维高斯(o s)求积公式式中n是积分(jfn)点的数目， 是积分(jfn)点i的局部坐标， 是加权系数。下表是n4的数值第59页/共65页第六十页，共65页。iiH31539598高斯(o s)求积公式中的积分点坐标与加权系数n213040.861,136,311,5940.339,981,043,5840.347,854,845

40、,1370.652,145,154,862第60页/共65页第六十一页，共65页。)(f)(f21mn)(f 积分点数目n的选取与被积函数 有关，当 是m次多项式时，则取 。当 不是多项式时，则需通过一些计算来判断(pndun)选取适当的n值；n不能取得过大，因为计算工作量将会随着积分点的增多而急剧地增加。逐次地利用(lyng)上式。则有下列公式 nimjlkkjikjinimjjijifHHHdddffHHddf111111111111111),(),(),(),(它们可以(ky)应用于二维和三维问题。第61页/共65页第六十二页，共65页。 具体的说，对于平面等参元，求积点的个数可选为n2

41、；对于空间(kngjin)等参元，求积点个数可选为n3。阶数n可按下表选取。4结点8结点20，21结点二维三维n=2n=3n=2n=3 应该指出，高斯(o s)求积法仍有一定误差，但这一误差在某种程度上可起到一定的“补偿”作用。因为，用有限单元法求解结构的应力与变形问题时，往往使结构变得“过刚”，所得有限元法的近似解的变形能总小于精确解的变形能，因而使所求得的应变及应力在总刚上来说是偏小的；高斯(o s)积分的误差对有限元法近似解一般可以起到适当增大变形能的作用，所以说起到一定的“补偿”作用。第62页/共65页第六十三页，共65页。 ek66 ek1616 k1234567891012345

42、第三次作业第三次作业(zuy)： 6已知图示所示的梁已知图示所示的梁,三角形单元三角形单元 和八节点等参元和八节点等参元。求7求图示5节点(ji din)单元的形函数。第63页/共65页第六十四页，共65页。xN1yN2)21,21(xy1(40,50)2(5,40)3(10,10)4(30,20)5678Q 8图示为二次四边形单元，试计算 和 在自然坐标为 的点Q的数值（提示：因为单元的边是直线，可用4个节点定义(dngy)单元的几何形状）。 9二维8节点等参元，在x, y坐标中单元各边与坐标轴x, y平行，边长为a, b，确定下列载荷情况下的节点荷载(hzi)。 （a）在x正方向有一分布荷载(hzi)作用在=1的边上，荷载(hzi)在= 1为0；在=1为q0，呈直线变化； （b）在=1的边上作用有均布荷载(hzi)q0，方向压向单元； （c）在y正方向上作用有均匀的体积力b。第64页/共65页第六十五页，共65页。