abaqus子程序Dload的主要作用：

（1）可用于定義作為位置、時間、單元編號、被加載積分點數量等的函數分布載荷大小的變化。

（2）在應力分析期間，將在每個基于單元或基于表面的非均勻分布載荷定義的載荷積分點處調用；

（3）將在每個積分點調用，以計算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；

（4）不能在基于模態的程序中用于描述負載的時間變化；并且忽略可能與相關聯的階躍定義或非均勻分布負載定義一起出現的任何幅度參考。

子程序接口界面

      SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,
     1 COORDS,JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION TIME(2), COORDS (3)
      CHARACTER*80 SNAME

      user coding to define F


      RETURN
      END

待定義變量

F：分布載荷的大小。表面載荷的單位為FL?2，體力的單位為FL?3。F將作為基于單元或基于表面的分布式載荷定義的一部分指定的載荷大小傳遞到例程中。如果未定義大小，F將作為零傳入。對于使用修正Riks法（靜態應力分析）的靜態分析，F必須定義為荷載比例系數λ的函數。分布式負載大小不可用于輸出目的。

用于傳遞信息的變量

KSTEP：Step 編號

KINC：增量數

TIME(1)：當前分析步對應的當前時間

TIME(2)：所有分析步對應的當前時間

NOEL：單元編號

NPT：根據荷載類型，構件內或構件表面上的荷載積分點編號。

LAYER：層編號

KSPT：當前層截面點編號

COORDS：包含被加載積分點坐標的數組。如果在步驟中考慮幾何非線性，則這些是當前坐標（參見定義分析）；否則，數組包含點的原始坐標。對于允許非軸對稱變形的軸對稱單元，COORDS（3）是積分點的角位置，單位為度

JLTYP：標識要調用DLOAD的加載類型。載荷類型可以是體力、基于表面的載荷或基于單元的表面載荷。對于基于單元的表面加載，此變量標識要調用DLOAD的元素面。當多個不同的非均勻分布載荷同時施加在一個元件上時，該信息非常有用。有關圖元面標識，請參見關于圖元庫。關鍵如下：

SNAME:基于曲面的載荷定義的曲面名稱（JLTYP=0）。對于體力或基于元素的表面載荷，表面名稱將作為空白傳遞。

案例演示

三維六面體中心圓形面載荷加載模擬

幾何模型20*20*1mm

材料采用純彈性的鋼參數，楊氏模量210Gpa，泊松比0.33，使用線性靜力分析，平面四周被完全固定，中心半徑為2.0mm區域施加100Mpa的表面力，并應用Dload子程序進行實現

結果與后處理

應力分布情況：

位移分布情況：

類似的可以生成其他形式的載荷，如隨著時間移動的載荷，如鐵軌上的火車的移動等