基于有限元法的空調(diào)換熱管的塑性失穩(wěn)分析

岳建軍

（英格索蘭亞太工程技術(shù)中心，江蘇 太倉215400）

0 引 言

在工程設備和結(jié)構(gòu)中存在著許多受壓構(gòu)件，對這些受壓構(gòu)件的設計校核不僅要進行強度（應力）和剛度（變形）的計算，往往還要進行穩(wěn)定性計算（也稱為屈曲分析）。為了節(jié)省材料和減輕質(zhì)量，寬長薄壁結(jié)構(gòu)和細長桿狀結(jié)構(gòu)在工程設計中得以廣泛應用，在受壓狀態(tài)下，這類結(jié)構(gòu)經(jīng)常因為結(jié)構(gòu)失穩(wěn)而導致失效或破壞，并非因為結(jié)構(gòu)強度不夠或者剛度不足。歷史上也存在不少由于結(jié)構(gòu)失穩(wěn)而導致的重大事故：1907年發(fā)生的加拿大魁北克省勞倫斯河上的魁北克橋的倒塌事故就是由于該橋梁的懸臂下弦在受壓作用失穩(wěn)所致；1987年發(fā)生在美國哈特福德市的一座體育館網(wǎng)架屋頂坍塌事故則是由于鋼屋架結(jié)構(gòu)發(fā)生壓桿屈曲而導致的。結(jié)構(gòu)失穩(wěn)會導致構(gòu)件突然喪失承壓能力而發(fā)生整體或局部垮塌，這種破壞有可能是災難性的。因此，研究結(jié)構(gòu)的穩(wěn)定性問題，對受壓構(gòu)件進行屈曲分析是工程結(jié)構(gòu)設計校核的主要任務之一。

工程中常見的結(jié)構(gòu)失穩(wěn)包括壓桿失穩(wěn)[1]，以及薄板、薄殼的受壓失穩(wěn)[2]。例如在內(nèi)壓作用下的圓柱形薄殼，壁內(nèi)應力為拉應力，這就是一個強度問題，比如商用空調(diào)機組的管殼式換熱器（如圖1）的圓筒形殼體的強度校核。但如果圓柱形薄殼受到了均勻外壓的作用，此時壁內(nèi)應力變成了壓應力，則當外壓達到某一臨界值時，薄殼的圓形平衡就變?yōu)椴环€(wěn)定，這就是外壓作用下的薄殼失穩(wěn)問題，比如管殼式換熱器中的換熱管會受到制冷劑的外壓作用，在設計時候就需要考慮其穩(wěn)定性問題[3]。

圖1 管殼式換熱器

在工程計算時，要評估結(jié)構(gòu)在載荷作用下的失穩(wěn)風險，可以先進行線性屈曲分析（或稱特征值屈曲分析），線性屈曲分析可以預測結(jié)構(gòu)在失穩(wěn)或塌陷之前所能承受的最大荷載（或稱臨界載荷）。基于有限元模型的特征值分析不僅可以得到屈曲載荷，還可以顯示相應的失穩(wěn)模態(tài)，該方法分析簡單，計算速度快。但特征值屈曲分析基于線彈性理論，其往往產(chǎn)生非保守的結(jié)果，計算得到的臨界載荷要高于實際情況。因此，在實際工程結(jié)構(gòu)分析中一般不采用特征值屈曲分析，但特征值屈曲分析可以作為非線性屈曲分析的初步評估。

在實際的工程應用中，結(jié)構(gòu)在外力作用下往往會經(jīng)歷材料的非線性或邊界的非線性過程，而結(jié)構(gòu)本身也可能存在初始幾何缺陷，因此要研究此情況下的結(jié)構(gòu)失穩(wěn)問題，就需要進行非線性屈曲分析。非線性屈曲分析基于非線性靜力學分析，在計算時隨著荷載的逐漸增加來尋找結(jié)構(gòu)失穩(wěn)的荷載水平，直到求解過程開始趨于發(fā)散為止[4]，如圖2所示。計算模型中可以包括結(jié)構(gòu)的初始幾何缺陷、材料的塑性行為、邊界條件的非線性及大撓度響應等特征[5-7]。基于彈塑性材料的非線性屈曲分析還可以模擬結(jié)構(gòu)的塑性失穩(wěn)，即構(gòu)件在載荷作用下在整體上或某一局部的全厚度上發(fā)生過度塑性變形而不能保持原有形狀的一種失穩(wěn)的失效形式。例如受外壓的薄殼（如前文提到的管殼式換熱器中的換熱銅管）在外壓較大時會被壓潰。

圖2 非線性屈曲分析載荷-位移（撓度）曲線

1 基于有限元法的非線性屈曲分析

1.1 算例

本文以某一長度為200 mm、外徑為26 mm、內(nèi)徑為24 mm、壁厚為1 mm的薄壁圓管為例（如圖3），假設的邊界條件為兩端剛性固定，受均勻外壓作用，并給定材料的彈性模量為2×105MPa，泊松比為0.3。運用大型商用有限元分析軟件ANSYS Workbench進行非線性屈曲分析，通過給定不同的幾何初始缺陷、不同的材料屈服強度及彈塑性材料模型來分析研究幾何初始缺陷、材料的屈服強度及塑性硬化對結(jié)構(gòu)失穩(wěn)的臨界載荷的影響。

圖3 薄壁圓管算例

1.2 結(jié)果分析

1.2.1 初始幾何缺陷對屈曲臨界載荷的影響

在該算例中，分析時采用雙線性隨動硬化的彈塑性材料模型, 并假定材料的屈曲應力（以Ys表示）為200 MPa，塑性硬化階段材料的切線模量（以Et表示）為1500 MPa。需要說明的是，初始幾何缺陷系數(shù)的施加需基于ANSYS Workbench特征值屈曲分析的計算結(jié)果，特征值屈曲分析結(jié)果的第一階屈曲模態(tài)如圖4所示。

圖4 線性屈曲分析結(jié)果（第一階屈曲模態(tài)）

由圖4可知，該屈曲模態(tài)的最大位移約為1.12 mm，圓管的初始幾何缺陷可以通過缺陷系數(shù)乘以線性屈曲模態(tài)位移來施加。比如定義缺陷系數(shù)為0.1，其對應的最大初始缺陷為0.112 mm。

計算不同的幾何缺陷系數(shù)（Defect Factor，以DF表示）下的圓管在均勻外壓作用下的載荷位移曲線，如圖5所示。該圓管屈曲后的結(jié)構(gòu)變形如圖6所示。從該曲線可以看出，隨著載荷的增加，屈曲載荷位移曲線在到達某一臨界點后出現(xiàn)突變，該突變處所對應的載荷即為結(jié)構(gòu)的失穩(wěn)臨界載荷，在該臨界點之后，隨著載荷的微小增加，位移呈現(xiàn)急劇上升的趨勢。在初始幾何缺陷系數(shù)分別為0.100、0.010、0.001的情況下，該圓管的屈曲臨界載荷分別為12.4、18.8、20.5 MPa。可見，改善初始幾何缺陷可以有效提高結(jié)構(gòu)失穩(wěn)的臨界載荷。同時也可以看到，當幾何缺陷系數(shù)從0.01降低到0.001時，屈曲臨界載荷僅從18.6 MPa增加到20.5 MPa，這說明當初始幾何缺陷能夠控制在足夠小的范圍內(nèi)時，幾何缺陷的進一步降低對結(jié)構(gòu)失穩(wěn)臨界載荷的改善呈下降趨勢。

圖5 不同幾何缺陷系數(shù)DF下的載荷位移曲線

1.2.2 材料屈服強度對屈曲臨界載荷的影響

在該算例中，假定圓管的初始幾何缺陷系數(shù)（以DF表示）為0.01，分析時仍采用雙線性隨動硬化的彈塑性材料模型，并假定塑性硬化階段材料的切線模量Et為1500 MPa，通過給定三種不同的屈服應力水平（200、225、250 MPa）來研究材料的屈服強度對屈曲臨界載荷的影響。

計算不同的材料屈服強度Ys下的圓管在均勻外壓作用下的載荷位移曲線如圖7所示。

圖6 圓管屈曲后的結(jié)構(gòu)變形

圖7 不同材料屈服強度Ys下的載荷位移曲線

從圖7可以看出，在材料的屈服強度分別為200、225、250 MPa的情況下，該圓管的屈曲臨界載荷分別為18.8、20.8、22.5 MPa。由此可見，由于材料屈服強度的提高，使得結(jié)構(gòu)的承載能力更強，結(jié)構(gòu)失穩(wěn)的臨界載荷也會相應增大。

1.2.3 材料應變硬化對屈曲臨界載荷的影響

塑性材料在超過屈服階段之后，此時材料應力將會隨著應變的增大而不斷加大，這說明材料又恢復了抵抗變形的能力，要使其繼續(xù)變形，就必須增加載荷，這一現(xiàn)象稱為材料的強化（或稱應變硬化）[1]。材料經(jīng)過應變硬化階段后其屈服強度得以提高，根據(jù)1.2.2節(jié)的計算結(jié)果可以看出，經(jīng)過材料應變硬化過程的結(jié)構(gòu)，其抵抗失穩(wěn)的能力也會得以提高。

材料在屈服之后，應力應變曲線呈現(xiàn)非線性。在進行工程分析時，通常需要對應力應變曲線進行簡化，如前述的雙線性隨動硬化的彈塑性模型，將塑性應變硬化階段用直線代替，其對應的斜率以切線模量Et來表征，這樣與彈性階段（對應的斜率為材料的彈性模量E）共同組成一個雙線性的材料本構(gòu)關(guān)系，如圖8所示。

在實際工程應用中，某些材料在屈服后有明顯的屈服流動階段，當屈服流動階段比較長或者硬化程度比較小時可以忽略，這時可以簡化為理想彈塑性材料模型[8]，即應力達到屈服點之前，應力應變呈線性關(guān)系，應力達到屈服點之后，應力保持為常數(shù)σs。此時對應的切線模量Et=0。在本算例中，假定圓管的初始幾何缺陷系數(shù)（以DF表示）為0.01，并假定材料的屈服強度（以Ys表示）為200 MPa，通過給定三種不同的切線模量（0、1500、3000 MPa）來研究材料的塑性應變硬化階段的材料切線模量對屈曲臨界載荷的影響。

計算不同的材料塑性應變硬化階段的切線模量Et下的圓管在均勻外壓作用下的載荷位移曲線如圖9所示。

圖8 材料拉伸應力應變簡化曲線

圖9 不同材料切線模量Et 下的載荷位移曲線

從圖9可以看出，在材料的切線模量分別為0（即理想彈塑性材料）、1500、3000 MPa的情況下，該圓管的屈曲臨界載荷分別為18.6、18.8、19.1 MPa。由此可見，由于材料切線模量的增大（通常體現(xiàn)在材料具有更高的抗拉強度或更小的延伸率和韌性），材料呈現(xiàn)更短的屈服流動階段及更大的塑性硬化程度，結(jié)構(gòu)失穩(wěn)的臨界載荷也會相應增大。

1.3 實際案例

結(jié)合以上分析，對某實際空調(diào)產(chǎn)品中的換熱管設計進行了有限元仿真計算。在實際應用中，為了增加換熱管的換熱面積，提高換熱效果，通常采用的是螺旋翅片管（包括光管段、過渡段及翅片段），如圖10所示。

盡管在翅片加工時由于材料的擠壓變形，導致翅片段的最小壁厚小于光管段，但在有限元分析時候，考慮到翅片加工時材料的應變硬化行為，翅片段材料的屈服強度得以提高，其抵抗失穩(wěn)的載荷反而要高于光管段和過渡段，這也與實驗結(jié)果相吻合。實際的高效換熱管的樣件測試結(jié)果表明，塌陷（失穩(wěn)）通常發(fā)生在光管段和過渡段。而基于有限元分析得到結(jié)構(gòu)發(fā)生失穩(wěn)時的臨界載荷與試驗結(jié)果也十分接近（如表1）。

圖10 螺旋翅片換熱管

表1 實驗失穩(wěn)壓力與FEA計算結(jié)果對比

2 結(jié) 論

相比于線性特征值屈曲分析，非線性屈曲分析可以在計算模型中考慮結(jié)構(gòu)的幾何初始缺陷、材料的塑性行為、邊界條件的非線性及大撓度響應等特征，能較好地得到結(jié)構(gòu)的屈曲特性，在工程分析中越來越受到重視。考慮彈塑性材料的非線性屈曲分析還可以模擬結(jié)構(gòu)的塑性失穩(wěn)，從而可以進一步提高工程分析的精度。本文結(jié)合基于彈塑性材料的薄壁圓管非線性屈曲分析算例，研究了幾何初始缺陷、材料的屈服強度及應變硬化對結(jié)構(gòu)失穩(wěn)臨界載荷的影響，為空調(diào)管殼式換熱器中換熱管的設計提供了借鑒與指導。計算結(jié)果表明，在進行空調(diào)管殼式換熱器中換熱管的失穩(wěn)分析時，不僅要考慮結(jié)構(gòu)的初始幾何缺陷，必要時還需考慮材料的彈塑性力學行為（比如翅片換熱管加工時所發(fā)生的材料應變硬化行為）。在條件允許的情況下，建議通過拉伸試驗得到材料的實際應力應變曲線作為有限元分析的輸入，從而可以更準確地預測結(jié)構(gòu)失穩(wěn)的臨界載荷，達到降低實驗成本、優(yōu)化結(jié)構(gòu)設計的目的。