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

    基于可靠度的輸電塔抗風優(yōu)化研究

    2014-09-27 13:52:39鄭敏梁樞果熊鐵華
    湖南大學學報·自然科學版 2014年6期
    關(guān)鍵詞:失效模式優(yōu)化

    鄭敏+梁樞果+熊鐵華

    文章編號:16742974(2014)06003508

    收稿日期:20130716

    基金項目:國家自然科學基金資助項目(51078296)

    作者簡介:鄭 敏(1981-),女,浙江金華人,武漢大學講師,博士

    通訊聯(lián)系人,E-mail:liangsg@public.wh.hb.cn

    摘 要:提出一種基于可靠度約束的輸電塔實用抗風優(yōu)化設(shè)計方法.首先,識別出體系的失效模式并找出關(guān)鍵桿件同時提取各失效模式對應(yīng)的安全余量方程,計算各失效模式的可靠指標.然后,對可靠度不滿足的失效模式通過加強關(guān)鍵桿件的方法進行優(yōu)化.最后,對各失效模式進行綜合求出體系可靠指標,若體系可靠度滿足要求則優(yōu)化結(jié)束,若體系可靠度不滿足則通過優(yōu)化最小可靠指標的失效模式來提高體系可靠度直到滿足為止.編制程序?qū)λ苓M行了若干個風速的優(yōu)化,分析結(jié)果表明:1)在完整失效路徑上存在關(guān)鍵的破壞桿件,對關(guān)鍵桿件進行優(yōu)化可以有效地提高體系可靠度;2) 這是一種直接消除結(jié)構(gòu)薄弱點的優(yōu)化,優(yōu)化效果好; 3)隨優(yōu)化風速的增大,優(yōu)化的效率將變低;4)方法計算簡單,結(jié)果合理.

    關(guān)鍵詞:輸電塔;可靠指標;失效模式;優(yōu)化

    中圖分類號:TU312 文獻標識碼:A

    Study on the Reliabilitybased Wind Resistant 

    Optimization of Transmission Line Towers

    

    ZHENG Min, LIANG Shuguo, XIONG Tiehua

    (School of Civil and Architectural Engineering, Wuhan Univ, Wuhan, Hubei 430072,China),

    Abstract:A practical design method was proposed for the reliabilitybased wind resistant optimization of transmission line towers. Firstly, the failure modes and their key elements were identified, and their safety margin equations were established. Second, the failure modes with a reliability index of less than the target reliability index were optimized by strengthening its key element. Lastly, the system reliability index was calculated, and the failure mode with the least reliability index would be optimized if the system reliability index was less than the target reliability. An optimization example of a tower under several wind speeds was given and the results show: 1) there were key elements in the whole failure path, and the structure would be well strengthened if the key element was strengthened;2) it was a method in which the weak points were eliminated directly to strengthen the whole structure and it was very effective; 3) as the wind speed increased, it would be less effective; and 4) the method was simple and its results were reasonable.

    Key words: towers; reliability; failure modes; optimization

    

    塔架結(jié)構(gòu)廣泛應(yīng)用于廣播、電視、通信、輸電等領(lǐng)域,隨著建筑規(guī)范的更新,已建的很多塔架已不能滿足現(xiàn)行規(guī)范的要求,急需對其進行加固;在塔架設(shè)計中,為了更好地利用資源也存在優(yōu)化的問題.格構(gòu)式塔架的質(zhì)量較一般高層建筑結(jié)構(gòu)要輕得多,所以一般而言,風荷載對格構(gòu)式塔架的影響更大,是格構(gòu)式高聳結(jié)構(gòu)的主要側(cè)向荷載.由于塔架所受風荷載的復(fù)雜性以及塔線耦聯(lián)效應(yīng)的復(fù)雜性,以往塔架優(yōu)化中所用的風荷載都是一些經(jīng)驗公式[1-3],降低了優(yōu)化結(jié)果的可靠性.結(jié)構(gòu)設(shè)計中存在諸多不確定性,如材料強度不確定性、荷載不確定性等等,對其進行基于概率統(tǒng)計的可靠度分析更符合實際,是結(jié)構(gòu)設(shè)計理論的研究方向.目前,基于可靠度的結(jié)構(gòu)優(yōu)化取得了一定的成果,不過大多數(shù)只是在構(gòu)件的層次上[4-8],認為各個構(gòu)件同等重要, 這種處理對于靜定結(jié)構(gòu)當然沒有問題,但對于高次超靜定的塔架,把所有構(gòu)件對結(jié)構(gòu)的安全視為同等重要顯然是不合理的,將造成材料的浪費.基于體系可靠度的結(jié)構(gòu)優(yōu)化設(shè)計[9-11]理論比較全面、合理, 但由于體系可靠度計算的復(fù)雜性提高了問題的難度,對于實際結(jié)構(gòu),通常要進行大量的計算,缺乏實用性.

    針對目前輸電塔優(yōu)化存在的不足,本文主要從兩方面作了改進.一方面,本文采用的風荷載為由完全氣彈模型風洞試驗處理得到的等效靜力風荷載,更加符合實際.另一方面,由可靠度理論可知,體系可靠指標小于各種失效模式的可靠指標,且體系可靠指標對可靠指標小的失效模式更加敏感.本文提出一種計算量相對較少的可應(yīng)用于工程實際的基于體系可靠度優(yōu)化設(shè)計法.

    1 等效靜力風荷載

    1.1 風洞試驗

    本次風洞試驗?zāi)P蜑橐粋€輸電塔線體系完全氣動彈性模型,依據(jù)相似理論進行設(shè)計.本次模型設(shè)計除了滿足幾何相似以外,還滿足了氣動彈性模型最重要的相似參數(shù):Strouhal數(shù)相似、弗勞德數(shù)相似、彈性參數(shù)和慣性參數(shù)相似,氣彈模型各相似比見表1.本次風洞試驗?zāi)P偷脑褪呛舾邽?5 m跨度為750 m的酒杯塔線體系,塔構(gòu)件為不同截面尺寸的角鋼.

    考慮到該塔處于城市郊區(qū),試驗中利用塔尖和粗糙單元來模擬B類風場.本試驗在西南交通大學XNJD3風洞進行,風洞截面尺寸22.5 m×4.5 m,試驗段長度36 m,最高實驗風速可達17 m/s.本次試驗的采樣頻率為256 Hz,采樣時間為2 min.本試驗總共布置了8個激光位移測點.風洞試驗照片以及風向、測點的布置如圖1所示,測點所測位移方向垂直于紙面.

    表1 輸電塔的氣彈模型相似比(n=1∶40)

    Tab.1 The similar ratios

    相似比

    長度

    λL

    風速

    λU

    頻率

    λn

    拉伸剛度

    λEA

    質(zhì)量

    λm

    公式

    n

    n0.5

    n-0.5

    n3

    n3

    取值

    1∶40

    1∶400.5

    1∶40-0.5

    1∶403

    1∶403

    圖1 風洞試驗

    Fig.1 The wind tunnel test

    進行了單塔和塔線體系的多個風速下的風洞試驗.圖1(d)為90°風向角,也進行了0°,22.5°,45°,60°風向角的風洞試驗.由于該輸電塔線體系中線上的風荷載比塔直接所受風荷載大得多,而90°風向角時由于線的迎風面最大線所受的風荷載也最大,所以90°風向角時的風荷載是塔線體系的控制風荷載.下面的等效靜力風荷載是針對90°風荷載情況進行的.

    1.2 等效靜力風荷載

    有關(guān)等效風荷載的計算,詳見文獻[12],這里僅給出計算結(jié)果.這里采用的等效靜力風荷載分為平均風荷載和擬靜力脈動風荷載兩部分.

    1.2.1 平均風荷載

    1) 輸電塔直接所受的平均風荷載: 輸電塔直接所受的平均風荷載,參照《規(guī)范》[13],采用式(1)的形式.

    W1=w10?μz?μs1?Af1. (1)

    式中:W1為塔架風荷載平均值,N;w10為按10 m高的平均風速算得的風壓,w10=0.5?ρ?v210,N/m2;μz為風壓高度變化系數(shù),μz=z/10α ,α由地面粗糙度決定,這里按B類地貌,取α=0.15;μs1為體型系數(shù),為待擬合參數(shù),與風速無關(guān),z取塔架上各節(jié)點的標高,m;Af1為構(gòu)件承受風壓投影面積計算值,m2.

    根據(jù)試驗得到4個順風向測點位移平均值,按最小二乘法擬合得到μs1=2.24.

    2)線傳到塔上的平均風荷載:線傳到塔上的平均風荷載只考慮順風向,且平均地分配到導(dǎo)線的掛線處.線傳到塔上的平均風荷載參照《規(guī)范》[13],采用式(2)的形式.

    W2=w10?μz?μs2?Af2.(2) 

    式中:W2為線傳給塔架的風荷載平均值,N;μs2為體型系數(shù),為待擬合參數(shù);Af2為線承受風壓投影面積計算值,m2;μz為風壓高度變化系數(shù),μz=z/10α,α由地面粗糙度決定,這里按B類地貌,取α=0.15,z取塔架上掛線點的標高,m.

    由特定風速下輸電塔掛線與不掛線時同一測點平均位移差計算出該風速下的μs2,最后擬合得到:

    μs2=1.087 8+0.019 2?v.(3)

    式中:v為實際結(jié)構(gòu)10 m高的風速, m/s.

    1.2.2 擬靜力脈動風荷載

    擬靜力脈動風荷載采用模態(tài)風振力為基本標架,以廣義位移方差為標架值的形式.確定風振響應(yīng)均方根時先在模型上施加各模態(tài)風振力,再由振型分解法的思想按式(8)把各模態(tài)風振力下響應(yīng)匯總得到風振響應(yīng)均方根.

    1)模態(tài)選取: 為了得到風荷載動力響應(yīng)中起控制作用的模態(tài),在模型上施加按準定常得到的風壓譜.經(jīng)過計算,得到對于輸電塔部分而言,起控制作用的模態(tài)有兩階,振型如圖2所示,分別為順風向、橫風向以塔變形為主的第一階模態(tài).

    2)模態(tài)風振力: 單位廣義位移下模態(tài)風振力按式(4)計算.

    Fi=mi?ω2j?φij.(4)

    其中Fi為單位廣義位移下的模態(tài)風振力,量綱為[F]/[L];mi為節(jié)點i的質(zhì)量;φij為第j階模態(tài)振型在第i節(jié)點的值;n為塔的節(jié)點數(shù);ωj為第j階模態(tài)的自振頻率.

    3)廣義位移方差:根據(jù)廣義位移(協(xié))方差與測點位移(協(xié))方差的關(guān)系[式(5)],可由測點位移(協(xié))方差Ryij擬合廣義位移(協(xié))方差Rqmn.考慮到第一階主要是順風向變形,第二階主要是橫風向變形,可不考慮其相關(guān)性,兩階可分開擬合.

    Ryij=∑km=1∑kn=1φim?φjn?Rqmn. (5)

    式中:Ryij為測點位移(協(xié))方差,由風洞試驗數(shù)據(jù)處理得到;φim,φjn為第m和n階模態(tài)振型在i和j測點處的取值,由ANSYS模態(tài)提取;k為所選取的模態(tài)數(shù);Rqmn為第m和n階模態(tài)廣義位移(協(xié))方差,為待擬合的參數(shù).

    (a) 第1階

    (b)第2階

    圖2 所選模態(tài)

    Fig.2 The selected modes

    

    運用最小二乘法,分別由順風向、橫風向的測點位移(協(xié))方差擬合出第一階、第二階模態(tài)廣義位移方差Rq11和Rq22,最后由不同風速下的值擬合得到隨風速變化的公式(6)和(7),以便設(shè)計中使用.

    Rq11=3.844 3e-4?v2.74147, (6) 

    Rq22=7.540 7e-4?v2.04401.(7)

    式中: v為實際結(jié)構(gòu)10 m高的風速,m/s;Rq11和Rq22分別為第一階、第二階模態(tài)廣義位移方差, m2.

    4)風振響應(yīng)均方根: 有了廣義位移方差,按式(8)可以得到各種響應(yīng)的均方根.

    σi=φ2i1?Rq11+φ2i2?Rq2212. (8)

    式中:φij為第j模態(tài)響應(yīng)函數(shù)在i桿件的值,φij通過模態(tài)風振力施加到有限元模型上得到,可以是軸力、剪力、彎矩等內(nèi)力,也可以是軸向應(yīng)力、剪應(yīng)力等應(yīng)力,或者是其組合.

    2 失效模式識別

    2.1 構(gòu)件的破壞準則

    由于實際塔架的節(jié)點一般采用多個螺栓連接或焊接,這里構(gòu)件采用梁單元,考慮構(gòu)件的軸向力和彎矩.參照規(guī)范[13],構(gòu)件的破壞準則如下.

    1) 拉彎單元:

    Nm?An+MmM?W≤f. (9)

    式中:N,M為桿件的軸力和彎矩;m為構(gòu)件強度折減系數(shù),對于單肢連接的角鋼構(gòu)件(肢寬>40 mm)取0.70;An為構(gòu)件凈截面面積;mM為受彎構(gòu)件穩(wěn)定強度折減系數(shù),對應(yīng)角鋼,取1.0;W為桿件截面抗彎抵抗矩;f為鋼材屈服強度.

    2) 壓彎單元: 

    Nφ?mN?A+MmM?W≤f. (10)

    式中:φ為桿件軸心受壓穩(wěn)定系數(shù);mN為壓桿構(gòu)件穩(wěn)定系數(shù),角鋼根據(jù)翼緣板自由外伸寬度與厚度比值確定,其他符號同式(9).各符號的含義及取值詳見文獻[13].式(9)和(10)統(tǒng)一寫為式(11).

    σ≤f.(11)

    失效模式識別中,采用的是荷載及強度的平均值,荷載同時考慮自重和風荷載,自重視作確定值.式(11)可寫成式(12)的形式.

    σ平均風±σ自重≤fy. (12)

    式中:σ自重和σ平均風分別為由自重和平均風引起的應(yīng)力,按式(9)或(10)計算,其中內(nèi)力N和M取平均值;fy為鋼材的屈服強度平均值.

    某單元的承載力因子可取為

    re=fyσ自重σ平均風.(13)

    可看出,在均值意義上,re>1時桿件安全,re<1 時桿件破壞,re=1為臨界狀態(tài).所以在尋找失效模式時,以此為標準來判斷單元的失效情況.

    2.2 失效模式識別

    輸電塔屬于超靜定結(jié)構(gòu),某一個單元的失效并不意味著結(jié)構(gòu)的失效,只有失效的單元達到一定的數(shù)量,形成某一完整的失效模式,結(jié)構(gòu)成為機構(gòu),結(jié)構(gòu)才算失效.失效模式的識別參照董聰?shù)娜峙R界強度分枝約界法[14]的基本思想進行.

    記ck為階段約界系數(shù),cs為體系約界系數(shù),rmin為各個桿件承載力因子中最小值,rs為最小的體系承載力因子.X為失效模式集合,Xwz為完整的失效模式集合.以下失效模式的識別都是在均值下進行的,為了進行可靠度的計算,在每一步中都計算出順風向、橫風向第一模態(tài)風振力下的應(yīng)力φij1和φij2,以及自重應(yīng)力σ自重,平均風應(yīng)力σ平均風,并記錄到失效路徑X和Xwz中,其中i為破壞階段號,j為失效模式號.失效模式識別的流程如下:

    ①置完整失效路徑集合Xwz=φ,置體系承載力因子Rs=

    SymboleB@

     ,計算初始狀態(tài)各個桿件的承載力因子并找出其中最小值rmin,把承載力因子re

    ②在失效路徑集合X中,檢查是否形成機構(gòu).

    ③如果形成機構(gòu),把該失效路徑加入到完整失效路徑集合Xwz中.找出該路徑上各個桿件各自破壞時的承載力因子的最大值rmax ,rmax 就是這種失效模式的承載力因子R,而rmax 對應(yīng)的桿件即為該失效模式的關(guān)鍵桿件,判斷該桿件破壞時的階段即為該失效模式破壞的關(guān)鍵階段.如果該R值比體系承載力因子Rs小,更新體系承載力Rs為R,同時把該路徑從失效路徑X中刪除.

    ④如果未形成機構(gòu),把失效路徑上的單元從結(jié)構(gòu)中刪除,對于受拉破壞的桿件,刪除桿件的同時,在桿件兩端的節(jié)點處加上沿桿件長度的力以模擬受拉破壞后的屈服后拉力.計算剩余各個桿件的承載力因子re,找出其中的最小值rmin,把re<min ck?rmin ,cs?Rs的桿件加到失效路徑上并加入到失效路徑X中,同時把原來的失效路徑從失效路徑集合中刪除.

    ⑤若X=φ則搜索結(jié)束,若X≠φ轉(zhuǎn)入②.

    找到所有的完整失效路徑Xwz就是各種失效模式,且Xwz中還包含失效路徑上破壞單元對應(yīng)的σ自重, σ平均風, φij1, φij2等信息,為可靠度計算做好準備. 

    3 可靠度計算

    可靠度是指可靠的程度,用可靠指標β表示,它與失效概率pf的關(guān)系如式(14)所示.

    pf=φ-β.(14)

    后面的計算假設(shè)隨機變量均為正態(tài)分布,可靠度均用可靠指標β來表示.

    3.1 破壞階段安全余量方程

    3.1.1 安全余量方程提取

    各失效模式各破壞階段的安全余量方程如下:

    mij=fy-σ自重-σ平均風-φij1?q1-φij2?q2. (15)

    式中:j為失效模式號;i為失效路徑上破壞階段號;fy為鋼材屈服強度,為一隨機變量;σ自重為自重應(yīng)力,為一確定變量,由2.2部分得到;σ平均風為平均風應(yīng)力,為一確定變量,由2.2部分得到;φij1, φij2分別

    為順風向、橫風向第一模態(tài)應(yīng)力振型在破壞單元的取值,由2.2部分得到;q1,q2分別為順風向、橫風向第一模態(tài)廣義位移,為隨機變量.

    關(guān)于fy,以Q345為例,Q345表示屈服強度標準值fybz為345 N/mm2,是具有95%保證率的屈服強度的較低值,則屈服強度平均值可按式(16)得到:

    fy=fybz+1.645σfy. (16)

    記屈服強度變異系數(shù)為δ,用σfy=δ?fy 代入 式(16)整理可得:

    fy=fybz1-1.645δ. (17)

    q1,q2為順風向、橫風向第一模態(tài)廣義位移,均值為0,由文獻得到了順風向、橫風向第一模態(tài)廣義位移的均方差Rq11和Rq22,且不考慮兩模態(tài)的相關(guān)性,故q1和q2為相互獨立的隨機變量,均方根σq1和σq2可由式(18)和(19)確定.

    σq1=3.844 3e-4?v2.741 47,(18)

    σq2=7.540 7e-4?v2.044 01.(19)

    3.1.2 安全余量方程標準化

    式(15)中用U1,U2,U3代替fy,q1,q2,U1=fy-fy/σfy,U2=q1/σq1,U3=q2/σq2,U1,U2,U3為相互獨立的標準正態(tài)分布隨機變量,對系數(shù)做歸一化處理,即系數(shù)滿足式(21).得到各失效模式對應(yīng)的各階段的標準形式的安全余量方程(20).式(20)中右邊的常數(shù)項βij就是第j失效模態(tài)第i破壞階段對應(yīng)的可靠指標. 

    Mij=aij1?U1-aij2?U2-aij3?U3+βij, (20)

    ∑3k=1a2ijk=1. (21)

    3.2 單個失效模態(tài)可靠指標

    某失效模式的各個階段屬于并聯(lián),只有所有階段都破壞了才形成失效模式,即結(jié)構(gòu)按某失效模式失效的條件是該失效模式各個階段的安全余量方程均小于0,發(fā)生第j失效模式應(yīng)滿足式(22).

    a1j1?U1-a1j2?U2-a1j3?U3+β1j<0,

    a2j1?U1-a2j2?U2-a2j3?U3+β2j<0,

    

    amj1?U1-amj2?U2-amj3?U3+βmj<0.(22) 

    式中:m為失效階段數(shù).可以用逐步等效線性化Johnson求交法來得到第j失效模式的等效安全余量方程.關(guān)于該求解方法詳見文獻[15].計算得到的第j失效模式等效安全余量方程如式(23),式中系數(shù)滿足式(24).

    Mj=aj1?U1-aj2?U2-aj3?U3+βj,(23) 

    ∑3k=1a2jk=1.(24)

    式中:βj為第j失效模式的可靠度指標.

    3.3 體系可靠指標

    體系的各種失效模式屬于串聯(lián),只要有一種失效模式形成,結(jié)構(gòu)即告破壞.不過如果從結(jié)構(gòu)安全角度來看,各種失效模式也是并聯(lián),只有所有的失效模式都安全,結(jié)構(gòu)才安全,即結(jié)構(gòu)安全應(yīng)滿足式(25).

    a11?U1-a12?U2-a13?U3+β1>0,

    a21?U1-a22?U2-a23?U3+β2>0,

    am1?U1-am2?U2-am3?U3+βm>0.25 

    為了能用逐步等效線性化Johnson求交法來計算等效安全余量方程,式兩邊同乘以(-1),得式(26).

    -a11?U1-a12?U2-a13?U3+β1<0,

    -a21?U1-a22?U2-a23?U3+β2<0,

    -am1?U1-am2?U2-am3?U3+βm<0. (26)

    利用逐步等效線性化Johnson求交法來計算得到式(26)的等效的不等式(27).

    -a1?U1-a2?U2-a3?U3+β<0.(27)

    式(27)為結(jié)構(gòu)安全應(yīng)滿足的不等式,結(jié)構(gòu)失效事件與結(jié)構(gòu)安全事件互補,故結(jié)構(gòu)失效的等效不等式如式(28).

    a1?U1-a2?U2-a3?U3+β<0.(28)

    故結(jié)構(gòu)體系的安全余量方程如式(29)所示.

    M=a1?U1-a2?U2-a3?U3+β. (29)

    式中:β即為體系可靠度指標.

    4 塔架優(yōu)化

    4.1 優(yōu)化模型

    這里做體系可靠指標滿足要求的結(jié)構(gòu)質(zhì)量增量最小的優(yōu)化,即

    find Aii=1,2,…,m, 

    min ΔWA=∑i=mi=1CiAi-Ai0,

    st. β≥βu,

    AiεB1,B2,B3,…,BB.

    其中:Ai為需要優(yōu)化的失效模式對應(yīng)的關(guān)鍵桿件的截面; m為需要優(yōu)化的關(guān)鍵桿件個數(shù); Ai必須在B1,B2,B3,…,BB中選取,B為可選取的截面類型數(shù)目;ΔWA為結(jié)構(gòu)的總質(zhì)量的增量;Ci 為與材料、桿件長度有關(guān)的系數(shù);Ai0為Ai優(yōu)化前的截面; β為體系可靠指標,按第3部分計算;βu為目標可靠指標,由結(jié)構(gòu)的安全等級以及破壞的性質(zhì)確定,具體取值可參見規(guī)范[16].

    4.2 優(yōu)化流程

    根據(jù)體系可靠指標的兩個特點,采用準則法進行優(yōu)化.第一,由于體系可靠指標一定小于任一失效模式的可靠指標,提出首先把各個失效模式優(yōu)化到滿足目標可靠指標的要求.第二,因為體系可靠指標對可靠指標小的失效模式更敏感,提出在各個失效模式單個可靠指標滿足要求的前提下若體系可靠指標不滿足則提高最小可靠指標的那個失效模式的可靠指標來優(yōu)化體系,直到體系可靠指標滿足要求為止.其優(yōu)化流程如下:

    ①尋找失效模式并提取各種失效模式失效路徑上破壞單元的σ自重,σ平均風,φij1,φij2,按2.2部分進行計算.

    ②計算各失效模式各自的可靠指標βj,按3.1和3.2部分進行.若各失效模式βj均大于目標可靠指標,轉(zhuǎn)入④步.若存在βj≤βu 的失效模式,轉(zhuǎn)入③步.

    ③對于βj≤βu的失效模式,優(yōu)化關(guān)鍵桿件截面使該失效模式可靠指標βj大于目標可靠指標.

    ④計算體系可靠指標β,按3.3部分進行.若β≥βu轉(zhuǎn)入⑥步,若β<βu,轉(zhuǎn)入第⑤步.

    ⑤找出最小可靠指標βmin的那種失效模式,把對應(yīng)的關(guān)鍵桿件截面加大一號來提高體系可靠度β,轉(zhuǎn)入第④步.

    ⑥優(yōu)化結(jié)束.

    5 算例

    一呼高為75 m的塔架,圖1所示為其風洞試驗的照片,風場為B類風場,對其在10 m高風速分別為19 m/s,20 m/s,21 m/s,22 m/s,23 m/s,24 m/s的情況下進行可靠度下的質(zhì)量最小化的優(yōu)化,桿件截面在表2中選取.目標可靠指標取3.7.

    表2 截面型號

    Tab.2 Sections

    序號

    1

    2

    3

    4

    5

    截面

    L40X4

    L45X4

    L50X4

    L45X5

    L56X4

    序號

    6

    7

    8

    9

    10

    截面

    L50X5

    L56X5

    L63X5

    L70X5

    L75X5

    序號

    11

    12

    13

    14

    15

    截面

    L70X6

    L75X6

    L80X6

    L80X7

    L90X7

    序號

    16

    17

    18

    19

    20

    截面

    L100X7

    L110X8

    L125X8

    L110X10

    L125X10

    序號

    21

    22

    23

    24

    25

    截面

    L140X10

    L140X12

    L160X12

    L180X14

    L180X16

    序號

    26

    27

    28

    29

    30

    截面

    L200X16

    L200X18

    L200X20

    L200X24

    2L200X18

    在風速24 m/s范圍內(nèi)進行各風速下基于可靠度的優(yōu)化,通過數(shù)值計算,找到了9種典型的失效模式,如圖3所示.

    圖3 前9種失效模式

    Fig.3 The first ninefailure modes

    

    各失效模式破壞風速、破壞起因詳見表3,表中破壞風速對應(yīng)實際結(jié)構(gòu)10 m高的風速.一種典型的失效模式包含多種起因的失效模式,表3中列出

    表3 各失效模式情況一覽表

    Tab.3 The failure modes

    模式

    破壞風速

    /(m?s-1)

    破壞起因

    失效路徑

    1

    19

    塔身一級腹桿失穩(wěn)破壞

    770~1 895

    2

    21

    塔身一級腹桿失穩(wěn)破壞

    1 859~1 904

    3

    21

    塔身一級腹桿失穩(wěn)破壞

    1 039~509

    4

    23

    塔身立柱受拉強度破壞

    1 947

    5

    23

    塔身立柱受拉強度破壞

    1 960

    6

    23

    塔身立柱受拉強度破壞

    1 954

    7

    24

    塔腳立柱受壓強度破壞

    832

    8

    24

    上屈臂立柱受拉強度破壞

    386

    9

    24

    塔身立柱受拉強度破壞

    1 940

    的為破壞風速較低的一種情況,而且在失效路徑中只給出了主要的失效單元.失效路徑所對應(yīng)的桿件號如圖4所示.

    (a) 塔頭桿件號

    (b) 塔身塔腳桿件號

    圖4桿件號示意圖

    Fig.4 Element numbers

    

    對于不同的風速要求,對可靠度不滿足要求的失效模式進行優(yōu)化,優(yōu)化通過加大關(guān)鍵桿件的截面來實現(xiàn).由于塔架是一個高次超靜定結(jié)構(gòu),優(yōu)化過程中會產(chǎn)生內(nèi)力重分布,消除某種失效模式往往有多種途徑,而且?guī)追N需要優(yōu)化的失效模式的關(guān)鍵桿件之間往往互相影響,優(yōu)化時必須綜合考慮,在保證承載力的條件下找出使結(jié)構(gòu)質(zhì)量增加最小的一種優(yōu)化方案.

    表4為各種風速的優(yōu)化結(jié)果,限于篇幅, 24 m/s的計算結(jié)果未列出.表5列出優(yōu)化的效果,圖5為質(zhì)量增加百分比隨優(yōu)化風速的變化圖.可以看出,優(yōu)化的效果較好,在質(zhì)量增加不多的情況下,抗風承載力得到了很大提高.也可以看出,隨優(yōu)化風速的增大,質(zhì)量增加百分比斜率總的趨勢是增加的,也就是意味著隨著優(yōu)化風速的增大優(yōu)化的效果是在變差的,這主要是因為優(yōu)化風速增加,結(jié)構(gòu)趨于等強度,同樣的風速增加時產(chǎn)生的失效模式變多,而且同一種失效模式需要優(yōu)化的設(shè)計變量變多.由表5可以看出,除了因為桿件截面不連續(xù)引起的可靠指標稍大一點外(優(yōu)化風速為19 m/s),優(yōu)化后體系可靠指標均略大于目標可靠指標3.7,優(yōu)化比較合理.

    表4 優(yōu)化結(jié)果表

    Tab.4 The optimized results

    風速

    / (m?s-1)

    優(yōu)化

    變量

    所處位置

    截面序號

    優(yōu)化前

    優(yōu)化后

    失效

    模式

    19

    63

    塔身腹桿

    18

    20

    20

    63

    塔身腹桿

    18

    20

    21

    63

    塔身腹桿

    18

    20

    67

    塔身腹桿

    16

    17

    140

    下屈臂腹桿

    9

    10

    138

    下屈臂腹桿

    8

    9

    22

    63

    塔身腹桿

    18

    21

    67

    塔身腹桿

    16

    17

    125

    下屈臂立桿

    20

    21

    140

    下屈臂腹桿

    9

    12

    138

    下屈臂腹桿

    8

    9

    23

    63

    塔身腹桿

    18

    21

    67

    塔身腹桿

    16

    17

    56

    塔身立柱

    23

    24

    140

    下屈臂腹桿

    9

    12

    121

    塔頸立柱

    20

    21

    125

    下屈臂立桿

    20

    23

    123

    下屈臂立桿

    20

    22

    148

    塔頸腹桿

    15

    16

    138

    下屈臂腹桿

    8

    9

    53

    塔身立柱

    24

    25

    64

    塔身腹桿

    17

    18

    57

    塔身立柱

    23

    24

    54

    塔身立柱

    24

    25

    六

    表5 優(yōu)化效果

    Tab.5 The optimization effect

    序號

    優(yōu)化風速/

    (m?s-1)

    可靠指標

    優(yōu)化前

    優(yōu)化后

    承載力

    增量/%

    質(zhì)量

    增量/%

    1

    18

    3.9

    3.9

    -

    -

    2

    19

    3.1

    4.4

    11

    0.8

    3

    20

    2.5

    4.0

    23

    0.8

    4

    21

    1.9

    3.9

    36

    1.3

    5

    22

    1.3

    3.7

    49

    1.9

    6

    23

    0.4

    3.8

    63

    3.8

    7

    24

    -0.4

    3.7

    78

    7.0 

    破壞風速/(m?s-1)

    圖5 優(yōu)化效果曲線

    Fig.5The optimization effect curve

    6 結(jié) 論

    通過本文分析,得到以下結(jié)論:

    1)在完整失效路徑上存在關(guān)鍵的破壞桿件,對關(guān)鍵桿件進行優(yōu)化可以有效地提高體系可靠度.

    2)這是一種直接消除結(jié)構(gòu)薄弱點的優(yōu)化,優(yōu)化效果好.

    3)隨優(yōu)化風速的增大,優(yōu)化的效率將變低.這主要是因為隨著優(yōu)化的進行,結(jié)構(gòu)趨于等強度,同樣的風速增加時需要優(yōu)化的失效模式變多,而且同一種失效模式需要優(yōu)化的設(shè)計變量變多.

    4)本文提出的基于可靠度的優(yōu)化方法計算簡單,結(jié)果合理.

    參考文獻

    [1] 田琪凌,伋雨林,陳振,等.500 kV輸電塔承載能力分析及優(yōu)化設(shè)計[J]. 華中科技大學學報:城市科學版,2010,27(1):64-69.

    TIAN Qiling , JI Yulin , CHEN Zhen, et al. Analysis on the bearing capacity of 500 kV overhead iced transmission tower and structural optimization[J]. Journal of HUST: Urban Science, 2010,27(1):64-69. (In Chinese)

    [2] 熊鐵華,梁樞果,鄒良浩.考慮斷線時輸電鐵塔的失效模式及其極限荷載[J]. 土木工程學報,2009,42(11):86-90.

    XIONG Tiehua, LIANG Shuguo , ZOU Lianghao. Dominant failure modes of transmission towers and the limit loads with wire break off[J]. China Civil Engineering Journal, 2009, 42(11): 86-90. (In Chinese)

    [3] 熊鐵華,梁樞果,吳海洋.某輸電線路鐵塔覆冰條件下的失效模式分析[J].計算力學學報,2011, 28(3): 468-478.

    XIONG Tiehua, LIANG Shuguo, WU Haiyang. Failure modes analysis of a broken down transmission tower under ice loads[J]. Chinese Journal of Computational Mechanics, 2011, 28(3): 468-478. (In Chinese)

    [4] NOILUBLAO N, BUREERAT S. Simultaneous topology, shape and sizing optimization of a threedimensional slender truss tower using multiobjective evolutionary algorithms[J]. Computers and Structures,2011,89(23/24):2531-2538.

    [5] RAO G W. Optimum designs for transmission line towers[J]. Computer and Structreres, 1995, 57(1):81-92. 

    [6] 方從嚴,劉寧.基于可靠度的非線性隨機優(yōu)化[J].水利水電科技進展,2005,25(5):31-32.

    FANG Congyan, LIU Ning. Reliabilitybased structural nonlinear stochastic optimization[J]. Advances in Science and Technology of Water Resources,2005,25(5):31-32.(In Chinese)

    [7] 王光遠,吳波.基于可靠度最優(yōu)分布的桁架優(yōu)化設(shè)計[J]. 哈爾濱建筑工程學院學報,1990,23(2):10-23.

    WANG Guangyuan, WU Bo. Design of trusses based on optimal reliability distribution[J]. Journal of Harbin Architecture & Civil Engineering Institute, 1990, 23(2):10-23. (In Chinese)

    [8] 黃炳生,王文濤.工程系統(tǒng)可靠度優(yōu)化的多階段決策算子法[J]. 工程力學,1997,14(1):69-78.

    HUANG Bingsheng, WANG Wentao. Multistage decision operator technique for reliability optimization of engineering system[J]. Engineering Mechanics,1997,14(1):69-78.(In Chinese)

    [9] NATARARAJAN K, SANTHAKUMAR A R. Reliabilitybased optimization of transmission line towers[J]. Computers and Structures, 1995,55(3):387-403.

    [10]方從嚴, 劉寧,卓家壽.結(jié)構(gòu)可靠度優(yōu)化的敏感性分析[J].應(yīng)用力學學報,2005,22(1):63-66.

    FANG Congyan, LIU Ning, ZHUO Jiashou. Structural reliability optimization with sensitivity analysis[J]. Chinese Journal of Applied Mechanics,2005,22(1):63-66. (In Chinese)

    [11]許林. 基于結(jié)構(gòu)可靠度的優(yōu)化設(shè)計研究[D].大連:大連理工大學,2004:27-114.

    XU Lin. A study on reliability based structural optimization[D]. Dalian: Dalian University of Technology,2004:27-114. (In Chinese)

    [12]鄭敏,梁樞果,熊鐵華.基于塔線體系氣彈模型風洞試驗的輸電塔風致響應(yīng)計算[J].武漢理工大學學報,2013,35(11): 124-131.

    ZHENG Min, LIANG Shuguo, XIONG Tiehua. Windinduced response calculation of the tower in transmission towerline system based on wind tunnel tests[J]. Journal of Wuhan University of Technology, 2013, 35(11): 124-131. (In Chinese)

    [13]DL/T 5154—2002 架空送電線路桿塔結(jié)構(gòu)設(shè)計技術(shù)規(guī)定[S]. 北京: 中國電力出版社, 2002:9-30.

    DL/T 5154—2002 Technical regulation of design for tower and pole structures of overhead transmission line[S]. Beijing: China Electric Power Press,2002:9-30. (In Chinese)

    [14]董聰.現(xiàn)代結(jié)構(gòu)系統(tǒng)可靠性理論及其應(yīng)用[M].北京:科學出版社,2001:62-225.

    DONG Cong.Reliability theory and its application in the modern structural systems[M].Beijing: Science Press,2001:62-225. (In Chinese)

    [15]劉寧, 李同春. 用三維隨機有限元尋求結(jié)構(gòu)的最大可能失效模式[J]. 水利學報, 1996(3): 36-43.

    LIU Ning, LI Tongchun. A method for searching the mostprobable failure mode by on 3D SFEM[J]. Journal of Hydraulic Engineering, 1996(3): 36-43. (In Chinese)

    [16]GB 50068—2001建筑結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標準[S].北京:中國建筑工業(yè)出版社,2001:1-19.

    GB 50068—2001 Unified standard for reliability design of building structures[S]. Beijing: China Architecture & Building Press, 2001:1-19. (In Chinese)

    WANG Guangyuan, WU Bo. Design of trusses based on optimal reliability distribution[J]. Journal of Harbin Architecture & Civil Engineering Institute, 1990, 23(2):10-23. (In Chinese)

    [8] 黃炳生,王文濤.工程系統(tǒng)可靠度優(yōu)化的多階段決策算子法[J]. 工程力學,1997,14(1):69-78.

    HUANG Bingsheng, WANG Wentao. Multistage decision operator technique for reliability optimization of engineering system[J]. Engineering Mechanics,1997,14(1):69-78.(In Chinese)

    [9] NATARARAJAN K, SANTHAKUMAR A R. Reliabilitybased optimization of transmission line towers[J]. Computers and Structures, 1995,55(3):387-403.

    [10]方從嚴, 劉寧,卓家壽.結(jié)構(gòu)可靠度優(yōu)化的敏感性分析[J].應(yīng)用力學學報,2005,22(1):63-66.

    FANG Congyan, LIU Ning, ZHUO Jiashou. Structural reliability optimization with sensitivity analysis[J]. Chinese Journal of Applied Mechanics,2005,22(1):63-66. (In Chinese)

    [11]許林. 基于結(jié)構(gòu)可靠度的優(yōu)化設(shè)計研究[D].大連:大連理工大學,2004:27-114.

    XU Lin. A study on reliability based structural optimization[D]. Dalian: Dalian University of Technology,2004:27-114. (In Chinese)

    [12]鄭敏,梁樞果,熊鐵華.基于塔線體系氣彈模型風洞試驗的輸電塔風致響應(yīng)計算[J].武漢理工大學學報,2013,35(11): 124-131.

    ZHENG Min, LIANG Shuguo, XIONG Tiehua. Windinduced response calculation of the tower in transmission towerline system based on wind tunnel tests[J]. Journal of Wuhan University of Technology, 2013, 35(11): 124-131. (In Chinese)

    [13]DL/T 5154—2002 架空送電線路桿塔結(jié)構(gòu)設(shè)計技術(shù)規(guī)定[S]. 北京: 中國電力出版社, 2002:9-30.

    DL/T 5154—2002 Technical regulation of design for tower and pole structures of overhead transmission line[S]. Beijing: China Electric Power Press,2002:9-30. (In Chinese)

    [14]董聰.現(xiàn)代結(jié)構(gòu)系統(tǒng)可靠性理論及其應(yīng)用[M].北京:科學出版社,2001:62-225.

    DONG Cong.Reliability theory and its application in the modern structural systems[M].Beijing: Science Press,2001:62-225. (In Chinese)

    [15]劉寧, 李同春. 用三維隨機有限元尋求結(jié)構(gòu)的最大可能失效模式[J]. 水利學報, 1996(3): 36-43.

    LIU Ning, LI Tongchun. A method for searching the mostprobable failure mode by on 3D SFEM[J]. Journal of Hydraulic Engineering, 1996(3): 36-43. (In Chinese)

    [16]GB 50068—2001建筑結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標準[S].北京:中國建筑工業(yè)出版社,2001:1-19.

    GB 50068—2001 Unified standard for reliability design of building structures[S]. Beijing: China Architecture & Building Press, 2001:1-19. (In Chinese)

    WANG Guangyuan, WU Bo. Design of trusses based on optimal reliability distribution[J]. Journal of Harbin Architecture & Civil Engineering Institute, 1990, 23(2):10-23. (In Chinese)

    [8] 黃炳生,王文濤.工程系統(tǒng)可靠度優(yōu)化的多階段決策算子法[J]. 工程力學,1997,14(1):69-78.

    HUANG Bingsheng, WANG Wentao. Multistage decision operator technique for reliability optimization of engineering system[J]. Engineering Mechanics,1997,14(1):69-78.(In Chinese)

    [9] NATARARAJAN K, SANTHAKUMAR A R. Reliabilitybased optimization of transmission line towers[J]. Computers and Structures, 1995,55(3):387-403.

    [10]方從嚴, 劉寧,卓家壽.結(jié)構(gòu)可靠度優(yōu)化的敏感性分析[J].應(yīng)用力學學報,2005,22(1):63-66.

    FANG Congyan, LIU Ning, ZHUO Jiashou. Structural reliability optimization with sensitivity analysis[J]. Chinese Journal of Applied Mechanics,2005,22(1):63-66. (In Chinese)

    [11]許林. 基于結(jié)構(gòu)可靠度的優(yōu)化設(shè)計研究[D].大連:大連理工大學,2004:27-114.

    XU Lin. A study on reliability based structural optimization[D]. Dalian: Dalian University of Technology,2004:27-114. (In Chinese)

    [12]鄭敏,梁樞果,熊鐵華.基于塔線體系氣彈模型風洞試驗的輸電塔風致響應(yīng)計算[J].武漢理工大學學報,2013,35(11): 124-131.

    ZHENG Min, LIANG Shuguo, XIONG Tiehua. Windinduced response calculation of the tower in transmission towerline system based on wind tunnel tests[J]. Journal of Wuhan University of Technology, 2013, 35(11): 124-131. (In Chinese)

    [13]DL/T 5154—2002 架空送電線路桿塔結(jié)構(gòu)設(shè)計技術(shù)規(guī)定[S]. 北京: 中國電力出版社, 2002:9-30.

    DL/T 5154—2002 Technical regulation of design for tower and pole structures of overhead transmission line[S]. Beijing: China Electric Power Press,2002:9-30. (In Chinese)

    [14]董聰.現(xiàn)代結(jié)構(gòu)系統(tǒng)可靠性理論及其應(yīng)用[M].北京:科學出版社,2001:62-225.

    DONG Cong.Reliability theory and its application in the modern structural systems[M].Beijing: Science Press,2001:62-225. (In Chinese)

    [15]劉寧, 李同春. 用三維隨機有限元尋求結(jié)構(gòu)的最大可能失效模式[J]. 水利學報, 1996(3): 36-43.

    LIU Ning, LI Tongchun. A method for searching the mostprobable failure mode by on 3D SFEM[J]. Journal of Hydraulic Engineering, 1996(3): 36-43. (In Chinese)

    [16]GB 50068—2001建筑結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標準[S].北京:中國建筑工業(yè)出版社,2001:1-19.

    GB 50068—2001 Unified standard for reliability design of building structures[S]. Beijing: China Architecture & Building Press, 2001:1-19. (In Chinese)

    猜你喜歡
    失效模式優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
    長輸管道失效模式診斷技術(shù)綜述
    電動汽車動力電池系統(tǒng)故障診斷研究
    覆冰荷載作用下高壓輸電塔可靠性研究
    科技視界(2016年24期)2016-10-11 18:43:38
    醫(yī)療失效模式與效應(yīng)分析
    A公司新產(chǎn)品項目開發(fā)過程中的技術(shù)風險管理研究
    黄色毛片三级朝国网站| 狠狠狠狠99中文字幕| 不卡一级毛片| 一本大道久久a久久精品| 午夜日韩欧美国产| 亚洲伊人久久精品综合| 操美女的视频在线观看| 国产欧美日韩综合在线一区二区| 成人三级做爰电影| 国产黄频视频在线观看| 亚洲av国产av综合av卡| 黄频高清免费视频| 亚洲国产日韩一区二区| 亚洲男人天堂网一区| 国产精品二区激情视频| 亚洲欧美激情在线| 黑人猛操日本美女一级片| 亚洲人成电影免费在线| 99香蕉大伊视频| 午夜福利乱码中文字幕| 新久久久久国产一级毛片| 久久国产精品男人的天堂亚洲| 曰老女人黄片| 中文字幕精品免费在线观看视频| 国产1区2区3区精品| 亚洲精品国产色婷婷电影| 久久久精品国产亚洲av高清涩受| 国产国语露脸激情在线看| 欧美日韩福利视频一区二区| 汤姆久久久久久久影院中文字幕| 人成视频在线观看免费观看| 亚洲,欧美精品.| 亚洲国产看品久久| 丁香六月欧美| 黑人欧美特级aaaaaa片| 多毛熟女@视频| 国产区一区二久久| www.熟女人妻精品国产| 欧美日韩亚洲国产一区二区在线观看 | 精品熟女少妇八av免费久了| 久久精品成人免费网站| 麻豆av在线久日| 欧美在线一区亚洲| 亚洲色图 男人天堂 中文字幕| 日韩有码中文字幕| 三级毛片av免费| 首页视频小说图片口味搜索| 国产欧美日韩一区二区精品| 夜夜夜夜夜久久久久| 国产高清国产精品国产三级| 国产激情久久老熟女| 久久国产精品人妻蜜桃| 亚洲av电影在线观看一区二区三区| 国产成人精品久久二区二区免费| 老司机在亚洲福利影院| 欧美日韩国产mv在线观看视频| 热re99久久国产66热| 99久久综合免费| 91精品伊人久久大香线蕉| 一本久久精品| 欧美 亚洲 国产 日韩一| 一级黄色大片毛片| 亚洲欧美精品综合一区二区三区| 美女国产高潮福利片在线看| 亚洲三区欧美一区| 日韩制服丝袜自拍偷拍| 热99久久久久精品小说推荐| 在线精品无人区一区二区三| 欧美午夜高清在线| 18禁国产床啪视频网站| 欧美xxⅹ黑人| 91国产中文字幕| 成人黄色视频免费在线看| 热99re8久久精品国产| av网站免费在线观看视频| kizo精华| 亚洲,欧美精品.| 99国产综合亚洲精品| www.999成人在线观看| 成人免费观看视频高清| 免费女性裸体啪啪无遮挡网站| 久久性视频一级片| av又黄又爽大尺度在线免费看| 欧美精品av麻豆av| 精品一区二区三区四区五区乱码| 国产精品秋霞免费鲁丝片| 成人国语在线视频| 777米奇影视久久| 中文字幕av电影在线播放| 五月天丁香电影| 纯流量卡能插随身wifi吗| 亚洲精品一区蜜桃| 纯流量卡能插随身wifi吗| 久久久国产欧美日韩av| h视频一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产91精品成人一区二区三区 | 蜜桃在线观看..| 日韩欧美一区二区三区在线观看 | 欧美乱码精品一区二区三区| 亚洲欧美激情在线| 人人妻,人人澡人人爽秒播| 久久天堂一区二区三区四区| 高清视频免费观看一区二区| 欧美xxⅹ黑人| 色94色欧美一区二区| 欧美日韩视频精品一区| 国产成人免费观看mmmm| 久久人人爽av亚洲精品天堂| 国产精品久久久久成人av| 亚洲三区欧美一区| 久久久久久久大尺度免费视频| 国产在线观看jvid| 国产一区二区三区综合在线观看| 一个人免费在线观看的高清视频 | 久久久久久久精品精品| 一级a爱视频在线免费观看| 韩国高清视频一区二区三区| 久久狼人影院| 中文字幕人妻丝袜一区二区| 黄色视频不卡| 一级片免费观看大全| 欧美日韩视频精品一区| 久久久久久久久久久久大奶| 一级,二级,三级黄色视频| 在线观看免费日韩欧美大片| 中国美女看黄片| 亚洲色图 男人天堂 中文字幕| 亚洲精华国产精华精| 黄色片一级片一级黄色片| 99香蕉大伊视频| 亚洲欧洲精品一区二区精品久久久| 精品少妇内射三级| 少妇粗大呻吟视频| 日韩 亚洲 欧美在线| 丁香六月天网| 岛国毛片在线播放| 国产精品 国内视频| 免费观看av网站的网址| 中文字幕人妻丝袜制服| 老汉色av国产亚洲站长工具| 51午夜福利影视在线观看| 丝袜人妻中文字幕| 日韩有码中文字幕| 一级片免费观看大全| 成人亚洲精品一区在线观看| 丝瓜视频免费看黄片| 三上悠亚av全集在线观看| 亚洲国产看品久久| 麻豆av在线久日| 伦理电影免费视频| 中国国产av一级| 国产精品国产三级国产专区5o| 精品亚洲乱码少妇综合久久| 亚洲成人国产一区在线观看| 99精国产麻豆久久婷婷| 亚洲天堂av无毛| 中文字幕最新亚洲高清| 国产亚洲精品一区二区www | 最黄视频免费看| 亚洲人成电影免费在线| 三级毛片av免费| 亚洲avbb在线观看| 日韩一卡2卡3卡4卡2021年| 一级片'在线观看视频| 精品少妇内射三级| 九色亚洲精品在线播放| 啦啦啦啦在线视频资源| 国产高清视频在线播放一区 | 国产一区二区三区在线臀色熟女 | 9色porny在线观看| 日韩制服丝袜自拍偷拍| 在线观看免费日韩欧美大片| 日韩免费高清中文字幕av| 国产精品1区2区在线观看. | 日韩制服丝袜自拍偷拍| 黄片大片在线免费观看| 亚洲精品国产av成人精品| 国产不卡av网站在线观看| 波多野结衣一区麻豆| 91九色精品人成在线观看| 日韩 欧美 亚洲 中文字幕| 一区二区三区精品91| 在线观看免费高清a一片| 十八禁高潮呻吟视频| 午夜福利,免费看| bbb黄色大片| 日本av手机在线免费观看| 麻豆乱淫一区二区| 日韩欧美免费精品| 精品国产乱子伦一区二区三区 | 丰满人妻熟妇乱又伦精品不卡| 亚洲专区国产一区二区| 最近中文字幕2019免费版| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 久久精品国产亚洲av高清一级| 欧美精品av麻豆av| 亚洲第一av免费看| 欧美日韩一级在线毛片| 久久国产精品影院| 黄色视频,在线免费观看| 欧美国产精品一级二级三级| 免费黄频网站在线观看国产| 欧美亚洲 丝袜 人妻 在线| 热re99久久国产66热| 欧美变态另类bdsm刘玥| 国产国语露脸激情在线看| 欧美日韩精品网址| 久久人妻熟女aⅴ| 成人18禁高潮啪啪吃奶动态图| 丝瓜视频免费看黄片| 男女边摸边吃奶| 久久国产亚洲av麻豆专区| 美女中出高潮动态图| 一区在线观看完整版| 免费在线观看视频国产中文字幕亚洲 | 国产精品成人在线| 黄片播放在线免费| 久久久久久久大尺度免费视频| 久久女婷五月综合色啪小说| 少妇的丰满在线观看| 狂野欧美激情性bbbbbb| 亚洲国产毛片av蜜桃av| 亚洲欧美精品自产自拍| www.精华液| 黄色 视频免费看| 国产av一区二区精品久久| 亚洲人成电影免费在线| 在线十欧美十亚洲十日本专区| 精品国产一区二区三区四区第35| 成人手机av| 国产成人精品久久二区二区免费| 涩涩av久久男人的天堂| 狠狠狠狠99中文字幕| 男女国产视频网站| 欧美日韩精品网址| 欧美变态另类bdsm刘玥| 欧美成狂野欧美在线观看| 9色porny在线观看| 妹子高潮喷水视频| 欧美黑人精品巨大| 欧美午夜高清在线| 亚洲一码二码三码区别大吗| 欧美人与性动交α欧美精品济南到| 欧美日韩一级在线毛片| 国产高清国产精品国产三级| 国产老妇伦熟女老妇高清| 97在线人人人人妻| 岛国毛片在线播放| 日韩制服骚丝袜av| 精品国产乱子伦一区二区三区 | 一级片免费观看大全| 制服诱惑二区| 欧美国产精品一级二级三级| 欧美人与性动交α欧美软件| 蜜桃国产av成人99| 久久中文看片网| 日韩中文字幕视频在线看片| 久久久久久久大尺度免费视频| 男女边摸边吃奶| 中文字幕精品免费在线观看视频| 黑丝袜美女国产一区| 国产日韩一区二区三区精品不卡| 国产免费av片在线观看野外av| 日韩制服骚丝袜av| 每晚都被弄得嗷嗷叫到高潮| 日韩人妻精品一区2区三区| 一区二区三区四区激情视频| 18禁国产床啪视频网站| 国产片内射在线| 国产亚洲精品第一综合不卡| 国产av又大| 久久久久国产精品人妻一区二区| 亚洲精品中文字幕一二三四区 | 狠狠精品人妻久久久久久综合| 18在线观看网站| 中文字幕高清在线视频| 国产亚洲欧美在线一区二区| 久久久久久久精品精品| 中文字幕人妻丝袜一区二区| 黑人欧美特级aaaaaa片| 久久精品国产亚洲av香蕉五月 | 免费在线观看黄色视频的| 99热全是精品| 99国产精品一区二区蜜桃av | 一级毛片女人18水好多| 国产一卡二卡三卡精品| 国产男女内射视频| 狠狠精品人妻久久久久久综合| 亚洲精品中文字幕一二三四区 | 在线十欧美十亚洲十日本专区| 又紧又爽又黄一区二区| 亚洲 国产 在线| 欧美日本中文国产一区发布| 国产有黄有色有爽视频| 99久久国产精品久久久| 青青草视频在线视频观看| 999久久久精品免费观看国产| 国产日韩欧美在线精品| 精品久久久久久久毛片微露脸 | 亚洲情色 制服丝袜| 超碰97精品在线观看| 18禁国产床啪视频网站| 永久免费av网站大全| 18禁裸乳无遮挡动漫免费视频| 日韩免费高清中文字幕av| 中文字幕制服av| 桃红色精品国产亚洲av| 黑人猛操日本美女一级片| 九色亚洲精品在线播放| 国产在线视频一区二区| 国产精品成人在线| 日日爽夜夜爽网站| 成年动漫av网址| 亚洲欧美精品自产自拍| 国产不卡av网站在线观看| 欧美精品人与动牲交sv欧美| 一本—道久久a久久精品蜜桃钙片| 窝窝影院91人妻| 国产日韩欧美亚洲二区| 精品久久久久久久毛片微露脸 | 久久久久久亚洲精品国产蜜桃av| 真人做人爱边吃奶动态| 国产在线观看jvid| 一级,二级,三级黄色视频| 50天的宝宝边吃奶边哭怎么回事| 999久久久国产精品视频| 国产亚洲精品久久久久5区| 午夜成年电影在线免费观看| 999久久久精品免费观看国产| 国产一区二区三区在线臀色熟女 | www日本在线高清视频| 两个人免费观看高清视频| 青草久久国产| 亚洲欧美一区二区三区久久| 岛国毛片在线播放| 男男h啪啪无遮挡| 久久久久精品人妻al黑| 亚洲情色 制服丝袜| 亚洲成人国产一区在线观看| 人人澡人人妻人| 国产又色又爽无遮挡免| 久久久久精品国产欧美久久久 | 国产区一区二久久| 亚洲欧美色中文字幕在线| 十分钟在线观看高清视频www| 亚洲少妇的诱惑av| 国产区一区二久久| 巨乳人妻的诱惑在线观看| 亚洲免费av在线视频| 飞空精品影院首页| 国产成人精品无人区| e午夜精品久久久久久久| 亚洲精品美女久久av网站| 国产欧美日韩一区二区精品| 两个人看的免费小视频| 在线观看人妻少妇| 成人黄色视频免费在线看| 在线 av 中文字幕| 久久久水蜜桃国产精品网| 精品欧美一区二区三区在线| 国产精品久久久久久精品古装| 久久精品国产亚洲av高清一级| 国精品久久久久久国模美| av在线app专区| 亚洲人成电影观看| 国产老妇伦熟女老妇高清| 国产欧美亚洲国产| 亚洲精品国产一区二区精华液| 久久久久国产精品人妻一区二区| 国产一区二区激情短视频 | 国产精品99久久99久久久不卡| 人人妻人人澡人人爽人人夜夜| 国产成人av激情在线播放| 亚洲第一av免费看| 亚洲欧美成人综合另类久久久| 欧美日韩精品网址| 亚洲中文av在线| 老司机午夜福利在线观看视频 | 操出白浆在线播放| 老司机深夜福利视频在线观看 | 亚洲国产精品一区三区| 亚洲专区中文字幕在线| 亚洲av日韩在线播放| 欧美精品av麻豆av| 99国产极品粉嫩在线观看| 久久精品成人免费网站| 亚洲专区字幕在线| 一二三四社区在线视频社区8| 777米奇影视久久| 啦啦啦啦在线视频资源| 欧美精品人与动牲交sv欧美| 男男h啪啪无遮挡| 欧美精品高潮呻吟av久久| 99久久国产精品久久久| 黄网站色视频无遮挡免费观看| 欧美在线一区亚洲| 国产亚洲精品一区二区www | 精品福利观看| 一二三四社区在线视频社区8| 欧美日韩视频精品一区| 午夜福利视频精品| 丁香六月天网| 亚洲一区二区三区欧美精品| 热re99久久国产66热| 国产一卡二卡三卡精品| 中文字幕人妻熟女乱码| 精品福利观看| 久久精品人人爽人人爽视色| 国产av又大| 老司机影院成人| 亚洲精品美女久久av网站| 国产精品免费大片| 成人国产一区最新在线观看| 国产精品 国内视频| 在线观看免费高清a一片| 考比视频在线观看| 国产区一区二久久| 一二三四在线观看免费中文在| 久久女婷五月综合色啪小说| 午夜成年电影在线免费观看| 在线亚洲精品国产二区图片欧美| 伊人亚洲综合成人网| 咕卡用的链子| 黄频高清免费视频| 国产精品99久久99久久久不卡| 亚洲人成电影观看| 男人舔女人的私密视频| 十分钟在线观看高清视频www| 国产真人三级小视频在线观看| 国产男女内射视频| 亚洲成av片中文字幕在线观看| 久久国产亚洲av麻豆专区| 久久中文字幕一级| 色婷婷久久久亚洲欧美| 十八禁人妻一区二区| 12—13女人毛片做爰片一| 国产国语露脸激情在线看| 亚洲成人免费av在线播放| 大码成人一级视频| 亚洲成人国产一区在线观看| 久久人人97超碰香蕉20202| 妹子高潮喷水视频| 国产成人精品久久二区二区91| 少妇粗大呻吟视频| 亚洲av国产av综合av卡| 黑人巨大精品欧美一区二区mp4| 丝袜美腿诱惑在线| 精品国内亚洲2022精品成人 | 别揉我奶头~嗯~啊~动态视频 | 国产色视频综合| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看| 国产亚洲精品一区二区www | 精品国产超薄肉色丝袜足j| 伊人亚洲综合成人网| 国产精品久久久久成人av| 老司机影院毛片| 亚洲色图综合在线观看| 国产又色又爽无遮挡免| www.精华液| 美女中出高潮动态图| 亚洲第一欧美日韩一区二区三区 | 蜜桃国产av成人99| 亚洲五月婷婷丁香| 黑人巨大精品欧美一区二区mp4| av网站免费在线观看视频| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲综合色网址| 日韩欧美一区视频在线观看| 我的亚洲天堂| 国产精品自产拍在线观看55亚洲 | 一本大道久久a久久精品| 最近最新免费中文字幕在线| 亚洲国产成人一精品久久久| 首页视频小说图片口味搜索| 亚洲精品一二三| 精品国产乱子伦一区二区三区 | 亚洲精品成人av观看孕妇| 欧美日韩视频精品一区| av在线老鸭窝| 国产人伦9x9x在线观看| 一个人免费在线观看的高清视频 | 香蕉丝袜av| 97精品久久久久久久久久精品| 别揉我奶头~嗯~啊~动态视频 | 欧美久久黑人一区二区| 国产精品九九99| xxxhd国产人妻xxx| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| 9色porny在线观看| 精品一区二区三卡| 国产97色在线日韩免费| 99国产极品粉嫩在线观看| 99精品久久久久人妻精品| 一本大道久久a久久精品| 亚洲欧美精品自产自拍| 两性午夜刺激爽爽歪歪视频在线观看 | 成年人免费黄色播放视频| 可以免费在线观看a视频的电影网站| 少妇粗大呻吟视频| 国产国语露脸激情在线看| 999久久久国产精品视频| 欧美人与性动交α欧美精品济南到| 久久久久久免费高清国产稀缺| 欧美日韩视频精品一区| 啦啦啦啦在线视频资源| 欧美成狂野欧美在线观看| 免费在线观看影片大全网站| 亚洲五月婷婷丁香| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频| av免费在线观看网站| 丝袜人妻中文字幕| 国产亚洲一区二区精品| 嫁个100分男人电影在线观看| 色老头精品视频在线观看| 午夜福利在线免费观看网站| 黑人操中国人逼视频| 国产成人精品久久二区二区91| 亚洲专区国产一区二区| 亚洲激情五月婷婷啪啪| 亚洲午夜精品一区,二区,三区| 天天操日日干夜夜撸| 99久久99久久久精品蜜桃| 精品国产一区二区三区久久久樱花| 亚洲自偷自拍图片 自拍| 肉色欧美久久久久久久蜜桃| 一本大道久久a久久精品| 国产成人精品久久二区二区91| a级毛片黄视频| 丝袜人妻中文字幕| 久久久久久久国产电影| 免费在线观看黄色视频的| 99久久人妻综合| 日韩一卡2卡3卡4卡2021年| 日韩制服丝袜自拍偷拍| 中文精品一卡2卡3卡4更新| 中文字幕制服av| 国产麻豆69| 美女午夜性视频免费| av网站在线播放免费| 人人妻,人人澡人人爽秒播| 国产男女内射视频| av视频免费观看在线观看| av一本久久久久| 欧美黄色片欧美黄色片| 亚洲精品久久午夜乱码| 免费在线观看完整版高清| 亚洲av成人一区二区三| 亚洲人成77777在线视频| 天天躁日日躁夜夜躁夜夜| 成年美女黄网站色视频大全免费| 夜夜骑夜夜射夜夜干| 大型av网站在线播放| 不卡一级毛片| 老司机午夜十八禁免费视频| 久久人人97超碰香蕉20202| av又黄又爽大尺度在线免费看| 国产在线观看jvid| 免费av中文字幕在线| 十八禁网站网址无遮挡| 午夜福利影视在线免费观看| 成年美女黄网站色视频大全免费| a级毛片黄视频| 久久精品成人免费网站| 精品少妇一区二区三区视频日本电影| 一进一出抽搐动态| 天天影视国产精品| 少妇的丰满在线观看| 久久天堂一区二区三区四区| 欧美久久黑人一区二区| www.av在线官网国产| 午夜福利影视在线免费观看| 丝袜人妻中文字幕| 不卡一级毛片| 亚洲国产av新网站| 丁香六月天网| av网站在线播放免费| 黄片小视频在线播放| 精品国产国语对白av| 男女免费视频国产| 国产精品1区2区在线观看. | 日韩中文字幕视频在线看片| 亚洲自偷自拍图片 自拍| 久久天躁狠狠躁夜夜2o2o| 手机成人av网站| 国产欧美日韩一区二区精品| 美女视频免费永久观看网站| 欧美日韩视频精品一区| 亚洲av成人不卡在线观看播放网 | 天天躁日日躁夜夜躁夜夜| 国产亚洲精品久久久久5区| 丝袜脚勾引网站| 制服诱惑二区| 日韩欧美国产一区二区入口| 各种免费的搞黄视频| 另类亚洲欧美激情| 制服人妻中文乱码| www.999成人在线观看| 9色porny在线观看| 欧美激情高清一区二区三区| 国产精品1区2区在线观看. | 又紧又爽又黄一区二区| 中文字幕人妻熟女乱码| 日本91视频免费播放| 成人影院久久| 国产一区二区激情短视频 | 精品人妻熟女毛片av久久网站| 色播在线永久视频| 欧美在线一区亚洲| 性色av一级| 国产麻豆69|