金属成型中CEL与各类网格重划结果的对比分析

1)CEL
詳細可以參考6.8 EF版 Abaqus Example Problems Manual 2.3.1 Revit forming.僅能於explicit中求解,由於CEL只能使用3D element EC3D8R,因此在Z軸方向設定對稱邊界來模擬平面應變,
 
由於Eulerian material的邊界並不是依據離散的element邊線,
 
而是依據每個element所佔material的體積率作平均及插值得到material的邊界,
 

因此在後處理時會因為數值平均得到近似值.
 
而在偵測接觸時也會受其影響而失準(如下圖變形體之右上角及右下內轉角所示).

CEL penetration.jpg

因此在許可的狀況下,建議在Lagrangian element與Eulerian element material的接觸邊界能夠預留一層Eulerian element,或者細化網格以減少誤差.
 

不過使用不同版本計算結果差異蠻大的,可能修改了程式,不同版本之比較如下圖所示,左圖是6.8.1之運算結果(看似較合理),右圖是6.8 EF之運算結果(ALE及map solution皆有此現象發生),至於哪個版本準確可能需要配合實驗作評斷.

CEL.gif


2)ALE

詳細可以參考Abaqus Analysis User's Manual 12.2 ALE adaptive meshing,
 
ALE在explicit使用的效果較為顯著,適用元素種類為1階減積分實體元素( 4-node quadrilaterals, 3-node triangles, 8-node hexahedra, and 4-node tetrahedra),,
 
本模型使用CPE4R element,在step module/ others中設定參數,
 
調試了幾次參數仍在成型後期產生單元零體積現象而終止計算(請高手賜教能否解決),
 
在這個成型模型中似不及CEL優越.
 
未設定ALE,部份element設定ALE及全部element設定ALE之動畫如下圖所示.

ALE.gif
3)remesh
詳細可以參考Abaqus Analysis User's Manual 12.3 Adaptive remeshing, remesh目前只能在standard中使用,而且需有幾何外型,無法使用導入網格,並且必需搭配ABAQUS/CAE使用, 元素種類適用於平面實體及殼元素( 4-node quadrilaterals, 3-node triangles)和Tetrahedrals實體元素, 本模型使用2D plane strain quad element,因此在切網格時需特別選擇mesh control為Quad-dominated/Advancing front,在mesh module/ adaptive中設定參數,在job module/ adaptive process manager中執行, ABAQUS/CAE會根據計算結果重切網格再做疊代計算.
 
使用remesh在這個成型模型中基本上難以收斂,element也扭曲得很厲害而使得結果不合理(請高手賜教能否解決),似不及CEL優越.
 
上傳各次疊代之動畫(remesh.gif)供參考(其中第3次疊代時,內部之element已扭曲超過表面之element而穿透剛體).

remesh.gif


4)Map solution

又稱為mesh-to-mesh solution mapping,詳細可以參考Abaqus Analysis User's Manual 12.4.1 Mesh-to-mesh solution mapping, 目前只能於standard中使用,而且在ABAQUS/CAE中不支援edit keyword,因此只能在inp文件中直接加指令,於command中執行.

目前適用元素種類為實體元素. 本模型使用2D plane strain quad element,初始網格在ABAQUS/CAE中切割, 使用map solution需設定restart request,當分析結果元素畸變太嚴重時,個人作法如下(如有錯誤或有更好的作法請高手指正)
 

  a.新建model, import/ part/ 選擇odb/ 選擇instance/ 勾選import deformed configuration/ 選擇step及frame
 

  b.ABAQUS/CAE似乎不支援無幾何之網格重劃,因此寫出此model之inp文件至hypermesh做網格重劃再import回ABAQUS/CAE
 

  特別注意的是,和import一樣,map solution的model,所有的interaction及load module都需要重新定義,
 

  step time和Amplitude要調整為剩下的分析時間. 位移荷載也要調整為剩下的分析位移量.
 

  最後寫出inp文件,在文件中加入map solution(詳細請參考help及附件)
 

  在ABAQUS command window中執行此inp,程式即會自動詢問所需的old job name.
 

  在使用上似乎若存在接觸,map solution就會遇到standard的收斂問題,而且不好解決.
 

  下圖在0.035 s導入map solution解決左半中段元素扭曲的問題,但很明顯做1次map並不足夠,
 

  有興趣者可以嘗試再接著做map解決到成型完成所有的網格畸變問題.
 

  Ps.建議使用**MAP SOLUTION, UNBALANCED STRESS=STEP指令,才能得到正確之反力,但似乎較難收斂.

map-0_1_join.gif

CEL的例子.
 

免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空