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

    基于ANSYS優(yōu)化技術(shù)的非對(duì)稱管板的分析與設(shè)計(jì)

    2017-07-03 14:58:16孫偉明石秀真李劍虹
    關(guān)鍵詞:管區(qū)管板非對(duì)稱

    孫偉明,石秀真,周 文,常 波,李劍虹

    (1.浙江工業(yè)大學(xué) 化工機(jī)械設(shè)計(jì)研究所,浙江 杭州 310014;2.衢州市特種設(shè)備檢驗(yàn)中心,浙江 衢州 324000)

    基于ANSYS優(yōu)化技術(shù)的非對(duì)稱管板的分析與設(shè)計(jì)

    孫偉明1,石秀真1,周 文2,常 波1,李劍虹1

    (1.浙江工業(yè)大學(xué) 化工機(jī)械設(shè)計(jì)研究所,浙江 杭州 310014;2.衢州市特種設(shè)備檢驗(yàn)中心,浙江 衢州 324000)

    非對(duì)稱管殼式蒸發(fā)器是低溫余熱發(fā)電的核心設(shè)備,在近年來應(yīng)用越來越廣泛.非對(duì)稱管板即管板的管束布置呈非對(duì)稱狀,而管束的特殊分布使得管板的應(yīng)力情況復(fù)雜,目前沒有非對(duì)稱管板板厚的具體計(jì)算公式,只能借助數(shù)值求解的方法來計(jì)算.從安全生產(chǎn)的角度出發(fā),為了在滿足結(jié)構(gòu)強(qiáng)度的要求下達(dá)到成本最小的目的,利用ANSYA有限元軟件對(duì)構(gòu)件進(jìn)行了應(yīng)力分析,通過ANSYS優(yōu)化技術(shù),選取管板厚度為設(shè)計(jì)變量,蒸發(fā)器質(zhì)量為目標(biāo)函數(shù),對(duì)管板參數(shù)進(jìn)行優(yōu)化,最終得出了合理的結(jié)果.

    非對(duì)稱管板;有限元;優(yōu)化設(shè)計(jì)

    膨脹機(jī)蒸發(fā)器是一種有特殊用途的管殼式換熱器,應(yīng)用于石油相關(guān)的氣體加工、能量回收、蒸汽與各種氣體的減壓等場(chǎng)合[1].膨脹機(jī)蒸發(fā)器在設(shè)計(jì)時(shí)為給殼程一個(gè)氣相膨脹空間,設(shè)計(jì)成上下非對(duì)稱管板,換熱管布置在管板下半部分.由于蒸發(fā)器管板的形狀和受力復(fù)雜,所以在結(jié)構(gòu)設(shè)計(jì)時(shí)不能依據(jù)常規(guī)設(shè)計(jì)標(biāo)準(zhǔn)熱交換器GB/T 151—2014[2].隨著有限元分析技術(shù)的出現(xiàn),人們?cè)絹碓蕉嗟貙⑦@種技術(shù)應(yīng)用于非對(duì)稱管殼式換熱器的管板結(jié)構(gòu)設(shè)計(jì).有限元設(shè)計(jì)方法能滿足結(jié)構(gòu)強(qiáng)度的要求,但不一定是結(jié)構(gòu)設(shè)計(jì)中的最優(yōu)方法[3-10].有限元計(jì)算和優(yōu)化方案相結(jié)合,是尋找最優(yōu)方案的一種有效方法.利用ANSYS提供的參數(shù)優(yōu)化功能,提出了非對(duì)稱式管殼式換熱器管板的優(yōu)化設(shè)計(jì)的方法.

    1 蒸發(fā)器結(jié)構(gòu)和參數(shù)

    蒸發(fā)器為固定管板式換熱器,管板下半部分的換熱管采用對(duì)稱式布置,共有換熱管1 850根,管長(zhǎng)9 000 mm,換熱管與管板之間采用脹接方式連接,其管板結(jié)構(gòu)如圖1所示.

    這里重點(diǎn)研究的是蒸發(fā)器的管板,管板初始厚度為79 mm,假設(shè)為平蓋厚度按照壓力容器GB 150—2011[11]確定.管板的主要性能參數(shù)如表1所示.

    表1 管板主要參數(shù)

    圖1 管板結(jié)構(gòu)及布管方式Fig.1 The tube sheet and piping manner

    2 有限元分析

    2.1 有限元模型

    如圖1所示,蒸發(fā)器為左右對(duì)稱結(jié)構(gòu),為了減少計(jì)算的量,建模時(shí)取其1/2對(duì)稱模型.由于本次計(jì)算主要是針對(duì)管板的厚度,在計(jì)算過程中只考慮與管板有關(guān)的法蘭、殼體和換熱管而忽略了殼體內(nèi)部的其他非承壓部件.模型采用solid186單元對(duì)模型進(jìn)行網(wǎng)格劃分.模型中,共有384 780個(gè)單元,609 524個(gè)節(jié)點(diǎn),其模型如圖2所示.

    圖2 有限元模型圖Fig.2 Finite element model

    2.2 約束和載荷

    管板在工作中承受管程壓力、殼程壓力、墊片壓緊力和螺栓預(yù)緊力的作用.其中管程壓力和殼程壓力為已知條件,墊片壓緊力和螺栓預(yù)緊力根據(jù)GB 150—2011進(jìn)行計(jì)算,計(jì)算公式為

    Fp=6.28DGbmpc

    (1)

    (2)

    式中:Fp為墊片壓緊力,N;DG為墊片壓緊力作用中心圓直徑,mm;b為墊片有效密封寬度,mm;m為墊片系數(shù);pc為計(jì)算壓力,MPa;Wp為螺栓預(yù)緊力,N.蒸發(fā)器中共有螺栓56個(gè),通過計(jì)算得到墊片壓緊力為4.92N,單個(gè)螺栓的預(yù)緊力為13 338N.

    從工作參數(shù)可以看到蒸發(fā)器的殼程設(shè)計(jì)溫度為144 ℃,管程設(shè)計(jì)溫度為150 ℃,由于法蘭、管板、殼體和換熱管的溫度在144~150 ℃范圍內(nèi),溫差較小,可以忽略不計(jì),本次有限元分析過程中未考慮溫度場(chǎng)對(duì)管板的影響.

    2.3 結(jié)果分析

    根據(jù)設(shè)計(jì)條件,對(duì)蒸發(fā)器管板在管程與殼程同時(shí)承受壓力、管程單獨(dú)受壓以及殼程單獨(dú)受壓等三種情況下進(jìn)行應(yīng)力計(jì)算分析,得到管板最大應(yīng)力位置處應(yīng)力情況如表2所示.由表2可知:當(dāng)殼程單獨(dú)受壓時(shí),產(chǎn)生的Stressintensity應(yīng)力強(qiáng)度、一次薄膜應(yīng)力值和彎曲應(yīng)力值等皆大于管程單獨(dú)受壓、管程與殼程同時(shí)承受壓力時(shí)的應(yīng)力值,所以,以下優(yōu)化計(jì)算中以殼程單獨(dú)受壓為準(zhǔn).

    表2 管板應(yīng)力

    從上述計(jì)算結(jié)果可以看出:三種加載情況下最大應(yīng)力值都出現(xiàn)在管板與筒體連接的焊縫上,如圖3(a)所示.根據(jù)平板的受力情況可知,筒體內(nèi)應(yīng)力均勻分布,最大應(yīng)力出現(xiàn)在管板與筒體連接處.而均勻布管管板,由管板彎矩分布及彎曲應(yīng)力分布情況得出,管板最大徑向彎矩及彎曲應(yīng)力出現(xiàn)在布管區(qū)邊界位置.而非對(duì)稱管板布管在下半部分,管板上半部分為平板,管板承受內(nèi)壓,布管區(qū)有換熱管支撐,所以最大應(yīng)力發(fā)生在上半部分,又焊縫處應(yīng)力集中,所以非對(duì)稱管板應(yīng)力分布介于平板和均勻布管管板之間.

    1—焊縫/最大應(yīng)力位置;2—管板中心;3—邊緣;4—布管邊界圖3 管板應(yīng)力分布云圖Fig.3 Stress contour plot of tube sheet

    管板的安全分析不只取決于最大應(yīng)力值發(fā)生位置,由JB/T 4732—2005[12]可知:平板的危險(xiǎn)位置有板邊緣處、板中心處以及板與筒體連接處,而均勻布管管板需要進(jìn)行應(yīng)力校核的地方有管板中心、管板布管區(qū)與非布管區(qū)邊界處以及管板邊緣處.因?yàn)榉菍?duì)稱管板的特殊性,結(jié)合平板與管板,對(duì)其板中心、布管區(qū)與非布管區(qū)邊界處、管板與筒體連接的焊縫上以及邊緣處等進(jìn)行應(yīng)力分析.表3為對(duì)各處的應(yīng)力進(jìn)行分類得出的結(jié)果.

    表3 各危險(xiǎn)界面應(yīng)力分類

    2.4 對(duì)比分析

    目前,對(duì)于非對(duì)稱管板的厚度計(jì)算以及應(yīng)力校核都沒有具體公式,所以針對(duì)這類管板的設(shè)計(jì)與應(yīng)力校核都存在一定困難.由上文可知:此類管板的應(yīng)力分布介于平板與均勻布管管板,在此根據(jù)JB/T 4732—2005計(jì)算平板邊緣處、板中心處以及板與筒體連接處的最大應(yīng)力強(qiáng)度,根據(jù)GB/T 151—2014計(jì)算均勻布管管板中心、管板布管區(qū)與非布管區(qū)邊界處以及管板邊緣處的最大應(yīng)力強(qiáng)度.

    表4為通過JB/T 4732—2005,GB/T 151—2014中的公式對(duì)平板及均勻布管管板的邊緣、中心、板與筒體連接處以及布管區(qū)與非布管區(qū)邊界處等危險(xiǎn)界面進(jìn)行應(yīng)力計(jì)算得出的結(jié)果,與表3中非對(duì)稱管板的應(yīng)力強(qiáng)度進(jìn)行對(duì)比.

    由上述對(duì)比結(jié)果可知:非對(duì)稱管板的最大應(yīng)力粗略校核可參考相同工況下的平板與管板.具體危險(xiǎn)界面的參考如下:管板中心及布管邊界可參考相同工況下的均勻布管管板,焊縫處及管板邊緣可參考相同工況下的平板.

    表4 平板與管板各危險(xiǎn)界面的應(yīng)力

    2.5 基于ANSYS進(jìn)行優(yōu)化

    由上述分析可知:最危險(xiǎn)的位置位于管板與筒體連接的焊縫上,以下優(yōu)化部分則針對(duì)此位置進(jìn)行.

    管板與筒體的焊縫上,存在一次應(yīng)力和二次應(yīng)力.JB/T 4732—2005中有關(guān)應(yīng)力強(qiáng)度的許用極限的規(guī)定:一次薄膜應(yīng)力加一次彎曲應(yīng)力小于等于1.5倍許用應(yīng)力,一次應(yīng)力加二次應(yīng)力小于等于3倍許用應(yīng)力.而GB/T 151—2014中規(guī)定,僅殼程壓力作用,不計(jì)膨脹差時(shí),最大應(yīng)力強(qiáng)度應(yīng)小于1.5倍許用應(yīng)力;計(jì)入膨脹差時(shí),最大應(yīng)力強(qiáng)度應(yīng)小于3倍許用應(yīng)力.本次計(jì)算中,忽略熱膨脹,綜合考慮JB/T 4732—2005和GB/T 151—2014,選取最大應(yīng)力強(qiáng)度σmax<1.5[σ],一次薄膜應(yīng)力與彎曲應(yīng)力的和σn<1.5[σ].由計(jì)算結(jié)果可知σmax=125.225 MPa.在最大應(yīng)力發(fā)生位置進(jìn)行應(yīng)力分類,σn=157.8 MPa.σn,σmax都遠(yuǎn)小于材料的1.5[σ]=250 MPa,雖然滿足了安全要求,但是材料利用率太小,導(dǎo)致經(jīng)濟(jì)成本升高.

    優(yōu)化設(shè)計(jì)是一種尋找最優(yōu)設(shè)計(jì)方案的設(shè)計(jì)方法,而最優(yōu)設(shè)計(jì)方案是既能夠滿足所有設(shè)計(jì)要求的方案,又能保證支出(如體積、重量和費(fèi)用等)最小,是最有效率的方案.ANSYS有限元軟件提供了適用于優(yōu)化分析的一階優(yōu)化方法[13].

    由于優(yōu)化計(jì)算是針對(duì)管板強(qiáng)度的安全設(shè)計(jì),所以在保證安全的情況下,是以確定用材最少的設(shè)計(jì)方案為計(jì)算目的的.因此,除板厚之外其他尺寸不得改變.為使結(jié)構(gòu)能滿足設(shè)計(jì)要求,選擇Stress intensity應(yīng)力強(qiáng)度為狀態(tài)變量,確定管板厚度為設(shè)計(jì)變量.因?yàn)椴牧厦芏葹槌?shù),尋找滿足安全要求的最小管板厚度,同時(shí)也是構(gòu)件的最小質(zhì)量,因此可以選擇構(gòu)件的質(zhì)量作為目標(biāo)函數(shù).狀態(tài)變量為,假設(shè)管板厚度為T,變化范圍為20~100 mm,那么該問題的數(shù)學(xué)模型[14]可表示為

    設(shè)計(jì)變量:T

    目標(biāo)函數(shù):minWT(x)

    當(dāng)20 mm

    按照優(yōu)化過程的步驟,以原有計(jì)算模型為基礎(chǔ),根據(jù)上述優(yōu)化數(shù)學(xué)模型,采用參數(shù)化設(shè)計(jì)語(yǔ)言[15],由ANSYS軟件的一階優(yōu)化方法進(jìn)行運(yùn)算,在第8個(gè)循環(huán)步結(jié)束時(shí)得到了最優(yōu)結(jié)果,計(jì)算終止.表5給出了每個(gè)循環(huán)步的計(jì)算結(jié)果,從表5中的數(shù)據(jù)可以看出:蒸發(fā)器非對(duì)稱管板經(jīng)ANSYS軟件優(yōu)化,最大應(yīng)力強(qiáng)度已經(jīng)從125.225 MPa升高到了250.99 MPa,管板厚度變?yōu)?3.594 mm,對(duì)比初始值79 mm降低了43%,構(gòu)件質(zhì)量從2 630.8 kg降低到了2 325.5 kg,且滿足安全要求,大大提高了材料的利用率,降低了產(chǎn)品成本.因此最終優(yōu)化結(jié)果應(yīng)選擇表5中第8步的計(jì)算結(jié)果.

    表5 優(yōu)化數(shù)據(jù)的變化情況

    注:1) 第8步為最優(yōu)設(shè)計(jì).

    對(duì)設(shè)計(jì)變量?jī)?yōu)化結(jié)果T進(jìn)行圓整,為方便加工,取管板厚度T=45 mm.修改計(jì)算模型,重新計(jì)算,得到管板應(yīng)力云圖,其中最大應(yīng)力強(qiáng)度σmax=249.30 MPa,一次薄膜應(yīng)力加一次彎曲應(yīng)力σn=204.43 MPa,構(gòu)件應(yīng)力強(qiáng)度云圖如圖3(b)所示.由圖3(b)看出:管板厚度發(fā)生變化后,最大應(yīng)力位置沒有改變,且對(duì)比其他危險(xiǎn)截面,焊縫處應(yīng)力明顯大于其他危險(xiǎn)截面,故在對(duì)管板進(jìn)行優(yōu)化設(shè)計(jì)時(shí),僅考慮焊縫處.

    根據(jù)GB/T 151—2014規(guī)定,對(duì)帶法蘭的管板的厚度進(jìn)行設(shè)計(jì)時(shí),不僅要保證管板的強(qiáng)度,也要保證管板的剛度以及換熱管的軸向應(yīng)力σt,根據(jù)GB 150—2011,帶法蘭管板剛度指數(shù)J及換熱管軸向應(yīng)力的計(jì)算式分別為

    (3)

    (4)

    式中:VI為整體法蘭系數(shù);M0為法蘭設(shè)計(jì)力矩,N·mm;λ為系數(shù);δ0位法蘭有效厚度,mm;KI為剛度系數(shù),取0.3;h0為參數(shù);a為1根換熱管管壁金屬的橫截面積,mm2;d為換熱管外徑,mm;σt為換熱管軸向應(yīng)力,MPa;[σ]t為設(shè)計(jì)溫度下?lián)Q熱管材料的許用應(yīng)力,MPa.換熱管選用材料為20號(hào)鋼,外徑為19 mm,壁厚為2 mm.

    通過式(3,4)計(jì)算得出管板剛度指數(shù)J=6.72×10-5<1,軸向應(yīng)力|σt|=2.92 MPa<[σ]t=147 MPa.計(jì)算結(jié)果滿足標(biāo)準(zhǔn)要求,可以保證管板優(yōu)化后的結(jié)構(gòu)安全.

    3 結(jié) 論

    借助ANSYS有限元分析軟件對(duì)蒸發(fā)器非對(duì)稱管板進(jìn)行應(yīng)力校核,并提供了一種簡(jiǎn)略的校核方法.利用有限元優(yōu)化模塊,在滿足強(qiáng)度要求的前提下建立數(shù)學(xué)模型,對(duì)管板的結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化,可得出最優(yōu)設(shè)計(jì)方案.利用ANSYS有限元優(yōu)化技術(shù)求解結(jié)構(gòu)優(yōu)化問題,最突出的優(yōu)點(diǎn)就是可避免煩瑣的建模、計(jì)算工作,從而可以提高設(shè)計(jì)效率,同時(shí)也滿足了機(jī)械設(shè)計(jì)中安全性、經(jīng)濟(jì)性的需要.

    [1] 陸征,蔣國(guó).螺桿膨脹機(jī)的發(fā)展及應(yīng)用[C]//2010年設(shè)備、機(jī)泵、攪拌年會(huì)會(huì)刊.大連:全國(guó)化工設(shè)備設(shè)計(jì)技術(shù)中心,2010.

    [2] 中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).熱交換器: GB/T 151—2014[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2015:2.

    [3] 徐學(xué)真,黃建冰,孫偉明.帶金屬O型環(huán)法蘭的密封性能有限元分析[J].浙江工業(yè)大學(xué)學(xué)報(bào),2014,42(6):686-689.

    [4] 劉天豐,陳建良,林興華.非對(duì)稱管殼式換熱器的三維有限元結(jié)構(gòu)分析[J].壓力容器,2009(8):27-31.

    [5] LIU Jingcheng, ZHANG Shuyou, ZHAO Xinyue, et al. Influence of fin arrangement on fluid flow and heat transfer in the inlet of a plate-fin heat exchanger[J]. Journal of zhejiang university science a,2015,16(4):279-294.

    [6] 劉天豐.非對(duì)稱管殼式換熱器結(jié)構(gòu)分析及改進(jìn)中的問題研究[D].杭州:浙江大學(xué),2005.

    [7] 祝鈴鈺,姜波,任奕丞.基于集合算法的復(fù)雜精餾建模與求解[J].浙江工業(yè)大學(xué)學(xué)報(bào),2014,42(2):152-156.

    [8] 譚蔚,楊星,楊向濤.高參數(shù)換熱器管板熱應(yīng)力分析模型的研究[J].壓力容器,2011,28(2):44-50.

    [9] 徐璋,霍玉雷,陳勇,等.新型組合式垂直軸風(fēng)機(jī)的優(yōu)化設(shè)計(jì)與性能研究[J].浙江工業(yè)大學(xué)學(xué)報(bào),2015,43(3):261-264.

    [10] 郭奇凡.非對(duì)稱管殼式換熱器結(jié)構(gòu)分析及改進(jìn)中的問題研究[J].中國(guó)化工貿(mào)易,2013(4):90-90.

    [11] 中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).壓力容器:GB 150.1-150.4—2011[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2011.

    [12] 原全國(guó)壓力容器標(biāo)準(zhǔn)化技術(shù)委員會(huì)(SAC/TC 262).鋼制壓力容器—分析設(shè)計(jì)標(biāo)準(zhǔn):JB/T 4732—2005[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2005.

    [13] 王金龍,王清明,王偉章.ANSYS12.0有限元分析與范例解析[M].北京:機(jī)械工業(yè)出版社,2010:203-241.

    [14] 陳水生,屠立群,孫偉明.基于ANSYS優(yōu)化技術(shù)的不等截面塔連接部位強(qiáng)度設(shè)計(jì)[J].化工裝備技術(shù),2006(4):56-59.

    [15] 龔曙光,謝桂蘭,黃云清.ANSYS操作命令與參數(shù)化編程[M].北京:機(jī)械工業(yè)出版社,2004:343-355.

    (責(zé)任編輯:陳石平)

    Strength analysis and design of asymmetric tube sheet based on ANSYS optimization techniques

    SUN Weiming1, SHI Xiuzhen1, ZHOU Wen2, CHANG Bo1, LI Jianhong1

    (1.Institute of Process Equipment and Control Engineering, Zhejiang University of Technology, Hangzhou 310014, China;2.Quzhou Special Equipment Inspection Center, Quzhou 324000, China)

    Unsymmetrical shell and tube evaporator is the core equipment, and it is used in the low temperature waste heat power generation process in recent years. Because of the tube bundle is asymmetric so the stress of the unsymmetrical tube sheet is very complicated, and there is no specific formula to calculate the thickness of the unsymmetrical tube sheet. From the safety purpose and the minimum cost, the tube sheet thickness is selected as design variable, component quality is selected as objective function. The stress of tube sheet is analyzed by ANSYS finite element method, and the thickness of the tube sheet is optimized and the reasonable result is got finally.

    unsymmetrical tube sheet; finite element method; design optimization

    2016-11-06

    國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局科技計(jì)劃項(xiàng)目(2014QK195)

    孫偉明(1964—),男,浙江杭州人,教授,博士,研究方向?yàn)檫^程裝備的結(jié)構(gòu)完整性技術(shù),E-mail:wmsun@zjut.edu.cn.

    TB653

    A

    1006-4303(2017)04-0366-04

    猜你喜歡
    管區(qū)管板非對(duì)稱
    非正式治理的制度化運(yùn)作:管區(qū)制與鄉(xiāng)鎮(zhèn)治理現(xiàn)代化的路徑選擇
    20世紀(jì)30年代警管區(qū)制在無錫的實(shí)踐*
    國(guó)產(chǎn)690TT合金U形管彎管區(qū)性能測(cè)定
    行政包干制:鄉(xiāng)鎮(zhèn)“管區(qū)”治理的邏輯與機(jī)制
    理論與改革(2021年5期)2021-09-16 11:56:36
    非對(duì)稱Orlicz差體
    大型管板鉆銑一體化裝備的開發(fā)
    點(diǎn)數(shù)不超過20的旗傳遞非對(duì)稱2-設(shè)計(jì)
    管板堆焊鎳基合金625焊接工藝
    焊接(2016年6期)2016-02-27 13:04:57
    非對(duì)稱負(fù)載下矩陣變換器改進(jìn)型PI重復(fù)控制
    HVM膜過濾器管板改造總結(jié)
    精品久久久久久久末码| 亚洲熟妇中文字幕五十中出| 美女高潮到喷水免费观看| 一级片免费观看大全| 日本一区二区免费在线视频| 精品一区二区三区av网在线观看| 精品国产超薄肉色丝袜足j| 国产一区在线观看成人免费| 淫秽高清视频在线观看| 热99re8久久精品国产| 黄色丝袜av网址大全| 亚洲最大成人中文| www日本黄色视频网| 亚洲av电影在线进入| 国产精华一区二区三区| 亚洲欧洲精品一区二区精品久久久| 老司机深夜福利视频在线观看| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久毛片微露脸| 一本一本综合久久| 免费在线观看日本一区| 欧美亚洲日本最大视频资源| 黄色视频不卡| tocl精华| 亚洲 欧美一区二区三区| 亚洲成人免费电影在线观看| 99精品欧美一区二区三区四区| 69av精品久久久久久| 国产区一区二久久| 一本大道久久a久久精品| 91字幕亚洲| 国产欧美日韩精品亚洲av| 母亲3免费完整高清在线观看| 国产亚洲欧美98| 精品高清国产在线一区| 日本五十路高清| 国产精品日韩av在线免费观看| 给我免费播放毛片高清在线观看| 国产高清videossex| 亚洲一区二区三区不卡视频| 韩国精品一区二区三区| 视频区欧美日本亚洲| 青草久久国产| 观看免费一级毛片| 麻豆成人av在线观看| 一a级毛片在线观看| 一夜夜www| 国产区一区二久久| 中文资源天堂在线| 美女大奶头视频| 叶爱在线成人免费视频播放| xxxwww97欧美| 日本熟妇午夜| 婷婷六月久久综合丁香| 色精品久久人妻99蜜桃| 亚洲第一欧美日韩一区二区三区| 成在线人永久免费视频| 国产精品,欧美在线| 超碰成人久久| 日韩精品中文字幕看吧| 亚洲精品在线美女| 黄色视频不卡| 亚洲免费av在线视频| 久久久久久久久中文| 日韩免费av在线播放| 十八禁网站免费在线| 嫁个100分男人电影在线观看| 国内久久婷婷六月综合欲色啪| 手机成人av网站| 啦啦啦免费观看视频1| 国产亚洲精品久久久久5区| 亚洲成人久久爱视频| 免费在线观看日本一区| 日韩中文字幕欧美一区二区| 欧美日韩精品网址| 在线观看午夜福利视频| 久久中文字幕人妻熟女| 两个人视频免费观看高清| а√天堂www在线а√下载| 国产亚洲欧美98| 香蕉久久夜色| 欧美黑人精品巨大| 久久天躁狠狠躁夜夜2o2o| 久久伊人香网站| 亚洲一区中文字幕在线| 亚洲,欧美精品.| 久久久久九九精品影院| 成人欧美大片| 亚洲国产精品999在线| 欧美日韩乱码在线| 色av中文字幕| 侵犯人妻中文字幕一二三四区| 久久精品国产亚洲av香蕉五月| 亚洲黑人精品在线| 人妻丰满熟妇av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久毛片微露脸| av片东京热男人的天堂| 欧美不卡视频在线免费观看 | 国产精品久久久av美女十八| 亚洲一码二码三码区别大吗| 国产又黄又爽又无遮挡在线| 一边摸一边做爽爽视频免费| 亚洲av成人不卡在线观看播放网| 亚洲 欧美一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 国产1区2区3区精品| 不卡av一区二区三区| 12—13女人毛片做爰片一| 丰满的人妻完整版| 一夜夜www| 十分钟在线观看高清视频www| 看免费av毛片| 12—13女人毛片做爰片一| 岛国视频午夜一区免费看| 1024视频免费在线观看| 18禁黄网站禁片午夜丰满| 欧美久久黑人一区二区| 欧美激情久久久久久爽电影| 亚洲自拍偷在线| 听说在线观看完整版免费高清| 国产成人一区二区三区免费视频网站| 国产又黄又爽又无遮挡在线| 国产亚洲av高清不卡| 人人妻人人澡欧美一区二区| a级毛片a级免费在线| 久久这里只有精品19| 可以在线观看毛片的网站| 久久香蕉精品热| 97人妻精品一区二区三区麻豆 | 满18在线观看网站| 午夜激情福利司机影院| 日日夜夜操网爽| 麻豆久久精品国产亚洲av| 一a级毛片在线观看| 国产亚洲欧美在线一区二区| 成人手机av| 欧美另类亚洲清纯唯美| av电影中文网址| 女人高潮潮喷娇喘18禁视频| 国产精品精品国产色婷婷| 18禁黄网站禁片午夜丰满| 男女那种视频在线观看| 变态另类成人亚洲欧美熟女| av欧美777| 成人午夜高清在线视频 | 不卡一级毛片| 精品乱码久久久久久99久播| 欧美在线一区亚洲| 免费在线观看亚洲国产| 啦啦啦韩国在线观看视频| 女同久久另类99精品国产91| 美女大奶头视频| 久久香蕉国产精品| 老鸭窝网址在线观看| 性欧美人与动物交配| 午夜日韩欧美国产| 亚洲专区中文字幕在线| 18禁国产床啪视频网站| 亚洲国产精品成人综合色| 两性夫妻黄色片| 亚洲成a人片在线一区二区| 国产亚洲精品一区二区www| 欧美在线黄色| 亚洲成av人片免费观看| 免费一级毛片在线播放高清视频| 男女下面进入的视频免费午夜 | 国产在线精品亚洲第一网站| 国产精品久久久久久亚洲av鲁大| 国内毛片毛片毛片毛片毛片| 亚洲成人久久爱视频| 亚洲男人天堂网一区| 色综合站精品国产| 男人的好看免费观看在线视频 | 男女之事视频高清在线观看| 给我免费播放毛片高清在线观看| 国产精品亚洲av一区麻豆| 国产人伦9x9x在线观看| 757午夜福利合集在线观看| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕在线视频| 亚洲欧美精品综合一区二区三区| 大型黄色视频在线免费观看| 国产成年人精品一区二区| 亚洲精品美女久久av网站| 欧美激情高清一区二区三区| 国内毛片毛片毛片毛片毛片| 欧美日韩亚洲国产一区二区在线观看| 999久久久精品免费观看国产| 亚洲久久久国产精品| 夜夜夜夜夜久久久久| av免费在线观看网站| 国产黄a三级三级三级人| 女人被狂操c到高潮| 老司机在亚洲福利影院| 男女下面进入的视频免费午夜 | 一级毛片高清免费大全| www日本在线高清视频| 91字幕亚洲| 国产免费av片在线观看野外av| 18禁美女被吸乳视频| 国内精品久久久久久久电影| 免费在线观看日本一区| 亚洲国产高清在线一区二区三 | e午夜精品久久久久久久| 婷婷精品国产亚洲av| 黄片小视频在线播放| 国产成人一区二区三区免费视频网站| 在线观看免费视频日本深夜| 免费女性裸体啪啪无遮挡网站| 级片在线观看| 精品久久久久久,| 一本久久中文字幕| 久久久久国产精品人妻aⅴ院| 欧美黄色片欧美黄色片| 欧美另类亚洲清纯唯美| 波多野结衣高清无吗| 最近最新免费中文字幕在线| 黑人操中国人逼视频| 久久99热这里只有精品18| 亚洲人成电影免费在线| 日本成人三级电影网站| 俄罗斯特黄特色一大片| 精品第一国产精品| 国产又色又爽无遮挡免费看| 欧美色欧美亚洲另类二区| 在线观看66精品国产| 露出奶头的视频| 国产片内射在线| 人人妻,人人澡人人爽秒播| 国产精品 国内视频| 久久国产亚洲av麻豆专区| 国内揄拍国产精品人妻在线 | 亚洲三区欧美一区| 69av精品久久久久久| 窝窝影院91人妻| 一a级毛片在线观看| 亚洲人成电影免费在线| 日本黄色视频三级网站网址| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频| 国产精品免费一区二区三区在线| 国产97色在线日韩免费| 少妇粗大呻吟视频| 男人的好看免费观看在线视频 | 最好的美女福利视频网| 久9热在线精品视频| 美女国产高潮福利片在线看| videosex国产| 怎么达到女性高潮| 少妇被粗大的猛进出69影院| 久久精品国产99精品国产亚洲性色| 午夜激情av网站| 欧美激情极品国产一区二区三区| 欧美久久黑人一区二区| 久久婷婷人人爽人人干人人爱| 无人区码免费观看不卡| 妹子高潮喷水视频| www.www免费av| av电影中文网址| 亚洲国产欧洲综合997久久, | 欧美午夜高清在线| 美女扒开内裤让男人捅视频| 欧美日本亚洲视频在线播放| bbb黄色大片| 久热爱精品视频在线9| 淫秽高清视频在线观看| 国产99久久九九免费精品| av欧美777| 国内精品久久久久久久电影| 欧美一区二区精品小视频在线| 日本 欧美在线| 欧美日本视频| 亚洲无线在线观看| 啦啦啦韩国在线观看视频| 婷婷六月久久综合丁香| 国产亚洲欧美精品永久| 亚洲自拍偷在线| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产亚洲在线| 免费人成视频x8x8入口观看| 美国免费a级毛片| 亚洲熟妇熟女久久| 淫妇啪啪啪对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产视频内射| 久久久久亚洲av毛片大全| 欧美乱妇无乱码| 亚洲专区字幕在线| 在线观看日韩欧美| 一区二区三区国产精品乱码| 国产精品综合久久久久久久免费| 亚洲精品色激情综合| 精品无人区乱码1区二区| 熟妇人妻久久中文字幕3abv| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 午夜久久久久精精品| 日本黄色视频三级网站网址| 午夜a级毛片| 色精品久久人妻99蜜桃| 欧美又色又爽又黄视频| 免费电影在线观看免费观看| 好男人在线观看高清免费视频 | 欧美日本亚洲视频在线播放| cao死你这个sao货| 欧美黑人精品巨大| a在线观看视频网站| 成人18禁在线播放| 一级a爱视频在线免费观看| 男女床上黄色一级片免费看| 亚洲人成伊人成综合网2020| 欧美精品啪啪一区二区三区| 亚洲第一av免费看| 青草久久国产| 精品高清国产在线一区| 午夜影院日韩av| 亚洲精品中文字幕一二三四区| 日本一本二区三区精品| 亚洲全国av大片| 少妇裸体淫交视频免费看高清 | av视频在线观看入口| 国产片内射在线| 夜夜看夜夜爽夜夜摸| 亚洲av电影在线进入| 亚洲男人的天堂狠狠| 麻豆久久精品国产亚洲av| 日韩欧美免费精品| 男女那种视频在线观看| 美女大奶头视频| 俺也久久电影网| 久久精品国产综合久久久| 欧美午夜高清在线| 三级毛片av免费| 一本大道久久a久久精品| 亚洲精品美女久久久久99蜜臀| 亚洲第一欧美日韩一区二区三区| 欧美亚洲日本最大视频资源| 国产精品九九99| 久久久精品国产亚洲av高清涩受| 免费在线观看视频国产中文字幕亚洲| 精品国产美女av久久久久小说| 麻豆国产av国片精品| 女人爽到高潮嗷嗷叫在线视频| 18禁国产床啪视频网站| 欧美成人免费av一区二区三区| 亚洲av成人不卡在线观看播放网| 国产精华一区二区三区| 色老头精品视频在线观看| 露出奶头的视频| 淫妇啪啪啪对白视频| 久久婷婷人人爽人人干人人爱| 操出白浆在线播放| 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 日本一本二区三区精品| 在线观看一区二区三区| 中文字幕精品亚洲无线码一区 | 午夜福利免费观看在线| 久久伊人香网站| 99久久久亚洲精品蜜臀av| 99久久无色码亚洲精品果冻| 在线播放国产精品三级| 美女免费视频网站| a在线观看视频网站| 久久国产精品人妻蜜桃| 午夜福利高清视频| 久久久久久久精品吃奶| 50天的宝宝边吃奶边哭怎么回事| 午夜视频精品福利| 欧美黄色片欧美黄色片| 老司机福利观看| 亚洲精品在线观看二区| 这个男人来自地球电影免费观看| 久久国产亚洲av麻豆专区| 在线观看一区二区三区| 热99re8久久精品国产| 欧美三级亚洲精品| 久久久久久久午夜电影| 精品久久久久久久久久免费视频| 亚洲国产精品sss在线观看| 国产精品美女特级片免费视频播放器 | 午夜a级毛片| www国产在线视频色| 久久婷婷成人综合色麻豆| 国产精品精品国产色婷婷| 在线播放国产精品三级| 午夜福利欧美成人| 麻豆成人午夜福利视频| 曰老女人黄片| 美女高潮到喷水免费观看| 午夜免费鲁丝| 久久精品夜夜夜夜夜久久蜜豆 | videosex国产| 国产精品一区二区免费欧美| 日韩欧美一区二区三区在线观看| 97超级碰碰碰精品色视频在线观看| 成年版毛片免费区| 露出奶头的视频| 国产高清videossex| 亚洲真实伦在线观看| www.www免费av| 美国免费a级毛片| 十八禁网站免费在线| 亚洲av电影不卡..在线观看| 国产高清视频在线播放一区| 人成视频在线观看免费观看| 男男h啪啪无遮挡| 欧美午夜高清在线| 成年免费大片在线观看| 一进一出抽搐gif免费好疼| 免费看十八禁软件| 欧美三级亚洲精品| 免费在线观看黄色视频的| 久久精品人妻少妇| 国产日本99.免费观看| 精品久久久久久,| 91麻豆av在线| 一进一出好大好爽视频| 免费在线观看日本一区| 嫩草影院精品99| 一级毛片高清免费大全| 久久久久久久久久黄片| 精品久久久久久成人av| 久久国产精品男人的天堂亚洲| 精品一区二区三区四区五区乱码| 精品一区二区三区视频在线观看免费| 久久香蕉国产精品| 色尼玛亚洲综合影院| 亚洲精品一区av在线观看| 成年人黄色毛片网站| 日韩精品中文字幕看吧| 亚洲av电影在线进入| 黄色a级毛片大全视频| 在线天堂中文资源库| 一级毛片高清免费大全| 一夜夜www| 国产乱人伦免费视频| 99riav亚洲国产免费| 999久久久精品免费观看国产| 精品熟女少妇八av免费久了| 亚洲美女黄片视频| 欧美国产日韩亚洲一区| 一二三四社区在线视频社区8| 国产1区2区3区精品| 99国产综合亚洲精品| 岛国在线观看网站| 桃红色精品国产亚洲av| 亚洲 国产 在线| 夜夜夜夜夜久久久久| 久久久久久久精品吃奶| 免费人成视频x8x8入口观看| 热re99久久国产66热| 黄色视频不卡| 啦啦啦观看免费观看视频高清| 国产精品久久久久久精品电影 | 国产99白浆流出| 麻豆av在线久日| 国产精品1区2区在线观看.| 老鸭窝网址在线观看| 欧美国产精品va在线观看不卡| 中文字幕另类日韩欧美亚洲嫩草| 女人爽到高潮嗷嗷叫在线视频| 女生性感内裤真人,穿戴方法视频| 国产又黄又爽又无遮挡在线| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 精品不卡国产一区二区三区| 国产黄色小视频在线观看| 看片在线看免费视频| 十八禁网站免费在线| 91大片在线观看| 国产主播在线观看一区二区| 亚洲av电影不卡..在线观看| 亚洲精品一区av在线观看| 亚洲免费av在线视频| 国内揄拍国产精品人妻在线 | 国产成人系列免费观看| 中文字幕最新亚洲高清| 亚洲av日韩精品久久久久久密| 免费在线观看视频国产中文字幕亚洲| 日韩高清综合在线| 少妇粗大呻吟视频| 久久久久国产一级毛片高清牌| 欧美午夜高清在线| 男男h啪啪无遮挡| 日日干狠狠操夜夜爽| 亚洲熟女毛片儿| 国产成人av教育| 国产在线观看jvid| 久久久久久久久免费视频了| 极品教师在线免费播放| 日本成人三级电影网站| 亚洲成人久久爱视频| 亚洲国产毛片av蜜桃av| 久久午夜综合久久蜜桃| 成人亚洲精品一区在线观看| 欧美激情久久久久久爽电影| 午夜精品在线福利| 久久人人精品亚洲av| 久久伊人香网站| 中文字幕人妻丝袜一区二区| 亚洲国产欧美网| 69av精品久久久久久| 动漫黄色视频在线观看| 非洲黑人性xxxx精品又粗又长| 美女午夜性视频免费| 一本精品99久久精品77| 亚洲精品美女久久久久99蜜臀| 黑人操中国人逼视频| 一个人观看的视频www高清免费观看 | 欧美乱妇无乱码| 婷婷精品国产亚洲av在线| 亚洲精品国产一区二区精华液| 十八禁网站免费在线| 无限看片的www在线观看| 国产成人影院久久av| 波多野结衣高清无吗| 国产成人啪精品午夜网站| 啦啦啦 在线观看视频| 在线十欧美十亚洲十日本专区| 嫁个100分男人电影在线观看| 亚洲精品在线观看二区| 正在播放国产对白刺激| 天天添夜夜摸| 变态另类成人亚洲欧美熟女| 午夜激情福利司机影院| 国产蜜桃级精品一区二区三区| 精品久久久久久久人妻蜜臀av| 国产成人精品无人区| 欧美日韩黄片免| 性欧美人与动物交配| 99热只有精品国产| 国产伦一二天堂av在线观看| 久久性视频一级片| 成人三级黄色视频| 亚洲在线自拍视频| 国产伦人伦偷精品视频| 动漫黄色视频在线观看| 可以免费在线观看a视频的电影网站| 日韩欧美三级三区| 国产免费男女视频| 91国产中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 50天的宝宝边吃奶边哭怎么回事| 亚洲av第一区精品v没综合| 身体一侧抽搐| 国产亚洲欧美精品永久| 亚洲精品国产精品久久久不卡| 欧美激情久久久久久爽电影| 亚洲人成电影免费在线| 欧美黑人精品巨大| 亚洲国产欧美网| 18禁裸乳无遮挡免费网站照片 | 国产熟女午夜一区二区三区| 国产又色又爽无遮挡免费看| 国产激情偷乱视频一区二区| 亚洲精华国产精华精| 久久欧美精品欧美久久欧美| 国产精品永久免费网站| or卡值多少钱| 在线观看66精品国产| 欧美色欧美亚洲另类二区| 国产精品免费一区二区三区在线| 国产成+人综合+亚洲专区| 香蕉国产在线看| 精品国产一区二区三区四区第35| 啦啦啦韩国在线观看视频| 亚洲黑人精品在线| 波多野结衣高清无吗| 国产aⅴ精品一区二区三区波| 精品国产乱码久久久久久男人| 国产一区二区三区视频了| 少妇的丰满在线观看| 国产av不卡久久| 777久久人妻少妇嫩草av网站| 99热6这里只有精品| 国产精品,欧美在线| 亚洲一区二区三区不卡视频| 亚洲色图 男人天堂 中文字幕| 午夜福利欧美成人| 中文字幕高清在线视频| 99riav亚洲国产免费| 亚洲国产毛片av蜜桃av| 窝窝影院91人妻| 久久精品91无色码中文字幕| 亚洲精品一区av在线观看| 淫妇啪啪啪对白视频| 久久久久久久精品吃奶| 久久欧美精品欧美久久欧美| xxx96com| 亚洲aⅴ乱码一区二区在线播放 | 国产精品自产拍在线观看55亚洲| 久久精品91蜜桃| 中出人妻视频一区二区| 琪琪午夜伦伦电影理论片6080| 国产精品亚洲一级av第二区| 少妇的丰满在线观看| 亚洲精品粉嫩美女一区| 久热爱精品视频在线9| 成人国语在线视频| 亚洲精品粉嫩美女一区| 女同久久另类99精品国产91| 精品久久蜜臀av无| 色综合站精品国产| 成人特级黄色片久久久久久久| 精品午夜福利视频在线观看一区| 亚洲av日韩精品久久久久久密| 亚洲中文字幕日韩| 午夜福利成人在线免费观看| 成人亚洲精品一区在线观看| 在线观看www视频免费| 99热这里只有精品一区 | 村上凉子中文字幕在线| 99精品在免费线老司机午夜|