高斯积分法及其应用.doc

《高斯积分法及其应用.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、§4-4高斯积分法及其应用l由§4-3知，在计算空间等参数单元的载荷列阵及刚度矩阵时，需用到如下形式的定积分：；其中被积分函数f(ξ，η，ζ)一般是很复杂的，即使能够得出它的显式，其积分也是很繁的。因此，一般用数值积分来代替函数的定积分。l数值积分：在积分区域内按一定规则选出一些点，称为积分点，算出被积函数f(ξ，η，ζ)在这些积分点处的值，然后再乘以相应的加权系数并求和，作为近似的积分值。l数值积分的方法有多种，其中高斯积分法可以用相同的积分点数达到较高的精度，或者说用较少的积分数达到同样的

2、精度。一、高斯积分法1．一维积分的高斯公式l一维积分的高斯公式(4-47)其中f(ξi)是被积函数在积分点ξi处的数值，Hi为加数系数，n为积分点数目。l可以证明，²对于n个积分点，只要选取适当的加数系数及积分点位置，能够使(4-47)式在被积分函数为不超过(2n-1)次多项式时精确成立。²由于多数函数可表示成多项式形式，这种积分适应于大多数函数。l例如，²n=1时(a)不论f(ξ)的次数是0还是1，只需取H１=2，ξ１＝０，上式均是精确成立的。因为(b)(c)²当n=2时，能保证(4-47)式精

3、确成立所允许的多项式的最高次数是3，此时，f(ξ)的通式为(d)其精确积分为(e)数值积分为(f)为了在C0～C3取任意值(包括取零值在内)时公式(f)是精确的，显然应有，，所以，应取²同样，对于不超过五次的多项式，只要取n=3即可保证得到精确的积分值。²对于n=1到n=6时，高斯求积公式中积分点坐标ξi及加权系Hi的数值列于表4-1。表3-1高斯求积公式积分点坐标与加权系数点数n坐标ξi加权系数Hi1234560.±0.577350269189626±0.7745966692414831.000

4、000000000000±0.861136311594053±0.339981043584856±0.906179845938664±0.5384693101056830.000000000000000±0.932469514203152±0.661209386466256±0.2386191860831972.1.0000000000000000.5555555555555550.8888888888888880.3478548451374540.6521451548625460.2369268

5、850561890.4786286704993660.5688888888888890.1713244923791790.3607615730481390.4679139345726912．二维积分的高斯公式l以一维高斯积分公式(4-47)为基础，极易导出二维及三维公式。求二维重积分的数值时，可以先对ξ进行积分，这时将η当作常量，于是由(4-47)式得到(g)再对η进行积分，得出(h)将式(g)代入，即得或改写为(4-48)这就是二维的高斯积分公式。3．三维积分的高斯公式同样，可以求得三维高斯积分

6、公式：(4-49)式(4-48)及(4-49)中的n，m，l是分别关于变量ξ，η，ζ的积分点数目。l各个维数上的积分点数目由各个自变量在被积函数中可能出现的最高次数分别决定，一般并不要求相同。但为应用方便，常常在各个方向取相同的积分数，即统一为最高值(4-48a)(4-49a)由前面的推导可见，当在每个方向取n个积分点时，只要多项式被积函数中自变量的次数m≤2n-1，则用高斯求积公式求得的积分值是完全精确的。反过来，对于m次多项式的被积函数，为了积分值完全精确，积分点的数目必须取。二、合适积分点数

7、的确定在空间等参数单元分析中，计算节点载荷及刚度矩阵时要用到高斯积分公式，下面以20节点的等参数单元（二次）为例，来探讨一下相应公式中积分点的数目。这直接关系到收敛性、计算精度和计算时长。1．分布体力★分布体力的等效节点载荷列阵计算见公式(4-35)，即(4-35)l由式(4-16)可见，形函数Ni对每个局部坐标来说，一般都是二次式，因而在元素中，每个局部坐标都可能以二次幂出现，同时，由(4-15)式又知，在整体坐标的表达式中，每个局部坐标也都可能以二次幂出现。l根据(4-21)式，在中每个局部坐

8、标可能出现的次数是5。l据此，当体力分量为常量时(例如体力为均质单元自重)，(4-35)式中的被积函数对于每个局部坐标来说，可能出现的最高幂次是m=5+2=7。为了使积分值完全精确，每个坐标方向积分点的数目应为。对于三维积分，其积分点总数应为43=64个。2．分布面力★一般分布面力的节点载荷列阵为(4-36)在这里，由于被积函数中包含有根式，不能化为多项式，因而无法精确判明所需积分点的数目，只有结合面力的分布规律，依次增加积分点数目，进行若干次试算，才能大致判明。l当分布面力为静水