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

    峽谷和埡口地形條件下輸電線路風(fēng)偏特性

    2021-05-17 05:32:10劉春城孫紅運(yùn)
    振動(dòng)與沖擊 2021年9期
    關(guān)鍵詞:風(fēng)偏埡口山脈

    劉春城,孫紅運(yùn)

    (東北電力大學(xué) 建筑工程學(xué)院,吉林 吉林 132012)

    我國(guó)是一個(gè)多山國(guó)家,兩條山脈并排形成峽谷和埡口是常見(jiàn)的山地地形。峽谷和埡口會(huì)形成的天然通道,受狹管效應(yīng)[1]的影響,風(fēng)速顯著增大。對(duì)于跨越峽谷和埡口的輸電線路而言,風(fēng)的加速效應(yīng)會(huì)導(dǎo)致導(dǎo)線風(fēng)荷載加大,對(duì)線路風(fēng)偏產(chǎn)生不利影響。目前針對(duì)山地地形中輸電線路抗風(fēng)偏設(shè)計(jì),我國(guó)規(guī)范相關(guān)規(guī)定比較簡(jiǎn)單,沒(méi)有詳細(xì)考慮到一些地形因素(如山脈長(zhǎng)度,山脈坡度,山脈間距)的影響。僅規(guī)定線路跨越峽谷、河道或位于暴露的山脊、頂峰等線路區(qū)段,其風(fēng)速值應(yīng)適當(dāng)增加,如無(wú)可靠資料時(shí),一般可按附近平地風(fēng)速增加10%[2]來(lái)考慮。這種過(guò)于籠統(tǒng)的規(guī)定會(huì)導(dǎo)致風(fēng)偏計(jì)算結(jié)果與實(shí)際相差較大,不利于線路安全運(yùn)行。

    國(guó)內(nèi)外學(xué)者針對(duì)山地風(fēng)場(chǎng)下輸電線路風(fēng)偏研究很少有涉及,大多數(shù)關(guān)于輸電線路風(fēng)偏響應(yīng)的研究都基于平坦地貌邊界層風(fēng)場(chǎng)[3-8]。Bullard等[9]研究了峽谷走向?qū)饬鬟\(yùn)動(dòng)趨勢(shì)的影響。Bowen等[10-12]研究了山體坡度,風(fēng)向角,山體粗糙度等因素對(duì)加速效應(yīng)的影響,但山體模型都是過(guò)于簡(jiǎn)單的山丘,并未考慮山體長(zhǎng)度。孫毅等[13]提出了新的脈動(dòng)風(fēng)及平均風(fēng)分布模型。李正良等[14]提出了復(fù)雜山地風(fēng)場(chǎng)平均及脈動(dòng)風(fēng)速均方根計(jì)算修正模型。李天昊等[15-19]較系統(tǒng)地研究并分析了峽谷長(zhǎng)度、山頂間距、山脈坡度3種地貌因素對(duì)平均風(fēng)加速效應(yīng)的影響。總的來(lái)看,目前對(duì)于山地線路風(fēng)偏的研究過(guò)于簡(jiǎn)單籠統(tǒng),并未考慮地形因素。因此獲得不同地形因素下峽谷和埡口輸電線路風(fēng)速特性,準(zhǔn)確計(jì)算線路的風(fēng)偏響應(yīng),對(duì)于架設(shè)在山地中輸電線路抗風(fēng)偏設(shè)計(jì)具有重要的借鑒和指導(dǎo)意義。

    本文建立了不同地形因素下峽谷和埡口風(fēng)場(chǎng),得到線路各點(diǎn)的平均橫向和豎向風(fēng)速隨山脈地形因素的變化規(guī)律。計(jì)算了線路在山地風(fēng)速和平地風(fēng)速下的風(fēng)偏響應(yīng),通過(guò)風(fēng)偏角增大百分比來(lái)量化山地地形對(duì)線路風(fēng)偏的影響。最后,用響應(yīng)面方法得到山脈長(zhǎng)度、坡度、間距三種因素和風(fēng)偏角增大百分比的二次回歸方程,并得到最不利地形的參數(shù)。研究結(jié)果為指導(dǎo)山地地形下輸電線路抗風(fēng)偏設(shè)計(jì)提供技術(shù)參考。

    1 峽谷和埡口地形風(fēng)場(chǎng)的數(shù)值模擬

    1.1 峽谷和埡口地形模型

    本文山體橫斷面選用余弦型山體。峽谷和埡口的三維模型及線路走向如圖1和圖2所示。其中:H表示山體高度(本文中山體高度統(tǒng)一取100 m),D表示山脈底部直徑,則山脈的坡度可表示為2H/D。

    圖1 峽谷模型Fig.1 Canyon model

    圖2 埡口模型Fig.2 Puerto model

    1.2 CFD數(shù)值模擬

    1.2.1 計(jì)算域及網(wǎng)格劃分

    利用ANSYS軟件中的SPACECLAIM建立峽谷和埡口模型,并用ICEM劃分網(wǎng)格。其中計(jì)算域高度為6倍模型高度,計(jì)算域長(zhǎng)度:山前取3倍模型長(zhǎng)度,山后取5倍模型長(zhǎng)度,寬度為3倍模型迎風(fēng)寬度[20-21]。計(jì)算域采用結(jié)構(gòu)化六面體網(wǎng)格,山體表面最大網(wǎng)格尺寸為5 m,增長(zhǎng)率為1.1,豎向網(wǎng)格首層高度取為1 m,增長(zhǎng)率為1.05。計(jì)算域模型及網(wǎng)格劃分如圖3所示。

    (a) 計(jì)算域模型

    1.2.2 入口邊界條件

    入口剖面包括平均風(fēng)剖面和湍流風(fēng)剖面,參照我國(guó)GB 50009—2012《建筑結(jié)構(gòu)荷載規(guī)范》[22]中使用的指數(shù)率模型,平均風(fēng)剖面的計(jì)算公式為

    (1)

    式中:U(z)為離地高度z處的平均風(fēng)速;Ur為參考高度zr處的平均風(fēng)速;α為地面粗糙度指數(shù)。

    湍流風(fēng)剖面計(jì)算包括湍動(dòng)能和耗散率兩個(gè)參數(shù),為實(shí)現(xiàn)入口剖面的自保持性,對(duì)湍流動(dòng)能k和耗散率ε分別采用式(2)和(3)來(lái)計(jì)算

    k(z)=0.3[I(z)×U(z)]2

    (2)

    (3)

    式中:k=0.42;Cμ=0.09;I(z)為z高度處湍流度;Lu為湍流積分尺度。

    這兩個(gè)參數(shù)均采用日本規(guī)范[23]中的公式:(Zg為梯度風(fēng)高度,B類(lèi)地貌下Zb=5 m)

    (4)

    (5)

    1.2.3 其余參數(shù)

    湍流模型采用Realizablek-ε湍流模型,選用非平衡壁面函數(shù),計(jì)算域入口為速度入口,出口為自由出口,兩側(cè)和頂部為對(duì)稱(chēng)邊界,平地和山體表面為固定壁面。其中Fluent軟件中粗糙度高度(roughness height)ks與空氣動(dòng)力學(xué)中粗糙度長(zhǎng)度Z0的換算關(guān)系為ks=20Z0,這里取山體和平地的粗糙高度分別為1 m和0.5 m。

    2 輸電線路參數(shù)和計(jì)算工況

    2.1 輸電線路參數(shù)

    線路走向如圖1或圖2中所示,采用“耐張塔-直線塔-直線塔-耐張塔”組成三跨線路,跨內(nèi)無(wú)轉(zhuǎn)角。線路電壓等級(jí)為500 kV,B類(lèi)地區(qū),基本風(fēng)速v0為30 m/s,耐張塔呼稱(chēng)高36 m,直線塔呼稱(chēng)高39 m,導(dǎo)線型號(hào)4×LGJ-400/35,絕緣子型號(hào)XWP2-160。為便于比較不同地形導(dǎo)線各點(diǎn)風(fēng)速,每檔導(dǎo)線均劃分45個(gè)單元,整條線路節(jié)點(diǎn)個(gè)數(shù)為136。

    2.2 導(dǎo)線風(fēng)荷載及風(fēng)偏響應(yīng)模擬

    在有限元軟件中模擬導(dǎo)線風(fēng)偏時(shí),將4分裂導(dǎo)線等效為一根導(dǎo)線,并將截面積、風(fēng)荷載等參數(shù)做等效化處理。將導(dǎo)線離散成多段索單元,計(jì)算每個(gè)單元所承受的風(fēng)荷載,作為靜力施加到節(jié)點(diǎn)上。每個(gè)單元所受風(fēng)荷載計(jì)算公式

    (6)

    式中:Fx為單元所受風(fēng)荷載;vh為單元所在高度h處的風(fēng)速;μsc為導(dǎo)線體型系數(shù),取1.1;βc為風(fēng)荷載系數(shù),計(jì)算風(fēng)偏時(shí)取1;d為導(dǎo)線直徑;l為單元長(zhǎng)度。

    2.3 計(jì)算工況

    為定量分析山地風(fēng)場(chǎng)下加速效應(yīng)對(duì)絕緣子風(fēng)偏位移的影響,本文設(shè)置2種計(jì)算工況。

    工況二:風(fēng)速vh取峽谷和埡口地形的橫向風(fēng)速。

    基于不同工況下直線塔絕緣子串風(fēng)偏角,定義風(fēng)偏角增大百分比為

    η=(θ2-θ1)/θ1

    (7)

    式中,θ1,θ2分別為工況1,2的直線塔絕緣子風(fēng)偏角計(jì)算結(jié)果。

    3 峽谷地形輸電線路風(fēng)偏特性

    3.1 模型工況

    峽谷指的是山頂間距(W)大于或等于山脈底部直徑(D)的地形。此種地形下研究的地貌因素包括山脈長(zhǎng)度,山頂距離,山脈坡度。所有涉及到的數(shù)值模擬地形工況如表1所示。

    表1 峽谷模型參數(shù)Tab.1 Canyon model parameters

    3.2 線路架設(shè)位置

    輸電線路在橫跨峽谷時(shí),可能架設(shè)在山脈中的任意位置。當(dāng)風(fēng)向?yàn)轫槏{谷時(shí),由于峽谷中風(fēng)場(chǎng)的復(fù)雜性,不同位置處的線路所承受的風(fēng)速也不同。因此,在山脈長(zhǎng)度300 m的峽谷中每隔50 m架設(shè)一條線路,共七條線路,如圖4所示。每條線路所在位置橫向風(fēng)速由圖5所示。圖中邊跨和靠近山頂處的線路風(fēng)速隨著到入口距離的增加先減小后增大。在距離入口200 m處達(dá)到最小,這是因?yàn)橛捎陲L(fēng)越過(guò)山頂后是一段平坦的山脈,且這兩部分的線路距離山體表面較近,空氣的流動(dòng)受山體表面摩擦的影響較大,導(dǎo)致線路上風(fēng)速的減小。但是在峽谷的后半段線路上風(fēng)速卻在增大,這是因?yàn)榭諝獯┻^(guò)峽谷后是開(kāi)闊的平地,有利于峽谷中空氣的流出,所以導(dǎo)致后半段的風(fēng)速會(huì)有小幅的提升。線路中跨中點(diǎn)處的風(fēng)速并非在入口處最大,而是在距入口50 m處最大,這是因?yàn)閸{谷中軸線上的風(fēng)速在入口處仍然處在一個(gè)加速的過(guò)程,在50 m處達(dá)到最大,而后由于山體的摩擦逐漸變小。綜合邊跨和中跨的風(fēng)速來(lái)看,在距離入口150~200 m處線路風(fēng)速最小,但是由于風(fēng)向的不確定性,建議在線路選線時(shí)在峽谷中部跨越峽谷更為保守。峽谷入口處風(fēng)的加速效應(yīng)最大,所以入口處線路是風(fēng)偏最不利的線路,也是受地形因素影響最大的位置。因此在研究峽谷地形因素的影響時(shí),線路的位置均架設(shè)在峽谷入口處。

    圖4 線路架設(shè)位置Fig.4 Line erection position

    3.3 山脈長(zhǎng)度的影響

    3.3.1 山脈長(zhǎng)度對(duì)風(fēng)場(chǎng)的影響

    不同的山脈長(zhǎng)度會(huì)對(duì)來(lái)流風(fēng)的加速效應(yīng)產(chǎn)生不同影響,處于峽谷入口的線路也會(huì)根據(jù)山脈長(zhǎng)度的不同而承受不同的風(fēng)速。不同山脈長(zhǎng)度入口處輸電線路所受橫向風(fēng)速,及平地風(fēng)速如圖6所示。由圖6可知,峽谷對(duì)風(fēng)的加速效應(yīng)很明顯,相同高度時(shí),峽谷中輸電線路所受風(fēng)速要高于平地風(fēng)速。隨著峽谷長(zhǎng)度的逐漸增大,峽谷對(duì)風(fēng)的加速效應(yīng)逐漸減小,線路所承受的風(fēng)速也逐漸變小。圖中山脈長(zhǎng)度在達(dá)到300 m后,間隔從100 m變成300 m時(shí)風(fēng)速的變化與之前相同,說(shuō)明山脈長(zhǎng)度的影響在逐漸減小,最后趨于穩(wěn)定。因此在山脈長(zhǎng)度達(dá)到一定程度后,再增加長(zhǎng)度也不會(huì)對(duì)風(fēng)速產(chǎn)生影響。不同山脈長(zhǎng)度輸電線路所受豎向風(fēng)速,及平地豎向風(fēng)速如圖7所示。由圖7可知,不管峽谷長(zhǎng)度多少,整條線路的豎向風(fēng)速都很小,最大處僅有3.8 m/s,這是因?yàn)楫?dāng)風(fēng)向?yàn)轫槏{谷時(shí),峽谷的受風(fēng)面積小,流動(dòng)的空氣可以順利從山脈兩側(cè)繞過(guò)去,不會(huì)產(chǎn)生大面積的爬升。因此,整條線路上的豎向風(fēng)加速效應(yīng)不顯著。

    圖6 不同山脈長(zhǎng)度橫向風(fēng)速Fig.6 Cross wind speeds of different mountain lengths

    圖7 不同山脈長(zhǎng)度豎向風(fēng)速Fig.7 Vertical wind speeds of different mountain lengths

    3.3.2 山脈長(zhǎng)度對(duì)線路風(fēng)偏的影響

    由圖6可知,峽谷風(fēng)速要大于平地風(fēng)速,在此基礎(chǔ)上,可以進(jìn)一步計(jì)算出不同山脈長(zhǎng)度下輸電線路風(fēng)偏角增大百分比的變化情況。如圖8所示,在山脈長(zhǎng)度為100 m時(shí),山頂直線塔的風(fēng)偏角要比平地大7.58%,當(dāng)山脈長(zhǎng)度增大時(shí),風(fēng)偏角增大百分比變小。說(shuō)明山脈長(zhǎng)度越長(zhǎng),峽谷地形對(duì)線路風(fēng)偏的影響越小。

    圖8 不同山脈長(zhǎng)度的風(fēng)偏角增大百分比ηFig.8 The η of different mountain lengths

    3.4 山頂間距的影響

    3.4.1 山頂間距對(duì)風(fēng)場(chǎng)的影響

    山頂間距也會(huì)對(duì)峽谷內(nèi)風(fēng)場(chǎng)產(chǎn)生影響,導(dǎo)致峽谷內(nèi)輸電線路所承受的風(fēng)速有所差異。不同山頂間距線路橫向風(fēng)速,及平地風(fēng)速如圖9所示。當(dāng)山頂間距變化時(shí),線路邊跨上的風(fēng)速不變,說(shuō)明峽谷中兩個(gè)山脈之間的距離并不影響山脈外側(cè)的風(fēng)場(chǎng)。而對(duì)于線路中跨,風(fēng)速隨著山頂間距的增大而減小。造成這種結(jié)果的原因一方面是由于山頂間距越大,線路檔距增大,導(dǎo)線弧垂越大,導(dǎo)線高度變小,風(fēng)速就會(huì)變小。另一方面是因?yàn)樯巾旈g距變大,導(dǎo)致峽谷中的空間更加開(kāi)闊,對(duì)風(fēng)的加速效應(yīng)減小,風(fēng)速變小。對(duì)比峽谷風(fēng)速和平地風(fēng)速可以發(fā)現(xiàn),不同山頂間距時(shí),數(shù)值模擬的風(fēng)速結(jié)果變化更大,印證了山脈間距會(huì)對(duì)峽谷內(nèi)風(fēng)場(chǎng)產(chǎn)生影響。不同山頂間距線路豎向風(fēng)速,及平地風(fēng)速如圖10所示。由圖可知,不論山頂間距多大,線路上的豎向風(fēng)速都很小,因此整條線路上的豎向風(fēng)加速效應(yīng)不顯著。

    圖9 不同山脈間距橫向風(fēng)速Fig.9 Cross wind speeds of different mountain spacing

    3.4.2 山頂間距對(duì)線路風(fēng)偏的影響

    進(jìn)一步計(jì)算山脈間距對(duì)線路風(fēng)偏的影響,如圖11所示。風(fēng)偏角增大百分比隨著間距的增大而減小,在兩條山脈山腳相接時(shí)最大。間距從300 m增加到500 m時(shí),η變化范圍并不大,說(shuō)明山脈間距對(duì)線路風(fēng)偏的影響較小。

    圖11 不同山脈間距的風(fēng)偏角增大百分比ηFig.11 The η of different mountain spacing

    3.5 山脈坡度的影響

    3.5.1 山脈坡度對(duì)風(fēng)場(chǎng)的影響

    本小節(jié)研究峽谷中不同山脈坡度對(duì)風(fēng)場(chǎng)的影響,本文是通過(guò)改變山脈底部直徑來(lái)改變坡度,而峽谷的定義是山頂間距W≥底部直徑D。因此,在改變山脈底部直徑的同時(shí)改變山頂間距,使W=D,即兩條山脈的山腳相接。不同山脈坡度線路橫向風(fēng)速,及平地風(fēng)速如圖12所示。線路邊跨的風(fēng)速隨著山脈坡度的減小而減小,這是因?yàn)槠露仍叫?,地?shì)越平緩,對(duì)風(fēng)的加速效應(yīng)越小。而山頂?shù)娘L(fēng)速隨著坡度的減小先增大后減小。出現(xiàn)這一結(jié)果是因?yàn)樯巾數(shù)娘L(fēng)速是由山脈體型和山脈坡度共同作用的。一方面坡度的變化會(huì)導(dǎo)致山脈體型的變化,坡度越大,山體越顯得細(xì)長(zhǎng),空氣就越容易在山的兩側(cè)產(chǎn)生繞流,而很難爬升到山頂。坡度越小,山體越顯得寬厚,空氣不易繞流,繼續(xù)爬升到山頂而導(dǎo)致山頂風(fēng)速變大,在底部直徑為400 m時(shí)達(dá)到最大。另一方面隨著坡度減小到一定程度,平緩的地勢(shì)對(duì)山頂?shù)娘L(fēng)速起到了主導(dǎo)作用,因此山頂風(fēng)速會(huì)越來(lái)越小。線路中跨在靠近山脈的部分與山頂風(fēng)速規(guī)律和原因基本相同。而線路中跨中間部分隨著山脈坡度的減小而減小,一方面由于山頂間距越大導(dǎo)致導(dǎo)線弧垂越大,導(dǎo)線高度越小,風(fēng)速越小。另一方面坡度越小,峽谷內(nèi)的空間越大,狹道的加速效應(yīng)越小導(dǎo)致風(fēng)速變小。不同山脈坡度線路豎向風(fēng)速,及平地風(fēng)速如圖13所示。不論山脈坡度多大,線路上的豎向風(fēng)速都很小,因此整條線路上的豎向風(fēng)加速效應(yīng)不顯著。

    圖12 不同山脈坡度下的橫向風(fēng)速Fig.12 Cross wind speed of different mountain slopes

    圖13 不同山脈坡度下的豎向風(fēng)速Fig.13 Vertical wind speed of different mountain slopes

    3.5.2 山脈坡度對(duì)線路風(fēng)偏的影響

    進(jìn)一步計(jì)算山脈坡度對(duì)線路風(fēng)偏的影響,如圖14所示。各個(gè)坡度的風(fēng)偏角增大百分比,隨著坡度的逐漸減小先變大后減小,在直徑為300 m時(shí)達(dá)到最大,說(shuō)明當(dāng)峽谷的山腳相接時(shí),山脈的坡度為2/3時(shí),對(duì)線路風(fēng)偏的影響最大。

    圖14 不同山脈坡度的風(fēng)偏角增大百分比ηFig.14 The η of different mountain slopes

    4 埡口地形輸電線路風(fēng)偏特性

    4.1 模型工況

    埡口地形指的是山頂距離W小于山底直徑D的地形,不同于峽谷之處是兩個(gè)山脈交疊在一起,在山脈之間形成一個(gè)低凹的地勢(shì)。此種地形下研究的地貌因素包括山脈長(zhǎng)度,山頂間距,山脈坡度。所有涉及到的數(shù)值模擬地形工況如表2所示。

    表2 埡口模型參數(shù)Tab.2 Puerto model parameters

    4.2 線路架設(shè)位置

    線路可在任意位置跨越埡口,風(fēng)向?yàn)轫槇嚎跁r(shí),埡口中不同位置的線路風(fēng)速也不同。在山脈長(zhǎng)300 m的埡口地形中,每隔50 m架設(shè)一條線路共7條線路,如圖15所示。每條線路所在位置橫向風(fēng)速如圖16所示。雖然埡口地形兩山間距相比峽谷變小了,但是圖中線路邊跨,山頂風(fēng)速規(guī)律與峽谷基本相同,都是先減小再增大,在距入口200 m時(shí)最小,其原因也同峽谷一樣,在此不再贅述。說(shuō)明兩山之間的間距并不影響山脈外側(cè)的風(fēng)場(chǎng)。線路中跨中間部分與峽谷不同,在入口處風(fēng)速最大,然后依次減小。原因是埡口地形兩條山脈距離很近,導(dǎo)致低凹處空間狹小。且導(dǎo)線架設(shè)處距離山體很近,使在兩側(cè)山體繞流的空氣更加迅速的吹向?qū)Ь€。來(lái)流風(fēng)吹向埡口時(shí),空氣迅速涌進(jìn)入口導(dǎo)致入口風(fēng)速突增,然后由于山體的摩擦平穩(wěn)減小。圖中可以看到中跨中間部分風(fēng)速會(huì)有一個(gè)小突起。這是由于山脈之間的交疊使低凹處形成一個(gè)夾角,由于此處空間過(guò)于狹小,使風(fēng)速略大于周?chē)以趭A角上方一定高度內(nèi)風(fēng)速都會(huì)受到影響。而且這個(gè)突起會(huì)隨著線路到入口距離的增加而變大,這是因?yàn)橥黄鹛帉?dǎo)線距離山體表面較遠(yuǎn),風(fēng)速受山體摩擦的影響沒(méi)有附近線路的大,導(dǎo)致風(fēng)速差距越來(lái)越大。由于邊跨和中跨導(dǎo)線到山體表面的距離相比于峽谷差別更小,所以風(fēng)速受山體表面摩擦作用相近,風(fēng)速變化的程度也想近。綜合來(lái)看,線路跨越埡口時(shí),在距入口200 m處的風(fēng)速最小,但是由于風(fēng)向的不確定性,架設(shè)在埡口中間更為保守。埡口在入口處風(fēng)的加速效應(yīng)最大,所以入口處線路是風(fēng)偏最不利的線路,也是受地形因素影響最大的位置。因此在研究埡口地形因素的影響時(shí),線路的位置均架設(shè)在埡口入口處。

    圖16 線路各點(diǎn)風(fēng)速Fig.16 Wind speed at each point of the line

    4.3 山脈長(zhǎng)度的影響

    4.3.1 山脈長(zhǎng)度對(duì)風(fēng)場(chǎng)的影響

    不同山脈長(zhǎng)度入口處輸電線路所受橫向風(fēng)速,及平地風(fēng)速如圖17所示。由圖可知,與峽谷地形風(fēng)速相比,峽谷的中跨風(fēng)速要比山頂略小,而埡口地形線路中跨風(fēng)速與山頂風(fēng)速基本相同。但是風(fēng)速變化規(guī)律和峽谷相同,都是隨著山脈長(zhǎng)度的增加而減小。圖中山脈長(zhǎng)度由100 m到200 m時(shí),風(fēng)速減小的程度略大于其他。說(shuō)明埡口的長(zhǎng)度不大時(shí),長(zhǎng)度的變化對(duì)風(fēng)的加速效應(yīng)更加敏感。不同山脈長(zhǎng)度入口處輸電線路所受豎向風(fēng)速,及平地風(fēng)速如圖18所示。埡口的豎向風(fēng)速相比于峽谷地形,雖然中跨導(dǎo)線更接近于山體,使豎向風(fēng)速略有增加,但是最大豎向風(fēng)速也僅有4.9 m/s,所以豎向風(fēng)加速效應(yīng)不顯著。

    圖17 不同山脈長(zhǎng)度橫向風(fēng)速Fig.17 Cross wind speeds of different mountain lengths

    圖18 不同山脈長(zhǎng)度豎向風(fēng)速Fig.18 Vertical wind speeds of different mountain lengths

    4.3.2 山脈長(zhǎng)度對(duì)線路風(fēng)偏的影響

    計(jì)算出埡口不同山脈長(zhǎng)度下輸電線路風(fēng)偏角增大百分比的變化情況。如圖19所示,在山脈長(zhǎng)度為100 m時(shí),山頂直線塔的風(fēng)偏角要比平地大7.92%,當(dāng)山脈長(zhǎng)度增大時(shí),風(fēng)偏角增大百分比變小。且曲線逐漸趨于平緩,隨著山脈長(zhǎng)度的增加,風(fēng)偏角增大百分比會(huì)趨近與一個(gè)固定值。

    圖19 不同山脈長(zhǎng)度的風(fēng)偏角增大百分比ηFig.19 The η of different mountain lengths

    4.4 山頂間距的影響

    4.4.1 山頂間距對(duì)風(fēng)場(chǎng)的影響

    不同山頂間距線路橫向風(fēng)速,及平地風(fēng)速如圖20所示。圖中可以看出線路邊跨的風(fēng)速并不隨著間距的變化而改變,說(shuō)明埡口地形山脈之間的間距不影響山脈外側(cè)的風(fēng)場(chǎng),這一點(diǎn)和峽谷地形相同。線路整個(gè)中跨的風(fēng)速隨著間距的增大而減小,出現(xiàn)這個(gè)結(jié)果的原因一部分是由于中跨的弧垂隨著檔距的增大而增大,導(dǎo)致導(dǎo)線高度減小風(fēng)速減小。但是對(duì)比平地風(fēng)速變化的程度,埡口的風(fēng)速變化范圍更大,這是由于地形變化導(dǎo)致的。間距的增大,導(dǎo)致埡口低凹處空間變開(kāi)闊,對(duì)風(fēng)的加速效應(yīng)減小。隨著埡口間距的逐漸變小,線路中跨的風(fēng)速越來(lái)越接近山頂?shù)娘L(fēng)速。這是因?yàn)殚g距的減小使埡口內(nèi)地表的高度不斷升高接近山頂,使埡口看起來(lái)更像是一個(gè)孤立山脈,所以導(dǎo)線風(fēng)速更接近于山頂。間距很小時(shí),導(dǎo)線中跨中間的風(fēng)速有一個(gè)小突起,但隨著間距的增大突起逐漸消失直至凹陷。這是由于山體二維輪廓形狀的特性導(dǎo)致的,間距越小,山體交叉點(diǎn)的高度越高,山體輪廓的斜率越大,導(dǎo)致低凹處的夾角越小。由4.2節(jié)可知,低凹處的夾角會(huì)導(dǎo)致其上方一定高度內(nèi)風(fēng)速略大于周?chē)L(fēng)速。但隨著間距的增大,山體輪廓逐漸平緩,交叉處的夾角也變大,且導(dǎo)線距離夾角的高度變高,所以,交叉處對(duì)于導(dǎo)線風(fēng)速的影響變小。不同山頂間距線路豎向風(fēng)速,及平地風(fēng)速如圖21所示。圖中可以看出山頂間距最小時(shí)線路豎向風(fēng)速最大,但僅有5.2 m/s,所以豎向風(fēng)加速效應(yīng)不顯著。

    圖20 不同山脈間距橫向風(fēng)速Fig.20 Cross wind speeds of different mountain spacing

    圖21 不同山脈間距豎向風(fēng)速Fig.21 Vertical wind speeds of different mountain spacing

    4.4.2 山頂間距對(duì)線路風(fēng)偏的影響

    不同山頂間距的風(fēng)偏角增大百分比如圖22,在山頂間距100 m時(shí)達(dá)到最大,隨著間距的增大,η逐漸變小。但是埡口的風(fēng)偏角增大百分比普遍大于峽谷地形,說(shuō)明埡口對(duì)線路風(fēng)偏的影響更加明顯。

    圖22 不同山脈間距的風(fēng)偏角增大百分比ηFig.22 The η of different mountain spacing

    4.5 山脈坡度的影響

    4.5.1 山脈坡度對(duì)風(fēng)場(chǎng)的影響

    對(duì)于埡口地形,山脈坡度的改變將影響整個(gè)埡口的形狀。因此為了排除不同山頂間距的影響,將每個(gè)模型的山頂間距取底面直徑的一半,即每個(gè)模型兩條山脈的交點(diǎn)都在50 m高度處。不同山脈坡度線路橫向風(fēng)速,及平地風(fēng)速如圖23所示。圖中橫坐標(biāo)采用和峽谷同樣的方法表示。線路邊跨風(fēng)速隨著坡度的減小而減小,這是因?yàn)槠露仍叫。貏?shì)越平緩,對(duì)風(fēng)的加速效應(yīng)越小。甚至在底面直徑900 m時(shí),有一部分線路的風(fēng)速要小于平地風(fēng)速。雖然埡口兩條山脈的距離相比峽谷來(lái)說(shuō)很近,但是在山頂和山頂附近風(fēng)速變化規(guī)律和峽谷相同,都是隨著坡度的減小先增大在減小,在底面直徑400 m時(shí)最大。說(shuō)明山脈距離的改變并不影響山脈頂部及附近的風(fēng)場(chǎng),這一點(diǎn)與3.4.1中的結(jié)論相同。但是對(duì)于線路中跨中間部分風(fēng)速,埡口和峽谷是有差異的。峽谷地形中此處風(fēng)速隨著坡度的增大而減小,埡口地形中是先增大再減小。這是因?yàn)閸{谷地形中兩條山脈并不相交,可以看成兩個(gè)孤立山脈。線路中跨中部可以看做處于兩條山脈中間,其下方和前方基本可以看做平地,所以此處的風(fēng)速只與狹道的大小有關(guān)。但是對(duì)于埡口,兩條山脈相互交疊,不可看做獨(dú)立的山脈,而是一個(gè)整體。線路中跨中部處于山體中央,此處的風(fēng)速受山脈坡度和山脈體型的共同影響。坡度越大,加速效應(yīng)越明顯,但是在山體坡度很大時(shí),山體的體型很瘦小,遮擋面積很小,所以風(fēng)速不會(huì)過(guò)大。隨著底面直徑的增加,遮擋面積變大,所以風(fēng)速普遍變大。但是隨著底面直徑不斷增大,山脈的地勢(shì)越來(lái)越平緩,接近平地時(shí),山脈坡度又起到了主導(dǎo)作用,風(fēng)速會(huì)越來(lái)越小。綜合山脈坡度和體型的影響,導(dǎo)致線路中跨中部風(fēng)速在底面直徑300~400 m時(shí)最大。不同山脈坡度線路豎向風(fēng)速,及平地風(fēng)速如圖24所示。圖中山脈坡度越大豎向風(fēng)速越大,但是最大也只有5.2 m/s,所以豎向風(fēng)加速效應(yīng)不顯著。

    圖23 不同山脈坡度下的橫向風(fēng)速Fig.23 Cross wind speed of different mountain slopes

    圖24 不同山脈坡度下的豎向風(fēng)速Fig.24 Vertical wind speed of different mountain slopes

    4.5.2 山脈坡度對(duì)線路風(fēng)偏的影響

    不同山脈坡度的風(fēng)偏增大百分比如圖25,隨著坡度的減小,η先增大再減小。當(dāng)兩條山脈山頂間距為底面直徑一半時(shí),山脈坡度為2/3時(shí),對(duì)線路風(fēng)偏的影響最大。

    圖25 不同山脈坡度的風(fēng)偏角增大百分比ηFig.25 The η of different mountain slopes

    5 地形因素與風(fēng)偏角增大百分比η的響應(yīng)面分析

    5.1 響應(yīng)面回歸方程的建立

    由于現(xiàn)實(shí)中山地地形復(fù)雜多變,無(wú)法通過(guò)數(shù)值模擬計(jì)算出每個(gè)地形的η。因此本文采用Box-Behnken方法進(jìn)行仿真工況設(shè)計(jì),以風(fēng)偏角增大百分比為響應(yīng)值(Y),選用山脈長(zhǎng)度(A),山脈坡度(B),山脈間距(C)為影響因素,其中山脈坡度用底面直徑表示。運(yùn)用Design-Expert軟件進(jìn)行多項(xiàng)擬合回歸分析,得到二次響應(yīng)面回歸方程為

    (8)

    5.2 仿真分析因素及水平及仿真結(jié)果分析

    仿真分析因素及水平見(jiàn)表3,響應(yīng)面分析仿真結(jié)果見(jiàn)表4,響應(yīng)面回歸模型方差分析結(jié)果見(jiàn)表5。

    表3 仿真分析因素及水平Tab.3 Simulation analysis factors and levels

    表4 仿真結(jié)果Tab.4 Simulation results

    表5 回歸模擬的方差分析Tab.5 Variance analysis of regression model

    通過(guò)模型顯著性檢驗(yàn)(P<0.05),表明該模型具有統(tǒng)計(jì)學(xué)意義。其自變量一次項(xiàng)A極顯著(P<0.000 1),C,A2,B2高度顯著(P<0.01)。失擬項(xiàng)用來(lái)表示所用模型與仿真擬合的程度,P值0.074 3>0.05,無(wú)失擬因素存在,因此可用該回歸方程代替仿真真實(shí)點(diǎn)對(duì)仿真結(jié)果進(jìn)行分析。相關(guān)系數(shù)R2=0.963 0,AdjR2=0.915 5,均大于0.8,說(shuō)明模型擬合優(yōu)度較好。

    5.3 各因素交叉相互作用分析

    圖26~28直觀顯示了各個(gè)因素交互作用的響應(yīng)值,圖中各因素的邊線越平坦則該因素對(duì)風(fēng)偏角增大百分比η的影響越小,邊線越陡則該因素對(duì)風(fēng)偏角增大百分比的影響越大[24]。圖中顯示,各因素的陡峭程度為A>C>B,與方差分析結(jié)果一致。

    5.4 最大風(fēng)偏角增大百分比

    結(jié)合數(shù)值模擬和響應(yīng)面分析結(jié)果,可以發(fā)現(xiàn)山脈長(zhǎng)度越短,山脈間距越近越不利于線路風(fēng)偏。由于本文設(shè)置了長(zhǎng)度和間距的下限,所以根據(jù)回歸模型得到最不利于線路風(fēng)偏的地形參數(shù)為:山脈長(zhǎng)度100 m,山脈坡度0.46,山脈間距100 m。最大風(fēng)偏角增大百分比為8.33%。

    6 結(jié) 論

    (1) 線路跨越峽谷和埡口時(shí),所受風(fēng)速隨著線路至山脈入口距離的增加,先減小后增大。所以山脈入口處風(fēng)加速效應(yīng)最為顯著,不利于線路安全運(yùn)行,建議線路從山脈中部跨越峽谷和埡口。

    (2) 峽谷和埡口中線路上的豎向風(fēng)速可以忽略不計(jì)。峽谷中η隨著山脈長(zhǎng)度,山脈間距的增大而減小,并且η會(huì)趨近與一個(gè)固定值。當(dāng)山脈重疊部分相同時(shí),η隨著山脈坡度的減小先增大再減小,在坡度為2/3時(shí)最大。

    (3) 埡口中η受地形因素的影響規(guī)律與峽谷相同,但是由于山脈間距更近,η值會(huì)普遍大于峽谷。

    (4) 山脈長(zhǎng)度、山脈坡度、山脈間距三種因素對(duì)風(fēng)偏角增大百分比影響的靈敏度從大到小依次為山脈長(zhǎng)度、山脈間距和山脈坡度。在本文研究的模型范圍內(nèi),山脈長(zhǎng)度100 m,山脈坡度為0.46,山脈間距100 m時(shí),風(fēng)偏角增大百分比達(dá)到8.33%。

    猜你喜歡
    風(fēng)偏埡口山脈
    亞丁大轉(zhuǎn)山 對(duì)世外桃源的最好詮釋
    它,就在那里
    500 kV緊湊型輸電線路覆冰厚度對(duì)導(dǎo)線布置的影響分析
    人,山脈和海洋
    贊馬關(guān)石埡口碉樓
    含笑花(2017年3期)2017-05-26 17:48:30
    500kV線路絕緣子串風(fēng)偏跳閘故障分析及對(duì)策
    廣西電力(2016年6期)2016-07-10 12:03:37
    往事
    連續(xù)多跨輸電線路動(dòng)態(tài)風(fēng)偏特征及計(jì)算模型
    冥王星上有一條山脈
    世界上海拔最高、最雄偉的山脈
    河北遙感(2014年4期)2014-07-10 13:54:58
    干丝袜人妻中文字幕| 又黄又爽又刺激的免费视频.| 真实男女啪啪啪动态图| 欧美日本视频| 亚洲精品自拍成人| 精品久久久久久久久av| 高清毛片免费观看视频网站| 久久精品国产亚洲av天美| 欧美激情国产日韩精品一区| 亚洲av中文av极速乱| 99在线人妻在线中文字幕| 最后的刺客免费高清国语| 91精品国产九色| 赤兔流量卡办理| 国产精品久久久久久久电影| 国产极品天堂在线| av.在线天堂| 亚洲av一区综合| 99热只有精品国产| 中文欧美无线码| 日本成人三级电影网站| 人妻制服诱惑在线中文字幕| 啦啦啦韩国在线观看视频| 观看免费一级毛片| 亚洲久久久久久中文字幕| av在线老鸭窝| 成人综合一区亚洲| 亚洲在线自拍视频| 国产成人午夜福利电影在线观看| 色5月婷婷丁香| 亚洲最大成人中文| 免费看光身美女| 欧美+日韩+精品| 少妇的逼水好多| 国产精品不卡视频一区二区| 久久久久久久久久成人| 亚洲精品乱码久久久久久按摩| 一边亲一边摸免费视频| 给我免费播放毛片高清在线观看| 亚洲不卡免费看| 搞女人的毛片| 欧美日本亚洲视频在线播放| 少妇裸体淫交视频免费看高清| 一进一出抽搐动态| 成人综合一区亚洲| 12—13女人毛片做爰片一| 高清日韩中文字幕在线| 插阴视频在线观看视频| 一本精品99久久精品77| 女同久久另类99精品国产91| 国产精品久久久久久精品电影| 午夜亚洲福利在线播放| 久久精品夜色国产| 精品久久久久久久人妻蜜臀av| 高清午夜精品一区二区三区 | 久久午夜亚洲精品久久| 直男gayav资源| 国产国拍精品亚洲av在线观看| 狠狠狠狠99中文字幕| 国产精品一二三区在线看| 欧美日本视频| 免费看a级黄色片| 国国产精品蜜臀av免费| 在线免费观看不下载黄p国产| 大又大粗又爽又黄少妇毛片口| 在线天堂最新版资源| 少妇被粗大猛烈的视频| 日韩成人av中文字幕在线观看| 欧美色欧美亚洲另类二区| 麻豆一二三区av精品| 亚洲国产高清在线一区二区三| 亚洲综合色惰| 日韩欧美精品免费久久| 直男gayav资源| 欧美最新免费一区二区三区| 免费电影在线观看免费观看| 18禁裸乳无遮挡免费网站照片| 国内久久婷婷六月综合欲色啪| 三级经典国产精品| 丰满的人妻完整版| 97热精品久久久久久| 亚洲人成网站在线播| 日韩国内少妇激情av| 欧美日本亚洲视频在线播放| 少妇猛男粗大的猛烈进出视频 | 22中文网久久字幕| 国产精品日韩av在线免费观看| 亚洲国产色片| 看黄色毛片网站| 日韩精品青青久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 18禁在线播放成人免费| or卡值多少钱| 国产亚洲精品久久久com| 能在线免费观看的黄片| 青春草视频在线免费观看| 级片在线观看| 久久精品国产亚洲av涩爱 | 成人国产麻豆网| 黄片无遮挡物在线观看| 九九在线视频观看精品| 全区人妻精品视频| 国产av在哪里看| 国产一级毛片七仙女欲春2| 波多野结衣高清无吗| 亚洲中文字幕日韩| 欧美不卡视频在线免费观看| 岛国毛片在线播放| 午夜福利在线观看免费完整高清在 | 婷婷六月久久综合丁香| 久久这里只有精品中国| 国产国拍精品亚洲av在线观看| 91在线精品国自产拍蜜月| 精品一区二区三区视频在线| 老司机福利观看| 尤物成人国产欧美一区二区三区| 亚洲人成网站高清观看| 国产精品.久久久| 亚洲国产精品sss在线观看| av又黄又爽大尺度在线免费看 | 丰满的人妻完整版| 一个人免费在线观看电影| 能在线免费观看的黄片| 1024手机看黄色片| 久久久成人免费电影| 国产一区二区三区在线臀色熟女| 美女内射精品一级片tv| 国产精品人妻久久久影院| 又爽又黄a免费视频| 中国国产av一级| 日本一二三区视频观看| 精品久久久久久久久久久久久| 又黄又爽又刺激的免费视频.| 九草在线视频观看| 欧美一区二区精品小视频在线| 色综合亚洲欧美另类图片| 亚洲精品影视一区二区三区av| 国产精品无大码| 国产精品精品国产色婷婷| 国产精品嫩草影院av在线观看| 久久精品国产鲁丝片午夜精品| 免费看光身美女| 日韩av不卡免费在线播放| 久久人人精品亚洲av| 国产高清视频在线观看网站| 床上黄色一级片| 床上黄色一级片| 国产亚洲精品久久久com| 国产91av在线免费观看| 啦啦啦韩国在线观看视频| 99riav亚洲国产免费| av视频在线观看入口| videossex国产| 中文字幕精品亚洲无线码一区| 国产91av在线免费观看| 真实男女啪啪啪动态图| 久久久久久久久中文| 91久久精品电影网| 成人国产麻豆网| 我要看日韩黄色一级片| 欧美日韩综合久久久久久| 九九久久精品国产亚洲av麻豆| 亚洲18禁久久av| 欧美激情在线99| a级毛色黄片| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 国产午夜精品久久久久久一区二区三区| 国产成人精品婷婷| 波多野结衣高清作品| 日韩av在线大香蕉| 国产中年淑女户外野战色| 免费黄网站久久成人精品| 国产高清三级在线| 欧美高清成人免费视频www| 搞女人的毛片| 99视频精品全部免费 在线| 亚洲成人久久性| 丰满的人妻完整版| 中国国产av一级| 可以在线观看的亚洲视频| 亚洲精品456在线播放app| 老女人水多毛片| 国产黄片美女视频| 日日啪夜夜撸| 成人国产麻豆网| 搡老妇女老女人老熟妇| 久久精品国产99精品国产亚洲性色| 国产在线男女| www.色视频.com| 美女 人体艺术 gogo| 99精品在免费线老司机午夜| 精品久久久久久成人av| 午夜免费激情av| 日韩一区二区三区影片| 乱码一卡2卡4卡精品| 国产免费男女视频| 亚洲第一区二区三区不卡| 精品一区二区免费观看| 一区二区三区四区激情视频 | 看免费成人av毛片| 亚洲真实伦在线观看| АⅤ资源中文在线天堂| 亚洲av免费高清在线观看| 深爱激情五月婷婷| 亚洲高清免费不卡视频| 亚洲成人久久爱视频| 99久久成人亚洲精品观看| 免费看光身美女| 日韩制服骚丝袜av| 天美传媒精品一区二区| 99热全是精品| 亚洲,欧美,日韩| 搡女人真爽免费视频火全软件| 一进一出抽搐动态| 在线观看免费视频日本深夜| 国产精品精品国产色婷婷| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 能在线免费观看的黄片| 国产熟女欧美一区二区| 日本爱情动作片www.在线观看| 日本一本二区三区精品| 日韩中字成人| 人妻系列 视频| 成人亚洲欧美一区二区av| 成人漫画全彩无遮挡| 中文在线观看免费www的网站| 欧美日本亚洲视频在线播放| 国产黄色视频一区二区在线观看 | 九九久久精品国产亚洲av麻豆| 99热6这里只有精品| 看十八女毛片水多多多| 久久亚洲精品不卡| 亚洲欧美清纯卡通| 久久久久久久久久成人| 成人特级av手机在线观看| 亚洲国产精品国产精品| 精华霜和精华液先用哪个| 国产精品久久久久久精品电影小说 | 欧美最新免费一区二区三区| 在线观看美女被高潮喷水网站| 久久久精品大字幕| 色噜噜av男人的天堂激情| 狂野欧美白嫩少妇大欣赏| 小说图片视频综合网站| 女人十人毛片免费观看3o分钟| 欧美变态另类bdsm刘玥| 亚洲国产欧洲综合997久久,| 中文字幕人妻熟人妻熟丝袜美| 日日啪夜夜撸| 熟女人妻精品中文字幕| 亚洲欧美精品自产自拍| 国产精品久久视频播放| 日韩,欧美,国产一区二区三区 | 热99re8久久精品国产| 国产白丝娇喘喷水9色精品| 国产精品爽爽va在线观看网站| 日本在线视频免费播放| 美女大奶头视频| 麻豆国产av国片精品| 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 老女人水多毛片| 国产一区二区在线观看日韩| 麻豆精品久久久久久蜜桃| 寂寞人妻少妇视频99o| 夜夜爽天天搞| 国产毛片a区久久久久| 久久久久久大精品| 国产 一区 欧美 日韩| 黄片wwwwww| 欧美xxxx性猛交bbbb| 久久精品91蜜桃| 日韩av在线大香蕉| 日本熟妇午夜| 日本av手机在线免费观看| 99热精品在线国产| av又黄又爽大尺度在线免费看 | 成人欧美大片| 一进一出抽搐gif免费好疼| 亚洲无线在线观看| 校园春色视频在线观看| 国产精品综合久久久久久久免费| 偷拍熟女少妇极品色| 中文在线观看免费www的网站| 国产精品久久久久久亚洲av鲁大| 91在线精品国自产拍蜜月| 十八禁国产超污无遮挡网站| 久久久久久九九精品二区国产| 有码 亚洲区| 51国产日韩欧美| 亚洲丝袜综合中文字幕| 国产午夜福利久久久久久| 精品午夜福利在线看| 99久久九九国产精品国产免费| 日产精品乱码卡一卡2卡三| 青春草亚洲视频在线观看| 秋霞在线观看毛片| 精品久久久久久久久久久久久| 身体一侧抽搐| 亚洲最大成人av| 91精品国产九色| 九九在线视频观看精品| videossex国产| 黄色配什么色好看| 最近中文字幕高清免费大全6| 看黄色毛片网站| 一边亲一边摸免费视频| 亚洲在线自拍视频| 国产高清视频在线观看网站| 天堂网av新在线| www.色视频.com| 校园春色视频在线观看| 久久精品久久久久久噜噜老黄 | 色5月婷婷丁香| 色综合站精品国产| 亚洲在线自拍视频| 亚洲av成人av| 成人一区二区视频在线观看| 长腿黑丝高跟| 亚洲欧美日韩高清在线视频| АⅤ资源中文在线天堂| 一级毛片aaaaaa免费看小| 久久精品国产清高在天天线| 别揉我奶头 嗯啊视频| 欧美bdsm另类| 高清毛片免费看| 欧美日韩精品成人综合77777| 久久精品久久久久久久性| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩高清专用| 边亲边吃奶的免费视频| 国产淫片久久久久久久久| 免费观看在线日韩| 精品欧美国产一区二区三| 女同久久另类99精品国产91| 亚洲精品国产av成人精品| 免费无遮挡裸体视频| 校园人妻丝袜中文字幕| 看免费成人av毛片| 少妇猛男粗大的猛烈进出视频 | 国产片特级美女逼逼视频| av.在线天堂| 一本久久中文字幕| 成人特级黄色片久久久久久久| 一区二区三区免费毛片| 99视频精品全部免费 在线| 欧美xxxx黑人xx丫x性爽| 欧美不卡视频在线免费观看| 精品久久久久久久末码| 插逼视频在线观看| 久久久色成人| 国产高潮美女av| 国产片特级美女逼逼视频| 国内精品美女久久久久久| 国内精品宾馆在线| 色哟哟哟哟哟哟| 久久午夜福利片| 亚洲经典国产精华液单| 18禁在线无遮挡免费观看视频| 99久久九九国产精品国产免费| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲av天美| 国产乱人视频| 国模一区二区三区四区视频| 亚洲无线观看免费| 亚洲av男天堂| 婷婷色av中文字幕| 成人亚洲欧美一区二区av| 九九热线精品视视频播放| 婷婷色av中文字幕| 亚洲一区二区三区色噜噜| 性插视频无遮挡在线免费观看| 少妇猛男粗大的猛烈进出视频 | 男女视频在线观看网站免费| 国产精品.久久久| 成人午夜高清在线视频| 日韩国内少妇激情av| 最近中文字幕高清免费大全6| 日本熟妇午夜| 六月丁香七月| 成人特级黄色片久久久久久久| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 男人和女人高潮做爰伦理| 两个人的视频大全免费| 精品一区二区三区视频在线| 亚洲精品日韩在线中文字幕 | 婷婷六月久久综合丁香| 真实男女啪啪啪动态图| 国产黄a三级三级三级人| 国产激情偷乱视频一区二区| 男女边吃奶边做爰视频| 中国美白少妇内射xxxbb| 综合色丁香网| 免费无遮挡裸体视频| 亚洲婷婷狠狠爱综合网| 成年版毛片免费区| 国产成人影院久久av| 久久精品91蜜桃| 久久久久网色| 神马国产精品三级电影在线观看| 男人狂女人下面高潮的视频| 久久这里只有精品中国| 国产成人a区在线观看| 黑人高潮一二区| 欧美潮喷喷水| 偷拍熟女少妇极品色| 伊人久久精品亚洲午夜| 一区二区三区免费毛片| 精品一区二区三区人妻视频| 亚洲成人久久爱视频| 2022亚洲国产成人精品| 女的被弄到高潮叫床怎么办| www日本黄色视频网| 亚洲图色成人| 国产男人的电影天堂91| 国产成人午夜福利电影在线观看| 69人妻影院| 人人妻人人澡人人爽人人夜夜 | 高清日韩中文字幕在线| 精品久久久噜噜| 床上黄色一级片| 18+在线观看网站| 18禁在线播放成人免费| 网址你懂的国产日韩在线| 性色avwww在线观看| 亚洲一区二区三区色噜噜| 可以在线观看的亚洲视频| 日产精品乱码卡一卡2卡三| 久久草成人影院| 精品久久久久久久久亚洲| 直男gayav资源| 深夜精品福利| 97超碰精品成人国产| 日本爱情动作片www.在线观看| 亚洲色图av天堂| 男插女下体视频免费在线播放| 日韩欧美 国产精品| 欧美xxxx性猛交bbbb| 中文欧美无线码| 内地一区二区视频在线| 亚洲自偷自拍三级| 97超碰精品成人国产| 一区二区三区四区激情视频 | 在线免费观看不下载黄p国产| 欧美日韩在线观看h| 亚洲不卡免费看| 神马国产精品三级电影在线观看| 久久草成人影院| 久久久成人免费电影| 国内精品久久久久精免费| 色综合站精品国产| 国内精品一区二区在线观看| 美女cb高潮喷水在线观看| 日本免费一区二区三区高清不卡| 国产亚洲欧美98| 久久亚洲精品不卡| 国产探花在线观看一区二区| 九九热线精品视视频播放| 日韩国内少妇激情av| 亚洲精品久久国产高清桃花| 日韩人妻高清精品专区| 女的被弄到高潮叫床怎么办| 黑人高潮一二区| 在线国产一区二区在线| 久久韩国三级中文字幕| 免费观看在线日韩| 麻豆精品久久久久久蜜桃| 午夜a级毛片| 亚洲第一电影网av| 国产精品三级大全| 亚洲人成网站高清观看| av在线观看视频网站免费| 一进一出抽搐动态| 悠悠久久av| 日韩高清综合在线| 波野结衣二区三区在线| 综合色丁香网| 亚洲无线观看免费| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 欧美一区二区国产精品久久精品| 国产精品电影一区二区三区| 国产成人精品婷婷| 国产色婷婷99| 欧美精品国产亚洲| 亚洲人成网站高清观看| 亚洲一级一片aⅴ在线观看| 精品久久久久久久久久久久久| 国产精品av视频在线免费观看| 午夜福利成人在线免费观看| av国产免费在线观看| 成人二区视频| 日韩精品有码人妻一区| 中文字幕免费在线视频6| 久久精品夜夜夜夜夜久久蜜豆| 国产精品永久免费网站| 久久婷婷人人爽人人干人人爱| 午夜精品国产一区二区电影 | 色噜噜av男人的天堂激情| 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| 精品无人区乱码1区二区| 久久欧美精品欧美久久欧美| 99热6这里只有精品| 波多野结衣高清无吗| 中文字幕精品亚洲无线码一区| 久久欧美精品欧美久久欧美| 日韩一区二区三区影片| av天堂在线播放| 久久99热6这里只有精品| 日日啪夜夜撸| 一级黄片播放器| 波多野结衣高清无吗| a级毛色黄片| 国产真实伦视频高清在线观看| 美女高潮的动态| av免费观看日本| 亚洲人成网站在线播| 国产精品,欧美在线| 国产亚洲av片在线观看秒播厂 | 欧美日韩综合久久久久久| 国产精品久久久久久精品电影小说 | 国产一级毛片七仙女欲春2| 国产亚洲欧美98| 久久久成人免费电影| 在线国产一区二区在线| 亚洲精品久久国产高清桃花| 日本五十路高清| 久久久久久久久久久丰满| 一级毛片电影观看 | 日本-黄色视频高清免费观看| 狂野欧美白嫩少妇大欣赏| 国产精品av视频在线免费观看| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久av不卡| 亚洲精品亚洲一区二区| 精品久久久久久久久久久久久| 中文字幕熟女人妻在线| 国产精品1区2区在线观看.| 午夜激情福利司机影院| 只有这里有精品99| 久久久a久久爽久久v久久| 国产精品女同一区二区软件| 一级黄色大片毛片| 精品无人区乱码1区二区| 免费av毛片视频| 中文资源天堂在线| 卡戴珊不雅视频在线播放| 日韩av不卡免费在线播放| 日本在线视频免费播放| a级毛片免费高清观看在线播放| .国产精品久久| 一本一本综合久久| 久久久久网色| 免费观看a级毛片全部| 啦啦啦观看免费观看视频高清| 女人十人毛片免费观看3o分钟| 午夜福利高清视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品国产成人久久av| 国产精品1区2区在线观看.| 99热全是精品| 国产精品一及| 国产伦理片在线播放av一区 | 伊人久久精品亚洲午夜| 简卡轻食公司| 少妇被粗大猛烈的视频| 亚洲丝袜综合中文字幕| 国产黄色小视频在线观看| 亚洲成人久久性| 欧洲精品卡2卡3卡4卡5卡区| 我要看日韩黄色一级片| 永久网站在线| 成年av动漫网址| 中文字幕熟女人妻在线| 日本-黄色视频高清免费观看| 女的被弄到高潮叫床怎么办| 一级毛片电影观看 | 日本黄大片高清| 国产91av在线免费观看| 久久久精品欧美日韩精品| 国产精品.久久久| 欧美xxxx性猛交bbbb| 色综合亚洲欧美另类图片| 麻豆成人av视频| 成人av在线播放网站| 精品一区二区免费观看| 国产伦一二天堂av在线观看| 国产精品久久久久久精品电影小说 | 中文字幕av成人在线电影| 精品人妻一区二区三区麻豆| 美女xxoo啪啪120秒动态图| 日本黄色视频三级网站网址| 晚上一个人看的免费电影| 精品久久国产蜜桃| 在线a可以看的网站| 亚洲成a人片在线一区二区| 国产成人精品婷婷| .国产精品久久| 婷婷色av中文字幕| 国产一区二区亚洲精品在线观看| 精品不卡国产一区二区三区| 中出人妻视频一区二区| 亚洲无线在线观看| 国产午夜福利久久久久久| 免费黄网站久久成人精品| 国产精品女同一区二区软件| 国产精品综合久久久久久久免费| 中文字幕人妻熟人妻熟丝袜美| 亚洲七黄色美女视频| 日日啪夜夜撸| 国产探花在线观看一区二区|