移動車輛沖擊荷載在ABAQUS中的模擬研究

胡興宇

(土木工程學院（重慶交通大學）， 重慶 400041)

0 引言

在模擬車輛荷載計算時，移動荷載功能是必不可少的。ABAQUS中用戶界面可以通過在結構表面通過partution命令創建移動荷載的軌跡，通過相鄰的矩形方塊組成，根據移動的速度和軌跡距離來激活相應的矩形方塊荷載，但是如果軌跡長度非常長，會導致工作相當繁瑣，分析步可能會成百上千，俗稱“荷載移動帶”。子程序可以完全解決這些問題，并且可以根據用戶需要的方式繼續模擬，大大減少了建模時間，計算結果十分精確，避免了多個分析步造成的錯誤。

1 子程序的編寫

模擬車輛荷載的作用首先需要確定作用面積和速度，然而在ABAQUS軟件中，作用面積是不能直接人為選定的，需要用戶在模擬荷載時簡單創建處荷載整個過程需要經過的軌跡，之后只需要在子程序中賦予每一個車輛的長度即可。

Abaqus已經可以通過python語言進行功能開發了，但本文采用通用數學Fortran語言進行編寫。在用戶手冊中，荷載類型主要用DLOAD和UTRACLOAD兩種，在車輛荷載模擬中，DLOAD主要用于模擬車輪豎向荷載，UTRACLOAD用于模擬車輪水平摩擦荷載。

DLOAD和UTRACLOAD編寫格式類似，本節以

DLOAD為例編寫。

DLOAD實現格式如下：

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’S CODES

RETURN

END

上述代碼中其實部分是子程序的固定模式，F代表的荷載值，KSTEP代表在第幾個分析步時被調用，KINC即是代表增量步，TIME(1)和TIME(2)分別代表當前時間和總時間，這對用戶控制變化荷載極為有用，NOEL和NPT分別代表單元號和積分點號，COORDS（I）用于激活移動荷載的方向，I=1、2、3分別代表下沿x、y、z方向與移動，用戶只用通過IF從句賦予F即可模擬出荷載移動。

2 沖擊荷載模擬

2.1 模型參數

車輛通過路面或者橋面時，遇到路面不平整將會產生跳車現象，本節根據這一現象模擬車輛后輪對橋面板的沖擊進行模擬。

模型參數:混凝土材料選用C40，橋面板長為10m，寬為4m，厚度為0.4m，鋼筋采用HRB335，分布在橋面板上層和下層，跳車車輛荷載峰值取30KN，采用文獻中公式推導曲線，作用時間為0.03s。

2.2 模擬沖擊荷載方法

由于跳車荷載時隨著時間不斷變化的，本文通過得到的沖擊荷載曲線進行公式擬合，得到相應的公式，再通過一下條件語句賦予:

IF((X_NOW.GE.ROLL_1_1).AND.(X_NOE.LE.ROLL_1_0))THEN IF(KSTEP.EQ.1)THEN F=F(t) END IF ELSE F=0.0 END IF

將由于擬合函數的自變量賦予成時間t，可以實現荷載隨著時間增長而不斷變化，再其余未賦予荷載的路徑上均為0。

由于車輛再行駛過程中沖擊時間很小，其余時間均屬于勻速狀態，作用在橋面板的荷載時不變動的，需要在REAL關鍵字賦予相應的時間TIME1、TIME2等控制荷載移動過程的時間控制，如下所示實現：

Real*6 TIME1、TIME2、TIME3、V_CAR..

V_CAR=10E3 …

IF((X_NOW.GE.ROLL_1_1).AND.(X_NOE.LE.ROLL_1_0))THEN

IF(KSTEP.EQ.1)THEN

F=0.1

ELSE IF(KSTEP.EQ.2)THEN

F=SIN(2*3.14/0.06*(TIME(2)-0.1)) ！擬合公式

ELSE IF(KSTEP.EQ.3)THEN

F=0.1

END IF ELSE F=0.0 END IF

上述可以實現車輛勻速過程中突然的跳車沖擊，但需要用戶在step中賦予三段分析步來描述不同荷載作用狀態。

3 結果分析

3.1 加速度分析

計算過程采用動力隱式模塊進行求解，三個分析步分別設置為0.1s、0.03s和0.87s，經計算得到跳車時橋面板跨中的加速度時程曲線。

根據結果得知，車輛荷載在開始到0.1s是勻速前進，到達0.1s時發生跳車現象，跨中豎向加速度增大，振動響應劇烈，0.13s之后，隨著車輛繼續行駛，加速度開始下降，振動趨勢穩定至0.5s，此時車輛荷載到達跨中，此處由于車輪駛過加速度曲線較為密集，0.5s后加速度值逐漸減小。可知子程序模擬移動沖擊荷載是有效可行的。

3.2 應力分析

從odb文件中提取得到橋面板的等效應力云圖，是可以明顯看出車輛沖擊的部位出現了明顯的壓應力，其值為0.2MPa左右。

4 小結

本文主要介紹了DLOAD子程序的編寫和運行過程，UTRACLOAD使用方法與其類似，ABAQUS平臺提供了非常人性化的界面與命令，用戶可以根據自己的需要對子程序進行編寫。通過全文的簡單介紹得到了一下幾點結論：

1）子程序的使用替代了傳統的“荷載移動帶”的方法，為設計節約了時間，同時減少了錯誤的發生，是種有效可行的方法。

2）通過Fortran語言可以解決大部分用戶的需求，本文只模擬了一對車輪荷載，如果需要模擬多對車輪荷載，需要對代碼進行補充，引入DISTANCE_ROLL參數來賦予前后車輪邊之間的距離，便可以實現多個荷載模擬。

3）本文編寫的子程序只對靜力通用和動力隱式分析有效，但對于動力顯示是無效的，由于隱式和顯示的求解方程不同，導致子程序格式有些差異。

綜上所述，ABAQUS平臺功能強大，需要用戶對編程語言較為熟悉，這樣才能將許多繁瑣的工作變得簡單。