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

    預(yù)應(yīng)力箱型梁橋遭受超高車(chē)輛在不同位置撞擊下的動(dòng)態(tài)響應(yīng)

    2016-02-22 00:54:44馮振寧

    田 力, 馮振寧

    (天津大學(xué) 建筑工程學(xué)院, 天津 300072)

    ?

    預(yù)應(yīng)力箱型梁橋遭受超高車(chē)輛在不同位置撞擊下的動(dòng)態(tài)響應(yīng)

    田力,馮振寧

    (天津大學(xué)建筑工程學(xué)院, 天津300072)

    摘要:隨著城鎮(zhèn)化的加快,超高車(chē)輛撞擊城市預(yù)應(yīng)力箱型橋梁上部結(jié)構(gòu)的事故時(shí)有發(fā)生。因此,有必要對(duì)此類(lèi)撞擊事故的動(dòng)態(tài)響應(yīng)進(jìn)行研究。運(yùn)用LS-DYNA建立了精細(xì)化T型梁橋上部結(jié)構(gòu)的三維分離式模型,與罐體模型進(jìn)行碰撞分析。與相應(yīng)試驗(yàn)結(jié)果進(jìn)行了對(duì)比,驗(yàn)證了LS-DYNA模擬此類(lèi)撞擊的正確性。運(yùn)用LS-DYNA建立了精細(xì)化預(yù)應(yīng)力箱型梁橋上部結(jié)構(gòu)的三維分離式模型,分別建立了三種工況下車(chē)輛與橋梁上部結(jié)構(gòu)碰撞的耦合模型,由計(jì)算結(jié)果可以看出三種工況下橋梁損傷的位置和狀況大致相似,撞擊位置距離支座越近,橋梁損傷就越嚴(yán)重。結(jié)果表明:撞擊位置距離支座越近,碰撞力越大,撞擊區(qū)域位移就越大。因此,超高車(chē)輛撞擊預(yù)應(yīng)力箱梁考慮最不利情況時(shí),應(yīng)研究車(chē)輛撞擊橋梁支座處的破壞情況。

    關(guān)鍵詞:箱型梁橋;汽車(chē)撞擊;分離式模型;動(dòng)態(tài)響應(yīng)

    在繁華的現(xiàn)代大都市中,運(yùn)輸繁忙的立交橋和其它一些跨線城市橋梁經(jīng)常遭受超高車(chē)輛撞擊,造成車(chē)輛和橋梁損傷。特別是最近二十年,隨著中國(guó)經(jīng)濟(jì)技術(shù)的迅猛發(fā)展,城鎮(zhèn)化進(jìn)程的加快,越來(lái)越多的人們涌入大城市中,使得城市立體交通體系成為城市交通的重要組成部分。這一系列的變化改善了城市交通,方便人們的出行,但也致使超高車(chē)輛撞擊城市橋梁上部結(jié)構(gòu)的事故發(fā)生更加頻繁。

    盡管并非每次超高車(chē)輛撞擊橋梁上部結(jié)構(gòu)會(huì)造成車(chē)毀人亡的重大損失,但是橋梁被撞會(huì)使得城市交通受阻、橋梁壽命減少,造成一定的經(jīng)濟(jì)損失。因此,此類(lèi)問(wèn)題的頻繁發(fā)生得到了社會(huì)多方面的重視。英國(guó)Arup公司曾經(jīng)做過(guò)車(chē)撞橋梁上部結(jié)構(gòu)的足尺實(shí)驗(yàn),但并未公布實(shí)驗(yàn)結(jié)果。歐洲規(guī)范EUROCODE 1[1]對(duì)車(chē)撞橋梁上部結(jié)構(gòu)提出了相關(guān)的撞擊力計(jì)算方法,但對(duì)橋梁的動(dòng)態(tài)響應(yīng)沒(méi)有具體研究。我國(guó)的JTG D60-2004《公路橋梁通用設(shè)計(jì)規(guī)范》[2]參考外國(guó)規(guī)范,簡(jiǎn)單規(guī)定汽車(chē)撞擊力的等效撞擊力大小取值(等效撞擊力在車(chē)輛行駛方向取1000 kN,垂直車(chē)輛行駛方向取500 kN,兩個(gè)方向不能同時(shí)考慮),但是對(duì)于不同橋型、不同撞擊位置,規(guī)范并沒(méi)有詳細(xì)劃分。重慶交通大學(xué)的畢晨華等[3]對(duì)T型梁橋遭受超高車(chē)輛撞擊的損傷進(jìn)行有限元分析。清華大學(xué)陸新征等[4~6]對(duì)超高車(chē)輛撞擊橋梁上部結(jié)構(gòu)進(jìn)行了一系列的數(shù)值模擬分析,研究了部分橋型的破壞過(guò)程和損壞機(jī)理,分析了各個(gè)橋梁遭受撞擊后不同的動(dòng)態(tài)響應(yīng)。然而,對(duì)于預(yù)應(yīng)力混凝土薄壁箱梁這一常見(jiàn)橋型的研究還沒(méi)有涉及。而且,由于橋梁上部結(jié)構(gòu)配筋的復(fù)雜性,目前研究橋梁上部結(jié)構(gòu)多采用整體式建模。顯然為了獲取更精確的結(jié)果,建立分離式模型模擬超高車(chē)輛撞擊橋梁上部結(jié)構(gòu)顯得十分必要。

    本文運(yùn)用有限元分析軟件LS-DYNA對(duì)超高車(chē)輛撞擊預(yù)應(yīng)力箱型梁橋的過(guò)程進(jìn)行了數(shù)值模擬。建立了精細(xì)化的預(yù)應(yīng)力混凝土箱梁分離式模型,驗(yàn)證了該軟件數(shù)值模擬的有效性,比較分析了汽車(chē)在不同撞擊位置下,預(yù)應(yīng)力混凝土橋梁的動(dòng)態(tài)響應(yīng)。

    1數(shù)值模型和材料

    1.1汽車(chē)數(shù)值模型

    作為美國(guó)政府部門(mén)汽車(chē)安全的最高主管部門(mén),美國(guó)高速公路安全管理局(NHTSA)在汽車(chē)碰撞實(shí)驗(yàn)分析和汽車(chē)安全性能分析領(lǐng)域具有話語(yǔ)權(quán)。與其合作的“國(guó)家碰撞分析中心”(National Crash Analysis Center, NCAC)擁有世界上最準(zhǔn)確的汽車(chē)碰撞試驗(yàn)資料庫(kù),其官方網(wǎng)站為研究人員免費(fèi)提供可靠的汽車(chē)有限元模型。本文使用的汽車(chē)有限元模型即是NCAC提供的標(biāo)準(zhǔn)雙軸卡車(chē)有限元模型。該車(chē)模型原重為7.17 t,筆者通過(guò)改變相關(guān)參數(shù)將汽車(chē)質(zhì)量提高到8 t,以符合我國(guó)常見(jiàn)的雙軸卡車(chē)重量,有限元模型如圖1所示。

    圖1 汽車(chē)數(shù)值模型

    1.2箱型梁橋數(shù)值模型

    參考文獻(xiàn)[7]建立了典型的預(yù)應(yīng)力混凝土箱型梁橋上部結(jié)構(gòu)數(shù)值模型。該橋梁由5個(gè)箱型梁組成,其中單個(gè)箱梁橫截面尺寸及配筋圖如圖2所示。橋梁標(biāo)準(zhǔn)跨徑為30 m,橋面寬度為15 m,可供四車(chē)道車(chē)輛同時(shí)行駛。橋梁邊界條件采用板式橡膠支座,上部結(jié)構(gòu)與支座定義為自動(dòng)面面接觸,摩擦系數(shù)為0.3[2]。

    圖2 箱梁橫截面尺寸及配筋/mm

    采用分離式建模,鋼筋采用梁?jiǎn)卧獎(jiǎng)澐郑炷敛捎昧骟w實(shí)體單元?jiǎng)澐?。鋼筋和混凝土采用共?jié)點(diǎn)處理,不考慮它們之間的粘結(jié)滑移效應(yīng)[8]。其中,混凝土強(qiáng)度等級(jí)取為C50,單根箱梁預(yù)應(yīng)力鋼筋面積為7848 mm2,預(yù)應(yīng)力值為1200 MPa。梁箍筋每隔100 mm配452 mm2的HRB300鋼筋。其它配筋皆按照交通部專(zhuān)家委員會(huì)編制的《公路橋涵通用圖》布置,并皆為二級(jí)鋼筋,橋梁上部結(jié)構(gòu)有限元模型如圖3所示。

    圖3 單根箱梁有限元模型

    1.3汽車(chē)-箱型梁橋耦合模型

    本文運(yùn)用后處理軟件LS-PREPOST將汽車(chē)和預(yù)應(yīng)力混凝土箱型梁橋的數(shù)值模型耦合在一起。計(jì)算中,汽車(chē)行駛速度取為80 km/h,汽車(chē)車(chē)廂與橋梁上部結(jié)構(gòu)接觸定義為面面接觸。參考?xì)W洲規(guī)范EUROCODE 1[1]的規(guī)定,本文將汽車(chē)超高值取為250 mm(車(chē)輛凈高與橋下凈空的差值),其余設(shè)置皆按照文獻(xiàn)[9]設(shè)置。汽車(chē)與預(yù)應(yīng)力箱型梁橋的耦合模型如圖4所示。

    圖4 汽車(chē)-橋梁上部結(jié)構(gòu)的耦合模型

    1.4材料選擇

    本文普通鋼筋材料采用LS-DYNA中*MAT_PLASTIC_KINEMATIC(MAT3)[10]。鋼筋材料模型的的屈服應(yīng)力可以表示為

    (1)

    采用的普通鋼筋材料參數(shù)如表1。

    表1 鋼筋材料參數(shù)

    預(yù)應(yīng)力鋼筋材料采用LS-DYNA中*MAT_ELASTIC_PLASTIC_THERMAL(MAT4)[10]。利用此種材料“熱脹冷縮”特性,當(dāng)溫度下降時(shí),鋼筋收縮,從而壓縮混凝土,混凝土中產(chǎn)生預(yù)壓應(yīng)力。利用*LOAD_THERMAL_LOAD_CURVE關(guān)鍵字加載溫度荷載。

    混凝土材料采用*MAT_CONCRETE_DAMAGE_REAL3(MAT72)[10]。由于在沖擊荷載作用下混凝土的抗壓強(qiáng)度與靜載作用時(shí)會(huì)有不同程度的提高,所以有必要考慮混凝土的應(yīng)變率效應(yīng)?;炷敛牧蠎?yīng)變率效應(yīng)應(yīng)根據(jù)CEB模型的建議,通常用材料強(qiáng)度動(dòng)力增大系數(shù)DIF(Dynamic Increase Factor)來(lái)表示,DIF為在某一應(yīng)變率下材料的動(dòng)力強(qiáng)度與靜力強(qiáng)度之比。

    混凝土受拉時(shí)按以下公式計(jì)算強(qiáng)度增大系數(shù)TDIF:

    (2)

    混凝土受壓時(shí)按以下公式計(jì)算強(qiáng)度增大系數(shù)CDIF:

    (3)

    2數(shù)值模型計(jì)算驗(yàn)證

    圖5 試驗(yàn)T梁模型

    圖6 試驗(yàn)有限元模型

    圖7 T梁有限元模型及配筋/mm

    本文建立了文獻(xiàn)[11]所做試驗(yàn)對(duì)應(yīng)的有限元模型,其試驗(yàn)?zāi)P秃陀邢拊P头謩e如圖5和圖6所示。該試驗(yàn)是以縮尺試驗(yàn)的方法,用罐體模型撞擊鋼筋混凝土T梁模型,以模擬現(xiàn)實(shí)中罐車(chē)撞擊T型梁橋上部結(jié)構(gòu)的情形。由圖5可以看出,該試驗(yàn)采用的是反位試驗(yàn)。其中,罐體通過(guò)腹部若干點(diǎn)由纜繩約束,抬升至一定高度使之落下,獲得的初速度為2.855 m/s,罐體與T梁的撞擊區(qū)高度(罐體頂面到橋梁底面)為50 mm。在對(duì)應(yīng)的有限元模型中,罐體采用殼單元來(lái)進(jìn)行劃分,其材料模型采用上文提到的*MAT_PLASTIC_KINEMATIC(MAT3),支座采用實(shí)體單元自由劃分,支座材料采用*MAT_RIGID(MAT20)[10]剛體材料,其它材料的選擇與上文一致。采用精細(xì)化建模,建立了與試驗(yàn)?zāi)P鸵恢碌匿摻罨\,鋼筋單元和混凝土單元采用共節(jié)點(diǎn),如圖7所示,材料強(qiáng)度取值與試驗(yàn)一致。T梁與支座、鋼筋混凝土梁接觸皆采用自動(dòng)面面接觸,接觸面間的摩擦系數(shù)取為0.5。將試驗(yàn)通過(guò)加速度積分算出的T梁受撞擊部位最大點(diǎn)的位移時(shí)程與有限元模擬相應(yīng)部位最大點(diǎn)位移時(shí)程結(jié)果對(duì)比,如圖8所示。

    圖8 跨中位移時(shí)程比較

    從圖8可以看出,本文采用LS-DYNA建立的數(shù)值模型計(jì)算結(jié)果與試驗(yàn)測(cè)試結(jié)果吻合較好。無(wú)論T梁頂面還是底面,最大位移值相差不大,位移發(fā)展基本趨勢(shì)大體相似。因此,可以斷定,用有限元軟件LS-DYNA以及采用上文提到的模型建立和材料參數(shù)設(shè)置來(lái)進(jìn)行模擬超高車(chē)輛撞擊橋梁上部結(jié)構(gòu)動(dòng)態(tài)響應(yīng)分析是真實(shí)可靠的。

    3不同撞擊位置的橋梁損傷對(duì)比

    3.1數(shù)值分析初始條件

    本文將計(jì)算過(guò)程分為兩個(gè)階段,第一階段為初應(yīng)力加載階段,第二階段為車(chē)-橋碰撞階段。

    在第一階段中,為減少計(jì)算單元,只建立橋梁上部結(jié)構(gòu)模型,對(duì)橋梁上部結(jié)構(gòu)施加重力、橋面上部荷載以及預(yù)應(yīng)力荷載。為了讓結(jié)構(gòu)能夠快速穩(wěn)定下來(lái),通過(guò)試算,估計(jì)出結(jié)構(gòu)臨界阻尼,加到橋梁模型上。在施加重力荷載時(shí),分兩步加載:荷載首先由零開(kāi)始線性增加至重力加速度大小,加載時(shí)長(zhǎng)為0.5 s,待橋梁上部結(jié)構(gòu)穩(wěn)定之后,保持荷載不變,直到第二階段碰撞結(jié)束。

    本文利用了LS-DYNA的重啟動(dòng)功能,在第二階段中撤掉結(jié)構(gòu)上的阻尼,避免阻尼對(duì)碰撞結(jié)果產(chǎn)生影響。同時(shí),加入車(chē)輛模型,給予車(chē)輛80 km/h的初速度,分別撞擊橋梁跨中處,距橋梁端部10 m處,距橋梁端部5 m處,本文分別稱(chēng)為工況一、工況二、工況三,如圖9所示。

    圖9 三種工況

    3.2不同撞擊位置橋梁損傷對(duì)比

    圖10是三種工況下模型計(jì)算至1.2 s時(shí)的橋梁損傷圖,可以看出,橋梁的局部損傷主要是在橋梁底部和箱梁被撞面另一側(cè)。其中,底部混凝土大面積脫落、鋼筋屈服,這是因?yàn)檐?chē)廂變形后對(duì)箱梁底的沖切所致。而箱梁另一側(cè)的混凝土脫落是因?yàn)橄淞菏軟_擊使得另一面混凝土局部受壓破壞。通過(guò)對(duì)比三種工況下橋梁損傷情況可以發(fā)現(xiàn),工況三混凝土脫落面積最大,工況二次之,工況一最小,說(shuō)明超高車(chē)輛撞擊橋梁的位置不同,對(duì)橋梁造成的損傷也會(huì)不同。撞擊位置離支座距離越近,橋梁受到的破壞就越大,原因是離支座越近,橋梁受到支座的約束越強(qiáng),剛度相對(duì)就越大,使得撞擊更為猛烈。

    圖10 三種工況下橋梁的損傷

    3.3不同撞擊位置下點(diǎn)位移及撞擊力對(duì)比

    在被撞箱梁另一側(cè)取三個(gè)碰撞區(qū)域關(guān)聯(lián)點(diǎn),分別命名為點(diǎn)A、B和C,如圖11所示。其中,A點(diǎn)為工況一箱梁被撞區(qū)域?qū)?yīng)另一側(cè)跨中一點(diǎn),其距離箱梁底面高度為250 mm,距離梁端距離為15 m;B點(diǎn)為為工況二箱梁被撞區(qū)域?qū)?yīng)另一側(cè)跨中一點(diǎn),其距離箱梁底面高度為250 mm,距離較近梁端為10 m;C點(diǎn)為為工況三箱梁被撞區(qū)域?qū)?yīng)另一側(cè)跨中一點(diǎn),其距離箱梁底面高度為250 mm,距離較近梁端為5 m。圖12是三種工況下,三個(gè)點(diǎn)的X方向位移時(shí)程曲線圖,可以看出,三種工況各自在其撞擊區(qū)域?qū)?yīng)的那一點(diǎn)位移較大,其它兩點(diǎn)位移較小。而且,工況三對(duì)應(yīng)點(diǎn)的最大位移﹥工況二對(duì)應(yīng)點(diǎn)最大位移﹥工況一的對(duì)應(yīng)點(diǎn)最大位移。產(chǎn)生這種現(xiàn)象的原因是由于碰撞區(qū)域離端部越近,橋梁受到端部支座的約束越強(qiáng),使得碰撞更為猛烈,致使局部損傷更為嚴(yán)重,表現(xiàn)為位移越大。另外,通過(guò)圖12可以看出,工況一和工況二對(duì)應(yīng)的最大位移時(shí)程曲線出現(xiàn)了反彎的現(xiàn)象,工況三沒(méi)有反彎現(xiàn)象,而且隨著撞擊位置離端部越近,反彎現(xiàn)象越不明顯,原因是因?yàn)榕鲎矃^(qū)域離端部越遠(yuǎn),受端部支座的約束越弱,使得剛度相對(duì)較弱,相對(duì)“柔一點(diǎn)”,因此會(huì)有反彎現(xiàn)象。圖13給出了三種工況下碰撞力的時(shí)程對(duì)比,可以看出三種工況的碰撞力峰值由大到小依次為為工況三、工況二、工況一。這進(jìn)一步的驗(yàn)證了上文提到的,撞擊區(qū)域離端部越近,撞擊更為猛烈。

    圖11 數(shù)據(jù)采集點(diǎn)

    圖12 不同工況下三點(diǎn)位移時(shí)程曲線

    圖13 不同工況下碰撞力時(shí)程曲線

    3.4撞擊過(guò)程的能量轉(zhuǎn)化

    使用有限元單點(diǎn)積分方法可以大幅度降低計(jì)算成本,但關(guān)鍵要看能不能控制沙漏能的大小。因此,查看沙漏能(Hourglass Energy)與總能量(Total Energy)隨時(shí)間變化的數(shù)值,判斷其能量變化是否合理,顯得非常必要。由于三種工況數(shù)值模擬方法一致,故單獨(dú)拿出工況一的能量變化進(jìn)行分析,如圖14所示。該圖顯示了車(chē)橋碰撞過(guò)程中各種能量之間的動(dòng)態(tài)變化走向。由此圖可以看出,總能量基本保持不變,90%以上的動(dòng)能轉(zhuǎn)化為了內(nèi)能,沙漏能占總能量的比例不到10%。因此,可以斷定這次碰撞過(guò)程能量變化是合理的,也證明了這次數(shù)值模擬的正確性。

    圖14 能量變化

    4結(jié)論

    (1)根據(jù)相應(yīng)試驗(yàn)建立的有限元模型數(shù)值分析結(jié)果與試驗(yàn)結(jié)果對(duì)比發(fā)現(xiàn):利用LS-DYNA軟件進(jìn)行超高車(chē)輛撞擊橋梁上部結(jié)構(gòu)的模擬分析方法是可信的。

    (2)不同撞擊位置下,超高車(chē)輛以80 km/h的速度撞擊預(yù)應(yīng)力箱梁上部結(jié)構(gòu)可以造成橋梁底部混凝土大面積脫落、鋼筋屈服,箱梁被撞區(qū)域?qū)?yīng)另一側(cè)混凝土局部受壓破壞,混凝土脫落。

    (3)撞擊位置越是靠近橋梁端部,車(chē)橋碰撞產(chǎn)生的碰撞力越大,局部損傷也就越大,表現(xiàn)為局部位移較大,局部混凝土脫落面積較大。因此,超高車(chē)輛撞擊預(yù)應(yīng)力箱梁考慮最不利情況應(yīng)研究車(chē)輛撞擊橋梁支座處的破壞情況。

    [1]EN 1991-1-4-2005, EUROCODE 1: Actions on Structures[S].

    [2]JTG D60-2004, 公路橋涵設(shè)計(jì)通用規(guī)范[S].

    [3]畢晨華, 張永水. 超高車(chē)輛碰撞橋梁的有限元損傷分析[J]. 重慶交通大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 30(6): 1294-1297.

    [4]陸新征, 張炎圣, 葉列平, 等. 超高車(chē)輛-橋梁上部結(jié)構(gòu)碰撞的破壞模式與荷載計(jì)算[J]. 中國(guó)公路學(xué)報(bào), 2009, 22(5): 60-67.

    [5]陸新征, 張炎圣, 寧?kù)o, 等. 超高車(chē)輛與立交橋梁碰撞的高精度非線性有限元仿真[J].石家莊鐵道學(xué)院學(xué)報(bào), 2007, 20(1): 29-34.

    [6]陸新征,何水濤,黃盛楠. 超高車(chē)輛撞擊橋梁上部結(jié)構(gòu)研究——破壞機(jī)理、設(shè)計(jì)方法和防護(hù)研究[M]. 北京:中國(guó)建筑工業(yè)出版社, 2011.

    [7]閆志剛.鋼筋混凝土及預(yù)應(yīng)力混凝土簡(jiǎn)支梁橋結(jié)構(gòu)設(shè)計(jì)[M]. 北京:機(jī)械工業(yè)出版社, 2009.

    [8]Krauthammer T, Shahriar S, Shanaa H M. Response of reinforced concrete element severe impulsive loads[J]. Journal of Structural Engineering New York N Y, 1990, 116(4): 1061-1079.

    [9]趙海鷗. LS-DYNA動(dòng)力分析指南[M]. 北京: 兵器工業(yè)出版社, 2003.

    [10]LS-DYNA. Keyword User’s Manual[Z]. Livermore, California:Livermore Software Technology Cooporation, 2006.

    [11]何水濤, 陸新征, 盧 嘯, 等. 超高車(chē)輛撞擊鋼筋混凝土T梁橋主梁試驗(yàn)研究[J]. 蘭州交通大學(xué)學(xué)報(bào), 2011, 30(6): 20-25.

    Dynamic Response of the Pre-stressed Box Girder Bridge to

    Over-high Truck Impact Under Different Impact Position

    TIANLi,FENGZhen-ning

    (School of Civil Engineering, Tianjin University, Tianjin 300072, China)

    Abstract:With the acceleration of the urbanization, the accidents of the over-high truck impacting the pre-stressed box girder bridge superstructure have occurred. Therefore, it is necessary to study the dynamic response of this kind of the impact. A three-dimensional separation model of the fine T bridge superstructure is established by LS-DYNA, and the collision analysis of tank and the T bridge is made. The accuracy of the LS-DYNA simulation was verified by comparison with the experiments. The three-dimensional separation model of the fine pre-stressed box girder bridge superstructure is established by LS-DYNA, and the coupling models of the vehicle and bridge superstructure under three working conditions are established respectively. Through the calculation results can be seen in the bridge damage location and condition are roughly similar. The more close the impact position to the bearing, the more serious damage to the bridge. The results show that the more close the impact position to the bearing, the greater the impact force and the greater the displacement of the impact area. It can be Concluded that the impact of the high vehicle impact on the pre-stressed concrete box girder, the most unfavorable situation should be considered to study the failure of the vehicle impact bridge bearing.

    Key words:box girder bridge; vehicle impact; discrete model; dynamic response

    中圖分類(lèi)號(hào):U447; U441+.3

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):2095-0985(2016)01-0018-06

    基金項(xiàng)目:國(guó)家自然科學(xué)基金 (51178310;50678116)

    作者簡(jiǎn)介:田力(1970-),男,天津人,副教授,博士,研究方向?yàn)楣こ探Y(jié)構(gòu)抗爆、抗沖擊理論(Email:ltian@tju.edu.cn)

    收稿日期:2015-08-03修回日期: 2015-10-29

    久久精品国产综合久久久| 9热在线视频观看99| 午夜视频精品福利| 亚洲久久久国产精品| 欧美成人性av电影在线观看| 一区在线观看完整版| 亚洲avbb在线观看| 一边摸一边抽搐一进一出视频| 亚洲少妇的诱惑av| 色综合站精品国产| 很黄的视频免费| 国产高清国产精品国产三级| 一区二区三区激情视频| 在线免费观看的www视频| 精品国产乱码久久久久久男人| 波多野结衣一区麻豆| av天堂久久9| 欧美不卡视频在线免费观看 | 热re99久久国产66热| 久久人人爽av亚洲精品天堂| 美女扒开内裤让男人捅视频| 女警被强在线播放| 中文字幕人妻熟女乱码| 亚洲色图 男人天堂 中文字幕| 一级片免费观看大全| 国产成人精品久久二区二区91| a在线观看视频网站| 久久久久九九精品影院| 免费人成视频x8x8入口观看| www.www免费av| 老鸭窝网址在线观看| 丝袜人妻中文字幕| 1024视频免费在线观看| 国产国语露脸激情在线看| 久久国产精品人妻蜜桃| 丝袜在线中文字幕| 久久久水蜜桃国产精品网| 两性夫妻黄色片| 欧美大码av| 变态另类成人亚洲欧美熟女 | 久久久久久久久久久久大奶| 夫妻午夜视频| 少妇 在线观看| 国产高清激情床上av| 国产成人精品久久二区二区91| 中文字幕色久视频| 美女国产高潮福利片在线看| 国产一区在线观看成人免费| 久久精品91无色码中文字幕| 淫妇啪啪啪对白视频| 精品国产亚洲在线| 国产一区二区在线av高清观看| 久9热在线精品视频| 夜夜躁狠狠躁天天躁| 动漫黄色视频在线观看| 国产三级黄色录像| 国产区一区二久久| 激情在线观看视频在线高清| 国产又色又爽无遮挡免费看| 这个男人来自地球电影免费观看| 精品第一国产精品| 老熟妇仑乱视频hdxx| 天堂√8在线中文| 男人操女人黄网站| 免费在线观看日本一区| 人人妻,人人澡人人爽秒播| 色综合站精品国产| 中出人妻视频一区二区| 又紧又爽又黄一区二区| www.999成人在线观看| 久久久国产成人免费| 看免费av毛片| 国产欧美日韩一区二区三区在线| 久久中文字幕人妻熟女| 精品国产一区二区久久| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品sss在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 在线观看一区二区三区激情| 久久天躁狠狠躁夜夜2o2o| www.熟女人妻精品国产| 在线观看www视频免费| 岛国视频午夜一区免费看| 在线观看免费日韩欧美大片| 国产精品综合久久久久久久免费 | 热re99久久国产66热| 亚洲精品国产区一区二| 国产精品二区激情视频| 极品人妻少妇av视频| 一区二区三区国产精品乱码| 精品久久久久久,| 日日干狠狠操夜夜爽| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 另类亚洲欧美激情| 视频区图区小说| 午夜日韩欧美国产| 天天影视国产精品| 黄片大片在线免费观看| 一区二区三区国产精品乱码| 午夜福利欧美成人| 在线观看免费日韩欧美大片| 精品国产国语对白av| 日韩成人在线观看一区二区三区| 最近最新免费中文字幕在线| 欧美最黄视频在线播放免费 | 91国产中文字幕| 我的亚洲天堂| 午夜影院日韩av| 免费女性裸体啪啪无遮挡网站| 精品欧美一区二区三区在线| 亚洲五月色婷婷综合| 亚洲第一青青草原| 欧美日韩视频精品一区| 久久久久久久久中文| 村上凉子中文字幕在线| 国产精品久久视频播放| 久久人人爽av亚洲精品天堂| 最好的美女福利视频网| 久久精品91无色码中文字幕| 麻豆一二三区av精品| 国产av在哪里看| 日日干狠狠操夜夜爽| 国内久久婷婷六月综合欲色啪| 亚洲熟女毛片儿| 天堂动漫精品| 亚洲一区高清亚洲精品| 男女做爰动态图高潮gif福利片 | 咕卡用的链子| 亚洲国产精品合色在线| 国产亚洲av高清不卡| www.自偷自拍.com| 国产又色又爽无遮挡免费看| 一a级毛片在线观看| 久久精品影院6| 成人免费观看视频高清| 亚洲七黄色美女视频| 日韩大码丰满熟妇| 久久久国产一区二区| 国产精品一区二区免费欧美| 精品福利观看| 亚洲国产欧美日韩在线播放| 黄色女人牲交| 视频区欧美日本亚洲| 欧美在线一区亚洲| 久久久精品国产亚洲av高清涩受| 国产精品免费视频内射| 色综合站精品国产| 色尼玛亚洲综合影院| 热re99久久精品国产66热6| 国产成人精品在线电影| 成人18禁在线播放| 午夜免费成人在线视频| 一边摸一边抽搐一进一小说| 99国产精品99久久久久| 久久亚洲精品不卡| 香蕉国产在线看| 国产成人精品久久二区二区免费| 国产欧美日韩一区二区三| 波多野结衣一区麻豆| 精品久久久精品久久久| 久久天堂一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费 | 亚洲精品久久成人aⅴ小说| 桃红色精品国产亚洲av| 久久精品亚洲精品国产色婷小说| 亚洲激情在线av| 欧美日韩一级在线毛片| 黄色成人免费大全| 亚洲狠狠婷婷综合久久图片| 国产精品 欧美亚洲| 久久婷婷成人综合色麻豆| 亚洲性夜色夜夜综合| 国产激情久久老熟女| 丰满饥渴人妻一区二区三| 91成年电影在线观看| 婷婷丁香在线五月| 好男人电影高清在线观看| 亚洲精品中文字幕一二三四区| 国产片内射在线| 久久久久久人人人人人| 最近最新中文字幕大全免费视频| 国产激情欧美一区二区| 婷婷六月久久综合丁香| 99热国产这里只有精品6| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成国产人片在线观看| av免费在线观看网站| 一二三四社区在线视频社区8| e午夜精品久久久久久久| 麻豆成人av在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲第一欧美日韩一区二区三区| 久久香蕉精品热| 亚洲性夜色夜夜综合| 丰满人妻熟妇乱又伦精品不卡| 精品电影一区二区在线| 国产精品久久视频播放| 亚洲情色 制服丝袜| 一二三四在线观看免费中文在| 欧美日本中文国产一区发布| bbb黄色大片| 欧美最黄视频在线播放免费 | aaaaa片日本免费| 在线观看66精品国产| 色综合站精品国产| 久久人人爽av亚洲精品天堂| 日本 av在线| 午夜精品在线福利| 亚洲色图av天堂| 久久久久久久午夜电影 | 啦啦啦 在线观看视频| 中出人妻视频一区二区| www.自偷自拍.com| 国产精华一区二区三区| 欧美乱色亚洲激情| 久热这里只有精品99| 国产精华一区二区三区| 日韩免费高清中文字幕av| 欧美老熟妇乱子伦牲交| 成人三级做爰电影| 侵犯人妻中文字幕一二三四区| 悠悠久久av| 国产精品二区激情视频| 桃色一区二区三区在线观看| 老司机午夜十八禁免费视频| 欧美精品亚洲一区二区| 亚洲男人天堂网一区| 最新在线观看一区二区三区| 欧美日韩黄片免| 日韩欧美免费精品| 久久热在线av| 天天影视国产精品| 人人妻人人澡人人看| 一边摸一边抽搐一进一出视频| 天堂√8在线中文| 成人国语在线视频| 久久久久久免费高清国产稀缺| 国产精品av久久久久免费| 久久久久精品国产欧美久久久| 久久人人97超碰香蕉20202| 精品福利永久在线观看| 91在线观看av| 午夜精品在线福利| 欧美乱色亚洲激情| 久久香蕉国产精品| 国产aⅴ精品一区二区三区波| 亚洲精品在线观看二区| 午夜精品在线福利| 精品少妇一区二区三区视频日本电影| 性色av乱码一区二区三区2| 精品久久久久久久久久免费视频 | 国产麻豆69| 国产精品一区二区精品视频观看| 香蕉丝袜av| 中文字幕色久视频| 50天的宝宝边吃奶边哭怎么回事| 国产熟女xx| 欧美性长视频在线观看| 亚洲国产精品合色在线| 一级黄色大片毛片| 亚洲av成人不卡在线观看播放网| 亚洲欧美日韩无卡精品| 99国产精品免费福利视频| 国产精品自产拍在线观看55亚洲| 午夜福利在线观看吧| 无人区码免费观看不卡| 在线播放国产精品三级| 99国产精品99久久久久| 午夜激情av网站| 大型黄色视频在线免费观看| 国产成年人精品一区二区 | 午夜日韩欧美国产| 日本精品一区二区三区蜜桃| 亚洲欧洲精品一区二区精品久久久| 婷婷精品国产亚洲av在线| 香蕉国产在线看| 757午夜福利合集在线观看| 国产黄a三级三级三级人| 日韩欧美在线二视频| 99久久人妻综合| 在线观看免费日韩欧美大片| 欧美精品亚洲一区二区| 99国产精品一区二区蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产综合亚洲精品| 黑人欧美特级aaaaaa片| 国产av一区二区精品久久| 99久久精品国产亚洲精品| 欧美亚洲日本最大视频资源| 成人手机av| 国产精品99久久99久久久不卡| 大陆偷拍与自拍| 我的亚洲天堂| 国产91精品成人一区二区三区| 国产单亲对白刺激| 99国产极品粉嫩在线观看| 中文欧美无线码| a级毛片黄视频| 天天添夜夜摸| 亚洲av成人一区二区三| 国产午夜精品久久久久久| 高清在线国产一区| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 91精品三级在线观看| 精品久久久久久成人av| 天天影视国产精品| 中出人妻视频一区二区| 搡老岳熟女国产| 美女 人体艺术 gogo| 丝袜在线中文字幕| 成人手机av| 亚洲五月色婷婷综合| 18禁美女被吸乳视频| 黄色毛片三级朝国网站| 久久精品影院6| 国产片内射在线| 亚洲第一av免费看| 亚洲成国产人片在线观看| 久久欧美精品欧美久久欧美| 国产精品美女特级片免费视频播放器 | 不卡一级毛片| 亚洲色图 男人天堂 中文字幕| 国产精品99久久99久久久不卡| 亚洲av第一区精品v没综合| 国产成人免费无遮挡视频| 免费在线观看完整版高清| www.熟女人妻精品国产| 国产欧美日韩一区二区三| 免费久久久久久久精品成人欧美视频| avwww免费| 国产精品亚洲av一区麻豆| 日韩免费av在线播放| 久久九九热精品免费| 九色亚洲精品在线播放| 久久亚洲真实| 亚洲色图av天堂| 久久人妻熟女aⅴ| 精品卡一卡二卡四卡免费| 国产三级黄色录像| 99国产精品99久久久久| 午夜福利在线免费观看网站| 丝袜美足系列| 淫秽高清视频在线观看| 国产成人精品久久二区二区免费| 成人免费观看视频高清| 国产1区2区3区精品| 中文字幕人妻丝袜一区二区| 看免费av毛片| 久久久水蜜桃国产精品网| 黑人欧美特级aaaaaa片| 亚洲熟妇中文字幕五十中出 | 露出奶头的视频| 国产99久久九九免费精品| 欧美激情 高清一区二区三区| 亚洲欧洲精品一区二区精品久久久| 两个人免费观看高清视频| 欧美另类亚洲清纯唯美| 精品一品国产午夜福利视频| 欧美中文综合在线视频| 婷婷六月久久综合丁香| 高清在线国产一区| 精品乱码久久久久久99久播| 亚洲欧美日韩高清在线视频| 十八禁网站免费在线| 欧美久久黑人一区二区| 国产精品二区激情视频| 一级毛片精品| 欧美日韩中文字幕国产精品一区二区三区 | 精品日产1卡2卡| 亚洲欧美日韩无卡精品| 成人亚洲精品av一区二区 | 50天的宝宝边吃奶边哭怎么回事| 老汉色av国产亚洲站长工具| 在线观看免费高清a一片| 在线天堂中文资源库| 黄色a级毛片大全视频| 欧美黄色片欧美黄色片| 精品一区二区三区av网在线观看| 亚洲久久久国产精品| 亚洲情色 制服丝袜| 亚洲第一青青草原| av国产精品久久久久影院| 一二三四社区在线视频社区8| 免费观看精品视频网站| 国产成人精品在线电影| 国产熟女xx| 日本免费a在线| 国产成+人综合+亚洲专区| 99国产综合亚洲精品| 国产伦一二天堂av在线观看| 一夜夜www| 中文字幕精品免费在线观看视频| 69av精品久久久久久| 亚洲一区中文字幕在线| 国产黄色免费在线视频| 黄色 视频免费看| 一区在线观看完整版| 悠悠久久av| 黄色片一级片一级黄色片| 老司机靠b影院| 亚洲人成电影免费在线| 国产野战对白在线观看| 精品日产1卡2卡| 国产主播在线观看一区二区| 久久人人爽av亚洲精品天堂| 日韩三级视频一区二区三区| 国产熟女xx| 亚洲国产精品一区二区三区在线| 少妇 在线观看| 日韩大码丰满熟妇| 一进一出抽搐gif免费好疼 | 久久精品国产亚洲av香蕉五月| 国产高清国产精品国产三级| 亚洲成人国产一区在线观看| 脱女人内裤的视频| 天天添夜夜摸| 亚洲午夜精品一区,二区,三区| 精品久久久久久久毛片微露脸| 黄色a级毛片大全视频| 丁香欧美五月| 国产欧美日韩一区二区精品| 手机成人av网站| 女性生殖器流出的白浆| 精品无人区乱码1区二区| 欧美成人午夜精品| 狠狠狠狠99中文字幕| 窝窝影院91人妻| 男女下面进入的视频免费午夜 | 国产欧美日韩一区二区三区在线| 一个人观看的视频www高清免费观看 | 另类亚洲欧美激情| 首页视频小说图片口味搜索| 天堂俺去俺来也www色官网| 亚洲黑人精品在线| 欧美不卡视频在线免费观看 | 亚洲五月天丁香| 黄色毛片三级朝国网站| 午夜a级毛片| 欧美精品一区二区免费开放| 亚洲一区二区三区不卡视频| 精品人妻在线不人妻| 久久香蕉激情| 亚洲午夜精品一区,二区,三区| 精品久久久精品久久久| 国产精品98久久久久久宅男小说| 精品国产乱码久久久久久男人| 国产99久久九九免费精品| 国产91精品成人一区二区三区| 午夜久久久在线观看| 乱人伦中国视频| 成人黄色视频免费在线看| 亚洲专区国产一区二区| 高清黄色对白视频在线免费看| 国产亚洲精品一区二区www| 女生性感内裤真人,穿戴方法视频| 91老司机精品| 99精品欧美一区二区三区四区| 日韩欧美在线二视频| 国产又色又爽无遮挡免费看| 可以免费在线观看a视频的电影网站| 亚洲成人免费av在线播放| 色婷婷av一区二区三区视频| 久久人人97超碰香蕉20202| 热99re8久久精品国产| 黄色 视频免费看| 日本一区二区免费在线视频| 日韩免费av在线播放| 日本一区二区免费在线视频| 成人国产一区最新在线观看| 日本一区二区免费在线视频| 成人国产一区最新在线观看| 久久久精品国产亚洲av高清涩受| 久久久久亚洲av毛片大全| 99热国产这里只有精品6| 每晚都被弄得嗷嗷叫到高潮| 午夜精品久久久久久毛片777| 欧美亚洲日本最大视频资源| 亚洲色图 男人天堂 中文字幕| 国产无遮挡羞羞视频在线观看| 又大又爽又粗| a级毛片黄视频| 51午夜福利影视在线观看| 久久精品亚洲精品国产色婷小说| 大陆偷拍与自拍| 欧美日韩国产mv在线观看视频| av有码第一页| 免费在线观看日本一区| 欧美黄色片欧美黄色片| 制服人妻中文乱码| 80岁老熟妇乱子伦牲交| 免费av中文字幕在线| 欧美黑人精品巨大| 人妻久久中文字幕网| 在线视频色国产色| av片东京热男人的天堂| 韩国精品一区二区三区| 18禁黄网站禁片午夜丰满| 18禁观看日本| 欧美日韩一级在线毛片| 久久人人97超碰香蕉20202| 三级毛片av免费| 婷婷精品国产亚洲av在线| 极品人妻少妇av视频| 精品人妻1区二区| 在线av久久热| 亚洲国产看品久久| 欧美最黄视频在线播放免费 | 人妻久久中文字幕网| 亚洲久久久国产精品| 一本大道久久a久久精品| 视频在线观看一区二区三区| 18禁观看日本| 午夜成年电影在线免费观看| 日本黄色日本黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜91福利影院| 夜夜爽天天搞| 成人国产一区最新在线观看| 成人18禁在线播放| 国产精品爽爽va在线观看网站 | aaaaa片日本免费| 新久久久久国产一级毛片| 精品国产一区二区久久| 久久久久久人人人人人| 久久性视频一级片| 亚洲专区国产一区二区| 老鸭窝网址在线观看| 香蕉久久夜色| 成年人免费黄色播放视频| 熟女少妇亚洲综合色aaa.| 午夜影院日韩av| 精品一区二区三区四区五区乱码| 亚洲一区二区三区色噜噜 | 日本wwww免费看| 国产亚洲av高清不卡| 精品国产乱码久久久久久男人| 一区二区三区国产精品乱码| 在线天堂中文资源库| 可以在线观看毛片的网站| 99热只有精品国产| 国产aⅴ精品一区二区三区波| 国产av一区二区精品久久| 国产高清国产精品国产三级| 亚洲第一av免费看| 黄色片一级片一级黄色片| 三上悠亚av全集在线观看| 一边摸一边抽搐一进一小说| 国产精品 国内视频| 中文亚洲av片在线观看爽| 久久性视频一级片| 女人被躁到高潮嗷嗷叫费观| 人人妻人人澡人人看| 无遮挡黄片免费观看| 久久草成人影院| 黄色女人牲交| 搡老乐熟女国产| 亚洲男人天堂网一区| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 在线观看舔阴道视频| 日本免费a在线| 九色亚洲精品在线播放| 国产精华一区二区三区| 老汉色∧v一级毛片| 两人在一起打扑克的视频| 天天影视国产精品| 国产成人啪精品午夜网站| 日韩欧美三级三区| av国产精品久久久久影院| 日韩视频一区二区在线观看| 欧美日韩瑟瑟在线播放| 亚洲aⅴ乱码一区二区在线播放 | 涩涩av久久男人的天堂| www国产在线视频色| av欧美777| 亚洲精品国产色婷婷电影| 在线永久观看黄色视频| 国产片内射在线| 99精品欧美一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 美女高潮喷水抽搐中文字幕| 国产aⅴ精品一区二区三区波| 午夜福利一区二区在线看| 亚洲avbb在线观看| 久久亚洲精品不卡| 成熟少妇高潮喷水视频| 亚洲片人在线观看| 亚洲欧洲精品一区二区精品久久久| 国产蜜桃级精品一区二区三区| 色尼玛亚洲综合影院| 亚洲精品一二三| 黄色丝袜av网址大全| 国产蜜桃级精品一区二区三区| 久久午夜综合久久蜜桃| 啦啦啦 在线观看视频| 亚洲午夜理论影院| 亚洲精品中文字幕一二三四区| 亚洲精品一卡2卡三卡4卡5卡| 黄片大片在线免费观看| 国产成人欧美| 精品一区二区三区四区五区乱码| 日本欧美视频一区| 久久人妻av系列| 国产一区二区在线av高清观看| 精品国产国语对白av| 亚洲一区二区三区欧美精品| 国产欧美日韩一区二区三区在线| 欧美中文日本在线观看视频| 午夜福利一区二区在线看| 脱女人内裤的视频| 黄频高清免费视频| 丰满饥渴人妻一区二区三|