基于Fourier譜元法的拱壩結構靜動力分析探究

唐波

簡介：將譜方法和有限元理論結合形成的譜元法，同時具有譜方法任意階收斂和有限元可操作性的優點。將Fourier積分函數運用于譜元法對拱壩結構計算中，具有快速、高效、精確的求解優勢，對提高數值計算精度和譜元法在工程中的實用性，具有較高的參考價值。

關鍵詞：拱壩結構靜動力分析、譜元法、Fourier積分函數

Research on static and dynamic analysis of dam structure

based on Fourier spectral element method

Tang Bo

School of Civil Engineering, Tong ji University, Shang Hai, China, 200092

Abstract：The spectral element method originated from the finite element method (FEM) and the spectral method has the advantages of excellent operability and convergence of arbitrary order. The dam structural calculation using spectral element method with Fourier integral function can solve equations fast, efficiently and accurately. Therefore, it has a large value in improving the numerical accuracy and the practicality of spectral method in practical engineering.

Keywards：Static and dynamic analysis of dam structure, spectral element method, Fourier integral function

1、引言

結構靜、動力分析的高速、有效和準確性，一直以來都是研究人員關注的問題，在結構計算理論也不斷發展的基礎上，新的計算方法也不斷被提出，例如數值差分法、有限元方法等。本文在前人的基礎上，利用譜方法的優點，對結構各種荷載激勵下的靜、動力計算分析方法做出了初步探究。

2、基本理論

我們知道，一般靜、動力問題的基本理論方程為：

幾何方程： (1)

為應變向量，為微分算子，為形函數，為幾何矩陣，為節點位移向量，為拱壩結構位移矩陣。

物理方程：(2)

為結構應力向量，為剛度矩陣。

力平衡方程：(3)

是微分算子，是該點在單位體積的體積力向量，和分別為質量密度和阻尼系數。

力的邊界條件：（4-1）

為邊界上已知彈性體單位面積上作用的面積力量；為邊界外法線方向余弦形成的列向量；

幾何邊界條件：（4-2）

為邊界上已知的位移向量。

3、基本動力方程

由有限元法可知，由達朗貝爾原理可知結構的基本動力方程為：

（5）

單元剛度矩陣：（5-1）

單元質量矩陣： （5-2）

單元阻尼矩陣：（5-3）

矩陣形式表達為：

上式中，為單元荷載向量。

則對整個結構可知：（6）

分別為結構剛度矩陣，阻尼矩陣，質量矩陣，荷載矩陣。

4、利用Fourier譜元法形成計算形函數矩陣

對于任意函數，采用fourier逼近則有：

由最佳平方逼近可知得：

對于在離散點 給定值的函數，可以類似的求出在離散點的離散Fourier系數和，即得到如下表達式：

（）（7）

利用三角變換可得：（8）

其中：

在實際運用中，采用如下點集： ，則代入上式可得：

（9）

其中：

則可以得到：

對于函數導數的計算可以通過對的導數計算，配置點處的導數為：

（10）

對三維函數，在標準的正方體單元內，分別定義三個方向上的節點系，其中：

則有：

則的傅里葉譜展開式為：

（11）

其中：，，

顯然：

若共i個點，則（13）式可表示為：

得到位移形函數再代入下面的計算中。

5、剛度矩陣、質量（阻尼）矩陣、荷載向量推導

類似有限元方法，推導在變換矩陣和矩陣

同樣按照有限元的推導方法，得到：

由得到單元剛度分塊矩陣

其中： （12）

的表達形式類似于有限元法的表達形式，不過形函數形式如前所推導。

由得到單元質量分塊矩陣，

其中： （13）

由得到荷載列向量

其中: 體力荷載（14-1）

面力荷載（14-2）

得到積分形式之后，我們可以參考有限元等參單元和高斯積分將積分形式化為代數乘積表達式，從而得到相應矩陣。

6、編程進行結構靜、動力計算分析

拱壩結構受力分析計算一般計算按照如下步驟進行，模型單元劃分——→設定基本參數及規定邊界條件和初始條件——→讀取荷載——→計算單元應力、應變——→結果顯示

7、結論

以上的理論推導，筆者認為：譜元法運用于拱壩結構的靜、動力分析具有精度高、快速收斂的特點，在實際工程中是可行的，并且由于其獨特的優勢，將來的運用會越來越廣泛。當然，由于譜元法還不夠成熟，成熟的運用于實際工程還有很大困難，需要更多的科研人員努力。

參考文獻（References）：

[1]: 趙廷剛.若干發展方程的譜方法和譜元法:[博士學位論文].上海：上海大學數學系，2007

[2]：Usik Lee, Joohong Kim, and Andrew Y. T. Leung，The Spectral Element Method in Structural Dynamics. The Shock and Vibration Digest,2000 Vol.32:415～465

[3]：M. Krawczuk, M. Palaczb, W. Ostachowiczb, The dynamic analysis of a cracked Timoshenko beam by the spectral element method. Journal of Sound and Vibration,2003:1139～1153

[4]：王海波，李德玉.拱壩抗震設計理論與實踐. 北京：中國水利水電出版社，2006

[5]：朱伯芳.有限單元法原理與應用. 北京：中國水利水電出版社，2009

[6]：向新民.譜方法的數值分析.北京：科學出版社，2000

[7]：John P.Boyd . Chebyshev and Fourier Spectral Methods . University of Michigan , 2000,19-80

[8]：王勖成.有限單元法. 北京：清華大學出版社，2003

注：文章內所有公式及圖表請以PDF形式查看。