第四章彈塑性有限元法基本理論與模擬方法演示文稿

第四章彈塑性有限元法基本理論與模擬方法演示文稿目前一頁\總數(shù)一百零五頁\編于十八點(diǎn)優(yōu)選第四章彈塑性有限元法基本理論與模擬方法目前二頁\總數(shù)一百零五頁\編于十八點(diǎn)非線性問題可以分為三類：材料非線性：體系的非線性由材料的應(yīng)力應(yīng)變關(guān)系的非線性引起。如金屬變形彈塑性行為、橡膠的超彈性行為等幾何非線性：結(jié)構(gòu)的位移使體系的受力狀態(tài)發(fā)生了顯著的變化。如板殼的大撓度問題——平衡方程必須建立于變形后的狀態(tài)接觸非線性：接觸狀態(tài)的變化所引起。如金屬成形、跌落試驗(yàn)、多零件裝配體等目前三頁\總數(shù)一百零五頁\編于十八點(diǎn)碰到障礙物的懸臂梁（端部碰到障礙物時(shí)，梁端部的邊界條件發(fā)生了突然變化，阻止了進(jìn)一步的豎向撓度。）板料的沖壓成形接觸非線性例子目前四頁\總數(shù)一百零五頁\編于十八點(diǎn)隨著有限元算法理論、計(jì)算機(jī)硬件和軟件技術(shù)的進(jìn)步及實(shí)際工業(yè)的需求，CAE技術(shù)的應(yīng)用逐步由線性模擬為主向非線性模擬為主快速發(fā)展。1969年，第一個(gè)商業(yè)非線性有限元程序——Marc誕生。目前幾乎所有的商業(yè)有限元軟件都具備較強(qiáng)的非線性問題的分析求解能力。非線性求解技術(shù)的先進(jìn)性與穩(wěn)健性已經(jīng)成為衡量一個(gè)結(jié)構(gòu)分析程序優(yōu)劣的標(biāo)準(zhǔn)。目前五頁\總數(shù)一百零五頁\編于十八點(diǎn)非線性問題的有限元求解方法非線性方程(組)的求解方法直接迭代法Newton-Raphson迭代法修正的Newton-Raphson迭代法非線性問題通常采用增量法求解（追蹤加載過程中應(yīng)力和變形的演變歷史。）每個(gè)增量步采用Newton-Raphson迭代法非線性問題有限元控制方程：目前六頁\總數(shù)一百零五頁\編于十八點(diǎn)非線性方程的迭代求解方法直接迭代法Newton-Raphson迭代修正的N-R迭代目前七頁\總數(shù)一百零五頁\編于十八點(diǎn)非線性方程組的迭代求解方法目前八頁\總數(shù)一百零五頁\編于十八點(diǎn)直接迭代法N-R迭代修正的N-R迭代第四章彈塑性有限元法基本理論與模擬方法目前九頁\總數(shù)一百零五頁\編于十八點(diǎn)非線性問題的增量法求解過程(1)將總的外力載荷分為一系列載荷段(2)在每一載荷段中進(jìn)行迭代，直至收斂(3)所有載荷段循環(huán)，并將結(jié)果進(jìn)行累加目前十頁\總數(shù)一百零五頁\編于十八點(diǎn)(1)將總的外力載荷分為一系列載荷段目前十一頁\總數(shù)一百零五頁\編于十八點(diǎn)(2)在每一載荷段中進(jìn)行迭代，直至收斂N-R迭代:目前十二頁\總數(shù)一百零五頁\編于十八點(diǎn)(3)所有載荷段循環(huán)，并將結(jié)果進(jìn)行累加目前十三頁\總數(shù)一百零五頁\編于十八點(diǎn)4.2材料非線性問題及分類概念：由于材料的應(yīng)力應(yīng)變非線性關(guān)系引起的非線性。分類：不依賴時(shí)間的彈、塑性問題非線性彈性——橡膠彈塑性——沖壓成形依賴于時(shí)間的粘（彈、塑）性問題蠕變——載荷不變，變形隨時(shí)間繼續(xù)變化松弛——變形不變，應(yīng)力隨時(shí)間衰減目前十四頁\總數(shù)一百零五頁\編于十八點(diǎn)非線性彈性材料行為橡膠應(yīng)力應(yīng)變關(guān)系曲線目前十五頁\總數(shù)一百零五頁\編于十八點(diǎn)彈塑性材料進(jìn)入塑性的特征：載荷卸去后存在不可恢復(fù)的永久變形。應(yīng)力應(yīng)變之間不是單值對應(yīng)關(guān)系，與加載歷史有關(guān)。目前十六頁\總數(shù)一百零五頁\編于十八點(diǎn)單軸應(yīng)力狀態(tài)下彈塑性材料行為單軸（一維）應(yīng)力狀態(tài)下材料的應(yīng)力應(yīng)變行為可以從拉伸試驗(yàn)中獲得。目前十七頁\總數(shù)一百零五頁\編于十八點(diǎn)目前十八頁\總數(shù)一百零五頁\編于十八點(diǎn)單調(diào)加載硬化塑性理想彈塑性目前十九頁\總數(shù)一百零五頁\編于十八點(diǎn)各向同性硬化:運(yùn)動(dòng)硬化:混合硬化:反向加載運(yùn)動(dòng)硬化各向同性硬化混合硬化目前二十頁\總數(shù)一百零五頁\編于十八點(diǎn)

在簡單拉伸的情況下，當(dāng)材料發(fā)生塑性變形后卸載，此后再重新加載，則應(yīng)力和應(yīng)變的變化仍服從彈性關(guān)系，直至應(yīng)力到達(dá)卸載前曾經(jīng)達(dá)到的最高應(yīng)力點(diǎn)時(shí)，材料才再次屈服（后繼屈服）。這個(gè)最高應(yīng)力點(diǎn)的應(yīng)力就是材料在經(jīng)歷了塑性變形后的新的屈服應(yīng)力。由于材料的強(qiáng)化特性，它比初始屈服應(yīng)力大。第四章彈塑性有限元法基本理論與模擬方法目前二十一頁\總數(shù)一百零五頁\編于十八點(diǎn)為了與初始屈服應(yīng)力相區(qū)別，我們稱之為后繼屈服應(yīng)力。與初始屈服應(yīng)力不同，它不是一個(gè)材料常數(shù)，而是依賴于塑性變形的大小和歷史。后繼屈服應(yīng)力是在簡單拉伸下，材料在經(jīng)歷一定塑性變形后再次加載時(shí)，變形是按彈性還是塑性規(guī)律變化的界限。第四章彈塑性有限元法基本理論與模擬方法目前二十二頁\總數(shù)一百零五頁\編于十八點(diǎn)

和簡單應(yīng)力狀態(tài)相似，材料在復(fù)雜應(yīng)力狀態(tài)下同樣存在初始屈服和后繼屈服的問題。

材料在復(fù)雜應(yīng)力狀態(tài)下，在經(jīng)歷初始屈服和發(fā)生塑性變形后，此時(shí)卸載，將再次進(jìn)入彈性狀態(tài)（稱為后繼彈性狀態(tài)）。第四章彈塑性有限元法基本理論與模擬方法目前二十三頁\總數(shù)一百零五頁\編于十八點(diǎn)

把復(fù)雜應(yīng)力狀態(tài)下，確定材料后繼彈性狀態(tài)的界限的準(zhǔn)則就稱為后繼屈服條件，又稱為加載條件。問題:

當(dāng)材料處于后繼彈性狀態(tài)而繼續(xù)加載時(shí)，應(yīng)力（或變形）發(fā)展到什么程度材料再一次開始屈服呢？第四章彈塑性有限元法基本理論與模擬方法目前二十四頁\總數(shù)一百零五頁\編于十八點(diǎn)一般應(yīng)力狀態(tài)下彈塑性材料行為屈服準(zhǔn)則（初始屈服條件）硬化法則（后繼屈服函數(shù)、加載函數(shù)、加載曲面）流動(dòng)法則加載、卸載準(zhǔn)則目前二十五頁\總數(shù)一百零五頁\編于十八點(diǎn)屈服準(zhǔn)則（初始屈服條件）在單向受力情況下，當(dāng)應(yīng)力達(dá)到材料的屈服強(qiáng)度時(shí)材料開始產(chǎn)生塑性變形。對于一般復(fù)雜的應(yīng)力狀態(tài)，應(yīng)力狀態(tài)由六個(gè)應(yīng)力分量決定時(shí)，顯然不能根據(jù)某個(gè)單獨(dú)應(yīng)力分量的數(shù)值作為判斷材料是否進(jìn)入塑性變形的標(biāo)準(zhǔn)。為此，引入以應(yīng)力分量為坐標(biāo)的應(yīng)力空間，根據(jù)代表不同應(yīng)力路徑的實(shí)驗(yàn)結(jié)果，可以定出從彈性階段進(jìn)入塑性階段的各個(gè)界限，即屈服應(yīng)力點(diǎn)。在應(yīng)力空間中，這些屈服應(yīng)力點(diǎn)形成一個(gè)區(qū)分彈性和塑性的分界面——屈服面。描述這個(gè)屈服面的數(shù)學(xué)表達(dá)式就是我們所要尋求的一般應(yīng)力狀態(tài)下的屈服準(zhǔn)則。目前二十六頁\總數(shù)一百零五頁\編于十八點(diǎn)常用的各向同性Von-Mises屈服準(zhǔn)則：各向同性屈服準(zhǔn)則：各個(gè)方向屈服應(yīng)力相同各向異性屈服準(zhǔn)則：不同方向屈服應(yīng)力有差異目前二十七頁\總數(shù)一百零五頁\編于十八點(diǎn)三維主應(yīng)力空間π平面上的屈服軌跡σ3=0平面上的屈服軌跡第四章彈塑性有限元法基本理論與模擬方法目前二十八頁\總數(shù)一百零五頁\編于十八點(diǎn)硬化法則塑性硬化法則規(guī)定了材料進(jìn)入塑性變形后的后繼屈服函數(shù)（又稱加載函數(shù)或加載曲面）各向同性硬化運(yùn)動(dòng)硬化混合硬化目前二十九頁\總數(shù)一百零五頁\編于十八點(diǎn)運(yùn)動(dòng)硬化：該模型假設(shè)材料隨塑性變形發(fā)展時(shí)，屈服面的大小和形狀不變，僅是整體在應(yīng)力空間作平動(dòng)。

各向同性硬化：材料進(jìn)入塑性變形以后，屈服面在各方向均勻地向外擴(kuò)張，其形狀、中心及其在應(yīng)力空間的方位均保持不變。

材料的強(qiáng)化只與總的塑性變形功有關(guān)而與加載路徑無關(guān)。

應(yīng)力有反復(fù)變化時(shí)，等向強(qiáng)化模型與實(shí)驗(yàn)結(jié)果不相符合。

目前三十頁\總數(shù)一百零五頁\編于十八點(diǎn)混合硬化：其實(shí)質(zhì)就是將隨動(dòng)強(qiáng)化模型和等向強(qiáng)化模型結(jié)合起來，即認(rèn)為后繼屈服面的形狀、大小和位置一起隨塑性變形的發(fā)展而變化。該模型能夠更好的反映材料的Bauschinger效應(yīng)。各向同性硬化運(yùn)動(dòng)硬化第四章彈塑性有限元法基本理論與模擬方法目前三十一頁\總數(shù)一百零五頁\編于十八點(diǎn)流動(dòng)法則塑性應(yīng)變增量和應(yīng)力分量的關(guān)系：塑性應(yīng)變沿后繼屈服面F=0的法線方向——是一正的待定系數(shù)，其具體數(shù)值和材料硬化準(zhǔn)則有關(guān)目前三十二頁\總數(shù)一百零五頁\編于十八點(diǎn)加載、卸載準(zhǔn)則對于硬化材料（當(dāng)材料處于某一塑性狀態(tài)）：目前三十三頁\總數(shù)一百零五頁\編于十八點(diǎn)4.3幾何非線性問題及分類概念：由于大位移、大轉(zhuǎn)動(dòng)而引起的非線性。分類：大位移、大轉(zhuǎn)動(dòng)、小應(yīng)變問題

——板殼的大撓度和后屈曲大位移、大轉(zhuǎn)動(dòng)、大應(yīng)變問題

——薄板成形、彈性材料的受力目前三十四頁\總數(shù)一百零五頁\編于十八點(diǎn)比較：線彈性—幾何非線性線彈性：小變形假設(shè)——假定物體發(fā)生的位移遠(yuǎn)小于物體本身的幾何尺寸，應(yīng)變遠(yuǎn)小于1。建立平衡方程時(shí)不考慮物體位置和形狀的變化。幾何非線性：物體發(fā)生有限變形——大位移、大轉(zhuǎn)動(dòng)的情況。建立平衡方程時(shí)必須考慮物體位置和形狀的變化。目前三十五頁\總數(shù)一百零五頁\編于十八點(diǎn)4.4彈塑性矩陣

應(yīng)力與應(yīng)變的關(guān)系有各種不同的近似表達(dá)式和簡化式。根據(jù)普蘭特爾—羅伊斯（Prandtl-Reuss）假設(shè)和密賽斯屈服準(zhǔn)則，當(dāng)外作用力較小時(shí)，變形體內(nèi)的等效應(yīng)力小于屈服極限時(shí)為彈性狀態(tài)。當(dāng)外力增大到某一值，等效應(yīng)力達(dá)到屈服應(yīng)力，材料進(jìn)入塑性狀態(tài)，這時(shí)變形包括彈性變形和塑性變形兩部分，即：

式中下腳e、p分別表示彈、塑性狀態(tài)。

(4-10)在彈性階段，應(yīng)力與應(yīng)變關(guān)系符合虎克定律。進(jìn)入塑性狀態(tài)后，符合Prandtl-Reuss假設(shè)。目前三十六頁\總數(shù)一百零五頁\編于十八點(diǎn)4.4.1彈性階段

在彈性階段，應(yīng)力和應(yīng)變的關(guān)系是線性的，應(yīng)變僅取決于最后的應(yīng)力狀態(tài)，并且一一對應(yīng)，而與變形過程無關(guān)，有下列全量形式：

式中為彈性矩陣。

(4-11)目前三十七頁\總數(shù)一百零五頁\編于十八點(diǎn)對于各向同性材料，由廣義虎克定律可得：

或：(4-12)式中：是材料的彈性模量，是泊松比。目前三十八頁\總數(shù)一百零五頁\編于十八點(diǎn)對于各向同性材料，廣義虎克定律：

是材料的剪切彈性模量式中：是材料的彈性模量，是泊松比。公式（4-12）的具體推導(dǎo)：

目前三十九頁\總數(shù)一百零五頁\編于十八點(diǎn)（2）+（3）有：將其帶入（1）得：將其帶入（1）得：目前四十頁\總數(shù)一百零五頁\編于十八點(diǎn)同理可推得得表達(dá)式，寫成矩陣的形式，就是：目前四十一頁\總數(shù)一百零五頁\編于十八點(diǎn)4.4.2彈塑性階段

當(dāng)材料所受外力達(dá)到一定值時(shí)，等效應(yīng)力達(dá)到屈服極限，應(yīng)力應(yīng)變關(guān)系曲線由彈塑性矩陣決定，現(xiàn)推導(dǎo)彈塑性矩陣。

等效應(yīng)力為：

對應(yīng)力求導(dǎo)得：

式中為應(yīng)力偏量，

(4-13)(4-14)目前四十二頁\總數(shù)一百零五頁\編于十八點(diǎn)公式（4-14）對應(yīng)力求導(dǎo)的具體推導(dǎo)：

目前四十三頁\總數(shù)一百零五頁\編于十八點(diǎn)目前四十四頁\總數(shù)一百零五頁\編于十八點(diǎn)

由普蘭特爾—羅伊斯關(guān)系有：將式4-14代入式4-15得：寫成矩陣形式為：(4-15)(4-16)(4-17)目前四十五頁\總數(shù)一百零五頁\編于十八點(diǎn)公式（4-16）的具體推導(dǎo)：

目前四十六頁\總數(shù)一百零五頁\編于十八點(diǎn)式中:

又因有：

寫成矩陣乘積的形式為：(4-18)(4-19)目前四十七頁\總數(shù)一百零五頁\編于十八點(diǎn)設(shè)為硬化曲線上任一點(diǎn)的斜率，即將式4-20代入式4-19中得：

將式4-11寫成增量形式為：

(4-20)(4-21)

(4-22)

目前四十八頁\總數(shù)一百零五頁\編于十八點(diǎn)再利用式4-10就可得到：

兩邊同乘以后可得:利用式4-17和式4-21，可將上式寫成：(4-23)(4-24)(4-25)目前四十九頁\總數(shù)一百零五頁\編于十八點(diǎn)由此得：將式4-26代入式4-17得：(4-27)

(4-26)目前五十頁\總數(shù)一百零五頁\編于十八點(diǎn)將式4-27代入式4-23得：(4-28)

目前五十一頁\總數(shù)一百零五頁\編于十八點(diǎn)由式4-14有：則:(4-29)(4-30)

目前五十二頁\總數(shù)一百零五頁\編于十八點(diǎn)將式4-12代入式4-30，并注意到：得：因有：

(4-31)

(4-32)

(4-33)

目前五十三頁\總數(shù)一百零五頁\編于十八點(diǎn)因?yàn)椋核裕汗剑?-31）的具體推導(dǎo)：

目前五十四頁\總數(shù)一百零五頁\編于十八點(diǎn)故：注意：目前五十五頁\總數(shù)一百零五頁\編于十八點(diǎn)因?yàn)椋核裕汗剑?-33）的具體推導(dǎo)：

即：那么：目前五十六頁\總數(shù)一百零五頁\編于十八點(diǎn)或者：目前五十七頁\總數(shù)一百零五頁\編于十八點(diǎn)故令：式4-28可寫成：(4-34)

(4-35)

目前五十八頁\總數(shù)一百零五頁\編于十八點(diǎn)利用上述關(guān)系式可將式4-34表示成顯式，即：(4-36)目前五十九頁\總數(shù)一百零五頁\編于十八點(diǎn)(4-31)公式（4-36）的具體推導(dǎo)：

目前六十頁\總數(shù)一百零五頁\編于十八點(diǎn)目前六十一頁\總數(shù)一百零五頁\編于十八點(diǎn)對于平面應(yīng)力狀態(tài)，，則有：(4-38)(4-37)

目前六十二頁\總數(shù)一百零五頁\編于十八點(diǎn)公式（4-37）的具體推導(dǎo)：

對于平面應(yīng)力狀態(tài)，，則有：(1)(2)得：目前六十三頁\總數(shù)一百零五頁\編于十八點(diǎn)同理可推得的表達(dá)式:寫成矩陣的形式，就是：目前六十四頁\總數(shù)一百零五頁\編于十八點(diǎn)按照上述同樣方法可得：式中：在塑性區(qū)：(4-39)目前六十五頁\總數(shù)一百零五頁\編于十八點(diǎn)公式（4-39）的具體推導(dǎo)：

所以：

目前六十六頁\總數(shù)一百零五頁\編于十八點(diǎn)

由普蘭特爾—羅伊斯關(guān)系有：寫成矩陣形式為：目前六十七頁\總數(shù)一百零五頁\編于十八點(diǎn)仿照前面，不難推得：而：目前六十八頁\總數(shù)一百零五頁\編于十八點(diǎn)目前六十九頁\總數(shù)一百零五頁\編于十八點(diǎn)目前七十頁\總數(shù)一百零五頁\編于十八點(diǎn)其中：目前七十一頁\總數(shù)一百零五頁\編于十八點(diǎn)對于平面應(yīng)變問題，有

只需從式4-12和式4-36中消去上述為零的分量，就可得到下列各式。(4-40)目前七十二頁\總數(shù)一百零五頁\編于十八點(diǎn)(4-42)(4-41)目前七十三頁\總數(shù)一百零五頁\編于十八點(diǎn)4.5變剛度法

變剛度法又稱切線剛度法，它所采用的應(yīng)力與應(yīng)變關(guān)系見圖3-1。在等效應(yīng)力達(dá)到屈服極限后，應(yīng)力與應(yīng)變不再是線性關(guān)系，而是由下列關(guān)系式所確定。(4-43)目前七十四頁\總數(shù)一百零五頁\編于十八點(diǎn)

彈塑性矩陣[D]ep中含有應(yīng)力，它是加載過程的函數(shù)。直接求解是困難的，通常采用增量形式來近似代替微分形式，這樣使求解成為可能。計(jì)算中由于[D]ep在?{σ}范圍內(nèi)變化不大，因此可假設(shè)在每一加載步中是一個(gè)常數(shù)，并以該加載步前的應(yīng)力狀態(tài)近似計(jì)算出[D]ep，即：

單元?jiǎng)偠染仃嘯k]在一個(gè)加載步中也同樣取作常數(shù)，即：

(4-44)(4-45)目前七十五頁\總數(shù)一百零五頁\編于十八點(diǎn)

在一個(gè)變形體中，不僅各點(diǎn)的應(yīng)力狀態(tài)是不相同的，而且隨著加載而變化著，通常變形體受外力作用時(shí)，從一個(gè)區(qū)域到另一個(gè)區(qū)域，等效應(yīng)力是逐漸地達(dá)到屈服極限，即進(jìn)入塑性（彈塑性）狀態(tài)。為了簡化，本章所指進(jìn)入塑性即為彈塑性狀態(tài)，這就是說在變形體中，各單元的應(yīng)力和應(yīng)變狀態(tài)是不一樣的，隨著加載又是變化的，且各有各的變化規(guī)律。變形體內(nèi)的單元按狀態(tài)可分為三類:

彈性單元塑性單元過渡單元各類單元有不同的本構(gòu)關(guān)系和單元?jiǎng)偠染仃?。目前七十六頁\總數(shù)一百零五頁\編于十八點(diǎn)在加載過程中，各單元的狀態(tài)是變化的，為此[K]也是變化的。在計(jì)算中，每增加一個(gè)載荷增量，就得重新計(jì)算一次整體剛度矩陣[K]，這也就是變剛度法名稱的由來。式中[K]整體剛度矩陣；

n1、n2、n3分別為彈性單元、塑性單元和過渡單元的數(shù)量；

[k]e、[k]ep、[k]g

分別為彈性單元、塑性單元和過渡單元的單元?jiǎng)偠染仃嚒τ谡w來說，可用下列關(guān)系式表示：(4-46)目前七十七頁\總數(shù)一百零五頁\編于十八點(diǎn)整體剛度矩陣求得之后，就可根據(jù)下列載荷和位移的線性方程組求解出未知的節(jié)點(diǎn)位移增量。有了節(jié)點(diǎn)位移增量就能求得各單元的應(yīng)變及應(yīng)力增量。(4-47)(4-49)(4-48)目前七十八頁\總數(shù)一百零五頁\編于十八點(diǎn)4.5.1定加載法

定加載法又稱等量加載法。它每次的加載量是預(yù)先給定的。這種加載法的加載量一般較大。由于每次加載量較大，每次加載中由彈性單元轉(zhuǎn)變?yōu)閺椝苄詥卧倪^渡單元較多。過渡單元在加載步中達(dá)到屈服，式中m為加權(quán)系數(shù)，0≤m≤1，采用不同的加載方法，過渡單元的處理也有所不同。下面介紹幾種加載方法。(4-50)目前七十九頁\總數(shù)一百零五頁\編于十八點(diǎn)

m的取值需要進(jìn)行迭代來逼近，收斂性一般都很好，只需進(jìn)行2～3次迭代就能達(dá)到滿意的精確度。目前八十頁\總數(shù)一百零五頁\編于十八點(diǎn)定加載法計(jì)算程序框圖目前八十一頁\總數(shù)一百零五頁\編于十八點(diǎn)4.5.2變加載法這種方法又被稱作r因子法。用這種方法計(jì)算，每次加載量是變化的，其大小是由計(jì)算結(jié)果來確定。計(jì)算開始時(shí)預(yù)先施加一個(gè)單位載荷增量，然后求出各單元在施加單位載荷增量后的等效應(yīng)力增加量。根據(jù)這個(gè)增加量求出各彈性單元當(dāng)達(dá)到屈服時(shí)所需要施加的增量值，最后取這些增量值中最小的一個(gè)增量值作為本次加載的加載量。目前八十二頁\總數(shù)一百零五頁\編于十八點(diǎn)各彈性單元的加載因子按下式進(jìn)行計(jì)算：式中,i單元前次加載后的等效應(yīng)力；

,i單元本次施加單位載荷增量后的等效應(yīng)力；

,i單元達(dá)到屈服所需施加單位加載量的倍數(shù)。(4-51)目前八十三頁\總數(shù)一百零五頁\編于十八點(diǎn)為了加快計(jì)算步伐，常假設(shè)單元的等效應(yīng)力接近屈服極限時(shí)就由彈性單元轉(zhuǎn)變?yōu)閺椝苄詥卧Ｒ话憧扇卧刃?yīng)力，在下一次施加載荷增量的計(jì)算中，這單元就按彈塑性單元處理。采用這種處理方法，能保證每次加載后彈性單元中等效應(yīng)力的最大者正好達(dá)到屈服。在下一次加載中該單元按彈塑性單元處理。這種方法能避免在每個(gè)加載步中單元由彈性轉(zhuǎn)變?yōu)閺椝苄运枰?jì)算m因子的過程，還能保證足夠好的計(jì)算精度。目前八十四頁\總數(shù)一百零五頁\編于十八點(diǎn)變加載法的計(jì)算程序框圖目前八十五頁\總數(shù)一百零五頁\編于十八點(diǎn)4.5.3位移法

對于有些問題如圓柱體鐓粗，每次施加的增量不是控制加載力，而是控制壓下量。假設(shè)工具為剛性體，在工具與坯料的接觸面上，各節(jié)點(diǎn)的位移相同。而接觸面上的壓力分布是未知的。這樣在計(jì)算中需對與接觸表面上節(jié)點(diǎn)有關(guān)的方程進(jìn)行處理。計(jì)算這類問題時(shí)，一般先假設(shè)接觸表面上有一個(gè)已知的軸向位移增量，根據(jù)這個(gè)已知位移增量求出各單元的應(yīng)力和應(yīng)變增量。這種施加位移增量的方法又分為定位移增量法和變位移增量法。目前八十六頁\總數(shù)一百零五頁\編于十八點(diǎn)4.5.3.1定位移增量法定位移增量法每次施加的位移增量相同。計(jì)算過程先算出各單元的應(yīng)力和應(yīng)變增量，有了應(yīng)力增量，累加后可得全應(yīng)力和計(jì)算得到等效應(yīng)力。根據(jù)等效應(yīng)力的大小將單元分成彈性單元、塑性單元和過渡單元。對于過渡單元的處理與定加載法相同，為此這種方法在求取m因子時(shí)也需要進(jìn)行迭代。

目前八十七頁\總數(shù)一百零五頁\編于十八點(diǎn)4.5.3.2變位移增量法

這種方法與變加載法相似，每次施加的位移增量由計(jì)算結(jié)果來決定。每步計(jì)算也是先施加一個(gè)單位位移增量，然后根據(jù)這位移增量計(jì)算應(yīng)力增量，從而找出彈性單元中等效應(yīng)力增長最快而又最先達(dá)到屈服極限的單元，以這個(gè)單元達(dá)到屈服所需要的位移增量為本次施加的位移增量。達(dá)到屈服極限的彈性單元在下次計(jì)算中按彈塑性單元處理。

目前八十八頁\總數(shù)一百零五頁\編于十八點(diǎn)4.6初載荷法

初載荷法是將塑性變形問題試圖轉(zhuǎn)化為彈性問題來求解，它把塑性變形部分視作初應(yīng)力或初應(yīng)變來處理。

在彈性有限元中，當(dāng)彈性體的單元中存在初應(yīng)變{ε}，如因溫度而引起的應(yīng)變，或有初應(yīng)力，如殘余應(yīng)力，則應(yīng)力和應(yīng)變關(guān)系分別為：(4-52)(4-53)目前八十九頁\總數(shù)一百零五頁\編于十八點(diǎn)這樣可得有初應(yīng)變時(shí)的變形能為：即：

(4-54)(4-55)目前九十頁\總數(shù)一百零五頁\編于十八點(diǎn)因：得：上式中最后一項(xiàng)為與節(jié)點(diǎn)位移無關(guān)的變形能，由此可得到與節(jié)點(diǎn)位移有關(guān)的變形能為：式中∑是對所有單元求和。

(4-56)(4-57)目前九十一頁\總數(shù)一百零五頁\編于十八點(diǎn)通過變分可得到初載荷時(shí)的有限元公式為：

同樣可導(dǎo)出有初應(yīng)力時(shí)的變形能為：這樣，基本方程就比無初應(yīng)變或無初應(yīng)力時(shí)多一項(xiàng){R}，{R}是作為載荷存在于基本方程中，稱為載荷向量。它是由于有初應(yīng)變或初應(yīng)力而引起的，下面分別導(dǎo)出{R}的表達(dá)式。(4-58)(4-59)目前九十二頁\總數(shù)一百零五頁\編于十八點(diǎn)4.6.1初應(yīng)力法

在小位移的彈塑性問題中，應(yīng)力應(yīng)變關(guān)系為：

為使問題線性化，對于區(qū)域中已達(dá)到屈服的單元，采用逐次加載法。如每次加載的載荷較小，可將上式的微分形式近似地寫作增量形式，即：

(4-60)(4-61)目前九十三頁\總數(shù)一百零五頁\編于十八點(diǎn)得：

因有：將式:

代入上式，并將右邊表示成矩陣乘積形式，變分后得到：

令：由前面可知：(4-64)(4-63)(4-65)(4-62)目前九十四頁\總數(shù)一百零五頁\編于十八點(diǎn)上式要表示成（*）的形式，必有：

寫成增量形式為：式中：(4-66)(4-68)(4-67)目前九十五頁\總數(shù)一百零五頁\編于十八點(diǎn)由上式可以看出，初載荷{R}不僅與加載前的應(yīng)力有關(guān)，而且與該次加載引起的應(yīng)變增量有關(guān)，即式：的兩邊都含未知數(shù)。因此每此加載時(shí)，必須利用迭代法來求解。具體迭代過程：可先取，則得，求得初載矢量，在載荷下，解方程式：此時(shí)方程可寫成：(4-69)目前九十六頁\總數(shù)一百零五頁\編于十八點(diǎn)解方程：得到后，再照上述方法依次迭代計(jì)算下去。寫作一般迭代式為：

當(dāng)相鄰兩次迭代所得初應(yīng)力相差甚小時(shí)，可認(rèn)為迭代結(jié)束。

求得后，再由求和初載荷矢量，然后進(jìn)行第二次迭代計(jì)算。(4-70)(4-71)目前九十七頁\總數(shù)一百零五頁\編于十八點(diǎn)為了考慮過渡單元，運(yùn)用前面引入的加權(quán)系數(shù)m，這樣初載荷矢量為：當(dāng)單元為彈性狀態(tài)時(shí)，取m=1單元為塑性狀態(tài)時(shí)，取m=0若為過渡單元，則取計(jì)算出的m值。(4-72)目前九十八頁\總數(shù)一百零五頁\編于十八點(diǎn)

計(jì)算步驟如下：

⑴施加全部載荷{P}

于結(jié)構(gòu)，按線彈性計(jì)算；⑵算出各單元的等效應(yīng)力，并取出最大值σmax

。若σmax≤σs

，則材料尚未達(dá)到或剛巧達(dá)到屈服，所計(jì)算結(jié)果就是最終結(jié)果。否則，令β=σs/σmax

，則β{P}是