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

    熱環(huán)境下功能梯度環(huán)板的譜幾何法自由振動(dòng)解

    2022-05-30 10:55:16石先杰
    振動(dòng)與沖擊 2022年10期
    關(guān)鍵詞:環(huán)板固有頻率梯度

    石先杰, 左 朋,2

    (1. 中國工程物理研究院 總體工程研究所,四川 綿陽 621999;2. 中國科學(xué)技術(shù)大學(xué) 近代力學(xué)系,合肥 230026)

    功能梯度材料[1-3]作為一種新興材料,它克服了傳統(tǒng)材料在過熱環(huán)境下機(jī)械性能下降的缺陷,具有極強(qiáng)的耐熱性,而由其所制成的環(huán)板結(jié)構(gòu)有著質(zhì)量輕、結(jié)構(gòu)剛度大、力學(xué)性能好、承載能力強(qiáng)等優(yōu)點(diǎn),在航空航天工程、飛機(jī)工程、核反應(yīng)堆等領(lǐng)域的發(fā)動(dòng)機(jī)端蓋上有著廣泛應(yīng)用。因此,這類結(jié)構(gòu)在實(shí)際工程應(yīng)用中通常會(huì)受到復(fù)雜工作環(huán)境的影響而產(chǎn)生較大的振動(dòng)甚至引起疲勞損傷。因此,準(zhǔn)確分析預(yù)測(cè)功能梯度環(huán)板結(jié)構(gòu)振動(dòng)特性一直是國內(nèi)外專家學(xué)者的研究熱點(diǎn)。為此,國內(nèi)外學(xué)者已開展了一系列數(shù)值求解方法研究,如Rayleigh-Ritz法[4-6]、微分求積法[7-10]、擬格林函數(shù)法[11]和有限環(huán)形棱柱法[12]等。

    然而,上述研究均是基于常溫環(huán)境,在高溫環(huán)境下,結(jié)構(gòu)的材料屬性會(huì)隨著溫度的變化而改變,從而導(dǎo)致其自由振動(dòng)的固有頻率相比于常溫下變得明顯不同。對(duì)于熱環(huán)境下功能梯度結(jié)構(gòu)的振動(dòng)特性,國內(nèi)外學(xué)者進(jìn)行了一些研究分析。假設(shè)材料特性依賴于溫度,并在板的厚度方向上逐漸變化,Haddadpour等[13]采用伽遼金法求解運(yùn)動(dòng)方程,對(duì)簡支功能梯度圓柱殼在熱環(huán)境中的自由振動(dòng)特性進(jìn)行了分析。Mirtalaie[14]根據(jù)微分求積法,求解了熱環(huán)境下功能梯度薄環(huán)形扇形板自由振動(dòng)頻率。在三維彈性理論的框架下,Yang等[15]分析了功能梯度環(huán)板的熱振響應(yīng)。在一階剪切變形理論的框架內(nèi)應(yīng)用微分求積法,Malekzadeh等[16-17]研究了旋轉(zhuǎn)功能梯度圓柱殼和帶有彈性支撐的功能梯度環(huán)板在熱場(chǎng)中的自由振動(dòng)特性。用傅里葉級(jí)數(shù)結(jié)合輔助函數(shù)來描述板的位移容許函數(shù),Zhou等[18]研究了彈性地基上具有溫度梯度的超音速多孔功能梯度材料板的振動(dòng)和顫振特性。Li等[19]采用特征正交多項(xiàng)式詳細(xì)研究了功能梯度多孔階梯圓柱殼在熱場(chǎng)中的自由振動(dòng)、穩(wěn)態(tài)和瞬態(tài)響應(yīng)特性。蒲育等[20-21]根據(jù)二維彈性理論,采用微分求積法分析了材料屬性徑向變化的功能梯度薄環(huán)板在熱環(huán)境下的面內(nèi)自由振動(dòng)特性。同樣是考慮材料屬性沿徑向變化功能梯度薄環(huán)板,呂朋等[22]采用改進(jìn)傅里葉級(jí)數(shù)建立了其在熱場(chǎng)中受彈性約束的面內(nèi)振動(dòng)分析模型。而對(duì)于沿厚度方向變化的功能梯度環(huán)板,石鵬等[23]應(yīng)用三維Chebyshev-Ritz法對(duì)其在熱環(huán)境下的自由振動(dòng)特性進(jìn)行了求解。

    綜上所述,目前關(guān)于功能梯度環(huán)板在熱場(chǎng)中的自由振動(dòng)分析雖然有一定的研究成果,但大部分研究工作集中在材料屬性沿徑向變化的薄環(huán)板面內(nèi)振動(dòng)問題。而對(duì)于厚環(huán)板問題,往往是采用三維彈性理論求解,存在著計(jì)算效率低等問題。因此,針對(duì)工程中更為常見的中厚環(huán)板,考慮材料特性沿厚度方向的梯度變化,本文以一階剪切變形理論為基礎(chǔ),通過沿邊界設(shè)置具有可變剛度的邊界約束彈簧的方式來考慮任意邊界約束。同時(shí),采用譜幾何法建立功能梯度環(huán)板熱振分析模型。在譜幾何法求解框架下,板的位移函數(shù)總是表示為標(biāo)準(zhǔn)的傅里葉余弦級(jí)數(shù),并輔以四個(gè)正弦函數(shù)項(xiàng)來消除邊界處位移函數(shù)的導(dǎo)數(shù)不連續(xù)性問題。這種級(jí)數(shù)展開方法理論上能夠表征任意有限長結(jié)構(gòu)的位移變量,且函數(shù)形式統(tǒng)一,易于參數(shù)化分析。在數(shù)值算例中,通過與文獻(xiàn)解、有限元解的對(duì)比研究,驗(yàn)證了該模型的正確性,進(jìn)而分析了邊界條件、環(huán)境溫度、梯度指數(shù)等對(duì)功能梯度環(huán)板固有頻率的影響。

    1 理論方法

    1.1 模型描述

    考慮如圖1所示的功能梯度環(huán)板,圓柱坐標(biāo)系(r,θ,z)位于結(jié)構(gòu)的幾何中面上。圖1中:R0和R1分別為環(huán)板沿徑向方向r的內(nèi)徑和外徑;θ為圓周方向;h為環(huán)板的厚度。文中研究的功能梯度材料由陶瓷和金屬材料兩種成分組成,其有效材料參數(shù)P(包括彈性模量E、密度ρ、泊松比ν和熱膨脹系數(shù)α等)沿環(huán)板厚度方向呈梯度變化,可以被描述為[24]

    (1)

    式中:Pt和Pu分別為環(huán)板上、下表面的材料參數(shù);-h/2≤z≤h/2;p為功能梯度指數(shù)??紤]熱環(huán)境的影響,功能梯度材料的有效材料參數(shù)Pξ(ξ=t,u)在溫度值T的變化影響下有著如式(2)所示的性質(zhì)

    Pξ(z,T)=P0ξ(P-1ξT-1+1+P1ξT+P2ξT2+P3ξT3)

    (2)

    式中,P-1ξ,P0ξ,P1ξ,P2ξ和P3ξ為溫度相關(guān)系數(shù)。文中研究的是均勻溫度分布的情形[25],T0為參考溫度,其溫度變化為ΔT=T-T0。

    圖1 功能梯度環(huán)板的幾何形狀以及坐標(biāo)系Fig.1 Geometric shape and coordinate system of functionally graded annular plate

    1.2 能量方程

    根據(jù)一階剪切變形理論,環(huán)板上任一點(diǎn)沿r,θ和z方向的位移場(chǎng)分量u,v,w的表達(dá)式可以寫成

    u(r,θ,z,t)=u0(r,θ,t)+zφr(r,θ,t),
    v(r,θ,z,t)=v0(r,θ,t)+zφθ(r,θ,t),
    w(r,θ,z,t)=w0(r,θ,t)

    (3)

    式中:u0,v0,w0分別為環(huán)板在幾何中面上沿r,θ和z方向的位移分量;φr和φθ分別為在rz和θz平面上的轉(zhuǎn)角位移分量。

    根據(jù)中厚板理論,環(huán)板結(jié)構(gòu)中的應(yīng)變和位移之間存在線性關(guān)系式

    (4)

    (5)

    ε=[εr,εθ,γrθ,γrz,γθz]T,
    σ=[σr,σθ,τrθ,τrz,τθz]T,
    σ=Q(5×5)ε

    (6)

    式中,Q(5×5)為與溫度值T和厚度坐標(biāo)z有關(guān)的彈性常數(shù)所組成的矩陣(參考2011年Malekzadeh等的研究),其元素的具體值為

    (7)

    對(duì)應(yīng)力分量σ在z方向積分,引入剪切修正系數(shù)κ=5/6, 得到環(huán)板的本構(gòu)方程表達(dá)式

    N=Aε0+Bχ,M=Bε0+Dχ,Qκ=Aκγ

    (8)

    式中:N=[Nr,Nθ,Nrθ]T和Q=[Qr,Qθ]T分別為環(huán)板的合力和橫向剪應(yīng)力分量;M=[Mr,Mθ,Mrθ]T為合彎矩分量;A(Aκ),B和D分別為以結(jié)構(gòu)的拉伸剛度、拉伸-彎曲耦合剛度和彎曲剛度為元素組成的矩陣,其具體表達(dá)式可見2011年Malekzadeh等的研究。

    根據(jù)建立的本構(gòu)方程表達(dá)式,結(jié)合2011年Malekzadeh等給出的熱應(yīng)力描述,得到環(huán)板的應(yīng)變能表達(dá)式

    (9)

    代入位移場(chǎng)分量可以得到環(huán)板的動(dòng)能表達(dá)式

    (10)

    文中通過在環(huán)板的內(nèi)部和外部放置邊界彈簧實(shí)現(xiàn)邊界條件的模擬。邊界約束彈簧所儲(chǔ)存的能量可以寫成

    (11)

    式中,上標(biāo)0和1分別為設(shè)置在環(huán)板內(nèi)部和外部的邊界彈簧剛度值。

    1.3 位移容許函數(shù)以及方程的求解

    文中分別采用了譜幾何法和傅里葉正余弦函數(shù)來表示環(huán)板沿徑向方向和周向方向的位移容許函數(shù)

    (12)

    如上所述,熱環(huán)境下功能梯度環(huán)板的拉格朗日能量泛函可以寫成

    (13)

    根據(jù)1.2節(jié)的能量方程表達(dá)式,結(jié)合相應(yīng)的位移容許函數(shù)進(jìn)行求解,并采用Rayleigh-Ritz法對(duì)未知展開系數(shù)求偏導(dǎo),可以得到熱環(huán)境下功能梯度環(huán)板自由振動(dòng)特征方程

    (K-ω2M)H=0

    (14)

    式中:K和M分別為環(huán)板的剛度矩陣和質(zhì)量矩陣;H為未知系數(shù)向量,其中溫度變化會(huì)影響K。

    2 數(shù)值算例與分析

    根據(jù)上述分析獲得的自由振動(dòng)特征方程,本章將分析熱環(huán)境下功能梯度環(huán)板自由振動(dòng)特性。根據(jù)航空航天飛行器的工程實(shí)際應(yīng)用,在接下來的數(shù)值算例中選擇了幾何尺寸為R0=0.5 m,R1=2.0 m,h=0.10 m的環(huán)板結(jié)構(gòu)作為研究對(duì)象。通過設(shè)置邊界約束彈簧的剛度值,文中研究考慮了包括經(jīng)典和彈性邊界在內(nèi)的10種任意邊界條件,并且采用了不同字母來表示,其中: C,SS,SD和F分別為固支、簡支、剪切-膜片邊界和自由邊界; E1,E2,E3,E4,E5和E6為彈性邊界。它們的彈簧剛度值按Su等的研究設(shè)置。此外,位移容許函數(shù)在徑向和周向方向的級(jí)數(shù)在理論上可以展開為無窮項(xiàng)。但是考慮到實(shí)際數(shù)值計(jì)算需要,文中對(duì)其截?cái)嘀劣邢揄?xiàng)。參考文獻(xiàn)[26]對(duì)譜幾何法收斂性分析結(jié)果,文中研究將徑向和周向的級(jí)數(shù)均截?cái)酁镸=N=14。

    2.1 數(shù)值驗(yàn)證

    不同邊界條件和梯度指數(shù)p下功能梯度環(huán)板在不考慮熱環(huán)境影響情況下的第1階固有頻率,如表1所示。算例包括了兩種板的厚度h=0.05 m和h=0.10 m, 其余幾何尺寸與默認(rèn)參數(shù)保持一致。環(huán)板的上、下表面材料分別為鋁和氧化鋯,材料參數(shù)為:Et=168 GPa,νt=0.3,ρt=5 700 kg/m3,Eu=70 GPa,νu=0.3,ρu=2 707 kg/m3。表1中將譜幾何法模型計(jì)算結(jié)果與Su等的研究結(jié)果進(jìn)行了對(duì)比。從表1可以看出,文中結(jié)果與Su等的解非常接近,存在的偏差主要來源于徑向所采用的位移函數(shù)差異。這表明文中模型邊界約束彈簧剛度以及截?cái)鄶?shù)的設(shè)置是合理的,可以有效求解具有不同梯度指數(shù)p的功能梯度環(huán)板的自由振動(dòng)特性。

    在開展熱環(huán)境下功能梯度環(huán)板結(jié)構(gòu)自由振動(dòng)分析時(shí),由于缺乏相關(guān)文獻(xiàn)數(shù)據(jù)作為參考。因此,文中采用有限元法(finite element method,F(xiàn)EM)計(jì)算結(jié)果作為對(duì)比數(shù)據(jù)。環(huán)板上、下表面的材料分別為氮化硅(Si3N4)和304不銹鋼(SUS304),它們?cè)跓釄?chǎng)中的材料參數(shù)如表2所示。環(huán)板的幾何尺寸與默認(rèn)參數(shù)保持一致,參考溫度取為T0=300 K。同樣是考慮了C-C,SS-SS和C-SS三種經(jīng)典邊界條件。不同溫度變化ΔT下環(huán)板前5階固有頻率對(duì)比情況(梯度指數(shù)p=1.0),如表3所示。不同梯度指數(shù)p下的前9階固有頻率對(duì)比(溫度變化ΔT=150 K),如表4所示。兩種方法得到的頻率結(jié)果可以較好地匹配,最大偏差在0.8%以內(nèi)。偏差的來源主要是文中方法和FEM分別采用了不同的板殼理論和求解方法體系。此外,F(xiàn)EM計(jì)算結(jié)果的求解精度依賴于網(wǎng)格密度,需要不斷地細(xì)化網(wǎng)格才能獲得更加準(zhǔn)確的結(jié)果,但同時(shí)也會(huì)帶來計(jì)算成本的增加。兩種方法求解獲得的模態(tài)振型對(duì)比情況,如圖2所示。其中溫度變化ΔT=150 K, 梯度指數(shù)p=0.5, 邊界條件為C-C。從圖2可以看出,基于兩種方法獲得的模態(tài)振型一致性較好。綜上所述,文中方法可以有效地分析熱環(huán)境下功能梯度環(huán)板自由振動(dòng)特性。

    表1 具有不同邊界條件和梯度指數(shù)的功能梯度環(huán)板的第1階固有頻率Tab.1 The first natural frequency of functionally graded annular plate with different boundary conditions and gradient indexes

    2.2 參數(shù)化分析

    在驗(yàn)證文中構(gòu)建模型準(zhǔn)確性和通用性的基礎(chǔ)上,針對(duì)航空航天飛行器工程運(yùn)用中的各種復(fù)雜工況,本節(jié)對(duì)功能梯度環(huán)板在熱環(huán)境下的振動(dòng)特性開展參數(shù)化分析。結(jié)構(gòu)的幾何尺寸以及材料的選擇與表3、表4保持一致。

    首先,當(dāng)溫度變化ΔT=50 K時(shí),各種復(fù)雜邊界條件下功能梯度環(huán)板固有頻率,如表5所示,材料梯度指數(shù)p=1。從表5可以看出,對(duì)于經(jīng)典邊界,C-SD邊界下的固有頻率最高,SS-SD和SD-SD邊界下的頻率結(jié)果非常接近,而F邊界的存在會(huì)很明顯地降低結(jié)構(gòu)的固有頻率。而對(duì)于彈性邊界,在E1-E1和E4-E4邊界條下環(huán)板分別有著最高和最低的固有頻率。

    溫度變化ΔT對(duì)功能梯度環(huán)板固有頻率(n=1,m=1)的影響情況,如圖3所示。圖3中,材料的梯度指數(shù)p=1。不難看出,在經(jīng)典、經(jīng)典-彈性以及彈性邊界條件下,隨著ΔT增大,環(huán)板的固有頻率會(huì)不斷降低。這是因?yàn)棣變大降低應(yīng)變能導(dǎo)致了功能梯度環(huán)板剛度降低,從而使結(jié)構(gòu)固有頻率減小。

    梯度指數(shù)p以及周向波數(shù)n的變化對(duì)環(huán)板固有頻率的影響,如圖4所示,溫度變化ΔT=50 K。從圖4可以發(fā)現(xiàn),結(jié)構(gòu)固有頻率會(huì)隨著周向波數(shù)n的增加而不斷增加。而隨著梯度指數(shù)p的增加,功能梯度環(huán)板的固有頻率先是迅速下降,達(dá)到一定值后,影響逐漸不明顯。因?yàn)殡S著p的增加,環(huán)板材料屬性更加接近于下表面的金屬組分,從而降低了結(jié)構(gòu)的固有頻率。

    圖4 周向波數(shù)n和梯度指數(shù)的變化對(duì)功能梯度環(huán)板固有頻率的影響Fig.4 The influence of changes in circumferential wavenumber n and gradient indexes on the natural frequencies of functionally graded annular plates

    3 結(jié) 論

    綜合考慮功能梯度環(huán)板材料屬性溫度相關(guān)性及其在厚度方向上梯度連續(xù)變化,本文在一階剪切變形理論框架下推導(dǎo)出功能梯度環(huán)板振動(dòng)能量方程,并采用譜幾何法求解得到功能梯度環(huán)板在熱場(chǎng)中的自由振動(dòng)特性。通過將文中計(jì)算結(jié)果與相關(guān)文獻(xiàn)和有限元法結(jié)果進(jìn)行對(duì)比,驗(yàn)證了所構(gòu)建分析模型的正確性和可靠性。在此基礎(chǔ)上,研究了邊界約束條件、溫度變化和梯度指數(shù)等參數(shù)對(duì)熱環(huán)境下功能梯度環(huán)板自由振動(dòng)特性的影響,獲得以下結(jié)論:

    (1) 文中所構(gòu)建的計(jì)算分析模型可以有效分析預(yù)示熱環(huán)境下功能梯度環(huán)板自由振動(dòng)特性。

    (2) 對(duì)于經(jīng)典邊界,SS和SD邊界對(duì)功能梯度環(huán)板自由振動(dòng)固有頻率的影響類似,F(xiàn)邊界則會(huì)很明顯地降低結(jié)構(gòu)固有頻率。對(duì)于彈性邊界,E1和E4邊界下的功能梯度環(huán)板分別有著最高和最低的固有頻率。

    (3) 功能梯度環(huán)板固有頻率隨著溫度變化ΔT的增加而降低;而梯度指數(shù)p的增加,則會(huì)導(dǎo)致功能梯度環(huán)板的固有頻率呈現(xiàn)先迅速下降后逐漸穩(wěn)定的趨勢(shì)。

    猜你喜歡
    環(huán)板固有頻率梯度
    一個(gè)改進(jìn)的WYL型三項(xiàng)共軛梯度法
    現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
    湖北某電站發(fā)電機(jī)轉(zhuǎn)子支架環(huán)板變形原因剖析及處理
    一種自適應(yīng)Dai-Liao共軛梯度法
    敦化水輪機(jī)座環(huán)開口高度的控制措施
    一類扭積形式的梯度近Ricci孤立子
    變厚度圓板、環(huán)板振動(dòng)分析的傳遞矩陣法*
    總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
    EPR核電站RIS系統(tǒng)貫穿件與膨脹節(jié)環(huán)板焊接變形控制
    電焊機(jī)(2015年8期)2015-01-09 09:44:02
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    满18在线观看网站| 人妻久久中文字幕网| 精品国产乱子伦一区二区三区 | 亚洲av欧美aⅴ国产| 最新在线观看一区二区三区| 美女福利国产在线| 成人国产一区最新在线观看| 一二三四社区在线视频社区8| 欧美日韩亚洲综合一区二区三区_| 制服诱惑二区| 啦啦啦中文免费视频观看日本| 一边摸一边做爽爽视频免费| 日韩 欧美 亚洲 中文字幕| 美女高潮到喷水免费观看| 国产色视频综合| 法律面前人人平等表现在哪些方面 | 乱人伦中国视频| 精品人妻熟女毛片av久久网站| 视频在线观看一区二区三区| √禁漫天堂资源中文www| 欧美精品啪啪一区二区三区 | 久久精品成人免费网站| 首页视频小说图片口味搜索| 热re99久久国产66热| 免费在线观看影片大全网站| 91麻豆精品激情在线观看国产 | 国产视频一区二区在线看| av在线老鸭窝| 亚洲黑人精品在线| 一级毛片电影观看| 久久久久国产精品人妻一区二区| 日本欧美视频一区| 电影成人av| 大片电影免费在线观看免费| 老司机亚洲免费影院| 精品亚洲成国产av| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 亚洲精品国产av成人精品| 亚洲av成人不卡在线观看播放网 | 亚洲午夜精品一区,二区,三区| 日韩视频在线欧美| 啪啪无遮挡十八禁网站| av一本久久久久| 国产精品久久久久久精品电影小说| 欧美精品高潮呻吟av久久| 一个人免费在线观看的高清视频 | 亚洲精品在线美女| 日韩三级视频一区二区三区| 69精品国产乱码久久久| 亚洲av成人一区二区三| 午夜福利视频在线观看免费| 亚洲欧美一区二区三区黑人| 国产熟女午夜一区二区三区| 亚洲色图综合在线观看| 99久久人妻综合| 国产亚洲欧美在线一区二区| 精品一区二区三区四区五区乱码| 99久久99久久久精品蜜桃| 精品亚洲乱码少妇综合久久| 午夜福利在线免费观看网站| 人人妻人人澡人人看| 日日摸夜夜添夜夜添小说| 欧美日韩亚洲高清精品| 99精国产麻豆久久婷婷| 老司机福利观看| 欧美国产精品va在线观看不卡| 亚洲国产毛片av蜜桃av| 久久免费观看电影| svipshipincom国产片| 日韩 亚洲 欧美在线| 精品一品国产午夜福利视频| a级毛片在线看网站| 黄色 视频免费看| 午夜免费观看性视频| 国产成人精品无人区| 高清av免费在线| 男人爽女人下面视频在线观看| 国产精品香港三级国产av潘金莲| 久久精品国产亚洲av高清一级| 精品欧美一区二区三区在线| 久久精品亚洲熟妇少妇任你| 在线观看免费视频网站a站| 国产精品久久久久久人妻精品电影 | 男女国产视频网站| 国产精品国产三级国产专区5o| 亚洲成人国产一区在线观看| 亚洲精华国产精华精| 国产片内射在线| 午夜两性在线视频| 日韩,欧美,国产一区二区三区| av福利片在线| 亚洲av片天天在线观看| 巨乳人妻的诱惑在线观看| 天堂中文最新版在线下载| 大片电影免费在线观看免费| 亚洲欧美日韩另类电影网站| 亚洲国产精品999| 日本av免费视频播放| 精品国内亚洲2022精品成人 | 亚洲精品av麻豆狂野| 狠狠精品人妻久久久久久综合| 女人久久www免费人成看片| 久久久精品国产亚洲av高清涩受| xxxhd国产人妻xxx| 国精品久久久久久国模美| 久久久国产一区二区| 在线av久久热| 老司机在亚洲福利影院| 国产在线视频一区二区| 日韩大片免费观看网站| 岛国在线观看网站| 免费在线观看完整版高清| 亚洲熟女精品中文字幕| 亚洲av成人一区二区三| 国产一区有黄有色的免费视频| 满18在线观看网站| 精品国内亚洲2022精品成人 | 国产成人免费无遮挡视频| 亚洲 国产 在线| 99精国产麻豆久久婷婷| 欧美精品av麻豆av| 亚洲全国av大片| 国产精品自产拍在线观看55亚洲 | 人人妻人人澡人人看| 亚洲av成人一区二区三| 久久毛片免费看一区二区三区| 91国产中文字幕| 啦啦啦在线免费观看视频4| 亚洲全国av大片| 91av网站免费观看| 国产高清国产精品国产三级| 欧美日韩亚洲高清精品| 欧美精品一区二区大全| 男女国产视频网站| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 精品第一国产精品| 久久久国产欧美日韩av| 中文欧美无线码| 国产99久久九九免费精品| 丝瓜视频免费看黄片| 精品免费久久久久久久清纯 | 各种免费的搞黄视频| 下体分泌物呈黄色| 欧美黑人欧美精品刺激| 亚洲国产av新网站| 黄色视频在线播放观看不卡| 一区福利在线观看| 中文字幕人妻熟女乱码| 黄色视频在线播放观看不卡| 国产在线视频一区二区| 欧美另类亚洲清纯唯美| 日韩制服丝袜自拍偷拍| 1024视频免费在线观看| netflix在线观看网站| 精品国产一区二区三区四区第35| 免费人妻精品一区二区三区视频| 热re99久久国产66热| 日韩 欧美 亚洲 中文字幕| 免费人妻精品一区二区三区视频| 黄色毛片三级朝国网站| 午夜福利乱码中文字幕| 亚洲人成电影免费在线| 亚洲国产精品成人久久小说| 亚洲精品美女久久久久99蜜臀| 精品国产一区二区久久| 曰老女人黄片| 日韩熟女老妇一区二区性免费视频| 一本—道久久a久久精品蜜桃钙片| 亚洲国产精品成人久久小说| 男女床上黄色一级片免费看| 亚洲av男天堂| 热re99久久精品国产66热6| 国产精品.久久久| 欧美日韩亚洲国产一区二区在线观看 | 国产精品成人在线| 高清欧美精品videossex| 精品一区二区三卡| 精品国产一区二区久久| 国产精品国产三级国产专区5o| 成人亚洲精品一区在线观看| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三区在线| 大型av网站在线播放| 亚洲情色 制服丝袜| 国产野战对白在线观看| 日本av免费视频播放| 一本综合久久免费| 在线精品无人区一区二区三| 狂野欧美激情性xxxx| 国产精品国产三级国产专区5o| 人人妻人人澡人人爽人人夜夜| 日本猛色少妇xxxxx猛交久久| 黄色 视频免费看| 日本a在线网址| 9色porny在线观看| 日韩 欧美 亚洲 中文字幕| 制服诱惑二区| 欧美精品啪啪一区二区三区 | 国产高清videossex| 每晚都被弄得嗷嗷叫到高潮| 大型av网站在线播放| 国产99久久九九免费精品| 国产真人三级小视频在线观看| 欧美成狂野欧美在线观看| 日韩免费高清中文字幕av| 亚洲视频免费观看视频| 精品久久蜜臀av无| 男人舔女人的私密视频| 国产精品1区2区在线观看. | 亚洲欧美日韩另类电影网站| 90打野战视频偷拍视频| 国产一区二区 视频在线| 日本91视频免费播放| 中国美女看黄片| 精品卡一卡二卡四卡免费| 国产免费av片在线观看野外av| 精品国产一区二区久久| 999精品在线视频| 亚洲av日韩在线播放| 制服人妻中文乱码| 欧美在线黄色| 免费不卡黄色视频| 黄色 视频免费看| 天天操日日干夜夜撸| 啦啦啦中文免费视频观看日本| 日韩电影二区| 亚洲精品乱久久久久久| 亚洲av成人不卡在线观看播放网 | 亚洲av成人不卡在线观看播放网 | av线在线观看网站| 日本av免费视频播放| 老司机亚洲免费影院| 欧美激情极品国产一区二区三区| 国产精品 欧美亚洲| 人人妻,人人澡人人爽秒播| 中亚洲国语对白在线视频| netflix在线观看网站| 免费在线观看完整版高清| 久久久久国产一级毛片高清牌| 婷婷色av中文字幕| 十八禁人妻一区二区| 91九色精品人成在线观看| 欧美激情 高清一区二区三区| 十分钟在线观看高清视频www| 亚洲精品国产区一区二| 亚洲一区二区三区欧美精品| 精品国产国语对白av| 成人国产一区最新在线观看| 涩涩av久久男人的天堂| 亚洲视频免费观看视频| 亚洲国产欧美一区二区综合| av片东京热男人的天堂| 国产亚洲精品一区二区www | 十分钟在线观看高清视频www| 精品卡一卡二卡四卡免费| 嫩草影视91久久| 美女高潮到喷水免费观看| 国产极品粉嫩免费观看在线| 天堂8中文在线网| 欧美日韩黄片免| 十八禁网站网址无遮挡| 亚洲人成电影免费在线| 看免费av毛片| 新久久久久国产一级毛片| 亚洲九九香蕉| 亚洲精品国产一区二区精华液| 自线自在国产av| 十八禁人妻一区二区| 欧美精品啪啪一区二区三区 | 精品第一国产精品| 热99久久久久精品小说推荐| 中文字幕色久视频| 99国产极品粉嫩在线观看| 嫁个100分男人电影在线观看| 欧美 亚洲 国产 日韩一| 国产亚洲欧美在线一区二区| 天堂8中文在线网| 精品人妻在线不人妻| 亚洲精品乱久久久久久| 丁香六月天网| 亚洲一卡2卡3卡4卡5卡精品中文| 黄网站色视频无遮挡免费观看| 欧美xxⅹ黑人| 老司机深夜福利视频在线观看 | 精品久久久精品久久久| 国产一区二区三区在线臀色熟女 | 国内毛片毛片毛片毛片毛片| 国产免费视频播放在线视频| 热99久久久久精品小说推荐| 免费久久久久久久精品成人欧美视频| 国产真人三级小视频在线观看| 性少妇av在线| 热re99久久精品国产66热6| 一级片免费观看大全| 亚洲自偷自拍图片 自拍| www.av在线官网国产| 99国产精品一区二区三区| 在线天堂中文资源库| 久久久久精品国产欧美久久久 | 操出白浆在线播放| 热re99久久精品国产66热6| 美女大奶头黄色视频| 夜夜骑夜夜射夜夜干| 五月开心婷婷网| 制服人妻中文乱码| 亚洲,欧美精品.| 中文字幕色久视频| 黄频高清免费视频| 淫妇啪啪啪对白视频 | 亚洲国产av影院在线观看| 午夜激情av网站| 久久人人爽人人片av| 国产99久久九九免费精品| 777久久人妻少妇嫩草av网站| 少妇被粗大的猛进出69影院| 色精品久久人妻99蜜桃| 久久青草综合色| 日韩一区二区三区影片| 久久精品aⅴ一区二区三区四区| 亚洲黑人精品在线| 淫妇啪啪啪对白视频 | 亚洲国产毛片av蜜桃av| 永久免费av网站大全| 久久午夜综合久久蜜桃| 亚洲精品国产av蜜桃| 美女脱内裤让男人舔精品视频| 亚洲国产中文字幕在线视频| 国产一级毛片在线| 老司机影院毛片| 中文字幕色久视频| 激情视频va一区二区三区| 老司机影院毛片| 伦理电影免费视频| 999久久久精品免费观看国产| 欧美激情极品国产一区二区三区| 超色免费av| 亚洲情色 制服丝袜| 国产欧美日韩一区二区三区在线| 国产在线免费精品| 中文字幕av电影在线播放| 欧美黄色片欧美黄色片| 50天的宝宝边吃奶边哭怎么回事| 深夜精品福利| 啦啦啦中文免费视频观看日本| 欧美黑人欧美精品刺激| 免费在线观看日本一区| 啦啦啦 在线观看视频| 午夜免费鲁丝| www.av在线官网国产| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品大桥未久av| 欧美日韩成人在线一区二区| 黄片播放在线免费| 亚洲精品中文字幕在线视频| 日本欧美视频一区| 69精品国产乱码久久久| 少妇粗大呻吟视频| 亚洲人成电影免费在线| av线在线观看网站| 成人国语在线视频| 精品国产一区二区三区四区第35| 熟女少妇亚洲综合色aaa.| 欧美在线黄色| 一进一出抽搐动态| 亚洲熟女毛片儿| 视频区欧美日本亚洲| 国产男女内射视频| 国产精品国产三级国产专区5o| 亚洲精品粉嫩美女一区| 日韩制服骚丝袜av| 欧美性长视频在线观看| 国产xxxxx性猛交| 交换朋友夫妻互换小说| 老熟妇乱子伦视频在线观看 | 久久热在线av| 在线天堂中文资源库| 窝窝影院91人妻| 日韩,欧美,国产一区二区三区| av国产精品久久久久影院| 久久久久精品人妻al黑| 午夜福利视频精品| 黑丝袜美女国产一区| 啦啦啦视频在线资源免费观看| 久久久国产欧美日韩av| 亚洲黑人精品在线| 久久国产精品人妻蜜桃| 亚洲欧美日韩另类电影网站| 亚洲国产精品999| 一本—道久久a久久精品蜜桃钙片| 亚洲精品中文字幕在线视频| 在线看a的网站| 国产99久久九九免费精品| 午夜免费观看性视频| 少妇人妻久久综合中文| 精品国产超薄肉色丝袜足j| 高清黄色对白视频在线免费看| 99国产精品99久久久久| 1024香蕉在线观看| 成年女人毛片免费观看观看9 | 午夜福利视频精品| 国产成人欧美| 一个人免费在线观看的高清视频 | 天堂俺去俺来也www色官网| 国产亚洲精品久久久久5区| 高潮久久久久久久久久久不卡| 精品一区二区三卡| 精品国产一区二区久久| 热re99久久精品国产66热6| 国产福利在线免费观看视频| 色播在线永久视频| 男女床上黄色一级片免费看| 亚洲欧洲日产国产| 在线观看人妻少妇| 久久久国产成人免费| 亚洲精品第二区| 精品国产乱码久久久久久男人| 日本av手机在线免费观看| 欧美黑人精品巨大| 国产成人欧美在线观看 | 色婷婷av一区二区三区视频| 国产老妇伦熟女老妇高清| 国产日韩欧美在线精品| cao死你这个sao货| 脱女人内裤的视频| 少妇被粗大的猛进出69影院| 夜夜夜夜夜久久久久| 飞空精品影院首页| 桃红色精品国产亚洲av| 美女国产高潮福利片在线看| 亚洲av电影在线进入| 淫妇啪啪啪对白视频 | 成年人黄色毛片网站| 久久久久国产精品人妻一区二区| 黄频高清免费视频| 久久精品国产综合久久久| 老汉色av国产亚洲站长工具| 成人三级做爰电影| 巨乳人妻的诱惑在线观看| 欧美少妇被猛烈插入视频| 欧美午夜高清在线| 黄色毛片三级朝国网站| 天堂中文最新版在线下载| 亚洲精品一卡2卡三卡4卡5卡 | 97人妻天天添夜夜摸| 久久性视频一级片| 午夜免费成人在线视频| 老司机午夜福利在线观看视频 | 国产日韩欧美亚洲二区| 免费人妻精品一区二区三区视频| 高清欧美精品videossex| 在线精品无人区一区二区三| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲全国av大片| 亚洲色图 男人天堂 中文字幕| 亚洲五月色婷婷综合| 亚洲av日韩在线播放| 欧美乱码精品一区二区三区| 大香蕉久久网| 青青草视频在线视频观看| 久久亚洲精品不卡| 亚洲第一欧美日韩一区二区三区 | 丝袜脚勾引网站| 亚洲七黄色美女视频| 国产在视频线精品| 一个人免费在线观看的高清视频 | 国产亚洲精品第一综合不卡| 久久精品久久久久久噜噜老黄| 韩国高清视频一区二区三区| 美女主播在线视频| 在线观看一区二区三区激情| 久久午夜综合久久蜜桃| 国产成人精品在线电影| a 毛片基地| 亚洲熟女精品中文字幕| 免费黄频网站在线观看国产| 国产片内射在线| 欧美激情极品国产一区二区三区| 91精品伊人久久大香线蕉| 日本精品一区二区三区蜜桃| 成人18禁高潮啪啪吃奶动态图| 人妻一区二区av| av天堂久久9| 久久中文看片网| av一本久久久久| 亚洲男人天堂网一区| 国产免费视频播放在线视频| 69av精品久久久久久 | 亚洲精品美女久久av网站| netflix在线观看网站| 伊人亚洲综合成人网| 老司机靠b影院| 在线av久久热| 国产伦理片在线播放av一区| 色播在线永久视频| 色精品久久人妻99蜜桃| 成人国产一区最新在线观看| 法律面前人人平等表现在哪些方面 | 日本一区二区免费在线视频| 亚洲国产毛片av蜜桃av| 一本久久精品| 最近最新免费中文字幕在线| 亚洲全国av大片| 高清av免费在线| 老汉色av国产亚洲站长工具| av不卡在线播放| 中文欧美无线码| 亚洲成人免费av在线播放| 久久ye,这里只有精品| 色94色欧美一区二区| 美女国产高潮福利片在线看| 99香蕉大伊视频| 国产日韩欧美亚洲二区| 制服人妻中文乱码| 美女脱内裤让男人舔精品视频| 又紧又爽又黄一区二区| 亚洲伊人久久精品综合| 午夜福利一区二区在线看| 国产1区2区3区精品| 天天操日日干夜夜撸| 女人被躁到高潮嗷嗷叫费观| 日韩一卡2卡3卡4卡2021年| 精品亚洲乱码少妇综合久久| 国产成人免费观看mmmm| 精品亚洲成国产av| 色婷婷av一区二区三区视频| 老熟女久久久| 免费在线观看黄色视频的| 12—13女人毛片做爰片一| 69精品国产乱码久久久| 欧美97在线视频| 9色porny在线观看| 国产亚洲精品久久久久5区| 亚洲精品av麻豆狂野| 高清黄色对白视频在线免费看| 亚洲精品中文字幕在线视频| 777米奇影视久久| 亚洲国产日韩一区二区| 9热在线视频观看99| 美女主播在线视频| h视频一区二区三区| 啦啦啦在线免费观看视频4| 纵有疾风起免费观看全集完整版| 亚洲精品国产av蜜桃| 国产成人啪精品午夜网站| 国产一区二区激情短视频 | 亚洲精品日韩在线中文字幕| 久久精品国产a三级三级三级| 欧美国产精品一级二级三级| 国产一区二区三区在线臀色熟女 | 日韩一卡2卡3卡4卡2021年| 欧美午夜高清在线| 一本久久精品| 欧美日韩视频精品一区| 在线观看一区二区三区激情| av在线播放精品| 伊人亚洲综合成人网| 欧美+亚洲+日韩+国产| 一级毛片电影观看| 黄色视频,在线免费观看| 精品少妇一区二区三区视频日本电影| 夜夜夜夜夜久久久久| 男女床上黄色一级片免费看| 在线 av 中文字幕| 欧美日韩一级在线毛片| 午夜免费观看性视频| 大片免费播放器 马上看| 汤姆久久久久久久影院中文字幕| 波多野结衣av一区二区av| 黄色怎么调成土黄色| 人成视频在线观看免费观看| 日韩制服丝袜自拍偷拍| 青春草亚洲视频在线观看| 精品免费久久久久久久清纯 | 亚洲成人免费电影在线观看| 最新在线观看一区二区三区| 十八禁网站网址无遮挡| 97人妻天天添夜夜摸| 日本撒尿小便嘘嘘汇集6| 亚洲男人天堂网一区| 欧美精品啪啪一区二区三区 | 国产无遮挡羞羞视频在线观看| 一级片'在线观看视频| 女人久久www免费人成看片| 黑人巨大精品欧美一区二区mp4| 国产成人免费观看mmmm| 老汉色av国产亚洲站长工具| 欧美精品亚洲一区二区| 少妇人妻久久综合中文| 中文字幕人妻丝袜一区二区| 亚洲午夜精品一区,二区,三区| 亚洲国产av新网站| 男女床上黄色一级片免费看| 曰老女人黄片| 99热网站在线观看| 日本精品一区二区三区蜜桃| a 毛片基地| 亚洲免费av在线视频| 欧美日本中文国产一区发布| 久久久久国产精品人妻一区二区| 国产精品亚洲av一区麻豆| 岛国毛片在线播放| 精品少妇一区二区三区视频日本电影| 成年女人毛片免费观看观看9 | 国产在线一区二区三区精| 亚洲av电影在线观看一区二区三区| 日本a在线网址| 亚洲av成人不卡在线观看播放网 | 久久中文字幕一级| 国产高清视频在线播放一区 | 一个人免费在线观看的高清视频 | 中国美女看黄片| 搡老岳熟女国产| 美女福利国产在线|