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

    半封閉氣墊式調(diào)壓室模型優(yōu)化與水力特性研究

    2021-12-09 03:07:25歐傳奇劉德有周領(lǐng)
    人民長江 2021年10期

    歐傳奇 劉德有 周領(lǐng)

    摘要:為掌握半封閉氣墊式調(diào)壓室的水力特性,通過分析各類調(diào)壓室共性特征、改進已有模型氣體熱力學(xué)過程的基本假定、統(tǒng)一方程形式和提升離散精度,推導(dǎo)并建立了完全由顯式方程組成、易于編程電算、兼適用開敞式與氣墊式的半封閉氣墊式調(diào)壓室水力計算數(shù)學(xué)模型。在此基礎(chǔ)上,通過對主要計算參數(shù)取值、大波動水力性能、水位波動穩(wěn)定性能、結(jié)構(gòu)及運行控制開展研究和分析,系統(tǒng)地探討了半封閉氣墊式調(diào)壓室的水力特性。研究結(jié)果表明:該類調(diào)壓室能自適應(yīng)運行,性能良好,結(jié)構(gòu)簡單,改造方便,兼有開敞式和氣墊式調(diào)壓室的優(yōu)點。

    關(guān)鍵詞:半封閉氣墊式調(diào)壓室; 空氣阻抗調(diào)壓室; 水力特性; 調(diào)壓室水力計算模型

    中圖法分類號: TV732.5+6

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2021.10.036

    0引 言

    半封閉氣墊式調(diào)壓室亦稱空氣阻抗調(diào)壓室[1]、空氣閥式調(diào)壓室[2]、半壓氣式調(diào)壓室[3]、半氣墊式調(diào)壓室[4-5]、制氣式調(diào)壓室[6],是在封閉氣室頂部設(shè)置通氣閥(孔)與外部大氣相連的一種調(diào)壓設(shè)施。當室內(nèi)水位緩慢波動時,氣體能夠自由出入,不會引起室內(nèi)氣壓顯著變化;當水位劇烈波動時,空氣來不及充分流出氣室,室內(nèi)氣壓將大幅變化,從而對調(diào)壓室水位波動起到制動作用。在某些情況下,采用開敞式或氣墊式調(diào)壓室并不經(jīng)濟[4,7],若改用半封閉氣墊式調(diào)壓室可能會起到意想不到的效果。

    迄今為止,人們對半封閉氣墊式調(diào)壓室仍缺乏足夠的了解,研究投入比較缺乏,推廣應(yīng)用尚不重視,行業(yè)技術(shù)標準[8-9]也未曾關(guān)注。在水力計算數(shù)學(xué)模型方面,楊開林[5]曾在日本高佃政信的研究成果上作了一些改進,給出了相對完整的建模過程,但在基本假定和某些參數(shù)取值上仍有不足。鄭大瓊等[6]針對計算數(shù)學(xué)模型進行了推導(dǎo),也存在類似的問題,之后相關(guān)研究幾乎毫無進展。關(guān)于水力特性及適應(yīng)條件,早在1991年季奎等[4]就提出了半氣墊式調(diào)壓室方案,但直到2001年,鄭大瓊等[6]才將其與開敞式調(diào)壓室在長尾水系統(tǒng)中應(yīng)用的可行性進行了對比研究,并指出其部分性能特征可應(yīng)用于上游調(diào)壓室,這為后續(xù)研究提供了思路。不過,相關(guān)成果并未涉及到體型參數(shù)優(yōu)化、大小波動穩(wěn)定性能、結(jié)構(gòu)及運行控制特性,而且這些至關(guān)重要的內(nèi)容至今仍少有研究跟進。

    本文在分析已有典型模型的基礎(chǔ)上,通過改進基本假定、方程表達形式以及離散精度,研究完善該類調(diào)壓室的水力計算數(shù)學(xué)模型,以增加理論上的合理性、提高模型的通用性以及電算的方便性;并通過對主要參數(shù)取值、大波動水力性能、小波動穩(wěn)定性以及結(jié)構(gòu)和運行控制等方面的分析,系統(tǒng)探討該類調(diào)壓室的水力特性及其運行特性,以便為研究和推廣這類半封閉氣墊式調(diào)壓室提供參考依據(jù)。

    1計算數(shù)學(xué)模型及其改進

    1.1現(xiàn)有模型的主要缺陷

    以文獻[5-6]中的模型為代表,現(xiàn)有數(shù)學(xué)計算模型存在以下主要缺陷。

    (1) 未嚴格區(qū)分2個性質(zhì)不同的熱力學(xué)過程,理論欠嚴謹。假定氣體服從可逆的多邊關(guān)系,卻將室內(nèi)、外氣體交換過程與室內(nèi)氣體多變過程的氣體多變指數(shù)視為同一個量,未區(qū)分氣體流入、流出調(diào)壓室以及調(diào)壓室內(nèi)部變化2個不同的氣體熱力學(xué)過程。

    (2) 基于氣體多變方程推導(dǎo),氣體壓力計算式卻不連續(xù)。根據(jù)多變可逆的假定,氣體多變指數(shù)在閉區(qū)間[1.0,1.4]上取值(對應(yīng)熱力學(xué)過程為從等溫-等熵)都應(yīng)是適用的,但壓力計算式中,當氣體多變指數(shù)取1時,則出現(xiàn)0/0的情況。

    (3) 氣體多變指數(shù)取值簡單化,依據(jù)不足欠妥當。文中采用區(qū)間[1.0,1.4]上的平均值1.2進行計算,而實際氣體多變指數(shù)與之存在著較大差別[10-11],且取平均值的做法不能確保計算結(jié)果的安全。例如,室內(nèi)氣體壓力及水位控制值分別對應(yīng)其可能的最大值和最小值,取平均值計算不能兩者兼顧。

    (4) 需迭代求解隱式方程,而質(zhì)量流量對相對壓力的導(dǎo)數(shù)不連續(xù)。文獻[6]及經(jīng)典空氣閥等類似涉及進排氣的相關(guān)模型[5],通常采用牛頓等算法求解隱式方程。由于質(zhì)量流量對相對壓力的導(dǎo)數(shù)在相對壓力pr=1附近不連續(xù)(趨于負無窮大,而且通氣閥(孔)的CA值越大變化越劇烈,見圖1),模型求解穩(wěn)定性不好,計算速度也受影響。

    (5) 采用通氣孔面積計算質(zhì)量流量,未考慮其阻力特性。由于頂部通氣閥(孔)孔口尺寸一般較小,過流時存在著一定的阻力,而且在實際工程中,為了增加阻抗效果、減小調(diào)壓室的體積,通常會將通氣閥(孔)設(shè)置成進出阻力不同的情況,對于模擬氣體流入、流出不同的阻力特性具有現(xiàn)實需求和意義。

    1.2計算數(shù)學(xué)模型的改進

    為了增強模型的通用性,對于3種調(diào)壓室展開了分析(見圖2)。正如圖2所示,若不考慮水氣交界面以上的氣體熱力學(xué)效應(yīng),將水氣交界面處的氣壓看作已知量,則各類調(diào)壓室在控制方程的型式上沒有差別,借助于特征線法,可得出通用的幾個控制方程。

    1.3改進的計算模型主要特點

    (1) 基于最切合實際的假設(shè)推導(dǎo)出顯式氣體熱力學(xué)補充方程。對于第1個假設(shè),因調(diào)壓室與外界進行氣體交換流速大,氣體與外界來不及熱交換,最接近等熵過程,而且空氣熵流入、流出的過程完全不需要引入等溫假設(shè)[14]。對于第2個假設(shè),因?qū)嶋H系統(tǒng)氣體的熱力學(xué)過程可能發(fā)生部分熱交換而介于等溫和等熵之間,假定室內(nèi)氣體熱力學(xué)變化符合多變過程代替過去常用等溫的假定,不僅更符合實際,也考慮到了計算模型的通用性及實用性。在對實際工程進行設(shè)計時,可就氣體多變指數(shù)可能的上下限2種極端情況分別進行計算,從而獲得相對安全的設(shè)計參數(shù);必要時,也可根據(jù)一個波動周期內(nèi)氣體多變指數(shù)的變化,進行更切合實際的模擬。

    (2) 計算模型的通用性和適應(yīng)性顯著增強。① 在通用性方面,通過引入全局和局部2套坐標系,統(tǒng)一了正、負特征線方程形式,底部節(jié)點方程考慮了多管匯聚,兼顧底部各種復(fù)雜系統(tǒng)的編程計算[15];通過統(tǒng)一方程形式及設(shè)定CA值適配各類調(diào)壓室,當CA值取零時,適用于氣墊式調(diào)壓室,取較大CA值時,模型適用于開敞式調(diào)壓室,略作修改還可用于空氣閥。② 在適應(yīng)性方面,補充的熱力學(xué)方程為顯式表達式,并取二階近似保障精度,有利于改善算法的精度和穩(wěn)定性,提升計算速度。同時統(tǒng)一簡化了質(zhì)量流量計算表達式,有助于編程的模塊化和提高程序執(zhí)行維護效率。經(jīng)測試,采用Fortran編程,解決上述問題代碼可減少15%,且不存在迭代收斂問題。

    2計算參數(shù)取值敏感性分析

    根據(jù)計算數(shù)學(xué)模型,影響半封閉氣墊式調(diào)壓室的水力性能的主要參數(shù)有調(diào)壓室體型、通氣閥(孔)CA值、大氣溫度及氣體多變指數(shù)等。以四川省平武縣境內(nèi)陰平水電站為例[16],該水電站的尾水調(diào)壓室斷面面積為140.8 m2,高為14.8 m。

    2.1體型結(jié)構(gòu)尺寸的影響

    若調(diào)壓室體積相同但體型不同,則初始穩(wěn)定時刻室內(nèi)氣體總量也不同,由此影響水力瞬變的整個過程。為了便于分析,記斷高比為

    R=S/h(21)

    式中:S為調(diào)壓室斷面面積,m2;h為調(diào)壓室高度,m。

    取一系列斷高比進行過渡過程計算,計算結(jié)果如圖4所示。顯然,斷高比取值越大,對各控制參數(shù)越有利,而且當斷高比較小時,增加斷高比效果比較明顯。因此,在滿足布置要求的情況下,宜取“矮胖”體型,并可結(jié)合布置要求,通過該項分析來選擇最佳的體型。

    2.2通氣閥(孔)CA的取值

    通氣閥(孔)CA值指氣體流經(jīng)通氣閥(孔)的過流系數(shù)與氣體過流面積的乘積,包括流入時的CAin及流出時的CAout。CA值直接影響氣室與外界的氣體交換的速度,從而影響到調(diào)壓室內(nèi)的水位變化過程。氣體流入、流出的CA值是可以不同的,限于篇幅僅考慮CAin=CAout的情況,并取一系列的值進行過渡過程計算,計算結(jié)果如圖5所示。這也表明,調(diào)壓室內(nèi)最大、最小水深與尾水管進口的最小壓力對CA值均比較敏感。其中,調(diào)壓室內(nèi)最小水深隨CA值增大而降低,而對于調(diào)壓室內(nèi)最大水深,存在某一CA值(記CA1,約為調(diào)壓室斷面面積的1‰)可最大程度地降低最高水位;對于尾水管進口最小壓力,亦存在某一CA值(記CA2,約為調(diào)壓室斷面面積的3‰)可最大程度地改善此負壓值,但CA1不等于CA2。因此,CA的取值應(yīng)綜合考慮對這幾個參數(shù)的影響效果。此外,為對最小水深產(chǎn)生更為有利的影響,流入、流出的CA值可分別進行優(yōu)化。

    2.3大氣溫度變化與氣體多變指數(shù)取值

    由數(shù)學(xué)模型可知,大氣溫度對流入及流出的氣體質(zhì)量流量有影響,繼而影響整個過渡過程。不過,主要控制參數(shù)對大氣溫度的變化不敏感,詳細情況如表1所列。氣體多變指數(shù)m對整個過渡過程的影響體現(xiàn)在式(16)中。由表1可以看出:尾水管進口最小壓力隨m的增大而減小,即增大m對此不利;而室內(nèi)最低水位隨m的增大而升高,最高水位隨m的增大而降低,即水位波幅隨m的增大而減小,因此,m的增大對減小調(diào)壓室水位波幅是有利的。不過在等溫過程與等熵過程之間,主要控制參數(shù)對m取值并不敏感。

    3水力及運行特性分析

    3.1大波動水力性能

    結(jié)合設(shè)置條件,針對四川省平武縣境內(nèi)平陰水電站分別設(shè)置開敞式調(diào)壓室和半封閉氣墊式調(diào)壓室的情況(調(diào)壓室總體積相等,阻抗相同)進行了計算,計算結(jié)果如表2所列。計算結(jié)果表明:相同條件下,半封閉氣墊式調(diào)壓室的大波動水力性能較開敞式具有明顯的優(yōu)勢,不僅能夠改善尾水管進口真空度,還能有效地限制調(diào)壓室水位波幅,增加調(diào)壓室的安全水深,且當采用較優(yōu)的體型時,這種優(yōu)勢更加明顯。

    3.2水位波動穩(wěn)定性能

    水位波動穩(wěn)定性包括大波動穩(wěn)定性和小波動穩(wěn)定性。

    (1) 大波動穩(wěn)定性較開敞式調(diào)壓室要好。當半封閉氣墊式調(diào)壓室水面作大幅度波動時,其室內(nèi)氣壓將發(fā)生一定程度的變化,因氣室與外界進行了氣體交換,室內(nèi)氣壓變化程度較氣墊式調(diào)壓室相對要弱,因此,其大波動穩(wěn)定性理應(yīng)介于開敞式和氣墊式調(diào)壓室之間??紤]到任意調(diào)壓室水位大波動過程,若其衰減則必然進入小波動范疇,因此任何類型的調(diào)壓室波動越大則穩(wěn)定性越難以保障[17-19]。氣墊式、半封閉氣墊式調(diào)壓室涌波幅值較之開敞式調(diào)壓室小,其大、小波動穩(wěn)定體型尺寸差別較開敞式更小,因此半封閉氣墊式調(diào)壓室較開敞式調(diào)壓室具有更好的大波動穩(wěn)定性,在適用的低壓系統(tǒng)中兩者具有相近的穩(wěn)定性。

    (2) 小波動穩(wěn)定性與開敞式調(diào)壓室相當。由于半封閉氣墊式調(diào)壓室與外界大氣直接相連,在抑制水位波動的振幅的同時并不改變調(diào)壓室的穩(wěn)定水位,當半封閉氣墊式調(diào)壓室水面緩慢波動,調(diào)壓室內(nèi)外氣體能通過通氣孔自由出入,不會顯著影響氣室內(nèi)的氣壓,其小波動穩(wěn)定性與開敞式相當。對于氣墊式調(diào)壓室,普遍認為氣墊式穩(wěn)定性要差,但這是基于臨界穩(wěn)定斷面的比較結(jié)果,而氣墊式調(diào)壓室的穩(wěn)定性主要取決于氣體體積[20],穩(wěn)定斷面與穩(wěn)定氣體體積不具有可比性。

    3.3結(jié)構(gòu)及運行控制特性

    (1) 體型更優(yōu),布置靈活、施工方便。較開敞式調(diào)壓室,其水位波幅較小,同等水力約束條件下可減小調(diào)壓室體積。在保證相同安全水深的前提下,可適當抬高底板高程以減少基礎(chǔ)開挖,顯著降低工程投入的同時,也使調(diào)壓室布置更靈活,在某些情況下,可以減少尾水隧洞、廠房及調(diào)壓室之間的施工干擾。較氣墊式調(diào)壓室而言,其無需深埋地下。

    (2) 結(jié)構(gòu)要求相對寬松,閉氣要求易于實現(xiàn)。氣墊式調(diào)壓室利用氣墊抑制水位波幅的同時,必然會對室內(nèi)氣壓產(chǎn)生影響,從而導(dǎo)致對調(diào)壓室結(jié)構(gòu)強度及密閉性提出特殊要求。對于半封閉氣墊式調(diào)壓室,其室內(nèi)初始壓力僅為當?shù)卮髿鈮?,與氣墊式調(diào)壓室相比要小得多,而且在整個過渡過程中,氣體質(zhì)量的交換大大緩解了室內(nèi)氣體壓力的變化劇烈程度,因而不會產(chǎn)生過大的壓力變化,從而結(jié)構(gòu)上的要求會比氣墊式調(diào)壓室寬松得多,一般與開敞式接近,無需特殊處理。

    (3) 在改善阻抗設(shè)置上具有明顯的優(yōu)勢。就結(jié)構(gòu)特點而言,開敞式調(diào)壓室、半封閉氣墊式調(diào)壓室、氣墊式調(diào)壓室是調(diào)壓室頂部與氣體接觸面積逐漸減小至零的一個過程。從阻抗理論上講,相當于調(diào)壓室底部阻抗部分轉(zhuǎn)移至頂部氣體阻抗的過程,且因氣體阻抗引起的能量損失更小,對于一些阻抗設(shè)置困難的情況(阻抗孔兼做閘門孔尺寸受限,或底部阻抗孔面積過小流態(tài)不穩(wěn)定),是一種較好的選擇,較氣墊式調(diào)壓室而言,又不會顯得過“硬”。

    (4) 改造及運行方便,適宜增容改造工程。將開敞式調(diào)壓室改造成半封閉氣墊式調(diào)壓室,僅需頂部封閉并設(shè)置通氣孔或空氣閥即可,結(jié)構(gòu)尺寸不變性能可增加;較氣墊式調(diào)壓室而言,相當于自適應(yīng)等壓模式運行控制的氣墊式調(diào)壓室,無需補、排氣設(shè)備及相應(yīng)的監(jiān)測設(shè)施,地質(zhì)和防滲要求也大為降低,即降低了固定投資成本,又省去了運行成本。

    綜上所述,半封閉氣墊式調(diào)壓室利用低壓氣墊自適應(yīng)運行,結(jié)構(gòu)簡單,改造方便,兼得開敞式和氣墊式調(diào)壓室的優(yōu)點,無論是在大波動水力性能還是在穩(wěn)定性能方面均有良好表現(xiàn),較之設(shè)置條件相近的開敞式調(diào)壓室,在低壓系統(tǒng)中應(yīng)具有更好、更廣泛的適應(yīng)性。若設(shè)置為尾水調(diào)壓室,可利用尾水管進口最低負壓與調(diào)壓室最低水位出現(xiàn)的時間差,使得在不顯著降低尾水管進口最低負壓的同時,又能提高調(diào)壓室最低水位,從而減小工程量。

    4結(jié) 論

    (1) 改進了現(xiàn)有半封閉氣墊式調(diào)壓室數(shù)學(xué)模型,改進模型是基于更契合實際的假定,由顯式方程組成,求解方便,易于編程電算,兼適用開敞式與氣墊式調(diào)壓室,略作修改還可用于空氣閥。

    (2) 半封閉氣墊式調(diào)壓室水力性能主要受其體型結(jié)構(gòu)尺寸、通氣閥(孔)CA值影響,對大氣溫度、氣體多變指數(shù)取值(可能的變化范圍內(nèi))不敏感。在結(jié)構(gòu)尺寸方面,總體積不變時斷高比取值越大越有利,而且當斷高比小于1時,效果明顯,因此在滿足布置要求的情況下,宜選“矮胖”體型。通氣閥(孔)CA的值決定著利用氣墊的程度,當CA值小于5‰調(diào)壓室斷面面積時,對調(diào)壓室最高、最低水位與尾水管進口的最小壓力等主要參數(shù)影響明顯,因此可根據(jù)參數(shù)控制要求找出最佳CA值。

    (3) 半封閉氣墊式調(diào)壓室兼有良好的大波動水力性能和穩(wěn)定性能。相同條件下,其大波動水力性能較開敞式優(yōu)勢明顯,不僅能夠改善尾水管進口的真空度,還能有效地限制調(diào)壓室水位波幅,增加調(diào)壓室安全水深,而且當采用較優(yōu)體型時,其優(yōu)勢更加突出。就水位波動穩(wěn)定性而言,其大波動穩(wěn)定性較開敞式調(diào)壓室要好,小波動穩(wěn)定性與開敞式調(diào)壓室相當。

    (4) 半封閉氣墊式調(diào)壓室可利用低壓氣墊自適應(yīng)運行,結(jié)構(gòu)簡單,改造方便,兼得開敞式和氣墊式調(diào)壓室的優(yōu)點,在阻抗設(shè)置困難和增容改造需新增調(diào)節(jié)能力時優(yōu)勢明顯,較之設(shè)置條件相近的開敞式調(diào)壓室具有更好的適用性。

    參考文獻:

    [1]劉啟釗,陳家遠.空氣制動調(diào)壓室[J].華東水利學(xué)院學(xué)報,1980(2):57-76.

    [2]潘家錚,傅華.水工隧洞和調(diào)壓室:調(diào)壓室部分[M].北京:中國水利水電出版社,1981.

    [3]大連工學(xué)院,武漢水利電力學(xué)院.水電站建筑物[M].北京:水利出版社,1982.

    [4]季奎,劉文立.關(guān)于設(shè)置尾水調(diào)壓室條件的商討[J].鄭州工學(xué)院學(xué)報,1991,12(4):17-21.

    [5]楊開林.電站與泵站中的水力瞬變及調(diào)節(jié)[M].北京:中國水利水電出版社,2000.

    [6]鄭大瓊,王念慎.制氣式調(diào)壓室在帶有長尾水隧洞的水電站中的應(yīng)用[J].水利學(xué)報,2001(3):87-91.

    [7]劉德有,張健,索麗生.氣墊調(diào)壓室研究進展[J].水電能源科學(xué),2000(4):1-5.

    [8]國家能源局.水電站調(diào)壓室設(shè)計規(guī)范:NB/T 35021-2014[S].北京:中國電力出版社,2014.

    [9]中華人民共和國水利部.水利水電工程調(diào)壓室設(shè)計規(guī)范:SL 655-2014[S].北京:中國水利水電出版社,2014.

    [10]張健,索麗生,鄭源,等.氣墊調(diào)壓室內(nèi)氣體動態(tài)特性研究[J].水力發(fā)電學(xué)報,2004,23(4):97-101.

    [11]歐傳奇,劉德有,陳廣志.氣墊式調(diào)壓室氣體多變指數(shù)合理取值分析[J].人民黃河,2009,31(7):71-73.

    [12]歐傳奇,劉德有.復(fù)雜管系的水力計算[J].人民長江,2008,39(24):64-67.

    [13]WYLIE E B,STREETER V L,SUO L S.Fluid transients in systems[M].Nanjing:Prentice Hall,1993.

    [14]劉德有.輸水管道系統(tǒng)含大體積氣團瞬變流理論及應(yīng)用研究[D].南京:河海大學(xué),2008.

    [15]JAEGER C.水力不穩(wěn)定流在水力發(fā)電工程中的應(yīng)用[M].王樹人,譯.大連:大連工學(xué)院出版社,1987.

    [16]河海大學(xué).四川陰平水電站水力過渡過程計算分析[R].南京:河海大學(xué),2006.

    [17]季奎,馬躍先,王世強.調(diào)壓室大波動穩(wěn)定斷面研究[J].水利學(xué)報,1990(5):50-54.

    [18]馬躍先,季奎.氣壓式調(diào)壓室大波動穩(wěn)定斷面分析[J].水利學(xué)報,1992(12):41-47.

    [19]劉啟釗,彭守拙.水電站調(diào)壓室[M].北京:水利電力出版社,1993.

    [20]歐傳奇.氣墊式調(diào)壓室水力性能及運行控制研究[D].南京:河海大學(xué),2010.

    (編輯:趙秋云)

    Abstract:In order to master the hydraulic characteristics of semi-closed air cushion surge chamber(SACSC),a hydraulic calculation mathematical model of SACSC composed of explicit equations,easy to programming and suitable for both open and air cushion surge chamber,is derived and established by analyzing the common characteristics of all kinds of surge chambers.In the study,we improve the basic assumptions of gas thermodynamic process of existing models,unify the equation form and improve the discrete accuracy.On this basis,through the research and analysis of the main calculation parameters,large fluctuation hydraulic performance,water level fluctuation stability performance,structure and operation control,the hydraulic characteristics of SACSC are systematically discussed.The research shows that SACSC can operate adaptively,has good performance,simple structure and convenient reformation potential,possessing the advantages of open surge chamber and air cushion surge chamber.

    Key words:semi-closed air cushion surge chamber;air cushion surge chamber;hydraulics mathematical model of surges chamber;hydraulic characteristics

    麻豆成人午夜福利视频| 免费观看人在逋| 久久久久久国产a免费观看| 成年女人看的毛片在线观看| 免费一级毛片在线播放高清视频| 一进一出抽搐动态| 叶爱在线成人免费视频播放| av欧美777| 丰满的人妻完整版| 757午夜福利合集在线观看| 亚洲国产中文字幕在线视频| 毛片女人毛片| 亚洲国产欧美人成| 97人妻精品一区二区三区麻豆| 国产精品亚洲一级av第二区| 亚洲人成网站高清观看| 中文字幕精品亚洲无线码一区| 欧美激情在线99| 97超视频在线观看视频| 12—13女人毛片做爰片一| 日韩国内少妇激情av| 最新美女视频免费是黄的| 可以在线观看的亚洲视频| 欧美性猛交╳xxx乱大交人| 国产免费av片在线观看野外av| 午夜两性在线视频| 免费av观看视频| 欧美国产日韩亚洲一区| 十八禁人妻一区二区| 乱人视频在线观看| 国产高清激情床上av| 日韩人妻高清精品专区| 757午夜福利合集在线观看| 又爽又黄无遮挡网站| www.www免费av| 麻豆国产97在线/欧美| 亚洲精品一区av在线观看| 欧美最新免费一区二区三区 | 丰满的人妻完整版| 两个人的视频大全免费| 久久精品国产亚洲av涩爱 | 亚洲美女视频黄频| 成人av一区二区三区在线看| 内射极品少妇av片p| 国产成人a区在线观看| xxxwww97欧美| 国产亚洲av嫩草精品影院| 叶爱在线成人免费视频播放| 88av欧美| 男女做爰动态图高潮gif福利片| 一进一出抽搐gif免费好疼| 日韩国内少妇激情av| 日韩大尺度精品在线看网址| 色av中文字幕| 三级毛片av免费| 亚洲国产精品成人综合色| 久久久久久国产a免费观看| 国产免费av片在线观看野外av| 亚洲精品一区av在线观看| 久久精品91蜜桃| 久久久久国产精品人妻aⅴ院| 国产成人a区在线观看| 精品午夜福利视频在线观看一区| 嫁个100分男人电影在线观看| 淫秽高清视频在线观看| 在线a可以看的网站| 亚洲天堂国产精品一区在线| 中国美女看黄片| 免费看美女性在线毛片视频| 变态另类丝袜制服| 午夜免费观看网址| av国产免费在线观看| 嫩草影视91久久| 日本黄大片高清| 黄色成人免费大全| 国产免费男女视频| 中文字幕高清在线视频| 国产免费男女视频| 免费在线观看影片大全网站| 2021天堂中文幕一二区在线观| av欧美777| 精品电影一区二区在线| 亚洲人成网站在线播放欧美日韩| 熟妇人妻久久中文字幕3abv| 国产伦一二天堂av在线观看| 国产av麻豆久久久久久久| 成人18禁在线播放| 精品午夜福利视频在线观看一区| 免费在线观看亚洲国产| 国产精品美女特级片免费视频播放器| 久久性视频一级片| 成人无遮挡网站| 精品欧美国产一区二区三| 国产一区二区三区在线臀色熟女| 国产高清视频在线观看网站| 在线免费观看的www视频| 3wmmmm亚洲av在线观看| 精品人妻1区二区| 又粗又爽又猛毛片免费看| 精品不卡国产一区二区三区| 成年版毛片免费区| 婷婷精品国产亚洲av| 99久久九九国产精品国产免费| 特级一级黄色大片| 亚洲一区二区三区不卡视频| 在线观看舔阴道视频| 可以在线观看毛片的网站| 国产色婷婷99| 观看免费一级毛片| 日韩欧美国产在线观看| 欧美国产日韩亚洲一区| 午夜福利成人在线免费观看| 精品99又大又爽又粗少妇毛片 | 亚洲一区二区三区色噜噜| 久久香蕉精品热| 变态另类成人亚洲欧美熟女| 国产精品国产高清国产av| 成人永久免费在线观看视频| 在线国产一区二区在线| 国产精品99久久久久久久久| 欧美一区二区精品小视频在线| 精品国产三级普通话版| 亚洲成人中文字幕在线播放| 亚洲精品一区av在线观看| 2021天堂中文幕一二区在线观| 久久欧美精品欧美久久欧美| 国产高清视频在线播放一区| 1000部很黄的大片| 精品国产亚洲在线| 午夜两性在线视频| 国产视频一区二区在线看| 黄色成人免费大全| 麻豆成人av在线观看| 国产精品亚洲一级av第二区| 亚洲电影在线观看av| 日本黄色视频三级网站网址| 成年女人看的毛片在线观看| 精品国产亚洲在线| 亚洲一区二区三区色噜噜| 国产99白浆流出| 国产老妇女一区| 女同久久另类99精品国产91| 国产精品亚洲一级av第二区| 夜夜看夜夜爽夜夜摸| 国产极品精品免费视频能看的| 亚洲中文字幕一区二区三区有码在线看| 亚洲电影在线观看av| 我的老师免费观看完整版| 国产精品综合久久久久久久免费| 男人舔奶头视频| 日韩成人在线观看一区二区三区| 国产一区在线观看成人免费| 亚洲美女视频黄频| 高清日韩中文字幕在线| 国产高清视频在线观看网站| 19禁男女啪啪无遮挡网站| 成年女人永久免费观看视频| av在线天堂中文字幕| 国产高清有码在线观看视频| 母亲3免费完整高清在线观看| 18美女黄网站色大片免费观看| 男女视频在线观看网站免费| 精品久久久久久久末码| 桃红色精品国产亚洲av| 午夜精品久久久久久毛片777| 熟女人妻精品中文字幕| 久久久久性生活片| 91麻豆av在线| 亚洲国产精品久久男人天堂| 久久久久精品国产欧美久久久| 99热精品在线国产| 一本综合久久免费| 亚洲av二区三区四区| av在线蜜桃| 国产精品影院久久| 午夜福利在线在线| 91av网一区二区| 婷婷六月久久综合丁香| 国产精品免费一区二区三区在线| 亚洲精品国产精品久久久不卡| 国产一区二区三区在线臀色熟女| 在线观看美女被高潮喷水网站 | 夜夜看夜夜爽夜夜摸| 中文字幕人成人乱码亚洲影| 亚洲精品色激情综合| 波多野结衣高清作品| 哪里可以看免费的av片| 精品久久久久久久毛片微露脸| 免费大片18禁| 夜夜爽天天搞| 久久久久九九精品影院| 久久久久久国产a免费观看| 久久中文看片网| 哪里可以看免费的av片| 嫩草影院精品99| 18禁美女被吸乳视频| 国产精品久久久久久人妻精品电影| 亚洲五月天丁香| 毛片女人毛片| 国产视频内射| 亚洲avbb在线观看| 中文在线观看免费www的网站| 欧美一级a爱片免费观看看| 99久久99久久久精品蜜桃| 成人高潮视频无遮挡免费网站| 日本成人三级电影网站| 一区二区三区激情视频| ponron亚洲| 亚洲精品日韩av片在线观看 | 天堂√8在线中文| 18禁在线播放成人免费| 日本a在线网址| 99久久无色码亚洲精品果冻| a在线观看视频网站| 欧美成狂野欧美在线观看| 日韩有码中文字幕| 在线免费观看不下载黄p国产 | 成人国产一区最新在线观看| 男人舔女人下体高潮全视频| 88av欧美| 极品教师在线免费播放| 日本一二三区视频观看| 国产97色在线日韩免费| 日韩欧美国产在线观看| 久久精品国产亚洲av涩爱 | 麻豆成人午夜福利视频| 国产日本99.免费观看| 九色国产91popny在线| 香蕉久久夜色| 国产免费一级a男人的天堂| 亚洲成人精品中文字幕电影| 国产精品一及| av在线天堂中文字幕| 久久九九热精品免费| 国产精品久久久久久久久免 | 国产黄a三级三级三级人| 狠狠狠狠99中文字幕| 国产淫片久久久久久久久 | 免费看美女性在线毛片视频| 18美女黄网站色大片免费观看| 久99久视频精品免费| 高清在线国产一区| 中文字幕高清在线视频| av福利片在线观看| av天堂在线播放| 99久久精品热视频| 女生性感内裤真人,穿戴方法视频| 欧美在线一区亚洲| 美女 人体艺术 gogo| 日本三级黄在线观看| 国产亚洲精品综合一区在线观看| 免费无遮挡裸体视频| 大型黄色视频在线免费观看| 香蕉av资源在线| 两个人视频免费观看高清| 日韩亚洲欧美综合| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 免费无遮挡裸体视频| 欧美又色又爽又黄视频| 99热只有精品国产| 我要搜黄色片| 亚洲人成电影免费在线| 欧美中文综合在线视频| 搡老熟女国产l中国老女人| 三级国产精品欧美在线观看| 久久久久久久久大av| 精品免费久久久久久久清纯| 麻豆国产av国片精品| 国产久久久一区二区三区| 欧美乱码精品一区二区三区| 国产高潮美女av| 香蕉av资源在线| 狠狠狠狠99中文字幕| 久久午夜亚洲精品久久| 色尼玛亚洲综合影院| 国产亚洲精品一区二区www| 12—13女人毛片做爰片一| 久久亚洲精品不卡| 嫩草影院精品99| 亚洲人成电影免费在线| 99热这里只有精品一区| 在线观看舔阴道视频| 身体一侧抽搐| 男女做爰动态图高潮gif福利片| 18禁在线播放成人免费| 1024手机看黄色片| 国产综合懂色| 欧美激情久久久久久爽电影| 美女大奶头视频| 国产一区二区三区在线臀色熟女| 欧美日韩福利视频一区二区| 国内久久婷婷六月综合欲色啪| 久久精品国产综合久久久| 在线a可以看的网站| 一级a爱片免费观看的视频| 变态另类丝袜制服| 中国美女看黄片| 亚洲第一欧美日韩一区二区三区| 俄罗斯特黄特色一大片| 别揉我奶头~嗯~啊~动态视频| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站| 午夜精品一区二区三区免费看| 国产免费av片在线观看野外av| 亚洲欧美精品综合久久99| av中文乱码字幕在线| 99精品久久久久人妻精品| 成人特级黄色片久久久久久久| 嫩草影院入口| 免费av不卡在线播放| 男女之事视频高清在线观看| 小蜜桃在线观看免费完整版高清| 国产亚洲欧美在线一区二区| 亚洲在线自拍视频| 中文字幕人妻丝袜一区二区| 尤物成人国产欧美一区二区三区| 国产精品久久久人人做人人爽| 特大巨黑吊av在线直播| 亚洲国产精品999在线| 老司机午夜十八禁免费视频| 亚洲真实伦在线观看| 国产麻豆成人av免费视频| 我的老师免费观看完整版| 桃红色精品国产亚洲av| 成年女人看的毛片在线观看| 亚洲最大成人手机在线| 身体一侧抽搐| 麻豆成人午夜福利视频| 国产精品综合久久久久久久免费| 日韩有码中文字幕| 成人国产一区最新在线观看| 精品无人区乱码1区二区| 日韩免费av在线播放| 欧美色视频一区免费| 欧美三级亚洲精品| 成年人黄色毛片网站| 熟女电影av网| 亚洲最大成人手机在线| 亚洲一区二区三区色噜噜| 男人的好看免费观看在线视频| 免费av不卡在线播放| 国产精品久久久久久人妻精品电影| 伊人久久精品亚洲午夜| 亚洲七黄色美女视频| 熟妇人妻久久中文字幕3abv| 91在线精品国自产拍蜜月 | 午夜福利免费观看在线| 免费av观看视频| 九九热线精品视视频播放| 国产伦精品一区二区三区四那| 老汉色av国产亚洲站长工具| 国产欧美日韩精品一区二区| 久久久久久大精品| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 香蕉久久夜色| 日本a在线网址| 天天一区二区日本电影三级| 色播亚洲综合网| 麻豆一二三区av精品| avwww免费| 欧美中文日本在线观看视频| 亚洲自拍偷在线| 禁无遮挡网站| 一个人看的www免费观看视频| 91久久精品国产一区二区成人 | 乱人视频在线观看| 欧美+日韩+精品| 日本免费a在线| 内射极品少妇av片p| 日本免费a在线| 亚洲va日本ⅴa欧美va伊人久久| 精品免费久久久久久久清纯| 99久久九九国产精品国产免费| 欧美激情在线99| 手机成人av网站| 免费搜索国产男女视频| 国产视频一区二区在线看| 国产成人影院久久av| 熟女电影av网| 国产淫片久久久久久久久 | 国产成人系列免费观看| 午夜福利成人在线免费观看| avwww免费| 亚洲精品乱码久久久v下载方式 | 9191精品国产免费久久| 天美传媒精品一区二区| 老司机在亚洲福利影院| 国产精品久久视频播放| 亚洲中文字幕一区二区三区有码在线看| 亚洲成av人片免费观看| 国产乱人视频| 成人鲁丝片一二三区免费| 91av网一区二区| 成人亚洲精品av一区二区| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 欧美黑人欧美精品刺激| 真人做人爱边吃奶动态| 香蕉丝袜av| 亚洲欧美日韩东京热| 色播亚洲综合网| 在线a可以看的网站| 桃红色精品国产亚洲av| 老熟妇仑乱视频hdxx| 一进一出抽搐动态| 精品久久久久久,| 人妻久久中文字幕网| 国产高清激情床上av| 在线播放国产精品三级| 亚洲av成人av| 亚洲成av人片在线播放无| 我要搜黄色片| 国产精品98久久久久久宅男小说| 亚洲最大成人手机在线| 亚洲精品亚洲一区二区| 成年女人毛片免费观看观看9| 少妇丰满av| 精品人妻一区二区三区麻豆 | 久久伊人香网站| 丰满的人妻完整版| 99热精品在线国产| 18禁在线播放成人免费| 精品福利观看| 国产精品嫩草影院av在线观看 | 免费观看的影片在线观看| 1000部很黄的大片| av片东京热男人的天堂| 少妇丰满av| 黄色丝袜av网址大全| 深夜精品福利| 国产亚洲精品综合一区在线观看| 狠狠狠狠99中文字幕| av天堂中文字幕网| 十八禁人妻一区二区| 99久久九九国产精品国产免费| 亚洲内射少妇av| 淫妇啪啪啪对白视频| 99久久精品热视频| 最新美女视频免费是黄的| 一二三四社区在线视频社区8| 成人一区二区视频在线观看| 18禁在线播放成人免费| 国语自产精品视频在线第100页| 一本久久中文字幕| 老熟妇仑乱视频hdxx| 国产日本99.免费观看| 啦啦啦免费观看视频1| 黄色丝袜av网址大全| 长腿黑丝高跟| 国产精品亚洲美女久久久| 色播亚洲综合网| 日本黄色视频三级网站网址| 中文字幕精品亚洲无线码一区| 美女大奶头视频| 久久久久精品国产欧美久久久| 亚洲熟妇中文字幕五十中出| 香蕉久久夜色| 99国产精品一区二区蜜桃av| 色老头精品视频在线观看| 99riav亚洲国产免费| 久久久久久人人人人人| 免费看日本二区| 色综合亚洲欧美另类图片| 国产成人av激情在线播放| 国产精品爽爽va在线观看网站| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 午夜福利高清视频| 亚洲欧美精品综合久久99| 高清日韩中文字幕在线| 中文字幕av在线有码专区| 91在线精品国自产拍蜜月 | 热99re8久久精品国产| 美女大奶头视频| 中文字幕人妻熟人妻熟丝袜美 | 亚洲成av人片在线播放无| 51午夜福利影视在线观看| 一个人看的www免费观看视频| 精品欧美国产一区二区三| 又粗又爽又猛毛片免费看| 首页视频小说图片口味搜索| 高清日韩中文字幕在线| 午夜日韩欧美国产| 男女视频在线观看网站免费| 每晚都被弄得嗷嗷叫到高潮| 成年女人毛片免费观看观看9| 精品国产美女av久久久久小说| 国产精品野战在线观看| av视频在线观看入口| 动漫黄色视频在线观看| 成年免费大片在线观看| 高清日韩中文字幕在线| 日韩欧美在线乱码| 日本 欧美在线| 欧美一区二区精品小视频在线| 久久草成人影院| 久久中文看片网| 中文字幕久久专区| 免费观看人在逋| 久久香蕉精品热| 国内精品久久久久久久电影| 国产精品av视频在线免费观看| 欧美最黄视频在线播放免费| 亚洲七黄色美女视频| 成年女人看的毛片在线观看| 亚洲一区高清亚洲精品| 亚洲,欧美精品.| 中文资源天堂在线| 日韩高清综合在线| 午夜免费成人在线视频| 国产高清videossex| 久久久久久大精品| 婷婷精品国产亚洲av| 熟女人妻精品中文字幕| 亚洲片人在线观看| 看片在线看免费视频| 91久久精品国产一区二区成人 | 最后的刺客免费高清国语| 97碰自拍视频| 高清在线国产一区| 亚洲av电影在线进入| 亚洲专区国产一区二区| 亚洲人成伊人成综合网2020| 丰满的人妻完整版| 在线播放无遮挡| 一区二区三区免费毛片| 亚洲精品在线观看二区| 国产高清videossex| 亚洲av电影在线进入| 变态另类成人亚洲欧美熟女| x7x7x7水蜜桃| 91麻豆精品激情在线观看国产| 免费av观看视频| 一边摸一边抽搐一进一小说| 日韩欧美国产一区二区入口| 尤物成人国产欧美一区二区三区| 嫩草影院精品99| 在线观看免费视频日本深夜| 桃色一区二区三区在线观看| 国产私拍福利视频在线观看| 香蕉久久夜色| 久久午夜亚洲精品久久| 欧美在线黄色| 丁香六月欧美| 免费搜索国产男女视频| 欧美+亚洲+日韩+国产| 欧美国产日韩亚洲一区| 亚洲狠狠婷婷综合久久图片| 全区人妻精品视频| 女人被狂操c到高潮| 天堂动漫精品| 成人国产综合亚洲| 国产成人av教育| 一进一出好大好爽视频| 欧美日韩国产亚洲二区| 欧美成狂野欧美在线观看| 亚洲精品成人久久久久久| 免费人成在线观看视频色| 免费搜索国产男女视频| 久久欧美精品欧美久久欧美| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 观看美女的网站| av在线蜜桃| 亚洲男人的天堂狠狠| 九九热线精品视视频播放| 小说图片视频综合网站| xxxwww97欧美| 中文字幕人妻丝袜一区二区| 国产毛片a区久久久久| 成人性生交大片免费视频hd| 九九久久精品国产亚洲av麻豆| 在线播放无遮挡| 国产亚洲av嫩草精品影院| 18禁黄网站禁片免费观看直播| 窝窝影院91人妻| 久久天躁狠狠躁夜夜2o2o| 女人高潮潮喷娇喘18禁视频| 亚洲在线观看片| 亚洲黑人精品在线| 99国产极品粉嫩在线观看| 久久久久久久亚洲中文字幕 | 精品久久久久久成人av| 久久久国产精品麻豆| 欧美日韩综合久久久久久 | 欧美激情在线99| 嫩草影视91久久| 国产 一区 欧美 日韩| 91在线精品国自产拍蜜月 | 亚洲人成网站在线播| 国产精品爽爽va在线观看网站| 亚洲国产高清在线一区二区三| 欧美一级毛片孕妇| 成人特级av手机在线观看| 国产精品久久久久久久电影 | 中文字幕人妻丝袜一区二区| 少妇高潮的动态图| 亚洲欧美激情综合另类| 亚洲国产欧洲综合997久久,| 网址你懂的国产日韩在线| 成年免费大片在线观看| 国产国拍精品亚洲av在线观看 | 最近最新中文字幕大全免费视频| 亚洲 国产 在线| 国产精品 欧美亚洲| 香蕉丝袜av| 精品人妻偷拍中文字幕| 国产又黄又爽又无遮挡在线| av视频在线观看入口| 亚洲精品在线观看二区| 啦啦啦观看免费观看视频高清| 最近最新免费中文字幕在线| 亚洲欧美日韩无卡精品| 一区二区三区国产精品乱码|