在工程實際應用中，我們經常會遇到移動載荷的例子，如車轍實驗，汽車過橋等。今天就介紹一下用abaqus子程序Dload實現移動載荷——“CAE仿真實驗室”出品。

（1）Dolad簡介

     Dload 函數體和各參數物理意義如下:  SUBROUTINEDLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT, COORDS,JLTYP,SNAME) 

F	需要自己定義的載荷
KSTEP	載荷步數(如第一個載荷步、第二個載荷步)
KINC	增量步數
TIME	Time(1)載荷步時間、Time(2)總時間
NOEL	單元編號
NPT	單元積分點編號
LAYER、KSPT	和LAYER SOILD相關
COORDS	積分點坐標
JLTYP	載荷類型
SNAME	名稱

（2）案例

    如圖所示，物體表面作用一個壓力載荷P(t)，同事這個載荷在面內以一定速度V移動

Dload.for子程序源碼

SUBROUTINEDLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,

    1 COORDS,JLTYP,SNAME)

     INCLUDE 'ABA_PARAM.INC'

     DIMENSION TIME(2), COORDS (3)

     CHARACTER*80 SNAME

     PARAMETER(P0=5e4,Dm0=250,Dmiu=0.002,speed=80/3.6,DL=6.0,width=0.2)

     Pi=3.1415926535897932384626433832795

     omega=2.0*Pi*speed/DL !角速度 speed為移動速度

     Disp=TIME(1)*speed  !移動距離

     xMin=Disp+0.002   !0.002 設置誤差范圍，積分點坐標和節點坐標不一樣

     xMax=Disp+width-0.002  ! width 加載寬度

     IF(COORDS(1).ge.xMin.and.COORDS(1).lt.xMax)THEN  !判斷坐標位置

     F=P0+Dm0*Dmiu*omega*omega*sin(omega*TIME(1))      !定義載荷大小

     ELSE

     F=0

     END IF

     RETURN

     END 

（3）Abaqus操作

    1、用general static 分析步  設置成fixed 增量步  這樣可以等間隔的時間輸出結果

 

      2、定義載荷 load的時候  選 pressure   在distribution中選  user-defined

    3、建立job 在general選項中 添加Dload.for 文件目錄

4、簡單結果

（5）總結

   Abaqus的Dload子程序讓用戶可以自定義各種復雜的非均布載荷，比如和坐標或者時間相關的載荷，其實只要知道了Dload各參數的意思，用起來也比較得心應手。最后小伙伴們可以長按掃描下面二維碼關注我的公眾號‘CAE仿真實驗室’哦，我們將定期給大家分享有限元仿真干貨,謝謝~