Abaqus中施加移动车辆荷载（待整理）
在工程實(shí)際應(yīng)用中，我們經(jīng)常會(huì)遇到移動(dòng)載荷的例子，如車(chē)轍實(shí)驗(yàn)，汽車(chē)過(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 | 載荷類(lèi)型 |
| SNAME | 名稱(chēng) |
（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子程序讓用戶(hù)可以自定義各種復(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í)，通常使用用戶(hù)子程序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)包含不均勻的加載類(lèi)型標(biāo)簽時(shí)，子程序被調(diào)用。
）
*DLOAD 1, P4NU（Apply nonuniform DLOAD to face 4 of element 1，）
（
“ABAQUS /Standard”手冊(cè)中給出了可用于任何特定單元的非均勻分布載荷類(lèi)型列表。
）
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ù)載類(lèi)型，在元素或元素表面加載積分點(diǎn)數(shù)。
LAYER——層數(shù)
KSPT——當(dāng)前圖層中的Section point
COORDS——積分點(diǎn)坐標(biāo)（數(shù)組）
JLTYP——DLOAD調(diào)用的載荷類(lèi)型的標(biāo)識(shí)
?SNAME——表面加載定義的表面名稱(chēng)（JLTYP = 0）。 對(duì)于身體力或基于單元的表面載荷，表面名稱(chēng)以空白形式傳遞
供參考的子程序
————————————————————————————————————?
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?————————————————————————————————————?
?（
用戶(hù)只能定義變量F 這是分布式負(fù)載的大小，單位取決于分布式負(fù)載的類(lèi)型：
FL–1 對(duì)于沿一維（梁）單元的線(xiàn)載荷，
FL–2 表面負(fù)荷（例如壓力）和
FL–3 體力（例如重力，向心力，加速度）。
）
?
轉(zhuǎn)載于:https://www.cnblogs.com/jiangkejie/p/10231775.html
總結(jié)
以上是生活随笔為你收集整理的Abaqus中施加移动车辆荷载（待整理）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇： VS Code插件之Cordova To
- 下一篇： aws相关文档