• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-12-18 13:24:01徐司慧王炳龍周順華楊新文李堯臣
    力學學報 2017年6期
    關(guān)鍵詞:合板層狀剪切

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

    *(同濟大學道路與交通工程教育部重點實驗室,上海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

    猜你喜歡
    合板層狀剪切
    軋制復合制備TA1/AZ31B/TA1層狀復合材料組織與性能研究
    重型機械(2019年3期)2019-08-27 00:58:44
    寬厚板剪切線控制系統(tǒng)改進
    山東冶金(2018年5期)2018-11-22 05:12:46
    混凝土短梁斜向開裂后的有效剪切剛度與變形
    層合板上層建筑側(cè)壁抗空爆性能研究
    兩級結(jié)構(gòu)層狀Ti-TiBw/Ti復合材料擴散焊及其拉伸行為
    焊接(2016年9期)2016-02-27 13:05:22
    土-混凝土接觸面剪切破壞模式分析
    高韌性抗層狀撕裂Q345FTE-Z35鋼板開發(fā)
    新疆鋼鐵(2015年2期)2015-11-07 03:27:52
    基于玻璃纖維增強隔音復合材料的層合板的隔音性能
    濕熱環(huán)境對CCF300復合材料層合板的載荷放大系數(shù)影響
    單釘機械連接孔邊應(yīng)力及失效分析
    教練機(2014年3期)2014-12-02 06:16:02
    老司机亚洲免费影院| 视频区欧美日本亚洲| tube8黄色片| 老熟妇仑乱视频hdxx| av免费在线观看网站| 午夜91福利影院| av天堂久久9| 日韩一卡2卡3卡4卡2021年| 国产麻豆69| 交换朋友夫妻互换小说| 性色av乱码一区二区三区2| 国产一区二区三区在线臀色熟女 | 国产激情欧美一区二区| 夜夜夜夜夜久久久久| 精品国产一区二区三区久久久樱花| 亚洲成人国产一区在线观看| 久久狼人影院| 亚洲欧美色中文字幕在线| 黑人巨大精品欧美一区二区蜜桃| 成人av一区二区三区在线看| 精品一区二区三区视频在线观看免费 | 久久精品国产综合久久久| 久久人人爽av亚洲精品天堂| 免费久久久久久久精品成人欧美视频| 高清av免费在线| 日韩一卡2卡3卡4卡2021年| 男人操女人黄网站| 国产熟女午夜一区二区三区| 69精品国产乱码久久久| 久久精品亚洲av国产电影网| 久久人妻熟女aⅴ| 亚洲精品中文字幕一二三四区| 国产一区二区三区在线臀色熟女 | 午夜福利乱码中文字幕| 国产精品香港三级国产av潘金莲| 两个人免费观看高清视频| 国产一区有黄有色的免费视频| 极品少妇高潮喷水抽搐| 在线视频色国产色| 久久久精品免费免费高清| 人妻 亚洲 视频| 老熟女久久久| 国产深夜福利视频在线观看| 午夜免费成人在线视频| 欧美av亚洲av综合av国产av| 精品一区二区三卡| 国产又色又爽无遮挡免费看| 国产主播在线观看一区二区| 日韩欧美一区二区三区在线观看 | 人妻丰满熟妇av一区二区三区 | 欧美不卡视频在线免费观看 | 国产亚洲精品久久久久5区| 国产精品一区二区在线观看99| 香蕉丝袜av| 国内毛片毛片毛片毛片毛片| 久久中文字幕一级| 国产成人精品在线电影| 亚洲性夜色夜夜综合| 色94色欧美一区二区| 亚洲av成人av| 人人妻人人澡人人爽人人夜夜| 日韩三级视频一区二区三区| 日韩三级视频一区二区三区| av有码第一页| 欧美日韩亚洲综合一区二区三区_| 99re6热这里在线精品视频| 一a级毛片在线观看| 亚洲片人在线观看| 欧美午夜高清在线| 欧美日本中文国产一区发布| 免费在线观看完整版高清| 99国产精品免费福利视频| 人人妻人人澡人人爽人人夜夜| 免费人成视频x8x8入口观看| 欧美成狂野欧美在线观看| 成人国语在线视频| 欧美最黄视频在线播放免费 | 交换朋友夫妻互换小说| 午夜福利在线观看吧| 亚洲人成电影观看| 精品亚洲成a人片在线观看| 成人av一区二区三区在线看| 久久人人97超碰香蕉20202| 十八禁高潮呻吟视频| av在线播放免费不卡| 国产精品国产av在线观看| 国产精品免费一区二区三区在线 | 人人妻人人澡人人爽人人夜夜| 美国免费a级毛片| 亚洲中文日韩欧美视频| 亚洲成人手机| 丁香六月欧美| 国产在视频线精品| 久久久久国产精品人妻aⅴ院 | 精品国产乱码久久久久久男人| 久久人妻福利社区极品人妻图片| 欧美激情高清一区二区三区| 高清黄色对白视频在线免费看| 久久精品成人免费网站| 麻豆国产av国片精品| 操美女的视频在线观看| 亚洲精品国产一区二区精华液| 高清在线国产一区| 成人18禁在线播放| 国产aⅴ精品一区二区三区波| 9色porny在线观看| 啦啦啦视频在线资源免费观看| 9191精品国产免费久久| 国产精品98久久久久久宅男小说| 欧美日韩黄片免| 新久久久久国产一级毛片| 巨乳人妻的诱惑在线观看| avwww免费| 欧美午夜高清在线| 精品人妻熟女毛片av久久网站| 亚洲av日韩精品久久久久久密| 午夜亚洲福利在线播放| 一区二区三区精品91| 国产精品秋霞免费鲁丝片| 久久草成人影院| 亚洲人成电影免费在线| 久99久视频精品免费| 18禁美女被吸乳视频| 高清毛片免费观看视频网站 | 亚洲人成77777在线视频| 久久99一区二区三区| 男人的好看免费观看在线视频 | 视频区欧美日本亚洲| 久久精品91无色码中文字幕| 精品熟女少妇八av免费久了| 国产极品粉嫩免费观看在线| 在线av久久热| 国产高清国产精品国产三级| xxxhd国产人妻xxx| 一本大道久久a久久精品| 日本黄色视频三级网站网址 | 女人被狂操c到高潮| 久久影院123| 天天躁日日躁夜夜躁夜夜| 色老头精品视频在线观看| 激情在线观看视频在线高清 | 免费看十八禁软件| bbb黄色大片| 狂野欧美激情性xxxx| av天堂久久9| 精品卡一卡二卡四卡免费| 97人妻天天添夜夜摸| 纯流量卡能插随身wifi吗| 久久国产精品男人的天堂亚洲| 日韩三级视频一区二区三区| 久久久久久久精品吃奶| 制服诱惑二区| 大片电影免费在线观看免费| 亚洲av片天天在线观看| 久久国产精品影院| 国产男靠女视频免费网站| 久久天堂一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 中国美女看黄片| 免费在线观看亚洲国产| 国产高清国产精品国产三级| 搡老熟女国产l中国老女人| 国产亚洲一区二区精品| 亚洲第一av免费看| 国产午夜精品久久久久久| 天天操日日干夜夜撸| 免费观看a级毛片全部| 午夜久久久在线观看| 91精品国产国语对白视频| 一级a爱视频在线免费观看| 99国产极品粉嫩在线观看| 人妻 亚洲 视频| 香蕉丝袜av| 麻豆国产av国片精品| 国产精品欧美亚洲77777| 国产成+人综合+亚洲专区| 精品一区二区三区视频在线观看免费 | 老司机午夜十八禁免费视频| www.精华液| 精品高清国产在线一区| 在线观看www视频免费| 久久亚洲精品不卡| 国产精品av久久久久免费| 免费观看a级毛片全部| 国产精品欧美亚洲77777| 香蕉久久夜色| 亚洲欧美日韩另类电影网站| 精品亚洲成国产av| 男女床上黄色一级片免费看| 满18在线观看网站| 99热国产这里只有精品6| 国产激情久久老熟女| 免费一级毛片在线播放高清视频 | 最新的欧美精品一区二区| 一级作爱视频免费观看| 成熟少妇高潮喷水视频| 中文字幕高清在线视频| 日韩中文字幕欧美一区二区| 国产亚洲精品第一综合不卡| 亚洲精品粉嫩美女一区| 亚洲国产欧美一区二区综合| 99久久精品国产亚洲精品| av国产精品久久久久影院| 下体分泌物呈黄色| 伊人久久大香线蕉亚洲五| 另类亚洲欧美激情| 精品国产乱码久久久久久男人| 麻豆国产av国片精品| 法律面前人人平等表现在哪些方面| 欧美日韩国产mv在线观看视频| 婷婷成人精品国产| 国产欧美日韩一区二区三| 日本a在线网址| 国产成人av教育| 久久热在线av| a级毛片黄视频| 欧美日韩亚洲高清精品| 国产免费av片在线观看野外av| 国产av精品麻豆| 国产熟女午夜一区二区三区| 中文字幕最新亚洲高清| 欧美不卡视频在线免费观看 | 丰满人妻熟妇乱又伦精品不卡| 国产成人精品久久二区二区免费| 亚洲五月色婷婷综合| 精品久久久久久,| videosex国产| 久久精品91无色码中文字幕| 亚洲精品乱久久久久久| 亚洲av片天天在线观看| 亚洲国产欧美日韩在线播放| 欧美老熟妇乱子伦牲交| 一边摸一边抽搐一进一小说 | 亚洲伊人色综图| 日本精品一区二区三区蜜桃| 亚洲精华国产精华精| 欧美一级毛片孕妇| 欧美性长视频在线观看| 欧美+亚洲+日韩+国产| 国产精品二区激情视频| 亚洲aⅴ乱码一区二区在线播放 | 久久久精品免费免费高清| 欧美激情高清一区二区三区| 免费在线观看亚洲国产| 韩国av一区二区三区四区| 999久久久国产精品视频| 十八禁网站免费在线| 两性夫妻黄色片| 正在播放国产对白刺激| 亚洲午夜理论影院| 飞空精品影院首页| 人妻一区二区av| 久久精品91无色码中文字幕| 日本一区二区免费在线视频| 久久久久国内视频| 欧美黄色片欧美黄色片| 热re99久久精品国产66热6| 国产精品电影一区二区三区 | 啦啦啦在线免费观看视频4| 777米奇影视久久| 日韩免费高清中文字幕av| 丝袜在线中文字幕| 免费在线观看完整版高清| 亚洲国产欧美一区二区综合| 久久久久国产一级毛片高清牌| 精品高清国产在线一区| 建设人人有责人人尽责人人享有的| 国产精品国产av在线观看| 电影成人av| √禁漫天堂资源中文www| 搡老熟女国产l中国老女人| 老司机靠b影院| 亚洲熟妇中文字幕五十中出 | 在线观看午夜福利视频| www.精华液| 国产精品1区2区在线观看. | 黑丝袜美女国产一区| 69精品国产乱码久久久| 国内久久婷婷六月综合欲色啪| 美国免费a级毛片| 国产男靠女视频免费网站| 成人免费观看视频高清| 欧美成狂野欧美在线观看| www.精华液| 欧美激情久久久久久爽电影 | 两人在一起打扑克的视频| 啦啦啦视频在线资源免费观看| 韩国精品一区二区三区| 色94色欧美一区二区| 国产麻豆69| 9191精品国产免费久久| 久久久久久久国产电影| 99在线人妻在线中文字幕 | 久热爱精品视频在线9| 亚洲国产精品sss在线观看 | 亚洲少妇的诱惑av| 丝袜美腿诱惑在线| 欧美一级毛片孕妇| 一本大道久久a久久精品| 国产真人三级小视频在线观看| 成在线人永久免费视频| 国产97色在线日韩免费| 欧美 亚洲 国产 日韩一| 国产不卡一卡二| 国产精品 国内视频| 午夜久久久在线观看| 变态另类成人亚洲欧美熟女 | 大片电影免费在线观看免费| 欧美 亚洲 国产 日韩一| 精品亚洲成国产av| 天堂中文最新版在线下载| 99久久综合精品五月天人人| 热99国产精品久久久久久7| 久久国产精品男人的天堂亚洲| 中文字幕人妻丝袜制服| 大码成人一级视频| 亚洲色图av天堂| www.999成人在线观看| 麻豆av在线久日| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲一区二区精品| 亚洲成av片中文字幕在线观看| 麻豆国产av国片精品| 色综合婷婷激情| 免费看a级黄色片| 狠狠狠狠99中文字幕| 大片电影免费在线观看免费| 国产亚洲av高清不卡| 久久精品国产清高在天天线| 成人黄色视频免费在线看| 看片在线看免费视频| 精品亚洲成a人片在线观看| 久久精品国产亚洲av高清一级| 久久草成人影院| 国产精品免费视频内射| 99热只有精品国产| 久久性视频一级片| 99精品在免费线老司机午夜| 乱人伦中国视频| 久久人妻熟女aⅴ| 日韩欧美在线二视频 | а√天堂www在线а√下载 | 在线观看一区二区三区激情| а√天堂www在线а√下载 | 黑人巨大精品欧美一区二区蜜桃| 大陆偷拍与自拍| 18禁裸乳无遮挡动漫免费视频| 老司机深夜福利视频在线观看| 国产又色又爽无遮挡免费看| 国产色视频综合| 热re99久久国产66热| 亚洲国产毛片av蜜桃av| 亚洲五月天丁香| av网站在线播放免费| 免费在线观看视频国产中文字幕亚洲| 亚洲国产看品久久| av超薄肉色丝袜交足视频| 丰满人妻熟妇乱又伦精品不卡| 国产成人欧美| 99国产综合亚洲精品| 黑人猛操日本美女一级片| 日本黄色视频三级网站网址 | 精品第一国产精品| 成人av一区二区三区在线看| 欧美日韩亚洲综合一区二区三区_| 日韩欧美在线二视频 | 精品亚洲成a人片在线观看| 久久国产精品大桥未久av| 久久草成人影院| 亚洲精品成人av观看孕妇| 亚洲精品美女久久av网站| 国产成人av教育| 成年版毛片免费区| 国产欧美日韩一区二区三| 亚洲人成电影免费在线| 久久精品国产综合久久久| 成年动漫av网址| 国产欧美日韩一区二区精品| 啦啦啦在线免费观看视频4| 新久久久久国产一级毛片| 一级片免费观看大全| 脱女人内裤的视频| 免费看a级黄色片| 老鸭窝网址在线观看| 最新美女视频免费是黄的| av天堂在线播放| 亚洲一码二码三码区别大吗| 精品第一国产精品| 午夜精品国产一区二区电影| 少妇裸体淫交视频免费看高清 | 国产亚洲一区二区精品| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站| 成年人午夜在线观看视频| 精品人妻熟女毛片av久久网站| 免费在线观看亚洲国产| 国产乱人伦免费视频| 久久精品人人爽人人爽视色| 多毛熟女@视频| 18禁观看日本| 亚洲av日韩在线播放| 一级毛片女人18水好多| 久久国产精品男人的天堂亚洲| 如日韩欧美国产精品一区二区三区| 久久影院123| 日韩中文字幕欧美一区二区| 久久精品人人爽人人爽视色| 国产成人啪精品午夜网站| 国产成人欧美| 99精国产麻豆久久婷婷| 美女扒开内裤让男人捅视频| 日韩中文字幕欧美一区二区| 香蕉丝袜av| 亚洲精品久久午夜乱码| 亚洲欧美色中文字幕在线| 亚洲熟女精品中文字幕| 亚洲午夜精品一区,二区,三区| 成年人免费黄色播放视频| 日韩欧美在线二视频 | 嫁个100分男人电影在线观看| 美国免费a级毛片| 久久国产乱子伦精品免费另类| 欧美日韩亚洲综合一区二区三区_| 国内久久婷婷六月综合欲色啪| 激情在线观看视频在线高清 | 首页视频小说图片口味搜索| 天堂俺去俺来也www色官网| 黄频高清免费视频| 淫妇啪啪啪对白视频| 精品亚洲成国产av| 欧美人与性动交α欧美软件| 精品福利观看| 美女高潮到喷水免费观看| 欧美另类亚洲清纯唯美| 国产欧美日韩一区二区三区在线| 国产精品久久视频播放| 精品卡一卡二卡四卡免费| 国产蜜桃级精品一区二区三区 | 午夜视频精品福利| 亚洲精华国产精华精| 好看av亚洲va欧美ⅴa在| 亚洲男人天堂网一区| 人妻久久中文字幕网| 日韩免费av在线播放| 亚洲av熟女| 一级毛片高清免费大全| 午夜免费观看网址| 色综合欧美亚洲国产小说| 18禁观看日本| 久久久久国内视频| 日韩免费高清中文字幕av| 亚洲精品国产精品久久久不卡| 夫妻午夜视频| 亚洲精品在线观看二区| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 岛国毛片在线播放| 欧美精品av麻豆av| 天天影视国产精品| 国产精品秋霞免费鲁丝片| 身体一侧抽搐| 久久精品成人免费网站| 又黄又粗又硬又大视频| 国产一区二区三区综合在线观看| 免费久久久久久久精品成人欧美视频| 久久久久久久国产电影| 又大又爽又粗| 久久精品亚洲精品国产色婷小说| 在线天堂中文资源库| 国产aⅴ精品一区二区三区波| 国产真人三级小视频在线观看| 视频区图区小说| 亚洲国产欧美日韩在线播放| 亚洲av熟女| 亚洲一区二区三区不卡视频| 亚洲精品美女久久久久99蜜臀| 久久久国产一区二区| 99久久人妻综合| 三级毛片av免费| 亚洲欧美激情在线| 高清黄色对白视频在线免费看| 色综合婷婷激情| 99re6热这里在线精品视频| 精品午夜福利视频在线观看一区| 国产精品影院久久| 操出白浆在线播放| 丝袜在线中文字幕| 成人手机av| 亚洲五月色婷婷综合| 亚洲黑人精品在线| 国产在视频线精品| 久久久水蜜桃国产精品网| 久久午夜综合久久蜜桃| 80岁老熟妇乱子伦牲交| 中出人妻视频一区二区| 国产日韩欧美亚洲二区| 国产国语露脸激情在线看| 丝袜在线中文字幕| 亚洲熟妇中文字幕五十中出 | 多毛熟女@视频| 搡老岳熟女国产| 午夜精品久久久久久毛片777| 建设人人有责人人尽责人人享有的| 国产精品 国内视频| 热99国产精品久久久久久7| 亚洲九九香蕉| 欧美最黄视频在线播放免费 | 国产亚洲精品久久久久久毛片 | 国产精品免费视频内射| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频 | 欧美日韩国产mv在线观看视频| 亚洲成av片中文字幕在线观看| 国产精品久久电影中文字幕 | cao死你这个sao货| 精品卡一卡二卡四卡免费| 色94色欧美一区二区| 国产激情欧美一区二区| 国产欧美亚洲国产| 亚洲精品久久午夜乱码| 亚洲熟妇中文字幕五十中出 | 亚洲在线自拍视频| 人妻一区二区av| 久久午夜综合久久蜜桃| 黄片小视频在线播放| 久久影院123| 亚洲,欧美精品.| 搡老乐熟女国产| 久久精品国产清高在天天线| 国产精品久久电影中文字幕 | 最近最新中文字幕大全免费视频| 欧美性长视频在线观看| 午夜亚洲福利在线播放| 高清欧美精品videossex| 精品久久久久久久毛片微露脸| 99久久99久久久精品蜜桃| 天天添夜夜摸| 国产成人精品久久二区二区免费| 成人免费观看视频高清| 精品国产国语对白av| 极品教师在线免费播放| 久久精品国产综合久久久| 久久中文字幕人妻熟女| 欧美乱码精品一区二区三区| 精品久久久久久久毛片微露脸| 成人精品一区二区免费| 午夜福利免费观看在线| 王馨瑶露胸无遮挡在线观看| 亚洲九九香蕉| 757午夜福利合集在线观看| 欧美在线一区亚洲| 97人妻天天添夜夜摸| 亚洲精品自拍成人| 亚洲中文字幕日韩| 两性夫妻黄色片| 精品一区二区三区视频在线观看免费 | 国产精品1区2区在线观看. | svipshipincom国产片| 另类亚洲欧美激情| 久久午夜综合久久蜜桃| 黄片小视频在线播放| 人人澡人人妻人| 19禁男女啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 老司机影院毛片| 一本大道久久a久久精品| 亚洲欧美一区二区三区黑人| 欧美激情久久久久久爽电影 | 亚洲熟女毛片儿| 嫁个100分男人电影在线观看| 国产无遮挡羞羞视频在线观看| 国产欧美日韩一区二区精品| 久久久久久久国产电影| 热99re8久久精品国产| 十八禁高潮呻吟视频| 久久久国产成人精品二区 | 国产蜜桃级精品一区二区三区 | 精品一区二区三卡| 色综合欧美亚洲国产小说| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人免费av一区二区三区 | 欧美精品人与动牲交sv欧美| av免费在线观看网站| 久久精品亚洲av国产电影网| 美女视频免费永久观看网站| 成人黄色视频免费在线看| 欧美日韩亚洲高清精品| 亚洲人成伊人成综合网2020| 国产成人系列免费观看| 老司机福利观看| 亚洲,欧美精品.| 日日摸夜夜添夜夜添小说| 亚洲九九香蕉| 777久久人妻少妇嫩草av网站| 久久精品成人免费网站| 国产一区二区三区视频了| 亚洲熟女精品中文字幕| 亚洲精品在线观看二区| av片东京热男人的天堂| 国产xxxxx性猛交| 成人18禁高潮啪啪吃奶动态图| 亚洲熟妇中文字幕五十中出 | 在线观看舔阴道视频| 久久中文字幕一级| 九色亚洲精品在线播放| 超色免费av| 国产99白浆流出| 色婷婷av一区二区三区视频| 在线观看免费视频网站a站| 淫妇啪啪啪对白视频| 老司机影院毛片| 黄色成人免费大全| 欧美日韩一级在线毛片| 国产在线观看jvid| 国产精品一区二区在线不卡|