今天來解讀一個在ABAQUS中使用USDFLD子程序實現損傷計算的例子。

USDFLD子程序的界面如下：

需要定義的變量

數組FIELD（NFIELD）包含當前材料（積分）點的場變量fi。

——在當前增量步的末尾，這些傳入了由節點插值得到的值，由*INITIAL CONDITIONS選項或*FIELD選項定義。

——更新的fi用于計算材料屬性值，其作為場變量的函數。更新的fi傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點可以被調用。

可能需要定義的變量

• 包含依賴于解的狀態變量的數組STATEV（NSTATV），能夠在USDFLD中定義。

——這些在增量步的開始作為值傳入。

——在該子程序中STATEV在所有情況下都能夠更新，所有更新的值被傳遞到其他用戶子程序中（CREEP、HETVAL、UEXPAN，UHARD、UHYPEL、UMAT、UMATHT和UTRS），在該材料點可以被調用。

——與該材料點相關的狀態變量的數量由*DEPVAR選項定義。

• 建議的新時間步長與現在正在使用的時間步長（DTIME）的比例PNEWDT，可以給定。

——該變量允許用戶在ABAQUS中的自動時間增量步長算法中提供輸入。

只提供信息的變量

• NFIELD 在該點存在的場變量個數

• DIRECT 與當前積分點相關的材料方向在全局坐標系下的方向余弦

——DIRECT(#，1)定義第一個材料方向

• T 殼或膜材料方向相對于單元法向的任意轉動的方向余弦

等等…..

例子：層壓復合板的失效

該問題模擬了中心有孔的層壓復合板在受到面內壓縮時發生的損傷。

板由石墨-環氧樹脂層組成，纖維方向為（-45/45）疊層。

該問題使用了四分之一對稱有限元模型，使用兩層CPS4單元，而不是用殼單元對復合板進行建模，因為板的厚度足夠大，面外位移應該最小。