嵌入式共固化復合材料夾芯板結構的自由振動

翟彥春, 梁 森, 付小靜

(1. 濰坊科技學院 機械工程學院， 山東 壽光 262700；2. 青島理工大學 機械工程學院， 山東 青島 266000)

嵌入式共固化復合材料夾芯板結構的自由振動

翟彥春1，2, 梁 森2, 付小靜2

(1. 濰坊科技學院 機械工程學院， 山東 壽光 262700；2. 青島理工大學 機械工程學院， 山東 青島 266000)

基于一階ZIG-ZAG理論和夾芯層的常復數模型，采用Navier解形式對鋪層方式為(0°/90°/0°)的薄的四邊簡支三層夾芯板的自由振動進行分析。研究了在夾芯板總厚度不變的前提下，固有頻率隨著表層厚度的增加呈現(xiàn)先減小后增大的趨勢，呈上凹拋物線狀；損耗因子值隨著表層厚度的增加呈現(xiàn)先增大后降低的趨勢，分布呈下凸拋物線狀，最大值出現(xiàn)在表層和底層厚度相等時，因此對稱結構夾芯板為最優(yōu)結構設計；在夾芯層厚度和位置一定的前提下，固有頻率的值隨著表層和底層厚度的增加而增大，損耗因子值隨著表層和底層厚度的增加呈現(xiàn)先增加后減小的趨勢，損耗因子存在最大值，為夾芯板結構的最優(yōu)設計提供參考。

復合材料；夾芯板；自由振動；固有頻率；損耗因子

嵌入式共固化復合材料(Embedded Co-cured Composite Damping Structure，ECCDS)不但保持了原有的材料動力學特性又提高了構件的阻尼性能，與傳統(tǒng)材料相比具有良好的整體性能，如：熱穩(wěn)定性好、高比強、高比模、高阻尼及可裁剪性等等，因此在先進工業(yè)領域得到廣泛應用[1-2]。在實際工程中，大量的復雜幾何問題需要用板殼理論來解釋，復合材料層合板和夾芯板是目前研究的熱點方向。

文獻[3]分析了復合材料層合板的拉伸疲勞特性，為疲勞損傷演化的分析提供了依據；Yang等[4]在一階剪切變形理論基礎上，提出一種適用于任何邊界條件的傅里葉-里茲解(Fourier-Ritz Solution)，通過修正傅里葉級數與瑞利里茲法相結合，得到了較為精確的固有頻率和損耗因子解。文獻[5]基于一階ZIG-ZAG理論和Kelvin黏彈性芯材本構特征模型，采用Navier解形式分析了四邊簡支的三層夾芯板的自由振動，探討了阻尼層參數變化對結構固有自由振動特性的影響。文獻[6]基于三階等效單層板理論，采用有限元模型對層合板和夾芯板的彎曲和振動進行了研究。文獻[7]基于Reddy分層理論分別對約束阻尼層合板的層間應力和自由振動進行了研究；文獻[8]基于分層理論建立有限元模型，對層合板的自由振動和阻尼優(yōu)化進行分析。文獻[9]提出面向任意邊界條件和任意非均勻載荷下的角鋪設層合板的彎曲和應力分析的等效單層板理論。文獻[10]對三層復合材料夾芯板的振動和阻尼進行分析，用虛位移原理得到了自由振動的關于截面力和力矩的本構方程，以及相關的邊界條件，使用廣義微分積分法求解上述特征方程，得到頻率和損耗因子值。

目前四邊簡支三層夾芯板的自由振動研究中，重點在夾芯板平衡方程的理論模型推導、求解平衡方程和材料物理參數對自由振動的影響等等，鮮有相關文獻探究薄的四邊簡支夾芯板的最優(yōu)結構布置形式。由于目前復合材料夾芯板在實際工程中的大量應用，因而有必要探究薄夾芯板的最優(yōu)結構設計。

為了探究薄夾芯板的最優(yōu)結構設計，本文基于一階ZIG-ZAG理論和夾芯層常復數模型，采用Navier解形式對鋪層方式為(0°/90°/0°)的薄四邊簡支三層夾芯板的自由振動進行分析，得到的結論可為三層薄夾芯板的最優(yōu)結構設計提供參考。研究了在夾層板總厚度不變的前提下，芯層不同位置對固有頻率和損耗因子的影響，以及在夾芯層厚度和位置一定的前提下，研究了表層和底層不同厚度對固有頻率和損耗因子的影響。

1 夾芯板控制方程的推導及求解

1.1 基本假設

為了推導鋪層為(0°/90°/0°)的薄的四邊簡支三層夾芯板的控制方程及各運動關系，我們做出如下假設：層間界面無滑移；表層和底層板均為正交對稱各向異性層合板，位移符合基爾霍夫假設；垂直于板面方向的變形忽略不計，即εz=0；黏彈層板為正交鋪設層合板，只考慮剪切變形且剪切模量為常復數形式，泊松比為常數。

為了建立夾芯板的位移坐標系，規(guī)定各層板的Z向坐標以各層板的中面為基準。復合材料夾芯板的幾何結構及尺寸，如圖1所示；位移坐標及面內位移關系圖，如圖2所示。

圖1 夾芯板的幾何結構及尺寸Fig.1 Geometry and configuration of the sandwich plate

圖2 夾芯板的面內位移Fig.2 In-plane displacement of sandwich plate

1.2 位移-應變關系

根據所做出的假設和坐標系規(guī)定，表層和底層板的位移表達式為

(1)

根據層間界面無滑移的假設條件，夾芯層中面位移表達式為

(2)

式中：u2，v2分別為夾芯層中面位移；α2，β2分別為夾芯層法線相對于x、y軸的轉角。

進而可以推導出夾芯層的位移表達式為

(3)

第i層的位移-應變關系表達式為

(4)

式中，i=1,2,3。

1.3 截面法向力、剪力、彎矩的表達式

阻尼層的剪力表達式為

(5)

則表層和底層板上的法向力和彎矩的表達式為

(6)

(7)

1.4 復合材料夾芯板控制方程的推導

(8)

1.5 復合材料夾芯板控制方程的求解

四邊簡支復合材料夾芯板的自由振動邊界條件為

x=0,a:

項目施工實際安全保障水平是項目安全狀態(tài)在考核期內的量化體現(xiàn)。安全保障水平越高，項目越安全，發(fā)生安全事故的可能性就越小。以施工安全保障水平指標評價體系為基礎，對各項指標進行定量和定性分析，獲得項目施工實際安全保障水平。

y=0,b:

i=1,3。

采用Navier型解法，滿足四邊簡支邊界條件的位移函數表達式為

u1(x,y,t)=U1cos(nπx/a)sin(mπy/b)eiω*t

v1(x,y,t)=V1sin(nπx/a)cos(mπy/b)eiω*t

w(x,y,t)=Wsin(nπx/a)sin(mπy/b)eiω*t

u3(x,y,t)=U3cos(nπx/a)sin(mπy/b)eiω*t

v3(x，y,t)=V3sin(nπx/a)cos(mπy/b)eiω*t

(9)

M[-(ω*)2]X=(K+jC)X

(10)

式中：M為質量矩陣;(K+jC)為剛度矩陣;X=(U1,U3,V1,V3,W)T。

板的圓頻率和損耗因子按式(11)計算[11]

(11)

2 實例分析

2.1 數值計算

為了檢驗所推導的公式的正確性,分別計算文獻[11]中的(0/90/0)鋪設的四邊簡支夾芯板的振動頻率和損耗因子，材料參數分別見表1和表2，并與文獻值對照。

表1 材料參數

表2 材料參數

表3和表4列出了本文計算的頻率和損耗因子與文獻[11]中的數據的對比，可以看出固有頻率的誤差在7%以內，1～3階損耗因子的誤差在5%以內，4階損耗因子誤差9.8%。表4中的誤差較大是因為式(8)基于經典層合板理論，忽略了剪切變形的影響，而表4中的數據是基于高階剪切變形理論，考慮了剪切變形對頻率和損耗因子的影響，從整體看計算結果說明所推導的公式是正確和可靠的。

表3 固有頻率和損耗因子計算結果

表4 固有頻率和損耗因子計算結果

2.2 表層和底層厚度和一定條件下，改變夾芯層中面位置

從圖3中可以看出，不同模態(tài)下的損耗因子值隨著夾芯層位置下移(表層厚度的增加)呈現(xiàn)先增大后減小的趨勢，呈上凸拋物線狀，最大值出現(xiàn)在表層和底層厚度相等時，即h1=h3，夾芯板為對稱結構，之后逐漸降低，這是因為在對稱結構下，夾芯層會產生最大的剪應力來消耗振動，因此在夾芯板的設計中對稱結構為最佳設計，損耗因子值為最大值。從圖4中可以看出，固有頻率隨著夾芯層位置下移(表層厚度的增加)呈現(xiàn)先減小后增大的趨勢，呈上凹拋物線狀，最小值出現(xiàn)在表層和底層厚度相等時，即h1=h3。綜合說明對稱結構夾芯板在幾和尺寸不變的前提下具有最大損耗因子和最小固有頻率，為最佳設計。

圖3 損耗因子隨夾芯層中面位置變化(h1+h3為定值)Fig.3 Variation of the loss factor with mid-plane of h2(constant h1+h3)

圖4 固有頻率隨夾芯層中面位置變化(h1+h3為定值)Fig.4 Variation of the frequency with mid-plane of h2(constant h1+h3)

2.3 夾芯層厚度和位置一定的前提下，同時改變表層和底層厚度

從圖5和圖6中可以看出，在夾芯層厚度和位置一定的前提下，固有頻率的值隨著表層和底層厚度的增加而增大，1階損耗因子的值隨著表層和底層厚度的增加而增加，達到最大值后基本保持不變，但是2～4階損耗因子的值隨著表層和底層厚度的增加呈現(xiàn)非常明顯的先增加后減小的趨勢，1～4階損耗因子都存在最大值。

圖5 損耗因子隨表層厚度h1的變化(h1=h3)Fig.5 Variation of the loss factor with h1(h1=h3)

圖6 固有頻率隨表層厚度h1的變化(h1=h3)Fig.6 Variation of the frequency with h1(h1=h3)

3 結 論

本文基于一階ZIG-ZAG理論和夾芯層的常復數模型，采用Navier解形式對鋪層方式為(0°/90°/0°)的薄的四邊簡支三層夾芯板的自由振動進行分析。結果表明：

(1)在夾芯板厚度不變的前提下，固有頻率隨著表層厚度的增加呈現(xiàn)先減小后增大的趨勢，呈上凹拋物線狀，最小值出現(xiàn)在表層和底層厚度相等時；不同頻率下的損耗因子的變化情況具有很高的一致性，損耗因子值隨著表層厚度的增加呈現(xiàn)先增大后降低的趨勢，分布呈拋物線狀，最大值出現(xiàn)在表層和底層厚度相等時，因此對稱結構夾芯板的損耗因子最大。

(2)在夾芯層厚度和位置都不變的前提下，固有頻率的值隨著表層和底層厚度的增加而增大，一階損耗因子的值隨著表層和底層厚度的增加而增加，達到最大值后基本保持不變，但是2～4階損耗因子值隨著表層和底層厚度的增加呈現(xiàn)非常明顯的先增加后減小的趨勢，1～4階損耗因子都存在最大值。

[ 1 ] 梁森，梁磊，米鵬. 嵌入式共固化復合材料阻尼結構的新進展[J].應用力學學報, 2010, 27(4):767-771.

LIANG Sen, LIANG Lei,MI Peng. New development of the embedded and co-cured composite damping structres[J]. Chinese Journal of Applied Mechanics, 2010, 27(4):767-771.

[ 2 ] 許廣舉,陳慶樟,李學智,等.復合阻尼材料降低柴油機噪聲的試驗研究[J]. 機械設計與制造, 2015(5):160-163.

XU Guangju,CHEN Qingzhang,LI Xuezhi, et al. Experimental study on reducing noise of diesel engine by composite damping material[J].Mechinery Design & Manufacture, 2015(5):160-163.

[ 3 ] 李祚軍, 張娟, 王佩艷,等. T700/3234碳纖維層合板的拉伸疲勞特性分析[J]. 機械設計與制造, 2013(10):4-6.

LI Zuojun，ZHANG Juan，WANG Peiyan, et al.Tensile and fatigue properities of T700/3234 carbon fiber reinforced composite laminates[J]. Mechinery Design & Manufacture, 2013(10):4-6.

[ 4 ] YANG C, JIN G, YE X, et al. A modified Fourier-Ritz solution for vibration and damping analysis of sandwich plates with viscoelastic and functionally graded materials[J]. International Journal of Mechanical Sciences, 2016, 106:1-18.

[ 5 ] 楊坤, 梅志遠,李華東.黏彈性芯材復合材料夾芯板的自由振動分析[J].海軍工程大學學報, 2012, 24(5):6-11.

YANG Kun,MEI Zhiyuan, LI Huadong. Analysis of free vibration of composite sandwich with viscoelastic core layer[J]. Journal of Naval University of Engineering, 2012, 24(5): 6-11.

[ 6 ] TU T M, LE N T, QUOC T H. Finite element modeling for bending and vibration analysis of laminated and sandwich composite plates based on higher-order theory[J]. Computational Materials Science, 2010, 49(4):390-394.

[ 7 ] 胡明勇, 王安穩(wěn). 纖維增強黏彈性復合材料層合板的自由振動和應力分析[J]. 工程力學, 2010, 27(8):10-14.

HU Mingyong,WANG Anwen.Free vibration and stresses analysis of fiber-reinforced viscoelastic composite laminated plates[J]. Engineering Mechanics, 2010, 27(8):10-14.

[ 9 ] ALIPOUR M M. An analytical approach for bending and stress analysis of cross/angle-ply laminated composite plates under arbitrary non-uniform loads and elastic foundations[J]. Archives of Civil & Mechanical Engineering, 2016, 16(2):193-210.

[10] ARIKOGLU A, OZKOL I. Vibration analysis of composite sandwich plates by the generalized differential quadrature method[J]. Aiaa Journal, 2012, 50(3):620-630.

[11] CUPIAL P, NIZIOL J. Vibration and damping analysis of a three-layered composite plate with a viscoelastic mid-layer[J]. Journal of Sound & Vibration, 1995, 183(1):99-114.

Free vibration of a embedded co-cured composite sandwich plate

ZHAI Yanchun1,2, LIANG Sen2, FU Xiaojing2

(1. College of Mechanical Engineering, Weifang University of Science and Technology, Shouguang 262700,China;2. College of Mechanical Engineering, Qingdao Technological University, Qingdao 266000, China)

Analysis of free vibration of a thin composite sandwich plate with four simply supported edges and (0°/90°/0°) laying-up was conducted using the Navier solution based on the first-order ZIG-ZAG theory and a constant complex model. On the condition that the thickness of the total plate remained unchanged, nature frequency that tended to go from descent to ascent with the increase of the thickness of surface layer presented a concave parabolic distribution; loss factor that tended to go from increase to decrease with the ascent of the thickness of surface layer presented a concave parabolic distribution. So the symmetrical structure sandwich plate is the optimal structure design. On the condition that thickness and location of viscoelastic layer remains as constants, the natural frequencies can be improved with the increase of thickness of the surface layer; loss factor tended to go from increase to decrease with the ascent of the thickness of the surface layer and presented a parabolic distribution. All of this can provide the reference for the optimal design of sandwich plate structure.

embedded co-cured composite; sandwich plate; free vibration; natural frequency; loss factor

國家自然科學基金(51375248)

2016-05-13 修改稿收到日期： 2016-09-21

翟彥春 男，博士生，1984年11月生

梁森 男，博士，教授，1962年6月生

TH212；TH213.3

A

10.13465/j.cnki.jvs.2017.16.036