
注冊
/
登錄是這樣的,我想用UEL子程序來仿真斷裂失效,目前的計劃是利用abaqus設置狀態變量來實現單元刪除,但是我看好像只有UMAT或者VUMAT等少數子程序能夠通過關鍵字*Depvar來定義狀態變量實現單元刪除,沒有看到UEL相關的設置。所以我后來想能不能通過同時調用兩個子程序,用UMAT讀取我UEL程序輸出的狀態變量值,然后再umat里面判斷改單元是否失效。
各位老師好,我在進行abaqus用戶自定義子程序uel的開發,單元有兩個節點1、2,需要用到節點2的外力判斷是否發生滑移,也就是節點2的水平向的外力fx與μ×豎向外力fy之間的關系,因為這個單元也是應用在一個大模型中的,所以節點的外力是不斷變化的,所以我就想獲取到節點2處的外力,進行一個判斷,還請各位老師進行指導,萬分感謝!





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


按理來說除了約束住位移的節點之外,其他節點的z方向位移應該都指向z軸反方向且由大變小,abaqus自帶單元算的結果正是如此;但是我編寫的UEL子程序,除了約束住位移的節點之外,其他節點的z方向位移指向z軸反方向和正方向都存在,且節點位移也不是由大到小,很亂沒有規律。我的問題就是以上內容,望有大佬能指點一二。萬分感謝!!!!

請問我開發了UEL單元做了如圖所示的接觸模型,邊界條件上端施加位移,下端固定,結果只有壓頭有應力數值,下端的模型如圖發生了變形,卻沒有應力,是input文件的問題還是for文件的問題呢?又或者是其他方面的問題,還請老師不吝賜教






跳至頁
TOP