劉澤榮，龍仁榮，張慶明，陳 利

（北京理工大學(xué)爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100081）

多層板結(jié)構(gòu)廣泛存在于航天器、艦船等結(jié)構(gòu)中[1]，然而航天器會(huì)遭受在軌超高速空間碎片的撞擊，艦船等目標(biāo)會(huì)遭受超高速動(dòng)能武器的撞擊[2]。在超高速?gòu)椡枳矒粝?，多層板結(jié)構(gòu)中的第1 層板形成穿孔，第2 層及其后各層板則承受碎片云的沖擊，即分布式載荷。碎片云的破壞形式與超高速?gòu)椡杳黠@不同，其破壞特征與尺寸是毀傷評(píng)估的重要參數(shù)，研究碎片云撞擊下薄板的變形與破壞十分必要。

碎片云碰撞下薄板的變形與破壞可以看作薄板在具有空間分布的瞬態(tài)載荷下發(fā)生塑性動(dòng)態(tài)大變形問(wèn)題。薄壁結(jié)構(gòu)是一種常見(jiàn)的結(jié)構(gòu)，薄板變形理論由來(lái)已久，其中薄板小變形理論研究相對(duì)完善[3]，而薄板動(dòng)態(tài)塑性大變形理論仍在發(fā)展中。在研究薄板的大變形過(guò)程時(shí)，一方面需要同時(shí)考慮彎矩和膜力的作用，所建模型相對(duì)復(fù)雜，另一方面其幾何方程具有非線性項(xiàng)，使得位移求解非常困難。Craggs[4]提出用薄膜理論分析薄板大變形，在模型中不考慮彎矩作用而只考慮膜力，這雖然可以簡(jiǎn)化理論模型，但是卻只適用于很薄且強(qiáng)度較低的薄板大變形撓曲問(wèn)題，難以描述常用金屬薄板大變形問(wèn)題[5]。Reissner[6]首先提出板殼大變形基本理論，隨后Jones[7–8]應(yīng)用此理論求解了理想剛塑性板在動(dòng)載荷下的撓曲變形，通過(guò)忽略面內(nèi)變形，從而簡(jiǎn)化方程得到理論解。Yu 等[9]引入膜力因子，將膜力對(duì)板變形的影響通過(guò)膜力因子乘以平衡方程中的彎矩體現(xiàn)出來(lái)，由于方程中不包含膜力項(xiàng)，求解薄板撓曲變形方程的難度大大降低。然而，這些模型忽略或簡(jiǎn)化了膜力作用，難以實(shí)現(xiàn)大變形高精度的計(jì)算。

本研究利用二級(jí)輕氣炮開(kāi)展球形彈丸超高速正撞擊多層板實(shí)驗(yàn)，分析總結(jié)薄板在碎片云撞擊下的變形與破壞特征，在此基礎(chǔ)上建立考慮彎矩和膜力作用的環(huán)板變形模型，將大變形的薄板材料運(yùn)動(dòng)分解為橫向的外凸運(yùn)動(dòng)和環(huán)向的膨脹運(yùn)動(dòng)，并結(jié)合Grady 破碎理論，建立環(huán)向膨脹花瓣撕裂模型。

1 超高速撞擊多層板實(shí)驗(yàn)及結(jié)果分析

采用二級(jí)輕氣炮將直徑dp=9.13 mm 的LY12 鋁合金球形彈丸加載至速度v0=5 km/s 后撞擊多層板防護(hù)結(jié)構(gòu)。多層板防護(hù)結(jié)構(gòu)由4 塊尺寸為100 mm×100 mm 的LY12 鋁合金薄板組成，相鄰兩板間距為：S1=17 mm，S2=50 mm，S3=38 mm；第1 層至第4 層板的厚度分別為：hb1=1.0 mm，hb2=2.0 mm，hb3=2.0 mm，hb4=1.5 mm，如圖1 所示。

圖1 彈丸超高速碰撞多層板結(jié)構(gòu)Fig. 1 Multi-shock shield structure subjected to hypervelocity impact

實(shí)驗(yàn)后各層板的照片如圖2 所示，第1 層板只發(fā)生穿孔破壞，靶板中心有一個(gè)規(guī)整圓孔。第1 層與第2 層板的間距較小，彈丸與第1 層板碰撞形成的一次碎片云沒(méi)能充分膨脹，碎片云載荷依然較為集中，因此第2 層板中央形成了一個(gè)大的穿孔，孔四周出現(xiàn)一個(gè)由小碎片撞擊形成的環(huán)狀小坑區(qū)域。第3 層板與第4 層板的破壞特征相似，板中間有穿孔破壞，孔周邊材料發(fā)生拉伸完全變形，并且發(fā)生撕裂破壞，如圖2(c)、圖2(d)、圖2(e)所示。第3 層板彎曲撕裂至靶板邊界，邊沿夾持的卡具也發(fā)生了彎曲大變形。由于第3 層板的阻擋，作用到第4 層板的碎片云能量大幅減小，第4 層板的彎曲變形主要集中在中部區(qū)域，邊沿和卡具均沒(méi)有明顯的彎曲變形。由于裂紋擴(kuò)展所需要的能量較小，第4 層板的裂紋依然擴(kuò)展至靶板邊沿。第4 層板表面覆蓋一層白色物體，彈坑也相對(duì)較少，經(jīng)分析，白色物是熔融鋁冷卻凝固物，熔融態(tài)的碎片云對(duì)第4 層板的作用更趨于整體作用，而不是大尺寸固態(tài)碎片的局部作用。整個(gè)多層板結(jié)構(gòu)在超高速?gòu)椡枳矒粝碌淖冃翁卣魅鐖D1 所示。

圖2 高速?zèng)_擊后各層板的破壞照片F(xiàn)ig. 2 Damage results of each layer after hypervelocity impact

實(shí)驗(yàn)結(jié)果表明，彈丸與集中分布的碎片云撞擊薄板會(huì)造成穿孔破壞，隨著撞擊次數(shù)的增加，彈靶材料破碎得更細(xì)，且熔融態(tài)的比例增加，碎片云載荷分布越分散，薄板的變形與破壞模式轉(zhuǎn)變?yōu)橹行拇┛准碍h(huán)孔區(qū)域大凹陷變形并伴有撕裂。第3 層板阻擋了大部分碎片云能量，第4 層板的彎曲變形比第3 層板小很多，其中心穿孔直徑的測(cè)量平均值為46 mm，凹陷變形區(qū)域直徑為穿孔直徑的2～3 倍，凹陷板與靶板原平面的夾角較小，約為12°，說(shuō)明其凹陷變形的速度較小。薄板穿孔四周材料的凹陷變形可簡(jiǎn)單類(lèi)比成圓環(huán)的膨脹變形，當(dāng)膨脹速度大于某值后，圓環(huán)在膨脹過(guò)程會(huì)斷裂成多段。材料起裂后裂紋在板中傳播形成長(zhǎng)裂紋，與實(shí)驗(yàn)中薄板撕裂成多塊的結(jié)果一致。根據(jù)前述分析，第4 層板環(huán)孔凹陷區(qū)域的變形速度較低，其環(huán)向膨脹更趨近于薄板起裂的臨界速度。

2 碎片云作用下環(huán)形薄板變形模型

由實(shí)驗(yàn)結(jié)果總結(jié)得到在分散碎片云載荷下薄板的響應(yīng)特征如下：薄板中間區(qū)域材料在碎片云中部高幅值沖擊載荷下發(fā)生穿孔，環(huán)孔區(qū)域的環(huán)向板承受的碎片云載荷稍微小一點(diǎn)，環(huán)板材料發(fā)生凹陷變形，孔邊材料在凹陷膨脹過(guò)程中破碎成幾段，裂紋擴(kuò)展直至薄板裂成幾塊。碎片云速度分布與環(huán)板變形如圖3 所示。碎片云撞擊薄板穿孔尺寸的計(jì)算已有研究[10]，在此不做考慮，本研究側(cè)重碎片云載荷下環(huán)形薄板的變形建模。

圖3 環(huán)板變形側(cè)視圖Fig. 3 Lateral schematic of annular plate deformation

2.1 薄板變形假設(shè)與一般方程

根據(jù)實(shí)驗(yàn)結(jié)果可知，碎片云分布載荷作用下的變形屬于典型的薄板大變形問(wèn)題。薄板在變形過(guò)程中發(fā)生彎曲，需要考慮彎矩作用，另外薄板面內(nèi)還有膜應(yīng)力。要建立薄板變形的一般方程，需作如下基本假設(shè)：

(1) 考慮到薄板撓度是板厚的幾倍甚至幾十倍，為此忽略板的彈性小變形，采用理想剛塑性體模型描述薄板變形；

(2) 薄板變形服從Love-Kirchhoff 假設(shè)[11]，即中面法線變形后仍保持為中面法線，厚度變化忽略不計(jì)，同時(shí)略去中面法向應(yīng)力；

(3) 不區(qū)分變形前后參考位形的差異，采用變形前的位形作為參考基準(zhǔn)建立平衡方程；

(4) 薄板內(nèi)力除徑向彎矩Mr和環(huán)向彎矩Mθ外，還需要考慮徑向膜力Nr和環(huán)向膜力Nθ。

基于以上假設(shè)，可以建立軸對(duì)稱(chēng)薄板變形過(guò)程的橫向、環(huán)向和徑向動(dòng)力學(xué)方程[12]

式中：r為徑向坐標(biāo)，θ 為環(huán)向坐標(biāo)，u為徑向位移，w為撓度，Qr為剪力，μ為板的面密度，p為軸對(duì)稱(chēng)的橫向分布載荷。

考慮面內(nèi)變形與出面變形的薄板軸對(duì)稱(chēng)幾何方程，其內(nèi)力（或稱(chēng)廣義應(yīng)力[3]）徑向彎矩Mr、環(huán)向彎矩Mθ、徑向膜力Nr和環(huán)向膜力Nθ對(duì)應(yīng)的廣義應(yīng)變分別為徑向膜應(yīng)變?chǔ)舝、環(huán)向膜應(yīng)變?chǔ)纽?、徑向曲率κr和環(huán)向曲率κθ。幾何方程如下

采用理想剛塑性體模型，材料發(fā)生剛性位移或塑性變形，則薄板發(fā)生塑性變形時(shí)膜力Nr、Nθ達(dá)到N0，彎矩Mr、Mθ達(dá)到M0

式中：h為板厚度，σy為屈服強(qiáng)度。

薄板材料是理想剛塑性體模型，薄板的應(yīng)力狀態(tài)位于屈服面上，為此采用Hodge[13]提出的包含彎矩、膜力作用的近似屈服準(zhǔn)則，如圖4 所示?？梢钥闯?，屈服面上的膜力和彎矩獨(dú)立作用，各自位于Tresca 屈服面上。這樣的近似處理使求解過(guò)程變得非常簡(jiǎn)便。另外，橫向剪力Qr出現(xiàn)在平衡方程中，但是由于不考慮其產(chǎn)生的剪切變形，因此不出現(xiàn)在屈服條件中，并且后面計(jì)算中可以看到，化簡(jiǎn)的平衡方程不包含剪力。

圖4 中mr、mθ、nr、nθ為無(wú)量綱內(nèi)力，mr=Mr/M0，mθ=Mθ/M0，nr=Nr/N0，nθ=Nθ/N0；A、B、C、D、E、F是屈服條件上的角點(diǎn)，各線段代表材料的某一內(nèi)力狀態(tài)。正交法則規(guī)定一點(diǎn)廣義塑性應(yīng)變?cè)隽康姆较蚺c屈服面外法線方向相同，如圖4 所示。對(duì)于具體的問(wèn)題，在約束允許的機(jī)動(dòng)場(chǎng)中，根據(jù)廣義塑性應(yīng)變?cè)隽浚ɑ驈V義塑性應(yīng)變率）的方向，可以確定塑性變形由哪一部分屈服面控制，從而為具體問(wèn)題的求解確定內(nèi)力大小。

圖4 Tresca 屈服條件與正交法則Fig. 4 Tresca condition and flow rule

2.2 環(huán)板變形方程簡(jiǎn)化

根據(jù)薄板變形方程，結(jié)合屈服條件和正交法則，就可以求解薄板變形的初邊值問(wèn)題，從而得到環(huán)板的變形。然而多個(gè)偏微分方程的相互耦合以及非線性項(xiàng)的引入，使方程很難得到解析解，需要根據(jù)本問(wèn)題的特征對(duì)方程進(jìn)行簡(jiǎn)化。理想剛塑性體的使用、碎片云撞擊載荷的特性以及環(huán)板構(gòu)型的特殊性會(huì)使方程變量減少，復(fù)雜程度降低。

如圖5 所示，薄板中部發(fā)生穿孔（0 l）在本問(wèn)題中則可認(rèn)為是未變形區(qū)?？梢哉J(rèn)為碎片云沖擊載荷為瞬態(tài)載荷作用于環(huán)板面上。根據(jù)已有的碎片云模型[10]，得到作用在環(huán)板上的沖量面密度分布i=i(r)，碰撞過(guò)程瞬時(shí)完成，將沖量載荷轉(zhuǎn)化為環(huán)板橫向速度初始條件v(r) =i/μ(t= 0)，如圖3 所示。這樣將碎片云碰撞薄板問(wèn)題簡(jiǎn)化為圓環(huán)板在初始橫向速度條件下的變形問(wèn)題。

圖5 環(huán)板變形俯視圖Fig. 5 Top schematic of annular plate deformation

剛塑性環(huán)板撓曲變形采用錐形機(jī)動(dòng)場(chǎng)[14]，即a≤r

代入式(14)，得

2.3 動(dòng)力學(xué)方程求解

由于只有單一未知量，使用能量方程微分形式求解更方便，即塑性鉸上塑性功率與動(dòng)能變化率平衡

對(duì)于a

代入式(34)得到環(huán)板凹陷變形終態(tài)角度θt

由式(33)和式(34)即可求出環(huán)板變形最終撓度wta和終態(tài)角度θt，由式(29)可得到環(huán)板變形與速度的關(guān)系式。下面將計(jì)算環(huán)板在碎片云撞擊下的具體變形問(wèn)題，得出環(huán)板徑向速度，進(jìn)而計(jì)算環(huán)板花瓣型撕裂。

3 計(jì)算驗(yàn)證

圖6 撓度與徑向關(guān)系示意圖Fig. 6 Schematic of relationship between deflection and radial displacement

圖7 沖量載荷下環(huán)板變形示意圖Fig. 7 Deflection of annular plate subjected to impulse

為計(jì)算環(huán)板撕裂起始時(shí)的變形，首先計(jì)算出環(huán)板在碎片云撞擊下達(dá)到的最大徑向速度。為得到環(huán)板變形時(shí)材料受力狀態(tài)，將板厚h= 1.5 mm、材料屈服強(qiáng)度 σy= 265 MPa 代入式(8)，求得剛塑性彎矩M0和剛塑性膜力N0。根據(jù)多次碎片云模型[10]，可以計(jì)算本研究實(shí)驗(yàn)工況下3 次碎片云對(duì)第4 層環(huán)板施加的沖量面密度i（見(jiàn)圖8）以及碎片云作用范圍的半徑l= 47.2 mm，再通過(guò)式(30)～式(32)求得碎片云撞擊環(huán)板時(shí)環(huán)板的初始動(dòng)能Ed0= 97.5 J。由式(29)和式(17)得到內(nèi)環(huán)邊徑向速度u˙a隨內(nèi)環(huán)邊撓度wa變化曲線，如圖9 所示，得到最大內(nèi)環(huán)邊徑向速度u˙ma為25.99 m/s。

圖8 沖量面密度載荷分布曲線Fig. 8 Distribution of the impulse per area

圖9 環(huán)板內(nèi)環(huán)邊徑向速度與撓度的關(guān)系Fig. 9 Relationship between inner edge radial velocity and deflection

上述模型可以用于描述碎片云載荷分布較分散的環(huán)形薄板變形，其環(huán)板凹陷變形區(qū)域集中在碎片云載荷分布區(qū)域內(nèi)，由于載荷不太集中，環(huán)板的變形不太快，變形區(qū)與未變形區(qū)分界的塑性鉸采用駐定塑性鉸，其計(jì)算結(jié)果與第4 層板的實(shí)驗(yàn)結(jié)果相符合。實(shí)驗(yàn)中二次碎片云撞擊第3 層板時(shí)的沖量面密度相對(duì)較高，第3 層板凹陷變形的速度較快，拉著未變形區(qū)域發(fā)生凹陷，其變形區(qū)與未變形區(qū)的塑性鉸是移動(dòng)的，因此第3 層板的變形需采用移行的塑性鉸圓來(lái)描述。對(duì)于移行鉸，l=l(t)，l˙ ≠0，由圖6 可知，這時(shí)內(nèi)環(huán)邊撓度wa和徑向位移ua為

4 結(jié) 論

通過(guò)超高速撞擊實(shí)驗(yàn)總結(jié)了薄板在軸對(duì)稱(chēng)碎片云載荷作用下的變形與破壞模式，并對(duì)此完成了建模研究，結(jié)論如下：

(1) 實(shí)驗(yàn)中薄板在分散的碎片云載荷下經(jīng)歷中部穿孔，形成環(huán)板，環(huán)板在剩余碎片云載荷下發(fā)生凹陷，孔邊材料繼而發(fā)生拉伸破碎，形成花瓣型撕裂；

(2) 基于實(shí)驗(yàn)現(xiàn)象與經(jīng)典理論，建立了環(huán)形薄板在碎片云沖擊下的剛塑性變形模型，該模型同時(shí)考慮彎矩和膜力，通過(guò)簡(jiǎn)化得到了環(huán)板變形微分方程，可以求解得到板橫向速度場(chǎng)和徑向速度場(chǎng)；

(3)將孔周邊的薄板材料膨脹視為圓環(huán)的膨脹，根據(jù)模型計(jì)算得到環(huán)板內(nèi)環(huán)最大環(huán)向應(yīng)變率，代入Grady 破碎模型，計(jì)算實(shí)驗(yàn)中第4 層板撕裂的花瓣數(shù)，另外根據(jù)模型還計(jì)算得到環(huán)板凹陷變形終態(tài)角度，計(jì)算結(jié)果均與實(shí)驗(yàn)結(jié)果吻合；

(4) 對(duì)于高幅值碎片載荷的環(huán)形薄板變形，宜采用移行塑性鉸描述變形區(qū)域的擴(kuò)展過(guò)程，其變形的理論建模和求解還需要在本研究的基礎(chǔ)上進(jìn)一步開(kāi)展。