蔡志東， 鐘銳， 劉濤， 王青山

(1.中南大學 機電工程學院，湖南 長沙 410083； 2.中南大學 高性能綜合制造國家重點實驗室， 湖南 長沙 410083)

復合材料結構由于其高強度、高模量和低密度的特點，越來越被工程領域所重視，尤其在土木工程、車輛工程、海洋工程以及航空工程等領域被廣泛運用[1-5]，而復合材料矩形疊層板結構作為基礎構件，在工程領域的應用十分廣泛。在實際應用過程中，矩形疊層板結構需要適應各種復雜工況，因此其在前期設計過程中的振動特性預測對于設計人員的工況匹配研究至關重要。目前，矩形疊層板結構自由振動特性的相關研究較多，但是大部分都是研究橫向彎曲振動，對于高頻的面內振動特性研究相對較少。然而在工程應用中，矩形板的面內振動也起著極為重要的作用，如直線型超聲電機的研究[6]、超高速飛行器的研制、耦合板結構振動能量傳輸研究[7]等。因此，建立一種適用于一般邊界條件下矩形疊層板結構的面內振動特性分析模型顯得尤為迫切。

Bardell等[8]整理了板結構面內振動的研究文獻，總結了幾類經典邊界條件下板結構的面內振動模態(tài)，并獲得了精確的振動數(shù)據；Farag等[9]研究了對邊固支邊界條件下的板結構面內自由振動特性；Gorman[10-12]采用疊加方法對自由邊界條件下彈性板面內振動固有特性進行了建模分析，并進一步對簡支、固支以及彈性邊界條件下板結構面內自由振動特性研究進行了拓展。史冬巖等[13]采用改進傅里葉級數(shù)法(improved fourier series method,IFSM)對矩形板在任意邊界下的面內自由振動特性進行了研究，結合Rayleigh-Ritz法將矩形板的面內問題轉變成了一個標準特征值求解問題。杜敬濤等[14]采用二維改進傅里葉級數(shù)方法對彈性板結構面內振動耦合位移場函數(shù)進行構建，并開展了相關試驗測量進行驗證分析。目前，對于矩形疊層板的自由振動分析也取得了較多的成果：費雷拉等[15]通過多二次徑向基函數(shù)法和一階剪切變形理論(first-order shear deformation theory，F(xiàn)SDT)求解運動方程，研究了對稱層壓復合板的自由振動。Ngo-Cong等[16]基于FSDT提出了一種徑向基函數(shù)配置方法。Ye等[17]在FSDT的基礎上，開發(fā)了一種修正的傅里葉解決方案，用于分析中等厚度的復合層板的自由振動；Qin等[18]提出了一種Jacobi-Ritz分析方法，用于在任意邊界條件下對復合矩形疊層板的自由振動進行分析。

綜上，目前對于矩形疊層板的面振動特性研究對于橫向振動而言還相對匱乏，且已有的面內振動研究多集中于理論分析層面，相關實驗研究較為少見。為此，本文基于能量法推導了矩形疊層板結構的控制微分方程，結合勒讓德多項式和人工彈簧技術，建立矩形疊層板的數(shù)值分析模型。將求解結果與有限元法的結果進行對比，驗證本文方法的正確性。同時，開展矩形板的面內振動試驗，采用試驗數(shù)據進一步驗證本文方法的有效性。此外，本文針對矩形疊層板的面內振動特性展開參數(shù)化研究，分析比較纖維鋪設角和纖維疊層數(shù)等參數(shù)對面內振動特性的影響，總結了各個參數(shù)對矩形疊層板的面內振動特性的影響規(guī)律。

1 矩形疊層板的模型建立

1.1 模型描述

圖1所示為本文建立的一般邊界條件下矩形疊層板面內振動模型，其長度為a，寬度為b，厚度為h，纖維層數(shù)為Nk，建立坐標系的x軸、y軸和z軸分別與該矩形疊層板的長、寬、高方向相重合。此外，θ表示矩形疊層板第k層的纖維基本方向與x軸之間的夾角。為了將矩形疊層板的面內振動模型統(tǒng)一化，采用人工彈簧技術，也就是圖示中的2組線性彈簧系數(shù)(kt(t=u，v))對邊界條件進行等效模擬，通過設置各邊界約束彈簧的剛度值即可輕松得到不同的邊界條件。

圖1 彈性邊界條件下矩形疊層板面內振動模型Fig.1 In-plane vibration model of laminated rectangular plate with elastic boundary conditions

1.2 模型描述

根據線彈性理論，矩形疊層板在任意點的法向應變和剪切應變可以表示為：

(1)

式中u(x,y)和v(x,y)分別為矩形疊層板在x方向和y方向的面內位移。在胡克定律的基礎上，矩形疊層板的應力可以表示為：

(2)

(3)

(4)

(5)

式中：Nxx和Nyy表示法向力；Nxy表示剪切力。

(6)

式中Aij表示延伸剛度，根據層彈性系數(shù)進行確定。

矩形疊層板的應變能可以表示為：

(7)

在本研究中，采用人工彈簧技術來等效處理一般邊界條件。因此，與邊界彈簧有關的勢能為：

(ku,xLu2+kv,xLv2)|x=a]dy+

(ku,yLu2+kv,yLv2)|y=b]dx

(8)

式中：kt,x0、kt,xL、kt,y0和kt,yL(t=u，v)分別是邊界x=0，x=a，y=0和y=b邊界處的彈簧剛度。矩形疊層板面內振動引起的總動能為：

(9)

1.3 振動求解

在建立矩形疊層板的結構勢能、結構動能以及邊界勢能的基礎上，矩形疊層板的拉格朗日能量函數(shù)可以表示為：

L=T-U+V

(10)

在能量法求解過程中，位移容許函數(shù)對于結構振動特性的求解至關重要。常見的用于構建結構位移場的函數(shù)主要有改進傅里葉級數(shù)與正交多項式(即第一類和第二類切比雪夫多項式、勒讓德多項式、特征正交多項式等)[19]，Qin等[19]在石墨烯圓柱殼的自由振動特性的研究中發(fā)現(xiàn)：正交多項式擁有相似的計算效率和收斂性，且都優(yōu)于改進傅里葉級數(shù)。因此，對于矩形疊層板的面內振動的位移變量u(x,y)和v(x,y)，本文選取勒讓德多項式進行表示：

(11)

(12)

式中：Umn和Vmn是未知的位移未知擴展系數(shù)；Pm(x)和Pn(y)分別是x方向上的階數(shù)m和y方向上的階數(shù)n的勒讓德多項式；M和N分別代表Pm(x)和Pn(y)階數(shù)的最大值。勒讓德多項式的具體表達式為：

(13)

式中：i≥2；χ∈[-1,1]。

特別注意的是，由于勒讓德多項式的求解區(qū)域與結構求解域不一致，因此在代入求解之前需要進行x=ξxχx+ζx和y=ξyχy+ζy的線性坐標變換，其中ξx=ζx=a/2，ξy=ζy=b/2。

在結構位移容許函數(shù)建立的基礎上，將其代入方程(10)，采用里茲法位移未知擴展系數(shù)Umn和Vmn求極值即可獲得矩形疊層板的面內振動特性求解方程：

(K-Mω2)A=0

(14)

式中：K、M和A分別表示剛度矩陣、質量矩陣和未知系數(shù)的向量。通過求解方程(14)，可以獲得矩形疊層板的自然頻率和相關特征向量，然后代入方程(10)即可獲得相對應的面內振動模態(tài)振型。

2 數(shù)值與試驗結果分析

本文已經對矩形疊層板的面內自由振動特性求解方程進行了描述，在獲得理論推導的基礎上，本文首先進行模型收斂性研究，確定邊界剛度值與多項式截斷值，隨后開展數(shù)值仿真研究，并將計算結果與ABAQUS算得的結果和矩形板的面內振動試驗結果進行對比。在此基礎上展開參數(shù)化研究，討論各參數(shù)對矩形疊層板面內振動特性的影響規(guī)律。

2.1 模型收斂性研究

正如前面的建模過程可知，對于矩形疊層板的邊界條件模擬是通過在邊界位置處引入了線性彈簧來進行等效處理。但是對于固定邊界的等效需要設置一個合理大小的彈簧剛度值來進行近似處理，所以有必要開展邊界剛度的收斂性研究，相關的計算結果如圖2所示。在圖2中，給出了2種纖維敷設方案的矩形疊層板，方案分別為[0°/30°/0°]和[45°/-45°/45°/-45°]。結構參數(shù)為1 m×1 m×0.05 m，材料參數(shù)為：E1=150 GPa；E2=10 GPa；G12=5 GPa；G13=5 GPa；G23=6 GPa。ku和kv分別代表矩形疊層板沿x軸與y軸方向的邊界剛度值，坐標值分別取以10為底ku和kv的對數(shù)值，采用等高線云圖表現(xiàn)頻率的收斂性，剛度值分別從104變化到1014。

圖2 邊界剛度值收斂性分析Fig.2 Convergence analysis of boundary stiffness value

從圖2中可以看出，當邊界彈簧ku以及kv的剛度值在107～1010N/m內變化時，頻率參數(shù)隨之變化(圖中綠色區(qū)域)，此時可將其定義為彈性邊界(E)；而當結構的邊界剛度值在1010～1014N/m變化時，頻率參數(shù)已基本收斂(圖中紅色區(qū)域)，對應的邊界可視為固定支撐邊界(C)；而剛度值為0時(圖中藍色區(qū)域)，可將邊界定義為自由邊界(F)。

接下來，開展勒讓德多項式收斂性研究，將截斷值M和N從2取值到12得到的收斂性結果如圖3所示，結構邊界條件分別為C-C-C-C、F-F-F-F和C-F-C-F邊界，結構參數(shù)和材料參數(shù)與圖2保持一致，纖維敷設角度為[45°/-45°/45°/-45°]。從圖中可以看出，隨著截斷值M和N的增加，計算結果迅速收斂，當截斷值達到12時，與截斷值為10時的誤差均小于0.5‰。雖然從理論上講，當截斷值無限大時，計算結果可以完全等于精確分析結果。但是，考慮到實際的硬件計算資源和計算效率，通常需要選擇一個合理的截斷值。因此，在后續(xù)算例分析中，勒讓德多項式的截斷值選定為M=N=12。

圖3 勒讓德多項式截斷值收斂性分析Fig.3 Convergence analysis of truncation value of Legendre polynomial

2.2 模型正確性研究

表2給出了矩形疊層板在本文求解分析模型和有限元分析模型下所獲得的前6階自然固有頻率，結構參數(shù)和材料參數(shù)與圖2保持一致。對于結構振動特性而言，模態(tài)振型可以直觀反映出結構的固有特性，因此圖4還給出了本文方法與有限元法所得到的模態(tài)振型對比圖。通過表1數(shù)據和圖4可知，采用本文方法計算所得的數(shù)值與ABAQUS仿真得到的結果吻合較好，數(shù)值誤差較小。因此，可以驗證本文所提出的方法對矩形疊層板面內自由振動計算的有效性。

為了進一步驗證該模型的正確性，本文開展了實驗對比分析，搭建了如圖5中所示的振動試驗系統(tǒng)。試驗件采用鋁板材激光加工而成，尺寸為a×b×h=1 250 mm×1 127 mm×10 mm，彈性模量為E=70 GPa，泊松比μ=0.33。在鋁板頂邊兩端加工有2個吊裝孔，通過彈性繩將矩形鋁板吊裝在臺架上，從而實現(xiàn)了試驗件的自由約束邊界條件設置。測試方法選用單點面內激勵多點拾振，激勵點選擇為1號測點的x方向，將矩形鋁板的表面被均勻劃分為10×10個網格，共計121個測點，在測點處布置三向加速度計，按照編號依次挪動傳感器位置測量各測點的響應數(shù)據。

對于每次測量采用力錘進行敲擊，拾取被測點x、y方向上的振動響應信號，并保存在數(shù)據與采集分析系統(tǒng)中，力錘敲擊的激勵力信號數(shù)據也被同步保存在系統(tǒng)中。每個測點測量5次數(shù)據，系統(tǒng)自動求取平均值進行數(shù)據處理，以減少試驗誤差。測試完成后，通過東華DHDAS軟件輔助工具處理試驗中產生的激勵與響應數(shù)據，最終可以得到該矩形鋁板的面內模態(tài)試驗結果。

表1 固支邊界下矩形疊層板固有頻率參數(shù)對比

采用本文方法預測與試驗測試的得到前8階面內模態(tài)頻率對比如表2所示。通過表2可知，采用本文方法得到預測值也可以很好地與試驗測量值吻合，前8階面內模態(tài)頻率的誤差普遍在1%以內，最大誤差不超過2%。綜合分析試驗結果，可以進一步論證本文所建立的矩形疊層板面內自由振動模型的準確性。

圖4 本文方法與有限元法得到的前3階陣型圖對比Fig.4 Comparison of the first three-order array diagram obtained by present method and the finite element method

圖5 矩形板面內振動模態(tài)測試系統(tǒng)Fig.5 In-plane vibration mode test system of rectangular plate

表2 前8階面內模態(tài)頻率預測與測量結果對比

2.3 參數(shù)化研究

通過現(xiàn)有文獻調研可知，對于矩形疊層板振動特性研究目前大多數(shù)都是研究低頻的彎曲振動，對于高頻振動的面內自由振動的研究還十分鮮有。因此為了豐富現(xiàn)有研究成果，本小節(jié)在上述收斂性和有效性研究的基礎上，進一步開展參數(shù)化研究，主要包括纖維鋪設角和纖維疊層數(shù)等因素。

為研究材料纖維鋪設角的變化對矩形疊層板振動特性的影響，選取了2組纖維鋪設方案：[0°/α°/0°]和[0°/α°/0°/α°]，其中敷設角度α°從0°到90°變化。對于每種敷設方案，選取了3種不同的邊界條件CCCC、FFFF和CFCF，幾何和材料屬性與表1相同，選取前3階模態(tài)頻率進行參數(shù)化結果展示，如圖6所示。從圖中可知，各邊界下2種不同鋪設方案的鋪設角變化對結構面內振動頻率的影響具有相似的規(guī)律，且對于CCCC及CFCF邊界，頻率隨著鋪設角的增大，總體上呈現(xiàn)出增大的趨勢，而對于FFFF邊界，其變化大體是先增大后減少的趨勢，這主要是由于不同邊界條件對于結構的剛度影響存在差別引起的。

為研究纖維層數(shù)對矩形疊層板面內自由振動的影響，對不同纖維層數(shù)下的面內自由振動頻率進行分析計算。如圖7所示，矩形疊層板采用奇數(shù)層纖維鋪設角度為0°，偶數(shù)層層數(shù)的纖維鋪設角為90°的纖維鋪設方案，隨著纖維層數(shù)Nk取值的不同，纖維敷設層數(shù)也隨著發(fā)生相應的變化，結構邊界條件分別為CCCC、FFFF和CFCF，選取前3階模態(tài)頻率進行參數(shù)化結果展示。從圖7中可以看出，隨著纖維層數(shù)的增大，偶數(shù)層的模態(tài)頻率保持一致，這是因為剛度系數(shù)不隨偶數(shù)層數(shù)目Nk的變化而改變。而奇數(shù)層的改變，會導致剛度系數(shù)隨之發(fā)生變化，因此其對應的頻率參數(shù)也會發(fā)生相應的波動。

圖6 不同纖維敷設方案前3階面內頻率隨鋪設角變化規(guī)律Fig.6 Variation of frequency parameters with laying angle under different lamination schemes

3 結論

1)對比分析3種方法得到的數(shù)據結果，證明了本文計算方法的正確性和有效性。

2)由于不同邊界條件對于結構的剛度影響，矩形疊層板的纖維鋪設角變化對結構面內振動頻率有較大影響。

3)矩形疊層板的頻率參數(shù)會隨著鋪設層數(shù)的變化而發(fā)生波動。

4)通過本文的研究可以對未來進一步的更復雜的復合材料結構在高頻區(qū)域的振動特性研究提供一定的理論基礎。