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

    結(jié)構(gòu)風(fēng)管通風(fēng)分析與優(yōu)化研究

    2016-12-28 07:56:56姚正鋼倪崇本
    造船技術(shù) 2016年5期
    關(guān)鍵詞:機(jī)艙風(fēng)管管內(nèi)

    姚正鋼, 倪崇本, 姚 雯

    (1.江蘇科技大學(xué) 船舶與海洋工程學(xué)院, 江蘇 鎮(zhèn)江 212003; 2.上海船舶工藝研究所, 上海 200032)

    ?

    結(jié)構(gòu)風(fēng)管通風(fēng)分析與優(yōu)化研究

    姚正鋼1, 倪崇本2, 姚 雯2

    (1.江蘇科技大學(xué) 船舶與海洋工程學(xué)院, 江蘇 鎮(zhèn)江 212003; 2.上海船舶工藝研究所, 上海 200032)

    主要選取海洋工程平臺(tái)內(nèi)機(jī)艙和變壓器間的結(jié)構(gòu)風(fēng)管作為載體,運(yùn)用HyperMesh和Fluent數(shù)值模擬軟件作為工具,對(duì)通風(fēng)管道系統(tǒng)進(jìn)行流場(chǎng)分析,繪制管道阻力曲線和風(fēng)管內(nèi)速度矢量,對(duì)結(jié)構(gòu)風(fēng)管進(jìn)行優(yōu)化,并與優(yōu)化前的管道通風(fēng)系統(tǒng)進(jìn)行對(duì)比,以此來(lái)提升通風(fēng)效率。通過(guò)數(shù)值模擬可以為通風(fēng)系統(tǒng)的設(shè)計(jì)提供依據(jù)和參考。

    計(jì)算流體力學(xué) 氣流組織 數(shù)值模擬 優(yōu)化

    0 引言

    通風(fēng)管道就是把符合衛(wèi)生標(biāo)準(zhǔn)的新鮮空氣輸送到各艙室需要的地點(diǎn),把艙室內(nèi)局部地區(qū)或設(shè)備散發(fā)的污濁、有害氣體直接排送到艙外或經(jīng)凈化處理后排送到艙外的管道。管道系統(tǒng)包括通風(fēng)除塵管道、空調(diào)管道等,作用是把通風(fēng)進(jìn)氣口、空氣的熱、濕及凈化處理設(shè)備、送(排)氣口、部件和風(fēng)機(jī)連成一個(gè)整體,使之有效運(yùn)轉(zhuǎn)[1]。通風(fēng)管道系統(tǒng)設(shè)計(jì)的目的是,在滿足規(guī)范要求和保證使用效果的前提下,合理組織氣體流動(dòng),實(shí)現(xiàn)管道阻力最小,使系統(tǒng)制造和日常運(yùn)行維護(hù)費(fèi)用最優(yōu)。

    船舶、海工結(jié)構(gòu)風(fēng)管,是一種附著在船舶結(jié)構(gòu)上的通風(fēng)管道,結(jié)構(gòu)風(fēng)管作為通風(fēng)系統(tǒng)的主通道,不但可以節(jié)省空間,而且加強(qiáng)了結(jié)構(gòu)強(qiáng)度。結(jié)構(gòu)風(fēng)管大部分為矩形,較方便制作,不同的區(qū)域?qū)τ诮Y(jié)構(gòu)風(fēng)管的形狀、厚度有不同的要求,它的作用主要是通風(fēng)換氣。船舶機(jī)艙通風(fēng)的目的是降低機(jī)艙溫度,排除機(jī)艙內(nèi)的油氣和水蒸氣,向機(jī)艙內(nèi)供應(yīng)充足的新鮮空氣,從而保證動(dòng)力裝置的正常工作以及改善機(jī)艙管理人員的工作條件和衛(wèi)生條件。計(jì)算流體力學(xué)(簡(jiǎn)稱CFD),是用電子計(jì)算機(jī)和離散化的數(shù)值方法對(duì)流體力學(xué)問(wèn)題進(jìn)行數(shù)值模擬和分析的一個(gè)新分支,應(yīng)用領(lǐng)域非常廣泛。隨著CFD技術(shù)的發(fā)展,已經(jīng)在模擬通風(fēng)分析方面實(shí)現(xiàn)了應(yīng)用,通過(guò)對(duì)船舶的通風(fēng)系統(tǒng)進(jìn)行數(shù)值模擬,分析結(jié)構(gòu)風(fēng)管內(nèi)的氣流組織形式,提出氣流組織的優(yōu)化方案。對(duì)結(jié)構(gòu)風(fēng)管建立流體數(shù)學(xué)模型和設(shè)置邊界條件,通過(guò)模擬的方法來(lái)實(shí)現(xiàn)對(duì)通風(fēng)系統(tǒng)的分析,為通風(fēng)系統(tǒng)的設(shè)計(jì)、優(yōu)化提供支持[2]。

    1 風(fēng)管流場(chǎng)建模技術(shù)分析

    結(jié)構(gòu)風(fēng)管是通風(fēng)系統(tǒng)的主通道,由風(fēng)機(jī)將外界新鮮空氣送到各艙室內(nèi),與艙室內(nèi)氣體進(jìn)行交換,排除有毒有害氣體,同時(shí)降低艙室內(nèi)的溫度[3]。模擬結(jié)構(gòu)風(fēng)管內(nèi)氣流組織,首先要建立結(jié)構(gòu)風(fēng)管的三維CAD模型。文中通過(guò)二維CAD圖紙,根據(jù)各種型材、板材之間的相對(duì)布置,在AUTOCAD軟件中直接三維建模(見(jiàn)圖1),建立了結(jié)構(gòu)風(fēng)管的簡(jiǎn)化三維模型。

    圖1 結(jié)構(gòu)風(fēng)管三維模型

    在此基礎(chǔ)上,根據(jù)劃分網(wǎng)格軟件HyperMesh要求對(duì)模型簡(jiǎn)化處理,根據(jù)以往劃分網(wǎng)格的經(jīng)驗(yàn),先選擇合適的網(wǎng)格類型、尺寸、數(shù)量等,然后對(duì)模型的出入口條件進(jìn)行設(shè)置,導(dǎo)入Fluent模擬計(jì)算軟件格式的計(jì)算模型,下面簡(jiǎn)單介紹劃網(wǎng)格過(guò)程。

    1.1 幾何導(dǎo)入、修復(fù)與清理

    從已有的機(jī)艙與變壓器間布置圖模型文件中抽取機(jī)艙的結(jié)構(gòu)風(fēng)管模型,對(duì)機(jī)艙風(fēng)管布置圖中多余的模型進(jìn)行刪除,導(dǎo)入到HyperMesh軟件當(dāng)中,在通常情況下我們需要簡(jiǎn)化模型,這樣對(duì)于得到更好的網(wǎng)格質(zhì)量十分重要。

    1.2 面網(wǎng)格劃分

    HyperMesh軟件的強(qiáng)大功能就是自動(dòng)劃分網(wǎng)格,通過(guò)設(shè)置網(wǎng)格的尺寸、類型能夠快速地在一個(gè)面或者多個(gè)面自動(dòng)生成網(wǎng)格。我們?cè)O(shè)置容差來(lái)檢查出不滿足最小尺寸的網(wǎng)格單元,觀察網(wǎng)格排列是否整齊,不規(guī)則的地方可以重新劃分。網(wǎng)格劃分的思路是由面網(wǎng)格引出體網(wǎng)格,由于面網(wǎng)格的質(zhì)量決定了體網(wǎng)格的質(zhì)量,因此面網(wǎng)格的質(zhì)量決定了總體網(wǎng)格的質(zhì)量。

    1.3 體網(wǎng)格劃分

    對(duì)于體網(wǎng)格的形式,我們采用的是四面體網(wǎng)格,原因是:對(duì)于復(fù)雜的模型,劃分體網(wǎng)格都是依靠自動(dòng)或半自動(dòng)網(wǎng)格生成器產(chǎn)生,一般的三維域不能總是被分解為規(guī)則的六面體網(wǎng)格或其它類型網(wǎng)格,但它可以隨時(shí)被劃分為規(guī)則的四面體網(wǎng)格;對(duì)于比較大的模型,四面體網(wǎng)格計(jì)算速度最快,減少了計(jì)算時(shí)間且存儲(chǔ)空間顯著;考慮到相等的精確性和時(shí)間等因素,我們選用四面體網(wǎng)格。

    1.4 網(wǎng)格區(qū)域的映射性

    我們創(chuàng)建實(shí)體單元,可以從節(jié)點(diǎn)、曲面、實(shí)體幾何創(chuàng)建,在實(shí)體網(wǎng)格劃分過(guò)程中,經(jīng)過(guò)編輯的實(shí)體需在某個(gè)方向或者三個(gè)方向上可映射,我們稱這種屬性為可映射性,即要保持網(wǎng)格之間的連續(xù)性。區(qū)域劃分的目的就是使用邊界線和切割線分割實(shí)體,使每一塊實(shí)體都成為可映射形式,在平面的二維網(wǎng)格指引下,生成三維實(shí)體網(wǎng)格,然后由一個(gè)實(shí)體映射到另一個(gè)實(shí)體,最后生成完整的有限元模型。

    1.5 單元質(zhì)量和檢查

    單元的質(zhì)量決定了計(jì)算結(jié)果的精度,通過(guò)一些指標(biāo)可以清楚地顯示出單元質(zhì)量的好壞。常用的單元質(zhì)量指標(biāo)為偏斜度、長(zhǎng)寬比、翹曲度、雅可比等。對(duì)于劃分的網(wǎng)格,大部分都是連續(xù)的,我們?cè)O(shè)定壁面條件為流場(chǎng)邊界。為了模擬流場(chǎng)的實(shí)際情況,設(shè)置擋板處和小肋板處的網(wǎng)格不連續(xù),為插入邊界,如圖2所示。

    圖2 網(wǎng)格壁面條件圖

    2 風(fēng)管數(shù)值模擬仿真與優(yōu)化

    結(jié)構(gòu)風(fēng)管的計(jì)算模型導(dǎo)入軟件后,對(duì)計(jì)算結(jié)果進(jìn)行處理,即可得到通風(fēng)管道內(nèi)的壓力、速度矢量圖。通過(guò)對(duì)優(yōu)化前后整個(gè)結(jié)構(gòu)風(fēng)管內(nèi)的速度、壓力分析,可以得到管道內(nèi)通風(fēng)量的變化,驗(yàn)證是否滿足各艙室設(shè)計(jì)要求的通風(fēng)量,為管道的通風(fēng)設(shè)計(jì)與優(yōu)化提供依據(jù)和參考。

    2.1 結(jié)構(gòu)風(fēng)管的模擬仿真

    2.1.1 創(chuàng)建計(jì)算模型

    通過(guò)提供的機(jī)艙與變壓器間布置圖與三維模型文件(dxf格式)為基礎(chǔ),從模型文件中抽取機(jī)艙與變壓器間的風(fēng)管模型,導(dǎo)入HyperMesh軟件中,并對(duì)風(fēng)管內(nèi)部的流場(chǎng)區(qū)域進(jìn)行網(wǎng)格劃分,如圖3所示。

    圖3 風(fēng)管內(nèi)流場(chǎng)網(wǎng)格模型

    2.1.2 設(shè)定物理模型

    風(fēng)機(jī)風(fēng)量采用的是360 000 m3/h,按照矩形風(fēng)管的空氣流速公式:

    V=L/AB

    式中:V為流速;L為流量;A、B分別為長(zhǎng)和寬。對(duì)風(fēng)管截面面積估算后,風(fēng)管內(nèi)最大風(fēng)速約為50 m/s,約為0.16馬赫。按歐拉方程分析氣體流動(dòng)中密度變化:

    式中:Ma為馬赫數(shù);p、v分別是密度和速度;若Ma<0.3,密度變化可以忽略,屬于不可壓縮流動(dòng)范疇,視為無(wú)壓縮氣體流動(dòng)。氣體在管道內(nèi)流動(dòng)時(shí),流動(dòng)狀態(tài)分為層流和湍流,流體質(zhì)點(diǎn)的軌跡中規(guī)則的光滑曲線為層流,不規(guī)則運(yùn)動(dòng)的是湍流。根據(jù)雷諾數(shù)公式:

    Re=ρvd/η

    式中:ρ為空氣密度;v為流速;d為管道等效直徑;η為空氣運(yùn)動(dòng)粘性,取0.133×10-5Pa·s,通過(guò)計(jì)算得出雷諾數(shù)為3.6241×107,臨界雷諾數(shù)一般取值為2 000~4 000,當(dāng)Re遠(yuǎn)大于4 000,認(rèn)為氣流為充分發(fā)展的湍流運(yùn)動(dòng),而且需要考慮產(chǎn)生的渦漩,我們采用目前廣泛使用的Realizable k-ε湍流模型。

    +Gk-Gb+ρ ε-Ym+Sk

    式中:Gk和Gb分別是由平均速度和浮力所引起的湍流動(dòng)能k產(chǎn)生項(xiàng)。在船舶與海工結(jié)構(gòu)風(fēng)管中不存在隨時(shí)間變化邊界條件等因素,因此選用定常流模型。

    2.1.3 設(shè)置邊界條件

    將網(wǎng)格轉(zhuǎn)換為計(jì)算流體力學(xué)軟件Fluent的輸入格式后,對(duì)風(fēng)管內(nèi)流場(chǎng)進(jìn)行分析。風(fēng)管內(nèi)流場(chǎng)邊界條件如圖4所示,進(jìn)風(fēng)口為質(zhì)量流入口,出風(fēng)口為壓力出口。根據(jù)設(shè)計(jì)信息,變壓器間排風(fēng)扇靜壓為300 Pa,設(shè)定變壓器間出風(fēng)口與機(jī)艙出風(fēng)口壓差為300 Pa,即

    P機(jī)艙-P變壓器間= 300 Pa

    圖4 計(jì)算模型邊界條件

    設(shè)置風(fēng)管管壁作為壁面,取常規(guī)鍍鋅板的粗糙高度0.2 mm。風(fēng)機(jī)排風(fēng)口為速度邊界,取額定風(fēng)量360 000 m3/h,將速度方向設(shè)定為Z軸正方向流出。進(jìn)氣口邊界壓強(qiáng)為大氣壓,由于邊壁或流量的變化,均勻流在局部地區(qū)遇到干擾,使氣流的大小、方向和分布產(chǎn)生變化,甚至產(chǎn)生渦流,由此產(chǎn)生局部損失[5]。其局部流動(dòng)阻力表述為

    式中:ζ為局部阻力系數(shù);v為局部空氣流速,m/s。根據(jù)廠商提供的阻力系數(shù);ζ排風(fēng)=2.0,空氣=1.24 kg/m3,擬合后的局部流動(dòng)阻力為

    z=1.24·v2

    減小局部阻力著眼于防止或推遲氣流與壁面的分離,避免漩渦區(qū)的產(chǎn)生或減小漩渦區(qū)的大小和強(qiáng)度[4]。把計(jì)算的阻力結(jié)果輸入邊界條件,不對(duì)流速或壓力進(jìn)行假定,可以更加真實(shí)地反映實(shí)際流動(dòng)情況??刂品匠滩捎肧IMPLE方法求解,壓力項(xiàng)采用標(biāo)準(zhǔn)離散格式,動(dòng)量項(xiàng)、對(duì)流項(xiàng)和擴(kuò)散項(xiàng)均使用二階迎風(fēng)格式離散。進(jìn)行迭代計(jì)算后流場(chǎng)趨于穩(wěn)定,然后對(duì)流場(chǎng)內(nèi)流速較高的網(wǎng)格,以及速度梯度較大的網(wǎng)格進(jìn)行自適應(yīng)細(xì)分,以反映流場(chǎng)的細(xì)節(jié)。網(wǎng)格自適應(yīng)后再次進(jìn)行迭代計(jì)算,直至計(jì)算收斂。

    2.1.4 管路阻力與實(shí)際風(fēng)量計(jì)算

    按照計(jì)算實(shí)際風(fēng)量的需求,計(jì)算管路阻力以繪制管路阻力ptf曲線。共取4個(gè)工況,分別為通風(fēng)量30×104m3/h、32×104m3/h、34×104m3/h和36×104m3/h,對(duì)應(yīng)的管路阻力計(jì)算結(jié)果見(jiàn)表1。

    表1 管路阻力計(jì)算結(jié)果

    兩個(gè)進(jìn)風(fēng)口各配一臺(tái)風(fēng)機(jī) (Novenco S31 A-B),廠家提供的風(fēng)機(jī)曲線如圖5所示。

    圖5 Novenco S31 A-B風(fēng)機(jī)曲線

    風(fēng)機(jī)葉片安裝角為58°,擬合風(fēng)機(jī)曲線和管路阻力曲線后得出單臺(tái)風(fēng)機(jī)實(shí)際氣體流量應(yīng)為qv= 47.4 m3/s,管路阻力ptf = 1 820 Pa,如圖6所示。

    圖6 管路阻力與實(shí)際風(fēng)量計(jì)算

    運(yùn)轉(zhuǎn)點(diǎn)就是性能曲線和阻力曲線的交點(diǎn),因?yàn)檫@一點(diǎn)上的動(dòng)力矩和阻力矩相等,可以使風(fēng)機(jī)可持續(xù)勻速運(yùn)轉(zhuǎn)。如果阻力大于額定動(dòng)力,則風(fēng)機(jī)功率偏小,導(dǎo)致帶不動(dòng)或超載運(yùn)行而燒壞電機(jī),如果阻力小于額定動(dòng)力,則風(fēng)機(jī)輕載運(yùn)行,浪費(fèi)了一部分功率。通過(guò)圖6我們可以看出,風(fēng)機(jī)的運(yùn)轉(zhuǎn)點(diǎn)就是流量為qv= 47.4 m3/s所對(duì)應(yīng)的點(diǎn),折算后機(jī)艙與變壓器間實(shí)際通風(fēng)量約為341 280 m3/h。

    2.1.5 風(fēng)量分配計(jì)算

    按照實(shí)際通風(fēng)量設(shè)置流場(chǎng)邊界條件,利用Fluent軟件計(jì)算機(jī)艙與變壓器間各出風(fēng)口的送風(fēng)量分配,結(jié)果見(jiàn)表2。

    表2 機(jī)艙與變壓器間送風(fēng)量分配

    根據(jù)計(jì)算結(jié)果,大約僅10%的通風(fēng)量分配到變壓器間。對(duì)風(fēng)管內(nèi)流場(chǎng)進(jìn)行分析,風(fēng)管內(nèi)壓力分布如圖7所示。圖7中右側(cè)風(fēng)管的直角彎處壓力損失較大(見(jiàn)畫(huà)圈處)。

    圖7 風(fēng)管內(nèi)壓力分布

    圖8(a)為結(jié)構(gòu)風(fēng)管與變壓器間風(fēng)管連接處速度矢量,顯示大部分氣體直接從出風(fēng)口進(jìn)入機(jī)艙,僅少量氣體進(jìn)入變壓器間風(fēng)管??梢钥紤]加裝導(dǎo)流板以改善變壓器間通風(fēng)狀況。圖8(b)為變壓器間風(fēng)管內(nèi)速度矢量圖,顯示變壓器間風(fēng)管中間段內(nèi)基本無(wú)氣體流動(dòng),是左右兩側(cè)風(fēng)機(jī)通風(fēng)的分界點(diǎn)。該位置靠近左側(cè)風(fēng)機(jī),右側(cè)風(fēng)機(jī)送風(fēng)路程較長(zhǎng),沿程損耗較大。

    圖8 結(jié)構(gòu)風(fēng)管內(nèi)速度矢量圖

    速度矢量圖可以直觀地表示速度的大小、方向和渦流情況,模擬了氣體在管道中的流動(dòng)狀態(tài)。網(wǎng)格單元是最小的分析實(shí)體,在網(wǎng)格的中心用剪頭表示矢量,用顏色和長(zhǎng)度表示大小。因此選用速度矢量圖來(lái)對(duì)結(jié)構(gòu)風(fēng)管內(nèi)氣流組織進(jìn)行分析。通過(guò)機(jī)艙與變壓器間模型和風(fēng)機(jī)參數(shù)為基礎(chǔ),進(jìn)行通風(fēng)流場(chǎng)的計(jì)算與分析,得出機(jī)艙與變壓器間實(shí)際通風(fēng)量計(jì)算結(jié)果為341 280 m3/h,其中變壓器間30 498 m3/h,機(jī)艙310 782 m3/h。我們分析認(rèn)為變壓器間與機(jī)艙通風(fēng)量分配不均問(wèn)題較為突出,可采取的改進(jìn)措施如下:

    (1) 加裝導(dǎo)流板,增加變壓器間通風(fēng)量[6];

    (2) 選用其他型號(hào)的變壓器間排風(fēng)扇,加大機(jī)艙與變壓器間壓差(當(dāng)前為300 Pa);

    (3) 縮短變壓器間送風(fēng)長(zhǎng)度,可考慮將直角彎改為兩個(gè)45°彎、風(fēng)口位置前移縮短與風(fēng)機(jī)距離。

    2.2 結(jié)構(gòu)風(fēng)管的優(yōu)化

    2.2.1 安裝導(dǎo)流板

    導(dǎo)流板設(shè)計(jì)是為了減小空氣阻力??諝饬鲃?dòng)的阻力有兩種,一種是因?yàn)榭諝獗旧淼恼硿耘c管壁的摩擦而產(chǎn)生的阻力,稱為摩擦阻力;另一種是空氣流經(jīng)管件及擋板處,因?yàn)榱魉俚拇笮『头较蜃兓约爱a(chǎn)生渦流而產(chǎn)生的阻力,為局部阻力。如圖9所示,在擋板處氣流速度太大,擋板區(qū)域速度達(dá)到了40 m/s,而其它區(qū)域?yàn)?0 m/s,根據(jù)管道阻力公式,為了減小摩擦阻力,我們需要降低氣體的流速。

    式中:l為局部阻力當(dāng)量長(zhǎng)度,m;ζ為局部阻力系數(shù);λ為摩擦系數(shù);Δpf為阻力損失,Pa。

    在常壓下,氣體管道的流速一般為10 m/s~30 m/s,流速越大,阻力越大。另一方面,從圖9可以看出,在擋板處有小渦漩,產(chǎn)生了繞流,對(duì)氣體的通過(guò)產(chǎn)生了阻力,需要消除。

    圖9 結(jié)構(gòu)風(fēng)管內(nèi)繞流圖

    通過(guò)上述分析,為了減小風(fēng)管阻力,需要進(jìn)行導(dǎo)流板的設(shè)計(jì),以現(xiàn)有的模型為基礎(chǔ),根據(jù)先前流場(chǎng)分析結(jié)果,提出機(jī)艙與變壓器間的風(fēng)管內(nèi)導(dǎo)流板安裝方案(見(jiàn)圖10)。風(fēng)管內(nèi)共安裝四塊導(dǎo)流板,導(dǎo)流板A位于進(jìn)風(fēng)口下方遮住突出的結(jié)構(gòu),導(dǎo)流板B位于結(jié)構(gòu)風(fēng)管內(nèi)遮住小肋骨,導(dǎo)流板C、D分別位于左右節(jié)流孔下方遮住凹陷的結(jié)構(gòu)。

    圖10 導(dǎo)流板安裝方案圖

    2.2.2 優(yōu)化后模擬仿真

    對(duì)安裝導(dǎo)流板后風(fēng)管內(nèi)流場(chǎng)進(jìn)行建模,重新用Hypermesh軟件劃分網(wǎng)格,輸入Fluent軟件進(jìn)行流場(chǎng)分析。風(fēng)管內(nèi)流場(chǎng)邊界條件設(shè)置進(jìn)風(fēng)口為質(zhì)量流入口,出風(fēng)口為壓力出口,設(shè)定變壓器間出風(fēng)口與機(jī)艙出風(fēng)口壓差為300 Pa。

    2.2.3 優(yōu)化后效果分析

    通過(guò)優(yōu)化前后的速度矢量對(duì)比,得出通風(fēng)量的變化??芍?,優(yōu)化后減小了管道的阻力,有更多新鮮的空氣進(jìn)入各個(gè)艙室,且優(yōu)化后的系統(tǒng)可以合理地分配各出風(fēng)口的風(fēng)量,使管道內(nèi)氣體分布更合理。對(duì)模擬數(shù)據(jù)進(jìn)行處理后,繪制管路阻力ptf曲線,并與未安裝導(dǎo)流板的管路阻力進(jìn)行對(duì)比。對(duì)比結(jié)果如表3所示。

    表3 管路阻力對(duì)比

    擬合風(fēng)機(jī)曲線和管路阻力曲線后對(duì)比導(dǎo)流板的效果,未裝導(dǎo)流板時(shí),單臺(tái)風(fēng)機(jī)氣體流量估算值為qv0= 47.4 m3/s,管路阻力ptf0= 1 820 Pa;安裝導(dǎo)流板后,單臺(tái)風(fēng)機(jī)氣體流量估算值qv= 54.0 m3/s,管路阻力ptf = 1 707 Pa,如圖11所示。根據(jù)分析結(jié)果,在不改變其他配置情況下,安裝導(dǎo)流板大約可以提升約1/8的風(fēng)量,折算安裝導(dǎo)流板后機(jī)艙與變壓器間風(fēng)量預(yù)計(jì)能達(dá)到388 000 m3/h。

    圖11 風(fēng)量與管路阻力對(duì)比

    同樣通風(fēng)量的工況下(360 000 m3/h)風(fēng)管內(nèi)流場(chǎng)進(jìn)行對(duì)比,未裝導(dǎo)流板風(fēng)管內(nèi)局部最大風(fēng)速為62.2 m/s,安裝導(dǎo)流板后風(fēng)管內(nèi)局部最大風(fēng)速將為52.0 m/s,說(shuō)明導(dǎo)流板可以降低最大風(fēng)速以減少管路內(nèi)流動(dòng)阻力。圖12、圖13為導(dǎo)流板前后風(fēng)管內(nèi)速度矢量對(duì)比,從圖上可以看出,導(dǎo)流板對(duì)通風(fēng)口下方擋板處氣流改善效果明顯。

    圖12 風(fēng)管內(nèi)速度矢量分布對(duì)比

    圖13 擋板處速度矢量分布對(duì)比

    通過(guò)機(jī)艙與變壓器間風(fēng)管模型和風(fēng)機(jī)參數(shù)為基礎(chǔ),進(jìn)行導(dǎo)流板方案設(shè)計(jì)與分析。CFD對(duì)比分析表明,安裝導(dǎo)流板大約可以提升約1/8的風(fēng)量,安裝導(dǎo)流板后機(jī)艙與變壓器間風(fēng)量估算值為388 000 m3/h,通過(guò)對(duì)結(jié)構(gòu)風(fēng)管的優(yōu)化,減小了管道內(nèi)的阻力,使出風(fēng)量增加,管道內(nèi)的氣流組織分布更加合理、均勻、穩(wěn)定,同時(shí)也能夠降低船體的震動(dòng)和噪音。

    對(duì)于該方案的計(jì)算與優(yōu)化,得到船廠與船東的認(rèn)可,并已采用該方案。

    3 結(jié)論

    本文采用CFD技術(shù)對(duì)船舶與海工產(chǎn)品的機(jī)艙、變壓器間的結(jié)構(gòu)風(fēng)管進(jìn)行通風(fēng)計(jì)算,本文主要針對(duì)船舶空通風(fēng)專業(yè)進(jìn)行CFD模型的研究與選取,風(fēng)管計(jì)算工作的網(wǎng)格可靠性驗(yàn)證,管路性能與風(fēng)機(jī)性能的匹配,風(fēng)管優(yōu)化方案設(shè)計(jì)等進(jìn)行研究,通過(guò)對(duì)通風(fēng)管道系統(tǒng)的設(shè)計(jì)計(jì)算,借鑒有限體積法,運(yùn)動(dòng)代數(shù)法生成矩形風(fēng)管計(jì)算的結(jié)構(gòu)化網(wǎng)格,建立控制方程并

    [][]進(jìn)行離散,對(duì)k-ε湍流求解直至計(jì)算結(jié)果收斂,對(duì)優(yōu)化后的結(jié)構(gòu)風(fēng)管質(zhì)量流分析,得到了結(jié)構(gòu)風(fēng)管通風(fēng)的優(yōu)化方案,取得了階段性的進(jìn)展,為后續(xù)船舶與海工平臺(tái)艙室通風(fēng)的設(shè)計(jì)提供了指導(dǎo)性的參考和技術(shù)支持。

    [1] 金小闖,蔣合葉,黃林峰.船舶艙室通風(fēng)設(shè)計(jì)與分析[J].機(jī)電信息,2013(24):154-155.

    [2] 倪崇本,朱航,丁金鴻,等.LNG船機(jī)艙內(nèi)通風(fēng)氣流組織的數(shù)值分析[J].船舶工程,2014(3):57-61.

    [3] 郭昂,封海寶,王馳明.某海監(jiān)船機(jī)艙通風(fēng)系統(tǒng)的設(shè)計(jì)和優(yōu)化方案[J].江蘇船舶,2014(3):19-21.

    [4] 中國(guó)船舶工業(yè)集團(tuán)公司.船舶設(shè)計(jì)使用手冊(cè)-輪機(jī)分冊(cè)[M].北京:國(guó)防工業(yè)出版社,2013.

    [5] 羅繼杰,張蔚東.《通風(fēng)管道沿程阻力計(jì)算選用表》編制及設(shè)計(jì)應(yīng)用[J].暖通空調(diào),2009(9):97-99.

    [6] 張衛(wèi)東,王世忠,王新海.船舶艙室空調(diào)送風(fēng)優(yōu)化技術(shù)[J].艦船科學(xué)技術(shù),2011, 33(2):121-125.

    Research on Air Distribution and Optimization of Structural Air Duct

    YAO Zheng-gang1, NI Chong-ben2, YAO Wen2

    (1.School of Naval Architecture and Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang Jiangsu 212003, China;2.Shanghai Shipbuilding Technology Research Institute, Shanghai 200032, China)

    The marine platform structural air duct between engine room and transformer room was mainly selected as a carrier, and the simulation software of HyperMesh and Fluent as a tool was used to analyze ventilation duct system flow field. Getting the pipe resistance curve and wind tube velocity vector. And then, the structure of the wind tube was optimized and ventilation systems were compared in order to improve the ventilation efficiency. It can provide analysis of the numerical simulation for the design of the ventilation system.

    CFD Air distribution Numerical simulation Optimization

    姚正鋼(1987-),男,工程師,研究方向?yàn)榇棒秆b設(shè)計(jì)應(yīng)用。

    P75

    A

    猜你喜歡
    機(jī)艙風(fēng)管管內(nèi)
    一種橡膠煉膠設(shè)備的冷卻裝置
    船舶機(jī)艙火災(zāi)的原因分析及預(yù)防
    水上消防(2022年2期)2022-07-22 08:43:56
    船舶機(jī)艙通風(fēng)相關(guān)要求及常見(jiàn)關(guān)閉裝置分析
    水上消防(2021年3期)2021-08-21 03:12:20
    船舶機(jī)艙常見(jiàn)消防隱患及防控
    水上消防(2020年5期)2020-12-14 07:16:20
    12000t/d水泥熟料生產(chǎn)線三次風(fēng)管的優(yōu)化改造*
    水泥工程(2020年3期)2020-10-28 01:45:26
    上傾管內(nèi)油水兩相流流型實(shí)驗(yàn)研究
    揚(yáng)礦管內(nèi)高速螺旋流的數(shù)值模擬與仿真
    機(jī)艙污水井應(yīng)急除油設(shè)計(jì)
    腕管內(nèi)原發(fā)性腱鞘結(jié)核誤診1例
    微細(xì)管內(nèi)CO2流動(dòng)沸騰換熱特性研究
    大香蕉久久成人网| 国产一区二区三区综合在线观看| 午夜福利在线免费观看网站| 无限看片的www在线观看| svipshipincom国产片| 日韩中文字幕欧美一区二区| 十八禁人妻一区二区| 在线亚洲精品国产二区图片欧美| 成年人午夜在线观看视频| 脱女人内裤的视频| 性少妇av在线| 国产xxxxx性猛交| 亚洲中文日韩欧美视频| 18禁黄网站禁片午夜丰满| 国产精品久久久久久人妻精品电影 | 丝瓜视频免费看黄片| 久久精品国产综合久久久| 高清在线国产一区| 亚洲伊人久久精品综合| 桃红色精品国产亚洲av| 久久精品人人爽人人爽视色| 久久久欧美国产精品| 国产亚洲av片在线观看秒播厂| 超碰成人久久| 成人国产av品久久久| 久久久久精品国产欧美久久久 | 人人妻人人爽人人添夜夜欢视频| 黄片小视频在线播放| 男女午夜视频在线观看| 国产高清videossex| 亚洲精品日韩在线中文字幕| 欧美黑人欧美精品刺激| 亚洲自偷自拍图片 自拍| 可以免费在线观看a视频的电影网站| 国产精品成人在线| av有码第一页| av在线老鸭窝| 999久久久精品免费观看国产| 欧美 亚洲 国产 日韩一| 日本av手机在线免费观看| 男女高潮啪啪啪动态图| 最新在线观看一区二区三区| 国产av又大| 母亲3免费完整高清在线观看| 高清黄色对白视频在线免费看| 超碰97精品在线观看| 50天的宝宝边吃奶边哭怎么回事| 高清黄色对白视频在线免费看| 天天影视国产精品| 国产精品熟女久久久久浪| 中文字幕最新亚洲高清| 热99re8久久精品国产| 国产精品久久久久久人妻精品电影 | 国产高清videossex| 99香蕉大伊视频| 国内毛片毛片毛片毛片毛片| 日韩欧美国产一区二区入口| 一级a爱视频在线免费观看| 纵有疾风起免费观看全集完整版| 国产在线观看jvid| 精品亚洲成a人片在线观看| 亚洲国产欧美日韩在线播放| 在线观看www视频免费| 久久国产精品影院| 少妇 在线观看| 日韩 亚洲 欧美在线| 精品一区在线观看国产| 免费人妻精品一区二区三区视频| 国产精品国产av在线观看| 日韩精品免费视频一区二区三区| 亚洲国产看品久久| 又紧又爽又黄一区二区| 91精品国产国语对白视频| 亚洲av男天堂| 中亚洲国语对白在线视频| 啪啪无遮挡十八禁网站| 亚洲中文av在线| 在线观看免费视频网站a站| 国产一级毛片在线| 五月天丁香电影| 女性被躁到高潮视频| 热re99久久国产66热| 啦啦啦免费观看视频1| 亚洲av美国av| 亚洲中文av在线| 亚洲男人天堂网一区| 国产无遮挡羞羞视频在线观看| 久久精品亚洲熟妇少妇任你| 黑人猛操日本美女一级片| 老鸭窝网址在线观看| 国产精品成人在线| 午夜福利一区二区在线看| 一个人免费看片子| 啦啦啦视频在线资源免费观看| 久久综合国产亚洲精品| 亚洲精品久久午夜乱码| 一级毛片女人18水好多| 啦啦啦视频在线资源免费观看| 国产成人啪精品午夜网站| 国产区一区二久久| 少妇人妻久久综合中文| 一区二区av电影网| 成人亚洲精品一区在线观看| 天天添夜夜摸| 老鸭窝网址在线观看| 交换朋友夫妻互换小说| 日韩电影二区| 老司机靠b影院| 国产成人欧美在线观看 | 精品免费久久久久久久清纯 | 精品熟女少妇八av免费久了| 久久精品国产亚洲av高清一级| 午夜福利免费观看在线| 日本精品一区二区三区蜜桃| 国产av一区二区精品久久| 亚洲欧美色中文字幕在线| 18禁观看日本| 老司机深夜福利视频在线观看 | 成人av一区二区三区在线看 | 一本一本久久a久久精品综合妖精| 欧美日韩精品网址| 不卡av一区二区三区| 咕卡用的链子| 咕卡用的链子| av电影中文网址| 亚洲专区字幕在线| 天堂8中文在线网| 久久久久视频综合| 伊人亚洲综合成人网| 亚洲七黄色美女视频| av超薄肉色丝袜交足视频| 99国产精品一区二区蜜桃av | www.精华液| 久久影院123| 人人澡人人妻人| 久久久久久久精品精品| 大片免费播放器 马上看| 丰满少妇做爰视频| 精品福利观看| 真人做人爱边吃奶动态| 国产成人a∨麻豆精品| 搡老熟女国产l中国老女人| 国产老妇伦熟女老妇高清| 成人三级做爰电影| 亚洲欧洲日产国产| 最新的欧美精品一区二区| 操出白浆在线播放| 19禁男女啪啪无遮挡网站| 成年人午夜在线观看视频| 男人操女人黄网站| 悠悠久久av| 真人做人爱边吃奶动态| 精品福利永久在线观看| 午夜免费观看性视频| 精品国产超薄肉色丝袜足j| 久久影院123| 免费高清在线观看日韩| 午夜影院在线不卡| 国产精品久久久av美女十八| 久久精品熟女亚洲av麻豆精品| 国产日韩欧美亚洲二区| 大片免费播放器 马上看| 法律面前人人平等表现在哪些方面 | 青春草亚洲视频在线观看| 久久九九热精品免费| 精品免费久久久久久久清纯 | 99久久人妻综合| 国产成人欧美| 午夜两性在线视频| 国产精品1区2区在线观看. | 国产熟女午夜一区二区三区| 18在线观看网站| 免费在线观看黄色视频的| 亚洲国产欧美在线一区| 国产日韩一区二区三区精品不卡| 王馨瑶露胸无遮挡在线观看| 欧美日韩精品网址| 成人av一区二区三区在线看 | 一区在线观看完整版| 久久天躁狠狠躁夜夜2o2o| 女人久久www免费人成看片| 午夜福利,免费看| 视频区欧美日本亚洲| 麻豆av在线久日| 成人黄色视频免费在线看| 日韩制服骚丝袜av| 久久久国产一区二区| 老司机福利观看| 老司机午夜十八禁免费视频| 一本—道久久a久久精品蜜桃钙片| www.精华液| 亚洲色图综合在线观看| 人妻人人澡人人爽人人| 韩国高清视频一区二区三区| 国产高清视频在线播放一区 | 丝袜美腿诱惑在线| 人成视频在线观看免费观看| 国产黄色免费在线视频| 国产成人一区二区三区免费视频网站| 夜夜夜夜夜久久久久| 亚洲伊人久久精品综合| 纵有疾风起免费观看全集完整版| www.av在线官网国产| 宅男免费午夜| 超碰成人久久| 亚洲av片天天在线观看| av超薄肉色丝袜交足视频| 亚洲精品粉嫩美女一区| 中文字幕精品免费在线观看视频| 亚洲精品久久久久久婷婷小说| videosex国产| 久久久久久久大尺度免费视频| 一级毛片女人18水好多| 老司机福利观看| 精品熟女少妇八av免费久了| 欧美亚洲 丝袜 人妻 在线| 肉色欧美久久久久久久蜜桃| 久久国产精品男人的天堂亚洲| 啦啦啦啦在线视频资源| 国产精品久久久久久精品电影小说| 王馨瑶露胸无遮挡在线观看| 青青草视频在线视频观看| 亚洲专区字幕在线| 亚洲 国产 在线| 美女午夜性视频免费| 波多野结衣av一区二区av| 亚洲成人免费av在线播放| 日本91视频免费播放| 国产欧美日韩一区二区三 | 99re6热这里在线精品视频| 亚洲精品美女久久av网站| 美女高潮到喷水免费观看| 午夜福利免费观看在线| 国产福利在线免费观看视频| 亚洲人成电影观看| 日韩制服骚丝袜av| 亚洲欧洲日产国产| av国产精品久久久久影院| 女人爽到高潮嗷嗷叫在线视频| 免费在线观看视频国产中文字幕亚洲 | 一本大道久久a久久精品| 欧美精品啪啪一区二区三区 | 亚洲欧美精品综合一区二区三区| 男人爽女人下面视频在线观看| 韩国精品一区二区三区| 日韩欧美一区二区三区在线观看 | 久久久久久久久久久久大奶| 99久久国产精品久久久| 天堂8中文在线网| 国精品久久久久久国模美| 国产成人免费无遮挡视频| 欧美日韩亚洲国产一区二区在线观看 | bbb黄色大片| 99香蕉大伊视频| 99精品久久久久人妻精品| 欧美精品av麻豆av| av不卡在线播放| 纯流量卡能插随身wifi吗| 亚洲视频免费观看视频| 亚洲七黄色美女视频| 午夜老司机福利片| 精品一区二区三区av网在线观看 | 女警被强在线播放| 色视频在线一区二区三区| 老司机亚洲免费影院| 亚洲伊人色综图| 欧美黑人精品巨大| 男女无遮挡免费网站观看| 国产伦人伦偷精品视频| 久久人妻福利社区极品人妻图片| 久9热在线精品视频| 亚洲熟女精品中文字幕| 黑人欧美特级aaaaaa片| 亚洲 国产 在线| 老熟女久久久| 欧美黑人欧美精品刺激| 免费在线观看完整版高清| 国产成人精品久久二区二区91| 午夜激情久久久久久久| 岛国在线观看网站| 国产亚洲精品一区二区www | 欧美激情极品国产一区二区三区| 亚洲人成电影观看| 伦理电影免费视频| 女性被躁到高潮视频| 亚洲第一av免费看| 少妇裸体淫交视频免费看高清 | 正在播放国产对白刺激| 男女高潮啪啪啪动态图| 免费在线观看影片大全网站| 亚洲精品美女久久av网站| 欧美黄色淫秽网站| 国产1区2区3区精品| 中文字幕高清在线视频| 久久久久久免费高清国产稀缺| 97精品久久久久久久久久精品| 热re99久久国产66热| 中文字幕高清在线视频| 亚洲熟女毛片儿| 国产日韩欧美在线精品| 美女高潮到喷水免费观看| 丝袜人妻中文字幕| 大码成人一级视频| 国产成人免费无遮挡视频| 麻豆乱淫一区二区| 精品亚洲成a人片在线观看| 老汉色∧v一级毛片| 在线亚洲精品国产二区图片欧美| 欧美精品高潮呻吟av久久| 亚洲av日韩精品久久久久久密| 狠狠婷婷综合久久久久久88av| 后天国语完整版免费观看| 欧美日韩视频精品一区| 午夜成年电影在线免费观看| av不卡在线播放| av在线老鸭窝| 午夜福利在线免费观看网站| 日本wwww免费看| 精品卡一卡二卡四卡免费| 亚洲专区中文字幕在线| 女人久久www免费人成看片| 国产欧美日韩综合在线一区二区| 久久久精品94久久精品| 久久精品国产综合久久久| 国产精品久久久久久精品古装| 精品一区二区三卡| 免费观看av网站的网址| 欧美黄色淫秽网站| 亚洲精品国产av成人精品| 日本91视频免费播放| 亚洲成av片中文字幕在线观看| 老鸭窝网址在线观看| 欧美中文综合在线视频| 日韩 亚洲 欧美在线| h视频一区二区三区| 人成视频在线观看免费观看| 一级片'在线观看视频| 热99国产精品久久久久久7| 国产精品久久久av美女十八| 亚洲欧美一区二区三区黑人| av又黄又爽大尺度在线免费看| 热99国产精品久久久久久7| 美女午夜性视频免费| 交换朋友夫妻互换小说| 亚洲欧洲日产国产| 国产精品久久久久成人av| 国产日韩欧美视频二区| 国产一区有黄有色的免费视频| 大型av网站在线播放| 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| 国产老妇伦熟女老妇高清| 丁香六月天网| 大香蕉久久网| 免费看十八禁软件| 一级黄色大片毛片| 亚洲男人天堂网一区| 香蕉国产在线看| 99久久综合免费| 国产成人系列免费观看| 99国产精品一区二区蜜桃av | 午夜视频精品福利| 久久精品人人爽人人爽视色| 狂野欧美激情性xxxx| 丝袜脚勾引网站| 亚洲国产日韩一区二区| 午夜激情久久久久久久| 中文字幕制服av| 51午夜福利影视在线观看| 免费在线观看日本一区| 亚洲全国av大片| 亚洲精品粉嫩美女一区| 亚洲人成电影免费在线| 国产99久久九九免费精品| 亚洲国产成人一精品久久久| videosex国产| 视频区图区小说| 啦啦啦视频在线资源免费观看| 欧美日韩黄片免| 欧美黑人欧美精品刺激| 后天国语完整版免费观看| 免费在线观看日本一区| 久久久久久久大尺度免费视频| 日韩三级视频一区二区三区| 麻豆国产av国片精品| 精品少妇内射三级| 成人国产一区最新在线观看| 成在线人永久免费视频| av一本久久久久| 亚洲成人免费电影在线观看| 久热这里只有精品99| 18在线观看网站| 国产成人免费无遮挡视频| 精品亚洲乱码少妇综合久久| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面 | 黄色怎么调成土黄色| 国产深夜福利视频在线观看| 日韩欧美免费精品| 欧美国产精品一级二级三级| 婷婷色av中文字幕| 亚洲性夜色夜夜综合| 精品欧美一区二区三区在线| 又紧又爽又黄一区二区| 在线精品无人区一区二区三| 18禁裸乳无遮挡动漫免费视频| 久久国产精品影院| 国产在线视频一区二区| 女人精品久久久久毛片| 两性午夜刺激爽爽歪歪视频在线观看 | 波多野结衣av一区二区av| 香蕉丝袜av| 欧美少妇被猛烈插入视频| 日韩中文字幕欧美一区二区| 成人亚洲精品一区在线观看| 可以免费在线观看a视频的电影网站| 亚洲精品一卡2卡三卡4卡5卡 | 久久久久久久久久久久大奶| 成人18禁高潮啪啪吃奶动态图| 成人国产一区最新在线观看| 黄频高清免费视频| 多毛熟女@视频| 中文字幕制服av| 性色av乱码一区二区三区2| 国产精品一区二区在线不卡| 老熟女久久久| av有码第一页| 免费女性裸体啪啪无遮挡网站| 视频在线观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 侵犯人妻中文字幕一二三四区| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄色免费在线视频| 国产亚洲午夜精品一区二区久久| 黄色 视频免费看| 亚洲九九香蕉| 日本猛色少妇xxxxx猛交久久| 啦啦啦啦在线视频资源| 香蕉丝袜av| 久久影院123| 自线自在国产av| 性色av乱码一区二区三区2| av天堂久久9| 国产深夜福利视频在线观看| 美女大奶头黄色视频| 久久久久久久精品精品| 亚洲精品国产精品久久久不卡| 50天的宝宝边吃奶边哭怎么回事| 一级黄色大片毛片| 欧美激情极品国产一区二区三区| 首页视频小说图片口味搜索| 在线观看www视频免费| 黄色视频,在线免费观看| 免费久久久久久久精品成人欧美视频| 婷婷色av中文字幕| 成在线人永久免费视频| 脱女人内裤的视频| 亚洲视频免费观看视频| 一级片'在线观看视频| 另类精品久久| 女性被躁到高潮视频| 丁香六月天网| 99久久精品国产亚洲精品| 天堂俺去俺来也www色官网| 一级,二级,三级黄色视频| 一边摸一边抽搐一进一出视频| 久久女婷五月综合色啪小说| 高潮久久久久久久久久久不卡| 国产av一区二区精品久久| 日韩人妻精品一区2区三区| 女警被强在线播放| 亚洲一区中文字幕在线| 国产在线免费精品| 中文精品一卡2卡3卡4更新| 亚洲专区中文字幕在线| 视频在线观看一区二区三区| 午夜91福利影院| 中文字幕人妻丝袜一区二区| 一本色道久久久久久精品综合| 久久久久久人人人人人| 老司机影院毛片| 久久午夜综合久久蜜桃| 午夜福利视频在线观看免费| 亚洲色图 男人天堂 中文字幕| 欧美性长视频在线观看| 美女高潮喷水抽搐中文字幕| 国产精品亚洲av一区麻豆| 国产视频一区二区在线看| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 国产成人欧美| 国产成人精品久久二区二区91| 91精品伊人久久大香线蕉| 精品乱码久久久久久99久播| 日本猛色少妇xxxxx猛交久久| 国产一区有黄有色的免费视频| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美日韩高清在线视频 | 精品国产一区二区久久| 一级黄色大片毛片| 99国产综合亚洲精品| 亚洲国产精品一区三区| 老熟女久久久| 国产欧美日韩一区二区精品| 老司机在亚洲福利影院| 午夜福利免费观看在线| 久久ye,这里只有精品| 国产亚洲精品第一综合不卡| 最黄视频免费看| 18禁裸乳无遮挡动漫免费视频| 精品久久久久久电影网| 午夜激情久久久久久久| 欧美精品av麻豆av| 又黄又粗又硬又大视频| 新久久久久国产一级毛片| 成年人免费黄色播放视频| 国产精品自产拍在线观看55亚洲 | 波多野结衣av一区二区av| 国产av一区二区精品久久| 亚洲成人国产一区在线观看| 国产精品自产拍在线观看55亚洲 | 国产亚洲欧美在线一区二区| 亚洲精品国产色婷婷电影| 少妇精品久久久久久久| 18禁黄网站禁片午夜丰满| 久久人人爽人人片av| 亚洲中文av在线| 久久久久精品国产欧美久久久 | 亚洲 欧美一区二区三区| 建设人人有责人人尽责人人享有的| 国产精品麻豆人妻色哟哟久久| 久久国产精品大桥未久av| 欧美久久黑人一区二区| 性色av一级| 国内毛片毛片毛片毛片毛片| a级毛片在线看网站| 99久久国产精品久久久| 欧美另类一区| netflix在线观看网站| 久久国产精品男人的天堂亚洲| 久久九九热精品免费| 丰满人妻熟妇乱又伦精品不卡| 日韩中文字幕欧美一区二区| 电影成人av| 熟女少妇亚洲综合色aaa.| 亚洲va日本ⅴa欧美va伊人久久 | www.自偷自拍.com| 亚洲国产欧美一区二区综合| 大码成人一级视频| 久久久久网色| 午夜影院在线不卡| 精品国产乱码久久久久久男人| 欧美国产精品一级二级三级| videos熟女内射| 国产精品偷伦视频观看了| 热99国产精品久久久久久7| 国产主播在线观看一区二区| 久久久久久亚洲精品国产蜜桃av| 日韩欧美国产一区二区入口| 在线观看免费午夜福利视频| 不卡一级毛片| 欧美在线黄色| 久久精品国产亚洲av高清一级| 午夜福利视频精品| 波多野结衣av一区二区av| 在线永久观看黄色视频| 女人精品久久久久毛片| 无遮挡黄片免费观看| 老司机福利观看| 国产野战对白在线观看| av片东京热男人的天堂| 在线观看人妻少妇| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| 国产激情久久老熟女| 国产精品国产av在线观看| 国产人伦9x9x在线观看| 欧美黑人欧美精品刺激| 看免费av毛片| 性高湖久久久久久久久免费观看| 午夜久久久在线观看| 国产精品国产av在线观看| 亚洲av电影在线观看一区二区三区| 人人妻人人爽人人添夜夜欢视频| 高清在线国产一区| 国产一卡二卡三卡精品| 亚洲av电影在线观看一区二区三区| 精品人妻在线不人妻| 91大片在线观看| 青春草视频在线免费观看| 国产片内射在线| www.999成人在线观看| 亚洲情色 制服丝袜| 19禁男女啪啪无遮挡网站| 午夜91福利影院| 亚洲九九香蕉| 亚洲精品国产av成人精品| 日韩欧美免费精品| 老司机影院成人| 午夜福利在线观看吧| 黄色a级毛片大全视频| 日本91视频免费播放| 美女扒开内裤让男人捅视频| 1024香蕉在线观看| xxxhd国产人妻xxx| 视频区欧美日本亚洲| 午夜福利影视在线免费观看| 男女免费视频国产| 亚洲精品一卡2卡三卡4卡5卡 | 国产国语露脸激情在线看| 日本五十路高清| 女人精品久久久久毛片| 日本欧美视频一区| 亚洲精品国产色婷婷电影| 欧美日本中文国产一区发布|