高斯积分点以及有限元中应用学习教案

1、会计学1高斯积分点以及有限元中应用高斯积分点以及有限元中应用第一页，编辑于星期三：七点 二十八分。 在计算空间等参数单元的载荷列阵及刚度矩阵时，需用到如下形式的定积分： 其中被积分函数f(，)一般是很复杂的，即使能够得出它的显式，其积分也是很繁的。因此，一般用数值积分来代替函数的定积分。dddf 111111),(ddf 1111),(第1页/共19页第二页，编辑于星期三：七点 二十八分。 数值积分数值积分：在积分区域内按一定规则选出一些点，称为积分点，算出被积函数f(，)在这些积分点处的值，然后再乘以相应的加权系数并求和，作为近似的积分值。 数值积分数值积分的方法的方法有多种，其中高斯积分法

2、可以用相同的积分点数达到较高的精度，或者说用较少的积分数达到同样的精度。第2页/共19页第三页，编辑于星期三：七点 二十八分。一、一维积分的高斯公式一、一维积分的高斯公式 其中f(i)是被积函数在积分点i处的数值，Hi为加数系数，n为积分点数目。 对于n个积分点，只要选取适当的加数系数及积分点位置，能够使式在被积分函数为不超过(2n-1)次多项式时精确成立。 由于多数函数可表示成多项式形式，这种积分适应于大多数函数。 niiifHdf111)()(第3页/共19页第四页，编辑于星期三：七点 二十八分。例如，n=1时 不论f()的次数是0还是1，只需取H=2，1，上式均是精确成立的。因为 )()

3、(1111fHdfI10)(CCf101( )22(0)IfdCf第4页/共19页第五页，编辑于星期三：七点 二十八分。332210)(CCCCf2011322)(CCdfI)()()()()(323222102313212101221121CCCCHCCCCHfHfHfHIiii第5页/共19页第六页，编辑于星期三：七点 二十八分。221 HH02211HH32222211HH0322311HH2 ,269,350,577. 031210 ,000,000,000. 121 HH， ， 第6页/共19页第七页，编辑于星期三：七点 二十八分。ln n个插值结点非等距分布个插值结点非等距分布l结

4、点和积分权系数可以查表结点和积分权系数可以查表111()(niiiAfdf第7页/共19页第八页，编辑于星期三：七点 二十八分。二维积分的高斯公式二维积分的高斯公式以一维高斯积分公式为基础，导出二维及三维公式。求二维重积分 的数值时，可以先对、进行积分， 或改写成 这就是二维的高斯积分公式。ddf 1111),()(),(),(111niiifHdfmjjjHd111)()( mjnijiijfHHddf111111),(),( nimjjijifHHddf111111),(),(第8页/共19页第九页，编辑于星期三：七点 二十八分。 nimjlkkjikjifHHHdddf111111111

5、),(),( ninjjijifHHddf111111),(),( ninjnkkjikjifHHHdddf111111111),(),(第9页/共19页第十页，编辑于星期三：七点 二十八分。 由前面的推导可见，当在每个方向取由前面的推导可见，当在每个方向取n n个积分点时，个积分点时，只要多项式被积函数中自变量的次数只要多项式被积函数中自变量的次数m2n-1m2n-1，则用高斯求积，则用高斯求积公式求得的积分值是完全精确的。公式求得的积分值是完全精确的。 反过来，对于反过来，对于m m次多项式的被积函数，为了积分值完全次多项式的被积函数，为了积分值完全精确，积分点的数目必须取精确，积分点的数

6、目必须取 。第10页/共19页第十一页，编辑于星期三：七点 二十八分。l高斯积分方法预先定义了积分点和相应的加权系数，求出被积分的函数在指定积分点上的数值，加权后求和，就得到了该函数的积分。l高斯积分方法具有最高的计算精度。采用n个积分点的高斯积分可以达到2n-1阶的精度，也就是说，如果被积分的函数是2n-1次多项式，用n个积分点的高斯积分可以得到精确的积分结果。 第11页/共19页第十二页，编辑于星期三：七点 二十八分。l积分阶次的选择直接影响计算的精度和计算工作量。l积分阶次的选择必须保证积分的精度。（完全精确积分）l很多情况下，实际选取的高斯积分点数低于精确积分的要求，往往可以取得较完全

7、精确积分更好的精度。（减缩积分）第12页/共19页第十三页，编辑于星期三：七点 二十八分。线性单元完全精确积分 二次单元减缩积分第13页/共19页第十四页，编辑于星期三：七点 二十八分。有限元分析主要步骤 我们知道，经过单元方程的组装以后，结构静力学有限元方程如下F=KU 其中，F-节点载荷向量；K-总体刚度矩阵；U-节点位移向量 在引入边界条件以后，解上述方程组，就可以得到节点位移向量U.这是求解结构静力学方程组所得到的第一组解，它是最精确的。 得到节点的位移解后，下面是求取应变解和应力解。与位移解不同，它们并不是直接在节点上获得，而是首先在积分点上获得的。 第14页/共19页第十五页，编辑

8、于星期三：七点 二十八分。有限元分析主要步骤 所谓积分点是指，在对单元建立方程时，例如刚度矩阵是需要通过积分而得到的，而积分时为了能够方便计算，大多数有限元软件采用了所谓高斯积分的方式，即在单元内分布一些高斯点 这样，有限元软件会首先获得这些高斯点的应力和应变，其方法如下：l在高斯积分点上，依据几何方程:=BUl计算出高斯积分点上的应变:l然后基于虎克定律及几何方程推导的结果来计算高斯积分点的应力。:=DBU第15页/共19页第十六页，编辑于星期三：七点 二十八分。有限元分析主要步骤 可见，在应变和应力计算方面，高斯积分点的应变和应力是最最准确的。 利用特定单元的形函数以及高斯点的应力，应变值

9、，将这些值外推到该单元的节点上，就得到了单元上节点的应力应变值。 显然，不同的单元会共用一些节点，而从不同单元内的积分点外推到这些公共节点的应变值和应力值一般不相同,将一个公共节点的多个应力进行平均，以代表该节点的应力值。第16页/共19页第十七页，编辑于星期三：七点 二十八分。有限元分析主要步骤总之，求解节点应力的步骤是：（1）根据总体方程，得到节点的位移解。（2）根据几何方程，得到单元高斯点的应变解。（3）根据物理方程，得到单元高斯点的应力解。（4）在某一个单元内，基于形函数，将高斯点的应力外 推到该单元的所有节点。（5）对于某一个公共节点，将该节点关联的所有单元所推出的该节点的应力解进行平均，最终得到该节点的应力解。第17页/共19页第十八页，编辑于星期三：七点 二十八分。积分点与节点的关系 我们需要对应变在单元内的面积上进行积分时，因为节点的应力、位移显然与x,y无关，我们只需要考虑对形函数积分。 采用Gauss-Legendre多项式计算积分时，我们只需要计算根据特定积分点的值（在自然坐标系下是固定的，可以查手册，这些点也叫高