
注冊
/
登錄

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




在ABAQUS中UEL隱式動力學(xué)編寫時,出現(xiàn)了網(wǎng)格劃分疏散,計算結(jié)果與abaqus原來計算結(jié)果一致,但是單元質(zhì)量矩陣輸出不同,如圖(左圖為ABAQUS自己計算,右圖為我編寫的UEL子程序計算結(jié)果) 若網(wǎng)格劃分密集,出現(xiàn)了計算結(jié)果完全不同的情況,如圖(左圖為ABAQUS自己計算,右圖為我編寫的UEL子程序計算結(jié)果) 對于單元剛度矩陣,也是如此

用相場法分析復(fù)合材料單胞的裂紋擴展,基體用各向同性相場,纖維用各向異性,umat里面可以用cmname關(guān)鍵詞在一個子程序里編輯多個材料本構(gòu),uel里面能不能也這樣,在一個uel子程序里面,編輯兩種相場模型,編輯兩種單元,用什么樣的關(guān)鍵詞呢
在兩個子程序中,duncunEB模型的UMAT子程序計算是沒問題的,但是加上goodman單元后就計算不了了,給出的警告是這個問題,請各位大佬幫忙看看到底是子程序還是哪里出了問題。 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)變不為零。自己的子程序用三維單元驗證是沒問題。請問大佬這是哪里出問題?



或者說子程序怎么提取單元的體積
我在用Cohesive PPR本構(gòu)進行計算的時候,發(fā)現(xiàn)了一個奇怪問題?問題如下:我看了PPR本構(gòu)中需要9個參數(shù),于是我創(chuàng)建好模型生成inp文件后,我更改UEL PROPERTY中輸入?yún)?shù) ,當(dāng)這九個參數(shù)同時放在一行時,程序就會出問題,而當(dāng)把最后一個參數(shù)放到第二行時,程序就沒問題。這個我很費解。難道inp文件每行最大只能輸入8個參數(shù),如果大于8個就需要換行?


跳至頁
TOP