第一作者繆馥星女，博士，碩士生導(dǎo)師，1980年2月生

郵箱：aifxmiao@gmail.com

非對稱層合梁瞬態(tài)響應(yīng)回傳射線矩陣法解

繆馥星1，孫國鈞2

(1.寧波大學(xué)機(jī)械工程與力學(xué)學(xué)院，浙江寧波315211； 2.上海交通大學(xué)船舶海洋與建筑工程學(xué)院，上海200240)

摘要：用回傳射線矩陣法分析計(jì)及一階剪切變形效應(yīng)的層合復(fù)合材料梁瞬態(tài)動力響應(yīng)；給出速度瞬態(tài)響應(yīng)的回傳射線矩陣法計(jì)算公式；計(jì)算分析軸向-彎曲-剪切耦合效應(yīng)對層合梁瞬態(tài)動力響應(yīng)影響。結(jié)果表明，回傳射線矩陣法可求解計(jì)及一階剪切變形效應(yīng)時(shí)層合梁瞬態(tài)動力響應(yīng)，且與有限元法解一致；橫向(軸向)加載可引起明顯的軸向(橫向)響應(yīng)，隨鋪層角增大軸向-彎曲-剪切耦合效應(yīng)增加，沿不同方向加載下軸向、橫向速度響應(yīng)幅值明顯增加。

關(guān)鍵詞：一階剪切變形效應(yīng)；軸向-彎曲-剪切耦合效應(yīng)；瞬態(tài)動力響應(yīng)；回傳射線矩陣法；層合梁

基金項(xiàng)目：國家自然科學(xué)青年基金(11302109)；爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室開放基金(KFJJ13-7M)；教育部博士點(diǎn)基金新教師類(20133305120001)；浙江省教育廳項(xiàng)目(Y201327004)；寧波市自然科學(xué)基金(2013A610052)；浙江省重中之重學(xué)科項(xiàng)目(zj1221)

收稿日期：2014-06-03修改稿收到日期：2014-09-25

中圖分類號:O.347.4

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2015.20.001

Abstract:The method of reverberation-ray matrix (MRRM) was extended to the transient dynamic response analysis of laminated composite beams based on the first shear deformation theory (FSDT). A theoretical model of laminated composite beam based on the MRRM and first shear deformation theory was established. The velocity responses of a laminated composite cantilever beam under smoothed triangular pulse force load were investigated. The results show their good agreement with those obtained by finite element method (FEM). The effect of ply-stacking sequence on transient dynamic response was also illustrated. It is concluded that the MRRM can solve accurately the dynamic response of laminated composite beams. The axial-bending-shear coupling effects in asymmetrically laminated beams were analyzed by the MRRM. It is found that the axial (transverse) load can cause transverse (axial) dynamic response with the increase of the axial-bending-shear coupling coefficient and the amplitude of the coupled response also increases with the increase of coupling coefficient apparently.

Transient dynamic response analysis of asymmetric laminated composite beam by using the reverberation-ray matrix method

MIAOFu-xing1,SUNGuo-jun2(1. Department of Mechanical Engineering and Mechanics Science, Ningbo University, Ningbo 315211, China；2. Department of Engineering Mechanics, Shanghai Jiaotong University, Shanghai 200240, China)

Key words:first-order shear deformation effect; axial-bending-shear coupling effect; transient dynamic response; reverberation-ray matrix method; laminated composite beam

復(fù)合材料懸臂構(gòu)件大量用于轉(zhuǎn)子系統(tǒng)、嵌入有微型傳動裝置、傳感器及電子器件智能復(fù)合材料裝置。對各向異性梁及結(jié)構(gòu)而言，因其拉伸模量與橫向剪切模量之比較高，剪切效應(yīng)明顯存在，橫向剪切變形不可忽略；另外，鋪層順序也會改變其行為，當(dāng)鋪層關(guān)于中面非對稱時(shí)軸向-彎曲-剪切耦合效應(yīng)明顯存在。因此，用經(jīng)典層合梁理論 (CLBT)[1]進(jìn)行瞬態(tài)動力響應(yīng)分析不夠精確，尤其高階頻率下[2-3]。

計(jì)及剪切變形研究成果較多，主要以有限元法分析為主。Chen等[4]用有限元法分析考慮剪切變形時(shí)層合梁中自由振動頻率。計(jì)及層合復(fù)合材料梁一階剪切變形效應(yīng)及慣性矩時(shí)，Chandrashekhara等[5]計(jì)算分析剪切變形影響依賴梁的長細(xì)比(l/h)表明，對細(xì)梁(l/h>100)，考慮一階剪切變形及旋轉(zhuǎn)矩影響時(shí)動力響應(yīng)與經(jīng)典層合梁理論結(jié)果一致；Abramovich等[6-7]對基于Timoshenko梁方程的對稱層合梁考慮軸向、橫向位移、剪切變形及慣性矩影響，分析認(rèn)為耦合效應(yīng)在復(fù)合材料結(jié)構(gòu)動力分析中不可忽略，并進(jìn)一步討論非對稱鋪層復(fù)合材料梁中的耦合效應(yīng)。Mahapatra等[8]用譜單元法分析探討層合梁中軸向-彎曲-剪切耦合效應(yīng)對其波傳播影響。然而，無論有限元法或譜單元法，在求解系統(tǒng)矩陣方程時(shí)均存在奇異性問題。Pao等[9-11]的回傳射線矩陣分析法用Neumann級數(shù)展開，能有效消除矩陣方程求解中的奇異性問題，且已用實(shí)驗(yàn)證明能精確計(jì)算框架結(jié)構(gòu)的瞬態(tài)動力響應(yīng)?？婐バ堑萚12-14]已將此法擴(kuò)展到基于經(jīng)典層合梁理論的層合復(fù)合材料梁及框架瞬態(tài)動力響應(yīng)分析。本文用此法分析計(jì)及一階剪切變形效應(yīng)的層合復(fù)合材料梁瞬態(tài)動力響應(yīng)。

基于一階剪切變形理論(First-order Shear Deformation Theory, FSDT)建立層合復(fù)合材料梁瞬態(tài)動力響應(yīng)的回傳射線矩陣分析方法，分析非對稱層合復(fù)合材料梁在沖擊載荷作用下的速度瞬態(tài)響應(yīng)，建立層合復(fù)合材料梁的散射、回傳射線矩陣；計(jì)算光滑三角脈沖力作用下懸臂梁速度響應(yīng)，并用文獻(xiàn)[8]有限元結(jié)果進(jìn)行驗(yàn)證；通過計(jì)算不同鋪層角的層合懸臂梁在半正弦波力脈沖下的速度響應(yīng)，分析梁中軸向-彎曲-剪切耦合效應(yīng)及長細(xì)比l/h=100懸臂梁在懸壁端沿軸向、橫向分別施加半正弦脈沖力時(shí)速度響應(yīng)，并與經(jīng)典層合梁理論的軸向、橫向速度瞬態(tài)響應(yīng)分析對比。

1理論分析

1.1控制方程

對圖1層合復(fù)合材料梁，坐標(biāo)系選擇及標(biāo)記方式同文獻(xiàn)[8]。設(shè)層合梁每層均關(guān)于梁材料主軸正交各向異性，據(jù)一階剪切變形層合梁理論，層合梁軸向、橫向位移場可寫為

(1)

線性應(yīng)變?yōu)?/p>

(2)

式中：U，W分別為復(fù)合材料纖維層任一點(diǎn)處軸向、橫向位移；u0，w，z分別為x-y-z坐標(biāo)系下沿中性面軸向、橫向位移及z向坐標(biāo)值；〈·〉,x為x偏微分。

設(shè)層合梁共有n層，第k層厚度用hk(k=1,2,…,n)表示，則層合梁中第k纖維層向本構(gòu)關(guān)系為

(3)

聯(lián)立式(2)、(3)得

(4)

(5)

圖1　層合梁單元 Fig.1 Laminated composite beam element

本文主要分析平面應(yīng)變假設(shè)情形。據(jù)Hamilton原理，計(jì)及一階剪切變形層合梁頻域內(nèi)控制微分方程為

(6)

式中：A11,B11,D11,A55為剛度系數(shù)；I0,I1,I2為慣性矩項(xiàng)相系數(shù)。

頻域內(nèi)軸力、剪力、彎矩可表示為

(7)

與剛度系數(shù)有關(guān)的軸向-彎曲耦合因子r與材料非對稱性相關(guān)的軸向-彎曲耦合因子s2及彎曲-剪切耦合因子s1分別為

(8)

頻域內(nèi)非耦合軸向波、彎曲波、剪切波波數(shù)分別為

(9)

(10)

式中：ca,cb，cs為軸向波、彎曲波、剪切波波速。

1.2一階剪切變形層合梁回傳射線矩陣法公式

(11)

(12)

顯然，αj，βj為剛度系數(shù)、波數(shù)及頻率的函數(shù)。

聯(lián)立式(7)、(11)，軸力、剪力、彎矩可寫為

(13)

式中：Sij(i=1,2,3,j=1,2,…,6, )見文獻(xiàn)[16]。

層合材料懸臂梁見圖2，固定端頻域內(nèi)位移為

(14)

圖2　層合懸臂梁與局部坐標(biāo)系示意圖 Fig.2 The sketch of lanminated composite beam and its local coordinate

(15)

將每個(gè)節(jié)點(diǎn)方程組裝后可得整個(gè)結(jié)構(gòu)的控制方程，推導(dǎo)求解后獲得出射波波幅矢量，表示為

{d(ω)}=[I-R(ω)]-1{s(ω)}

(16)

式中：[R(ω)]為一階剪切變形層合復(fù)合材料梁回傳射線矩陣；s(ω)為源矢量。

將{I-R(ω)}-1展開成Neumann級數(shù)形式為

{I-R}-1=I+R+R2+…+Rn+…

(17)

求出整體系的出射波矢量8ey00ek及入射波矢量{a}后代入控制方程及內(nèi)力方程，可得頻域內(nèi)位移、內(nèi)力等矢量；進(jìn)而求得層合梁相關(guān)物理量的瞬態(tài)響應(yīng)。再用Fourier逆變換獲得其在時(shí)域內(nèi)響應(yīng)。懸臂梁任意位置處橫向加載時(shí)懸臂端頻域軸向速度與橫向速度響應(yīng)分別為

(18)

(19)

2算例

2.1層合懸臂梁速度響應(yīng)

為檢驗(yàn)考慮一階剪切變形效應(yīng)層合梁瞬態(tài)動力響應(yīng)分析的回傳射線矩陣法理論，以AS/3501-6石墨-環(huán)氧纖維層合復(fù)合材料懸臂梁為例,材料特性為：E1= 144.48 GPa，E2=9.632GPa，G12=4.128GPa，ν12=0.3，ρ=1 389 kg/m3；梁長l=1.0 m,度b=0.01m,厚h=0.01 m。選時(shí)間單位t0=13.385 80 μs，長度單位l0=0.1 m。則層合梁無量綱化總厚度、寬度均為0.1，長度為10。梁長細(xì)比l/h=100，計(jì)算非對稱鋪層[05/905](軸向-彎曲-剪切耦合因子r=0.757,s1=3.338×10-4)。將整個(gè)層合梁視為具有2個(gè)節(jié)點(diǎn)、2個(gè)回傳射線矩陣單元、在懸臂端分別沿軸向、橫向施加光滑三角脈沖荷載，見圖3，其中嵌套圖為頻譜圖。懸臂端軸向、橫向加載時(shí)回傳射線矩陣法與有限元法速度瞬態(tài)響應(yīng)比較見圖4。由圖4看出，回傳射線矩陣法(MRRM)軸向、橫向速度瞬態(tài)響應(yīng)與有限元法吻合良好。

圖3 光滑三角力脈沖示意圖Fig.3Impactloadhistoryofthesmoothedtriangularforcepulse.Frequencyspectrumofloadisshownintheinset圖4 一階剪切層合懸臂梁在光滑三角脈沖力作用下響應(yīng)比較Fig.4ComparisonofFSDTlaminacompositecantileverbeambyMRRMandFEM,undersmoothedtriangularforcepulse

圖5　回傳射線矩陣法懸臂端分別在CLBT及FSDT下響應(yīng)比較 Fig.5 Comparison of lamina cantilever plate beam predicted respectively based on CLBT and FSDT by MRRM,under semi-sin force pulse

圖5為平面應(yīng)變假設(shè)下長細(xì)比l/h=100懸臂梁在懸臂端沿軸向、橫向分別施加半正弦脈沖力時(shí)，分別在經(jīng)典層合梁理論與一階剪切變形理論下回傳射線矩陣法兩向速度響應(yīng)比較。

2.2軸向-彎曲-剪切耦合效應(yīng)

為進(jìn)一步分析軸向-彎曲-剪切耦合效應(yīng)，計(jì)算半正弦波力脈沖作用情形。以T300/5208材料層合懸臂梁為例，材料特性見文獻(xiàn)[16]。選三種鋪層[010]、[03/±304/453]、[05/905]，分別在平面應(yīng)變、應(yīng)力假設(shè)下軸向-彎曲耦合因子見表1。3種鋪層下軸向波波速分別為10 659.83 m/s、8 488.68 m/s、7 749.14 m/s。在懸臂端分別沿軸向、橫向施加半正弦波力脈沖，所得懸臂梁軸向、耦合速度瞬態(tài)響應(yīng)及橫向速度瞬態(tài)響應(yīng)分別見圖6。

圖6　一階剪切變形層合懸臂梁懸臂端速度瞬態(tài)響應(yīng) Fig.6 Lamina cantilever plate beam based on FSDT at the cantilever tip due to axial and transverse tip load of semi-sin force pulse, respectively

3結(jié)果分析

由圖4看出，計(jì)及一階剪切變形效應(yīng)時(shí)，層合復(fù)合材料懸臂梁橫向速度響應(yīng)結(jié)果與有限元法一致，尤其軸向速度在響應(yīng)初期幾乎完全一致。一次反射后回傳射線矩陣法幅值結(jié)果大于有限元法。由此，回傳射線矩陣法亦可用于預(yù)測一階剪切變形層合復(fù)合材料梁的瞬態(tài)動力響應(yīng)。

由圖5看出，分別在CLBT及FSDT下的軸向速度響應(yīng)、耦合速度瞬態(tài)響應(yīng)幾乎完全一致，與理論上梁的長細(xì)比較大結(jié)論相符合；而橫向速度瞬態(tài)響應(yīng)初期完全一致，速度振蕩現(xiàn)象在FSDT下出現(xiàn)較晚。反應(yīng)出剪切變形對橫向速度瞬態(tài)響應(yīng)影響較大，計(jì)及剪切變形與否的橫向速度瞬態(tài)響應(yīng)必有差別。即不同理論對層合復(fù)合材料梁瞬態(tài)動力響應(yīng)預(yù)測亦不同。

表1　T300/5208三種鋪層軸向-彎曲-剪切耦合因子

由圖6看出，①隨鋪層角增大、軸向-彎曲-剪切耦合因子增加，不同方向加載下軸向、橫向速度響應(yīng)幅值均增加。對端部，力脈沖作用時(shí)產(chǎn)生明顯響應(yīng)；當(dāng)波傳播到懸臂端時(shí)則發(fā)生反射、散射。此時(shí)端部質(zhì)點(diǎn)仍有速度響應(yīng)，加之反射波引起的響應(yīng)，端部響應(yīng)幅值較初始響應(yīng)會略有增加；但隨多次往復(fù)反射，外激勵(lì)撤消后短時(shí)間內(nèi)響應(yīng)將趨于動態(tài)平衡，幅值不再增加。另外，速度反向時(shí)間增長，尤其軸向加載的速度響應(yīng)，特性更明顯。此因軸向波速隨軸向-彎曲-剪切耦合因子增加明顯減小，即軸向波在梁中傳播所需時(shí)間增加，致使端點(diǎn)處速度反射時(shí)間延長。②由圖6(b)知，當(dāng)軸向-彎曲-剪切耦合因子明顯存在時(shí)，軸向(橫向)加載可引起明顯的橫向(軸向)速度響應(yīng)，且耦合因子越大響應(yīng)幅值越大。③由圖6 (c)知，響應(yīng)早期，橫向加載下速度響應(yīng)隨耦合因子增加響應(yīng)幅值明顯增加，但約0.9 ms后響應(yīng)幅值增加較少。此因相同頻率下隨軸向-彎曲-剪切耦合因子增加彎曲波速增加較少所致。④相同耦合因子下橫向加載的速度響應(yīng)幅值較軸向約大一個(gè)量級。因彎曲波、剪切波的彌散性引起，在響應(yīng)時(shí)間內(nèi)軸向波速遠(yuǎn)大于彎曲波速。⑤對不同力脈沖，在軸向加載下速度響應(yīng)與橫向速度響應(yīng)均明顯不同，橫向加載的軸向速度響應(yīng)亦不同，但其橫向速度響應(yīng)較類似。此因低頻的彎曲波速遠(yuǎn)小于軸向波速，軸向波較彎曲波對外激勵(lì)更敏感。

4結(jié)論

(1)本文用回傳射線矩陣法分析計(jì)及一階剪切變形效應(yīng)的層合復(fù)合材料梁瞬態(tài)動力響應(yīng)。建立基于一階剪切變形理論的回傳射線矩陣法瞬態(tài)動力響應(yīng)計(jì)算公式。用Neumann級數(shù)展開法有效解決層合復(fù)合材料梁瞬態(tài)動力響應(yīng)的頻域矩陣方程求解中奇異性問題。

(2)層合懸臂梁速度瞬態(tài)動力響應(yīng)計(jì)算結(jié)果與有限元法一致。通過分析軸向-彎曲-剪切耦合效應(yīng)知，橫向(軸向)加載可引起明顯的軸向(橫向)響應(yīng)，鋪層角增大軸向-彎曲-剪切耦合效應(yīng)增加，沿不同方向加載的軸向、橫向速度響應(yīng)幅值均增加。

參考文獻(xiàn)

[1]Abramovich H, Livshits A. Free vibration of non-symmetric cross-ply laminated composite beams[J]. Journal of Sound and Vibration,1994,176(5): 597-612.

[2]Obst A W, Kapanin R K. Nonlinear static and transient analysis of laminated composite beams[J]. Journal of Composites Engineering,1992, 2(5/6/7):375-390.

[3]Bhimaraddi A. Computation of ply thickness of laminated beam for which Bernoulli-Euler theory is adequate[J]. Composite Structures,1994, 29(4):415-420.

[4]Chen A T,Yang T Y. Static and dynamic formulation of symmetrically laminated beam finite element for a microcomputer[J]. Journal of Composite Materials,1985, 19(5): 459-475.

[5]Chandrashekhara K, Krishnamurthy K, Roy S. Free vibration of composite beams including rotary inertia and shear deformation[J]. Composite Structures, 1990,14(4): 269-279.

[6]Abramovich H. Shear deformation and rotary inertia effects of vibrating composite beams[J]. Composite Structures, 1992, 20(3):165-173.

[7]Abramovich H,Eisenberger M,Shulepov O. Vibration and bucking of cross-ply nonsymmetric laminated composite beams [J]. Journal of AIAA,1996, 34(5):1064-1069.

[8]Mahapatra D R,Gopalakrishnan S. A spectral finite element modal for analysis of axial-flexural-shear coupled wave propagation in laminated composite beams[J]. Composite Structures, 2003,59(1):67-88.

[9]Pao Y H, Keh D C, Howard S M. Dynamic response and wave propagation in plane trusses and frames [J]. Journal of AIAA,1999, 37(5):594-603.

[10]Pao Y H, Sun G J. Dynamic bending strains in planar trusses with pinned or rigid joints[J]. Journal of Engineering Mechanics (ASCE),2003, 129(3):324-332.

[11]Su X Y, Tian J Y, Pao Y H. Application of the reverberation-ray matrix to the propagation of elastic waves in a layered solid [J]. International Journal of Solids and Structures,2002, 39(21/22):5447-5463.

[12]繆馥星,孫國鈞.非對稱層合復(fù)合材料框架的動態(tài)響應(yīng)[J].上海交通大學(xué)學(xué)報(bào):自然科學(xué)版,2009, 43(10): 1658-1661.

MIAO Fu-xing,SUN Guo-jun.Dynamics response of asymmetrically laminated composite frames[J]. Journal of Shanghai Jiaotong University Science,2009,43(10):1658-1661.

[13]Miao Fu-xing, Sun Guo-jun. Dynamics analysis of asymmetrically laminated composite frames[J]. International Conference on Computational & Experimental Engineering & Sciences (ICCES),2009, 309(1):1-2.

[14]Miao F X, Sun G J,Chen K F.Transient response analysis of balanced laminated composite beams by the method of reverberation-ray matrix[J]. International Journal of Mechanics Science, 2013, 77:121-129.

[15]Reddy J N. Mechanics of laminated composite plates[M]. USA:CRC Press,1997.

[16]繆馥星.框架振動和瞬態(tài)響應(yīng)的回傳波射矩陣法研究[D].上海:上海交通大學(xué),2009.