• <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
    男女边摸边吃奶| 蜜桃在线观看..| 91国产中文字幕| 国产精品久久久久成人av| 久久精品久久精品一区二区三区| 亚洲精品日韩在线中文字幕| 99久久精品一区二区三区| 春色校园在线视频观看| 三级国产精品欧美在线观看| 青春草亚洲视频在线观看| 午夜福利在线观看免费完整高清在| 中文字幕人妻熟人妻熟丝袜美| 国产成人免费无遮挡视频| 欧美精品人与动牲交sv欧美| 黄色欧美视频在线观看| 黄色怎么调成土黄色| 日韩熟女老妇一区二区性免费视频| 中国国产av一级| 男女边摸边吃奶| 国产不卡av网站在线观看| 黄色视频在线播放观看不卡| 丰满乱子伦码专区| 亚洲欧美成人综合另类久久久| 成年美女黄网站色视频大全免费 | av在线观看视频网站免费| 女性生殖器流出的白浆| 制服诱惑二区| 人人妻人人添人人爽欧美一区卜| 三级国产精品片| 欧美一级a爱片免费观看看| 久久女婷五月综合色啪小说| 国产精品秋霞免费鲁丝片| 丰满迷人的少妇在线观看| a级毛片黄视频| 久久久国产欧美日韩av| 肉色欧美久久久久久久蜜桃| 一本色道久久久久久精品综合| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 美女中出高潮动态图| 亚洲成人av在线免费| 亚洲色图 男人天堂 中文字幕 | 国产亚洲最大av| 大香蕉97超碰在线| 精品久久久噜噜| 国产男女超爽视频在线观看| 精品一区二区免费观看| 好男人视频免费观看在线| 晚上一个人看的免费电影| 在线播放无遮挡| 交换朋友夫妻互换小说| 一级a做视频免费观看| 免费观看的影片在线观看| av线在线观看网站| 久久久久久久久久人人人人人人| 日日撸夜夜添| av不卡在线播放| 高清毛片免费看| 99re6热这里在线精品视频| 亚洲av国产av综合av卡| 99九九线精品视频在线观看视频| 天堂俺去俺来也www色官网| 秋霞伦理黄片| 午夜日本视频在线| 老女人水多毛片| 中文天堂在线官网| 一边摸一边做爽爽视频免费| 一级毛片我不卡| 免费不卡的大黄色大毛片视频在线观看| 天天影视国产精品| 亚洲欧美一区二区三区国产| 亚洲成人手机| 成人漫画全彩无遮挡| 国产成人精品无人区| av不卡在线播放| 午夜久久久在线观看| av福利片在线| 日本黄大片高清| 制服诱惑二区| 国产黄色免费在线视频| 国产精品人妻久久久久久| 十八禁网站网址无遮挡| 亚洲成色77777| 国产精品偷伦视频观看了| 国产精品国产三级国产av玫瑰| 国产日韩一区二区三区精品不卡 | 日日摸夜夜添夜夜爱| a级毛片黄视频| 亚洲av在线观看美女高潮| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产女主播在线喷水免费视频网站| 视频在线观看一区二区三区| 极品人妻少妇av视频| 亚洲综合色网址| 最近中文字幕高清免费大全6| 美女cb高潮喷水在线观看| 建设人人有责人人尽责人人享有的| 亚洲美女搞黄在线观看| 久久婷婷青草| 日本黄色日本黄色录像| 寂寞人妻少妇视频99o| 欧美日韩视频精品一区| 高清av免费在线| 丝袜美足系列| a级毛色黄片| 久久久久久久亚洲中文字幕| 亚洲精品美女久久av网站| 在线精品无人区一区二区三| 一级a做视频免费观看| 我要看黄色一级片免费的| 日日啪夜夜爽| 成年人免费黄色播放视频| 国产精品一区二区三区四区免费观看| 蜜臀久久99精品久久宅男| 国产亚洲午夜精品一区二区久久| 国产精品久久久久成人av| 中文字幕免费在线视频6| 精品人妻在线不人妻| 国产在线视频一区二区| 满18在线观看网站| 亚洲精品乱码久久久久久按摩| 久久久午夜欧美精品| 国产精品 国内视频| 26uuu在线亚洲综合色| 亚洲欧洲日产国产| 插阴视频在线观看视频| 另类精品久久| 日本欧美视频一区| 日韩熟女老妇一区二区性免费视频| 母亲3免费完整高清在线观看 | 伦理电影免费视频| 日韩免费高清中文字幕av| 女的被弄到高潮叫床怎么办| 蜜桃在线观看..| 日本黄大片高清| 最近最新中文字幕免费大全7| 婷婷色综合www| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| 伦理电影免费视频| 日本欧美国产在线视频| 国产一级毛片在线| 少妇的逼水好多| 亚洲人成77777在线视频| 99热国产这里只有精品6| 亚洲精品av麻豆狂野| 亚洲精品第二区| 18+在线观看网站| 国产女主播在线喷水免费视频网站| 最黄视频免费看| 黑人猛操日本美女一级片| 免费观看的影片在线观看| 国产爽快片一区二区三区| 色5月婷婷丁香| 久久久午夜欧美精品| 亚洲欧洲日产国产| 亚洲四区av| 乱人伦中国视频| 亚洲国产日韩一区二区| 亚洲色图 男人天堂 中文字幕 | 欧美精品一区二区免费开放| 夫妻午夜视频| 亚洲国产色片| av视频免费观看在线观看| 日韩熟女老妇一区二区性免费视频| 人人妻人人澡人人爽人人夜夜| 曰老女人黄片| 中文乱码字字幕精品一区二区三区| 丝袜喷水一区| 男女免费视频国产| 午夜激情久久久久久久| 嫩草影院入口| 精品亚洲乱码少妇综合久久| 日本黄色日本黄色录像| 亚洲av中文av极速乱| 国产成人精品婷婷| 免费观看的影片在线观看| 国产乱人偷精品视频| 大香蕉久久网| 欧美日韩成人在线一区二区| 一本久久精品| 人妻制服诱惑在线中文字幕| 日本vs欧美在线观看视频| 午夜免费鲁丝| 欧美一级a爱片免费观看看| 人人妻人人澡人人爽人人夜夜| 亚洲美女视频黄频| a 毛片基地| 日韩,欧美,国产一区二区三区| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| 熟女电影av网| 亚洲,一卡二卡三卡| 中文字幕亚洲精品专区| 午夜精品国产一区二区电影| 亚洲精品日本国产第一区| 中文乱码字字幕精品一区二区三区| 日本黄色片子视频| 夜夜爽夜夜爽视频| 亚洲经典国产精华液单| av福利片在线| 青青草视频在线视频观看| 菩萨蛮人人尽说江南好唐韦庄| av网站免费在线观看视频| 少妇人妻 视频| 成人漫画全彩无遮挡| 午夜激情av网站| 国产老妇伦熟女老妇高清| 亚洲人成网站在线播| a级片在线免费高清观看视频| 国产永久视频网站| 国产亚洲午夜精品一区二区久久| av女优亚洲男人天堂| 在线观看免费视频网站a站| 91久久精品电影网| 亚洲精品色激情综合| 国产69精品久久久久777片| 熟女电影av网| 久久99热这里只频精品6学生| 亚洲欧美色中文字幕在线| 男人操女人黄网站| 国产成人aa在线观看| 性高湖久久久久久久久免费观看| 久久久久人妻精品一区果冻| 黄色欧美视频在线观看| 少妇人妻久久综合中文| 久久精品国产亚洲网站| 永久免费av网站大全| 国产av国产精品国产| 日韩精品免费视频一区二区三区 | 日本免费在线观看一区| 久久 成人 亚洲| 日韩欧美一区视频在线观看| 极品人妻少妇av视频| 国产日韩欧美视频二区| av网站免费在线观看视频| 男男h啪啪无遮挡| 国产成人午夜福利电影在线观看| 亚洲精品色激情综合| 如何舔出高潮| 欧美日韩成人在线一区二区| 国产精品一区www在线观看| 国产免费现黄频在线看| 欧美日韩视频高清一区二区三区二| 日本91视频免费播放| 哪个播放器可以免费观看大片| 日韩亚洲欧美综合| 精品国产露脸久久av麻豆| 国产综合精华液| 男人添女人高潮全过程视频| 日韩伦理黄色片| 51国产日韩欧美| .国产精品久久| 青春草亚洲视频在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲精品日韩av片在线观看| 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 久久99热这里只频精品6学生| 亚洲丝袜综合中文字幕| 亚洲国产av新网站| 中文字幕亚洲精品专区| 国产精品一国产av| 国产日韩欧美视频二区| 亚洲欧美一区二区三区黑人 | 国产高清国产精品国产三级| 久久精品熟女亚洲av麻豆精品| 免费播放大片免费观看视频在线观看| 国产一级毛片在线| 日本与韩国留学比较| a级毛片黄视频| 欧美最新免费一区二区三区| av国产精品久久久久影院| 少妇被粗大的猛进出69影院 | √禁漫天堂资源中文www| 精品亚洲乱码少妇综合久久| 欧美日韩视频精品一区| 久久久精品区二区三区| 国产黄色免费在线视频| 九色亚洲精品在线播放| 天天操日日干夜夜撸| 99九九在线精品视频| 亚洲久久久国产精品| 国产精品欧美亚洲77777| 国产乱人偷精品视频| 男人操女人黄网站| 亚洲精品成人av观看孕妇| 成年美女黄网站色视频大全免费 | 国产成人精品无人区| 91精品伊人久久大香线蕉| 人妻一区二区av| 2018国产大陆天天弄谢| 成人综合一区亚洲| 各种免费的搞黄视频| 精品99又大又爽又粗少妇毛片| 只有这里有精品99| 简卡轻食公司| 天天操日日干夜夜撸| 日韩免费高清中文字幕av| 少妇高潮的动态图| 中文字幕av电影在线播放| 欧美老熟妇乱子伦牲交| 国产精品欧美亚洲77777| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 亚洲成人一二三区av| 亚洲精品乱码久久久久久按摩| 亚洲四区av| 久久久国产精品麻豆| 欧美日韩视频高清一区二区三区二| 80岁老熟妇乱子伦牲交| 国产精品欧美亚洲77777| 亚洲av成人精品一区久久| 精品国产乱码久久久久久小说| 日本黄色片子视频| 在线观看免费日韩欧美大片 | 欧美人与性动交α欧美精品济南到 | tube8黄色片| 国产伦精品一区二区三区视频9| 欧美 日韩 精品 国产| 欧美精品人与动牲交sv欧美| 黄片播放在线免费| 热re99久久国产66热| 午夜91福利影院| 亚洲国产欧美日韩在线播放| 国产精品国产av在线观看| 日本色播在线视频| 亚洲成人手机| 午夜精品国产一区二区电影| 91精品三级在线观看| 少妇人妻久久综合中文| 亚洲国产精品成人久久小说| 91精品三级在线观看| 亚洲av男天堂| 大香蕉97超碰在线| 人妻一区二区av| 色吧在线观看| 18禁在线无遮挡免费观看视频| 国产成人精品一,二区| 母亲3免费完整高清在线观看 | 日日摸夜夜添夜夜爱| 美女中出高潮动态图| 精品久久久噜噜| 亚洲成人av在线免费| 国产熟女午夜一区二区三区 | 免费看av在线观看网站| 亚洲第一av免费看| 99热国产这里只有精品6| 少妇人妻精品综合一区二区| 欧美最新免费一区二区三区| 久久人人爽av亚洲精品天堂| av线在线观看网站| 成人毛片a级毛片在线播放| 天美传媒精品一区二区| 午夜精品国产一区二区电影| 熟女电影av网| 中文字幕最新亚洲高清| 97精品久久久久久久久久精品| 中文字幕人妻丝袜制服| 男女国产视频网站| 国产精品蜜桃在线观看| 亚洲激情五月婷婷啪啪| 热99国产精品久久久久久7| 岛国毛片在线播放| 亚洲av.av天堂| 超碰97精品在线观看| 高清黄色对白视频在线免费看| 视频中文字幕在线观看| 18禁观看日本| 纵有疾风起免费观看全集完整版| 校园人妻丝袜中文字幕| 亚洲人与动物交配视频| 纯流量卡能插随身wifi吗| 男女免费视频国产| 久久综合国产亚洲精品| 国产一区有黄有色的免费视频| av福利片在线| 最近2019中文字幕mv第一页| 成年人免费黄色播放视频| 国产国拍精品亚洲av在线观看| 久久精品国产a三级三级三级| 亚洲精品视频女| 人妻人人澡人人爽人人| 日本黄色日本黄色录像| 丰满乱子伦码专区| 午夜福利视频在线观看免费| 亚洲美女黄色视频免费看| 国产在视频线精品| 亚洲国产精品999| 3wmmmm亚洲av在线观看| 成人无遮挡网站| h视频一区二区三区| 成人黄色视频免费在线看| 色婷婷久久久亚洲欧美| 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 国产日韩欧美亚洲二区| 成人手机av| 美女国产高潮福利片在线看| 黄色一级大片看看| 黑人猛操日本美女一级片| av在线app专区| 欧美日韩在线观看h| 国产片特级美女逼逼视频| 亚洲经典国产精华液单| 亚洲国产精品一区三区| 欧美97在线视频| 国产国语露脸激情在线看| 一边摸一边做爽爽视频免费| 国产免费福利视频在线观看| 午夜福利视频精品| 中国美白少妇内射xxxbb| 大码成人一级视频| 中文字幕亚洲精品专区| 极品人妻少妇av视频| 日本91视频免费播放| 大码成人一级视频| 午夜视频国产福利| 色视频在线一区二区三区| 国产成人精品福利久久| av国产久精品久网站免费入址| 国产精品久久久久久久电影| 精品人妻熟女av久视频| 一区二区三区精品91| 婷婷色av中文字幕| av福利片在线| 91成人精品电影| 午夜激情久久久久久久| 欧美老熟妇乱子伦牲交| 久久 成人 亚洲| 精品国产一区二区三区久久久樱花| 大香蕉97超碰在线| 母亲3免费完整高清在线观看 | 欧美日韩精品成人综合77777| 纯流量卡能插随身wifi吗| 日韩制服骚丝袜av| 大又大粗又爽又黄少妇毛片口| 一级a做视频免费观看| 91精品一卡2卡3卡4卡| 青春草亚洲视频在线观看| 日本黄色日本黄色录像| 大片免费播放器 马上看| 狠狠婷婷综合久久久久久88av| a级毛片黄视频| 国产极品天堂在线| 中文字幕制服av| 亚洲av男天堂| 久热这里只有精品99| 国产精品女同一区二区软件| 日韩熟女老妇一区二区性免费视频| 插阴视频在线观看视频| 亚洲精品国产色婷婷电影| a级毛片在线看网站| 大香蕉97超碰在线| 伊人久久精品亚洲午夜| 少妇精品久久久久久久| 国产精品一国产av| 欧美3d第一页| 欧美日韩在线观看h| 一边摸一边做爽爽视频免费| 午夜福利网站1000一区二区三区| 国产成人a∨麻豆精品| 午夜福利,免费看| 丰满乱子伦码专区| 日韩电影二区| 精品午夜福利在线看| 日本91视频免费播放| 色吧在线观看| 国产亚洲精品第一综合不卡 | 成人国产av品久久久| 天堂中文最新版在线下载| 国产一区二区三区av在线| 日本爱情动作片www.在线观看| 国产亚洲精品第一综合不卡 | av线在线观看网站| av在线观看视频网站免费| 亚洲精品日韩av片在线观看| 亚洲国产精品成人久久小说| 一级毛片黄色毛片免费观看视频| 亚洲精品av麻豆狂野| 亚洲无线观看免费| av天堂久久9| 另类亚洲欧美激情| 日本猛色少妇xxxxx猛交久久| 欧美精品高潮呻吟av久久| 色吧在线观看| 国产一区有黄有色的免费视频| 精品一区在线观看国产| 成人免费观看视频高清| 五月天丁香电影| 亚洲美女黄色视频免费看| 十八禁高潮呻吟视频| 男女边吃奶边做爰视频| av有码第一页| 欧美日韩av久久| 哪个播放器可以免费观看大片| 国产亚洲av片在线观看秒播厂| 国产片内射在线| 国产精品.久久久| 国产片特级美女逼逼视频| 国产精品无大码| 新久久久久国产一级毛片| 黑人猛操日本美女一级片| 岛国毛片在线播放| 亚洲国产精品成人久久小说| 五月伊人婷婷丁香| 日韩视频在线欧美| av国产精品久久久久影院| 欧美xxⅹ黑人| 国产一级毛片在线| 亚洲国产欧美在线一区| 十八禁网站网址无遮挡| 国产精品熟女久久久久浪| 中文欧美无线码| 老女人水多毛片| 国产日韩一区二区三区精品不卡 | 综合色丁香网| 日本-黄色视频高清免费观看| 日本黄色片子视频| 97超视频在线观看视频| 亚洲,欧美,日韩| 高清欧美精品videossex| 在线观看免费日韩欧美大片 | 狂野欧美激情性xxxx在线观看| 国产免费现黄频在线看| 狂野欧美白嫩少妇大欣赏| 少妇人妻 视频| 国产av一区二区精品久久| 91精品伊人久久大香线蕉| 久久热精品热| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品成人av观看孕妇| 亚洲怡红院男人天堂| av视频免费观看在线观看| 嫩草影院入口| 最新的欧美精品一区二区| 精品亚洲乱码少妇综合久久| 国产精品一区二区在线不卡| 九九在线视频观看精品| 日本黄色片子视频| av又黄又爽大尺度在线免费看| 欧美+日韩+精品| 一本一本综合久久| 美女主播在线视频| av专区在线播放| 男男h啪啪无遮挡| 男女边摸边吃奶| 精品亚洲乱码少妇综合久久| 青春草亚洲视频在线观看| 在线天堂最新版资源| 啦啦啦啦在线视频资源| 自拍欧美九色日韩亚洲蝌蚪91| 成年人免费黄色播放视频| 男人爽女人下面视频在线观看| 欧美bdsm另类| 赤兔流量卡办理| 欧美日韩国产mv在线观看视频| 久久99精品国语久久久| 午夜影院在线不卡| 高清不卡的av网站| 日韩av在线免费看完整版不卡| 国产乱人偷精品视频| 一级二级三级毛片免费看| 欧美三级亚洲精品| 美女视频免费永久观看网站| 精品久久久久久久久av| 亚洲人成77777在线视频| a级毛色黄片| 国产一区亚洲一区在线观看| 成年人免费黄色播放视频| 午夜激情av网站| 日本vs欧美在线观看视频| 人妻一区二区av| 国产老妇伦熟女老妇高清| 成人国产麻豆网| 97超碰精品成人国产| 我要看黄色一级片免费的| 3wmmmm亚洲av在线观看| 亚洲人成网站在线观看播放| 日韩免费高清中文字幕av| 一区在线观看完整版| 久久这里有精品视频免费| 好男人视频免费观看在线| 91国产中文字幕| 大又大粗又爽又黄少妇毛片口| 热re99久久国产66热| 国产精品偷伦视频观看了| 成人综合一区亚洲| 嫩草影院入口| 亚洲美女黄色视频免费看| 成年女人在线观看亚洲视频| 国产免费一区二区三区四区乱码| 街头女战士在线观看网站| 久久精品熟女亚洲av麻豆精品| 精品卡一卡二卡四卡免费| 国产精品无大码| 99久久精品一区二区三区| 精品久久久精品久久久| 欧美人与性动交α欧美精品济南到 | 亚洲欧美成人综合另类久久久| 午夜福利视频在线观看免费| 亚洲欧洲日产国产| 久久青草综合色| 日韩熟女老妇一区二区性免费视频| 99久久综合免费| 久久久久久久久久久免费av| 亚洲国产色片| av福利片在线| 亚洲丝袜综合中文字幕| 久久久久久久精品精品| 久久久国产精品麻豆| 五月天丁香电影| 日本与韩国留学比较| 18禁动态无遮挡网站| 精品午夜福利在线看|