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

    任意邊界條件下旋轉功能梯度錐-柱連接殼行波模態(tài)分析

    2024-01-19 06:57:08張宇航劉文光
    船舶力學 2024年1期
    關鍵詞:環(huán)向行波邊界條件

    張宇航,劉文光,劉 超

    (1.南昌航空大學航空制造工程學院,南昌 330063;2.哈爾濱工業(yè)大學深圳理學院,廣東深圳 150001)

    0 引 言

    作為一種新型復合材料,功能梯度材料(functionally graded materials,F(xiàn)GMs)有望應用于船舶推進系統(tǒng)中連接殼的設計。但是連接殼通常服役于復雜的邊界條件環(huán)境中并伴隨著旋轉運動,導致連接殼結構出現(xiàn)行波模態(tài),所以研究任意邊界下旋轉FGMs連接殼的行波模態(tài)對于推進FGMs在船舶推進結構的動力學設計應用具有重要價值。

    近些年來,關于FGMs 殼的振動問題引起了研究者的廣泛關注?;赗ayleigh-Ritz 法,劉超等[1-2]推導了內環(huán)向加筋FGMs圓柱殼和含孔隙FGMs圓柱殼的模態(tài)頻率方程,分析了加筋方式、數(shù)量、位置及孔隙率對FGMs 圓柱殼模態(tài)頻率的影響。引入彈性支撐,黃小林等[3]采用應力函數(shù)法推導了FGMs圓錐薄殼的運動學方程,分析了含孔隙的FGMs 圓錐薄殼的振動響應。受實際因素的影響,經典邊界通常難以完全實現(xiàn),眾多學者開始致力于研究任意邊界條件下殼體結構的振動問題[4-7]。引入人工彈簧技術,Qin 等[7]以改進Fourier 級數(shù)、正交多項式和Chebyshev 多項式構造位移函數(shù)求解了圓柱殼的自由振動,比較了三者的收斂速度和計算效率,發(fā)現(xiàn)Chebyshev 多項式具有較高的計算速度??紤]旋轉運動時殼體產生的科氏力和離心力,Bryan[8]研究了自旋轉圓柱殼的動力學行為。結合改進的Fourier級數(shù)法和Rayleigh-Ritz 法,李文達等[9-10]分析了彈性邊界下旋轉薄壁圓柱殼的行波模態(tài)特性。利用廣義微分求積方法(GDQM)[11],Han等[12]分析了旋轉FGMs圓錐殼的動力學特性;Shakouri[13]研究了熱環(huán)境下旋轉FGMs 圓錐殼的自由振動。采用Chebyshev 多項式描述位移容許函數(shù),劉超等[14]分析了任意邊界條件下旋轉FGMs層合圓柱殼的行波模態(tài)特性,分析表明,旋轉運動對殼體結構的模態(tài)有很大影響,科氏力導致旋轉殼行波模態(tài)出現(xiàn)分叉現(xiàn)象。

    工程實際中,殼體結構往往以連接殼的形式出現(xiàn)。因此,連接殼結構的振動問題也得到了研究者的關注。Irie等[15]采用傳遞矩陣法,建立了圓錐-圓柱連接殼體的模態(tài)頻率方程,求解了不同幾何參數(shù)下連接殼的模態(tài)頻率。Bagheri等[16-19]研究了不同邊界條件下FGMs圓柱殼、圓錐殼和球殼各種組合殼的模態(tài)頻率。采用分區(qū)廣義變分和最小二乘加權殘值法,瞿葉高等[20-21]將振動問題轉化為滿足條件下的無約束泛函變分問題,分析了以圓錐殼、圓柱殼和球殼為單元的單殼與連接殼的自由振動問題。采用改進的Fourier 級數(shù)法,張帥等[22]研究了錐-柱-球組合連接殼的振動行為??紤]石墨烯增強復合材料和碳納米管增強材料對殼體結構的作用,研究者探討了增強復合材料對殼體結構動力學性能的影響。Soureshjani 等[23]分析了熱效應對碳納米管增強圓錐-圓錐連接殼動力學行為的影響;Damercheloo等[24]研究了不同邊界條件下石墨烯增強復合材料圓錐-圓錐連接殼的自由振動特性。

    雖然研究者對連接殼的振動問題進行了廣泛的研究,但大多數(shù)研究通常關注靜態(tài)結構下的駐波模態(tài)。由于旋轉運動下,離心力會產生環(huán)向應力,提高結構的剛度;同時,科氏力會導致頻率發(fā)生分叉現(xiàn)象,產生前后行波模態(tài),并且由于前后行波頻率不同,無法再疊加為常規(guī)的駐波模態(tài),因此研究旋轉結構的振動問題,只能從行波模態(tài)的角度出發(fā)。本文以旋轉FGMs錐-柱連接殼為對象,采用彈簧模擬結構的任意邊界條件和接觸邊界條件,推導連接殼的模態(tài)頻率方程,探討不同參數(shù)對旋轉連接殼行波模態(tài)的影響。研究結果可以為FGMs連接殼的動力學設計提供理論支撐。

    1 FGMs連接殼模型

    1.1 幾何模型

    如圖1 所示,F(xiàn)GMs 連接殼由圓錐殼和圓柱殼組成。假設兩個柱坐標系(xi,θi,zi)分別建立在圓錐殼和圓柱殼中面上,殼體的厚度均為h,并且殼體繞軸線以恒定角速度Ω旋轉。殼體的軸向、環(huán)向和法向的位移分別為ui、vi、wi,其下標con、cy分別表示圓錐殼和圓柱殼。以該柱坐標系為參照,圓錐殼小端中面半徑為R1,大端中面半徑為R2,長度為Lcon,半錐角為α0。相應的圓柱殼中面半徑為Rcy=R2,長度為Lcy。

    圖1 旋轉FGMs圓錐-圓柱連接殼的幾何模型Fig.1 Geometric model of a rotating FGMs joined conical-cylindrical shell

    1.2 材料屬性模型

    圓錐殼和圓柱殼均采用金屬陶瓷材料作為梯度組分,采用冪律函數(shù)來描述陶瓷材料沿厚度方向上的體積分數(shù)[13]為

    采用Voigt模型,F(xiàn)GMs沿殼體厚度方向的有效參數(shù)可以定義為

    式中,N表示陶瓷體積分數(shù)指數(shù),Vmi、Vci分別為金屬和陶瓷的體積分數(shù),Pmi、Pci分別為金屬和陶瓷的材料屬性,Pi為FGMs圓錐殼和圓柱殼的有效材料屬性。

    因此,F(xiàn)GMs連接殼的材料屬性如彈性模量E、密度ρ、泊松比v可分別表示為

    式中,下標m、c分別表示金屬材料和陶瓷材料。

    1.3 邊界連續(xù)性模型

    如圖2所示,在結構的兩端與圓柱-圓錐殼的連接處分別定義了4組連續(xù)分布的彈簧,以模擬任意邊界條件和連接處的連續(xù)性。在xcon=0和xcy=Lcy處設置軸向彈簧ki

    圖2 邊界約束彈簧和連接彈簧Fig.2 Boundary constraint springs and connecting springs

    1、環(huán)向剪切彈簧ki2、徑向剪切彈簧ki3與扭轉彈簧ki4。同樣在xcon=Lcon和xcy=0 處設置軸向彈簧kcon-cy1、環(huán)向剪切彈簧kcon-cy2、徑向剪切彈簧kcon-cy3、扭轉彈簧kcon-cy4。連接處的位移所滿足的連續(xù)性條件可表示為

    2 模態(tài)頻率方程

    2.1 應力應變方程

    根據(jù)Love薄殼理論,連接殼的應變分量與中面位移的關系[22]可定義為

    2.2 能量方程

    根據(jù)彈性力學,F(xiàn)GMs連接殼的應變能Us可以表示為

    殼體旋轉時產生的動能K由方程(17)和(18)表示,環(huán)向應力產生的初始應變能Uh由方程(19)和(20)表達。根據(jù)1.3 節(jié)中連續(xù)性條件,儲存在邊界彈簧和接觸彈簧中的彈性勢能Ue可由方程(21)~(23)表達。

    式中,上標“·”表示對時間的一階導數(shù);Ω的0次項表示因變形產生的動能;Ω的1次項表示因科氏力產生的動能;Ω的2次項表示因離心力所產生的動能;Ii

    1表示慣性矩,表達式為

    需要指出的是,當令轉速為0 時,結構簡化為靜止結構,其同樣適用于本文的方法進行模態(tài)分析。

    2.3 位移方程

    假設環(huán)向波數(shù)取n時的位移方程為

    用Chebyshev多項式展開振型函數(shù)以模擬任意邊界條件和接觸連續(xù)性條件:

    2.4 頻率方程

    構造拉格朗日能量函數(shù)

    式中,M1、M2、K分別表示質量陣和剛度陣。

    求解式(32)可得到旋轉連接殼的前后行波模態(tài)頻率。與旋轉方向一致且小于0為前行波模態(tài)頻率,相反為后行波模態(tài)頻率。為了更好地比較前后行波模態(tài),取前后行波模態(tài)頻率的絕對值分析。定義無量綱轉速和無量綱頻率:

    式中,ωb和ωf分別表示連接殼的后行波模態(tài)頻率和前行波模態(tài)頻率。

    3 模型驗證

    3.1 收斂性分析

    假設FGMs 連接殼的幾何參數(shù)分別?。篟1/Rcy=0.4,Lcy/Rcy=2.5,Rcy/h=100,α0=30°。材料參數(shù)分別取Em=2.08×1011Pa,Ec=3.22×1011Pa,ρm=8166 kg/m3,ρc=2370 kg/m3,vm=0.3177,vc=0.3。軸向波數(shù)m取1。如無特殊說明,以下數(shù)值計算取的參數(shù)不發(fā)生改變。

    因模態(tài)頻率與Chebyshev 多項式的截斷項數(shù)以及彈簧剛度取值有關,圖3 首先分析了Ω=0 時,連接彈簧和邊界彈簧剛度值的收斂情況。結果表明,無論是連接彈簧還是邊界彈簧,彈簧剛度值大于1012時,模態(tài)頻率收斂。

    圖3 彈簧剛度收斂性分析(A=6)Fig.3 Convergence analysis of spring stiffness(A=6)

    本文連接彈簧剛度值取1013,而邊界彈簧剛度值為0時視作自由邊界;邊界彈簧剛度為1013時視作固支邊界。因此,改變彈簧的剛度值可模擬不同的邊界條件。表1所示是自由(F)、簡支(S)、固支(C)各種經典邊界對應的彈簧剛度取值。表2分析了自由邊界條件下連接殼(R1/Rcy=0.4226,Rcy/Lcy=1,Rcy/h=100,α0=30°,Em=2.11×1011Pa,ρm=7800 kg/m3,vm=0.3)的無量綱模態(tài)頻率?。

    表1 經典邊界下彈簧剛度值Tab.1 Spring stiffness values with classical boundary conditions

    表2 自由邊界條件下連接殼頻率收斂性分析Tab.2 Convergence analysis of joined shell with free boundary conditions

    結果表明,殼的模態(tài)頻率隨截斷項數(shù)的增加逐漸趨于收斂。所以,后續(xù)分析過程中取截斷項值Ai=10。

    3.2 有效性分析

    將連接殼簡化為靜態(tài)均質圓柱殼、旋轉均質圓柱殼和旋轉FGMs圓錐殼以驗證模型的有效性。表3 對比了不同邊界均質圓柱殼(Lcy/Rcy=2,Rcy/h=500,Em=7.102×1010Pa,vm=0.3,ρm=2796 kg/m3)的模態(tài)頻率。表4 顯示了簡支邊界下旋轉均質圓柱殼(Lcy/Rcy=5,R2/h=500,Ωd=0.0013,Em=1.6806×1011Pa,vm=0.3,ρm=3000 kg/m3)的前行波模態(tài)頻率變化情況。通過與文獻[12]微分求積法(DQM)和文獻[28]廣義微分求積法(GDQM)所得結果的對比,表5 分析了旋轉FGMs圓錐殼(Lcon/R1=2.5,R1/h=20,Ωd=0.025,n=1,α0=45°,N=1)的后行波模態(tài)頻率隨圓錐角的變化情況。

    表3 均質圓柱殼模態(tài)頻率對比Tab.3 Comparison of frequencies of an isotropic cylindrical shell

    表4 旋轉均質圓柱殼前行波頻率對比Tab.4 Comparison of frequencies of forward wave of a rotating isotropic cylindrical shell

    表5 旋轉FGMs圓錐殼后行波頻率對比Tab.5 Comparison of frequencies of backward wave of a rotating FGMs conical shell

    結果表明,本文計算結果和文獻吻合性良好,說明了理論模型的合理性。同時,對于不同邊界條件只需要改變彈簧的剛度取值,就可大大降低構造滿足邊界條件的位移函數(shù)的復雜性。

    4 行波模態(tài)頻率分析

    4.1 環(huán)向波數(shù)的影響

    圖4 比較了Ωd=0.003 時不同邊界條件下旋轉FGMs連接殼行波模態(tài)隨環(huán)向波數(shù)的變化。結果表明,后行波頻率始終大于前行波頻率,兩種頻率均隨著波數(shù)的增大呈先減小后上升趨勢。在C-C 和S-S 邊界下,結構最低的頻率都出現(xiàn)在環(huán)向波數(shù)取5 的情況,而對于F-C邊界,則出現(xiàn)在取值為3時。后續(xù)研究中,為了避免結構發(fā)生低頻共振,選用對應最小頻率下的波數(shù)進行分析。

    圖4 不同的環(huán)向波數(shù)下的行波模態(tài)頻率Fig.4 Traveling wave frequencies with different circumferential wave numbers

    4.2 體積分數(shù)指數(shù)的影響

    圖5分別比較了Ωd=0.001、0.003、0.005時,體積分數(shù)指數(shù)對FGMs錐-柱連接殼前后行波模態(tài)頻率的影響。結果表明,前后行波模態(tài)頻率都隨體積分數(shù)指數(shù)N的增大呈下降趨勢。這是由于N越大,殼體中的陶瓷含量降低,結構剛度下降。隨著轉速的增大,后行波模態(tài)頻率上升,前行波模態(tài)頻率下降,前后行波的分叉越發(fā)明顯。在F-C邊界條件下,前后行波模態(tài)頻率隨轉速和體積分布指數(shù)變化最小??梢园l(fā)現(xiàn),當N>5 后,頻率變化對陶瓷體積分數(shù)的敏感度進一步降低,即對于旋轉FGM 結構,通過調節(jié)N的大小可以在不改變結構模態(tài)頻率的基礎上,得到工程中所需要的材料性能,如陶瓷含量增加所提高的耐熱性、金屬含量增加所提高的耐腐蝕性。

    圖5 不同體積分數(shù)下的行波模態(tài)頻率Fig.5 Traveling wave frequencies with different volume fraction exponents

    4.3 旋轉速度的影響

    圖6研究了圓錐角α0=30°、45°、60°時轉速對前后行波模態(tài)頻率的影響。結果表明,不同邊界條件下,模態(tài)頻率對于錐角的敏感程度為:F-C>S-S>C-C。在F-C 邊界條件下,前行波的模態(tài)頻率隨轉速的增大而下降,并將出現(xiàn)零頻率。因此,實際中考慮F-C邊界條件下,轉速不宜取過大,從而避免結構失穩(wěn)。不同于F-C邊界,C-C和S-S邊界條件下,前行波模態(tài)頻率隨轉速增大先下降后上升。

    圖6 不同旋轉速度下的行波模態(tài)頻率Fig.6 Traveling wave frequencies with different rotational speeds

    4.4 邊界條件的影響

    工程實際中,結構所處的邊界非常復雜,因此接下來分析任意邊界下連接殼的行波模態(tài)。取Ωd=0.003,圖7 和圖8 分別研究了k1i=k4i=0 時,由F-F 邊界變化到S-S 邊界連接殼前后行波模態(tài)頻率的變化和k2i=k3i=1013時,由S-S 邊界變化到C-C 邊界連接殼前后行波模態(tài)頻率的變化。結果表明,環(huán)向彈簧k2i和徑向彈簧k3i對旋轉FGMs連接殼的行波模態(tài)影響類似。當彈簧剛度取10-4~106時,結構的行波模態(tài)幾乎不受影響。而當剛度在106~1012時,前后行波模態(tài)頻率呈線性上升的趨勢,最后收斂在1013處。行波模態(tài)在k1i取10-4~108時,k4i取10-4~104時,軸向彈簧k1i和扭轉彈簧k4i的共同作用幾乎不受影響。在k1i取108~1012以及k4i取104~108時,結構的行波模態(tài)頻率快速上升,且最終分別收斂于1013、109,因此軸向彈簧的收斂速度要遠低于扭轉彈簧,其對結構行波模態(tài)的影響遠大于扭轉彈簧??梢钥闯?,在彈性邊界時,結構行波模態(tài)變化明顯,因此采用彈簧模擬邊界條件是必要的。

    圖7 兩組彈簧剛度(ki2,ki3)對旋轉FGMs錐-柱連接殼行波模態(tài)的影響Fig.7 Effects of two kinds of boundary spring stiffness coefficients(ki2,ki3)on the traveling wave mode for a rotating FGMs joined conical-cylindrical shell

    圖8 兩組彈簧剛度(ki1,ki4)對旋轉FGMs錐-柱連接殼行波模態(tài)的影響Fig.8 Effects of two kinds of boundary spring stiffness coefficients(ki1,ki4)on the traveling wave mode for a rotating FGMs joined conical-cylindrical shell

    由圖6可知,轉速對殼體的行波模態(tài)影響明顯,因此圖9探究了轉速和單組彈簧剛度共同影響下,旋轉FGMs錐-柱連接殼的行波模態(tài)變化。除所討論的邊界彈簧剛度外,其他邊界彈簧剛度值取1013。

    圖9 單組彈簧剛度和旋轉速度對旋轉FGMs錐-柱連接殼行波模態(tài)的影響Fig.9 Effects of single boundary spring stiffness coefficients and rotational speed on the traveling wave mode for a rotating FGMs joined conical-cylindrical shell

    分析發(fā)現(xiàn),軸向彈簧(k1i)對連接殼的行波模態(tài)影響最大,環(huán)向彈簧(k2i)對連接殼的行波模態(tài)作用最小,而徑向彈簧(k3i)和扭轉彈簧(k4i)所產生的效果類似。相比后行波模態(tài),前行波模態(tài)對邊界彈簧剛度更為敏感,在彈簧剛度取106~1012時,模態(tài)頻率呈現(xiàn)顯著上升的趨勢。隨著轉速的提高,后行波模態(tài)一直呈現(xiàn)上升態(tài)勢。對于前行波模態(tài),頻率總是先減小,而當轉速高于0.004 時呈現(xiàn)增大趨勢。在改變環(huán)向、徑向以及扭轉彈簧剛度時,旋轉FGMs錐-柱連接殼行波模態(tài)頻率主要受轉速的影響。而在軸向彈簧與轉速共同作用時,軸向彈簧對連接殼的行波模態(tài)起主要作用。

    5 結 論

    本文采用彈簧模擬旋轉FGMs 錐-柱連接殼的任意邊界條件和連接條件,推導了旋轉連接殼在考慮科氏力和離心力作用下任意邊界條件的理論模型,利用Chebyshev多項式和Rayleigh-Ritz 法求解了結構的行波模態(tài)頻率,分析了各種參數(shù)對行波模態(tài)的影響。主要結論如下:

    (1)采用人工彈簧技術,有效實現(xiàn)了實際工程中的彈性邊界,在計算求解過程中,相比于傳統(tǒng)能量法,減少了大量重復的計算,提高了計算效率。

    (2)隨著陶瓷體積分數(shù)的增大,旋轉FGMs錐-柱連接殼前后行波模態(tài)頻率呈現(xiàn)下降趨勢,當N>5后,行波模態(tài)頻率的變化趨勢逐漸平緩;工程實際中,可以通過合理設計N以突出梯度材料組分中某一材料的優(yōu)越性能。

    (3)由F-F 邊界過渡到S-S 邊界,環(huán)向彈簧剛度(ki2)和徑向彈簧剛度(ki3)對結構行波模態(tài)的影響趨勢基本一致,且行波模態(tài)收斂所對應的彈簧剛度都為1013;由S-S邊界過渡到C-C邊界,軸向彈簧剛度(ki1)相對扭轉彈簧剛度(ki4)對結構行波模態(tài)的影響更明顯,彈簧剛度分別取ki1>1013、ki4>109后,行波模態(tài)頻率趨于收斂。

    (4)針對短薄連接殼結構,軸向彈簧剛度對其行波模態(tài)影響最為顯著,而軸向彈簧相較于轉速對行波模態(tài)又有更顯著的影響;在臨界轉速之前,前行波模態(tài)頻率呈下降趨勢,之后上升;當結構在臨界轉速工況下,前行波模態(tài)頻率最小,結構更容易發(fā)生失穩(wěn),實際中應該避免臨界轉速的工況。

    猜你喜歡
    環(huán)向行波邊界條件
    金屬管材力學性能多角度環(huán)向拉伸測試方法研究
    自承式鋼管跨越結構鞍式支承處管壁環(huán)向彎曲應力分析
    特種結構(2022年6期)2023-01-12 09:17:02
    一類非局部擴散的SIR模型的行波解
    一類帶有Stieltjes積分邊界條件的分數(shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動邊值問題的漸近解
    城市供水管網中鋼筋混凝土岔管受力分析
    Joseph-Egri方程行波解的分岔
    英國MACAW公司依據(jù)CEPA 2015提出管道環(huán)向應力腐蝕開裂預防處理改進方法
    Kolmogorov-Petrovskii-Piskunov方程和Zhiber-Shabat方程的行波解
    帶Robin邊界條件的2維隨機Ginzburg-Landau方程的吸引子
    国产精品国产av在线观看| 视频区欧美日本亚洲| 国产有黄有色有爽视频| 老司机午夜十八禁免费视频| av国产久精品久网站免费入址| 亚洲av综合色区一区| 另类亚洲欧美激情| 婷婷色麻豆天堂久久| e午夜精品久久久久久久| 黄色a级毛片大全视频| 久热爱精品视频在线9| 男女午夜视频在线观看| 精品久久久久久久毛片微露脸 | 人人妻人人澡人人爽人人夜夜| 天天躁夜夜躁狠狠久久av| 国产一区亚洲一区在线观看| 老司机亚洲免费影院| 女人高潮潮喷娇喘18禁视频| 99精国产麻豆久久婷婷| 亚洲五月婷婷丁香| 久久午夜综合久久蜜桃| 国产亚洲午夜精品一区二区久久| 精品久久久精品久久久| 午夜老司机福利片| 精品国产乱码久久久久久小说| 欧美日韩一级在线毛片| 视频区图区小说| 免费久久久久久久精品成人欧美视频| 欧美少妇被猛烈插入视频| 一级毛片电影观看| 大香蕉久久成人网| 亚洲九九香蕉| 天天添夜夜摸| 18禁国产床啪视频网站| 黄片小视频在线播放| 亚洲精品美女久久久久99蜜臀 | 久久九九热精品免费| 免费人妻精品一区二区三区视频| 亚洲国产精品一区二区三区在线| 美女高潮到喷水免费观看| 国产三级黄色录像| 中文字幕制服av| 欧美亚洲日本最大视频资源| 免费不卡黄色视频| 黄频高清免费视频| 老司机午夜十八禁免费视频| 最新在线观看一区二区三区 | 久久久久久人人人人人| 少妇人妻久久综合中文| 久久性视频一级片| 性色av乱码一区二区三区2| 免费观看人在逋| 国产福利在线免费观看视频| 18在线观看网站| 精品福利永久在线观看| 欧美xxⅹ黑人| 亚洲中文av在线| 免费高清在线观看日韩| 精品久久久久久久毛片微露脸 | 日本91视频免费播放| 精品一区二区三卡| 少妇被粗大的猛进出69影院| 亚洲成人国产一区在线观看 | 乱人伦中国视频| 大香蕉久久网| 午夜老司机福利片| 国产老妇伦熟女老妇高清| av在线老鸭窝| 亚洲欧美精品自产自拍| 免费av中文字幕在线| 日韩 亚洲 欧美在线| 午夜福利一区二区在线看| 天天躁夜夜躁狠狠久久av| 一级毛片电影观看| 热re99久久精品国产66热6| 人妻一区二区av| 中文欧美无线码| av在线播放精品| 国产一级毛片在线| 国产高清不卡午夜福利| av又黄又爽大尺度在线免费看| 一本久久精品| 久热这里只有精品99| 美女大奶头黄色视频| 免费av中文字幕在线| 欧美成人午夜精品| 午夜福利,免费看| 国产色视频综合| av国产精品久久久久影院| 悠悠久久av| www.av在线官网国产| 又大又黄又爽视频免费| 2021少妇久久久久久久久久久| 在线观看一区二区三区激情| 亚洲av国产av综合av卡| av一本久久久久| 又紧又爽又黄一区二区| 亚洲成人国产一区在线观看 | 亚洲伊人色综图| 黑丝袜美女国产一区| 日本wwww免费看| 精品国产国语对白av| av在线老鸭窝| 国产精品偷伦视频观看了| 高清av免费在线| av电影中文网址| 精品福利观看| 久久精品国产亚洲av涩爱| 精品久久久久久久毛片微露脸 | 亚洲国产中文字幕在线视频| 狂野欧美激情性bbbbbb| 国产欧美亚洲国产| 中文精品一卡2卡3卡4更新| 中文乱码字字幕精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲七黄色美女视频| 午夜激情久久久久久久| 亚洲精品在线美女| 黄色视频不卡| 亚洲国产欧美在线一区| 国产亚洲一区二区精品| 国产欧美日韩精品亚洲av| 亚洲人成电影免费在线| 亚洲午夜精品一区,二区,三区| 免费看不卡的av| 亚洲美女黄色视频免费看| 久久久久久久精品精品| 少妇精品久久久久久久| 欧美亚洲 丝袜 人妻 在线| 777米奇影视久久| 亚洲欧美一区二区三区国产| 国产精品九九99| 亚洲伊人久久精品综合| 亚洲欧美精品综合一区二区三区| 黄色一级大片看看| 成年女人毛片免费观看观看9 | 免费观看av网站的网址| 欧美久久黑人一区二区| 精品人妻1区二区| 亚洲欧美一区二区三区黑人| 久久午夜综合久久蜜桃| 丝袜喷水一区| 国产黄频视频在线观看| 一本综合久久免费| 国产精品国产av在线观看| 69精品国产乱码久久久| 国产片特级美女逼逼视频| 亚洲,欧美精品.| 麻豆国产av国片精品| 欧美乱码精品一区二区三区| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 亚洲五月色婷婷综合| www.av在线官网国产| 国产成人欧美在线观看 | 日本午夜av视频| 国语对白做爰xxxⅹ性视频网站| 精品人妻1区二区| 欧美精品人与动牲交sv欧美| 国产免费视频播放在线视频| 久热这里只有精品99| 麻豆国产av国片精品| av电影中文网址| 另类亚洲欧美激情| 日本wwww免费看| 曰老女人黄片| 看十八女毛片水多多多| 欧美大码av| 免费看av在线观看网站| 国产免费福利视频在线观看| 日本av手机在线免费观看| 亚洲伊人色综图| 一级片免费观看大全| av在线老鸭窝| 免费人妻精品一区二区三区视频| 热99国产精品久久久久久7| 久久久欧美国产精品| av福利片在线| 国产视频一区二区在线看| 国产欧美亚洲国产| 久久久精品94久久精品| 亚洲少妇的诱惑av| av网站免费在线观看视频| 亚洲av成人精品一二三区| 亚洲中文日韩欧美视频| 9191精品国产免费久久| 黄色片一级片一级黄色片| 国产黄频视频在线观看| 国产成人啪精品午夜网站| 国语对白做爰xxxⅹ性视频网站| 欧美 日韩 精品 国产| a级片在线免费高清观看视频| 久久鲁丝午夜福利片| 永久免费av网站大全| 久久99一区二区三区| 久久九九热精品免费| 视频区图区小说| 一本—道久久a久久精品蜜桃钙片| 免费av中文字幕在线| 欧美中文综合在线视频| 精品久久久精品久久久| 黄色 视频免费看| h视频一区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲精品久久久久久婷婷小说| 国产成人av教育| 亚洲精品国产色婷婷电影| 色播在线永久视频| 欧美亚洲日本最大视频资源| 亚洲第一青青草原| 亚洲国产精品一区二区三区在线| 日日爽夜夜爽网站| 久久久欧美国产精品| 黄色片一级片一级黄色片| 亚洲欧美一区二区三区国产| 黄频高清免费视频| 大片免费播放器 马上看| 另类亚洲欧美激情| 一本—道久久a久久精品蜜桃钙片| 性色av一级| 老司机影院毛片| 一级毛片电影观看| 国产又爽黄色视频| 99久久精品国产亚洲精品| 一二三四在线观看免费中文在| 91麻豆精品激情在线观看国产 | www.av在线官网国产| 熟女少妇亚洲综合色aaa.| 丝袜在线中文字幕| 亚洲精品自拍成人| www日本在线高清视频| 精品少妇黑人巨大在线播放| 多毛熟女@视频| 午夜免费成人在线视频| 丝袜人妻中文字幕| www日本在线高清视频| 岛国毛片在线播放| 男女边吃奶边做爰视频| 免费少妇av软件| 亚洲人成电影免费在线| 性色av一级| 男女免费视频国产| netflix在线观看网站| 国产一卡二卡三卡精品| 国产亚洲欧美精品永久| 十八禁高潮呻吟视频| 久久九九热精品免费| 五月天丁香电影| 高清视频免费观看一区二区| 在线av久久热| 欧美日韩黄片免| 男女免费视频国产| 夜夜骑夜夜射夜夜干| 99re6热这里在线精品视频| 婷婷丁香在线五月| 乱人伦中国视频| 中文欧美无线码| 亚洲午夜精品一区,二区,三区| 久久精品人人爽人人爽视色| 50天的宝宝边吃奶边哭怎么回事| 两个人免费观看高清视频| 桃花免费在线播放| 少妇粗大呻吟视频| 一区二区三区精品91| 我要看黄色一级片免费的| 久久久久精品人妻al黑| 国产精品一区二区免费欧美 | 满18在线观看网站| 丰满少妇做爰视频| 国产免费现黄频在线看| 亚洲精品久久午夜乱码| 热99久久久久精品小说推荐| 成人三级做爰电影| 国产在线观看jvid| 一本色道久久久久久精品综合| 老汉色∧v一级毛片| 精品福利永久在线观看| 五月天丁香电影| 久久国产精品大桥未久av| 午夜免费鲁丝| 午夜福利,免费看| 日本猛色少妇xxxxx猛交久久| 精品亚洲成a人片在线观看| 日本午夜av视频| 国产亚洲av高清不卡| 性色av乱码一区二区三区2| 脱女人内裤的视频| 97在线人人人人妻| 99国产精品一区二区蜜桃av | 久久久国产欧美日韩av| 亚洲欧美成人综合另类久久久| 精品国产超薄肉色丝袜足j| cao死你这个sao货| 欧美性长视频在线观看| 少妇裸体淫交视频免费看高清 | 老司机影院成人| videosex国产| 高清不卡的av网站| 亚洲精品美女久久久久99蜜臀 | 人人妻人人澡人人爽人人夜夜| 人人澡人人妻人| 国产主播在线观看一区二区 | 免费在线观看完整版高清| 久久久久国产一级毛片高清牌| 成年av动漫网址| 超色免费av| 曰老女人黄片| 男的添女的下面高潮视频| 精品久久久精品久久久| 人人妻人人澡人人爽人人夜夜| 久久99热这里只频精品6学生| 国产欧美日韩一区二区三 | 脱女人内裤的视频| 国产欧美日韩精品亚洲av| 高清黄色对白视频在线免费看| 菩萨蛮人人尽说江南好唐韦庄| 老汉色∧v一级毛片| 亚洲人成电影观看| 建设人人有责人人尽责人人享有的| 国产成人精品无人区| 久久精品国产综合久久久| www.av在线官网国产| 黄色a级毛片大全视频| 国产欧美日韩一区二区三 | 午夜免费鲁丝| 一二三四社区在线视频社区8| 制服人妻中文乱码| 蜜桃在线观看..| 一级片'在线观看视频| 精品人妻一区二区三区麻豆| 国产一区二区在线观看av| 久久久久久人人人人人| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 国产精品一国产av| 又紧又爽又黄一区二区| 99国产精品一区二区蜜桃av | 超色免费av| 黄色视频在线播放观看不卡| 美女大奶头黄色视频| 不卡av一区二区三区| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| 大香蕉久久网| 脱女人内裤的视频| 国产一区二区三区综合在线观看| 日本av手机在线免费观看| 亚洲男人天堂网一区| 久久午夜综合久久蜜桃| 亚洲人成网站在线观看播放| 国产片内射在线| 日韩精品免费视频一区二区三区| 一边亲一边摸免费视频| 两个人免费观看高清视频| 日本猛色少妇xxxxx猛交久久| 精品国产一区二区三区四区第35| 日本猛色少妇xxxxx猛交久久| 男人添女人高潮全过程视频| √禁漫天堂资源中文www| 国产淫语在线视频| www.精华液| 国产精品香港三级国产av潘金莲 | 亚洲精品美女久久久久99蜜臀 | 黄频高清免费视频| 在线 av 中文字幕| 天天操日日干夜夜撸| 日本欧美国产在线视频| 满18在线观看网站| 日韩 亚洲 欧美在线| 久久天躁狠狠躁夜夜2o2o | 国产成人av教育| 久久人人97超碰香蕉20202| 日本午夜av视频| 少妇粗大呻吟视频| 一区二区av电影网| 免费黄频网站在线观看国产| 啦啦啦在线观看免费高清www| 亚洲国产毛片av蜜桃av| 国产成人欧美在线观看 | 国产精品三级大全| 国产91精品成人一区二区三区 | 久久天躁狠狠躁夜夜2o2o | 中文欧美无线码| 美女福利国产在线| 久久久欧美国产精品| 国产视频首页在线观看| 亚洲国产精品999| 欧美日韩视频高清一区二区三区二| 亚洲第一青青草原| 亚洲 欧美一区二区三区| 国产一区二区激情短视频 | 超碰97精品在线观看| 只有这里有精品99| 天堂俺去俺来也www色官网| 最新的欧美精品一区二区| 亚洲一区中文字幕在线| 日韩一本色道免费dvd| 一级黄片播放器| 丝袜人妻中文字幕| 亚洲国产欧美日韩在线播放| 又大又黄又爽视频免费| 亚洲成人免费av在线播放| 成人国产av品久久久| 最近手机中文字幕大全| 国产免费福利视频在线观看| 丝袜脚勾引网站| 纵有疾风起免费观看全集完整版| 一级a爱视频在线免费观看| 亚洲七黄色美女视频| 又黄又粗又硬又大视频| 99国产精品一区二区蜜桃av | 午夜av观看不卡| 亚洲欧美日韩另类电影网站| 日本色播在线视频| 亚洲精品国产区一区二| 操出白浆在线播放| 午夜福利视频在线观看免费| 久久久国产一区二区| 在线观看免费视频网站a站| 久久久久精品国产欧美久久久 | 91精品三级在线观看| 我的亚洲天堂| 女性生殖器流出的白浆| av电影中文网址| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| 观看av在线不卡| 日韩制服丝袜自拍偷拍| 国产成人欧美在线观看 | 99精品久久久久人妻精品| 三上悠亚av全集在线观看| 亚洲人成电影免费在线| 午夜福利视频在线观看免费| 欧美少妇被猛烈插入视频| 国产成人一区二区在线| 日日爽夜夜爽网站| 国产亚洲精品久久久久5区| 9热在线视频观看99| 一级片免费观看大全| 亚洲中文字幕日韩| 亚洲欧美成人综合另类久久久| 丁香六月天网| 久久久精品94久久精品| 日韩,欧美,国产一区二区三区| 亚洲国产欧美网| 国产男女超爽视频在线观看| 欧美精品啪啪一区二区三区 | 国产麻豆69| 精品人妻1区二区| 中文精品一卡2卡3卡4更新| 国产精品久久久人人做人人爽| 色视频在线一区二区三区| 午夜福利免费观看在线| 国产片特级美女逼逼视频| 国产成人影院久久av| 老汉色∧v一级毛片| 少妇人妻久久综合中文| 国产主播在线观看一区二区 | 久久人妻福利社区极品人妻图片 | 如日韩欧美国产精品一区二区三区| 91字幕亚洲| 日韩中文字幕视频在线看片| kizo精华| 50天的宝宝边吃奶边哭怎么回事| 最近中文字幕2019免费版| 欧美日韩av久久| 欧美精品一区二区免费开放| 精品人妻一区二区三区麻豆| svipshipincom国产片| 午夜福利,免费看| 热re99久久国产66热| 久久99一区二区三区| 日韩中文字幕欧美一区二区 | 在线观看免费视频网站a站| 国产精品一区二区免费欧美 | e午夜精品久久久久久久| 成年动漫av网址| 人人妻人人添人人爽欧美一区卜| 人体艺术视频欧美日本| 国产精品免费大片| 欧美老熟妇乱子伦牲交| 中文字幕av电影在线播放| 丝袜美腿诱惑在线| 精品少妇一区二区三区视频日本电影| 久久精品亚洲av国产电影网| 国产视频一区二区在线看| 十八禁人妻一区二区| 国产精品国产三级国产专区5o| 亚洲av日韩精品久久久久久密 | 爱豆传媒免费全集在线观看| 巨乳人妻的诱惑在线观看| 2018国产大陆天天弄谢| 秋霞在线观看毛片| 电影成人av| 久久精品久久久久久久性| 黄色视频不卡| 丁香六月欧美| 久久精品人人爽人人爽视色| 高潮久久久久久久久久久不卡| 91精品三级在线观看| 亚洲精品av麻豆狂野| 久久久欧美国产精品| 精品一区二区三区四区五区乱码 | 国产精品99久久99久久久不卡| 久久99精品国语久久久| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久毛片微露脸 | 首页视频小说图片口味搜索 | 18在线观看网站| 老司机亚洲免费影院| 欧美日本中文国产一区发布| 青春草亚洲视频在线观看| 亚洲第一青青草原| 国产片内射在线| 宅男免费午夜| 伦理电影免费视频| 欧美激情 高清一区二区三区| 国产视频首页在线观看| 宅男免费午夜| 黄片播放在线免费| 久久久久精品人妻al黑| 国产av国产精品国产| 亚洲精品美女久久av网站| 久久国产精品人妻蜜桃| 男女床上黄色一级片免费看| 啦啦啦中文免费视频观看日本| 啦啦啦 在线观看视频| 国产成人免费无遮挡视频| 午夜福利免费观看在线| 日日爽夜夜爽网站| 天天添夜夜摸| 欧美性长视频在线观看| 国产精品久久久久久人妻精品电影 | 岛国毛片在线播放| 下体分泌物呈黄色| 99精品久久久久人妻精品| 久久狼人影院| 校园人妻丝袜中文字幕| 狠狠婷婷综合久久久久久88av| 国产亚洲av高清不卡| 国产成人91sexporn| 亚洲成人国产一区在线观看 | 麻豆av在线久日| 亚洲av在线观看美女高潮| 久久ye,这里只有精品| 波野结衣二区三区在线| 久久久久国产精品人妻一区二区| 视频区图区小说| 亚洲人成网站在线观看播放| 亚洲欧美一区二区三区久久| 91国产中文字幕| 男女午夜视频在线观看| 欧美日韩视频高清一区二区三区二| 欧美xxⅹ黑人| 免费在线观看影片大全网站 | 三上悠亚av全集在线观看| 爱豆传媒免费全集在线观看| 捣出白浆h1v1| 97人妻天天添夜夜摸| 亚洲国产av影院在线观看| 国产91精品成人一区二区三区 | 国产99久久九九免费精品| 纵有疾风起免费观看全集完整版| 好男人电影高清在线观看| 久久久久精品国产欧美久久久 | 欧美亚洲 丝袜 人妻 在线| 精品亚洲成国产av| 国产人伦9x9x在线观看| 香蕉丝袜av| 免费在线观看影片大全网站 | 亚洲精品一区蜜桃| 国产人伦9x9x在线观看| 国产av精品麻豆| 日韩中文字幕欧美一区二区 | 国产精品二区激情视频| 2021少妇久久久久久久久久久| 只有这里有精品99| 亚洲少妇的诱惑av| 高潮久久久久久久久久久不卡| 欧美在线一区亚洲| 人人妻,人人澡人人爽秒播 | 在线观看免费视频网站a站| 中文字幕高清在线视频| 一区二区三区精品91| 美女午夜性视频免费| 一区二区三区四区激情视频| 久久久久久久久久久久大奶| 极品少妇高潮喷水抽搐| 国产片内射在线| 人体艺术视频欧美日本| 黄色 视频免费看| 国产精品 欧美亚洲| 亚洲成色77777| 在线观看www视频免费| 一级片'在线观看视频| 亚洲av电影在线进入| 久久久久久免费高清国产稀缺| 美女扒开内裤让男人捅视频| 黄色一级大片看看| 麻豆乱淫一区二区| 啦啦啦在线观看免费高清www| 久久久国产精品麻豆| 欧美日韩av久久| 女人精品久久久久毛片| 看免费成人av毛片| 亚洲精品成人av观看孕妇| 首页视频小说图片口味搜索 | 热99国产精品久久久久久7| 黄色片一级片一级黄色片| 亚洲国产精品999| 国产精品久久久久久精品电影小说| 极品人妻少妇av视频| 捣出白浆h1v1| 99热全是精品| av网站在线播放免费|