劉 賽， 張偉貴， 肖 凱， 蘇 玲， 王 悅

(1. 中國(guó)運(yùn)載火箭技術(shù)研究院， 北京 100076； 2. 中國(guó)科學(xué)院空間應(yīng)用工程與技術(shù)中心， 北京 100094)

梁桿結(jié)構(gòu)的動(dòng)態(tài)屈曲準(zhǔn)則是進(jìn)行工程設(shè)計(jì)的必要依據(jù)，現(xiàn)有梁桿結(jié)構(gòu)的動(dòng)態(tài)屈曲準(zhǔn)則可以分為三類：一類是以分叉理論作為基礎(chǔ)的屈曲準(zhǔn)則[1-4]，僅用于某些沖擊速度較高的能列出解析方程求解的簡(jiǎn)單工況，而且一般分叉解距離失穩(wěn)解較遠(yuǎn)，不滿足工程需要；第二類是基于放大函數(shù)法的失穩(wěn)準(zhǔn)則[5-7]，但得到的臨界載荷多依賴人為主觀判斷；第三類是基于載荷-響應(yīng)曲線的B-R(Budiansky-Roth)運(yùn)動(dòng)準(zhǔn)則[8-10]，此類準(zhǔn)則適于工程應(yīng)用，但需要測(cè)量的物理量種類較多。因此亟待有一種不僅滿足梁桿結(jié)構(gòu)工程需要，而且涉及的物理量測(cè)量簡(jiǎn)便的動(dòng)態(tài)屈曲準(zhǔn)則。

本文針對(duì)文獻(xiàn)[11]的試驗(yàn)研究對(duì)象，采用文獻(xiàn)[12]中仿真精度較高的有限元模型，進(jìn)行不同重物沖擊速度下扁長(zhǎng)桿的動(dòng)態(tài)響應(yīng)分析，研究桿發(fā)生彈塑性動(dòng)態(tài)屈曲的準(zhǔn)則，以及初始條件參數(shù)(材料屈服應(yīng)力、預(yù)制弓形幅值和截面尺寸等)對(duì)臨界沖擊速度的影響。

1 扁長(zhǎng)桿沖擊屈曲的仿真分析模型

在文獻(xiàn)[11]的試驗(yàn)中，矩形截面扁長(zhǎng)桿的一端鉸支，另一端為加載端，如圖1所示。一質(zhì)量為369 g、速度為2.45 m/s的重物撞擊扁長(zhǎng)桿的加載端，重物的側(cè)向運(yùn)動(dòng)受到約束，使得扁長(zhǎng)桿加載端運(yùn)動(dòng)等效于滑動(dòng)鉸支移動(dòng)。扁長(zhǎng)桿長(zhǎng)度、截面寬度和厚度分別為193 mm、8.7 mm和0.63 mm，扁長(zhǎng)桿預(yù)制為微弓形，其定義見(jiàn)表1。材料為Ni-Cr合金鋼，密度為7.8 g/cm3，彈性模量為206 GPa，屈服強(qiáng)度為1 470 MPa。

圖1 受重物沖擊的微弓形扁長(zhǎng)桿[13]

顯式動(dòng)力學(xué)有限元仿真分析模型網(wǎng)格如圖2～圖4所示。參考文獻(xiàn)[12]對(duì)單元類型和單元尺寸的探討，采用單點(diǎn)積分8節(jié)點(diǎn)實(shí)體單元建模，扁長(zhǎng)桿的實(shí)體單元尺寸為0.16 mm。有限元模型的邊界條件參數(shù)取文獻(xiàn)[12]辨識(shí)出的改進(jìn)邊界條件參數(shù)。扁長(zhǎng)桿的材料模型為線性隨動(dòng)強(qiáng)化彈塑性模型，應(yīng)變率效應(yīng)使用Cowper-Symonds模型，參數(shù)如表2所示。其他部件均使用線彈性材料模型，參數(shù)如表3所示。仿真分析結(jié)果達(dá)到了較高的精度，如圖5和圖6所示。

表1 扁長(zhǎng)桿的預(yù)制微弓形

表2 扁長(zhǎng)桿的材料模型參數(shù)

表3 其他部件的材料模型參數(shù)

2 軸向沖擊作用下桿的動(dòng)態(tài)響應(yīng)

圖5 t=0.75 ms扁長(zhǎng)桿的屈曲變形曲線

圖6 扁長(zhǎng)桿處軸向力

圖7 重物沖擊速度為2.45 m/s時(shí)重物的速度、位移和扁長(zhǎng)桿加載端位移

Fig.7 Velocity and displacement of the weight and displacement of the loaded end under the impact velocity 2.45 m/s

圖10 重物沖擊速度為0.7 m/s時(shí)重物的速度、位移和扁長(zhǎng)桿加載端位移

Fig.10 Velocity and displacement of the weight and displacement of the loaded end under the impact velocity 0.7 m/s

圖11 重物沖擊速度為0.7 m/s時(shí)扁長(zhǎng)桿的能量和處軸向力

根據(jù)B-R運(yùn)動(dòng)準(zhǔn)則，通過(guò)多次改變仿真分析模型的沖擊速度，得到扁長(zhǎng)桿的臨界沖擊速度為1.3 m/s，臨界狀態(tài)見(jiàn)圖9，圖12和圖13，扁長(zhǎng)桿的軸向力下降，重物短時(shí)靜止，內(nèi)能下降，扁長(zhǎng)桿的加載端與重物發(fā)生明顯分離——這些現(xiàn)象同時(shí)發(fā)生。

圖12 重物沖擊速度為1.3 m/s時(shí)重物的速度、位移和扁長(zhǎng)桿加載端位移

Fig.12 Velocity and displacement of the weight and displacement of the loaded end under the impact velocity 1.3 m/s

3 軸向沖擊作用下桿的彈塑性動(dòng)態(tài)屈曲準(zhǔn)則

隨著重物沖擊速度的增加，扁長(zhǎng)桿依次出現(xiàn)彈性振動(dòng)和彈性動(dòng)態(tài)屈曲——兩種不同的動(dòng)態(tài)響應(yīng)。比較兩種動(dòng)態(tài)響應(yīng)的統(tǒng)計(jì)量特點(diǎn)表明，存在三種等價(jià)的判斷準(zhǔn)則：① 桿在軸向加載過(guò)程中，軸向力突然下降(B-R準(zhǔn)則)；② 桿在軸向加載過(guò)程中，內(nèi)能突然下降或動(dòng)能突然上升；③ 桿在軸向加載過(guò)程中，桿的加載端與重物發(fā)生明顯的分離。

圖13 重物沖擊速度為1.3 m/s時(shí)扁長(zhǎng)桿的能量和處軸向力

后兩者存在一定的內(nèi)在關(guān)聯(lián)，扁長(zhǎng)桿的內(nèi)能突然下降時(shí)，其動(dòng)能突然上升，表明扁長(zhǎng)桿橫向的屈曲速度和軸向的運(yùn)動(dòng)速度突然增大，導(dǎo)致了扁長(zhǎng)桿的加載端與重物發(fā)生明顯的分離。第③個(gè)準(zhǔn)則較第②個(gè)更適合工程應(yīng)用，可將其定義為加載分離準(zhǔn)則，其數(shù)學(xué)表達(dá)式如下(參考坐標(biāo)系見(jiàn)圖1)：

若桿的加載端與重物未發(fā)生明顯分離，則發(fā)生彈性振動(dòng)。雖然加載分離準(zhǔn)則只是通過(guò)數(shù)值試驗(yàn)總結(jié)而來(lái)，尚未有較深刻的力學(xué)理論支撐，但數(shù)值試驗(yàn)表明，加載分離準(zhǔn)則和B-R運(yùn)動(dòng)準(zhǔn)則識(shí)別的扁長(zhǎng)桿動(dòng)態(tài)屈曲的臨界沖擊速度是一致的，均為1.3 m/s。B-R運(yùn)動(dòng)準(zhǔn)則需要重物位移和軸向力兩個(gè)物理量，屬于兩類物理量，試驗(yàn)中需要分別使用位移測(cè)量系統(tǒng)和應(yīng)變測(cè)量系統(tǒng)測(cè)得，而加載分離準(zhǔn)則需要加載端位移和重物位移兩個(gè)物理量，屬于同一類物理量，試驗(yàn)中僅需使用位移測(cè)量系統(tǒng)測(cè)得。加載分離準(zhǔn)則減少了需要測(cè)量的物理量種類，應(yīng)用更加簡(jiǎn)捷。因此從判別方法的實(shí)用性和臨界沖擊速度的識(shí)別精度而言，加載分離準(zhǔn)則可以滿足工程需要。加載分離準(zhǔn)則和B-R運(yùn)動(dòng)準(zhǔn)則都不能判斷桿發(fā)生的是彈性動(dòng)態(tài)屈曲還是塑性動(dòng)態(tài)屈曲。

4 扁長(zhǎng)桿的初始條件參數(shù)對(duì)臨界沖擊速度的影響

對(duì)于工程實(shí)際問(wèn)題——如何提高一定長(zhǎng)度桿件的動(dòng)態(tài)屈曲強(qiáng)度，一般的途徑都是更換比較強(qiáng)的材料、提高加工精度或增大截面尺寸，本文從三方面著手：① 將第1節(jié)扁長(zhǎng)桿的初始條件參數(shù)(材料屈服應(yīng)力、預(yù)制弓形幅值和截面尺寸等)作為歸一化的參考值，通過(guò)大量數(shù)值計(jì)算，得到扁長(zhǎng)桿的初始條件參數(shù)歸一化值與臨界沖擊速度的關(guān)系曲線，如圖14所示。當(dāng)一個(gè)參數(shù)變化時(shí)，其余參數(shù)是不變的，仍舊為1；② 大量數(shù)值計(jì)算表明，加載分離準(zhǔn)則和B-R運(yùn)動(dòng)準(zhǔn)則識(shí)別的扁長(zhǎng)桿動(dòng)態(tài)屈曲的臨界沖擊速度是一致的，這也說(shuō)明了加載分離準(zhǔn)則的一般性和有效性；③ 圖14表明，屈服應(yīng)力僅在低于350 MPa時(shí)對(duì)臨界沖擊速度的影響顯著，因?yàn)楸忾L(zhǎng)桿此時(shí)發(fā)生塑性動(dòng)態(tài)屈曲，高于此屈服應(yīng)力時(shí)發(fā)生彈性動(dòng)態(tài)屈曲；微小的預(yù)制弓形可明顯降低扁長(zhǎng)桿的動(dòng)態(tài)屈曲強(qiáng)度；增大截面尺寸是提高扁長(zhǎng)桿動(dòng)態(tài)屈曲強(qiáng)度的較有效途徑。

圖14 扁長(zhǎng)桿的初始條件參數(shù)歸一化值與臨界沖擊速度關(guān)系

5 結(jié) 論

(1) 通過(guò)不同重物沖擊速度下扁長(zhǎng)桿動(dòng)態(tài)響應(yīng)的有限元分析，提出了桿在軸向沖擊下的動(dòng)態(tài)屈曲準(zhǔn)則——加載分離準(zhǔn)則：重物對(duì)桿的軸向加載過(guò)程中，桿的加載端與重物發(fā)生明顯的分離。

(2) 雖然目前沒(méi)有較深刻的力學(xué)理論支撐，但與B-R運(yùn)動(dòng)準(zhǔn)則相比，加載分離準(zhǔn)則同樣可以滿足工程需要，但涉及的物理量測(cè)量更為簡(jiǎn)便。

(3) 進(jìn)一步分析了扁長(zhǎng)桿的初始條件參數(shù)(材料屈服應(yīng)力、預(yù)制弓形幅值和截面尺寸等)對(duì)臨界沖擊速度的影響規(guī)律，驗(yàn)證了此準(zhǔn)則的一般性和有效性，總結(jié)了提高此種桿件動(dòng)態(tài)屈曲強(qiáng)度的有效措施：與提高扁長(zhǎng)桿的材料強(qiáng)度和加工精度相比，增大截面尺寸是提高其動(dòng)態(tài)屈曲強(qiáng)度的較有效方法。

本文的研究成果具有一般性意義，適用于其它梁桿結(jié)構(gòu)的動(dòng)態(tài)屈曲問(wèn)題，也有助于其它結(jié)構(gòu)的動(dòng)態(tài)屈曲準(zhǔn)則的發(fā)展。