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

    超高層混合結(jié)構(gòu)增量動(dòng)力分析的地震動(dòng)強(qiáng)度指標(biāo)

    2021-07-11 17:13:42鄧暉付志旭傅冠翔潘建榮
    土木建筑與環(huán)境工程 2021年3期
    關(guān)鍵詞:調(diào)幅位數(shù)震動(dòng)

    鄧暉 付志旭 傅冠翔 潘建榮

    摘 要:在進(jìn)行超高層結(jié)構(gòu)設(shè)計(jì)時(shí),地震動(dòng)強(qiáng)度指標(biāo)是指導(dǎo)結(jié)構(gòu)抗震設(shè)計(jì)的重要依據(jù),對(duì)結(jié)構(gòu)的安全至關(guān)重要,但目前并沒(méi)有統(tǒng)一的結(jié)論,因此,對(duì)超高層結(jié)構(gòu)進(jìn)行地震動(dòng)強(qiáng)度指標(biāo)的研究具有重要意義。在綜述建筑結(jié)構(gòu)增量動(dòng)力分析地震動(dòng)強(qiáng)度指標(biāo)有效性研究現(xiàn)狀的基礎(chǔ)上,提出適用于超高層建筑結(jié)構(gòu)增量動(dòng)力分析的地震動(dòng)強(qiáng)度指標(biāo)Sa%,地震動(dòng)強(qiáng)度指標(biāo)Sa%是考慮結(jié)構(gòu)前a%振型質(zhì)量參與系數(shù)的自振周期所對(duì)應(yīng)譜加速度的幾何平均值;為了驗(yàn)證地震動(dòng)強(qiáng)度指標(biāo)Sa%的有效性,以實(shí)際超高層為例,使用PERFORM-3D彈塑性分析軟件建立非線性模型并進(jìn)行增量動(dòng)力分析,根據(jù)各種地震動(dòng)強(qiáng)度指標(biāo)的θmax和Vmax結(jié)果進(jìn)行離散性分析,發(fā)現(xiàn)地震動(dòng)強(qiáng)度指標(biāo)Sa%用于超高層結(jié)構(gòu)增量動(dòng)力分析時(shí)能有效減少θmax和Vmax數(shù)據(jù)的離散性,從而提高計(jì)算數(shù)據(jù)的有效性,為實(shí)際超高層結(jié)構(gòu)的應(yīng)用提供參考。

    關(guān)鍵詞:增量動(dòng)力分析;超高層;地震動(dòng)強(qiáng)度指標(biāo);離散性

    中圖分類號(hào):TU973 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2096-6717(2021)03-0109-10

    Abstract: In design of super high-rise structure, the ground motion intensity index guides the seismic design of the structure and is very important for safety of the structure,but there is no unified conclusion on the ground motion intensity index of super high-rise mixed structure. Therefore, it is of great significance to study the ground motion intensity index of the super high-rise structure. Based on the review of the current research status of the effectiveness of seismic intensity index in the incremental dynamic analysis of building structures, this article proposes a seismic intensity index a%, which is suitable for the incremental dynamic analysis of super high-rise building structure. The ground motion intensity index Sa% is the geometric average of the spectral acceleration corresponding to the natural vibration period of a% mode mass participation coefficient before considering the translational direction of the structure. In order to verify the effectiveness of the seismic intensity index Sa%, PERFORM-3D elastoplastic analysis software was used to establish a nonlinear model of an actual super high-rise building and perform incremental dynamic analysis. Based on the discrete analysis of θmax and Vmax results of various ground motion intensity indexes, it is concluded that Sa% of ground motion intensity index proposed in this paper can effectively reduce the discreteness of θmax and Vmax when applied to incremental dynamic analysis of super-high-rise structures, which improves the validity of calculated data and provides reference for the application of actual super high-rise structures.

    Keywords: incremental dynamic analysis; super high-rise building; seismic intensity index; dispersion

    地震動(dòng)強(qiáng)度指標(biāo)一直是眾多學(xué)者研究的熱點(diǎn)之一,劉恢先[1]認(rèn)為選擇地震動(dòng)強(qiáng)度指標(biāo)的方向主要有兩個(gè):一個(gè)是地震動(dòng)地面運(yùn)動(dòng)指標(biāo),如PGA、PGV等;另一個(gè)是采取地震動(dòng)的反應(yīng)譜值,如結(jié)構(gòu)第一自振周期對(duì)應(yīng)的譜加速度。眾多學(xué)者通過(guò)研究提出的地震動(dòng)強(qiáng)度指標(biāo)目前已有30種之多[2]。增量動(dòng)力分析[3](Incremental Dynamic Analysis, IDA)是近年來(lái)受到眾多學(xué)者關(guān)注的方法,該方法相比于靜力彈塑性方法和動(dòng)力彈塑性方法具有明顯的優(yōu)勢(shì)。但因?yàn)槠溆?jì)算量大和研究積累的成果較少,所以,還未能成為抗震分析的主流方法,目前主要被廣泛應(yīng)用于結(jié)構(gòu)基于性能的抗震評(píng)估中[4-6]。但隨著計(jì)算機(jī)科學(xué)技術(shù)的飛速發(fā)展,IDA將是未來(lái)建筑結(jié)構(gòu)彈塑性分析的重要發(fā)展方向。

    目前眾多學(xué)者在進(jìn)行增量動(dòng)力分析時(shí),都選用PGA或者Sa(T1,ξ)作為地震動(dòng)強(qiáng)度指標(biāo)進(jìn)行調(diào)幅,因?yàn)檫@兩種指標(biāo)概念清晰,且計(jì)算簡(jiǎn)單,滿足了比例調(diào)整魯棒性、充分性和地震危險(xiǎn)性可計(jì)算性的要求,但在有效性方面,兩種指標(biāo)的適用性一直存在較大爭(zhēng)議。對(duì)于以第一振型為主導(dǎo)的低層、多層結(jié)構(gòu),用Sa(T1,ξ)進(jìn)行調(diào)幅計(jì)算得到的結(jié)果離散性較小[7],這一結(jié)論已獲得了較一致的認(rèn)同。但對(duì)于高層或者超高層建筑則尚未有一致的結(jié)論。Guan等[8]以框架核心筒為研究對(duì)象,對(duì)地震動(dòng)強(qiáng)度指標(biāo)進(jìn)行研究,結(jié)果表明,PGV作為長(zhǎng)周期地震動(dòng)強(qiáng)度指標(biāo)時(shí)效果較好,而PGA不適用于長(zhǎng)周期結(jié)構(gòu)。盧嘯等[9]考慮到超高層結(jié)構(gòu)高階振型對(duì)結(jié)構(gòu)有較大影響的特點(diǎn),基于彎剪耦合簡(jiǎn)化模型,提出適用于超高層結(jié)構(gòu)的地震動(dòng)強(qiáng)度指標(biāo)Sa,并通過(guò)兩棟高層結(jié)構(gòu)的倒塌實(shí)例分析驗(yàn)證了該指標(biāo)的相關(guān)性和離散程度。周穎等[10]對(duì)某型鋼混凝土框架鋼筋混凝土核心筒超高層建筑進(jìn)行了IDA分析,對(duì)比不同地震動(dòng)參數(shù)IDA曲線的條件對(duì)數(shù)標(biāo)準(zhǔn)差平均值,得到對(duì)于高層建筑,考慮前3階譜值參數(shù)的S123離散型最小、PGV次之、PGA離散型較大的結(jié)論。候喆辰[11]對(duì)高層剪力墻結(jié)構(gòu)實(shí)際工程進(jìn)行增量動(dòng)力分析,對(duì)比分析結(jié)果的標(biāo)準(zhǔn)差,得到選用PGA作為地震動(dòng)強(qiáng)度指標(biāo)時(shí)計(jì)算結(jié)果的離散性比選用PGV要小的結(jié)論。郝永蕾[12]對(duì)某30層型鋼混凝土框架混凝土核心筒混合結(jié)構(gòu)進(jìn)行IDA分析,選擇分位數(shù)曲線法和條件對(duì)數(shù)標(biāo)準(zhǔn)差對(duì)3個(gè)參數(shù)所得結(jié)果進(jìn)行離散性對(duì)比分析,得到選用PGA作為地震動(dòng)參數(shù)時(shí)離散性最小,PGV離散性稍大,Sa(T1,ξ)離散性較大的結(jié)論。由此可見(jiàn),目前對(duì)于高層和超高層建筑的地震動(dòng)強(qiáng)度指標(biāo)并沒(méi)有較一致的看法,筆者提出適用于超高層建筑結(jié)構(gòu)增量動(dòng)力分析的地震動(dòng)強(qiáng)度指標(biāo)Sa%,并通過(guò)實(shí)際工程證明其指標(biāo)的有效性。

    1 地震動(dòng)強(qiáng)度指標(biāo)及Sa%的提出

    由于超高層建筑結(jié)構(gòu)在設(shè)計(jì)時(shí)大多采用基于性能的抗震設(shè)計(jì)方法,具有較高的安全儲(chǔ)備,在大震作用下,結(jié)構(gòu)的主要抗側(cè)力構(gòu)件基本保持彈性或進(jìn)入塑性的程度不高[13-14],非彈性反應(yīng)譜是地震動(dòng)特性和結(jié)構(gòu)材料非線性相互耦合作用的結(jié)果,其中混合了結(jié)構(gòu)的非線性影響的成分[10],故筆者暫不考慮采用非彈性譜作為地震動(dòng)強(qiáng)度指標(biāo)?;谀壳暗卣饎?dòng)強(qiáng)度指標(biāo)的研究現(xiàn)狀和地震動(dòng)強(qiáng)度指標(biāo)在進(jìn)行增量動(dòng)力分析時(shí)需要滿足的特性,納入幾個(gè)目前常用的地震動(dòng)強(qiáng)度指標(biāo),并納入概念簡(jiǎn)單、使用方便、且能夠考慮高階振型影響從而適用于超高層結(jié)構(gòu)的地震動(dòng)強(qiáng)度指標(biāo)。選用已有的地震動(dòng)強(qiáng)度指標(biāo)如下:

    1)地震動(dòng)的峰值加速度PGA。

    2)地震動(dòng)的峰值速度PGV。

    3)阻尼比為4%時(shí),結(jié)構(gòu)第一自振周期對(duì)應(yīng)的譜加速度Sa(T1,4%)。

    4)阻尼比為4%時(shí),地震動(dòng)加速度譜峰值PSA。

    5)阻尼比為4%時(shí),地震動(dòng)速度譜峰值PSV。

    6)文獻(xiàn)[3]中提出的考慮高階振型參與的冪函數(shù)乘積形式的地震動(dòng)強(qiáng)度指標(biāo)。它以結(jié)構(gòu)主要周期對(duì)應(yīng)的彈性加速度反應(yīng)譜值為底數(shù),以各周期的振型參與質(zhì)量系數(shù)比為指數(shù)的各振型結(jié)果的乘積。

    當(dāng)考慮前2階振型時(shí),其表達(dá)式為

    當(dāng)考慮前3階振型時(shí),其表達(dá)式為

    式中:Sa(Ti)為第i周期對(duì)應(yīng)的譜加速度;n為所考慮的結(jié)構(gòu)平動(dòng)振型數(shù),與結(jié)構(gòu)的基本周期相關(guān);Sa的物理含義是結(jié)構(gòu)前n階平動(dòng)自振周期對(duì)應(yīng)的譜加速度的幾何平均值,反映了結(jié)構(gòu)地震響應(yīng)中參與顯著的結(jié)構(gòu)振型數(shù)量。

    在第6)、7)地震動(dòng)強(qiáng)度指標(biāo)的基礎(chǔ)上,提出適用于超高層結(jié)構(gòu)的地震動(dòng)強(qiáng)度指標(biāo)Sa%。該指標(biāo)物理含義為考慮結(jié)構(gòu)前a%振型質(zhì)量參與系數(shù)的自振周期對(duì)應(yīng)譜加速度的幾何平均值。

    提出此地震動(dòng)強(qiáng)度指標(biāo)的理由如下:

    1)S12和S123考慮高階振型參與影響的能力有限。對(duì)于許多復(fù)雜高層結(jié)構(gòu),若考慮更高階振型參與的指標(biāo)形式S1234和S12345,會(huì)因?yàn)楦唠A振型的質(zhì)量參與系數(shù)占比不高使該指標(biāo)數(shù)值變化很小,這會(huì)在一定程度上弱化高階振型對(duì)結(jié)構(gòu)的影響。但其依據(jù)質(zhì)量參與系數(shù)來(lái)確定參與振型數(shù)量的方法取得了一定的效果,所以,Sa%也以平動(dòng)方向的質(zhì)量參與系數(shù)來(lái)確定振型的參與數(shù)量。

    2)Sa考慮的結(jié)構(gòu)前n階振型主要由經(jīng)驗(yàn)公式所確定,該經(jīng)驗(yàn)公式由彎曲剪切耦合的二維連續(xù)化模型進(jìn)行時(shí)程分析得到,并非三維的超高層模型。且用兩個(gè)超高層模型進(jìn)行驗(yàn)證,模型A的驗(yàn)證效果較好,模型B中該指標(biāo)的效果不如PGV。綜上來(lái)看,該指標(biāo)確定振型數(shù)量的經(jīng)驗(yàn)公式存在局限性,因?yàn)樵摻?jīng)驗(yàn)公式只與結(jié)構(gòu)的第一自振周期有關(guān)。而即使是第一自振周期相近的超高層結(jié)構(gòu),高階振型的參與程度也因受到諸多因素的影響而有較大差別。但該指標(biāo)計(jì)算各自振周期對(duì)應(yīng)譜加速度的幾何平均值的計(jì)算方法相比于S12和S123能更有效地反映出結(jié)構(gòu)高階振型的參與。

    在分析上述兩種指標(biāo)優(yōu)點(diǎn)和不足的基礎(chǔ)上,提出以平動(dòng)方向前a%振型質(zhì)量參與系數(shù)為依據(jù),確定參與的自振周期數(shù),并計(jì)算對(duì)應(yīng)譜加速度的幾何平均值。

    2 地震動(dòng)強(qiáng)度指標(biāo)離散性分析方法

    筆者研究的內(nèi)容是使用不同地震動(dòng)強(qiáng)度指標(biāo)進(jìn)行IDA調(diào)幅計(jì)算,將計(jì)算結(jié)果得到的數(shù)據(jù)進(jìn)行離散性分析,此項(xiàng)研究主要涉及兩個(gè)問(wèn)題:一是如何將由PGA調(diào)幅計(jì)算得到的數(shù)據(jù)轉(zhuǎn)化成為各項(xiàng)地震動(dòng)指標(biāo)調(diào)幅計(jì)算下得到的結(jié)果,二是采用何種方法來(lái)評(píng)定數(shù)據(jù)的離散性。

    由于IDA計(jì)算需要大量的非線性時(shí)程分析,算例結(jié)構(gòu)復(fù)雜,每次時(shí)程分析耗時(shí)較長(zhǎng),且考慮的地震動(dòng)強(qiáng)度指標(biāo)較多,但地震動(dòng)強(qiáng)度指標(biāo)都具有比例調(diào)整魯棒性,可以將PGA進(jìn)行調(diào)幅后計(jì)算得到的結(jié)果轉(zhuǎn)化成等效于用其他指標(biāo)調(diào)幅后計(jì)算所得到的結(jié)果,PGA調(diào)幅值與各地震動(dòng)指標(biāo)下相對(duì)應(yīng)的調(diào)幅值的具體轉(zhuǎn)換步驟如下:1)計(jì)算各條地震動(dòng)作用下結(jié)構(gòu)各個(gè)地震動(dòng)指標(biāo)的強(qiáng)度值;2)計(jì)算出各條地震動(dòng)作用下其他地震動(dòng)強(qiáng)度指標(biāo)與PGA的比值&;3)由于各條地震動(dòng)PGA的調(diào)幅值是按照50、100、200、300、400、600、800、1 000 cm/s2進(jìn)行,所以,將&與PGA對(duì)應(yīng)的各個(gè)調(diào)幅值(50、100、200、300、400、600、800、1 000 cm/s2)相乘,最終得到該條地震動(dòng)作用下PGA調(diào)幅值相對(duì)應(yīng)的各個(gè)地震動(dòng)指標(biāo)的強(qiáng)度值。由上述內(nèi)容可知,只要計(jì)算出每條地震動(dòng)作用下各地震動(dòng)指標(biāo)與PGA的比值,將該比值與PGA相應(yīng)調(diào)幅值相乘,就可以把PGA調(diào)幅后計(jì)算得到的響應(yīng)結(jié)果相應(yīng)地分配給其他地震動(dòng)強(qiáng)度指標(biāo)。

    為了評(píng)定數(shù)據(jù)的離散型,采用分位數(shù)曲線法判斷IDA計(jì)算數(shù)據(jù)的離散程度,然后再計(jì)算IDA計(jì)算數(shù)據(jù)對(duì)于地震動(dòng)強(qiáng)度指標(biāo)的條件對(duì)數(shù)標(biāo)準(zhǔn)差平均值來(lái)直觀判斷數(shù)據(jù)的離散程度。分位數(shù)曲線法假定每條IDA曲線均服從對(duì)數(shù)正態(tài)分布,然后通過(guò)對(duì)數(shù)據(jù)百分位數(shù)的判定來(lái)衡量數(shù)據(jù)所在的相對(duì)位置。計(jì)算方法為計(jì)算不同DM值相對(duì)應(yīng)的不同IM值的均值μM和不同IM值的對(duì)數(shù)標(biāo)準(zhǔn)差δM,得到(DM,μM)、(DM,μM×exp(+δM))、(DM,μM×exp(-δM))3條曲線,分別為50%、84%、16%比例曲線,可知不同地震動(dòng)強(qiáng)度指標(biāo)的離散程度。為了更加直觀地看出數(shù)據(jù)的離散程度,計(jì)算每個(gè)θmax和Vmax的分位數(shù)曲線離散程度的誤差

    式中:Δ84%DM、Δ16%DM分別為84%分位數(shù)曲線和16%分位數(shù)曲線相對(duì)于50%分位數(shù)曲線的誤差;IM84%、IM50%、IM16%分別為每個(gè)DM值對(duì)應(yīng)下84%、50%、16%分位數(shù)曲線上對(duì)應(yīng)的IM值。計(jì)算出某地震動(dòng)指標(biāo)相對(duì)于θmax和Vmax的分位數(shù)曲線離散程度誤差后,取平均值即可得到離散程度平均誤差。

    同時(shí)也計(jì)算出每個(gè)地震動(dòng)指標(biāo)下IM對(duì)于DM的條件對(duì)數(shù)標(biāo)準(zhǔn)差,并計(jì)算出條件對(duì)數(shù)標(biāo)準(zhǔn)差的平均值,平均值越小代表離散程度越小,計(jì)算數(shù)據(jù)的有效性越好。

    3 算例

    3.1 結(jié)構(gòu)基本信息

    深圳匯德大廈項(xiàng)目位于深圳北站附近,包含2#、4#裙房和1#、3#主樓。其中,1#主樓結(jié)構(gòu)高度為248.75 m,地上總共58層,地下3層。項(xiàng)目所在地抗震設(shè)防烈度為7度,場(chǎng)地類別是II類,設(shè)計(jì)地震分組為第一組,場(chǎng)地特征周期是0.35 s,設(shè)計(jì)基本地震加速度值為0.10g。1#主樓是鋼框架鋼筋混凝土核心筒混合結(jié)構(gòu)體系,在12~16層、24~28層、35~39層、46~50層設(shè)置斜柱實(shí)現(xiàn)外框架相對(duì)于核心筒斜切角式的縮進(jìn),使建筑表現(xiàn)出旋轉(zhuǎn)向上的效果。

    匯德大廈1#主樓第一轉(zhuǎn)換區(qū)以下的外框架平面尺寸為45.1 m×45.1 m的正方形,核心筒尺寸為23.1 m×23.7 m的矩形,外框架與核心筒的邊界并不平行,通過(guò)4個(gè)轉(zhuǎn)換區(qū)外框斜切角的轉(zhuǎn)換后,外框平面尺寸變?yōu)?8.4 m×38.4 m,且外框架邊與核心筒邊平行,部分樓層結(jié)構(gòu)平面布置如圖1所示。核心筒的厚度和平面布置變化較多,布置方式如圖2所示,混凝土強(qiáng)度等級(jí)與厚度變化如表1所示。外框架采用鋼管混凝土柱,主要鋼材牌號(hào)為Q345B,斜柱轉(zhuǎn)換層附近的部分柱牌號(hào)采用了Q390GJC。梁構(gòu)件由于截面尺寸眾多,不再列舉。

    3.2 模型的建立

    綜合考慮軟件的計(jì)算精度和效率,最終選擇PERFORM-3D對(duì)結(jié)構(gòu)進(jìn)行增量動(dòng)力分析。對(duì)于梁柱單元的模擬,PERFORM-3D中可采用弦轉(zhuǎn)角模型、塑性鉸模型和塑性區(qū)模型來(lái)模擬梁柱單元。弦轉(zhuǎn)角模型雖然計(jì)算簡(jiǎn)單,但由于其反彎點(diǎn)在跨中,所以只適用于對(duì)稱的結(jié)構(gòu)。塑性鉸模型由彈性段和塑性鉸組成,塑性鉸集中于一點(diǎn),無(wú)法模擬塑性區(qū)的長(zhǎng)度,研究人員需對(duì)塑性鉸的特點(diǎn)有較深入的了解。塑性區(qū)模型將構(gòu)件分為彈性段和塑性段,可以模擬塑性區(qū)段的剛度退化,更加符合實(shí)際結(jié)構(gòu)的破壞狀態(tài),PERFORM-3D中使用纖維單元對(duì)塑性區(qū)段進(jìn)行模擬,模型將混凝土梁劃分為10條混凝土纖維和2條鋼筋纖維;將鋼梁劃分為12條纖維,腹板10條纖維,兩翼緣各1條纖維;鋼管混凝土柱劃分為48條纖維,混凝土劃分為36條纖維,鋼管劃分為12條纖維,參考文獻(xiàn)[15]的建議,取截面高度的0.5倍作為塑性區(qū)的長(zhǎng)度。對(duì)于剪力墻,剪力墻的宏觀模型為等代桿系模型和纖維截面模型,對(duì)于高層建筑這類具有大量剪力墻構(gòu)件的結(jié)構(gòu),使用此類宏觀模型可以兼顧計(jì)算精度與效率,剪力墻構(gòu)件采用shear wall宏觀單元模擬,沿截面高度把混凝土和鋼筋各分配8條纖維,沿墻肢高度每層劃分一個(gè)單元。鋼材選取非屈曲型鋼材的材料本構(gòu),骨架曲線擬采用三折線本構(gòu)來(lái)模擬鋼筋與鋼材材料。混凝土可以分為非約束混凝土、箍筋約束混凝土和鋼管混凝土。一般混凝土構(gòu)件中非箍筋加密區(qū)部分的混凝土為非約束混凝土,本構(gòu)關(guān)系按照《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》(GB 50010—2010)附錄C中的混凝土本構(gòu)關(guān)系確定;箍筋約束混凝土采用Kent-Park模型;由于目前沒(méi)有非常權(quán)威的鋼管混凝土本構(gòu)關(guān)系,也采用Kent-Park模型來(lái)定義,相較于箍筋約束混凝土乘以1.1的強(qiáng)度提高系數(shù)。

    為了保證PERFORM-3D軟件建立結(jié)構(gòu)彈塑性分析模型的準(zhǔn)確性,根據(jù)項(xiàng)目的實(shí)際施工圖,先使用YJK軟件進(jìn)行建模并施加相應(yīng)的荷載,通過(guò)計(jì)算得到結(jié)構(gòu)的配筋和模態(tài),并根據(jù)施工圖修改剪力墻等主要構(gòu)件的配筋,為PERFORM-3D模型的建立與校核做好準(zhǔn)備。算例中結(jié)構(gòu)的YJK模型和PERFORM-3D模型如圖3所示。

    PERFORM-3D模型與YJK模型的荷載信息、配筋信息、幾何信息以及質(zhì)量信息都應(yīng)該相接近,PERFORM-3D模型在墻單元中設(shè)置了內(nèi)嵌梁來(lái)模擬連梁和墻的剛接狀態(tài),故需要將兩個(gè)模型的模態(tài)結(jié)果進(jìn)行對(duì)比來(lái)驗(yàn)證PERFORM-3D模型是否與YJK模型存在較大偏差。兩種軟件所建模型的前3階周期對(duì)比結(jié)果見(jiàn)表2,從表2可以看出,兩個(gè)模型的自振周期十分接近,說(shuō)明PERFORM-3D模型相比于YJK模型沒(méi)有存在較大的偏差,此外,還對(duì)兩個(gè)模型的能量誤差以及多遇與罕遇地震作用下基底剪力進(jìn)行校核,校核結(jié)果表明結(jié)構(gòu)整體的各項(xiàng)信息在合理范圍內(nèi),PERFORM-3D模型具有較高的準(zhǔn)確性,可以用于結(jié)構(gòu)的增量動(dòng)力分析。

    3.3 地震波選擇

    通過(guò)分析大量的強(qiáng)震觀測(cè)記錄,《抗規(guī)》規(guī)定設(shè)計(jì)地震動(dòng)用彈性加速度反應(yīng)譜表示,設(shè)計(jì)特征周期定義為“抗震設(shè)計(jì)用的地震影響系數(shù)曲線中,反映地震震級(jí)、震中距和場(chǎng)地類別等因素的下降段起始點(diǎn)對(duì)應(yīng)的周期值”[16]。當(dāng)選擇多條地震波對(duì)結(jié)構(gòu)進(jìn)行抗震分析時(shí),根據(jù)場(chǎng)地特征周期與地震動(dòng)的特征周期接近這一原則進(jìn)行選波,可以從平均意義上得到比較符合工程所在場(chǎng)地特性的地震波。

    綜合上述分析,所選地震波的特征周期應(yīng)與場(chǎng)地特征周期相接近,同時(shí),選擇的地震波還要滿足《高規(guī)》中有效時(shí)長(zhǎng)的要求。為保證結(jié)構(gòu)進(jìn)行增量動(dòng)力分析時(shí)的隨機(jī)性,不選用來(lái)自不同測(cè)站但是同一次地震的記錄。選取12條天然地震波并生成3條人工地震波進(jìn)行結(jié)構(gòu)的增量動(dòng)力分析,所選地震波信息如表3所示,各條地震波與規(guī)范的歸一化加速度反應(yīng)譜,即地震動(dòng)力系數(shù)曲線如圖4所示。

    3.4 增量動(dòng)力分析

    增量動(dòng)力分析(IDA)方法首先應(yīng)確定地震動(dòng)的強(qiáng)度指標(biāo),不同的地震動(dòng)強(qiáng)度指標(biāo)進(jìn)行IDA計(jì)算的結(jié)果具有不同的離散性。選定地震動(dòng)強(qiáng)度指標(biāo)后就要進(jìn)行結(jié)構(gòu)性能參數(shù)的選取,選取的參數(shù)有結(jié)構(gòu)最大基底剪力、結(jié)構(gòu)最大層間位移角。在確定地震動(dòng)的強(qiáng)度指標(biāo)和結(jié)構(gòu)性能參數(shù)后,就應(yīng)確定地震動(dòng)的調(diào)幅方法。

    目前,PGA和Sa(T1,ξ) 是最常用的地震動(dòng)強(qiáng)度指標(biāo),而Sa(T1,ξ)對(duì)于長(zhǎng)周期超高層結(jié)構(gòu)并不適用。由于PGA概念簡(jiǎn)單,獲取地震動(dòng)時(shí)程記錄后就可以直接獲得,中國(guó)規(guī)范目前選定PGA作為地震動(dòng)強(qiáng)度指標(biāo),先采用PGA進(jìn)行調(diào)幅計(jì)算,然后再對(duì)常用的地震動(dòng)強(qiáng)度指標(biāo)和新提出的Sa%進(jìn)行調(diào)幅計(jì)算,Sa%是結(jié)構(gòu)前a%振型質(zhì)量參與系數(shù)的自振周期所對(duì)應(yīng)譜加速度的幾何平均值,其首先要計(jì)算各階振型所對(duì)應(yīng)的質(zhì)量參與系數(shù),從而得到Sa%所對(duì)應(yīng)的振型,Sa%最終的計(jì)算結(jié)果如表4所示。對(duì)于PGA的調(diào)幅,由于中國(guó)《抗規(guī)》中9度烈度區(qū)罕遇地震加速度時(shí)程的最大值為620 cm/s2,且有發(fā)生“特大地震”的可能,所以,IDA分析中PGA上限值取1 000 cm/s2,為了保證結(jié)構(gòu)獲得較精確的計(jì)算結(jié)果并具有高效的計(jì)算效率,PGA調(diào)幅按照50、100、200、300、400、600、800、1 000 cm/s2進(jìn)行。PGA在400 cm/s2之前,考慮到結(jié)構(gòu)由彈性逐漸進(jìn)入彈塑性階段,構(gòu)件的剛度和強(qiáng)度均有不同程度的退化,如果步長(zhǎng)選取過(guò)長(zhǎng),會(huì)導(dǎo)致這一階段計(jì)算結(jié)果的準(zhǔn)確度不高,所以,這一階段調(diào)幅步長(zhǎng)選取相對(duì)較小。PGA在400 cm/s2之后已屬于罕遇地震,在結(jié)構(gòu)設(shè)計(jì)基準(zhǔn)期內(nèi)發(fā)生的概率非常小,并且當(dāng)發(fā)生這類大地震時(shí),結(jié)構(gòu)進(jìn)入塑性的程度較高,每次彈塑性時(shí)程分析都需要大量的時(shí)間,綜合考慮這兩點(diǎn)因素,故適當(dāng)增加了步長(zhǎng)。對(duì)于其他地震動(dòng)指標(biāo)相對(duì)于PGA的調(diào)幅值,可按前述方法進(jìn)行計(jì)算,在地震動(dòng)Y1作用下,PGA與各個(gè)地震動(dòng)強(qiáng)度指標(biāo)相對(duì)應(yīng)的調(diào)幅值如表5所示。為使用分位數(shù)曲線法,保證每條地震動(dòng)最終都能將θmax計(jì)算至0.01以上,Vmax都能計(jì)算至100 000 kN以上。

    在PERFORM-3D模型中,瑞利阻尼與模態(tài)阻尼是疊加的,結(jié)構(gòu)的阻尼采用統(tǒng)一的模態(tài)阻尼加較小的瑞利阻尼來(lái)模擬,模態(tài)阻尼比統(tǒng)一取0.04。每條地震動(dòng)按規(guī)范建議采取雙向輸入,X向和Y向的地震波峰值取為1∶0.85。

    4 地震動(dòng)強(qiáng)度指標(biāo)離散性分析

    主要研究結(jié)構(gòu)X向的地震動(dòng)響應(yīng),以下結(jié)構(gòu)的響應(yīng)(基底剪力、層間位移角)都是X方向。依據(jù)IDA計(jì)算結(jié)果,可以得到地震波在各地震動(dòng)強(qiáng)度指標(biāo)下的基底剪力和最大層間位移角的IDA曲線,限于篇幅,僅給出地震波在PGA下的基底剪力和最大層間位移角的IDA曲線,如圖5和圖6所示。根據(jù)IDA計(jì)算的結(jié)果和前述方法,可以得到各地震動(dòng)強(qiáng)度指標(biāo)對(duì)于θmax和Vmax的分位數(shù)曲線,限于篇幅,只給出PGA和PGV對(duì)于θmax的分位數(shù)曲線,如圖7和圖8所示。基于θmax計(jì)算得到的分位數(shù)曲線離散程度平均誤差如表6所示,基于Vmax計(jì)算得到的分位數(shù)曲線離散程度平均誤差如表7所示。計(jì)算各地震動(dòng)強(qiáng)度指標(biāo)相對(duì)于θmax和Vmax的條件對(duì)數(shù)標(biāo)準(zhǔn)差平均值如表8所示。

    從表6、表7和表8中可以得到以下結(jié)論:

    1)在最常用的3個(gè)地震動(dòng)強(qiáng)度指標(biāo)PGA、PGV、Sa(T1,4%)中,PGV得到的分位數(shù)曲線離散程度最小,PGA次之,Sa(T1,4%)最大,這也與眾多學(xué)者認(rèn)為Sa(T1,ξ)只適用于第一振型反應(yīng)為主導(dǎo)的低層、規(guī)則結(jié)構(gòu)的結(jié)論保持一致。對(duì)于Vmax這一結(jié)構(gòu)性能指標(biāo),PGV基于Vmax得到的分位數(shù)曲線離散程度相較于θmax有所增大,與PGA接近。PSA和PSV的概念雖然簡(jiǎn)單,但離散程度大于PGA和PGV,實(shí)用價(jià)值不高。

    2)S12、S123、S12345得到的分位數(shù)曲線離散程度也較大,相較于Sa(T1,4%)未能取得良好的效果,這是因?yàn)樗憷Y(jié)構(gòu)較復(fù)雜,高階振型的質(zhì)量參與系數(shù)很小,所以,從一定程度上弱化了高階振型對(duì)結(jié)構(gòu)的影響。S3(S70%)得到的分位數(shù)曲線離散程度小于S123(S70%)和S12345(S80%)也驗(yàn)證了這一點(diǎn)。

    3)隨著所考慮質(zhì)量參與系數(shù)的百分比不斷增加,Sa%離散程度也不斷減小。S90%得到的分位數(shù)曲線離散程度相比于S80%要小,這也說(shuō)明,若以S90%作為地震動(dòng)強(qiáng)度指標(biāo)調(diào)幅進(jìn)行IDA計(jì)算,得到的θmax和Vmax結(jié)果相對(duì)文中的其他地震動(dòng)指標(biāo)的有效性較好。

    通過(guò)以上研究可以看出,提出的Sa%在考慮結(jié)構(gòu)X方向平動(dòng)前90%振型質(zhì)量參與系數(shù)的自振周期對(duì)應(yīng)的譜加速度的幾何平均值S90%作為地震動(dòng)強(qiáng)度指標(biāo)調(diào)幅進(jìn)行IDA計(jì)算時(shí),可以有效減少θmax和Vmax數(shù)據(jù)的離散性,提高計(jì)算數(shù)據(jù)的有效性,從而適當(dāng)減少時(shí)程分析的次數(shù)。低層和多層等較規(guī)則結(jié)構(gòu)由于第一平動(dòng)振型所占的質(zhì)量參與系數(shù)較大,往往達(dá)到70%~80%甚至更多,所以,Sa(T1,ξ)等效于S70%或S80%,從而具有優(yōu)于PGA的有效性,但對(duì)于超高層結(jié)構(gòu)則不適用??紤]高階振型參與的冪函數(shù)乘積形式的地震動(dòng)強(qiáng)度指標(biāo)S12、S123和S12345的效果也不理想。PGV提高數(shù)據(jù)有效性的效果僅次于S90%,其相比于PGA具有一定優(yōu)勢(shì),也值得考慮使用。

    5 結(jié)論

    在適用于超高層建筑結(jié)構(gòu)的地震動(dòng)強(qiáng)度指標(biāo)S123和Sa的基礎(chǔ)上,提出考慮高階振型影響并適用于超高層結(jié)構(gòu)的地震動(dòng)強(qiáng)度指標(biāo)Sa%;根據(jù)某實(shí)際超高層建筑結(jié)構(gòu)的IDA計(jì)算結(jié)果,使用分位數(shù)曲線法,并以分位數(shù)離散程度平均誤差和條件對(duì)數(shù)標(biāo)準(zhǔn)差均值最小為依據(jù),驗(yàn)證提出的地震動(dòng)強(qiáng)度指標(biāo)S90%相對(duì)于其他地震動(dòng)指標(biāo)的有效性最好。得到的主要結(jié)論如下:

    1)對(duì)于最常用的幾個(gè)地震動(dòng)強(qiáng)度指標(biāo),PGV得到的分位數(shù)曲線離散程度最小,PGA次之,Sa(T1,4%)、PSA和PSV的離散程度較大,不適用于高階振型影響較大的復(fù)雜超高層結(jié)構(gòu)。

    2)相比于現(xiàn)有的若干個(gè)常用的和適用于超高層結(jié)構(gòu)的地震動(dòng)強(qiáng)度指標(biāo),在已有研究的基礎(chǔ)上提出的S90%在用于超高層混合結(jié)構(gòu)的增量動(dòng)力分析時(shí),能更有效減少θmax和Vmax數(shù)據(jù)的離散性,從而提高計(jì)算數(shù)據(jù)的有效性。

    參考文獻(xiàn):

    [1] 劉恢先.論地震力[J].土木工程學(xué)報(bào), 1958, 5(2): 86-106.

    LIU H X. On seismic force [J]. Journal of Civil Engineering, 1958, 5(2): 86-106. (in Chinese)

    [2] 楊迪雄, 龍海豐. 近斷層地震動(dòng)強(qiáng)度參數(shù)與雙線型SDOF體系地震反應(yīng)的回歸分析[J]. 地震工程與工程振動(dòng), 2014, 34(3): 9-17.

    YANG D X, LONG H F. Regressive analysis of intensity parameters of near-fault ground motions and seismic responses of bilinear SDOF systems [J]. Earthquake Engineering and Engineering Dynamics, 2014, 34(3): 9-17. (in Chinese)

    [3] BERTERO V V. Strength and deformation capacities of buildings under extreme environments [J]. Structural Engineering and Structural Mechanics, 1977, 53(1): 29-79.

    [4] ASGARIAN B, JALAEEFAR A. Incremental dynamic analysis of steel braced frames designed based on the first, second and third editions of the Iranian seismic code (standard no. 2800) [J]. The Structural Design of Tall and Special Buildings, 2011, 20(2): 190-207.

    [5] 周穎, 呂西林, 卜一. 增量動(dòng)力分析法在高層混合結(jié)構(gòu)性能評(píng)估中的應(yīng)用[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 38(2): 183-187, 193.

    ZHOU Y, LU X L, BU Y. Application of incremental dynamic analysis to seismic evaluation of hybrid structure [J]. Journal of Tongji University (Natural Science), 2010, 38(2): 183-187, 193. (in Chinese)

    [6] 馬千里, 葉列平, 陸新征, 等. 采用逐步增量彈塑性時(shí)程方法對(duì)RC框架結(jié)構(gòu)推覆分析側(cè)力模式的研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2008, 29(2): 132-140.

    MA Q L, YE L P, LU X Z, et al. Study on lateral load patterns of pushover analysis using incremental dynamical analysis for RC frame structures [J]. Journal of Building Structures, 2008, 29(2): 132-140. (in Chinese)

    [7] VAMVATSIKOS D. Seismic performance, capacity and reliability of structures as seen through incremental dynamic analysis [D]. Stanford: Stanford University, 2002.

    [8] GUAN M S, DU H B, CUI J, et al. Optimal ground motion intensity measure for long-period structures [J]. Measurement Science and Technology, 2015, 26(10): 105001.

    [9] 盧嘯, 陸新征, 葉列平, 等. 適用于超高層建筑的改進(jìn)地震動(dòng)強(qiáng)度指標(biāo)[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2014, 35(2): 15-21.

    LU X, LU X Z, YE L P, et al. Development of an improved ground motion intensity measure for super high-rise buildings [J]. Journal of Building Structures, 2014, 35(2): 15-21. (in Chinese)

    [10] 周穎, 蘇寧粉, 呂西林. 高層建筑結(jié)構(gòu)增量動(dòng)力分析的地震動(dòng)強(qiáng)度參數(shù)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2013, 34(2): 53-60.

    ZHOU Y, SU N F, LU X L. Study on intensity measure of incremental dynamic analysis for high-rise structures [J]. Journal of Building Structures, 2013, 34(2): 53-60.(in Chinese)

    [11] 侯喆辰. 高層住宅剪力墻結(jié)構(gòu)的性能設(shè)計(jì)方法和地震易損性分析研究[D]. 廣州: 華南理工大學(xué), 2017.

    HOU Z C. Study on performance-based design and seismic fragility analysis of high-rise residential shear wall structure [D]. Guangzhou: South China University of Technology, 2017. (in Chinese)

    [12] 郝永蕾. 基于IDA的高層鋼混凝土混合結(jié)構(gòu)地震易損性分析[D]. 西安: 西安科技大學(xué), 2018.

    HAO Y L. Seismic vulnerability analysis of SRC frame-RC core wall hybrid structure in high-rise buildings based on the IDA method [D]. Xian: Xian University of Science and Technology, 2018. (in Chinese)

    [13] 蔣歡軍, 和留生, 呂西林, 等. 上海中心大廈抗震性能分析和振動(dòng)臺(tái)試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2011, 32(11): 55-63.

    JIANG H J, HE L S, LU X L, et al. Analysis of seismic performance and shaking table tests of the Shanghai Tower [J]. Journal of Building Structures, 2011, 32(11): 55-63.(in Chinese)

    [14] FAN H, LI Q S, TUAN A Y, et al. Seismic analysis of the world′s tallest building [J]. Journal of Constructional Steel Research, 2009, 65(5): 1206-1215.

    [15] CALVI G M, PRIESTLEY M J N, KOWALSKY M J. Displacement-based seismic design of structures[J]. Earthquake Spectra, 2007, 24(2): 1-24.

    [16] 陳鵬, 劉文鋒, 付興潘. 關(guān)于場(chǎng)地卓越周期和特征周期的若干討論[J]. 青島理工大學(xué)學(xué)報(bào), 2009, 30(6): 30-35.

    CHEN P, LIU W F, FU X P. Discussions on site predominant period and characteristic period [J]. Journal of Qingdao Technological University, 2009, 30(6): 30-35. (in Chinese)

    (編輯 王秀玲)

    猜你喜歡
    調(diào)幅位數(shù)震動(dòng)
    五次完全冪的少位數(shù)三進(jìn)制展開(kāi)
    震動(dòng)減脂儀可以減肥?
    基于MATLAB調(diào)幅包絡(luò)檢波和相干解調(diào)性能設(shè)計(jì)與比較
    電子制作(2019年9期)2019-05-30 09:42:12
    一種用于調(diào)幅接收機(jī)AGC的設(shè)計(jì)與實(shí)現(xiàn)
    水電工程場(chǎng)地地震動(dòng)確定方法
    振動(dòng)攪拌 震動(dòng)創(chuàng)新
    調(diào)幅翻轉(zhuǎn)式超深翻犁的研制及推廣應(yīng)用
    調(diào)頻引信中噪聲調(diào)幅干擾的自適應(yīng)抑制
    人工合成最不利地震動(dòng)
    遙感衛(wèi)星CCD相機(jī)量化位數(shù)的選擇
    欧美最黄视频在线播放免费| 国产av麻豆久久久久久久| 久久亚洲真实| 制服丝袜大香蕉在线| 综合色av麻豆| 国产精品女同一区二区软件 | 我的女老师完整版在线观看| avwww免费| 97超视频在线观看视频| 99久久无色码亚洲精品果冻| 欧洲精品卡2卡3卡4卡5卡区| av天堂中文字幕网| 美女免费视频网站| 国产高清有码在线观看视频| 欧美3d第一页| 国产精品亚洲美女久久久| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 欧美乱妇无乱码| 丰满人妻熟妇乱又伦精品不卡| 好看av亚洲va欧美ⅴa在| .国产精品久久| 日日夜夜操网爽| 国产午夜福利久久久久久| 国产精品99久久久久久久久| 中文字幕高清在线视频| 可以在线观看毛片的网站| 在线观看舔阴道视频| 国语自产精品视频在线第100页| 亚洲一区二区三区色噜噜| 亚洲欧美日韩无卡精品| 欧美一区二区亚洲| 性色avwww在线观看| 可以在线观看的亚洲视频| 国产亚洲精品久久久久久毛片| 如何舔出高潮| 国模一区二区三区四区视频| 757午夜福利合集在线观看| 搞女人的毛片| 欧美又色又爽又黄视频| 亚洲av电影在线进入| 国产真实乱freesex| 久久国产精品人妻蜜桃| av在线蜜桃| 宅男免费午夜| 1024手机看黄色片| 欧美潮喷喷水| 国产精品综合久久久久久久免费| 国产精品日韩av在线免费观看| 久9热在线精品视频| 午夜精品一区二区三区免费看| 欧美区成人在线视频| 搡女人真爽免费视频火全软件 | 三级国产精品欧美在线观看| 久久久久久大精品| 久久精品国产99精品国产亚洲性色| 国产熟女xx| 自拍偷自拍亚洲精品老妇| 国产精品自产拍在线观看55亚洲| 俺也久久电影网| 乱码一卡2卡4卡精品| 在线a可以看的网站| 久久99热这里只有精品18| 丝袜美腿在线中文| 老女人水多毛片| 欧美高清成人免费视频www| 亚洲天堂国产精品一区在线| 十八禁网站免费在线| 午夜福利免费观看在线| 国产精品自产拍在线观看55亚洲| 久久6这里有精品| 嫁个100分男人电影在线观看| 国产淫片久久久久久久久 | 国产av麻豆久久久久久久| av天堂中文字幕网| 国产av在哪里看| 亚洲美女视频黄频| 久久精品国产亚洲av香蕉五月| 天堂动漫精品| 在线看三级毛片| 在线播放无遮挡| 午夜免费男女啪啪视频观看 | 少妇熟女aⅴ在线视频| 丁香欧美五月| 久久久久九九精品影院| 免费一级毛片在线播放高清视频| 不卡一级毛片| 国产老妇女一区| 国产69精品久久久久777片| 无人区码免费观看不卡| 精品久久久久久久久久久久久| 免费在线观看日本一区| 国产伦在线观看视频一区| 久久国产精品影院| 熟女电影av网| 此物有八面人人有两片| av欧美777| 哪里可以看免费的av片| 69人妻影院| 国产美女午夜福利| 人妻久久中文字幕网| 久久久久精品国产欧美久久久| 国产综合懂色| 麻豆一二三区av精品| 欧美+日韩+精品| 在线a可以看的网站| 中文字幕av成人在线电影| av欧美777| 狠狠狠狠99中文字幕| 成人性生交大片免费视频hd| 国内久久婷婷六月综合欲色啪| 国产av麻豆久久久久久久| 少妇高潮的动态图| 老熟妇仑乱视频hdxx| 国产私拍福利视频在线观看| or卡值多少钱| 国产精品嫩草影院av在线观看 | 男女做爰动态图高潮gif福利片| 免费av毛片视频| 深夜精品福利| 观看美女的网站| 最后的刺客免费高清国语| 首页视频小说图片口味搜索| 亚洲欧美日韩高清在线视频| 国产成+人综合+亚洲专区| 99久久久亚洲精品蜜臀av| 国产一区二区在线观看日韩| 99国产综合亚洲精品| 俺也久久电影网| 亚洲美女搞黄在线观看 | 狠狠狠狠99中文字幕| 老司机午夜十八禁免费视频| 亚洲黑人精品在线| 精品久久久久久久久av| 舔av片在线| 99在线视频只有这里精品首页| 男人狂女人下面高潮的视频| 男女那种视频在线观看| 熟妇人妻久久中文字幕3abv| 国产极品精品免费视频能看的| 最好的美女福利视频网| 有码 亚洲区| av视频在线观看入口| av女优亚洲男人天堂| 99久国产av精品| 人人妻人人看人人澡| 精品福利观看| 神马国产精品三级电影在线观看| 我的女老师完整版在线观看| 亚洲国产欧美人成| 变态另类丝袜制服| 长腿黑丝高跟| 琪琪午夜伦伦电影理论片6080| a级毛片a级免费在线| 亚洲成av人片免费观看| 欧美国产日韩亚洲一区| 波多野结衣巨乳人妻| 我的老师免费观看完整版| 一a级毛片在线观看| 床上黄色一级片| 18禁在线播放成人免费| 国产精品久久久久久人妻精品电影| 成年人黄色毛片网站| 很黄的视频免费| 人妻丰满熟妇av一区二区三区| 国产伦一二天堂av在线观看| 久久欧美精品欧美久久欧美| 色综合站精品国产| 能在线免费观看的黄片| 国内精品久久久久精免费| 日本黄色视频三级网站网址| 亚洲,欧美精品.| 最后的刺客免费高清国语| 黄色丝袜av网址大全| 婷婷亚洲欧美| 午夜福利欧美成人| 69人妻影院| 99久国产av精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美日韩无卡精品| or卡值多少钱| 亚洲人成电影免费在线| 十八禁国产超污无遮挡网站| 国产乱人伦免费视频| 亚洲成人精品中文字幕电影| 亚洲熟妇中文字幕五十中出| 69av精品久久久久久| 久久久国产成人精品二区| 少妇的逼水好多| av天堂在线播放| 久久精品国产99精品国产亚洲性色| 国产伦精品一区二区三区四那| av在线蜜桃| 男女床上黄色一级片免费看| 久久亚洲真实| 亚洲第一电影网av| 精品福利观看| 亚洲精品成人久久久久久| 亚洲精品粉嫩美女一区| 欧美乱妇无乱码| 色尼玛亚洲综合影院| 白带黄色成豆腐渣| 精华霜和精华液先用哪个| 亚洲av一区综合| 日本撒尿小便嘘嘘汇集6| 日日摸夜夜添夜夜添小说| 国产免费一级a男人的天堂| 深夜精品福利| 在线免费观看不下载黄p国产 | av天堂在线播放| 18禁裸乳无遮挡免费网站照片| 色综合婷婷激情| 久久婷婷人人爽人人干人人爱| 俺也久久电影网| 国产伦在线观看视频一区| 亚洲五月天丁香| 观看免费一级毛片| av欧美777| 51国产日韩欧美| 国产黄片美女视频| 一级毛片久久久久久久久女| 国产高清三级在线| 亚洲一区二区三区色噜噜| 亚洲精品色激情综合| 午夜福利视频1000在线观看| 直男gayav资源| 亚洲国产精品合色在线| 日韩成人在线观看一区二区三区| 18禁在线播放成人免费| 欧美日韩福利视频一区二区| 免费在线观看日本一区| 成年女人看的毛片在线观看| 欧美精品国产亚洲| 麻豆国产av国片精品| 怎么达到女性高潮| 身体一侧抽搐| 国产精品日韩av在线免费观看| 桃红色精品国产亚洲av| 国产大屁股一区二区在线视频| 国产视频一区二区在线看| 国产免费av片在线观看野外av| 国模一区二区三区四区视频| 日韩欧美在线乱码| 国产91精品成人一区二区三区| 中文字幕精品亚洲无线码一区| 三级男女做爰猛烈吃奶摸视频| 亚洲狠狠婷婷综合久久图片| 国产精品一区二区免费欧美| 搡女人真爽免费视频火全软件 | 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 欧美bdsm另类| 首页视频小说图片口味搜索| 中文字幕av成人在线电影| 国模一区二区三区四区视频| 亚洲av五月六月丁香网| 我要搜黄色片| 国产精品人妻久久久久久| 搡女人真爽免费视频火全软件 | 国产在线男女| 日韩免费av在线播放| 免费高清视频大片| 给我免费播放毛片高清在线观看| 乱码一卡2卡4卡精品| 一个人看的www免费观看视频| 国产蜜桃级精品一区二区三区| 天天躁日日操中文字幕| 中文在线观看免费www的网站| 亚洲第一区二区三区不卡| 我要看日韩黄色一级片| 真人做人爱边吃奶动态| 少妇的逼好多水| 久久久久精品国产欧美久久久| 色在线成人网| 久久久久性生活片| 国产一区二区三区视频了| 丝袜美腿在线中文| 亚洲一区二区三区不卡视频| a级毛片a级免费在线| 亚洲在线观看片| 国产主播在线观看一区二区| 精品日产1卡2卡| 此物有八面人人有两片| 一进一出抽搐动态| 国产爱豆传媒在线观看| 99riav亚洲国产免费| 亚洲欧美日韩东京热| 一个人免费在线观看的高清视频| 成人精品一区二区免费| 成年女人永久免费观看视频| 日韩欧美在线二视频| 最好的美女福利视频网| 精品熟女少妇八av免费久了| 成人av一区二区三区在线看| 给我免费播放毛片高清在线观看| 亚洲中文字幕日韩| 97热精品久久久久久| 久久精品国产亚洲av涩爱 | 国产精品一区二区三区四区免费观看 | 国产成+人综合+亚洲专区| 亚洲精品粉嫩美女一区| 可以在线观看毛片的网站| 精品乱码久久久久久99久播| 精品午夜福利视频在线观看一区| 国产极品精品免费视频能看的| 国产精品三级大全| 99热6这里只有精品| 特级一级黄色大片| 中文字幕人成人乱码亚洲影| 天天躁日日操中文字幕| 久久久精品大字幕| 男女下面进入的视频免费午夜| 两人在一起打扑克的视频| 动漫黄色视频在线观看| 欧美成人a在线观看| 欧美+日韩+精品| or卡值多少钱| 亚洲午夜理论影院| 国产精品爽爽va在线观看网站| 欧美高清性xxxxhd video| 精品国产亚洲在线| 综合色av麻豆| 亚洲男人的天堂狠狠| 亚洲电影在线观看av| 国产一区二区激情短视频| 欧美黄色淫秽网站| 国产免费男女视频| 草草在线视频免费看| 中文字幕精品亚洲无线码一区| 欧美激情久久久久久爽电影| 亚洲无线在线观看| 看黄色毛片网站| 久久精品影院6| 久久久久久久亚洲中文字幕 | 久久人妻av系列| 他把我摸到了高潮在线观看| 国产在视频线在精品| 男女做爰动态图高潮gif福利片| 人妻久久中文字幕网| 亚洲精品日韩av片在线观看| 国产成人欧美在线观看| 村上凉子中文字幕在线| 又爽又黄无遮挡网站| 免费看a级黄色片| 国内精品久久久久久久电影| 亚洲专区中文字幕在线| 一进一出抽搐动态| 人人妻,人人澡人人爽秒播| 嫩草影院入口| 国产探花在线观看一区二区| 美女大奶头视频| 亚洲五月婷婷丁香| 日韩欧美在线乱码| 少妇丰满av| 毛片女人毛片| 亚洲中文字幕一区二区三区有码在线看| 最后的刺客免费高清国语| 黄色日韩在线| 日韩高清综合在线| 999久久久精品免费观看国产| 一本久久中文字幕| 色哟哟哟哟哟哟| 久久九九热精品免费| 亚洲电影在线观看av| 午夜福利视频1000在线观看| 悠悠久久av| 久9热在线精品视频| 12—13女人毛片做爰片一| 精品一区二区三区视频在线| 国产爱豆传媒在线观看| 成人性生交大片免费视频hd| 波多野结衣高清无吗| 51国产日韩欧美| 嫩草影视91久久| 99热6这里只有精品| 人人妻,人人澡人人爽秒播| 免费在线观看成人毛片| 精品人妻视频免费看| 男人狂女人下面高潮的视频| avwww免费| 日韩欧美在线二视频| 国产伦精品一区二区三区视频9| 亚洲男人的天堂狠狠| 热99在线观看视频| 老女人水多毛片| 日本精品一区二区三区蜜桃| 欧美最黄视频在线播放免费| 一个人免费在线观看的高清视频| 特级一级黄色大片| 99久久无色码亚洲精品果冻| 亚洲人成网站在线播放欧美日韩| 久久久久九九精品影院| 午夜福利在线观看免费完整高清在 | 亚洲无线观看免费| 亚洲七黄色美女视频| 国产精品一区二区三区四区免费观看 | av国产免费在线观看| 特大巨黑吊av在线直播| 久久精品国产亚洲av香蕉五月| 精品久久久久久久久久久久久| 日本与韩国留学比较| 免费搜索国产男女视频| 国产精品三级大全| 亚洲精品456在线播放app | 十八禁国产超污无遮挡网站| 久久国产乱子伦精品免费另类| 给我免费播放毛片高清在线观看| 韩国av一区二区三区四区| 97超视频在线观看视频| 老司机午夜福利在线观看视频| 不卡一级毛片| 国产精品久久久久久人妻精品电影| 亚洲午夜理论影院| 波多野结衣高清作品| 久久人妻av系列| 亚洲av成人av| 波多野结衣巨乳人妻| 俄罗斯特黄特色一大片| 白带黄色成豆腐渣| 国产亚洲精品久久久com| 99久久久亚洲精品蜜臀av| 51国产日韩欧美| 最好的美女福利视频网| 99热这里只有精品一区| 亚洲av美国av| 中文字幕高清在线视频| 亚洲七黄色美女视频| 国产精品98久久久久久宅男小说| 毛片女人毛片| 国产精品久久久久久亚洲av鲁大| 欧美成人免费av一区二区三区| 久久精品国产亚洲av香蕉五月| 久久国产精品人妻蜜桃| 久久久久久久午夜电影| 伊人久久精品亚洲午夜| www.www免费av| av在线蜜桃| 精品无人区乱码1区二区| 91av网一区二区| 亚洲精品乱码久久久v下载方式| 久久中文看片网| 成人欧美大片| 男女之事视频高清在线观看| 国产精品亚洲一级av第二区| 精品人妻熟女av久视频| 精品人妻视频免费看| 亚洲不卡免费看| 宅男免费午夜| 看黄色毛片网站| 变态另类丝袜制服| 99久久精品国产亚洲精品| 亚洲国产精品久久男人天堂| 好男人电影高清在线观看| 性色avwww在线观看| 伊人久久精品亚洲午夜| 亚洲av美国av| 亚洲人成网站在线播| 精品人妻1区二区| 成人性生交大片免费视频hd| 成人一区二区视频在线观看| 国内少妇人妻偷人精品xxx网站| 俺也久久电影网| 国产毛片a区久久久久| 极品教师在线免费播放| 国产伦精品一区二区三区四那| 亚州av有码| 久久人妻av系列| 欧美午夜高清在线| 国产成人欧美在线观看| 欧美一区二区精品小视频在线| 18禁裸乳无遮挡免费网站照片| 白带黄色成豆腐渣| 亚洲av一区综合| 窝窝影院91人妻| 国产免费av片在线观看野外av| 中出人妻视频一区二区| 欧美午夜高清在线| 色综合婷婷激情| 嫩草影院新地址| 亚洲精品一卡2卡三卡4卡5卡| 欧美乱妇无乱码| 婷婷精品国产亚洲av在线| 中文字幕免费在线视频6| 美女高潮喷水抽搐中文字幕| av专区在线播放| 国产精品一区二区三区四区久久| 精品99又大又爽又粗少妇毛片 | 国产一区二区三区视频了| 久久午夜福利片| 悠悠久久av| 性欧美人与动物交配| 欧美午夜高清在线| 在线国产一区二区在线| 欧美在线一区亚洲| 身体一侧抽搐| 老司机午夜福利在线观看视频| 国产精品,欧美在线| 看片在线看免费视频| 亚洲七黄色美女视频| 久久精品91蜜桃| 三级毛片av免费| 午夜a级毛片| 啦啦啦韩国在线观看视频| 91狼人影院| 日韩欧美精品免费久久 | 91麻豆av在线| www.999成人在线观看| 亚洲无线观看免费| 国产成人影院久久av| 99国产极品粉嫩在线观看| 欧美国产日韩亚洲一区| 色综合站精品国产| 亚洲av二区三区四区| 国产真实乱freesex| 美女免费视频网站| 久久久精品大字幕| 日本a在线网址| 动漫黄色视频在线观看| 男人狂女人下面高潮的视频| 男女床上黄色一级片免费看| 国产探花极品一区二区| 男人和女人高潮做爰伦理| 亚洲真实伦在线观看| 中文字幕免费在线视频6| av欧美777| 一本精品99久久精品77| 九九在线视频观看精品| 久久中文看片网| 精品国内亚洲2022精品成人| 国产一区二区三区在线臀色熟女| 熟女人妻精品中文字幕| 日韩欧美精品v在线| 久久草成人影院| 少妇丰满av| 午夜久久久久精精品| 精品一区二区免费观看| 老女人水多毛片| 哪里可以看免费的av片| 亚洲av中文字字幕乱码综合| 国产亚洲欧美在线一区二区| 国产在线男女| 18禁在线播放成人免费| 欧美日韩国产亚洲二区| 少妇裸体淫交视频免费看高清| 国产免费男女视频| 精品免费久久久久久久清纯| 欧美三级亚洲精品| 亚洲第一欧美日韩一区二区三区| 国产av不卡久久| 人妻制服诱惑在线中文字幕| 亚洲av免费高清在线观看| 亚洲美女黄片视频| 国内精品久久久久精免费| av视频在线观看入口| 日韩免费av在线播放| 国产精品永久免费网站| 色综合站精品国产| 欧美区成人在线视频| 欧美3d第一页| 天堂√8在线中文| 少妇熟女aⅴ在线视频| 成人av在线播放网站| 国产高清激情床上av| 麻豆一二三区av精品| 中国美女看黄片| 午夜影院日韩av| av在线老鸭窝| 国产高清有码在线观看视频| 老熟妇仑乱视频hdxx| 色尼玛亚洲综合影院| 99国产综合亚洲精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲第一欧美日韩一区二区三区| 精品久久久久久久久av| 亚洲专区国产一区二区| 欧美丝袜亚洲另类 | 欧美性猛交黑人性爽| 欧美区成人在线视频| 又爽又黄无遮挡网站| 九色成人免费人妻av| 波野结衣二区三区在线| 色视频www国产| 国产成年人精品一区二区| 国产伦一二天堂av在线观看| 国产欧美日韩一区二区精品| 国产主播在线观看一区二区| 亚洲国产日韩欧美精品在线观看| 成年人黄色毛片网站| 最新中文字幕久久久久| 国产亚洲精品综合一区在线观看| 亚洲欧美清纯卡通| eeuss影院久久| 18+在线观看网站| 国产亚洲精品久久久com| 国产不卡一卡二| 免费av观看视频| 在现免费观看毛片| 一a级毛片在线观看| 免费av观看视频| 久久午夜亚洲精品久久| 亚洲美女视频黄频| 国产成+人综合+亚洲专区| 99精品久久久久人妻精品| 精品一区二区免费观看| 老司机深夜福利视频在线观看| 精品人妻熟女av久视频| 国产精品免费一区二区三区在线| 欧美日韩瑟瑟在线播放| bbb黄色大片| 国产91精品成人一区二区三区| 在线十欧美十亚洲十日本专区| 精品熟女少妇八av免费久了| 久久精品国产清高在天天线| 日本熟妇午夜| 亚洲人成伊人成综合网2020|