殼到實體子模型流程

1、殼到實體子模型流程雷挺 由殼到實體的子模型技術(shù)在以下情況下不能使用，SHELL91 (KEYOPT(11) 0), SHELL99 (KEYOPT(11) 0), SHELL181 （截面偏移）,或者SHELL281 （截面偏移）。 以車體結(jié)構(gòu)由殼到實體的子模型的流程為例，車體結(jié)構(gòu)總體模型由184488個殼單元組成，節(jié)點數(shù)為153722個，單元類型為SHELL181單元。門角子模型為solid185單元，共有單元數(shù)357088個，單元節(jié)點數(shù)量為775425。殼和實體結(jié)構(gòu)都采用全積分模式進行計算，殼選擇，實體單元在單元類型中選擇k2=0; k3=0; k6=0; k10=0; 所有仿真使用的材料

2、參數(shù)相同，子模型的空間坐標位于總體模型相應(yīng)位置處。子模型的邊界在相應(yīng)的總體模型中也應(yīng)該存在，即子模型的邊界為一些平面，這些平面大致垂直于總體模型的殼結(jié)構(gòu)，同時總體模型在邊界處也必須有相應(yīng)直線分割總模型的平面，以確保子模型的邊界能從總體模型中正確的讀取邊界條件。仿真流程 打開總體模型，對總體模型進行加載，進行計算，得出計算結(jié)果。圖1：總體有限元模型圖2：總體模型計算結(jié)果 打開事先建立好的子模型，選出邊界點，建立邊界點文件。圖3：局部子模型圖4：選出子模型的邊界點圖5：建立點文件圖6：點文件全部建立 打開總體模型，進入通用后處理界面，讀入總體模型的計算結(jié)果。圖7：打開總體模型圖8：顯示計算結(jié)果圖9

3、：子模型在總體模型中的位置 如圖9所示，總體模型使用的帶殼厚度特征的顯示。現(xiàn)在通用后處理菜單中出現(xiàn)了submodeling選項，這時候就可以通過submodeling選項下的interpolate DOF和interp Body Forc讀取子模型的邊界自由度和邊界載荷力了。首先讀取邊界節(jié)點自由度，然后使用相同方法讀取邊界載荷力。圖10：從總體模型中讀取截取點的位移自由度圖11：位移自由度讀取成功使用相同的方法讀取邊界載荷力，讀取邊界條件時，需要確保材料類型，子模型的空間坐標都與總體模型一致。然后進行邊界載荷文件編輯。首先編輯.CBDO文件。查找CB1,把CB1前面的/EOF刪除，保存更改文件

4、。打開.BFIN文件，如圖13可以看到里面只有一個/EOF，而且所有邊界載荷力為0，在本例中可以不用讀取.BFIN。 進入子模型文件中，用input方式讀入剛才的.BFIN，.CBDO文件，在讀入.BFIN文件時，沒有任何變化，讀入編輯后的.CBDO文件時，子模型邊界載荷自動添加。圖12：編輯.cbdo文件圖13：.BFIN文件的載荷力值為0圖14：讀入邊界條件圖15：讀入邊界條件后的模型圖16：計算結(jié)果邊界上的節(jié)點位移使用插值函數(shù)的方法進行計算，計算結(jié)果的精確度誤差較大，所以只選擇我們關(guān)心的園角處的應(yīng)力。圖17：圓角處子模型計算結(jié)果應(yīng)力值101MPa 從圖形結(jié)果顯示可以看出，結(jié)構(gòu)的應(yīng)力梯度顯示的也很平滑。圖18：總體模型中計算結(jié)果應(yīng)力值為90MPa圖19：打開殼厚度特征的計算結(jié)果顯示圖20：單元應(yīng)力顯示結(jié)果為99MPa總體模