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

數(shù)值分析范文第1篇

關(guān)鍵詞：數(shù)值分析；數(shù)值實(shí)驗(yàn)；數(shù)學(xué)建模

數(shù)值分析是一門(mén)與計(jì)算機(jī)使用密切結(jié)合的、實(shí)用性很強(qiáng)的課程。它內(nèi)容豐富，涉及數(shù)學(xué)分析、代數(shù)、方程和泛函分析等諸多學(xué)科，研究方法深刻，有自身嚴(yán)密的科學(xué)系統(tǒng)。科學(xué)與工程中的數(shù)值計(jì)算已經(jīng)成為各門(mén)自然科學(xué)和工程技術(shù)科學(xué)的一種重要手段，成為實(shí)驗(yàn)和理論并列的一個(gè)不可缺少的環(huán)節(jié)［1］。所以數(shù)值分析既是一個(gè)基礎(chǔ)性的，同時(shí)也是一個(gè)應(yīng)用性的數(shù)學(xué)學(xué)科，與其他學(xué)科的聯(lián)系十分緊密。那么在平時(shí)的教學(xué)中，如何取得良好的教學(xué)效果呢？本文從以下幾個(gè)方面進(jìn)行探討。

一、數(shù)值分析課程的教學(xué)特點(diǎn)

與其它純數(shù)學(xué)理論課程相比，數(shù)值分析除了具備數(shù)學(xué)的高度抽象性與嚴(yán)密科學(xué)性的特點(diǎn)之外，又有應(yīng)用的廣泛性與實(shí)際試驗(yàn)的高度技術(shù)性的特點(diǎn)。具體來(lái)說(shuō)，這門(mén)課程具有以下的教學(xué)特點(diǎn)：

1.知識(shí)面跨度大［2］

數(shù)值分析是數(shù)學(xué)與應(yīng)用數(shù)學(xué)、信息與計(jì)算科學(xué)和統(tǒng)計(jì)學(xué)專(zhuān)業(yè)的必修課程，它廣泛運(yùn)用多門(mén)數(shù)學(xué)學(xué)科的知識(shí)，內(nèi)容包括數(shù)值逼近、數(shù)值積分、線(xiàn)性代數(shù)方程組的直接解法和迭代方法、非線(xiàn)性方程組的計(jì)算方法、矩陣特征值與特征向量的計(jì)算、常微分方程數(shù)值計(jì)算等，涉及數(shù)學(xué)分析、代數(shù)學(xué)、微分方程、泛函分析等眾多數(shù)學(xué)理論。

2.有可靠的理論分析［2］

能任意逼近并達(dá)到精度要求，對(duì)近似算法要保證收斂性和數(shù)值穩(wěn)定性，還要對(duì)誤差進(jìn)行分析。

3.注重理論與應(yīng)用的結(jié)合

與傳統(tǒng)數(shù)學(xué)課程強(qiáng)調(diào)理論分析和邏輯推導(dǎo)不同，數(shù)值分析課程更注重運(yùn)用這些理論構(gòu)造適合計(jì)算機(jī)執(zhí)行的數(shù)值方法，要根據(jù)計(jì)算機(jī)特點(diǎn)提供實(shí)際可行的有效算法。數(shù)值分析主要研究那些在理論上有解而用手工無(wú)法計(jì)算、必需借助計(jì)算機(jī)求解的數(shù)學(xué)問(wèn)題。它的許多理論與方法本身并不是數(shù)學(xué)學(xué)科的產(chǎn)物，而是以“計(jì)算”為目標(biāo)發(fā)展起來(lái)的。

二、教學(xué)體會(huì)

針對(duì)數(shù)值分析課程的特點(diǎn)，筆者認(rèn)為在教學(xué)中應(yīng)注重以下幾個(gè)方面：

1.教學(xué)方法上注重?cái)?shù)值思想的傳授

計(jì)算方法這門(mén)課程最主要闡述的思想就是“近似計(jì)算”的思想。在實(shí)際的計(jì)算過(guò)程中，有許多問(wèn)題的計(jì)算量非常龐大，簡(jiǎn)單的筆算費(fèi)時(shí)費(fèi)力，借助計(jì)算機(jī)可以快速解決這些問(wèn)題。但由于計(jì)算機(jī)本身位數(shù)的限制，以及其它誤差影響，只能進(jìn)行近似計(jì)算。

（1）“誤差分析”思想。由于是近似計(jì)算，那么就存在一定的誤差，所以在計(jì)算過(guò)程中要分析誤差、控制誤差和比較誤差，只有控制好誤差才能找到好的近似值。誤差是衡量近似計(jì)算結(jié)果好壞的一個(gè)標(biāo)準(zhǔn)，例如，在求解線(xiàn)性方程組直接法時(shí)，通過(guò)誤差分析可以確定方程組是病態(tài)的還是良態(tài)的，只有良態(tài)的方程組才能保證解的準(zhǔn)確性。通過(guò)分析誤差可以判斷算法的穩(wěn)定性、收斂性及收斂速度。由此可見(jiàn)誤差分析是非常重要的。

（2）逼近和近似思想。函數(shù)逼近是數(shù)值分析方法中的主要內(nèi)容之一，許多數(shù)值方法都依賴(lài)于函數(shù)逼近的思想。如，各種插值方法、數(shù)值微分和數(shù)值積分、微分方程數(shù)值解等等。函數(shù)逼近中常常采取的各種近似，利用插值函數(shù)對(duì)數(shù)值近似處理，讓學(xué)生意識(shí)到數(shù)值分析課程不是在簡(jiǎn)單地做數(shù)學(xué)練習(xí)，而是在訓(xùn)練通過(guò)對(duì)原問(wèn)題的分析，如何利用已有的數(shù)學(xué)知識(shí)和工具去逼近和近似原來(lái)問(wèn)題的解。逼近和近似思想作為一種全新的思維方式，它使學(xué)生認(rèn)識(shí)到：不能解析或精確求解問(wèn)題并不可怕，可怕的是不會(huì)和不敢利用已學(xué)數(shù)學(xué)知識(shí)去近似、簡(jiǎn)化原來(lái)的問(wèn)題，從而獲得原來(lái)問(wèn)題的近似解答。

（3）“離散化”思想［6］。把求連續(xù)變量問(wèn)題轉(zhuǎn)化為求離散變量問(wèn)題，稱(chēng)為“離散化”。一個(gè)連續(xù)的數(shù)學(xué)問(wèn)題要實(shí)現(xiàn)上機(jī)計(jì)算，必須先進(jìn)行離散化。在工程計(jì)算中，常常需要求解連續(xù)性問(wèn)題，比如求微分方程的解。一般而言，微分方程很難找到解析解，所以數(shù)值求解微分方程是計(jì)算方法中的一個(gè)重要的內(nèi)容。數(shù)值求解微分方程并不是依靠計(jì)算機(jī)給出微分方程的解析形式，而是依靠它近似給出微分方程在指定點(diǎn)的函數(shù)值。在引人離散化思想對(duì)問(wèn)題離散后，可以采用各種數(shù)值方法來(lái)求解各點(diǎn)函數(shù)的值。通過(guò)離散化思想，原來(lái)的連續(xù)性問(wèn)題變成了一個(gè)離散問(wèn)題。離散化思想是數(shù)值計(jì)算的一個(gè)基本思想，現(xiàn)有的數(shù)值計(jì)算，幾乎完全依賴(lài)于對(duì)問(wèn)題的離散化解決。離散方法一直是數(shù)值分析研究中一個(gè)很重要的方面。

（4）“迭代”思想［5］。迭代是計(jì)算機(jī)中重要的概念，也是數(shù)值分析方法中的重要的概念。在數(shù)學(xué)建模過(guò)程中，對(duì)結(jié)果可能性的猜測(cè)可以在很大程度上幫助我們?cè)诮７较蛏线M(jìn)行選擇，使我們少走許多彎路。由于迭代方法大都只有有限的收斂區(qū)間，所以如何利用已有的信息對(duì)解進(jìn)行猜測(cè)是很重要的一點(diǎn)，這依賴(lài)于學(xué)生在實(shí)踐中能夠綜合運(yùn)用數(shù)學(xué)分析理論和各種方法的經(jīng)驗(yàn)。許多連續(xù)問(wèn)題在轉(zhuǎn)化為離散問(wèn)題后，利用迭代法可以求解離散問(wèn)題。

2.多媒體課件與板書(shū)相結(jié)合的教學(xué)手段［3］

使用多媒體教學(xué)方法，能增大教學(xué)容量，提高教學(xué)效率，有利于解決重點(diǎn)和難點(diǎn)問(wèn)題。多媒體教學(xué)可以在一定程度上突破時(shí)間和空間的限制，充實(shí)直觀內(nèi)容，能夠較徹底地分解知識(shí)技能信息的復(fù)雜度，減少信息在大腦中從形象到抽象，再由抽象到形象的加工轉(zhuǎn)換過(guò)程，充分傳達(dá)教學(xué)意圖，并可以通過(guò)計(jì)算機(jī)的豐富表現(xiàn)手段突出教學(xué)重點(diǎn)。如，龍格現(xiàn)象可以用屏幕動(dòng)態(tài)的顯示在哪個(gè)區(qū)間收斂，使用多媒體教學(xué)可以幫助教師在課堂上根據(jù)學(xué)生的信息反饋，進(jìn)行現(xiàn)場(chǎng)分析和答疑，以人機(jī)對(duì)話(huà)方式靈活方便地進(jìn)行啟發(fā)式教學(xué)。同時(shí)，精彩的多媒體課件也能激發(fā)學(xué)生的興趣，提高學(xué)生的主動(dòng)性。

數(shù)值分析范文第2篇

關(guān)鍵詞：數(shù)值分析 課程改革 教與學(xué)換位思考

中圖分類(lèi)號(hào)：G642 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1673-9795（2013）08（a）-0027-01

1 明確研究對(duì)象和任務(wù)，確保教與學(xué)有的放矢

數(shù)值分析是研究現(xiàn)代科學(xué)與工程計(jì)算中各種數(shù)學(xué)問(wèn)題求解的數(shù)值計(jì)算方法，關(guān)于數(shù)值分析的國(guó)內(nèi)外教材中涉及到的數(shù)學(xué)問(wèn)題大致相同，主要包含線(xiàn)性代數(shù)的方程組求解，矩陣特征值特征向量的計(jì)算以及微積分和微分方程求解等，至于章節(jié)的順序安排主要有兩種，詳見(jiàn)參考文獻(xiàn)[1～2]。區(qū)別于解析解法，數(shù)值計(jì)算方法是結(jié)合計(jì)算機(jī)求解數(shù)值問(wèn)題近似解的方法。所謂的數(shù)值問(wèn)題[2]是指有限個(gè)輸入數(shù)據(jù)（問(wèn)題的自變量、原始數(shù)據(jù)）與有限個(gè)輸出數(shù)據(jù)（待求解數(shù)據(jù)）之間函數(shù)關(guān)系的一個(gè)明確無(wú)歧義的描述。上述的數(shù)學(xué)問(wèn)題不一定是數(shù)值問(wèn)題，如求解一階常微分方程初值問(wèn)題：

，

要求得到定義于區(qū)間[0，1]的函數(shù)解析表達(dá)式，這實(shí)際上要求無(wú)窮多個(gè)輸出，因而它不是數(shù)值問(wèn)題。但利用歐拉方法選取合適的步長(zhǎng)，求得個(gè)點(diǎn)處的函數(shù)值的近似值時(shí)，便成為了一個(gè)數(shù)值問(wèn)題。

數(shù)值分析的研究對(duì)象就是數(shù)值問(wèn)題，其主要任務(wù)是構(gòu)造求解數(shù)值問(wèn)題近似解的方法，并且利用計(jì)算機(jī)編程實(shí)現(xiàn)算法，最后通過(guò)誤差分析判定算法的收斂性和穩(wěn)定性。

2 合理利用計(jì)算機(jī)軟件，讓數(shù)值分析不再抽象

目前有很多《數(shù)值分析》[1]教材嘗試了和計(jì)算機(jī)應(yīng)用的結(jié)合，加入了計(jì)算機(jī)編程內(nèi)容，融進(jìn)了利用C語(yǔ)言或MATLAB軟件實(shí)踐計(jì)算方法理論的內(nèi)容，也有不少同仁提出了如何搞好數(shù)值實(shí)驗(yàn)教學(xué)的建議[4]，對(duì)數(shù)值分析理論與實(shí)踐并重的改革起到了推進(jìn)作用。但是不同專(zhuān)業(yè)、不同層次的學(xué)生對(duì)數(shù)值分析的教學(xué)要求也不同。一般而言，學(xué)習(xí)數(shù)值分析這門(mén)課程的學(xué)生，可以區(qū)分為兩部分：其一是數(shù)學(xué)和計(jì)算機(jī)專(zhuān)業(yè)的，他們學(xué)習(xí)的主要目標(biāo)是研究和創(chuàng)造算法；而另一部分是工科專(zhuān)業(yè)的，他們學(xué)習(xí)的目的主要是使用算法。由于學(xué)習(xí)目標(biāo)不同，所以教材的內(nèi)容、體系及側(cè)重點(diǎn)也應(yīng)該有相應(yīng)的差別。對(duì)于數(shù)學(xué)和計(jì)算機(jī)專(zhuān)業(yè)的學(xué)生，大都已經(jīng)系統(tǒng)的學(xué)習(xí)過(guò)C語(yǔ)言，具備一定的獨(dú)立編程能力，因此，他們的學(xué)習(xí)目標(biāo)應(yīng)該放在對(duì)數(shù)學(xué)算法的理解和優(yōu)化上，如果只是讓他們通過(guò)簡(jiǎn)單的調(diào)用MATLAB軟件中現(xiàn)成的函數(shù)命令來(lái)實(shí)現(xiàn)算法，對(duì)他們的數(shù)學(xué)理解能力和編程水平的提高沒(méi)有多大幫助；而對(duì)于其他大多數(shù)工科專(zhuān)業(yè)的學(xué)生來(lái)說(shuō)，不具備一定的編程能力，若是讓他們用C語(yǔ)言來(lái)實(shí)現(xiàn)算法有一定的難度，此時(shí)MATLAB卻可以發(fā)揮其處理矩陣輕松，編程簡(jiǎn)潔以及容易入手等優(yōu)勢(shì)。這也是MATLAB適用于科學(xué)和工程計(jì)算的原因。同時(shí)，MATLAB軟件還具備強(qiáng)大的繪圖和可視化功能，可以讓學(xué)生根據(jù)幾何圖形更直觀的觀察計(jì)算結(jié)果?？梢赃_(dá)到“一圖勝千言”的效果，讓課堂教學(xué)不再只是冗長(zhǎng)的數(shù)學(xué)公式和抽象的推理證明，進(jìn)而調(diào)動(dòng)了學(xué)生的學(xué)習(xí)興趣。

搞好數(shù)值實(shí)驗(yàn)教學(xué)的重要意義在于詳細(xì)地講解數(shù)值計(jì)算方法的基本理論后，利用數(shù)學(xué)軟件將實(shí)際算例的幾種不同算法，符號(hào)解法及其圖形和數(shù)表進(jìn)行實(shí)際計(jì)算效果對(duì)比，讓學(xué)生從理論和幾何直觀上觀察同一種算法取不同參數(shù)時(shí)對(duì)計(jì)算結(jié)果的影響，比較和分析不同算法在計(jì)算同一問(wèn)題時(shí)的誤差大小，從而篩選出最佳算法。通過(guò)這樣相對(duì)完整的科學(xué)計(jì)算過(guò)程，使學(xué)生切實(shí)體會(huì)到學(xué)有所用的樂(lè)趣。通過(guò)教學(xué)實(shí)踐驗(yàn)證效果確實(shí)不錯(cuò)，越來(lái)越多的學(xué)生開(kāi)始喜歡數(shù)值分析這門(mén)課程了。

3 活躍課堂教學(xué)氣氛，培養(yǎng)學(xué)生學(xué)習(xí)主動(dòng)性

許多大學(xué)生談數(shù)學(xué)色變，在數(shù)學(xué)課堂上也不是很積極，通過(guò)觀察及與學(xué)生的交流，我發(fā)現(xiàn)學(xué)生對(duì)于學(xué)習(xí)目標(biāo)要求很低，考試能及格就滿(mǎn)足了。還有的學(xué)生是認(rèn)為數(shù)值分析這門(mén)課程很難，遇到聽(tīng)不懂的便望而卻步了。

對(duì)于上述學(xué)生的困惑，結(jié)合自身的學(xué)習(xí)和工作經(jīng)歷，我給學(xué)生反饋的建議是：教與學(xué)的換位思考。教師講課是積極主動(dòng)的，而不少學(xué)生聽(tīng)課是消極被動(dòng)的。學(xué)生聽(tīng)課只要聽(tīng)明白就可以了，甚至當(dāng)堂聽(tīng)不懂也沒(méi)關(guān)系，課下還有自學(xué)跟復(fù)習(xí)的機(jī)會(huì)。老師講課不可以講不明白，在站上講臺(tái)之前就必須備好課，不光要清楚所講的知識(shí)點(diǎn)，還得清楚如何講才能讓學(xué)生更容易理解。這就需要老師具備更寬廣的知識(shí)面，只有這樣才能做到融會(huì)貫通，深入淺出。實(shí)際上，老師的備課過(guò)程也是一個(gè)學(xué)習(xí)過(guò)程，只不過(guò)較之學(xué)生要更積極更主動(dòng)些，并且學(xué)習(xí)目標(biāo)的定位也更高一些。綜上所述，如果學(xué)生能以老師的標(biāo)準(zhǔn)來(lái)要求自己，又怎么會(huì)學(xué)不好呢！我做過(guò)這樣的教學(xué)換位試驗(yàn)，讓一位曾經(jīng)覺(jué)得數(shù)值分析很難的學(xué)生提前準(zhǔn)備一周，然后到講臺(tái)上講解一個(gè)知識(shí)點(diǎn)。雖然當(dāng)他第一次站上講臺(tái)，面對(duì)那么多同學(xué)說(shuō)話(huà)的時(shí)候有些緊張，個(gè)別地方語(yǔ)言表達(dá)也不夠嚴(yán)密，但是可以看出來(lái)他上課前的確下功夫準(zhǔn)備了，最起碼他自己對(duì)該知識(shí)點(diǎn)基本上掌握了，再結(jié)合老師的解釋?zhuān)蜁?huì)理解的更深刻。正如《為學(xué)》中所說(shuō)：“學(xué)之，則難者亦易矣；不學(xué)，則易者亦難矣”。從此該學(xué)生聽(tīng)課也認(rèn)真了，還經(jīng)常跟老師探討問(wèn)題，期末考試成績(jī)也很理想。我相信對(duì)他來(lái)說(shuō)，這肯定是一次受益匪淺終生難忘的學(xué)習(xí)經(jīng)歷。由此可見(jiàn)找回自信端正態(tài)度對(duì)學(xué)習(xí)乃至工作和生活的重要性。

為師者，不僅要傳其道還要解其惑。因此，我經(jīng)常跟學(xué)生探討學(xué)習(xí)心得與體會(huì)，及時(shí)了解他們?cè)趯W(xué)習(xí)過(guò)程中的需求與困惑，聽(tīng)取他們關(guān)于教與學(xué)的意見(jiàn)與建議。同時(shí)我也會(huì)結(jié)合自身的體會(huì)以及諸多同仁的經(jīng)驗(yàn)給予學(xué)生中肯的建議和滿(mǎn)意的答復(fù)。要想讓學(xué)生喜歡老師所講的課，首先要讓學(xué)生喜歡講課的老師。所以，要想成為一名受學(xué)生愛(ài)戴的優(yōu)秀老師，不光要提升自身的教學(xué)科研水平，還要站在學(xué)生的角度設(shè)身處地的為學(xué)生考慮，做學(xué)生的朋友，盡可能讓枯燥的課堂教學(xué)變得如好友聊天般輕松愉快！

4 結(jié)語(yǔ)

數(shù)值分析作為數(shù)學(xué)基礎(chǔ)課程的一部分，它與高等數(shù)學(xué)、線(xiàn)性代數(shù)等大部分?jǐn)?shù)學(xué)課程又有很大的不同。數(shù)值分析既注重理論分析又注重實(shí)踐應(yīng)用，它強(qiáng)調(diào)結(jié)合計(jì)算機(jī)解決實(shí)際工程計(jì)算問(wèn)題，所以它對(duì)數(shù)學(xué)建模課程及競(jìng)賽，數(shù)學(xué)軟件和計(jì)算機(jī)編程語(yǔ)言的學(xué)習(xí)與應(yīng)用，乃至對(duì)學(xué)生將來(lái)的工作和科研都有著重要而又深遠(yuǎn)的影響，因此，數(shù)值分析的課程建設(shè)和教學(xué)改革是一項(xiàng)值得長(zhǎng)期關(guān)注和重視的項(xiàng)目。

參考文獻(xiàn)

數(shù)值分析范文第3篇

關(guān)鍵詞：空心墩；OpenSees；纖維梁柱單元；縱筋拔出變形；滯回曲線(xiàn)

中圖分類(lèi)號(hào)：TV331文獻(xiàn)標(biāo)示碼：A

Numerical seismic analysis model for hollow reinforced concrete bridge piers

YANG Chun-xi

Tianjincommunicationarchitecturedesigninstitute

Abstract: Hollow reinforced concrete bridge piers are extensively used in highway and railway bridges. Simulation of the seismic behavior of hollow bridge piers is important to ensure the seismic safety of bridges. Based on fiber beam-column element and zero length rotation spring element, three hollow reinforced concrete bridge piers are modeled by using OpenSees software, in which the flexural and bond-slip deformation of the piers are considered. The simulated hysteretic curves are compared with test results. The results show that the simulated results agreed well with test results.

Key Words: Hollow bridge pier; OpenSees; fiber beam-column element; bond-slip deformation; hysteretic curves

鋼筋混凝土空心墩被廣泛應(yīng)用于我國(guó)鐵路、公路橋梁中，而國(guó)內(nèi)外對(duì)空心墩抗震性能認(rèn)識(shí)比較缺乏，開(kāi)展鋼筋混凝土空心墩的抗震數(shù)值模擬工作具有重要的工程意義[1]。孫治國(guó)基于纖維梁柱單元建立了鋼筋混凝土空心墩滯回分析模型，詳細(xì)討論了縱筋配筋、壁厚、混凝土強(qiáng)度、剪跨比等因素對(duì)空心墩變形能力的影響。在此基礎(chǔ)上設(shè)計(jì)了2個(gè)矩形薄壁空心墩試件，分別進(jìn)行定軸力和變軸力下的擬靜力試驗(yàn)，并指出使用修正的壓力場(chǎng)理論(Modified Compression Field Theory, MCFT)計(jì)算的薄壁空心墩抗剪強(qiáng)度最為準(zhǔn)確[2-3]。同時(shí)，同濟(jì)大學(xué)[4]、東南大學(xué)[5]、北京工業(yè)大學(xué)[6]、長(zhǎng)安大學(xué)[7]等單位也針對(duì)空心墩抗震問(wèn)題開(kāi)展了試驗(yàn)研究工作，對(duì)認(rèn)識(shí)空心墩抗震能力提供了基礎(chǔ)。

注意到目前國(guó)內(nèi)外對(duì)空心墩抗震開(kāi)展的數(shù)值模擬工作較少，本文基于OpenSees數(shù)值分析平臺(tái)，考慮空心墩的彎曲變形與縱筋拔出變形，建立了3個(gè)空心墩數(shù)值分析模型，并將模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行了對(duì)比，驗(yàn)證模型正確性。

1試驗(yàn)介紹

選擇3個(gè)空心墩抗震擬靜力試驗(yàn)結(jié)果，并以此為依據(jù)，建立了空心墩抗震數(shù)值分析模型，并通過(guò)與試驗(yàn)結(jié)果的對(duì)比驗(yàn)證模型準(zhǔn)確性。

第1個(gè)試件選取Zahn等[8]完成的空心墩抗震擬靜力試驗(yàn)中的Unit11試件，高度為1600mm，外徑為400mm，內(nèi)徑為212mm，剪跨比為4.0，軸壓比為0.08，縱筋配筋率為3.56%。第2個(gè)試件選取孫治國(guó)等[3]完成的薄壁空心墩試驗(yàn)中的定軸力試件，橋墩高度為4000mm，截面尺寸為1000×890mm，空心部分截面尺寸為860×750mm，軸壓比為0.2，縱筋配筋率為1%。第3個(gè)試件選取Pinto等[9]完成的大比例尺矩形空心墩中的矮墩試件，截面尺寸為2740 mm×1020 mm，空心部分截面尺寸為2320 mm×860mm，軸壓比為0.09，縱筋配筋率為0.4％。需要強(qiáng)調(diào)，所有橋墩試件最終均發(fā)生彎曲破壞，圖1為各試件截面及配筋形式。

(a) Unit 11試件 (b) 定軸力試件

(c) 矮墩試件

圖1 空心墩截面及配筋形式

2 模型建立

2.1 混凝土本構(gòu)及鋼筋材料模型

混凝土數(shù)值模型采用OpenSees中的Concrete01，該材料基于Kent-Scott-Park混凝土單軸受壓應(yīng)力-應(yīng)變關(guān)系，如圖2所示，k為約束效應(yīng)系數(shù)，ε0為峰值應(yīng)變，εu為極限應(yīng)變，fc’為混凝土抗壓強(qiáng)度。鋼筋材料采用OpenSees中的Steel02鋼筋模型，其應(yīng)力應(yīng)變模型是基于Giuffre-Menegotto-Pinto模型，骨架曲線(xiàn)為雙折線(xiàn)形式。

圖2 混凝土應(yīng)力-應(yīng)變關(guān)系

縱筋在底座中的拔出采用Zhao Jian提出的Bond_SP01[10]材料模擬，骨架曲線(xiàn)如圖3所示，其中E為鋼筋彈性模量，fy為鋼筋屈服應(yīng)力，Sy為屈服滑移量，fu 為極限應(yīng)力，Su為極限滑移量，b為剛度折減系數(shù)。Sy計(jì)算公式如下：

(1)

式中，db為鋼筋直徑，α是局部粘結(jié)滑移參數(shù)，取0.4。fc’為混凝土強(qiáng)度。另根據(jù)經(jīng)驗(yàn)計(jì)算可得，Su=(30~40)Sy，b取(0.3~0.5)，R取(0.5~1.0)。

圖3 Bond_SP01鋼筋應(yīng)力-滑移骨架曲線(xiàn)

2.2 數(shù)值分析模型

基于OpenSees中的纖維梁柱單元和零長(zhǎng)度轉(zhuǎn)動(dòng)彈簧單元建立數(shù)值分析模型，如圖4所示。纖維梁柱單元用于模擬橋墩的非線(xiàn)性彎曲變形，將Bond_SP01材料賦予零長(zhǎng)度轉(zhuǎn)動(dòng)彈簧單元，用于模擬底部縱筋拔出變形。纖維梁柱單元與零長(zhǎng)度轉(zhuǎn)動(dòng)彈簧單元基于相同的纖維劃分，唯一的區(qū)別是非線(xiàn)性梁柱單元截面內(nèi)的鋼筋材料使用steel02，而零長(zhǎng)度轉(zhuǎn)動(dòng)彈簧單元截面內(nèi)的鋼筋材料使用Bond_SP01。

圖4 數(shù)值分析模型

3 滯回曲線(xiàn)對(duì)比

數(shù)值模型考慮了彎曲變形和底部縱筋的拔出變形，圖5為模擬得到的Unit11試件、定軸力試件、矮墩試件墩頂滯回曲線(xiàn)以及與試驗(yàn)結(jié)果的對(duì)比情況。

首先分析各試件承載力的對(duì)比情況，可以看出，模擬的空心墩極限荷載與試驗(yàn)結(jié)果吻合很好。即數(shù)值模型在承載力角度對(duì)空心墩抗震試驗(yàn)進(jìn)行了很好的模擬。

然后分析模擬得到的滯回曲線(xiàn)初始剛度、卸載剛度及與試驗(yàn)結(jié)果的對(duì)比情況?？砂l(fā)現(xiàn)所建模型也很好的模擬了各空心墩的加載和卸載段的剛度。

最后考慮各模型模擬得到的殘余位移及與試驗(yàn)結(jié)果的對(duì)比，可發(fā)現(xiàn)Unit11試件、定軸力試件模擬的殘余位移與試驗(yàn)結(jié)果基本吻合。矮墩試件模擬的滯回曲線(xiàn)殘余位移小于試驗(yàn)結(jié)果，但總體在可接受范圍內(nèi)。

綜上可以看出，本文所建數(shù)值模型準(zhǔn)確，可對(duì)空心墩地震反應(yīng)進(jìn)行較為準(zhǔn)確的模擬分析。

(a) Unit11試件

(b) 定軸力試件

(c) 矮墩試件

圖5 模擬與試驗(yàn)滯回曲線(xiàn)的對(duì)比

4 結(jié)論

基于OpenSees平臺(tái)建立了3個(gè)考慮彎曲變形和縱筋拔出變形的空心橋墩抗震數(shù)值分析模型，并與試驗(yàn)滯回曲線(xiàn)進(jìn)行了對(duì)比?？傮w來(lái)看，模擬滯回曲線(xiàn)與試驗(yàn)結(jié)果吻合較好，表明模型建立正確，并具有較高的模擬精度。

參考文獻(xiàn):

[1] 孫治國(guó),王東升,李宏男,等.鋼筋混凝土空心橋墩應(yīng)用及抗震性能研究綜述[J].交通運(yùn)輸工程學(xué)報(bào),2013,13(3):22-32.

[2] 孫治國(guó),郭迅,王東升,等.鋼筋混凝土空心墩延性變形能力分析[J].鐵道學(xué)報(bào).2012,34(1):91-96.

[3] 孫治國(guó),王東升,郭迅,等.地震作用下RC薄壁空心墩抗剪強(qiáng)度比較研究[J].土木工程學(xué)報(bào).2013,46(12):81-89.

[4] 羅征,李建中.低周往復(fù)荷載下空心矩形墩抗震性能試驗(yàn)研究[J].振動(dòng)與沖擊,2013,32(8):183-188.

[5] 宗周紅,夏堅(jiān),徐綽然.橋梁高墩抗震研究現(xiàn)狀及展望[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,43(2):445-452.

[6] 杜修力,陳明琦,韓強(qiáng).鋼筋混凝土空心橋墩抗震性能試驗(yàn)研究[J].振動(dòng)與沖擊,2011,30(11):254-259.

[7] 崔海琴,賀拴海,宋一凡.空心矩形薄壁墩延性抗震性能試驗(yàn)[J].公路交通科技,2010,27(6):58-63.

[8] ZAHN F A. Design of reinforced concrete bridge columns for strength and ductility [D]. Christchurch; University of Canterbury, 1985.

數(shù)值分析范文第4篇

【關(guān)鍵詞】巖土工程；連續(xù)介質(zhì)力學(xué)；數(shù)值分析；本構(gòu)理論；本構(gòu)模型

1 引 言

本文首先介紹筆者對(duì)我國(guó)巖土工程數(shù)值分析現(xiàn)狀的調(diào)查結(jié)果，然后就巖土工程分析中的關(guān)鍵問(wèn)題，如何發(fā)展巖土本構(gòu)理論和數(shù)值分析在巖土工程分析中的地位這 3 個(gè)問(wèn)題提出粗淺的看法。

2巖土工程分析中的關(guān)鍵問(wèn)題

巖土工程分析中人們常常將用簡(jiǎn)化的物理模型去描述復(fù)雜的工程問(wèn)題，再將其轉(zhuǎn)化為數(shù)學(xué)問(wèn)題并用數(shù)學(xué)方法求解。一個(gè)很典型的例子是，飽和軟黏土地基大面積堆載作用下的沉降問(wèn)題被簡(jiǎn)化為T(mén)erzaghi 一維固結(jié)物理模型，再轉(zhuǎn)化為 Terzaghi 固結(jié)方程求解。采用連續(xù)介質(zhì)力學(xué)模型求解工程問(wèn)題一般包括

下述方程：①運(yùn)動(dòng)微分方程式（包括動(dòng)力和靜力分析兩大類(lèi)）；②幾何方程（包括小應(yīng)變分析和大應(yīng)變分析兩大類(lèi)）；③本構(gòu)方程（即力學(xué)本構(gòu)方程）。對(duì)一具體工程問(wèn)題，根據(jù)具體的邊界條件和初始條件求解上述方程即可得到解答，對(duì)復(fù)雜的工程問(wèn)題，一般需采用數(shù)值分析法求解。對(duì)不同的工程問(wèn)題采用連續(xù)介質(zhì)力學(xué)模型求解，所用的運(yùn)動(dòng)微分方程式和幾何方程是相同的，不同的是本構(gòu)方程、邊界條件和初始條件。當(dāng)材料為線(xiàn)性彈性體，本構(gòu)方程為廣義虎克定律。將巖土材料視為多相體，采用連續(xù)介質(zhì)力學(xué)模型分析巖土工程問(wèn)題一般包括下述方程[2]：①運(yùn)動(dòng)微分方程式（包括動(dòng)力和靜力分析兩大類(lèi)）；②總應(yīng)力= 有效應(yīng)力+ 孔隙壓力（有效應(yīng)力原理）；③連續(xù)方程（總體積變化為各相體積變化之和）；④幾何方程，包括小應(yīng)變分析和大應(yīng)變分析兩大類(lèi)；⑤本構(gòu)方程，即力學(xué)和滲流本構(gòu)方程。

將多相體與單相體比較，基本方程多了 2 個(gè)，即有效應(yīng)力原理和連續(xù)方程，且本構(gòu)方程中多了滲流本構(gòu)方程。對(duì)不同的巖土工程問(wèn)題，基本方程中運(yùn)動(dòng)微分方程式、有效應(yīng)力原理、連續(xù)方程和幾何方程的表達(dá)式是相同的，不同的是本構(gòu)方程。對(duì)一具體巖土工程問(wèn)題，根據(jù)具體的邊界條件和初始條件求解上述方程即可得到解答，一般需采用數(shù)值分析法求解。從上面分析可知，采用連續(xù)介質(zhì)力學(xué)模型分析不同的巖土工程問(wèn)題時(shí)，不同的是本構(gòu)模型、邊界條件和初始條件。對(duì)一個(gè)具體的巖土工程問(wèn)題，邊界條件和初始條件是容易確定的，而巖土的應(yīng)力-應(yīng)變關(guān)系十分復(fù)雜，采用的本構(gòu)模型及參數(shù)對(duì)計(jì)算結(jié)果影響極大。

采用連續(xù)介質(zhì)力學(xué)模型分析巖土工程問(wèn)題一般需采用數(shù)值分析法求解，有限單元法對(duì)各種邊界條件和初始條件，采用的各類(lèi)本構(gòu)方程都有較大的適應(yīng)性。土的應(yīng)力-應(yīng)變關(guān)系十分復(fù)雜，自 Roscoe 和他的學(xué)生建立劍橋模型至今已近半個(gè)世紀(jì)，理論上已提出數(shù)百個(gè)本構(gòu)方程，但得到工程應(yīng)用認(rèn)可的極少，或者說(shuō)還沒(méi)有。從這個(gè)角度講，采用連續(xù)介質(zhì)力學(xué)模型求解巖土工程問(wèn)題的關(guān)鍵問(wèn)題是如何建立巖土材料的工程實(shí)用本構(gòu)方程。

3 如何發(fā)展巖土本構(gòu)理論的思考

Janbu 認(rèn)為，反映作用與效應(yīng)之間的關(guān)系稱(chēng)為本構(gòu)關(guān)系，力學(xué)中的虎克定律、電學(xué)中的歐姆定律、滲流學(xué)中的達(dá)西定律等反映的都是最簡(jiǎn)單的本構(gòu)關(guān)系。巖土是自然、歷史的產(chǎn)物，具有下述特性：土體性質(zhì)區(qū)域性強(qiáng)，即使同一場(chǎng)地同一層土，沿深度和水平方向變化也很復(fù)雜；巖土體中的初始應(yīng)力場(chǎng)復(fù)雜且難以測(cè)定；土是多相體，一般由固相、液相和氣相三相組成，土體中的三相有時(shí)很難區(qū)分，而且處不同狀態(tài)時(shí)，土的三相之間可以相互轉(zhuǎn)化。土中水的狀態(tài)又十分復(fù)雜；土體具有結(jié)構(gòu)性，與土的礦物成分、形成歷史、應(yīng)力歷史和環(huán)境條件等因素有關(guān)，十分復(fù)雜；土的強(qiáng)度、變形和滲透特性測(cè)定困難。巖土的應(yīng)力-應(yīng)變關(guān)系與應(yīng)力路徑、加荷速率、應(yīng)力水平、成分、結(jié)構(gòu)、狀態(tài)等有關(guān)，土還具有剪脹性、各向異性等，因此，巖土體的本構(gòu)關(guān)系十分復(fù)雜。至今人們建立的土體的本構(gòu)模型類(lèi)別有彈性模型、剛塑性模型、非線(xiàn)性彈性模型、彈塑性模型、黏彈性模型、黏彈塑性模型、邊界面模型、內(nèi)時(shí)模型、多重屈服面模型、損傷模型、結(jié)構(gòu)性模型等等。已建立的本構(gòu)模型多達(dá)數(shù)百個(gè)，但得到工程師認(rèn)可的極少，或者說(shuō)還沒(méi)有。怎么走出困境？這是我們必須面對(duì)的難題。筆者認(rèn)為，對(duì)土體本構(gòu)模型研究應(yīng)分為兩大類(lèi)，科學(xué)型模型的研究和工程實(shí)用性模型的研究。科學(xué)型模型重在揭示、反映某些特殊規(guī)律，如土的剪脹性、主應(yīng)力軸旋轉(zhuǎn)的影響等。該類(lèi)模型也不能求全面，一個(gè)模型能反映一個(gè)或幾個(gè)特殊規(guī)律即為好模型。從事科學(xué)型模型研究是少數(shù)人，是科學(xué)家。工程實(shí)用性模型更不能求全面、通用，工程實(shí)用性模型應(yīng)簡(jiǎn)單、實(shí)用，參數(shù)少且易測(cè)定，能反映主要規(guī)律，能抓住主要矛盾，參數(shù)少且易測(cè)定即為好模型。工程實(shí)用性模型重在能夠應(yīng)用于具體工程分析，多數(shù)人應(yīng)從事工程實(shí)用性模型研究。研究中應(yīng)重視工程類(lèi)別（基坑工程、路堤工程、建筑工程等）、土類(lèi)（黏性土、砂土和黃土等）和區(qū)域性（上海黏土、杭州黏土和湛江黏土等）的特性的影響，如建立適用于基坑工程分析的杭州黏土本構(gòu)模型，適用于道路工程沉降分析的西黃土本構(gòu)模型和適用建筑工程沉降分析的上海黏土本構(gòu)模型等。工程實(shí)用性模型研究還要重視地區(qū)經(jīng)驗(yàn)的積累。采用考慮工程類(lèi)別、土類(lèi)和區(qū)域性特性影響的工程實(shí)用本構(gòu)模型，應(yīng)用連續(xù)介質(zhì)力學(xué)理論，并結(jié)合地區(qū)經(jīng)驗(yàn)進(jìn)行巖土工程數(shù)值分析可能是發(fā)展方向。

4數(shù)值分析在巖土工程分析中的地位

下面從巖土材料特性、巖土工程與結(jié)構(gòu)工程有限元分析誤差來(lái)源分析比較和巖土工程分析方法三方面來(lái)分析數(shù)值分析在巖土工程分析中的地位。前面已經(jīng)提到巖土材料是自然、歷史的產(chǎn)物，工程特性區(qū)域性強(qiáng)，巖土體中的初始應(yīng)力場(chǎng)復(fù)雜且難以測(cè)定，土是多相體，土體中的三相有時(shí)很難區(qū)分，土中水的狀態(tài)又十分復(fù)雜。巖土的應(yīng)力-應(yīng)變關(guān)系與應(yīng)力路徑、加荷速率、應(yīng)力水平、成份、結(jié)構(gòu)、狀態(tài)等有關(guān)，巖土體的本構(gòu)關(guān)系十分復(fù)雜。至今尚無(wú)工程師普遍認(rèn)可的工程實(shí)用的本構(gòu)模型，而采用連續(xù)介質(zhì)力學(xué)模型求解巖土工程問(wèn)題的關(guān)鍵問(wèn)題是

如何建立工程實(shí)用的巖土本構(gòu)方程，這是應(yīng)面對(duì)的現(xiàn)狀，也是考慮數(shù)值分析在巖土工程分析中的地位時(shí)必須重視的現(xiàn)實(shí)情況。

5結(jié) 論

通過(guò)對(duì)我國(guó)巖土工程數(shù)值分析現(xiàn)狀的調(diào)查研究和上述分析，筆者對(duì)巖土工程數(shù)值和解析分析的思考意見(jiàn)如下：

（1）基于對(duì)巖土工程分析對(duì)象――巖土材料特性的分析，并考慮巖土工程初始條件和邊界條件的復(fù)雜性，巖土工程分析很少能得到解析解，而目前巖土工程數(shù)值分析只能用于定性分析。所以，巖土工程設(shè)計(jì)要重視概念設(shè)計(jì)，重視巖土工程師的綜合判斷。巖土工程數(shù)值分析結(jié)果是巖土工程師在巖土工程分析過(guò)程中進(jìn)行綜合判斷的重要依據(jù)之一。

（2）自 Roscoe 和他的學(xué)生建立劍橋模型至今已近半個(gè)世紀(jì)，各國(guó)學(xué)者已提出數(shù)百個(gè)本構(gòu)方程，但得到工程應(yīng)用認(rèn)可的極少，或者說(shuō)還沒(méi)有。從這個(gè)角度講采用連續(xù)介質(zhì)力學(xué)模型求解巖土工程問(wèn)題的關(guān)鍵問(wèn)題是如何建立巖土的工程實(shí)用本構(gòu)方程。

3）巖土工程師在充分掌握分析工程地質(zhì)資料、了解土的工程性質(zhì)基礎(chǔ)上，采用合理的物理數(shù)學(xué)模型，通過(guò)多種方法進(jìn)行計(jì)算分析，然后結(jié)合工程經(jīng)驗(yàn)進(jìn)行綜合判斷，提出設(shè)計(jì)依據(jù)。在巖土工程計(jì)算分析中應(yīng)堅(jiān)持因地制宜、抓主要矛盾、宜粗不宜細(xì)、宜簡(jiǎn)不宜繁的原則。

作者簡(jiǎn)介：

數(shù)值分析范文第5篇

關(guān)鍵詞地基承載力，塑性力學(xué)上限，最優(yōu)化方法。

1前 言

地基承載力、土壓力和邊坡穩(wěn)定是土力學(xué)的3個(gè)重要領(lǐng)域。這 3個(gè)問(wèn)題都基于共同的極限平衡分析原理 ,可以采用相同的分析方法。但是 ,在長(zhǎng)期的實(shí)踐中 ,這 3個(gè)領(lǐng)域各形成了自己的體系。在地基承載力領(lǐng)域 ,目前常用的計(jì)算方法仍然是基于 Prandtl解的各種經(jīng)驗(yàn)修正公式。應(yīng)用塑性力學(xué)上下限原理 ,在建立地基承載力、土壓力和邊坡穩(wěn)定分析統(tǒng)一的理論和方法方面作了大量的工作 ,但其有關(guān)的研究一直是在變分原理基礎(chǔ)上進(jìn)行的 ,因此 ,難以擴(kuò)展到具有復(fù)雜邊界和分層土體的實(shí)際工程問(wèn)題中。曾提出一個(gè)基于滑楔破壞模式的分析方法 ,其普遍適用性還有待進(jìn)一步論證。顯然 ,只有開(kāi)發(fā)數(shù)值分析的方法 ,方可使大部分實(shí)際問(wèn)題方便地獲得解答。

近期 ,筆者在二維領(lǐng)域應(yīng)用塑性力學(xué)上限定理進(jìn)行邊坡穩(wěn)定的理論研究[4 ,5]。該方法從變形協(xié)調(diào)出發(fā) ,對(duì)于一個(gè)設(shè)定的滑裂面和斜分條模式 ,建立協(xié)調(diào)的速度場(chǎng) ,根據(jù)外力功和內(nèi)能耗散相平衡的原理確定相應(yīng)的安全系數(shù)或加載系數(shù) ,然后應(yīng)用最優(yōu)化方法 ,確定對(duì)應(yīng)于最小安全系數(shù)的那個(gè)臨界滑裂面和斜分條模式 (以下簡(jiǎn)稱(chēng)能量法 )。這一方法在精確地確定邊坡穩(wěn)定安全系數(shù)方面獲得了成功。由于地基實(shí)際上是一個(gè)坡度為零的邊坡 ,將該成果推廣到地基承載力 ,是一個(gè)十分具有理論和實(shí)用價(jià)值的課題。

2極限分析法的理論基礎(chǔ)和計(jì)算步驟

2.1上限定理的基本命題

在邊坡穩(wěn)定和地基承載力分析領(lǐng)域 ,對(duì)上限定理的描述可以用下面的命題表達(dá)(圖 1):

在塑性區(qū) Ω*,給出一個(gè)機(jī)動(dòng)可能的應(yīng)變場(chǎng)εij* ,并在滑裂面 Γ* 上給出一個(gè)相應(yīng)的速度場(chǎng) V*,那么,按照下式計(jì)算獲得的外荷載T* 將比一個(gè)包含有真實(shí)的塑性區(qū) Ω和真實(shí)的滑裂面 Γ的臨界荷載 T大或與其相等。

∫Ωσi*jεi*j dΩ +∫vdDs* = WV* + T* V*(1)

上式左邊的第一、第二項(xiàng)分別為塑性區(qū)內(nèi)和滑裂面的內(nèi)能耗散；W為塑性區(qū)土體重。因此 ,在諸多協(xié)調(diào)的位移場(chǎng)中給出最小的 T* 的那個(gè)一定最接近真實(shí)的臨界荷載 T。

在地基承載力問(wèn)題中 ,通常定義加載系數(shù) η* 為

η* = (T * - T0)/ T0 (2)

其中T0為地基的實(shí)際承受的外荷載 ,那么上限圖定理的命題具體化為尋找一個(gè)使 η* 獲得最小值 η的應(yīng)變場(chǎng)和速度場(chǎng)。

2．2計(jì)算內(nèi)能耗散

如果材料遵守莫爾-庫(kù)侖破壞準(zhǔn)則和相關(guān)聯(lián)的流動(dòng)法則，則可確認(rèn)沿滑面的速度V與滑面夾角為摩擦角φ，單位面積內(nèi)能耗散為（圖2）：

d D = (c cosφ u sinφ)v(3)

其中c為凝聚力；u為孔隙壓力；V為滑塊沿滑面的在單位荷載增量下產(chǎn)生的相對(duì)位移 ,通常稱(chēng)變形速率。

2.3計(jì)算多塊體破壞模式協(xié)調(diào)的速度場(chǎng)

對(duì)某一邊坡的塑性區(qū) ,將其用一系列傾斜的線(xiàn)分成若干楔塊 ,每一楔塊都視為剛體 ,其變形速率為 V。圖 3示出 3個(gè)塊體的系統(tǒng)。 V與滑面夾角為 φ,與右邊相鄰塊體的相對(duì)速度為 Vj ,V j與該兩塊體交界面的夾角為 φj。內(nèi)能耗散發(fā)生于該楔塊的底面和楔塊間的界面 ,在剛體內(nèi)為零。

根據(jù)位移協(xié)調(diào)要求 ,可以得到

1 (4)

1 (5)

其中Vl 和 Vr 分別為左側(cè)和右側(cè)條塊的速度 ;θj =π/2-δ+ φj, θl=π+αl-φl(shuí),θr =π+αr -φr;α為底面與 x軸正向夾角 ;δ為側(cè)面與 y軸正向夾角 ;θ為速度與正 x軸的夾角。如果將條塊的寬度取為無(wú)限小 (圖 4),還可通過(guò)積分獲得滑面上坐標(biāo)為 x的條塊的絕對(duì)速度和相對(duì)速度。

V0為左端點(diǎn)（x=x0）的速度。在滑裂面上第k個(gè)α或φ 發(fā)生突變。上標(biāo)l和代表該不連續(xù)點(diǎn)左和右的物理量。計(jì)算從第一個(gè)界面開(kāi)始，到第 n-1個(gè)界面終止。這樣，對(duì)滑面上橫坐標(biāo)為x的任意一點(diǎn)，其條塊絕對(duì)速度V和條塊側(cè)面的相對(duì)速度Vj 都可表達(dá)為滑面左端點(diǎn)x=x0處的速度V0的函數(shù)。將獲得各條塊的絕對(duì)速度和相對(duì)速度代入式 (3)再代入式 (1),其中式 (1)左側(cè)第一項(xiàng)可通過(guò)將Vj替代式（3）中 V獲得。消去左右側(cè) V0 ,就可求解加載系數(shù) η*如下。

定義

其中d W =土條重量 ;T0x, T0y分別為 T0在 x和 y軸的分量 ;L為土條側(cè)面長(zhǎng)度 ;η′為水平地震力系數(shù)。式 (10)最后一項(xiàng)計(jì)及滑面上 (n -1)個(gè) α或φ的不連續(xù)點(diǎn)相應(yīng)的界面上的內(nèi)能耗散。由式 (2)定義的加載系數(shù)可通過(guò)下式計(jì)算 :

η* = G/Gb(12)

2.4求解臨界滑動(dòng)模式

陳祖煜和邵長(zhǎng)明曾詳細(xì)介紹了對(duì)傳統(tǒng)的極限平衡法計(jì)算最小安全系數(shù)和臨界滑裂面的數(shù)值分析方法。最優(yōu)化方法為使用計(jì)算機(jī)搜索臨界滑動(dòng)模式創(chuàng)造了條件 ,這些研究成果可以方便地推廣到本文介紹的極限分析方法中。所不同的是 ,滑動(dòng)模式和垂直條分法相比 ,增加了一個(gè)土條界面傾角 δ。每一條塊的 δ也將成為自由度。最優(yōu)化方法將最終找到相應(yīng)最小加載系數(shù)的滑裂面和斜分條模式。具體計(jì)算步驟通過(guò)下節(jié) [例 2 ]介紹。

3驗(yàn)證

為了驗(yàn)證上述推導(dǎo)的正確性 ,下面通過(guò)兩個(gè)例題進(jìn)行分析探討。

[例 1 ]對(duì)具有垂直表面荷載的例題 (圖 5),索科洛夫斯基 ( Sokolovski , 1954)給出的臨界垂直荷載 q的計(jì)算公式如下 :

的 q=c ctgφexp-2x)tgφ-1

其中 χ為邊坡斜面相對(duì)水平線(xiàn)的夾角。相應(yīng)的臨界滑裂面由三段組成 ,AB , CD為直線(xiàn) ,分別與邊坡線(xiàn)和坡頂線(xiàn)夾角為 μ。

μ=

BC為一對(duì)數(shù)螺旋線(xiàn) ,其左右邊界線(xiàn) BO和 CO分別與邊坡線(xiàn)和 y軸線(xiàn)夾角為 μ。

當(dāng)邊坡處于極限狀態(tài)時(shí) ,加載系數(shù) η=0。對(duì) AB ,BC和 CD段分別進(jìn)行積分 ,按式 (13)確

定的 q將使按式 (10)確定的 G為零。

這一實(shí)例說(shuō)明 ,本文提出的上限定理的命題可通過(guò)解析解獲得印證。

[例 2 ]某一坡角為 35°的均質(zhì)邊坡 ,其水平頂面上作用一均布荷載 ,荷載方向相對(duì)鉛直

線(xiàn)夾角為 δ′(圖 6)。根據(jù)索科洛夫斯基 (1954)提出的滑移線(xiàn)方法 ,此題的理論破壞面由直線(xiàn) AB , CD和對(duì)數(shù)螺旋線(xiàn) BC組成 , CD和 CO分別相對(duì)鉛直線(xiàn)夾角 μ+ρ和 -(μ+ρ)。其中 μ = 45°-φ/ 2 ;ρ為大主應(yīng)力相對(duì)鉛直線(xiàn)的夾角。

主要參數(shù) :c = 720kPa ,φ = 37°,χ = 35°,δ′= 24°,理論解提供的解答是 q =6 . 228 MPa , ρ= 28. 4°。理論的滑裂面和土條側(cè)面示于圖 6線(xiàn) 4?；衙嫱ㄟ^(guò)聯(lián)結(jié) 4個(gè)點(diǎn)的樣條函數(shù)形成。對(duì)設(shè)定的初始滑裂面 1和相應(yīng)的斜分條模式使用式 (12)求得 η3 =0127。從這個(gè)滑裂面開(kāi)始 ,進(jìn)行最優(yōu)化方法計(jì)算最終得臨界滑裂面和條間界面 (滑裂面 3 ,虛線(xiàn) ),相應(yīng) η=01019?；衙?2是優(yōu)化計(jì)算過(guò)程中通過(guò)隨機(jī)搜索獲得的滑裂面。如果用 5個(gè)點(diǎn)來(lái)模擬該滑裂面 ,則可得到 η=010028。與理論解相比 ,無(wú)論是最小加載系數(shù) ,還是臨界滑裂面和臨界條間界面均十分接近。

通過(guò) [例 2 ]說(shuō)明 ,應(yīng)用最優(yōu)化方法可以自動(dòng)找到相應(yīng)最小加載系數(shù) η的臨界滑裂面和相應(yīng)的斜分條模式。

4能量法在地基極限承載力計(jì)算中的推廣

4．1傳統(tǒng)的承載力計(jì)算方法

地基極限承載力的計(jì)算包括兩個(gè)方面，一方面是允許位移的校核 ,另一方面是極限承載力的計(jì)算。對(duì)于后者， Prandtl于 1920年根據(jù)塑性力學(xué)理論導(dǎo)出了剛性基礎(chǔ)壓入無(wú)重量土中滑裂面的形狀及其相應(yīng)的極限承載力計(jì)算公式。由于數(shù)學(xué)上的嚴(yán)格解答在大部分的實(shí)際問(wèn)題中是不可能得到的 , Terzaghi ,Meyerhof ,Vesic等眾多學(xué)者在 Prandtl解的基礎(chǔ)上對(duì)承載力理論進(jìn)行了研究和發(fā)展 ,最終形成地基極限承載力的近似解答。這一近似解答的一般表達(dá)式為

Nγ為一半經(jīng)驗(yàn)數(shù)據(jù) ,可從地基規(guī)范承載力表中查取或用半經(jīng)驗(yàn)公式 (表 1)計(jì)算 ;B為基礎(chǔ)寬度 ;D為基礎(chǔ)埋深 ;γ為土容量 ; qu為地基的極限承載力 ,即 T在單位寬度上的強(qiáng)度。

4.2能量法在地基極限承載力計(jì)算中的應(yīng)用

選取寬度 B = 17m的條形基礎(chǔ)進(jìn)行計(jì)算分析 ,相應(yīng)參數(shù)為 :c = 144. 5kPa ,γ=0. 0kN/ m3。這個(gè)例子針對(duì)土的不同內(nèi)摩擦角 φ值進(jìn)行計(jì)算。對(duì)于具有理論解的實(shí)例 ,使用式 (13)獲得的 qu應(yīng)保證使用式 (12)獲得的 η的最小值為零。圖 7示出 φ =0°和φ = 20°兩種情況。使用同樣的初始破壞模式如圖 7 (a),應(yīng)用最優(yōu)化方法獲得的臨界破壞模式分別如圖 7 ( b)和 7 (c)所示 ,η分別為 01004和 01008。計(jì)算機(jī)在搜索最小 η值時(shí) ,準(zhǔn)確地將中部的土條側(cè)面收斂于地基的左側(cè)點(diǎn) ,由此將滑裂面分為 3個(gè)區(qū)域 :荷載作用下面的三角形區(qū)域 ,對(duì)應(yīng)于理論上的主動(dòng) Rankine區(qū) ;條間界面一端收斂于一點(diǎn)所對(duì)應(yīng)的放射形區(qū)域 ,對(duì)應(yīng)于理論上的 Prandtl區(qū);放射形區(qū)域另一端的另一個(gè)三角形區(qū)域?qū)?yīng)理論上的被動(dòng) Rankine區(qū)。 φ=0°時(shí) ,滑裂面形狀接近于圓弧 ;φ≠0°時(shí) ,滑裂面形狀為兩段直線(xiàn)接一段對(duì)數(shù)螺旋線(xiàn)。這就說(shuō)明 ,采用建立在塑性力學(xué)上限 解基礎(chǔ)上的地基極限承載力數(shù)值分析方法直接獲得了理論上嚴(yán)格的地基極限承載力解答。

表 2將一系列φ值的計(jì)算成果與理論解對(duì)比，可見(jiàn)成果的準(zhǔn)確性相當(dāng)穩(wěn)定。所得的η值與理論值的誤差均在1%左右，而且自動(dòng)搜索得到的臨界滑裂面形狀也與理論解一致。

表2 上限數(shù)值解和理論解成果對(duì)比

注：qn和 qp分別為根據(jù)數(shù)值解和理論解獲得的極限承載力 ;η= ( qn/ qp -1);q的單位為 kN/ m。

表 3比較了容重不為零情況下表 1所列的各種經(jīng)驗(yàn)方法的準(zhǔn)確性 ,并示于圖 8。可見(jiàn) W1F Chen的公式計(jì)算結(jié)果與采用 Prandtl的經(jīng)驗(yàn)公式求得的結(jié)果最為接近 ,但是在 φ值超過(guò) 25°后 ,η值為負(fù) ,意味著高估了地基的承載力。而 Mayerhof以及 Terzaghi的方法則偏于保守。

表 3對(duì) γ≠0情況各種不同的經(jīng)驗(yàn)公式和數(shù)值解成果對(duì)比

注 :下標(biāo) n,v, m,t,c分別代表數(shù)值解、采用VesicMeyerhofTerzaghi 和 W. F. Chen方法的計(jì)算成果，q的單位為 kN/ m。，

在有容重且有埋深的條件下 ,式 (15)的經(jīng)驗(yàn)公式將基礎(chǔ)兩側(cè)埋置深度以?xún)?nèi)的土重以連續(xù)均勻分布的荷載考慮 ,未能計(jì)及這部分土體的抗剪作用 ,因此不可避免地存在著一定的誤差。