薄壁彎管在內(nèi)壓和彎矩作用下的彈塑性坍塌分析

案例 Elastic-plastic collapse of a thin-walled elbow under in-plane bending and internal pressure

問題描述與目標(biāo)

本案例旨在研究一個薄壁90度不銹鋼彎管及其相鄰直管段，在面內(nèi)彎矩（張開和閉合彎矩）與內(nèi)部壓力共同作用下的彈塑性響應(yīng)直至結(jié)構(gòu)坍塌。通過此分析，可以深入理解彎管的橢圓化變形、幾何非線性效應(yīng)、內(nèi)部壓力的影響以及彎管與直管段之間的相互作用。

它是一個薄壁彎頭，彎頭系數(shù)，半徑比 R/r = 3.07

建模詳細(xì)步驟

步驟 1：模型規(guī)劃與假設(shè)

分析類型： 靜態(tài)、通用分析步（Static），考慮幾何非線性（NLGEOM）。

對稱性利用： 由于模型和載荷的對稱性，可以只建立四分之一的模型以顯著減少計算量。

單元選擇： 主體使用減縮積分殼單元（如S8R5或S4R）來模擬薄壁管道結(jié)構(gòu)。

步驟 2：創(chuàng)建部件

分別創(chuàng)建代表直管段和90度彎管段的殼體部件。或一體的彎管部件。

幾何尺寸需嚴(yán)格按照案例提供的圖紙進(jìn)行。關(guān)鍵尺寸包括：管道外徑、壁厚、彎管中心線半徑等。

步驟 3：定義材料屬性

材料模型： 定義彈塑性材料。楊氏模量 E = 193MPa，泊松比u= 0.264。名稱為PIPE。

塑性數(shù)據(jù)：Mechanical--plastic 輸入真實應(yīng)力-塑性應(yīng)變曲線數(shù)據(jù)如下。該曲線顯示材料具有明顯的屈服平臺和硬化行為。

截面屬性Section-LEG,直管段：壁厚0.37

截面屬性Section-LBOW，彎管段：壁厚0.41

步驟 4：裝配與定位

1.將直管彎管部件實例化裝配，直管和彎管間tie（若彎管用ELBOW31B單元，直管用B31單元）。---對應(yīng)elbowcollapse_elbow31b_b31.inp

2.若用S8R5單元則只有一個零件，不需裝配。

---對應(yīng)elbowcollapse_s8r5.inp和elbowcollapse_s8r5_fine.inp

步驟 5：設(shè)置分析步與輸出請求

創(chuàng)建分析步： 創(chuàng)建static Riks分析步，用于施加彎矩和壓力載荷。

控制參數(shù)：打開大變形NLgeom:On。由于涉及坍塌（極值點失穩(wěn)），通常需要使用弧長法（Riks） 或設(shè)置非常小的初始增量步0.05來控制求解過程。

場輸出請求： 確保輸出應(yīng)力（S）、應(yīng)變（E）、位移（U）等。

增加輸出請求： 輸出Nout點集合的施加彎矩一端的反作用力矩（RM）和轉(zhuǎn)角（UR），用于繪制力矩-轉(zhuǎn)角曲線、橢圓變形等。

步驟 6：定義相互作用

綁定約束： 使用“Tie”約束將彎管段與直管段的端面完全連接在一起。若用S8R5單元則只有一個零件，不需要。

對稱邊界條件： 在對稱面上定義對稱邊界條件（XSYMM）。一端全約束。

步驟 7：施加載荷與邊界條件

固定端： 約束遠(yuǎn)離彎管的直管段末端的全部自由度（ENCASTRE）。

加載端： 在另一個直管段的末端，創(chuàng)建一個參考點（RP），并將該端面的所有節(jié)點與RP進(jìn)行運動耦合約束（Kinematic Coupling），以模擬剛性端蓋。

載荷施加：

內(nèi)部壓力： 在分析步中，作為表面壓力載荷施加在所有管道的內(nèi)表面上，3.45 MPa。

面內(nèi)彎矩： 通過在加載端的參考點（RP）上施加轉(zhuǎn)角（Rotation） 來間接實現(xiàn)彎矩加載。分別施加正向和反向的轉(zhuǎn)角來模擬“張開”和“閉合”彎矩。

步驟 8：劃分網(wǎng)格

對殼體部件進(jìn)行網(wǎng)格劃分。為確保精度，特別是在彎管曲率較大的區(qū)域，需要進(jìn)行適當(dāng)?shù)木W(wǎng)格細(xì)化。

案例中使用了S8R5（8節(jié)點四邊形殼）單元，并進(jìn)行了網(wǎng)格收斂性研究。

步驟 9：提交作業(yè)與計算

創(chuàng)建作業(yè)并提交計算。監(jiān)控求解過程（.sta文件），注意是否有收斂困難。

后處理與結(jié)果分析

步驟 10：驗證與結(jié)果提取

力矩-轉(zhuǎn)角曲線： 繪制加載端參考點的反作用力矩（RM）與轉(zhuǎn)角（UR）的關(guān)系曲線。這是評估結(jié)構(gòu)剛度和預(yù)測坍塌彎矩的關(guān)鍵結(jié)果。

結(jié)果對比： 將不同內(nèi)部壓力下的張開、閉合彎矩曲線進(jìn)行對比，觀察壓力對承載能力和坍塌行為的影響。

應(yīng)力應(yīng)變云圖： 查看在坍塌時刻，管道外表面的環(huán)向應(yīng)變和軸向應(yīng)變分布，觀察橢圓化變形的模式和應(yīng)變集中區(qū)域。

坍塌彎矩預(yù)測： 從力矩-轉(zhuǎn)角曲線的峰值點確定坍塌彎矩，并與理論解或文獻(xiàn)中的實驗結(jié)果進(jìn)行對比驗證。

關(guān)鍵要點總結(jié)

幾何非線性至關(guān)重要： 即使轉(zhuǎn)角不大，由于彎管截面橢圓化與整體彎曲的耦合，也必須打開NLGEOM選項。

彎管-直管相互作用： 直管段會約束彎管的自由橢圓化，從而影響其柔性和坍塌行為，建模時必須包含足夠長的直管段。

內(nèi)部壓力的雙重作用： 壓力會產(chǎn)生環(huán)向拉應(yīng)力，一方面可能延緩塑性屈服，另一方面又可能促進(jìn)屈曲，需要綜合分析。

材料屬性的敏感性： 分析結(jié)果對輸入的材料應(yīng)力-應(yīng)變曲線非常敏感，微小的差異可能導(dǎo)致預(yù)測響應(yīng)的顯著偏差。

S8R5單元ELBOW31B和B31的對比

追求精度，驗證機理：應(yīng)使用S8R5或其他殼單元建立詳細(xì)的有限元模型。這是研究彎管復(fù)雜非線性行為（如彈塑性坍塌）最可靠的方法。

工程分析，平衡效率：對于包含彎頭的管道系統(tǒng)進(jìn)行系統(tǒng)級分析時，推薦采用 混合建模：彎頭部分使用 ELBOW31B。它能以遠(yuǎn)低于殼單元的成本，提供比標(biāo)準(zhǔn)梁單元合理得多的預(yù)測結(jié)果。直管段部分：使用 B31。在此區(qū)域，它的假設(shè)是合理的，且能最大程度節(jié)省計算資源。

重要提示：案例分析表明，直管段（B31）與彎管（ELBOW31B）之間的相互作用對坍塌載荷有顯著影響，不可忽略。因此，即使使用簡化單元，也應(yīng)建立包含相鄰直管段的完整模型。

初步設(shè)計，快速估算：在對精度要求不高的初期，或進(jìn)行大量工況篩查時，可全部使用 B31 單元進(jìn)行極度簡化的分析，但必須清楚其結(jié)果可能是保守或不保守的，不能用于最終設(shè)計。