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

    正壓防護(hù)服用供氣背包設(shè)計及其數(shù)值模擬與實(shí)驗(yàn)

    2025-07-18 00:00:00侯雨雷李希萌索槐李岳曾達(dá)幸
    中國機(jī)械工程 2025年6期

    中圖分類號:TP202 DOI:10.3969/j.issn.1004-132X.2025.06.023 開放科學(xué)(資源服務(wù))標(biāo)識碼(OSID):

    Design and Numerical Simulation and Experiments of Gas Supply Backpacks for Positive Pressure Protective Clothing

    HOU Yulei1 LI Ximeng1SUO Huai2LI YuelZENG Daxing' * 1.School of Mechanical Engineering,Yanshan University,Qinhuangdao,Hebei, 2.Qinhuangdao Gaha Medical Electronic Instrument Co.,Ltd.,Qinhuangdao,Hebei, 3.School of Mechanical Engineering,Dongguan University of Technology,Dongguan, Guangdong,

    Abstract:At present,the function of positive pressure protective clothing is relatively single,the high internal temperature and difficult pressure control resulted in low safety and comfort. Therefore, a supporting device integrating supply/exhaust, filtration and cooling functions was proposed. The data required to maintain a slightly positive pressure state within the protective clothing was tested,and the structure of the air-supply backpack was designed in modules. A heat transfer analysis was conducted on the semiconductor refrigeration module,and a simulation model was established to analyze the impacts of profile dimensions on refrigeration performance based on CFD (computational fluid dynamics). The parameters of temperature and pressure were stabilized by fuzzy PID control,the transfer function was solved,the controller was designed and verified by Simulink simulation. The FLUENT simulation model of micro positive pressure protective clothing was built to analyze the distribution law of flow field and temperature field in protective clothing and the influences of air supply temperature. The prototype of air supply backpacks was developed and the temperature and humidity of protective clothing were tested. The results show that the air supply backpacks may dehumidify to normal ambient humidity within 3O s and reduce 5°C compared with room temperature within 100s

    Key words: positive pressure protective clothing; gas supply backpack; semiconductor refrigera-tion;FLUENT simulation; temperature and humidity testing

    0 引言

    防護(hù)服是各種污染環(huán)境中進(jìn)行檢測、救援等任務(wù)的基本防護(hù)裝備,但目前傳染病防控及生化實(shí)驗(yàn)室人員所用防護(hù)裝備生理舒適性較差[1],正壓防護(hù)服若在保證用戶安全的同時能夠?yàn)槠涮峁崈艨諝?,并持續(xù)排出身體熱量和產(chǎn)生的濁氣,則可有效提高穿戴的舒適性[2]。防護(hù)服密封性強(qiáng)、透氣性差,尤其在天氣炎熱的夏季,易造成醫(yī)護(hù)人員中暑3,故保持防護(hù)服內(nèi)溫度和濕度值處于合適范圍內(nèi)十分必要。

    供風(fēng)裝置作為正壓防護(hù)服的核心組件[4],需要考慮供風(fēng)量、正壓值、噪聲大小、使用壽命以及異常報警等多方面。李國棟等[5設(shè)計了一款防止放射性污染的正壓防護(hù)服動力送風(fēng)單元,內(nèi)置高效過濾器及微型風(fēng)機(jī),設(shè)有流量、溫度反饋控制功能。目前,國內(nèi)外暫無正壓防護(hù)服統(tǒng)一標(biāo)準(zhǔn)6吳金輝等[通過氣密性壓力測試得出靜止時正壓防護(hù)服內(nèi)部正壓維持在 100~110Pa 之間,行走時維持在 70~200Pa 之間。劉倩靜8應(yīng)用人體工效學(xué)對正壓防護(hù)設(shè)備進(jìn)行綜合效能評估,結(jié)果表明,正壓防護(hù)設(shè)備進(jìn)風(fēng)量達(dá) 75~125L/min 時,可滿足負(fù)荷勞動下的生理需求。

    在保證防護(hù)安全的前提下,減小氣壓可增強(qiáng)穿戴的便利性,微正壓防護(hù)服內(nèi)氣壓微大于防護(hù)服外氣壓,屬于正壓防護(hù)服范疇,但當(dāng)前對此缺乏明確的界定。正壓防護(hù)服的氣壓控制[9要求較高,李鑫等[10提出一套微正壓防護(hù)服系統(tǒng),設(shè)定氣流量范圍為 60~130L/min ,可排除汗液,但未實(shí)現(xiàn)自動調(diào)壓。供氣型正壓防護(hù)服的研究還比較初步,其應(yīng)用仍需進(jìn)一步解決包括安全、舒適性等在內(nèi)的諸多問題[11]。

    相較于液體降溫及相變降溫服,半導(dǎo)體制冷[12]具有效果好、無噪聲及易攜帶等優(yōu)點(diǎn)。丁喜梅[13]設(shè)計了一款半導(dǎo)體制冷與水冷散熱結(jié)合的裝置,降溫效果好,但混合降溫服結(jié)構(gòu)較為復(fù)雜,操作不便。LOU等[14]利用半導(dǎo)體制冷片研發(fā)了一種可降溫除濕的穿戴式防護(hù)裝備,每小時可從防護(hù)服內(nèi)散發(fā)51.7W的熱量和 26.3g 的水分。SU等[15]開發(fā)了一種便攜式冷卻裝備,在 18~20 °C 空氣溫度下, 50L/min 的冷氣流量在頭部位置可營造一種舒適環(huán)境,但送風(fēng)難以流向腿部。

    半導(dǎo)體熱電制冷模塊的性能關(guān)乎降溫量和除濕率,通??捎脤?shí)驗(yàn)或FLUENT仿真的方法予以驗(yàn)證[16]。ELTAWEEL等[17]搭建實(shí)驗(yàn)平臺,通過改變半導(dǎo)體制冷片的輸人電壓及型材旋轉(zhuǎn)角度來研究熱電制冷模塊發(fā)揮出最佳性能時的條件參數(shù)。劉尚杰[18]建立了散熱型材模型并通過FLU-ENT軟件分別驗(yàn)證了型材長度、翅片間距及翅片厚度對除濕效果的影響。呂明杰[19]將熱電制冷模塊分為進(jìn)口域、出口域、內(nèi)部流體域及固體域,研究了導(dǎo)冷型材翅高、翅厚及翅間距對導(dǎo)熱性能的影響。當(dāng)前研究通常為單獨(dú)分析散熱端或?qū)Ю涠?,忽略了兩端之間的熱量傳遞及其相互影響,未將兩者結(jié)合一同分析,也未曾考慮型材底板厚度對導(dǎo)熱性能的影響。

    此外,通過FLUENT軟件還可分析防護(hù)服內(nèi)部空氣流動狀態(tài)及溫濕度分布。XU等[20]設(shè)計了一款新型通風(fēng)防護(hù)服,建立了簡化人體模型,并進(jìn)行了溫度分布仿真分析。SHEN等[21]研究了環(huán)境溫度及服裝材質(zhì)對傳熱的影響,得到了人體表面的氣體流動路徑和速度云圖。陳慧臻22采用空心圓柱模擬防護(hù)服軀干,并建立了曲面人體模型,分別研究了防護(hù)服厚度、流速等參數(shù)對防護(hù)服導(dǎo)熱性能的影響。供風(fēng)量大小、環(huán)境溫度等參數(shù)對正壓防護(hù)服的氣壓、溫度分布均有較重要的影響,但目前相關(guān)研究仍欠充分。

    綜上所述,當(dāng)前正壓防護(hù)服尚未有統(tǒng)一的微正壓范圍標(biāo)準(zhǔn),壓力控制還不夠精確,穿著后防護(hù)服內(nèi)部溫度較高,舒適感較低,為此,本文提出一種集供/排氣、多級過濾、制冷功能于一體的供氣背包,進(jìn)行相關(guān)理論計算、仿真分析及實(shí)驗(yàn)測試,以期提高正壓防護(hù)服穿戴舒適性。

    1供氣背包總體結(jié)構(gòu)設(shè)計

    1.1 正壓防護(hù)服微氣壓測試

    為測試防護(hù)服微正壓數(shù)值,制作供氣背包如圖1所示,以風(fēng)機(jī)為動力來源,氣體由過濾模塊進(jìn)入,通過供風(fēng)接頭送至防護(hù)服,防護(hù)服內(nèi)的濁氣由排風(fēng)接頭連接排至外界。

    圖1氣壓測試背包樣機(jī)Fig.1 Prototypeofairpressure testbackpack

    采用數(shù)字壓力表分別測試防護(hù)服達(dá)到微鼓狀態(tài)、蹲下和站立瞬間頭部和腹部的氣壓值。防護(hù)服尺碼越大,體積總量越大,其氣壓變化程度越大,蹲下及站立瞬間防護(hù)服內(nèi)氣體瞬時聚集或外排,氣壓迅速增大或減小。選擇160、170、180等多種型號防護(hù)服,分別記錄蹲下瞬時最大氣壓值和站立瞬時最小氣壓值。以170型號為例,記錄頭部、腹部不同狀態(tài)下氣壓測量值分別如表1、表2所示。

    同一試驗(yàn)條件下,腹部的氣壓值一般較頭部小,考慮拆裝方便,將氣壓傳感器放置于腹部,且壓差范圍設(shè)置為 5~20Pa ,則無論何種動作狀態(tài)下,頭部均能保持正壓狀態(tài),此狀態(tài)可阻止外界病毒侵入,符合安全要求。

    表1170型號防護(hù)服頭部不同狀態(tài)下的氣壓測量值 Tab.1Air pressure measurement value of 170 type protective clothing head under different conditionsPa
    表2170型號防護(hù)服腹部不同狀態(tài)下的氣壓測量值 Tab.2Air pressure measurement value of 170 type protective clothing under different abdominal conditions

    1.2 供氣背包結(jié)構(gòu)設(shè)計

    供氣背包總體結(jié)構(gòu)可分為過濾模塊、背包殼體、電源模塊、供氣及連接模塊、制冷模塊和控制模塊。從功能上分為普通供氣型和制冷型兩種,設(shè)計參數(shù)要求如表3所示。其中,出口降溫量指制冷型背包開啟制冷功能后,防護(hù)服進(jìn)風(fēng)口空氣的溫度相較室溫的降低量。

    表3供氣背包設(shè)計參數(shù)要求

    PaTab.3 Requirementsfordesignparametersofair supplybackpack

    過濾裝置決定著防護(hù)病毒等級、過濾效率和風(fēng)阻大小??紤]實(shí)現(xiàn) 99.99% 的過濾效率, 24h 的過濾壽命[23],濾芯材質(zhì)選用H14熔噴布,適量活性炭,濾芯進(jìn)風(fēng)口表面、活性炭進(jìn)/出風(fēng)口表面分別加一層F9等級的無紡布,過濾模塊結(jié)構(gòu)如圖2所示。

    過濾模塊設(shè)有壽命計時器,當(dāng)壽命不足時,系統(tǒng)會適時報警提醒。過濾外罩兩側(cè)與過濾上殼體之間留有間隙,即過濾裝置進(jìn)風(fēng)口,過濾外罩還可以阻止大面積灰塵、雜質(zhì)進(jìn)入過濾裝置,以延長過濾壽命。

    供氣連接模塊包括供氣背包進(jìn)出氣口、供風(fēng)燕尾槽 過濾殼體過濾外罩 出氣口底層無紡布活性炭 初層無紡布密封膠 隔板濾芯 濾芯支撐架 進(jìn)氣口卡緊結(jié)構(gòu) 壽命計時器(a)過濾模塊外殼體 (b)過濾模塊內(nèi)部布局管道和風(fēng)機(jī),風(fēng)機(jī)又分為供風(fēng)風(fēng)機(jī)、備用風(fēng)機(jī)、回 抽風(fēng)機(jī)和散熱風(fēng)機(jī)。供風(fēng)管道采用PE材質(zhì),供 氣背包與過濾裝置間的通氣接頭以雙O形密封 圈連接,O形圈選用硅膠材質(zhì),線徑為 1.9mm ,外 徑為 34mm ,壓縮率為 21.05% ;風(fēng)機(jī)按照各自功 能,考慮用電量、噪聲及供風(fēng)量等參數(shù)選取相應(yīng)型 號。通氣接頭結(jié)構(gòu)剖視圖見圖3。

    圖2過濾模塊結(jié)構(gòu)示意圖Fig.2 Schematic diagram of filtermodule structure圖3通氣接頭結(jié)構(gòu)剖視圖Fig.3 Sectionviewofventilationjointstructure

    電源模塊包括鋰電池及其外殼、電池倉,其中采用六、九節(jié)鋰電池分別作為普通供氣型和制冷型背包的供電來源,計算公式如下:

    式中: W 為用電功率; Ub 為用電電壓; Ib 為用電電流; E 為電池能量; V 為電池額定電壓; Ac 為電池容量; tc 為電池可用電時間。

    根據(jù)供電時間要求,經(jīng)計算,六、九節(jié)鋰電池容量分別為 Ac?8.125A?h 和 Ac?15A?h 結(jié)合電池規(guī)格,分別選取 10A?h 和 15A?h 容量的電池組。

    背帶用以承載背包,遵循人體工程學(xué)設(shè)計,左側(cè)肩帶夾層內(nèi)留直徑 4mm 線束通道,用來通信走線。背包殼體承載著所有零部件,需有一定的硬度,且希望密度小,能夠抗紫外線,故選用ASA材質(zhì)。殼體模塊如圖4所示。

    制冷型背包采用二級制冷模塊裝置,每一級制冷模塊均設(shè)有兩個12704半導(dǎo)體制冷片,兩級制冷模塊插接而成,用戶可根據(jù)實(shí)際需求選擇一級制冷還是二級制冷,一級制冷模塊模型見圖5。

    普通供氣背包可實(shí)現(xiàn)防護(hù)服內(nèi)空氣過濾、供風(fēng)和除濕;制冷型背包除前述功能外,還可實(shí)現(xiàn)降溫功能,整體結(jié)構(gòu)如圖6所示。

    圖4供氣背包殼體模塊
    Fig.4Gas supply backpack housing module
    Fig.5Schematicdiagramofrefrigeration圖6供氣背包整體示意圖Fig.6 Overall schematic diagram of gas supply backpack

    2 制冷模塊數(shù)值模擬分析

    2.1 制冷模塊參數(shù)設(shè)置

    制冷模塊由半導(dǎo)體制冷片、散熱型材、導(dǎo)冷型材、風(fēng)機(jī)、通道及保溫棉組成。進(jìn)行傳熱分析如下:假設(shè)熱電制冷系統(tǒng)傳熱只在垂直方向發(fā)生,外加電場在熱電偶上均勻一致。以一對熱電偶為例,熱電偶對長為 l ,面積為 A ,通過的電流為 I ,電阻為 R ,塞貝克系數(shù)為 s ,導(dǎo)熱參數(shù) λ 為常數(shù),熱電偶內(nèi)部溫度為 Td ,熱電制冷系統(tǒng)冷端溫度為Tc 、熱端溫度為 Th ,在第三類邊界條件 x=0 情況下,流體溫度為 Tfl ,對流傳熱系數(shù)為 α1:x=l 條件下流體溫度為 Tf2 ,對流傳熱系數(shù)為 α2°qV 為均勻內(nèi)熱源。熱電制冷系統(tǒng)工作穩(wěn)定時,熱傳遞方程[24]為

    d2Td/dx2+qV/λ=0

    第三類邊界條件:

    -λAdTd/dx=sITh1A(Th-Tfl)x=0

    -λAdTd/dx=sITc2A(Tf2-Tc)x=l

    綜合式 (4)~ 式(6),得到熱電偶內(nèi)部的溫度分布為

    0?x?l

    制冷量 Qg

    熱端散熱量為

    Qh=sITh+λA(dTd/dx)x=01A(Th-Tf1

    系統(tǒng)輸入功率為

    P=Qh-Qg=I2R+sI(Th-Tc

    系統(tǒng)輸入電壓為

    U=P/I=IR+s(Th-Tc

    系統(tǒng)制冷系數(shù)為

    制冷型供氣背包采用TEC1-12704型半導(dǎo)體制冷片制冷,參數(shù) l,A,s,λ 均為常數(shù),提高熱電系統(tǒng)的制冷性能,需從散熱性能、導(dǎo)冷方式和制冷模塊結(jié)構(gòu)來考慮。由式(10)可知,散熱量 Qh 與制冷量 Qg 成正相關(guān),若要提高 Qg ,需要良好的散熱結(jié)構(gòu)。由式(8)可知, Qg 還與對流傳熱系數(shù)、冷熱端溫差大小有關(guān),因此應(yīng)增加散熱量和導(dǎo)冷量以縮小半導(dǎo)體制冷片冷熱端的溫差,減少其間的熱量傳遞。

    利用FLUENT軟件分別對制冷模塊的散熱端和導(dǎo)冷端進(jìn)行仿真,以分析不同結(jié)構(gòu)參數(shù)和環(huán)境變量下的制冷性能。建立仿真模型并劃分流體域結(jié)果如圖7所示。

    圖5制冷模塊結(jié)構(gòu)示意圖圖7 制冷模塊流體域劃分Fig.7 Fluiddomaindivisionofrefrigerationmodule

    制冷模塊內(nèi)嵌在封閉外殼內(nèi),不與外界環(huán)境發(fā)生熱交換,是一個涉及固體域和流體域?qū)α?、?dǎo)熱互相耦合的過程,隨著半導(dǎo)體制冷片、熱端、冷端和流體域之間的相互影響,系統(tǒng)最終會達(dá)到穩(wěn)定狀態(tài)。

    無論是散熱端還是導(dǎo)冷端,空氣流動均由風(fēng)機(jī)提供動力,可視為強(qiáng)制對流換熱,由于需要研究散熱端和導(dǎo)冷端溫度隨時間的變化趨勢,因此進(jìn)行三維瞬態(tài)計算,其控制方程如下:

    式中: u,v,w 分別為 x,y,z 方向的分速度; τ 為半導(dǎo)體制冷片工作時間; ρ 為空氣密度; p 為大氣壓; cp 為空氣質(zhì)量熱容; T 為型材表面溫度: 為空氣熱導(dǎo)率; μ 為空氣黏度。

    使用雷諾數(shù)判斷空氣流動狀態(tài):

    M=3600FUv

    式中: Uv 為空氣流動速度; L 為特征長度; u1 為流體運(yùn)動黏度; M 為風(fēng)量, m3/h;F 為管道通風(fēng)面積, m2 。

    所選兩款風(fēng)機(jī)風(fēng)量分別為 7.5m3/h 和9.6m3/h 。以室溫 26°C 下的空氣運(yùn)動黏度 u1 為例,求得導(dǎo)冷通道的空氣流動速度為 2.95m/s ,雷諾數(shù)為5898,散熱通道的空氣流動速度為3.78m/s ,雷諾數(shù)為7550,雷諾數(shù)均大于2000,故兩個流體域均視為湍流[25]。

    采用標(biāo)準(zhǔn) k-ε 模型進(jìn)行求解,兩個流體域的入口邊界設(shè)置為速度入口,出口邊界設(shè)置為壓力出口,溫度為常數(shù)。流體域與固體交界處,速度無突變,溫度一致。

    FLUENT仿真計算中,為消除網(wǎng)格數(shù)量對計算結(jié)果的影響,有必要進(jìn)行網(wǎng)格的獨(dú)立性驗(yàn)證。選取導(dǎo)冷通道流體域進(jìn)出口熱交換效率為評估指標(biāo),計算結(jié)果誤差大小,并逐漸增加網(wǎng)格數(shù)量。當(dāng)兩個網(wǎng)格下的熱交換效率計算結(jié)果相對誤差足夠小時,認(rèn)為可以驗(yàn)證網(wǎng)格獨(dú)立性。熱交換效率計算公式如下:

    式中: Tin?Tout 分別為冷端進(jìn)出氣口溫度。

    相對誤差計算公式如下:

    式中: N 為網(wǎng)格劃分次數(shù); φN 為第 N 次劃分網(wǎng)格時的誤差; ηN,ηN-1 分別為第 N,N-1 次劃分網(wǎng)格的熱交換效率。

    由表4可知,隨著網(wǎng)格數(shù)量的增加,誤差逐漸減小,且第四、第五次的誤差接近,故綜合考慮仿真精度和計算時間成本,設(shè)置該模型的網(wǎng)格數(shù)量為238萬左右。

    表4流體域進(jìn)出口熱交換效率隨網(wǎng)格數(shù)的變化

    Tab.4Variation of heat exchange efficiency between inletandoutletof fluiddomainwithgrid number

    2.2 制冷模塊仿真分析

    分析散熱型材底板厚度、型材長度、齒片齒高以及導(dǎo)冷型材底板厚對制冷性能的影響,計算模型如圖8所示。其中, d. ,為型材底板厚度, dw 為型材齒寬,h為齒片齒高。

    圖8 計算模型圖Fig.8 Calculation model

    為測試不同底板厚度型材的散熱效果,對模型進(jìn)行瞬態(tài)仿真,分析不同底板厚度下散熱型材吸熱速率的變化規(guī)律。底板厚度分別為 3.0mm 、4.5mm 及 6.0mm 時底板內(nèi)側(cè)表面平均溫度隨時間的變化如圖9所示??梢?3mm 底板厚度型材表面溫度最高、吸熱速率最快,其次是 4.5mm ,較快的吸熱速率有助于型材快速散掉半導(dǎo)體制冷片熱端的熱量,進(jìn)而提高制冷能力,故選取 3mm 底板厚度散熱型材是較合適的。

    為在有限的空間內(nèi)實(shí)現(xiàn)最佳的散熱性能,進(jìn)一步探究型材尺寸對熱量分布的影響。選取四組不同長 × 寬 x 高的散熱型材對模型進(jìn)行瞬態(tài)仿真,分析型材長度、高度對散熱性能的影響。不同尺寸散熱型材散熱通道出風(fēng)口平均溫度隨時間變化如圖10所示??梢姽╋L(fēng)開始 30s 內(nèi), 120mm× 50mm×28mm 型材溫度上升最快,其次是160mm×50mm×28mm,200mm×50mm×28 mm 和 160mm×50mm×50mm 。由三種不同長度型材的對比結(jié)果可知,型材長度越小吸熱越快,三種型材表面最終穩(wěn)定溫度相近,故宜采用質(zhì)量較小的 120mm×50mm×28mm 型材;由齒高 28mm 和齒高 50mm 兩種型材的對比結(jié)果可知,型材高度越小吸熱越快,最終穩(wěn)定溫度也較高,有利于熱量快速散出。

    圖9 熱端不同底板厚度下內(nèi)表面溫度隨時間變化Fig.9 Variation diagramof innersurface temperaturewith time at different base plate thickness at hot end

    壓降用來表征型材的風(fēng)阻大小,風(fēng)阻越大,空氣流通受阻,不利于空氣帶走散熱型材導(dǎo)出的熱量,故壓降越小越好。四組型材通道內(nèi)的壓降 ?in 變化如圖11所示。結(jié)果可見, 200mm×50mm× 28mm 型材通道內(nèi)阻力最大,其次是 160mm× 50mm×50mm,160mm×50mm×28mm,120 mm×50mm×28mm 。分析可得,型材越長、高度越高,管道內(nèi)流動空氣的阻力越大,因此當(dāng)其他條件一致時,為保證通暢的空氣流動,宜選擇高度較小、長度較小的散熱型材。

    隨時間變化情況如圖12所示。結(jié)果可見,1.0mm 底板厚度型材表面溫度最低、導(dǎo)冷速率最快,其次是 2.5mm,5.0mm ,更快的導(dǎo)冷速率有助于快速傳遞半導(dǎo)體制冷片冷端的制冷量,進(jìn)而提高制冷能力,但底板過薄不易加工,加工成本高,綜合考慮選取 2.5mm 底板厚散熱型材是較合適的。

    圖10 不同型材散熱通道出風(fēng)口平均溫度隨時間變化Fig.10 The average temperature of the outlet of differentprofiles cooling channels changes with time圖11 不同尺寸型材通道內(nèi)的壓降曲線Fig.11 Pressuredropcurvesindifferentsizeprofilechannels圖12冷端不同底板厚內(nèi)表面溫度隨時間變化圖Fig.12Variation diagram of inner surface temperatureofcold endwithdifferent baseplate thicknesswith time

    3 模糊PID控制開發(fā)

    3.1 控制模型建立

    為保障穿戴的安全性和舒適性,防護(hù)服內(nèi)需維持穩(wěn)定、適宜的氣壓和溫度,本節(jié)根據(jù)防護(hù)服內(nèi)溫度、氣壓兩個參數(shù)的動態(tài)特性建立數(shù)學(xué)模型。

    防護(hù)服內(nèi)部為相對密閉空間,主要考慮由供風(fēng)風(fēng)機(jī)、回抽風(fēng)機(jī)進(jìn)行的空氣循環(huán)和熱交換過程。以降溫為例,防護(hù)服內(nèi)熱量與供氣背包制冷模塊提供的制冷量在時間序列上不斷變化,其表達(dá)式為

    式中: C 為防護(hù)服的容量系數(shù); θa 為防護(hù)服內(nèi)溫度; θb 為室溫; Ff 為防護(hù)服表面積; K 為防護(hù)服的熱導(dǎo)率。

    當(dāng)防護(hù)服內(nèi)溫度達(dá)到穩(wěn)定時,溫度變化率

    0,制冷模塊初始制冷量 Qg0

    Qg0=FfK(θa0b0

    式中: θa0、θb0 分別為防護(hù)服內(nèi)和室內(nèi)的初始溫度。

    當(dāng)防護(hù)服內(nèi)溫度處于動態(tài)時,則有

    其中, Δθa,ΔQg,Δθb 為各環(huán)節(jié)增量。

    將式(24)代入式(22),則有

    Qg0+ΔQg+FfK(θb0+Δθb

    將式(23)代人式(25),有

    由于室溫不會突變, Δθb=0 ,進(jìn)而有

    進(jìn)行拉氏變化,得初始傳遞函數(shù)為

    式中: k1 為靜態(tài)放大增益; T1 為時間常數(shù)。

    防護(hù)服內(nèi)的溫度變化存在滯后性,因此增加滯后環(huán)節(jié),得到溫度控制傳遞函數(shù)如下:

    式中: τ1 為滯后時間參數(shù)。

    同理,根據(jù)質(zhì)量守恒定律,防護(hù)服內(nèi)氣壓平衡方程為

    式中: Vd 為防護(hù)服內(nèi)體積; 為增壓量; ?s 為防護(hù)服內(nèi)氣壓; ΔG 為送風(fēng)量與抽風(fēng)量之差,即凈供風(fēng)量。

    同理,經(jīng)拉氏變化并增加滯后環(huán)節(jié),得到防護(hù)服氣壓控制傳遞函數(shù)如下:

    式中: Ψ:k2 為靜態(tài)放大增益; T2 為時間常數(shù); τ2 為滯后時間參數(shù)。

    3.2 傳遞函數(shù)求解

    式(30)式(32)被視為防護(hù)服內(nèi)溫度、氣壓模型的一階慣性環(huán)節(jié),未知參數(shù)通過階躍響應(yīng)法得出。以模擬蹲起時的增壓過程為例,供風(fēng)風(fēng)機(jī)和回抽風(fēng)機(jī)的風(fēng)量直接影響防護(hù)服內(nèi)氣壓的變化,實(shí)驗(yàn)測得防護(hù)服內(nèi)氣壓隨時間變化見圖13。由圖13可見, t=90 s時,防護(hù)服內(nèi)氣壓基本穩(wěn)定,供風(fēng)風(fēng)機(jī)與回抽風(fēng)機(jī)供氣量達(dá)到平衡,有 。系統(tǒng)給定的階躍電壓 u0 為 12V ,靜態(tài)放大增益 k2

    圖13 防護(hù)服內(nèi)氣壓變化曲線

    將階躍響應(yīng) y(t) 轉(zhuǎn)換為[0,1]內(nèi)量綱一形式,得

    在階躍響應(yīng)作用下, y*(t) 化簡為

    將 y?(t) 兩個數(shù)據(jù)組(30,0.41),(40,0.92)代人式(34),求得 T2=5,τ2=33 ,再代人式(32)得

    同理,對防護(hù)服內(nèi)溫度進(jìn)行采集,室溫 35°C 下采用四個制冷片以最大制冷功率進(jìn)行降溫,利用MATLAB軟件繪制溫度變化曲線,如圖14所示。

    Fig.13Air pressure change curve in protective clothing圖14 防護(hù)服內(nèi)溫度變化曲線Fig.14Temperature change curve in protective clothing

    由圖14可見,在 250s 左右防護(hù)服溫度基本達(dá)到平衡,其中 y(0)=35,y(∞)=22.6 ,求得

    將兩個數(shù)據(jù)組(90,0.52),(180,0.95)代入y?(t) 曲線中,求得 T1=40,τ1=61 ,進(jìn)而可得溫度控制傳遞函數(shù)為

    3.3 模糊PID控制器設(shè)計

    防護(hù)服內(nèi)溫度、氣壓參數(shù)存在時滯性和非線性,考慮傳統(tǒng)PID難以保證精度,故將模糊控制引入,其原理如圖15所示。圖中, Ii??I2 代表防護(hù)服內(nèi)溫度、氣壓目標(biāo)值, Y1,Y2 代表溫度、氣壓實(shí)際值, e1e2 代表實(shí)際值與目標(biāo)值的偏差, ec1?ec2 代表偏差變化率, O1…O2 代表模糊器輸出值, Δkp ?Δki?Δkd 為輸出參數(shù)增量。

    圖15溫度和氣壓模糊PID控制原理圖Fig.15 Schematic diagram of temperature andpressure fuzzy PID control

    此系統(tǒng)是兩參數(shù)二輸入三輸出結(jié)構(gòu),可隨工況實(shí)時調(diào)整,在降溫過程中,模糊PID控制輸出0~12V 電壓,以改變制冷片輸人功率;調(diào)壓過程中模糊PID控制輸出 0/12V 電壓,以改變回抽風(fēng)機(jī)PWM占空比,傳感器實(shí)時檢測,直至達(dá)到穩(wěn)定狀態(tài)。

    以溫度控制為例設(shè)計模糊控制器。首先根據(jù)測量數(shù)據(jù)取值范圍確定模糊論域,溫度誤差基本論域?yàn)閇-2.7,2.7],模糊論域?yàn)閇-2,2],量化因子為0.7,偏差率基本論域?yàn)椋邸?.4,1.4],模糊論域?yàn)閇—2,2],量化因子為1.4。比例增量模糊論域?yàn)閇一3,3],量化因子為0.04;積分增量模糊論域?yàn)閇一3,3],量化因子為0.000004;微分增量模糊論域?yàn)閇一3,3],量化因子為8。防護(hù)服內(nèi)溫度不會突變,輸入、輸出隸屬函數(shù)取三角形,模糊子集取{NB,NM,NS,ZO,PS,PM,PB},利用Mamdani模糊推理法建立模糊規(guī)則。通過重心法進(jìn)行模糊數(shù)據(jù)清晰化(即解模糊)對PID三個輸出參數(shù)進(jìn)行修正,從而獲得更加穩(wěn)定的性能。

    設(shè)置防護(hù)服初始溫度為 ,目標(biāo)值為22°C ,初始?xì)鈮褐禐?-15Pa ,目標(biāo)值為 20Pa ,利用Simulink創(chuàng)建傳統(tǒng)PID與模糊PID溫度控制仿真模型,如圖16所示,相應(yīng)仿真曲線見圖17。

    同理,氣壓控制仿真模型如圖18所示,仿真曲線見圖19。

    圖17溫度控制仿真曲線
    Fig.17Temperature control simulation curve
    圖16傳統(tǒng)PID與模糊PID溫度控制仿真模型 Fig.16 Simulation model of traditional PID and fuzzy PID temperature control圖18傳統(tǒng)PID與模糊PID氣壓控制仿真模型 Fig.18 Simulation model of traditional PID and fuzzy PID air pressure control圖19氣壓控制仿真曲線Fig.19 Pressure control simulation curve

    可見,模糊PID控制較傳統(tǒng)PID更快地實(shí)現(xiàn)了防護(hù)服內(nèi)的加壓和降溫目標(biāo),精確性也更高,防護(hù)服內(nèi)溫度、氣壓可快速準(zhǔn)確地達(dá)到設(shè)定值,有助于提高防護(hù)可靠性及穿戴者的舒適性。

    4正壓防護(hù)服仿真與實(shí)驗(yàn)

    4.1 正壓防護(hù)服流體仿真

    利用FLUENT軟件進(jìn)行正壓防護(hù)服仿真,需先搭建三維人體和防護(hù)服模型,并將兩者一體化處理。取人體與防護(hù)服間距為 5mm 來模擬微正壓狀態(tài),進(jìn)出氣口以圓柱孔替代,依次進(jìn)行抽取流體域、網(wǎng)格劃分、網(wǎng)格無關(guān)性驗(yàn)證和求解計算等,最終確定網(wǎng)格數(shù)量為 80.5×104 個,計算求解模型如圖20所示。

    參照實(shí)驗(yàn)條件,仿真參數(shù)取室溫 20.5°C ,穿著防護(hù)服時人體舒適的平均表面溫度取 33°C ,防護(hù)服供風(fēng)口設(shè)置為速度人口,大小為 1.5m/s ,排風(fēng)口設(shè)置為速度出口,大小為 -1.5m/s ,以殘差大小為 10-9 作為收斂條件,系統(tǒng)穩(wěn)定后,防護(hù)服內(nèi)部溫度流場分布如圖21所示。可見,靠近進(jìn)出氣口的背側(cè)溫度更低,這是因?yàn)楸硞?cè)空氣流動速度快,熱量可以及時排出,而手臂處空氣不易排出,故溫度相對較高。

    圖20 FLUENT仿真及求解模型
    Fig.20 FLUENTsimulation and solutionmodel

    在室溫 20.5°C 情況下,分別取供風(fēng)溫度 θ 為17.5°C?14.5°C?11.5°C 和 8.5°C ,其他條件不變,分析供風(fēng)溫度對防護(hù)服熱舒適性的影響,得到流體域背側(cè)溫度分布如圖22所示。可見,不同供風(fēng)溫度下,均出現(xiàn)了不同程度的熱量分布不均的情況,溫度越高,不均勻情況越明顯。以脖頸下方為測溫點(diǎn),對應(yīng)不同供風(fēng)溫度下的測量值依次為19.5°C?18.2°C?16.6°C 和 14.7°C 。相較于普通供風(fēng),制冷供風(fēng)起到了較為明顯的降溫作用。以低于室溫 12°C 的空氣供氣,測溫點(diǎn)較室溫降低了5.8°C ,大部分位置仍小于 3°C ,故要取得明顯的降溫效果,制冷型供氣背包應(yīng)具備 12°C 以上的制冷能力。

    圖21 FLUENT防護(hù)服內(nèi)流場及溫度場分布 Fig.21 Distribution of flow field and temperature field in FLUENT protective clothing圖22不同供風(fēng)溫度下背側(cè)溫度分布云圖
    Fig.22Nephogram of temperaturedistribution at the backunder different air supply temperatures

    4.2 正壓防護(hù)服溫度濕度測試

    防護(hù)服內(nèi)的溫度濕度值是決定用戶熱舒適性的關(guān)鍵因素,為此,搭建正壓防護(hù)服供氣套裝,如圖23所示。其中,普通供氣背包總體尺寸為340mm×330mm×120mm ,質(zhì)量為 2.4kg ;制冷型供氣背包總體尺寸為 340mm×330mm×205 mm ,質(zhì)量為 4.1kg 。

    圖23防護(hù)服供氣套裝測試工裝

    將溫濕度傳感器放置于脖頸下方,在溫度20.5°C 、濕度 31.7% 的環(huán)境下,分別測試無供風(fēng)、普通供風(fēng)及制冷供風(fēng)三種情況下防護(hù)服內(nèi)的溫度濕度變化。如圖24所示,無供風(fēng)時,防護(hù)服內(nèi)溫度濕度迅速提高,180s時濕度 RH 已達(dá)到90.5% ,溫度 30.5°C ,出現(xiàn)嚴(yán)重憋悶感;開啟普通供風(fēng)后,防護(hù)服內(nèi)溫濕度值先升高后降低直至平緩,升高是因?yàn)槿梭w呼吸所致,降低是因?yàn)楣╋L(fēng)加快了空氣循環(huán);普通供風(fēng)較于無供風(fēng)有較為明顯的降溫除濕功能,但增大風(fēng)擋對降溫除濕效果的提升有限;普通供風(fēng)情況下最終防護(hù)服內(nèi)濕度與外界接近,溫度高于外界溫度。

    Fig.23Protective clothingair supply suit test tooling圖24普通供風(fēng)時防護(hù)服內(nèi)的溫度濕度變化 Fig.24Changes of temperature and humidity in protectiveclothingduringnormal airsupply

    同樣在室溫 20.5°C 情況下,分別取供風(fēng)溫度為 17.5°C?14.5°C?11.5°C 和 8.5°C 進(jìn)行制冷供風(fēng)實(shí)驗(yàn),防護(hù)服內(nèi)溫濕度變化如圖25所示。由圖25可知,開啟制冷供風(fēng)后,防護(hù)服內(nèi)溫度濕度值變化趨勢仍為先升高后降低直至平緩,相較于普通供風(fēng),除濕降溫幅度更大,隨著供風(fēng)溫度的降低,濕度值也在逐漸下降,但效果并非特別顯著;防護(hù)服內(nèi)溫度隨著供風(fēng)溫度降低而減小,當(dāng)供風(fēng)溫度為 8.5°C 時,防護(hù)服溫度降至 15.2°C ,較室溫降低了 5.3°C ,較普通供風(fēng)降低了 8.6°C 。

    圖25制冷供風(fēng)時防護(hù)服內(nèi)的溫度濕度變化 Fig.25Temperature and humidity changes in protective clothing during cooling and air supply

    綜合可知,濕度方面,普通供氣背包通風(fēng)30s左右,防護(hù)服內(nèi)的濕度值已降低至與環(huán)境基本一致并逐漸達(dá)到穩(wěn)定,該款即可滿足除濕的功能需求;溫度方面,實(shí)驗(yàn)與仿真數(shù)據(jù)規(guī)律一致,普通供風(fēng)降溫效果有限,而制冷供風(fēng)效果較為顯著,制冷模塊所產(chǎn)生的熱量在通風(fēng) 100s 時可實(shí)現(xiàn)防護(hù)服內(nèi)溫度較環(huán)境溫度降低 5°C 左右。

    5結(jié)論

    1)完成了多個尺碼防護(hù)服氣壓測試,確定腹部測量位置較佳的微正壓壓差范圍為 5~20Pa ·設(shè)計了供氣背包各個模塊結(jié)構(gòu),對制冷模塊傳熱分析發(fā)現(xiàn),增強(qiáng)制冷模塊性能可從減少散熱和導(dǎo)冷兩端溫差入手,即加快熱端散熱和冷端導(dǎo)冷速度;仿真結(jié)果表明,散熱型材高度對散熱效果影響較明顯,型材底板厚度對導(dǎo)熱性能影響較小。

    2)提出基于溫度、氣壓雙參數(shù)的模糊PID自動控制策略,設(shè)計控制器并進(jìn)行Simulink仿真,結(jié)果表明,相較于傳統(tǒng)PID,采用本文控制策略可提前10s左右達(dá)到氣壓目標(biāo)設(shè)定值,溫控準(zhǔn)確度也更高,有助于保障正壓防護(hù)的可靠性及穿戴者的舒適性。

    3)基于FLUENT軟件完成防護(hù)服內(nèi)流場和溫度分布分析,試制了供氣背包樣機(jī),開展了供氣量和溫濕度測試。仿真與實(shí)驗(yàn)結(jié)果表明,所設(shè)計供氣背包實(shí)現(xiàn)了預(yù)期的降溫、除濕功能,能夠較好地提高穿戴者的舒適性。

    參考文獻(xiàn):

    [1] 周珣,吳洪兵,鄔佩瑤.疫情背景下供氣型正壓醫(yī)用防護(hù)服的研究與設(shè)計[J].黑龍江科學(xué),2022,8(13):22-24.ZHOU Xun,WU Hongbing,WU Peiyao.Researchand Design of Gas Supply Positive Pressure MedicalProtective Clothing under the Background of theEpidemic[J]. Heilongjiang Science,2022,8(13):22-24.

    [2] 劉曉敏,楊棟,陳亮,等.融合RE與BID的產(chǎn)品創(chuàng)新方法研究及應(yīng)用[J].機(jī)械工程學(xué)報,2022,58(23):227-239.LIUXiaomin,YANG Dong,CHEN Liang,etal.Research and Application of Product InnovationMethods Integrating RE and BID[J]. Journal of Me-chanical Engineering,2022,58(23):227-239.

    [3] YIW,CHAN A P.Optimizing Work-rest Schedulefor Construction Rebar Workers in Hot and HumidEnvironment[J]. Building and Environment,2013(6):104-113.

    [4] HAOL,QIJ,WANG R,etal.Protective Per-formance Testand Safety Risk Evaluation ofaPow-ered Air-purifying Suit[J].Biosafety and Health,2019,1(2):91-97.

    [5] 李國棟,陳法國,韓毅,等.放射性污染通風(fēng)防護(hù)服動力送風(fēng)單元設(shè)計[J].核電子學(xué)與探測技術(shù),2020,40(4):633-637.LI Guodong,CHEN Faguo,HAN Yi,et al. Designof Power Supply Air Unit for Radioactive PollutionVentilation and Protection Clothing[J]. NuclearElectronics and Detection Technology, 2020, 40(4) :633-637.

    [6]羅慶祥,吳碧君.醫(yī)用正壓防護(hù)服結(jié)構(gòu)設(shè)計及主要性能指標(biāo)淺析[J].中國醫(yī)療器械信息,2023,29(11):154-156.LUO Qingxiang,WU Bijun. Analysis of StructuralDesign and Main Performance Indicators of MedicalPositive Pressure Protective Clothing [J] ChineseMedical Device Information,2023,29(11):154-156.

    [7]吳金輝,衣穎,楊荊泉,等.全身正壓生物防護(hù)服防護(hù)性能研究[J].中國安全科學(xué)學(xué)報,2011,21(9):125-130.WU Jinhui,YI Ying,YANG Jingquan,et al. Re-search on the Protective Performance of Full BodyPositive Pressure Biological Protective Clothing[J]Chinese Journal of Safety Sciences,2011,21(9):125-130.

    [8]劉倩靜.正壓生物防護(hù)頭罩綜合效能評估研究[D].北京:軍事科學(xué)院,2022:22-30.LIU Qianjing. Research on the Comprehensive Efi-ciency Evaluation of Positive Pressure BiologicalProtective Headgear[D]. Beijing:Academy of Mili-tary Sciences,2022:22-30.

    [9]楊望明,楊鵬明,楊明亮.一種適用于正壓防護(hù)服的控制系統(tǒng):CN112631181A[P].2021-04-09.YANG Wangming, YANG Pengming, YANGMingliang. AControl System Suitable for PositivePressure Protective Clothing: CN112631181A[P].2021-04-09.

    [10]李鑫,李碩,顧陽光.采用兩級過濾的微正壓防護(hù)服系統(tǒng):CN113413687A[P].2021-09-21.LI Xin,LI Shuo,GU Yangguang. A Micro Posi-tive Pressure Protective Suit System Using Two-stage Filtration:CN113413687A[P]. 2021-09-21.

    [11]韓玲.改善醫(yī)用一次性防護(hù)服熱濕舒適性的研究進(jìn)展[J].棉紡織技術(shù),2020,48(4):75-78.HAN Ling. Research Progress on Improving theThermal and Humidity Comfort of DisposableMedical Protective Clothing[J]. Cotton TextileTechnology,2020,48 (4):75-78.

    [12]蘇宇,何寧,李亮,等.新型切削用冷風(fēng)發(fā)生裝置的研制及性能測試[J].中國機(jī)械工程,2007,18(10):1208-1211.SU Yu,HE Ning,LI Liang,et al. Developmentand Performance Testing of a New Type of ColdAir Generation Device for Cutting[J]. China Mechanical Engineering,2007,18(10):1208-1211.

    [13]丁喜梅.新型半導(dǎo)體降溫防護(hù)服研究與設(shè)計[D].西安:西安科技大學(xué),2018:35-60.DING Ximei. Research and Design of New Semi-conductor Cooling Protective Clothing[D]. Xian:Xian University of Science and Technology,2018:35-60.

    [14]LOU L, ZHOU Y Y,YAN Y S,et al. Wearable-cooling and Dehumidifying System for PersonalProtective Equipment[J]. Energy Buildings,2022(5):276-280.

    [15]SU X,TIAN SC,LI H,et al. Thermal and Humid Environment Improvement of the ProtectiveClothing for Medical Use with a Portable CoolingDevice:Analysis of Air Supply Parameters[J]. En-ergy Buildings,2021,24:1109-1112.

    [16]程顯耀.半導(dǎo)體制冷熱端散熱器傳熱特性研究[D].濟(jì)南:山東大學(xué),2016:50-56.CHENG Xianyao. Research on Heat TransferCharacteristics of Semiconductor Refrigeration HotEnd Heat Sink[D]. Jinan: Shandong University,2016:50-56.

    [17]ELTAWEEL M,HEGGY A H,YASEEN Z M,et al. Application of the ANOVA Method in theOptimization of a Thermoelectric Cooler-based De-humidification System[J]. Energy Reports,2022,81:1050-1054.

    [18]劉尚杰.基于熱電制冷的溫濕度調(diào)節(jié)機(jī)整體設(shè)計與關(guān)鍵技術(shù)研究[D].鄭州:鄭州大學(xué),2021:55-60.LIU Shangjie. The Overall Design and Key Tech-nology Research of Temperature and HumidityRegulation Machine Based on Thermoelectric Re-frigeration[D]. Zhengzhou: Zhengzhou University,2021:55-60.

    [19]呂明杰.半導(dǎo)體制冷耦合蒸氣壓縮制冷試驗(yàn)與數(shù)值模擬研究[D].長春:吉林大學(xué),2023:19-30.LYU Mingjie. Experimental and Numerical Simu-lation Research on Semiconductor RefrigerationCoupled with Vapor Compression Refrigeration[D]. Changchun:Jilin University,2023:19-30.

    [20]XU X,ZHANG L,DEYU M,et al. Research onthe Novel Medical Protective Clothing for COVID-19[J].Heliyon,2023,9(2):1337-1344.

    [21]SHEN H,ANYY,ZHANG HM,et al.3D Nu-merical Investigation of the Heat and Flow Trans-fer through Cold Protective Clothing Based on CFD[J]. International Journal of Heat and Mass Transfer,2021,175:188-197.

    [22]陳慧臻.基于CFD的人體-防化服-環(huán)境系統(tǒng)傳熱研究[D].蘇州:蘇州大學(xué),2022:20-60.

    日韩人妻高清精品专区| 免费看美女性在线毛片视频| 我要看日韩黄色一级片| 成人av在线播放网站| 久久精品国产99精品国产亚洲性色| 青春草视频在线免费观看| 岛国毛片在线播放| 国产白丝娇喘喷水9色精品| 老司机福利观看| 免费观看人在逋| 卡戴珊不雅视频在线播放| 一本久久中文字幕| 国产麻豆成人av免费视频| 99热全是精品| 99热这里只有精品一区| 久久久精品大字幕| 日日摸夜夜添夜夜添av毛片| 国产毛片a区久久久久| 成人特级黄色片久久久久久久| 国产精品综合久久久久久久免费| 国产精品无大码| 免费电影在线观看免费观看| 国产精品国产三级国产av玫瑰| 国产精品爽爽va在线观看网站| 久久6这里有精品| 久久九九热精品免费| 久久精品夜色国产| 欧美成人精品欧美一级黄| 高清在线视频一区二区三区 | 麻豆国产av国片精品| 中文字幕精品亚洲无线码一区| 晚上一个人看的免费电影| 国产综合懂色| 午夜视频国产福利| 岛国毛片在线播放| 两个人的视频大全免费| 村上凉子中文字幕在线| 久久99热这里只有精品18| 日本五十路高清| 在线天堂最新版资源| 搞女人的毛片| 黄色欧美视频在线观看| 国产91av在线免费观看| 久久午夜亚洲精品久久| 美女xxoo啪啪120秒动态图| 国产高清激情床上av| 麻豆久久精品国产亚洲av| 国产午夜精品论理片| 成人亚洲精品av一区二区| 99在线视频只有这里精品首页| 日韩成人av中文字幕在线观看| av免费在线看不卡| 人人妻人人澡欧美一区二区| 夫妻性生交免费视频一级片| 91久久精品国产一区二区成人| 男人舔女人下体高潮全视频| 国内精品美女久久久久久| 国内精品一区二区在线观看| 欧美三级亚洲精品| 蜜臀久久99精品久久宅男| 午夜a级毛片| 亚洲国产欧洲综合997久久,| 一级二级三级毛片免费看| av卡一久久| 亚洲国产色片| 久久草成人影院| 亚洲欧美日韩高清在线视频| 亚洲无线观看免费| 变态另类成人亚洲欧美熟女| 欧美精品一区二区大全| 日韩欧美精品v在线| 色吧在线观看| 中国国产av一级| 午夜福利在线观看吧| 国产成人aa在线观看| 国产久久久一区二区三区| 成人一区二区视频在线观看| 日韩国内少妇激情av| 亚洲中文字幕日韩| 蜜臀久久99精品久久宅男| 久久久成人免费电影| 久久久久久国产a免费观看| 99国产精品一区二区蜜桃av| 在线免费十八禁| 久久精品国产99精品国产亚洲性色| 老司机福利观看| 日韩欧美国产在线观看| 午夜视频国产福利| 一个人看的www免费观看视频| 国产高潮美女av| 五月伊人婷婷丁香| 久久人人爽人人爽人人片va| 日本免费a在线| 亚洲精品国产av成人精品| 日韩一区二区三区影片| 久久久久久久久久久丰满| 赤兔流量卡办理| av福利片在线观看| 日韩精品青青久久久久久| 国产伦精品一区二区三区视频9| 国产av不卡久久| 亚洲七黄色美女视频| 99在线视频只有这里精品首页| 特大巨黑吊av在线直播| 成人一区二区视频在线观看| 老熟妇乱子伦视频在线观看| 一进一出抽搐gif免费好疼| 久久亚洲精品不卡| 97在线视频观看| av卡一久久| 国产日本99.免费观看| 国产综合懂色| 少妇熟女aⅴ在线视频| 国产女主播在线喷水免费视频网站 | 久久久成人免费电影| 深爱激情五月婷婷| 在现免费观看毛片| 免费黄网站久久成人精品| 啦啦啦韩国在线观看视频| 黄色一级大片看看| 22中文网久久字幕| 男女做爰动态图高潮gif福利片| 国产精品日韩av在线免费观看| av卡一久久| 久久99热6这里只有精品| 国产高清三级在线| 岛国毛片在线播放| 国产三级中文精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人福利小说| 亚洲高清免费不卡视频| 91aial.com中文字幕在线观看| 赤兔流量卡办理| 99视频精品全部免费 在线| 自拍偷自拍亚洲精品老妇| 色吧在线观看| 亚洲七黄色美女视频| 51国产日韩欧美| 啦啦啦韩国在线观看视频| 人人妻人人澡人人爽人人夜夜 | av视频在线观看入口| 国产黄色视频一区二区在线观看 | 国产精品永久免费网站| 亚洲久久久久久中文字幕| 国产精品嫩草影院av在线观看| 国产黄色小视频在线观看| 久久久精品欧美日韩精品| 天堂√8在线中文| 亚洲18禁久久av| 久久久久久伊人网av| 中文欧美无线码| 男女啪啪激烈高潮av片| 日韩欧美 国产精品| 春色校园在线视频观看| avwww免费| 亚洲欧美成人综合另类久久久 | 国产久久久一区二区三区| 日韩 亚洲 欧美在线| 韩国av在线不卡| 免费看av在线观看网站| 国产高清三级在线| 伊人久久精品亚洲午夜| 婷婷精品国产亚洲av| 国产三级中文精品| 爱豆传媒免费全集在线观看| 狠狠狠狠99中文字幕| 99热这里只有是精品在线观看| 岛国毛片在线播放| 亚洲无线在线观看| 午夜久久久久精精品| 欧美一区二区精品小视频在线| 欧美日韩乱码在线| 又黄又爽又刺激的免费视频.| 波多野结衣高清无吗| 两个人视频免费观看高清| 国产在线精品亚洲第一网站| 国产三级中文精品| 一个人观看的视频www高清免费观看| 国产精品麻豆人妻色哟哟久久 | 成人国产麻豆网| 深夜a级毛片| 九九热线精品视视频播放| 18禁在线无遮挡免费观看视频| 久久草成人影院| 欧美日本视频| 日韩欧美 国产精品| 国产精品人妻久久久影院| 一本久久中文字幕| 成年av动漫网址| 91狼人影院| 神马国产精品三级电影在线观看| 国产高清视频在线观看网站| 我要搜黄色片| 国产亚洲av片在线观看秒播厂 | 成人综合一区亚洲| 亚洲成人中文字幕在线播放| 少妇人妻精品综合一区二区 | 免费大片18禁| 天堂网av新在线| 亚洲精品久久久久久婷婷小说 | 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 久久精品综合一区二区三区| 国产综合懂色| 午夜a级毛片| 熟女人妻精品中文字幕| 网址你懂的国产日韩在线| 黄片无遮挡物在线观看| a级毛片免费高清观看在线播放| 99在线视频只有这里精品首页| av在线老鸭窝| 熟女电影av网| 国产成人福利小说| 国产精品一二三区在线看| 免费人成在线观看视频色| 亚洲国产高清在线一区二区三| 久久精品久久久久久噜噜老黄 | 成人无遮挡网站| 男人狂女人下面高潮的视频| 色综合色国产| 免费av毛片视频| 色播亚洲综合网| 久久久国产成人精品二区| 国产精品爽爽va在线观看网站| 春色校园在线视频观看| 国内久久婷婷六月综合欲色啪| 少妇熟女aⅴ在线视频| 国产午夜精品久久久久久一区二区三区| 一级毛片电影观看 | 亚洲av男天堂| 99久久精品国产国产毛片| 又爽又黄a免费视频| 国产白丝娇喘喷水9色精品| av国产免费在线观看| 久久韩国三级中文字幕| 欧美+日韩+精品| 国产高清视频在线观看网站| 免费看av在线观看网站| 欧美+亚洲+日韩+国产| 欧美xxxx黑人xx丫x性爽| 18+在线观看网站| 亚洲av.av天堂| a级毛色黄片| 全区人妻精品视频| 一个人观看的视频www高清免费观看| 成人午夜精彩视频在线观看| 亚洲一级一片aⅴ在线观看| 日本黄色片子视频| 欧美+亚洲+日韩+国产| 最后的刺客免费高清国语| 91狼人影院| 国产亚洲91精品色在线| 夜夜夜夜夜久久久久| 亚洲成av人片在线播放无| 成人特级黄色片久久久久久久| 亚洲欧美成人精品一区二区| 欧美日韩综合久久久久久| 国产 一区精品| 免费大片18禁| 美女 人体艺术 gogo| 国产日韩欧美在线精品| 婷婷精品国产亚洲av| 亚洲欧美成人综合另类久久久 | 欧美日韩在线观看h| 少妇高潮的动态图| 久久人人精品亚洲av| 亚洲最大成人av| 国产亚洲欧美98| 亚洲精华国产精华液的使用体验 | 午夜福利在线观看吧| 在线播放无遮挡| av又黄又爽大尺度在线免费看 | 99热这里只有是精品50| 久久久成人免费电影| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 夜夜爽天天搞| 亚洲最大成人中文| 综合色丁香网| 国产成人一区二区在线| 国产精品爽爽va在线观看网站| 国产精品麻豆人妻色哟哟久久 | 久久欧美精品欧美久久欧美| 国产老妇伦熟女老妇高清| 国语自产精品视频在线第100页| 听说在线观看完整版免费高清| 亚洲av成人精品一区久久| 国产一区二区在线观看日韩| 久久久久性生活片| 三级国产精品欧美在线观看| 国产亚洲精品久久久com| 国产精品国产高清国产av| 听说在线观看完整版免费高清| 日韩大尺度精品在线看网址| 国产精品,欧美在线| 成人亚洲精品av一区二区| 只有这里有精品99| 免费在线观看成人毛片| 亚洲av中文av极速乱| 男人和女人高潮做爰伦理| 亚洲色图av天堂| 偷拍熟女少妇极品色| 国产激情偷乱视频一区二区| 亚洲av成人av| 久久久久久久亚洲中文字幕| 亚洲一级一片aⅴ在线观看| 精品久久久久久成人av| 黄片wwwwww| 麻豆乱淫一区二区| 日产精品乱码卡一卡2卡三| 亚洲国产色片| 国产成人精品久久久久久| 18禁在线无遮挡免费观看视频| 国产精品乱码一区二三区的特点| 搡老妇女老女人老熟妇| 91狼人影院| 亚洲精品国产成人久久av| 国产色婷婷99| 国产亚洲91精品色在线| 免费观看精品视频网站| 爱豆传媒免费全集在线观看| 亚洲av不卡在线观看| 变态另类成人亚洲欧美熟女| 精品久久久久久久久av| 午夜精品在线福利| 美女xxoo啪啪120秒动态图| 国产在视频线在精品| 国产精品一区二区在线观看99 | 久久亚洲国产成人精品v| 色哟哟哟哟哟哟| av免费在线看不卡| av在线播放精品| 国产一级毛片在线| 亚洲,欧美,日韩| 又粗又爽又猛毛片免费看| 国产精品人妻久久久影院| 亚洲国产色片| 成人美女网站在线观看视频| 久久久欧美国产精品| 永久网站在线| 亚洲精品自拍成人| av又黄又爽大尺度在线免费看 | 色播亚洲综合网| 亚洲三级黄色毛片| 精品久久久久久久久亚洲| 少妇人妻一区二区三区视频| 一本一本综合久久| 美女脱内裤让男人舔精品视频 | 国产亚洲av片在线观看秒播厂 | 中国美女看黄片| 亚洲国产精品合色在线| 丝袜美腿在线中文| 国产毛片a区久久久久| 日日干狠狠操夜夜爽| 国产午夜精品论理片| 欧美性感艳星| 亚洲欧美日韩东京热| 偷拍熟女少妇极品色| 亚洲无线观看免费| 91在线精品国自产拍蜜月| 天天一区二区日本电影三级| 亚洲七黄色美女视频| 亚洲最大成人中文| 欧美成人a在线观看| 少妇高潮的动态图| 一边摸一边抽搐一进一小说| 亚洲中文字幕日韩| 精品一区二区三区视频在线| 亚洲激情五月婷婷啪啪| 夜夜夜夜夜久久久久| 亚洲乱码一区二区免费版| 日韩一区二区三区影片| 日本色播在线视频| 蜜桃久久精品国产亚洲av| 波多野结衣高清无吗| 18禁在线无遮挡免费观看视频| 亚洲久久久久久中文字幕| 最近的中文字幕免费完整| 欧美日韩国产亚洲二区| 男女下面进入的视频免费午夜| 99在线视频只有这里精品首页| 人妻制服诱惑在线中文字幕| 亚洲国产欧洲综合997久久,| 特大巨黑吊av在线直播| 午夜福利高清视频| 一个人观看的视频www高清免费观看| 天堂av国产一区二区熟女人妻| 热99在线观看视频| 国产黄色视频一区二区在线观看 | 两个人视频免费观看高清| 十八禁国产超污无遮挡网站| 最好的美女福利视频网| 最近手机中文字幕大全| av在线亚洲专区| 赤兔流量卡办理| 久久欧美精品欧美久久欧美| 亚洲国产精品成人综合色| 日本一本二区三区精品| 91av网一区二区| 成人av在线播放网站| 精品不卡国产一区二区三区| 一级二级三级毛片免费看| 黄色欧美视频在线观看| 我的女老师完整版在线观看| 中文字幕制服av| 午夜福利在线在线| 嘟嘟电影网在线观看| 欧美日韩精品成人综合77777| 成人亚洲精品av一区二区| 亚洲欧美成人综合另类久久久 | 国产精品不卡视频一区二区| 一夜夜www| 黄色配什么色好看| 亚洲va在线va天堂va国产| 国产成人精品久久久久久| 能在线免费观看的黄片| 欧美日本亚洲视频在线播放| 亚洲人成网站在线播| 国产成人91sexporn| 久久人妻av系列| 欧美精品国产亚洲| 亚洲欧美成人精品一区二区| 搡老妇女老女人老熟妇| 在线观看66精品国产| 国产成人影院久久av| 精品一区二区三区人妻视频| 国产日本99.免费观看| 日本色播在线视频| 亚洲天堂国产精品一区在线| 长腿黑丝高跟| 亚洲国产日韩欧美精品在线观看| 最近中文字幕高清免费大全6| 黄色日韩在线| 好男人在线观看高清免费视频| 成人毛片a级毛片在线播放| 久久久久性生活片| 日韩精品有码人妻一区| 亚洲在线观看片| 亚洲无线在线观看| 欧美激情久久久久久爽电影| 午夜福利视频1000在线观看| 精品人妻偷拍中文字幕| 成人鲁丝片一二三区免费| 91aial.com中文字幕在线观看| 成年av动漫网址| 国产精品免费一区二区三区在线| 可以在线观看毛片的网站| 青春草亚洲视频在线观看| 国产免费男女视频| 午夜福利在线观看吧| 身体一侧抽搐| 国产亚洲欧美98| 成人国产麻豆网| 国产成人a区在线观看| 亚洲无线在线观看| 少妇高潮的动态图| 老司机影院成人| 国产精品久久久久久亚洲av鲁大| 美女国产视频在线观看| 99热这里只有是精品50| 乱人视频在线观看| 黄色日韩在线| 精品人妻熟女av久视频| 97人妻精品一区二区三区麻豆| 亚洲美女搞黄在线观看| 一边亲一边摸免费视频| 国产伦理片在线播放av一区 | 欧美另类亚洲清纯唯美| 美女大奶头视频| 搡老妇女老女人老熟妇| 日本黄色片子视频| 国产三级在线视频| 如何舔出高潮| 久久久午夜欧美精品| 免费人成视频x8x8入口观看| 欧美潮喷喷水| 综合色av麻豆| 尾随美女入室| 听说在线观看完整版免费高清| 亚洲成av人片在线播放无| 青青草视频在线视频观看| 国产精品不卡视频一区二区| 一级毛片电影观看 | 在现免费观看毛片| 高清午夜精品一区二区三区 | 在线观看66精品国产| 色播亚洲综合网| 黄色一级大片看看| 日本av手机在线免费观看| 久久亚洲国产成人精品v| 一边亲一边摸免费视频| 午夜精品一区二区三区免费看| 国产精品,欧美在线| 尤物成人国产欧美一区二区三区| 在线观看一区二区三区| 一级av片app| 日本色播在线视频| 男人狂女人下面高潮的视频| 亚洲欧美日韩无卡精品| 久久99精品国语久久久| 美女 人体艺术 gogo| ponron亚洲| 国产片特级美女逼逼视频| 不卡一级毛片| 国产淫片久久久久久久久| 一区二区三区高清视频在线| 一本久久精品| 女同久久另类99精品国产91| 国产精品久久视频播放| 日韩国内少妇激情av| 亚洲精品国产av成人精品| 国产av在哪里看| av黄色大香蕉| 嫩草影院精品99| 亚洲经典国产精华液单| 日韩 亚洲 欧美在线| 最新中文字幕久久久久| 99久久精品一区二区三区| 欧美性猛交╳xxx乱大交人| 在线观看免费视频日本深夜| 国产日本99.免费观看| 国产精品久久久久久精品电影| 波多野结衣巨乳人妻| 天天躁日日操中文字幕| 国产三级中文精品| 99久久成人亚洲精品观看| 狂野欧美激情性xxxx在线观看| 黄色日韩在线| 韩国av在线不卡| 日日撸夜夜添| 麻豆成人av视频| 日本一二三区视频观看| 人人妻人人澡欧美一区二区| 高清在线视频一区二区三区 | 亚洲五月天丁香| 又粗又爽又猛毛片免费看| 国产久久久一区二区三区| 九色成人免费人妻av| 看免费成人av毛片| 能在线免费观看的黄片| 搡女人真爽免费视频火全软件| 欧美变态另类bdsm刘玥| 国模一区二区三区四区视频| 国产亚洲5aaaaa淫片| 中文字幕免费在线视频6| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av涩爱 | 国产一区二区三区av在线 | 观看美女的网站| 日韩国内少妇激情av| 丝袜美腿在线中文| 国产在视频线在精品| 国产熟女欧美一区二区| 午夜老司机福利剧场| 国产高清不卡午夜福利| 91午夜精品亚洲一区二区三区| 简卡轻食公司| 国产精品久久久久久亚洲av鲁大| 69av精品久久久久久| 国产高清激情床上av| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区久久| 黄色欧美视频在线观看| 色哟哟·www| 中出人妻视频一区二区| 老熟妇乱子伦视频在线观看| 少妇猛男粗大的猛烈进出视频 | 深爱激情五月婷婷| 亚洲无线观看免费| 综合色丁香网| 变态另类成人亚洲欧美熟女| 久久国内精品自在自线图片| 天堂中文最新版在线下载 | 波野结衣二区三区在线| 亚洲婷婷狠狠爱综合网| 欧美日韩综合久久久久久| 99riav亚洲国产免费| 亚洲中文字幕日韩| 亚洲四区av| 人妻夜夜爽99麻豆av| 亚洲国产精品成人久久小说 | 欧美成人一区二区免费高清观看| 在线观看66精品国产| 插阴视频在线观看视频| 你懂的网址亚洲精品在线观看 | 成人亚洲精品av一区二区| 欧美精品国产亚洲| 天天一区二区日本电影三级| 99riav亚洲国产免费| 国产亚洲5aaaaa淫片| 精品人妻一区二区三区麻豆| 久久久精品欧美日韩精品| 一级黄色大片毛片| 少妇熟女欧美另类| 丰满的人妻完整版| 日韩在线高清观看一区二区三区| 中国美女看黄片| 亚洲精品久久久久久婷婷小说 | 可以在线观看的亚洲视频| avwww免费| 天堂网av新在线| 91久久精品国产一区二区成人| 一区二区三区高清视频在线| 国产久久久一区二区三区| 成人特级av手机在线观看| av福利片在线观看| 国产中年淑女户外野战色| 丰满的人妻完整版| 国产精品不卡视频一区二区| 国产成人影院久久av| 国产精品久久久久久亚洲av鲁大| 啦啦啦啦在线视频资源| 国产午夜精品久久久久久一区二区三区| 天堂中文最新版在线下载 | 国产 一区 欧美 日韩|