一、寫(xiě)在文前

滾動(dòng)輪胎模態(tài)仿真實(shí)際上是在輪荷加載的基礎(chǔ)之上的重啟動(dòng)分析。輪胎在穩(wěn)態(tài)滾動(dòng)過(guò)程中,會(huì)受到預(yù)加載荷、慣性力以及輪胎和地面的摩擦力的影響,這些力會(huì)對(duì)整個(gè)系統(tǒng)的剛度矩陣和阻尼矩陣產(chǎn)生影響,導(dǎo)致非對(duì)稱(chēng)性。故不能采用常規(guī)方法對(duì)動(dòng)力學(xué)方程進(jìn)行解耦,必須用復(fù)模態(tài)來(lái)解耦，所以滾動(dòng)輪胎的模態(tài)仿真其實(shí)是復(fù)模態(tài)的的提取。

在abaqus的穩(wěn)態(tài)滾動(dòng)中,輪胎實(shí)際上是不滾動(dòng)的,只是內(nèi)部材料的流動(dòng)（歐拉-拉格朗日法）。在提取復(fù)模態(tài)之前,必須保證輪胎滾動(dòng)的轉(zhuǎn)速和線速度相匹配，故需先進(jìn)行roll tire仿真調(diào)試：

roll tire計(jì)算中，先給定輪胎線速度、轉(zhuǎn)動(dòng)角速度然后提取輪胎輪心的力矩M,當(dāng)輪胎穩(wěn)態(tài)滾動(dòng)的時(shí)候, 輪胎輪心的力矩M應(yīng)該為0。在實(shí)際操作中,需要不斷的調(diào)節(jié)定義的ω值,使最終繞Y向的力矩M在[-10,10]之內(nèi)。

進(jìn)行roll tire計(jì)算時(shí)，首先進(jìn)行step1二維輪胎充氣仿真計(jì)算，然后進(jìn)行step2rev旋轉(zhuǎn)3D輪胎生成及輪荷加載計(jì)算，在此基礎(chǔ)上進(jìn)行Free roll計(jì)算Inp文件的編寫(xiě)，進(jìn)行計(jì)算，查看輪胎輪心的力矩M判斷輪胎是否處于穩(wěn)態(tài)滾動(dòng)狀態(tài)。下圖為step1.inp以及step2rev.inp運(yùn)行結(jié)果圖：





二、rolltire計(jì)算

step3roll_tire.inp文件編寫(xiě)如下:

*HEADING（重啟動(dòng)分析） Step3: roll_tire from step2rev.inp*RESTART,READ*STEP,INC=500,NLGEOM=YES,UNSYMM=YES（幾何非線性打開(kāi)；非對(duì)稱(chēng)打開(kāi)）4: roll_tire at 60km/h*STEADY STATE TRANSPORT,LONG TERM,INERTIA=YES（動(dòng)態(tài)分析，不考慮粘塑性影響，慣性打開(kāi)） 0.5, 1.0, 1E-6, 1.0*TRANSPORT VELOCITY（初始角速度=V/r） TIRE-RIM,57.0776*MOTION,TYPE=VELOCITY,TRANSLATION（線速度V） NODE_TIRE,1, （方向?yàn)?）,16666.67（60公里每小時(shí)化為毫米/S）*CHANGE FRICTION,INTERACTION=SRIGID（改變摩擦系數(shù)）*FRICTION,SLIP=0.02 1.0*NODE PRINT,FREQ=1*EL PRINT,FREQ=1*OUTPUT,HISTORY,FREQ=1*NODE PRINT,FREQUENCY=999,NSET=TIRE-RIM RF, （輸出輪心的力矩）*END STEP

abaqus job=step3roll_tire  oldjob=step2rev int cpus=4

點(diǎn)擊回車(chē)運(yùn)行

?

計(jì)算完成后可出現(xiàn)以下幾個(gè)文件：

?

其中： 提交計(jì)算文件的語(yǔ)法錯(cuò)誤信息在.dat文件中查看； 提交計(jì)算文件的計(jì)算步長(zhǎng)在.sta文件中查看； 所得計(jì)算結(jié)果的.odb文件可在Hyper view或Abaqus中查看；

因我們需要的是輪胎輪心的力矩M：

?

故用UE打開(kāi)step3roll_tire.dat文件查看NODE OUTPUT： 因我們需要的是繞Y向的力矩M（RM2）在[-10,10]之內(nèi)，這時(shí)才算穩(wěn)態(tài)滾動(dòng)，而輸出值RM2為-9.5877E 05故需調(diào)整ω值,使最終繞Y向的力矩M在[-10,10]之內(nèi)。在此大家自己調(diào)試即可。

?

三、滾動(dòng)輪胎模態(tài)提取

在ABAQUS/STANDARD中,復(fù)模態(tài)的提取采用子空間投影法,而且在復(fù)模態(tài)提取之前,一般先進(jìn)行具有對(duì)稱(chēng)剛度矩陣的無(wú)阻尼系統(tǒng)的特征值和特征模態(tài)計(jì)算。

滾動(dòng)輪胎模態(tài)提取inp文件編寫(xiě)如下：

*HEADING
 Step4:Caculation complex mode from step3roll_tire.inp
*RESTART,READ
*STEP,PERTURBATION
5:static modes analysis
*FREQUENCY,EIGENSOLVER=LANCZOS,NORMALIZATION=DISPLACEMENT（實(shí)模態(tài)提?。?,0.1,300.,,,
*ACOUSTIC FLOW VELOCITY,ROTATION（定義流體聲腔的角速度）
 NODE_AIR,57.0776
*BOUNDARY,OP=NEW
 ROAD,1,6（路面全約束）
*END STEP
**
*STEP,NLGEOM=YES
 6:Complex modes analysis,spindle free
*COMPLEX FREQUENCY,FRICTION DAMPING=YES,PROPERTY EVALUATION,NORMALIZATION=DISPLACEMENT（復(fù)模態(tài)提?。? ,
*ACOUSTIC FLOW VELOCITY,ROTATION
 NODE_AIR,57.0776
*END STEP

abaqus job=step4rolltire_mode oldjob=step3roll_tire int cpus=4

提交計(jì)算文件的語(yǔ)法錯(cuò)誤信息在.dat文件中查看；

提交計(jì)算文件的計(jì)算步長(zhǎng)在.sta文件中查看；

所得計(jì)算結(jié)果的.odb文件可在Hyper view或Abaqus中查看；

現(xiàn)用UE打開(kāi)step4rolltire_mode.dat文件可得滾動(dòng)輪胎的實(shí)模態(tài)頻率如下：

?

復(fù)模態(tài)頻率如下：

?

（注：step4rolltire_mode.inp是基于上文step3roll_tire.inp進(jìn)行運(yùn)算的，沒(méi)有進(jìn)行ω值的調(diào)試，故上復(fù)模態(tài)頻率不作為參考依據(jù)，下只呈現(xiàn)實(shí)模態(tài)振型云圖，大家可后續(xù)進(jìn)行ω值的,使繞Y向的力矩M在[-10,10]之內(nèi)再進(jìn)行復(fù)模態(tài)頻率的提?。?/p>

在Abaqus中打開(kāi)step4rolltire_mode的ODB文件可得：

審核編輯：劉清