康 煌，王 舒，閆文敏，崔海林，郭香華，張慶明

（1. 兵器工業(yè)第208 研究所瞬時(shí)沖擊國(guó)防重點(diǎn)實(shí)驗(yàn)室，北京102202；2. 北京理工大學(xué)爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100081）

飛機(jī)起降階段，起落架、輪胎等結(jié)構(gòu)承受?chē)?yán)酷的沖擊載荷作用，其支撐部件將會(huì)出現(xiàn)嚴(yán)重的變形，甚至可能出現(xiàn)破壞現(xiàn)象。據(jù)有關(guān)數(shù)據(jù)統(tǒng)計(jì)，飛機(jī)輪胎爆炸是最頻發(fā)且最危險(xiǎn)的事故之一。在極端的壓力載荷作用下充滿高壓氣體的輪胎將產(chǎn)生高速破片，并沖擊機(jī)身蒙皮結(jié)構(gòu)，導(dǎo)致內(nèi)部乘員的生命安全受到威脅。但是，在飛機(jī)的結(jié)構(gòu)設(shè)計(jì)中蒙皮外表面通常不采用任何防護(hù)結(jié)構(gòu)。因此，開(kāi)展機(jī)身蒙皮結(jié)構(gòu)抗輪胎破片沖擊性能的實(shí)驗(yàn)研究顯得尤為重要。然而，全尺寸機(jī)身結(jié)構(gòu)生產(chǎn)制造成本高昂，在實(shí)驗(yàn)研究中采用縮比模型替代全尺寸機(jī)身結(jié)構(gòu)是行之有效的辦法。

許多學(xué)者已開(kāi)展了結(jié)構(gòu)沖擊的相似理論和實(shí)驗(yàn)研究。Mines 等通過(guò)實(shí)驗(yàn)手段研究不同沖擊速度、不同沖擊角度下橡膠輪胎破片對(duì)鋁合金板的沖擊行為，測(cè)試了輪胎橡膠靜態(tài)和動(dòng)態(tài)的拉伸、壓縮材料力學(xué)性能，實(shí)驗(yàn)結(jié)果表明這一沖擊行為并不會(huì)引起機(jī)身結(jié)構(gòu)的破裂而會(huì)產(chǎn)生大變形。Karagiozova 等采用靜、動(dòng)態(tài)的實(shí)驗(yàn)數(shù)據(jù)校核了橡膠輪胎破片的材料模型，對(duì)30°入射角下不同初始速度的橡膠破片沖擊鋁合金板進(jìn)行了數(shù)值模擬，發(fā)現(xiàn)所采用的數(shù)值模型可以較好地預(yù)測(cè)全尺寸原型的沖擊實(shí)驗(yàn)結(jié)果。Jia 等采用數(shù)值模擬技術(shù)分析了輪胎破片沖擊機(jī)身的前艙壁和主艙段的動(dòng)態(tài)響應(yīng)，獲取了不同材質(zhì)機(jī)身結(jié)構(gòu)的極限毀傷速度；對(duì)比分析發(fā)現(xiàn)，復(fù)合材料機(jī)身的最大接觸壓力約是鋁合金機(jī)身的3.3 倍，并且鋁合金前艙壁能夠承受100 m/s 的破片沖擊，而復(fù)合材料前艙壁則不能。楊娜娜等基于近場(chǎng)動(dòng)力學(xué)方法分析了破片沖擊作用下艦船復(fù)合材料的結(jié)構(gòu)損傷，發(fā)現(xiàn)增加少量質(zhì)量的加筋板的抗沖擊性能優(yōu)于復(fù)合材料層合板，加筋板的筋條尺寸和破片沖擊位置對(duì)加筋板的損傷具有明顯的影響。陳映秋等采用相似理論的基本原理，詳細(xì)地分析了模型實(shí)驗(yàn)結(jié)果預(yù)報(bào)實(shí)船結(jié)構(gòu)應(yīng)力的可行性，考慮了大尺度縮放的模型與實(shí)船的幾何尺度不能用同一比例表達(dá)而產(chǎn)生的模型畸變，提出了能夠修正模型畸變的有限元預(yù)測(cè)系數(shù)法，通過(guò)該方法實(shí)現(xiàn)了模型實(shí)驗(yàn)預(yù)測(cè)實(shí)船的結(jié)構(gòu)應(yīng)力。秦健等采用π 定律導(dǎo)出了水下爆炸沖擊不同材料加筋板的相似方法，實(shí)現(xiàn)了通過(guò)普通鋼加筋板模型預(yù)報(bào)在水下爆炸沖擊作用下船體鋼加筋板原型的動(dòng)態(tài)響應(yīng)。沈雁鳴等采用量綱分析的方法，推導(dǎo)了超高速撞擊模擬中的相似率表達(dá)式，提出了適用于超高速撞擊模擬的相似準(zhǔn)則，采用光滑粒子流體動(dòng)力學(xué)數(shù)值模擬技術(shù)，通過(guò)計(jì)算厚板開(kāi)坑、薄板穿透2 種超高速撞擊現(xiàn)象，驗(yàn)證了相似準(zhǔn)則的有效性。楊亞?wèn)|等開(kāi)展了鋼筋混凝土結(jié)構(gòu)內(nèi)爆炸的實(shí)驗(yàn)研究，發(fā)現(xiàn)比例模型之間的超壓峰值和壓力波形時(shí)程曲線符合相似律，證明了內(nèi)爆炸特征參量之間存在相似規(guī)律。Oshiro 等推導(dǎo)了無(wú)量綱相似理論模型，考慮了應(yīng)變率效應(yīng)和扭曲幾何引起的不相似行為，通過(guò)調(diào)整外載荷速度的方式，修正了大尺度縮放模型與原型之間的不相似行為，利用3 種典型的結(jié)構(gòu)沖擊模型驗(yàn)證了該方法的有效性。在此基礎(chǔ)上，Kong 等開(kāi)展了半封閉區(qū)域爆炸加載柵格板的相似性研究，實(shí)現(xiàn)了采用扭曲厚度的相似模型預(yù)測(cè)原型的動(dòng)態(tài)響應(yīng)。Kang 等采用量綱分析的方法建立了無(wú)量綱相似模型，考慮了應(yīng)變率效應(yīng)、表面效應(yīng)和扭曲幾何對(duì)淺埋地雷爆炸沖擊圓板、方板和V 形板相似行為的影響，發(fā)現(xiàn)采用修正的非等比例相似模型能夠準(zhǔn)確地預(yù)測(cè)原型的結(jié)構(gòu)響應(yīng)，可以有效地指導(dǎo)裝甲車(chē)抗沖擊性能模型實(shí)驗(yàn)設(shè)計(jì)。Mazzariol 等探究了動(dòng)載荷作用下不同材料縮比模型與原型結(jié)構(gòu)之間的相似關(guān)系，考慮了尺度效應(yīng)、材料屬性對(duì)相似行為的影響，通過(guò)理論分析認(rèn)為結(jié)構(gòu)材料的密度是影響相似行為的關(guān)鍵因素，建立了一個(gè)能夠有效預(yù)測(cè)原型動(dòng)態(tài)響應(yīng)的相似模型，通過(guò)模擬在沖擊載荷作用下梁、板結(jié)構(gòu)的動(dòng)態(tài)力學(xué)行為，驗(yàn)證了該相似方法的準(zhǔn)確性。徐坤利用質(zhì)量、速度及應(yīng)力將其他力學(xué)參量無(wú)量綱化，推導(dǎo)了縮比結(jié)構(gòu)的相似方程，采用修正速度的方法，校正了應(yīng)變率效應(yīng)引起的縮比結(jié)構(gòu)不相似行為，通過(guò)反彈道Taylor 實(shí)驗(yàn)和自由梁橫向撞擊實(shí)驗(yàn)驗(yàn)證了該方法的有效性。蘇子星等研究了在沖擊載荷作用下材料的應(yīng)變率效應(yīng)對(duì)相似行為的影響，引入Cowper-Symonds 方程推導(dǎo)了修正外載系數(shù)的相似方法，并采用45 鋼和2024 鋁合金驗(yàn)證了該修正方法的精確性。劉源等基于侵徹阻力模型得到侵徹貫穿混凝土靶板實(shí)驗(yàn)的剛體過(guò)載相似準(zhǔn)則，提出了非等比例縮比彈體的設(shè)計(jì)方法，并對(duì)縮比模型進(jìn)行了數(shù)值模擬驗(yàn)證，發(fā)現(xiàn)非等比例縮比彈體過(guò)載的脈寬、幅值均可實(shí)現(xiàn)與原型一致的加載條件。陳材等采用量綱分析的方法分析了圓柱形彈藥空氣中爆炸的相似關(guān)系，得到了模型與原型峰值壓力相似的基本條件，通過(guò)實(shí)驗(yàn)和數(shù)值模擬驗(yàn)證了縮比模型同彈藥原型的峰值壓力誤差在允許的范圍內(nèi)，彈藥的縮比模型能夠應(yīng)用于爆炸損傷實(shí)驗(yàn)。綜上，眾多的學(xué)者采用相似理論研究了沖擊載荷作用下結(jié)構(gòu)沖擊模型的縮比規(guī)律，但并未針對(duì)輪胎破片沖擊機(jī)身結(jié)構(gòu)的相似規(guī)律進(jìn)行深入研究。因此，開(kāi)展橡膠輪胎破片沖擊機(jī)身結(jié)構(gòu)的相似模型研究是有必要的，并且具有可行性，能夠節(jié)約實(shí)驗(yàn)成本，降低不必要的資源耗費(fèi)，對(duì)指導(dǎo)工程設(shè)計(jì)具有重要的實(shí)際意義。

本文中，采用數(shù)值模擬的手段優(yōu)化機(jī)身材料及結(jié)構(gòu)尺寸，分析橡膠輪胎破片的沖擊角度和著靶姿態(tài)對(duì)機(jī)身抗沖擊性能的影響，選擇最優(yōu)的機(jī)身結(jié)構(gòu)和最嚴(yán)苛的沖擊角度作為相似原型。采用量綱分析的方法，并基于線性應(yīng)變率函數(shù)，建立無(wú)量綱相似關(guān)系，通過(guò)改變模型破片的加載速度，修正應(yīng)變率效應(yīng)和扭曲幾何造成的非等比例機(jī)身模型與機(jī)身原型的不相似行為，結(jié)合數(shù)值模擬的方法驗(yàn)證縮比模型設(shè)計(jì)方法的合理性。

1 機(jī)身結(jié)構(gòu)優(yōu)化設(shè)計(jì)

1.1 材料屬性

1.1.1 機(jī)身的材料參數(shù)

對(duì)機(jī)身的金屬材料采用Johnson-Cook (J-C)模型，將應(yīng)力、應(yīng)變率和溫度等影響因素表示為乘積關(guān)系，通過(guò)對(duì)材料參數(shù)的合理取值，能較好地描述金屬材料的動(dòng)態(tài)力學(xué)行為，相關(guān)表達(dá)式為：

本文中對(duì)鋁合金和鈦合金均采用J-C 材料模型，選用文獻(xiàn)[3-4]中已測(cè)定的材料力學(xué)性能參數(shù)，見(jiàn)表1，表中ρ 為材料初始密度，為彈性模量，為泊松比。

表1 機(jī)身材料參數(shù)[3-4]Table 1 Material parameters of fuselage[3-4]

1.1.2 輪胎的材料參數(shù)

對(duì)輪胎破片采用超彈性體本構(gòu)模型，該模型基于連續(xù)介質(zhì)力學(xué)的唯象理論模型，以材料的宏觀測(cè)試特性為依據(jù)，建立可以描述變形性質(zhì)的數(shù)學(xué)關(guān)系：

1.2 有限元模型的檢驗(yàn)

本文中，建立了長(zhǎng)×寬×厚為60 mm×30 mm×15 mm 的輪胎破片有限元模型和長(zhǎng)×寬×厚為300 mm×300 mm×1.6 mm 的鋁合金靶板有限元模型，有限元模型中靶板的固定方式與實(shí)驗(yàn)中的保持一致，用螺栓將方形剛性框架固定在鋁合金靶板的四周，其實(shí)際受載面的長(zhǎng)×寬×厚為260 mm×260 mm×1.6 mm；通過(guò)LS-DYNA 動(dòng)力學(xué)模擬軟件，對(duì)初始沖擊速度為135 m/s 的輪胎破片分別以30°和80°的角度沖擊鋁合金板的工況進(jìn)行了數(shù)值模擬。

通過(guò)對(duì)比實(shí)驗(yàn)數(shù)據(jù)，具體為：靶板中心的最終撓度、應(yīng)變時(shí)間歷程曲線、變形輪廓和橡膠破片的撞擊變形形態(tài)等，從而校驗(yàn)了數(shù)值模型的有效性。圖1 中對(duì)比了30°沖擊角度β 下不同時(shí)刻橡膠輪胎破片變形形貌的高速攝影和本文的數(shù)值計(jì)算結(jié)果，圖中沖擊角度β 為速度矢量與靶板的夾角。從圖1 可以看出，橡膠作為軟體侵徹物在沖擊過(guò)程中出現(xiàn)明顯的壓彎和翻滾現(xiàn)象，從沖擊點(diǎn)開(kāi)始一直貼附于靶板運(yùn)動(dòng)，直到彈飛出靶板的邊緣，這表明有限元模型可以較好地模擬橡膠的沖擊運(yùn)動(dòng)過(guò)程，與實(shí)驗(yàn)觀測(cè)到的結(jié)果相吻合。

圖1 高速攝影[2]和數(shù)值計(jì)算分別得到的橡膠輪胎破片以30°沖擊角度和135 m/s 的沖擊速度撞擊鋁合金靶板后不同時(shí)刻的變形形貌Fig. 1 Deformation morphologies of the rubber tire fragment with the initial impact velocity of 135 m/s at the impact angle of 30°after its impacting on an aluminum alloy target at different times obtained by high-speed photography[2] and simulation

圖2 中對(duì)比了實(shí)驗(yàn)和本文模擬的在輪胎破片以30°和80°的角度沖擊下鋁合金靶板的空間變形輪廓，圖3 中對(duì)比了實(shí)驗(yàn)和本文模擬的在輪胎破片以30°的角度沖擊下鋁合金靶板中心點(diǎn)的應(yīng)變時(shí)間歷程曲線?？梢钥闯觯瑪?shù)值模擬得到的受沖擊靶板變形輪廓基本與實(shí)驗(yàn)結(jié)果相同，中心監(jiān)測(cè)點(diǎn)的應(yīng)變時(shí)間歷程曲線和實(shí)驗(yàn)結(jié)果基本一致。

圖2 不同沖擊角度下靶板的空間變形輪廓Fig. 2 Spatial deformation profiles of target plates at different impact angles

圖3 中心監(jiān)測(cè)點(diǎn)的應(yīng)變時(shí)間歷程曲線Fig. 3 Time history of strain at the central point

表2 對(duì)比了靶板中心點(diǎn)的最終撓度的實(shí)驗(yàn)結(jié)果和模擬結(jié)果，2 個(gè)不同角度沖擊作用下的靶板分別表現(xiàn)出彈性和塑性主導(dǎo)的變形行為，模擬值與實(shí)驗(yàn)值的相對(duì)偏差分別為18.6%和17.4%。這也進(jìn)一步說(shuō)明，本文中建立的有限元模型能夠有效地模擬橡膠沖擊金屬靶板的物理過(guò)程。

表2 靶板中心點(diǎn)最終撓度的實(shí)驗(yàn)結(jié)果[2]與模擬結(jié)果的對(duì)比Table 2 Comparison of the residual deformations at the centers of the target plates between experiment[2] and simulation

1.3 破片偏航姿態(tài)的影響

為了研究破片偏航姿態(tài)對(duì)靶板沖擊響應(yīng)的影響，即破片的沖擊角度和著靶角度的共同作用，著靶角度為輪胎破片中心軸線與靶板的夾角 θ 。當(dāng) θ ≤90， β=θ ；當(dāng) 9 0＜θ＜180， β =θ-90，圖4 為破片以不同姿態(tài)沖擊靶板的空間位置示意圖。因此，設(shè)置沖擊角度 β 分別為30°、60°和90°；破片的著靶角度θ 為30°、120°，60°、150°、90°和180°，沖擊角度和著靶角度兩兩相互對(duì)應(yīng)。采用數(shù)值軟件LS-DYNA 模擬分析相同沖擊速度下不同著靶角度對(duì)靶板沖擊響應(yīng)的影響，具體模擬參數(shù)及結(jié)果如表3 所示。

圖4 不同姿態(tài)的破片沖擊靶板的空間位置Fig. 4 Spatial positions of the fragments with different attitudes impacting into target plates

表3 不同沖擊條件下方形鋁靶板中心的響應(yīng)參數(shù)Table 3 Response parameters at the centers of square aluminum plates under different impact conditions

表3 展示了不同沖擊條件下方形鋁靶板中心的響應(yīng)參數(shù)，包括中心瞬時(shí)最大撓度、中心最終撓度和脈寬，脈寬為撓曲線的第1 個(gè)波谷到第1 個(gè)波峰所對(duì)應(yīng)的2 個(gè)時(shí)刻的時(shí)間差。圖5 為不同著靶角度下靶板的最終變形輪廓曲線，圖6 為不同著靶角度下靶板中心點(diǎn)的撓度時(shí)間歷程曲線。由表3 和圖5～6 可知：靶板中心撓度（代指瞬時(shí)最大撓度和最終撓度）的幅值大小主要取決于沖擊角度，而破片的偏航姿態(tài)（著靶角度）對(duì)靶板中心撓度幅值的影響次之；靶板撓曲線的脈寬隨著靶角度的增大而增大；在相同沖擊角度的情況下，鈍角著靶比銳角著靶對(duì)靶板造成的沖擊作用更嚴(yán)重，表現(xiàn)為鈍角著靶下靶板中心撓度更大；然而，在30°沖擊角度下，靶板中心撓度隨著靶姿態(tài)的變化出現(xiàn)異常現(xiàn)象，即120°著靶下靶板中心的瞬時(shí)最大撓度比30°著靶下的小，而靶板中心的最終撓度更大。對(duì)于著靶角為銳角的情況，靶板的最終變形輪廓明顯不對(duì)稱(chēng)。

圖5 靶板的最終變形輪廓曲線Fig. 5 Residual deformation profiles of target plates

圖6 靶板中心撓度的時(shí)間歷程曲線Fig. 6 Deflection time history at the center points of target plates

針對(duì)上述發(fā)現(xiàn)的現(xiàn)象及規(guī)律，分析認(rèn)為：沖擊角度決定了靶板豎直方向（向）所受的沖擊速度分量，著靶角度決定了靶板中心區(qū)域的瞬時(shí)受載面積，以及破片初次沖擊靶板后的運(yùn)動(dòng)行為，即銳角著靶姿態(tài)出現(xiàn)反向回彈，鈍角著靶姿態(tài)出現(xiàn)沿靶板的橫向移動(dòng)，并可能發(fā)生二次碰撞。對(duì)于銳角著靶姿態(tài)下靶板輪廓的不對(duì)稱(chēng)現(xiàn)象，是由于銳角著靶所產(chǎn)生的橫向沖擊速度分量和較小的瞬時(shí)受載面積，二者共同作用造成了靶板變形輪廓的不對(duì)稱(chēng)。對(duì)于30°沖擊角度下靶板撓度的異?，F(xiàn)象，是由于靶板中心的瞬時(shí)最大撓度強(qiáng)烈依賴于靶板變形的局部化程度，而最終撓度取決于靶板的受載面積、塑性變形程度和回彈響應(yīng)；這就解釋了表3 中30°沖擊角度的2 種不同著靶姿態(tài)下出現(xiàn)的異?，F(xiàn)象，即在相同的瞬時(shí)沖擊速度下，銳角著靶的局部化現(xiàn)象導(dǎo)致靶板中心的瞬時(shí)撓度更大，而局部化區(qū)域以外并未產(chǎn)生足夠的塑性變形，因此銳角著靶姿態(tài)下靶板的回彈變形幅值更大（回彈量=瞬時(shí)最大變形-最終變形），而靶板中心的最終撓度更小。

1.4 機(jī)身材料及厚度的影響

機(jī)身結(jié)構(gòu)設(shè)計(jì)需要考慮整機(jī)質(zhì)量以及耐熱、強(qiáng)度等因素，通常情況下，機(jī)身不會(huì)配備質(zhì)量過(guò)大的裝甲防護(hù)結(jié)構(gòu)，這就限制了機(jī)身的抗沖擊性能。在輕量化的要求下，機(jī)身抵御外來(lái)破片沖擊的能力仍需提升。建立真實(shí)的機(jī)身結(jié)構(gòu)抗輪胎破片垂直侵徹的有限元模型，如圖7 所示。模型幾何尺寸參數(shù)為：機(jī)身直徑，3 856 mm；寬度，1 200 mm；內(nèi)部加筋間距，500 mm；主肋骨間距，600 mm；輪胎破片尺寸，400 mm×100 mm×28 mm。

圖7 輪胎破片正侵徹機(jī)身結(jié)構(gòu)的有限元模型Fig. 7 The finite element model for a tire fragment impacting the fuselage

機(jī)身結(jié)構(gòu)相對(duì)復(fù)雜造成數(shù)值模擬的結(jié)構(gòu)響應(yīng)時(shí)間過(guò)長(zhǎng)，而機(jī)身結(jié)構(gòu)產(chǎn)生的最大瞬時(shí)變形撓度是造成乘員及結(jié)構(gòu)損傷的重要因素，因此將其作為抗沖擊性能的衡量指標(biāo)，將靶板的最大瞬時(shí)撓度簡(jiǎn)稱(chēng)為靶板撓度。選用鋁合金和鈦合金2 種材料，常用的機(jī)身厚度為3～6 mm，鋁合金靶板厚度設(shè)置為4.0、5.0 和6.0 mm，鈦合金靶板厚度設(shè)置為3.0、3.5 和4.0 mm。選用最嚴(yán)苛的沖擊條件，破片沖擊速度為150 m/s，著靶角度為180°，沖擊角度為90°。機(jī)身幾何結(jié)構(gòu)及響應(yīng)參數(shù)如表4 所示，機(jī)身中心點(diǎn)的瞬時(shí)最大撓度的時(shí)間歷程曲線如圖8 所示。

表4 機(jī)身的幾何結(jié)構(gòu)及響應(yīng)參數(shù)Table 4 Geometric structure and response parameters of the fuselage

表4 和圖8 對(duì)比了不同材料和結(jié)構(gòu)厚度對(duì)靶板沖擊響應(yīng)的影響，可以看出：3.5 mm 厚鈦合金靶板、4.0 mm 厚鈦合金靶板和6.0 mm 厚鋁合金靶板三者中心的瞬時(shí)最大撓度較小，分別為7.28、6.24 和6.77 cm；而三者的質(zhì)量由大到小依次是4.0 mm 厚鈦合金靶板的質(zhì)量最大，6.0 mm厚鋁合金靶板的質(zhì)量次之，3.5mm 鈦合金靶板的質(zhì)量最小，三者占用的空間體積由大到小依次是6.0 mm 厚鋁合金靶板占用的空間體積最大，4.0 mm厚鈦合金靶板占用的空間體積次之，3.5 mm 厚鈦合金靶板占用的空間體積最??；6.0 mm 厚鋁合金靶板的變形撓度最小，但其厚度過(guò)大造成其空間占用體積幾乎是3.5 mm 厚鈦合金靶板的2 倍；同樣，3.5 mm 厚鈦合金靶板的質(zhì)量比4.0 mm 厚鈦合金靶板的小18.9 kg，盡管后者中心的瞬時(shí)最大撓度要小1.04 cm。在輕量化要求的前提條件，3.5 mm 厚鈦合金結(jié)構(gòu)兼顧空間占用體積小、輕量化和抗沖擊性能好等多種因素，是機(jī)身結(jié)構(gòu)和材料的最佳選擇之一。

圖8 機(jī)身中心點(diǎn)的瞬時(shí)最大撓度的時(shí)間歷程曲線Fig. 8 Time histories of the maximum instantaneous deflections at the centers of the fuselages

圖9 對(duì)比了靶板中心的瞬時(shí)最大撓度時(shí)刻機(jī)身結(jié)構(gòu)的等效應(yīng)力和有效塑性應(yīng)變?cè)茍D，可以看出，鈦合金靶板中心區(qū)域產(chǎn)生的平均應(yīng)力約為540 MPa，約為鋁合金靶板中心區(qū)域平均應(yīng)力的2 倍，并且鈦合金靶板的應(yīng)力分布區(qū)域更廣。鋁合金靶板中心區(qū)域的平均塑性應(yīng)變約為0.035，約為鈦合金靶板中心區(qū)域平均塑性應(yīng)變的2 倍，鋁合金靶板的有效塑性應(yīng)變分布區(qū)域更大。這表明，兩者抵抗輪胎破片沖擊的機(jī)制不同，鈦合金靶板形成全局的沖擊應(yīng)力，而鋁合金靶板以局部變形吸能的方式耗散沖擊動(dòng)能。

圖9 瞬時(shí)最大撓度時(shí)刻靶板中心的等效應(yīng)力和有效塑性應(yīng)變?cè)茍DFig. 9 Equivalent stress and plastic strain of the fuselage structures at the instant of the maximum transient deformation

2 無(wú)量綱修正關(guān)系

金屬結(jié)構(gòu)受到強(qiáng)沖擊載荷作用會(huì)產(chǎn)生明顯的應(yīng)變率效應(yīng)，而大尺度縮比后應(yīng)變率效應(yīng)是導(dǎo)致幾何不相似的關(guān)鍵因素，許多研究結(jié)論證實(shí)了這一觀點(diǎn)。因此，Oshiro 等開(kāi)展了結(jié)構(gòu)沖擊的相似研究，嚴(yán)格地證明了應(yīng)變率效應(yīng)的影響，并提出了合理的修正方法。在此基礎(chǔ)上，本文中采用該修正方法分析應(yīng)變率效應(yīng)對(duì)輪胎橡膠破片沖擊機(jī)身結(jié)構(gòu)的相似行為的影響，并考慮大尺度縮放扭曲厚度模型的相似行為，建立相應(yīng)的修正關(guān)系，實(shí)現(xiàn)扭曲厚度的非等比例模型預(yù)測(cè)原型響應(yīng)。

2.1 應(yīng)變率效應(yīng)的修正關(guān)系

2.2 扭曲厚度的修正關(guān)系

薄壁結(jié)構(gòu)的大尺度縮放會(huì)導(dǎo)致其厚度方向的尺寸過(guò)小，引起材料的剛度和強(qiáng)度變化，無(wú)法制造出符合材料性能要求的靶板。因此，提出扭曲厚度的縮放模型，以便滿足實(shí)際材料狀態(tài)的需求。厚度的變化必然會(huì)加劇相似模型的失真程度，為解決厚度畸變導(dǎo)致的不相似行為，提出扭曲比例系數(shù) α來(lái)表征相似模型厚度方向上的單獨(dú)變化，表示扭曲模型厚度的幾何尺寸。扭曲厚度模型與等比例相似模型的幾何比例因子：

圖10 扭曲厚度模型的修正函數(shù) f2 與厚度變化因子 α x/α 的理想關(guān)系曲線[10]Fig. 10 The ideal curve of function f2 and the thickness factor α x/α [10]

3 驗(yàn)證相似模型

本文中提出一種能夠快速獲取工程縮比模型的方法。采用有限元軟件LS-DYNA，分析應(yīng)變率效應(yīng)和扭曲厚度對(duì)相似行為的影響，并通過(guò)數(shù)值模擬驗(yàn)證了縮比模型與機(jī)身結(jié)構(gòu)原型的相似性，提出了快速獲取縮比模型實(shí)驗(yàn)的設(shè)計(jì)方法。工程設(shè)計(jì)過(guò)程中，首先，對(duì)比選擇最嚴(yán)酷的沖擊條件；然后，以此為基礎(chǔ)綜合分析多種限制因素并優(yōu)選出最佳的結(jié)構(gòu)模型。接著，選用最佳的結(jié)構(gòu)和最嚴(yán)酷的沖擊條件作為原型工況，設(shè)計(jì)等比例模型和幾何尺寸扭曲的非等比例模型，采用相似修正方程獲取對(duì)應(yīng)的修正參數(shù)。最后，基于數(shù)值模擬驗(yàn)證修正參數(shù)的精確性，設(shè)計(jì)修正后的模型作為實(shí)驗(yàn)用的縮比模型，相似模型設(shè)計(jì)流程如圖11 所示。

圖11 相似模型設(shè)計(jì)流程Fig. 11 Design process for the scaled-down model

3.1 應(yīng)變率效應(yīng)的影響

在第2 節(jié)的基礎(chǔ)上，選用輪胎破片對(duì)原型機(jī)身沖擊實(shí)驗(yàn)的具體參數(shù)為：橡膠輪胎破片沖擊速度150 m/s 條件下，沖擊角度90°、著靶角度為180°沖擊3.5 mm 厚鈦合金靶板結(jié)構(gòu)。機(jī)身機(jī)構(gòu)原型與等比例模型的比例因子分別為1、1/10、1/8、1/5 和1/2；采用數(shù)值軟件LS-DYNA 模擬150 m/s 輪胎破片沖擊未修正的等比例縮放模型，獲1 取相應(yīng)模型機(jī)身結(jié)構(gòu)中心的應(yīng)變率分別為392、3 856、3 185、1 963 和794 s，代入式(3)便可獲取相應(yīng)的修正速度，如表5 所示。

表5 對(duì)比模型與原型的瞬時(shí)最大撓度Table 5 Comparison of the maximum displacement between the model and the prototype

其中，應(yīng)變率近似定義為受沖擊靶板中心位置處有效塑性應(yīng)變隨時(shí)間的變化率，通過(guò)模擬數(shù)值模型中的監(jiān)測(cè)點(diǎn)獲取靶板中心位置處的應(yīng)變-時(shí)間曲線，并通過(guò)公式計(jì)算獲取相應(yīng)的應(yīng)變率，即 ε ˙=(ε-ε)/(-) 。未修正的完全相似模型的中心變形時(shí)間歷程曲線如圖12所示。

由圖12 和表5 可以看出，未修正的完全相似模型與原型結(jié)果的偏差較小，可以較好地預(yù)測(cè)原型的結(jié)構(gòu)變形。這說(shuō)明應(yīng)變率效應(yīng)對(duì)150 m/s輪胎破片沖擊機(jī)身結(jié)構(gòu)的相似行為幾乎沒(méi)有影響，未修正的幾何完全相似模型就可以用于預(yù)測(cè)原型結(jié)構(gòu)的動(dòng)態(tài)響應(yīng)。

圖12 對(duì)比未修正的等比例模型和原型的中心變形的時(shí)間歷程曲線Fig. 12 Comparison of deformation-time curves between the uncorrected scaled models and the prototype

3.2 扭曲厚度的影響

在3.1 節(jié)中，大尺度縮放模型的靶板厚度分別為0.350、0.438 和0.700 mm。由于制造技術(shù)的限制，難以生產(chǎn)出符合材料強(qiáng)度要求的大尺度縮放模型實(shí)物。因此，引入扭曲厚度的非等比例模型處理這一問(wèn)題，仍將考慮應(yīng)變率效應(yīng)微小變化對(duì)非等比例修正模型的影響。在本節(jié)中，采用等比例縮放模型的比例因子 α =1/8，模型的機(jī)身結(jié)構(gòu)厚度=0.438 mm。

基于2.2 節(jié)中提出的扭曲模型的修正方法，為獲取修正函數(shù)的指數(shù)，設(shè)定2 個(gè)不同的扭曲模型系數(shù)α=2.0 和α=3.5 ，對(duì)應(yīng)的機(jī)身厚度為0.87 和1.52 mm。設(shè)定4 組不同的指數(shù)，分別為1.00、1.05、1.12、1.15，通過(guò)數(shù)字計(jì)算獲取2 個(gè)扭曲模型靶板的變形撓度 δ ；當(dāng) δ≈δ，便可獲取對(duì)應(yīng)的扭曲函數(shù)的指數(shù)，不同指數(shù)對(duì)應(yīng)的沖擊速度和撓度如表6 所示。最后，將求解得到的指數(shù)代入=(α/α)，便可計(jì)算扭曲厚度非等比例模型的修正速度。

表6 不同指數(shù)對(duì)應(yīng)的沖擊速度和撓度Table 6 Impact velocities and central deflections in relation to different exponents

由表6 可以看出，隨著指數(shù)的增大，2 個(gè)扭曲模型撓度的相對(duì)偏差減?。划?dāng)=1.15 時(shí)，2 個(gè)扭曲模型撓度的偏差僅為3%；當(dāng)＞1.15 時(shí)，存在某一極值點(diǎn)使兩者的撓度偏差接近零；而進(jìn)一步增大指數(shù)，兩者的撓度偏差會(huì)逐漸增大，這是由修正函數(shù)的性質(zhì)所決定的。由模擬結(jié)果發(fā)現(xiàn)，當(dāng)＞1.15 時(shí)，機(jī)身結(jié)構(gòu)會(huì)出現(xiàn)明顯的破壞，3%的撓度相對(duì)偏差足以保證非等比例相似模型的修正精度。因此，當(dāng)=1.15 時(shí)，獲得扭曲厚度模型的修正方程為：

為了驗(yàn)證扭曲厚度模型的有效性，假定扭曲因子分別為2.29、2.74 和3.43，將求解的參數(shù)和分別代入式(3)和(8)，計(jì)算獲得扭曲厚度模型的修正速度，將參數(shù)代入數(shù)值模擬軟件DYNA 獲取未修正模型和修正結(jié)果。表7 對(duì)比原型與模型的靶板瞬時(shí)最大撓度，可以看出，修正方程能夠解決扭曲厚度所導(dǎo)致的偏差，修正后的3 個(gè)不同厚度的扭曲模型的預(yù)測(cè)值最大偏差不超過(guò)5.1%。

表7 對(duì)比原型與模型的靶板瞬時(shí)最大撓度Table 7 Comparison of the maximum deflection between the prototype and models

圖13 對(duì)比了扭曲厚度的非等比例模型與原型的撓度時(shí)間歷程曲線，可以看出，修正方法可以彌補(bǔ)因板厚扭曲而導(dǎo)致的撓度峰值和靶板響應(yīng)模型不一致的問(wèn)題。然而，在時(shí)間域上，3 個(gè)修正的扭曲厚度模型的撓曲線的分布與原型的相差較大。對(duì)比發(fā)現(xiàn)，原型撓度曲線的第1 個(gè)脈寬遠(yuǎn)大于其他修正模型的，未修正扭曲模型的撓度曲線的脈寬具有較好的一致性；此外，隨著扭曲板厚的增大，修正模型的撓度脈寬逐漸減小。分析發(fā)現(xiàn)，在時(shí)間域上，扭曲模型產(chǎn)生的行為不一致現(xiàn)象，是由橡膠破片的材料性質(zhì)所決定的；這是因?yàn)樵撔拚椒ㄍㄟ^(guò)調(diào)整載荷速度從而實(shí)現(xiàn)扭曲厚度模型行為的相似性，而橡膠破片沖擊靶板會(huì)出現(xiàn)明顯的“彈飛”現(xiàn)象，導(dǎo)致修正速度越高的扭曲模型，橡膠破片沖擊作用靶板的時(shí)間越短。因此，橡膠的材料屬性是引起修正扭曲模型響應(yīng)曲線出現(xiàn)不相似行為的關(guān)鍵因素。

圖13 對(duì)比非等比例模型與原型的撓度時(shí)間歷程曲線Fig. 13 Comparison of deformation-time curves for the center points of the plates

圖14 對(duì)比了非等比例模型與原型靶板的最大變形時(shí)刻的輪廓圖。可以看出，不同扭曲厚度模型的靶板撓度輪廓變形與原型結(jié)果基本一致。在空間域上，靶板的撓度與原型的變形曲線基本一致。因此，采用該修正方法仍然能夠準(zhǔn)確地預(yù)測(cè)原型板的最大撓度，能夠有效地指導(dǎo)工程設(shè)計(jì)。

圖14 非等比例模型與原型靶板的最大變形時(shí)刻的輪廓圖Fig. 14 The maximum deformation profiles of the plates for the incomplete scaling model and prototype

4 結(jié) 束 語(yǔ)

飛機(jī)起降階段產(chǎn)生的不可控因素會(huì)使起落架和輪胎承受巨大的壓力，導(dǎo)致飛機(jī)輪胎發(fā)生爆炸，產(chǎn)生四散的破片沖擊飛機(jī)四周及威脅艙內(nèi)乘員的安全。因此，為了解決全尺寸原型抗沖擊性能實(shí)驗(yàn)耗費(fèi)巨大的問(wèn)題，開(kāi)展輪胎破片沖擊機(jī)身結(jié)構(gòu)的相似性研究，建立相似模型設(shè)計(jì)方法，以便實(shí)現(xiàn)快速獲取相似模型的設(shè)計(jì)方案。

(1)基于測(cè)試實(shí)驗(yàn)數(shù)據(jù)，驗(yàn)證了橡膠輪胎破片沖擊鋁合金靶板有限元模型的有效性。為建立受沖擊機(jī)身的原型，綜合分析了多種因素，結(jié)果表明：150 m/s 輪胎破片沖擊速度下，90°沖擊角度和180°著靶角度為輪胎破片最嚴(yán)苛的沖擊條件；3.5 mm 厚鈦合金板是機(jī)身的最佳選擇。

(2)在沖擊角度為30°的情況下，銳角著靶比鈍角著靶的受沖擊靶板的瞬時(shí)最大撓度更大，而最終撓度更小。因?yàn)闆_擊角度較小導(dǎo)致受載靶板的局部化塑性變形程度更高，所以銳角著靶的瞬時(shí)最大撓度更大；然而，局部化區(qū)域以外受彈性變形主導(dǎo)，導(dǎo)致靶板的回彈幅值更大，所以最終撓度更小。

(3)基于量綱分析方法，建立了應(yīng)變率效應(yīng)和扭曲幾何的無(wú)量綱修正公式。為驗(yàn)證該相似方法，建立了不同縮放系數(shù)的等比例和非等比例相似模型。未修正模型的結(jié)果表明，破片沖擊機(jī)身結(jié)構(gòu)的應(yīng)變率效應(yīng)不會(huì)影響等比例縮放模型的相似性，而扭曲厚度會(huì)導(dǎo)致非等比例縮放模型產(chǎn)生嚴(yán)重的不相似行為。修正模型的結(jié)果表明：在時(shí)間域上，修正模型的撓度曲線的脈寬明顯小于原型的，這是由橡膠材料的特殊性質(zhì)所決定的；然而，在空間域上，修正模型的變形輪廓曲線和瞬時(shí)最大撓度的最大偏差不超過(guò)5.1%。這表明采用本文中所提出的修正方法，可以快速建立輪胎破片沖擊機(jī)身結(jié)構(gòu)的非等比例相似模型，實(shí)現(xiàn)準(zhǔn)確預(yù)測(cè)原型機(jī)身結(jié)構(gòu)的沖擊響應(yīng)。