Abaqus用Dload子程序?qū)崿F(xiàn)移動載荷
2017-03-15 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
在工程實(shí)際應(yīng)用中,我們經(jīng)常會遇到移動載荷的例子,如車轍實(shí)驗(yàn),汽車過橋等。
今天就介紹一下Abaqus用Dload子程序實(shí)現(xiàn)移動載荷——“CAE仿真實(shí)驗(yàn)室”出品,歡迎轉(zhuǎn)載。
(1)Dload簡介
Dload 函數(shù)體和各參數(shù)物理意義如下: SUBROUTINEDLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT, COORDS,JLTYP,SNAME)
F |
需要自己定義的載荷 |
KSTEP |
載荷步數(shù)(如第一個載荷步、第二個載荷步) |
KINC |
增量步數(shù) |
TIME |
Time(1)載荷步時間、Time(2)總時間 |
NOEL |
單元編號 |
NPT |
單元積分點(diǎn)編號 |
LAYER、KSPT |
和LAYER SOILD相關(guān) |
COORDS |
積分點(diǎn)坐標(biāo) |
JLTYP |
載荷類型 |
SNAME |
名稱 |
(2)案例
如圖所示,物體表面作用一個壓力載荷P(t),同事這個載荷在面內(nèi)以一定速度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 設(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 增量步 這樣可以等間隔的時間輸出結(jié)果
2、定義載荷 load的時候 選 pressure 在distribution中選 user-defined
3、建立job 在general選項(xiàng)中 添加Dload.for 文件目錄
(4)、簡單結(jié)果
(5)總結(jié)
Abaqus的Dload子程序讓用戶可以自定義各種復(fù)雜的非均布載荷,比如和坐標(biāo)或者時間相關(guān)的載荷,其實(shí)只要知道了Dload各參數(shù)的意思,用起來也比較得心應(yīng)手。
相關(guān)標(biāo)簽搜索:Abaqus用Dload子程序?qū)崿F(xiàn)移動載荷 abaqus分析培訓(xùn) abaqus技術(shù)教程 abaqus巖土分析 鋼筋混凝土仿真 abaqus分析理論 abaqus軟件下載 abaqus umat用戶子程序編程 Abaqus代做 Abaqus基礎(chǔ)知識 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn)