魏曉，張嵐，尤琳

（山東豪邁機械制造有限公司，山東 青島 266061）

關鍵字：彈塑性分析；棘輪分析；有限元

對于壓力容器受壓元件，通常會受到內壓和沿壁厚疊加的循環(huán)熱應力綜合作用。熱應力通常會在瞬態(tài)工況時產生，例如容器的開車和停工。加熱溫度或冷卻溫度的瞬態(tài)改變會形成沿壁厚方向的溫度梯度。對于設有保溫的容器，熱量的損耗非常小，在經歷溫度瞬變后，溫度沿壁厚會再次變?yōu)榫?。因此，對容器或管道進行彈性分析時，由開車的瞬態(tài)溫度變化產生的應力會在后續(xù)的操作中及停車后的一段時間后消失。

當溫度載荷為循環(huán)工況時，結構可能會發(fā)生棘輪失效。棘輪現象是元件受交變機械應力、溫度應力或兩者同時作用時出現的漸增的塑性變形或應變，而熱應力棘輪是由部分或整體的溫差應力引起的。棘輪現象的變形特點是每次加載循環(huán)的前半周和后半周在結構的不同部位輪流產生塑性變形，最終導致壓力容器因漸增性變形而垮塌[1]。

ASME 規(guī)范給出了兩種分析棘輪現象方法，即彈性分析法和彈塑性分析法，彈性分析是一種近似方法，在大部分情況下偏于安全和保守，對棘輪進行嚴密的評定應采用彈-塑性分析方法[2]。

在分析設計中，往往通過有限元程序對結構進行熱-結構耦合分析對棘輪失效模型進行校核。如果滿足以下任何一個條件，則棘輪準則滿足。如果不滿足下述的棘輪準則，則應修正元件的結構或降低作用的載荷并重新進行分析。①塑性應變?yōu)榱?，即表現為完全彈性行為或彈性安定，該準則最為保守；②彈性核準則，即出現塑性安定；③總體尺寸無永久性改變，即結構無明顯漸增性塑性垮塌[3]。

本文對一典型壓力容器結構，由所作用載荷的施加、移去和再次施加，采用彈塑性分析以防止棘輪失效。分析結果表明彈性核是存在的，這意味著該結構處于安定狀態(tài)，不會發(fā)生塑性棘輪失效。

除此之外，為了進一步展示彈性核消失和塑性變形遞增的現象，將25 mm 的球形封頭替換為50 mm 厚的平蓋結構。同時為了與彈性棘輪分析方法對比，本文依據規(guī)范對模型的關鍵部位進行線性化路徑選取，用彈性分析方法對防止棘輪失效進行評定。結果表明用彈性分析方法與彈塑性分析方法所得結論一致。

該結構包括一圓柱殼和球形封頭，三維模型如圖1 所示。建模厚度采用去除腐蝕余量并考慮成型減薄量的有效厚度，圓柱殼外徑為φ700 mm，壁厚為50 mm，球形封頭的厚度為25 mm。

圖1 三維模型Fig.1 3D model

設計壓力：25 MPa

溫度循環(huán)范圍：22 ～ 300 ℃

1 有限元分析

1.1 模型建立

為了加快求解速度和提高求解精度，根據結構與載荷的對稱性，取1/4 模型來進行有限元分析（該模型為軸對稱模型，亦可用2D 模型分析，為了更加直觀，這里筆者采用了3D 模型）。

1.2 材料定義

采用彈性-理想塑性材料模型，假定材料為某碳鋼，屈服強度設置為235 MPa，選擇VON-MISES 屈服函數和與之相關的流動準則，應力應變曲線如圖2所示，未考慮循環(huán)載荷作用下材料的強化行為。另外，還需定義材料的彈性模量、泊松比、線膨脹系數等參 數。

圖2 材料定義Fig.2 Material definition

1.3 劃分網格

本模型主要采用掃掠網格劃分的方法，單元類型為默認的SOLID186 單元。網格劃分完成后節(jié)點數990 036，單元數為362 538，圖3 為有限元網格模型。

圖3 劃分網格后的模型Fig.3 Model after meshing

1.4 耦合分析

圖4 為耦合分析的示意圖，先對結構進行多載荷步瞬態(tài)溫度分析，得到各時刻的溫度分布，再將溫度分布結果作為載荷條件用于靜態(tài)結構分析，需要注意的是二者的載荷步和時間需要一一對應，且導入溫度載荷時要選擇所有時刻的溫度計算結果。同時在求解設置中打開大變形開關，考慮幾何非線性的影響。

圖4 耦合分析Fig.4 Coupled analysis

對于瞬態(tài)熱分析，內表面溫度循環(huán)歷史如圖5所示，考慮前3.5 個循環(huán)，分14 個載荷步進行定義。首先因為熱介質的加入使設備內壁面升溫，升溫到300 ℃以后保載一段時間，此時溫度沿壁厚的分布逐漸趨于均勻，然后開始卸載所加溫度載荷，在卸載的初始時刻，外表面溫度仍保持不變，內表面溫度瞬間降低到環(huán)境溫度22 ℃，使得內外壁面產生較大的溫度梯度。

圖5 溫度載荷Fig.5 Temperature loads

位移邊界條件如圖6 所示，在兩個對稱面上施加對稱約束，并限制底部端面的軸向位移。在垂直容器內壁面上施加均勻的內壓25 MPa，該內壓會在殼體中產生平均環(huán)向應力，如圖7 所示。

圖6 約束Fig.6 Constraints

圖7 機械載荷（內壓）Fig.7 Mechanical load （internal pressure）

3 結果分析

內外壁溫度隨時間沿壁厚的分布曲線如圖8 所示，圖9 為溫度梯度示意圖，可見沿壁厚方向形成了明顯的溫度梯度。

圖8 溫度分布曲線（內壁與外壁）Fig.8 Curve of the temperature distribution （inner wall and outer wall ）

圖9 溫度梯度Fig.9 The temperature gradient

圖10 所示為最大塑性應變隨時間變化曲線，可見在第一循環(huán)發(fā)生初始的塑性應變后，后續(xù)循環(huán)沒有產生新的塑性應變，三個循環(huán)中未出現塑性變形的累積。圖11 為彈性核示意圖，通過在后處理結果中查看等效塑性應變，并通過Capped Iso-surface 命令僅顯示塑性應變?yōu)榱愕膮^(qū)域，便可獲得彈性核的云圖。可見彈性核在圓柱殼壁厚和球封頭處都是連續(xù)的，塑性應變僅出現在圓柱殼的內壁以及圓柱殼與球封頭連接處的局部區(qū)域。故所分析的結構處于安定狀態(tài)，也就是說在殼體的沿厚度方向的任一截面上都沒有發(fā)生棘輪現象[4]。

圖10 最大塑性應變隨時間變化曲線圖Fig.10 Curve of the maximum plastic strain along with time

圖11 彈性核示意圖Fig.11 The elastic core

圖12 顯示了球形封頭頂點處的最大位移在三個循環(huán)中隨時間變化的曲線，可看出位移增量在第二個循環(huán)后變得穩(wěn)定，沒有發(fā)生明顯漸增性塑性變形。

圖12 三個循環(huán)中最大位移隨時間變化曲線Fig.12 Curve of the maximum deformation of the three cycles along with time

4 與平蓋結構的結果對比

為了直觀地感受棘輪失效現象，筆者將25mm厚的球形封頭更改為50mm 厚的平蓋結構，保持圓柱殼尺寸不變，結構如圖13 所示。施加相同的溫度和壓力載荷進行棘輪分析，所得塑性應變隨時間的變化曲線和彈性核示意圖分別如圖14 和圖15 所示，平蓋中心處最大位移隨時間變化曲線如圖16 所示。可見在平蓋與筒體的連接處及平蓋厚度方向上彈性核消失，同時最大塑性應變隨時間產生累積，且平蓋中心處的最大位移也隨著每個循環(huán)遞增，發(fā)生漸增性塑性垮塌，故可知該結構是不安定的，會發(fā)生棘輪失效。

圖13 平蓋結構Fig.13 Flat cover structure

圖14 最大塑性應變隨時間變化曲線Fig.14 Curve of the maximum plastic strain along with time

圖15 彈性核示意圖Fig.15 Diagram of the elastic core

圖16 三個循環(huán)中最大位移隨時間變化曲線圖Fig.16 Curve of the maximum deformation of the three cycles along with time

5 與彈性棘輪分析方法的對比

為將上述彈塑性棘輪分析的結果與彈性棘輪分析進行對比，根據ASME 的規(guī)定，對兩種結構用彈性分析的方法進行防止棘輪現象的評定，則應滿足ΔSn，k≤SPS，其中一次加二次當量應力范圍ΔSn，k是由越過截面厚度最高值導得的一次薄膜應力加一次彎曲應力加二次應力（PL+Pb+Q）組合得到的當量應力范圍。SPS為一次加二次當量應力范圍的極限，取max[3S， 2Sy]；但當屈服強度與最大拉伸強度之比超過0.7 或S為與時間相關的性能決定時，SPS取3S[2]。

同樣先進行瞬態(tài)溫度分析，再進行靜態(tài)結構分析，時間和載荷步設置依然對應，但彈性棘輪分析方法不再考慮材料的非線性，僅定義材料的彈性模量、泊松比、線膨脹系數等參數。

兩種結構的線性化路徑如圖17 所示。提取各條路徑的一次應力加上二次應力ΔSn，k，并比較ΔSn，k和SPS的大小。計算結果和評定結果見表1，評定結果顯示球封頭結構不會發(fā)生棘輪，而平蓋結構會發(fā)生棘輪失效，這與彈塑性棘輪分析的結果是一致的。

表1 應力評定（棘輪評定）Tab.1 Stress evaluation （ratcheting assessment）

圖17 路徑定義Fig.17 Definition of paths

6 結論

本文從材料定義、求解設置和結果分析幾個方面，對彈塑性棘輪分析進行了詳細的描述。除了本文述及的圓筒、球形封頭、平蓋等結構外，壓力容器還會涉及更多復雜的結構和載荷工況，結構通常是超靜定的。對于類似相關結構，可以參照該方法進行有限元分析，并結合規(guī)范中的彈性核準則，完成棘輪評定。規(guī)范規(guī)定，在作用最少為三個完整的循環(huán)以后，按棘輪準則進行評定。為證實其收斂性，可能需要增加附加的循環(huán)次數。但在大多數情況下，通過三個循環(huán)就能夠確定彈性核是否存在及其分布區(qū)域[5]。

通過本例還可以看出，棘輪評定的三個準則間是既有區(qū)別又有聯系的[6]。零塑性應變準則實際上是彈性核準則的一個特例，即沒有發(fā)生任何塑性應變，整個結構就是一個大的彈性核，該方法對有些材料來說過于保守。在進行棘輪評定時，總體變形準則與彈性核準則是一致的，當基于彈性核準則結果評定為棘輪時，總體變形就會持續(xù)地增大；而當基于彈性核準則結果評定為安定時，總體變形也會很快趨于穩(wěn)定。

本文只針對棘輪失效模式進行評定，分析過程未考慮計算的完整性，工程中應根據實際情況并考慮多種失效模式進行全面的校核。