abaqus子程序Dload的主要作用：

（1）可用于定義作為位置、時(shí)間、單元編號(hào)、被加載積分點(diǎn)數(shù)量等的函數(shù)分布載荷大小的變化。

（2）在應(yīng)力分析期間，將在每個(gè)基于單元或基于表面的非均勻分布載荷定義的載荷積分點(diǎn)處調(diào)用；

（3）將在每個(gè)積分點(diǎn)調(diào)用，以計(jì)算承受不均勻荷載類型PENU和PINU的管道元件的有效軸向力ESF1；

（4）不能在基于模態(tài)的程序中用于描述負(fù)載的時(shí)間變化；并且忽略可能與相關(guān)聯(lián)的階躍定義或非均勻分布負(fù)載定義一起出現(xiàn)的任何幅度參考。

子程序接口界面

      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(xiàn)將作為零傳入。對(duì)于使用修正Riks法（靜態(tài)應(yīng)力分析）的靜態(tài)分析，F(xiàn)必須定義為荷載比例系數(shù)λ的函數(shù)。分布式負(fù)載大小不可用于輸出目的。

用于傳遞信息的變量

KSTEP：Step 編號(hào)

KINC：增量數(shù)

TIME(1)：當(dāng)前分析步對(duì)應(yīng)的當(dāng)前時(shí)間

TIME(2)：所有分析步對(duì)應(yīng)的當(dāng)前時(shí)間

NOEL：?jiǎn)卧幪?hào)

NPT：根據(jù)荷載類型，構(gòu)件內(nèi)或構(gòu)件表面上的荷載積分點(diǎn)編號(hào)。

LAYER：層編號(hào)

KSPT：當(dāng)前層截面點(diǎn)編號(hào)

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

JLTYP：標(biāo)識(shí)要調(diào)用DLOAD的加載類型。載荷類型可以是體力、基于表面的載荷或基于單元的表面載荷。對(duì)于基于單元的表面加載，此變量標(biāo)識(shí)要調(diào)用DLOAD的元素面。當(dāng)多個(gè)不同的非均勻分布載荷同時(shí)施加在一個(gè)元件上時(shí)，該信息非常有用。有關(guān)圖元面標(biāo)識(shí)，請(qǐng)參見關(guān)于圖元庫(kù)。關(guān)鍵如下：

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

案例演示

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

幾何模型20*20*1mm

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

結(jié)果與后處理

應(yīng)力分布情況：

位移分布情況：

類似的可以生成其他形式的載荷，如隨著時(shí)間移動(dòng)的載荷，如鐵軌上的火車的移動(dòng)等