王宇 郝育新

(北京信息科技大學(xué)機電學(xué)院,北京 100192)

引言

TC4鈦合金材料的組成為Ti-6Al-4V,具有優(yōu)良的耐腐蝕性,密度較小而比強度較高,同時具有良好的韌性和焊接性能等一系列優(yōu)點.該材料存在于汽車船舶、航空航天、石油化工、醫(yī)藥等各行業(yè)領(lǐng)域中.

懸臂板結(jié)構(gòu)是在工程中廣泛應(yīng)用的一種結(jié)構(gòu),也是動力學(xué)研究中頗為常見的一種研究模型.Tao[1]研究了Mindlin懸臂板的振動控制,提出了獨立模態(tài)空間控制和波控制,并應(yīng)用于抑制振動.呂書鋒[2]研究了橫向氣動載荷和參數(shù)激勵聯(lián)合作用下復(fù)合材料懸臂外伸矩形板在伸出過程中的非線性動力學(xué)問題.Li[3]采用Theodorsen機翼理論研究了懸臂板的模態(tài)耦合,并探究了各種實際因素對模態(tài)的細微影響.孫曉婷[4]研究了懸臂板的動力學(xué)特性,采用康特洛維奇法求解了懸臂板的固有特性.劉利軍[5]對轉(zhuǎn)動懸臂板結(jié)構(gòu)的動力學(xué)建模及發(fā)展過程進行了綜述總結(jié),研究了繞定軸轉(zhuǎn)動的剛—柔耦合懸臂板結(jié)構(gòu)的動態(tài)及振動特性.杜長城[6]研究了功能梯度矩形薄板的非線性動力學(xué)問題,分析了動力學(xué)系統(tǒng)的周期振動狀態(tài).

在板材加工和運輸過程中會因各種情況產(chǎn)生初始幾何缺陷,使得懸臂板產(chǎn)生各種細微的變化,但會對其動力學(xué)行為產(chǎn)生較為明顯的影響和變化,因此國內(nèi)外學(xué)者對于模型含幾何缺陷的動力學(xué)研究也比較多.Ankit[7]基于高階剪切變形理論對含有初始幾何缺陷梯度板進行了數(shù)值研究.Rafee[8]以正、余弦以及雙曲函數(shù)模擬了初始幾何缺陷,研究了含初始幾何缺陷的壓電功能梯度材料碳納米管增強層合板的非線性動力學(xué)穩(wěn)定性.Bruno[9]采用一種均一法,對含初始缺陷的夾層板進行了幾何非線性分析.Li[10]研究了含初始幾何缺陷的考慮剪切變形的各向異性梁的屈曲與后屈曲行為.Kirichenko[11]研究了含初始缺陷的熱彈性多層淺圓柱殼在非經(jīng)典數(shù)學(xué)模型下的耦合問題.王珂晟[12]研究了含初始缺陷的夾層圓柱殼的非線性屈曲問題,得到了屈曲載荷、缺陷幅值、缺陷波數(shù)、夾心模量等參量之間的關(guān)系.Chen[13]研究了含初始應(yīng)力和幾何缺陷的功能梯度板的非線性振動問題,得到了初始應(yīng)力、幾何缺陷和體積分數(shù)指數(shù)對非線性振動響應(yīng)的變化影響.

在諸多有關(guān)幾何缺陷的動力學(xué)研究中,對于固有頻率的研究少有涉及.基于經(jīng)典板殼理論,預(yù)設(shè)懸臂板的模態(tài)函數(shù)[14],采用Rayleigh-Ritz法求得懸臂板固有頻率.根據(jù)模擬實際情況中產(chǎn)生的缺陷所得到的擬合缺陷函數(shù)和三角函數(shù)來表示缺陷形式,加入動力學(xué)模型中,得到不同類型缺陷的板頻率.

1 動力學(xué)模型

1.1 動力學(xué)方程

如圖1所示的均質(zhì)懸臂薄版模型,本文研究該模型無阻尼自由振動情況下的前4階頻率.

圖1 均質(zhì)懸臂薄板Fig.1 Homogeneous cantilever thin plate

本文模型為各向同性的均質(zhì)材料懸臂板,采用彈性薄板小撓度理論基本假設(shè),撓度(w)沿著板的厚度方向的分量可以略去,同厚度中各處撓度等于中面的撓度.動力學(xué)方程具體如下：

位移分量[14]：

(1)

式中w2為缺陷項.

應(yīng)變分量[14]：

(2)

本構(gòu)關(guān)系[14]：

(3)

(4a～c)

式中E為彈性模量,μ為泊松比.

能量方程[14]：

變形能U、動能T分別如下：

(5a)

(5b)

1.2 缺陷形式

(A)全局型缺陷：

(6)

(B)局部缺陷：

(7a)

(7b)

公式(7a)表示自由-自由方向幾何缺陷,與x坐標相關(guān),(7b)表示懸臂方向幾何缺陷,與y坐標相關(guān).這兩種缺陷函數(shù)由實際可能出現(xiàn)的缺陷形狀擬合所得.

η為缺陷尺度的系數(shù),控制缺陷整體大小程度；δn控制半波數(shù)或缺陷的密集程度；ψn控制缺陷位置.全局型缺陷均布于整個板面,故不做不同ψn值的對比.

懸臂方向的局部缺陷與y相關(guān),自由方向的局部缺陷與x相關(guān).ψ1和ψ2分別控制自由、懸臂方向局部缺陷的中心位置,δ1和δ2控制局部缺陷的集中程度和影響范圍.

2 固有頻率求解

固有頻率也稱自然頻率(natural frequency),即物體自由振動時的振動頻率.對于板殼等結(jié)構(gòu)的固有頻率進行的研究有著廣泛而深遠的實用性.一方面,對于結(jié)構(gòu)安全、結(jié)構(gòu)使用壽命等實際應(yīng)用方面有著重要的意義,而且對含幾何缺陷的結(jié)構(gòu)進行研究更是具有實用意義；另一方面,對板殼的模態(tài)分析的成果和結(jié)論,可以為進一步分析求解其更多動力學(xué)特性提供有意義的參照.根據(jù)所研究對象的特點,本文采用Rayleigh-Ritz法對懸臂薄板的各階固有頻率的進行求解.

2.1 振型函數(shù)

此處w1(x,y,t)為假設(shè)的模態(tài)函數(shù),滿足懸臂邊界條件.Cmn為待定系數(shù),通過選擇該系數(shù)使得Umax-Tmax的值最小.

w1(x,y,t)=w0(x,y)sin(ωt)

(8)

(9)

根據(jù)懸臂邊界條件,即：x=0固支,x=a和y=0,b自由,振型函數(shù)形式如下：

Xi(x)=sinλix-sinhλix+αi(coshλix-cosλix)

(10)

Y1=1

(11a)

(11b)

Yj(y)= sinμjy+sinhμjy-βj(coshμjy+cosμjy)

(j>2)

(11c)

cosλia·coshλia+1=0

cosμjb·coshμjb-1=0

(11d)

式中：

λ1=1.875,λ2=4.694,λ3=7.854…

μ1=4.730,μ2=7.853,μ3=10.995…

由方程(1)得到加入缺陷項的位移形式,再將方程(10)、(11)帶入方程(9)中即可得到具體的振型函數(shù)；再將所得函數(shù)帶入(2)得到應(yīng)變方程；根據(jù)本關(guān)系(3)、(4)得到應(yīng)力方程.

2.2 求解

由所得應(yīng)力方程、應(yīng)變方程和位移形式(1)帶入方程(5a)、(5b)可得勢能、動能最大值表達式分別為：

(13a)

(w0dt)2]dV

(13b)

圓頻率的求解：

(14)

將方程(9)帶入方程(13)中,得到勢能、動能最大值的展開式；再將方程(13)帶入方程(14)中,得到含m×n個關(guān)于Cmn和的齊次方程的方程組,將Cmn作為自變量提取所得到的方程組的系數(shù)矩陣Amn(只包含常數(shù)和Cmn及w),求解|Amn|=0,即可得到各階圓頻率的值.

3 對比驗證與結(jié)果分析

3.1 對比結(jié)果

為了驗證本文算法的可靠性和準確性,取孫曉婷[4]中的算例進行對比驗證,并用Ansys進行有限元分析,其幾何物理參數(shù)如下：

a= 0.133 m,b=0.109 m,h=0.00135 m,

E=1.138×1011Pa,μ= 0.31,ρ= 4420 kg/m3

結(jié)果如表1所示,前兩階固有頻率值優(yōu)于參考文獻[4],三、四階誤差率略大于前兩階頻率.

表1 固有頻率與其他理論及ANSYS有限元分析的對比Table 1 Comparisons of natural frequencies of cantilever platefrom other theories and ANSYS finite element analysis

另外將本文模型利用ANSYS有限元分析結(jié)果求得前四階固有頻率,與本文算法另做對比,本文所用參數(shù)尺寸：

a= 0.266 m,b=0.209 m,h=0.002 m,

E=1.138×1011Pa,μ= 0.31,ρ= 4420 kg/m3

結(jié)果如表2所示由上可知本文算法比較準確可靠.

表2 懸臂板固有頻率與ANSYS有限元分析的對比Table 2 Comparisons of natural frequencies of cantilever plate from R-R method and ANSYS finite element analysis

3.2 固有頻率

(1)全局型缺陷(η=0.1,ψ0=0.5)

各階固有頻率隨著δ0值的增大而明顯變大,如圖2(a)、(b),δ0分別為1、3,可見隨著半波數(shù)δ0值的增大,缺陷越來越密集.

圖2 全局型幾何缺陷Fig.2 Global geometrical imperfections

如圖3所示.由此可知,懸臂板受缺陷影響隨著缺陷密度集中程度的增長而變大,即缺陷越為集中,對于振動中的板材在偏離平衡位置時的回彈起到積極作用,故頻率變大.

圖3 含全局型幾何缺陷的懸臂板的固有頻率Fig.3 Natural frequencies of cantilever plate with global geometrical imperfection

(2)局部缺陷(η=0.1)

自由-自由方向缺陷如圖4所示,(a)為ψ1=0.133,δ1=0.02；(b)為ψ1=0.133,δ1=0.01.

懸臂方向缺陷如圖5所示,(a)為ψ2=0.133,δ2=0.01；(b)為ψ2=0.133,δ2=0.01.

圖5 懸臂方向局部幾何缺陷Fig.5 Localized geometrical imperfections on the cantilever direction

同全局型缺陷一樣,各階固有頻率隨著δ1值的減小而明顯增大,如表3所示.故可認為缺陷越為集中,懸臂板的固有頻率數(shù)值越大.

另一方面,表3對比了不同方向各階固有頻率數(shù)值,由表3可見在懸臂方向除基礎(chǔ)頻率外各階頻率都大于自由-自由方向.綜合圖6所示的前四階模態(tài)的振型：圖(a)、(b)、(c)、(d)分別為無缺陷懸臂薄版的前四階固有頻率.

表3 含局部幾何缺陷的懸臂板的固有頻率(ψ1、ψ2=0.133)Table 3 Natural frequencies of cantilever plate with localized geometrical imperfections (ψ1、ψ2=0.133)

圖6 無幾何缺陷懸臂薄版ANSYS有限元分析云圖Fig.6 Stress and deformation nephogram of cantilever thin plate withoutgeometrical imperfection through ANSYS finite element analysis

由無缺陷懸臂板的振型可看出,在存在懸臂方向缺陷時,包括彎曲變形的自由振動(即二、三、四階振型)對應(yīng)的固有頻率的變化相對自由-自由方向缺陷影響更大；另一方面,自由方向的缺陷對基礎(chǔ)自由振動的影響比懸臂方向缺陷的影響要大很多.

(3)局部缺陷的位置(η=0.1,δ1=δ2=0.02)

如圖7,(a)、(b)為不同位置的自由-自由方向的幾何缺陷；(c)、(d)為不同位置的懸臂方向的幾何缺陷.

圖7 不同位置的幾何缺陷Fig.7 Geometrical imperfections on different locations

由表4可以看出,懸臂板局部缺陷的位置對固有頻率的影響主要體現(xiàn)在懸臂方向的缺陷上,自由方向的缺陷位置對頻率的影響較??；從懸臂方向不同位置缺陷所得固有頻率可以看出,位置對頻率的影響趨勢也隨振型相關(guān).

表4 含局部不同位置的幾何缺陷的懸臂板的固有頻率Table 4 Natural frequencies of cantilever plate with localized geometrical imperfections on different locations

(4)缺陷尺寸

由圖8數(shù)據(jù)可以明顯看出,懸臂板的初始幾何缺陷尺寸越大(即缺陷程度越大),對于固有頻率的影響也越大.

圖8 不同缺陷尺寸(η=0.3/0.5/0.7)的懸臂板的固有頻率Fig.8 Natural frequencies of cantilever plate with different sizes of imperfections (η=0.3/0.5/0.7)

4 結(jié)論

本文基于經(jīng)典板殼理論,彈性薄板小撓度理論基本假設(shè),預(yù)設(shè)懸臂板的振型函數(shù),采用Rayleigh-Ritz法求得懸臂薄板前四階固有頻率；采用余弦函數(shù)、自擬合函數(shù)描述不同類型缺陷帶入懸臂板的動力學(xué)模型求解,得到不同缺陷影響下的懸臂薄板的前四階固有頻率.從計算結(jié)果得到以下結(jié)論:

(1)初始幾何缺陷的存在會增大均質(zhì)材料懸臂薄板的固有頻率.

(2)懸臂板受缺陷影響隨著缺陷密度集中程度的增長而變大,即缺陷越為集中,對于振動中的矩形板結(jié)構(gòu)在偏離平衡位置時的回彈起到的積極作用越明顯,故頻率變大.

(3)懸臂板的初始幾何缺陷尺寸越大(即缺陷程度越大),其各階固有頻率也越大,缺陷的幾何尺寸對自由振動中的板的回彈同樣起到積極作用.

(4)局部缺陷的不同位置、不同方向?qū)冶郯宓淖杂烧駝拥挠绊懪c其各階具體振型有關(guān):

a.缺陷方向：自由方向缺陷對基礎(chǔ)頻率影響較大,而懸臂方向缺陷對高階頻率影響較大.

b.缺陷位置：自由方向缺陷的不同位置對頻率影響較小,懸臂方向的不同位置對頻率影響較大.

1Chen T, Hu C, Huang W H. Vibration control of cantilevered Mindlin-type plates.JournalofSoundandVibration, 2008.07.028

2呂書鋒,張偉. 復(fù)合材料懸臂外伸板的費線性動力學(xué)建模及數(shù)值研究. 動力學(xué)與控制學(xué)報, 2015,13(4):288～292 (Lu S F, Zhang W. Nonlinear analysis of deploying laminated composite cantilever plates.JournalofDynamicsandControl, 2015,13(4):288～292 (in Chinese))

3Huang L X, Zhang C. Modal analysis of cantilever plate flutter.JournalofFluidsandStructures, 2013,38(1):273～289

4孫曉婷. 懸臂板的動力學(xué)特性分析及減震研究[碩士學(xué)位論文]. 大連：大連理工大學(xué), 2015:09～14 (Sun X T. Analysis of dynamic characteristics an research on vibration suppression of cantilever plate[Master Thesis]. Dalian: Dalian University of Technology, 2015:09～14 (in Chinese))

5劉利軍,張志誼,杜冬. 作定軸轉(zhuǎn)動懸臂板的動力學(xué)建模及其動態(tài)特性分析. 機械科學(xué)與技術(shù), 2009,28(9):1153～1156 (Liu L J, Zhang Z Y, Du D. Dynamic modeling and modal analysis of a cantilever plate rotating around a rigid hub.MechanicalScienceandTechnologyforAerospaceEngineering, 2009,28(9):1153～1156 (in Chinese))

6杜長誠,李映輝. 功能梯度矩形板的非線性自由振動. 力學(xué)季刊, 2010,31(2):250～255 (Du C C, Li Y H. Nonlinear free vibration of functionally graded rectangular plates.ChineseQuarterlyofMechanics, 2010,31(2):250～255 (in Chinese))

7Gupta A, Talha M. An assessment of a non-polynomial based higher order shear and normal deformation theory for vibration response of gradient plates with initial geometric imperfections.CompositesPartB, 2016,107(1):141～161

8Rafidee M, He X Q, Liew K M. Non-linear dynamic stability of piezoelectric functionally graded carbon nanotube-reinforced composite plates with initial geometric imperfection.InternationalJournalofNon-LinearMechanics, 2014,59(1):37～51

9Goncalves B R, Jelovica J, Romanoff J. A homogenization method for geometric nonlinear analysis of sandwich structures with initial imperfections.InternationalJournalofSolidsandStructures, 2016,87(1):194～205

10 Li Z M, Q P Z. Buckling and postbuckling behavior of shear deformable anisotropic laminated beams with initial geometric imperfections subjected to axial compression.EngineeringStructures, 2015,85(1):277～292

11 Kirichenko V F, Awrejcewicz J, Kirichenko A V, et al. On the non-classical mathematical models of coupled problems of thermo-elasticity formulti-layer shallow shells with initial imperfections.InternationalJournalofNon-LinearMechanics, 2015,74(1):51～72

12 王珂晟,唐國金. 含軸對稱初始缺陷的具有正交各向異性表層的夾層圓柱殼的非線性屈曲. 力學(xué)季刊, 2003,24(4):560～566 (Wang K S, Tang G J. Nonlinear buckling of sandwich cylindrical shell with orthotropic facings containing axisymmetric shape initial imperfection.ChineseQuarterlyofMechanics, 2003,24(4):560～566 (in Chinese))

13 Chen C S, Tan A H. Imperfection sensitivity in the nonlinear vibration of initially stresses functionally graded plates.CompositeStructures, 2007,78(1):529～536

14 Reddy J N. Mechanics of laminated composite plates and shells. the United States: Crc Press, 2004