是這樣的，我想用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，題主目前需要利用dload子程序施加載荷，但是每個節(jié)點的載荷與此節(jié)點的相鄰節(jié)點的空間坐標(biāo)有關(guān)，題主嘗試過在dload中直接call GETVRN，但報錯jvgblock，后續(xù)發(fā)現(xiàn)GETVRN好像只能與UMESHMOTION配合使用，不知道吧里有沒有前輩做過相關(guān)的研究，或者其他的方式進行。金幣懸賞隨便填的，若解決了問題，可有償感謝！
各位老師好，我在進行abaqus用戶自定義子程序uel的開發(fā)，單元有兩個節(jié)點1、2，需要用到節(jié)點2的外力判斷是否發(fā)生滑移，也就是節(jié)點2的水平向的外力fx與μ×豎向外力fy之間的關(guān)系，因為這個單元也是應(yīng)用在一個大模型中的，所以節(jié)點的外力是不斷變化的，所以我就想獲取到節(jié)點2處的外力，進行一個判斷，還請各位老師進行指導(dǎo)，萬分感謝！
自己編寫材料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)變不為零。自己的子程序用三維單元驗證是沒問題。請問大佬這是哪里出問題？
在ABAQUS中UEL隱式動力學(xué)編寫時，出現(xiàn)了網(wǎng)格劃分疏散，計算結(jié)果與abaqus原來計算結(jié)果一致，但是單元質(zhì)量矩陣輸出不同，如圖（左圖為ABAQUS自己計算，右圖為我編寫的UEL子程序計算結(jié)果） 若網(wǎng)格劃分密集，出現(xiàn)了計算結(jié)果完全不同的情況，如圖（左圖為ABAQUS自己計算，右圖為我編寫的UEL子程序計算結(jié)果） 對于單元剛度矩陣，也是如此
26 跳至頁
工程師必備
- 項目客服
- 培訓(xùn)客服
- 平臺客服
TOP
