杜 杰，吳 振(沈陽(yáng)航空航天大學(xué) 航空航天工程學(xué)部(院)，沈陽(yáng) 110136)

機(jī)械與材料工程

熱環(huán)境下復(fù)合材料層合梁自由振動(dòng)分析

杜 杰，吳 振
(沈陽(yáng)航空航天大學(xué) 航空航天工程學(xué)部(院)，沈陽(yáng) 110136)

溫度升高導(dǎo)致彈性常數(shù)下降?；赗eddy理論，研究熱環(huán)境下材料常數(shù)變化對(duì)復(fù)合材料層合梁振動(dòng)特性影響?；诠茴D原理，推導(dǎo)了Reddy層合梁的動(dòng)力平衡方程，得到了層合梁的固有頻率，分析了長(zhǎng)細(xì)比對(duì)固有頻率的影響，對(duì)比了熱環(huán)境下Reddy梁及Euler梁的固有頻率。數(shù)值結(jié)果表明隨溫度升高，彈性常數(shù)下降，Reddy梁對(duì)應(yīng)的固有頻率相應(yīng)降低，而Euler梁的固有頻率無(wú)明顯變化。

熱環(huán)境；材料常數(shù)；Reddy理論；層合梁；自由振動(dòng)

由于復(fù)合材料具有高比強(qiáng)度、高比剛度、低密度及良好的抗疲勞特性，復(fù)合材料層合結(jié)構(gòu)已普遍應(yīng)用于航空航天工程中。航空航天工程中的復(fù)合材料通常處于變化溫度環(huán)境，研究其動(dòng)力特性須考慮到溫度對(duì)材料彈性常數(shù)影響。

基于一階剪切變形理論，Kahya等[1]使用五節(jié)點(diǎn)單元分析了功能梯度梁的屈曲及振動(dòng)特性。Sahu等基于一階剪切變形理論[2]分析了濕熱對(duì)編織纖維層合板振動(dòng)特性的影響。Thai等[3]指出一階理論適用于中厚度梁自由振動(dòng)分析，然而分析厚梁自由振動(dòng)問(wèn)題精度較低。Reddy[4-5]發(fā)展了三階理論并建立了各向同性梁的振動(dòng)方程，此理論不需要使用剪切修正系數(shù)且未知變量個(gè)數(shù)與一階剪切變形理論的未知變量個(gè)數(shù)相同。Jin等[6]基于Reddy理論分析了粘彈性?shī)A層梁的阻尼及振動(dòng)特性。Giunta等[7]分別基于經(jīng)典理論、一階理論及高階理論分析了層合梁的振動(dòng)特性。Nguyen等[8]基于增強(qiáng)型高階理論分析了功能梯度夾層板的彎曲、振動(dòng)及屈曲特性。Ebrahimi等[9]在濕熱環(huán)境下基于非局部應(yīng)變梯度理論分析了粘彈性功能梯度納米梁的振動(dòng)特性。Singh等[10]基于C0型高階鋸齒理論，使用九節(jié)點(diǎn)單元分析了熱機(jī)載荷作用下復(fù)合材料層合/夾層板的屈曲及振動(dòng)特性。李映輝等[11-12]在濕熱環(huán)境下分別分析了熱變形及彈性常數(shù)的變化對(duì)復(fù)合材料薄壁梁及旋轉(zhuǎn)梁振動(dòng)特性影響。Verette等[13]分析了溫度對(duì)碳/環(huán)氧材料彈性常數(shù)的影響。然而，前期工作基于Reddy理論研究復(fù)合材料層合結(jié)構(gòu)動(dòng)力問(wèn)題并沒(méi)有考慮不同溫度下彈性常數(shù)的改變對(duì)振動(dòng)頻率的影響。本文對(duì)不同溫度下Reddy層合簡(jiǎn)支梁的自由振動(dòng)進(jìn)行分析，結(jié)果表明隨著溫度升高材料的彈性常數(shù)下降，Reddy梁固有頻率會(huì)隨之降低。隨著頻率階數(shù)升高，頻率降低會(huì)更為明顯,然而Euler梁固有頻率無(wú)明顯變化。

1 理論公式

1.1 幾何方程

本文模型采用Reddy層合梁理論，其位移場(chǎng)為

u(x,t)=u0(x,t)+φ1u1(x,t)+φ2u3(x,t)

(1)

其中，u1(x,t)=-Ψ(x,t)，u3(x,t)=?w/?x，φ1=z-cz3，φ2=-cz3。式中Ψ為梁繞Y軸截面的轉(zhuǎn)角。其中c=4μ/3h2，當(dāng)μ=1時(shí)，該理論即為Reddy梁理論，當(dāng)μ=0時(shí)則退化為僅考慮一階剪切變形的Timoshenko梁理論，并且當(dāng)Ψ=-?w/?x時(shí)，該理論繼續(xù)退化為Euler-Bernoulli梁理論。

該理論的應(yīng)變表達(dá)式為

(2)

(3)

1.2 本構(gòu)方程

(4)

Q11=cos4φC11+sin4φC22+2cos2φsin2φ(C12+2C55)

(5)

Q44=cos2φC44+sin2φC55+2cos2φsin2φ(C12+2C55)

(6)

其中，φ為鋪設(shè)角。

1.3 動(dòng)力平衡方程和邊界條件

本文基于哈密頓原理推導(dǎo)了Reddy簡(jiǎn)支層合梁的動(dòng)力平衡方程，該原理可表述為

(7)

其中U為應(yīng)變能，W為外力做的功，K為動(dòng)能。其中，W的一階變分為

(8)

f和q分別為沿梁X軸和Y軸的體力。

U的一階變分為

(9)

K的一階變分為

(10)

將式 (8)、(9)、(10)帶入式 (7)即得Reddy梁的動(dòng)力平衡方程

(11)

(12)

(13)

滿足邊界條件的位移試函數(shù)為

(14)

其中u0n、u1n、wn為傅里葉系數(shù)，n為頻率的階數(shù)。

將式(11)、(12)、(13)帶入式(14)中可得

(15)

其中

2 算例

2.1 Reddy梁理論有效性驗(yàn)證

表1 梁無(wú)量綱化頻率比較

將本文結(jié)果與精確解(Exact)[14]，整體-局部高階理論有限元計(jì)算結(jié)果(HGLT)[15]及一階理論計(jì)算結(jié)果(FSDT)[14]進(jìn)行比較。表1表明，本文結(jié)果與精確解及整體-局部高階理論有限元計(jì)算結(jié)果符合良好，然而，一階理論精度較低。

2.2 考慮溫度變化的Reddy梁振動(dòng)分析

為分析溫度升高導(dǎo)致的彈性常數(shù)下降對(duì)層合梁振動(dòng)頻率的影響，以簡(jiǎn)支梁為例，分別基于Reddy梁理論及Euler梁理論，分析不同溫度下梁的前4階固有頻率，并討論長(zhǎng)細(xì)比在溫度變化時(shí)對(duì)梁固有頻率的影響。

層合梁彈性常數(shù)隨溫度的變化[11]如表2所示。

其中，G12=G13;G23=0.5G12;v12=0.3;ρ=1 600 kg/m3。

2.2.1 不同溫度下層合梁的固有頻率

表2 不同溫度下的彈性常數(shù)

表3 不同溫度下梁的前4階頻率

表3表明溫度升高時(shí)，Reddy梁各階頻率均會(huì)下降。圖1表明溫度升高時(shí)，高階振動(dòng)比低階振動(dòng)的頻率下降更快，即彈性常數(shù)下降對(duì)高階振動(dòng)影響更大，然而，基于Euler理論，溫度變化對(duì)頻率的影響并不明顯。

圖1 不同階數(shù)下頻率下降的程度

2.2.2 長(zhǎng)細(xì)比在溫度升高時(shí)對(duì)梁固有頻率的影響

表4表明Reddy梁頻率隨長(zhǎng)細(xì)比增加而增加,隨溫度升高而降低。圖2表明基于Reddy理論，當(dāng)長(zhǎng)細(xì)比較小時(shí)，溫度升高引起的彈性常數(shù)下降對(duì)頻率影響較大；隨著長(zhǎng)細(xì)比增加，彈性常數(shù)的下降對(duì)頻率影響逐漸減弱，然而，基于Euler理論，溫度變化對(duì)頻率的影響并不明顯。

表4 不同長(zhǎng)細(xì)比的梁在不同溫度下的固有頻率

圖2 長(zhǎng)細(xì)比對(duì)頻率下降程度的影響

3 結(jié)論

本文通過(guò)哈密頓原理，建立了Reddy梁的動(dòng)力平衡方程，得出了不同溫度下Reddy梁的固有頻率并與Euler梁的固有頻率進(jìn)行了比較。隨著溫度升高，材料彈性常數(shù)下降。數(shù)值結(jié)果表明：

(1)彈性常數(shù)下降會(huì)導(dǎo)致Reddy梁固有頻率降低，同時(shí)，隨著模態(tài)階數(shù)增加，其影響會(huì)更為明顯。然而，溫度升高導(dǎo)致的彈性常數(shù)下降對(duì)Euler梁固有頻率無(wú)明顯影響。

(2)長(zhǎng)細(xì)比越大，Reddy梁與Euler梁的固有頻率就越接近，即溫度引起的材料彈性常數(shù)改變對(duì)Reddy梁固有頻率的影響在長(zhǎng)細(xì)比較小時(shí)會(huì)更為明顯，當(dāng)長(zhǎng)細(xì)比接近40時(shí)，彈性常數(shù)的改變對(duì)固有頻率影響并不明顯。

[1] KAHYA V,TURAN M.Finite element model for vibration and buckling of functionally graded beams based on the first-order shear deformation theory[J].Composites Part B Engineering,2017,109(10):108-115.

[2] PANDA H S,SAHU S K,PARHI P K.Hygrothermal effects on free vibration of delaminated woven fiber composite plates-Numerical and experimental results[J].Composite Structures,2013,96(1):502-513.

[3] THAI C H,TRAN L V,TRAN D T,et al.Analysis of laminated composite plates using higher-order shear deformation plate theory and node-based smoothed discrete shear gap method[J].Applied Mathematical Modelling,2012,36(11):5657-5677.

[4] REDDY J N.A Simple Higher-order theory for laminated composite Plates[J].Appl Mech,1984(51):745-752.

[5] REDDY J N，HEYLIGER P R.A higher order beam finite element for bending and vibration problems[J].Journal of Sound and Vibration,1988,126(2):309-326.

[6] JIN G，YANG C，LIU Z.Vibration and damping analysis of sandwich viscoelastic-core beam using Reddy'shigher-order theory[J].Composite Structures,2016,140(1):390-409.

[7] GIUNTA G,BISCANI F,BELOUETTAR S,et al.Free vibration analysis of composite beams via refined theories[J].Composites Part B Engineering,2013,44(1):540-552.

[8] NGUYEN K T,TAI H T,VO T P.A refined higher-order shear deformation theory for bending,vibration and buckling analysis of functionally graded sandwich plates[J].Steel and Composite Structures,2015,18(1):91-120.

[9] EBRAHIMI F，BARATI MR.Hygrothermal effects on vibration characteristics of viscoelastic FG nanobeams based on nonlocal strain gradient theory[J].Composite Structures,2017,159(9):433-444.

[10]SINGH S K,CHAKRABARTI A.Static,vibration and buckling analysis of skew composite and sandwich plates under thermo mechanical loading.International Journal of Applied Mechanics and Engineering,2013,18(3):887-898.

[11]馬艷龍,李映輝.濕熱環(huán)境下復(fù)合材料薄壁梁振動(dòng)特性研究[J].振動(dòng)與沖擊,2016,35(15):154-160.

[12]蔣寶坤,李映輝.濕熱環(huán)境對(duì)旋轉(zhuǎn)復(fù)合材料梁擺振特性的影響[J].復(fù)合材料學(xué)報(bào),2015,32(2):579-585.

[13]VERETTE R M.Temperature/Humidity Effects on the Strength of Graphite/Epoxy Laminates[J].Journal of Aircraft,2015,14(1):90-94.

[14]KAPURIA S，DUMIR P C,JAIN N K.Assessment of zigzag theory for static loading,buckling,free and forced response of composite and sandwich beams[J].Composite Structures,2004,64(3-4):317-327.

[15]WU Z,CHEN W J.Free vibration of laminated composite and sandwich plates using global-local higher-order theory[J].Journal of Sound and Vibration,2006,298(1-2):333-349.

Freevibrationanalysisoflaminatedcompositebeamunderthermalenvironment

DU Jie,WU Zhen
(Faculty of Aerospace Engineering,Shenyang Aerospace University,Shenyang 110136,China)

The increase of temperature causes the decrease of material constants.Based on the Reddy's theory,the influences of material constants under thermal environment on vibration characteristics of a laminated composite beam were investigated.The equations of motion were deduced using Hamilton′s principle to obtain natural frequency of the beam and to analyze the effects of length-thickness ratios.The natural frequency of Reddy beam was compared with that of Euler beam under the thermal environment.Numerical results show that elastic constant decrease and natural frequencies obtained from the Reddy′s model also obviously decrease with the increase of temperature,while the natural frequencies calculated from Euler's theory slightly decrease.

thermal environment;material constants;Reddy′s theory;laminated composite beam;free vibration

2017-09-26

國(guó)家自然科學(xué)基金(項(xiàng)目編號(hào)：11402152)

杜 杰(1990-)，男，山西晉中人，碩士研究生，主要研究方向：層合梁熱振動(dòng)，E-mail：1165923722@qq.com;吳 振(1977-)，男，黑龍江佳木斯人，教授，博士，主要研究方向：復(fù)合材料力學(xué)研究，E-mail：wuzhenhk@163.com。

2095-1248(2017)06-0055-05

O326

A

10.3969/j.issn.2095-1248.2017.06.009

吳萍 英文審校：趙歡)