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

    增量動力分析地震動強度參數(shù)研究綜述

    2014-03-21 09:07:00勵勐劼
    結(jié)構(gòu)工程師 2014年6期
    關(guān)鍵詞:標量震動振型

    周 穎 勵勐劼

    (同濟大學土木工程防災(zāi)國家重點實驗室,上海200092)

    1 引言

    增量動力分析(Incremental Dynamic Analysis,IDA)被稱作動力推覆分析,是目前基于性能的地震工程中最具有發(fā)展前景的方法之一,已被美國聯(lián)邦緊急救援署(Federal Emergency Management Agency,F(xiàn)EMA)所采用,作為分析鋼框架結(jié)構(gòu)整體抗倒塌能力的方法[1]。

    增量動力分析結(jié)果是通過地震動強度參數(shù)(Intensity Measure,IM)與工程需求參數(shù)(Engineering demand parameter,EDP)關(guān)系來表示,可以反映結(jié)構(gòu)歷經(jīng)初始彈性、逐步退化直至整體動力失穩(wěn)的全過程響應(yīng)。目前常用的工程需求參數(shù),例如最大層間位移角θmax能較好反映結(jié)構(gòu)破壞程度[2],具有一定的合理性[3]。而選用不同地震動強度參數(shù)將獲得不同的分析結(jié)果,是影響結(jié)構(gòu)抗震性能評估結(jié)果的關(guān)鍵因素之一[4]。近年來各國學者從不同角度對地震動強度參數(shù)進行了研究,本文從標量型地震動強度參數(shù)和向量型地震動強度參數(shù)展開分析。

    2 地震動強度參數(shù)基本性質(zhì)

    地震動強度參數(shù)是表征地震動強度的指標,其應(yīng)具有以下性質(zhì),并可作為判斷地震動強度參數(shù)優(yōu)劣的標準。

    (1)有效性(Efficiency):有效性可降低需求中位值估計誤差,減少分析所需地震動數(shù)量,可通過條件對數(shù)標準差來判別。

    (2)充分性(Sufficiency):指地震動強度參數(shù)在統(tǒng)計意義上對地震特征參數(shù)(震級M,震中距R等)的條件獨立性,通常可用P值判斷。

    (3)相關(guān)性(Correlation):反映地震動強度參數(shù)與工程需求參數(shù)的關(guān)聯(lián)程度,可用相關(guān)系數(shù)的平方來判斷。

    (4)實用性(Practicality):指地震動強度參數(shù)與工程需求參數(shù)的關(guān)聯(lián)敏感性。地震動強度參數(shù)實用性不強,則工程需求參數(shù)很少或者不依賴于地震動強度參數(shù)的變化。實用性通??捎镁€性回歸系數(shù)b來判斷。

    (5)適用性(Proficiency):綜合考慮有效性和實用性的性質(zhì)。通常用條件對數(shù)標準差與線性回歸系數(shù)b的比值來判斷。

    (6)比例魯棒性(Scaling Robustness):指地震動強度參數(shù)按比例調(diào)整后,結(jié)構(gòu)的響應(yīng)結(jié)果具有的無偏性質(zhì)。

    (7)危險性可計算性(Hazard Computability):地震動強度參數(shù)在計算地震危險性時所應(yīng)具有的簡易特性。

    3 標量型地震動強度參數(shù)

    峰值加速度(Peak Ground Acceleration,PGA)概念直觀簡單,是最早使用也是目前使用最多的地震動強度參數(shù),大多國家都采用這一指標。

    在1998年,Bazzurro等[5]提出了用結(jié)構(gòu)彈性基本周期對應(yīng)的加速度反應(yīng)譜譜值Sa(T1)作為地震動強度參數(shù)。隨后Vamvatasikos等[6]在2002年對這兩個參數(shù)進行了比較,認為Sa(T1)可以明顯降低增量動力分析結(jié)果的離散程度。Sa(T1)也逐漸發(fā)展為目前最為常用的地震動強度參數(shù)之一。但隨著研究的不斷深入,發(fā)現(xiàn)Sa(T1)作為地震動強度參數(shù)并不完善,例如不適用于長周期結(jié)構(gòu)等,從而衍生了許多改進地震動強度參數(shù)或者通過其他角度切入的地震動強度參數(shù)。

    3.1 基于加速度反應(yīng)譜譜形的地震動強度參數(shù)

    當阻尼比確定時,地震動對應(yīng)的加速度彈性反應(yīng)譜的形狀即確定?;谧V形的地震動強度參數(shù)實質(zhì)就是將反應(yīng)譜與結(jié)構(gòu)關(guān)聯(lián)起來,通過譜值、斜率等參數(shù)來反映地震動與結(jié)構(gòu)的信息。目前的相關(guān)研究主要基于加速度反應(yīng)譜Sa,原因有兩點[7]:災(zāi)害危險性分布圖(Hazard Maps)大多基于加速度反應(yīng)譜;許多研究已經(jīng)討論了以加速度反應(yīng)譜作為地震動強度參數(shù)的優(yōu)缺點。將反應(yīng)譜與結(jié)構(gòu)特性相結(jié)合的地震動強度參數(shù)是研究與應(yīng)用的主要方向。

    3.1.1 基于結(jié)構(gòu)非線性的地震動強度參數(shù)

    Cordova等[8]在2000年提出了一個考慮非線性結(jié)構(gòu)周期延長的地震動強度參數(shù):

    式中,Tf是一個大于T1的值,一般情況下可采用α=0.5和Tf/T1=2,該參數(shù)可以降低增量動力分析曲線的整體離散性。

    在此基礎(chǔ)上,Mehanny等[9]在 2009年提出了:

    楊成等[10]在2010年做了兩點修正:將Tf/T1的比值由2改為由靜力推覆分析確定;將α取0.5改為隨增量動力分析加載級別控制的漸變函數(shù):

    式中,Sa,i(T1)表示第i級增量動力分析地震輸入下的結(jié)構(gòu)基本周期對應(yīng)的彈性反應(yīng)譜值;Sa,475y(T1)表示50年超越概率10%水平下對應(yīng)的規(guī)范設(shè)計反應(yīng)譜值;Sa,max(T1)表示動力失穩(wěn)前增量動力分析最大響應(yīng)譜值。

    該函數(shù)對增量動力分析曲線收斂性有一定的提高。

    上述學者在考慮結(jié)構(gòu)非線性周期延長均采用了延長單個周期所對應(yīng)的反應(yīng)譜譜值。對此,Bojórquez[11]在 2010 年提出了用多個譜值點反映周期延長段譜形的地震動強度參數(shù):

    式中,

    式中,Tn取2~2.5倍的基本周期,α在0~1之間取值。

    該地震動強度參數(shù)對增量動力分析曲線離散程度的降低明顯優(yōu)于前者提出的地震動強度參數(shù)。

    3.1.2 基于高階振型的地震動強度參數(shù)

    隨著結(jié)構(gòu)高度的增加,周期隨之增加,高階振型影響不可忽略。Asgarian等[12]在2012年提出了基于前三階振型的地震動強度參數(shù):

    式中,τa,τb,τc為結(jié)構(gòu)一個方向的前三階周期;β=γ=1/3。

    周穎等[4]在2013年提出了更進一層次的地震動強度參數(shù):

    與Asgarian等提出的地震動強度參數(shù)相比,α,β,γ用振型質(zhì)量參與系數(shù)的權(quán)重表示,反映出各階振型的影響。此外,Lu等[13]在2013年也提出了基于多階振型的地震動強度參數(shù):

    推導了階數(shù)n與基本周期T1的關(guān)系式,等效地對待了各階振型。

    上述基于高階振型影響的地震動強度參數(shù)在對中長周期結(jié)構(gòu)的增量動力分析曲線離散程度減小方面均得到了很好的驗證。

    3.1.3 基于非線性與高階振型的地震動強度參數(shù)

    除了單方面基于非線性影響和高階振型影響的地震動強度參數(shù)外,有學者開始提出了綜合考慮兩者的地震動強度參數(shù)。2009年,Bianchini等[14]提出了:

    與 Bojórquez 提出的 INp中的 Sa,avg(T1,…,Tn)不同,式(8a)中的Ti均不代表結(jié)構(gòu)的i階自振周期,包括T1,并非結(jié)構(gòu)基本周期。該組周期點只是普通的一組周期點,但范圍涵蓋了結(jié)構(gòu)高階周期點與結(jié)構(gòu)非線性周期延長的周期點。通過證明,對于多自由度系統(tǒng)Sa,avg可用10個周期點,并用系數(shù)與基本周期T(1)乘積表示,可將式(8a)寫為:

    式中,Kl考慮了高階振型,為選取周期點的下限;Ku考慮了非線性周期延長,為選取周期點的上限,該系數(shù)的值可根據(jù)結(jié)構(gòu)延性系數(shù)μ確定。

    經(jīng)過分析,Sa,avg的有效性與充分性均優(yōu)于Sa(T1)。

    Lin等[2]在2010年提出了針對考慮非線性和高階振型影響的兩個地震動強度參數(shù),通過判斷結(jié)構(gòu)基本周期T1的范圍來選擇:

    式中,SN1考慮了結(jié)構(gòu)非線性,適用于T1<1.5 s的框架,SN2考慮了高階振型,適用于T1>1.5 s的框架;C為考慮結(jié)構(gòu)非線性的周期延長系數(shù),根據(jù)美國土木工程學會標準7荷載規(guī)范ASCE/SEI 7 -05 建議取為 1.5[15]。

    Adeli等[16]在2012年提出了一種用加速度反應(yīng)譜面積表示的地震動強度參數(shù):區(qū)間范圍在[1.2Tm,1.5T1]的譜面積。其中,Tm為結(jié)構(gòu)振型質(zhì)量參與系數(shù)超過95%對應(yīng)的周期,該區(qū)間包含了高階振型以及非線性周期延長段。反應(yīng)譜的面積可認為是一種譜強度,面積也可看作是多個譜值的疊加。

    從上述研究成果中可以發(fā)現(xiàn),目前基于加速度反應(yīng)譜譜值的地震動強度參數(shù)研究開展的較為廣泛,并逐步從單一譜值點的形式發(fā)展至多個譜值點。一般可以認為,隨著地震動強度參數(shù)表達形式中反應(yīng)譜譜值點的增加,所描述的反應(yīng)譜譜形更為準確,同時反映出的反應(yīng)譜與結(jié)構(gòu)的信息也更加全面,其計算結(jié)果將更為合理。

    盡管基于反應(yīng)譜譜形相關(guān)研究較多,但目前仍缺乏全面反映結(jié)構(gòu)特性與地震動信息的參數(shù)。在考慮結(jié)構(gòu)非線性影響時,也僅從結(jié)構(gòu)剛度降低、周期延長的角度出發(fā),并未考慮阻尼等因素。這仍將會成為今后研究的主要方向。

    3.2 基于位移的地震動強度參數(shù)

    Luco等[17]在2007年提出了6個基于彈性位移反應(yīng)譜和非彈性位移反應(yīng)譜的地震動強度參數(shù):

    此外,2013年Lu等[18]在不同周期結(jié)構(gòu)下對峰值位移(Peak Ground Displacement,PGD)以及位移反應(yīng)譜在結(jié)構(gòu)基本周期的譜值Sd(T1)進行了比較,研究發(fā)現(xiàn)Sd(T1)與Sa(T1)相比在離散程度上并沒有降低,反而在低周期結(jié)構(gòu)效果不如Sa(T1)。峰值位移在長周期結(jié)構(gòu)有效性與IM1E&2E相近,優(yōu)于Sa(T1),但同樣也不適用短周期結(jié)構(gòu)。

    3.3 基于速度的地震動強度參數(shù)

    目前基于速度的地震動強度參數(shù)研究主要是針對已有速度型地震動強度參數(shù)進行對比。Wang等[19]在2012年通過比較,建議用速度譜強度(Velocity Spectrum Intensity,VSI)與峰值速度(Peak Ground Velocity,PGV)作為橋基耦合系統(tǒng)的地震動強度參數(shù)。Asgarian等[12]在2012年對速度反應(yīng)譜在結(jié)構(gòu)基本周期的譜值Sv(T1)進行了比較,認為近場地震作用最合適的地震動強度參數(shù)為 Sv(T1)。此外 Lu 等[13,18]在 2013 年也對地震動峰值速度作了對比研究,發(fā)現(xiàn)峰值速度對長周期結(jié)構(gòu)相關(guān)性與離散程度均表現(xiàn)優(yōu)異。

    從上述研究表明,基于速度的地震動強度參數(shù)相對于傳統(tǒng)的基于加速度的地震動強度參數(shù)有一定的優(yōu)勢,且日本也采用峰值速度作為烈度的物理標準。但我國抗震設(shè)計采用峰值加速度作為主要參數(shù)指標,相應(yīng)峰值速度的研究尚存空缺??傮w來說,基于位移和速度的地震動強度參數(shù)相對基于加速度反應(yīng)譜的地震動強度參數(shù)缺乏系統(tǒng)的研究,應(yīng)值得學者關(guān)注。

    3.4 基于彈塑性反應(yīng)譜的地震動強度參數(shù)

    楊成等[20]在2008年用彈塑性反應(yīng)譜譜值NSa(T1,5%)來代替Sa(T1)作為地震動強度參數(shù),提出改進的MIDA(Modified IDA)方法。楊成采用了Buffalo大學的NSPECTRA2.0程序計算彈塑性反應(yīng)譜,并定義了彈塑性反應(yīng)譜計算參數(shù)ξy。通過算例發(fā)現(xiàn),用NSa(T1,5%)作為地震動強度參數(shù)得到的增量動力分析曲線離散性明顯小于采用彈性反應(yīng)譜Sa(T1)。

    目前研究難點主要在于彈塑性反應(yīng)譜的確定,基于彈塑性反應(yīng)譜的地震動強度參數(shù)的研究并未得到很好的開展。

    4 向量型地震動強度參數(shù)

    由于一維標量型地震動強度參數(shù)存在著局限性,Baker[21]在2005 年將地震動強度參數(shù)從一維發(fā)展到了二維,提出了一種向量型地震動強度參數(shù) <Sa(T1),ε>。其中,ε為地震動加速度反應(yīng)譜譜值與衰減模型預(yù)測得到的加速度反應(yīng)譜均值之間的差值,同樣也是一種反映反應(yīng)譜譜形的參數(shù)。通過與震級M和震中距R進行對比,ε對分析結(jié)果的影響要更大。Baker還建議在選取地震動進行分析時要考慮ε。

    隨后,2008 年 Baker[22]在 Cordova 等[8]基礎(chǔ)上提出了向量型地震動強度參數(shù) <Sa(T1),RT1,T2>,其中,RT1,T2為T1和T2周期對應(yīng)的反應(yīng)譜譜值的比值。研究表明,傳統(tǒng)的Sa(T1)不能很好描述脈沖型地震動,而第二參數(shù)RT1,T2的加入使得地震動強度參數(shù)針對普通型和脈沖型地震動進行結(jié)構(gòu)響應(yīng)預(yù)測時沒有明顯區(qū)別,使地震動強度參數(shù)更具有充分性。

    2011 年 Bojórquez[23]在標量型 INp[11]基礎(chǔ)上也提出了新的向量型地震動強度參數(shù)<Sa(T1),Np>,其中,Np表達式與3.1.1節(jié)中式(4)相同。研究證明,在相同標準誤差e之下,向量型地震動強度參數(shù)所需的地震動數(shù)量n明顯小于標量型地震動強度參數(shù),其中近場脈沖型地震動的數(shù)量向量型僅為標量型的10%。由此可見,選用向量型地震動強度參數(shù)可提高計算分析的效率。

    隨后在2012 年,Bojórquez等[7]又作了進一步的探討,將向量型地震動強度參數(shù)劃分為峰值型(< Sa(T1),PGA > 、< Sa(T1),PGV >)、累積破壞型(<Sa(T1),tD>、<Sa(T1),ID>)和基于反應(yīng)譜譜形 (<Sa(T1),RT1,T2> 、<Sa(T1),Np>)三大類進行研究比較。向量型地震動強度參數(shù)第一參數(shù)均為Sa(T1),tD為地震動持時,ID=,tE為地震動總時長。結(jié)果表明基于反應(yīng)譜譜形的向量型地震動強度參數(shù) < Sa(T1),RT1,T2> 和 < Sa(T1),Np> 有效性較高。

    近年有學者在原有的標量型地震動強度參數(shù)基礎(chǔ)上進行組合。Jalayer等[24]在2012年對向量型地 震 動 強 度 參 數(shù) <Sa(T1),Sa(T2) > 、<Sa(T1),ε>、<PGA,M >進行了對比研究,其中考慮兩階振型的 <Sa(T1),Sa(T2)>充分性較好。

    Gehl等[25]在2013年研究了50個標量型地震動強度參數(shù)值(如峰值加速度、峰值速度、Arias強度等),對結(jié)構(gòu)屈服與倒塌兩種狀況進行比較,結(jié)果未發(fā)現(xiàn)有標量型地震動強度參數(shù)同時屬于這兩種狀況的最佳選擇。通過用向量型地震動強度參數(shù)將原有標量型地震動強度參數(shù)進行組合,汲取單個標量型地震動強度參數(shù)優(yōu)勢。研究表明,二維的向量型地震動強度參數(shù)有效性均要高于一維標量型地震動強度參數(shù),其可降低地震危險性分析的不確定性。

    向量型地震動強度參數(shù)相比標量型地震動強度參數(shù)可提供更多的信息,并可將各影響因素分開考慮,相對于綜合型的單一標量型地震動強度參數(shù)其意義更為明確。向量型地震動強度參數(shù)的充分性和有效性已經(jīng)在研究中得到了驗證,其具有一定的優(yōu)勢。目前向量型地震動強度研究剛剛開展,研究空間較大。

    5 結(jié)語

    降低增量動力分析曲線離散性,可減少地震動樣本數(shù)量;提高地震動強度參數(shù)充分性,可簡化地震危險性分析的積分式[22]。因此深入地震動強度參數(shù)研究對基于性能的地震工程分析的效率提高具有重要意義。目前無論是標量型地震動強度參數(shù)還是向量型地震動強度參數(shù),基于反應(yīng)譜譜形的地震動強度參數(shù)仍是研究的趨勢。隨著結(jié)構(gòu)高度增加、體系復雜化,選取地震動強度參數(shù)時需要考慮更多的影響因素。而向量型地震動強度參數(shù)相比標量型地震動強度參數(shù)可提供的信息更多,具有更大的發(fā)展前景。此外,由于目前的研究采用的地震動輸入均為單向,并不能完全滿足工程要求,因此多向地震動輸入的地震動強度參數(shù)研究也將是今后應(yīng)重視的研究內(nèi)容。

    [1] Federal Emergency Management Agency.Recommended seismic design criteria for new steel momentframe buildings[R].Washington DC:2000.

    [2] Lin L,Naumoski N,Saatcioglu M,et al.Improved intensity measures for probabilistic seismic demand analysis part 1:development of improved intensity measures[J].Canadian Journal of Civil Engineering,2010,38(1):79-88.

    [3] 周穎,呂西林,卜一.增量動力分析法在高層混合結(jié)構(gòu)性能評估中的應(yīng)用.同濟大學學報(自然科學版),2010,38(2):183-193.Zhou Ying,Lu Xilin,Bu Yi.Application of incremental dynamic analysis to seismic evaluation of hybrid structure[J].Journal of Tongji University(Natural Science Edition),2010,38(2):183-193.(in Chinese)

    [4] 周穎,蘇寧粉,呂西林.高層建筑結(jié)構(gòu)增量動力分析的地震動強度參數(shù)研究[J].建筑結(jié)構(gòu)學報,2013,34(2):53-60.Zhou Ying,Su Ningfen,Lu Xilin.Study on intensity measure of incremental dynamic analysis for high-rise structures[J].Journal of Building Structures,2013,34(2):53-60.(in Chinese)

    [5] Bazzurro P,Cornell C A,Shome N,et al.Three proposals for characterizing MDOF nonlinear seismic response[J].Journal of Structural Engineering,1998,124(11):1281-1289.

    [6] Vamvatsikos D,Cornell C A.Incremental dynamic analysis[J].Earthquake Engineering and Structural Dynamics,2002,31(3):491-514.

    [7] Bojórquez E,Iervolino I,Reyes-Salazar A,et al.Comparing vector-valued intensity measures for fragility analysis of steel frames in the case of narrow-band ground motions[J].Engineering Structures,2012,45:472-480.

    [8] Cordova P P,Deierlein G G,Mehanny S S F,et al.Development of a two-parameter seismic intensity measure and probabilistic assessment procedure[C].Proceedings of the 2nd US-Japan Workshop on Performance-Based Seismic Design Methodology for Reinforced Concrete Building Structures,2000.

    [9] Mehanny S S F.A broad-range power-law form scalar-based seismic intensity measure[J].Engineering Structures,2009,31(7):1354-1368.

    [10] 楊成,潘毅,趙世春,等.烈度指標對IDA曲線離散性的影響[J].工程力學,2010,27(增刊 I):68-72.Yang Cheng,Pan Yi,Zhao Shichun,et al.Effect of seismic intensity measure function to the IDA curves dispesion[J].Engineering Mechanics,2010,27(S I):68-72.(in Chinese)

    [11] Bojórquez E,Iervolino I.A Spectral shape-based scalar ground motion intensity measure for maximum and cumulative structural demands[C].Proceedings.of 14th European Conference on Earthquake Engineering,2010.

    [12] Asgarian B,Khazaee H,Mirtaheri M.Performance evaluation of different types of steel moment resisting frames subjected to strong ground motion through incremental dynamic analysis[J].International Journal of Steel Structures,2012,12(3):363-379.

    [13] Lu X,Ye L P,Lu X Z,et al.An improved ground motion intensity measure for super high-rise buildings[J].Science China Technological Sciences,2013,56(6):if25-if33.

    [14] Bianchini M,Diotallevi P P,Baker J W.Prediction of inelastic structural response using an average of spectral accelerations[C].10th International Conference on Structural Safety and Reliability(ICOSSAR09),Osaka,Japan.2009.

    [15] ASCE/SEI 7-05.Minimum design loads for buildings and other structures[C].American Society of Civil Engineers,2006.

    [16] Adeli M M,Banazadeh M,Deylami A,et al.Introducing a new spectral intensity measure parameter to estimate the seismic demand of steel moment-resisting frames using bayesian statistics[J].Advances in Structural Engineering,2012,15(2):231-248.

    [17] Luco N,Cornell C A.Structure-specific scalar intensity measures for near-source and ordinary earthquake ground motions[J].Earthquake Spectra,2007,23(2):357-392.

    [18] Lu X,Lu X Z,Guan H,et al.Comparison and selection of ground motion intensity measures for seismic design of super high-rise buildings[J].Advances in Structural Engineering,2013,16(7):1249-1262.

    [19] Wang Z,Osorio L,Padgett J E.Optimal intensity measures for probabilistic seismic response analysis of bridges on liquefiable and non-liquefiable soils[C].Proceedings.of the 2012 Structures Congress,2012.

    [20] 楊成,徐騰飛,李英民,等.應(yīng)用彈塑性反應(yīng)譜對IDA方法的改進研究[J].地震工程與工程震動,2008,28(4):64-69.Yang Cheng,Xu Tengfei,Li Yingmin,et al.Research on modified IDA method by inelastic response spectrum[J].Journal of Earthquake Engineering and Engineering Vibration,2008,28(4):64-69.(in Chinese)

    [21] Baker J W,Cornell C A.A vector-valued ground motion intensity measure consisting of spectral acceleration and epsilon[J].Earthquake Engineering and Structural Dynamics,2005,34(10):1193-1217.

    [22] Baker J W,Cornell C A.Vector-valued intensity measures for pulse-like near-fault ground motions[J].Engineering Structures,2008,30(4):1048-1057.

    [23] Bojórquez E,Iervolino I.Spectral shape proxies and nonlinear structural response[J].Soil Dynamics and Earthquake Engineering,2011,31(7):996-1008.

    [24] Jalayer F,Beck J L,Zareian F.Analyzing the sufficiency of alternative scalar and vector intensity measures of ground shaking based on information theory[J].Journal of Engineering Mechanics,2011,138(3):307-316.

    [25] Gehl P,Seyedi D M,Douglas J.Vector-valued fragility functions for seismic risk evaluation[J].Bulletin of Earthquake Engineering,2013,11(2):365-384.

    猜你喜歡
    標量震動振型
    關(guān)于模態(tài)綜合法的注記
    力學與實踐(2022年5期)2022-10-21 08:10:34
    縱向激勵下大跨鋼桁拱橋高階振型效應(yīng)分析
    震動減脂儀可以減肥?
    科學大眾(2020年12期)2020-08-13 03:22:32
    一種高效的橢圓曲線密碼標量乘算法及其實現(xiàn)
    塔腿加過渡段輸電塔動力特性分析
    水電工程場地地震動確定方法
    一種靈活的橢圓曲線密碼并行化方法
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    結(jié)構(gòu)振型幾何辨識及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    人工合成最不利地震動
    婷婷色综合www| 日韩一卡2卡3卡4卡2021年| 午夜免费观看性视频| 一级爰片在线观看| 建设人人有责人人尽责人人享有的| 王馨瑶露胸无遮挡在线观看| 精品卡一卡二卡四卡免费| 卡戴珊不雅视频在线播放| 精品少妇久久久久久888优播| 中国国产av一级| 欧美日韩视频高清一区二区三区二| 观看美女的网站| 国产成人a∨麻豆精品| 国精品久久久久久国模美| 免费在线观看视频国产中文字幕亚洲 | 观看美女的网站| 免费黄色在线免费观看| 一区二区日韩欧美中文字幕| 亚洲av国产av综合av卡| 黑人欧美特级aaaaaa片| 人妻系列 视频| 成年女人在线观看亚洲视频| 国产女主播在线喷水免费视频网站| 成年女人毛片免费观看观看9 | 性少妇av在线| 亚洲欧美清纯卡通| 日本欧美视频一区| 日韩一本色道免费dvd| 熟女少妇亚洲综合色aaa.| 国产精品久久久久成人av| 久久精品国产鲁丝片午夜精品| 亚洲精品美女久久久久99蜜臀 | 国产成人精品久久二区二区91 | 国产综合精华液| 国产日韩欧美视频二区| 亚洲在久久综合| 免费人妻精品一区二区三区视频| 在线观看美女被高潮喷水网站| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 超色免费av| 国产视频首页在线观看| 国产精品久久久久久精品古装| 欧美黄色片欧美黄色片| 黄片无遮挡物在线观看| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 成人免费观看视频高清| 女人精品久久久久毛片| 日本欧美国产在线视频| 日韩欧美一区视频在线观看| 亚洲人成电影观看| 亚洲色图综合在线观看| 男女午夜视频在线观看| 国产女主播在线喷水免费视频网站| 777米奇影视久久| 亚洲av免费高清在线观看| 久久久久久久亚洲中文字幕| 日韩av在线免费看完整版不卡| 日韩av免费高清视频| 在现免费观看毛片| 18+在线观看网站| 亚洲精品久久久久久婷婷小说| 国产人伦9x9x在线观看 | 国产国语露脸激情在线看| 久久 成人 亚洲| 狠狠精品人妻久久久久久综合| 在线天堂中文资源库| 欧美日韩一级在线毛片| av电影中文网址| 免费黄频网站在线观看国产| 精品人妻偷拍中文字幕| 精品少妇久久久久久888优播| 国产精品蜜桃在线观看| 国产av精品麻豆| 十分钟在线观看高清视频www| 久久韩国三级中文字幕| 日本91视频免费播放| av.在线天堂| 午夜激情久久久久久久| 国产黄色视频一区二区在线观看| 99re6热这里在线精品视频| 免费观看a级毛片全部| 日日爽夜夜爽网站| 精品亚洲乱码少妇综合久久| 国产成人精品无人区| 夫妻午夜视频| 99热全是精品| 性高湖久久久久久久久免费观看| 久久久久视频综合| 男人爽女人下面视频在线观看| 在线观看免费视频网站a站| 日韩av不卡免费在线播放| 国产男人的电影天堂91| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放| 少妇人妻久久综合中文| 最近手机中文字幕大全| 亚洲三级黄色毛片| 欧美日韩成人在线一区二区| 国产精品久久久久久av不卡| 新久久久久国产一级毛片| 久久狼人影院| 欧美少妇被猛烈插入视频| 亚洲熟女精品中文字幕| 亚洲精品一二三| 一级毛片黄色毛片免费观看视频| 少妇人妻 视频| 午夜免费鲁丝| 婷婷成人精品国产| 国产av一区二区精品久久| 成人亚洲欧美一区二区av| 老鸭窝网址在线观看| 久久久亚洲精品成人影院| 国产成人精品久久二区二区91 | 亚洲美女视频黄频| 欧美精品亚洲一区二区| 尾随美女入室| 日本vs欧美在线观看视频| 最近手机中文字幕大全| 国产成人午夜福利电影在线观看| 亚洲一区二区三区欧美精品| 亚洲国产欧美在线一区| 91精品国产国语对白视频| 男人爽女人下面视频在线观看| 少妇熟女欧美另类| 国产精品蜜桃在线观看| 亚洲国产精品一区二区三区在线| 成人影院久久| 男女午夜视频在线观看| 亚洲精品日韩在线中文字幕| 久久久久久久精品精品| 一级爰片在线观看| 国产精品一国产av| 国产精品久久久久成人av| 午夜免费观看性视频| 街头女战士在线观看网站| 高清不卡的av网站| 日韩精品免费视频一区二区三区| 久久精品国产自在天天线| 午夜91福利影院| 国产免费又黄又爽又色| 91国产中文字幕| 久久韩国三级中文字幕| 久久久欧美国产精品| 美女脱内裤让男人舔精品视频| 欧美老熟妇乱子伦牲交| 人人妻人人爽人人添夜夜欢视频| 黄频高清免费视频| 亚洲av电影在线进入| 如日韩欧美国产精品一区二区三区| 精品福利永久在线观看| 美女主播在线视频| 日韩不卡一区二区三区视频在线| 少妇被粗大猛烈的视频| 亚洲一级一片aⅴ在线观看| 天美传媒精品一区二区| 99久久综合免费| 咕卡用的链子| 国产精品久久久久久精品古装| 美女视频免费永久观看网站| av线在线观看网站| 丝袜喷水一区| 伊人久久大香线蕉亚洲五| 亚洲欧美色中文字幕在线| 国产精品国产三级国产专区5o| 成年av动漫网址| 成年女人在线观看亚洲视频| 最近最新中文字幕大全免费视频 | 亚洲成人一二三区av| 高清在线视频一区二区三区| 美女中出高潮动态图| 成人免费观看视频高清| 久久99精品国语久久久| 一边亲一边摸免费视频| 日日摸夜夜添夜夜爱| 一区二区三区激情视频| 精品国产一区二区三区四区第35| 99国产综合亚洲精品| 26uuu在线亚洲综合色| 亚洲国产精品一区三区| 97人妻天天添夜夜摸| 男女无遮挡免费网站观看| 制服丝袜香蕉在线| 一级片'在线观看视频| 国产精品嫩草影院av在线观看| 秋霞伦理黄片| 色婷婷av一区二区三区视频| 超碰97精品在线观看| 国产成人精品在线电影| 成人漫画全彩无遮挡| 久久精品国产鲁丝片午夜精品| 永久免费av网站大全| 精品少妇久久久久久888优播| 伊人久久国产一区二区| 一级,二级,三级黄色视频| 高清不卡的av网站| 国产精品国产av在线观看| 免费播放大片免费观看视频在线观看| 九色亚洲精品在线播放| 久久久国产精品麻豆| 亚洲内射少妇av| 国产精品一二三区在线看| 国产成人精品在线电影| 国产片内射在线| 国产爽快片一区二区三区| 精品福利永久在线观看| 国精品久久久久久国模美| 在线观看一区二区三区激情| 国产欧美日韩综合在线一区二区| 狂野欧美激情性bbbbbb| 你懂的网址亚洲精品在线观看| 在线观看国产h片| 久久精品国产a三级三级三级| 国产白丝娇喘喷水9色精品| 亚洲一级一片aⅴ在线观看| 三级国产精品片| 99精国产麻豆久久婷婷| 免费观看性生交大片5| 香蕉精品网在线| a级片在线免费高清观看视频| 欧美日韩综合久久久久久| 汤姆久久久久久久影院中文字幕| av免费在线看不卡| 成人影院久久| 我要看黄色一级片免费的| 99久久人妻综合| www.av在线官网国产| 日本午夜av视频| 丰满迷人的少妇在线观看| 久久久久国产精品人妻一区二区| 日韩电影二区| av.在线天堂| 色婷婷av一区二区三区视频| 国产1区2区3区精品| 国产成人欧美| www.av在线官网国产| 婷婷色综合大香蕉| 高清不卡的av网站| 欧美成人午夜免费资源| 精品人妻偷拍中文字幕| 国产精品久久久久成人av| 99re6热这里在线精品视频| 一级片'在线观看视频| 女的被弄到高潮叫床怎么办| 国产黄色免费在线视频| 久久久久网色| 国产精品av久久久久免费| 亚洲情色 制服丝袜| 免费黄频网站在线观看国产| 国产国语露脸激情在线看| 久久久久国产精品人妻一区二区| 制服丝袜香蕉在线| 午夜老司机福利剧场| 亚洲男人天堂网一区| 精品午夜福利在线看| 久久久久久人人人人人| 国产极品天堂在线| 亚洲欧美中文字幕日韩二区| 黄频高清免费视频| 亚洲国产最新在线播放| 久久久久国产网址| 人人澡人人妻人| 男女无遮挡免费网站观看| 国产国语露脸激情在线看| 国产精品99久久99久久久不卡 | 国产熟女午夜一区二区三区| 成年人午夜在线观看视频| 看免费成人av毛片| 国产一区二区 视频在线| 国产成人精品一,二区| 99热国产这里只有精品6| 少妇 在线观看| 色网站视频免费| 久久久久人妻精品一区果冻| 久久久久久久亚洲中文字幕| 性少妇av在线| 乱人伦中国视频| 国产黄频视频在线观看| 国产黄色视频一区二区在线观看| 大话2 男鬼变身卡| 亚洲精品国产av成人精品| 久久久久国产网址| 观看美女的网站| 成年动漫av网址| 菩萨蛮人人尽说江南好唐韦庄| 性色av一级| 国产福利在线免费观看视频| 一本大道久久a久久精品| 国产一区二区三区综合在线观看| 青春草亚洲视频在线观看| tube8黄色片| 免费在线观看黄色视频的| 日韩一区二区三区影片| 丝袜人妻中文字幕| 在线 av 中文字幕| 亚洲第一区二区三区不卡| 黄色毛片三级朝国网站| 免费大片黄手机在线观看| 久久ye,这里只有精品| 免费少妇av软件| 国产麻豆69| 亚洲精品在线美女| av在线观看视频网站免费| 国产精品久久久久久精品电影小说| 亚洲经典国产精华液单| 美女福利国产在线| 精品人妻熟女毛片av久久网站| 男女边吃奶边做爰视频| 日韩欧美一区视频在线观看| 高清av免费在线| 丁香六月天网| 桃花免费在线播放| 中文字幕制服av| 麻豆av在线久日| 国产精品嫩草影院av在线观看| 80岁老熟妇乱子伦牲交| 99香蕉大伊视频| 99热全是精品| 人体艺术视频欧美日本| 亚洲综合色网址| 国产在线一区二区三区精| 欧美精品一区二区免费开放| 建设人人有责人人尽责人人享有的| 久久免费观看电影| 免费日韩欧美在线观看| 一级毛片黄色毛片免费观看视频| 精品久久久精品久久久| 午夜福利在线免费观看网站| 亚洲欧美一区二区三区久久| 欧美最新免费一区二区三区| 午夜日韩欧美国产| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 午夜福利在线免费观看网站| 精品少妇内射三级| 色网站视频免费| 精品少妇内射三级| 夫妻性生交免费视频一级片| 在线观看免费高清a一片| 日本黄色日本黄色录像| 久久午夜综合久久蜜桃| a级片在线免费高清观看视频| 国产精品 欧美亚洲| 日韩视频在线欧美| 一边亲一边摸免费视频| 成年动漫av网址| 777米奇影视久久| 波多野结衣av一区二区av| 国产成人91sexporn| 九九爱精品视频在线观看| 日本av免费视频播放| 国产日韩欧美视频二区| 91国产中文字幕| 国产日韩欧美视频二区| 免费人妻精品一区二区三区视频| 建设人人有责人人尽责人人享有的| 亚洲精品日韩在线中文字幕| 国产乱来视频区| 电影成人av| 国产精品久久久久久精品电影小说| 亚洲欧美一区二区三区黑人 | 婷婷色麻豆天堂久久| 亚洲欧美一区二区三区国产| 天天躁夜夜躁狠狠久久av| 久久久久精品性色| 精品人妻偷拍中文字幕| 午夜激情av网站| 国产女主播在线喷水免费视频网站| 国产精品熟女久久久久浪| 寂寞人妻少妇视频99o| 欧美 亚洲 国产 日韩一| 国产有黄有色有爽视频| 99久国产av精品国产电影| 国产有黄有色有爽视频| 成年人午夜在线观看视频| 最近最新中文字幕大全免费视频 | 美女福利国产在线| 大陆偷拍与自拍| 一本久久精品| 亚洲av日韩在线播放| 99久国产av精品国产电影| 久久精品国产亚洲av高清一级| 亚洲av免费高清在线观看| 亚洲人成77777在线视频| 午夜福利视频精品| 精品国产一区二区三区久久久樱花| 一级毛片我不卡| 五月开心婷婷网| 国产成人免费无遮挡视频| 黄色一级大片看看| 啦啦啦视频在线资源免费观看| 国产在线免费精品| 久久精品国产综合久久久| 国产高清不卡午夜福利| 我要看黄色一级片免费的| 久久精品国产亚洲av高清一级| 欧美国产精品va在线观看不卡| 久久人妻熟女aⅴ| 国产精品熟女久久久久浪| 高清黄色对白视频在线免费看| 久久久久网色| 18在线观看网站| 一区二区三区激情视频| 最黄视频免费看| 熟妇人妻不卡中文字幕| 桃花免费在线播放| 久久狼人影院| 国产日韩一区二区三区精品不卡| 人人妻人人爽人人添夜夜欢视频| 人人澡人人妻人| 少妇被粗大的猛进出69影院| 波多野结衣一区麻豆| 高清黄色对白视频在线免费看| 男女高潮啪啪啪动态图| 免费看av在线观看网站| 最近最新中文字幕大全免费视频 | 精品人妻熟女毛片av久久网站| 久久精品久久精品一区二区三区| 国产在线免费精品| 9191精品国产免费久久| 制服丝袜香蕉在线| 黄色视频在线播放观看不卡| 欧美亚洲 丝袜 人妻 在线| 在线观看www视频免费| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 成人国产av品久久久| 一级爰片在线观看| 黑人猛操日本美女一级片| 久久av网站| 综合色丁香网| 五月开心婷婷网| 午夜老司机福利剧场| 男女免费视频国产| 久久久久网色| 天天躁夜夜躁狠狠久久av| 建设人人有责人人尽责人人享有的| 午夜av观看不卡| 国产精品蜜桃在线观看| 色婷婷av一区二区三区视频| 亚洲国产看品久久| 国产视频首页在线观看| 国产在线免费精品| 亚洲四区av| 1024视频免费在线观看| 中文字幕制服av| 性少妇av在线| 精品一品国产午夜福利视频| 综合色丁香网| 亚洲av电影在线进入| 2021少妇久久久久久久久久久| 国产成人精品久久二区二区91 | 久久人妻熟女aⅴ| 午夜日韩欧美国产| 高清黄色对白视频在线免费看| 国产精品免费视频内射| 亚洲精品,欧美精品| 久久久久久久国产电影| 18禁裸乳无遮挡动漫免费视频| 亚洲国产看品久久| 丰满迷人的少妇在线观看| 女人久久www免费人成看片| 精品久久久精品久久久| 老鸭窝网址在线观看| 韩国高清视频一区二区三区| 久久这里只有精品19| 久久久国产欧美日韩av| 色哟哟·www| 免费在线观看完整版高清| 中文字幕av电影在线播放| 黑丝袜美女国产一区| 老司机影院成人| 2022亚洲国产成人精品| 国产爽快片一区二区三区| 色吧在线观看| 少妇的逼水好多| 观看美女的网站| 最近中文字幕2019免费版| 国产片特级美女逼逼视频| av在线观看视频网站免费| 精品99又大又爽又粗少妇毛片| 爱豆传媒免费全集在线观看| 亚洲人成网站在线观看播放| 亚洲,一卡二卡三卡| 精品国产乱码久久久久久男人| 最近2019中文字幕mv第一页| 99香蕉大伊视频| 精品亚洲成国产av| 久久久精品94久久精品| 五月开心婷婷网| 不卡av一区二区三区| 国产精品国产三级专区第一集| 亚洲美女搞黄在线观看| 亚洲综合色网址| 欧美黄色片欧美黄色片| 一级毛片我不卡| 亚洲一区二区三区欧美精品| 久久久久久久国产电影| 欧美日韩亚洲国产一区二区在线观看 | 成人影院久久| 色94色欧美一区二区| av女优亚洲男人天堂| 大香蕉久久成人网| 十八禁高潮呻吟视频| 免费av中文字幕在线| 日本黄色日本黄色录像| 超碰成人久久| 亚洲色图 男人天堂 中文字幕| 精品亚洲成a人片在线观看| 色哟哟·www| 欧美日本中文国产一区发布| 五月天丁香电影| 国产精品一国产av| 午夜福利视频在线观看免费| 美女xxoo啪啪120秒动态图| 色94色欧美一区二区| 亚洲一区中文字幕在线| 日本色播在线视频| 日本欧美视频一区| 亚洲精品久久午夜乱码| 久久人妻熟女aⅴ| 妹子高潮喷水视频| 美女国产视频在线观看| 热99久久久久精品小说推荐| 午夜日韩欧美国产| freevideosex欧美| 亚洲av在线观看美女高潮| 国产男女内射视频| 国产白丝娇喘喷水9色精品| a级毛片黄视频| 五月伊人婷婷丁香| h视频一区二区三区| 亚洲精品久久午夜乱码| 天堂中文最新版在线下载| 韩国高清视频一区二区三区| 桃花免费在线播放| 午夜福利,免费看| 春色校园在线视频观看| 成人手机av| 精品亚洲成a人片在线观看| 成人漫画全彩无遮挡| 叶爱在线成人免费视频播放| 亚洲三区欧美一区| 最近最新中文字幕免费大全7| 久久亚洲国产成人精品v| 午夜91福利影院| 亚洲av国产av综合av卡| 极品少妇高潮喷水抽搐| 丰满迷人的少妇在线观看| 在线天堂最新版资源| 亚洲成人一二三区av| 久久久久国产一级毛片高清牌| 亚洲在久久综合| www日本在线高清视频| 国产精品一国产av| 青春草视频在线免费观看| 午夜免费鲁丝| 欧美精品一区二区大全| 中文字幕制服av| 欧美国产精品一级二级三级| 午夜福利在线免费观看网站| 亚洲三区欧美一区| 亚洲精品乱久久久久久| 国产一区二区三区综合在线观看| 亚洲国产精品成人久久小说| 亚洲美女视频黄频| 1024香蕉在线观看| 99久久精品国产国产毛片| 日本色播在线视频| 成年女人在线观看亚洲视频| 国产成人免费无遮挡视频| 亚洲av欧美aⅴ国产| 久久人人97超碰香蕉20202| 成人影院久久| av.在线天堂| 丝袜在线中文字幕| 午夜久久久在线观看| 性高湖久久久久久久久免费观看| 在线免费观看不下载黄p国产| 精品一区二区免费观看| 女的被弄到高潮叫床怎么办| 欧美中文综合在线视频| 亚洲av电影在线进入| 三级国产精品片| 免费不卡的大黄色大毛片视频在线观看| 精品人妻熟女毛片av久久网站| www日本在线高清视频| 最新中文字幕久久久久| 天天影视国产精品| 观看美女的网站| 国产有黄有色有爽视频| 黄色 视频免费看| 亚洲欧洲日产国产| 国产高清不卡午夜福利| 青春草国产在线视频| 成年女人毛片免费观看观看9 | 麻豆乱淫一区二区| 国产97色在线日韩免费| 99re6热这里在线精品视频| 国产片特级美女逼逼视频| 亚洲欧洲国产日韩| 亚洲,欧美精品.| 在线 av 中文字幕| 久久久久精品久久久久真实原创| 亚洲精品久久成人aⅴ小说| 久久99热这里只频精品6学生| 亚洲欧美精品自产自拍| 韩国精品一区二区三区| 国产高清国产精品国产三级| 中文字幕人妻丝袜制服| 亚洲成色77777| 日韩 亚洲 欧美在线| 少妇人妻久久综合中文| 亚洲成av片中文字幕在线观看 |