
注冊(cè)
/
登錄是這樣的,我想用UEL子程序來(lái)仿真斷裂失效,目前的計(jì)劃是利用abaqus設(shè)置狀態(tài)變量來(lái)實(shí)現(xiàn)單元刪除,但是我看好像只有UMAT或者VUMAT等少數(shù)子程序能夠通過(guò)關(guān)鍵字*Depvar來(lái)定義狀態(tài)變量實(shí)現(xiàn)單元刪除,沒(méi)有看到UEL相關(guān)的設(shè)置。所以我后來(lái)想能不能通過(guò)同時(shí)調(diào)用兩個(gè)子程序,用UMAT讀取我UEL程序輸出的狀態(tài)變量值,然后再umat里面判斷改單元是否失效。






在ABAQUS子程序文件.for中,使用了兩個(gè)子程序umeshmotion和fric,想實(shí)現(xiàn)在umeshtion中的數(shù)據(jù)也夠在計(jì)算中在fric程序中使用。在單個(gè)子程序中,使數(shù)據(jù)進(jìn)行保留,能夠實(shí)現(xiàn)在不同增量步中進(jìn)行數(shù)據(jù)傳遞是使用common對(duì)變量進(jìn)行定義就行。但我想將umeshmotion的數(shù)據(jù)在同一增量步時(shí)能用于fric中,因此可使用什么方法?
在ABAQUS中UEL隱式動(dòng)力學(xué)編寫時(shí),出現(xiàn)了網(wǎng)格劃分疏散,計(jì)算結(jié)果與abaqus原來(lái)計(jì)算結(jié)果一致,但是單元質(zhì)量矩陣輸出不同,如圖(左圖為ABAQUS自己計(jì)算,右圖為我編寫的UEL子程序計(jì)算結(jié)果) 若網(wǎng)格劃分密集,出現(xiàn)了計(jì)算結(jié)果完全不同的情況,如圖(左圖為ABAQUS自己計(jì)算,右圖為我編寫的UEL子程序計(jì)算結(jié)果) 對(duì)于單元剛度矩陣,也是如此



用相場(chǎng)法分析復(fù)合材料單胞的裂紋擴(kuò)展,基體用各向同性相場(chǎng),纖維用各向異性,umat里面可以用cmname關(guān)鍵詞在一個(gè)子程序里編輯多個(gè)材料本構(gòu),uel里面能不能也這樣,在一個(gè)uel子程序里面,編輯兩種相場(chǎng)模型,編輯兩種單元,用什么樣的關(guān)鍵詞呢
在兩個(gè)子程序中,duncunEB模型的UMAT子程序計(jì)算是沒(méi)問(wèn)題的,但是加上goodman單元后就計(jì)算不了了,給出的警告是這個(gè)問(wèn)題,請(qǐng)各位大佬幫忙看看到底是子程序還是哪里出了問(wèn)題。 dc_eb+goodman.for
自己編寫材料vumat,利用平面應(yīng)力單元模擬,但通過(guò)輸出應(yīng)變?cè)隽堪l(fā)現(xiàn)strainInc(k,3)在過(guò)程中始終等于零,所以最后3方向的應(yīng)變也始終為零。但是利用abaqus自帶材料模型模擬平面應(yīng)力狀態(tài)會(huì)發(fā)現(xiàn)3的應(yīng)變不為零。自己的子程序用三維單元驗(yàn)證是沒(méi)問(wèn)題。請(qǐng)問(wèn)大佬這是哪里出問(wèn)題?




重啟動(dòng)NROPT,FULLEKILL,ALL若按上述命令流執(zhí)行,重啟動(dòng)后在/POST1中輸入的命令將失效,無(wú)法選擇指定的單元上述命令是參照ansys Advanced Analysis Guide 中生死單元使用例子書寫的,原文如下 請(qǐng)問(wèn)哪里出現(xiàn)了問(wèn)題,或者能否利用其它方法在重啟動(dòng)后不進(jìn)入/POST1的情況下殺死溫度超過(guò)2000的單元,謝謝!
跳至頁(yè)
TOP