概述：采用UEL接口二次開發實現八節點單元，考慮BBAR修正，避免體積自鎖，對標ABAQUS自帶的C3D8單元，計算的剛度矩陣、質量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進行應力和應變數據的可視化，計算的應力應變數據同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動力簡諧荷載為例，采用上述程序，應用動力隱式計算分析步，最終計算的位移、應變等時程曲線均與ABAQUS保持一致。

（）模型信息

模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：

一面的所有節點均固定。另一面所有的節點施加簡諧荷載，簡諧荷載曲線為（詳細的參數見附件）：

計算的網格圖為：

網格尺寸為1，共計10x10x=1000個單元。總的計算時長為1s，這只固定增量步長為0.01s，所以總增量步數為100。

（）計算結果

以上面網格圖中中間角點為例，提取加載向位移時程曲線如下圖：

0.89s時刻x向應力云圖為：

注：左一為ABAQUS計算結果，中間為umat實現線彈性各向同性本構計算結果，右一為UEL+UMAT計算結果，后面的云圖也是這樣排布。

0.89s時刻y向應力云圖為：

寫文字好麻煩，不寫了！放個動圖算了，有興趣的移步附件。

加載向位移時程動圖為：

附件：

-----------------------
-----------------------
ABAQUS
  JOB-1.INP
  RUN.BAT
  UM.CAE
  UM.JNL
UEL-UMAT
  BBAR.OBJ
  JOB-1.INP
  RUN.BAT
UMAT
  UM.CAE
  JOB-1.INP
  UM.JNL
  UMAT.FOE
  RUN.BAT
-----------------------
-----------------------