黏彈性層狀周期板動力計算的近似理論與解答

徐司慧 王炳龍 周順華 楊新文 李堯臣

*(同濟大學道路與交通工程教育部重點實驗室，上海201804)

?(同濟大學航空航天與力學學院，上海200092)

黏彈性層狀周期板動力計算的近似理論與解答

徐司慧*,1)王炳龍*周順華*楊新文*李堯臣?

*(同濟大學道路與交通工程教育部重點實驗室，上海201804)

?(同濟大學航空航天與力學學院，上海200092)

由于周期性隔振結(jié)構(gòu)動力計算中較少考慮軌道交通載荷及材料黏彈性,因此，本文以黏彈性層狀周期板為研究對象，提出了垂向移動簡諧載荷下，可以考慮材料黏彈性及板內(nèi)橫向剪切變形的黏彈性層狀周期板動力計算近似理論并給出解析解答.設(shè)板中性面的橫向剪切變形為橫截面的整體剪切變形，利用Reissner-Mindlin假設(shè)及提出的剪切變形補充計算條件，得到了中性面法線轉(zhuǎn)角與中性面剪應(yīng)力的關(guān)系.基于平衡方程和應(yīng)力連續(xù)條件，建立了黏彈性層狀周期板振動控制方程，推導了對邊簡支對邊自由條件下，板垂向位移的簡化Fourier級數(shù)形式解.與經(jīng)典層合板模型和有限元計算結(jié)果進行了比較，驗證了本文解答的有效性.結(jié)果表明：(1)黏彈性層狀周期板可以顯著降低單一材料板在自振頻率處的振動響應(yīng)，但會引起局部低頻頻段的振動放大；(2)板的垂向位移隨著載荷速度的增大而增大，當載荷速度超過300km/h后，其對板振動響應(yīng)的影響減弱；(3)黏彈性層剪切模量存在最佳設(shè)計值，可使結(jié)構(gòu)的隔振性能最佳；(4)黏彈性層的阻尼特性在低頻范圍內(nèi)對結(jié)構(gòu)振動影響較?。?5)可在滿足工程實際的情況下適當增加板長，以提高結(jié)構(gòu)的隔振性能.

周期性結(jié)構(gòu)，軌道交通載荷，Reissner-Mindlin假設(shè)，解析解，振動響應(yīng)

引言

周期性結(jié)構(gòu)是指由基本周期單元按照相同的排列方式在單方向或多方向上堆積而成的結(jié)構(gòu)[1-2].固體物理學研究表明，周期性結(jié)構(gòu)具有獨特的頻率衰減域特性，即處于某些頻段(衰減域)的波不能順利透過該結(jié)構(gòu)，導致振動減弱[3-5].受此啟發(fā)，石志飛等[6-8]于2007年將周期性的結(jié)構(gòu)引入土木工程領(lǐng)域，提出了周期性隔震基礎(chǔ)和周期性排樁等，隔震/振效果顯著，為周期性結(jié)構(gòu)在實際工程中的應(yīng)用打下了基礎(chǔ).

目前，研究較廣泛的周期性隔震/振基礎(chǔ)主要為在單方向上呈周期性排列的層狀周期板，這也是周期性結(jié)構(gòu)的基本形式.研究內(nèi)容集中于地震載荷下的建筑物隔震技術(shù)，采用的理論和方法主要有傳遞矩陣法、有限元法、模型試驗等.石志飛等[9]利用傳遞矩陣法研究了地震載荷作用下層狀周期板基礎(chǔ)的動力衰減域特性，得到了衰減域的主要控制因素；黃建坤[10]采用平面波展開法研究了周期性排樁和波屏障在環(huán)境減振中的應(yīng)用，并設(shè)計出了具有隔振效果的波屏障；陳志寶等[11-12]通過振動臺試驗證明了具有黏彈性的層狀周期結(jié)構(gòu)對外部激勵具有較好的隔離或減弱效果.上述研究雖然取得了一系列有益成果，但是載荷形式較單一，多為地震載荷，較少涉及工程中廣泛存在的軌道交通載荷，難以體現(xiàn)載荷的作用范圍以及移動速度.

軌道交通載荷是軌道交通結(jié)構(gòu)振動的主要振源，研究軌道交通載荷下的結(jié)構(gòu)減振技術(shù)對軌道交通領(lǐng)域的環(huán)境振動控制意義重大.本文將周期性結(jié)構(gòu)應(yīng)用于軌道交通減振領(lǐng)域，采用黏彈性層狀周期板作為軌道交通領(lǐng)域板結(jié)構(gòu)的基本形式.為了考慮軌道交通載荷及材料的黏彈性，采用經(jīng)典軌道動力學中輪軌垂向作用力的基本形式，即垂向移動簡諧載荷[13-14]，并將具有黏彈性的層狀周期板視為復合材料層合板結(jié)構(gòu)，利用層合板理論進行動力計算.從既有文獻的調(diào)研結(jié)果來看，求解層合板問題常采用基于Kirchho ff假設(shè)的經(jīng)典板理論[15]，即將每層板視為不考慮橫向剪應(yīng)力的薄板，然后聯(lián)立求解.然而，由于忽視了板內(nèi)橫向剪應(yīng)力，經(jīng)典層合板理論無法考慮板的面內(nèi)位移在板厚方向上的zig-zag效應(yīng)以及層間剪應(yīng)力連續(xù)條件[16-18]，因此無法準確計算板面的撓度分布，難以預測層狀周期板的減振效果.

基于經(jīng)典層合板理論，學者從優(yōu)化剪切函數(shù)、細化應(yīng)力場分布等角度提出了多個復合材料層合板理論[19-21].Chen等[22-23]采用三維彈性理論得到了對邊簡支層合板的自由振動半解析解.Yang等[24]利用高階變形理論對夾芯板的受迫振動進行了分析.Faruk Firat Calim[25]在拉氏域內(nèi)分析了非均勻復合材料梁的自由振動和強迫振動，研究了非均勻參數(shù)對梁動力特性的影響.Ganapathi等[26]引入材料特性剪切修正因子，結(jié)合一階剪切變形理論建立了層合板的振動方程.李堯臣等[27-28]利用簡化理論對壓電功能梯度圓板和矩形板進行了分析并得到了近似解析解.Thai等[29]基于一階剪切變形板理論來計算夾芯板的力學特性并得到了很好的結(jié)果.Pandey和Pradyumna[30]利用Layerwise理論對功能梯度夾芯板的自由振動問題進行分析.一般認為，復合材料層合板理論可以得到簡單邊界條件下較為精確的結(jié)構(gòu)整體響應(yīng)，比如撓度、基頻振動頻率等.然而，由于獨立場變量的個數(shù)隨層合板鋪層數(shù)的增加而增加，導致該類方法計算量較大，不適用于大規(guī)模工程問題的計算.因此，針對黏彈性層狀周期板，探索表達形式和算法相對簡潔，同時精度又能滿足要求的近似理論很有必要.

本文將周期性結(jié)構(gòu)應(yīng)用于軌道交通減振領(lǐng)域，采用黏彈性層狀周期板作為軌道交通領(lǐng)域板結(jié)構(gòu)的基本形式，并考慮軌道交通結(jié)構(gòu)振動計算的特點.提出了垂向移動簡諧載荷作用下，對邊自由對邊簡支板動力計算的近似理論與解答.首先，引入了板理論的Reissner-Mindlin假設(shè)，根據(jù)等效截面法確定層狀周期板中性軸位置，并結(jié)合應(yīng)力邊界條件計算出中性面上橫向剪應(yīng)力；然后，對Reissner-Mindlin假設(shè)補充計算條件，以中性面上的橫向剪切變形作為板橫截面上的整體剪切變形，建立了中性面法線轉(zhuǎn)角與中性面上剪應(yīng)力的關(guān)系；最后，結(jié)合平衡方程和應(yīng)力連續(xù)條件，推導了黏彈性層狀周期板振動控制方程，提出了考慮橫向剪切變形的黏彈性層狀周期板近似理論及移動簡諧載荷下黏彈性層狀周期板振動位移簡化Fourier級數(shù)形式解.分別與傳統(tǒng)方法和有限元計算結(jié)果進行了比較，該解在形式和計算上都較為簡潔，且精度滿足要求，可為黏彈性層狀周期板在軌道交通減振領(lǐng)域的設(shè)計及推廣提供一定理論基礎(chǔ).

1 黏彈性層狀周期板近似理論及求解

1.1 基本方程

黏彈性層狀周期板沿板厚方向設(shè)置周期子結(jié)構(gòu)層，如圖1所示.圖1中，層A和層B為一個周期單元：層A為彈性層，層高為h，長度l，寬度b；層B為黏彈性層，尺寸同層A.黏彈性層B在提供一定黏彈性的同時也可作為結(jié)構(gòu)承重層，可選擇的黏彈性層B的材料有瀝青混凝土[31]、環(huán)氧瀝青混凝土[32]等.

對圖1所示的黏彈性周期板做垂向簡諧載荷作用下的動力分析，為考慮板內(nèi)橫向剪應(yīng)力及層間剪切變形，引入Reissner-Mindlin假設(shè)，即：

(1)變形前的中性面法線在變形后仍保持為直線，但不一定與中性面垂直；

(2)與面內(nèi)應(yīng)力σx,σy,τxy相比，另外3個應(yīng)力分量σz,τzx,τzy很小，在計算應(yīng)變時可以忽略不計.

考慮材料為各向同性，黏彈性周期板的基本方程如下：

圖1 黏彈性層狀周期板示意圖Fig.1 Viscoelastic layered periodic plate

應(yīng)力平衡方程(不計體力)

幾何方程

黏彈性層物理方程

彈性層物理方程

式中，σx,σy,σz,τzy,τzx,τxy為應(yīng)力；εx,εy,εz,γyz,γzx,γxy為應(yīng)變；u,v,w分別為x,y,z方向的位移；Y(1)(t)和Y(2)(t)為松弛模量，分別描述應(yīng)力的偏量部分和球量部分的松弛特性[33]；*表示斯蒂爾吉斯卷積；C11，C12，C13，C33，C44，C66為彈性系數(shù).

由假設(shè)(2)可知σz=0，代入式(3a)和式(3b)并進行Laplace變換，可得

對彈性層A(各向同性材料)，有

對黏彈性層B，有

1.2 移動簡諧載荷下板振動位移解析求解

計算黏彈性層狀周期板中性面的位置，設(shè)中性面到中面的距離為z0(如圖2所示)，則中性面的位置可由下式[34]確定

式中上標(1)～(4)表示板每一層編號.

圖2 中性面位置示意圖Fig.2 Location of the neutral plane

令中性面上z=0，由假設(shè)(1)可得板內(nèi)位移為

式中，βx和βy分別為中性面法線繞x和y軸的轉(zhuǎn)角，正負號按右螺旋法則確定；u0和v0為中面上的面內(nèi)位移.

將式(1)和式(7)進行Laplace變換得

假設(shè)初始位移和初始速率為0，將式(4)代入式(8a)和式(8b)，并利用式(9)和幾何方程式(2)，再從0到z積分，考慮到初值定理，可得

由圖2可知，板的上表面為ztop=2h?z0，下表面為zbot=?(2h+z0).將ztop及zbot表達式代入式(10)和式(11)可得

式中，標號(A)和(B)分別代表板一個周期單元內(nèi)，彈性層A和黏彈性層B的材料參數(shù).

令板上下表面剪應(yīng)力為0，將式(13)和式(14)代入式(12),并從?z0?2h到2h?z0對z分段積分，得到

式中

當板上的載荷為沿板中心線勻速移動的單位簡諧載荷時

式中，x0和y0為載荷初始位置的坐標，v為載荷的移動速度，wf為載荷的角頻率，Kd為板下支承層剛度.

式(15)包含未知量w,βx,βy，想要求解垂向位移w，根據(jù)本構(gòu)關(guān)系(3a)和(3b)和變形幾何條件(2)，對Reissner-Mindlin假設(shè)引入補充條件：

(1)以中性面上的橫向剪切變形作為板橫截面上的整體剪切變形；

(2)變形前的中性面法線在變形后與變形后的中性面法線的夾角等于該處中性面上的剪應(yīng)變.于是有

將式(13)和式(14)代入式(19)并進行Laplace變換，得

到此，聯(lián)立式 (15)，式 (17)，式(18)，式 (20)和式(21)，并根據(jù)邊界條件，即可求解單位垂向簡諧載荷下層狀周期板的振動位移.

設(shè)板在x=0和x=l的邊界上簡支，即

將式(22a)～式 (22d)代入式(15)，式(17)，式(18)，式(20)和式(21)可得求解級數(shù)系數(shù)的方程

式中，ri為級數(shù)系數(shù)方程對應(yīng)齊次方程的特征根；為與ri有關(guān)的系數(shù)；為待定常數(shù)，由y=0,y=b兩邊的邊界條件以及板在載荷作用點兩側(cè)的位移、轉(zhuǎn)角、彎矩、扭矩連續(xù)和剪力的平衡條件確定.兩邊自由時，滿足以下邊界條件

式中，My，Qy，Mxy分別為zoy平面內(nèi)的彎矩、剪力、扭矩，按傳統(tǒng)方法定義.

為了得到時域中板的動力響應(yīng)，需對Laplace變換域中相應(yīng)的物理量進行逆變換.由于象函數(shù)復雜，且精度要求較高，本文采用Durbin改進后的Laplace逆變換法[35]

取cT=5,T=20,k=2000.

2 模型驗證

采用本文模型計算黏彈性層狀周期板的振動位移并與既有計算方法進行比較.計算周期板的幾何尺寸如圖3所示.由圖3，板長度l=4.962m，寬度b=2.55m，厚度4h=0.2m.彈性層選用混凝土，黏彈性層選用黏彈性橡膠材料，相關(guān)計算參數(shù)如表1所示.

圖3 板幾何尺寸Fig.3 Geometric size of the plate

表1 驗證模型計算參數(shù)Table 1 Parameters of model

分別采用經(jīng)典層合板模型 (基于 Kirchho ff假設(shè))、本文模型及有限元法計算某時刻板中心線上的垂向位移w，如圖4所示.

圖4 本文模型解答與既有計算方法比較(t=0.2s)Fig.4 Comparison of displacement curves by presented and other solutions(t=0.2s)

經(jīng)典層合板模型是指將式(19)中略去剪切變形項，即

式(31)表示變形前的中性面法線在變形后仍與中面垂直，即經(jīng)典層合板模型采用的Kirchho ff假設(shè).將式(31)代入式(15)可得只含板振動位移w的控制方程，即

將式(22a)代入式(32)，并結(jié)合邊界條件式(33)即可求得板垂向位移w

有限元解采用有限元軟件ABAQUS通過三維建模得到，選用六面體八節(jié)點縮減積分單元(C3D8R)劃分實體單元網(wǎng)格，并在載荷區(qū)域進行網(wǎng)格加密，單元數(shù)為80500.

由圖 4可知，本文模型解與有限元解吻合較好，可整體反應(yīng)出板在載荷作用下的瞬態(tài)變形，經(jīng)典模型解與有限元解有約12%的誤差，明顯大于本文模型解與有限元解的差值.究其原因，數(shù)值解是基于有限元軟件計算得到的三維數(shù)值結(jié)果，單元數(shù)量足夠多，可以認為該解答的精度最高，與實際情況最接近.經(jīng)典模型忽略了板內(nèi)的橫向剪切變形，即無法考慮剪變形和擠壓變形引起的位移[37]，故計算結(jié)果小于有限元解.本文模型在經(jīng)典模型的基礎(chǔ)上考慮了板的橫向剪切變形，即增加了剪變形引起的位移項，因此比經(jīng)典模型解更接近于有限元解.但是本文模型仍然對板的內(nèi)部變形做了約束，人為定義了板內(nèi)的橫向剪應(yīng)力分布，導致板相較真實情況剛度增大，所以本文模型解仍小于有限元解，但誤差在合理范圍內(nèi).

表2為板上4個測點的振動位移在不同計算方法下的數(shù)值.由表2可知，測點1與載荷作用點距離大，振動位移幅值小，且解析解與有限單元解差別較小.測點4比測點1更接近載荷作用點，振動位移幅值大，且振動位移的解析解與有限單元解的差值增大.在本文計算條件下，本文模型與有限元解的最大差值在合理范圍內(nèi)，可以認為計算模型是有效的.

表2 點1～4垂向位移的解析解與有限單元解比較(y=1.275m)Table 2 Comparison of displacement by analytical and FEM solutions at some points(y=1.275m)

進一步對比不同時刻下板的變形規(guī)律，如圖5所示.圖5為t=0.1s和0.2s時刻下板的整體位移響應(yīng)，由圖可看出，隨著計算時間增加，載荷沿著板中心線移動，板位移的峰值也相繼發(fā)生變化.圖5(a)和圖5(b)中板位移峰值分別為x=1.1m和x=2.2m，計算結(jié)果吻合類似文獻[38]和工程實踐的認識，可以認為本文所得模型是合理的.

圖5 不同時刻下板的振動響應(yīng)Fig.5 Vibration response of plate at di ff erent time

3 計算結(jié)果分析

3.1 導納特性分析

本文從頻域角度出發(fā)，以頻率響應(yīng)函數(shù)(導納)為評價指標，研究板的垂向振動特性.圖6為板的垂向位移導納，計算參數(shù)見表1.由圖6可知，黏彈性層狀周期板結(jié)構(gòu)及材料組成多樣，振動模態(tài)豐富，垂向位移導納呈現(xiàn)多個峰值.位移導納在達到第5階共振頻率前，隨頻率的增大而增大；超過第5階共振頻率后，隨著頻率的增大而減小.當激勵頻率在50Hz附近時，板位移導納達到最大，周期板在共振頻率下的垂向位移導納見表3.

圖6 板垂向位移導納Fig.6 Vertical displacement admittance of plate

表3 板位移導納峰值Table 3 Vertical displacement admittance peak of plate

圖7 兩種板位移導納的對比Fig.7 Comparison of vertical displacement admittance by di ff erent plates

為衡量黏彈性層狀周期板相較于單一材料板的減振能力，圖7對比了兩種板的垂向位移導納值，其中，單一材料板的材料參數(shù)同黏彈性層狀周期板的彈性層，具體數(shù)值見表1.由圖7可知，單一材料板自振頻率成分單一，位移導納極值為4.6×10?7m/N，對應(yīng)頻率為51Hz；黏彈性層狀周期板導納極大值為6.6×10?8m/N，對應(yīng)頻率為48Hz.由此可見，采用黏彈性層狀周期板，可以顯著降低單一材料板在自振頻率處的振動響應(yīng)，最大減少量達到3.9×10?7m/N，占整體位移響應(yīng)的85%，且兩種板位移導納極值對應(yīng)的自振頻率接近.但是，在30～40Hz范圍內(nèi)，黏彈性層狀周期板的振動響應(yīng)大于單一材料板，最大放大值可達到1.25×10?8m/N，占位移響應(yīng)的28%；60Hz以上的頻段，兩種板振動響應(yīng)相差不大.

3.2 載荷移動速度的影響

取載荷速度分別為40km/h,100km/h,200km/h,300km/h,400km/h,500km/h,600km/h，計算黏彈性層狀周期板垂向振動位移的極大值，計算結(jié)果如圖8所示.由圖8可知，在本文計算條件下，隨著載荷移動速度的增加，板垂向振動位移增大.載荷速度小于300km/h時，板振動位移增大幅度顯著，尤其是100km/h～200km/h區(qū)間，位移增大值達 2.5×10?8m.當載荷速度超過300km/h，板振動位移增大幅度減緩：300km/h～500km/h區(qū)間，位移增大0.84×10?8m，500km/h～600km/h區(qū)間，位移增大 0.08×10?8m.

圖8 載荷移動速度的影響Fig.8 The in fl uence of load velocity

3.3 材料參數(shù)及幾何尺寸的影響

探討?zhàn)椥詫硬牧蠀?shù)及板幾何尺寸對隔振性能的影響，采用力傳遞率評價板的隔振性能[39]，定義為

式中，F(xiàn)stran是板傳給板下基礎(chǔ)的力的幅值和，F(xiàn)stran=|Fz|=Kd|w|；F是外激勵的幅值.

圖9為穩(wěn)態(tài)振動下，黏彈性層剪切模量GB=10?6GA～10?2GA時，板力傳遞率的變化曲線.由圖8可看出，黏彈性層剪切模量較小時，結(jié)構(gòu)力傳遞率較大，達到0.88；增加黏彈性層的剪切模量，結(jié)構(gòu)力傳遞率先快速減小，從0.88降到了0.22，減小量達75%；之后，當黏彈性層與彈性層剪切模量比值超過10?4后，結(jié)構(gòu)力傳遞率趨近穩(wěn)定，呈緩幅增大，從0.22增加到0.26.由此可見，黏彈性層剪切模量存在最佳設(shè)計值，可使結(jié)構(gòu)的力傳遞率達到最小，這給黏彈性層狀周期板的材料參數(shù)設(shè)計提供有益思路.

圖9 力傳遞率隨著剪切模量的變化曲線Fig.9 Force transfer rate varies with shear modulus

圖10為穩(wěn)態(tài)振動下，黏彈性層損耗因子的變化對板力傳遞率的影響.由圖10可知，黏彈性層的阻尼特性對結(jié)構(gòu)振動響應(yīng)具有一定抑制作用.當結(jié)構(gòu)沒有黏性，即損耗因子取0時，力傳遞率為0.598.之后損耗因子增至0.5，力傳遞率緩幅減小，但衰減幅度低，力傳遞率僅降低0.02.究其原因，激勵頻率較低，為20Hz，黏彈性層的阻尼特性在低頻范圍內(nèi)對結(jié)構(gòu)振動響應(yīng)的影響較小.

圖10 力傳遞率隨損耗因子的變化曲線Fig.10 Force transfer rate varies with loss factor

圖11為穩(wěn)態(tài)振動下，板的長度變化對隔振性能的影響曲線，評價指標取單位長度力傳遞率，即力傳遞率T與板長l的商.由圖11可知，隨著板的長度增加，結(jié)構(gòu)參振質(zhì)量增大，力傳遞率減小，隔振性能增加.當板的長度由5m增加至11m時，單位長度力傳遞率衰減顯著，從0.197降到了0.077，衰減度達76%；之后繼續(xù)增加板長，結(jié)構(gòu)單位長度力傳遞率平緩降低，從0.077減小至0.039.對于軌道交通領(lǐng)域的板結(jié)構(gòu)，板的長度受限于施工條件等因素，因此可考慮適當增加板長，以提高結(jié)構(gòu)的隔振性能.

圖11 板長對隔振性能的影響Fig.11 Vibration isolation performance varies with length

4 結(jié)論

針對單位移動簡諧載荷下，對邊簡支對邊自由黏彈性層狀周期板動力計算問題，本文提出了相應(yīng)的近似理論與解析解答.模型計算形式簡潔，計算效率高，可以考慮板內(nèi)橫向剪應(yīng)力的影響.通過與經(jīng)典層合板模型和三維有限元計算結(jié)果的比較，驗證了解答的有效性.本文模型可為黏彈性層狀周期板在軌道交通減振領(lǐng)域的設(shè)計及推廣提供一定理論基礎(chǔ).

利用本文計算方法研究了某黏彈性層狀周期板在頻域內(nèi)的動力響應(yīng)，得到如下結(jié)論：

(1)采用黏彈性層狀周期板可以顯著降低單一材料板在自振頻率處的振動響應(yīng)，但會引起局部低頻頻段的振動放大；

(2)本文計算條件下，板的垂向位移隨著載荷速度的增加而增大，超過300km/h后，載荷速度對板振動響應(yīng)的影響減弱；

(3)黏彈性層剪切模量存在最佳設(shè)計值，可使結(jié)構(gòu)的隔振性能最佳；

(4)黏彈性層的阻尼特性在低頻范圍內(nèi)對結(jié)構(gòu)振動的影響較小;

(5)可在滿足工程實際的情況下適當增加板長，以提高結(jié)構(gòu)的隔振性能.

1 Kushwaha MS,Halevi P,Mart′?nez G,et al.Theory of acoustic band structure of periodic elastic composites.Phys Rev B Condens Matter,1994,49(4):2313-2322

2 Brillouin L.Wave Propagation in Periodic Structures:Electric Filters and Crystal Lattics.Dover Publications,Inc,2003

3 溫激鴻，王剛，郁殿龍等.聲子晶體振動帶隙及減振特性研究.中國科學,2007,37(9):1126-1139(Wen Jihong,Wang Gang,Yu Dianlong,et al.Vibration attenuation and band gap characteristics of phononic crystals.Scientia Sinica Technologica,2007,37(9):1126-1139(in Chinese))

4 黃毓.彈性波帶隙材料/結(jié)構(gòu)優(yōu)化設(shè)計.[博士論文].大連：大連理工大學,2014(Huang Yu.Optimization design of elastic wave bandgap material/Structures.[PhD Thesis].Dalian:Dalian University of Technology,2014(in Chinese))

5 朱曉輝，李隆球，張廣玉等.基于層狀周期性結(jié)構(gòu)的聲波調(diào)控技術(shù)研究.機械工程學報,2017,53(6):10-15(Zhu Xiaohui,Li Longqiu,Zhang Guangyu,et al.Investigation of acoustic manipulation by layered periodic composites.Journal of Mechanical Engineering,2017,53(6):10-15(in Chinese))

6 Shi Z,Huang J.Feasibility of reducing three-dimensional wave energy by introducing periodic foundations.Soil Dynamics&Earthquake Engineering,2013,50(1):204-212

7 Yan Y,Laskar A,Cheng Z,et al.Seismic isolation of two dimensional periodic foundations.Journal of Applied Physics,2014,116(4):1-12

8 Huang J,Shi Z,Huang W,et al.A periodic foundation with rotational oscillators for extremely low-frequency seismic isolation:analysis and experimental veri fi cation.Smart Material Structures,2017,26

9 程志寶，石志飛，向宏軍.層狀周期結(jié)構(gòu)動力衰減域特性研究.振動與沖擊,2013,32(9):178-182(Cheng Zhibao,Shi Zhifei,Xang Hongjun.Vibration attenuation zones of a laminated periodic structure.Journal of Vibration and Shock,2013,32(9):178-182(in Chinese))

10 黃建坤.周期性排樁和波屏障在土木工程減振中的應(yīng)用研究.[博士論文].北京:北京交通大學,2014(Huang Jiankun.Application of periodic pile barriers and wave barriers to reduce vibrations in civil engineering.[PhD Thesis].Beijing:Beijing Jiaotong University,2014(in Chinese))

11 程志寶.周期性結(jié)構(gòu)及周期性隔震基礎(chǔ).[博士論文].北京：北京交通大學,2014(Cheng Zhibao.Periodic structures and periodic seismic isolation foundations.[PhD Thesis].Beijing:Beijing Jiaotong University,2014(in Chinese))

12 Cheng ZB,Yan YQ,Menq FY,et al.3D periodic foundation-based structural vibration isolation.Proceedings of the World Congress on Engineering&Computer Scie,2013

13 翟婉明.車輛--軌道耦合動力學.北京：科學出版社,2007(Zhai Wanming.Vehicle-track Coupling Dynamics.Beijing:Science Press,2007(in Chinese))

14 Knothe K.Gleisdynamik.Berlin:Ernst&Sohn,2001

15 葉敏，呂敬，丁千等.復合材料層合板1:1參數(shù)共振的分岔研究.力學學報,2004,36(1):64-71(Ye Min,LJing,Ding Qian,et al.The bifurcation analysis of the laminated composite plate with 1:1 parametrically resonance.Chinese Journal of Theoretical and Applied Mechanics,2004,36(1):64-71(in Chinese))

16 王曉丹.滿足層間連續(xù)條件的層合板理論及其準確和高效的板單元.[博士論文].天津：天津大學,2015(Wang Xiaodan.Laminated composite plate theories accounting for interlaminar continuity and its accurate,efficient plate elements.[PhD Thesis].Tianjin:Tianjin University,2015(in Chinese))

17 Duan TC,Li LX.The uni fi ed solution for a beam of rectangular cross-section with di ff erent higher-order shear deformation models.Latin American Journal of Solids and Structures,2016,13:1716-1737

18 段鐵城，李錄賢.厚板的高階剪切變形理論研究.力學學報,2016,48(5):1096-1113(Duan Tiecheng,Li Luxian.Study on higherorder shear deformation theories of thick-plate.Chinese Journal of Theoretical and Applied Mechanics,2016,48(5):1096-1113(in Chinese))

19 李念,陳普會.復合材料層合板低速沖擊損傷分析的連續(xù)介質(zhì)損傷力學模型.力學學報,2015,47(3):458-470(Li Nian,Chen Puhui.Continuum damage mechanics model for low-velocity impact damage analysis of composite laminates.Chinese Journal of Theoretical and Applied Mechanics,2015,47(3):458-470(in Chinese))

20 徐巍,王立峰,蔣經(jīng)農(nóng).基于應(yīng)變梯度中厚板單元的石墨烯振動研究.力學學報,2015,47(5):751-761(Xu Wei,Wang Lifeng,Jiang Jingnong.Finite element analysis of strain gradient middle thick plate model on the vibration of graphene sheets.Chinese Journal of Theoretical and Applied Mechanics,2015,47(5):751-761(in Chinese))

21 陳建霖,勵爭,儲鵬程.大開口復合材料層合板強度破壞研究.力學學報,2016,48(6):1326-1333(Chen Jianlin Li Zheng,Chu Pengcheng.Strength analysis of fi ber reinforced composite laminates with big cutouts.Chinese Journal of Theoretical and Applied Mechanics,2016,48(6):1326-1333(in Chinese))

22 Chen WQ,Ding HJ.On free vibration of a functionally graded piezoelectric rectangular plate.Acta Mechanica,2002,153(3-4):207-216

24 Yang M,Qiao P.Higher-order impact modeling of sandwich structures with fl exible core.International Journal of Solids&Structures,2005,42(20):5460-5490

25 C?alim FF.Free and forced vibrations of non-uniform composite beams.Composite Structures,2009,88(3):413-423

26 Ganapathi M,Kalyani A,Mondal B,et al.Free vibration analysis of simply supported composite laminated panels.Composite Structures,2009,90(1):100-103

27 李堯臣，亓峰，仲政.功能梯度矩形板的近似理論與解析解.力學學報,2010,42(4):670-681(Li Yaochen,Qi Feng,Zhong Zheng.Approximate theory and analytical solution for graded piezoelectric rectangular plates.Chinese Journal of Theoretical and Applied Mechanics,2010,42(4):670-681(in Chinese))

28 李堯臣，聶國雋，楊昌錦.面內(nèi)功能梯度矩形板的近似理論與解答.力學學報,2013,45(4):560-567(Li Yaochen,Nie Guojun,Yang Changjin.Approximate theory and analytical solution for rectangular plates with in-plane sti ff ness gradient.Chinese Journal of Theoretical and Applied Mechanics,2013,45(4):560-567(in Chinese))

29 Thai HT,Nguyen TK,Vo TP,et al.Analysis of functionally graded sandwich plates using a new fi rst-order shear deformation theory.European Journal of Mechanics-A/Solids,2014,45(2):211-225

30 Pandey S,Pradyumna S.Free vibration of functionally graded sandwich plates in thermal environment using a layerwise theory.European Journal of Mechanics-A/Solids,2015,51:55-66

31 方明鏡.高速鐵路瀝青混凝土軌下基礎(chǔ)結(jié)構(gòu)行為與材料設(shè)計.[博士論文].成都：西南交通大學,2012(Fang MingJing.Structural behavior and mix design for asphalt concrete substructures in highspeed rail.[PhD Thesis].Chengdu:Southwest Jiaotong University,2012(in Chinese))

32 錢振東，劉陽，劉昕依.高速鐵路環(huán)氧瀝青混凝土基礎(chǔ)結(jié)構(gòu)動力響應(yīng)及設(shè)計.東南大學學報(自然科學版),2017,47(1):118-123(Qian Zhendong,Liu Yang,Liu Xinyi.Dynamic response and design of epoxy asphalt concrete substructure for high-speed railway.Journal of Southeast University(Natural Science Edition),2017,47(1):118-123(in Chinese))

33 郭大智，任瑞波.層狀黏彈性體系力學.哈爾濱：哈爾濱工業(yè)大學出版社,2001(Guo Dazhi,Ren Ruibo.Mechanics of Layered Viscoelastic Systems.Harbin:Harbin Engineering University Press,2001(in Chinese))

34 孫訓方,方孝淑，關(guān)來泰.材料力學(第五版).北京：高等教育出版社,2009(Sun Xunfang,Fang Xiaoshu,Guan Laitai.al.Mechanics of Materials(5th Edn).Beijing:Higher Education Press,2009(in Chinese))

35 Durbin F.Numerical inversion of laplace transforms:An efficient improvement to dubner and abate’s method.Computer Journal,2013,17(4):371-376

36 楊坤，梅志遠，李華東.黏彈性復合材料夾芯板穩(wěn)態(tài)響應(yīng)分析.振動與沖擊,2013,32(7):88-92(Yang Kun,Mei Zhiyuan,Li Huadong.Steady response analysis for a composite sandwich plate with viscoelastic core layer based on Kelvin model.Journal of Vibration and Shock,2013,32(7):88-92(in Chinese))

37 曹志遠. 板殼振動理論.北京：中國鐵道出版社,1989(Cao Zhiyuan.Vibration Theory of Plates and Shells.Beijing:China Railway Publishing House,1988(in Chinese))

38 王春玲，黃必成，曹彩芹.移動載荷下彈性半空間地基上矩形薄板的動力響應(yīng)解析研究.土木工程學報,2016,49(1):116-121(Wang Chunling,Huang Bicheng,Cao Caiqin.Analytical study on dynamic response of rectangular thin plate on an elsticc half-space under moving loads.China Civil Engineering Journal,2016,49(1):116-121(in Chinese))

39 楊新文，楊建近.MTMDs對浮置板軌道結(jié)構(gòu)隔振性能的影響分析.鐵道學報,2015(4):87-93(Yang Xinwen,Yang Jianjin.E ff ect of MTMDS on the vibration isolation of fl oating slab track.Journal of The China Railway Society,2015(4):87-93(in Chinese))

APPROXIMATE THEORY AND ANALYTICAL SOLUTION FOR DYNAMIC CALCULATION OF VISCOELASTIC LAYERED PERIODIC PLATE

Xu Sihui*,1)Wang Binglong*Zhou Shunhua*Yang Xinwen*Li Yaochen?
*(Key Laboratory of Road and Traffic Engineering of the Ministry of Education,Tongji University,Shanghai201804,China)
?(School of Aerospace Engineering and Applied Mechanics,Tongji University,Shanghai200092,China)

Rail transit loading and viscoelastic of material are mostly ignored in the previous dynamic calculations of the periodic vibration isolation structures.Approximate theory and analytical solution for viscoelastic layered periodic plate subjected to vertical moving harmonic loading is established,and the viscoelastic of material and the transverse shear deformation are considered.In this theory,Reissner-Mindlin assumption and additional equation of shear deformation are introduced,and the relation between the normal rotation and the shear stress of neutral plane is obtained on the assumption that the transverse shear deformation of the plate’s neutral plane is the overall shear deformation of the cross section.Vibrationgoverningequationofviscoelasticlayeredperiodicplateisproposedaccordingtoequilibriumequations and stress continuity conditions,and vertical displacement in Fourier series is derived as well.The model is validated by the good agreement with solution of the classical laminate model and the fi nite element method(FEM).The results show that:(1)Vibration response at the natural frequency of plate can be signi fi cantly reduced by substituting viscoelastic layered periodic plate for homogeneous one,but vibration ampli fi cation in local low frequency band is aroused as well.(2)The vertical displacement of the plate increases with the increment of the loading velocity,and increase trend slows down once velocity above 300km/h.(3)Shear modulus of viscoelastic layer can be designed to achieve the optimal vibration isolation characteristic.(4)Vibration response is not susceptible to damping characteristic in low frequency band.(5)It’s appropriate to increase the plate length,within the engineering requirement,to improve the vibration isolation performance.

periodic structure,rail transit loading,Reissner-Mindlin assumption,analytical solution,vibration response

O328

A doi：10.6052/0459-1879-17-248

2017–07–07 收稿，2017–09–13 錄用,2017–09–15 網(wǎng)絡(luò)版發(fā)表.

1)徐司慧，博士研究生，主要研究方向：軌道結(jié)構(gòu)減振.E-mail:xsyh1234@163.com

徐司慧,王炳龍,周順華,楊新文,李堯臣.黏彈性層狀周期板動力計算的近似理論與解答.力學學報,2017,49(6):1348-1359

Xu Sihui,Wang Binglong,Zhou Shunhua,Yang Xinwen,Li Yaochen.Approximate theory and analytical solution for dynamic calculation of viscoelastic layered periodic plate.Chinese Journal of Theoretical and Applied Mechanics,2017,49(6):1348-1359