
注冊
/
登錄是這樣的,我想用UEL子程序來仿真斷裂失效,目前的計(jì)劃是利用abaqus設(shè)置狀態(tài)變量來實(shí)現(xiàn)單元刪除,但是我看好像只有UMAT或者VUMAT等少數(shù)子程序能夠通過關(guān)鍵字*Depvar來定義狀態(tài)變量實(shí)現(xiàn)單元刪除,沒有看到UEL相關(guān)的設(shè)置。所以我后來想能不能通過同時(shí)調(diào)用兩個(gè)子程序,用UMAT讀取我UEL程序輸出的狀態(tài)變量值,然后再umat里面判斷改單元是否失效。
在ABAQUS中UEL隱式動力學(xué)編寫時(shí),出現(xiàn)了網(wǎng)格劃分疏散,計(jì)算結(jié)果與abaqus原來計(jì)算結(jié)果一致,但是單元質(zhì)量矩陣輸出不同,如圖(左圖為ABAQUS自己計(jì)算,右圖為我編寫的UEL子程序計(jì)算結(jié)果) 若網(wǎng)格劃分密集,出現(xiàn)了計(jì)算結(jié)果完全不同的情況,如圖(左圖為ABAQUS自己計(jì)算,右圖為我編寫的UEL子程序計(jì)算結(jié)果) 對于單元剛度矩陣,也是如此


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


或者說子程序怎么提取單元的體積









跳至頁
TOP