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

    熱環(huán)境中的功能梯度扁球殼非線性穩(wěn)定性分析

    2020-08-28 02:30:26趙偉東高士武黃永玉
    工程力學(xué) 2020年8期
    關(guān)鍵詞:球殼構(gòu)形變溫

    趙偉東,高士武,黃永玉

    (1. 青海大學(xué)土木工程學(xué)院,西寧 810016;2. 青海大學(xué)機(jī)械工程學(xué)院,西寧 810016)

    薄壁彈性扁球殼結(jié)構(gòu)在核反應(yīng)、石油工業(yè)、艦船、航空航天、精密儀器等領(lǐng)域有著廣泛的應(yīng)用。為使此類結(jié)構(gòu)能夠可靠的服役,需要對其屈曲行為進(jìn)行精確分析。早在20 世紀(jì)30 年代,von Karman 等[1]首次用非線性分析方法分析了扁球殼的屈曲行為。在此之后,用非線性分析方法分析薄殼的屈曲行為吸引了諸多學(xué)者的注意。Huang[2]研究了均質(zhì)扁球殼在外壓力作用下的非對稱屈曲問題,并指出非對稱變形可能是軸對稱屈曲理論與試驗(yàn)臨界壓力不一致的原因。葉開源等[3]提出了修正迭代法,并用該方法成功地研究了均質(zhì)扁球殼的非線性穩(wěn)定問題。Uemura[4]采用兩項(xiàng)撓度近似法處理了均布外壓作用的周邊夾緊帶缺陷均質(zhì)材料扁球殼的軸對稱屈曲問題,并利用總勢能的二次變分評估了球殼矢高對解的非唯一性的影響。Yan[5]采用分步加載法研究了均布外壓作用下變厚度開頂扁球殼的非線性穩(wěn)定性。Nie[6]提出了一種漸近迭代法,并用該方法分析了外側(cè)壓力作用的均勻材料扁球殼的非線性屈曲問題,分析過程考慮了缺陷、邊界彈性約束和彈性地基的影響。Shahsiah 等[7]運(yùn)用解析方法,考慮熱-機(jī)荷載聯(lián)合作用,研究了均質(zhì)扁球殼線性屈曲行為,并考察了缺陷對解的的影響。Li 等[8]基于一階剪切變形殼理論,運(yùn)用修正迭代法研究了均勻外壓下均質(zhì)扁薄球殼的屈曲行為??紤]橫向剪切變形效應(yīng),Zhu 等[9]考察了對稱層合圓柱正交各向異性扁球殼的熱屈曲行為。Xu 等[10]采用修正迭代法研究了雙層網(wǎng)狀圓形扁球殼在均布外壓作用下的屈曲行為。運(yùn)用有限元方法,張平等[11]對用于高速戰(zhàn)斗機(jī)自適應(yīng)進(jìn)氣道的均質(zhì)扁球殼的雙穩(wěn)態(tài)特性進(jìn)行了討論。李忱等[12]針對均布外壓與溫度耦合作用的均質(zhì)薄球殼的屈曲行為開展了研究。趙偉東等[13]基于數(shù)值打靶法分析了熱環(huán)境中的周邊簡支均質(zhì)扁球殼在外側(cè)均布壓力作用下的穩(wěn)定性問題。

    自從功能梯度材料[14]問世以來,對功能梯度材料板殼結(jié)構(gòu)力學(xué)行為的研究吸引了很多研究人員的關(guān)注。Woo 等[15]給出了在壓力和溫度場作用下的功能梯度材料板和淺圓柱殼的后屈曲行為的解析解。并指出熱-機(jī)耦合效應(yīng)和邊界條件對功能梯度板/殼在邊緣壓力作用下的響應(yīng)中起著重要作用。Shahsiah 等[16]基于Sanders 應(yīng)力-應(yīng)變關(guān)系和變分法研究了功能梯度扁球殼的熱彈性失穩(wěn)。Prakash 等[17]和Ganapathi[18]基于一階剪切變形理論,用有限元方法研究了固支邊功能梯度扁球殼的非線性軸對稱動力穩(wěn)定性問題。Bich 等[19]利用經(jīng)典殼理論研究了功能梯度扁球殼在均勻外壓(包括溫度效應(yīng))作用下的非線性軸對稱響應(yīng),研究過程考慮了球殼缺陷對其響應(yīng)的影響。Boroujerdy等[20 ? 21]研究了壓電-功能梯度扁球殼的非線性軸對稱熱-機(jī)械響應(yīng)和熱屈曲。Mao 等[22]考察了功能梯度扁球殼在低速沖擊下的非線性動力響應(yīng)和損傷行為。Fu 等[23]對功能梯度扁球殼在隨時間變化的熱機(jī)械荷載作用下的非線性瞬態(tài)響應(yīng)進(jìn)行了分析。

    用有限單元法分析功能梯度材料板殼結(jié)構(gòu),為了能夠反映組分材料特性,通常需要沿厚度方向劃分單元,這使求解過程變得復(fù)雜。另外,要想用有限單元法得到殼體屈曲過程的完整的荷載-撓度曲線,其計(jì)算工作量非常繁重。邵玉龍等[24]的工作表明:用二階一致無網(wǎng)格法對功能梯度材料進(jìn)行數(shù)值分析,有較高的計(jì)算效率。張鵬飛等[25]用有限質(zhì)點(diǎn)法分析了薄殼屈曲問題,該方法有望成為分析復(fù)雜邊界功能梯度扁球殼屈曲完整過程的有效數(shù)值方法??紤]到工程中較為常見的周邊固支或簡支功能梯度圓底扁球殼是軸對稱結(jié)構(gòu),若荷載和溫度場也是軸對稱的(事實(shí)上這種情況較常見,如外側(cè)均布壓力和沿厚度均勻或一維熱傳導(dǎo)溫度場),在這種情況下,對于滿足一定幾何條件的扁球殼,其彎曲乃至屈曲變形經(jīng)常是軸對稱的?;诖?,如果利用球殼中曲面的橫向和徑向位移函數(shù)建立數(shù)學(xué)模型,則扁球殼軸對稱彎曲/屈曲控制方程和邊界條件為典型的常微分方程兩點(diǎn)邊值問題。趙偉東等[26]曾經(jīng)運(yùn)用變分原理推導(dǎo)了均勻溫度場中的不可移夾緊邊功能梯度圓底扁球殼在外側(cè)均布壓力作用下的軸對稱彎曲/屈曲幾何非線性常微分控制方程及其邊界條件,并借助打靶法得到了一些對工程設(shè)計(jì)有幫助的數(shù)值結(jié)果。根據(jù)對文獻(xiàn)的調(diào)閱可以發(fā)現(xiàn),到目前為止,可以用來參考的簡支邊功能梯度圓底扁球殼彎曲/屈曲的研究數(shù)據(jù)還比較稀缺。因此,本文在文獻(xiàn)[26]工作的基礎(chǔ)上,基于經(jīng)典殼理論和von Karman 幾何非線性理論,借用球殼平衡微分方程,導(dǎo)出了均勻外壓與均勻溫度場聯(lián)合作用下的功能梯度圓底扁球殼軸對稱變形的位移型幾何非線性控制方程,推導(dǎo)了功能梯度圓底扁球殼位移型簡支邊界條件。借用數(shù)值打靶法求解了該兩點(diǎn)邊值問題,分析了簡支邊功能梯度圓底扁球殼的幾何非線性力學(xué)行為??紤]到工程界更關(guān)注扁球殼的承載力(即上臨界載荷),文中提供了此類結(jié)構(gòu)承載力的兩個數(shù)表和一些有益的數(shù)值曲線,供設(shè)計(jì)人員參考。

    1 問題的基本方程

    如圖1 所示簡支邊功能梯度圓底扁球殼,假設(shè)其材料性質(zhì)沿厚度方向連續(xù)變化。其中曲面曲率半徑、殼體厚度和底圓半徑分別為R、h 和a。球坐標(biāo)系(φ,θ,z)分別為子午、圓周和厚度坐標(biāo),z(?h/2≤z≤h/2)坐標(biāo)軸正向與球殼中曲面的外法線方向相同,坐標(biāo)原點(diǎn)位于中曲面的中心。

    1.1 幾何方程

    對于扁球殼,考慮球殼底圓徑向坐標(biāo)r(r=Rsinφ)及近似關(guān)系dr=Rdφ(cosφ≈1)會給方程的推導(dǎo)和求解帶來方便。在軸對稱變形和基于上述坐標(biāo)系的情況下,扁球殼中曲面上任意點(diǎn)A(φ,θ,0)的橫向和徑向位移W 和U 僅是半徑r 的一元函數(shù)。根據(jù)經(jīng)典殼理論與von Karman 幾何非線性理論,殼體上任意點(diǎn)A(φ,θ,z)處的應(yīng)變分量與殼體中曲面位移分量和橫向坐標(biāo)z 的關(guān)系為[20]:

    式中,( ·),r表示( ·)對r 的一階導(dǎo)數(shù)。

    圖1 扁球殼的幾何形狀和坐標(biāo)系Fig.1 Geometry and the coordinate system of a shallow spherical shell

    1.2 物理方程

    計(jì)及溫度應(yīng)力的球殼線彈性應(yīng)力-應(yīng)變關(guān)系為:

    式中:T 為球殼相對初始無應(yīng)力狀態(tài)的均勻變溫場;彈性模量E、熱膨脹系數(shù)α 和泊松比ν 為坐標(biāo)z 的連續(xù)函數(shù)。通常功能梯度板的泊松比變化不大,本文將其假設(shè)為常數(shù)。

    球殼的內(nèi)力分量與應(yīng)力分量及橫向坐標(biāo)的關(guān)系為:

    式中:Nr和Nθ分別指殼體徑向和環(huán)向截面內(nèi)單位寬度上的薄膜力;Mr和Mθ分別指殼體徑向和環(huán)向截面內(nèi)單位寬度上的彎矩。將式(1)代入式(2),再將式(2)代入式(3),沿厚度方向積分可得用位移表示的殼體薄膜力和彎矩表達(dá)式:

    其中,剛度系數(shù)定義分別為:

    上述剛度系數(shù)可分別表示為:

    式(4)~式(7)中的熱薄膜力S3和熱彎矩S4可分別表示為:

    其中:

    式中:Am為參考均勻金屬球殼的拉伸剛度; ?i為 無量綱系數(shù)。

    球殼的彈性模量及熱膨脹系數(shù)可表示為:

    式中:ξ=z/h; Em、 αm分別為參考均勻金屬球殼的彈性模量及熱膨脹系數(shù);ψi(ξ)為坐標(biāo) ξ的連續(xù)函數(shù)。于是可給出無量綱系數(shù)?i的計(jì)算公式:

    考慮材料性質(zhì)在橫向按冪率函數(shù)變化(需要指出,如果材料性質(zhì)按其他函數(shù)形式變化,只要函數(shù)僅是厚度坐標(biāo)z 或ξ 的一元函數(shù),在軸對稱熱機(jī)荷載的情況下,本文的求解方法依然適用),則函數(shù)ψi(ξ)可表示為[27]:

    式中:η1=(Ec?Em)/Em;η2=(αc?αm)/αm;0≤k<∞; Ec=E(h/2);αc=α(h/2)。

    將式(14)代入式(13),可得無量綱系數(shù)?i的 解析表達(dá)式如下:

    1.3 平衡方程

    在軸對稱的情況下,殼體平衡微分方程可以由文獻(xiàn)[20]給出的相應(yīng)方程退化為如下形式:

    式中,q 為橫向均勻分布荷載,沿z 坐標(biāo)軸方向?yàn)檎?/p>

    2 兩點(diǎn)邊值問題

    2.1 位移型控制方程

    為便于討論,采用如下無量綱變換:

    式中:x 為無量綱的徑向坐標(biāo);w 和u 分別為球殼中曲面無量綱的橫向和徑向位移; δ為球殼底圓半徑a 與球殼厚度h 的比; γ為球殼底圓半徑a 與球殼中曲面曲率半徑R 的比;μ為球殼相對初始無應(yīng)力狀態(tài)的無量綱均勻變溫場;Q 為無量綱的法向均勻分布力集度。

    將式(9)~式(11)代入式(4)~式(7),再將結(jié)果代入式(20)~式(21),考慮無量綱變換式(22),得到一組無量綱形式的位移型幾何非線性常微分控制方程:

    式中,?2(·)=(·),x/x+(·),xx,(·),x為 (·)對x 的一階導(dǎo)數(shù)。

    2.2 位移型簡支邊界條件

    在簡單支撐邊界面上,彎矩Mr為0。將Si(i=1,2,4)代入式(6),考慮無量綱變換式(22),得到x=1 處的一個位移型邊界條件如下:

    易知,在x=1 處還有兩個位移型邊界條件如下:

    扁球殼在軸對稱變形情況下有如下平衡微分方程:

    式中,Qr為殼體單位寬度徑向截面上的剪力。

    在均布法向荷載作用下,據(jù)內(nèi)力對稱特性,在r=0 處有:

    聯(lián)立式(27)~式(28),在r=0 處有:

    將式(6)~式(7)代入式(29),并考慮Si(i=1,2,4)的結(jié)果及無量綱變換式(22),得到在x=0 處的一個位移型中心對稱條件如下:

    易知,在x=0 處還有兩個位移型中心對稱條件如下:

    3 打靶法求解說明

    式(23)~式(24)、邊界條件式(25)~式(26)及中心對稱條件式(30)~式(31)為熱環(huán)境中的簡支邊功能梯度圓底扁球殼在沿球面法向均布荷載作用下的軸對稱彎曲變形的常微分控制方程及其邊界條件,屬于典型的位移型兩點(diǎn)邊值問題??捎么虬蟹ㄈ菀椎孬@得該兩點(diǎn)邊值問題的足夠精確的數(shù)值解。

    對上述兩點(diǎn)邊值問題用打靶法求解時,可設(shè):

    式(23)、式(24)及式(30)在x=0 處存在奇異性,為使數(shù)值解收斂,需要將解域定義在區(qū)間[Δx ,1]上,其中,Δx 為能夠滿足解的精度要求的足夠小的正數(shù)。則原來的兩點(diǎn)邊值問題轉(zhuǎn)化為如下初值問題:

    式中,φ1和φ2分別由式(23)和式(24)得到。因?yàn)槲灰坪瘮?shù)在x=0 處連續(xù)可導(dǎo),當(dāng)Δx 足夠小時,可以用Δx 處的值代替x=0 處的邊界條件:

    式中:ζ 為x=Δx 處的撓度初值(賦值);V1、V2和V3為待定參數(shù),需要通過x=1 處的邊界條件確定。

    4 數(shù)值結(jié)果分析

    由于第1 節(jié)定義的橫向撓度函數(shù)W(r)及均布荷載q 均以凸面?zhèn)韧夥ň€方向一致為正,而本文重點(diǎn)考察扁球殼在凸面?zhèn)葔毫ψ饔孟碌那袨?,在這種情況下,撓度W(r)及荷載q 的數(shù)值結(jié)果均為負(fù)值。為了方便描述和容易理解,下文給出的橫向撓度W(r)和荷載集度q 在數(shù)值計(jì)算時添加了負(fù)號(也即向球殼內(nèi)凹的一側(cè)為正)。

    4.1 簡化模型數(shù)值結(jié)果對比

    由于功能梯度簡支邊圓底扁球殼屈曲數(shù)值解很難找到??紤]到當(dāng)k=0 或趨于+∞時,功能梯度材料扁球殼分別退化為彈性模量為Ec或Em的均勻材料扁球殼。需要指出的是,當(dāng)k 趨于+∞時,與Ec=Em,αc=αm及取任意的k 是等效的。現(xiàn)在考慮泊松比ν=0.3、E=Ec=Em=70 GPa 的均質(zhì)扁球殼,取球殼中曲面曲率半徑R=1600 mm、底圓半徑a=160 mm 和厚度h=5 mm,表1 給出了該均質(zhì)扁球殼在多個外側(cè)壓力集度q 作用下的球殼中心有量綱撓度W(0)的打靶法解和有限元解(表中加粗的數(shù)據(jù)為殼體上臨界荷載qu)??梢钥闯?,當(dāng)扁球殼接近臨界失穩(wěn)狀態(tài)時,打靶法解與有限元解有一定程度偏離,其余情況下,兩者一致性較好。需要指出的是,對于均勻材料扁球殼,用打靶法得到的球殼上臨界荷載(指球殼S 形荷載-撓度曲線上極大值點(diǎn)對應(yīng)的荷載。當(dāng)外側(cè)壓力從零逐漸增加并達(dá)到該值時,球殼的平衡形態(tài)會由整體外凸形態(tài)突然改變?yōu)檎w或局部內(nèi)凹形態(tài),該過程是以“跳躍”的形式實(shí)現(xiàn)的,并伴隨明顯的聲音信號)與三次樣條函數(shù)法得到的相應(yīng)解一致性很好(見文獻(xiàn)[28])。據(jù)此可以判斷,用打靶法求解扁球殼屈曲問題,解的精度是可以保證的。

    表1 均質(zhì)扁球殼在外側(cè)均布壓力q 作用下的中心撓度W(0)Table1 Uniform external pressure q & central deflection W(0) for an homogeneous shallow spherical shell

    4.2 功能梯度扁球殼屈曲平衡路徑

    為了揭示不變溫情況下簡支邊功能梯度圓底扁球殼屈曲平衡路徑,針對陶瓷及金屬組分材料分別為黏土(Ec=380 GPa)和鋁(Em=70 GPa)、泊松比均為ν=0.3,在指定指數(shù)k=1、圓底半徑a=160 mm和厚度h=5 mm 的情況下,圖2 給出了中曲面曲率半徑R=1600 mm、1750 mm、2000 mm、2500 mm、3400 mm、9000 mm 和1010mm 分別對應(yīng)的扁球殼型結(jié)構(gòu)的屈曲平衡路徑。計(jì)算表明:當(dāng)R≥3400 mm時,撓度與荷載一一對應(yīng),球殼可視為具有一定初曲率的圓板。當(dāng)R<3400 mm 時,屈曲平衡路徑出現(xiàn)了迂回。迂回特性是淺殼從凸面?zhèn)瘸惺軌毫Χ钠胶饴窂降牡湫吞卣鳎祟惽€明確地示出了殼體的上臨界荷載及下臨界荷載(指球殼S 形荷載-撓度曲線上極小值點(diǎn)對應(yīng)的荷載。若規(guī)定凸面?zhèn)葔毫檎?,?dāng)該值大于零時,其指要維持后屈曲狀態(tài)的最小壓力,換句話說,在后屈曲狀態(tài),如果凸面?zhèn)葔毫π∮谠撝?,球殼將會回跳到屈曲前狀態(tài);當(dāng)該值小于零時,表示對于屈曲后的球殼,即使沒有凸面?zhèn)葔毫ψ饔茫驓ひ膊荒茏詣踊靥角盃顟B(tài),該值指使球殼能夠回跳的最小反向壓力)。從圖2 可以看出:隨著R 的減小,球殼上臨界荷載單調(diào)增加,這是扁球殼承載力變化規(guī)律的普遍特征,不必贅述。計(jì)算表明:對于簡支邊功能梯度圓底扁球殼(不考慮變溫),在指數(shù)k、圓底半徑a、厚度h 給定的情況下,隨中曲面曲率半徑R 逐漸減小(球殼深度逐漸增加),其下臨界荷載變化規(guī)律是復(fù)雜的。圖2 所涉及的功能梯度扁球殼,其材料彈性模量從外側(cè)向內(nèi)側(cè)逐漸減小,從圖2 可以看出,隨R 逐漸減小,球殼下臨界荷載略有增加。為了對比,在球殼幾何參數(shù)、梯度指數(shù)以及材料泊松比與圖2 中的相應(yīng)參數(shù)完全相同的情況下,圖3 給出了組分材料為氧化鋯(Ec=151 GPa)和鋼(Em=206 GPa)的球殼(其材料彈性模量從外側(cè)向內(nèi)側(cè)逐漸增加)的屈曲平衡路徑,從圖3 可以看出,隨R 逐漸減小,下臨界荷載是逐漸減小的。上述結(jié)果表明:對于簡支邊功能梯度扁球殼,在梯度指數(shù)、球殼圓底半徑以及厚度給定的情況下,隨球殼深度增加(中曲面曲率半徑R 減小),球殼下臨界荷載的變化規(guī)律與殼體橫向材料彈性模量梯度特性有關(guān)。

    圖2 FGM 扁球殼曲率半徑R 對屈曲平衡路徑的影響Fig.2 Effects of curvature radius R on buckling equilibrium paths of FGM shallow spherical shells

    圖3 FGM 扁球殼曲率半徑R 對屈曲平衡路徑的影響Fig.3 Effects of curvature radius R on buckling equilibrium paths of FGM shallow spherical shells

    在不變溫的情況下,為了揭示材料體積分?jǐn)?shù)指數(shù)k 對FGM 扁球殼屈曲平衡路徑的影響,現(xiàn)在考慮幾何參數(shù)a/R=0.1、R/h=320 的球殼,圖4 給出了多個梯度指數(shù)k 對應(yīng)的球殼的屈曲平衡路徑。從圖4 易見,當(dāng)指數(shù)k 增加時,球殼上臨界荷載是單調(diào)減小的,k 越大,影響程度越小。指數(shù)k 對下臨界荷載的影響規(guī)律比較復(fù)雜,從圖4可以看出,當(dāng)k=0 時,下臨界荷載最小(指代數(shù)量,而不是力的絕對大小,下同);對于圖4 中給定的一些離散的k 值,當(dāng)k 從0 逐漸增加到1 時,下臨界荷載逐漸增加并達(dá)到最大值,在此之后,隨k 增加,下臨界荷載逐漸減小。當(dāng)k 趨于正無窮(球殼已過渡為均勻材料的金屬殼)時,其下臨界荷載仍大于均質(zhì)陶瓷球殼(k=0)的下臨界荷載。

    圖4 體積分?jǐn)?shù)指數(shù)k 對FGM 扁球殼屈曲平衡路徑的影響Fig.4 Effects of volume fraction index k on buckling equilibrium paths of FGM shallow spherical shells

    為了考察組分材料彈性模量對功能梯度扁球殼屈曲平衡路徑的的影響,在其他參數(shù)不變的情況下,將圖4 中的彈性模量為Ec=380 GPa 的陶瓷材料(黏土)用彈性模量為Ec=151 GPa 的陶瓷材料(氧化鋯)替換,并由圖5 給出該功能梯度扁球殼的屈曲平衡路徑。比較圖4 與圖5 不難發(fā)現(xiàn),在金屬組分材料彈性模量不變的情況下,當(dāng)陶瓷組分材料彈性模量顯著減小時,功能梯度扁球殼上臨界荷載會明顯降低。對于圖5 中給定的一些離散的k 值,當(dāng)k=2 時,下臨界荷載達(dá)到最大(其值與k=1 對應(yīng)的球殼下臨界荷載相差很小),其他變化規(guī)律與圖4 相似。需要指出:雖然圖4、圖5 豎坐標(biāo)是有量綱荷載,但只要球殼幾何參數(shù)(R,a,h)滿足圖中比例關(guān)系,其有量綱荷載不會改變,該性質(zhì)對工程設(shè)計(jì)是有意義的。

    為了揭示沿殼體厚度均勻溫度改變對功能梯度扁球殼屈曲平衡路徑的影響,在梯度指數(shù)k=1時,選取球殼幾何參數(shù)a/R=0.1 和R/h=320,圖6給出了均勻變溫T=60 ℃、40 ℃、20 ℃、0 ℃、?20 ℃和?40 ℃6 個變溫參數(shù)時的該功能梯度扁球殼的屈曲平衡路徑。從圖6 易見,溫度變化對球殼屈曲平衡路徑有明顯的影響,溫度升高時,其上臨界荷載明顯增加,下臨界荷載明顯減小。據(jù)此推測,球殼試驗(yàn)環(huán)境溫度相對初始無應(yīng)力狀態(tài)的殼體溫度變化可能是屈曲理論解與試驗(yàn)臨界壓力不一致的一個重要原因。因此,對不涉及變溫的球殼進(jìn)行穩(wěn)定性試驗(yàn),建議在嚴(yán)格恒溫環(huán)境中制作被試驗(yàn)球殼,并在相同環(huán)境溫度下進(jìn)行穩(wěn)定性試驗(yàn)。若要測試變溫對球殼臨界壓力的影響,則是一個比較困難的問題。因?yàn)榍驓さ某o定特性,在球殼溫度相對初始無應(yīng)力狀態(tài)的殼體溫度發(fā)生改變時,如何保證支座的溫度不發(fā)生變化是困難的。另外易見,溫度升高,屈曲過程跳躍幅度增加。

    圖5 體積分?jǐn)?shù)指數(shù)k 對FGM 扁球殼屈曲平衡路徑的影響Fig.5 Effects of volume fraction index k on buckling equilibrium paths of FGM shallow spherical shells

    圖6 變溫場對FGM 扁球殼屈曲平衡路徑的影響Fig.6 Effects of temperature field on buckling equilibrium paths of FGM shallow spherical shells

    4.3 功能梯度扁球殼雙穩(wěn)態(tài)構(gòu)形

    在球殼工作環(huán)境溫度相對初始無應(yīng)力狀態(tài)溫度不發(fā)生改變的情況下,取外側(cè)均布壓力q=1 MPa,圖7 給出了梯度指數(shù)k=0.2 和1.0 對應(yīng)的扁球殼屈曲前和屈曲后的穩(wěn)態(tài)構(gòu)形。為了對比,還給出了扁球殼的初始構(gòu)形(即扁球殼未發(fā)生變形的初始幾何形狀)。對于屈曲前的構(gòu)形,從圖7 可以看出,k 越大,其變形越大。這是因?yàn)殡S指數(shù)k 增加,球殼彎曲剛度減小,進(jìn)而引起殼體變形增加,與預(yù)期結(jié)果一致。對于屈曲后的情況,此類球殼(球殼材料彈性模量從外側(cè)向內(nèi)側(cè)逐漸減小),屈曲后邊界并未充分翻轉(zhuǎn),進(jìn)而使屈曲后構(gòu)形呈現(xiàn)了較復(fù)雜的特性。在外側(cè)均布壓力q=1 MPa 和梯度指數(shù)k=1 時,圖8 給出了三個變溫值T=20 ℃、0 ℃和?20 ℃對應(yīng)扁球殼屈曲前和屈曲后的穩(wěn)態(tài)構(gòu)形。為了對比,還給出了扁球殼的初始構(gòu)形。對于屈曲前的構(gòu)形,從圖8 可以看出,隨溫度升高,球殼變形顯著減小,反之則顯著增大。這是因?yàn)椋呵驓で埃S溫度升高,熱膨脹引起的球殼橫向變形沿凸面外法線方向,與荷載引起的向內(nèi)側(cè)的撓度相抵消;另外,熱膨脹引起球殼深度增加,導(dǎo)致球殼承載力(即上臨界荷載)顯著增加,對應(yīng)于給定荷載的變形(指屈曲前)會減小。還容易看出,溫度升高會導(dǎo)致球殼屈曲過程的跳躍幅度明顯增加。另外,球殼后屈曲構(gòu)形在邊界處的形狀再次顯示了此類球殼(其材料彈性模量從外側(cè)向內(nèi)側(cè)逐漸減小)屈曲后其邊界并不能充分翻轉(zhuǎn)。圖9 給出了材料特性參數(shù)Em=206 GPa,Ec=151 GPa,ν=0.3,幾何參數(shù)R=1600 mm,a=160 mm,h=5 mm的扁球殼在荷載q=0.6 MPa 和未變溫情況下的屈曲前和屈曲后穩(wěn)態(tài)構(gòu)形,同時給出了球殼的初始構(gòu)形。從圖9 可以看出,此類功能梯度球殼(其材料彈性模量從外側(cè)向內(nèi)側(cè)逐漸增加)后屈曲狀態(tài)下邊界發(fā)生了充分的翻轉(zhuǎn)。

    圖7 給定荷載和指數(shù)k 對應(yīng)的FGM 扁球殼的雙穩(wěn)態(tài)構(gòu)形Fig.7 Bi-stable configuration of FGM shallow spherical shells for given load and index k

    圖8 給定荷載和變溫場對應(yīng)的FGM 扁球殼的雙穩(wěn)態(tài)構(gòu)形Fig.8 Bi-stable configuration of FGM shallow spherical shells for given load and temperature field

    圖9 給定荷載對應(yīng)的FGM 扁球殼的雙穩(wěn)態(tài)構(gòu)形Fig.9 Bi-stable configuration of FGM shallow spherical shells for given load

    5 實(shí)用數(shù)表和曲線

    分析本文無量綱兩點(diǎn)邊值問題可知:在球殼工作環(huán)境溫度相對初始無應(yīng)力狀態(tài)的溫度不發(fā)生改變的情況下,對給定幾何參數(shù)a2/(Rh)、指數(shù)k 和組分材料泊松比ν 的扁球殼,其無量綱荷載與球殼中心無量綱撓度的關(guān)系與球殼組分材料彈性模量的絕對值無關(guān),僅取決于其彈性模量比??紤]到工程界更關(guān)注殼體上臨界荷載,即承載力。因此,在未變溫的情況下,針對工程界常見的鋁/黏土(Em/Ec=7/38,ν=0.3)和鋁/氧化鋯(Em/Ec=70/151,ν=0.3)兩種組分材料組合,考慮一組幾何參數(shù)a2/(Rh)和一組梯度指數(shù)k,通過表2 和表3 分別給出了對應(yīng)的扁球殼的無量綱上臨界荷載Qu。對表中的數(shù)據(jù)進(jìn)行換算,能夠方便的得到具有實(shí)際參數(shù)簡支邊功能梯度圓底扁球殼的有量綱上臨界荷載。

    考慮到扁球殼工作環(huán)境溫度會發(fā)生變化,因此,在計(jì)及變溫的情況下,制作類似表2 與表3的數(shù)表能夠方便地指導(dǎo)結(jié)構(gòu)設(shè)計(jì)工作。遺憾的是,計(jì)及變溫時,即使在組分材料常數(shù)給定的情況下,球殼的幾何參數(shù)、材料體積分?jǐn)?shù)指數(shù)和無量綱變溫的組合種類仍然很多,用數(shù)表無法表達(dá)。若用平面曲線來表達(dá),也需要非常多的曲線圖形。限于篇幅,在給定幾何參數(shù)a2/(Rh)=3.2 和a/h=32(球殼有量綱溫度不但與a2/(Rh)的值有關(guān),還與a/h 的值有關(guān))和一組體積分?jǐn)?shù)指數(shù)k(見圖10、圖11)的情況下,針對組分材料組合鋁/黏土和鋁/氧化鋯的簡支邊功能梯度圓底扁球殼,分別由圖10和圖11 給出了球殼的有量綱上臨界荷載qu隨均勻變溫T 的變化曲線,供設(shè)計(jì)人員參考。

    表2 不同幾何參數(shù)a2/(Rh)與材料體積分?jǐn)?shù)指數(shù)k 對應(yīng)的扁球殼無量綱上臨界荷載QuTable2 Dimensionless upper critical load Qu of shallow spherical shells for different geometric parameters a2/(Rh) and volume fraction index k

    表3 不同幾何參數(shù)a2/(Rh)與材料體積分?jǐn)?shù)指數(shù)k 對應(yīng)的扁球殼無量綱上臨界荷載QuTable3 Dimensionless upper critical load Qu of shallow spherical shells for different geometrical parameters a2/(Rh) and volume fraction index k

    圖10 體積分?jǐn)?shù)指數(shù)k 對變溫場中的FGM扁球殼的上臨界荷載的影響Fig.10 Effects of volume fraction index k on upper critical load of FGM shallow spherical shells for different temperature fields

    圖11 體積分?jǐn)?shù)指數(shù)k 對變溫場中的FGM扁球殼的上臨界荷載的影響Fig.11 Effects of volume fraction index k on upper critical load of FGM shallow spherical shells for different temperature fields

    6 結(jié)論

    本文基于經(jīng)典殼理論和von Karman 幾何非線性理論,考慮材料性質(zhì)沿厚度按冪率函數(shù)變化,建立了簡支邊功能梯度材料圓底扁球殼在熱機(jī)械荷載作用下的彈性穩(wěn)定性問題的位移型控制方程及其邊界條件。采用數(shù)值打靶法求得了殼體的平衡路徑和平衡構(gòu)形。分析了簡支邊功能梯度圓底扁球殼在熱機(jī)載荷共同作用下的跳躍屈曲行為。得出如下主要結(jié)論:

    (1) 材料體積分?jǐn)?shù)指數(shù)增加時,殼體上臨界荷載顯著減小,殼體下臨界荷載變化規(guī)律比較復(fù)雜;

    (2) 組分材料彈性模量減小時,殼體上臨界荷載顯著減小,殼體下臨界荷載變化不明顯;

    (3) 溫度均勻升高時,殼體上臨界荷載顯著增加,殼體下臨界荷載顯著減?。?/p>

    (4) 橫向梯度特性對球殼屈曲平衡路徑和后屈曲穩(wěn)態(tài)構(gòu)形有顯著影響。

    猜你喜歡
    球殼構(gòu)形變溫
    空間分布系列球殼聲場耦合與輻射特性研究
    雙星跟飛立體成像的構(gòu)形保持控制
    通有構(gòu)形的特征多項(xiàng)式
    氯乙烯生產(chǎn)中變溫吸附脫水工藝的使用及改進(jìn)
    點(diǎn)電荷和介質(zhì)球殼系統(tǒng)的電勢分布
    焊接殘余應(yīng)力對深潛器耐壓球殼承載能力的影響
    對一個幾何構(gòu)形的探究
    薄壁球殼受剛性平面壓縮的連續(xù)變形測量
    凍融處理對甘薯變溫壓差膨化干燥動力學(xué)的影響
    非共面四頻激光陀螺變溫零偏周期性波動
    大陆偷拍与自拍| 一级毛片我不卡| 在线观看www视频免费| 国产免费又黄又爽又色| 美女脱内裤让男人舔精品视频| 日本午夜av视频| 人妻人人澡人人爽人人| 香蕉精品网在线| 久热这里只有精品99| av.在线天堂| 日韩制服骚丝袜av| 在线观看免费高清a一片| 妹子高潮喷水视频| 久热这里只有精品99| 国产福利在线免费观看视频| 日韩不卡一区二区三区视频在线| 黄片播放在线免费| 97精品久久久久久久久久精品| 丁香六月天网| 99视频精品全部免费 在线| 亚洲熟女精品中文字幕| 久久久久人妻精品一区果冻| 9热在线视频观看99| 国产色爽女视频免费观看| 免费av不卡在线播放| 免费在线观看黄色视频的| 国产成人av激情在线播放| 一级黄片播放器| h视频一区二区三区| 国产免费福利视频在线观看| 各种免费的搞黄视频| 毛片一级片免费看久久久久| 成年av动漫网址| 成年人免费黄色播放视频| 久久99蜜桃精品久久| 伦理电影免费视频| 在线免费观看不下载黄p国产| 青春草视频在线免费观看| 中文字幕人妻丝袜制服| 99久久中文字幕三级久久日本| 亚洲婷婷狠狠爱综合网| 久久精品国产a三级三级三级| 看免费成人av毛片| 欧美人与性动交α欧美精品济南到 | 美女主播在线视频| 久久99热这里只频精品6学生| 建设人人有责人人尽责人人享有的| 赤兔流量卡办理| 22中文网久久字幕| 国产激情久久老熟女| 在线 av 中文字幕| 精品视频人人做人人爽| 天天操日日干夜夜撸| 精品午夜福利在线看| 高清欧美精品videossex| 国产日韩欧美亚洲二区| 熟女电影av网| 各种免费的搞黄视频| 高清欧美精品videossex| tube8黄色片| 国产av码专区亚洲av| 亚洲国产精品专区欧美| 久久精品国产综合久久久 | 国产亚洲欧美精品永久| 美国免费a级毛片| 亚洲图色成人| av视频免费观看在线观看| 少妇的逼水好多| 男女国产视频网站| 成人无遮挡网站| 亚洲成人一二三区av| 日本欧美视频一区| 大香蕉久久网| 成人毛片a级毛片在线播放| 成人国产麻豆网| 成人影院久久| 老司机亚洲免费影院| 国产精品一二三区在线看| 久久久a久久爽久久v久久| 午夜视频国产福利| 亚洲在久久综合| 成人亚洲精品一区在线观看| 菩萨蛮人人尽说江南好唐韦庄| 最后的刺客免费高清国语| 极品少妇高潮喷水抽搐| 成年美女黄网站色视频大全免费| 国产国拍精品亚洲av在线观看| 成年人午夜在线观看视频| 成年动漫av网址| 一级,二级,三级黄色视频| av一本久久久久| 国产精品国产三级国产专区5o| 尾随美女入室| 又大又黄又爽视频免费| 丰满少妇做爰视频| 国产精品麻豆人妻色哟哟久久| 满18在线观看网站| 黑人猛操日本美女一级片| 亚洲精品第二区| 中国美白少妇内射xxxbb| 在线观看免费日韩欧美大片| 免费观看无遮挡的男女| 久久精品国产a三级三级三级| 国产av码专区亚洲av| 国产精品人妻久久久影院| 亚洲国产精品成人久久小说| 另类亚洲欧美激情| av免费观看日本| 久久韩国三级中文字幕| tube8黄色片| 欧美精品人与动牲交sv欧美| 亚洲情色 制服丝袜| 韩国av在线不卡| 亚洲在久久综合| 亚洲欧洲国产日韩| 久久国内精品自在自线图片| 久久精品aⅴ一区二区三区四区 | 亚洲激情五月婷婷啪啪| 在线观看美女被高潮喷水网站| 香蕉精品网在线| 色94色欧美一区二区| 成人国产av品久久久| 久久ye,这里只有精品| 啦啦啦啦在线视频资源| 精品酒店卫生间| 99热网站在线观看| a级毛片黄视频| 亚洲欧美成人精品一区二区| 精品国产乱码久久久久久小说| 亚洲高清免费不卡视频| 精品国产一区二区三区四区第35| 婷婷色麻豆天堂久久| 久久精品aⅴ一区二区三区四区 | 汤姆久久久久久久影院中文字幕| 看免费av毛片| 国产亚洲av片在线观看秒播厂| 男女免费视频国产| 久久午夜综合久久蜜桃| 免费黄网站久久成人精品| 国产有黄有色有爽视频| 久久久精品94久久精品| 999精品在线视频| 亚洲av电影在线进入| 亚洲精品色激情综合| 丝袜脚勾引网站| 亚洲欧美成人精品一区二区| 日本wwww免费看| 我要看黄色一级片免费的| 久久精品aⅴ一区二区三区四区 | 亚洲国产欧美在线一区| 少妇熟女欧美另类| 少妇高潮的动态图| xxx大片免费视频| 秋霞在线观看毛片| 久久久国产精品麻豆| 亚洲精品乱久久久久久| 亚洲欧美日韩另类电影网站| 欧美日本中文国产一区发布| 最近最新中文字幕免费大全7| 久久久久精品性色| 久久精品国产亚洲av涩爱| 国产精品秋霞免费鲁丝片| 在线 av 中文字幕| 99国产精品免费福利视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品成人久久小说| 看免费成人av毛片| 亚洲av福利一区| 人妻系列 视频| 成人亚洲精品一区在线观看| 国国产精品蜜臀av免费| 18禁观看日本| 亚洲精品aⅴ在线观看| 国产精品女同一区二区软件| 丝袜喷水一区| 国产黄频视频在线观看| 日韩欧美精品免费久久| 80岁老熟妇乱子伦牲交| 26uuu在线亚洲综合色| 春色校园在线视频观看| 欧美老熟妇乱子伦牲交| www.色视频.com| 国产精品一区www在线观看| 欧美+日韩+精品| 亚洲美女视频黄频| a 毛片基地| 99re6热这里在线精品视频| 亚洲综合色惰| 麻豆精品久久久久久蜜桃| 国产成人精品在线电影| 男女国产视频网站| 伊人久久国产一区二区| 18在线观看网站| 大片免费播放器 马上看| 日韩精品免费视频一区二区三区 | 欧美xxxx性猛交bbbb| 少妇高潮的动态图| 国产精品欧美亚洲77777| 中文字幕亚洲精品专区| 母亲3免费完整高清在线观看 | 在现免费观看毛片| 欧美激情国产日韩精品一区| 亚洲成国产人片在线观看| 韩国av在线不卡| 成人亚洲欧美一区二区av| av一本久久久久| 亚洲精品一区蜜桃| 视频在线观看一区二区三区| 精品久久久精品久久久| 熟女av电影| 精品久久久久久电影网| 一边摸一边做爽爽视频免费| 国产男女超爽视频在线观看| 丰满迷人的少妇在线观看| 极品少妇高潮喷水抽搐| 街头女战士在线观看网站| 亚洲国产最新在线播放| 宅男免费午夜| 国产男人的电影天堂91| 国产精品蜜桃在线观看| 一级片免费观看大全| 卡戴珊不雅视频在线播放| 美女中出高潮动态图| 欧美bdsm另类| 国产黄色视频一区二区在线观看| 成人亚洲欧美一区二区av| av天堂久久9| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 人妻少妇偷人精品九色| 国产免费现黄频在线看| 久久午夜综合久久蜜桃| 在线观看一区二区三区激情| 亚洲色图 男人天堂 中文字幕 | 国产精品一区二区在线观看99| 丰满乱子伦码专区| 免费黄频网站在线观看国产| 久久精品国产自在天天线| 黄色视频在线播放观看不卡| 99视频精品全部免费 在线| 丰满迷人的少妇在线观看| tube8黄色片| 观看av在线不卡| 人妻少妇偷人精品九色| 人体艺术视频欧美日本| 91在线精品国自产拍蜜月| 精品国产一区二区三区四区第35| 国产一区二区在线观看av| 国产精品一二三区在线看| 午夜福利在线观看免费完整高清在| 欧美性感艳星| 午夜福利乱码中文字幕| 十分钟在线观看高清视频www| 蜜桃国产av成人99| 国产精品 国内视频| 成人无遮挡网站| 欧美+日韩+精品| 日韩av免费高清视频| 国产深夜福利视频在线观看| 亚洲熟女精品中文字幕| 狂野欧美激情性bbbbbb| 精品熟女少妇av免费看| 午夜福利视频精品| 99国产精品免费福利视频| 免费久久久久久久精品成人欧美视频 | 亚洲在久久综合| 国产精品免费大片| 免费黄频网站在线观看国产| 亚洲av成人精品一二三区| 麻豆精品久久久久久蜜桃| 日本黄大片高清| 嫩草影院入口| 成年人免费黄色播放视频| www.熟女人妻精品国产 | 国产视频首页在线观看| 国产xxxxx性猛交| 国产麻豆69| 国产亚洲一区二区精品| 日韩av不卡免费在线播放| 一本久久精品| 久久精品人人爽人人爽视色| 观看美女的网站| 国产av一区二区精品久久| 18禁裸乳无遮挡动漫免费视频| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 老司机亚洲免费影院| 午夜免费观看性视频| 亚洲av男天堂| 亚洲精品国产av蜜桃| 亚洲精品,欧美精品| 国产精品久久久久久精品电影小说| videosex国产| 啦啦啦视频在线资源免费观看| 日韩在线高清观看一区二区三区| 精品久久国产蜜桃| 巨乳人妻的诱惑在线观看| 国产一区二区在线观看av| 亚洲国产欧美日韩在线播放| av播播在线观看一区| 国产精品一区www在线观看| 你懂的网址亚洲精品在线观看| 国产1区2区3区精品| 国产精品不卡视频一区二区| www.熟女人妻精品国产 | 激情五月婷婷亚洲| 久久婷婷青草| 亚洲综合色网址| 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 丝袜在线中文字幕| 国产黄频视频在线观看| 秋霞在线观看毛片| 男女边摸边吃奶| 精品一区二区三区视频在线| 制服丝袜香蕉在线| 高清欧美精品videossex| 国产精品一区二区在线观看99| 最近中文字幕高清免费大全6| 欧美老熟妇乱子伦牲交| 久久ye,这里只有精品| 久久精品熟女亚洲av麻豆精品| 七月丁香在线播放| 乱人伦中国视频| av免费观看日本| 久久精品久久久久久久性| 午夜日本视频在线| 免费黄色在线免费观看| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| av.在线天堂| 日韩欧美精品免费久久| 成人亚洲欧美一区二区av| 久久国产亚洲av麻豆专区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久影院123| h视频一区二区三区| 欧美精品国产亚洲| 国产精品久久久av美女十八| 精品午夜福利在线看| av在线老鸭窝| 色吧在线观看| 国产一区有黄有色的免费视频| 日本wwww免费看| 久久久久久人妻| 美女主播在线视频| 中文字幕制服av| 日本色播在线视频| 国产精品久久久久久久久免| 五月伊人婷婷丁香| 亚洲国产成人一精品久久久| 99热国产这里只有精品6| 国产精品久久久久久久久免| 欧美bdsm另类| 丝瓜视频免费看黄片| 免费播放大片免费观看视频在线观看| 欧美成人午夜免费资源| 精品国产一区二区久久| 高清在线视频一区二区三区| 十分钟在线观看高清视频www| 极品少妇高潮喷水抽搐| 成人国产av品久久久| 一二三四在线观看免费中文在 | 一本大道久久a久久精品| 最近中文字幕2019免费版| av国产久精品久网站免费入址| videosex国产| 好男人视频免费观看在线| 国产免费现黄频在线看| 欧美日韩国产mv在线观看视频| 成人亚洲欧美一区二区av| 好男人视频免费观看在线| 国产精品一区二区在线观看99| 亚洲成人手机| 一级毛片黄色毛片免费观看视频| 日本欧美国产在线视频| 纵有疾风起免费观看全集完整版| 久久久久国产精品人妻一区二区| 超色免费av| 国产精品一二三区在线看| 国产黄频视频在线观看| 国产精品久久久久久久久免| 26uuu在线亚洲综合色| 精品99又大又爽又粗少妇毛片| 日韩伦理黄色片| 久久久久久久大尺度免费视频| 欧美日韩视频精品一区| 国产爽快片一区二区三区| 日韩不卡一区二区三区视频在线| 涩涩av久久男人的天堂| xxx大片免费视频| 精品视频人人做人人爽| 国产亚洲最大av| 宅男免费午夜| 国产av国产精品国产| 亚洲中文av在线| 国产精品久久久久久av不卡| 亚洲成人av在线免费| 国产在视频线精品| 中文天堂在线官网| 三级国产精品片| 中文欧美无线码| 亚洲成人一二三区av| 欧美日韩精品成人综合77777| 久久久久久久国产电影| 日韩制服骚丝袜av| 免费少妇av软件| 国产男女内射视频| 尾随美女入室| 91精品国产国语对白视频| 成人免费观看视频高清| 午夜福利影视在线免费观看| 欧美精品一区二区大全| 99热6这里只有精品| 黄片无遮挡物在线观看| 亚洲国产成人一精品久久久| 桃花免费在线播放| 免费人成在线观看视频色| 国国产精品蜜臀av免费| 欧美激情国产日韩精品一区| 99久久综合免费| 日本免费在线观看一区| 亚洲欧美一区二区三区国产| 亚洲色图 男人天堂 中文字幕 | 亚洲精品,欧美精品| 国产一级毛片在线| 精品久久久久久电影网| 美女中出高潮动态图| 亚洲欧洲精品一区二区精品久久久 | 在线观看免费高清a一片| 成人亚洲欧美一区二区av| 亚洲四区av| 22中文网久久字幕| 999精品在线视频| 免费播放大片免费观看视频在线观看| 国产高清三级在线| 国产精品欧美亚洲77777| 色5月婷婷丁香| 久久青草综合色| 亚洲人成77777在线视频| 99久久精品国产国产毛片| 久久久久久久久久久久大奶| 久久99热6这里只有精品| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 夫妻午夜视频| 日本黄色日本黄色录像| 午夜老司机福利剧场| 亚洲第一区二区三区不卡| 国产精品一区www在线观看| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| 最近最新中文字幕免费大全7| 久久精品熟女亚洲av麻豆精品| 久久午夜综合久久蜜桃| 国产白丝娇喘喷水9色精品| 国产精品久久久久成人av| 欧美最新免费一区二区三区| 久久 成人 亚洲| 美女福利国产在线| 国产日韩欧美亚洲二区| 香蕉国产在线看| 97精品久久久久久久久久精品| 最近最新中文字幕大全免费视频 | 少妇人妻精品综合一区二区| 午夜激情久久久久久久| 亚洲国产精品一区二区三区在线| 国产免费一区二区三区四区乱码| 亚洲精品美女久久av网站| 午夜激情久久久久久久| 一个人免费看片子| 精品一区二区免费观看| 丰满饥渴人妻一区二区三| 欧美xxⅹ黑人| 又黄又粗又硬又大视频| 国产av码专区亚洲av| 国语对白做爰xxxⅹ性视频网站| 午夜福利网站1000一区二区三区| 超碰97精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久久狼人影院| 久久99热这里只频精品6学生| 精品国产一区二区三区久久久樱花| 欧美日韩av久久| 波多野结衣一区麻豆| 97在线人人人人妻| 伦理电影免费视频| 久久久国产精品麻豆| 人妻人人澡人人爽人人| 亚洲精品日韩在线中文字幕| 高清毛片免费看| 亚洲av福利一区| 国产成人一区二区在线| 丰满少妇做爰视频| 国产黄色免费在线视频| 99久久人妻综合| 一级a做视频免费观看| 国产一区二区三区综合在线观看 | 欧美精品一区二区大全| 999精品在线视频| 中国三级夫妇交换| 免费看不卡的av| 在线观看免费高清a一片| 国产福利在线免费观看视频| 丝袜在线中文字幕| 亚洲国产欧美在线一区| 精品人妻熟女毛片av久久网站| 婷婷色综合www| 欧美亚洲 丝袜 人妻 在线| 老熟女久久久| 激情五月婷婷亚洲| 免费高清在线观看日韩| 97人妻天天添夜夜摸| 男人操女人黄网站| 精品亚洲成国产av| 女人精品久久久久毛片| 国产精品人妻久久久久久| 久久韩国三级中文字幕| 亚洲av欧美aⅴ国产| 亚洲国产成人一精品久久久| 一边亲一边摸免费视频| 国产成人精品福利久久| 亚洲国产精品成人久久小说| 久久精品国产自在天天线| 亚洲国产看品久久| 日日爽夜夜爽网站| 美女大奶头黄色视频| 中文字幕制服av| 日韩av免费高清视频| 免费大片黄手机在线观看| 久久久久网色| 少妇被粗大猛烈的视频| 日韩一区二区视频免费看| 久久97久久精品| 人成视频在线观看免费观看| 欧美成人精品欧美一级黄| 午夜激情久久久久久久| 岛国毛片在线播放| 水蜜桃什么品种好| 精品一区二区三区四区五区乱码 | 日韩大片免费观看网站| 欧美老熟妇乱子伦牲交| 婷婷色综合www| 91aial.com中文字幕在线观看| 午夜福利视频在线观看免费| 久久午夜综合久久蜜桃| 美女脱内裤让男人舔精品视频| 亚洲精品第二区| 久久精品久久久久久噜噜老黄| 午夜福利乱码中文字幕| 大片免费播放器 马上看| 丝袜人妻中文字幕| 97在线人人人人妻| 久久综合国产亚洲精品| 肉色欧美久久久久久久蜜桃| 美女视频免费永久观看网站| 久久人人爽人人爽人人片va| 伦理电影大哥的女人| 51国产日韩欧美| 99热网站在线观看| 欧美+日韩+精品| 在现免费观看毛片| 久久久久精品久久久久真实原创| 成人午夜精彩视频在线观看| 毛片一级片免费看久久久久| 亚洲精品久久久久久婷婷小说| 久久久国产欧美日韩av| 亚洲精品久久午夜乱码| 亚洲av欧美aⅴ国产| 亚洲图色成人| 欧美性感艳星| 日本黄大片高清| 狠狠精品人妻久久久久久综合| 欧美少妇被猛烈插入视频| 精品国产一区二区三区久久久樱花| 99热6这里只有精品| 最近最新中文字幕免费大全7| 99re6热这里在线精品视频| 国产片特级美女逼逼视频| 99香蕉大伊视频| 一区二区日韩欧美中文字幕 | 在线精品无人区一区二区三| 精品一区二区三区视频在线| 国产精品秋霞免费鲁丝片| a级片在线免费高清观看视频| 大香蕉97超碰在线| av福利片在线| 免费看不卡的av| 国产亚洲精品久久久com| 看非洲黑人一级黄片| 国产精品国产av在线观看| 99九九在线精品视频| 制服丝袜香蕉在线| 内地一区二区视频在线| 亚洲国产精品一区三区| 国产免费一级a男人的天堂| 亚洲美女搞黄在线观看| 亚洲国产精品一区三区| 国产成人精品无人区| 国产在线一区二区三区精| 人人澡人人妻人| av网站免费在线观看视频| 亚洲伊人色综图| 精品少妇久久久久久888优播| 国产永久视频网站| 在线观看人妻少妇| 高清不卡的av网站| 欧美亚洲 丝袜 人妻 在线| 夜夜爽夜夜爽视频| 黄色毛片三级朝国网站| 欧美日韩成人在线一区二区| 美女中出高潮动态图| 午夜91福利影院| 五月伊人婷婷丁香| 日韩在线高清观看一区二区三区|