高斯積分點以及有限元中應用學習教案

1、會計學1高斯積分點以及有限元中應用高斯積分點以及有限元中應用第一頁，編輯于星期三：七點 二十八分。 在計算空間等參數單元的載荷列陣及剛度矩陣時，需用到如下形式的定積分： 其中被積分函數f(，)一般是很復雜的，即使能夠得出它的顯式，其積分也是很繁的。因此，一般用數值積分來代替函數的定積分。dddf 111111),(ddf 1111),(第1頁/共19頁第二頁，編輯于星期三：七點 二十八分。 數值積分數值積分：在積分區(qū)域內按一定規(guī)則選出一些點，稱為積分點，算出被積函數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-節(jié)點載荷向量；K-總體剛度矩陣；U-節(jié)點位移向量 在引入邊界條件以后，解上述方程組，就可以得到節(jié)點位移向量U.這是求解結構靜力學方程組所得到的第一組解，它是最精確的。 得到節(jié)點的位移解后，下面是求取應變解和應力解。與位移解不同，它們并不是直接在節(jié)點上獲得，而是首先在積分點上獲得的。 第14頁/共19頁第十五頁，編輯

8、于星期三：七點 二十八分。有限元分析主要步驟 所謂積分點是指，在對單元建立方程時，例如剛度矩陣是需要通過積分而得到的，而積分時為了能夠方便計算，大多數有限元軟件采用了所謂高斯積分的方式，即在單元內分布一些高斯點 這樣，有限元軟件會首先獲得這些高斯點的應力和應變，其方法如下：l在高斯積分點上，依據幾何方程:=BUl計算出高斯積分點上的應變:l然后基于虎克定律及幾何方程推導的結果來計算高斯積分點的應力。:=DBU第15頁/共19頁第十六頁，編輯于星期三：七點 二十八分。有限元分析主要步驟 可見，在應變和應力計算方面，高斯積分點的應變和應力是最最準確的。 利用特定單元的形函數以及高斯點的應力，應變值

9、，將這些值外推到該單元的節(jié)點上，就得到了單元上節(jié)點的應力應變值。 顯然，不同的單元會共用一些節(jié)點，而從不同單元內的積分點外推到這些公共節(jié)點的應變值和應力值一般不相同,將一個公共節(jié)點的多個應力進行平均，以代表該節(jié)點的應力值。第16頁/共19頁第十七頁，編輯于星期三：七點 二十八分。有限元分析主要步驟總之，求解節(jié)點應力的步驟是：（1）根據總體方程，得到節(jié)點的位移解。（2）根據幾何方程，得到單元高斯點的應變解。（3）根據物理方程，得到單元高斯點的應力解。（4）在某一個單元內，基于形函數，將高斯點的應力外 推到該單元的所有節(jié)點。（5）對于某一個公共節(jié)點，將該節(jié)點關聯的所有單元所推出的該節(jié)點的應力解進行平均，最終得到該節(jié)點的應力解。第17頁/共19頁第十八頁，編輯于星期三：七點 二十八分。積分點與節(jié)點的關系 我們需要對應變在單元內的面積上進行積分時，因為節(jié)點的應力、位移顯然與x,y無關，我們只需要考慮對形函數積分。 采用Gauss-Legendre多項式計算積分時，我們只需要計算根據特定積分點的值（在自然坐標系下是固定的，可以查手冊，這些點也叫高