韓大偉，王安穩(wěn)

（海軍工程大學(xué)理學(xué)院力學(xué)系，武漢 430033）

矩形薄板在面內(nèi)撞擊下動(dòng)力屈曲的實(shí)驗(yàn)研究

韓大偉，王安穩(wěn)

（海軍工程大學(xué)理學(xué)院力學(xué)系，武漢 430033）

利用Hopkinson壓桿系統(tǒng)對(duì)矩形薄板在面內(nèi)撞擊下的動(dòng)力屈曲進(jìn)行了實(shí)驗(yàn)研究。結(jié)果表明矩形薄板存在屈曲后強(qiáng)度，在屈曲發(fā)生后面內(nèi)載荷可以繼續(xù)增加，通過(guò)記錄屈曲發(fā)生的時(shí)間求得臨界屈曲長(zhǎng)度，并對(duì)實(shí)驗(yàn)?zāi)Ｐ瓦M(jìn)行了理論分析。對(duì)比分析結(jié)果表明，當(dāng)壓應(yīng)力波向前傳播一定距離后矩形薄板發(fā)生局部屈曲，得到的臨界屈曲長(zhǎng)度與理論計(jì)算結(jié)果吻合良好。

動(dòng)力屈曲；矩形薄板；臨界屈曲長(zhǎng)度；Hopkinson壓桿；實(shí)驗(yàn)研究

近年來(lái)，利用夾心構(gòu)件的失穩(wěn)吸能來(lái)提高結(jié)構(gòu)的抗沖擊性能的研究越來(lái)越受到國(guó)內(nèi)外學(xué)者的重視，對(duì)夾心構(gòu)件動(dòng)力失穩(wěn)機(jī)理的研究逐漸成為固體力學(xué)研究領(lǐng)域的熱點(diǎn)之一。結(jié)構(gòu)在沖擊載荷作用下的動(dòng)力屈曲有顯著不同于靜力屈曲的特性，必須考慮橫向慣性效應(yīng)的影響。文獻(xiàn)［1－4］從理論上研究了幾何完善結(jié)構(gòu)的動(dòng)力屈曲，分析結(jié)果表明，結(jié)構(gòu)在沖擊載荷作用下會(huì)發(fā)生局部屈曲，而且動(dòng)力屈曲的臨界載荷要比靜力屈曲的情況大很多。文獻(xiàn)［5］用數(shù)值方法研究了中等速度撞擊下有初始缺陷矩形薄板的動(dòng)力屈曲，結(jié)果表明初始缺陷和載荷持續(xù)時(shí)間對(duì)板的動(dòng)力屈曲有重要影響，初始缺陷和載荷持續(xù)時(shí)間越小板動(dòng)力屈曲的臨界載荷越大。文獻(xiàn)［6］用解析和數(shù)值計(jì)算方法研究了蜂窩夾心板面板和心板的質(zhì)量分布對(duì)抗沖擊性能的影響。實(shí)驗(yàn)研究方面，以往文獻(xiàn)主要集中在夾心板在爆炸載荷下整體抗沖擊變形吸能的研究［7－9］，而對(duì)單個(gè)夾心基本構(gòu)件動(dòng)力屈曲方面的實(shí)驗(yàn)研究很少。

利用Hopkinson壓桿系統(tǒng)對(duì)矩形薄板在面內(nèi)撞擊下的動(dòng)力屈曲問(wèn)題進(jìn)行實(shí)驗(yàn)研究，取單個(gè)夾心矩形薄板加以分析研究，將實(shí)驗(yàn)結(jié)果與理論計(jì)算結(jié)果進(jìn)行對(duì)照比較，旨在探討矩形薄板在沖擊載荷作用下動(dòng)力屈曲的機(jī)理，對(duì)板的后屈曲吸能方面的研究奠定基礎(chǔ)。實(shí)驗(yàn)結(jié)果表明矩形薄板存在屈曲后強(qiáng)度，在屈曲發(fā)生后面內(nèi)載荷仍然可以繼續(xù)增加。理論與實(shí)驗(yàn)對(duì)比分析結(jié)果表明，當(dāng)壓應(yīng)力波向前傳播一定距離后矩形薄板會(huì)發(fā)生局部屈曲，得到的臨界屈曲長(zhǎng)度與理論計(jì)算結(jié)果吻合良好。

1 實(shí)驗(yàn)裝置及應(yīng)變片位置

實(shí)驗(yàn)主要設(shè)備包括空氣炮、Hopkinson壓桿系統(tǒng)（SHPB）、試件、動(dòng)態(tài)應(yīng)變儀、數(shù)據(jù)采集系統(tǒng)。通過(guò)壓縮空氣推動(dòng)子彈撞擊輸入桿產(chǎn)生一脈沖載荷施加到試件上（見圖1）。子彈、入射桿、透射桿和吸收桿為同一材質(zhì)直徑相同的圓截面鋼桿。

圖1 實(shí)驗(yàn)裝置圖Fig．1 Equipments of the experiment

試件為6061鋁合金夾芯板，沿子彈撞擊方向?yàn)閮蓧K對(duì)稱布置的矩形薄板，為保證入射桿上傳來(lái)的脈沖載荷作用到芯層薄板上為面內(nèi)載荷，薄板兩端垂直嵌于兩矩形厚板上，同時(shí)保證了薄板受載邊和對(duì)邊為夾支。應(yīng)變片貼于芯層矩形板沿寬度方向的正中間，1號(hào)～8號(hào)應(yīng)變片具體布置位置見圖2。

圖2 應(yīng)變片位置圖Fig．2 Locations of the strain gages

圖2中兩厚板之間的距離H＝90．64 mm，夾心矩形板厚度h＝1．53 mm，寬度B＝44．30 mm。1號(hào)和2號(hào)、3號(hào)和4號(hào)、5號(hào)和6號(hào)應(yīng)變片分別沿薄板兩側(cè)對(duì)稱布置，7號(hào)和8號(hào)應(yīng)變片分別貼于入射桿和透射桿上。

2 實(shí)驗(yàn)原理及測(cè)量方案

高壓空氣推動(dòng)子彈以一定初速度與入射桿進(jìn)行高速同軸撞擊，在入射桿中產(chǎn)生彈性壓應(yīng)力波并向前傳播，當(dāng)入射桿中的壓應(yīng)力波傳到試件部位時(shí)一部分被反射回入射桿形成反射拉伸波，其余部分作用到試件上。Hopkinson壓桿實(shí)驗(yàn)在進(jìn)行數(shù)據(jù)處理時(shí)有兩點(diǎn)基本的假設(shè)：①一維應(yīng)力波理論；②假定試件中的應(yīng)力是均勻的。通過(guò)應(yīng)變片1號(hào)～6號(hào)記錄試件中矩形薄板內(nèi)傳播的應(yīng)力波，通過(guò)應(yīng)變片7號(hào)和8號(hào)分別記錄入射桿和透射桿中傳播的應(yīng)力波。

應(yīng)變片1號(hào)～6號(hào)在矩形薄板上沿長(zhǎng)度方向布置，兩兩沿薄板厚度方向?qū)ΨQ分布。當(dāng)矩形薄板無(wú)彎曲變形時(shí)同一測(cè)量點(diǎn)的兩應(yīng)變片信號(hào)是重合的；當(dāng)屈曲發(fā)生時(shí)產(chǎn)生彎曲變形，兩應(yīng)變片中分別增加了拉、壓信號(hào)，使得兩波形分離（分叉），可以記錄臨界屈曲時(shí)間和屈曲時(shí)刻作用到薄板上的應(yīng)力幅值。

3 實(shí)驗(yàn)結(jié)果及理論分析

6061鋁合金材料的應(yīng)力－應(yīng)變曲線見圖3。其中彈性模量Ee＝68．5 GPa，塑性模量Et＝1．24 GPa，彈性極限σs＝310 MPa，密度ρ＝2 700 kg／m3，泊松比υ＝0．33。

圖3 材料的應(yīng)力－應(yīng)變曲線Fig．3 The stress-strain curve of thematerial

實(shí)驗(yàn)中應(yīng)力波在試件中由碰撞端沿長(zhǎng)度方向傳播，動(dòng)態(tài)應(yīng)變儀記錄了試件中各應(yīng)變片位置處薄板的應(yīng)變變化歷程。圖4為應(yīng)變片1號(hào)、2號(hào)測(cè)得的應(yīng)變隨時(shí)間的變化曲線，圖4中時(shí)間軸以應(yīng)力波傳到應(yīng)變片1號(hào)、2號(hào)處的一刻為原點(diǎn)。圖5為應(yīng)變片1、2處應(yīng)力隨時(shí)間變化的曲線。

圖4 應(yīng)變片1號(hào)、2號(hào)測(cè)得的應(yīng)變隨時(shí)間變化的曲線Fig．4 The strain-time curves obtained by strain gage 1 and strain gage 2

圖5 應(yīng)變片1號(hào)、2號(hào)處應(yīng)力隨時(shí)間變化的曲線Fig．5 The stress-time curves in strain gage 1 and strain gage 2

從圖5可以看出矩形薄板在12μs時(shí)刻發(fā)生屈曲，屈曲時(shí)應(yīng)變片1號(hào)、2號(hào)位置處應(yīng)力大小為208．5 MPa，處于彈性范圍之內(nèi)，根據(jù)圖4曲線可以將加載函數(shù)近似為

式中，k＝17．375 MPa／μs。

根據(jù)應(yīng)力波傳播速度以及加載函數(shù)（1）可以近似確定臨界屈曲時(shí)刻試件心板中的應(yīng)力分布見圖7。

圖6 應(yīng)力波在試件中的傳播Fig．6 The propagation of the stress waves

圖7 屈曲時(shí)刻矩形薄板長(zhǎng)度方向應(yīng)力波分布Fig．7 Stress distribution in the thin rectangular plate

從圖5可知，屈曲發(fā)生后面內(nèi)應(yīng)力仍然可以繼續(xù)增加，表明矩形薄板存在屈曲后強(qiáng)度。在應(yīng)力波向前傳播的波前并非突加階躍載荷，而是短時(shí)內(nèi)逐漸加載的，為便于下面用解析方法對(duì)板的屈曲進(jìn)行分析，對(duì)屈曲前的動(dòng)載荷按沖量相等的原則等效簡(jiǎn)化為矩形脈沖載荷，等效簡(jiǎn)化后的載荷幅值為310 MPa，計(jì)算可得等效簡(jiǎn)化后的臨界屈曲長(zhǎng)度為L(zhǎng)cr＝60．4 mm。

下面用解析方法對(duì)夾心矩形板的動(dòng)力屈曲進(jìn)行分析。由圖7可知，塑性應(yīng)力波的幅值不是很大，而且屈曲時(shí)刻在板中傳播的距離相當(dāng)短，分析時(shí)可簡(jiǎn)化為彈性屈曲并假定彈性應(yīng)力波在矩形薄板內(nèi)的傳播是一維的。夾心矩形板在矩型脈沖載荷作用下的模型見圖8，在其中面建立如圖所示直角坐標(biāo)系。

屈曲瞬間對(duì)板微元進(jìn)行動(dòng)力平衡分析，忽略轉(zhuǎn)動(dòng)慣性項(xiàng)，得板的動(dòng)力屈曲平衡方程

圖8 矩形板受力模型示意圖Fig．8 Model of force in the plate

圖9為模型最終變形圖，這里要說(shuō)明一下，塑性大變形主要是由于屈曲發(fā)生后面內(nèi)應(yīng)力和彎曲應(yīng)力的疊加作用導(dǎo)致的，對(duì)于屈曲后變形的發(fā)展，本文不作分析，本文研究的只是臨界屈曲時(shí)刻的彈性屈曲。從圖9可知夾心矩形板屈曲變形模態(tài)沿y方向近似無(wú)變化而只與x方向有關(guān)，故方程（2）中可設(shè)

圖9 模型最終變形圖Fig．9 Final deformation of themodel

對(duì)（2）式分離變量得

夾心矩形板沿x方向在受載邊為夾支，存在如下邊界條件：

臨界屈曲時(shí)刻在壓應(yīng)力波未傳到的部分未受擾動(dòng)，存在如下邊界條件：

式中：L為板沿x方向的臨界屈曲長(zhǎng)度，L＝cetcr，tcr表示臨界屈曲時(shí)間。由屈曲瞬間能量轉(zhuǎn)化率守恒條件可導(dǎo)出壓縮波前附加約束條件：

當(dāng)α2＞2時(shí)，方程（4）有滿足邊界條件和波前約束條件的解

將式（11）代入邊界條件式（6）、式（7）、式（8）、式（9）和波前附加約束條件式（10），得出

由式（15）得到β1和β2，然后將β1，β2和式（16）代入式（11），得

對(duì)于動(dòng)力失穩(wěn)，應(yīng)有λ＝ω2＞0，方程（5）的解為

將式（17）、式（18）代入式（3），得到屈曲位移w為

式中：C為無(wú)限小積分常數(shù)且C＝－C0C2。

式（4）中α2＝Fx／D且α2＝β21+β22，將由式（15）得到β1和β2代入聯(lián)立求解可得臨界屈曲長(zhǎng)度

當(dāng)n＝1時(shí)，L取得最小值，對(duì)應(yīng)本實(shí)驗(yàn)中σx＝310 MPa，由式（19）計(jì)算得到的臨界屈曲長(zhǎng)度為L(zhǎng)＝67．9 mm。理論計(jì)算與實(shí)驗(yàn)結(jié)果之間的誤差為

4 產(chǎn)生誤差的原因及結(jié)論

實(shí)驗(yàn)得到的矩形薄板臨界屈曲長(zhǎng)度比理論計(jì)算的值偏小，導(dǎo)致誤差的原因主要有以下幾點(diǎn)：

（1）矩形薄板兩端達(dá)不到理想的夾支約束，導(dǎo)致實(shí)驗(yàn)中的臨界屈曲長(zhǎng)度偏小。

（2）在沖擊載荷作等效簡(jiǎn)化處理時(shí)應(yīng)力幅值的選取當(dāng)然會(huì)對(duì)臨界屈曲長(zhǎng)度的確定產(chǎn)生影響，應(yīng)力幅值選的越大，按沖量相等的原則得出的臨界屈曲長(zhǎng)度就會(huì)越小。

（3）實(shí)驗(yàn)中沖擊載荷并非理想的矩形脈沖階躍載荷，在應(yīng)力波傳播的波前是短時(shí)內(nèi)逐漸加載的，在確定實(shí)驗(yàn)中臨界屈曲長(zhǎng)度以及做理論計(jì)算時(shí)固定應(yīng)力幅值不變，對(duì)屈曲前的動(dòng)載荷按沖量相等的原則對(duì)此做了等效簡(jiǎn)化處理，也是導(dǎo)致誤差存在的原因之一，更加精確的理論計(jì)算應(yīng)考慮時(shí)變沖擊載荷作用下時(shí)間內(nèi)應(yīng)力分布不均勻的情況，這種情況有待于后續(xù)的研究。

本次實(shí)驗(yàn)利用Hopkinson壓桿系統(tǒng)對(duì)矩形薄板在面內(nèi)撞擊下的彈性動(dòng)力屈曲問(wèn)題進(jìn)行了實(shí)驗(yàn)研究，實(shí)驗(yàn)結(jié)果表明矩形薄板存在屈曲后強(qiáng)度，在屈曲發(fā)生后面內(nèi)載荷仍然可以繼續(xù)增加。通過(guò)記錄屈曲發(fā)生的時(shí)間得到臨界屈曲長(zhǎng)度，并對(duì)實(shí)驗(yàn)?zāi)Ｐ瓦M(jìn)行了理論分析。理論與實(shí)驗(yàn)對(duì)比分析結(jié)果表明，當(dāng)壓應(yīng)力波向前傳播一定距離后矩形薄板發(fā)生局部屈曲，得到的臨界屈曲長(zhǎng)度與理論計(jì)算結(jié)果吻合良好。

［1］王安穩(wěn)．彈性壓應(yīng)力波下直桿動(dòng)力失穩(wěn)的機(jī)理與判據(jù)［J］．力學(xué)學(xué)報(bào)，2001，33（6）：812－820．

WANG An-wen．Mechanism and criterion for dynamic intability of bars under elastic compression wave［J］．Acta Mechanica Sinic，2001，33（6）：812－820．

［2］王安穩(wěn)．軸向壓應(yīng)力波下圓柱殼彈性動(dòng)力失穩(wěn)的判據(jù)與機(jī)理［J］．固體力學(xué)學(xué)報(bào)，2001，22（2）：185－191．

WANG An-wen．Dynamic instability mechanism of cylindrical shells under axial compression wave［J］．Acta Mechanica Solida Sinic，2001，22（2）：185－191．

［3］施連會(huì)，王安穩(wěn)．面內(nèi)階躍載荷下薄板彈性動(dòng)力屈曲差分解［J］．海軍工程大學(xué)學(xué)報(bào)，2008，20（4）：25－29．

SHI Lian-hui，WANG An-wen．Difference solution to elastic dynamic buckling of thin plates subjected to in-plane step load［J］．Journal of Naval University of Engineering，2008，20（4）：25－29．

［4］韓大偉，王安穩(wěn)．三邊簡(jiǎn)支一邊固支矩形薄板動(dòng)力屈曲解析解［J］．振動(dòng)與沖擊，2013，32（2）：140－142．

HAN Da-wei，WANG An-wen．Analytical resolution to dynamic buckling of rectangular thin plates with one edge clamped and the other three edges simply supported［J］．Journal of Vibration and Shock，2013，32（2）：140－142．

［5］Cui Shi-jie，Hao Hong，Cheong H K．Numerical analysis of dynamic buckling of rectangular plates subjected to intermediate-velocity impact［J］．International Journal of Impact Engineering，2001，25：147－167．

［6］Joseph A M，George AG．Uniaxial crushing of sandwich plates under air blast：Influence ofmass distribution［J］．International Journal of Solids and Structures，2008，45：2297－2321．

［7］Nasser S N，Kang W J，McGee J D，et al．Experimental investigation ofenergy-absorption chacteristics of components of sandwich structures［J］．International Journal of Impact Engineering，2007，34：1119－1146．

［8］Zhu Feng，Zhao Long-mao，Lu Guo-xing，et al．Deformation and failure of blast-loaded metallic sandwich panels－Experimental investigations［J］．International Journalof Impact Engineering，2008，35：937－951．

［9］Dharmasena K P，Wadley H N G，Zhenyu X，et al．Mechanical response of metallic honeycomb sandwich panel structures to high-intensity dynamic loading［J］．International Journal of Impact Engineering，2008，35：1063－1074．

Experimental investigation on dynam ic buck ling of thin rectangular plates subjected to in-plane impact

HAN Da-wei，WANG An-wen
（College of Science，Naval University．of Engineering，Wuhan 430033，China）

The dynamic buckling of thin rectangular plates subjected to in-plane impact was investigated experimentally．The result indicates that after buckling the plate can be in-plane loaded further and is of postbuckling strength．The critical buckling length of the thin rectangular plate was obtained by recording the buckling time．The experimentalmodelwas also analysed by using theoreticalmethod．By comparing the results it is indicated that partial buckling occurs in the thin rectangular plate when the compression wave propagates a certain distance．The critical buckling length obtained in the experiment is in reasonable agreementwith the analytical calculation result．

dynamic buckling；thin rectangular plate；critical buckling length；Hopkinson bar；experimental investigation

O347

A

10．13465／j．cnki．jvs．2015．12．016

國(guó)家自然科學(xué)技術(shù)基金（11172330）；海軍工程大學(xué)博士創(chuàng)新基金（HGYJSCXJJ2012005）

2014－03－12 修改稿收到日期：2014－06－19

韓大偉男，博士，1986年生

王安穩(wěn) 男，教授，博士生導(dǎo)師，1945年生