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

    基于EC法的風(fēng)浪聯(lián)合作用主塔-基礎(chǔ)體系極限荷載效應(yīng)

    2017-11-06 02:29:41涂志斌黃銘楓樓文娟
    振動(dòng)與沖擊 2017年19期
    關(guān)鍵詞:有限元效應(yīng)

    涂志斌,黃銘楓,樓文娟

    (浙江大學(xué) 建筑工程學(xué)院結(jié)構(gòu)工程研究所,杭州 310058)

    基于EC法的風(fēng)浪聯(lián)合作用主塔-基礎(chǔ)體系極限荷載效應(yīng)

    涂志斌,黃銘楓,樓文娟

    (浙江大學(xué) 建筑工程學(xué)院結(jié)構(gòu)工程研究所,杭州 310058)

    針對(duì)基于Rosenblatt映射變換的極限狀態(tài)曲線(xiàn)計(jì)算困難的現(xiàn)狀,提出了多維隨機(jī)變量極限狀態(tài)曲線(xiàn)的簡(jiǎn)化算法。該算法通過(guò)尋找同時(shí)滿(mǎn)足邊緣分布、聯(lián)合分布和可靠指標(biāo)的變量估計(jì)值,將隨機(jī)變量條件分布函數(shù)及其逆函數(shù)的求解轉(zhuǎn)化為邊緣分布函數(shù)的求逆,從而達(dá)到簡(jiǎn)化極限狀態(tài)曲線(xiàn)計(jì)算的目的。以某跨海大橋橋塔-基礎(chǔ)體系為工程實(shí)例,以潿洲島海洋站風(fēng)浪同步觀(guān)測(cè)資料為環(huán)境變量統(tǒng)計(jì)資料,通過(guò)基于Copula函數(shù)的聯(lián)合分布模型構(gòu)造了風(fēng)浪聯(lián)合分布函數(shù),根據(jù)簡(jiǎn)化算法計(jì)算了風(fēng)浪極限狀態(tài)曲線(xiàn),利用EC法估計(jì)了橋塔-基礎(chǔ)體系的基底剪力極限荷載效應(yīng),并與外推法的估計(jì)結(jié)果進(jìn)行了比較。結(jié)果表明,基于簡(jiǎn)化算法和EC法的跨海大橋橋塔-基礎(chǔ)體系基底剪力極限荷載效應(yīng)具有較高的準(zhǔn)確性。

    極限荷載效應(yīng);極限狀態(tài)曲線(xiàn);風(fēng)浪聯(lián)合作用;EC法;Copula函數(shù)

    跨海大橋在施工和運(yùn)營(yíng)過(guò)程中將承受復(fù)雜的、隨時(shí)間變化的隨機(jī)環(huán)境作用,主要包括風(fēng)和波浪。為使結(jié)構(gòu)在使用期內(nèi)正常發(fā)揮預(yù)定功能,采用極限荷載效應(yīng)進(jìn)行結(jié)構(gòu)設(shè)計(jì)十分必要[1-2]。

    在極限荷載效應(yīng)估計(jì)中,構(gòu)造聯(lián)合分布函數(shù)是考慮多維隨機(jī)環(huán)境變量聯(lián)合作用的有效方法。多位學(xué)者采用傳統(tǒng)聯(lián)合分布模型來(lái)構(gòu)造多維隨機(jī)環(huán)境變量的聯(lián)合分布函數(shù)[3-7]。然而該模型要求變量服從同類(lèi)邊緣分布,相關(guān)關(guān)系采用Pearson線(xiàn)性相關(guān)系數(shù)表達(dá),適用范圍有限?;贑opula函數(shù)的聯(lián)合分布模型彌補(bǔ)了傳統(tǒng)聯(lián)合分布模型的不足。該模型將邊緣分布和相關(guān)關(guān)系分開(kāi)考慮,不要求變量服從同類(lèi)邊緣分布,可描述任意相關(guān)關(guān)系,具有極強(qiáng)的靈活性,在多維變量聯(lián)合分布函數(shù)構(gòu)造中的應(yīng)用十分廣泛[8-15]。

    在基于概率的結(jié)構(gòu)設(shè)計(jì)方法中,外推法是多維隨機(jī)環(huán)境變量聯(lián)合作用下極限荷載效應(yīng)分析的有效方法,具有較高的可信度和適應(yīng)性,但需在環(huán)境變量分布范圍內(nèi)進(jìn)行結(jié)構(gòu)動(dòng)力有限元分析,計(jì)算效率較低[16]。為提高極限荷載效應(yīng)的估計(jì)效率,Winterstein等[17-18]在外推法的基礎(chǔ)上提出了EC(Environmental Contour)法。EC法的核心思想是在具有指定超越概率的環(huán)境變量極限狀態(tài)曲面上搜尋極限荷載效應(yīng)及其對(duì)應(yīng)的環(huán)境變量組合點(diǎn),結(jié)構(gòu)動(dòng)力有限元分析只需在極限狀態(tài)曲面上完成,效率大大提高。Haver等[19-23]利用EC法估計(jì)了多維隨機(jī)環(huán)境變量聯(lián)合作用下工程結(jié)構(gòu)的極限荷載效應(yīng)。然而在EC法中,極限狀態(tài)曲面的估計(jì)依賴(lài)于環(huán)境變量條件分布函數(shù)及其逆函數(shù)的求解。當(dāng)聯(lián)合分布函數(shù)較為復(fù)雜時(shí),條件分布函數(shù)及其逆函數(shù)的求解十分困難,極限狀態(tài)曲面的計(jì)算難以實(shí)現(xiàn),EC法的適應(yīng)性較低。Montes-Iturrizaga欲根據(jù)Copula函數(shù)的性質(zhì)提出極限狀態(tài)曲面的簡(jiǎn)化算法[24],但在推導(dǎo)過(guò)程中不恰當(dāng)?shù)乩昧寺?lián)合分布函數(shù)與Copula函數(shù)的關(guān)系。

    為提高EC法的適應(yīng)性,本文提出了極限狀態(tài)曲面的簡(jiǎn)化算法。以某跨海大橋橋塔-基礎(chǔ)體系為工程實(shí)例,以潿洲島海洋站的風(fēng)浪同步觀(guān)測(cè)資料為環(huán)境變量統(tǒng)計(jì)資料,通過(guò)Copula函數(shù)構(gòu)造了風(fēng)浪聯(lián)合分布函數(shù),根據(jù)簡(jiǎn)化算法計(jì)算了風(fēng)浪極限狀態(tài)曲線(xiàn),利用EC法估計(jì)了橋塔-基礎(chǔ)體系的基底剪力極限荷載效應(yīng),并與外推法的估計(jì)結(jié)果進(jìn)行了比較,討論了EC法的適用范圍。

    1 基于Copula函數(shù)的聯(lián)合分布模型

    Copula函數(shù)是一種構(gòu)造聯(lián)合分布的工具,其優(yōu)勢(shì)是將變量間的相關(guān)關(guān)系與邊緣分布分開(kāi)考慮,能適應(yīng)任意多維隨機(jī)變量的聯(lián)合分布分析。二維隨機(jī)變量聯(lián)合分布函數(shù)與Copula函數(shù)的關(guān)系為

    F(s1,s2)=C[F1(s1),F2(s2),θ]

    (1)

    式中:F(s1,s2)為隨機(jī)變量s1、s2的聯(lián)合分布函數(shù);C[·]為Copula函數(shù);F1(s1)、F2(s2)為s1、s2的邊緣分布函數(shù);θ為二者的相關(guān)系數(shù),能衡量隨機(jī)變量間的任意相關(guān)關(guān)系。對(duì)式(1)等號(hào)兩側(cè)同時(shí)求二階混合偏導(dǎo),可得到聯(lián)合概率密度函數(shù)f(s1,s2)與Copula函數(shù)的關(guān)系

    (2)

    在眾多Copula函數(shù)中,同時(shí)滿(mǎn)足以下兩個(gè)條件的Copula函數(shù)稱(chēng)為最優(yōu)Copula:①v1、v2能準(zhǔn)確描述s1、s2的邊緣分布特性;②θ能準(zhǔn)確描述各隨機(jī)變量間的相關(guān)關(guān)系。目前常用的最優(yōu)Copula評(píng)價(jià)準(zhǔn)則為AIC準(zhǔn)則和BIC準(zhǔn)則[26-27]。二者則均建立在Copula函數(shù)參數(shù)估計(jì)的基礎(chǔ)上,表達(dá)式分別為

    (3)

    (4)

    表1 常用的Copula函數(shù)Tab.1 Common Copula functions

    2 基于EC法的極限荷載效應(yīng)

    2.1EC法的基本原理

    EC法的核心思想是在超越概率為pE的環(huán)境變量極限狀態(tài)曲面上搜尋極限荷載效應(yīng)及其對(duì)應(yīng)的環(huán)境變量組合點(diǎn)。與外推法在整個(gè)環(huán)境變量分布范圍內(nèi)完成結(jié)構(gòu)動(dòng)力有限元分析不同,EC法只需要在極限狀態(tài)曲面上完成分析,計(jì)算效率大大提高。EC法有兩個(gè)基本假定:①荷載效應(yīng)極值的變異性可以忽略,②環(huán)境變量極限狀態(tài)曲面為線(xiàn)性凸曲面。根據(jù)假定,超越概率為pE的極限荷載效應(yīng)是相同超越概率下極限狀態(tài)曲面上的最大荷載效應(yīng)極值中位值。

    本節(jié)仍以二維隨機(jī)變量來(lái)說(shuō)明EC法的基本原理和實(shí)施步驟,此時(shí)極限狀態(tài)曲面退化為極限狀態(tài)曲線(xiàn)。對(duì)于具有相關(guān)性的隨機(jī)變量s1、s2,直接在物理空間中計(jì)算二者的極限狀態(tài)曲線(xiàn)并不容易。根據(jù)Rosenblatt映射變換,具有任意相關(guān)關(guān)系的隨機(jī)變量都可轉(zhuǎn)換為標(biāo)準(zhǔn)正態(tài)空間中的獨(dú)立隨機(jī)變量,即

    Φ1(u1)=F1(s1),Φ2(u2)=F2|1(s2)

    (5)

    (6)

    式中:u1、u2為標(biāo)準(zhǔn)正態(tài)空間中相互獨(dú)立的隨機(jī)變量,服從標(biāo)準(zhǔn)正態(tài)分布;F2|1(s2)為s2的條件分布函數(shù),與s2的邊緣分布特性及s1和s2的相關(guān)關(guān)系有關(guān)。根據(jù)式(5)和式(6),標(biāo)準(zhǔn)正態(tài)空間中的極限狀態(tài)曲線(xiàn)u1-u2可轉(zhuǎn)換為物理空間中的極限狀態(tài)曲線(xiàn)s1-s2。若將s2看作主變量,s1看作條件變量,式(6)可改寫(xiě)為

    (7)

    理論上式(6)和式(7)是等價(jià)的,以s2為主變量、s1為條件變量不會(huì)改變極限狀態(tài)曲線(xiàn)的估計(jì)結(jié)果。因此本節(jié)分析在式(6)的基礎(chǔ)上展開(kāi)。在標(biāo)準(zhǔn)正態(tài)空間中,極限狀態(tài)曲線(xiàn)u1-u2為半徑等于可靠指標(biāo)βT的圓

    (8)

    βT=-Φ-1(PE)

    (9)

    式中,ω為角坐標(biāo)。式(9)說(shuō)明了可靠指標(biāo)βT與超越概率PE的關(guān)系。

    圖1詳細(xì)說(shuō)明了風(fēng)、浪聯(lián)合作用下基于EC法的極限荷載效應(yīng)lT的計(jì)算流程,包含4個(gè)主要步驟,如4個(gè)虛線(xiàn)框所示。①根據(jù)重現(xiàn)期T計(jì)算超越概率PE和可靠指標(biāo)βT。②在[0,2π)范圍內(nèi)離散ω,離散點(diǎn)數(shù)為m;根據(jù)式(8)計(jì)算標(biāo)準(zhǔn)正態(tài)空間中的極限狀態(tài)曲線(xiàn)離散點(diǎn)(u1,i,u2,i);根據(jù)式(6)轉(zhuǎn)換為物理空間中的離散點(diǎn)(ui,hi)。③調(diào)整隨機(jī)數(shù)cj來(lái)完成離散點(diǎn)(ui,hi)處的隨機(jī)流場(chǎng)同步模擬,包括脈動(dòng)風(fēng)速和隨機(jī)波浪;生成隨機(jī)荷載時(shí)程,完成結(jié)構(gòu)動(dòng)力有限元計(jì)算,并根據(jù)POT法提取荷載效應(yīng)極值Lj。當(dāng)k次結(jié)構(gòu)動(dòng)力有限元計(jì)算完成時(shí),擬合Lj的分布函數(shù),取分布函數(shù)值為0.5時(shí)的荷載效應(yīng)極值為中位值Lm,i。④重復(fù)③直至i=m,在m個(gè)中位值Lm,i中尋找最大值,作為超越概率為PE的極限荷載效應(yīng)lT。在該流程中離散點(diǎn)(ui,hi)的計(jì)算是難點(diǎn)。在EC法中,極限狀態(tài)曲線(xiàn)上的結(jié)構(gòu)動(dòng)力有限元計(jì)算次數(shù)為k×m。為進(jìn)一步提高極限荷載效應(yīng)lT的估計(jì)效率并保證準(zhǔn)確性,極限狀態(tài)曲線(xiàn)的離散應(yīng)遵循以下原則:在中位值Lm,i較大的區(qū)域加密離散點(diǎn),在Lm,i較小的區(qū)域減少甚至不設(shè)置離散點(diǎn)。

    圖1 基于EC法的極限荷載效應(yīng)lT計(jì)算流程Fig.1 Flow chart of extreme load effect lT based on EC method

    一般而言,基于POT(Peak Over Threshold)的數(shù)據(jù)樣本服從廣義帕累托分布(Generalized Pareto Distribution,GPD)或三參數(shù)威布爾分布(Three-Parameter Weibull Distribution,W3P)[28],二者的表達(dá)式分別為

    GPD

    F0,POT(lT)=1-(1-kglT/ag)kg≠0

    F0,POT(lT)=1-exp(-lT/ag)kg=0

    (10)

    W3P

    (11)

    式中:ag、aw為尺度參數(shù);kg、kw為形狀參數(shù);μw為位置參數(shù)。荷載效應(yīng)極值分布函數(shù)F0(lT)與F0,POT(lT)的關(guān)系為

    F0(lT)=[F0,POT(lT)]nPOT

    (12)

    式中,nPOT為荷載效應(yīng)極值L的樣本容量。

    2.2極限狀態(tài)曲線(xiàn)的簡(jiǎn)化算法

    (13)

    結(jié)合式(1)和式(5)有

    (14)

    圖2 網(wǎng)格劃分示意圖
    Fig.2 Grid meshing

    3 工程背景

    3.1風(fēng)浪同步觀(guān)測(cè)資料

    (15)

    式中:s=u、h;Fs(s)為邊緣分布函數(shù);μs為位置參數(shù);σs為尺度參數(shù)。子樣1、子樣2的邊緣分布參數(shù)擬合結(jié)果見(jiàn)表2。

    表2 邊緣分布函數(shù)的參數(shù)擬合值Tab.2 Fitted values of marginal CDF

    從表1中選擇最優(yōu)Copula來(lái)構(gòu)造子樣1、子樣2的聯(lián)合分布函數(shù)。采用極大似然法估計(jì)Copula函數(shù)的參數(shù),采用AIC準(zhǔn)則和BIC準(zhǔn)則選擇最優(yōu)Copula,結(jié)果見(jiàn)表3和表4。對(duì)于子樣1、子樣2,兩種準(zhǔn)則的評(píng)價(jià)結(jié)果一致,即Gaussian Copula為最優(yōu)Copula,參數(shù)θ的擬合值分別為0.745 5、0.605 6。圖3和圖4為子樣1、子樣2的聯(lián)合概率密度,符號(hào)Y表示考慮風(fēng)浪相關(guān)性(θ≠0),N表示不考慮風(fēng)浪相關(guān)性(θ=0)。由圖可知,子樣1/2-Y的聯(lián)合概率密度峰值大于子樣1/2-N。根據(jù)Turkstra法則及其拓展[29-30],子樣1、子樣2為所有樣本的上下邊界。因此基于子樣1、子樣2的極限荷載效應(yīng)是基于其他樣本的極限荷載效應(yīng)的上下邊界。

    表3Copula函數(shù)的擬合參數(shù)、AIC及BIC:子樣1
    Tab.3Fttedparameters,AICandBICofCopulafunctionsofsubsample1

    模型θ或θ/λAICBICGumbel-1.9982-13.8458-12.4785Frank-5.2697-11.2529-9.8856Clayton-2.0411-14.4579-13.0906Gaussian-0.7455-14.5385-13.1712t-0.6471/1.4007-13.8336-11.0990Galambos-1.2909-13.9756-12.6083HuslerReiss-1.8369-14.1997-12.8324

    3.2橋塔-基礎(chǔ)體系

    某跨海大橋橋塔-基礎(chǔ)為鋼筋混凝土結(jié)構(gòu),混凝土

    表4Copula函數(shù)的擬合參數(shù)、AIC及BIC:子樣2
    Tab.4Fittedparameters,AICandBICofCopulafunctionsofsubsample2

    模型θ或θ/λAICBICGumbel-1.6187-6.4222-5.0549Frank-4.1293-6.7862-5.4190Clayton-1.1386-6.2017-4.8344Gaussian-0.6056-6.8906-5.5233t-0.6030/6.1884-5.4983-2.7637Galambos-0.8928-6.1899-4.8226HuslerReiss-1.2865-5.7686-4.4013

    (a) 子樣1-Y

    (b) 子樣1-N圖3 子樣1的聯(lián)合概率密度函數(shù)Fig.3 JPDF of subsample 1

    (a) 子樣2-Y

    (b) 子樣2-N圖4 子樣2的聯(lián)合概率密度函數(shù)Fig.4 JPDF of subsample

    的強(qiáng)度等級(jí)為C40。橋塔為鉆石型,高460 m,位于水面以上;橋塔89 m處設(shè)有4道高9 m、寬4 m的橫梁;橋塔截面尺寸沿高度逐漸變化,其中x、y方向上塔底中心間距為28 m、40 m,塔底尺寸為20 m、16 m,塔頂尺寸為15 m、14 m;橫梁以下塔腿內(nèi)、外側(cè)x方向上的傾斜度為5.70∶1、8.39∶1,y方向上的傾斜度為4.78∶1、6.21∶1,橫梁以上塔腿內(nèi)、外側(cè)x方向上的傾斜度為15.46∶1、17.73∶1,y方向上的傾斜度為11.53∶1、10.95∶1?;A(chǔ)由沉井和承臺(tái)組成,位于水面以下。沉井為大直徑圓形沉井,直徑90 m,井壁厚2.5 m;沉井內(nèi)部x、y方向上等間距地設(shè)有5道厚度為1.5 m的隔墻;沉井頂部為承臺(tái),厚度7 m;泥面以上基礎(chǔ)的總高度為50 m。借助ANSYS軟件建立橋塔-基礎(chǔ)體系的有限元模型,如圖5。橋塔采用空間梁?jiǎn)卧M,沉井和承臺(tái)采用三維實(shí)體單元模擬,在橋塔和承臺(tái)的接觸面處建立剛性區(qū)域以傳遞自由度。沉井底部采用固定端約束。表5為橋塔-基礎(chǔ)體系的前三階模態(tài)信息。

    表5 橋塔-基礎(chǔ)體系的模態(tài)信息Tab.5 Mode shapes of bridge-tower system

    圖5 橋塔-基礎(chǔ)體系有限元模型Fig.5 Finite element model of tower-basement system

    4 基底剪力極限荷載效應(yīng)

    4.1極限狀態(tài)曲線(xiàn)

    (a) 子樣1

    (b) 子樣2圖6 子樣1、子樣2的極限狀態(tài)曲線(xiàn)Fig.6 Limit state lines of subsample 1 and 2

    4.2橋塔-基礎(chǔ)體系的輸入荷載

    (a) 主變量

    (b) 條件變量Hs圖7 子樣1平均風(fēng)速和有效波高Hs離散點(diǎn)Fig.7 Discrete points of and Hs for subsample 1

    (a) 條件變量

    (b) 主變量Hs圖8 子樣2平均風(fēng)速和有效波高Hs離散點(diǎn)Fig.8 Discrete points of and Hs for subsample 2

    圖9 脈動(dòng)風(fēng)速時(shí)程Fig.9 History of fluctuating wind speed: =47.4 m/s

    (a) 塔腿部分

    (b) 塔頂部分圖10 橋塔阻力系數(shù)CDFig.10 Drag force coefficients CD of bridge tower

    隨機(jī)波浪的模擬參數(shù)見(jiàn)表6,Ts為有效周期,γ為譜峰因子。圖11為模擬隨機(jī)波面時(shí)程,此時(shí)Hs=10.8 m(T=100a、ω=18.75°,考慮風(fēng)浪相關(guān)性)。作用在基礎(chǔ)上的隨機(jī)波浪壓力可根據(jù)繞射理論計(jì)算。

    表6 隨機(jī)波浪模擬參數(shù)Tab.6 Simulation parameters of random wave height

    在時(shí)域內(nèi)完成橋塔-基礎(chǔ)體系的動(dòng)力有限元分析,橋塔的輸入荷載為風(fēng)荷載時(shí)程,基礎(chǔ)的輸入荷載為隨機(jī)波浪壓力時(shí)程,分析方法為完全瞬態(tài)法。根據(jù)DNV(Det Norske Veritas)[32]的研究成果,風(fēng)浪聯(lián)合作用時(shí)風(fēng)向與波向的差異很小??绾4髽蛩诘氐臍庀筚Y料顯示,橋位處的常年風(fēng)向以NW為主,與模型坐標(biāo)y軸一致。因此風(fēng)攻角和波浪入射方向均與y軸一致。

    圖11 隨機(jī)波浪時(shí)程:Hs=10.8 mFig.11 History of random wave height: Hs=10.8 m

    4.3荷載效應(yīng)極值分布函數(shù)

    根據(jù)IEC 61400-1的建議,離散點(diǎn)處的結(jié)構(gòu)動(dòng)力有限元重復(fù)計(jì)算次數(shù)取k=6,閾值取6條荷載效應(yīng)時(shí)程的均值與1.4倍標(biāo)準(zhǔn)差之和。結(jié)合極限狀態(tài)曲線(xiàn)的離散點(diǎn)數(shù),本文共完成了k×m=150次有限元計(jì)算。

    圖12 基底剪力時(shí)程Fig.12 History of base shear force: =47.4 m/s、Hs=10.8 m

    圖13 荷載效應(yīng)分布函數(shù)Fig.13 CDF of base shear force: =47.4 m/s、Hs=10.8 m

    4.4極限荷載效應(yīng)

    圖14為各離散點(diǎn)處荷載效應(yīng)極值中位值Lm,i的變化規(guī)律。由圖可知:①隨離散點(diǎn)的變化,Lm,i無(wú)突變且變化軌跡為拋物線(xiàn);②考慮風(fēng)浪相關(guān)性時(shí),各重現(xiàn)期下子樣1、子樣2的max(Lm,i)分別出現(xiàn)在ω=22.5°、78.75°;③不考慮風(fēng)浪相關(guān)性時(shí),各重現(xiàn)期下子樣1、子樣2的max(Lm,i)分別出現(xiàn)在ω=78.75°、84.375°。這表明極限狀態(tài)曲線(xiàn)離散點(diǎn)的設(shè)置是合理的。為了進(jìn)一步減小計(jì)算量,提高極限荷載效應(yīng)lT的估計(jì)效率,離散點(diǎn)的分布范圍和數(shù)量可進(jìn)一步減少。

    (a) 子樣1

    (b) 子樣2圖14 各離散點(diǎn)處的荷載效應(yīng)極值中位值Lm,iFig.14 Medians Lm,i of peak load effect at discrete points

    基于EC法的極限荷載效應(yīng)lT將與基于外推法的lT進(jìn)行對(duì)比。由于不引入任何假定,基于外推法的計(jì)算結(jié)果具有較高的可信度,可作為基于EC法的計(jì)算結(jié)果的對(duì)比標(biāo)準(zhǔn)?;谕馔品ǖ臉蛩?基礎(chǔ)體系基底剪力極限荷載效應(yīng)lT的具體實(shí)現(xiàn)和計(jì)算結(jié)果參考涂志斌等[33]的文章。涂志斌等指出采用外推法估計(jì)lT時(shí),風(fēng)浪分布范圍內(nèi)的橋塔-基礎(chǔ)體系動(dòng)力有限元計(jì)算次數(shù)為432,而本文的計(jì)算次數(shù)為150。因此與外推法相比,EC法的計(jì)算效率顯著提高。事實(shí)上,通過(guò)合理地設(shè)置離散點(diǎn),橋塔-基礎(chǔ)體系的動(dòng)力有限元計(jì)算次數(shù)還可進(jìn)一步降低,EC法的計(jì)算效率也可進(jìn)一步提高。

    (a) 子樣1-N

    (b) 子樣1-Y圖15 荷載效應(yīng)極值中位值Lm,i的等值線(xiàn):子樣1Fig.15 Contours of peak load effect medians Lm,i: subsample 1

    (a) 子樣2-N

    表8為各重現(xiàn)期下基于子樣2、外推法和EC法的橋塔-基礎(chǔ)體系基底剪力極限荷載效應(yīng)lT,各符號(hào)的意義與表7相同。對(duì)表8進(jìn)行分析可得到與表7類(lèi)似的結(jié)論,不再贅述。對(duì)于子樣1、子樣2,極限荷載效應(yīng)lT的計(jì)算結(jié)果并不相同。在實(shí)際工程中為了保證結(jié)構(gòu)安全,應(yīng)選擇兩者中較大的極限荷載效應(yīng)進(jìn)行結(jié)構(gòu)設(shè)計(jì)。

    表7 子樣1的極限荷載效應(yīng)lTTab.7 Extreme load effects lT for subsample 1 ×105 kN

    表8 子樣2的極限荷載效應(yīng)lTTab.8 Extreme load effects lT for subsample 2 ×105 kN

    5 結(jié) 論

    針對(duì)基于Rosenblatt映射變換的極限狀態(tài)曲線(xiàn)計(jì)算困難的現(xiàn)狀,本文提出了極限狀態(tài)曲線(xiàn)的簡(jiǎn)化算法,提高了EC法的可行性。根據(jù)簡(jiǎn)化算法和EC法完成了某跨海大橋橋塔-基礎(chǔ)體系基底剪力極限荷載效應(yīng)估計(jì),并與外推法的計(jì)算結(jié)果進(jìn)行了對(duì)比,得到了以下結(jié)論:

    (1) 簡(jiǎn)化算法能有效的估計(jì)風(fēng)浪的極限狀態(tài)曲線(xiàn);是否考慮風(fēng)浪相關(guān)性對(duì)極限狀態(tài)曲線(xiàn)有顯著影響。

    (2) 忽略風(fēng)浪相關(guān)性使極限荷載效應(yīng)的估計(jì)值偏小。

    (3) 與外推法相比,EC法的計(jì)算效率顯著提高。EC法的計(jì)算效率與極限狀態(tài)曲線(xiàn)離散點(diǎn)的設(shè)置有關(guān)。為進(jìn)一步提高EC法的計(jì)算效率,極限狀態(tài)曲線(xiàn)的離散點(diǎn)可僅設(shè)置在使結(jié)構(gòu)荷載效應(yīng)較大的區(qū)域。

    (4) 基于EC法的極限荷載效應(yīng)略小于基于外推法的極限荷載效應(yīng),誤差出現(xiàn)的原因是忽略荷載效應(yīng)極值的變異性。為降低荷載效應(yīng)極值的變異性,進(jìn)一步減小基于EC法的極限荷載效應(yīng)的估計(jì)誤差,可適當(dāng)提高統(tǒng)計(jì)樣本的提取閾值。由于估計(jì)誤差較小,基于EC法的極限荷載效應(yīng)可直接用于工程結(jié)構(gòu)設(shè)計(jì)。

    [1] Design of offshore wind turbine structure: DNY-OS-J101[S]. Oslo: Det Norske Veritas, 2010: 9-22.

    [2] Wind turbines-Part 1: design requirements: IEC 61400-1[S]. Geneva: Electrotechnical Commission, 2005: 19-32.

    [3] 李鋒. 海洋工程雙變量環(huán)境條件設(shè)計(jì)參數(shù)估計(jì)[D]. 青島:中國(guó)海洋大學(xué), 2005:17-40.

    [4] 周道成, 段忠東. 耿貝爾邏輯模型在極值風(fēng)速和有效波高聯(lián)合概率分布中的應(yīng)用[J]. 海洋工程, 2003, 21(2): 45-51.

    ZHOU Daocheng, DUAN Zhongdong. The Gumbel-logistic model for joint probability distribution of extreme-value wind speeds and effective wave heights [J]. The Ocean Engineering, 2003, 21(2): 45-51.

    [5] 歐進(jìn)萍, 肖儀清, 段忠東, 等. 基于風(fēng)浪聯(lián)合概率模型的海洋平臺(tái)結(jié)構(gòu)系統(tǒng)可靠度分析[J]. 海洋工程, 2003, 21(4): 1-7.

    OU Jinping, XIAO Yiqing, DUAN Zhongdong, et al. Structual system reliability analysis for offshore platforms based on the joint probabilistic model of wind and wave[J]. The Ocean Engineering, 2003, 21(4): 1-7.

    [6] YUE S. The Gumbel logistic model for representing a multivariate storm event [J]. Advance in Water Resources, 2001, 24(2): 179-185.

    [7] 劉偉. 基于最大熵分布的海洋平臺(tái)環(huán)境條件聯(lián)合重現(xiàn)值推算[D]. 青島: 中國(guó)海洋大學(xué), 2011: 64-87.

    [8] GENES C, FAVRE A C. Everything you wanted to know about Copula modeling but were afraid to ask [J]. Journal of Hydrologic Engineering, 2007, 12(4):347-368.

    [9] WIST H T, MYRHAUG D, RUE H. Statistical properties of successive wave heights ans successive wave periods [J]. Applied Ocean Research, 2005, 26(3/4): 114-136.

    [10] SILVA-GONZLEZ F, HEREDIA-ZAVONI E, MONTES-ITURRIZAGA R. Development of environmental contours using Nataf distribution model [J]. Ocean Engineering, 2013, 58:27-34.

    [11] 陶山山. 多維最大熵模型及其在海岸及海洋工程中的應(yīng)用[D]. 青島:中國(guó)海洋大學(xué), 2013: 76-113.

    [12] 董勝, 翟金金, 陶山山. 基于Archimedean Copula函數(shù)的風(fēng)浪聯(lián)合統(tǒng)計(jì)分析[J]. 中國(guó)海洋大學(xué)學(xué)報(bào), 2014, 44(10): 134-141.

    DONG Sheng, ZHAI Jinjin, TAO Shanshan. The joint statistical analysis of wind and wave based on Archimedean Copula functions [J]. Periodical of Ocean University of China, 2014, 44(10): 134-141.

    [13] 陳子燊. 波高與風(fēng)速聯(lián)合概率分布研究[J]. 海洋通報(bào), 2011, 30(2): 158-163.

    CHEN Zishen. Study on joint probability distribution of wave height and wind velocity[J]. Marine Science Bulletin, 2011, 30(2): 158-163.

    [14] YANG X C, ZHANG Q H. Joint probability distribution of winds and waves from wave simulation of 20 years (1989-2008) in Bohai Bay [J]. Water Science and Engineering, 2013, 6(3): 296-307.

    [15] ZHANG Y, BEER M, QUEK S T. Long-term performance assessment and design of offshore structures [J]. Computers and Structures, 2015, 154:101-115.

    [16] Wind turbines-Part 3: design requirements for offshore wind turbines: IEC 61400-3[S]. [S.l.]: Electrotechnical Commission, TC88 WG3 Committee Draft, 2005: 18-21.

    [17] WINTERSTEIN S R, UDE T C, CORNELL C A,et al. Environmental parameters for extreme response: inverse FORM with omission factors: ICOSSAR-93[R]. Innsbruck: [s.n.], 1993.

    [18] WINTERSTEIN S R, ENGEBRETSEN K. Reliability-based prediction of design loads and responses for floating ocean structures [C]∥Appeared, Proceedings, 17th Intl. Conf. on Offshore Mechanics and Arctic Engineering. Lisbon: OMAR, 1998.

    [19] HAVER S, WINTERSTEIN S R. Environmental contour lines: a method for estimating long term extremes by a short term analysis[C]∥ Transactions-Society of Naval Architects and Marine Engineers.[S.l.]:[s.n.], 2008: 116.

    [20] SARANYASOONTORN K, MANUEL L. Design loads for wind turbines using the environmental contour method [J]. Journal of Solar Energy Engineering, 2006, 128(4):554-561.

    [21] AGARWAL P, MANUEL L. Simulation of offshore wind turbine response for long-term extreme load prediction[J]. Engineering Structures, 2009, 31(10):2236-2246.

    [22] HUSEBY A B, VANEM E, NATVIG B. A new approach to environmental contours for ocean engineering applications based on direct Monte Carlo simulations[J]. Ocean Engineering, 2013, 60(60):124-135.

    [23] HUSEBY A B, VANEM E, NATVIG B. Alternative environmental contours for structural reliability analysis[J]. Structural Safety, 2015, 54:32-45.

    [24] MONTES-ITURRIZAGA R, HEREDIA-ZAVONI E. Environmental contours using copulas[J]. Applied Ocean Research, 2015, 52:125-139.

    [25] GHORBEL A, TRABELSI A. Energy portfolio risk management using time-varying extreme value copula methods[J]. Economic Modelling, 2014, 38: 470-485.

    [26] TRIVEDI P K, ZIMMER D M. Copula modeling: an introduction for practitioners. Foundations and Trends in Econometrics [M]. [S.l.]: Now Publishers, 2007: 53-70.

    [27] SUNDARESAN A, VARSHNEY P K. Location estimation of a random signal source based on correlated sensor observations[J]. Signal Processing, IEEE Transactions on, 2011, 59(2): 787-799.

    [28] RAGAN P, MANUEL L. Statistical extrapolation methods for estimating wind turbine extreme loads[J]. Journal of Solar Energy Engineering, 2008, 130(3): 1-19.

    [29] TURKSTRA C J, MADSEN H O. Load combinations in condified structural design [J]. Journal of the Structural Division, 1980, 106(12): 2527-2543.

    [30] YEO D. Multiple points-in-time estimation of peak wind effects on structures [J]. Journal of Structural Engineering, 2013, 139(3): 462-471.

    [31] MYRHAUG D, ONG M C. Effect of wave age on wind gust spectra over wind waves[J]. Journal of Offshore Mechanics and Arctic Engineering, 2009, 131(3): 034501.

    [32] Environmental conditions and environmental loads Recommended Practice: DNV-RP-C205[S]. [S.l.]: DNV, 2012.

    [33] 涂志斌, 黃銘楓, 樓文娟. 風(fēng)浪耦合作用下橋塔-基礎(chǔ)體系的極限荷載效應(yīng)[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2016, 50(4): 601-610.

    TU Zhibin, HUANG Mingfeng, LOU Wenjuan. Extreme load effects on a bridge tower-basement system due to the joint actions of wind and wave [J]. Journal of Zhejinag University (Engineering Science), 2016, 50(4): 601-610.

    Extremeloadeffectsonabridgetower-basementsystemunderthejointactionofwindandwavebasedontheECmethod

    TU Zhibin,HUANG Mingfeng,LOU Wenjuan

    (Institute of Structural Engineering, Zhejiang University, Hangzhou 310058, China)

    A simplified algorithm was developed to construct the limit state line of multi random variables, aiming at overcoming the difficulty in the construction of limit state line by Rosenblatt transformation. By searching the expected values of variables which can satisfy the marginal distributions, joint cumulative distribution and reliability index simultaneously, the calculation of conditional distribution function and its inverse function was turned to find the inverse function of marginal distribution, and thus the construction of the limit state line was simplified. Taking some bridge tower-basement system as an engineering example and the wind-wave simultaneous observation data at the Weizhou marine station as statistic samples, the joint cumulative distribution function of wind and wave was estimated by an Copula function, the limit state line of wind and wave was constructed by the proposed algorithm, and the extreme load effects of the base shear force were estimated by the environmenlal contour (EC) method and compared with those estimated by a statistical extrapolation method. It is demonstrated that the extreme load effects estimated by the EC method together with the simplified algorithm is highly accurate.

    extreme load effect; limit state line; joint action of wind and wave; environmental contour method; copula function

    TU411; TU472.5

    A

    10.13465/j.cnki.jvs.2017.19.019

    交通運(yùn)輸部科技項(xiàng)目(2011318223170);國(guó)家自然科學(xué)基金資助項(xiàng)目(51578504)

    2016-04-19 修改稿收到日期:2016-08-04

    涂志斌 女,博士生,1988年9月生

    黃銘楓 男,博士,副教授,博士生導(dǎo)師,1976年10月生

    猜你喜歡
    有限元效應(yīng)
    鈾對(duì)大型溞的急性毒性效應(yīng)
    懶馬效應(yīng)
    場(chǎng)景效應(yīng)
    新型有機(jī)玻璃在站臺(tái)門(mén)的應(yīng)用及有限元分析
    基于有限元的深孔鏜削仿真及分析
    基于有限元模型對(duì)踝模擬扭傷機(jī)制的探討
    應(yīng)變效應(yīng)及其應(yīng)用
    磨削淬硬殘余應(yīng)力的有限元分析
    偶像效應(yīng)
    基于SolidWorks的吸嘴支撐臂有限元分析
    国产男靠女视频免费网站| 免费黄网站久久成人精品| 国产成人福利小说| 亚洲av日韩精品久久久久久密| 变态另类成人亚洲欧美熟女| 听说在线观看完整版免费高清| 精品一区二区免费观看| 最好的美女福利视频网| 成人av一区二区三区在线看| 免费搜索国产男女视频| 国产激情偷乱视频一区二区| 老熟妇仑乱视频hdxx| 久久久久久九九精品二区国产| 久久久久久九九精品二区国产| 成人国产一区最新在线观看| 欧美一级a爱片免费观看看| 午夜福利欧美成人| 免费人成在线观看视频色| 啪啪无遮挡十八禁网站| 日本-黄色视频高清免费观看| 嫁个100分男人电影在线观看| 免费av毛片视频| 午夜日韩欧美国产| 国产精品免费一区二区三区在线| 丝袜美腿在线中文| 国产精品乱码一区二三区的特点| 在线观看美女被高潮喷水网站| 高清在线国产一区| 久久久久久久久大av| 国产精品电影一区二区三区| 日韩高清综合在线| 最新中文字幕久久久久| 亚洲成人久久爱视频| 自拍偷自拍亚洲精品老妇| 成人鲁丝片一二三区免费| 三级男女做爰猛烈吃奶摸视频| 日本一二三区视频观看| 99久久精品一区二区三区| 午夜免费激情av| 蜜桃亚洲精品一区二区三区| 精品人妻1区二区| 亚洲av电影不卡..在线观看| 麻豆成人av在线观看| 欧美色视频一区免费| 一级黄色大片毛片| 成年人黄色毛片网站| 久久99热这里只有精品18| 免费大片18禁| 亚洲精品在线观看二区| 少妇的逼水好多| 深爱激情五月婷婷| 亚洲精品久久国产高清桃花| 自拍偷自拍亚洲精品老妇| 搡女人真爽免费视频火全软件 | 男女下面进入的视频免费午夜| 熟女人妻精品中文字幕| 婷婷精品国产亚洲av在线| 日本一本二区三区精品| 美女高潮喷水抽搐中文字幕| 成年人黄色毛片网站| 国产精品野战在线观看| 97超视频在线观看视频| 久久亚洲真实| 美女 人体艺术 gogo| 精品一区二区免费观看| 国产在线男女| 国产美女午夜福利| 少妇的逼好多水| 亚洲欧美日韩东京热| 亚洲av第一区精品v没综合| 小说图片视频综合网站| 免费在线观看成人毛片| 免费观看在线日韩| 日韩强制内射视频| 久久久久免费精品人妻一区二区| 91久久精品国产一区二区成人| 国产免费一级a男人的天堂| 1024手机看黄色片| 亚洲av成人av| 国产视频一区二区在线看| 亚洲男人的天堂狠狠| 嫩草影院新地址| 91久久精品电影网| 免费在线观看成人毛片| 俄罗斯特黄特色一大片| 国产探花在线观看一区二区| 露出奶头的视频| 伦理电影大哥的女人| 香蕉av资源在线| av女优亚洲男人天堂| 国产男靠女视频免费网站| 欧美最黄视频在线播放免费| 又爽又黄无遮挡网站| 色视频www国产| 免费人成视频x8x8入口观看| 国产91精品成人一区二区三区| 久久精品人妻少妇| 成年版毛片免费区| 高清在线国产一区| 欧美日韩中文字幕国产精品一区二区三区| 国产午夜福利久久久久久| 国产激情偷乱视频一区二区| 亚洲一区高清亚洲精品| 午夜老司机福利剧场| 亚洲自拍偷在线| 亚洲国产日韩欧美精品在线观看| 老司机午夜福利在线观看视频| 一本久久中文字幕| 麻豆国产97在线/欧美| 国产伦精品一区二区三区视频9| 美女黄网站色视频| 男人舔奶头视频| 国产伦在线观看视频一区| 免费看日本二区| 日韩国内少妇激情av| 亚洲国产日韩欧美精品在线观看| 国产美女午夜福利| 91久久精品国产一区二区三区| 少妇丰满av| 三级国产精品欧美在线观看| 国产精品国产高清国产av| 亚洲中文日韩欧美视频| 一本精品99久久精品77| 一级av片app| 午夜福利在线观看吧| 国产成人福利小说| 日韩欧美国产一区二区入口| 日韩在线高清观看一区二区三区 | 成人av在线播放网站| 亚洲av日韩精品久久久久久密| 国产伦精品一区二区三区视频9| 国产精品福利在线免费观看| 三级国产精品欧美在线观看| 一进一出抽搐gif免费好疼| 变态另类成人亚洲欧美熟女| 亚洲性久久影院| 国内精品一区二区在线观看| 白带黄色成豆腐渣| 亚洲天堂国产精品一区在线| 亚洲av电影不卡..在线观看| 亚洲精品456在线播放app | 成人av一区二区三区在线看| 国产熟女欧美一区二区| 亚洲,欧美,日韩| 在线a可以看的网站| 黄色日韩在线| 又黄又爽又刺激的免费视频.| 亚洲人成伊人成综合网2020| 无人区码免费观看不卡| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| 在线国产一区二区在线| 国产av麻豆久久久久久久| 国内精品宾馆在线| 国产精品女同一区二区软件 | 成人午夜高清在线视频| a级毛片a级免费在线| 久久久成人免费电影| 国产高清视频在线播放一区| 在线播放无遮挡| 蜜桃久久精品国产亚洲av| 日韩亚洲欧美综合| 精品乱码久久久久久99久播| 日本黄大片高清| 国产午夜精品论理片| 欧美国产日韩亚洲一区| 亚洲人与动物交配视频| 春色校园在线视频观看| 深爱激情五月婷婷| 国产亚洲91精品色在线| 免费av毛片视频| 大型黄色视频在线免费观看| 欧美精品啪啪一区二区三区| 午夜视频国产福利| 成人美女网站在线观看视频| 天堂影院成人在线观看| 亚洲第一电影网av| 精品午夜福利视频在线观看一区| 别揉我奶头~嗯~啊~动态视频| 久久久久久久精品吃奶| 亚洲狠狠婷婷综合久久图片| 欧美一区二区亚洲| 国产精品日韩av在线免费观看| 国产精品一及| 日本色播在线视频| 男女啪啪激烈高潮av片| h日本视频在线播放| 久久精品国产亚洲av涩爱 | 少妇的逼水好多| 精品一区二区三区av网在线观看| 亚洲,欧美,日韩| 精品不卡国产一区二区三区| 国产精品亚洲一级av第二区| 人妻夜夜爽99麻豆av| 国产精品一区二区性色av| 欧美最新免费一区二区三区| 成年女人看的毛片在线观看| 国产伦一二天堂av在线观看| 在线天堂最新版资源| 亚洲国产欧洲综合997久久,| 国产aⅴ精品一区二区三区波| 国产精品野战在线观看| 99视频精品全部免费 在线| 人妻夜夜爽99麻豆av| 亚洲成人精品中文字幕电影| 综合色av麻豆| 亚洲精华国产精华精| 无人区码免费观看不卡| xxxwww97欧美| 欧美另类亚洲清纯唯美| 极品教师在线视频| 1024手机看黄色片| 日韩欧美精品免费久久| 亚洲黑人精品在线| 午夜免费激情av| aaaaa片日本免费| 给我免费播放毛片高清在线观看| 国产成人av教育| 一级黄片播放器| 亚洲欧美日韩高清专用| 日本一二三区视频观看| 欧美成人一区二区免费高清观看| 综合色av麻豆| 国产精品国产三级国产av玫瑰| 久久中文看片网| 中文字幕熟女人妻在线| 中国美女看黄片| 精品人妻熟女av久视频| 午夜精品久久久久久毛片777| 国产成人影院久久av| 欧美色视频一区免费| 乱码一卡2卡4卡精品| 欧美潮喷喷水| 亚洲中文日韩欧美视频| netflix在线观看网站| 在线观看舔阴道视频| 免费观看在线日韩| 搡女人真爽免费视频火全软件 | 日日摸夜夜添夜夜添av毛片 | 欧美黑人巨大hd| xxxwww97欧美| 国产在视频线在精品| 小蜜桃在线观看免费完整版高清| 又黄又爽又免费观看的视频| 国产av一区在线观看免费| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 欧美日本视频| 国产探花在线观看一区二区| 国产亚洲精品av在线| 亚洲三级黄色毛片| 亚洲,欧美,日韩| 变态另类成人亚洲欧美熟女| 午夜影院日韩av| 国产伦精品一区二区三区视频9| 免费av毛片视频| 美女黄网站色视频| 日韩欧美精品免费久久| 久久精品国产99精品国产亚洲性色| 美女 人体艺术 gogo| 国产麻豆成人av免费视频| 国产欧美日韩一区二区精品| 国内精品宾馆在线| 男人的好看免费观看在线视频| 欧美xxxx黑人xx丫x性爽| 一进一出抽搐动态| 精品99又大又爽又粗少妇毛片 | 欧美三级亚洲精品| 日日撸夜夜添| 亚洲欧美日韩卡通动漫| 99热6这里只有精品| 欧美另类亚洲清纯唯美| 精品国内亚洲2022精品成人| 国产精品99久久久久久久久| 婷婷亚洲欧美| 国产一级毛片七仙女欲春2| 久久精品91蜜桃| 成年女人看的毛片在线观看| 成人永久免费在线观看视频| 久久精品国产亚洲av涩爱 | 久久热精品热| 又粗又爽又猛毛片免费看| 日韩人妻高清精品专区| 麻豆久久精品国产亚洲av| 国产精品98久久久久久宅男小说| av.在线天堂| 欧美精品啪啪一区二区三区| 成人永久免费在线观看视频| 国产毛片a区久久久久| 97碰自拍视频| 麻豆国产97在线/欧美| 午夜爱爱视频在线播放| 97热精品久久久久久| 久99久视频精品免费| 99久久精品一区二区三区| 精品一区二区三区人妻视频| 久久精品国产亚洲av香蕉五月| 免费观看人在逋| 熟女人妻精品中文字幕| 乱码一卡2卡4卡精品| 看黄色毛片网站| 在线a可以看的网站| 少妇丰满av| 国产一区二区三区在线臀色熟女| 女人十人毛片免费观看3o分钟| 亚洲成人久久爱视频| 日韩欧美国产在线观看| 婷婷亚洲欧美| 精品久久久久久,| or卡值多少钱| 中亚洲国语对白在线视频| 午夜爱爱视频在线播放| 日韩欧美国产在线观看| 亚洲成人免费电影在线观看| 亚洲国产精品sss在线观看| 亚洲av日韩精品久久久久久密| 久9热在线精品视频| 国产精品永久免费网站| 国内精品一区二区在线观看| 99国产极品粉嫩在线观看| 色综合站精品国产| 午夜爱爱视频在线播放| 黄色丝袜av网址大全| 欧美精品啪啪一区二区三区| 亚洲精品日韩av片在线观看| 舔av片在线| 久久99热6这里只有精品| 午夜视频国产福利| 久久久久久九九精品二区国产| 国产国拍精品亚洲av在线观看| 色播亚洲综合网| 别揉我奶头 嗯啊视频| 在线观看舔阴道视频| av在线蜜桃| 男女啪啪激烈高潮av片| 久久久久久大精品| 免费看美女性在线毛片视频| 此物有八面人人有两片| 男人舔女人下体高潮全视频| 天堂动漫精品| 成年女人永久免费观看视频| 国产人妻一区二区三区在| 久久久久精品国产欧美久久久| 日本成人三级电影网站| 日韩欧美三级三区| 日本爱情动作片www.在线观看 | 少妇被粗大猛烈的视频| 长腿黑丝高跟| 亚洲七黄色美女视频| 最新在线观看一区二区三区| 日韩欧美国产在线观看| 国产色婷婷99| 亚洲,欧美,日韩| 国产精品亚洲一级av第二区| 亚洲av二区三区四区| 男人舔奶头视频| 日韩欧美 国产精品| 国产高潮美女av| 久久这里只有精品中国| 淫秽高清视频在线观看| 国产69精品久久久久777片| 国产 一区 欧美 日韩| h日本视频在线播放| 日本黄大片高清| 欧美日韩精品成人综合77777| 色综合亚洲欧美另类图片| 男人舔女人下体高潮全视频| 日韩欧美精品免费久久| 国产麻豆成人av免费视频| 男人和女人高潮做爰伦理| 黄色丝袜av网址大全| 亚洲午夜理论影院| 嫩草影院精品99| 午夜精品一区二区三区免费看| 乱系列少妇在线播放| 亚洲精品久久国产高清桃花| 蜜桃久久精品国产亚洲av| 99视频精品全部免费 在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精华国产精华精| 18禁黄网站禁片免费观看直播| 亚洲精品亚洲一区二区| 成人永久免费在线观看视频| 久久天躁狠狠躁夜夜2o2o| 国产国拍精品亚洲av在线观看| 欧美高清成人免费视频www| 91精品国产九色| 成年女人毛片免费观看观看9| 午夜精品久久久久久毛片777| 人妻丰满熟妇av一区二区三区| 日韩精品中文字幕看吧| 久久久精品欧美日韩精品| 啦啦啦啦在线视频资源| 天堂网av新在线| 亚洲美女视频黄频| 一边摸一边抽搐一进一小说| 全区人妻精品视频| 性插视频无遮挡在线免费观看| 中文在线观看免费www的网站| 亚洲成av人片在线播放无| 国产不卡一卡二| 欧美成人a在线观看| 又黄又爽又免费观看的视频| 日本五十路高清| 国产在线男女| 免费在线观看成人毛片| 乱系列少妇在线播放| 老师上课跳d突然被开到最大视频| 亚洲av五月六月丁香网| 精品国产三级普通话版| 色精品久久人妻99蜜桃| 老熟妇乱子伦视频在线观看| 嫩草影院精品99| 午夜精品在线福利| 小说图片视频综合网站| 国产国拍精品亚洲av在线观看| 亚洲图色成人| 日本五十路高清| 一区二区三区免费毛片| 日韩精品中文字幕看吧| 国产一区二区激情短视频| 中文字幕精品亚洲无线码一区| 精品一区二区三区人妻视频| 欧美成人性av电影在线观看| 国产免费男女视频| 九九热线精品视视频播放| 亚洲成av人片在线播放无| 黄色视频,在线免费观看| 亚洲无线观看免费| 99国产精品一区二区蜜桃av| 成年版毛片免费区| 日韩高清综合在线| 日本黄大片高清| 久99久视频精品免费| 波多野结衣高清无吗| 乱系列少妇在线播放| av天堂中文字幕网| 亚洲精品久久国产高清桃花| 91狼人影院| 成人国产麻豆网| 老司机福利观看| 日本-黄色视频高清免费观看| 欧美在线一区亚洲| 国产亚洲精品av在线| 亚洲成人中文字幕在线播放| 中文字幕免费在线视频6| 美女大奶头视频| videossex国产| av福利片在线观看| 俄罗斯特黄特色一大片| 亚洲,欧美,日韩| 色噜噜av男人的天堂激情| 国产精品嫩草影院av在线观看 | 三级国产精品欧美在线观看| 亚洲欧美清纯卡通| 91狼人影院| 嫩草影院新地址| 黄色欧美视频在线观看| 精品一区二区三区av网在线观看| 日本免费一区二区三区高清不卡| 亚洲国产精品成人综合色| 日韩欧美国产在线观看| 国产精品av视频在线免费观看| 在线免费观看不下载黄p国产 | 欧美一区二区亚洲| 黄色欧美视频在线观看| 免费av毛片视频| 亚洲在线观看片| 琪琪午夜伦伦电影理论片6080| 日韩在线高清观看一区二区三区 | 午夜亚洲福利在线播放| 欧美激情久久久久久爽电影| 色哟哟·www| 人人妻人人看人人澡| 日日撸夜夜添| 99热精品在线国产| 久久这里只有精品中国| 精品一区二区三区av网在线观看| 18禁黄网站禁片免费观看直播| 中文字幕免费在线视频6| 亚洲精品亚洲一区二区| 夜夜夜夜夜久久久久| 日本黄色视频三级网站网址| 在现免费观看毛片| av福利片在线观看| 欧美性猛交╳xxx乱大交人| 又黄又爽又免费观看的视频| 亚洲精品国产成人久久av| 搞女人的毛片| 久久精品综合一区二区三区| 国产精品不卡视频一区二区| 国产精品,欧美在线| 日本一本二区三区精品| 精品乱码久久久久久99久播| 精品99又大又爽又粗少妇毛片 | 国产黄色小视频在线观看| 亚洲国产精品合色在线| 一边摸一边抽搐一进一小说| 欧美一区二区精品小视频在线| 日韩欧美精品免费久久| 欧美一区二区精品小视频在线| 成人无遮挡网站| 悠悠久久av| 夜夜爽天天搞| 国产熟女欧美一区二区| 欧美日韩乱码在线| 日本成人三级电影网站| 国产精品女同一区二区软件 | 18禁在线播放成人免费| 久久精品国产亚洲av涩爱 | 麻豆成人av在线观看| 国内少妇人妻偷人精品xxx网站| 色噜噜av男人的天堂激情| 欧美日韩中文字幕国产精品一区二区三区| 露出奶头的视频| 日韩欧美精品免费久久| 午夜久久久久精精品| 欧美激情久久久久久爽电影| 免费看a级黄色片| 成年女人永久免费观看视频| 午夜福利视频1000在线观看| 日本三级黄在线观看| 欧美绝顶高潮抽搐喷水| av在线观看视频网站免费| 色综合站精品国产| 欧美高清性xxxxhd video| 国产精品亚洲美女久久久| 男女做爰动态图高潮gif福利片| 成人av在线播放网站| 婷婷精品国产亚洲av在线| 真人一进一出gif抽搐免费| 午夜福利高清视频| 最后的刺客免费高清国语| 一区二区三区高清视频在线| 久久天躁狠狠躁夜夜2o2o| 国产精品综合久久久久久久免费| 一边摸一边抽搐一进一小说| 日本五十路高清| 午夜爱爱视频在线播放| 精品国产三级普通话版| 成年免费大片在线观看| 高清日韩中文字幕在线| 亚洲国产精品久久男人天堂| 国产成人福利小说| 国产 一区 欧美 日韩| 国产黄a三级三级三级人| 日本一本二区三区精品| 舔av片在线| 久久人妻av系列| 日韩强制内射视频| 大又大粗又爽又黄少妇毛片口| 欧美日韩中文字幕国产精品一区二区三区| 国内精品久久久久精免费| or卡值多少钱| 少妇的逼水好多| 琪琪午夜伦伦电影理论片6080| 不卡视频在线观看欧美| 亚洲七黄色美女视频| 亚洲电影在线观看av| 高清毛片免费观看视频网站| 午夜免费成人在线视频| 51国产日韩欧美| 日韩欧美精品v在线| 欧美日韩亚洲国产一区二区在线观看| 黄片wwwwww| 美女被艹到高潮喷水动态| 国产高清视频在线观看网站| 欧美日韩瑟瑟在线播放| 午夜亚洲福利在线播放| 成人特级黄色片久久久久久久| 欧美成人a在线观看| 91午夜精品亚洲一区二区三区 | 欧美一区二区精品小视频在线| 麻豆国产av国片精品| 日韩中字成人| 又爽又黄无遮挡网站| 欧美精品国产亚洲| 亚洲人成网站高清观看| 一区二区三区四区激情视频 | 亚洲avbb在线观看| 日本与韩国留学比较| 国产高清三级在线| 免费av毛片视频| 国产亚洲精品综合一区在线观看| 日韩强制内射视频| 成年女人看的毛片在线观看| 成人无遮挡网站| 亚洲国产精品sss在线观看| 亚洲不卡免费看| 丰满的人妻完整版| 久久久久久久久中文| 久久久久久大精品| 女同久久另类99精品国产91| 欧美在线一区亚洲| av中文乱码字幕在线| 三级毛片av免费| 国产精品,欧美在线| 日本一本二区三区精品| 亚洲va在线va天堂va国产| 国产三级在线视频| 中文在线观看免费www的网站| 日日撸夜夜添| 成人高潮视频无遮挡免费网站| 色哟哟·www| 免费av不卡在线播放| 欧美色欧美亚洲另类二区| 久久欧美精品欧美久久欧美| 淫妇啪啪啪对白视频| 国产成年人精品一区二区| 色5月婷婷丁香| 国产精品女同一区二区软件 | 久久久久久久午夜电影| 波多野结衣巨乳人妻| 乱人视频在线观看| 国产一区二区三区视频了|