• <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)展
    日本a在线网址| 国产成人影院久久av| 久久精品影院6| 欧美日本亚洲视频在线播放| 欧美性长视频在线观看| 岛国在线观看网站| 一二三四在线观看免费中文在| 老司机在亚洲福利影院| 午夜两性在线视频| 国产无遮挡羞羞视频在线观看| 纯流量卡能插随身wifi吗| 精品人妻1区二区| 精品国产一区二区三区四区第35| avwww免费| 国产精品免费一区二区三区在线| 国产黄a三级三级三级人| 又大又爽又粗| 久久国产精品人妻蜜桃| 男人的好看免费观看在线视频 | 中文字幕av电影在线播放| 一二三四在线观看免费中文在| 精品一区二区三卡| 免费搜索国产男女视频| 亚洲精品国产色婷婷电影| 欧美激情久久久久久爽电影 | 美女高潮到喷水免费观看| 久久天躁狠狠躁夜夜2o2o| 可以免费在线观看a视频的电影网站| 成人精品一区二区免费| 亚洲精品一区av在线观看| 久久香蕉国产精品| 99精品久久久久人妻精品| 久久热在线av| 少妇粗大呻吟视频| 亚洲第一青青草原| 欧美日韩视频精品一区| 91大片在线观看| 亚洲av五月六月丁香网| 在线观看一区二区三区激情| 99国产精品一区二区蜜桃av| 免费av中文字幕在线| 欧美日韩瑟瑟在线播放| 天天影视国产精品| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 国产乱人伦免费视频| 啦啦啦免费观看视频1| 中文字幕另类日韩欧美亚洲嫩草| av网站在线播放免费| 高清欧美精品videossex| 久久人人97超碰香蕉20202| 国产一区二区三区综合在线观看| 18禁黄网站禁片午夜丰满| 中文字幕最新亚洲高清| 国产精品免费视频内射| 久久精品aⅴ一区二区三区四区| 丝袜美腿诱惑在线| 老司机福利观看| 一本综合久久免费| av中文乱码字幕在线| 国产精品九九99| 香蕉国产在线看| 搡老乐熟女国产| 午夜激情av网站| 亚洲第一av免费看| 免费看十八禁软件| 亚洲精品一区av在线观看| 天堂俺去俺来也www色官网| 亚洲片人在线观看| 91在线观看av| 大码成人一级视频| 欧美激情高清一区二区三区| 好看av亚洲va欧美ⅴa在| 婷婷丁香在线五月| a在线观看视频网站| 午夜a级毛片| 亚洲精品成人av观看孕妇| 日日摸夜夜添夜夜添小说| 精品久久蜜臀av无| 国产极品粉嫩免费观看在线| 国产精品免费一区二区三区在线| 韩国精品一区二区三区| 亚洲中文字幕日韩| 中文字幕最新亚洲高清| 亚洲精品美女久久av网站| 大型黄色视频在线免费观看| 美女国产高潮福利片在线看| 亚洲精品久久午夜乱码| 老司机午夜福利在线观看视频| 电影成人av| 黄色视频不卡| 一区二区日韩欧美中文字幕| 成人av一区二区三区在线看| 97人妻天天添夜夜摸| 久久久久久久久中文| 成人18禁高潮啪啪吃奶动态图| 法律面前人人平等表现在哪些方面| 一级作爱视频免费观看| 国产精品久久视频播放| 亚洲国产精品一区二区三区在线| 国产欧美日韩精品亚洲av| 国产亚洲精品第一综合不卡| 18禁观看日本| 亚洲专区中文字幕在线| 丰满饥渴人妻一区二区三| 午夜两性在线视频| 国产精品1区2区在线观看.| 美女国产高潮福利片在线看| 亚洲色图av天堂| 日日摸夜夜添夜夜添小说| 少妇的丰满在线观看| 99久久精品国产亚洲精品| 久久午夜亚洲精品久久| 99riav亚洲国产免费| 桃红色精品国产亚洲av| 大陆偷拍与自拍| 亚洲精品美女久久久久99蜜臀| 日韩成人在线观看一区二区三区| 91成人精品电影| 亚洲一区二区三区色噜噜 | aaaaa片日本免费| av超薄肉色丝袜交足视频| 一区二区日韩欧美中文字幕| a级毛片黄视频| 国产成人av激情在线播放| 亚洲一区高清亚洲精品| 嫁个100分男人电影在线观看| 亚洲九九香蕉| 亚洲第一青青草原| 一进一出好大好爽视频| 成人国产一区最新在线观看| 久久婷婷成人综合色麻豆| 一级黄色大片毛片| 99国产精品一区二区三区| 老司机深夜福利视频在线观看| 超碰成人久久| 搡老乐熟女国产| 国产欧美日韩综合在线一区二区| 中亚洲国语对白在线视频| av在线天堂中文字幕 | 老鸭窝网址在线观看| 亚洲av第一区精品v没综合| 亚洲熟妇中文字幕五十中出 | 精品久久久久久久毛片微露脸| 国产成年人精品一区二区 | 久久精品亚洲av国产电影网| av超薄肉色丝袜交足视频| 亚洲国产欧美一区二区综合| 99精国产麻豆久久婷婷| 久久久精品欧美日韩精品| 国产又爽黄色视频| 亚洲美女黄片视频| 亚洲精品一区av在线观看| 日韩免费av在线播放| 麻豆久久精品国产亚洲av | 狠狠狠狠99中文字幕| 日韩人妻精品一区2区三区| 免费女性裸体啪啪无遮挡网站| 久久久久久大精品| 久久精品影院6| 久久国产精品男人的天堂亚洲| 亚洲国产毛片av蜜桃av| av天堂在线播放| 最好的美女福利视频网| 老熟妇仑乱视频hdxx| 波多野结衣高清无吗| 黑人巨大精品欧美一区二区蜜桃| 亚洲 国产 在线| 无遮挡黄片免费观看| 色在线成人网| av天堂在线播放| av天堂中文字幕网| av国产免费在线观看| 老司机福利观看| 久久久久久久久久成人| 精品欧美国产一区二区三| 两个人视频免费观看高清| 欧美日韩乱码在线| 欧美成人性av电影在线观看| 成年女人永久免费观看视频| 长腿黑丝高跟| 久久亚洲精品不卡| eeuss影院久久| 国产不卡一卡二| 一区福利在线观看| 哪里可以看免费的av片| 久久久久久国产a免费观看| 精品一区二区免费观看| 有码 亚洲区| 久久久久国内视频| 成人鲁丝片一二三区免费| 欧美一级a爱片免费观看看| 夜夜躁狠狠躁天天躁| 人人妻人人看人人澡| 人妻制服诱惑在线中文字幕| 一级黄色大片毛片| 91午夜精品亚洲一区二区三区 | 黄色一级大片看看| 麻豆成人av在线观看| 人妻夜夜爽99麻豆av| 久久国产乱子免费精品| 十八禁网站免费在线| 51午夜福利影视在线观看| 精品欧美国产一区二区三| 老司机午夜福利在线观看视频| 级片在线观看| 天堂网av新在线| 国产精品国产高清国产av| 色哟哟哟哟哟哟| www.色视频.com| 一夜夜www| 久久九九热精品免费| 女人被狂操c到高潮| 欧美国产日韩亚洲一区| 一本精品99久久精品77| 亚洲成av人片在线播放无| 久久香蕉精品热| 日本与韩国留学比较| 亚洲精品色激情综合| 91久久精品国产一区二区成人| 欧美色欧美亚洲另类二区| 人妻丰满熟妇av一区二区三区| 内射极品少妇av片p| 一个人看视频在线观看www免费| 麻豆久久精品国产亚洲av| 免费电影在线观看免费观看| 波多野结衣高清作品| 91在线精品国自产拍蜜月| 久久精品影院6| 欧美在线一区亚洲| 99在线视频只有这里精品首页| 国产精品亚洲av一区麻豆| 一级黄片播放器| 成人美女网站在线观看视频| 亚洲色图av天堂| 国产三级在线视频| 亚洲精品粉嫩美女一区| 日本免费a在线| 亚洲精华国产精华精| 成人性生交大片免费视频hd| 色av中文字幕| 高潮久久久久久久久久久不卡| 国产精品人妻久久久久久| 老熟妇乱子伦视频在线观看| 悠悠久久av| www日本黄色视频网| 国产一区二区在线av高清观看| 深爱激情五月婷婷| 日韩亚洲欧美综合| 国产单亲对白刺激| 又爽又黄a免费视频| 看黄色毛片网站| 免费看光身美女| 999久久久精品免费观看国产| 午夜福利18| 99热只有精品国产| 99国产综合亚洲精品| 女人十人毛片免费观看3o分钟| 禁无遮挡网站| 成人美女网站在线观看视频| 黄色一级大片看看| 两个人的视频大全免费| 亚洲国产精品999在线| 我的女老师完整版在线观看| 国产黄a三级三级三级人| 偷拍熟女少妇极品色| 99riav亚洲国产免费| 午夜激情福利司机影院| 国产成人啪精品午夜网站| 成人美女网站在线观看视频| 怎么达到女性高潮| 特级一级黄色大片| 全区人妻精品视频| 人人妻人人澡欧美一区二区| 亚洲av免费在线观看| 久久久久九九精品影院| 精品人妻一区二区三区麻豆 | 啦啦啦观看免费观看视频高清| 成年版毛片免费区| 久久人人精品亚洲av| 欧美另类亚洲清纯唯美| 国产精品人妻久久久久久| 国产精品影院久久| 国产成人福利小说| 国产高清有码在线观看视频| 1024手机看黄色片| 免费人成在线观看视频色| 日韩欧美 国产精品| 床上黄色一级片| 午夜福利在线观看吧| 日韩有码中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美国产在线观看| 熟妇人妻久久中文字幕3abv| 午夜影院日韩av| 一夜夜www| 亚洲人成网站在线播放欧美日韩| 91午夜精品亚洲一区二区三区 | 国产精品野战在线观看| 在线观看免费视频日本深夜| 久久精品国产清高在天天线| 在线观看美女被高潮喷水网站 | 亚洲国产精品久久男人天堂| 中文亚洲av片在线观看爽| 国产91精品成人一区二区三区| 亚洲国产精品成人综合色| 国产日本99.免费观看| 99精品在免费线老司机午夜| 99久久久亚洲精品蜜臀av| 久久精品夜夜夜夜夜久久蜜豆| 制服丝袜大香蕉在线| 成人鲁丝片一二三区免费| 日韩精品中文字幕看吧| 日韩欧美一区二区三区在线观看| 午夜免费激情av| 国产在线精品亚洲第一网站| 国产乱人伦免费视频| 男人舔奶头视频| 国产精品一区二区三区四区久久| 91久久精品国产一区二区成人| 国产成人啪精品午夜网站| 黄色一级大片看看| avwww免费| 夜夜爽天天搞| 日韩 亚洲 欧美在线| 国产 一区 欧美 日韩| 免费无遮挡裸体视频| 久久国产乱子伦精品免费另类| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 精品人妻一区二区三区麻豆 | 日韩成人在线观看一区二区三区| 天堂动漫精品| 中文字幕精品亚洲无线码一区| 18+在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品久久久久久一区二区三区 | 成人美女网站在线观看视频| 亚洲专区国产一区二区| 久久精品国产自在天天线| 黄色视频,在线免费观看| 国产人妻一区二区三区在| 尤物成人国产欧美一区二区三区| 一个人看的www免费观看视频| 在线免费观看的www视频| 禁无遮挡网站| 中文字幕高清在线视频| 久久精品人妻少妇| 久久精品国产自在天天线| 国产精品女同一区二区软件 | 男人和女人高潮做爰伦理| 成人欧美大片| 精品久久久久久,| 毛片一级片免费看久久久久 | 国产国拍精品亚洲av在线观看| 欧美日韩福利视频一区二区| 赤兔流量卡办理| 日韩欧美精品v在线| 亚洲精品影视一区二区三区av| 麻豆av噜噜一区二区三区| 久久久久久久亚洲中文字幕 | 亚洲性夜色夜夜综合| 亚洲五月婷婷丁香| 成人美女网站在线观看视频| 国产精品综合久久久久久久免费| 国产精品爽爽va在线观看网站| 免费看日本二区| 窝窝影院91人妻| 午夜两性在线视频| 国产精品一及| 亚洲精品久久国产高清桃花| 如何舔出高潮| 老熟妇乱子伦视频在线观看| 久久亚洲精品不卡| a在线观看视频网站| 校园春色视频在线观看| 国产毛片a区久久久久| 搡女人真爽免费视频火全软件 | 日本一二三区视频观看| 99热这里只有是精品50| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 亚洲一区高清亚洲精品| 成人高潮视频无遮挡免费网站| 国产伦精品一区二区三区视频9| 午夜精品久久久久久毛片777| 色综合婷婷激情| 99久久精品热视频| 成人性生交大片免费视频hd| 看片在线看免费视频| 久久久久久九九精品二区国产| 哪里可以看免费的av片| 国产亚洲精品久久久com| 一本综合久久免费| 成年女人永久免费观看视频| 熟女人妻精品中文字幕| 99国产精品一区二区三区| 国产欧美日韩一区二区三| 一夜夜www| 亚洲人成网站在线播| 精品不卡国产一区二区三区| 精品国产三级普通话版| 亚洲国产精品成人综合色| 国产大屁股一区二区在线视频| 每晚都被弄得嗷嗷叫到高潮| 我要看日韩黄色一级片| 亚洲avbb在线观看| 亚洲国产高清在线一区二区三| 国产欧美日韩精品亚洲av| 欧美色欧美亚洲另类二区| 亚洲av五月六月丁香网| 欧美一区二区国产精品久久精品| 成熟少妇高潮喷水视频| 欧美不卡视频在线免费观看| avwww免费| 国产精华一区二区三区| 怎么达到女性高潮| 国产三级中文精品| 国产黄色小视频在线观看| 色综合欧美亚洲国产小说| 欧美最新免费一区二区三区 | 我的女老师完整版在线观看| 午夜免费激情av| 免费黄网站久久成人精品 | 嫩草影院新地址| 老女人水多毛片| 女同久久另类99精品国产91| 亚洲18禁久久av| 精品人妻熟女av久视频| 伦理电影大哥的女人| 亚洲国产欧洲综合997久久,| 亚洲第一欧美日韩一区二区三区| 免费av不卡在线播放| 国产私拍福利视频在线观看| 禁无遮挡网站| 国产 一区 欧美 日韩| 国产极品精品免费视频能看的| 日本精品一区二区三区蜜桃| 成年免费大片在线观看| 网址你懂的国产日韩在线| 成年女人永久免费观看视频| 99热这里只有是精品在线观看 | 在线免费观看的www视频| 亚洲人成电影免费在线| 日本熟妇午夜| 国产三级黄色录像| 欧美区成人在线视频| 国产69精品久久久久777片| 亚洲第一区二区三区不卡| 最近中文字幕高清免费大全6 | 成人午夜高清在线视频| 亚洲在线观看片| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 亚洲久久久久久中文字幕| 中文字幕av成人在线电影| 在线观看午夜福利视频| 一个人看视频在线观看www免费| 久久久久久久午夜电影| 国产色爽女视频免费观看| 国产真实伦视频高清在线观看 | 欧美日韩乱码在线| 九色国产91popny在线| 人人妻,人人澡人人爽秒播| 男插女下体视频免费在线播放| 韩国av一区二区三区四区| 宅男免费午夜| 老司机福利观看| 国产高清视频在线播放一区| 欧美性感艳星| 国产精品永久免费网站| 在线观看av片永久免费下载| 99久久无色码亚洲精品果冻| 看片在线看免费视频| 婷婷六月久久综合丁香| 很黄的视频免费| 99视频精品全部免费 在线| 亚洲最大成人av| 老女人水多毛片| 日韩欧美精品v在线| 亚洲国产精品成人综合色| 美女被艹到高潮喷水动态| eeuss影院久久| 两人在一起打扑克的视频| 日韩中文字幕欧美一区二区| 99国产综合亚洲精品| 亚洲av成人精品一区久久| 自拍偷自拍亚洲精品老妇| 国产人妻一区二区三区在| 色综合婷婷激情| 亚洲激情在线av| 夜夜夜夜夜久久久久| 麻豆成人av在线观看| 亚洲狠狠婷婷综合久久图片| 国产久久久一区二区三区| 亚洲 国产 在线| 午夜福利在线观看吧| 看黄色毛片网站| 国产精品,欧美在线| aaaaa片日本免费| 日韩欧美 国产精品| 亚洲精华国产精华精| 赤兔流量卡办理| 欧美黄色片欧美黄色片| aaaaa片日本免费| 91久久精品电影网| 国产精品久久视频播放| 91狼人影院| 国产精品日韩av在线免费观看| 久久人妻av系列| 此物有八面人人有两片| 男女视频在线观看网站免费| 真人做人爱边吃奶动态| 美女大奶头视频| 真人一进一出gif抽搐免费| 免费av不卡在线播放| 最近中文字幕高清免费大全6 | 亚洲精品久久国产高清桃花| 午夜日韩欧美国产| 国产精品免费一区二区三区在线| 老司机午夜福利在线观看视频| 99久久久亚洲精品蜜臀av| 桃红色精品国产亚洲av| 成人高潮视频无遮挡免费网站| 亚洲av电影在线进入| 脱女人内裤的视频| 丁香欧美五月| 高清毛片免费观看视频网站| 国产精品电影一区二区三区| 少妇被粗大猛烈的视频| aaaaa片日本免费| 亚洲经典国产精华液单 | 在线播放国产精品三级| av在线观看视频网站免费| 精品熟女少妇八av免费久了| 精品一区二区三区视频在线观看免费| 国产成人影院久久av| 国产精品99久久久久久久久| 国产亚洲欧美在线一区二区| 国产 一区 欧美 日韩| www.色视频.com| 欧美极品一区二区三区四区| 人人妻人人澡欧美一区二区| 午夜精品久久久久久毛片777| 悠悠久久av| 亚洲精华国产精华精| 成人国产综合亚洲| 我要搜黄色片| 91九色精品人成在线观看| 舔av片在线| 欧美乱色亚洲激情| 亚洲av一区综合| 老司机深夜福利视频在线观看| 波多野结衣高清作品| 国产在视频线在精品| 久久这里只有精品中国| 99久久精品国产亚洲精品| 99国产极品粉嫩在线观看| 亚洲成人久久爱视频| 亚洲电影在线观看av| 18美女黄网站色大片免费观看| 黄色一级大片看看| 亚洲七黄色美女视频| 特级一级黄色大片| 51午夜福利影视在线观看| 欧美绝顶高潮抽搐喷水| 我要搜黄色片| 两性午夜刺激爽爽歪歪视频在线观看| 一级毛片久久久久久久久女| 欧美高清成人免费视频www| 露出奶头的视频| 日本一本二区三区精品| 国产精品久久久久久精品电影| 欧美高清性xxxxhd video| 3wmmmm亚洲av在线观看| 免费黄网站久久成人精品 | 午夜精品久久久久久毛片777| 国产高潮美女av| 可以在线观看的亚洲视频| 不卡一级毛片| 欧美成狂野欧美在线观看| 日本黄色视频三级网站网址| 男人的好看免费观看在线视频| 欧美日韩综合久久久久久 | 国产高清激情床上av| 亚洲美女搞黄在线观看 | 免费黄网站久久成人精品 | 国产成人欧美在线观看| 亚洲自拍偷在线| 亚洲精品456在线播放app | 亚洲成人精品中文字幕电影| 色视频www国产| 亚洲真实伦在线观看| 热99re8久久精品国产| 欧美+亚洲+日韩+国产| 无人区码免费观看不卡| 欧美+亚洲+日韩+国产| 中文字幕精品亚洲无线码一区| 国产乱人视频| 少妇人妻精品综合一区二区 | 欧美最黄视频在线播放免费| 亚洲熟妇熟女久久| 国产精品野战在线观看| 琪琪午夜伦伦电影理论片6080| 赤兔流量卡办理| 国产精品精品国产色婷婷| 一本综合久久免费| 亚洲av第一区精品v没综合| 成人永久免费在线观看视频| 99国产综合亚洲精品| 欧美乱色亚洲激情| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧|