尹濤 王震威

摘要： 脫層是復(fù)合材料層合結(jié)構(gòu)最主要的破壞形式之一。基于EulerBernoulli梁理論和譜有限元方法得到周期結(jié)構(gòu)中健康基本周期單元的動(dòng)剛度矩陣。再考慮脫層位于周期結(jié)構(gòu)的某一基本周期單元中，并假設(shè)脫層邊緣處橫截面變形服從平截面假定，忽略脫層前沿的應(yīng)力奇異，建立分層模型，根據(jù)分層邊緣處各子單元位移、轉(zhuǎn)角連續(xù)及分層界面處力平衡條件形成含脫層的周期層合梁超單元?jiǎng)觿偠染仃嚒；趥鬟f矩陣法得到健康及含脫層周期單元的傳遞矩陣，進(jìn)而形成含脫層周期層合梁結(jié)構(gòu)的傳遞矩陣及總體動(dòng)剛度矩陣，并通過(guò)傳遞矩陣計(jì)算獲得的波傳播常數(shù)來(lái)分析周期結(jié)構(gòu)的振動(dòng)與波傳播特性。分別對(duì)無(wú)脫層情況與脫層引起局部失諧情況下一周期簡(jiǎn)支梁結(jié)構(gòu)的波傳播特性進(jìn)行數(shù)值計(jì)算研究，并結(jié)合ANSYS有限元仿真分析，對(duì)方法進(jìn)行了驗(yàn)證。

關(guān)鍵詞： 周期結(jié)構(gòu)； 層合梁； 脫層； 譜有限元法； 傳播常數(shù)

中圖分類(lèi)號(hào)： O327； O343文獻(xiàn)標(biāo)志碼： A文章編號(hào)： 10044523（2016）04056809

DOI：10.16385/j.cnki.issn.10044523.2016.04.002

引言

近年來(lái)，復(fù)合材料由于其強(qiáng)度高、自重輕等優(yōu)點(diǎn)在航空航天、機(jī)械、船舶、生物工程以及土木工程等領(lǐng)域都得到越來(lái)越廣泛的應(yīng)用。而對(duì)于復(fù)合材料層合結(jié)構(gòu)而言，脫層的出現(xiàn)能夠明顯改變其局部以及整體結(jié)構(gòu)的動(dòng)力特性，對(duì)其進(jìn)行基于振動(dòng)與波傳播理論的動(dòng)力特性研究就顯得尤其重要，這同時(shí)也為脫層缺陷檢測(cè)提供理論依據(jù)。目前國(guó)內(nèi)外已有學(xué)者對(duì)含脫層復(fù)合材料層合結(jié)構(gòu)的動(dòng)力特性進(jìn)行研究。如，王德明和錢(qián)管良分別通過(guò)模態(tài)分析法研究了含裂縫各向同性梁的振動(dòng)問(wèn)題[12]；Wang等人對(duì)脫層梁的振動(dòng)問(wèn)題進(jìn)行了深入研究，并證明當(dāng)脫層較短時(shí)，梁的固有頻率不會(huì)發(fā)生顯著變化[3]；Zhang從應(yīng)力波的頻散及反射現(xiàn)象角度對(duì)結(jié)構(gòu)內(nèi)部缺陷檢測(cè)進(jìn)行了較詳細(xì)研究[4]；羅松南等對(duì)含有任意脫層復(fù)合材料梁中脫層位置及其長(zhǎng)度對(duì)復(fù)合材料梁動(dòng)力特性的影響進(jìn)行了研究[5] 。

另一方面，周期結(jié)構(gòu)在工程領(lǐng)域中獲得了廣泛的應(yīng)用，例如聲子晶體、超高層剪切型建筑、超長(zhǎng)石油天然氣管道、多跨高架橋、航天器太陽(yáng)能電池帆板、渦輪葉片以及復(fù)合材料層合結(jié)構(gòu)等。從20世紀(jì)60年代開(kāi)始，國(guó)內(nèi)外學(xué)者對(duì)周期結(jié)構(gòu)進(jìn)行了較深入地研究，并發(fā)現(xiàn)周期結(jié)構(gòu)不同于其他結(jié)構(gòu)的一些特性：周期結(jié)構(gòu)中會(huì)產(chǎn)生特有的帶隙（即頻率通帶和阻帶）現(xiàn)象[68]，其表現(xiàn)在當(dāng)波處在結(jié)構(gòu)的頻率通帶范圍內(nèi)時(shí)，波能夠傳遍整個(gè)結(jié)構(gòu)而不發(fā)生能量和振幅的衰減；但當(dāng)波處在結(jié)構(gòu)的頻率阻帶范圍時(shí)，波將發(fā)生能量與振幅衰減，而不能傳遍整個(gè)結(jié)構(gòu)。研究還發(fā)現(xiàn)失諧周期結(jié)構(gòu)會(huì)在單元的交界面處出現(xiàn)振動(dòng)及波動(dòng)的局部化現(xiàn)象[910]，導(dǎo)致波動(dòng)幅值沿失諧周期結(jié)構(gòu)以空間指數(shù)的形式衰減并產(chǎn)生能量集中現(xiàn)象，據(jù)此可以對(duì)振動(dòng)波在結(jié)構(gòu)中特定頻率范圍內(nèi)的傳播規(guī)律進(jìn)行深入研究，有利于從理論上指導(dǎo)結(jié)構(gòu)的振動(dòng)控制實(shí)施。

然而，目前對(duì)復(fù)合材料層合結(jié)構(gòu)脫層引起的結(jié)構(gòu)動(dòng)力特性改變問(wèn)題的研究基本都是針對(duì)非周期結(jié)構(gòu)開(kāi)展的，而對(duì)于周期性復(fù)合材料層合結(jié)構(gòu)由脫層引起結(jié)構(gòu)失諧并導(dǎo)致其出現(xiàn)振動(dòng)與波動(dòng)局部化現(xiàn)象的研究尚鮮見(jiàn)報(bào)道。深入掌握其規(guī)律有助于基于振動(dòng)和波傳播特性的周期復(fù)合材料層合結(jié)構(gòu)的脫層缺陷檢測(cè)。

本文對(duì)含局部脫層的周期簡(jiǎn)支層合梁結(jié)構(gòu)的動(dòng)力特性進(jìn)行研究。基于EulerBernoulli梁理論與譜有限元法[11]，分別建立健康與含橫向脫層基本周期單元的動(dòng)剛度矩陣與超單元?jiǎng)觿偠染仃?，并基于傳遞矩陣法分別獲得健康與含脫層基本周期單元傳遞矩陣，進(jìn)而形成含脫層周期梁結(jié)構(gòu)的傳遞矩陣。通過(guò)數(shù)值仿真分別研究基本周期單元長(zhǎng)度、脫層幾何尺寸及位置等參數(shù)改變對(duì)于健康與含脫層基本周期單元傳播常數(shù)的影響，并對(duì)振動(dòng)波在結(jié)構(gòu)中傳播時(shí)出現(xiàn)的幅值衰減和相位變化特征進(jìn)行分析。另通過(guò)單元組裝法建立周期結(jié)構(gòu)總體動(dòng)剛度矩陣，對(duì)基于傳遞矩陣法的計(jì)算結(jié)果進(jìn)行了對(duì)比驗(yàn)證。此外，基于ANSYS仿真計(jì)算，對(duì)健康和含脫層的周期結(jié)構(gòu)在頻率通帶與阻帶下諧響應(yīng)位移進(jìn)行計(jì)算，進(jìn)一步驗(yàn)證本文方法關(guān)于通阻帶特性的傳播常數(shù)分析結(jié)果。

Abstract： Delamination is one of the main types of damage for composite laminated structures. The dynamic stiffness matrix of the undamaged periodic cell of the periodic structure is firstly derived by employing both the EulerBernoulli beam theory and spectral finite element method. Then， considering the delamination is located in one periodic cell of the periodic beam and neglecting the effect of stress singularity at the delamination tip， a delamination model is developed with the plain section assumption for the crosssection of delamination edges. The superelement dynamic stiffness matrix of delaminated periodic cell is formed by utilizing the continuous condition of displacement and rotation at the edge of each subelement as well as the force balance condition at the delamination interface. After that， the transfer matrices of both healthy and delaminated periodic cells are formulated based on transfer matrix method， and then the global transfer matrix as well as the global dynamic stiffness matrix is formulated. The propagation constants obtained from the transfer matrices are employed to investigate the wave propagation properties of the whole periodic structure. In addition， for verification purpose， a set of numerical simulations are carried out for investigating the wave propagation characteristics of a periodically simplysupported beam under both the healthy and disordered conditions. Furthermore， the ANSYS software is also utilized to verify the proposed methodology.

Key words： periodic structure； laminated beam； delamination； spectral finite element method； propagation constant