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

    考慮矩震級(jí)影響的位移譜阻尼調(diào)整系數(shù)

    2014-06-07 07:15:48王國(guó)弢胡克旭
    結(jié)構(gòu)工程師 2014年5期
    關(guān)鍵詞:影響模型

    王國(guó)弢 胡克旭 雷 敏

    (同濟(jì)大學(xué)結(jié)構(gòu)工程與防災(zāi)研究所,上海200092)

    考慮矩震級(jí)影響的位移譜阻尼調(diào)整系數(shù)

    王國(guó)弢*胡克旭 雷 敏

    (同濟(jì)大學(xué)結(jié)構(gòu)工程與防災(zāi)研究所,上海200092)

    基于408條地震地面運(yùn)動(dòng)記錄統(tǒng)計(jì)分析,首先研究了矩震級(jí)(Mw)對(duì)位移譜阻尼調(diào)整系數(shù)(DSF)的影響。結(jié)果表明:Mw對(duì)DSF的影響與周期(T)和阻尼比(ξ)有關(guān);在T>0.6 s的范圍內(nèi),當(dāng)ξ<5%時(shí),DSF隨著Mw的增長(zhǎng)而增長(zhǎng);當(dāng)ξ>5%時(shí),DSF隨著Mw的增長(zhǎng)而減小。隨著T的增長(zhǎng),Mw對(duì)DSF的影響越顯著;ξ越遠(yuǎn)離5%,Mw對(duì)DSF的影響越顯著。在T<0.6 s的范圍內(nèi),Mw對(duì)DSF無(wú)顯著影響;當(dāng)Mw=5~6時(shí),DSF受T的影響顯著;隨著Mw的增大,在T>~0.6 s范圍內(nèi),DSF受T的影響逐漸減弱;當(dāng)Mw=7~8時(shí),在T>0.6 s范圍內(nèi),除阻尼比為0.5%和1%的DSF外,DSF基本不隨T而變化;DSF相對(duì)于ξ存在異方差性。然后提出了只考慮周期和阻尼比影響的DSF回歸方程,計(jì)算出在指定周期和阻尼比處的殘差,根據(jù)在長(zhǎng)周期處該回歸方程的殘差隨Mw顯直線分布,最后提出了包含線性矩震級(jí)項(xiàng)的DSF回歸模型并給出了模型的回歸系數(shù)。該模型的殘差隨Mw和相對(duì)能量持時(shí)(D5-95)的分布表明:模型能體現(xiàn)Mw對(duì)DSF的影響,能間接地體現(xiàn)D5-95對(duì)DSF的影響。模型的標(biāo)準(zhǔn)差能體現(xiàn)DSF相對(duì)于ξ的異方差性。

    矩震級(jí),阻尼調(diào)整系數(shù),周期,阻尼比,重要持時(shí)

    Keywordsmomentmagnitude,damping scaling factor,period,damping ratio,significant duration

    1 引 言

    在結(jié)構(gòu)的抗震設(shè)計(jì)中需要不同阻尼比下的彈性反應(yīng)譜。一般可采用兩種方法來(lái)獲得不同阻尼比下的彈性反應(yīng)譜:①直接形成不同阻尼比下的反應(yīng)譜預(yù)測(cè)方程;②采用阻尼調(diào)整系數(shù)將5%阻尼比的彈性反應(yīng)譜調(diào)整到不同于5%阻尼比的彈性反應(yīng)譜。各國(guó)抗震規(guī)范中普遍采用的是第二種方法。對(duì)于不同的反應(yīng)譜,有不同的阻尼調(diào)整系數(shù),本文所研究的阻尼調(diào)整系數(shù)為位移譜阻尼調(diào)整系數(shù)。近年來(lái),對(duì)阻尼調(diào)整系數(shù)的研究表明:矩震級(jí)對(duì)阻尼調(diào)整系數(shù)有顯著的影響;但中國(guó)和歐美現(xiàn)行抗震規(guī)范[1-6]中提供的阻尼調(diào)整系數(shù)只與阻尼比和結(jié)構(gòu)自振周期有關(guān)或僅與阻尼比有關(guān),不能反映矩震級(jí)對(duì)阻尼調(diào)整系數(shù)的影響。到目前為止,已有許多學(xué)者對(duì)阻尼調(diào)整系數(shù)進(jìn)行了研究。Ashour[7]、Tolis和Faccioli[8]及Priestley[9]提出了僅考慮阻尼比影響的阻尼調(diào)整系數(shù)回歸模型;Newmark和Hall[10]、Wu和Hanson[11]及Idriss[12]提出了考慮阻尼比和周期影響的阻尼調(diào)整系數(shù)回歸模型。Lin和Chang[13]研究了場(chǎng)地條件對(duì)阻尼調(diào)整系數(shù)的影響,并提出了考慮阻尼比、周期和場(chǎng)地條件影響的阻尼調(diào)整系數(shù)回歸模型。以上研究者所提出的回歸模型均沒(méi)有考慮矩震級(jí)對(duì)阻尼調(diào)整系數(shù)的影響。Bommer和Mendis[14]研究了矩震級(jí)對(duì)阻尼調(diào)整系數(shù)的影響,對(duì)阻尼比大于5%的情況,提出了阻尼調(diào)整系數(shù)隨矩震級(jí)和觀測(cè)點(diǎn)到斷層面最近距離的增大而減小的結(jié)論;Cameron和Green[15]認(rèn)為阻尼調(diào)整系數(shù)與地震動(dòng)的頻譜特性和持時(shí)有關(guān)并且他們將地面運(yùn)動(dòng)記錄分別按地質(zhì)構(gòu)造環(huán)境、場(chǎng)地條件和震級(jí)進(jìn)行分組,以表格的形式列出了在指定阻尼比和周期處的阻尼調(diào)整系數(shù)的中值和標(biāo)準(zhǔn)差,但他們均沒(méi)有提出一個(gè)簡(jiǎn)單統(tǒng)一的能反映矩震級(jí)影響的回歸模型。本文的研究目的是提出一個(gè)簡(jiǎn)單統(tǒng)一的阻尼調(diào)整系數(shù)回歸模型并在模型中包含矩震級(jí)這個(gè)變量,使該模型能體現(xiàn)矩震級(jí)對(duì)阻尼調(diào)整系數(shù)的影響。此外,根據(jù)Kempton和Stewart[16],持時(shí)與矩震級(jí)有正的強(qiáng)相關(guān)性,所以包含了矩震級(jí)的回歸模型還可以間接體現(xiàn)持時(shí)對(duì)阻尼調(diào)整系數(shù)的影響。本文所采用的持時(shí)為相對(duì)能量持時(shí)D[17]5-95。

    2 地面運(yùn)動(dòng)加速度時(shí)程記錄

    在本文的研究中使用了408條地震地面運(yùn)動(dòng)記錄,這些地面運(yùn)動(dòng)記錄下載于太平洋地震工程研究中心的強(qiáng)震數(shù)據(jù)庫(kù)(http://peer.berkeley.edu/peer ground motion database)。地面運(yùn)動(dòng)記錄的相對(duì)能量持時(shí)(D5-95)分布在1.4 s~80.3 s之間,矩震級(jí)(Mw)分布在5~7.62之間,觀測(cè)點(diǎn)到斷層破裂面的最近距離(Rrup)分布在0.1 km~100 km之間,地表厚度30 m內(nèi)平均剪切波速(Vs30)分布在116.3 m/s~2 016.1 m/s之間,峰值地面加速度(PGA)分布在0.004 7 g~1.434 5 g之間。圖1表明了本文所用記錄的Mw和Rrup的分布。

    圖1 Mw—Rrup分布Fig.1 Magnitude Vs.distance distribution

    3 矩震級(jí)對(duì)阻尼調(diào)整系數(shù)的影響

    式中,DSF代表位移譜阻尼調(diào)整系數(shù);ξ代表阻尼比;T代表單自由度體系的自振周期;SD(ξ,T)和PSA(ξ,T)分別是ξ≠5%,周期為T時(shí)彈性單自由度體系的位移譜值和擬加速度譜值;SD(5%,T)和PSA(5%,T)分別是ξ=5%,周期為T時(shí)彈性單自由度體系的位移譜值和擬加速度譜值。

    本文采用軟件SeismoSignal對(duì)每條輸入地震波分別計(jì)算了11個(gè)阻尼比下位移譜。11個(gè)阻尼比分別為0.5%、1%、2%、3%、5%、7%、10%、15%、20%、25%和30%。對(duì)每個(gè)阻尼比,單自由度體系取如下20個(gè)周期:0.02 s、0.04 s、0.06 s、0.08 s、0.1 s、0.14 s、0.2 s、0.24 s、0.3 s、0.4 s、0.5 s、0.6 s、0.8 s、1 s、2 s、3 s、4 s、5 s、7.5 s和10 s。然后再根據(jù)式(1),對(duì)每條輸入地震波可求

    位移譜阻尼調(diào)整系數(shù)被定義為得各阻尼比下的位移譜阻尼調(diào)整系數(shù)。

    為了研究Mw對(duì)DSF的影響,本文將地面運(yùn)動(dòng)記錄按Mw分為三組,分組情況見(jiàn)表1。

    表1 地面運(yùn)動(dòng)記錄分組Table 1 Groups of ground motion records

    圖2表明了ξ=0.5%、2%、7%和25%時(shí),Mw對(duì)DSF中值的影響;圖3表明了ξ=1%、3%、10%和30%時(shí),Mw對(duì)DSF中值的影響。由圖2和圖3可見(jiàn),Mw對(duì)DSF的影響可歸納為以下幾點(diǎn):

    (1)Mw對(duì)DSF的影響與T和ξ有關(guān);在T>0.6 s的范圍內(nèi),當(dāng)ξ<5%時(shí),DSF隨著Mw的增長(zhǎng)而增長(zhǎng);當(dāng)ξ>5%時(shí),DSF隨著Mw的增長(zhǎng)而減小.

    (2)隨著T的增長(zhǎng),Mw對(duì)DSF的影響越顯著;ξ越遠(yuǎn)離5%,Mw對(duì)DSF的影響越顯著。在T<0.6 s的范圍內(nèi),Mw對(duì)DSF無(wú)顯著影響.

    (3)當(dāng)Mw=5~6時(shí),DSF受T的影響顯著;隨著Mw的增大,在T>0.6 s范圍內(nèi),DSF受T的影響逐漸減弱;當(dāng)Mw=7~8時(shí),除阻尼比為0.5%和1%的DSF外,在T>0.6 s范圍內(nèi),DSF基本不隨T而變化。

    圖2 Mw對(duì)DSF中值的影響(ξ=0.5%、2%、7%和25%)Fig.2 Influence of Mwon themedian DSF(ξ=0.5%、2%、7%and 25%)

    圖3 Mw對(duì)DSF中值的影響(ξ=1%、3%、10%和30%)Fig.3 Influence of Mwon themedian DSF(ξ=1%、3%、10%和30%)

    圖4 給出了T=1 s時(shí),Mw與DSF之間的統(tǒng)計(jì)關(guān)系,圖中直線為Mw和DSF的線性擬合線。由圖4也可見(jiàn):Mw對(duì)DSF的影響與ξ有關(guān)。此外,圖4還表明了DSF相對(duì)于ξ的異方差性。隨著ξ越遠(yuǎn)離5%,DSF的離散程度越大。DSF相對(duì)于ξ的異方差性將在后面的方差模型中考慮。

    圖4 Mw對(duì)DSF的影響Fig.4 Influence of Mwon DSF

    綜上所述,Mw對(duì)DSF有顯著的影響,所以在DSF的回歸模型中應(yīng)考慮Mw的影響。

    4 DSF的概率分布

    通??梢哉J(rèn)為PSA服從對(duì)數(shù)正態(tài)分布,對(duì)式(1)兩邊取自然對(duì)數(shù)有:

    由式(2)可知,如果PSA(ξ,T)與PSA(5%,T)不相關(guān),則ln(DSF(ξ,T))服從正態(tài)分布,DSF服從對(duì)數(shù)正態(tài)分布。但是,PSA(ξ,T)與PSA(5%,T)是相關(guān)的,所以在理論上DSF并不嚴(yán)格服從對(duì)數(shù)正態(tài)分布。根據(jù)對(duì)本文數(shù)據(jù)的分析,在T∈[0.1,7.5]的范圍內(nèi),ln(DSF(ξ,T))與正態(tài)分布曲線擬合得較好。圖5給出了ξ=2%時(shí),在T=0.2 s和1 s處,ln(DSF)數(shù)據(jù)點(diǎn)的直方圖、累積分布函數(shù)曲線圖及相應(yīng)的正態(tài)分布函數(shù)的擬合結(jié)果,圖中CDF表示累積分布函數(shù)(cumulative distribution function),PDF表示概率密度函數(shù)(probability density function)。在此周期范圍外,擬合結(jié)果并不理想。

    綜上,本文假設(shè)在指定的周期和阻尼比處,DSF近似服從對(duì)數(shù)正態(tài)分布,并取回歸模型為

    ln(DSF)=μ(ξ,T,Mw)+?(3)式中,μ(ξ,T,Mw)為ln(DSF)的期望;?代表隨機(jī)誤差且?服從正態(tài)分布N(0,σ2);σ2為方差。

    圖5 ln(DSF)的分布和擬合曲線Fig.5 Distribution of ln(DSF)and fitted curve

    5 建立回歸模型

    5.1 回歸方程μ(ξ,T,Mw)

    首先設(shè)回歸方程僅是ξ和T的函數(shù),取回歸模型為

    在指定周期Ti和阻尼比ξj處,μ(ξj,Ti)的估計(jì)值為

    其中,ln[DSF(ξj,Ti)k]為在Ti和ξj處ln(DSF)的第k個(gè)數(shù)據(jù)值,n=408為數(shù)據(jù)總數(shù);在Ti和ξj處,第k個(gè)數(shù)據(jù)與回歸方程的殘差[Δ(ξj,Ti)]k為

    按式(6)可求出在ξj和Ti處每個(gè)數(shù)據(jù)與回歸方程的殘差,然后可作出殘差相對(duì)于Mw和D5-95的分布。圖6和圖7給出了在T=0.1 s和3 s處,ξ=2%和20%時(shí),殘差相對(duì)于Mw和D5-95的分布。分別將計(jì)算出的殘差按Mw和ln(D5-95)各分為9組。按Mw分組的情況為:5≤Mw<5.3、5.3≤Mw<5.6、5.6≤Mw<5.9、5.9≤Mw<6.2、6.2≤Mw<6.5、6.5≤Mw<6.8、6.8≤Mw<7.1、7.1≤Mw<7.4和Mw≥7.4。

    按ln(D5-95)分組的情況為ln(D5-95)<1.7、1.7<ln(D5-95)<2、2<ln(D5-95)<2.3、2.3<ln(D5-95)<2.6、2.6<ln(D5-95)<2.9、2.9<ln(D5-95)<3.2、3.2<ln(D5-95)<3.5、3.5<ln(D5-95)<3.8和ln(D5-95)>3.8。

    圖中黑色方框代表每組殘差的均值,黑線為均值的連線,表明了殘差隨Mw和D5-95的分布。根據(jù)殘差分析可知:在短周期處(如T=0.1 s),各阻尼比下的殘差基本隨機(jī)對(duì)稱的分布于零水平線的兩側(cè)(見(jiàn)圖6和圖7),這表明了在短周期處Mw及D5-95對(duì)DSF沒(méi)有顯著的影響;在長(zhǎng)周期處(如:T=3 s),各阻尼比下的殘差相對(duì)于Mw近似顯直線分布,相對(duì)于D5-95顯曲線分布,這表明了在長(zhǎng)周期處Mw及D5-95對(duì)DSF有顯著的影響。以上通過(guò)殘差分析所反映的Mw對(duì)DSF的影響是與前文的分析結(jié)果是相一致的。

    圖6 殘差相對(duì)于Mw的分布Fig.6 Residuals vs.Mw

    圖7 殘差相對(duì)于D5-95的分布Fig.7 Residuals vs.D5-95

    為了改善長(zhǎng)周期處殘差相對(duì)于Mw的分布,在指定周期點(diǎn)Ti處,本文將回歸方程取為Mw的線性函數(shù):

    在指定周期點(diǎn)Ti處,回歸過(guò)程分為以下兩步。

    第一步:在Ti處,用最小二乘法回歸出式(7)中在指定阻尼比ξj處,主系數(shù)的估計(jì)值(^c0(ξj,Ti)和^c1(ξj,Ti))。

    對(duì)11個(gè)阻尼比重復(fù)上述過(guò)程,可獲得所有阻尼比處主系數(shù)的估計(jì)值。

    第二步:在Ti處,采用y=a+b lnξ+c(lnξ)2的函數(shù)形式來(lái)擬合主系數(shù)的估計(jì)值與ξ之間的關(guān)系(ξ取為百分號(hào)中的數(shù)值,例如,當(dāng)阻尼比為5%時(shí),取ξ=5,以下各方程中均類同)。設(shè)

    以11個(gè)阻尼比處主系數(shù)的估計(jì)值為數(shù)據(jù)點(diǎn),分別對(duì)式(8)和式(9)執(zhí)行回歸分析,可求得式(8)和式(9)中次系數(shù)的估計(jì)值(如^b0(Ti)、^b1(Ti)和^b2(Ti)等)。圖8給出了在T=0.1 s和3 s處,式(8)和式(9)的擬合結(jié)果。對(duì)20個(gè)周期點(diǎn)重復(fù)上述回歸過(guò)程,可獲得所有周期點(diǎn)處的回歸方程。在Ti處的回歸方程如下式:

    表2給出了20個(gè)周期點(diǎn)處回歸系數(shù)的估計(jì)值。其余周期點(diǎn)處回歸系數(shù)的估計(jì)值可由線性插值確定。在Ti和ξj處,第k個(gè)數(shù)據(jù)與回歸方程的殘差為:

    按式(11)求出在ξj和Ti處每個(gè)數(shù)據(jù)與回歸方程的殘差,然后可作出殘差相對(duì)于Mw和D5-95的分布。圖9和圖10給出了在T=0.1 s和3 s處,ξ=2%和20%時(shí),殘差相對(duì)于Mw和D5-95的分布。根據(jù)殘差分析可知:在短周期處(如T=0.1 s),采用式(11)的回歸方程對(duì)殘差相對(duì)于Mw和D5-95的分布無(wú)明顯影響(將圖6、圖7與圖9、圖10對(duì)比可知),再次說(shuō)明在短周期處Mw和D5-95對(duì)DSF無(wú)顯著影響,各阻尼比的殘差隨機(jī)對(duì)稱的分布于零水平線的兩側(cè);在長(zhǎng)周期處(如T=3 s),采用式(10)的回歸方程使得殘差相對(duì)于Mw不再具有直線分布,各阻尼比下的殘差隨機(jī)對(duì)稱的分布于零水平線的兩側(cè)(見(jiàn)圖9),這說(shuō)明了本文回歸模型的合理性且能反映出長(zhǎng)周期處Mw對(duì)DSF的影響;此外,采用式(10)的回歸方程極大地改善了長(zhǎng)周期處殘差相對(duì)于D5-95的分布(對(duì)比圖7和圖10),間接地反映了D5-95對(duì)DSF的影響。圖11給出了在ξ=0.5%、2%、7%、15%和25%處,按本文模型(取Mw=6.5)所計(jì)算的DSF中值與本文數(shù)據(jù)(Mw=6~7)所計(jì)算的DSF中值的對(duì)比。由圖可見(jiàn),兩者吻合得較好。

    現(xiàn)以一算例進(jìn)一步說(shuō)明所建議模型的有效性及模型中考慮矩震級(jí)影響的必要性。對(duì)于基地隔震結(jié)構(gòu),當(dāng)上部結(jié)構(gòu)的抗側(cè)剛度遠(yuǎn)大于隔震層的抗側(cè)剛度時(shí),可將其簡(jiǎn)化為一單自由度體系。設(shè)有一基地隔震結(jié)構(gòu),上部結(jié)構(gòu)的抗側(cè)剛度遠(yuǎn)大于隔震層的抗側(cè)剛度,基本周期為3 s,阻尼比為30%。該結(jié)構(gòu)可簡(jiǎn)化為自振周期為3 s,阻尼比為30%的單自由度體系。當(dāng)該體系遭受表1中矩震級(jí)Mw=7~8的地面運(yùn)動(dòng)記錄激勵(lì)時(shí),通過(guò)分析可求得該體系精確的阻尼調(diào)整系數(shù)中值為0.488(圖3)。對(duì)式(5)取指數(shù),可求得不考慮矩震級(jí)影響的阻尼調(diào)整系數(shù)模型的阻尼調(diào)整系數(shù)中值為0.579。取Mw=7.5,采用本文所建議的考慮矩震級(jí)影響的阻尼調(diào)整系數(shù)模型所求得的阻尼調(diào)整系數(shù)中值為0.48。由上述計(jì)算可見(jiàn):采用不考慮矩震級(jí)影響的阻尼調(diào)整系數(shù)模型計(jì)算的阻尼調(diào)整系數(shù)中值所產(chǎn)生的誤差為+18.6%(“+”代表高估);采用本文建議的考慮矩震級(jí)影響的阻尼調(diào)整系數(shù)模型計(jì)算的阻尼調(diào)整系數(shù)中值所產(chǎn)生的誤差僅為-1.6%(“-”代表低估),在工程允許的誤差范圍內(nèi)。以上分析說(shuō)明了考慮矩震級(jí)對(duì)阻尼調(diào)整系數(shù)影響的必要性,尤其對(duì)于長(zhǎng)周期大阻尼比的結(jié)構(gòu),不考慮矩震級(jí)的影響將帶來(lái)較大的誤差,也說(shuō)明了本文模型的有效性和優(yōu)越性。

    圖8 ^c0,^c1和擬合曲線Fig.8 ^c0,^c1and fitted curve

    圖9 殘差相對(duì)于Mw的分布Fig.9 Residuals vs.Mw

    表2 模型回歸系數(shù)(×10-2)Table 2 Regression coefficients for themodel(×10-2)

    圖10 殘差相對(duì)于D5-95的分布Fig.10 Residuals vs.D5-95

    圖11 本文模型和數(shù)據(jù)的DSF中值對(duì)比Fig.11 Comparison on median DSF between the model and data

    5.2 標(biāo)準(zhǔn)差σln(DSF)

    對(duì)于一元線性方程,在Ti和ξj處,^σlnDSF(ξj,Ti)的無(wú)偏估計(jì)可用下式[18]計(jì)算:在指定的Ti處,用式(12)可求得11個(gè)阻尼比處

    標(biāo)準(zhǔn)差的估計(jì)值。在Ti處,本文采用的函數(shù)形式來(lái)擬合標(biāo)準(zhǔn)差的估計(jì)值與ξ之間的關(guān)系。?。?/p>

    圖12給出了在T=0.1 s和3 s處式(13)的擬合結(jié)果。由圖可見(jiàn),DSF相對(duì)于ξ存在異方差性。隨著ξ越遠(yuǎn)離5%,DSF的離散程度越大。這與前文的分析結(jié)果是一致的。各周期點(diǎn)處回歸系數(shù)的估計(jì)值列于表2。

    圖12 ^σln(DSF)和擬合曲線Fig.12 ^σln(DSF)and fitted curve

    6 結(jié) 論

    基于408條地震地面運(yùn)動(dòng)記錄,本文研究了Mw對(duì)DSF的影響,并在此基礎(chǔ)上提出了考慮Mw影響的回歸模型,得出以下結(jié)論:

    (1)Mw對(duì)DSF的影響與T和ξ有關(guān);在T>0.6 s的范圍內(nèi),當(dāng)ξ<5%時(shí),DSF隨著Mw的增長(zhǎng)而增長(zhǎng);當(dāng)ξ>5%時(shí),DSF隨著Mw的增長(zhǎng)而減小。隨著T的增長(zhǎng),Mw對(duì)DSF的影響越顯著;ξ越遠(yuǎn)離5%,Mw對(duì)DSF的影響越顯著。在T<0.6 s的范圍內(nèi),Mw對(duì)DSF無(wú)顯著影響;

    (2)當(dāng)Mw=5~6時(shí),DSF受T的影響顯著;隨著Mw的增大,在T>0.6 s范圍內(nèi),DSF受T的影響逐漸減弱;當(dāng)Mw=7~8時(shí),除ξ=0.5%和1%外,在T>0.6 s范圍內(nèi),DSF基本不隨T而變化。

    (3)DSF相對(duì)于ξ存在異方差性;ξ越遠(yuǎn)離5%,DSF的離散程度越大。

    (4)本文所提出的回歸模型能體現(xiàn)Mw對(duì)DSF的影響,及DSF相對(duì)于ξ的異方差性,也能間接地體現(xiàn)D5-95對(duì)DSF的影響。

    [1] 中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.GB 50011—2010建筑抗震設(shè)計(jì)規(guī)范[S].北京:中國(guó)建筑工業(yè)出版社,2010.Ministry of Housing and Urban-Rural Department of the People’s Republic of China.GB 50011—2010 Code for seismic design of buildings[S].Beijing:China Architecture and Building Press,2010.(in Chinese)

    [2] UBC.Uniform building code[S].Whittier,CA:International Conference of Building Officials,1994,1997.

    [3] IBC.International building code[S].Whittier,CA:International Conference of Building Officials,2000.

    [4] NEHRP.Recommended provisions for seismic regulations for new buildings[R].Washington,DC:Federal Emergency Management Agency,1994,1997,2000.

    [5] FEMA-273.NEHRPGuidelines for the seismic rehabilitation of buildings[R].Washington,DC:Federal Emergency Management Agency,1997.

    [6] Eurocode 8.Design of structures for earthquake resistance,part 1:general rules,seismic actions and rules for buildings[S].Brussels:European Committee for Standardization,2004.

    [7] Ashour SA.Elastic seismic response of buildingswith supplemental damping[D].Ann Arbor:University of Michigan,1987.

    [8] Tolis SV,F(xiàn)accioli E.Displacement design spectra[J].Journal of Structural Engineering,1999,3(1):107-125.

    [9] Priestley M.Myths and fallacies in earthquake engineering,revisited[M].Pavia(Italy):IUSS Press,2003.

    [10] Newmark N M,Hall W J.Earthquake spectra and design[M].Oakland:Earthquake Engineering Research Institute,1982.

    [11] Wu J,Hanson R D.Study of inelastic spectra with high damping[J].Journal of Structural Engineering,1989,115(6):1412-1431.

    [12] Idriss IM.Procedures for selecting earthquake ground motions at rock sites[R].Report NISTGCR 93-625,Washington,D.E:National Institute of Standards and Technology,1993.

    [13] Lin Y Y,Chang K C.Effects of site class on damping reduction factors[J].Journal of Structural Engineering,2004,130(11):1667-1675.

    [14] Bommer J J,Mendis R.Scaling of spectral displacement ordinates with damping ratios[J].Earthquake Engineering and Structural Dynamics,2005,34(2):145-165.

    [15] Cameron W I,Green I.Damping correction factors for horizontal ground-motion response spectra[J].Bulletin of the Seismological Society of America,2007,97(3):934-960.

    [16] Kempton J J,Stewart J P.Prediction equations for significant duration ofearthquake groundmotions considering site and near-source effects[J],Earthquake Spectra,2006,22(4):985-1013.

    [17] 謝禮立,張曉志.地震動(dòng)記錄持時(shí)和工程持時(shí)[J].地震工程與工程振動(dòng),1988,8(1):31-38.Xie Lili,Zhang Xiaozhi.Accelerogram-based duration and engineering duration of ground motion[J].Earthquake Engineering and Engineering Vibration.1988,8(1):31-38.(in Chinese)

    [18] Ang A H S,TangW H.Probability concepts in engineering planning and design volume1 basic principles[M].New York:John Wily and Sons,Inc,1975.

    Damping Scaling Factor w ith a Consideration of Influence of M oment M agnitude

    WANG Guotao*HU Kexu LEIMin
    (Resecnch Institute of Struetural Engineering and Disaster Reduction,Tongji University,Shanghai200092,China)

    The influence ofmomentmagnitude(Mw)on damping scaling factor(DSF)was firstly investigated based on the statistical analysis of408 earthquake ground motion records.It is showed that the effect of Mwon DSF depends on period of vibration(T)and damping ratio(ξ).The DSF increases as Mwincreases for the periods greater than about0.6 s ifξ<5%,but it decreases ifξ>5%.The pattern with Mwismuchmore significant as T increases orξdeviates from 5%.Almostno pattern with Mwis seen for the periods less than about 0.6 s.DSF is significantly affected by T when Mw=5~6.The influence of T on DSF becomes less significant with the increase of Mwfor the periods greater than about 0.6 s.When Mw=7~8 DSF does not basically change with T within the range of periods greater than 0.6 s except DSF corresponding to 0.5%and 1% damping.The heteroscedasticity in the DSF with respect toξexists.Then DSF regression equation was proposed that only depends on T andξand equation residualswas calculated at specified T andξ.Finally according to the residual plots showing an almost linear relationship with Mwat long periods,the regressionmodel that includes a linearmagnitude term was proposed and regression coefficients were tabulated.Residual analysis shows that the proposed model can reflect the influence of Mwon DSF and can indirectly reflect the influence of significant duration(D5-95)on DSF.The standard deviation of themodel can reflect heteroscedasticity in the DSF with respect toξ.

    2013-06-17

    國(guó)家科技支撐計(jì)劃重點(diǎn)項(xiàng)目(2009BAJ28B02)*聯(lián)系作者,Email:802212wgt@#edu.cn

    猜你喜歡
    影響模型
    一半模型
    是什么影響了滑動(dòng)摩擦力的大小
    哪些顧慮影響擔(dān)當(dāng)?
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    沒(méi)錯(cuò),痛經(jīng)有時(shí)也會(huì)影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    3D打印中的模型分割與打包
    擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
    基于Simulink的跟蹤干擾對(duì)跳頻通信的影響
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    69av精品久久久久久| 99精品久久久久人妻精品| 美女高潮的动态| 欧美最黄视频在线播放免费| av天堂中文字幕网| 色综合亚洲欧美另类图片| 国产伦人伦偷精品视频| 性欧美人与动物交配| 国产成人a区在线观看| 久久久久免费精品人妻一区二区| 国产av在哪里看| 国产精品三级大全| 可以在线观看的亚洲视频| 观看免费一级毛片| 特级一级黄色大片| 成人三级黄色视频| 欧美成人性av电影在线观看| 精品日产1卡2卡| 欧美高清成人免费视频www| 一区二区三区高清视频在线| 国产一区二区亚洲精品在线观看| 三级毛片av免费| 免费人成在线观看视频色| 亚洲国产高清在线一区二区三| 草草在线视频免费看| 亚洲欧美日韩卡通动漫| 国产精品亚洲av一区麻豆| 看黄色毛片网站| 欧美一区二区精品小视频在线| 亚州av有码| 97热精品久久久久久| 国产精品久久电影中文字幕| 亚洲专区国产一区二区| 亚洲av美国av| 午夜福利视频1000在线观看| 免费人成在线观看视频色| 不卡一级毛片| 床上黄色一级片| 真人一进一出gif抽搐免费| 最近最新中文字幕大全电影3| av福利片在线观看| 深夜精品福利| 精品欧美国产一区二区三| 男女那种视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲av日韩精品久久久久久密| 午夜精品久久久久久毛片777| 精品一区二区三区av网在线观看| 日本免费一区二区三区高清不卡| 中文字幕熟女人妻在线| 国产v大片淫在线免费观看| 欧美激情国产日韩精品一区| 久久精品久久久久久噜噜老黄 | 欧美zozozo另类| 又爽又黄无遮挡网站| 国产成+人综合+亚洲专区| 人人妻,人人澡人人爽秒播| a级毛片免费高清观看在线播放| 亚洲国产色片| 欧美高清成人免费视频www| 欧美又色又爽又黄视频| 51午夜福利影视在线观看| 亚洲人成伊人成综合网2020| bbb黄色大片| 禁无遮挡网站| 国产精品久久久久久人妻精品电影| 色综合站精品国产| 成人特级av手机在线观看| 午夜福利免费观看在线| 亚洲狠狠婷婷综合久久图片| 欧美xxxx黑人xx丫x性爽| 国产精品久久视频播放| 欧美成人一区二区免费高清观看| 一区二区三区免费毛片| 亚洲人与动物交配视频| 人妻久久中文字幕网| 又黄又爽又刺激的免费视频.| 人人妻人人看人人澡| 小说图片视频综合网站| 美女黄网站色视频| 久99久视频精品免费| 美女被艹到高潮喷水动态| 久久久久久九九精品二区国产| 有码 亚洲区| 狂野欧美白嫩少妇大欣赏| 波多野结衣巨乳人妻| 久久精品人妻少妇| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av天美| 国产精品1区2区在线观看.| 草草在线视频免费看| 十八禁人妻一区二区| 久久久精品大字幕| 色哟哟哟哟哟哟| 亚洲精品粉嫩美女一区| 日韩有码中文字幕| 亚洲最大成人中文| 亚洲av美国av| 两人在一起打扑克的视频| 搞女人的毛片| 精品乱码久久久久久99久播| 免费观看的影片在线观看| 99热精品在线国产| 亚洲av成人不卡在线观看播放网| 欧美成人性av电影在线观看| 久久久精品大字幕| 国产精品久久久久久亚洲av鲁大| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 老熟妇仑乱视频hdxx| 综合色av麻豆| 如何舔出高潮| av在线蜜桃| 国产一区二区在线av高清观看| 午夜激情欧美在线| 国语自产精品视频在线第100页| 久久99热这里只有精品18| 桃红色精品国产亚洲av| 91九色精品人成在线观看| 3wmmmm亚洲av在线观看| 久久久色成人| 国产精品亚洲美女久久久| 在线观看一区二区三区| 精品欧美国产一区二区三| 中文字幕高清在线视频| 精品国产三级普通话版| 国产淫片久久久久久久久 | 日本免费a在线| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区免费观看 | 国产视频内射| 国产精品1区2区在线观看.| 国产成年人精品一区二区| 国产淫片久久久久久久久 | 99在线人妻在线中文字幕| 久久久国产成人免费| 色综合婷婷激情| 日韩成人在线观看一区二区三区| 午夜免费激情av| 国产激情偷乱视频一区二区| 久久国产乱子免费精品| 91久久精品国产一区二区成人| 亚洲一区高清亚洲精品| 久久精品91蜜桃| 床上黄色一级片| 亚洲一区高清亚洲精品| 日本黄大片高清| 亚洲无线在线观看| 久久精品国产亚洲av天美| 国内精品美女久久久久久| 最近视频中文字幕2019在线8| 欧美成人免费av一区二区三区| 一级黄色大片毛片| av在线天堂中文字幕| 亚洲人成网站在线播| 999久久久精品免费观看国产| 国产精品三级大全| 丰满人妻熟妇乱又伦精品不卡| 欧美三级亚洲精品| 欧美一区二区国产精品久久精品| 精品久久久久久成人av| 亚洲av成人不卡在线观看播放网| 在线观看午夜福利视频| 中文字幕高清在线视频| 国产精品,欧美在线| 能在线免费观看的黄片| 特级一级黄色大片| 久久久久久久久久黄片| 性欧美人与动物交配| 99视频精品全部免费 在线| 免费av毛片视频| 国产精品自产拍在线观看55亚洲| 国产 一区 欧美 日韩| 黄色配什么色好看| 欧美午夜高清在线| 亚洲人成电影免费在线| 久久久久久久久久成人| 国产高清视频在线观看网站| 又爽又黄无遮挡网站| 午夜两性在线视频| 51午夜福利影视在线观看| 久久国产精品人妻蜜桃| 97人妻精品一区二区三区麻豆| 久久久久国内视频| 亚洲最大成人中文| 国产精品爽爽va在线观看网站| 日本三级黄在线观看| 1000部很黄的大片| 日本 欧美在线| 亚洲av中文字字幕乱码综合| 欧美午夜高清在线| 91麻豆av在线| 久久精品人妻少妇| 亚洲av免费高清在线观看| 韩国av一区二区三区四区| 亚洲综合色惰| 色综合欧美亚洲国产小说| 永久网站在线| 性色av乱码一区二区三区2| 乱码一卡2卡4卡精品| 精品人妻1区二区| 国产高清有码在线观看视频| 精品日产1卡2卡| 久久久久免费精品人妻一区二区| 变态另类成人亚洲欧美熟女| 亚洲精品日韩av片在线观看| 天堂√8在线中文| 毛片一级片免费看久久久久 | 亚洲av日韩精品久久久久久密| 日韩高清综合在线| 国产麻豆成人av免费视频| 窝窝影院91人妻| 一个人看的www免费观看视频| 他把我摸到了高潮在线观看| 午夜福利在线观看吧| 亚洲欧美精品综合久久99| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲专区中文字幕在线| 久久人人爽人人爽人人片va | 中文字幕av在线有码专区| 好男人电影高清在线观看| а√天堂www在线а√下载| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 亚洲最大成人手机在线| 国产黄片美女视频| 偷拍熟女少妇极品色| 精品久久久久久久久久免费视频| 在线观看一区二区三区| 首页视频小说图片口味搜索| 我的老师免费观看完整版| 日韩av在线大香蕉| 久久久久久久久久成人| 天堂√8在线中文| 国产精品久久电影中文字幕| 久久国产乱子免费精品| 日本撒尿小便嘘嘘汇集6| 动漫黄色视频在线观看| 啪啪无遮挡十八禁网站| 精品久久久久久久人妻蜜臀av| 国产主播在线观看一区二区| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 久久天躁狠狠躁夜夜2o2o| 9191精品国产免费久久| 久久久久久久久久黄片| 尤物成人国产欧美一区二区三区| 午夜福利18| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 欧美性猛交黑人性爽| 黄色丝袜av网址大全| 亚洲美女黄片视频| 无人区码免费观看不卡| 不卡一级毛片| 男插女下体视频免费在线播放| 久久久久久久久久成人| 最近视频中文字幕2019在线8| 亚洲在线自拍视频| 色吧在线观看| 亚洲美女搞黄在线观看 | 一进一出抽搐gif免费好疼| 久久亚洲真实| 日本熟妇午夜| 午夜福利免费观看在线| 国产激情偷乱视频一区二区| 国产av不卡久久| 国产久久久一区二区三区| 久久久久亚洲av毛片大全| 欧美乱妇无乱码| 757午夜福利合集在线观看| 亚洲av电影在线进入| 观看免费一级毛片| 久久精品国产清高在天天线| 日本三级黄在线观看| 国内精品久久久久精免费| av国产免费在线观看| 一进一出抽搐动态| 在线天堂最新版资源| 69人妻影院| 精品久久久久久久人妻蜜臀av| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 亚洲精品日韩av片在线观看| 成人鲁丝片一二三区免费| 91久久精品电影网| 成年人黄色毛片网站| av福利片在线观看| 最近最新免费中文字幕在线| h日本视频在线播放| 欧美日韩乱码在线| 欧美乱色亚洲激情| 舔av片在线| 99精品在免费线老司机午夜| 深夜精品福利| 国产成人a区在线观看| 色噜噜av男人的天堂激情| 色视频www国产| 亚洲 国产 在线| 久久久色成人| 欧美日韩黄片免| 老鸭窝网址在线观看| 日韩国内少妇激情av| or卡值多少钱| 欧美激情在线99| 国产爱豆传媒在线观看| a级毛片免费高清观看在线播放| 亚洲第一欧美日韩一区二区三区| av天堂中文字幕网| 91麻豆精品激情在线观看国产| 亚洲电影在线观看av| 久久久久九九精品影院| 亚洲在线自拍视频| 亚洲第一区二区三区不卡| 国产精品三级大全| 色av中文字幕| 永久网站在线| 一夜夜www| 中文字幕久久专区| 悠悠久久av| 一区二区三区四区激情视频 | 久久精品国产99精品国产亚洲性色| 国产成+人综合+亚洲专区| 日本精品一区二区三区蜜桃| 两人在一起打扑克的视频| 欧美不卡视频在线免费观看| 精品久久国产蜜桃| 蜜桃亚洲精品一区二区三区| 亚洲成人久久性| 日日干狠狠操夜夜爽| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 九色成人免费人妻av| 两人在一起打扑克的视频| 国产精品亚洲美女久久久| 国产精品久久久久久精品电影| 欧美黄色淫秽网站| 蜜桃亚洲精品一区二区三区| 夜夜夜夜夜久久久久| 亚洲中文字幕日韩| 色哟哟·www| 婷婷丁香在线五月| ponron亚洲| 亚洲av美国av| 国产69精品久久久久777片| 午夜福利在线观看免费完整高清在 | 变态另类丝袜制服| 亚洲精品456在线播放app | 国产极品精品免费视频能看的| 91久久精品国产一区二区成人| 人人妻人人看人人澡| 国产欧美日韩精品亚洲av| av黄色大香蕉| av天堂在线播放| 午夜老司机福利剧场| 亚洲最大成人中文| 舔av片在线| 亚洲人成网站高清观看| 此物有八面人人有两片| 国产精品久久视频播放| 少妇的逼水好多| 美女黄网站色视频| 成熟少妇高潮喷水视频| 一个人看的www免费观看视频| 国产高清有码在线观看视频| 成人美女网站在线观看视频| 久久亚洲精品不卡| 亚洲欧美激情综合另类| 精品一区二区三区视频在线观看免费| 日韩av在线大香蕉| 欧美绝顶高潮抽搐喷水| 乱人视频在线观看| 成人欧美大片| 欧美午夜高清在线| xxxwww97欧美| 老熟妇仑乱视频hdxx| 久久九九热精品免费| 久久国产精品影院| 色播亚洲综合网| 伊人久久精品亚洲午夜| 极品教师在线免费播放| 丁香六月欧美| 波多野结衣高清作品| 国产 一区 欧美 日韩| www.色视频.com| 亚洲人成网站在线播| 精华霜和精华液先用哪个| 日本黄色片子视频| 亚洲av成人不卡在线观看播放网| 99在线视频只有这里精品首页| 日韩中文字幕欧美一区二区| 夜夜爽天天搞| 日本熟妇午夜| 天美传媒精品一区二区| 午夜免费男女啪啪视频观看 | 最近最新中文字幕大全电影3| 欧美日韩综合久久久久久 | 国产中年淑女户外野战色| av天堂中文字幕网| 久久99热这里只有精品18| 三级国产精品欧美在线观看| 97超视频在线观看视频| av女优亚洲男人天堂| 亚洲av五月六月丁香网| 国产一区二区激情短视频| 国产高清三级在线| 久久精品国产自在天天线| 给我免费播放毛片高清在线观看| 亚洲最大成人中文| 久久久久久大精品| 中文字幕人妻熟人妻熟丝袜美| 欧美又色又爽又黄视频| a在线观看视频网站| 人妻制服诱惑在线中文字幕| 琪琪午夜伦伦电影理论片6080| 国产又黄又爽又无遮挡在线| 亚洲精品日韩av片在线观看| 成人性生交大片免费视频hd| 国语自产精品视频在线第100页| 美女xxoo啪啪120秒动态图 | 丁香欧美五月| 中文字幕精品亚洲无线码一区| 一区二区三区激情视频| 亚洲欧美日韩卡通动漫| 国产熟女xx| 哪里可以看免费的av片| 午夜久久久久精精品| avwww免费| 欧美黄色淫秽网站| 亚洲精品在线观看二区| 久久草成人影院| 国产伦精品一区二区三区视频9| 久久国产精品人妻蜜桃| 女生性感内裤真人,穿戴方法视频| 9191精品国产免费久久| 免费人成视频x8x8入口观看| 国产午夜福利久久久久久| 亚洲成人久久爱视频| 村上凉子中文字幕在线| 日本三级黄在线观看| 97超视频在线观看视频| 人人妻人人看人人澡| 精品免费久久久久久久清纯| 九九在线视频观看精品| 精品熟女少妇八av免费久了| 亚洲国产精品sss在线观看| 久久精品国产亚洲av天美| 美女高潮喷水抽搐中文字幕| 波多野结衣巨乳人妻| 制服丝袜大香蕉在线| 白带黄色成豆腐渣| 在现免费观看毛片| 日本与韩国留学比较| 欧美日韩亚洲国产一区二区在线观看| 97超视频在线观看视频| 国产精品人妻久久久久久| 性色av乱码一区二区三区2| 免费av不卡在线播放| 好男人在线观看高清免费视频| 精品人妻偷拍中文字幕| 特大巨黑吊av在线直播| 99久久无色码亚洲精品果冻| 麻豆成人av在线观看| 麻豆国产av国片精品| 精品人妻1区二区| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产清高在天天线| 少妇的逼水好多| а√天堂www在线а√下载| 欧美区成人在线视频| 最近在线观看免费完整版| 免费看a级黄色片| 亚洲av电影不卡..在线观看| 黄片小视频在线播放| 小蜜桃在线观看免费完整版高清| 成人午夜高清在线视频| 色哟哟哟哟哟哟| 夜夜夜夜夜久久久久| 久久婷婷人人爽人人干人人爱| 精品一区二区三区av网在线观看| 欧美性猛交╳xxx乱大交人| 国产精品国产高清国产av| 亚洲最大成人手机在线| 亚洲乱码一区二区免费版| 欧美国产日韩亚洲一区| 国产精品日韩av在线免费观看| 熟妇人妻久久中文字幕3abv| 日本一二三区视频观看| av在线观看视频网站免费| 欧洲精品卡2卡3卡4卡5卡区| 精品久久国产蜜桃| 黄色一级大片看看| 国产精品嫩草影院av在线观看 | www.色视频.com| 国产av麻豆久久久久久久| 日本a在线网址| 久久香蕉精品热| 亚洲人成伊人成综合网2020| avwww免费| 国产精品电影一区二区三区| 一区二区三区四区激情视频 | 亚洲一区高清亚洲精品| 色吧在线观看| 搡老熟女国产l中国老女人| 国产精品久久视频播放| 日韩精品青青久久久久久| 色综合欧美亚洲国产小说| 国产高清视频在线观看网站| 一个人看视频在线观看www免费| 亚洲欧美清纯卡通| 99国产极品粉嫩在线观看| 亚洲avbb在线观看| 国产一区二区激情短视频| 黄色一级大片看看| 简卡轻食公司| 亚洲片人在线观看| 午夜视频国产福利| 亚洲黑人精品在线| 亚洲 国产 在线| 亚洲男人的天堂狠狠| 十八禁网站免费在线| 老司机深夜福利视频在线观看| 俺也久久电影网| 最好的美女福利视频网| 直男gayav资源| 亚洲无线在线观看| 网址你懂的国产日韩在线| 一进一出抽搐动态| 九九久久精品国产亚洲av麻豆| av专区在线播放| 亚洲七黄色美女视频| 亚洲 欧美 日韩 在线 免费| 国产亚洲欧美98| 久久久久久久久久黄片| 久久精品国产亚洲av涩爱 | 婷婷色综合大香蕉| 看免费av毛片| 亚洲黑人精品在线| 此物有八面人人有两片| 天堂影院成人在线观看| 18禁在线播放成人免费| 夜夜看夜夜爽夜夜摸| 成年免费大片在线观看| 97超视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 真人一进一出gif抽搐免费| 亚洲精品成人久久久久久| 亚洲在线观看片| 嫩草影院新地址| 99久久久亚洲精品蜜臀av| 精品久久久久久成人av| 久久久久久大精品| 久久香蕉精品热| 国产精品亚洲av一区麻豆| 波多野结衣高清作品| 久久久久国产精品人妻aⅴ院| 欧美日韩福利视频一区二区| 国产伦精品一区二区三区视频9| 18+在线观看网站| 国产成人影院久久av| 村上凉子中文字幕在线| 人妻丰满熟妇av一区二区三区| 亚洲av熟女| 最新中文字幕久久久久| 日日夜夜操网爽| 亚洲精品影视一区二区三区av| 波多野结衣高清无吗| 国产aⅴ精品一区二区三区波| 少妇裸体淫交视频免费看高清| 国产精品一区二区性色av| 亚洲avbb在线观看| or卡值多少钱| 亚洲精品在线美女| 色av中文字幕| 综合色av麻豆| 动漫黄色视频在线观看| 亚洲av电影不卡..在线观看| 欧美性猛交黑人性爽| 久久香蕉精品热| 亚洲欧美激情综合另类| 欧美黑人欧美精品刺激| 国产精品一区二区三区四区免费观看 | 一二三四社区在线视频社区8| 国产精品一区二区三区四区免费观看 | 亚洲av日韩精品久久久久久密| 欧美区成人在线视频| 成年女人毛片免费观看观看9| 麻豆国产97在线/欧美| 又紧又爽又黄一区二区| 久久精品国产99精品国产亚洲性色| 亚洲熟妇熟女久久| 久久国产精品影院| 一级av片app| 成年女人永久免费观看视频| 国产视频一区二区在线看| 性欧美人与动物交配| 国产欧美日韩精品亚洲av| 波野结衣二区三区在线| 99热这里只有是精品在线观看 | 一个人看视频在线观看www免费| 精品熟女少妇八av免费久了| 亚洲av.av天堂| 欧美日韩瑟瑟在线播放| 免费在线观看亚洲国产| 欧美精品国产亚洲| 免费av不卡在线播放| 男人狂女人下面高潮的视频| 成年女人永久免费观看视频| 久久香蕉精品热| 亚洲av电影不卡..在线观看| 最近视频中文字幕2019在线8| av在线蜜桃| 精品人妻一区二区三区麻豆 | 久久婷婷人人爽人人干人人爱| 国产伦一二天堂av在线观看| 国产欧美日韩一区二区精品|