在工程實(shí)際應(yīng)用中，我們經(jīng)常會(huì)遇到移動(dòng)載荷的例子，如車轍實(shí)驗(yàn)，汽車過(guò)橋等。今天就介紹一下用abaqus子程序Dload實(shí)現(xiàn)移動(dòng)載荷。

Dolad：隱式求解算法

Vdload：顯示求解算法

（1）Dolad簡(jiǎn)介

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

F	需要自己定義的載荷
KSTEP	載荷步數(shù)(如第一個(gè)載荷步、第二個(gè)載荷步)
KINC	增量步數(shù)
TIME	Time(1)載荷步時(shí)間、Time(2)總時(shí)間
NOEL	單元編號(hào)
NPT	單元積分點(diǎn)編號(hào)
LAYER、KSPT	和LAYER SOILD相關(guān)
COORDS	積分點(diǎn)坐標(biāo)
JLTYP	載荷類型
SNAME	名稱

（2）案例

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

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為移動(dòng)速度

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

xMin=Disp+0.002 !0.002 設(shè)置誤差范圍，積分點(diǎn)坐標(biāo)和節(jié)點(diǎn)坐標(biāo)不一樣

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

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

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

ELSE

F=0

END IF

RETURN

END

（3）Abaqus操作

1、用general static?分析步設(shè)置成fixed?增量步這樣可以等間隔的時(shí)間輸出結(jié)果

2、定義載荷load的時(shí)候選pressure?在distribution中選user-defined

3、建立job 在general選項(xiàng)中 添加Dload.for 文件目錄

4、結(jié)果

（5）總結(jié)

Abaqus的Dload子程序讓用戶可以自定義各種復(fù)雜的非均布載荷，比如和坐標(biāo)或者時(shí)間相關(guān)的載荷，其實(shí)只要知道了Dload各參數(shù)的意思，用起來(lái)也比較得心應(yīng)手。

?

?

DLOAD子程序說(shuō)明

————————————————————————————————————?

Use this subroutine to define nonuniform, distributed
mechanical loads (pressures and body forces).

使用此子程序定義非均勻的分布式機(jī)械負(fù)載（壓力和體力）。

（

當(dāng)負(fù)載是時(shí)間和/或位置的復(fù)雜函數(shù)時(shí)，通常使用用戶子程序DLOAD。

通常可以用* AMPLITUDE選項(xiàng)建模簡(jiǎn)單的時(shí)間函數(shù)。

子程序也可以用來(lái)定義一個(gè)隨元素號(hào)和/或積分點(diǎn)數(shù)而變化的負(fù)載。

）

————————————————————————————————————?

一般用法，在INP文件中定義載荷施加的位置（關(guān)鍵字），求解時(shí)讓其與子DLOAD子程序關(guān)聯(lián)。

（

當(dāng)* DLOAD或* DSLOAD選項(xiàng)包含不均勻的加載類型標(biāo)簽時(shí)，子程序被調(diào)用。

）

*DLOAD 1, P4NU

（Apply nonuniform DLOAD to face 4 of element 1，）

（

“ABAQUS /Standard”手冊(cè)中給出了可用于任何特定單元的非均勻分布載荷類型列表。

）

DLOAD子程序接口

?

SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,1 COORDS,JLTYP,SNAME) CINCLUDE 'ABA_PARAM.INC' CDIMENSION TIME(2), COORDS (3)CHARACTER*80 SNAMEuser coding to define FRETURNEND

?

?

子程序中的各個(gè)參數(shù)的意義可以參考ABAQUS幫助文檔

F——分布負(fù)載的大小

KSTEP——分析步

TIME(1)——步進(jìn)時(shí)間的當(dāng)前值或負(fù)載比例因子的當(dāng)前值

TIME(2)——當(dāng)前總時(shí)間值

NOEL——單元編號(hào)

NPT——根據(jù)負(fù)載類型，在元素或元素表面加載積分點(diǎn)數(shù)。

LAYER——層數(shù)

KSPT——當(dāng)前圖層中的Section point

COORDS——積分點(diǎn)坐標(biāo)（數(shù)組）

JLTYP——DLOAD調(diào)用的載荷類型的標(biāo)識(shí)

?SNAME——表面加載定義的表面名稱（JLTYP = 0）。 對(duì)于身體力或基于單元的表面載荷，表面名稱以空白形式傳遞

供參考的子程序

————————————————————————————————————?

SUBROUTINE DLOAD(F, KSTEP, KINC, TIME ,NOEL, NPT,1 LAYER, KSPT, COORDS, JLTYP, SNAME) C C EXPONENTIAL PRESSURE LOAD CINCLUDE ‘ABA_PARAM.INC’ CDIMENSION COORDS(3),TIME(2)CHARACTER*80 SNAMEDATA TEN,ONE,CONST /10.,1.,-23.03/F=TEN*(ONE-(EXP(CONST*TIME(1))))IF(NPT.EQ.1) WRITE(6,*) ‘ LOAD APPLIED’,F,’AT TIME=’,TIME(1)RETURNEND

?————————————————————————————————————?

?（

用戶只能定義變量F 這是分布式負(fù)載的大小，單位取決于分布式負(fù)載的類型：

FL–1 對(duì)于沿一維（梁）單元的線載荷，

FL–2 表面負(fù)荷（例如壓力）和

FL–3 體力（例如重力，向心力，加速度）。

）

?

轉(zhuǎn)載于:https://www.cnblogs.com/jiangkejie/p/10231775.html

總結(jié)

以上是生活随笔為你收集整理的Abaqus中施加移动车辆荷载（待整理）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。