基于能量法研究剛性質(zhì)量塊軸向撞擊圓柱殼的屈曲

辛朝煒， 韓志軍

(1. 太原理工大學 機械與運載工程學院，太原 030024； 2. 太原理工大學 生物醫(yī)學工程學院，太原 030024)

隨著現(xiàn)代科學技術迅猛發(fā)展，結(jié)構屈曲的研究受到廣大學者的極大關注。目前靜力屈曲問題的研究已較為充分，而動力屈曲問題的研究仍充滿挑戰(zhàn)。作為結(jié)構基本原件的圓柱殼廣泛應用于實際生活和各種工程實際中，例如：各種飛行器以及運載火箭主體結(jié)構均采用殼體。所以圓柱殼動力屈曲的研究有著很強的工程應用價值、理論價值和重要意義。

Zimicik等[1]忽略應力波效應，用理論和實驗研究了滿載作用下，徑厚比不同的薄壁圓柱殼的軸向動力屈曲，指出圓柱殼的屈曲形式與初缺陷的分布有著密切的關系。Lo Lepik[2]考慮應力波效應，用準分叉法分析了圓柱殼的動力屈曲，用Runge-Kutta法給出其數(shù)值解。Bisagni[3-4]忽略應力波效應，使用Abaqus軟件對有初缺陷的圓柱殼在軸向沖擊載荷下的動力屈曲進行了數(shù)值分析，并用實驗驗證了數(shù)值模型，數(shù)值計算結(jié)果表明初缺陷和載荷持續(xù)時間對動力屈曲有影響。Karagiozova等[5-6]考慮應力波效應的影響，給出了圓柱殼發(fā)生動力屈曲的條件，用FE碼Abaqus/Standard對其進行了數(shù)值模擬，表明材料和幾何特性對屈曲是有影響的；其實驗研究，表明沖擊速度對圓柱殼動力屈曲有大的影響，而且材料特性和慣性效應對于屈曲模態(tài)的形成有重大影響。Xu等[7-8]通過將哈密頓系統(tǒng)引入圓柱殼的動態(tài)屈曲中，考慮應力波效應，得到了哈密頓正則方程的數(shù)值解。Fan等[9-10]不考慮應力波效應，導出了圓柱殼在滿載作用下屈曲臨界荷載隨時間變化的解析式，并討論了圓柱殼徑厚比和外壓加載速度對屈曲臨界載荷的影響。桂夷斐等[11]研究了在軸向階躍載荷作用下，應力波對圓柱殼動力屈曲的影響，通過數(shù)值計算，給出了屈曲臨界載荷與臨界長度、圓柱殼徑厚比的關系。上述研究表明，考慮應力波效應時一般采用階躍載荷加載，忽略應力波效應時一般以滿載方式加載，且研究問題的方法有理論分析、數(shù)值模擬或計算機仿真、實驗研究，這樣考慮問題的角度和所用方法不同，所得結(jié)論各異，且對具有工程應用背景的質(zhì)量塊撞擊圓柱殼的動力屈曲研究較少也不夠深入。Karagiozova等[12-13]考慮應力波效應，利用 Abaqus對在質(zhì)量塊沖擊作用下的圓柱殼進行模擬并得到不同的屈曲模態(tài)，結(jié)果表明，受到高速沖擊的殼體，其動力屈曲不僅取決于施加載荷的大小，還取決于初始動能的大小。顧紅軍等[14]采用落錘實驗對圓柱殼的動力屈曲行為進行研究，結(jié)果表明圓柱殼的動力屈曲與徑厚比有關。馬健等[15]對質(zhì)量塊撞擊圓柱殼進行實驗研究，運用LS-DYNA對其進行了仿真模擬，兩者得到的沖擊加速度響應曲線基本一致。上述質(zhì)量塊撞擊圓柱殼的屈曲表明，因應力波傳播比較復雜，這給理論研究帶來一定困難，目前主要集中于實驗和計算機仿真分析。

基于此，本文考慮應力波效應，利用能量法，建立系統(tǒng)的拉格朗日函數(shù)，將經(jīng)計算獲得的符合邊界條件的準試函數(shù)代入第二類拉格朗日方程，得到二階線性偏微分方程，分析方程解的特性，得到質(zhì)量塊撞擊圓柱殼屈曲臨界速度的解析表達式。算例分析，討論了臨界長度、質(zhì)量比、軸向模態(tài)數(shù)、環(huán)向模態(tài)數(shù)、徑厚比等參數(shù)對屈曲臨界速度的影響。

1 屈曲臨界速度解析表達式的推導

如圖1所示，均質(zhì)鋼圓柱殼左端夾支，右端固支，長度為L，半徑為R，厚度為h，彈性模量為E，密度為ρ，泊松比μ，橫截面積為A，選取柱坐標系x,y,z，其相應的位移為u,v,w。左端受質(zhì)量為M的剛性質(zhì)量塊撞擊。應力波沿x方向傳播，如圖2所示，當傳播到距夾支端為lcr時，圓柱殼受力段發(fā)生突然屈曲，且lcr為臨界長度，其中l(wèi)cr=tcr·c。

圖1 受剛性質(zhì)量塊撞擊的圓柱殼

圖2 應力波傳播示意圖

應力波在圓柱殼內(nèi)傳播[16]

(1)

于是圓柱殼內(nèi)各段軸力

(2)

當應力波傳到lcr時，系統(tǒng)的動能T、應變能U，以及外力功W，可表示為

(3)

Mθκθ+Mxθκxθ}Rdxdθ

(4)

(5)

在小撓度理論下圓柱殼的幾何關系為

(6)

(7)

圓柱殼的內(nèi)力

(8)

內(nèi)力矩為

(9)

忽略軸向慣性，整理系統(tǒng)動能T、系統(tǒng)應變能U、系統(tǒng)做功W的表達式得到

(10)

(11)

(12)

設系統(tǒng)的拉格朗日函數(shù)為L，有

L=T-U+W

(13)

把式(10)、式(11)、式(12)代入式(13)，得到

(14)

可設試函數(shù)為

w=f(x,θ)T(t)

(15)

其中,

(16)

式中：n為軸向模態(tài)數(shù)；m為環(huán)向模態(tài)數(shù)；α和β為相關待定系數(shù)。式(15)應滿足下述邊界條件。

(17)

將式(15)、式(16)代入式(17)中可得

(18)

若式(18)有解，則sinn1π與sinn2π一定同奇偶，即n1=n,n2=n+2。又因為α+β=0，故α=-β，若α=1，則β=-1。得到圓柱殼的有效試函數(shù)是

(n=0,1,2,…)

(19)

將試函數(shù)式(19)代入式(14)，整理得

(20)

根據(jù)第二類拉格朗日方程

(21)

把式(20)代入式(21)，得

(22)

將方程整理得

(23)

其中，

(24)

(25)

(26)

將式(16)分別代入式(24)、式(25)、式(26)得

C1=-ρhRπl(wèi)cr

(27)

(28)

(29)

其中,

(30)

分析式(23)解的性質(zhì)，則有

當C2-C3<0時，得到的解以指數(shù)形式隨時間擴散，表明圓柱殼已經(jīng)發(fā)生屈曲；

當C2-C3>0時，得到的解是隨時間穩(wěn)態(tài)的三角函數(shù)，表明圓柱殼不發(fā)生屈曲；

當C2-C3=0時，處于屈曲臨界狀態(tài)，這時可得到屈曲的臨界條件，其為

(31)

式(31)即為圓柱殼屈曲臨界速度解析表達式。

2 算例分析

在算例分析中，取圓柱殼長L=1 m，半徑R=0.1 m，厚度h=0.002 m，徑厚比k=(2R)/h，E=2×1011Pa，μ=0.3，ρ=7 800 kg/m3。

在m=1,k=200 ,質(zhì)量比為10，n分別為0，1，2的情況下，臨界速度隨臨界長度的變化關系，如圖3所示。

圖3 當n不同時，臨界速度隨臨界長度的變化關系圖

圖3表明：臨界速度隨臨界長度的增大而減小，這也表明應力波效應對動力屈曲有明顯影響；曲線可分為敏感區(qū)和非敏感區(qū)，其分界點約在0.4 m左右，曲線在敏感區(qū)變化比較陡，在非敏感區(qū)變化比較平緩；當臨界長度相同時，臨界速度隨軸向模態(tài)數(shù)增大而增大，說明高速沖擊越容易激發(fā)圓柱殼屈曲的高階模態(tài)。

圖4中的曲線可分為敏感區(qū)和非敏感區(qū)，曲線在敏感區(qū)變化比較陡，在非敏感區(qū)變化比較平緩；臨界速度隨臨界時間的增加而減小，表明應力波效應對動力屈曲有明顯影響。與Karagiozova等研究所得曲線變化規(guī)律、結(jié)論一致。

圖4 臨界速度與臨界時間的變化關系曲線

在n=0,k=200 ,質(zhì)量比為10，m分別為1,2,3的情況下，臨界速度隨臨界長度的變化關系,如圖5所示。

圖5 當m不同時，臨界速度隨臨界長度的變化關系圖

圖5表明：臨界速度隨臨界長度的增大而減小，這也表明應力波效應對動力屈曲有明顯影響；曲線敏感區(qū)和非敏感區(qū)的分界點約在0.3 m左右，結(jié)合圖3，表明：軸向模態(tài)對屈曲的影響比環(huán)向模態(tài)還要明顯；當臨界長度相同時，臨界速度隨環(huán)向模態(tài)數(shù)增大而增大，說明高速沖擊容易激發(fā)圓柱殼屈曲的環(huán)向高階模態(tài)；臨界長度約在0.1左右時，曲線幾乎重合，說明高速下屈曲模態(tài)躍遷不明顯，這與軸向模態(tài)對屈曲的影響有不同之處。

在m=1,k=200 ,lcr=0.4 m,n分別為0,1,2的情況下，臨界速度隨質(zhì)量比的變化關系,如圖6所示。

圖6 當n不同時，臨界速度隨質(zhì)量比的變化關系圖

圖6表明：臨界速度隨質(zhì)量比的增加而減小，表明初始動能對屈曲有明顯影響，這與Karagiozova等研究中數(shù)值模擬所得曲線變化規(guī)律、結(jié)論一致；曲線可分為敏感區(qū)和非敏感區(qū)，其分界點約在質(zhì)量比為5左右；在質(zhì)量比相同時，臨界速度隨軸向模態(tài)數(shù)增大而增大，初始沖擊動能也隨模態(tài)數(shù)增大而增大，進一步說明初始動能對屈曲有明顯影響，同時也說明高速沖擊或較大的初始動能容易激發(fā)圓柱殼屈曲的高階軸向模態(tài)。

在n=0,k=200 ,lcr=0.4 m,m分別為1，2，3的情況下，臨界速度隨質(zhì)量比的變化關系，如圖7所示。

圖7 當m不同時，臨界速度隨質(zhì)量比的變化關系圖

圖7表明:臨界速度隨質(zhì)量比的增大而減小；曲線可分為敏感區(qū)和非敏感區(qū)，其分界點約在質(zhì)量比為5左右；在質(zhì)量比相同時，臨界速度隨環(huán)向模態(tài)增大而增大，初始沖擊動能也隨模態(tài)數(shù)增大而增大，說明初始動能對屈曲有明顯影響，同時也說明高速沖擊或較大的初始動能容易激發(fā)圓柱殼屈曲的高階環(huán)向模態(tài)。

在n=0，m=1，質(zhì)量比為10，改變k，使其分別為100，150，200的情況下，臨界速度隨臨界長度的變化關系,如圖8所示。

圖8表明:臨界速度隨臨界長度的增大而減小，這也表明應力波效應對動力屈曲有明顯影響；曲線敏感區(qū)和非敏感區(qū)的分界點約在lcr=0.4 m左右；在臨界長度相同時，臨界速度的變化隨徑厚比的減小而增大，這說明徑厚比越大，高速沖擊越容易激發(fā)圓柱殼的屈曲。

在n=0，m=1,lcr=0.4 m，改變k,使其分別為100，150，200的情況下，臨界速度隨質(zhì)量比的變化關系,如圖9所示。

圖9 當k不同時，臨界速度隨質(zhì)量比的變化關系圖

圖9表明:臨界速度隨質(zhì)量比的增大而減小；曲線敏感區(qū)和非敏感區(qū)的分界點約在質(zhì)量比為7左右；在質(zhì)量比相同時，臨界速度隨徑厚比的減小而增大，初始沖擊動能也隨徑厚比的減小而增大，這說明徑厚比越大，高速沖擊越容易激發(fā)圓柱殼的屈曲。這與顧紅軍等所得的規(guī)律、結(jié)論一致。

3 結(jié) 論

本文通過理論研究和算例分析得出以下結(jié)論：

(1)基于能量法，考慮應力波效應，建立了拉格朗日函數(shù)，將獲得的滿足邊界條件的準試函數(shù)代入第二類拉格朗日方程，得到圓柱殼屈曲臨界速度的解析表達式。

(2)算例分析，討論了軸向模態(tài)數(shù)、環(huán)向模態(tài)數(shù)、徑厚比、臨界長度、質(zhì)量比對動力屈曲的影響，結(jié)果表明:

當模態(tài)數(shù)不同時，臨界速度和臨界長度的關系曲線可分為較陡的敏感區(qū)和變化平緩非敏感區(qū)，臨界速度隨臨界長度的增加而減小，表明應力波效應對圓柱殼的動力屈曲有明顯影響；臨界速度與質(zhì)量比的關系曲線也可分為敏感區(qū)和非敏感區(qū)，當質(zhì)量比相同時，初始沖擊動能隨模態(tài)數(shù)增大而增大，表明初始動能對屈曲有明顯影響；臨界速度隨模態(tài)數(shù)增大而增大，說明高速沖擊容易激發(fā)圓柱殼屈曲的高階模態(tài)。

當徑厚比不同時，臨界速度隨徑厚比的減小而增大，這說明徑厚比越大，高速沖擊越容易激發(fā)圓柱殼的屈曲。