• <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é)
    亚洲九九香蕉| 国产亚洲欧美98| 午夜久久久久精精品| netflix在线观看网站| 亚洲精品国产区一区二| 欧美激情 高清一区二区三区| 国产伦人伦偷精品视频| 精品一区二区三区四区五区乱码| 夜夜躁狠狠躁天天躁| 在线观看免费视频日本深夜| 宅男免费午夜| 久久性视频一级片| 国内精品久久久久精免费| 亚洲欧美日韩另类电影网站| 国产真人三级小视频在线观看| 国产精品精品国产色婷婷| 久99久视频精品免费| 50天的宝宝边吃奶边哭怎么回事| 老司机靠b影院| 一区福利在线观看| 午夜福利在线观看吧| 成人三级做爰电影| 亚洲午夜精品一区,二区,三区| 精品久久久久久,| 欧美国产精品va在线观看不卡| 自线自在国产av| 亚洲精品美女久久久久99蜜臀| 又大又爽又粗| 精品久久久久久久久久免费视频| 国产精品免费视频内射| 久久精品国产清高在天天线| a级毛片在线看网站| 老司机福利观看| 在线观看舔阴道视频| 大型黄色视频在线免费观看| 久久久久久大精品| 精品无人区乱码1区二区| 日日爽夜夜爽网站| 国产欧美日韩一区二区三区在线| 久久久国产精品麻豆| 久99久视频精品免费| 成人亚洲精品av一区二区| 欧美国产日韩亚洲一区| 国产精品亚洲一级av第二区| 久久中文看片网| 美女高潮到喷水免费观看| 69av精品久久久久久| 最近最新免费中文字幕在线| 久久久久久久精品吃奶| 大陆偷拍与自拍| 黑丝袜美女国产一区| 亚洲国产精品999在线| 国产黄a三级三级三级人| 成人av一区二区三区在线看| 午夜免费激情av| 亚洲 国产 在线| 丝袜人妻中文字幕| 亚洲第一av免费看| 无限看片的www在线观看| 大型av网站在线播放| 人人妻人人澡欧美一区二区 | 亚洲av片天天在线观看| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区三区四区久久 | 99国产精品一区二区蜜桃av| 波多野结衣高清无吗| netflix在线观看网站| 欧美日韩福利视频一区二区| 国产一区二区在线av高清观看| 亚洲欧美日韩另类电影网站| 国产高清videossex| 午夜福利视频1000在线观看 | 大码成人一级视频| 国产精品久久久久久人妻精品电影| 亚洲专区中文字幕在线| 久久青草综合色| 欧美人与性动交α欧美精品济南到| 亚洲激情在线av| 国产精品久久久av美女十八| www日本在线高清视频| 久久草成人影院| 男人舔女人下体高潮全视频| 看片在线看免费视频| 国产精品久久视频播放| 一级a爱片免费观看的视频| 可以在线观看毛片的网站| 又大又爽又粗| 亚洲国产精品999在线| 在线国产一区二区在线| 制服诱惑二区| 国产精品久久久久久精品电影 | 亚洲色图 男人天堂 中文字幕| 午夜a级毛片| 在线观看舔阴道视频| 欧美成狂野欧美在线观看| 久久狼人影院| 国产1区2区3区精品| 日韩av在线大香蕉| 老司机靠b影院| 欧美久久黑人一区二区| 人成视频在线观看免费观看| 999精品在线视频| 国产97色在线日韩免费| 国产xxxxx性猛交| 免费无遮挡裸体视频| 久久久国产成人精品二区| or卡值多少钱| 国产野战对白在线观看| 免费在线观看视频国产中文字幕亚洲| av片东京热男人的天堂| 俄罗斯特黄特色一大片| 国产av在哪里看| 熟妇人妻久久中文字幕3abv| 人妻久久中文字幕网| 一级,二级,三级黄色视频| 亚洲成av片中文字幕在线观看| 每晚都被弄得嗷嗷叫到高潮| 一级毛片女人18水好多| 视频在线观看一区二区三区| 精品人妻在线不人妻| 国产精品国产高清国产av| 久久久久久亚洲精品国产蜜桃av| 久久人妻熟女aⅴ| 搡老岳熟女国产| 村上凉子中文字幕在线| 国产成人系列免费观看| 九色亚洲精品在线播放| 制服人妻中文乱码| 亚洲国产精品成人综合色| 97人妻精品一区二区三区麻豆 | 精品久久久久久久久久免费视频| 亚洲精品久久国产高清桃花| 日本 av在线| 在线观看免费视频日本深夜| 久久热在线av| 久久九九热精品免费| 国产黄a三级三级三级人| 日本黄色视频三级网站网址| 国产人伦9x9x在线观看| 欧美日韩乱码在线| 99国产极品粉嫩在线观看| 国产成人啪精品午夜网站| 亚洲七黄色美女视频| av免费在线观看网站| 亚洲欧美激情综合另类| 国产免费av片在线观看野外av| 又黄又爽又免费观看的视频| 少妇的丰满在线观看| 亚洲av日韩精品久久久久久密| 精品国产一区二区三区四区第35| 亚洲专区字幕在线| 91大片在线观看| 免费在线观看黄色视频的| 亚洲av日韩精品久久久久久密| 男人舔女人下体高潮全视频| 亚洲 国产 在线| 欧美日韩瑟瑟在线播放| 韩国av一区二区三区四区| 麻豆一二三区av精品| 看黄色毛片网站| 久久人妻av系列| 窝窝影院91人妻| 欧美日韩亚洲国产一区二区在线观看| 大码成人一级视频| 色综合欧美亚洲国产小说| 久久久久国产一级毛片高清牌| 成人亚洲精品av一区二区| 校园春色视频在线观看| 亚洲成a人片在线一区二区| 精品国产一区二区久久| 免费在线观看黄色视频的| 老司机午夜十八禁免费视频| 97超级碰碰碰精品色视频在线观看| 精品无人区乱码1区二区| 亚洲欧美激情综合另类| 久久久久国产一级毛片高清牌| 少妇的丰满在线观看| 91麻豆av在线| 精品熟女少妇八av免费久了| 99久久久亚洲精品蜜臀av| 宅男免费午夜| 夜夜夜夜夜久久久久| 亚洲国产精品久久男人天堂| 香蕉国产在线看| 一个人观看的视频www高清免费观看 | 在线观看www视频免费| 亚洲自偷自拍图片 自拍| 91在线观看av| 又黄又粗又硬又大视频| 天天一区二区日本电影三级 | www.999成人在线观看| 99国产精品免费福利视频| 天堂动漫精品| АⅤ资源中文在线天堂| 久久久国产精品麻豆| 大陆偷拍与自拍| 大陆偷拍与自拍| 性色av乱码一区二区三区2| 亚洲黑人精品在线| 久久亚洲真实| 久久午夜亚洲精品久久| 久久香蕉激情| 国产欧美日韩综合在线一区二区| 久久性视频一级片| 欧美乱妇无乱码| 久久九九热精品免费| 99久久精品国产亚洲精品| 中文字幕av电影在线播放| 女性被躁到高潮视频| 欧美精品啪啪一区二区三区| 亚洲自偷自拍图片 自拍| 午夜成年电影在线免费观看| cao死你这个sao货| 国产不卡一卡二| 亚洲精品中文字幕在线视频| 成在线人永久免费视频| 国产高清视频在线播放一区| 国产在线精品亚洲第一网站| 国产亚洲av嫩草精品影院| 国产精品二区激情视频| 此物有八面人人有两片| 波多野结衣巨乳人妻| 9色porny在线观看| 午夜免费鲁丝| 久久精品国产综合久久久| 女性被躁到高潮视频| 久久精品91蜜桃| 国产高清视频在线播放一区| 色在线成人网| 亚洲av成人av| 欧美日韩乱码在线| 免费人成视频x8x8入口观看| 精品欧美国产一区二区三| 免费不卡黄色视频| 亚洲精品国产一区二区精华液| 久久午夜综合久久蜜桃| 国产精品国产高清国产av| 亚洲精品中文字幕在线视频| 国产真人三级小视频在线观看| 一进一出好大好爽视频| 视频区欧美日本亚洲| 日本 欧美在线| 可以在线观看的亚洲视频| 亚洲av电影不卡..在线观看| 亚洲色图综合在线观看| av电影中文网址| 一二三四在线观看免费中文在| 97碰自拍视频| 亚洲片人在线观看| videosex国产| 久久精品成人免费网站| x7x7x7水蜜桃| 国产av一区在线观看免费| 国产1区2区3区精品| 在线观看66精品国产| 久久久久久人人人人人| 国产av精品麻豆| 亚洲avbb在线观看| 国产精品一区二区精品视频观看| 老熟妇仑乱视频hdxx| 亚洲片人在线观看| 国产精品日韩av在线免费观看 | 国产片内射在线| 久久 成人 亚洲| 一边摸一边抽搐一进一出视频| 亚洲中文字幕日韩| 97碰自拍视频| 91精品三级在线观看| 老鸭窝网址在线观看| 操出白浆在线播放| 亚洲天堂国产精品一区在线| 一本大道久久a久久精品| 成人av一区二区三区在线看| 男女下面进入的视频免费午夜 | 色婷婷久久久亚洲欧美| 精品电影一区二区在线| 亚洲视频免费观看视频| 国产成人精品久久二区二区91| 久久这里只有精品19| 黄片播放在线免费| 9热在线视频观看99| 免费看a级黄色片| 天堂√8在线中文| 视频在线观看一区二区三区| 成人国语在线视频| 色综合婷婷激情| 国产高清视频在线播放一区| 久久精品国产清高在天天线| 日本 欧美在线| 91麻豆av在线| 成人精品一区二区免费| 黄色 视频免费看| 久久天堂一区二区三区四区| 日本a在线网址| 国产极品粉嫩免费观看在线| 美女扒开内裤让男人捅视频| 韩国av一区二区三区四区| ponron亚洲| 日韩欧美一区二区三区在线观看| 18禁美女被吸乳视频| 国产欧美日韩一区二区精品| 最新美女视频免费是黄的| 亚洲最大成人中文| 亚洲自拍偷在线| 黄色视频,在线免费观看| 一级片免费观看大全| 狂野欧美激情性xxxx| 热re99久久国产66热| 成人av一区二区三区在线看| 美女高潮喷水抽搐中文字幕| 国产精华一区二区三区| 午夜成年电影在线免费观看| 免费在线观看视频国产中文字幕亚洲| 国产高清视频在线播放一区| 国产精品免费一区二区三区在线| 亚洲aⅴ乱码一区二区在线播放 | 黄网站色视频无遮挡免费观看| 丰满的人妻完整版| 欧美人与性动交α欧美精品济南到| 精品欧美一区二区三区在线| 搡老熟女国产l中国老女人| 亚洲中文字幕日韩| 国产主播在线观看一区二区| 亚洲性夜色夜夜综合| 两性夫妻黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产综合久久久| 一区福利在线观看| 男人舔女人下体高潮全视频| 成在线人永久免费视频| 久久香蕉激情| 国产99久久九九免费精品| 在线观看免费视频日本深夜| 亚洲第一av免费看| 亚洲狠狠婷婷综合久久图片| 国产熟女午夜一区二区三区| 亚洲第一电影网av| 亚洲午夜精品一区,二区,三区| 99国产极品粉嫩在线观看| 久久精品影院6| 给我免费播放毛片高清在线观看| 在线播放国产精品三级| 两性午夜刺激爽爽歪歪视频在线观看 | 日本五十路高清| 国产又爽黄色视频| а√天堂www在线а√下载| 成年女人毛片免费观看观看9| 国产精品野战在线观看| 欧美+亚洲+日韩+国产| 90打野战视频偷拍视频| 国产极品粉嫩免费观看在线| 窝窝影院91人妻| 国产一级毛片七仙女欲春2 | 伊人久久大香线蕉亚洲五| 国产欧美日韩精品亚洲av| 波多野结衣一区麻豆| 人人澡人人妻人| 成年版毛片免费区| 男女午夜视频在线观看| 亚洲狠狠婷婷综合久久图片| x7x7x7水蜜桃| 99久久99久久久精品蜜桃| 欧美+亚洲+日韩+国产| 亚洲无线在线观看| 国产成人一区二区三区免费视频网站| 日日摸夜夜添夜夜添小说| 亚洲精品国产区一区二| 日本精品一区二区三区蜜桃| 香蕉国产在线看| 99精品在免费线老司机午夜| av天堂在线播放| videosex国产| 日本vs欧美在线观看视频| 在线国产一区二区在线| 亚洲av成人不卡在线观看播放网| 午夜福利免费观看在线| 成人国语在线视频| 亚洲自偷自拍图片 自拍| 亚洲国产毛片av蜜桃av| 国产高清有码在线观看视频 | 欧美一级a爱片免费观看看 | 久久精品成人免费网站| 午夜久久久在线观看| 日韩精品中文字幕看吧| av欧美777| 日本a在线网址| 午夜福利成人在线免费观看| 又黄又爽又免费观看的视频| 久久婷婷人人爽人人干人人爱 | 精品久久久久久久人妻蜜臀av | 国产精品一区二区三区四区久久 | 黄片小视频在线播放| 国产精华一区二区三区| 精品国内亚洲2022精品成人| 97碰自拍视频| 久久人妻福利社区极品人妻图片| 狂野欧美激情性xxxx| 国产精品一区二区精品视频观看| 国产成年人精品一区二区| 麻豆国产av国片精品| 精品国产国语对白av| 国产精品久久久av美女十八| 电影成人av| 两个人看的免费小视频| 亚洲精品中文字幕在线视频| 色综合站精品国产| 亚洲一区高清亚洲精品| 手机成人av网站| 久久中文看片网| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产欧美网| 久久亚洲精品不卡| 日本 欧美在线| 国产精品一区二区免费欧美| 精品电影一区二区在线| 亚洲狠狠婷婷综合久久图片| 久久久国产成人精品二区| 欧美精品亚洲一区二区| 国产欧美日韩一区二区三| 欧美色视频一区免费| 午夜日韩欧美国产| 欧美在线黄色| 日韩大码丰满熟妇| 9色porny在线观看| 久久人人97超碰香蕉20202| 高清黄色对白视频在线免费看| 黄色丝袜av网址大全| 少妇粗大呻吟视频| 国产一区二区三区在线臀色熟女| 国产成人欧美在线观看| 久热爱精品视频在线9| 亚洲熟妇熟女久久| 久久久国产成人精品二区| av免费在线观看网站| www.999成人在线观看| 欧美色欧美亚洲另类二区 | 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合久久99| 亚洲成人久久性| 人成视频在线观看免费观看| 欧美久久黑人一区二区| 男人舔女人的私密视频| 亚洲av第一区精品v没综合| 最新美女视频免费是黄的| 97碰自拍视频| 亚洲七黄色美女视频| 高潮久久久久久久久久久不卡| 午夜福利18| 久久精品成人免费网站| 国产激情欧美一区二区| 精品福利观看| 1024香蕉在线观看| 女警被强在线播放| 国产精华一区二区三区| 天堂影院成人在线观看| 波多野结衣高清无吗| 欧美成人一区二区免费高清观看 | 91精品国产国语对白视频| 男女下面进入的视频免费午夜 | 一二三四社区在线视频社区8| 在线av久久热| 一区二区三区激情视频| 亚洲九九香蕉| 黄色a级毛片大全视频| 村上凉子中文字幕在线| 禁无遮挡网站| 欧美 亚洲 国产 日韩一| 亚洲精品美女久久久久99蜜臀| 成年女人毛片免费观看观看9| 亚洲 欧美 日韩 在线 免费| 久久中文字幕一级| 人妻丰满熟妇av一区二区三区| 美女午夜性视频免费| 午夜精品在线福利| 欧美乱码精品一区二区三区| 亚洲视频免费观看视频| 久久草成人影院| 午夜视频精品福利| 日本 欧美在线| av电影中文网址| 免费av毛片视频| 国产亚洲欧美精品永久| 黄色 视频免费看| 国产精品久久久久久人妻精品电影| 一级毛片女人18水好多| 在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品合色在线| 亚洲久久久国产精品| 亚洲一区高清亚洲精品| 一级黄色大片毛片| 久久中文字幕人妻熟女| 两个人看的免费小视频| 美女大奶头视频| 在线观看日韩欧美| 国产精品亚洲一级av第二区| 淫秽高清视频在线观看| 国内久久婷婷六月综合欲色啪| 天天躁夜夜躁狠狠躁躁| 夜夜爽天天搞| 搡老熟女国产l中国老女人| 宅男免费午夜| 大码成人一级视频| 精品一区二区三区视频在线观看免费| 亚洲狠狠婷婷综合久久图片| 九色国产91popny在线| 巨乳人妻的诱惑在线观看| 99香蕉大伊视频| 97人妻天天添夜夜摸| 女生性感内裤真人,穿戴方法视频| 自拍欧美九色日韩亚洲蝌蚪91| 日本免费一区二区三区高清不卡 | 亚洲av成人一区二区三| 国产91精品成人一区二区三区| e午夜精品久久久久久久| 色哟哟哟哟哟哟| 一个人观看的视频www高清免费观看 | 欧美黄色片欧美黄色片| 狂野欧美激情性xxxx| 日韩精品青青久久久久久| 久久久久九九精品影院| 不卡av一区二区三区| 成人亚洲精品一区在线观看| 国产日韩一区二区三区精品不卡| 国产xxxxx性猛交| 夜夜夜夜夜久久久久| 欧美一区二区精品小视频在线| 日韩高清综合在线| 精品福利观看| 黄频高清免费视频| netflix在线观看网站| 国产欧美日韩精品亚洲av| 久久午夜综合久久蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| 麻豆av在线久日| 色婷婷久久久亚洲欧美| 黄色视频,在线免费观看| 精品国产乱子伦一区二区三区| 黄片小视频在线播放| √禁漫天堂资源中文www| 国产又爽黄色视频| 国产麻豆69| 欧美日韩黄片免| 桃红色精品国产亚洲av| 9热在线视频观看99| 叶爱在线成人免费视频播放| 久久久久久久久中文| 欧美黑人欧美精品刺激| 欧美中文综合在线视频| 国产三级在线视频| 很黄的视频免费| 精品国产亚洲在线| 18美女黄网站色大片免费观看| 少妇熟女aⅴ在线视频| 免费观看人在逋| 丝袜美足系列| 国产亚洲欧美在线一区二区| 欧美一级a爱片免费观看看 | 中文字幕最新亚洲高清| 一a级毛片在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲精品中文字幕在线视频| 亚洲人成电影观看| 欧美成人午夜精品| 久久精品91无色码中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 久久中文字幕一级| 国产亚洲精品综合一区在线观看 | 美女扒开内裤让男人捅视频| 精品一品国产午夜福利视频| 精品日产1卡2卡| 亚洲美女黄片视频| 中文字幕人妻熟女乱码| 亚洲精品中文字幕在线视频| 免费人成视频x8x8入口观看| av电影中文网址| 亚洲无线在线观看| 波多野结衣巨乳人妻| 亚洲精品美女久久av网站| 国产亚洲欧美精品永久| 色在线成人网| 欧美黄色淫秽网站| 波多野结衣高清无吗| 99香蕉大伊视频| 日本撒尿小便嘘嘘汇集6| 亚洲一区中文字幕在线| 丝袜美腿诱惑在线| 亚洲av片天天在线观看| 女性生殖器流出的白浆| 国产成人免费无遮挡视频| 老汉色av国产亚洲站长工具| 最近最新中文字幕大全电影3 | 欧美激情久久久久久爽电影 | 亚洲欧美日韩高清在线视频| 亚洲,欧美精品.| 色尼玛亚洲综合影院| 婷婷六月久久综合丁香| 午夜成年电影在线免费观看| 色在线成人网| 两人在一起打扑克的视频| 亚洲国产毛片av蜜桃av| 亚洲免费av在线视频| av天堂在线播放| x7x7x7水蜜桃| 此物有八面人人有两片| 三级毛片av免费| 97碰自拍视频| 久久久久久亚洲精品国产蜜桃av| 人人妻,人人澡人人爽秒播| 女警被强在线播放| 国产亚洲精品第一综合不卡| 久久草成人影院| 欧美在线黄色|