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

*(同濟(jì)大學(xué)道路與交通工程教育部重點(diǎn)實(shí)驗(yàn)室，上海201804)

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

黏彈性層狀周期板動(dòng)力計(jì)算的近似理論與解答

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

*(同濟(jì)大學(xué)道路與交通工程教育部重點(diǎn)實(shí)驗(yàn)室，上海201804)

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

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

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

引言

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

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

軌道交通載荷是軌道交通結(jié)構(gòu)振動(dòng)的主要振源，研究軌道交通載荷下的結(jié)構(gòu)減振技術(shù)對(duì)軌道交通領(lǐng)域的環(huán)境振動(dòng)控制意義重大.本文將周期性結(jié)構(gòu)應(yīng)用于軌道交通減振領(lǐng)域，采用黏彈性層狀周期板作為軌道交通領(lǐng)域板結(jié)構(gòu)的基本形式.為了考慮軌道交通載荷及材料的黏彈性，采用經(jīng)典軌道動(dòng)力學(xué)中輪軌垂向作用力的基本形式，即垂向移動(dòng)簡(jiǎn)諧載荷[13-14]，并將具有黏彈性的層狀周期板視為復(fù)合材料層合板結(jié)構(gòu)，利用層合板理論進(jìn)行動(dòng)力計(jì)算.從既有文獻(xiàn)的調(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]，因此無法準(zhǔn)確計(jì)算板面的撓度分布，難以預(yù)測(cè)層狀周期板的減振效果.

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

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

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

1.1 基本方程

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

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

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

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

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

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

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

幾何方程

黏彈性層物理方程

彈性層物理方程

式中，σ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)并進(jìn)行Laplace變換，可得

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

對(duì)黏彈性層B，有

1.2 移動(dòng)簡(jiǎn)諧載荷下板振動(dòng)位移解析求解

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

式中上標(biāo)(1)～(4)表示板每一層編號(hào).

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

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

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

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

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

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

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

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

式中

當(dāng)板上的載荷為沿板中心線勻速移動(dòng)的單位簡(jiǎn)諧載荷時(shí)

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

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

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

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

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

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

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

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

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

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

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

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

2 模型驗(yàn)證

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

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

表1 驗(yàn)證模型計(jì)算參數(shù)Table 1 Parameters of model

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

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

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

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

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

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

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

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

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

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

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

3 計(jì)算結(jié)果分析

3.1 導(dǎo)納特性分析

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

圖6 板垂向位移導(dǎo)納Fig.6 Vertical displacement admittance of plate

表3 板位移導(dǎo)納峰值Table 3 Vertical displacement admittance peak of plate

圖7 兩種板位移導(dǎo)納的對(duì)比Fig.7 Comparison of vertical displacement admittance by di ff erent plates

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

3.2 載荷移動(dòng)速度的影響

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

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

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

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

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

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

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

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

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

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

圖11 板長(zhǎng)對(duì)隔振性能的影響Fig.11 Vibration isolation performance varies with length

4 結(jié)論

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

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

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

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

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

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

(5)可在滿足工程實(shí)際的情況下適當(dāng)增加板長(zhǎng)，以提高結(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 溫激鴻，王剛，郁殿龍等.聲子晶體振動(dòng)帶隙及減振特性研究.中國(guó)科學(xué),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è)計(jì).[博士論文].大連：大連理工大學(xué),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ù)研究.機(jī)械工程學(xué)報(bào),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)動(dòng)力衰減域特性研究.振動(dòng)與沖擊,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)用研究.[博士論文].北京:北京交通大學(xué),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ǔ).[博士論文].北京：北京交通大學(xué),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 翟婉明.車輛--軌道耦合動(dòng)力學(xué).北京：科學(xué)出版社,2007(Zhai Wanming.Vehicle-track Coupling Dynamics.Beijing:Science Press,2007(in Chinese))

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

15 葉敏，呂敬，丁千等.復(fù)合材料層合板1:1參數(shù)共振的分岔研究.力學(xué)學(xué)報(bào),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ù)條件的層合板理論及其準(zhǔn)確和高效的板單元.[博士論文].天津：天津大學(xué),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 段鐵城，李錄賢.厚板的高階剪切變形理論研究.力學(xué)學(xué)報(bào),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 李念,陳普會(huì).復(fù)合材料層合板低速?zèng)_擊損傷分析的連續(xù)介質(zhì)損傷力學(xué)模型.力學(xué)學(xué)報(bào),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)變梯度中厚板單元的石墨烯振動(dòng)研究.力學(xué)學(xué)報(bào),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 陳建霖,勵(lì)爭(zhēng),儲(chǔ)鵬程.大開口復(fù)合材料層合板強(qiáng)度破壞研究.力學(xué)學(xué)報(bào),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 李堯臣，亓峰，仲政.功能梯度矩形板的近似理論與解析解.力學(xué)學(xué)報(bào),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 李堯臣，聶國(guó)雋，楊昌錦.面內(nèi)功能梯度矩形板的近似理論與解答.力學(xué)學(xué)報(bào),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è)計(jì).[博士論文].成都：西南交通大學(xué),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)動(dòng)力響應(yīng)及設(shè)計(jì).東南大學(xué)學(xué)報(bào)(自然科學(xué)版),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 郭大智，任瑞波.層狀黏彈性體系力學(xué).哈爾濱：哈爾濱工業(yè)大學(xué)出版社,2001(Guo Dazhi,Ren Ruibo.Mechanics of Layered Viscoelastic Systems.Harbin:Harbin Engineering University Press,2001(in Chinese))

34 孫訓(xùn)方,方孝淑，關(guān)來泰.材料力學(xué)(第五版).北京：高等教育出版社,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 楊坤，梅志遠(yuǎn)，李華東.黏彈性復(fù)合材料夾芯板穩(wěn)態(tài)響應(yīng)分析.振動(dò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 曹志遠(yuǎn). 板殼振動(dòng)理論.北京：中國(guó)鐵道出版社,1989(Cao Zhiyuan.Vibration Theory of Plates and Shells.Beijing:China Railway Publishing House,1988(in Chinese))

38 王春玲，黃必成，曹彩芹.移動(dòng)載荷下彈性半空間地基上矩形薄板的動(dòng)力響應(yīng)解析研究.土木工程學(xué)報(bào),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對(duì)浮置板軌道結(jié)構(gòu)隔振性能的影響分析.鐵道學(xué)報(bào),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

徐司慧,王炳龍,周順華,楊新文,李堯臣.黏彈性層狀周期板動(dòng)力計(jì)算的近似理論與解答.力學(xué)學(xué)報(bào),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