基于ABAQUS復合材料薄壁圓筒的屈曲分析

1、基于ABAQUS復合材料薄壁圓筒的屈曲分析由于玻璃鋼復合材料的薄壁圓筒結(jié)構(gòu)具有強度高、重量輕、剛度大、耐腐蝕，電絕緣及透微波等優(yōu)點，目前已廣泛應用于航空航天和民用領域中。工程中廣泛使用的這些薄壁圓筒，當它們受壓縮、剪切、彎曲和扭轉(zhuǎn)等荷載作用時，最常見的失效模式為屈曲。因此，為了保證結(jié)構(gòu)的安全，需要進行屈曲分析。對結(jié)構(gòu)進行屈曲分析，涉及到較復雜的彈(塑)性理論和數(shù)學計算，要通過求解高階偏微分方程組，才能求解失穩(wěn)臨界荷載，而且只有少數(shù)簡單結(jié)構(gòu)才能求得精確的解析解。因此，只能采用能量法、數(shù)值方法和有限元方法等近似的分析方法進行分析。近20年來，隨著計算機和有限元方法的迅猛發(fā)展，形成了許多的實用分析程

2、序，提高了對復雜結(jié)構(gòu)進行屈曲分析的能力和設計水平。ABAQUS就是其中的杰出代表。1.屈曲有限元理論有限元方法中，對結(jié)構(gòu)的屈曲失穩(wěn)問題的分析方法主要有兩類：一類是通過特征值分析計算屈曲載荷，另一類是利用結(jié)合NewtonRaphson迭代的弧長法來確定加載方向，追蹤失穩(wěn)路徑的幾何非線性分析方法，能有效分析高度非線性屈曲和后屈曲問題。1.1線性屈曲假設結(jié)構(gòu)受到的外載荷模式為P0。，幅值大小為，結(jié)構(gòu)內(nèi)力為Q，則靜力平衡方程應為P0=Q進一步考察結(jié)構(gòu)在(+)P0載荷作用下的平衡方程，得到KE+KSS+S+KGu+uu=P0由于結(jié)構(gòu)達到保持穩(wěn)定的臨界載荷時有，代入上式得KE+KS+KGuu=0該方程對應

3、的特征值問題為detKE+KS+KGu=0如果忽略幾何剛度增量的影響，屈曲分析的方程又可進一步簡化為detKE+KS=0該方程即為求解線性屈曲的特征值方程。為屈曲失穩(wěn)載荷因子，u為結(jié)構(gòu)失穩(wěn)形態(tài)的特征向量。1.2非線性屈曲非線性屈曲分析方法多采用弧長法進行分步迭代計算，在增量非線性有限元分析中，沿著平衡路徑迭代位移增量的大小(也叫弧長)和方向，確定載荷增量的自動加載方案，可用于高度非線性的屈曲失穩(wěn)問題。與提取特征值的線性屈曲分析相比，弧長法不僅考慮剛度奇異的失穩(wěn)點附近的平衡，而且通過追蹤整個失穩(wěn)過程中實際的載荷、位移關系，獲得結(jié)構(gòu)失穩(wěn)前后的全部信息，適合于高度非線性的屈曲失穩(wěn)問題。2.ABAQU

4、S的線性屈曲分析ABAQUS中提供兩種分析方法來確定結(jié)構(gòu)的臨界荷載和結(jié)構(gòu)發(fā)生屈曲響應的特征形狀：線性屈曲分析(特征值屈曲分析)、非線性屈曲分析。線性屈曲分析用于預測一個理想的彈性結(jié)構(gòu)的理論屈曲強度。它是預期的線性屈曲荷載的上限，可以作為非線性屈曲分析的給定荷載，在漸進加載達到此荷載前，非線性求解必然發(fā)散；它還可以作為施加初始缺陷或擾動荷載的依據(jù)。所以預先進行特征值屈曲分析有助于非線性屈曲分析，進行特征值屈曲分析是必要的。3.算例3.1問題概述圖3-1 實例模型如圖所示兩端開口的復合材料薄壁圓筒，底端固支，頂端作用有均勻分布的軸壓邊載。半徑R=152mm，高度300mm，厚度t=0.804mm，

5、對稱鋪層45,0s，單層厚度為0.134mm。復合材料圓筒的材料參數(shù)如下表：表1 AS4/3501-6 石墨/環(huán)氧的彈性參數(shù) 彈性模量參數(shù)E1/GPa142E2/GPa9.7G12/GPa6G13/GPa3.6G23/GPa3.60.33.2建模、劃分網(wǎng)格及分析過程3.2.1 ABAQUS進行前處理建立模型首先先繪制一個半徑為152mm的三維輪廓圖，然后進行拉伸300mm創(chuàng)建一個具有平面殼體單元的三維變形體。其次進入property模塊生成材料，參考上表編輯其彈性參數(shù)。編輯截面進行復合層板的鋪層完成材料屬性、截面屬性的定義。圖3-2 鋪層 接著生成裝配件，定義屈曲分析步。對所建模型施加邊界條件

6、和施加載荷。將圓筒的一端固支，另一端施加單位1的均布軸邊壓載。最后進行網(wǎng)格劃分?？紤]到在復合材料層合殼模型中剪切柔度的影響，可采用厚殼單元S4R來模擬它。所得到的模型如圖所示：圖3-3 模型網(wǎng)格圖計算結(jié)果分析在JOD模塊中建立屈曲分析模塊進行分析，可得到薄壁圓筒的六階屈曲失穩(wěn)載荷因子。表2 各階模態(tài)的屈曲載荷因子 模態(tài)屈曲載荷因子148.015248.015348.98448.981550.188650.188則取第一階模態(tài)的屈曲載荷因子計算臨界載荷。由屈曲載荷因子與臨界載荷的關系：Pcr=*p當對薄壁圓筒施加48.015N/mm的邊載時，圓筒將屈曲失穩(wěn)。對應的各階模態(tài)云圖如下所示：模態(tài)1位移

7、軸向和截面方向云圖模態(tài)2位移軸向和截面方向云圖模態(tài)3位移軸向和截面方向云圖模態(tài)4位移軸向和截面方向云圖模態(tài)5位移軸向和截面方向云圖模態(tài)6位移軸向和截面方向云圖圖3-4 各階模態(tài)云圖 由上圖可以看出當施加一階模態(tài)的載荷時，圓筒已經(jīng)屈曲失穩(wěn)。所以保留一階的臨界載荷更有現(xiàn)實意義。4.總結(jié)由算例的計算過程可以得到屈曲載荷因子的大小與劃分的網(wǎng)格數(shù)和網(wǎng)格類型有關，網(wǎng)格劃分過小會導致計算步驟過多影響計算速度甚至無法計算出結(jié)果，網(wǎng)格過大導致結(jié)果不精確。所以在選擇網(wǎng)格時應選擇適當選擇類型和大小。參考文獻1 呂桂萍、趙冰燃，水輪機圓筒閥筒體屈曲分析，閥門，2011年第5期.2 劉力濤，蘇文獻,內(nèi)壓和軸向載荷作用下開孔薄壁短圓筒屈曲的數(shù)值研究,石油化工設備，第37卷，第