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

    高速列車前端多胞吸能結(jié)構(gòu)的耐撞性優(yōu)化

    2017-06-19 19:35:43張秧聰鄧雯苑車全偉
    振動(dòng)與沖擊 2017年12期
    關(guān)鍵詞:撞擊力薄壁列車

    張秧聰, 許 平, 彭 勇, 鄧雯苑,2, 車全偉

    (1. 中南大學(xué) 軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙 410075; 2. 廣東技術(shù)師范學(xué)院 汽車學(xué)院,廣州 510000)

    高速列車前端多胞吸能結(jié)構(gòu)的耐撞性優(yōu)化

    張秧聰1, 許 平1, 彭 勇1, 鄧雯苑1,2, 車全偉1

    (1. 中南大學(xué) 軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,長(zhǎng)沙 410075; 2. 廣東技術(shù)師范學(xué)院 汽車學(xué)院,廣州 510000)

    為設(shè)計(jì)具有良好耐撞性能的高速列車前端多胞吸能結(jié)構(gòu),基于顯式動(dòng)力學(xué)有限元軟件LS-DYNA,建立此吸能結(jié)構(gòu)的有限元模型。通過(guò)臺(tái)車碰撞試驗(yàn)驗(yàn)證了有限元模型的準(zhǔn)確性,結(jié)合驗(yàn)證的有限元模型與全因子試驗(yàn)設(shè)計(jì),構(gòu)造了吸能結(jié)構(gòu)的比吸能SEA和撞擊平臺(tái)力關(guān)于設(shè)計(jì)參數(shù)單元胞邊長(zhǎng)和壁厚的Kriging代理模型,并進(jìn)行了誤差分析,采用多目標(biāo)粒子群優(yōu)化算法,對(duì)多胞吸能結(jié)構(gòu)的截面尺寸和厚度進(jìn)行了優(yōu)化設(shè)計(jì)。結(jié)果表明,壁厚比單元胞邊長(zhǎng)對(duì)多胞吸能結(jié)構(gòu)耐撞性影響更顯著,通過(guò)合理匹配壁厚和邊長(zhǎng),能有效提高撞擊平臺(tái)力和比吸能。

    固體力學(xué);多胞結(jié)構(gòu);數(shù)值仿真;碰撞試驗(yàn);多目標(biāo)優(yōu)化

    隨著高速列車運(yùn)行速度提高[1],一旦發(fā)生碰撞事故,將造成重大人員傷亡和財(cái)產(chǎn)損失。2011-07-23T20:30,甬溫線發(fā)生高速列車追尾事故,造成40多人死亡200多人受傷。為實(shí)現(xiàn)列車乘員安全保護(hù),減少事故損失,對(duì)列車端部吸能裝置的要求也越來(lái)越高,所以設(shè)計(jì)和優(yōu)化合理的列車端部吸能裝置已經(jīng)成為列車耐撞性研究中的熱點(diǎn)問(wèn)題。

    理想的列車吸能結(jié)構(gòu)應(yīng)當(dāng)位于車體的前后部分,在可控制的變形區(qū)域內(nèi)發(fā)生塑性變形,吸收撞擊動(dòng)能,同時(shí)保障乘客區(qū)域不發(fā)生嚴(yán)重破壞,并且在碰撞過(guò)程中不會(huì)產(chǎn)生過(guò)大的撞擊力峰值,使撞擊減速度在人體的承受范圍內(nèi)[2]。薄壁結(jié)構(gòu)因其吸能效率高、重量輕、成本低等特點(diǎn),被廣泛應(yīng)用于車輛、輪船、航空航天等領(lǐng)域[3],對(duì)高速列車來(lái)說(shuō),采用薄壁管作為吸能元件是比較合適的。國(guó)內(nèi)外學(xué)者已經(jīng)在列車吸能結(jié)構(gòu)方面做了很多工作并取得了一定成果[4-9],姚松等[10-14]對(duì)吸能結(jié)構(gòu)中常用的典型薄壁結(jié)構(gòu)碰撞吸能的影響因素如結(jié)構(gòu)橫截面形狀、結(jié)構(gòu)尺寸、材料特性參數(shù)、結(jié)構(gòu)預(yù)變形等進(jìn)行了研究,探討了不同的因素對(duì)結(jié)構(gòu)碰撞吸能特性的影響規(guī)律。

    過(guò)去的研究多集中于圓管、方管等單胞吸能結(jié)構(gòu),很少有關(guān)于多胞吸能結(jié)構(gòu)方面的研究。高廣軍等[15]對(duì)六邊形蜂窩吸能結(jié)構(gòu)進(jìn)行了撞擊數(shù)值仿真,發(fā)現(xiàn)組合構(gòu)型多胞薄壁管的耐撞性要優(yōu)于單胞吸能結(jié)構(gòu);文桂林等[16]基于多項(xiàng)式代理模型和多目標(biāo)粒子群優(yōu)化算法,對(duì)不同截面形狀的多胞結(jié)構(gòu)進(jìn)行了耐撞性優(yōu)化設(shè)計(jì);亓昌等[17]采用有限元仿真,以比能量SEA和撞擊峰值力Fp為評(píng)價(jià)指標(biāo),研究了不同沖擊角度下的錐形多胞薄壁方管的吸能特性。

    基于以上研究,可以看出多胞薄壁吸能結(jié)構(gòu)具有廣泛的應(yīng)用前景,本文以實(shí)際的工程項(xiàng)目為背景,以某型號(hào)高速動(dòng)車組前端多胞吸能結(jié)構(gòu)為研究對(duì)象,首先對(duì)吸能結(jié)構(gòu)進(jìn)行數(shù)值仿真研究其耐撞性,接著通過(guò)臺(tái)車沖擊試驗(yàn)驗(yàn)證數(shù)值仿真結(jié)果。在經(jīng)過(guò)驗(yàn)證的數(shù)值仿真的基礎(chǔ)上,結(jié)合Kriging近似模型和多目標(biāo)粒子群優(yōu)化算法,以最大比吸能SEA(Specific Energy Absorption)和撞擊平臺(tái)力Fm為優(yōu)化目標(biāo),對(duì)高速列車前端吸能結(jié)構(gòu)進(jìn)行耐撞性優(yōu)化,并確定合理的截面尺寸參數(shù)。

    1 模型的建立

    1.1 幾何模型

    高速動(dòng)車組的前端吸能結(jié)構(gòu),如圖1所示。是以八邊形管為基礎(chǔ),正六邊形管圍繞在其四周,形成的五孔組合蜂窩結(jié)構(gòu)形式[18]。正六邊形單元胞的邊長(zhǎng)a為56 mm,八邊形管上、下邊長(zhǎng)b為56 mm,側(cè)邊長(zhǎng)c為51 mm,斜邊與側(cè)邊的夾角為150°,上、下邊與斜邊的夾角為120°,橫截面的大小為280 mm×245 mm。五孔吸能結(jié)構(gòu)薄壁厚度為5 mm,長(zhǎng)度為500 mm。

    圖1 高速動(dòng)車組前端吸能結(jié)構(gòu)Fig.1 High-speed train front-end energy-absorbing structure

    1.2 有限元模型

    1.2.1 有限元模型的定義

    為了研究吸能結(jié)構(gòu)的耐撞性,本文使用非線性有限元軟件LS-DYNA對(duì)其進(jìn)行數(shù)值仿真分析。利用Hypermesh前處理軟件對(duì)吸能結(jié)構(gòu)的幾何模型進(jìn)行網(wǎng)格劃分,得到如圖2所示的有限元模型,包括固定剛性墻、吸能結(jié)構(gòu)、固定軌道、試驗(yàn)臺(tái)車四個(gè)部分。為了準(zhǔn)確模擬吸能結(jié)構(gòu)的變形,胞壁均采用5mm×5mm的Beltschko-Tsay四邊形殼單元,單元厚度方向采用3點(diǎn)積分,面內(nèi)采用單點(diǎn)積分。同時(shí)為了提高計(jì)算效率,試驗(yàn)臺(tái)車和軌道的單元尺寸設(shè)為20mm×20mm。

    在模擬碰撞時(shí),吸能結(jié)構(gòu)的自接觸采用“AUTOMATIC_SINGLE_SURFACE”接觸算法,吸能結(jié)構(gòu)與臺(tái)車、輪對(duì)與軌道之間采用“AUTOMATIC_SURFACE_TO_SURFACE”接觸算法。接觸的靜摩擦因數(shù)設(shè)為0.3,動(dòng)摩擦因數(shù)設(shè)為0.1。吸能結(jié)構(gòu)固定在質(zhì)量為3 000kg的臺(tái)車前端,臺(tái)車在固定軌道上以初始沖擊速度v0=15.51m/s撞擊固定剛性墻。

    圖2 有限元模型及邊界條件Fig.2 FEM and boundary conditions

    1.2.2 材料模型

    吸能結(jié)構(gòu)的材料為6008系鋁合金,為了獲得材料的力學(xué)性能,從相同的吸能結(jié)構(gòu)上截取一塊標(biāo)準(zhǔn)件,使用MTS647液壓萬(wàn)能材料試驗(yàn)機(jī)進(jìn)行準(zhǔn)靜態(tài)拉伸試驗(yàn),如圖3所示。因金屬材料在彈性階段應(yīng)力隨應(yīng)變呈線性變化,試驗(yàn)時(shí)采用恒定載荷增量加載以縮短加載時(shí)間,加載載荷增量為20kN/min;當(dāng)進(jìn)入初始塑性后,采用15mm/min的恒定速度加載。得到材料的應(yīng)力-應(yīng)變曲線,如圖4所示。材料的參數(shù),如表1所示。使用LS-DYNA材料庫(kù)中的Mat.024材料模型定義6008系鋁合金的材料屬性。為了減少計(jì)算時(shí)間,在計(jì)算精度不降低的情況下,臺(tái)車和鋼軌使用“Mat.020-MAT_RIGID”材料模型定義屬性。

    表1 吸能結(jié)構(gòu)材料參數(shù)Tab.1 Energy-absorbing structurematerial properties

    圖3 準(zhǔn)靜態(tài)拉伸試驗(yàn)Fig.3 Quasi-tensile test

    圖4 材料拉伸試驗(yàn)數(shù)據(jù)Fig.4 Tensile test data for the material

    1.2.3 耐撞性的評(píng)價(jià)指標(biāo)

    高速列車前端吸能結(jié)構(gòu)在碰撞發(fā)生時(shí),主要通過(guò)壓縮變形吸收沖擊動(dòng)能,耐撞性優(yōu)化目標(biāo)是得到具有最大吸能量和最小峰值力的結(jié)構(gòu)參數(shù)。撞擊峰值力(Fp),撞擊平臺(tái)力(Fm),變形位移(S),吸能量(EA)都是廣泛用于吸能結(jié)構(gòu)的耐撞性參數(shù)。

    吸能量(EA)是衡量吸能結(jié)構(gòu)吸收撞擊能量的一個(gè)重要指標(biāo),在數(shù)學(xué)上定義如下:

    (1)

    式中:F(s)為瞬時(shí)撞擊力,是位移S的函數(shù)。

    撞擊平臺(tái)力(Fm)為衡量吸能結(jié)構(gòu)緩沖變形的能力,可以通過(guò)吸能量(EA)除以壓縮位移S得到:

    (2)

    另外比吸能SEA用來(lái)表示吸能結(jié)構(gòu)單位質(zhì)量的吸能,可以通過(guò)吸能量(EA)除以質(zhì)量m得到:

    (3)

    2 有限元模型的試驗(yàn)驗(yàn)證

    2.1 試驗(yàn)設(shè)置

    為了觀察吸能結(jié)構(gòu)在高速?zèng)_擊下的撞擊特性和驗(yàn)證數(shù)值仿真模型的準(zhǔn)確性,對(duì)高速列車前端吸能結(jié)構(gòu)進(jìn)行全尺寸的動(dòng)態(tài)沖擊試驗(yàn),如圖5所示。電機(jī)通過(guò)鋼絲繩牽引3 000kg的臺(tái)車加速,當(dāng)速度達(dá)到15.51m/s時(shí),電機(jī)收到信號(hào)停止轉(zhuǎn)動(dòng),臺(tái)車依靠慣性沖擊剛性墻。在軌道的右側(cè)布置有NAC MEMRECAM HX-3高速攝像機(jī)(采樣頻率為50k fps),記錄臺(tái)車和吸能結(jié)構(gòu)的變形過(guò)程和位移。固定剛性墻安裝有測(cè)力傳感器,可以采集瞬時(shí)撞擊力。

    圖5 試驗(yàn)臺(tái)車與吸能結(jié)構(gòu)Fig.5 Crash trolley and energy-absorbing structure

    2.2 結(jié)果分析

    從圖6可知,臺(tái)車沖擊試驗(yàn)和有限元仿真中吸能結(jié)構(gòu)的變形基本一致,得到的撞擊力-時(shí)間曲線,如圖7所示。兩條曲線的趨勢(shì)基本一致。試驗(yàn)和仿真的撞擊力分別在0.001 s和0.000 9 s達(dá)到第一個(gè)峰值點(diǎn),大小分別為1 313.25 kN和1 277.11 kN,隨后急劇下降,然后在873.871 kN和890.385 kN上下波動(dòng)起伏,最后在0.032 s和0.035 s開(kāi)始降為0,試驗(yàn)和仿真的撞擊持續(xù)時(shí)間分別0.039 5 s和0.041 s。由于實(shí)際中試驗(yàn)條件要比仿真更加復(fù)雜,所以撞擊力-時(shí)間曲線的每個(gè)峰值點(diǎn)不能完全吻合,但是峰值點(diǎn)的數(shù)量和撞擊力的振幅基本一致。吸能結(jié)構(gòu)沖擊試驗(yàn)和仿真的結(jié)果對(duì)比與誤差見(jiàn)表2,撞擊峰值力Fp、撞擊平臺(tái)力Fm、比吸能SEA、撞擊持續(xù)時(shí)間誤差都很小,從而驗(yàn)證了有限元模型的準(zhǔn)確性。

    圖6 臺(tái)車沖擊試驗(yàn)和有限 元仿真變形對(duì)比圖Fig.6 Contrast diagram of car crash tests and finite element simulation deformation contrast diagram

    圖7 試驗(yàn)與仿真的 撞擊力-時(shí)間曲線Fig.7 Crash force-time curves for experimental and simulation表2 試驗(yàn)與數(shù)值仿真關(guān)鍵參數(shù)對(duì)比Tab.2 Comparison of critical indexes between the experimental and the numerical simulations

    參數(shù)峰值力Fp/kN平臺(tái)力Fm/kN持續(xù)時(shí)間/s比吸能SEA/(kJ·kg)試驗(yàn)1313.25873.8710.039529.482仿真1277.11890.3850.04131.103誤差/%-2.751.890.0385.5

    3 吸能結(jié)構(gòu)的耐撞性優(yōu)化

    3.1 優(yōu)化問(wèn)題的提出

    在實(shí)際工程應(yīng)用中,作為高速列車的主吸能元件,一方面要求吸能效率高,在一定范圍內(nèi),撞擊平臺(tái)力Fm越高越好,另一方面,高速列車前端吸能結(jié)構(gòu)還作為車輛承載結(jié)構(gòu),其撞擊平臺(tái)力Fm越高,承載能力就越強(qiáng)。所以本文以最大撞擊平臺(tái)力Fm和最大比吸能SEA為優(yōu)化目標(biāo),對(duì)吸能結(jié)構(gòu)進(jìn)行多目標(biāo)耐撞性優(yōu)化,選取壁厚t和正六邊形單元胞邊長(zhǎng)a作為設(shè)計(jì)變量,a和t的取值范圍分別是40 mm≤a≤60 mm和3 mm≤t≤5 mm,所以吸能結(jié)構(gòu)的耐撞性優(yōu)化用數(shù)學(xué)公式表達(dá)如下:

    (4)

    3.2 試驗(yàn)設(shè)計(jì)

    試驗(yàn)設(shè)計(jì)(DOE),就是構(gòu)建近似模型的采樣策略,決定了采樣點(diǎn)的數(shù)量和空間分布[19]。全因子設(shè)計(jì)、正交試驗(yàn)設(shè)計(jì)、中心復(fù)合設(shè)計(jì)(CCD)、Box-Behnken設(shè)計(jì)、拉丁超立方采樣、哈默斯雷采樣這些試驗(yàn)設(shè)計(jì)類型被廣泛應(yīng)用到DOE。吸能結(jié)構(gòu)具有壁厚t和正六邊形單元胞邊長(zhǎng)a兩個(gè)因子,兩個(gè)因子各有5個(gè)取值水平,考慮所有可能的因子組合及因子的每個(gè)取值水平,采用全因子設(shè)計(jì),所需要的試驗(yàn)次數(shù)為25次。

    3.3 Kriging近似模型與誤差分析

    以全因子試驗(yàn)設(shè)計(jì)的仿真結(jié)果為樣本點(diǎn),基于Kriging法構(gòu)建撞擊平臺(tái)力Fm和比吸能SEA的近似模型,如圖8所示。撞擊平臺(tái)力Fm隨著壁厚t的增大而增大,而且曲線變化趨勢(shì)明顯,說(shuō)明壁厚對(duì)平臺(tái)力Fm的影響顯著;撞擊平臺(tái)力Fm隨著正六邊形單元胞邊長(zhǎng)a增大而增大,但是變化趨勢(shì)不明顯,撞擊平臺(tái)力Fm在達(dá)到1 000 kN時(shí)開(kāi)始趨于平穩(wěn)。比吸能SEA隨壁厚t增大而增大,增大趨勢(shì)顯著;比吸能SEA隨正六邊形單元胞邊長(zhǎng)a增加而減小。

    圖8 撞擊平臺(tái)力Fm和比吸能SEA的響應(yīng)面Fig.8 Response surface of Fm and SEA

    為檢驗(yàn)Kriging代理模型的擬合精度,從25個(gè)樣本點(diǎn)中隨機(jī)抽取10個(gè)檢測(cè)樣本進(jìn)行誤差分析,擬合精度分析,見(jiàn)表3。相對(duì)誤差RE曲線,如圖9所示。相對(duì)誤差RE可以通過(guò)如下公式計(jì)算得到:

    (5)

    表3 擬合精度分析結(jié)果Tab.3 Fitting accuracy analysis results

    從圖7和表3可知,Kriging模型擬合精度很好,撞擊平臺(tái)力Fm擬合的相對(duì)誤差為0.087%~6.445%,比吸能SEA的擬合精度誤差為0.276%~1.806%,比吸能SEA擬合精度要比撞擊平臺(tái)力Fm高。

    3.4 多目標(biāo)粒子群優(yōu)化算法和流程

    基于Kriging法構(gòu)建撞擊平臺(tái)力Fm和比吸能SEA

    的近似模型,需要借助優(yōu)化算法對(duì)代理模型進(jìn)行尋優(yōu)。吸能結(jié)構(gòu)的耐撞性優(yōu)化是典型的多目標(biāo)優(yōu)化問(wèn)題,相比于傳統(tǒng)的多目標(biāo)遺傳算法NSGA,多目標(biāo)粒子群算法收斂速度快、產(chǎn)生Pareto前沿分布均勻、算法簡(jiǎn)單,被廣泛應(yīng)用于車輛設(shè)計(jì)中[20],所以本文采用多目標(biāo)粒子群算法得到目標(biāo)函數(shù)撞擊平臺(tái)力Fm和比吸能SEA的Pareto前沿。多目標(biāo)粒子群算法的流程,如圖10所示。

    圖9 相對(duì)誤差曲線Fig.9 Relative Error curves

    圖10 多目標(biāo)優(yōu)化算法流程圖Fig.10 Flow chart of multi-objective optimization algorithm

    3.5 優(yōu)化結(jié)果

    多目標(biāo)粒子群優(yōu)化得到的Pareto前沿,如圖11所示。理論上,Pareto前沿上的任意一點(diǎn)都可作為最佳點(diǎn),從圖中可知,提高撞擊平臺(tái)力Fm必然會(huì)減小比吸能SEA,高速列車前端吸能結(jié)構(gòu)作為承載式吸能結(jié)構(gòu),其承載能力在撞擊中尤為關(guān)鍵,設(shè)計(jì)要求撞擊平臺(tái)力Fm>1 000 kN,所以從符合條件的Pareto解集中選擇平臺(tái)力Fm最大的優(yōu)化結(jié)果,對(duì)應(yīng)的最優(yōu)結(jié)構(gòu)參數(shù)壁厚t和正六邊形單元胞邊長(zhǎng)a分別為5 mm和60 mm。優(yōu)化前后吸能結(jié)構(gòu)的撞擊力-時(shí)間曲線、優(yōu)化結(jié)果,如圖12和表4所示。優(yōu)化后撞擊平臺(tái)力從973.871 kN提高至1 004.24 kN,不僅滿足了設(shè)計(jì)要求的平臺(tái)力限制,在相同壓縮位移內(nèi)還提高了吸能量;比吸能從29.482提高至31.119,說(shuō)明優(yōu)化后的吸能結(jié)構(gòu)單位質(zhì)量吸能量更高。

    表4 吸能結(jié)構(gòu)優(yōu)化結(jié)果Tab.4 The optimization results of energy absorption

    圖11 Pareto前沿Fig.11 Pareto fronts

    圖12 優(yōu)化后撞擊力-時(shí)間曲線Fig.12 The optimized crash force-time curve

    4 結(jié) 論

    本文為研究和優(yōu)化高速列車前端結(jié)構(gòu)的耐撞特性,建立吸能結(jié)構(gòu)的有限元模型,基于LS-DYNA進(jìn)行數(shù)值仿真。為了驗(yàn)證有限元模型的準(zhǔn)確性,以15.51 m/s的速度進(jìn)行臺(tái)車沖擊試驗(yàn),試驗(yàn)和仿真的結(jié)果吻合得很好。經(jīng)驗(yàn)證的吸能結(jié)構(gòu)有限元模型,能有效提高分析效率。

    以最大化撞擊平臺(tái)力Fm和比吸能SEA為優(yōu)化目標(biāo),采用全因子試驗(yàn)設(shè)計(jì)得到25組樣本點(diǎn),基于Kriging構(gòu)件撞擊平臺(tái)力Fm和比吸能SEA兩個(gè)響應(yīng)的近似模型,從兩個(gè)響應(yīng)面模型可以看出壁厚t相比單元胞邊長(zhǎng)a對(duì)多胞吸能結(jié)構(gòu)的撞擊平臺(tái)力Fm和比吸能SEA影響更顯著。隨機(jī)抽取10組樣本進(jìn)行誤差分析,撞擊平臺(tái)力Fm擬合的相對(duì)誤差為0.087%~6.445%,比吸能SEA的擬合的相對(duì)誤差為0.276%~1.806%,驗(yàn)證了Kriging模型的擬合精度。

    在撞擊平臺(tái)力Fm和比吸能SEA兩個(gè)優(yōu)化目標(biāo)的Kriging模型基礎(chǔ)上,基于多目標(biāo)粒子群優(yōu)化算法對(duì)吸能結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì),得到最優(yōu)撞擊平臺(tái)力Fm和比吸能SEA的Pareto前沿,優(yōu)化后結(jié)構(gòu)的撞擊平臺(tái)力Fm達(dá)到1 004.24 kN,比吸能SEA達(dá)到31.119 kJ/kg。通過(guò)合理匹配壁厚和單元胞邊長(zhǎng),能有效提高撞擊平臺(tái)力和比吸能。

    [1] 張志新, 肖守訥, 陽(yáng)光武, 等. 高速列車乘員碰撞安全性研究[J]. 鐵道學(xué)報(bào), 2013, 35(10): 24-32. ZHANG Zhixin, XIAO Shoune, YANG Guangwu, et al. Research on collision safety of high-speed train crews & passengers[J]. Journal of the China Railway Society, 2013, 35(10): 24-32.

    [2] 劉建新, 趙華. 高速動(dòng)力車能量吸收裝置[J]. 鐵道學(xué)報(bào), 1997, 19(3): 32-36. LIU Jianxin, ZHAO Hua. Energy absorbing devices for high-speed locomotive under crash[J]. Journal of the China Railway Society, 1997, 19(3): 32-36.

    [3] 蔣致禹, 顧敏童, 趙永生. 一種薄壁吸能結(jié)構(gòu)的設(shè)計(jì)優(yōu)化[J]. 振動(dòng)與沖擊, 2010, 29(2): 111-116. JIANG Zhiyu, GU Mintong, ZHAO Yongsheng. A thin-walled energy-absorbing structure design optimization[J]. Journal of Vibration and Shock, 2010, 29(2): 111-116.

    [4] LEWIS J H, RASAIAH W G, SCHOLES A. Validation of measures to improve rail vehicle crashworthiness[J].Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 1996, 210(2): 73-85.

    [5] LU G. Energy absorption requirement for crashworthy vehicles[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2002, 216(1): 31-39.

    [6] 李松晏, 鄭志軍. 高速列車吸能結(jié)構(gòu)設(shè)計(jì)和耐撞性分析[J]. 爆炸與沖擊, 2015,35(2):164-170. LI Songyan, ZHENG Zhijun. Energy-absorbing structure design and crashworthiness analysis of high-speed trains[J]. Explosion and Shock Waves, 2015,35(2): 164-170.

    [7] 雷成, 肖守訥, 羅世輝. 基于顯式有限元的高速列車吸能裝置吸能原理研究[J]. 鐵道機(jī)車車輛, 2012, 32(2): 1-5. LEI Cheng, XIAO Shoune, LUO Shihui. Research on the energy-absorbing theory of high-speed train energy-absorbing component based on the explicit finite element[J]. Railway Locomotive & Car, 2012, 32(2): 1-5.

    [8] 陳秉智, 張向海, 馬紀(jì)軍, 等. 高速動(dòng)車組被動(dòng)安全性和耐撞性研究[J]. 計(jì)算力學(xué)學(xué)報(bào), 2011, 28(1): 152-158. CHEN Bingzhi, ZHANG Xianghai, MA Jijun, et al. High-speed EMUs passive safety and crashworthiness study[J]. Chinese Journal of Computational Mechanics, 2011, 28(1): 152-158.

    [9] 張志新. 高速列車耐撞性結(jié)構(gòu)及安全性研究[D]. 成都:西南交通大學(xué), 2012.

    [10] 姚松, 田紅旗. 車輛吸能部件的薄壁結(jié)構(gòu)碰撞研究[J]. 中國(guó)鐵道科學(xué), 2001, 22(2): 55-60. YAO Song, TIAN Hongqi. Research on vehicle energy absorption component of thin-walled structure under impact[J] China Railway Science, 2001, 22(2): 55-60.

    [11] 栗蔭帥. 車輛薄壁結(jié)構(gòu)碰撞吸能特性分析與改進(jìn)[D]. 大連:大連理工大學(xué), 2007.

    [12] 賈宇, 肖守訥. 耐碰撞車體吸能裝置的薄壁結(jié)構(gòu)研究[J]. 鐵道車輛, 2005, 43(5): 6-10. JIA Yu, XIAO Shoune. Research on the thin-wall structure of the energy absorption equipment on the collision resistant car-body[J]. Rolling Stock, 2005, 43(5): 6-10.

    [13] 謝素超. 耐沖擊地鐵車輛吸能結(jié)構(gòu)研究[D]. 長(zhǎng)沙:中南大學(xué), 2007.

    [14] 張平, 馬建, 那景新. 波紋管耐撞性的多目標(biāo)優(yōu)化[J]. 振動(dòng)與沖擊, 2015, 34(15): 12-16. ZHANG Ping, MA Jian, NA Jingxin. Multi-objective optimization on the crashworthiness of corrugated tube[J]. Journal of Vibration and Shock,2015, 34(15): 12-16.

    [15] 高廣軍, 姚松. 車輛薄壁結(jié)構(gòu)撞擊吸能特性研究[J]. 鐵道機(jī)車車輛, 2002 (2): 8-10. GAO Guangjun, YAO Song. Research on energy absorption attribute in collision of thin-walled structure between trains[J]. Railway Locomotive & Car, 2002 (2): 8-10.

    [16] 文桂林, 孔祥正, 尹漢鋒, 等. 泡沫填充夾芯墻多胞結(jié)構(gòu)的耐撞性多目標(biāo)優(yōu)化設(shè)計(jì)[J]. 振動(dòng)與沖擊, 2015, 34(5): 115-121. WEN Guilin, KONG Xiangzheng, YIN Hanfeng, et al. Multi-objective optimization design of form-filled sandwich wall multi-cell structure[J]. Journal of Vibration and Shock, 2015, 34(5): 115-121.

    [17] 亓昌, 董方亮, 楊姝, 等. 錐形多胞薄壁管斜向沖擊吸能特性仿真研究[J]. 振動(dòng)與沖擊, 2012, 31(24): 102-107. QI Chang, DONG Fangliang, YANG Shu, et al. Energy absorbing characteristics of a tapered multi-cell thin-walled tube under oblique impact[J]. Journal of Vibration and Shock,2012, 31(24): 102-107.

    [18] 張?jiān)谥? 姚曙光, 許平,等. 高速列車車體前端吸能結(jié)構(gòu)的碰撞仿真與試驗(yàn)[J]. 機(jī)車車輛工藝, 2015(3):33-34. ZHANG Zaizhong, YAO Shuguang, XU Ping, et al. High-speed train car body front-end energy-absorbing structure of crash simulation and experiment[J]. Locomotive & Rolling Stock Technology, 2015(3):33-34.

    [19] BENYOUNIS K Y, OLABI A G. Optimization of different welding processes using statistical and numerical approaches—A reference guide[J]. Advances in Engineering Software, 2008, 39(6): 483-496.

    [20] RAQUEL C R, NAVAL JR P C. An effective use of crowding distance in multiobjective particle swarm optimization[C]∥ Proceedings of the 7th Annual conference on Genetic and Evolutionary Computation. Washington, DC: ACM, 2005: 257-264.

    Crashworthiness optimization of high-speed train front multi-cell energy-absorbing structures

    ZHANG Yangcong1, XU Ping1,PENG Yong1,DENG Wenyuan1,2,CHE Quanwei1

    (1. Key Laboratory of Traffic Safety on Track of Ministry of Education, Central South University, Changsha 410075, China;2. School of Automotive Engineering, Guangdong Polytechnic Normal University, Guangzhou 510000, China)

    To design high-speed train front multi-cell energy-absorbing structures with excellent crashworthiness, a finite element model (FEM) of energy-absorbing structure was established using LS-DYNA based on explicit dynamic. The FEM was verified using trolley collision test. Combined with the validated FEM and full factorial experiment design, the Kriging surrogate models of specific energy absorption (SEA) and mean crushing force with respect to design parameters were respectively constructed, and error analysis was given. The section size and thickness of the energy-absorbing structure were optimized by using a multi-objective particle swarm optimization algorithm. The result shows that wall thickness has more obvious impact on crashworthiness of the multi-cell energy-absorbing structure than cell-wall length. The mean crushing force and SEA can be effectively improved by reasonable matching wall thickness and length.

    solid mechanics; multi-cellstructure; numerical simulation; impact experiments; multi-objective optimization

    國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(U1334208);國(guó)家支撐項(xiàng)目(2015BAG12B01);教育部重點(diǎn)項(xiàng)目(113051A)

    2016-01-11 修改稿收到日期: 2016-05-11

    張秧聰 男,碩士,1990年生

    彭勇 男,博士,教授,1984年生

    U270.2;O342

    A

    10.13465/j.cnki.jvs.2017.12.006

    猜你喜歡
    撞擊力薄壁列車
    登上末日列車
    關(guān)愛(ài)向列車下延伸
    磁流變液仿生薄壁吸能管及其耐撞性可控度的研究
    汽車工程(2021年12期)2021-03-08 02:34:06
    超小型薄壁線圈架注射模設(shè)計(jì)
    模具制造(2019年7期)2019-09-25 07:29:58
    穿越時(shí)空的列車
    接觸面對(duì)駁船撞擊橋墩動(dòng)力響應(yīng)的影響
    一種薄壁件的加工工藝介紹
    受撞橋梁結(jié)構(gòu)撞擊力仿真分析研究
    西去的列車
    樁基布置對(duì)高樁碼頭撞擊力分配的影響
    水道港口(2014年1期)2014-04-27 14:14:40
    欧美激情久久久久久爽电影| 91麻豆av在线| 国产高清有码在线观看视频| 黄片小视频在线播放| 丁香欧美五月| 午夜福利在线观看免费完整高清在 | 波多野结衣高清作品| 久久99热这里只有精品18| 国产不卡一卡二| 国产野战对白在线观看| ponron亚洲| 国产精品1区2区在线观看.| 在线播放无遮挡| 国产精品野战在线观看| 亚洲在线观看片| 欧美另类亚洲清纯唯美| 热99re8久久精品国产| 久99久视频精品免费| 麻豆成人午夜福利视频| av女优亚洲男人天堂| 欧美成狂野欧美在线观看| 99精品久久久久人妻精品| 性欧美人与动物交配| 性欧美人与动物交配| 狂野欧美激情性xxxx| 欧美另类亚洲清纯唯美| 夜夜夜夜夜久久久久| 国产精华一区二区三区| 一本一本综合久久| 欧美成狂野欧美在线观看| 99riav亚洲国产免费| 免费看十八禁软件| 99riav亚洲国产免费| 男人和女人高潮做爰伦理| 日韩亚洲欧美综合| tocl精华| 日本一二三区视频观看| 男人的好看免费观看在线视频| 一本久久中文字幕| 俺也久久电影网| 黄片大片在线免费观看| 免费看十八禁软件| 欧美成人一区二区免费高清观看| 亚洲欧美日韩卡通动漫| 色播亚洲综合网| 香蕉丝袜av| 一个人免费在线观看电影| 国产av麻豆久久久久久久| 18禁裸乳无遮挡免费网站照片| 免费人成视频x8x8入口观看| 成年女人永久免费观看视频| or卡值多少钱| 亚洲性夜色夜夜综合| 天美传媒精品一区二区| 国产精品久久视频播放| 免费看十八禁软件| 美女cb高潮喷水在线观看| 真人一进一出gif抽搐免费| 亚洲成人久久爱视频| 非洲黑人性xxxx精品又粗又长| 国产麻豆成人av免费视频| 欧美色欧美亚洲另类二区| 桃色一区二区三区在线观看| 国内精品久久久久久久电影| aaaaa片日本免费| 国产亚洲欧美98| 51国产日韩欧美| xxxwww97欧美| 精品午夜福利视频在线观看一区| 亚洲欧美激情综合另类| 美女 人体艺术 gogo| 国产单亲对白刺激| 亚洲av不卡在线观看| 婷婷丁香在线五月| 高清在线国产一区| 久久久久精品国产欧美久久久| 久久久久久久亚洲中文字幕 | 成年人黄色毛片网站| 久99久视频精品免费| 午夜激情福利司机影院| 一进一出抽搐动态| 99久久精品国产亚洲精品| 欧美日韩亚洲国产一区二区在线观看| svipshipincom国产片| 国产成人a区在线观看| 国产精品电影一区二区三区| 午夜免费男女啪啪视频观看 | 一级毛片高清免费大全| 日韩 欧美 亚洲 中文字幕| 少妇丰满av| 最新中文字幕久久久久| 国内揄拍国产精品人妻在线| 精品国产美女av久久久久小说| 丁香欧美五月| 久久久成人免费电影| 亚洲av免费在线观看| 亚洲国产欧美人成| 欧美成人性av电影在线观看| 国产精品久久视频播放| 精品欧美国产一区二区三| 男人和女人高潮做爰伦理| xxxwww97欧美| 国产精品 欧美亚洲| 淫妇啪啪啪对白视频| 淫妇啪啪啪对白视频| 99精品在免费线老司机午夜| 18+在线观看网站| 欧美乱码精品一区二区三区| 日韩国内少妇激情av| 欧美精品啪啪一区二区三区| 国产乱人伦免费视频| 日韩中文字幕欧美一区二区| 12—13女人毛片做爰片一| 亚洲美女黄片视频| 亚洲国产精品999在线| 免费av毛片视频| 免费高清视频大片| 欧美性猛交╳xxx乱大交人| 免费搜索国产男女视频| or卡值多少钱| 国产免费一级a男人的天堂| 91av网一区二区| 久久精品国产亚洲av涩爱 | 午夜亚洲福利在线播放| 老司机午夜福利在线观看视频| 久久久久精品国产欧美久久久| 免费人成在线观看视频色| 无限看片的www在线观看| 国产老妇女一区| 男女午夜视频在线观看| 高清日韩中文字幕在线| 99国产综合亚洲精品| 中亚洲国语对白在线视频| 亚洲成人久久爱视频| svipshipincom国产片| 天堂av国产一区二区熟女人妻| 无限看片的www在线观看| 国产一区二区三区视频了| av天堂在线播放| 97人妻精品一区二区三区麻豆| 香蕉av资源在线| 日本a在线网址| 男人的好看免费观看在线视频| 国产视频内射| 亚洲av成人精品一区久久| 伊人久久精品亚洲午夜| 日本黄色视频三级网站网址| 午夜免费成人在线视频| 九九久久精品国产亚洲av麻豆| 噜噜噜噜噜久久久久久91| 久久九九热精品免费| 国产精品亚洲av一区麻豆| 国产综合懂色| 无限看片的www在线观看| 动漫黄色视频在线观看| 国产午夜精品久久久久久一区二区三区 | 别揉我奶头~嗯~啊~动态视频| 欧美日韩亚洲国产一区二区在线观看| 国产激情偷乱视频一区二区| 老司机在亚洲福利影院| 日韩欧美 国产精品| 香蕉av资源在线| 丝袜美腿在线中文| 热99re8久久精品国产| 免费观看精品视频网站| 国产乱人伦免费视频| 国产男靠女视频免费网站| 国产成年人精品一区二区| 一夜夜www| 无人区码免费观看不卡| 一本精品99久久精品77| 岛国在线观看网站| 国产极品精品免费视频能看的| 一级a爱片免费观看的视频| 色吧在线观看| 亚洲av一区综合| 精品熟女少妇八av免费久了| 国产精品永久免费网站| 国产伦精品一区二区三区四那| 精品福利观看| 一进一出抽搐gif免费好疼| 草草在线视频免费看| 日本免费a在线| 日韩高清综合在线| e午夜精品久久久久久久| 黑人欧美特级aaaaaa片| 美女被艹到高潮喷水动态| 全区人妻精品视频| 一区福利在线观看| 亚洲国产欧美人成| 久久久久九九精品影院| 国产成人欧美在线观看| 成人欧美大片| 国产亚洲精品久久久久久毛片| 日本与韩国留学比较| 国产精品一区二区三区四区久久| 99国产精品一区二区蜜桃av| 在线观看免费午夜福利视频| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av涩爱 | 性欧美人与动物交配| 亚洲在线观看片| 小说图片视频综合网站| 日韩成人在线观看一区二区三区| 亚洲熟妇中文字幕五十中出| 最近最新免费中文字幕在线| 欧美丝袜亚洲另类 | 97人妻精品一区二区三区麻豆| 久久久久久久久久黄片| 久久亚洲真实| 国产精品亚洲av一区麻豆| 国产激情欧美一区二区| 精品电影一区二区在线| 18禁国产床啪视频网站| 51午夜福利影视在线观看| 国产欧美日韩一区二区精品| 国产高清三级在线| 欧美乱色亚洲激情| 熟女电影av网| 国产蜜桃级精品一区二区三区| 丰满的人妻完整版| 亚洲欧美日韩卡通动漫| 亚洲第一欧美日韩一区二区三区| 高潮久久久久久久久久久不卡| 又紧又爽又黄一区二区| 丁香欧美五月| 日韩欧美免费精品| 最近在线观看免费完整版| 中文字幕人妻熟人妻熟丝袜美 | 国产不卡一卡二| 国产一区二区三区在线臀色熟女| 最近在线观看免费完整版| 欧美一级a爱片免费观看看| 免费无遮挡裸体视频| 最新中文字幕久久久久| 中文字幕熟女人妻在线| 亚洲国产色片| 高清在线国产一区| 免费av不卡在线播放| 欧美bdsm另类| 免费大片18禁| 可以在线观看的亚洲视频| av黄色大香蕉| 欧美日韩精品网址| 国产真人三级小视频在线观看| 19禁男女啪啪无遮挡网站| 亚洲av熟女| 午夜激情欧美在线| 亚洲在线自拍视频| 精品不卡国产一区二区三区| 最新中文字幕久久久久| 亚洲狠狠婷婷综合久久图片| 熟女人妻精品中文字幕| 免费一级毛片在线播放高清视频| 一区二区三区激情视频| 国产伦精品一区二区三区视频9 | 亚洲avbb在线观看| 国产探花极品一区二区| 此物有八面人人有两片| or卡值多少钱| 97超级碰碰碰精品色视频在线观看| 窝窝影院91人妻| 国产乱人视频| 丰满人妻一区二区三区视频av | 在线观看日韩欧美| 怎么达到女性高潮| 99国产精品一区二区三区| 日本成人三级电影网站| 美女大奶头视频| 国产欧美日韩一区二区三| 一个人免费在线观看电影| 亚洲人与动物交配视频| 99久久精品一区二区三区| 精品电影一区二区在线| 亚洲熟妇熟女久久| 男人和女人高潮做爰伦理| 中文字幕久久专区| 国产高清视频在线播放一区| 欧美一区二区国产精品久久精品| 国产真实乱freesex| 免费在线观看日本一区| 国产精品永久免费网站| 高潮久久久久久久久久久不卡| 国产老妇女一区| 99在线人妻在线中文字幕| 中国美女看黄片| 黑人欧美特级aaaaaa片| av在线蜜桃| 99精品欧美一区二区三区四区| 青草久久国产| 精品久久久久久久人妻蜜臀av| 免费在线观看影片大全网站| 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 天堂影院成人在线观看| 欧美日韩精品网址| 丰满的人妻完整版| 十八禁人妻一区二区| 欧美一级a爱片免费观看看| 成人av一区二区三区在线看| 少妇的逼好多水| 久99久视频精品免费| 中文资源天堂在线| 美女大奶头视频| 亚洲国产精品成人综合色| www.熟女人妻精品国产| 国产老妇女一区| АⅤ资源中文在线天堂| 久久九九热精品免费| 国产高清videossex| 久久国产精品影院| 国产精品久久久久久精品电影| 岛国在线免费视频观看| 色av中文字幕| 欧美日韩精品网址| 精品国产美女av久久久久小说| 桃红色精品国产亚洲av| 免费人成在线观看视频色| 精品久久久久久,| 丝袜美腿在线中文| 国产 一区 欧美 日韩| 午夜福利在线观看免费完整高清在 | 超碰av人人做人人爽久久 | 天天一区二区日本电影三级| 一区福利在线观看| 欧美xxxx黑人xx丫x性爽| 一本精品99久久精品77| 在线观看免费视频日本深夜| 在线观看一区二区三区| 黄片大片在线免费观看| av专区在线播放| 成人三级黄色视频| 99国产精品一区二区蜜桃av| av黄色大香蕉| 国产精品99久久99久久久不卡| 国产精品免费一区二区三区在线| 欧美xxxx黑人xx丫x性爽| 精品久久久久久久久久免费视频| 日韩免费av在线播放| 变态另类丝袜制服| 成年免费大片在线观看| 日本熟妇午夜| 国产精品,欧美在线| 少妇人妻精品综合一区二区 | 99在线人妻在线中文字幕| 精品无人区乱码1区二区| 精品福利观看| АⅤ资源中文在线天堂| 亚洲精品456在线播放app | 欧美色欧美亚洲另类二区| 色老头精品视频在线观看| 午夜免费观看网址| 在线免费观看不下载黄p国产 | 99精品久久久久人妻精品| 人人妻人人看人人澡| 搡老岳熟女国产| 久久久久国内视频| 亚洲av中文字字幕乱码综合| 色在线成人网| 久久久久国内视频| 看黄色毛片网站| 最新中文字幕久久久久| 给我免费播放毛片高清在线观看| 91在线观看av| 精品国产三级普通话版| 日韩高清综合在线| 国产欧美日韩精品一区二区| 亚洲成人精品中文字幕电影| 桃红色精品国产亚洲av| 小蜜桃在线观看免费完整版高清| 一区福利在线观看| 亚洲人成电影免费在线| netflix在线观看网站| 国产一区在线观看成人免费| 青草久久国产| 欧美绝顶高潮抽搐喷水| 亚洲18禁久久av| 日韩有码中文字幕| 少妇高潮的动态图| 久久精品影院6| 亚洲欧美精品综合久久99| 亚洲电影在线观看av| 观看美女的网站| 国产久久久一区二区三区| 大型黄色视频在线免费观看| 精品免费久久久久久久清纯| 久久国产精品人妻蜜桃| 日韩免费av在线播放| 免费看美女性在线毛片视频| 日本在线视频免费播放| 日韩人妻高清精品专区| av在线天堂中文字幕| 最新美女视频免费是黄的| 亚洲真实伦在线观看| 午夜免费激情av| 搞女人的毛片| www.999成人在线观看| 久久草成人影院| 久久久久久九九精品二区国产| 日日干狠狠操夜夜爽| 高潮久久久久久久久久久不卡| 精品人妻1区二区| 性色avwww在线观看| 国产午夜精品论理片| 免费av毛片视频| 精品无人区乱码1区二区| 九九在线视频观看精品| 免费看美女性在线毛片视频| 黑人欧美特级aaaaaa片| 亚洲欧美日韩无卡精品| 高清毛片免费观看视频网站| 欧美绝顶高潮抽搐喷水| 99在线人妻在线中文字幕| 男人舔奶头视频| 国产精品 国内视频| 看片在线看免费视频| 国产91精品成人一区二区三区| 精品无人区乱码1区二区| 在线观看免费视频日本深夜| 99国产综合亚洲精品| 久久久精品大字幕| 欧美av亚洲av综合av国产av| 亚洲成av人片免费观看| 欧美+日韩+精品| 国产国拍精品亚洲av在线观看 | 久久精品国产综合久久久| 欧美午夜高清在线| 中出人妻视频一区二区| 精品国产超薄肉色丝袜足j| 69人妻影院| 日日摸夜夜添夜夜添小说| 婷婷六月久久综合丁香| 午夜激情福利司机影院| 欧美区成人在线视频| 国产黄色小视频在线观看| 国产精品野战在线观看| 他把我摸到了高潮在线观看| 久久久国产成人精品二区| 又粗又爽又猛毛片免费看| 十八禁人妻一区二区| 日韩国内少妇激情av| 国产精品日韩av在线免费观看| 国产在视频线在精品| 无遮挡黄片免费观看| 一进一出抽搐动态| 午夜两性在线视频| 俄罗斯特黄特色一大片| 欧美一级a爱片免费观看看| 精品国产三级普通话版| 无人区码免费观看不卡| 欧美区成人在线视频| 最近最新中文字幕大全免费视频| 看片在线看免费视频| 亚洲,欧美精品.| 熟女少妇亚洲综合色aaa.| 97碰自拍视频| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 香蕉丝袜av| 免费人成在线观看视频色| 精品乱码久久久久久99久播| 欧美成人免费av一区二区三区| 亚洲 国产 在线| 亚洲五月婷婷丁香| 久久久久久久精品吃奶| 国产精品嫩草影院av在线观看 | 少妇的逼好多水| 精品久久久久久成人av| 免费看a级黄色片| 久久久国产成人免费| 国产精品乱码一区二三区的特点| a级毛片a级免费在线| 天堂网av新在线| 国产精品久久久久久精品电影| 中出人妻视频一区二区| 人人妻人人澡欧美一区二区| 国产亚洲av嫩草精品影院| 网址你懂的国产日韩在线| 51午夜福利影视在线观看| 国产亚洲精品一区二区www| 久久久久久国产a免费观看| 色av中文字幕| 日本与韩国留学比较| 欧美日韩福利视频一区二区| 999久久久精品免费观看国产| 久久久久九九精品影院| 国产精品一区二区三区四区免费观看 | www.色视频.com| 狂野欧美激情性xxxx| 久久久久亚洲av毛片大全| 国产真实乱freesex| 老汉色∧v一级毛片| 天天躁日日操中文字幕| 伊人久久精品亚洲午夜| 极品教师在线免费播放| 国产探花在线观看一区二区| 久久中文看片网| 欧洲精品卡2卡3卡4卡5卡区| 搡老岳熟女国产| xxxwww97欧美| 男人舔奶头视频| 国产精品美女特级片免费视频播放器| 国产精品久久久久久久电影 | 黄色成人免费大全| 亚洲黑人精品在线| 好看av亚洲va欧美ⅴa在| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3| 他把我摸到了高潮在线观看| 欧洲精品卡2卡3卡4卡5卡区| 啦啦啦免费观看视频1| 久久久国产成人精品二区| 精品人妻1区二区| 日韩有码中文字幕| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 在线a可以看的网站| 欧美zozozo另类| 少妇熟女aⅴ在线视频| 国产精品女同一区二区软件 | 深爱激情五月婷婷| 真人做人爱边吃奶动态| 有码 亚洲区| 国产精品电影一区二区三区| 色吧在线观看| 亚洲天堂国产精品一区在线| avwww免费| 人妻久久中文字幕网| 97人妻精品一区二区三区麻豆| 免费在线观看成人毛片| 人人妻,人人澡人人爽秒播| 黄色片一级片一级黄色片| 亚洲 国产 在线| 成人av在线播放网站| 91久久精品电影网| 欧美色视频一区免费| 午夜福利在线观看免费完整高清在 | 亚洲最大成人中文| 91在线精品国自产拍蜜月 | 国产一区二区三区在线臀色熟女| 99久久精品一区二区三区| 国产高清视频在线观看网站| 高潮久久久久久久久久久不卡| 变态另类丝袜制服| 精品久久久久久久人妻蜜臀av| 变态另类丝袜制服| 看片在线看免费视频| 丰满乱子伦码专区| 18+在线观看网站| 国产一区二区亚洲精品在线观看| 日韩大尺度精品在线看网址| 天堂动漫精品| 日韩大尺度精品在线看网址| 19禁男女啪啪无遮挡网站| 丁香六月欧美| 久久精品综合一区二区三区| 精品久久久久久久人妻蜜臀av| 国产一区二区激情短视频| 日日夜夜操网爽| 久久精品综合一区二区三区| 国产伦精品一区二区三区四那| 欧美三级亚洲精品| 国产亚洲精品久久久com| 亚洲电影在线观看av| 亚洲五月婷婷丁香| 午夜两性在线视频| 毛片女人毛片| 一区二区三区高清视频在线| 在线观看免费午夜福利视频| 午夜福利在线在线| 老司机福利观看| 国产欧美日韩一区二区三| 亚洲在线观看片| 噜噜噜噜噜久久久久久91| 国产精品亚洲美女久久久| 一个人观看的视频www高清免费观看| 熟女电影av网| 亚洲成人久久爱视频| 午夜福利18| 亚洲国产中文字幕在线视频| 97超视频在线观看视频| 波多野结衣巨乳人妻| 亚洲av第一区精品v没综合| 99视频精品全部免费 在线| 国产欧美日韩一区二区精品| 亚洲精品亚洲一区二区| 亚洲成a人片在线一区二区| av国产免费在线观看| www.熟女人妻精品国产| 在线免费观看的www视频| 人妻丰满熟妇av一区二区三区| 久久久久久国产a免费观看| 国产亚洲精品综合一区在线观看| 香蕉丝袜av| 色老头精品视频在线观看| 亚洲精华国产精华精| 动漫黄色视频在线观看| 欧美日韩一级在线毛片| 亚洲精品一区av在线观看| 婷婷丁香在线五月| 99久久精品热视频| 亚洲人成网站高清观看| 国产一区二区亚洲精品在线观看| 搡老熟女国产l中国老女人| 两个人视频免费观看高清| av福利片在线观看| 国产亚洲av嫩草精品影院| 国产爱豆传媒在线观看| 国产日本99.免费观看| 国产不卡一卡二| 九色国产91popny在线| 欧美一区二区国产精品久久精品| 少妇的逼水好多| 亚洲精品国产精品久久久不卡| 看黄色毛片网站|