唐波

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

關(guān)鍵詞：拱壩結(jié)構(gòu)靜動力分析、譜元法、Fourier積分函數(shù)

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、引言

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

2、基本理論

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

幾何方程： (1)

為應(yīng)變向量，為微分算子，為形函數(shù)，為幾何矩陣，為節(jié)點位移向量，為拱壩結(jié)構(gòu)位移矩陣。

物理方程：(2)

為結(jié)構(gòu)應(yīng)力向量，為剛度矩陣。

力平衡方程：(3)

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

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

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

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

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

3、基本動力方程

由有限元法可知，由達(dá)朗貝爾原理可知結(jié)構(gòu)的基本動力方程為：

（5）

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

單元質(zhì)量矩陣： （5-2）

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

矩陣形式表達(dá)為：

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

則對整個結(jié)構(gòu)可知：（6）

分別為結(jié)構(gòu)剛度矩陣，阻尼矩陣，質(zhì)量矩陣，荷載矩陣。

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

對于任意函數(shù)，采用fourier逼近則有：

由最佳平方逼近可知得：

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

（）（7）

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

其中：

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

（9）

其中：

則可以得到：

對于函數(shù)導(dǎo)數(shù)的計算可以通過對的導(dǎo)數(shù)計算，配置點處的導(dǎo)數(shù)為：

（10）

對三維函數(shù)，在標(biāo)準(zhǔn)的正方體單元內(nèi)，分別定義三個方向上的節(jié)點系，其中：

則有：

則的傅里葉譜展開式為：

（11）

其中：，，

顯然：

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

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

5、剛度矩陣、質(zhì)量（阻尼）矩陣、荷載向量推導(dǎo)

類似有限元方法，推導(dǎo)在變換矩陣和矩陣

同樣按照有限元的推導(dǎo)方法，得到：

由得到單元剛度分塊矩陣

其中： （12）

的表達(dá)形式類似于有限元法的表達(dá)形式，不過形函數(shù)形式如前所推導(dǎo)。

由得到單元質(zhì)量分塊矩陣，

其中： （13）

由得到荷載列向量

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

面力荷載（14-2）

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

6、編程進(jìn)行結(jié)構(gòu)靜、動力計算分析

拱壩結(jié)構(gòu)受力分析計算一般計算按照如下步驟進(jìn)行，模型單元劃分——→設(shè)定基本參數(shù)及規(guī)定邊界條件和初始條件——→讀取荷載——→計算單元應(yīng)力、應(yīng)變——→結(jié)果顯示

7、結(jié)論

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

參考文獻(xiàn)（References）：

[1]: 趙廷剛.若干發(fā)展方程的譜方法和譜元法:[博士學(xué)位論文].上海：上海大學(xué)數(shù)學(xué)系，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]：王海波，李德玉.拱壩抗震設(shè)計理論與實踐. 北京：中國水利水電出版社，2006

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

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

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

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

注：文章內(nèi)所有公式及圖表請以PDF形式查看。