做了一年多的複合材料壓力容器的分析工作,也積累了一些分析經驗,到了總結的時候了,回想起來,總最初採用i-deas,到到最後選定ansys爲自己的分析工具,確實有一些東西值得和大家分享,與從事複合材料結構分析的朋友門共同探討。
一)概述篇
複合材料是由一種以上具有不同性質的材料構成,其主要優點是具有優異的材料性能,在工程應用中典型的一種複合材料爲纖維增強複合材料,這種材料的特性表現爲正交各向異性,對於這種材料的類比,很多的程式都提供了一些處理方法,在i-deas、nastran、ansys中都有相應的處理方法。筆者最初是用i-deas下建立各項異性材料結合三維實體結構單元來類比(由於研究物件是厚壁容器,不宜採用殼單元),分析結果還是非常好的,而且i-deas強大的建模功能,但由於課題要求要進行壓力容器的優化分析,而且必須要自己寫優化程式,i-deas的二次開發功能開放性不是很強,所以改爲提供了一種非常好的二次開發編程語言pcl(以後在msc的版中專門給大家貼出這部分內容),採用patran結合nastran的分析環境,建立了基於正交各項異性和各項異性兩種分析模型,但最終發現,在得到的最後結果中,複合材料層之間的應力結果始終不合理,而模型是沒有問題的(因爲在i-deas中,相同的模型結果是合理的),於是最後轉向ansys,剛開始接觸ansys,真有相見恨晚的感覺,豐富的單元庫,開放的二次開發環境(apdl語言),下面就重點寫ansys的內容。
在ansys程式中,可以通過各項異性單元(solid 64)來類比,另外還專門提供了一類層合單元(layer elements)來類比層合結構(shell 99, shell 91, shell 181, solid 46 和solid 191)的複合材料。
採用ansys程式對複合材料結構進行處理的主要問題如下:
1)選擇單元類型
針對不同的結構和輸出結果的要求,選用不同的單元類型。
shell 99 ——線性結構殼單元,用於較小或中等厚度複合材料板或殼結構,一般長度方向和厚度方向的比值大於10;
shell 91 ——非線性結構殼單元,這種單元支援材料的塑性和大應變行爲;
shell 181—— 有限應變殼單元,這種單元支援幾乎所有的包括大應變在內的材料的非線性行爲;
solid 46 ——三維實體結構單元,用於厚度較大的複合材料層合殼或實體結構;
solid 191—— 三維實體結構單元,高精度單元,不支援材料的非線性和大變形。
2) 定義層屬性配置
主要是定義單層的層屬性,對於纖維增強複合材料,在這裏可以定義單層厚度、纖維方向等。
3)定義失效準則
支援多種失效準則,不過我還是沒有用他,而是自己寫了通過應力結果採用二次蔡胡準則程式來判斷的。
4)其他的一些建模技巧和後處理指導
在我的分析工作中,主要採用了三維實體結構單元。
關於solid 46單元
1)solid 46是用於類比複合材料厚殼或實體的8節點三維層合結構單元,單元節點有x,y和z方向三個結構自由度,單元允許最多250層不同的材料;
2)這種單元的定義包括:8個節點、各層厚度、各層材料方向角和正交各項異性材料屬性,其中每層可以爲面內兩個方向雙線性的不等厚層; (3)在材料定義時,只需定義材料主方向和材料坐標系(單元坐標系)一致的材料參數,不一致的複合材料層通過定義材料方向角(該層材料主方向和材料坐標系所成的角度)由程式自動轉換;
4)通過選擇不同的層直接在單元座標下獲取單元應力,包括三個方向的應力和麵內剪切應力,而不需要通過應力應變的轉換來獲取;
複合材料是一種各向異性材料,對於纖維增強複合材料又是一種正交各向異性材料,因此,在進行複合材料結構建模的時候要特別注意的一個重要的問題,就是材料的方向性。下面,就我個人的分析經驗,對複合材料結構的建模作一個總結。
1.結構坐標系、單元坐標系、材料坐標系和結果坐標系
建立複合材料結構模型,存在一個結構坐標系,用於確定幾何元素的位置,這個座標可以是笛卡爾坐標系、柱坐標系或者是球坐標系;單元坐標系是每個單元的局部坐標系,一般用來描述整個單元;材料坐標系是確定材料屬性方向的坐標系,一般沒有專門建立的材料坐標系,而是參考其他坐標系,如整體結構坐標系,或單元坐標系,在ansys程式中,材料座標是由單元座標唯一確定的,要確定材料座標,只要確定單元座標就行了;結果坐標系是在進行結果輸出時所使用的坐標系,也是一般參考其他坐標系。在ansys程式中,關於坐標系有人做過專門的總結。見附件。
2.用於複合材料結構分析的單元
用於複合材料分析的單元主要有兩類,一類是層合單元,如shell 99, shell 91, shell 181, solid 46 和solid 191;另一類是各向異性單元,如solid64;這些材料都有不同的處理方法,層合單元,在一個單元內可以包含多層資訊,包括各層的材料、厚度和方向;各項各向異性單元,在一個單元內,只能包含一種材料資訊,而且所得到的計算結果還要進行一些處理,因此有一定的局限性。
3.單元座標的一致性問題
在進行複合材料結構建模的時候,有些時候結構幾何比較複雜,很難用統一的座標來確定單元坐標系,即使對一些規則的幾何(如圓桶),在用旋轉方法生成幾何時,不同的面法向也會帶來單元座標的不一致,這就使得材料輸入的時候存在問題並使計算結果錯誤,因此,在幾何建模時要特別注意這一問題,筆者也沒有得到一些複雜幾何進行單元劃分時保持單元一致的合適方法。
4.一個實例
下面的命令流顯示了不同的幾何生成方法會産生不同的單元座標方向:
prep7
***create material***
mptemp,,,
mptemp,1,0
mpdata,ex,1,,2.068e8
mpdata,prxy,1,,0.29
mptemp,,,
mptemp,1,0
mpdata,dens,1,,7.82e-6
create element type
et,1,solid95
keyopt,1,1,1
keyopt,1,5,0
keyopt,1,6,0
keyopt,1,11,0
csys,1
hs=80
**create two keypoints along axial
k,101,0,0,0,
k,102,0,0,400,
**create keypoints
k,1,61,0,0,
k,2,hs,0,0,
k,5,100,0,0,
k,11,61,0,178,
k,12,hs,0,178,
k,15,hs+10,0,178,
k,111,61,0,178,
k,112,hs,0,178,
k,115,hs+10,0,178,
k,21,61,0,2450,
k,22,hs-4,0,2450,
k,25,hs+6,0,2450,
**create areas by keypoints
flst,2,4,3
fitem,2,21
fitem,2,111
fitem,2,112
fitem,2,22
a,p51x
flst,2,4,3
fitem,2,22
fitem,2,112
fitem,2,115
fitem,2,25
a,p51x
flst,2,2,5,orde,2
fitem,2,1
fitem,2,-2
flst,8,2,3
fitem,8,101
fitem,8,102
vrotat,p51xp51x, ,90,1,
type, 1
mat, 1
real,
esys, 0
secnum,
mshape,0,3d
mshkey,1
flst,5,2,6,orde,2
fitem,5,1
fitem,5,-2
cm,_y,volu
vsel, ,p51x
cm,_y1,volu
chkmsh,'volu'
cmsel,s,_y
vmesh,_y1
cmdele,_y
cmdele,_y1
cmdele,_y2
運行上述命令流,查看一下單元座標,再把命令流中下列部分
flst,2,4,3
fitem,2,21
fitem,2,111
fitem,2,112
fitem,2,22
a,p51x
改爲: flst,2,4,3
fitem,2,22
fitem,2,21
fitem,2,111
fitem,2,112
a,p51x
再看一下單元座標。
ansys坐標系總結
工作平面(working plane)
工作平面是創建幾何模型的參考(x,y)平面,在前處理器中用來建模(幾何和網格)
總體坐標系
在每開始進行一個新的ansys分析時,已經有三個坐標系預先定義了。它們位於模型的總體原點。三種類型爲:
cs,0: 總體笛卡爾坐標系
cs,1: 總體柱坐標系
cs,2: 總體球坐標系
資料庫中節點座標總是以總體笛卡爾坐標系,無論節點是在什麽坐標系中創建的。
局部坐標系
局部坐標系是用戶定義的坐標系。局部坐標系可以通過功能表路徑workplane>local cs>create lc來創建。
啟動的坐標系是分析中特定時間的參考系。缺省爲總體笛卡爾坐標系。當創建了一個新的坐標系時,新坐標系變爲啟動坐標系。
這表明後面的啟動坐標系的命令。功能表中啟動坐標系的路徑 workplane>change active cs to>。
節點坐標系
每一個節點都有一個附著的坐標系。節點坐標系缺省總是笛卡爾坐標系並與總體笛卡爾坐標系平行。節點力和節點邊界條件(約束)指的是節點坐標系的方向。
時間歷程後處理器 /post26 中的結果資料是在節點坐標系下表達的。而通用後處理器/post1中的結果是按結果坐標系進行表達的。
陶瓷基复合材料的制备原理与工艺
材料制备原理课程 题目陶瓷基复合材料的制备原理与工艺。学院材料科学与工程学院 专业班级。学生姓名。2012 年 3 月 28日。前言 科学技术的发展对材料提出了越来越高的要求,陶瓷基复合材料由于在破坏过程中表现出非脆性断裂特性,具有高可靠性,在新能源 国防军工 航空航天 交通运输等领域具有广阔的应用...
ANSYS与结构分析
3 强大的非线性分析能力 可进行几何非线性 材料非线性 接触非线性和单元非线性分析。其中,材料非线性包括压电材料和形状记忆合金等。4 强大的网格划分能力 可智能网格划分,根据几何模型的特点自动生成有限元网格。也可根据用户的要求,实现多种网格划分。5 良好的优化能力 通过ansys的优化设计功能,确实...
ANSYS结构分析指南疲劳
疲劳是指结构在低于静态极限强度载荷的重复载荷作用下,出现断裂破坏的现象。例如一根能够承受 300 kn 拉力作用的钢杆,在 200 kn 循环载荷作用下,经历 1,000,000 次循环后亦会破坏。导致疲劳破坏的主要因素如下 载荷的循环次数 每一个循环的应力幅 每一个循环的平均应力 存在局部应力集中...