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

    橋梁地震易損性函數(shù)中不確定性參數(shù)的影響規(guī)律分析*

    2019-01-24 12:55:38吳文朋
    關(guān)鍵詞:易損性不確定性抗震

    吳文朋, 梁 鵬

    (1.湘潭大學(xué) 土木工程與力學(xué)學(xué)院,湖南 湘潭 411105;2.長(zhǎng)沙理工大學(xué) 南方地區(qū)橋梁長(zhǎng)期性能提升技術(shù)國(guó)家地方聯(lián)合工程實(shí)驗(yàn)室,湖南 長(zhǎng)沙 4100143.湘潭大學(xué) 工程結(jié)構(gòu)動(dòng)力學(xué)與可靠性分析湖南省高等學(xué)校重點(diǎn)實(shí)驗(yàn)室,湖南 湘潭 411105)

    相比傳統(tǒng)確定性的結(jié)構(gòu)抗震分析方法,地震易損性分析的優(yōu)點(diǎn)是能較好地考慮地震波和結(jié)構(gòu)模型等不確定性對(duì)地震響應(yīng)的影響.然而,如何量化各類(lèi)不確定性的影響卻又是易損性分析面臨的挑戰(zhàn)之一.結(jié)構(gòu)地震易損性分析中的不確定性反映了對(duì)結(jié)構(gòu)動(dòng)力響應(yīng)或損傷概率的預(yù)測(cè)精度.不確定性主要包括偶然不確定性和認(rèn)知不確定性?xún)纱箢?lèi)[1].偶然不確定性是事物或事件本身所固有的且不可預(yù)測(cè)的屬性.例如,地震事件的震級(jí)和震中距、鋼筋的彈性模量和混凝土抗壓強(qiáng)度等,因此,這類(lèi)不確定性是無(wú)法徹底消除的.認(rèn)知不確定性是由于人類(lèi)目前掌握的知識(shí)水平的不足所導(dǎo)致的.例如,在分析過(guò)程中采用的各種模型假定(能力模型、需求模型)或收集數(shù)據(jù)的不完整性等.因此,當(dāng)人類(lèi)的知識(shí)水平上升到較高的高度或收集的數(shù)據(jù)量足夠大時(shí),這類(lèi)不確定性在理論上是可以避免的.從概念上容易區(qū)分這兩類(lèi)不確定性的差別,但是在實(shí)際應(yīng)用中卻很難將兩者完全分離.

    Wen等[1-2]在對(duì)大量的文獻(xiàn)進(jìn)行綜述以后認(rèn)為,結(jié)構(gòu)抗震性能評(píng)估中的不確定性主要反映在結(jié)構(gòu)抗震能力和地震需求上.美國(guó)聯(lián)邦應(yīng)急管理局于2009年出版了《FEMA-P695研究報(bào)告》[3],該報(bào)告提出了一種用于量化建筑結(jié)構(gòu)系統(tǒng)可靠度的抗震性能評(píng)估方法(簡(jiǎn)稱(chēng)“P695方法”).P695方法中詳細(xì)闡述了在結(jié)構(gòu)倒塌易損性函數(shù)中的各種不確定性確定方法.2001年美國(guó)聯(lián)邦應(yīng)急管理局和美國(guó)加州應(yīng)用技術(shù)委員會(huì)展開(kāi)合作,一起實(shí)施了ATC-58研究計(jì)劃.該計(jì)劃第一個(gè)階段(2000-2006年)的任務(wù)是確定基于性能抗震設(shè)計(jì)的性能目標(biāo),并形成明確的研究計(jì)劃,其主要研究成果是提出了下一代基于性能的抗震設(shè)計(jì)指南FEMA-445[4].近年來(lái),國(guó)內(nèi)外大量有關(guān)基于性能結(jié)構(gòu)抗震設(shè)計(jì)的研究成果都是依據(jù)該設(shè)計(jì)指南展開(kāi)的.該計(jì)劃第二個(gè)階段(2006-2012年)在第一階段研究成果的基礎(chǔ)之上,進(jìn)一步發(fā)展了基于性能的抗震性能評(píng)估方法,其主要成果是2012年出版的一系列報(bào)告,其中FEMA-P58-1 Methodology[5](簡(jiǎn)稱(chēng)“P58方法”)是其基本理論研究成果.

    國(guó)內(nèi)關(guān)于結(jié)構(gòu)地震易損性的眾多研究中,部分研究都只考慮了地震動(dòng)不確定性的影響[6-7],也有部分研究同時(shí)考慮了多種不確定性因素的影響,例如,李立峰等[8]、陳力波等[9]在易損性研究中同時(shí)考慮了地震動(dòng)和結(jié)構(gòu)主要材料參數(shù)的不確定性.然而,以上研究并沒(méi)有對(duì)易損性函數(shù)中的不確定性進(jìn)行量化分析,與國(guó)外較成熟的研究成果相比還有很大差距.本文首先總結(jié)了結(jié)構(gòu)地震易損性函數(shù)中各類(lèi)不確定性的研究現(xiàn)狀;然后根據(jù)我國(guó)橋梁工程設(shè)計(jì)與施工過(guò)程中存在的實(shí)際情況,對(duì)橋梁地震易損性函數(shù)中需要考慮的各種不確定性進(jìn)行了分類(lèi);結(jié)合既往的研究中常采用的橋梁地震易損性分析函數(shù),討論了各種不確定性在既有地震易損性函數(shù)中的處理方法,并對(duì)其優(yōu)缺點(diǎn)進(jìn)行了比較;最后通過(guò)參數(shù)分析研究了不確定性對(duì)地震易損性函數(shù)的影響規(guī)律.

    1 易損性函數(shù)中不確定性的分類(lèi)

    易損性函數(shù)定義了在給定地震動(dòng)水平下結(jié)構(gòu)失效的概率,一般可用兩參數(shù)對(duì)數(shù)正態(tài)累積分布函數(shù)來(lái)表示結(jié)構(gòu)的地震易損性曲線(xiàn):

    (1)

    式中:Pf表示在IM=x的地震動(dòng)作用下結(jié)構(gòu)的需求響應(yīng)(D)超過(guò)其抗震能力(C)的概率;Φ[·]為標(biāo)準(zhǔn)正態(tài)累計(jì)分布函數(shù)(CDF);θ為易損性函數(shù)的均值;β為易損性函數(shù)的標(biāo)準(zhǔn)數(shù).實(shí)際上,結(jié)構(gòu)抗震性能評(píng)估中的不確定性主要反映在結(jié)構(gòu)抗震能力和地震需求上[1-2],抗震能力的不確定性主要反映在結(jié)構(gòu)材料特性和幾何尺寸的隨機(jī)性上,地震需求的不確定性主要反映在地震波庫(kù)的選擇以及地震動(dòng)參數(shù)選擇中.Wen和Ellingwood 等[1-2]提出了同時(shí)考慮多種不確定性的易損性函數(shù)計(jì)算式:

    (2)

    FEMA-P695報(bào)告[3]中給出了結(jié)構(gòu)倒塌易損性函數(shù)中的各種不確定性的確定方法.該報(bào)告將影響結(jié)構(gòu)倒塌易損性的重要不確定性因素分為四個(gè)方面,并認(rèn)為這四個(gè)方面的不確定性在統(tǒng)計(jì)意義上是相互獨(dú)立的,因此系統(tǒng)總的不確定性用對(duì)數(shù)標(biāo)準(zhǔn)差可表示為:

    (3)

    式中:βRTR表示地震波對(duì)地震波的不確定性,反映了地震波不同而導(dǎo)致結(jié)構(gòu)地震響應(yīng)的離散性,βRTR主要來(lái)源于兩個(gè)方面:(1) 工程場(chǎng)地危險(xiǎn)性導(dǎo)致的選波差異性;(2) 地震波本身頻譜和動(dòng)力特性的差異性.βDR表示與設(shè)計(jì)要求相關(guān)的不確定性,反映了結(jié)構(gòu)設(shè)計(jì)要求的完整性、魯棒性以及對(duì)結(jié)構(gòu)倒塌失效防范的可靠性.βTD表示與試驗(yàn)數(shù)據(jù)相關(guān)的不確定性,反映了用于定義結(jié)構(gòu)系統(tǒng)的試驗(yàn)數(shù)據(jù)的完整性和魯棒性.βMDL表示與結(jié)構(gòu)建模相關(guān)的不確定性,反映了結(jié)構(gòu)有限元模型的建模精度,用于衡量有限元模型是否能準(zhǔn)確表達(dá)真實(shí)結(jié)構(gòu)地震響應(yīng)的特征.盡管βTD與βMDL在一定程度上是相互關(guān)聯(lián)的,但P695方法中認(rèn)為二者是完全不相同的.

    與以往基于性能的抗震設(shè)計(jì)方法相比,P58方法[5]最大的不同是其采用了能讓業(yè)主或決策者容易理解的性能指標(biāo),同時(shí)在建立易損性函數(shù)過(guò)程中綜合考慮了多種不確定性因素.該方法從結(jié)構(gòu)抗震能力和地震需求兩個(gè)方面考慮各種不確定性的影響,并針對(duì)不同的結(jié)構(gòu)類(lèi)型給出了考慮這些不確定性因素的方法.例如,對(duì)于基于強(qiáng)度損傷的延性構(gòu)件(彎曲破壞)和脆性構(gòu)件(剪切破壞),其強(qiáng)度極限損傷狀態(tài)或能力的中值可分別表示為:

    θductile=Cqe(2.054β)φRn,

    (4)

    θbrittle=Cqe(2.81β)φRn.

    (5)

    式中:φRn表示設(shè)計(jì)規(guī)范中規(guī)定的構(gòu)件設(shè)計(jì)強(qiáng)度;Cq為考慮失效模式對(duì)施工質(zhì)量敏感性的調(diào)整系數(shù);β為損傷狀態(tài)總的不確定性,可采用下式計(jì)算:

    (6)

    式中:βDe表示與結(jié)構(gòu)設(shè)計(jì)公式相關(guān)的不確定性,對(duì)于延性構(gòu)件可取0.05或0.10,脆性構(gòu)件取0.25;βMa表示與材料強(qiáng)度相關(guān)的不確定性,根據(jù)結(jié)構(gòu)材料類(lèi)型不同而取不同的值,例如,木材、混凝土和鋼材的不確定性默認(rèn)值分別取為0.3、0.2和0.1;βCon表示與施工質(zhì)量相關(guān)的不確定性,在沒(méi)有特殊要求的情況下其默認(rèn)值可按P58方法的建議進(jìn)行取值.在地震需求分析方面,P58方法將影響結(jié)構(gòu)需求參數(shù)的不確定性分為:(1) 結(jié)構(gòu)模型的不確定性;(2) 地震波對(duì)地震波的不確定性;(3) 基于地區(qū)性衰減模型地震動(dòng)不確定性.其中后兩種不確定性與式(2)中的βRTR相當(dāng).而結(jié)構(gòu)建模的不確定性則采用下式計(jì)算:

    (7)

    式中:βc是與結(jié)構(gòu)建模過(guò)程中單元模型、阻尼參數(shù)以及質(zhì)量參數(shù)的選取相關(guān)的不確定性;βq是與非線(xiàn)性模型的建模質(zhì)量相關(guān)的不確定性.同樣,P58方法也給出了βc和βq的默認(rèn)取值.

    既往研究根據(jù)各自研究目標(biāo)的不同,在地震易損性分析中考慮的不確定性因素會(huì)有一定的差異.在建筑結(jié)構(gòu)倒塌易損性研究中,由于倒塌易損性直接與震后人員傷亡情況的評(píng)估有關(guān),且在結(jié)構(gòu)臨近倒塌時(shí)其地震響應(yīng)的離散性會(huì)更大,因此倒塌易損性函數(shù)中包括的不確定性最充分.需要特別指出的是,P695方法和P58方法主要針對(duì)結(jié)構(gòu)的倒塌地震易損性曲線(xiàn),其研究對(duì)象是各種不同類(lèi)型的建筑結(jié)構(gòu),其建議的各種不確定性的取值主要是針對(duì)大量不同建筑結(jié)構(gòu)的地震易損性函數(shù)統(tǒng)計(jì)出來(lái)的.FEMA方法主要適用于美國(guó)的建筑結(jié)構(gòu)工程.然而,一方面,以上方法是否適用于設(shè)計(jì)和施工完全不同的橋梁工程仍有待研究;另一方面,我國(guó)橋梁設(shè)計(jì)和施工的時(shí)代背景與其他國(guó)家有著很大的不同,目前在國(guó)內(nèi)的橋梁地震易損性分析中仍主要考慮輸入地震波的差異性導(dǎo)致的需求響應(yīng)的不確定性,對(duì)其他類(lèi)型不確定性因素的影響也還有待研究.因此,考慮到橋梁工程的結(jié)構(gòu)特點(diǎn)以及不確定性因素的復(fù)雜性,本文認(rèn)為橋梁地震易損性函數(shù)可以從以下幾個(gè)方面展開(kāi)不確定性的研究.

    1.1 橋梁建模相關(guān)的不確定性

    地震作用下橋梁非線(xiàn)性時(shí)程分析的結(jié)果很大程度上依賴(lài)于建立的有限元模型的準(zhǔn)確性.橋梁建模就是要從結(jié)構(gòu)體系的角度,根據(jù)橋梁幾何尺寸離散各構(gòu)件單元并精確模擬其力學(xué)特性,處理好各單元間的力學(xué)連接特性,使模型的數(shù)值分析結(jié)果盡可能準(zhǔn)確地預(yù)測(cè)實(shí)際橋梁在地震激勵(lì)下的動(dòng)力響應(yīng)[10].因此,橋梁建模相關(guān)的不確實(shí)性大小能夠衡量橋梁地震響應(yīng)預(yù)測(cè)結(jié)果的精度.橋梁建模相關(guān)的不確實(shí)性可以概括為以下兩個(gè)方面:(1) 橋梁結(jié)構(gòu)本身的不確實(shí)性.例如,設(shè)計(jì)師的判斷、施工控制的質(zhì)量、養(yǎng)護(hù)條件的好壞、橋梁設(shè)計(jì)尺寸以及工程材料的變異性等都會(huì)導(dǎo)致該類(lèi)不確定性.(2) 有限元建模技術(shù)的不確定性.例如,非線(xiàn)性單元類(lèi)型的選擇、土壤邊界和支座邊界等的模擬都會(huì)導(dǎo)致這一類(lèi)不確定性.其中,設(shè)計(jì)師的判斷、施工質(zhì)量養(yǎng)護(hù)條件等導(dǎo)致的不確定性需要通過(guò)大量調(diào)查研究得到.因此,具體研究中可分為以下三類(lèi)影響結(jié)構(gòu)建模的不確定性進(jìn)行展開(kāi):(a) 橋梁幾何尺寸變異性;(b) 邊界條件參數(shù)變異性;(c) 鋼筋和混凝土材料參數(shù)變異性.值得指出的是,由于我國(guó)目前還沒(méi)有可供研究人員使用的較完備的橋梁統(tǒng)計(jì)信息數(shù)據(jù)庫(kù),現(xiàn)有的研究仍主要偏向于特定橋梁的地震易損性分析,即研究的對(duì)象仍為單體橋梁.因此,上述橋梁幾何尺寸變異性主要是指,由施工等因素導(dǎo)致的實(shí)際構(gòu)件尺寸與原設(shè)計(jì)尺寸不完全一致,而并非區(qū)域性橋梁地震易損性分析中橋與橋之間的差異性[11].

    1.2 地震波相關(guān)的不確定性

    地震波相關(guān)的不確定性包括與場(chǎng)地類(lèi)型衰減模型相關(guān)的不確定性,以及地震波對(duì)地震波的不確定性.其中,與場(chǎng)地類(lèi)型衰減模型相關(guān)的不確定性涉及地震工程學(xué)的內(nèi)容,易損性分析中一般不予考慮,而主要通過(guò)收集不同地震波數(shù)據(jù)庫(kù)來(lái)反映地震波相關(guān)的不確定性.在橋梁非線(xiàn)性動(dòng)力分析中,地震波本身頻譜特性、持續(xù)時(shí)間等對(duì)橋梁地震響應(yīng)的影響很大,因此,地震波對(duì)地震波的變異性最終體現(xiàn)在對(duì)橋梁結(jié)構(gòu)或構(gòu)件最大響應(yīng)值預(yù)測(cè)的不確定性.地震波的不確定性是其本身的固有屬性,是人力所無(wú)法消除的.然而,在實(shí)際結(jié)構(gòu)概率地震需求分析中,可以通過(guò)一定的措施盡量減小這類(lèi)不確定性對(duì)分析結(jié)果的影響,例如,增加分析中使用的地震波數(shù)量、選擇合適的地震動(dòng)強(qiáng)度參數(shù)進(jìn)行后處理等.盡管既往橋梁地震易損性研究中引入了地震波不確定性的影響,但很少對(duì)其影響大小進(jìn)行定量的分析.橋梁地震易損性分析可以在以下幾個(gè)方面進(jìn)行:(a) 分別研究近場(chǎng)脈沖型地震波、遠(yuǎn)場(chǎng)無(wú)脈沖地震波、人工生成地震波三類(lèi)地震波庫(kù)對(duì)橋梁概率地震需求模型(PSDM)的影響.(b) 針對(duì)以上三類(lèi)地震波庫(kù),從“效率性”“實(shí)用性”“充分性”的角度,分別研究每類(lèi)地震波合理地震動(dòng)參數(shù)的選擇.

    1.3 損傷指標(biāo)的不確定性

    確定橋梁在不同損傷狀態(tài)下的損傷指標(biāo)是橋梁地震易損性分析的一個(gè)非常重要的環(huán)節(jié).橋梁地震易損性分析中常將結(jié)構(gòu)的損傷劃分為:無(wú)損傷(DS0)、輕微損傷(DS1)、中等損傷(DS2)、嚴(yán)重?fù)p傷(DS3)和完全破壞(DS4)五個(gè)等級(jí)[12].美國(guó)FEMA最新的災(zāi)害損失風(fēng)險(xiǎn)評(píng)估軟件HAZUS-MH-MR5中[13]對(duì)橋梁結(jié)構(gòu)后四種損傷狀態(tài)進(jìn)行了詳細(xì)的描述.對(duì)不同的損傷狀態(tài)有了定性的描述以后,還需要進(jìn)一步量化每一種構(gòu)件的損傷指標(biāo).損傷指標(biāo)不僅要與工程需求參數(shù)(EDP)形式相一致(如延性比、相對(duì)位移),還要能反映出對(duì)橋梁實(shí)際使用功能的影響程度.例如,以位移延性比(μd)為墩柱構(gòu)件的損傷指標(biāo),當(dāng)μEDP>μd時(shí)橋梁發(fā)生中等損傷,此時(shí)橋梁需要關(guān)閉15 d以上進(jìn)行抗震加固和維修.因此,在橋梁地震易損性分析中還將考慮不同損傷狀態(tài)下?lián)p傷指標(biāo)(極限狀態(tài))的不確定性的影響.對(duì)不同損傷狀態(tài)損傷指標(biāo)的量化通常采用兩種方法:基于計(jì)算或試驗(yàn)的規(guī)定性方法和基于統(tǒng)計(jì)調(diào)查的描述性方法[14].然而,無(wú)論采用何種方法得到的損傷極限指標(biāo)都不可能是完全確定的.首先損傷的定性描述存在主觀的不確定性,其次計(jì)算方法、試驗(yàn)過(guò)程和調(diào)查過(guò)程都存在不確定性.因此,通過(guò)調(diào)查和專(zhuān)家問(wèn)卷調(diào)研來(lái)收集信息,或者對(duì)國(guó)內(nèi)既有的較離散的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)、歸類(lèi),并建立區(qū)域性的實(shí)驗(yàn)數(shù)據(jù)庫(kù),是確定我國(guó)橋梁結(jié)構(gòu)損傷指標(biāo)不確定性的有效途徑.

    2 易損性函數(shù)中不確定性的處理方式

    目前基于非線(xiàn)性數(shù)值模擬結(jié)果的理論易損性函數(shù)主要包括“云圖法”“需求能力比對(duì)數(shù)回歸法”“縮放法”“極大似然估計(jì)法”四種類(lèi)型[15],文獻(xiàn)[15]對(duì)各種方法有詳細(xì)介紹,這里不再贅述.本文主要介紹各類(lèi)不確定性在四類(lèi)易損性函數(shù)中的處理方式并對(duì)其各自的優(yōu)缺點(diǎn)進(jìn)行比較.

    “云圖法”能夠同時(shí)考慮與橋梁建模、地震波和損傷指標(biāo)相關(guān)的不確定性.考慮橋梁建模和地震波的不確定性時(shí),所有的相關(guān)參數(shù)都不再是恒定值.有的參數(shù)被假定為服從某一分布特征的隨機(jī)變量,例如,假設(shè)墩柱混凝土的抗壓強(qiáng)度f(wàn)c服從正態(tài)分布,如圖1(a)所示,墩柱縱向鋼筋的屈服強(qiáng)度f(wàn)y服從對(duì)數(shù)正態(tài)分布,如圖1(b)所示.有的參數(shù)在一定范圍內(nèi)完全隨機(jī)分布,例如,圖1(c)所示,為地震波的震中距和震級(jí)分布情況.在假定了所有參數(shù)的隨機(jī)分布特征以后,進(jìn)行完全隨機(jī)組合或采用拉丁超立方抽樣方法進(jìn)行組合,可以生成大量的“橋梁模型-地震波”分析樣本,然后對(duì)每組樣本分別進(jìn)行非線(xiàn)性時(shí)程分析,以此來(lái)考慮結(jié)構(gòu)建模和地震波相關(guān)的不確定性.對(duì)于不同構(gòu)件損傷指標(biāo)的不確定性,則直接體現(xiàn)在易損性函數(shù)的對(duì)數(shù)標(biāo)準(zhǔn)差βc中.

    “需求能力比對(duì)數(shù)回歸法”同樣能考慮與橋梁建模、地震波和損傷指標(biāo)相關(guān)的不確定性因素.前兩類(lèi)不確定性體現(xiàn)在橋梁有限元建模和非線(xiàn)性動(dòng)力分析過(guò)程中.有所不同的是對(duì)于損傷指標(biāo)的不確定性處理.該方法假定不同損傷狀態(tài)的損傷指標(biāo)也為服從一定分布的隨機(jī)變量,并參與到拉丁超立方抽樣過(guò)程中.因此該方法最終生成一系列的“橋梁模型-地震波-損傷指標(biāo)”樣本,然后進(jìn)行非線(xiàn)性分析以及回歸分析,直接得到易損性函數(shù)的兩個(gè)估計(jì)參數(shù).然而,該方法的缺點(diǎn)是對(duì)于單個(gè)構(gòu)件的各個(gè)損傷極限狀態(tài)都得進(jìn)行對(duì)數(shù)回歸分析,當(dāng)考慮的橋梁失效模式過(guò)多時(shí),勢(shì)必增加了對(duì)數(shù)回歸分析的工作量.

    “縮放法”和 “極大似然估計(jì)法”對(duì)于三種不確定性的處理方式,與“云圖法”基本上是一致的,即在有限元建模和非線(xiàn)性時(shí)程分析的過(guò)程中通過(guò)隨機(jī)抽樣的方式引入建模相關(guān)和地震波相關(guān)的不確定性,然后通過(guò)對(duì)數(shù)標(biāo)準(zhǔn)差βc來(lái)體現(xiàn)損傷指標(biāo)的不確定性.然而,有所不同的是,“縮放法”和 “極大似然估計(jì)法”只能先建立結(jié)構(gòu)的地震易損性曲線(xiàn),再通過(guò)修正易損性函數(shù)的對(duì)數(shù)標(biāo)準(zhǔn)差來(lái)體現(xiàn)不確定性的影響.需要特別說(shuō)明的是,“縮放法”所采用的地震波輸入并非原始地震波,而需要經(jīng)過(guò)一系列的縮放處理[9].然而,一方面地震波的縮放過(guò)程本身又會(huì)引入新的不確定性;另一方面,如果考慮的不確定性因素太多,勢(shì)必導(dǎo)致非線(xiàn)性時(shí)程分析的計(jì)算量顯著增加.因此,“縮放法”不適合處理大量橋梁建模相關(guān)的不確定性.

    由此可見(jiàn),基于“云圖法”的地震易損性函數(shù)與另外三種方法相比,在同時(shí)處理多種不確定性時(shí)有一定的優(yōu)勢(shì).主要原因是該方法將地震易損性分析中需求分析過(guò)程和能力分析過(guò)程分開(kāi)進(jìn)行,也就是將影響結(jié)構(gòu)地震需求和抗震能力的不確定性分開(kāi)處理,從而簡(jiǎn)化結(jié)構(gòu)地震易損性分析的過(guò)程.

    3 參數(shù)對(duì)易損性函數(shù)的影響

    橋梁地震易損性分析往往主要考慮地震波不確定性的影響,其他不確定性的影響則是通過(guò)對(duì)易損性函數(shù)中的對(duì)數(shù)標(biāo)準(zhǔn)差進(jìn)行修正來(lái)考慮.然而,橋梁建模相關(guān)不確定性也會(huì)對(duì)概率地震需求分析有影響,這種影響不僅體現(xiàn)在對(duì)數(shù)標(biāo)準(zhǔn)差β上,還有可能體現(xiàn)在結(jié)構(gòu)地震響應(yīng)均值θ上.易損性函數(shù)中均值θ和對(duì)數(shù)標(biāo)準(zhǔn)差β的變化都會(huì)影響最終的結(jié)構(gòu)地震易損性曲線(xiàn).因此,本節(jié)將結(jié)合通常的易損性函數(shù)(式(1)),通過(guò)參數(shù)分析來(lái)研究參數(shù)變化對(duì)地震易損性曲線(xiàn)的影響,以掌握易損性函數(shù)兩個(gè)參數(shù)的敏感性變化規(guī)律.

    3.1 均值變化的影響

    首先考慮對(duì)數(shù)標(biāo)準(zhǔn)差β=0.5時(shí),研究均值θ的取值變化對(duì)損傷概率密度函數(shù)以及易損性曲線(xiàn)的影響.當(dāng)θ從0.2 g增大到1.0 g時(shí),對(duì)應(yīng)的易損性函數(shù)以及相應(yīng)的概率密度函數(shù)的變化情況分別如圖2(a)和圖2(b)所示.由圖可知,地震易損性函數(shù)中的均值θ的理論意義是:當(dāng)結(jié)構(gòu)或構(gòu)件的失效超越概率達(dá)到50%時(shí)所對(duì)應(yīng)的地震動(dòng)強(qiáng)度值.由圖2(a)可知,均值θ越大時(shí),在相同的地震動(dòng)作用下結(jié)構(gòu)或構(gòu)件失效的概率越小.正是因?yàn)槿绱?,Padgett等[16]在研究中直接采用均值θ來(lái)衡量不同抗震加固措施對(duì)橋梁地震易損性的影響大小.由圖2(b)可知,由于β=0.5保持不變,其概率密度分布形式幾乎不變,只是其峰值點(diǎn)對(duì)應(yīng)的PGA隨均值θ的增大而逐漸向右移動(dòng).

    3.2 對(duì)數(shù)標(biāo)準(zhǔn)差變化的影響

    保持均值θ=0.5 g不變以研究對(duì)數(shù)標(biāo)準(zhǔn)差β的變化對(duì)損傷概率密度函數(shù)以及易損性曲線(xiàn)的影響.當(dāng)對(duì)數(shù)標(biāo)準(zhǔn)差β從0.2增大到1.0時(shí),易損性函數(shù)以及相應(yīng)的概率密度函數(shù)的變化情況分別如圖3(a)和3(b)所示.

    由圖3可知,由于均值θ=0.5 g保持不變,結(jié)構(gòu)或構(gòu)件失效超越概率為50%時(shí)所對(duì)應(yīng)的PGA值以及概率密度函數(shù)的峰值對(duì)應(yīng)的PGA值仍為0.5 g不變.由圖3(a)可知,當(dāng)?shù)卣饎?dòng)強(qiáng)度小于0.5 g時(shí),失效概率隨β的增加而增大;當(dāng)?shù)卣饎?dòng)強(qiáng)度大于0.5 g時(shí),失效概率則隨β的增大而減小.總的來(lái)說(shuō),地震易損性曲線(xiàn)隨β的增加會(huì)繞50%的失效概率點(diǎn)發(fā)生順時(shí)針?lè)较蜣D(zhuǎn)動(dòng).由圖3(b)可知,隨著對(duì)數(shù)標(biāo)準(zhǔn)差β的增加,概率密度函數(shù)由尖聳逐漸變得平緩,換句話(huà)說(shuō),PGA的取值集中于均值θ=0.5 g的程度越來(lái)越差,離散性越來(lái)越大.如果以θ=0.5 g作為橋梁地震作用下易損性大小預(yù)測(cè)值,則β反映了這種預(yù)測(cè)的準(zhǔn)確性.因此對(duì)數(shù)標(biāo)準(zhǔn)差β的理論意義是:概率性地震分析中結(jié)構(gòu)或構(gòu)件地震響應(yīng)預(yù)測(cè)均值的不確定性.

    3.3 兩參數(shù)同時(shí)變化的影響

    以上單獨(dú)分析了單個(gè)參數(shù)變化對(duì)易損性函數(shù)的影響.在實(shí)際中某些不確定性參數(shù)的變化(如選擇地震波集的變化)會(huì)同時(shí)導(dǎo)致均值θ和對(duì)數(shù)標(biāo)準(zhǔn)差β發(fā)生變化.因此,研究?jī)蓚€(gè)參數(shù)同時(shí)變化對(duì)易損性曲線(xiàn)的影響更貼近實(shí)際情況.圖4(a)和4(b)為兩個(gè)參數(shù)同時(shí)變化對(duì)易損性曲線(xiàn)和概率密度函數(shù)的影響.

    如圖4所示,當(dāng)均值θ和對(duì)數(shù)標(biāo)準(zhǔn)差β同時(shí)變化時(shí),易損性曲線(xiàn)的變化呈現(xiàn)多種變化形式.但總的來(lái)說(shuō),均值θ的變化更容易影響易損性曲線(xiàn)失效概率的大小,直接反映結(jié)構(gòu)抗震性能評(píng)估的結(jié)果;而β的變化則主要反映了該性能評(píng)估結(jié)果的準(zhǔn)確性,β越小則結(jié)構(gòu)抗震性能預(yù)測(cè)的精度越高.

    4 結(jié)論

    (1) 既往研究中的不確定性取值主要是根據(jù)大量不同建筑結(jié)構(gòu)的地震易損性函數(shù)統(tǒng)計(jì)出來(lái)的,為促進(jìn)我國(guó)橋梁地震易損性分析的發(fā)展和應(yīng)用,需要對(duì)橋梁建模、地震波和損傷指標(biāo)三個(gè)方面的不確定性分別展開(kāi)研究.

    (2) 云圖法易損性函數(shù)能將地震需求分析和能力分析過(guò)程分開(kāi)進(jìn)行,從而將影響結(jié)構(gòu)地震需求和抗震能力的不確定性分開(kāi)處理,與其他地震易損性函數(shù)相比,該方法在同時(shí)處理多種不確定性時(shí)表現(xiàn)出一定的優(yōu)勢(shì).

    (3) 地震易損性函數(shù)中均值θ,表示結(jié)構(gòu)或構(gòu)件的失效超越概率達(dá)50%時(shí)所對(duì)應(yīng)的地震動(dòng)強(qiáng)度值,θ的變化對(duì)易損性失效概率的影響較大,能直接反映結(jié)構(gòu)抗震性能評(píng)估的結(jié)果;對(duì)數(shù)標(biāo)準(zhǔn)差β表示概率性地震分析中對(duì)結(jié)構(gòu)或構(gòu)件地震響應(yīng)預(yù)測(cè)均值的不確定性,能反映結(jié)構(gòu)抗震性能評(píng)估結(jié)果的準(zhǔn)確性,β越小則結(jié)構(gòu)抗震性能預(yù)測(cè)的精度越高。

    猜你喜歡
    易損性不確定性抗震
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    關(guān)于房建結(jié)構(gòu)抗震設(shè)計(jì)的思考
    基于IDA的預(yù)應(yīng)力混凝土連續(xù)梁橋易損性分析
    古建抗震有絕招
    英鎊或繼續(xù)面臨不確定性風(fēng)險(xiǎn)
    抗震建筑材料簡(jiǎn)介
    上海建材(2018年6期)2018-03-21 06:24:28
    基于PSDM和IDA法的深水隔震橋梁地震易損性分析比較
    具有不可測(cè)動(dòng)態(tài)不確定性非線(xiàn)性系統(tǒng)的控制
    他成了抗震搶險(xiǎn)的先鋒
    基于性能的FRP加固RC框架結(jié)構(gòu)地震易損性分析
    搡老岳熟女国产| www.色视频.com| 欧美日韩精品成人综合77777| 亚洲自拍偷在线| 草草在线视频免费看| 97超视频在线观看视频| 欧美日韩乱码在线| 麻豆国产av国片精品| 国产v大片淫在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 又爽又黄a免费视频| 亚洲丝袜综合中文字幕| 久久久久精品国产欧美久久久| 观看免费一级毛片| 看十八女毛片水多多多| 一进一出抽搐gif免费好疼| 精品人妻视频免费看| 国产精品爽爽va在线观看网站| 国产欧美日韩一区二区精品| 久久精品国产亚洲av天美| 国产精品乱码一区二三区的特点| 日韩精品有码人妻一区| а√天堂www在线а√下载| 色视频www国产| 亚洲丝袜综合中文字幕| 欧美一级a爱片免费观看看| 亚洲av成人精品一区久久| 欧美最黄视频在线播放免费| 国产伦一二天堂av在线观看| 欧美+亚洲+日韩+国产| 波多野结衣高清无吗| 国产午夜精品久久久久久一区二区三区 | 欧美性感艳星| 国产精华一区二区三区| 99热6这里只有精品| 国产亚洲精品久久久com| 非洲黑人性xxxx精品又粗又长| 中文字幕免费在线视频6| 国产色婷婷99| 婷婷亚洲欧美| 精品人妻视频免费看| 国产男人的电影天堂91| 97热精品久久久久久| 国产美女午夜福利| 国产av一区在线观看免费| 久久亚洲国产成人精品v| 日本黄色视频三级网站网址| 亚洲成人精品中文字幕电影| 亚洲18禁久久av| 91久久精品国产一区二区成人| 激情 狠狠 欧美| 一本久久中文字幕| 国产v大片淫在线免费观看| 欧美高清性xxxxhd video| 国产av一区在线观看免费| 小蜜桃在线观看免费完整版高清| 老司机午夜福利在线观看视频| 69av精品久久久久久| 免费高清视频大片| 3wmmmm亚洲av在线观看| av国产免费在线观看| 色噜噜av男人的天堂激情| 九九久久精品国产亚洲av麻豆| 久久久成人免费电影| 日产精品乱码卡一卡2卡三| 99国产极品粉嫩在线观看| 国产一区二区三区av在线 | 高清午夜精品一区二区三区 | 我的老师免费观看完整版| 亚洲成人中文字幕在线播放| 女生性感内裤真人,穿戴方法视频| 国模一区二区三区四区视频| 啦啦啦啦在线视频资源| 国产欧美日韩精品一区二区| 欧美最新免费一区二区三区| 午夜影院日韩av| 国产成年人精品一区二区| 亚洲熟妇熟女久久| 亚洲最大成人av| 给我免费播放毛片高清在线观看| 精品午夜福利视频在线观看一区| 天堂av国产一区二区熟女人妻| 在线a可以看的网站| 久久人人精品亚洲av| 亚洲av成人av| 在线观看午夜福利视频| 一本一本综合久久| 亚洲人成网站在线播| 久久久欧美国产精品| 国产一级毛片七仙女欲春2| 一区福利在线观看| 深爱激情五月婷婷| 老司机影院成人| 欧美另类亚洲清纯唯美| 亚洲人成网站在线观看播放| 小说图片视频综合网站| 美女大奶头视频| 九九热线精品视视频播放| 简卡轻食公司| av在线老鸭窝| 网址你懂的国产日韩在线| 国产精品福利在线免费观看| 91在线观看av| 两个人的视频大全免费| 久久久久久伊人网av| 亚洲av熟女| 97超级碰碰碰精品色视频在线观看| 联通29元200g的流量卡| 午夜福利在线观看免费完整高清在 | 久久久久免费精品人妻一区二区| 国产爱豆传媒在线观看| 97超视频在线观看视频| 欧美成人a在线观看| 成人二区视频| 免费av观看视频| ponron亚洲| 日韩欧美三级三区| 国产精品久久视频播放| 午夜福利成人在线免费观看| 伦精品一区二区三区| 日日啪夜夜撸| 国产精品三级大全| 国内精品宾馆在线| 国产 一区 欧美 日韩| 欧美最新免费一区二区三区| 在线观看一区二区三区| 综合色丁香网| 97热精品久久久久久| 亚洲国产精品久久男人天堂| 国产精品野战在线观看| 国内精品美女久久久久久| 久久精品国产清高在天天线| 欧美zozozo另类| 国产精品伦人一区二区| h日本视频在线播放| 久久久久久久久久黄片| 国产人妻一区二区三区在| 91在线观看av| 在线观看一区二区三区| 亚洲欧美清纯卡通| 成人欧美大片| 亚洲精品一区av在线观看| 欧美在线一区亚洲| 国产av一区在线观看免费| 色哟哟哟哟哟哟| or卡值多少钱| 亚洲欧美精品综合久久99| 变态另类成人亚洲欧美熟女| 国产av不卡久久| 麻豆av噜噜一区二区三区| 欧美不卡视频在线免费观看| 国产精品国产三级国产av玫瑰| 国产亚洲av嫩草精品影院| 女人被狂操c到高潮| 狠狠狠狠99中文字幕| 在线免费观看的www视频| 国产高清有码在线观看视频| 久久久精品欧美日韩精品| 一个人观看的视频www高清免费观看| av女优亚洲男人天堂| 国产单亲对白刺激| 日日摸夜夜添夜夜添小说| 伦理电影大哥的女人| 国产av在哪里看| 精品久久久久久久久av| 最近手机中文字幕大全| 日本在线视频免费播放| 伦理电影大哥的女人| 国产一区二区在线观看日韩| 校园人妻丝袜中文字幕| 久久99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 免费在线观看影片大全网站| 香蕉av资源在线| 欧美成人a在线观看| 成年av动漫网址| 亚洲欧美日韩东京热| 成人三级黄色视频| 少妇裸体淫交视频免费看高清| 亚洲精品456在线播放app| 国语自产精品视频在线第100页| 国产精品,欧美在线| 好男人在线观看高清免费视频| 久久久a久久爽久久v久久| 日韩欧美 国产精品| 免费观看人在逋| 级片在线观看| 国产av在哪里看| 久久久久国内视频| 久久久久久久久中文| 在线观看一区二区三区| 久久久国产成人免费| 特级一级黄色大片| 亚州av有码| 97超级碰碰碰精品色视频在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品三级大全| 成人亚洲精品av一区二区| 寂寞人妻少妇视频99o| 一级毛片久久久久久久久女| 丝袜喷水一区| 国产午夜精品久久久久久一区二区三区 | 一个人观看的视频www高清免费观看| 99精品在免费线老司机午夜| 精华霜和精华液先用哪个| 99久久精品国产国产毛片| 亚洲性久久影院| 久久99热这里只有精品18| 亚洲成人久久爱视频| 国语自产精品视频在线第100页| 久久久久久国产a免费观看| 免费黄网站久久成人精品| 一进一出抽搐gif免费好疼| 日日摸夜夜添夜夜爱| 大又大粗又爽又黄少妇毛片口| 精华霜和精华液先用哪个| 最近在线观看免费完整版| 免费看光身美女| 欧美精品国产亚洲| 国产精品伦人一区二区| 精品一区二区三区人妻视频| 日韩成人av中文字幕在线观看 | 成人特级av手机在线观看| 成年av动漫网址| 大型黄色视频在线免费观看| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 91av网一区二区| 好男人在线观看高清免费视频| 久久鲁丝午夜福利片| 精品久久国产蜜桃| 国内久久婷婷六月综合欲色啪| 97碰自拍视频| 一级毛片我不卡| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 成人二区视频| 成人毛片a级毛片在线播放| 国产 一区精品| 免费观看的影片在线观看| av在线亚洲专区| 国产精品人妻久久久影院| 久久久国产成人精品二区| av在线蜜桃| 日日干狠狠操夜夜爽| 免费一级毛片在线播放高清视频| 日本 av在线| 三级毛片av免费| 全区人妻精品视频| 综合色丁香网| 国产成人a区在线观看| 日本与韩国留学比较| 男人的好看免费观看在线视频| 欧美日韩国产亚洲二区| 国产精品一区二区三区四区久久| 欧美成人免费av一区二区三区| 国产高潮美女av| 欧美xxxx黑人xx丫x性爽| 久久久久免费精品人妻一区二区| 日本精品一区二区三区蜜桃| 晚上一个人看的免费电影| 日韩av不卡免费在线播放| 亚洲精品影视一区二区三区av| 免费av观看视频| 天堂动漫精品| 中文字幕免费在线视频6| 国产精品久久久久久亚洲av鲁大| 美女cb高潮喷水在线观看| 伊人久久精品亚洲午夜| 麻豆国产97在线/欧美| 久久久久久久久久久丰满| 国产一区二区三区av在线 | 国产av麻豆久久久久久久| 久久精品国产亚洲av香蕉五月| 最新在线观看一区二区三区| 老熟妇仑乱视频hdxx| 亚洲精品久久国产高清桃花| 久久久久久久久久成人| 床上黄色一级片| 老师上课跳d突然被开到最大视频| 亚洲精品成人久久久久久| 热99re8久久精品国产| 午夜福利在线观看免费完整高清在 | 在线观看午夜福利视频| 不卡视频在线观看欧美| 国产精品日韩av在线免费观看| 老女人水多毛片| 五月伊人婷婷丁香| 日本撒尿小便嘘嘘汇集6| 悠悠久久av| 国产精品久久视频播放| 人妻制服诱惑在线中文字幕| 黄色视频,在线免费观看| 可以在线观看毛片的网站| 免费电影在线观看免费观看| 国产91av在线免费观看| 人人妻人人澡欧美一区二区| 国产淫片久久久久久久久| 三级男女做爰猛烈吃奶摸视频| av在线播放精品| av在线亚洲专区| 国产黄a三级三级三级人| 精品人妻一区二区三区麻豆 | 麻豆乱淫一区二区| 日本黄大片高清| 日韩高清综合在线| 大香蕉久久网| 亚洲性夜色夜夜综合| 日韩大尺度精品在线看网址| 久久99热这里只有精品18| 成熟少妇高潮喷水视频| 国产一区二区三区av在线 | 国产精品电影一区二区三区| 久久99热这里只有精品18| 如何舔出高潮| 免费看a级黄色片| 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 久久久久九九精品影院| 婷婷亚洲欧美| 校园人妻丝袜中文字幕| 少妇的逼水好多| 免费电影在线观看免费观看| 欧美另类亚洲清纯唯美| 国产三级在线视频| 秋霞在线观看毛片| ponron亚洲| 国产片特级美女逼逼视频| 国产老妇女一区| 免费搜索国产男女视频| 少妇的逼水好多| 免费av毛片视频| 午夜亚洲福利在线播放| 国产欧美日韩精品亚洲av| 在现免费观看毛片| 人妻丰满熟妇av一区二区三区| 亚洲国产精品sss在线观看| 一本一本综合久久| 欧美人与善性xxx| 日韩精品有码人妻一区| 午夜福利视频1000在线观看| 亚洲中文日韩欧美视频| 日本色播在线视频| 亚洲va在线va天堂va国产| 18+在线观看网站| 成人国产麻豆网| 1024手机看黄色片| 熟妇人妻久久中文字幕3abv| 天堂动漫精品| 国产v大片淫在线免费观看| 亚洲最大成人中文| 久久久久久久久大av| 久久久精品94久久精品| 97热精品久久久久久| 日韩一区二区视频免费看| 麻豆久久精品国产亚洲av| a级毛片免费高清观看在线播放| 简卡轻食公司| 又黄又爽又刺激的免费视频.| 亚洲精品日韩在线中文字幕 | 99久久精品热视频| 亚洲国产精品成人久久小说 | 久久久久免费精品人妻一区二区| 男人舔奶头视频| 女人十人毛片免费观看3o分钟| 日本一本二区三区精品| 2021天堂中文幕一二区在线观| 中国美女看黄片| 免费av观看视频| 国产私拍福利视频在线观看| 麻豆乱淫一区二区| 在线观看免费视频日本深夜| 舔av片在线| 级片在线观看| 亚洲自拍偷在线| 欧美色视频一区免费| 黄色一级大片看看| 久久久久久久午夜电影| 日本三级黄在线观看| 亚洲精品在线观看二区| 午夜老司机福利剧场| 91狼人影院| 十八禁网站免费在线| 国产精品美女特级片免费视频播放器| 精品久久久久久久久久免费视频| 寂寞人妻少妇视频99o| 成年女人永久免费观看视频| 女的被弄到高潮叫床怎么办| 亚洲中文字幕日韩| 天美传媒精品一区二区| 亚洲欧美日韩高清在线视频| 亚洲精品成人久久久久久| 乱系列少妇在线播放| 国内精品久久久久精免费| 自拍偷自拍亚洲精品老妇| 一级毛片我不卡| 国产精品99久久久久久久久| 国产男人的电影天堂91| 午夜精品国产一区二区电影 | 免费观看的影片在线观看| 白带黄色成豆腐渣| 日本免费a在线| 中文字幕av在线有码专区| 亚洲一级一片aⅴ在线观看| 精品99又大又爽又粗少妇毛片| 最好的美女福利视频网| 可以在线观看的亚洲视频| 婷婷精品国产亚洲av| 亚洲七黄色美女视频| 中文资源天堂在线| 九色成人免费人妻av| 午夜久久久久精精品| 在线观看一区二区三区| 女人被狂操c到高潮| 国产三级中文精品| 国产亚洲av嫩草精品影院| 伦精品一区二区三区| 亚洲熟妇熟女久久| 亚洲美女搞黄在线观看 | 激情 狠狠 欧美| 国产乱人偷精品视频| 亚洲精品一区av在线观看| 中国美女看黄片| 午夜日韩欧美国产| 97超视频在线观看视频| 精品99又大又爽又粗少妇毛片| 99国产极品粉嫩在线观看| 午夜福利视频1000在线观看| 天天一区二区日本电影三级| 亚洲国产欧美人成| 国产 一区 欧美 日韩| avwww免费| 日本黄色片子视频| 一区二区三区高清视频在线| 一级毛片我不卡| 欧美不卡视频在线免费观看| 欧美性感艳星| 成人综合一区亚洲| 亚洲一级一片aⅴ在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲欧美日韩无卡精品| 韩国av在线不卡| 天天躁夜夜躁狠狠久久av| 99在线视频只有这里精品首页| 欧美最新免费一区二区三区| 夜夜爽天天搞| 春色校园在线视频观看| 国产精品一区二区三区四区免费观看 | 综合色丁香网| 国产又黄又爽又无遮挡在线| 婷婷六月久久综合丁香| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 乱系列少妇在线播放| 免费在线观看影片大全网站| 久久久久九九精品影院| 又黄又爽又免费观看的视频| 欧美一区二区亚洲| 成人欧美大片| av免费在线看不卡| 午夜精品国产一区二区电影 | 一本精品99久久精品77| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲欧美98| 最新中文字幕久久久久| 国产老妇女一区| 深爱激情五月婷婷| 小蜜桃在线观看免费完整版高清| 搡女人真爽免费视频火全软件 | 日本在线视频免费播放| 国产一区二区在线观看日韩| 非洲黑人性xxxx精品又粗又长| 国产极品精品免费视频能看的| 美女内射精品一级片tv| 欧美激情国产日韩精品一区| 亚洲精品国产成人久久av| 国产色婷婷99| 亚洲av免费在线观看| 一本精品99久久精品77| 久久99热这里只有精品18| 晚上一个人看的免费电影| av中文乱码字幕在线| 毛片一级片免费看久久久久| 一个人看视频在线观看www免费| 国产一区亚洲一区在线观看| 国产高清视频在线观看网站| 蜜桃久久精品国产亚洲av| 亚洲自拍偷在线| 国产成人aa在线观看| 中文亚洲av片在线观看爽| 日日干狠狠操夜夜爽| 国产伦一二天堂av在线观看| 国产黄a三级三级三级人| 淫妇啪啪啪对白视频| 亚洲综合色惰| 狂野欧美白嫩少妇大欣赏| 成人特级黄色片久久久久久久| 日韩欧美精品免费久久| 国产老妇女一区| 女人十人毛片免费观看3o分钟| 亚洲人成网站在线观看播放| 日本精品一区二区三区蜜桃| 内地一区二区视频在线| av在线观看视频网站免费| 99热全是精品| 欧美+亚洲+日韩+国产| 国产精品美女特级片免费视频播放器| 三级国产精品欧美在线观看| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 精品欧美国产一区二区三| 国产熟女欧美一区二区| 亚洲精品国产av成人精品 | 嫩草影院精品99| 尾随美女入室| 俄罗斯特黄特色一大片| 日韩欧美在线乱码| 男女边吃奶边做爰视频| 小蜜桃在线观看免费完整版高清| 日韩强制内射视频| 久久久久免费精品人妻一区二区| 国产精品女同一区二区软件| 熟妇人妻久久中文字幕3abv| 午夜亚洲福利在线播放| 婷婷亚洲欧美| 午夜福利在线观看吧| 在线观看免费视频日本深夜| 午夜老司机福利剧场| 高清午夜精品一区二区三区 | 久久久午夜欧美精品| 99热全是精品| 99久久久亚洲精品蜜臀av| 一进一出抽搐动态| 一个人看视频在线观看www免费| 国产成人精品久久久久久| 美女内射精品一级片tv| 亚洲第一区二区三区不卡| 亚洲人与动物交配视频| 无遮挡黄片免费观看| 成熟少妇高潮喷水视频| 日日撸夜夜添| 亚洲精品乱码久久久v下载方式| 99热网站在线观看| 大又大粗又爽又黄少妇毛片口| 综合色av麻豆| 伦理电影大哥的女人| 在线观看午夜福利视频| 国产一级毛片七仙女欲春2| 一进一出抽搐动态| 深夜a级毛片| 一卡2卡三卡四卡精品乱码亚洲| 搡老岳熟女国产| 晚上一个人看的免费电影| 成人国产麻豆网| 亚洲在线观看片| 一级毛片电影观看 | 深夜精品福利| 精品一区二区三区av网在线观看| 国产黄a三级三级三级人| 亚洲丝袜综合中文字幕| 九九热线精品视视频播放| 特级一级黄色大片| 少妇猛男粗大的猛烈进出视频 | aaaaa片日本免费| 国产精品福利在线免费观看| 女的被弄到高潮叫床怎么办| 最新中文字幕久久久久| 中文在线观看免费www的网站| 国产精品一及| 乱码一卡2卡4卡精品| 日韩三级伦理在线观看| 内地一区二区视频在线| 99精品在免费线老司机午夜| 久久精品国产清高在天天线| 国产 一区 欧美 日韩| 国产精品一区二区三区四区久久| 欧美最黄视频在线播放免费| 97在线视频观看| 亚洲国产日韩欧美精品在线观看| or卡值多少钱| 国产伦在线观看视频一区| 亚洲av第一区精品v没综合| 亚洲国产精品合色在线| 久久6这里有精品| 婷婷亚洲欧美| 美女cb高潮喷水在线观看| 日本免费一区二区三区高清不卡| 精品午夜福利视频在线观看一区| 国国产精品蜜臀av免费| 日韩精品中文字幕看吧| 最后的刺客免费高清国语| 日本黄大片高清| 国产一区二区激情短视频| 免费看av在线观看网站| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类| 国产精品福利在线免费观看| 亚洲久久久久久中文字幕| 国产精品国产三级国产av玫瑰| 天堂av国产一区二区熟女人妻| 日韩高清综合在线| av国产免费在线观看| 中文字幕av在线有码专区| 毛片女人毛片| 亚洲美女视频黄频| 中文资源天堂在线| 国产片特级美女逼逼视频| 菩萨蛮人人尽说江南好唐韦庄 | 九色成人免费人妻av| 免费看av在线观看网站| 婷婷精品国产亚洲av| 给我免费播放毛片高清在线观看| 精品熟女少妇av免费看| 99在线视频只有这里精品首页| 白带黄色成豆腐渣|