趙 凱，高光發(fā)，王肖鈞

(中國科學(xué)技術(shù)大學(xué) 近代力學(xué)系，合肥 230027)

柱殼是被廣泛應(yīng)用于大型油氣存儲(chǔ)工程，水下航行器，航空航天飛行器，各種殺傷性彈藥武器的典型結(jié)構(gòu)的代表，如何提高柱殼結(jié)構(gòu)的抗沖擊性能一直以來都是各國工程和防護(hù)領(lǐng)域科研工作者的重要研究課題。早期Rath等[1]基于薄殼理論研究了柱殼的振動(dòng)問題，Qatu[2-3]提出了中厚復(fù)合殼的低階近似振動(dòng)方程，薄壁圓柱殼在動(dòng)態(tài)和靜態(tài)載荷下的失穩(wěn)屈曲，一直以來都是學(xué)者們重點(diǎn)關(guān)注的課題[4-6]。近年來，Ebrahim等[7]采用各種三維剪切理論進(jìn)一步研究了厚殼的靜態(tài)和震動(dòng)問題，Liu等[8]基于廣義射線矩陣法和回傳射線矩陣法研究了沖擊載荷下層合柱殼的瞬態(tài)響應(yīng)問題，Hu等[9]將薄殼理論和流固耦合方法RAVS相結(jié)合，提出了一種研究水中爆炸載荷下層合柱殼的瞬態(tài)響應(yīng)計(jì)算近似方法，戴向勝等[10]用修正的Singace疊縮模型，研究了沖擊載荷下含應(yīng)變強(qiáng)化效應(yīng)、應(yīng)變率強(qiáng)化效應(yīng)和溫度效應(yīng)的圓柱殼吸能情況。然而，這些理論方法畢竟具有一定的應(yīng)用局限性。

在現(xiàn)代科學(xué)技術(shù)研究中，許多問題的解決都需要通過大量的實(shí)驗(yàn)才能實(shí)現(xiàn)。然而，實(shí)物的原型試驗(yàn)存在復(fù)雜性和難度很高、耗資巨大等種種問題，有時(shí)甚至是不可能實(shí)現(xiàn)的。因此在了解相應(yīng)物理現(xiàn)象的基本規(guī)律的前提下，需要通過所謂模擬實(shí)驗(yàn)(或數(shù)值模擬)取得所研究問題中的一些基本規(guī)律和特征量，用于實(shí)物原型上?；诹烤V分析中的相似律[11]，能讓我們?cè)谀Ｐ脱芯恐蝎@得某種規(guī)律或特征量后，以無量綱特征量的方式直接用于原型，因此在學(xué)術(shù)研究和工程實(shí)踐中都具有廣泛的應(yīng)用價(jià)值。本文以沖擊載荷下柱殼的瞬態(tài)響應(yīng)為背景，通過量綱分析研究了可簡化為平面應(yīng)變問題的柱殼在沖擊載荷下的最大等效應(yīng)力和最大等效應(yīng)變所依賴的主要無量綱自變量，并開展相應(yīng)的數(shù)值模擬，得到了不同相似比時(shí)的具體函數(shù)關(guān)系。

1 雙層柱殼結(jié)構(gòu)抗沖擊性能的量綱分析

考慮沖擊載荷作用下雙層圓柱殼的瞬態(tài)響應(yīng)問題。設(shè)給定一個(gè)長度為L、外半徑為R的圓柱殼，殼體由兩層不同材料組成(如圖1(a)所示)。假設(shè)兩層材料都滿足各向同性理想塑性材料的Mises屈服準(zhǔn)則。材料1殼環(huán)的厚度為H1，材料性能參數(shù)為ρ1(質(zhì)量密度)、E1(楊氏模量)、ν1(泊松比)、Y1(簡單拉壓屈服應(yīng)力)；材料2殼環(huán)的厚度為H2，材料性能參數(shù)為ρ2(質(zhì)量密度)、E2(楊氏模量)、ν2(泊松比)、Y2(簡單拉壓屈服應(yīng)力)。

殼體的任何響應(yīng)量都是外載荷的空間分布和歷時(shí)函數(shù)的泛函，我們考慮如下的外載：假設(shè)殼體在180°張角的外側(cè)沿軸向均勻分布、而沿θ角按余弦規(guī)律分布的壓力載荷p(t)cosθ的作用(如圖1(b)所示)，載荷在θ=0°處的壓力峰值峰值為p0、作用總時(shí)間為τ0、加卸載規(guī)律為三角形脈沖(如圖1(c)所示)。

圖1 柱殼上表面受均布載荷p的示意圖

上述分布規(guī)律的外載荷作用下的柱殼的動(dòng)力響應(yīng)顯然只依賴于載荷參數(shù)p0和τ0?？紤]柱殼受載區(qū)域的中心點(diǎn)A處的最大Mises等效應(yīng)力σmax和最大Mises等效應(yīng)變?chǔ)舖ax隨外載的變化規(guī)律。如果柱殼材料是完全理想彈塑性材料，則我們可以確定如下一些物理量是問題的主定量(自變量)：柱體長度L、外徑R、殼環(huán)厚H1和H2、材料性能參數(shù)ρ1,E1,ν1,Y1和ρ2,E2,ν2,Y2、外載荷參數(shù)p0和τ0，于是圓柱殼中任何點(diǎn)處(包括A點(diǎn))的最大Mises等效應(yīng)變?chǔ)舖ax、最大Mises等效應(yīng)力σmax應(yīng)是這些物理量的函數(shù)，即

(1)

這里f(x)、g(x)表示某種函數(shù)關(guān)系。采用LMT單位制分析，式(1)中，因變量應(yīng)變?chǔ)舖ax本身已是無量綱量；L、R、H1、H2的量綱均是長度L；ρ1、ρ2表示材料的質(zhì)量密度，其量綱可以表示為ML-3；因變量σmax和自變量E1、E2、Y1、Y2和p0都具有應(yīng)力的量綱，可以表示為ML-1T-2；自變量ν1和ν2本身已是無量綱量；自變量τ0的量綱為時(shí)間T。

上述2個(gè)函數(shù)關(guān)系式在LMT單位制里，共有三個(gè)獨(dú)立的量綱，可以取其中任何3個(gè)量綱彼此獨(dú)立的物理量作為其基本量。我們將基本量組取為：圓柱殼體的半徑R、材料1的質(zhì)量密度ρ1和材料1的楊氏模量E1，因?yàn)檫@3個(gè)量是獨(dú)立的、有較清晰的物理意義的量，且其它的因變量和自變量的量綱都可由它們的量綱表達(dá)出來。

(2)

在(2) 式中，對(duì)無量綱因變量起作用的無量綱自變量有11個(gè)：前3個(gè)是幾何長度之比，即幾何相似準(zhǔn)數(shù)，它們反映了對(duì)結(jié)構(gòu)幾何尺度的幾何相似要求；中間6個(gè)則是材料歸一化力學(xué)性能參數(shù)，即材料相似準(zhǔn)數(shù)，它們反映了對(duì)結(jié)構(gòu)材料相似性的要求；倒數(shù)第2個(gè)表示以第一種材料的楊氏模量為基準(zhǔn)而量度的外載荷歸一化強(qiáng)度，即外載的動(dòng)力學(xué)相似準(zhǔn)數(shù)，它反映了對(duì)外載荷強(qiáng)度的動(dòng)力相似要求；最后1個(gè)表示以第一種材料桿中彈性縱波波速傳過結(jié)構(gòu)特征尺度R所需時(shí)間為基準(zhǔn)而量度的外載荷歸一化歷時(shí)。

函數(shù)關(guān)系(2)式的具體形式并不能單純地由量綱分析方法而得出，而我們的任務(wù)就是通過實(shí)驗(yàn)、或者理論分析、或者模擬計(jì)算的方法確定上述多元函數(shù)的具體形式。

我們只討論柱殼長度L足夠大、而且整個(gè)殼的長度上都受有如上載荷的問題，則三維的柱體響應(yīng)問題可以視為平面應(yīng)變問題，此時(shí)L的大小將不影響柱殼中的歸一化最大Mises等效應(yīng)變和歸一化最大Mises等效應(yīng)力(嚴(yán)格而言當(dāng)殼體的長度L為有限時(shí)還應(yīng)引入殼兩端的邊界條件及相關(guān)的外載參數(shù))。因此(2)式便簡化為平面應(yīng)變問題的動(dòng)力響應(yīng)問題，即

(3)

如果我們要用模型實(shí)驗(yàn)或者數(shù)值模擬的方法來尋求無量綱函數(shù)的具體形式，則必須保證模型問題(用下標(biāo)m表述)和原型問題(用下標(biāo)p表示)的完全物理相似，即

(εmax)m=(εmax)p,(σmax/E1)m=(σmax/E1)p

(4)

(5)

如果在下述條件下開展模擬實(shí)驗(yàn)或者進(jìn)行模擬計(jì)算，可以進(jìn)一步作如下簡化：

①保持模型和原型滿足幾何相似條件，即

(6)

②假定模型和原型使用相同的材料，則自然就保證了材料相似，即有，

(7)

因而我們需要考慮的相似準(zhǔn)數(shù)便只剩下外載荷的動(dòng)力相似準(zhǔn)數(shù)，此時(shí)(3)式可簡化為

(8)

于是，在保證模型與原型幾何相似和材料完全相同的前提下，我們要研究結(jié)構(gòu)的響應(yīng)規(guī)律(如某點(diǎn)的無量綱化最大Mises等效應(yīng)力(σmax/E1)m和最大Mises等效應(yīng)變值(εmax)m，只需通過模擬實(shí)驗(yàn)或者模擬計(jì)算的方法求出二元函數(shù)(8)式的相應(yīng)數(shù)據(jù)，并通過擬合而求出這兩個(gè)二元函數(shù)所表達(dá)的曲面就行了，也就是說，二元函數(shù)(8)式即柱殼最大等效應(yīng)力和應(yīng)變的具體依賴函數(shù)關(guān)系。

需要指出的是，雖然(8)式中沒有顯含材料的屈服強(qiáng)度Y，但是，無論對(duì)彈性材料還是對(duì)理想塑性材料，甚至不難證明對(duì)有硬化效應(yīng)的彈塑性材料，在保持模型和原型幾何相似及材料相似的前提下，上述問題的無量綱因變量與自變量的依賴關(guān)系都可由二元函數(shù)(8)所表達(dá)。對(duì)彈性材料和彈塑性材料進(jìn)行模擬實(shí)驗(yàn)或模擬計(jì)算的時(shí)候，所得出的二元函數(shù)f或g的形式則必然是不同的，因?yàn)閷?duì)彈性問題屈服應(yīng)力Y不會(huì)發(fā)生作用，而對(duì)彈塑性問題要求材料進(jìn)入屈服，所以屈服應(yīng)力Y自然會(huì)發(fā)生作用。這里面隱含著另兩個(gè)無量綱參數(shù)Y1/E1和Y2/E2的作用，但是在保持模型和原型材料完全一樣的前提下，它們就是絕對(duì)常數(shù)，并不增加函數(shù)的復(fù)雜性，所以我們要尋求的仍然是二元函數(shù)。

2 雙層柱殼結(jié)構(gòu)抗沖擊性能二元函數(shù)的數(shù)值模擬

采用國際商用沖擊動(dòng)力學(xué)軟件LSDYNA，對(duì)上述問題進(jìn)行數(shù)值模擬。柱殼外層材料是鋼，內(nèi)層材料是鋁，采用線性硬化彈塑性本構(gòu)模型，其基本力學(xué)參數(shù)如表1所示。柱殼的原型尺寸為，半徑R=2 000 mm，H1=20 mm，H2=10 mm。原型上作用的脈沖載荷強(qiáng)度p0=500 MPa，τ0=20 μs。

表1 材料基本力學(xué)性能

通過上述25個(gè)算例，可以得到不同半徑、不同強(qiáng)度載荷下柱殼A點(diǎn)最大等效應(yīng)力、最大應(yīng)變的計(jì)算結(jié)果，如表2、表3所示。

表2 A點(diǎn)最大等效應(yīng)力(單位：MPa)

表3 A點(diǎn)最大等效應(yīng)變

圖2 σmax/E1 vs. p0 /E1

圖3 εmax vs. p0 /E1

圖4 無量綱最大等效應(yīng)力函數(shù)曲面

圖5 無量綱最大等效應(yīng)變函數(shù)曲面

(9)

可得，a1=-2.540 45，b1=0.747 07，c1=35.233 81，d1=3 513.485 71，f1=-180.818 56，g1=0.045 62；a2=-0.050 28，b2=1.240 49，c2=-0.014 62，d2=-469.714 29，f2=22.848 68，g2=0.002 13。

3 結(jié) 論

以雙層柱殼結(jié)構(gòu)在脈沖載荷作用下的瞬態(tài)響應(yīng)為例，通過量綱分析得出了可簡化為平面應(yīng)變問題的柱殼在沖擊載荷下某點(diǎn)最大等效應(yīng)力和等效應(yīng)變無量綱因變量對(duì)主要無量綱自變量的函數(shù)依賴關(guān)系，并在模型和原型幾何相似、所取材料相同的條件下，通過數(shù)值模擬得到了具體的二元函數(shù)表達(dá)，得出以下結(jié)論：

(2)數(shù)值模擬顯示，在載荷強(qiáng)度不大、模型尺寸較大時(shí)，柱殼的最大等效應(yīng)力和等效應(yīng)變的無量綱函數(shù)呈線性，而隨載荷強(qiáng)度的增大、模型的縮小，函數(shù)的非線性增長趨勢(shì)顯著。

[1]Rath B K, Das Y C. Vibration of layered shells[J]. Journal Sound and Vibration, 1973,28(4):737-757.

[2]Qatu M S. Accurate stress resultants equations for laminated composite deepthick shells[J]. Composites for the Pressure Vessel Industry, 1995(302):39-46.

[3]Qatu M S. Accurate equations for laminated composite deep thick shells[J]. Interational Jourmal of Solids and Structures,1999,36(19):2917-2941.

[4]Yamaki N．Elastic stability of circular cylindrical shell[M]．Amsterdam：North-Holland，1984．

[5]顧紅軍，趙國志，陸廷金，等.軸向沖擊下薄壁圓柱殼的屈曲行為的實(shí)驗(yàn)研究[J].振動(dòng)與沖擊，2004,23(4):58-63.

GU Hong-jun,ZHAO Guo-zhi,LU Ting-jin,et al.Buckling of thin-wall cylindrical shell under axial impact[J].Journai of Vibration and Shock,2004,23(4):58-63.

[6]張善元，張濤.圓柱殼的軸向動(dòng)力屈曲、參數(shù)共振與混沌運(yùn)動(dòng)[J].振動(dòng)與沖擊，2010,29(12)2:34-38.

ZHANG Shan-yuan,ZHANG Tao.Axial dynamic buckling,parametric resonance and chaotic motion of a closed cylindrical shell[J].Journai of Vibration and Shock,2010,29(12):34-38.

[7]Asadi E, Wang Wen-chao, Qatu M S. Static and vibration analyses of thick deep laminated cylindrical shells using 3D and various shear deformation theories[J]. Composite Structures, 2012,94(2): 494-500.

[8]LIU Chun-chuan, LI Feng-ming, HUANG Wen-hu. Transient wave propagation and early short time transient responses of laminated composite cylindrical shells[J]. Composite Structures, 2011,93(10): 2587-2597.

[9]HU Gang-yi, XIA Fei, LI Jun, The transient responses of two-layered cylindrical shells attacked by underwater explosive shock waves[J]. Composite Structures, 2010,92(7): 1551-1560.

[10]戴向勝,馬建敏.沖擊載荷作用下金屬圓柱殼能量吸收研究[J].振動(dòng)與沖擊,2012,31(6):100-103.

DAI Xiang-sheng,MA Jian-min.Energy absorbed by a metal tube under axial crush load[J].Journal of Vibration and Shock,2012,31(6):110-103.

[11]談慶明. 量綱分析[M]. 合肥：中國科學(xué)技術(shù)大學(xué)出版社, 2005.