abaqus隱式計(jì)算很慢的案例
ABAQUS 二次開發(fā)（UEL+UMAT C3D8 單元?jiǎng)恿?em>隱式計(jì)算）
概述：采用UEL接口二次開發(fā)實(shí)現(xiàn)八節(jié)點(diǎn)單元，考慮BBAR修正，避免體積自鎖，對(duì)標(biāo)ABAQUS自帶的C3D8單元，計(jì)算的剛度矩陣、質(zhì)量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進(jìn)行應(yīng)力和應(yīng)變數(shù)據(jù)的可視化，計(jì)算的應(yīng)力應(yīng)變數(shù)據(jù)同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動(dòng)力簡(jiǎn)諧荷載為例，采用上述程序，應(yīng)用動(dòng)力隱式計(jì)算分析步，最終計(jì)算的位移、應(yīng)變等時(shí)程曲線均與ABAQUS保持一致。
（）模型信息
模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：
一面的所有節(jié)點(diǎn)均固定。另一面所有的節(jié)點(diǎn)施加簡(jiǎn)諧荷載，簡(jiǎn)諧荷載曲線為（詳細(xì)的參數(shù)見附件）：
計(jì)算的網(wǎng)格圖為：
網(wǎng)格尺寸為1，共計(jì)10x10x=1000個(gè)單元。總的計(jì)算時(shí)長(zhǎng)為1s，這只固定增量步長(zhǎng)為0.01s，所以總增量步數(shù)為100。
（）計(jì)算結(jié)果
以上面網(wǎng)格圖中中間角點(diǎn)為例，提取加載向位移時(shí)程曲線如下圖：
0.89s時(shí)刻x向應(yīng)力云圖為：
注：左一為ABAQUS計(jì)算結(jié)果，中間為umat實(shí)現(xiàn)線彈性各向同性本構(gòu)計(jì)算結(jié)果，右一為UEL+UMAT計(jì)算結(jié)果，后面的云圖也是這樣排布。
0.89s時(shí)刻y向應(yīng)力云圖為：
寫文字好麻煩，不寫了！放個(gè)動(dòng)圖算了，有興趣的移步附件。
展開 