前言：想要寫出一篇令人眼前一亮的文章嗎？我們特意為您整理了5篇數(shù)值方法范文，相信會為您的寫作帶來幫助，發(fā)現(xiàn)更多的寫作思路和靈感。

數(shù)值方法范文第1篇

關(guān)鍵字：內(nèi)域；數(shù)值計算；有限元

中圖分類號：TU315.3

1.引言

雖然目前技術(shù)和計算設(shè)備發(fā)展十分迅速，計算能力不斷提高，一些大型通用有限元軟件已經(jīng)具備十分強大的分析功能，它解決了地震反應(yīng)的許多實際工程結(jié)構(gòu)分析的問題。但對于一些大型復(fù)雜體系而言，空間離散的自由度數(shù)目非常龐大，數(shù)值穩(wěn)定性的限制要求時間離散的步距也不能過大。這樣，開展結(jié)構(gòu)地震反應(yīng)分析時所需要完成的時空四維數(shù)值計算的工作量將變的很大。在工程設(shè)計中，需要分析各種工況下結(jié)構(gòu)的地震反應(yīng)行為，對比不同的設(shè)計方案并做出優(yōu)化決策，從而要求在設(shè)計期內(nèi)多次完成結(jié)構(gòu)的地震反應(yīng)計算。在數(shù)值精度的基礎(chǔ)上，保證系統(tǒng)的穩(wěn)定和提高結(jié)構(gòu)地震反應(yīng)分析方法的核心計算效率。研究高效率的大型復(fù)雜體系地震反應(yīng)數(shù)值分析方法，減少計算費用仍然是非常現(xiàn)實的考慮，它具有重要理論意義和實用價值。

地震反應(yīng)分析的數(shù)學(xué)物理模型就是波動方程。波動數(shù)值模擬包含兩個部分，一是對人工邊界的數(shù)值模擬，二是對內(nèi)域的數(shù)值模擬。這里僅討論對內(nèi)域的運動節(jié)點的數(shù)值模擬問題。

現(xiàn)有的內(nèi)域波動的數(shù)值模擬方法，按能量等效式，可劃分為三類。一類是空間有限元方法。所謂空間有限元方法是指的空間用有限元法進行離散，而使用一個直接的方式離散時間。第二類是時空有限元方法。這里是指對時間使用和對空間一樣的有限元法方式離散。時空有限元的時空域是空間有限元的空間域在時間域上增加了一個時間維度，區(qū)別只是處理時間域的先后上，空間有限元是先處理空間問題，然后處理時間問題，而時空有限元是同時處理時間和空間的問題。第三類是微分求積方法，是直接的方式對空間和時間的離散。下面分別闡述這三類方法的應(yīng)用發(fā)展過程――

2.空間有限元方法

在1956年，有限元的概念首次被Turner等人提出，最早應(yīng)用于彈性力學(xué)平面應(yīng)力問題上。1963年，Besseling、Melosh和Jones等人發(fā)現(xiàn)有限元法和基于變分原理的里茲法是等效的。有限元法在處理連續(xù)介質(zhì)問題上比普通里茲法更有優(yōu)勢。隨后幾十年，在解決復(fù)雜工程問題上，有限元法得到廣泛的應(yīng)用。

波動方程是時空耦合的，基于廣義Hamilton原理的波動有限元方法通常也是時空解耦的數(shù)值過程。傳統(tǒng)有限元方法的離散過程通常包含兩步，先進行空間離散，將微分方程轉(zhuǎn)化為常微分方程，然后對時間進行離散，即在時域?qū)ΤＮ⒎址匠踢M行數(shù)值積分。由于時空耦合的數(shù)值過程包含過多的自由度，求解這類方程在實際工程中很難實現(xiàn)，建立時空解耦的波動數(shù)值分析方法是這方面重要的工作。最直接的做法是實現(xiàn)空間及時間域的解耦，通常是只建立空間的有限元離散方法，而時間采用直接的假設(shè)，最常用的是采用逐步積分的方式進行離散。

逐步積分法簡單來講就是把最終速度和位移由它們的初始值和一個積分表達式來表示。加速度歷程的積分決定速度的變化，速度的積分決定位移的變化。換句話說，加速度控制了速度的變化，因而可以由這一步向前獲得下一時間步。解答這類問題，第一步先考慮時間步內(nèi)的加速度問題，假設(shè)加速度是如何變化的，依據(jù)加速度和位移的關(guān)系，得到關(guān)于時間步的遞推公式。所謂的Euler-Gauss法就是假設(shè)在時間間隔內(nèi)加速度為常數(shù)。而Newmark[1]法是加入系數(shù)從而可以改變初始和最終加速度的權(quán)重從而得到加速度的一種方法。Wilson- [2]方法是假定在時間步距內(nèi)加速度為線性加速度的一種數(shù)值方法，用內(nèi)插公式得到體系在下一刻的運動。α方法[3]是在Newmark方法的基礎(chǔ)上，通過修改結(jié)構(gòu)動力方程的時間離散形式得到的。Chung 和 Hulbert 發(fā)展了一種無條件穩(wěn)定的隱式廣義α法，它由三個參數(shù)控制數(shù)值損耗。Runge-Kutta[4， 5]方法是在一個時間步距中內(nèi)插若干計算點，利用這些計算點上函數(shù)值線性組合來代替函數(shù)的泰勒展開中的高階導(dǎo)數(shù)，從而提高精度階。不同于兩步信息預(yù)測，線性多步法[6]發(fā)展了多步信息來預(yù)測下一步，從而獲得了更高的精度。

逐步積分是最主要的時域積分方法，而它最常規(guī)的做法是差分法。時域有限差分法（Finite Difference Method， 簡稱FDM），是地震波傳播模擬最廣泛被使用的一類方法[kelly―Marfurt，1990][7-10]。有限元差分是將微分方程中的微分項用相應(yīng)的差商代替，從而將微分方程轉(zhuǎn)化為代數(shù)形式的差分方程。

由微商和差商的定義可以知道，微分的有限形式是差分，而導(dǎo)數(shù)的有限形式是差商。而微分和導(dǎo)數(shù)是以極限形式表示。數(shù)值計算方法導(dǎo)數(shù)可以用差商的自變量趨近于零來代替，換句話說，位移對時間的求導(dǎo)可以用有限差分的方式得到，位移的一階導(dǎo)數(shù)是速度，二階導(dǎo)數(shù)加速度。當(dāng)世界步距為等步長，得到中心差分。差商代替微分方程中的導(dǎo)數(shù)，就可以得到微分方程的有限差分形式。。較之傳統(tǒng)有限元法，雖然在定義幾何結(jié)構(gòu)上不夠靈活，但時域有限差分法具有顯式計算的優(yōu)勢，計算效率高，計算精度高于顯示有限元法。這些方法的缺點是精度不高，只有一階或者二階精度，難以模擬高頻問題，這類無法避免算法阻尼，從而形成較大的誤差。

為了克服低精度的問題，很多高精度的數(shù)值積分方法相繼提出。不僅僅四階，五階精度，甚至更高精度的數(shù)值積分方法處在發(fā)展之中。Golley[11]為了得到三階精度的格式，對時間域采用高斯點作為配置點。在哈密頓變分原理的基礎(chǔ)上，Riff和Barch[12] 采用3次多項式構(gòu)造插值函數(shù)，得到了有條件穩(wěn)定四階精度數(shù)值方法。Argyris[13]等在前人的基礎(chǔ)上，采用Hermitian插值，得到無條件穩(wěn)定四階精度數(shù)值方法。Kujawski和Gallagher[14]從另外一個角度，利用廣義最小二乘法，在無阻尼結(jié)構(gòu)中，構(gòu)造了一種四階精度的無條件穩(wěn)定積分格式。Tarnow和Simo[15]利用二階精度算法的結(jié)果，在此基礎(chǔ)上縮短時間間隔，構(gòu)造了一種四階精度算法。鐘萬勰[16-19]在1994年提出了精細時程積分法。在保守系統(tǒng)下，積分結(jié)果保持系統(tǒng)守恒量不變。1995年在之前工作之上，鐘萬勰提出了子域精細時程積分法，提高了計算效率，使工作量大大減小，存儲量大幅減小，為精細法應(yīng)用提供了基礎(chǔ)。2000年，顧元憲[20]提出了增維精細積分法，改進了矩陣的運算，提高了精度。但局限條件較多。2004年，汪夢甫[21]利用精細積分方法基本原理，采用高斯積分方法，建立了更新精細積分方法，這種精細方法適應(yīng)度高，為得到了廣泛的應(yīng)用提供了條件，并且提高了精度。理論上汪夢甫分析了精細方法得到任意精度的可能。2009年，富明慧[22-25]在汪夢甫研究的基礎(chǔ)上提出了高效高精度廣義精細積分法。

這類方法的困難在于不容易構(gòu)造精度較高的時間離散模式，并且空間有限元在時域上每個時間步逐步推進，因而會產(chǎn)生誤差累積[26-32]。

3.時空有限元方法

時空有限元最早由J.TOden，I.Fried和J.H.Argyris等人提出。利用哈密頓原理建立關(guān)于時間邊界的變分原理。幾十年來，在各個領(lǐng)域得到全面的發(fā)展。在波動問題，動力問題，流體問題等方向得到廣泛的應(yīng)用。

傳統(tǒng)的數(shù)值方法假定時間和空間是相互獨立的，這樣的假定廣泛應(yīng)用在實踐當(dāng)中并且在數(shù)學(xué)上很好理解。但是，使用上述技術(shù)同時也導(dǎo)致了數(shù)學(xué)上的困難。因為有用的信息可能通過速度在結(jié)構(gòu)傳播，而沒有分離的時空格式能規(guī)避這種類型的數(shù)值困難。這就要求對時間的離散和對空間離散一樣使用有限元。例如，當(dāng)結(jié)構(gòu)是三維時，這種格式需要四維的維度來表示。從而需要對時間和對空間使用相同的離散方式。空間有限元對時間和空間分別離散，空間節(jié)點形成的網(wǎng)格只能處在同一時間下，形成的是規(guī)則網(wǎng)格。時空有限元同時對時間和空間進行離散，理論上可以把網(wǎng)格劃成任意形狀，不必考慮相同的時間值，可以靈活的劃分。有限元方法區(qū)別于其它方法在于它利用能量等效原理將偏微分方程進行積分。得到方程的弱形式。恰當(dāng)?shù)淖兎中问绞怯邢拊欠衲軌蚶媚芰康刃У年P(guān)鍵。而時空有限元能否成功取決于能否找到對應(yīng)的變分原理。

R.Riff和M.Baruch等人建立了一種能同時求解所有變量的時間有限元，這種有限元格式借鑒了空間有限元，把整個時間域看作是空間的延續(xù)，從而能夠求解不同時刻的變量。馮康是提出 羅恩在馮康的基礎(chǔ)上，完善了哈密頓型變分原理，發(fā)展了非傳統(tǒng)哈密頓變分原理。羅恩對比等價的正則方程和相空間變分原理，認為即使是等價的，但由于形式的不同，產(chǎn)生的算法并不一定會有相同的效果。其結(jié)果就是相空間變分原理計算效率更高，也更接近物理問題的本質(zhì)特征。沿著這個思路，羅恩建立了非傳統(tǒng)相空間哈密頓變分原理。劉世奎對哈密頓原理進行了推廣，構(gòu)建時間邊界條件，建立了有兩個參數(shù)的廣義哈密頓變分原理，完成了從單一變量向多種變量的轉(zhuǎn)變。卓家壽等人在哈密頓體系下，推導(dǎo)了幾種變分原理的等價形式。羅恩在此基礎(chǔ)上推導(dǎo)了類變量廣義哈密頓變分原理，它包含了所有的條件。基于這種變分原理提出在時間子域上進行五次樣條插值的時間子域法。 2007年鐘萬勰[38]發(fā)現(xiàn)時空有限元可以更高效的解決邊界問題和多尺度剛度問題。2013年朱寶[39]在鐘萬勰的基礎(chǔ)上進一步研究多尺度和邊界問題，對其穩(wěn)定性進行了研究。

獲得高階精度，只需要提高多項式基函數(shù)的次數(shù)，從理論上來說，時空有限元可以獲得任意精度。空間域增加時間域之后，單元的幾何性質(zhì)簡單，避免了空間有限元的復(fù)雜邊界。讓傳統(tǒng)的有限元得到更廣的應(yīng)用。是一種有很大發(fā)展空間的數(shù)值方法。

4.微分求積方法

Bellman和Casti[40]在1971年提出微分求積法的基本原理。此后，微分求積法因為原理簡單，廣泛應(yīng)用在工程問題中，微分求積法得到快速發(fā)展。

微分求積法即DQM方法，本質(zhì)上來說，函數(shù)和它的導(dǎo)數(shù)在給點節(jié)點的值用全部節(jié)點的函數(shù)值乘以系數(shù)并求和來代替。從而讓微分方程可以轉(zhuǎn)化成關(guān)于節(jié)點函數(shù)值的一組代數(shù)方程組。由此可知，DQM是一種數(shù)值技術(shù)，它通常被用來解決初值和邊界問題。從本質(zhì)上來說，DQM是特殊的一種加權(quán)殘值法，而且是高階的有限差分法。DQM方法相對有限元方法而言，并不需要變分原理就可以求解微分方法。

從微分積分法的原理出發(fā)，可以發(fā)現(xiàn)影響數(shù)值精度主要由兩個因素構(gòu)成。一方面是權(quán)系數(shù)的值，另一個方面是選取合適的網(wǎng)格離散點。其中網(wǎng)格離散點的選取和假設(shè)試函數(shù)的模式可以確定權(quán)系數(shù)，從另一個角度來說，試函數(shù)的假設(shè)和網(wǎng)格點的選取是決定精度的關(guān)鍵因素。從而研究人員也沿著這個思路對微分求積方法進行了探索。利用多項式是有限元試函數(shù)選取的基本思路。Bert和Wang[41]為試函數(shù)來求權(quán)系數(shù)，此時構(gòu)成線性方程組的系數(shù)矩陣是勒得蒙矩陣。但由于當(dāng)離散點數(shù)目增多以后，勒得蒙矩陣會出現(xiàn)病態(tài)。所以出現(xiàn)很大的誤差。后來。Quan[42， 43]用采用了Lagrange插值，得到了微分積分法一階和二階精度的公式。Bert和Striz[41]建立了HDQ方法，采用用不同于多項式的諧函數(shù)作為試函數(shù)，開闊了研究思路。由上可知，試函數(shù)的選取并不是單一的，可以從多個角度來選取。不僅是諧函數(shù)或者多項式，甚至是指數(shù)函數(shù)，對數(shù)函數(shù)等初始函數(shù)都可以作為試函數(shù)進行研究。根據(jù)需要選擇混合初始函數(shù)來得到試函數(shù)是值得探索的方向。網(wǎng)格點的選擇方面，研究發(fā)現(xiàn)一些問題對節(jié)點選取是很敏感的，等距網(wǎng)點因為使用方便而被先采用，但是結(jié)果發(fā)現(xiàn)得到的解不夠理想。其實真實的狀況讓均勻網(wǎng)格模擬顯得不夠合理，發(fā)展非均勻網(wǎng)格更可能得到高精度方法。 Bellman[40]就用勒讓得多項式零點進行了研究，發(fā)現(xiàn)用其作為網(wǎng)格點提高了精度。在這啟發(fā)之下，Quan等研究了切比雪夫多項式零點，并且用之與其它正交多項式作了對比，發(fā)現(xiàn)切比雪夫多項式零點更有優(yōu)勢。此外，微分求積的研究的方向是更加具體的研究，Bellman[44]用微分求積法應(yīng)用到初值非線性偏微分方法得到高效精確的解法。 在多維問題上，微分積分法也得到了應(yīng)用，Civan[45]得到了多維積分微分方程。Bert[46]首先將這種方法用到結(jié)構(gòu)力學(xué)問題的求解。2001年在DQM法則的基礎(chǔ)上，F(xiàn)ung[47， 48]建立了一種不同于邊值問題的動力微分方程解法，解決了初值問題的動力微分方程。并研究時間網(wǎng)點選擇方式對數(shù)值精度和穩(wěn)定性的影響，提出了一種高精度的時間網(wǎng)格離散方法。

微分求積法雖然發(fā)展的歷史比較短，但是由于這種方法原理簡單，精度高，計算效率高，處理數(shù)據(jù)方便等優(yōu)點。在工程領(lǐng)域有廣泛的應(yīng)用。

5.結(jié)語

有限元方法一直在數(shù)值模擬中很占有重要地位，這種思想在理論研究和實際應(yīng)用中發(fā)揮著很重要的作用。以有限元原理為基礎(chǔ)，發(fā)展的新方法讓數(shù)值計算展現(xiàn)出新的活力。但是數(shù)值模擬是一門深奧的學(xué)問，在理論上和實際應(yīng)用中還有很多不完善的地方，需要克服的缺點還有很多，本文作者僅僅就自己所涉及的研究領(lǐng)域做了一些簡單的論述。

參 考 文 獻

[1]. Newmark， N.M. A method of computation for structural dynamics[C]. in Proc. ASCE. 1959.

[2]. Wilson E L. A computer program for the dynamic stress analysis of underground structures[R]. DTIC Document， 1968.

[3]. Trefethen L N. Finite difference and spectral methods for ordinary and partial differential equations[M]. Cornell University [Department of Computer Science and Center for Applied Mathematics]， 1996.

[4]. 易大義，陳道琦. 數(shù)值分析引論[M].杭州：浙江大學(xué)出版社，1998.

[5]. 胡健偉，湯懷民. 微分方程數(shù)值方法[M].北京：科學(xué)出版社，1999.

[6]. Chew W C. Waves and fields in inhomogeneous media[M]. IEEE press New York， 1995.

[7].Wilson E L， Bathe K. Numerical methods in finite element analysis[M]. Prentice-Hall， 1976.

[8]. Clough R W， Penzien J. Dynamics of structures[R]. 1975.

[9]. Wilson E L， Farhoomand I， Bathe K J. Nonlinear dynamic analysis of complex structures[J]. Earthquake Engineering & Structural Dynamics. 1972， 1（3）： 241-252.

[10]. 劉晶波，杜修力.結(jié)構(gòu)動力學(xué)[M]. 北京：機械工業(yè)出版社，2005.

[11]. Golley B W. A TIMESTEPPING PROCEDURE FOR STRUCTURAL DYNAMICS USING GAUSS POINT COLLOCATION[J]. International Journal for Numerical Methods in Engineering. 1996， 39（23）： 3985-3998.

[12]. Riff R， Baruch M. Time finite element discretization of Hamilton's law of varying action[J]. AIAA journal. 1984， 22（9）： 1310-1318.

[13]. Argyris J H， Vaz L E， Willam K J. Higher order methods for transient diffusion analysis[J]. Computer Methods in Applied Mechanics and Engineering. 1977， 12（2）： 243-278.

[14]. Kujawski J， Gallagher R H. A generalized leastsquares family of algorithms for transient dynamic analysis[J]. Earthquake engineering & structural dynamics. 1989， 18（4）： 539-550.

[15]. Tarnow N， Simo J C. How to render second order accurate time-stepping algorithms fourth order accurate while retaining the stability and conservation properties[J]. Computer Methods in Applied Mechanics and Engineering. 1994， 115（3）： 233-252.

[16]. 鐘萬勰. 結(jié)構(gòu)動力方程的精細時程積分法[J].大連理工大學(xué)學(xué)報.1994，34（2）：131-136.

[17]. 鐘萬勰. 計算結(jié)構(gòu)力學(xué)與最優(yōu)控制[M].大連：大連理工大學(xué)出版社，1993.

[18].鐘萬勰. 應(yīng)用力學(xué)的辛數(shù)學(xué)方法[M].北京：高等教育出版社，2006.

[19]. Zhong W X， Williams F W. A precise time step integration method[J]. Proceedings of the Institution of Mechanical Engineers， Part C： Journal of Mechanical Engineering Science. 1994， 208（6）： 427-430.

[20]. 顧元憲，陳飚松，張洪武.結(jié)構(gòu)動力方程的增維精細積分法[J].力學(xué)學(xué)報.2000，32（4）.

[21].汪夢甫，周錫元. 結(jié)構(gòu)動力方程的更新精細積分方法[J].力學(xué)學(xué)報.2004，36（2）：191-195.

[22]. 富明慧，劉祚秋. 固體力學(xué)中的變分差分方法[J].中山大學(xué)學(xué)報：自然科學(xué)版.2001，40（2）：13-15.

[23]. 富明慧，林敬華，劉祚秋.結(jié)構(gòu)動力分析的廣義精細積分法[Z].第九屆全國振動理論及應(yīng)用學(xué)術(shù)會議論文摘要集.2007.

[24]. 富明慧，林敬華. 精細積分法在非線性動力學(xué)問題中的應(yīng)用[J].中山大學(xué)學(xué)報：自然科學(xué)版.2008，47（3）： 1-5.

[25]. 富明慧，廖子菊，劉祚秋. 結(jié)構(gòu)動力方程的樣條精細積分法[J].計算力學(xué)學(xué)報.2009， 26（3）： 379-384.

[26]. 裘春航，呂和祥，鐘萬勰. 求解非線性動力學(xué)方程的分段直接積分法[J].力學(xué)學(xué)報.2002， 34（3）： 369-378.

[27]. 鄭兆昌，沈松，蘇志霄.非線性動力學(xué)常微分方程組高精度數(shù)值積分方法[J].力學(xué)學(xué)報.2003，35（3）： 284-295.

[28]. 胡海昌.彈性力學(xué)的變分原理及其應(yīng)用[M].北京：科學(xué)出版社，1981.

[29]. 趙秋玲.非線性動力學(xué)方程的精細積分算法[J].力學(xué)與實踐.1998，20（6）：24-26.

[30]. 王金東，高鵬. 波動方程的精細逐步積分法[J].力學(xué)季刊.2000，21（3）：316-321.

[31]. 張洪武.關(guān)于動力分析精細積分算法精度的討 ）[J]. 力學(xué)學(xué)報. 2001，33.

[32]. 梁立孚胡海昌.一般力學(xué)中三類變量的廣義變分原理[J].中國科學(xué)： A 輯. 2000，30（12）：1130-1135.

[33]. 羅恩.幾何非線性彈性動力學(xué)中廣義 Hamilton 型擬變分原理[J].中山大學(xué)學(xué)報（自然科學(xué)版）.1990， 29（2）：15-19.

[34]. 羅恩，黃偉江.相空間非傳統(tǒng) Hamilton 型變分原理與辛算法[J].中國科學(xué)：A輯.2002，32（12）： 1119-1126.

[35].羅恩，梁立孚，李緯華.分析力學(xué)的非傳統(tǒng) Hamilton 型變分原理[J].中國科學(xué)：G 輯.2007，36（6）： 633-643.

[36]. 梁立孚，羅恩，馮曉九.分析力學(xué)初值問題的一種變分原理形式[J].力學(xué)學(xué)報. 2007， 23（1）： 106-111.

[37].羅恩，朱慧堅，黃偉江.Hamilton 彈性動力學(xué)及其辛算法[J].中山大學(xué)學(xué)報（自然科學(xué)版）. 2003，42（5）： 131-132.

[38]. 鐘萬勰，高強. 時間-空間混和有限元[J].動力學(xué)與控制學(xué)報，2007，5（1）：1-7.

[39]. 朱寶，高強，鐘萬勰. 三維時間-空間混和有限元[J].計算力學(xué)學(xué)報，2013，30（3）：331-335.

[40]. Bellman R， Casti J. Differential quadrature and long-term integration[J]. Journal of Mathematical Analysis and Applications. 1971， 34（2）： 235-238.

[41]. Bert C W， Xinwei W， Striz A G. Differential quadrature for static and free vibration analyses of anisotropic plates[J]. International Journal of Solids and Structures. 1993， 30（13）： 1737-1744.

[42]. Quan J R， Chang C T. New insights in solving distributed system equations by the quadrature method―I. Analysis[J]. Computers & chemical engineering. 1989， 13（7）： 779-788.

[43]. Quan J R， Chang C. New insights in solving distributed system equations by the quadrature method―II. Numerical experiments[J]. Computers & chemical engineering. 1989， 13（9）： 1017-1024.

[44]. Bellman R， Kashef B G， Casti J. Differential quadrature： a technique for the rapid solution of nonlinear partial differential equations[J]. Journal of Computational Physics. 1972， 10（1）： 40-52.

[45]. Civan F， Sliepcevich C M. Differential quadrature for multi-dimensional problems[J]. Journal of Mathematical Analysis and Applications. 1984， 101（2）： 423-443.

[46]. Bert C， Jang S， Striz A. New methods for analyzing vibration of structural components[Z]. Structures， Structural Dynamics and Materials Conference， 28 th， Monterey， CA， Apr. 6-8， 1987 and AIAA Dynamics Specialists Conference， Monterey， CA. 1987，936-943.

[47]. Fung T C. Solving initial value problems by differential quadrature method―part 1： firstorder equations[J]. International Journal for Numerical Methods in Engineering. 2001， 50（6）： 1411-1427.

[48]. Fung T C. Solving initial value problems by differential quadrature method―part 2： secondand higherorder equations[J]. International Journal for Numerical Methods in Engineering. 2001， 50（6）： 1429-1454.

數(shù)值方法范文第2篇

Lighthill利用保角變換的方法首先提出了二維翼型的反設(shè)計方法，Hicks，Murman和Henne等人將此方法發(fā)展為可應(yīng)用機設(shè)計的工程設(shè)計方法。后Campbell等提出過一種帶約束的直接迭代的表面曲率（CDISC）方法，Yu將其與N-S解算器耦合形成了一種翼型和機翼的設(shè)計方法。波音公司則將此方法發(fā)展成工程應(yīng)用的設(shè)計方法，并廣泛地應(yīng)用于波音的B757，B777和B737NG等型號的設(shè)計過程，取得了很好的效果。例如在B777研制中由于使用了反設(shè)計方法，僅經(jīng)過三輪機翼的設(shè)計便取得了滿意的結(jié)果，使風(fēng)洞實驗的機翼模型大大少于過去B757和B767設(shè)計時的數(shù)目，充分表明了該設(shè)計工具的作用。可以說，反設(shè)計方法曾對民機設(shè)計起過革新性的推動作用；但反設(shè)計方法也有其固有的弱點（參見文獻[13]的附錄D）：首先，對于高度三維的流動要找到“好”的壓強分布很困難；其次，不能保證所得結(jié)果為最優(yōu)，即既具有高速巡航低阻的特性又在非設(shè)計條件下具有可接受的性能；最后，其他學(xué)科的約束會導(dǎo)致反復(fù)迭代。

低可信度CFD模型的數(shù)值優(yōu)化方法

隨著計算能力和數(shù)值優(yōu)化方法的快速發(fā)展，應(yīng)用基于CFD的數(shù)值優(yōu)化方法于民機設(shè)計得到了很大的發(fā)展。這一方法的應(yīng)用也從低可信度CFD模型開始，逐漸發(fā)展到采用先進的N-S方程解算器。波音公司發(fā)展了一種耦合TRANAIR[16]（一種全速勢方程的有限元方法，可參見文獻[13]附錄B）和梯度優(yōu)化方法的數(shù)值優(yōu)化氣動力設(shè)計方法，并在1992年形成了TRANAIR優(yōu)化器的雛形[17]。經(jīng)過近十年的改進，得到了一個適用于位勢流/邊界層耦合飛行條件的氣動力優(yōu)化設(shè)計工具[18-20]，具有多點優(yōu)化設(shè)計能力，可處理高達600個幾何自由度和45000個非線性不等式的約束條件（圖1表示了TRANAIR優(yōu)化過程示意圖）。作為一個例子，圖2給出了采用該軟件對機翼/發(fā)動機短艙設(shè)計計算前后壓強分布的對比，圖a和圖b分別表示了設(shè)計前后等馬赫數(shù)線的分布。可以看出圖a中掛架處出現(xiàn)激波；圖b中短艙附近的機翼表面上消除了由于短艙干擾形成的激波。算例結(jié)果表明該設(shè)計軟件可以處理很復(fù)雜的飛機/發(fā)動機綜合設(shè)計問題。

高可信度CFD模型的數(shù)值優(yōu)化方法現(xiàn)代優(yōu)化算法可以分為依賴和不依賴梯度的方法兩大類。

1.依賴梯度的優(yōu)化算法

目前可用的大多數(shù)依賴梯度的數(shù)值優(yōu)化方法都是從控制理論出發(fā)的，Jameson是此類方法的先驅(qū)者之一。盡管最初是由Pironneau提出利用控制理論進行橢圓方程系主控的外形優(yōu)化的[21-22]，但Jameson首先提出了通過控制理論自動進行外形優(yōu)化的伴隨方程方法[23]并應(yīng)用于跨聲速流動。后來，Jameson和他的合作者，還有其他研究者，大力發(fā)展此方法，從全位勢方程到Euler/N-S方程，從無粘設(shè)計到有粘設(shè)計，甚至從氣動設(shè)計到氣動/結(jié)構(gòu)的耦合設(shè)計，形成了大量文獻[24-36]。此方法不同于一般梯度優(yōu)化方法之處在于它將外形作為一個自由表面，促使流動解和最終優(yōu)化的外形同時趨于收斂，因而使優(yōu)化方法具有很高的效率（其基本思想可參見文獻[13]附錄D）。

2.不依賴梯度的優(yōu)化算法

最早無需梯度的優(yōu)化算法有Powell（共軛方向法）[37]和Nolder-Mead的單純形法[38]。最近Sturdza還應(yīng)用后者于空氣動力的設(shè)計[39]。近二十多年來人們更多地使用諸如模擬退火法[40]和遺傳算法(GeneticAlgorithm-GA)等的搜索方法，特別后者更為人們所關(guān)注。Holland利用進化理論創(chuàng)造了遺傳算法[41]（可參閱文獻[13]附錄D），即模仿生物的自然選擇進行搜索以尋求最優(yōu)解。與傳統(tǒng)的搜索和優(yōu)化方法相比，遺傳算法具有下述4個特點[42-45]：1）不是直接作用于參變量集本身，而是對參變量集的某種編碼運算。2）不是對單個點而是對多個點構(gòu)成的群體進行搜索。3）直接計算適應(yīng)值（函數(shù)），無需導(dǎo)數(shù)和其他輔助信息。4）利用概率轉(zhuǎn)移原則，而非傳統(tǒng)優(yōu)化方法中的確定性原則。已有愈來愈多的研究和民機研制機構(gòu)表現(xiàn)出了對這種隨機尋優(yōu)方法的濃厚興趣，也已出現(xiàn)了不少利用遺傳算法進行翼型或機翼優(yōu)化計算的文獻[46-56]。

3.對高可信度CFD模型數(shù)值優(yōu)化方法的要求

分析最近十余年中出現(xiàn)的大量基于Euler/N-S方程的數(shù)值優(yōu)化方法和文獻，可以看出多數(shù)仍表現(xiàn)為學(xué)院式的探討，提供可直接用于工程設(shè)計的方法和工具顯得尚很有限，盡管已開始向這方面努力。這可能是因為：1）只是近幾年來隨DPW研討會等的進行，數(shù)值模擬才可以比過去更正確地估算阻力值。2）工程界的空氣動力外形優(yōu)化需要在高維搜索空間中進行并存在大量的非線性約束，使優(yōu)化問題十分復(fù)雜且計算開銷巨大；3）巨大的計算量要求很豐富的計算資源和很長的計算時間，這與工程問題要求的迅速反饋相悖。

因此要使基于CFD的空氣動力優(yōu)化方法和軟件成為日常的工程設(shè)計手段和工具需解決如下技術(shù)關(guān)鍵：1）具有建立準(zhǔn)確計算諸如升力、阻力、力矩等敏感氣動特性的正確流動模型的能力。比較現(xiàn)有的氣動力優(yōu)化方法可知，大多數(shù)方法還在使用不完善的流動模型，如基于Euler方程，甚至全位勢方程等。雖然它們在一定條件下，如巡航小迎角飛行狀態(tài)，可以提供合理的結(jié)果，但工程應(yīng)用常要求準(zhǔn)確地估算出阻力、俯仰力矩等敏感的氣動特性，要求可計算整個飛行包線的飛行狀態(tài)以及不同的復(fù)雜的幾何外形等，這只能通過求解N-S方程來實現(xiàn)。順便指出，有些文獻（如文獻[28]）雖以N-S方程為主控方程，但優(yōu)化時的伴隨運算子卻是在沒有考慮粘性流動的假設(shè)下得出的（參見文獻[28]第6節(jié)）。為了提高計算準(zhǔn)確度,最好在離散N-S方程時使用高階的差分算子[53-54]。2）具有尋求全局最優(yōu)的能力。通常基于梯度的算法容易陷入局部最優(yōu)，而遺傳算法等隨機搜索的方法則具有取得總體最優(yōu)的優(yōu)點。3）能有效地處理大量幾何和氣動力的非線性約束。優(yōu)化問題的最優(yōu)解常常是位于不同維超曲面（hyper-surface）的交匯處，遺傳算法不同于基于梯度的方法，不限于目標(biāo)函數(shù)的光滑擴展，可應(yīng)用于多重約束的情況[53-54]。4）可應(yīng)用于不同的幾何外形和設(shè)計條件。5）掃描高維搜索空間的計算有效性高，以滿足設(shè)計周期和研制成本的要求。遺憾的是這正是遺傳算法的主要缺點，即估算適應(yīng)函數(shù)的高代價。可以采用多處理器上的有效并行計算來大大減少計算時間[57]，或在估算適應(yīng)函數(shù)值時采用近似模型，如降階模型[54,58]或響應(yīng)面模型[50]等。

數(shù)值優(yōu)化方法的發(fā)展現(xiàn)狀和驗證研究#p#分頁標(biāo)題#e#

1.空氣動力優(yōu)化設(shè)計計算的系列研討會

近年來CFD學(xué)術(shù)界和航空業(yè)界都十分關(guān)注計算阻力的精度問題，這也是CFD應(yīng)用于工程設(shè)計時所面臨的第一個具有挑戰(zhàn)性的計算。AIAA的應(yīng)用空氣動力學(xué)專業(yè)委員會在各方支持下，自2001年開始舉行了DPW（DragPredictionWorkshop）系列會議[59]，參與者都用N-S方程求解相同的幾何外形（翼/身組合體，翼/身/短艙/掛架的復(fù)雜組合體等），得到了一個巨大的計算結(jié)果數(shù)據(jù)集，可與已有的已經(jīng)過修正的風(fēng)洞試驗值比較。由于參與的計算者所采用的數(shù)值方法、湍流模型、計算網(wǎng)格形式及數(shù)目等各不相同，此數(shù)據(jù)集可用作分析和討論各種因素對CFD計算結(jié)果的影響。該系列會議至今已舉行了5屆，對推動和提高CFD計算阻力的精度很有意義。文獻[13]的附錄C中給出了前3屆結(jié)果的分析和討論。鑒于DPW系列會議的成功，AIAA應(yīng)用空氣動力學(xué)專業(yè)委員會針對CFD面臨的第二個挑戰(zhàn)---計算三維高升力外形的最大升力CLmax，于2009年發(fā)起并組織了類似的高升力計算研討會，其第一次會議(HiLiftPW-I)已于2010年6月在美國舉行,文獻[60]是該次會議的總結(jié)。在上述工作的基礎(chǔ)上，2013年1月AIAA又進一步在其ASM會議過程中形成了以加拿大McHill大學(xué)Nadarajah教授為首的空氣動力優(yōu)化設(shè)計討論組，作為空氣動力優(yōu)化設(shè)計計算系列研討會實際的組委會。討論組討論了：1）建立可供在一個有約束的設(shè)計空間中測試氣動優(yōu)化方法的一組標(biāo)準(zhǔn)算例。2）舉行研討會的時間。與會者一致認為，由于工業(yè)界對基于CFD的氣動外形數(shù)值優(yōu)化方法有強烈的需求，優(yōu)化方法和工具的研制也已有了相當(dāng)?shù)陌l(fā)展，可以以類似于DPW的研討會形式，通過對一系列復(fù)雜氣動外形的優(yōu)化，來評估現(xiàn)有的尋求最小阻力外形的各種優(yōu)化方法的能力，并將結(jié)果向工業(yè)界/研究機構(gòu)公布。與會者還認為第一次會議從二維和三維機翼外形開始是合適的，并請加拿大的與會者準(zhǔn)備標(biāo)準(zhǔn)算例。第一次會議擬于2013或2014年的AIAA應(yīng)用空氣動力會議期間舉行。

2.先導(dǎo)性的研究

事實上為準(zhǔn)備此研討會，波音的Vassberg，斯坦福的Jameson，以色列的Epstein及Peigin等三方從2007年起即開始了先導(dǎo)性的研究（pilotproject）,以積累經(jīng)驗和發(fā)現(xiàn)問題。三方用各自己開發(fā)的優(yōu)化軟件（MDOPT，SYN107，OPTIMAS）對第三屆DPW會議的測試機翼DPW-W1獨立地作優(yōu)化計算[61,62]。波音研制的MDOPT[63]（也可參見文獻[13]的1.7.3節(jié)）可使用響應(yīng)面模型（InterpolatedRe-sponseSurface—IRS）的數(shù)值優(yōu)化格式[64]，也可直接從流場解計算設(shè)計變量的靈敏度代替IRS模型完成優(yōu)化。其流場解軟件為TLNS3D[65]，計算網(wǎng)格點為3582225。Jameson開發(fā)的SYN107采用基于梯度的“連續(xù)”伴隨方程方法[23,31]，其流場解軟件為FLO107，計算網(wǎng)格點為818,547。

以色列航空公司開發(fā)的OPTIMAS采用降階模型的GA算法，流場解軟件為NES[66-68]，計算網(wǎng)格點為250,000。對三方獨立優(yōu)化后所得的外形再用不參與優(yōu)化的流場解軟件OVERFLOW[69]作評估計算，計算網(wǎng)格點數(shù)為4,000,000，以便能準(zhǔn)確地計算阻力。結(jié)果表明，4個分析軟件計算得到的阻力增量值的分散度在Ma=0.76時為5counts(1count=0.0001)，Ma=0.78時為10counts，因此很難確定哪個優(yōu)化后外形最優(yōu)。但從Ma=0.76，C=0.5單設(shè)計點的阻力改進結(jié)果（表1）[61]看，OPTIMAS優(yōu)化后的04外形明顯優(yōu)于MDOPT優(yōu)化后的M5和SYN107優(yōu)化后的S4。文獻[61]還討論了從比較中可吸取的經(jīng)驗和教訓(xùn)。

一種基于高可信度CFD模型的數(shù)值

優(yōu)化方法的構(gòu)造本節(jié)將以O(shè)PTIMAS為例對如何滿足可應(yīng)用于工程實踐的高可信度CFD模型數(shù)值優(yōu)化方法的要求做一說明。

1.優(yōu)化方法的構(gòu)造及其特點

OPTIMAS是將遺傳優(yōu)化算法和求解全N-S方程的分析算法相結(jié)合的一種有效并魯棒的三維機翼優(yōu)化方法。1)其全N-S方程的流場并行解算器NES[66-67]基于高階低耗散的ENO概念（適用于在多區(qū)點對點對接網(wǎng)格中的多重網(wǎng)格計算）[66,71]和通量插值技術(shù)相結(jié)合的數(shù)值格式，采用SA湍流模型，可快速準(zhǔn)確地完成氣動力計算，因此具有計算大量不同流動和幾何條件的魯棒性。作為例子圖3和4給出了ARA翼身組合體Ma=0.80，Re=13.110時的升阻極線和CL=0.40時的阻力發(fā)散曲線[68]，使用的網(wǎng)格點數(shù)分別為，細網(wǎng)格（3lev）：900,000，中等網(wǎng)格（2lev）：115,000。可見升阻極線直到大升力狀態(tài)的計算與實驗都很一致。對比圖中還給出的TLNS3D在細網(wǎng)格（2,000,000）中的計算值可見，無論升阻極線或阻力發(fā)散曲線NES的都更優(yōu)。作為數(shù)值優(yōu)化軟件的特點之一是其在流場解算器中首次使用了高精度格式。2)優(yōu)化計算的遺傳算法中采用了十進制編碼、聯(lián)賽選擇算子[42]、算術(shù)交叉算子、非均勻?qū)崝?shù)編碼變異算子[72]和最佳保留機制。為解決搜索時總體尋優(yōu)耗時大和求解N-S方程估算適應(yīng)函數(shù)代價高的問題，在尋優(yōu)過程中估算適應(yīng)函數(shù)時采用當(dāng)?shù)財?shù)據(jù)庫中的降階模型[54,58]獲取流場解（當(dāng)?shù)財?shù)據(jù)庫是在搜索空間中離散的基本點處求解全N-S方程建立的），并以多區(qū)預(yù)測-修正方法來彌補這種近似帶來的誤差。多區(qū)預(yù)測-修正方法即在搜索空間的多個區(qū)域并行搜索得到各區(qū)的優(yōu)化點，再通過求解全N-S方程的驗證取得最優(yōu)點。為保證優(yōu)化的收斂，尋優(yōu)過程采用了迭代方法。3)在整個空間構(gòu)筑尋優(yōu)路徑（圖5），擴大了搜索空間和估算適應(yīng)函數(shù)的區(qū)間[54]。4）為提高計算效率，OPTIMAS包含了五重并行計算：Level1并行地求解N-S方程Level2并行地掃描多個幾何區(qū)域，提供多個外形的適應(yīng)函數(shù)的計算（level1隱于level2中）。Level3并行的GA優(yōu)化過程（level3隱于level4中）。Level4并行地GA搜索多個空間。Level5并行地生成網(wǎng)格。5）采用單參數(shù)或雙參數(shù)的BezierSpline函數(shù)對搜索空間參數(shù)化;并基于優(yōu)化外形與原始外形的拓撲相似自動地實現(xiàn)空間網(wǎng)格的快速變換。

2.優(yōu)化設(shè)計的典型結(jié)果

文獻[53]~文獻[58]給出的大量算例充分表明了OPTIMAS優(yōu)化軟件的優(yōu)異性能。本文5.2中給出了其優(yōu)化三維機翼的性能，這里再補充兩例。1）翼身組合體整流（fairing）外形的優(yōu)化文獻[73]討論了某公務(wù)機翼身組合體機翼外形優(yōu)化的單點和多點設(shè)計兩者性能的比較。結(jié)果表明，多點優(yōu)化設(shè)計能同時保證設(shè)計的巡航狀態(tài)時，和高Ma數(shù)飛行，起飛等非設(shè)計狀態(tài)時的良好性能。文獻[74]進一步討論了翼身組合體整流外形的優(yōu)化設(shè)計。流動的復(fù)雜性（三維粘流/無粘流強相互作用的流動區(qū)域）和幾何的復(fù)雜性（三維非線性表面）使整流外形的設(shè)計經(jīng)歷了傳統(tǒng)的試湊法，基于Euler解的試湊法等，最后才發(fā)展為現(xiàn)代完全N-S解的數(shù)值優(yōu)化方法。文獻[74]采用了這種方法，先作機翼外形優(yōu)化，再作整流外形優(yōu)化，然后再作機翼優(yōu)化，整流外形優(yōu)化，……依次迭代，直至收斂。優(yōu)化中用雙參數(shù)的BezierSpline函數(shù)將整流外形參數(shù)化，所得搜索空間的維數(shù)ND=（2N-2）*（M-1）決定的參數(shù)化整流外形與實際外形的差別在M=10，N=4，ND=54時可準(zhǔn)確到0.3mm（滿足工程需求）。計算網(wǎng)格數(shù)為90萬。表3給出了設(shè)計條件和約束，表4給出了設(shè)計點的阻力值比較。由表4可知，GBJ2的減阻為16.7，50%DC，GBJFR1的減阻為10.7，32.1%DC，GBJFR2的減阻為5.9，兩次優(yōu)化機翼的減阻總計為22.6，67.9%DC，優(yōu)化機翼和優(yōu)化整流外形減阻作用分別約占2/3和1/3，可見整流外形的優(yōu)化也是十分重要的。約束則使減阻損失4.6（如GBJFR3-GBJFR1）。圖6至圖9分別為原始外形，GBJ2，GBJFR2和GBJFR4的整流處等壓線分布，可見整流外形的優(yōu)化消除了原始外形和GBJ2中存在的激波。圖10和圖11分別給出了Ma=0.8時升阻極線和CL=0.4時阻力發(fā)散曲線的比較，可見優(yōu)化設(shè)計不僅對設(shè)計點，對非設(shè)計狀態(tài)也都有好處。2）翼身融合體飛機氣動外形的優(yōu)化[75]優(yōu)化設(shè)計以英國克朗菲爾德大學(xué)設(shè)計的BWB外形[76]為出發(fā)外形，該外形的主要設(shè)計點為，。在數(shù)值優(yōu)化計算中還考慮了，的第二個設(shè)計點和，（起飛狀態(tài)）的第三個設(shè)計點。幾何約束有剖面相對厚度，前緣半徑，后緣角，每個剖面的樑處還附加兩個厚度約束，其中上標(biāo)b表示出發(fā)外形，*表示優(yōu)化外形，下標(biāo)i表示第i個剖面。附加的空氣動力約束為對俯仰力矩的規(guī)定。采用Bezier樣條描述幾何外形，總設(shè)計變量為93個。表5給出了設(shè)計計算各狀態(tài)的條件和約束，其中是權(quán)系數(shù)。表6給出了優(yōu)化計算結(jié)果。#p#分頁標(biāo)題#e#

單點優(yōu)化的BWB-1結(jié)果與文獻[77]的結(jié)果相比較可見，文獻[77]采用Euler方程的無黏優(yōu)化使阻力降低了26counts；而這里的BWB-1全N-S方程優(yōu)化使阻力降低了52counts，顯示了此黏性優(yōu)化方法的優(yōu)點。比較有、無俯仰力矩約束時優(yōu)化得到的BWB-2和BWB-1表明，盡管BWB-1阻力降低的效果突出，但其值過大，出于穩(wěn)定性考慮而不能接受；BWB-2的阻力雖比BWB-1大了1.9counts，卻滿足了力矩的要求。表6中的雙點優(yōu)化設(shè)計（BWB-4），使第三設(shè)計點（低速狀態(tài)）的達到1.671（消除了BWB-2達不到設(shè)計要求1.63的缺點），且基本保持了主設(shè)計點的阻力收益，為196.6。然而BWB-4在時的阻力達216.6，高于BWB-2的213.4，表明需要三個設(shè)計點的優(yōu)化設(shè)計（BWB-3）。BWB-3在時，為202.5（比兩點設(shè)計值減小了14.1），同時滿足了其它兩個設(shè)計點的性能要求。圖12至圖15給出了所有設(shè)計狀態(tài)和時的極曲線，時的阻力發(fā)散曲線和時的隨迎角α變化的曲線。由圖可見，時所有優(yōu)化設(shè)計的極曲線都非常接近，相比于原始外形的極曲線，性能有了很大改進；時也一樣，特別是三點優(yōu)化設(shè)計的BWB-3，優(yōu)點更明顯。阻力發(fā)散曲線也都有了很大改進，在前所有的總阻力基本保持常值，單點與兩點優(yōu)化的阻力發(fā)散點接近，而三點優(yōu)化的可達附近。由圖15可知，沒有考慮低速目標(biāo)的BWB-1和BWB-2具有較低的，將低速目標(biāo)計入設(shè)計狀態(tài)的BWB-3和BWB-4所得的皆優(yōu)于原始外形的。上述結(jié)果表明三點優(yōu)化設(shè)計具有最佳的優(yōu)化效果和總體最好的氣動性能。Fig.15LiftcoefficientCLvsangleofattackatMa=0.2最后，上述各優(yōu)化結(jié)果在（主設(shè)計點）時的阻力值基本相同，但幾何外形卻差別不小，由此可見，外形阻力優(yōu)化問題沒有唯一解[75]。上述計算是在具有456GBRAM，114MB二級高速緩存的機群環(huán)境下通過“過夜”方式完成單點優(yōu)化設(shè)計，在1.5-2天的計算時間內(nèi)完成三點優(yōu)化設(shè)計的，計算時間可滿足應(yīng)用于工程設(shè)計的要求[75]。

結(jié)束語

數(shù)值方法范文第3篇

1概述

隨著科學(xué)技術(shù)水平的不斷提高和工程建設(shè)規(guī)模的不斷擴大，在土木建筑、水利工程和路橋工程中樁基承載力、沉降量大小和堤壩穩(wěn)定性等力學(xué)問題變得十分復(fù)雜。而事實上，軟土路基呈現(xiàn)出空間非線性沉降變形規(guī)律，土體的變形協(xié)調(diào)條件和應(yīng)力平衡條件也十分復(fù)雜。這些問題已經(jīng)很少能用數(shù)學(xué)方法求得精確解或通過模擬試驗得到定量解，大多數(shù)課題必須借助于計算機和計算數(shù)學(xué)用數(shù)值分析的方法求出近似解。目前用于地基沉降量分析的數(shù)值方法主要是差分法、有限元法和邊界元法，其發(fā)展趨勢是有限元法與差分法或與邊界元法相結(jié)合解決課題，以發(fā)揮各種方法的優(yōu)越性。

2分析方法解析

2.1差分法

差分法的基本精神就是將研究區(qū)域用差分網(wǎng)格離散，對每一個節(jié)點通過差商代替導(dǎo)數(shù)把固結(jié)微分方程化成差分方程，然后結(jié)合初始條件和邊界條件，求解線性方程組得到數(shù)值解。

以平面問題為例，差分法可得到所研究平面內(nèi)在各個時間的孔隙壓力的分布，因此可以導(dǎo)出初始沉降Si與任何時間t的總沉降S或固結(jié)沉降量Sc。由于土的豎向應(yīng)變?yōu)?/p>

（2.1）

地基中沿著某一鉛垂線的沉降為：

（2.2）

H為有效壓縮層厚度。由于上式中采用不排水指標(biāo)，，，所以孔隙壓力，應(yīng)用總應(yīng)力來計算。總沉降為：

（2.3）

式中，E和隨有效應(yīng)力而變化，因此上式可算得任意時間的總沉降。任何時間的固結(jié)沉降為。當(dāng)計算最終總沉降時，式中的孔隙壓力。

2.1有限元法

有限元是地基和結(jié)構(gòu)作為一個整體來分析，將其劃分網(wǎng)格，形成離散體結(jié)構(gòu)，形成有限數(shù)目的區(qū)域單元，這些單元體只在結(jié)點處有力的聯(lián)系。材料的應(yīng)力－應(yīng)變關(guān)系可表示為

（2.4）

由虛位移原理可建立單元體的結(jié)點力與結(jié)點位移之間的關(guān)系，進而寫出總體平衡方程：

（2.5）

式中，分別為勁度矩陣，結(jié)點位移列陣和結(jié)點荷載列陣。然后結(jié)合初始和邊界條件求解線性方程組，在荷載作用下算得任一時刻地基和結(jié)構(gòu)各點的位移和應(yīng)力，得到問題的數(shù)值解。有限元法可以將地基作為二維甚至三維問題來考慮，反映了側(cè)向變形的影響[1]。

（1）彈性有限元法

土的彈性應(yīng)力——應(yīng)變數(shù)學(xué)模型包括線性和非線性彈性模型兩種。用線性彈性模型計算地基的位移和沉降，只適用于不排水加荷情況，并且對破壞要有較大的安全系數(shù)，一般不發(fā)生屈服的情況。實際上土體中的應(yīng)力狀態(tài)都可能發(fā)生屈服，其應(yīng)力——應(yīng)變的關(guān)系是非線性的。此外，除了理論建模還可通過試驗擬合的途徑，即根據(jù)土體的應(yīng)力應(yīng)變試驗曲線，用彈性系數(shù)的連續(xù)變化來逼近真實的試驗曲線建立起各種不同形式的非線性彈性模型[49]，第一種是以E（彈性模量）和μ（泊松比）兩個彈性常數(shù)表達的，稱E-μ彈性模型，第二種是以K（彈性體積模量）和G（剪切模量）兩個彈性常數(shù)表達的，稱K－G模型。另外一種是南京水利科學(xué)研究所采用的非線性的變彈性體模型，它的特點是不用常規(guī)的彈性常數(shù)，改用兩個非線性函數(shù)，來表達應(yīng)力與應(yīng)變之間的關(guān)系。

典型的E-μ彈性模型是Duncan-Chang的雙曲線模型。Duncan-Chang的雙曲線模型可以考慮應(yīng)力歷史對變形的影響，若應(yīng)力低于或高于前期固結(jié)應(yīng)力，則采用不同的彈性模量計算公式。它還可考慮土與結(jié)構(gòu)共同工作，考慮復(fù)雜的邊界條件，考慮施工逐級加荷，考慮土層的各向異性等。

土體在路堤荷載作用下的變形過程，伴隨著主應(yīng)力大小的不斷變化及主應(yīng)力方向的不斷偏轉(zhuǎn)，土的彈性模量及泊松比也隨之改變。路基真實應(yīng)力場為初始應(yīng)力場疊加上由路堤填筑荷載引起的附加應(yīng)力場。假定路堤填筑荷載引起的附加應(yīng)力場可近似采用同樣荷載作用下在線彈性半空間無限體所產(chǎn)生的附加應(yīng)力場。采用鄧肯－張非線性模型描述土體本構(gòu)關(guān)系，由于路基縱向應(yīng)變?yōu)榱悖瑒t其增量形式的本構(gòu)關(guān)系可表示為[2]：

（2.6）

式中：，分別為切線彈性模量和切線泊松比，可表示為：

鄧肯－張模型含有8個參數(shù)須三軸排水試驗確定。

K-G彈性模型用三軸等向壓縮試驗測p和體積應(yīng)變，由此建立K的公式。Naylor取切線體積模量為p的線性函數(shù)，即：

（3.34）

由三軸剪切試驗可建立切線模量的公式：

（3.35）

參數(shù)均由三軸試驗確定。

（2）彈塑性有限元法

土的彈塑性模型將土的應(yīng)變分為可以恢復(fù)的彈性應(yīng)變和不可恢復(fù)的塑性應(yīng)變兩部分。彈性應(yīng)變增量可以用彈性理論計算，塑性應(yīng)變增量可以用增量塑性理論求解。土的彈塑性計算模型一般分為理想塑性和硬（軟）化塑性模型2種。

彈性非線性模型是假定全部變形都是彈性的，用改變彈性常數(shù)的方法來反映非線性；彈塑性模型則把總的變形分成彈性變形和塑性變形兩部分，用虎克定律計算彈性變形部分，用塑性理論來求解塑性變形部分。典型的彈塑性應(yīng)力－應(yīng)變模型有：劍橋模型、修正劍橋模型、拉德模型、橢圓－拋物雙屈服面模型、沈珠江雙屈服面彈塑性模型、“空間準(zhǔn)滑面”模型（SMP模型）等[3]。

3結(jié)語

土的應(yīng)力－應(yīng)變關(guān)系非常復(fù)雜，任何模型都有它的局限性。過去利用E-μ彈性模型和K-G彈性模型，雖然對壩和地基的應(yīng)力應(yīng)變分析都曾作出有意義的貢獻，但它們都忽視了土的剪脹性和應(yīng)力路徑的影響。而劍橋模型只適用于只有剪縮而沒有剪脹的正常固結(jié)粘土與松砂。計算時要找出一個數(shù)學(xué)模型來全面正確地表達土的這種特性是難以想象的。

參 考 文 獻

1 錢家歡，殷宗澤. 土工原理與計算（第二版）. 北京：中國水利水電出版社，1996

2 費正華，鄧水明. 應(yīng)用鄧肯－張非線性模型近似計算路基沉降. 中南公路工程，2001

3 河海大學(xué)，江蘇寧滬高速公路股份有限公司. 交通土建軟土地基工程手冊. 北京：人民交通出版社，2001

4 Casagrande,A. Classification and Identification of soils,Trans.ASCE,113,1948:901-991

數(shù)值方法范文第4篇

關(guān)鍵詞： 特殊角三角函數(shù)值 數(shù)形結(jié)合 函數(shù)圖像 函數(shù)單調(diào)性

高一下學(xué)期一開始，教學(xué)內(nèi)容就進入了三角函數(shù)。這一節(jié)公式很多，需要記憶的東西很多，但是只要學(xué)生能夠每天定時定量地練習(xí)題目，公式自然能夠熟練應(yīng)用，而且爛熟于心。而且學(xué)生本身對公式也比較重視，因為公式的各種靈活運用，能夠激發(fā)學(xué)生的興趣。他們做完一道題目之后，會互相討論，看還有沒有其他方法。這源于筆者平時在教學(xué)過程中不斷地鼓勵學(xué)生去思考、去總結(jié)，不但要學(xué)會，而且要會學(xué)；把新課標(biāo)強調(diào)的“提高學(xué)生自主學(xué)習(xí)能力和探究學(xué)習(xí)能力”這一思想。盡管公式學(xué)生已經(jīng)很熟悉了，但是仍有學(xué)生會在三角函數(shù)的題目上卡住。為什么呢？因為這一節(jié)還出現(xiàn)了大量的特殊角，如30°，45°，120°，甚至還有75°。學(xué)生覺得特殊角不如公式靈活，只能去死記硬背。因為對特殊角不熟悉，導(dǎo)致他們看到，卻不知道這就是tan60°；看到cos120°，還要苦想該用哪個誘導(dǎo)公式來誘導(dǎo)。雖然他們不止一次地體會到特殊角的重要性，但是他們?nèi)圆荒芙邮苡脖尺@樣傳統(tǒng)的學(xué)習(xí)方法。隨著高一課程的結(jié)束，高二的解析幾何、立體幾何中仍舊會出現(xiàn)這些特殊角。現(xiàn)在學(xué)生若是沒有記住，到了高二的時候怎么辦？

針對這個問題，筆者查閱了很多資料，大概是這個問題基本都是靠硬背來解決，因此所能找到的資料甚少。一個偶然的機會，筆者看到學(xué)生在算sin30°的時候，畫了一個30°的直角三角形，很顯然這個方法不能解決sin210°，但是筆者還是表揚了這個學(xué)生，因為他在想辦法解決問題。這個發(fā)現(xiàn)使筆者體會到，通過高一上函數(shù)部分的強化，學(xué)生現(xiàn)在已經(jīng)有了畫圖解決問題的思想，能不能用數(shù)形結(jié)合的辦法來解決這個一直讓學(xué)生比較頭痛的問題呢？其實學(xué)生在特殊角這部分暴露的問題很明顯，對［0，90°］范圍內(nèi)的角度接觸時間較久，比較熟悉，只是對高中階段才推廣的“大”角比較陌生。通過跟學(xué)生共同探討，筆者發(fā)現(xiàn)以下幾個方法比較適用。

一、利用三角函數(shù)圖像

y=sinx， y=cosx， y=tanx的圖像，在教材里面有三節(jié)內(nèi)容，對它們的圖像和性質(zhì)研究是三角函數(shù)部分的重點內(nèi)容。因此，若學(xué)生產(chǎn)能夠畫出它們的圖像，不要說cos150°，哪怕是sin225°，或者是更大的角，也能夠一眼看出。但這種方法的前提條件是，學(xué)生必須得記住這三個三角函數(shù)圖像。

二、 利用直角坐標(biāo)系

以sin225°為例，在平面直角坐標(biāo)系中，畫出225°所在的終邊，再做出它的延長線，這樣在第一象限內(nèi)就出現(xiàn)了一個以它的延長線為終邊的角，而此時學(xué)生就可選取非常熟悉的45°為此延長線的代表角。接下來做的事情，只需在延長線上取點P（x，y）。由圖像可知，兩線關(guān)于原點對稱，故此，兩線上的點的縱橫坐標(biāo)均互為相反數(shù)，則在原線上可以作出P點關(guān)于原點的對稱點P′（x，y），由任意角三角函數(shù)的定義即可得出sin225°==- =- sin45°=- 。通過剛才的推導(dǎo)過程，也可以得出這樣的結(jié)論：若兩角的終邊關(guān)于原點對稱，則它們的正弦值（或余弦值）互為相反數(shù)。這個結(jié)論的得出，讓學(xué)生知道了第三象限的角與他們熟悉的［0，90°］的角的關(guān)系，自然他們想到了第二象限和第四象限。

通過剛才的推導(dǎo)可知：若兩角終邊關(guān)于y軸對稱，則它們的余弦值互為相反數(shù)，正弦值相等。從另一個方面來看，這兩個角為互補的關(guān)系，所以剛才得出的結(jié)論也可敘述為：互補的兩角正弦值相等，余弦值相反。第四象限的角推導(dǎo)過程與上述過程類似。通過作出其關(guān)于x軸的對稱軸可知：若兩角終邊關(guān)于x軸對稱，則它們的余弦值相等，正弦值互為相反數(shù)。

綜上可知，若要解決特殊角的三角函數(shù)值，只需要在坐標(biāo)系中，畫出它的終邊所在的位置，通過做它關(guān)于原點（或x軸、或y軸）對稱線，找出第一象限我們非常熟悉的角，判斷出符號即可。

三、利用函數(shù)單調(diào)性

對于某些連［0，90°］都記不住的學(xué)生，除了用本文一開始提出的畫特殊三角形以外，還可以利用三角函數(shù)本身的單調(diào)性。由于特殊角的三角函數(shù)值只有幾個數(shù)值：0，，，，1，聯(lián)系y=sinx在［0，90°］內(nèi)單調(diào)遞增，故對號入座，sin0°=0，sin30°=，sin45°=，sin60°=，sin90°=1，相應(yīng)余弦值則可通過直角三角形里得出的結(jié)論，互余的兩角正余弦值互換得到。對于數(shù)值比較麻煩的15°和75°，我們可以通過構(gòu)造成兩角和或兩角差的方法，快速算出它們對應(yīng)的三角函數(shù)值。

筆者提出了這幾個方法后，很多學(xué)生都不再覺得特殊角三角函數(shù)值很難背了。究其原因，是在筆者提出的方法的基礎(chǔ)上，他們非常熟練地運用函數(shù)圖像、函數(shù)單調(diào)性等函數(shù)知識。這些方法中所涉及的數(shù)形結(jié)合思想，鍛煉了他們的數(shù)學(xué)思維能力，記憶的過程也成為了他們思考問題的過程。學(xué)生覺得學(xué)有所得，學(xué)有所用，這些特殊角三角函數(shù)值的記憶過程不再是枯燥無趣的幾個數(shù)字，而是生動形象的圖像、函數(shù)知識。而且，在這一過程中所涉及的數(shù)形結(jié)合方法，其本身就是高中數(shù)學(xué)階段重要方法之一。

參考文獻：

數(shù)值方法范文第5篇

【關(guān)鍵詞】地鐵；火災(zāi)；數(shù)值模擬

隨著城市的發(fā)展，地鐵已經(jīng)成為城市交通的命脈。人們大多喜歡乘坐地鐵出行，高峰時期地鐵中人員密度非常大。由于地鐵火災(zāi)的特征不同于地面建筑，當(dāng)?shù)罔F發(fā)生火災(zāi)時及易造成嚴重的財產(chǎn)損失和人員傷亡。地鐵火災(zāi)實驗很難進行，為了便于問題的探討，常選擇數(shù)值模擬的方法來對其進行研究。[1]

1.地鐵車站模型選擇

為了便于問題的研究，地鐵站火災(zāi)模擬中選取沈陽地區(qū)典型淺埋島式地鐵車站為研究對象。車站公共區(qū)域通風(fēng)空調(diào)系統(tǒng)按站廳、站臺均勻送-回（排）風(fēng)設(shè)計（閉式系統(tǒng)）。站臺層設(shè)置軌底排風(fēng)風(fēng)道，排風(fēng)口與剎車電阻箱對齊；軌頂排風(fēng)風(fēng)道，排風(fēng)口與列車空調(diào)冷凝器對齊；站臺排風(fēng)風(fēng)道。車站空調(diào)回（排）風(fēng)機兼作車站公共區(qū)消防排煙風(fēng)機，回（排）風(fēng)風(fēng)道兼作排煙風(fēng)道。站廳層送風(fēng)風(fēng)量為66000 m3/h，站臺層排煙量為132000m3/h。地鐵車站公共區(qū)域剖面如圖1所示，車站主體兩層結(jié)構(gòu)均在地面以下。

2.地鐵車站計算區(qū)域

合理簡化計算區(qū)域，在不影響整體研究效果的情況下，適當(dāng)?shù)慕档陀嬎愠杀尽Ｕ緩d與站臺層之間的結(jié)構(gòu)不作為計算區(qū)域（即忽略兩層的導(dǎo)熱作用）。站臺層長120m，候車區(qū)計算層高3.5m，寬11.8m；隧道出入口4個，尺寸為3.75m×5.26m；軌頂排煙口24個，尺寸為1m×0.6m；軌底排煙口48個，尺寸為0.8m×0.3m；送回風(fēng)口72個，尺寸為0.5m×0.5m。站廳層長88m，寬19.3m，計算層高3.5m；送風(fēng)口30個，尺寸為0.5m×0.5m；車站出入口4個，尺寸為4.5m×2.5m。

3.火源邊界條件設(shè)置

世界各國對于各種可燃物的釋熱速率并沒有明確的界定，根據(jù)香港地鐵工程技術(shù)人員的保守估計，地鐵車站公共區(qū)火災(zāi)釋熱速率為2MW[2]；而英國Building Research Establish出版的報告中顯示，在人員聚集的公共場所可能的火災(zāi)釋熱速率為2～2.5MW。根據(jù)文獻[3]，本次研究地鐵車站公共區(qū)域火災(zāi)釋熱速率5分鐘后穩(wěn)定在5MW的火災(zāi)情景。

4.數(shù)值實驗方法驗證

數(shù)值模擬方法是否正確常常需要通過相應(yīng)的實驗來進行驗證，但是地鐵站火災(zāi)實驗很難進行。下面驗證的過程是用地鐵站火災(zāi)數(shù)值模擬實驗的方法，針對已有隧道火災(zāi)實驗數(shù)據(jù)來進行的對比，以表明火災(zāi)數(shù)值模擬實驗方法的正確。

4.1模型建立

Kumar曾對隧道火災(zāi)進行了實驗測試，隧道尺寸長390m、寬5m、高4m[4，5]。火源釋熱速率在10s內(nèi)，從0線性增長到7.225MW。根據(jù)上述條件建立相應(yīng)物理模型，在此基礎(chǔ)上采用與地鐵站火災(zāi)數(shù)值模擬實驗相同的方法來進行對比。總計算時間為120s，截取部分計算區(qū)域，對數(shù)值實驗結(jié)果進行分析。

4.2數(shù)值實驗結(jié)果

自然通風(fēng)情況下，120s時火源中心X=2.5m截面處，可以看到強烈對流產(chǎn)生的溫度分布和速度矢量分布，基本以火源為中心成對稱分布；熱煙氣向上運動，并在隧道頂部形成熱煙氣層向隧道兩端擴散，最大擴散速度為7m/s；由于煙羽流的作用，煙氣層下方的空氣向著火源方向運動，速度較慢小于2m/s，與文獻[5]的計算值相似。

4.3數(shù)值模擬實驗與實驗測試對比

圖2為距火源10m處，隧道中心線上的溫度分布和速度分布與文獻[5]的計算值和Kumar實測值的對比。速度分布相近，趨勢相同；溫度分布稍有偏差，主要是沒有考慮輻射傳熱而引起的溫度稍高，但總體趨勢相同。綜上所述，數(shù)值模擬方法可以較好的描述火災(zāi)的發(fā)展過程。

5.總結(jié)與展望

在地鐵站火災(zāi)數(shù)值模擬研究中，地鐵站物理模型的選擇，模擬時邊界條件的設(shè)置，數(shù)值模擬方法是否正確，都對實驗結(jié)果造成一定的影響。物理模型的選擇可以根據(jù)已有或在建的項目來構(gòu)建，邊界條件的設(shè)置可以根據(jù)相關(guān)研究來確定，數(shù)值模擬方法是否正確就要和實際燃燒實驗來進行對比。

研究地鐵站發(fā)生火災(zāi)時，不同防排煙系統(tǒng)模式下，溫度場、速度場、以及火災(zāi)煙氣蔓延和排煙效果，以及列車活塞風(fēng)對地鐵站火災(zāi)溫度場、速度場以及煙氣蔓延的影響，對火災(zāi)中人員安全疏散起到了指導(dǎo)作用。

參考文獻：

[1]徐碩.地鐵火災(zāi)煙氣蔓延研究初探.中國科技縱橫，2011年7月（下）總第122期.

[2]李啟榮，黎少其.地鐵火災(zāi)系統(tǒng)保障研究.香港：城市軌道交通研究，2001.

[3]楊昀，曹麗英.地鐵火災(zāi)場景設(shè)計探討.自然災(zāi)害學(xué)報，2006，Vol.15（4）：122～125.