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

    橋梁復(fù)合材料防車撞結(jié)構(gòu)的參數(shù)靈敏度分析及耐撞性優(yōu)化*

    2017-06-05 14:21:32吳亞鋒許明財(cái)
    關(guān)鍵詞:防撞橋墩靈敏度

    吳亞鋒 潘 晉 方 涵 許明財(cái)

    (北京汽車股份有限公司汽車研究院1) 北京 101300) (武漢理工大學(xué)交通學(xué)院2) 武漢 430063) (華中科技大學(xué)船舶與海洋工程學(xué)院3) 武漢 430074) (武漢力拓橋科防撞設(shè)施有限公司4) 武漢 430040)

    橋梁復(fù)合材料防車撞結(jié)構(gòu)的參數(shù)靈敏度分析及耐撞性優(yōu)化*

    吳亞鋒1)潘 晉2)方 涵1)許明財(cái)3,4)

    (北京汽車股份有限公司汽車研究院1)北京 101300) (武漢理工大學(xué)交通學(xué)院2)武漢 430063) (華中科技大學(xué)船舶與海洋工程學(xué)院3)武漢 430074) (武漢力拓橋科防撞設(shè)施有限公司4)武漢 430040)

    設(shè)計(jì)了以纖維增強(qiáng)復(fù)合材料和蜂窩形結(jié)構(gòu)形式為基礎(chǔ)的橋梁防車撞結(jié)構(gòu).對(duì)結(jié)構(gòu)設(shè)計(jì)變量進(jìn)行了靈敏度分析,獲得對(duì)結(jié)構(gòu)的耐撞特性影響最大的3個(gè)設(shè)計(jì)變量.以靈敏度分析結(jié)果中的3個(gè)設(shè)計(jì)變量為試驗(yàn)因子,采用均勻試驗(yàn)設(shè)計(jì)方法獲得了合理的樣本空間,并借助非線性數(shù)值仿真方法和響應(yīng)面法得到關(guān)于結(jié)構(gòu)設(shè)計(jì)變量和比吸能、碰撞力的響應(yīng)面代理模型,通過遺傳算法優(yōu)化出了比較合理的結(jié)構(gòu)尺寸.結(jié)果表明,響應(yīng)面的二次多項(xiàng)式線性回歸模型比線性回歸模型更準(zhǔn)確,優(yōu)化后的結(jié)構(gòu)尺寸不僅具有優(yōu)良的耐撞特性,可以滿足橋墩設(shè)計(jì)抗撞力的要求,而且自重減小,從而降低了防撞成本.

    GFRP;靈敏度;比吸能;響應(yīng)面法

    0 引 言

    復(fù)合材料以其比強(qiáng)度和比模量高、抗疲勞性能好、可設(shè)計(jì)性好和破損安全性好等優(yōu)良特性越來越多地被用到防撞結(jié)構(gòu)設(shè)計(jì)中[1-2].Ramakrish等[3]的研究表明,借助于合理的結(jié)構(gòu)設(shè)計(jì)和鋪層設(shè)計(jì),玻璃纖維增強(qiáng)復(fù)合材料(glassfiber reinforced plastic,GFRP)具有比一般金屬材料更為優(yōu)越的緩沖吸能特性.Shyr等[4]通過研究認(rèn)為復(fù)合材料層合板的層數(shù)是影響其吸能機(jī)理的重要因素.Hitchen等[5]通過進(jìn)一步的研究發(fā)現(xiàn)鋪層方式將影響復(fù)合材料結(jié)構(gòu)內(nèi)部的脫層損傷模式和吸能效果.Robinson等[6]的研究試驗(yàn)結(jié)果表明,沖擊損傷與沖擊能量和結(jié)構(gòu)尺寸存在復(fù)雜的非線性關(guān)系.在數(shù)值仿真方面,崔維成[7]對(duì)復(fù)合材料結(jié)構(gòu)的破壞過程進(jìn)行了數(shù)值計(jì)算,證明用計(jì)算機(jī)來模擬復(fù)合材料結(jié)構(gòu)的破壞過程是完全可行的.

    在結(jié)構(gòu)的優(yōu)化領(lǐng)域,基于響應(yīng)面法(response surface method,RSM)和遺傳算法的耐撞性優(yōu)化方法主要應(yīng)用在汽車碰撞和航空航天領(lǐng)域.如Kurtaran等[8]采用響應(yīng)面近似技術(shù),從內(nèi)能方面對(duì)圓柱殼、簡(jiǎn)化車體等結(jié)構(gòu)進(jìn)行了耐撞性約束優(yōu)化.張立新等[9]采用響應(yīng)面方法將結(jié)構(gòu)碰撞的響應(yīng)表示為設(shè)計(jì)變量的顯式函數(shù),構(gòu)造出代理模型,并利用Matlab中的優(yōu)化工具箱進(jìn)行優(yōu)化求解.Lanzi等[10]采用了響應(yīng)面法和遺傳算法相結(jié)合的方法,對(duì)復(fù)合材料薄壁元件的耐撞性能進(jìn)行了多目標(biāo)約束優(yōu)化.

    目前GFRP復(fù)合材料在橋梁防船撞領(lǐng)域已有應(yīng)用,但將其應(yīng)用于橋梁防車撞領(lǐng)域還不多見,而對(duì)橋梁復(fù)合材料防撞結(jié)構(gòu)進(jìn)行整體優(yōu)化的則更少[11-12].因此,文中以GFRP為材料,以蜂窩形為結(jié)構(gòu)形式,構(gòu)造了具有低波阻抗特性、可以增加應(yīng)力波傳播途徑的橋梁復(fù)合材料防車撞結(jié)構(gòu),同時(shí)在結(jié)構(gòu)中使用了具有良好吸能特性的夾層板結(jié)構(gòu),并對(duì)該防車撞結(jié)構(gòu)進(jìn)行了優(yōu)化.首先設(shè)計(jì)了復(fù)合材料防車撞結(jié)構(gòu)的基本尺寸,對(duì)其主要設(shè)計(jì)參數(shù)進(jìn)行了靈敏度分析;接著基于均勻試驗(yàn)設(shè)計(jì)方法構(gòu)造了重要結(jié)構(gòu)尺寸參數(shù)的樣本空間,并對(duì)結(jié)構(gòu)進(jìn)行了非線性數(shù)值仿真分析;在此基礎(chǔ)上結(jié)合響應(yīng)面法和遺傳算法對(duì)結(jié)構(gòu)進(jìn)行耐撞性優(yōu)化,以期在滿足耐撞性的前提下,減輕結(jié)構(gòu)的自重,提高結(jié)構(gòu)的安全可靠性和經(jīng)濟(jì)性.

    1 橋梁復(fù)合材料防車撞結(jié)構(gòu)設(shè)計(jì)及有限元模型

    1.1 結(jié)構(gòu)設(shè)計(jì)

    1) 結(jié)構(gòu)材料 主要考慮車橋正撞的工況,近似采用一維應(yīng)力波理論進(jìn)行分析[13],碰撞力F的計(jì)算見式(1),應(yīng)力波波速Cv的計(jì)算見式(2)

    (1)

    (2)

    式中:ρ0CvA為廣義波阻抗;σ′(ε)為線彈性材料的彈性模量;E1為低應(yīng)變率下的材料彈性模量;E2為高應(yīng)變率下的材料彈性模量;ρ為防車撞結(jié)構(gòu)所采用的材料的密度;A為車橋碰撞接觸面積;v為撞擊初始速度.

    由式(1)~(2)可知,碰撞力主要由撞擊初始速度和廣義波阻抗來決定,在碰撞速度一定的情況下減小碰撞力的最有效方法為減小廣義波阻抗.文中研究的車橋碰撞工況中撞擊接觸面積近似為不變量,因此,減小材料的密度和彈性模量是減小廣義波阻抗的最有效途徑.與普通鋼材相比,GFRP材料具有強(qiáng)度高、模量低、密度小、抗疲勞性能好、可設(shè)計(jì)性好等優(yōu)良特性,其密度約為普通鋼材的27%,主平面內(nèi)的彈性模量約普通鋼材的29%.

    同時(shí),GFRP材料的高粘彈性將在結(jié)構(gòu)中形成顯著的應(yīng)變滯后現(xiàn)象,粘性耗散能比較大,并使前一輪的反向應(yīng)力波的遲滯過程抵消后一輪正向應(yīng)力波的撞擊效應(yīng),如此反復(fù),既使碰撞過程得到柔性緩沖,又延長(zhǎng)了碰撞歷程,最終使經(jīng)防車撞結(jié)構(gòu)傳遞到橋墩的碰撞力減小.

    2) 結(jié)構(gòu)形式 車橋碰撞過程中的碰撞力和能量耗散不僅與材料特性有關(guān),還與結(jié)構(gòu)的具體形式密切相關(guān).

    增加結(jié)構(gòu)吸能最主要的是增加參與能量交換的質(zhì)量,主要有二種方法:①使用具有低模量常數(shù)的GFRP以增加應(yīng)力波傳播速度;②依靠合理的結(jié)構(gòu)形式,增加應(yīng)力波的傳播途徑.在防撞吸能結(jié)構(gòu)中具有廣泛應(yīng)用的蜂窩形可以有效地增加應(yīng)力波傳播途徑.除此之外,蜂窩形結(jié)構(gòu)柔度較大不僅能夠有效的保護(hù)橋梁和汽車,還可以增加碰撞時(shí)間歷程,有利于能量的吸收和轉(zhuǎn)換.其次,設(shè)置在蜂窩形結(jié)構(gòu)外的夾層板結(jié)構(gòu)能夠有效傳遞結(jié)構(gòu)變形,使更多的結(jié)構(gòu)材料參與到變形中,從而增強(qiáng)結(jié)構(gòu)吸能的能力,因此,文中采用蜂窩形為結(jié)構(gòu)的基本形式,見圖1.

    圖1 結(jié)構(gòu)形式圖

    1.2 有限元模型

    有限元模型包括防撞結(jié)構(gòu)和撞擊車輛二部分,見圖2,其中車輛模型取自美國(guó)國(guó)家碰撞中心(national crash analysis center,NCAC)網(wǎng)站發(fā)布的有限元模型,該模型經(jīng)過了系統(tǒng)的驗(yàn)證,具有較高的準(zhǔn)確性.

    圖2 有限元模型

    1) 單元和材料模型 GFRP蜂窩結(jié)構(gòu)和夾層板的鋼板部分均采用shell163單元進(jìn)行模擬,夾層板的夾層部分采用solid單元進(jìn)行模擬.考慮到GFRP在碰撞的過程中極易發(fā)生損傷破壞,造成強(qiáng)度和剛度的下降,文中選用可以考慮失效的復(fù)合材料模型.在文中該材料模型采用Hughes-Liu積分規(guī)則,主要輸入的材料性能見表1[14].

    表1 材料性能參數(shù)

    2) 網(wǎng)格控制 碰撞是一個(gè)動(dòng)態(tài)響應(yīng)過程,有限元模型必須嚴(yán)格控制網(wǎng)格密度,既要防止單元尺寸過大引起的主從面單元互相穿透,又要避免單元尺寸過小引起的計(jì)算時(shí)間過長(zhǎng).因此,在結(jié)構(gòu)正撞面,鋼板和蜂窩等結(jié)構(gòu)的單元尺寸為40 mm×40 mm,非迎撞面的單元為100 mm×100 mm.對(duì)于發(fā)生接觸的汽車首部和防撞結(jié)構(gòu)的正撞面劃分的單元大小相近,見圖2.

    3) 接觸、邊界和載荷 防撞結(jié)構(gòu)中各個(gè)蜂窩在碰撞過程中會(huì)發(fā)生相互接觸,設(shè)置自接觸.汽車與防撞結(jié)構(gòu)之間設(shè)置為面面接觸.根據(jù)防車撞結(jié)構(gòu)的實(shí)際約束情況,對(duì)其正撞部分內(nèi)表面節(jié)點(diǎn)的豎向和橫向位移進(jìn)行約束,同時(shí)約束其側(cè)面內(nèi)表面節(jié)點(diǎn)的豎向和縱向位移進(jìn)行約束,以模擬得到傳遞到橋墩的碰撞力.對(duì)汽車模型的轉(zhuǎn)動(dòng)和平動(dòng)不作約束,只施加初始速度V0=100 km/h.

    2 結(jié)構(gòu)設(shè)計(jì)變量靈敏度分析

    在實(shí)際的工程應(yīng)用中,為使橋梁防車撞結(jié)構(gòu)充分發(fā)揮其良好的耐撞特性,需要對(duì)其進(jìn)行合理、經(jīng)濟(jì)的優(yōu)化設(shè)計(jì),但是,設(shè)計(jì)變量過多會(huì)降低優(yōu)化設(shè)計(jì)的效率,因此,可采用靈敏度分析獲得最重要的設(shè)計(jì)變量,作為優(yōu)化設(shè)計(jì)的依據(jù).

    文中的復(fù)合材料防車撞結(jié)構(gòu)采用對(duì)稱均勻鋪層,結(jié)構(gòu)設(shè)計(jì)變量包括層合板厚度T1、夾層面板厚度T2、夾芯厚度T3、蜂窩形結(jié)構(gòu)直徑T4、單層板厚度T5和四種基本鋪層角度T6,T7,T8和T9,初始值的設(shè)計(jì)主要根據(jù)成本要求和生產(chǎn)工藝要求.根據(jù)橋梁防車撞結(jié)構(gòu)的設(shè)計(jì)依據(jù),其耐撞性能主要通過傳遞到橋墩的碰撞力、比吸能和結(jié)構(gòu)的撞深來衡量.文中對(duì)防車撞結(jié)構(gòu)的設(shè)計(jì)變量進(jìn)行了基于結(jié)構(gòu)耐撞性的靈敏度分析,分析工況包括設(shè)計(jì)變量取初始值和9組變化值時(shí)的10種工況.靈敏度計(jì)算結(jié)果見表2,其中SF,Sη和SH分別表示遞到橋墩的碰撞力、結(jié)構(gòu)比吸能和結(jié)構(gòu)撞深對(duì)設(shè)計(jì)變量的靈敏度.

    由表2的靈敏度計(jì)算結(jié)果分析可知:①影響傳遞到橋墩的碰撞力最大的2個(gè)變量分別為T3和T4.與此同時(shí),T3和T42個(gè)變量的靈敏度均為負(fù)值,說明與碰撞力呈負(fù)相關(guān)關(guān)系,因此為降低傳遞到橋墩的碰撞力最有效的辦法是減小它們的取值;②影響結(jié)構(gòu)比吸能和結(jié)構(gòu)撞深最大的3個(gè)變量均為T2,T3和T4.其中T2的靈敏度系數(shù)為負(fù)值,與比吸能和撞深均呈負(fù)相關(guān)關(guān)系,T3和T4的靈敏度系數(shù)均為正值,與比吸能和撞深呈正相關(guān)關(guān)系,因此為增加結(jié)構(gòu)的比吸能并減小結(jié)構(gòu)撞深最有效的方法是減小T2的取值,并增大T3和T4的取值.

    表2 靈敏度計(jì)算結(jié)果

    圖3 耐撞特性曲線變化圖

    圖3為靈敏度分析的10種工況下碰撞力、比吸能和撞深變化曲線.由圖3可知,結(jié)構(gòu)遞到橋墩的碰撞力、比吸能和結(jié)構(gòu)撞深之間并不是完全獨(dú)立的關(guān)系,在碰撞力較大的情況下,比吸能和撞深的值較??;在碰撞力較小的情況下,比吸能和撞深的值比較大.結(jié)構(gòu)撞深的范圍較為合理,因此結(jié)構(gòu)的耐撞性最主要通過碰撞力和比吸能來體現(xiàn).

    綜上分析可知,在進(jìn)行結(jié)構(gòu)耐撞性優(yōu)化的過程中可主要考慮設(shè)計(jì)變量T2,T3和T4,其他設(shè)計(jì)變量的影響較小可暫時(shí)忽略不計(jì).另外,在對(duì)結(jié)構(gòu)耐撞性優(yōu)化過程中,主要考慮設(shè)計(jì)變量變化對(duì)于碰撞力和比吸能的影響.

    3 車-防車撞結(jié)構(gòu)碰撞分析

    3.1 均勻試驗(yàn)設(shè)計(jì)

    以T2,T3和T4三個(gè)設(shè)計(jì)變量作為均勻試驗(yàn)設(shè)計(jì)的試驗(yàn)因子,考慮到復(fù)合材料的工藝水平和常用尺寸,T2取為2~8,T3取為80~320,T4取為8~32.其中各因子取為4水平,因此,文中設(shè)計(jì)的試驗(yàn)方案見表3.

    表3 均勻試驗(yàn)設(shè)計(jì)表

    3.2 有限元計(jì)算與分析

    在橋梁車撞的情況下,主要從以下2個(gè)方面評(píng)價(jià)防撞結(jié)構(gòu)的性能:①橋墩受到的撞擊力,文中以節(jié)點(diǎn)約束模擬橋墩的作用,因此節(jié)點(diǎn)的約束反力F即指橋墩受到的撞擊力;②防撞結(jié)構(gòu)的耐撞性,文中用比吸能(specific energy absorption,SEA)作為衡量結(jié)構(gòu)耐撞性的指標(biāo),用η來表示,指單位質(zhì)量吸收的能量,它值越大說明結(jié)構(gòu)耐撞性越好.分別對(duì)均勻試驗(yàn)設(shè)計(jì)中的各結(jié)構(gòu)參數(shù)組合進(jìn)行防車撞結(jié)構(gòu)的碰撞數(shù)值模擬,得到的計(jì)算結(jié)果見表4,其中M為結(jié)構(gòu)質(zhì)量,E為結(jié)構(gòu)在碰撞中吸收的總能量.

    4 基于響應(yīng)面法的耐撞性優(yōu)化

    在防車撞結(jié)構(gòu)的耐撞性優(yōu)化設(shè)計(jì)中,因每次調(diào)整結(jié)構(gòu)參數(shù)均要進(jìn)行相應(yīng)的碰撞計(jì)算,而碰撞計(jì)算需要的計(jì)算機(jī)時(shí)較大,嚴(yán)重影響了耐撞性優(yōu)化的分析進(jìn)度,因此,有必要采用代理模型技術(shù)對(duì)設(shè)計(jì)變量空間進(jìn)行描述.響應(yīng)面法對(duì)于接觸和碰撞這類具有較強(qiáng)非線性的問題進(jìn)行耐撞性優(yōu)化時(shí)快速且高效.因此,采用響應(yīng)面法對(duì)上節(jié)得到的結(jié)果進(jìn)行擬合,分別獲得以比吸能和碰撞力為因變量,設(shè)計(jì)變量為自變量的函數(shù)表達(dá)式,以簡(jiǎn)化耐撞性優(yōu)化的計(jì)算工作量.

    表4 有限元計(jì)算結(jié)果

    4.1 響應(yīng)面法

    以比吸能為目標(biāo)函數(shù)、碰撞力為約束函數(shù),采用線性回歸和二次多項(xiàng)式回歸方法來獲得響應(yīng)面模型.為了避免不同參數(shù)單位數(shù)量級(jí)差別太大引起的誤差,將表3中的輸入?yún)?shù)通過最大、最小值進(jìn)行無量綱化,即

    (3)

    式中:xij為輸入?yún)?shù)的量綱一的量值(xij∈[-1,1]);ξij是實(shí)際參數(shù)單位的數(shù)值.

    1) 線性回歸 采用線性回歸的前提是響應(yīng)值的分布大致成一條直線,通過最小二乘法計(jì)算得目標(biāo)函數(shù)和約束函數(shù)的回歸方程.

    η=0.534 1-0.398 1T1-0.015D+0.051 7T2

    (4)

    F=1.524 8+0.207 6T1-0.064 2D-0.101T2

    (5)

    2) 二次多項(xiàng)式回歸 為提高回歸的精度,同時(shí)采用了含交叉項(xiàng)的二次多項(xiàng)式,同樣通過最小二乘法計(jì)算可得目標(biāo)函數(shù)和約束函數(shù)的回歸方程.

    η= 0.366 4-0.407 2T1-0.015 6D+

    0.028 1T1T2+0.4079DT2

    (6)

    F= 1.723 6+0.220 4T1-0.065 1D-

    0.040 6T1T2-0.575 1DT2

    (7)

    4.2 優(yōu)化

    利用遺傳算法對(duì)上節(jié)中的二次多項(xiàng)式響應(yīng)面代理模型進(jìn)行優(yōu)化.根據(jù)AASHTO抗力設(shè)計(jì)法規(guī)范對(duì)碰撞力的規(guī)定,并考慮文中的防撞工況要求橋墩受到的碰撞力不大于1.8 MN.通過計(jì)算,優(yōu)化后的結(jié)構(gòu)尺寸參數(shù)為T2=2,T3=100,T4=32.為驗(yàn)證優(yōu)化效果,利用優(yōu)化后的結(jié)構(gòu)尺寸建立有限元模型并計(jì)算得到相關(guān)耐撞性參數(shù),將結(jié)果與表3中比吸能最大的樣本點(diǎn)2進(jìn)行比較.耐撞性參數(shù)對(duì)比結(jié)果見表5,圖4~6分別為碰撞力和比吸能時(shí)間歷程對(duì)比圖、結(jié)構(gòu)撞深變形圖和汽車應(yīng)力變形圖.

    表5 優(yōu)化結(jié)果

    圖4 碰撞力對(duì)比圖

    圖5 比吸能對(duì)比圖

    圖6 結(jié)構(gòu)位移變形圖

    通過表4、圖5~6可知,復(fù)合材料防車撞結(jié)構(gòu)在進(jìn)行耐撞性優(yōu)化后,其自重明顯減輕,碰撞力峰值明顯減小,結(jié)構(gòu)碰撞后的破壞程度減弱,且比吸能大于樣本點(diǎn)2,因此在總吸能較大、碰撞力滿足要求的情況下,優(yōu)化后的結(jié)構(gòu)能夠有效地降低橋梁的防撞成本.

    圖7為不同時(shí)刻汽車應(yīng)力變形圖,可以看出在汽車與優(yōu)化后結(jié)構(gòu)的碰撞前期,汽車的變形很小,這是因?yàn)榉雷步Y(jié)構(gòu)的柔度起到了較大的緩沖作用;而在碰撞計(jì)算結(jié)束時(shí)刻,汽車的變形也比較合理,能夠有效的保護(hù)乘員的生命安全.

    圖7 汽車應(yīng)力-變形圖

    5 結(jié) 論

    1) 橋梁復(fù)合材料防車撞結(jié)構(gòu)能夠有效增大應(yīng)力波的傳播途徑和參與碰撞的結(jié)構(gòu)質(zhì)量,因此能夠吸收較多的汽車初始動(dòng)能,使汽車的變形能較小,并減小傳遞到橋墩的碰撞力,結(jié)構(gòu)具備良好的耐撞性.

    2) 在復(fù)合材料防車撞結(jié)構(gòu)的優(yōu)化設(shè)計(jì)中,應(yīng)以對(duì)結(jié)構(gòu)耐撞性能影響較大的蜂窩尺寸、夾層結(jié)構(gòu)尺寸和GFRP的尺寸等設(shè)計(jì)變量為主要設(shè)計(jì)依據(jù).

    響應(yīng)面法能夠較準(zhǔn)確的描述本文結(jié)構(gòu)的設(shè)計(jì)空間,提高計(jì)算效率,其中響應(yīng)面的二次多項(xiàng)式回歸模型比線性回歸模型更準(zhǔn)確.最終優(yōu)化得到的結(jié)構(gòu)吸收的總能量比較大,傳遞到橋墩的撞擊力滿足規(guī)范要求,結(jié)構(gòu)的耐撞性能大大提高;優(yōu)化后的結(jié)構(gòu)自重較小,有效地降低了橋梁防撞成本,可以廣泛地應(yīng)用于橋梁防車撞的工程中.

    [1]國(guó)家發(fā)展和改革委員會(huì).中長(zhǎng)期鐵路網(wǎng)規(guī)劃方案(2008調(diào)整)[Z].北京:國(guó)務(wù)院國(guó)家發(fā)展和改革委員會(huì),2008.

    [2]尚順幫,陳豐蘭.中國(guó)高速鐵路橋梁建設(shè)的發(fā)展[J].價(jià)值工程,2013,315(32):87-88.

    [3]RAMAKRISHNA S, HAMADA H. Energy absorption characteristics of crash worthy structural composite materials[J]. Key Engineering Materials,1998(1):585-562.

    [4]SHYR T W, PAN Y H. Impact resistance and damage characteristics of composite laminates[J]. Composite Structures,2003,62(2):193-203.

    [5]HITCHEN S A, KEMP R M J. The effect of stacking sequence on impact damage in a carbon fiber/epoxy composite[J]. Composites,1995,26(3):207-214.

    [6]ROBINSON P, DAVIES G A O. Impactor mass and specimen geometry effects in low velocity impact of laminated composites[J]. International Journal of Impact Engineering,1992,12(2):189-207.

    [7]崔維成.復(fù)合材料結(jié)構(gòu)破壞過程的計(jì)算機(jī)模擬[J].復(fù)合材料學(xué)報(bào),1996,13(4):102-111.

    [8]KURTARAN, ESKANDARIAN A. MARZOUGUI D, et al. Crashworthiness design optimization using successive response surface approximations[J]. Computational Mechanics,2002,29(4-5):409-421.

    [9]張立新,隋允康,杜家政.基于響應(yīng)面方法的結(jié)構(gòu)耐撞性優(yōu)化[J].北京工業(yè)大學(xué)學(xué)報(bào),2007,33(2):129-133.

    [10]LANZI L, CASTELLETTI L M L, ANGHILERI M.Multi-objective optimization of composite absorber shape under crashworthiness requirements[J]. Composite Structures,2004,65(3):433-441.

    [11]金軒慧.復(fù)合材料防撞套箱的數(shù)值模擬研究[D].重慶:重慶交通大學(xué),2014.

    [12]姜華,耿波,張錫祥.橋墩新型防船撞裝置防撞性能研究[J].振動(dòng)與沖擊,2014,33(17):154-160.

    [13]王禮立.應(yīng)力波基礎(chǔ)[M].2版.北京:國(guó)防工業(yè)出版社,2005.

    [14]王耀光.復(fù)合材料力學(xué)與結(jié)構(gòu)設(shè)計(jì)[M].上海:華東理工大學(xué)出版社,2012.

    [15]KMIECIK M, SOARES C G. Response surface approach to the probability distribution of the strength of compressed plates[J]. Marine Structures,2002,15(2):139-156.

    [16]XU M C, TEIXEIRA A P, GUEDES S C. Polynomial based response surface probabilistic modelling of the ultimate strength of stiffened panels[J]. Journal of Marine Science and Technology,2013,18(1):115-132.

    Sensitivity Analysis and Crashworthiness Optimization of Composite Energy-absorbing Structure for Bridge-vehicle Crash

    WU Yafeng1)PAN Jin2)FANG Han2)XU Mingcai3,4)

    (BaicMotorCorporationLtd.R&DCenter,Beijing101300,China)1)(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)2)(SchoolofNavalArchitecture&OceanEngineering,HuazhongUniversityofScienceandTechnology,Wuhan430074,China)3)(WuhanLituoBridgeProtectionTechnologyCo.Ltd,Wuhan430040,China)4)

    Based on the impact dynamics theories, a fiber-reinforced composite plastics material and honeycomb-type structure for bridge-vehicle crash is designed in this paper. Sensitivity analysis method is adopted to determine the significance of the design variables on the structure crashworthiness, and three design variables are chosen as the basis of optimization design. Reasonable sample space is designed by using uniform design experimentation, and then nonlinear numerical simulation method and response surface method are adopted to obtain surrogate models about design variables, specific energy absorption and impact force. With Matlab software and genetic algorithm, more reasonable sizes of the structure are optimized. The results show that the quadratic polynomial regression model is more accurate than the linear regression model, and the optimized structure not only has excellent crashworthiness and meets the pier design requirements, but also reduces the crash costs.

    GFRP; sensitivity; specific energy absorption; response surface method

    2017-01-23

    *國(guó)家自然科學(xué)基金青年科學(xué)基金項(xiàng)目(51609192)、浙江省近海海洋工程技術(shù)重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(ZJOELAB-1602)、中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(2015MS103)資助

    U41

    10.3963/j.issn.2095-3844.2017.02.032

    吳亞鋒(1992—):男,碩士生,主要研究領(lǐng)域?yàn)檐囕v結(jié)構(gòu)分析及碰撞動(dòng)力學(xué)

    猜你喜歡
    防撞橋墩靈敏度
    橋墩加固方案分析研究
    導(dǎo)磁環(huán)對(duì)LVDT線性度和靈敏度的影響
    啄木鳥的防撞科學(xué)
    鋁合金在汽車前防撞梁輕量化中的應(yīng)用
    汽車工程師(2018年1期)2018-07-13 00:47:56
    軌道交通預(yù)制拼裝橋墩設(shè)計(jì)和施工
    上海公路(2018年1期)2018-06-26 08:37:36
    地下水非穩(wěn)定流的靈敏度分析
    霧天上高速 防撞有招數(shù)
    穿甲爆破彈引信對(duì)薄弱目標(biāo)的靈敏度分析
    空中交通防撞系統(tǒng)(TCAS Ⅱ)的排故與維護(hù)
    預(yù)制節(jié)段拼裝橋墩研究進(jìn)展
    久久久久国产网址| 久久99精品国语久久久| 亚洲国产av影院在线观看| 黑人欧美特级aaaaaa片| 国产乱来视频区| 一区福利在线观看| 男女免费视频国产| 18禁裸乳无遮挡动漫免费视频| 青青草视频在线视频观看| 亚洲国产精品成人久久小说| 久久人人97超碰香蕉20202| 欧美黄色片欧美黄色片| 日韩av不卡免费在线播放| 色吧在线观看| 性色avwww在线观看| 亚洲精品国产av蜜桃| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品人妻久久久影院| 亚洲国产欧美网| 日韩中文字幕欧美一区二区 | 精品亚洲成a人片在线观看| 欧美人与善性xxx| 这个男人来自地球电影免费观看 | 国产老妇伦熟女老妇高清| 日日爽夜夜爽网站| 亚洲图色成人| 欧美日韩一区二区视频在线观看视频在线| 欧美精品国产亚洲| 一本久久精品| 日本91视频免费播放| 久久精品人人爽人人爽视色| 极品人妻少妇av视频| 日日爽夜夜爽网站| 婷婷成人精品国产| 在线天堂中文资源库| 日本av免费视频播放| 精品第一国产精品| 晚上一个人看的免费电影| 国产精品二区激情视频| 久久精品国产自在天天线| 一级毛片电影观看| 国产精品一二三区在线看| 热99久久久久精品小说推荐| 欧美精品一区二区大全| 国产片内射在线| 亚洲精品国产av成人精品| 日日啪夜夜爽| 最近最新中文字幕大全免费视频 | 一区二区av电影网| 亚洲色图综合在线观看| 国产精品蜜桃在线观看| 欧美日本中文国产一区发布| 少妇人妻 视频| 91精品伊人久久大香线蕉| 亚洲精品国产av蜜桃| 国产精品.久久久| 国产极品粉嫩免费观看在线| av电影中文网址| 精品亚洲成a人片在线观看| 亚洲成人av在线免费| 国产一区二区在线观看av| av国产久精品久网站免费入址| 黑人巨大精品欧美一区二区蜜桃| 免费观看av网站的网址| 欧美精品国产亚洲| 国产日韩欧美视频二区| 亚洲国产日韩一区二区| 免费在线观看黄色视频的| 国精品久久久久久国模美| 亚洲国产欧美日韩在线播放| 色网站视频免费| 少妇人妻久久综合中文| 中文字幕精品免费在线观看视频| 精品人妻在线不人妻| 老鸭窝网址在线观看| 男女边吃奶边做爰视频| 超碰97精品在线观看| 亚洲,欧美,日韩| 亚洲色图 男人天堂 中文字幕| 色视频在线一区二区三区| 国产成人欧美| 国产精品av久久久久免费| 成人国产av品久久久| 久久久久久久大尺度免费视频| 日本黄色日本黄色录像| 亚洲国产av影院在线观看| 天天躁日日躁夜夜躁夜夜| 男女午夜视频在线观看| 永久网站在线| 丰满迷人的少妇在线观看| 久久久精品区二区三区| 免费高清在线观看日韩| 高清av免费在线| 久久久久国产精品人妻一区二区| 精品少妇黑人巨大在线播放| 极品人妻少妇av视频| 亚洲欧美清纯卡通| 久久久欧美国产精品| 大香蕉久久成人网| 欧美日韩国产mv在线观看视频| 日本91视频免费播放| 国产免费福利视频在线观看| 多毛熟女@视频| 建设人人有责人人尽责人人享有的| 尾随美女入室| 亚洲av男天堂| 麻豆乱淫一区二区| 黄色一级大片看看| 在线观看三级黄色| 久久精品久久久久久久性| 国产探花极品一区二区| 最近最新中文字幕免费大全7| 日本欧美国产在线视频| 国产又色又爽无遮挡免| 成人国产麻豆网| 午夜福利影视在线免费观看| 亚洲欧美成人精品一区二区| 久久久精品国产亚洲av高清涩受| 黄色视频在线播放观看不卡| 熟女av电影| 欧美日韩一级在线毛片| 国产老妇伦熟女老妇高清| 另类亚洲欧美激情| 少妇人妻 视频| 国产精品久久久久成人av| 男女无遮挡免费网站观看| 视频在线观看一区二区三区| 精品少妇内射三级| 国产日韩欧美在线精品| 青草久久国产| 高清不卡的av网站| 日韩电影二区| 婷婷色综合大香蕉| 亚洲一区中文字幕在线| 久久鲁丝午夜福利片| 久久久精品国产亚洲av高清涩受| 五月天丁香电影| 1024香蕉在线观看| 99热网站在线观看| 成年人午夜在线观看视频| 国产又色又爽无遮挡免| av片东京热男人的天堂| 最近2019中文字幕mv第一页| 久久久国产精品麻豆| 一本—道久久a久久精品蜜桃钙片| a级毛片黄视频| 1024香蕉在线观看| 一本—道久久a久久精品蜜桃钙片| 丝袜脚勾引网站| 欧美+日韩+精品| 久久久精品国产亚洲av高清涩受| 欧美日韩av久久| 亚洲男人天堂网一区| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 日韩熟女老妇一区二区性免费视频| 18禁裸乳无遮挡动漫免费视频| 一级a爱视频在线免费观看| 亚洲久久久国产精品| 国产一区二区在线观看av| 亚洲国产精品成人久久小说| 99re6热这里在线精品视频| 日本wwww免费看| www日本在线高清视频| 精品亚洲乱码少妇综合久久| 国产黄色免费在线视频| 一本久久精品| 十八禁网站网址无遮挡| 久久久久久伊人网av| a 毛片基地| 精品人妻熟女毛片av久久网站| 久久这里只有精品19| 亚洲欧美中文字幕日韩二区| 国产成人精品在线电影| 成人18禁高潮啪啪吃奶动态图| 久久精品国产亚洲av天美| 男女下面插进去视频免费观看| 青青草视频在线视频观看| 伊人久久国产一区二区| 久久免费观看电影| 午夜91福利影院| 两个人免费观看高清视频| 一区在线观看完整版| av电影中文网址| 蜜桃在线观看..| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 国产av精品麻豆| 91久久精品国产一区二区三区| 纯流量卡能插随身wifi吗| 天堂中文最新版在线下载| 久久热在线av| 一区福利在线观看| 少妇人妻 视频| 国产精品99久久99久久久不卡 | 两性夫妻黄色片| 咕卡用的链子| 亚洲精品一二三| 精品一区二区免费观看| 少妇人妻精品综合一区二区| 国产高清不卡午夜福利| 亚洲精品国产色婷婷电影| www.自偷自拍.com| 欧美精品一区二区免费开放| 我要看黄色一级片免费的| 人妻一区二区av| 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| 成人毛片60女人毛片免费| 欧美少妇被猛烈插入视频| 自线自在国产av| 久久免费观看电影| av不卡在线播放| 久久精品国产亚洲av涩爱| 亚洲精品乱久久久久久| 天堂俺去俺来也www色官网| 伦理电影大哥的女人| 免费久久久久久久精品成人欧美视频| 精品亚洲乱码少妇综合久久| av免费观看日本| 国产精品嫩草影院av在线观看| 日本黄色日本黄色录像| 大片电影免费在线观看免费| 极品少妇高潮喷水抽搐| 最黄视频免费看| 2021少妇久久久久久久久久久| 99re6热这里在线精品视频| 亚洲国产欧美在线一区| 乱人伦中国视频| 亚洲激情五月婷婷啪啪| 免费日韩欧美在线观看| 五月开心婷婷网| 欧美国产精品一级二级三级| 伊人久久国产一区二区| 99国产综合亚洲精品| 伦理电影免费视频| 精品亚洲乱码少妇综合久久| 免费在线观看黄色视频的| 日韩 亚洲 欧美在线| 老司机影院毛片| www.熟女人妻精品国产| 99久久人妻综合| 一区二区三区精品91| av不卡在线播放| www.av在线官网国产| 成年女人毛片免费观看观看9 | 黄片播放在线免费| 久久久精品免费免费高清| 少妇人妻久久综合中文| 丝袜在线中文字幕| 亚洲av欧美aⅴ国产| 丝袜美足系列| 亚洲国产精品999| 亚洲一级一片aⅴ在线观看| 成人二区视频| 亚洲精品日韩在线中文字幕| 9热在线视频观看99| 亚洲人成77777在线视频| 黄片播放在线免费| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看| 午夜精品国产一区二区电影| 街头女战士在线观看网站| 99久国产av精品国产电影| 大香蕉久久网| 成人18禁高潮啪啪吃奶动态图| 亚洲第一av免费看| 国产97色在线日韩免费| av线在线观看网站| 亚洲欧洲国产日韩| 在线亚洲精品国产二区图片欧美| 欧美精品一区二区免费开放| 黄网站色视频无遮挡免费观看| 亚洲av日韩在线播放| 国产1区2区3区精品| 国产精品国产av在线观看| 日韩大片免费观看网站| 亚洲伊人色综图| 成人黄色视频免费在线看| 97人妻天天添夜夜摸| 夫妻性生交免费视频一级片| 91精品国产国语对白视频| 午夜免费男女啪啪视频观看| 午夜久久久在线观看| 精品亚洲乱码少妇综合久久| 搡女人真爽免费视频火全软件| 一区二区三区激情视频| freevideosex欧美| 久久精品夜色国产| 老鸭窝网址在线观看| 婷婷色综合www| videosex国产| 日日爽夜夜爽网站| 成年人免费黄色播放视频| 久久久精品94久久精品| 久久婷婷青草| 日韩av免费高清视频| 女人高潮潮喷娇喘18禁视频| 日本wwww免费看| 飞空精品影院首页| 亚洲熟女精品中文字幕| 免费看av在线观看网站| 免费人妻精品一区二区三区视频| 一边摸一边做爽爽视频免费| 最新的欧美精品一区二区| av女优亚洲男人天堂| 熟女电影av网| 91精品伊人久久大香线蕉| 久久精品国产综合久久久| 女人高潮潮喷娇喘18禁视频| 久久精品国产自在天天线| 性色av一级| 日韩av在线免费看完整版不卡| 老司机亚洲免费影院| 夫妻性生交免费视频一级片| videossex国产| www.av在线官网国产| 国产综合精华液| 亚洲人成网站在线观看播放| 国产亚洲最大av| 亚洲国产看品久久| 自拍欧美九色日韩亚洲蝌蚪91| 日日啪夜夜爽| 七月丁香在线播放| 在线观看人妻少妇| 免费不卡的大黄色大毛片视频在线观看| 丝瓜视频免费看黄片| 亚洲天堂av无毛| 一区二区三区激情视频| 欧美精品一区二区免费开放| 在线观看人妻少妇| 亚洲精品乱久久久久久| 国产成人精品一,二区| 99久久中文字幕三级久久日本| 欧美成人午夜免费资源| 亚洲国产av影院在线观看| 狠狠精品人妻久久久久久综合| 两性夫妻黄色片| 国产日韩欧美视频二区| 久久99蜜桃精品久久| 亚洲精品成人av观看孕妇| 美女大奶头黄色视频| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 日韩 亚洲 欧美在线| 美女大奶头黄色视频| 国产老妇伦熟女老妇高清| 成年美女黄网站色视频大全免费| 国产精品嫩草影院av在线观看| 中文字幕人妻丝袜制服| 一区二区三区乱码不卡18| 亚洲成色77777| 久久 成人 亚洲| 国产一区二区三区综合在线观看| 亚洲 欧美一区二区三区| 精品少妇黑人巨大在线播放| 中文字幕av电影在线播放| 不卡视频在线观看欧美| 国产免费视频播放在线视频| 国产男女超爽视频在线观看| 一本久久精品| 最近最新中文字幕大全免费视频 | 精品国产国语对白av| 国产精品一区二区在线观看99| 国产成人免费无遮挡视频| 日本黄色日本黄色录像| 午夜福利视频精品| 大陆偷拍与自拍| 性色avwww在线观看| 国产精品av久久久久免费| 国语对白做爰xxxⅹ性视频网站| 老司机影院毛片| av国产久精品久网站免费入址| 夫妻性生交免费视频一级片| 又黄又粗又硬又大视频| 国产成人精品久久二区二区91 | 婷婷色av中文字幕| 观看av在线不卡| 飞空精品影院首页| 国产精品秋霞免费鲁丝片| 伊人久久大香线蕉亚洲五| 免费播放大片免费观看视频在线观看| 亚洲国产精品国产精品| 亚洲一级一片aⅴ在线观看| 国产精品麻豆人妻色哟哟久久| 可以免费在线观看a视频的电影网站 | 成年av动漫网址| 日韩精品有码人妻一区| 极品少妇高潮喷水抽搐| 国产男人的电影天堂91| 少妇的逼水好多| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 日本av手机在线免费观看| 日本黄色日本黄色录像| 久久午夜福利片| 人体艺术视频欧美日本| 亚洲国产色片| 丝袜人妻中文字幕| 香蕉精品网在线| 久久精品国产自在天天线| 999精品在线视频| 美女主播在线视频| 亚洲伊人色综图| 亚洲一级一片aⅴ在线观看| 国产精品久久久av美女十八| 日本猛色少妇xxxxx猛交久久| 午夜免费观看性视频| 精品一区在线观看国产| 99国产精品免费福利视频| 久久97久久精品| 久久女婷五月综合色啪小说| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品av麻豆狂野| 伦理电影大哥的女人| videossex国产| 亚洲第一区二区三区不卡| 久久久久网色| 亚洲精品,欧美精品| 欧美成人精品欧美一级黄| 日韩精品有码人妻一区| 国产精品麻豆人妻色哟哟久久| 2018国产大陆天天弄谢| 亚洲视频免费观看视频| 欧美中文综合在线视频| 亚洲第一av免费看| www.精华液| 亚洲视频免费观看视频| 欧美国产精品一级二级三级| 国产色婷婷99| 97在线视频观看| 亚洲国产精品999| 国精品久久久久久国模美| 大陆偷拍与自拍| 亚洲精品日本国产第一区| 午夜日本视频在线| 欧美 日韩 精品 国产| 久久久精品国产亚洲av高清涩受| 免费少妇av软件| 精品一区在线观看国产| www.av在线官网国产| 午夜老司机福利剧场| 久久99一区二区三区| 亚洲五月色婷婷综合| 一区在线观看完整版| 久久精品国产鲁丝片午夜精品| 中文字幕色久视频| 国产成人午夜福利电影在线观看| 啦啦啦在线免费观看视频4| 久久久久久久久久人人人人人人| 老司机亚洲免费影院| 一区二区三区乱码不卡18| 一级毛片黄色毛片免费观看视频| 黄色视频在线播放观看不卡| 男人操女人黄网站| 新久久久久国产一级毛片| 高清av免费在线| 国产伦理片在线播放av一区| 国产免费视频播放在线视频| 国产成人精品久久久久久| 97人妻天天添夜夜摸| 香蕉国产在线看| 亚洲图色成人| www.精华液| 99国产精品免费福利视频| 亚洲国产欧美网| 男女下面插进去视频免费观看| 国产精品国产三级国产专区5o| 中文字幕人妻熟女乱码| 国产精品麻豆人妻色哟哟久久| 国产片内射在线| 久久鲁丝午夜福利片| 一本久久精品| 丝袜人妻中文字幕| 丝瓜视频免费看黄片| 欧美精品国产亚洲| 99re6热这里在线精品视频| 国产免费一区二区三区四区乱码| 一级片免费观看大全| 黄色 视频免费看| 99久久中文字幕三级久久日本| 精品人妻在线不人妻| 看免费av毛片| 欧美xxⅹ黑人| 久久精品久久久久久久性| 久久这里只有精品19| 日韩精品免费视频一区二区三区| 日本爱情动作片www.在线观看| 叶爱在线成人免费视频播放| 亚洲欧美日韩另类电影网站| 日韩不卡一区二区三区视频在线| 精品一品国产午夜福利视频| 一区二区三区四区激情视频| 国产福利在线免费观看视频| 老汉色∧v一级毛片| 十分钟在线观看高清视频www| 考比视频在线观看| 在线观看免费高清a一片| 国产成人精品久久久久久| 丰满饥渴人妻一区二区三| 久久这里只有精品19| 午夜福利乱码中文字幕| 日日爽夜夜爽网站| 午夜福利网站1000一区二区三区| 亚洲成人一二三区av| 亚洲av.av天堂| 18+在线观看网站| 热re99久久国产66热| 亚洲av欧美aⅴ国产| 飞空精品影院首页| 中文精品一卡2卡3卡4更新| 亚洲国产精品999| 性色avwww在线观看| 婷婷色麻豆天堂久久| av线在线观看网站| 男女国产视频网站| 波野结衣二区三区在线| 97精品久久久久久久久久精品| 日韩不卡一区二区三区视频在线| 久久精品aⅴ一区二区三区四区 | a级片在线免费高清观看视频| 在线观看免费高清a一片| 亚洲av欧美aⅴ国产| 国产欧美日韩一区二区三区在线| 亚洲内射少妇av| 亚洲精品美女久久av网站| 女性被躁到高潮视频| 免费观看a级毛片全部| 国产日韩欧美在线精品| 久久久久久久精品精品| 久久人妻熟女aⅴ| 激情视频va一区二区三区| 国产97色在线日韩免费| 国产成人一区二区在线| 色哟哟·www| 九草在线视频观看| 国产 精品1| 午夜av观看不卡| 国产日韩欧美在线精品| 91在线精品国自产拍蜜月| 亚洲欧美精品自产自拍| 婷婷成人精品国产| 又粗又硬又长又爽又黄的视频| 亚洲精品美女久久av网站| 亚洲一区中文字幕在线| 免费在线观看黄色视频的| 日韩一区二区视频免费看| 国产1区2区3区精品| 精品人妻偷拍中文字幕| 亚洲伊人色综图| 在线观看国产h片| 国产av精品麻豆| 国产亚洲一区二区精品| 亚洲成国产人片在线观看| 亚洲国产精品999| 99久久精品国产国产毛片| 街头女战士在线观看网站| 天天影视国产精品| 亚洲av电影在线观看一区二区三区| 国产熟女欧美一区二区| 天堂俺去俺来也www色官网| 免费不卡的大黄色大毛片视频在线观看| 欧美+日韩+精品| 女性被躁到高潮视频| 久久久久精品久久久久真实原创| 国产午夜精品一二区理论片| 9191精品国产免费久久| 只有这里有精品99| 国产成人精品一,二区| 一区二区三区四区激情视频| 美女大奶头黄色视频| 国产在线一区二区三区精| 成人18禁高潮啪啪吃奶动态图| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久久久免| 色哟哟·www| 黑丝袜美女国产一区| 久久久久久免费高清国产稀缺| av又黄又爽大尺度在线免费看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久精品国产亚洲av涩爱| 亚洲国产色片| 精品一区二区免费观看| 亚洲,欧美,日韩| av女优亚洲男人天堂| 精品久久久久久电影网| 国产成人一区二区在线| 亚洲一区中文字幕在线| 亚洲国产av影院在线观看| 亚洲国产精品999| 五月天丁香电影| 一区二区av电影网| 国产毛片在线视频| 国产精品久久久久久精品电影小说| 亚洲国产av新网站| 久久这里有精品视频免费| av国产久精品久网站免费入址| 成人手机av| 亚洲国产看品久久| 国产免费一区二区三区四区乱码| 极品少妇高潮喷水抽搐| 亚洲av.av天堂| 看免费av毛片| 亚洲三区欧美一区| 大香蕉久久成人网| 侵犯人妻中文字幕一二三四区| 男女啪啪激烈高潮av片| 最近2019中文字幕mv第一页| 国产精品偷伦视频观看了| 妹子高潮喷水视频| 免费观看a级毛片全部| 蜜桃国产av成人99| 欧美日韩视频精品一区| 999久久久国产精品视频| 日韩不卡一区二区三区视频在线|