04-04 高斯积分法及其应用.docVIP

§4-4 高斯积分法及其应用 由§4-3知，在计算空间等参数单元的载荷列阵及刚度矩阵时，需用到如下形式的定积分： ； 其中被积分函数f(ξ，η，ζ)一般是很复杂的，即使能够得出它的显式，其积分也是很繁的。因此，一般用数值积分来代替函数的定积分。 数值积分：在积分区域内按一定规则选出一些点，称为积分点，算出被积函数f(ξ，η，ζ)在这些积分点处的值，然后再乘以相应的加权系数并求和，作为近似的积分值。 数值积分的方法有多种，其中高斯积分法可以用相同的积分点数达到较高的精度，或者说用较少的积分数达到同样的精度。 一、高斯积分法 1．一维积分的高斯公式 一维积分的高斯公式 (4-47) 其中f(ξi)是被积函数在积分点ξi处的数值，Hi为加数系数，n为积分点数目。 可以证明， 对于n个积分点，只要选取适当的加数系数及积分点位置，能够使(4-47)式在被积分函数为不超过(2n-1)次多项式时精确成立。 由于多数函数可表示成多项式形式，这种积分适应于大多数函数。 例如， n=1时 (a) 不论f(ξ)的次数是0还是1，只需取H１=2，ξ１＝０，上式均是精确成立的。因为 (b) (c) 当n=2时，能保证(4-47)式精确成立所允许的多项式的最高次数是3，此时，f(ξ)的通式为 (d) 其精确积分为 (e) 数值积分为 (f) 为了在C0～C3取任意值(包括取零值在内)时公式(f)是精确的，显然应有 ， ， 所以，应取 同样，对于不超过五次的多项式，只要取 n=3 即可保证得到精确的积分值。 对于n=1到n=6时，高斯求积公式中积分点坐标ξi及加权系 Hi的数值列于表4-1。 表3-1 高斯求积公式积分点坐标与加权系数 点数n 坐标ξi 加权系数Hi 1 2 3 4 5 6 0. ±0.57735 02691 89626 ±0.77459 66692 41483 1.00000 00000 00000 ±0.86113 63115 94053 ±0.33998 10435 84856 ±0.90617 98459 38664 ±0.53846 93101 05683 0.00000 00000 00000 ±0.93246 95142 03152 ±0.66120 93864 66256 ±0.23861 91860 83197 2. 1.00000 00000 00000 0.55555 55555 55555 0.88888 88888 88888 0.34785 48451 37454 0.65214 51548 62546 0.23692 68850 56189 0.47862 86704 99366 0.56888 88888 88889 0.17132 44923 79179 0.36076 15730 48139 0.46791 39345 72691 2．二维积分的高斯公式 以一维高斯积分公式(4-47)为基础，极易导出二维及三维公式。求二维重积分 的数值时，可以先对ξ进行积分，这时将η当作常量，于是由(4-47)式得到 (g) 再对η进行积分，得出 (h) 将式(g)代入，即得 或改写为 (4-48) 这就是二维的高斯积分公式。 3．三维积分的高斯公式 同样，可以求得三维高斯积分公式： (4-49) 式(4-48)及(4-49)中的n，m，l是分别关于变量ξ，η，ζ的积分点数目。 各个维数上的积分点数目由各个自变量在被积函数中可能出现的最高次数分别决定，一般并不要求相同。但为应用方便，常常在各个方向取相同的积分数，即统一为最高值 (4-48a) (4-49a) 由前面的推导可见，当在每个方向取n个积分点时，只要多项式被积函数中自变量的次数m≤2n-1，则用高斯求积公式求得的积分值是完全精确的。 反过来，对于m次多项式的被积函数，为了积分值完全精