• <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在线观看| 久久精品国产综合久久久| 好看av亚洲va欧美ⅴa在| 日本成人三级电影网站| 18禁观看日本| 美女午夜性视频免费| 亚洲久久久国产精品| 亚洲九九香蕉| 香蕉国产在线看| 性欧美人与动物交配| 久久青草综合色| 一级a爱视频在线免费观看| 村上凉子中文字幕在线| 亚洲久久久国产精品| 色精品久久人妻99蜜桃| 99热6这里只有精品| 黑人欧美特级aaaaaa片| 女性被躁到高潮视频| 12—13女人毛片做爰片一| 啦啦啦免费观看视频1| 淫秽高清视频在线观看| 18美女黄网站色大片免费观看| 亚洲成a人片在线一区二区| 国内精品久久久久精免费| 亚洲av五月六月丁香网| av天堂在线播放| 国产欧美日韩精品亚洲av| 精品久久久久久久末码| 国产亚洲欧美98| 亚洲精品av麻豆狂野| 一级黄色大片毛片| 午夜影院日韩av| 99在线人妻在线中文字幕| 午夜福利在线观看吧| 最新在线观看一区二区三区| 国产99久久九九免费精品| 欧美日韩精品网址| av天堂在线播放| 国产精品九九99| 久久九九热精品免费| 狂野欧美激情性xxxx| 日韩欧美在线二视频| 狠狠狠狠99中文字幕| cao死你这个sao货| 国产精品免费一区二区三区在线| 99re在线观看精品视频| 在线视频色国产色| 午夜影院日韩av| 好男人电影高清在线观看| 亚洲av五月六月丁香网| 国内精品久久久久精免费| 成人三级做爰电影| 亚洲五月婷婷丁香| 日本撒尿小便嘘嘘汇集6| 欧美日韩亚洲国产一区二区在线观看| 两个人看的免费小视频| 亚洲专区国产一区二区| 搡老岳熟女国产| 欧美日韩亚洲综合一区二区三区_| 男人操女人黄网站| 99久久综合精品五月天人人| netflix在线观看网站| 丝袜在线中文字幕| 亚洲精品色激情综合| 日韩大码丰满熟妇| 成人18禁在线播放| 亚洲午夜精品一区,二区,三区| 亚洲国产精品sss在线观看| 嫩草影院精品99| 三级毛片av免费| 成人一区二区视频在线观看| 国内少妇人妻偷人精品xxx网站 | 免费看美女性在线毛片视频| 好男人电影高清在线观看| 一区二区三区高清视频在线| 午夜福利高清视频| 久久99热这里只有精品18| 长腿黑丝高跟| 亚洲熟妇中文字幕五十中出| 一a级毛片在线观看| 午夜福利18| 久久久久久人人人人人| 好男人电影高清在线观看| 成人av一区二区三区在线看| 在线观看66精品国产| 黄色成人免费大全| 曰老女人黄片| 老司机在亚洲福利影院| 亚洲精品一卡2卡三卡4卡5卡| 一级a爱视频在线免费观看| 日韩 欧美 亚洲 中文字幕| 91av网站免费观看| 免费在线观看黄色视频的| 国内精品久久久久精免费| 亚洲精品在线美女| 日韩欧美三级三区| 高清毛片免费观看视频网站| 国产aⅴ精品一区二区三区波| 国产精品99久久99久久久不卡| 国产精品久久久av美女十八| 国产精品av久久久久免费| 男人操女人黄网站| 看免费av毛片| 免费在线观看影片大全网站| 国产午夜精品久久久久久| 国产极品粉嫩免费观看在线| 99国产精品99久久久久| 亚洲精品在线观看二区| 香蕉丝袜av| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 侵犯人妻中文字幕一二三四区| 久久久久久免费高清国产稀缺| 天天躁狠狠躁夜夜躁狠狠躁| 18禁美女被吸乳视频| 久久草成人影院| 超碰成人久久| 美女午夜性视频免费| 一个人观看的视频www高清免费观看 | 搞女人的毛片| 91字幕亚洲| 18禁国产床啪视频网站| 天堂√8在线中文| videosex国产| 久久久久九九精品影院| 国产久久久一区二区三区| 一本综合久久免费| 此物有八面人人有两片| 亚洲熟妇熟女久久| 久久性视频一级片| 色综合欧美亚洲国产小说| 国产精品野战在线观看| 黄色 视频免费看| 久久精品国产清高在天天线| 久久亚洲真实| e午夜精品久久久久久久| 精品高清国产在线一区| 777久久人妻少妇嫩草av网站| 嫩草影视91久久| 国产av不卡久久| 最近最新中文字幕大全免费视频| 色精品久久人妻99蜜桃| 国产国语露脸激情在线看| 亚洲第一欧美日韩一区二区三区| 午夜久久久在线观看| 国产成人av激情在线播放| 久久精品国产综合久久久| 国产精品98久久久久久宅男小说| 久久婷婷成人综合色麻豆| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区三| 淫秽高清视频在线观看| 大型av网站在线播放| 国产极品粉嫩免费观看在线| 久久久久久人人人人人| 九色国产91popny在线| 亚洲成国产人片在线观看| 午夜久久久在线观看| 啦啦啦韩国在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲午夜精品一区,二区,三区| 99久久久亚洲精品蜜臀av| 一区福利在线观看| 欧美性猛交黑人性爽| av超薄肉色丝袜交足视频| 老司机午夜十八禁免费视频| 亚洲国产中文字幕在线视频| 亚洲人成网站高清观看| 欧美日本视频| 老熟妇仑乱视频hdxx| 成在线人永久免费视频| 美女扒开内裤让男人捅视频| 国产视频一区二区在线看| 在线天堂中文资源库| 国产区一区二久久| 在线十欧美十亚洲十日本专区| 国产免费av片在线观看野外av| 后天国语完整版免费观看| 波多野结衣巨乳人妻| 美女免费视频网站| 久久久久九九精品影院| 精品高清国产在线一区| 国产精品久久久av美女十八| 69av精品久久久久久| www.熟女人妻精品国产| 老司机靠b影院| 色综合站精品国产| 成年女人毛片免费观看观看9| 在线观看66精品国产| 中文字幕人妻丝袜一区二区| 国产精品美女特级片免费视频播放器 | 啦啦啦 在线观看视频| 很黄的视频免费| 成人永久免费在线观看视频| 国产精品香港三级国产av潘金莲| 国产成人影院久久av| 一级毛片女人18水好多| 女人高潮潮喷娇喘18禁视频| 日韩三级视频一区二区三区| 亚洲av美国av| 每晚都被弄得嗷嗷叫到高潮| 精华霜和精华液先用哪个| av有码第一页| 久久性视频一级片| 日本三级黄在线观看| 亚洲欧美精品综合久久99| 丁香六月欧美| 国产黄片美女视频| 亚洲精品一区av在线观看| 久久中文看片网| 两个人视频免费观看高清| 欧美人与性动交α欧美精品济南到| 性欧美人与动物交配| 真人做人爱边吃奶动态| 久热爱精品视频在线9| 国产av在哪里看| 中文字幕人妻丝袜一区二区| 亚洲五月天丁香| 亚洲成a人片在线一区二区| 成人一区二区视频在线观看| 亚洲美女黄片视频| 18美女黄网站色大片免费观看| 可以在线观看毛片的网站| 一边摸一边做爽爽视频免费| 天堂动漫精品| 首页视频小说图片口味搜索| 久久香蕉精品热| 国产精品日韩av在线免费观看| 精品日产1卡2卡| 亚洲三区欧美一区| 免费女性裸体啪啪无遮挡网站| 亚洲欧美日韩高清在线视频| 久久狼人影院| 午夜激情福利司机影院| 国产亚洲精品久久久久久毛片| 成在线人永久免费视频| 91老司机精品| 97人妻精品一区二区三区麻豆 | 人妻丰满熟妇av一区二区三区| 亚洲美女黄片视频| 禁无遮挡网站| 一本大道久久a久久精品| 欧美激情 高清一区二区三区| 亚洲人成伊人成综合网2020| 亚洲精品一卡2卡三卡4卡5卡| a级毛片a级免费在线| √禁漫天堂资源中文www| 午夜视频精品福利| 大型av网站在线播放| 丁香六月欧美| 黄色女人牲交| 亚洲免费av在线视频| 国产视频一区二区在线看| 亚洲九九香蕉| 亚洲欧美日韩高清在线视频| 国产不卡一卡二| 18禁黄网站禁片午夜丰满| 亚洲欧美精品综合一区二区三区| 午夜福利免费观看在线| 日韩欧美国产在线观看| 亚洲欧美日韩无卡精品| 后天国语完整版免费观看| 美女国产高潮福利片在线看| 国产一区二区在线av高清观看| 日韩欧美在线二视频| 午夜福利在线在线| 午夜福利视频1000在线观看| 成人永久免费在线观看视频| 老汉色av国产亚洲站长工具| www.精华液| 午夜免费成人在线视频| 国产色视频综合| 色老头精品视频在线观看| 女生性感内裤真人,穿戴方法视频| svipshipincom国产片| 波多野结衣高清无吗| 日韩高清综合在线| 香蕉久久夜色| 国产欧美日韩精品亚洲av| www.精华液| 国产男靠女视频免费网站| 久久性视频一级片| 亚洲av五月六月丁香网| 国产精品av久久久久免费| 50天的宝宝边吃奶边哭怎么回事| 国产主播在线观看一区二区| 欧美黄色淫秽网站| 久久久久久九九精品二区国产 | 在线观看日韩欧美| 给我免费播放毛片高清在线观看| 国产高清视频在线播放一区| 免费看日本二区| 成人国语在线视频| 久久久久国产一级毛片高清牌| 国产精品久久电影中文字幕| 久久久久精品国产欧美久久久| 高潮久久久久久久久久久不卡| 日日爽夜夜爽网站| 午夜日韩欧美国产| 日韩国内少妇激情av| 又黄又爽又免费观看的视频| 亚洲专区国产一区二区| 精品一区二区三区视频在线观看免费| 亚洲熟女毛片儿| 日韩精品青青久久久久久| a级毛片a级免费在线| 国产区一区二久久| 欧美日韩精品网址| 美女高潮喷水抽搐中文字幕| 少妇裸体淫交视频免费看高清 | 久久热在线av| 亚洲国产欧美日韩在线播放| 亚洲狠狠婷婷综合久久图片| 免费观看精品视频网站| 一区福利在线观看| 国产v大片淫在线免费观看| 人人妻,人人澡人人爽秒播| 一级作爱视频免费观看| 国产午夜福利久久久久久| 51午夜福利影视在线观看| 亚洲国产日韩欧美精品在线观看 | 日韩国内少妇激情av| 淫秽高清视频在线观看| 免费一级毛片在线播放高清视频| 美女高潮喷水抽搐中文字幕| 午夜久久久久精精品| 久久狼人影院| 国产v大片淫在线免费观看| 亚洲精品一区av在线观看| 18禁黄网站禁片午夜丰满| 亚洲一区二区三区色噜噜| 丝袜人妻中文字幕| 亚洲国产精品合色在线| 少妇裸体淫交视频免费看高清 | 久久久久国产精品人妻aⅴ院| 成人亚洲精品一区在线观看| 国产1区2区3区精品| 日本精品一区二区三区蜜桃| 最近最新中文字幕大全电影3 | 免费看日本二区| 国产男靠女视频免费网站| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| 国产精品久久视频播放| 久久久久久久久中文| ponron亚洲| 黑人巨大精品欧美一区二区mp4| 亚洲激情在线av| 麻豆国产av国片精品| 亚洲va日本ⅴa欧美va伊人久久| 久热这里只有精品99| 成在线人永久免费视频| e午夜精品久久久久久久| 无限看片的www在线观看| 一进一出好大好爽视频| 日本一区二区免费在线视频| 国产久久久一区二区三区| 国产精品免费视频内射| 国产精品免费一区二区三区在线| 久久性视频一级片| 久久青草综合色| 国产精品九九99| 俺也久久电影网| 一级毛片女人18水好多| 中文亚洲av片在线观看爽| 国产三级黄色录像| 别揉我奶头~嗯~啊~动态视频| 男男h啪啪无遮挡| 国产亚洲av高清不卡| av有码第一页| 男男h啪啪无遮挡| videosex国产| 国产成人啪精品午夜网站| 无遮挡黄片免费观看| av片东京热男人的天堂| 国产亚洲精品久久久久久毛片| 久久久久国内视频| 一级a爱视频在线免费观看| 淫秽高清视频在线观看| 色播在线永久视频| 精品卡一卡二卡四卡免费| 动漫黄色视频在线观看| 日韩av在线大香蕉| 欧美精品亚洲一区二区| 美女免费视频网站| 熟女少妇亚洲综合色aaa.| 国产又爽黄色视频| 青草久久国产| 欧美 亚洲 国产 日韩一| 日韩精品青青久久久久久| 黄片小视频在线播放| 久9热在线精品视频| 精品国产国语对白av| 精品人妻1区二区| or卡值多少钱| 麻豆成人午夜福利视频| 国产精品国产高清国产av| 变态另类丝袜制服| 国内精品久久久久久久电影| 精品国产超薄肉色丝袜足j| 淫秽高清视频在线观看| 国产亚洲精品一区二区www| 免费在线观看完整版高清| 亚洲在线自拍视频| 国产主播在线观看一区二区| www.熟女人妻精品国产| 国产一区在线观看成人免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品中文字幕在线视频| 亚洲国产欧美一区二区综合| 亚洲七黄色美女视频| 久久香蕉国产精品| √禁漫天堂资源中文www| 久久 成人 亚洲| 欧美午夜高清在线| 18禁裸乳无遮挡免费网站照片 | 美国免费a级毛片| 久久精品影院6| 很黄的视频免费| 久久精品亚洲精品国产色婷小说| 一边摸一边做爽爽视频免费| 一个人免费在线观看的高清视频| 久久草成人影院| 亚洲午夜精品一区,二区,三区| 欧美久久黑人一区二区| 久久精品91蜜桃| 久久精品影院6| 夜夜爽天天搞| 国产亚洲精品一区二区www| 亚洲真实伦在线观看| 精品久久久久久久末码| 午夜福利成人在线免费观看| 亚洲午夜精品一区,二区,三区| 国产高清视频在线播放一区| 51午夜福利影视在线观看| 可以免费在线观看a视频的电影网站| 日韩欧美在线二视频| 午夜免费鲁丝| 日韩成人在线观看一区二区三区| 亚洲av第一区精品v没综合| 精品高清国产在线一区| 日本 av在线| 中国美女看黄片| 熟妇人妻久久中文字幕3abv| 黄色女人牲交| 国产成人啪精品午夜网站| 亚洲一区二区三区色噜噜| 精品久久久久久久毛片微露脸| 国产成人欧美在线观看| 男人舔奶头视频| 99久久无色码亚洲精品果冻| 一本综合久久免费| 91麻豆精品激情在线观看国产| 亚洲avbb在线观看| 999久久久国产精品视频| 亚洲成av人片免费观看| 国产精品一区二区免费欧美| 欧美三级亚洲精品| 我的亚洲天堂| 日本成人三级电影网站| 啦啦啦 在线观看视频| 制服诱惑二区| 男男h啪啪无遮挡| 此物有八面人人有两片| 男人操女人黄网站| 亚洲午夜理论影院| 亚洲国产欧洲综合997久久, | 久久久国产成人精品二区| 男人舔女人下体高潮全视频| 在线观看免费日韩欧美大片| 久久久久久九九精品二区国产 | 午夜福利18| 精品免费久久久久久久清纯| 可以免费在线观看a视频的电影网站| 久久久久久免费高清国产稀缺| 国产一区二区激情短视频| 国产成人欧美在线观看| 国产亚洲精品第一综合不卡| 69av精品久久久久久| 亚洲av成人一区二区三| 2021天堂中文幕一二区在线观 | 身体一侧抽搐| 免费看十八禁软件| 精品电影一区二区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产极品粉嫩在线观看| 欧美日韩福利视频一区二区| 亚洲国产精品999在线| 99久久国产精品久久久| 亚洲精品一卡2卡三卡4卡5卡| 麻豆成人午夜福利视频| 国产一区二区三区视频了| 国产视频内射| 日韩中文字幕欧美一区二区| 中文字幕精品亚洲无线码一区 | 国产欧美日韩精品亚洲av| 久久精品国产综合久久久| 成人免费观看视频高清| 女性被躁到高潮视频| АⅤ资源中文在线天堂| 久久精品国产综合久久久| 两个人免费观看高清视频| 久久热在线av| 中文字幕人妻丝袜一区二区| 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 一级a爱视频在线免费观看| 午夜福利免费观看在线| 亚洲精品一卡2卡三卡4卡5卡| 日韩有码中文字幕| 老司机靠b影院| 好男人在线观看高清免费视频 | e午夜精品久久久久久久| 色综合亚洲欧美另类图片| www国产在线视频色| 亚洲男人天堂网一区| 男女那种视频在线观看| av在线天堂中文字幕| 亚洲 欧美一区二区三区| 日日干狠狠操夜夜爽| 欧美亚洲日本最大视频资源| 欧美成狂野欧美在线观看| 亚洲成人精品中文字幕电影| 久久人妻av系列| 午夜影院日韩av| 在线av久久热| 村上凉子中文字幕在线| 日本一本二区三区精品| 欧美激情高清一区二区三区| 日韩精品中文字幕看吧| 99在线视频只有这里精品首页| 美女免费视频网站| 亚洲avbb在线观看| 琪琪午夜伦伦电影理论片6080| 一本大道久久a久久精品| 亚洲av熟女| 久久婷婷成人综合色麻豆| 性欧美人与动物交配| 午夜福利18| 在线观看午夜福利视频| 波多野结衣高清无吗| 级片在线观看| 中文字幕人妻丝袜一区二区| 少妇熟女aⅴ在线视频| 狠狠狠狠99中文字幕| 每晚都被弄得嗷嗷叫到高潮| 香蕉国产在线看| 波多野结衣av一区二区av| 最新美女视频免费是黄的| 久久婷婷人人爽人人干人人爱| 在线国产一区二区在线| 久久久精品欧美日韩精品| 欧美成人午夜精品| 丝袜人妻中文字幕| 久久国产乱子伦精品免费另类| 精品一区二区三区av网在线观看| 日韩三级视频一区二区三区| 亚洲av五月六月丁香网| 欧美在线一区亚洲| 成人国产一区最新在线观看| 免费在线观看成人毛片| 99热只有精品国产| 我的亚洲天堂| 国产精品久久久久久人妻精品电影| 男人舔女人的私密视频| 1024手机看黄色片| 亚洲 国产 在线| 级片在线观看| 长腿黑丝高跟| 美女国产高潮福利片在线看| 欧美黑人巨大hd| 此物有八面人人有两片| 十八禁人妻一区二区| 中文亚洲av片在线观看爽| 99精品欧美一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 亚洲国产中文字幕在线视频| 精品免费久久久久久久清纯| www.熟女人妻精品国产| 国产一卡二卡三卡精品| 亚洲人成伊人成综合网2020| 国产精品一区二区精品视频观看| 色精品久久人妻99蜜桃| 欧美绝顶高潮抽搐喷水| 精品久久久久久久毛片微露脸| 韩国av一区二区三区四区| АⅤ资源中文在线天堂| 国产片内射在线| 午夜福利18| 日韩有码中文字幕| 亚洲国产精品999在线| 亚洲avbb在线观看| 村上凉子中文字幕在线| 日本三级黄在线观看| 88av欧美| 午夜久久久久精精品| 亚洲精品国产精品久久久不卡| 91麻豆av在线| 婷婷六月久久综合丁香| 免费在线观看日本一区| 国产真人三级小视频在线观看| 免费在线观看成人毛片| 搡老岳熟女国产| 两人在一起打扑克的视频| 亚洲全国av大片| 国产成人系列免费观看| 1024视频免费在线观看| 色哟哟哟哟哟哟| 欧美成人一区二区免费高清观看 | 国产精品日韩av在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 欧美黄色片欧美黄色片| 伦理电影免费视频| 国产精品影院久久| 成人手机av|