多種材料的靜力學仿真

1. case描述:

兩個共面的立方體,大小為1*1*1.

第一個立方體E為1e8,P為0.3

第二個立方體E為1e7, P為0.4

一端四點固支,另一端一條邊上施加X=10000,Y=30000的分布荷載

2. 實現(xiàn)

2.1.FastModel中創(chuàng)建兩個立方體(圖1)

2.2.材料庫中創(chuàng)建兩種新材料

2.3.添加兩種3Dsolid單元,關(guān)聯(lián)新建的材料,分別賦給兩個立方體

2.4.選中一個立方體一端的四點,設(shè)置為約束為XYZ固定,選中另一個立方體的一邊,施加 X=10000,Y=30000的力

2.6.求解設(shè)置為static,結(jié)果輸出應(yīng)力和位移

2.7.分別導出Ansys/Patran/Abaqus格式

2.8.分別用Nastran/Radioss求解查看結(jié)果

圖1

Radioss 計算結(jié)果:

總結(jié):

1. Patran中對于多種材料的網(wǎng)格劃分比較麻煩,接觸部分需要通過 Equivalence,乃至手工編輯來保持共面的結(jié)點和單元的連續(xù)性,Fastmodel中則直接通過內(nèi)部運算保持了邊界的連續(xù)性,而且可以通過設(shè)置網(wǎng)格參數(shù),直接對不同的實體進行網(wǎng)格大小設(shè)置。

2. Nastran計算結(jié)果目測正確,而Radioss好像沒有把兩個實體當成一個,可能是bdf文件格式不識別,也可能需要設(shè)置額外的參數(shù)來支持不同材料的情況。

Stress計算結(jié)果 與 Ansys/Abaqus/Nastran 有少許誤差

Displacement 一樣

有可能是計算VonMises 時小數(shù)位數(shù)截斷誤差造成的