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

    正八邊形多胞薄壁管吸能特性仿真和優(yōu)化

    2015-11-30 14:56:47白中浩王飛虎郭厚銳
    關(guān)鍵詞:優(yōu)化設(shè)計(jì)

    白中浩 王飛虎+郭厚銳

    摘要:基于正八邊形多胞薄壁管相比其他截面形狀薄壁管具有結(jié)構(gòu)緊湊、比吸能高等特點(diǎn),采用有限元軟件Abaqus研究了其在準(zhǔn)靜態(tài)軸向壓縮下的吸能特性,分析了不同截面形狀和結(jié)構(gòu)參數(shù)對(duì)吸能特性的影響。根據(jù)吸能特性指標(biāo),運(yùn)用非支配遺傳算法對(duì)正八邊形多胞管邊對(duì)邊(S2S)結(jié)構(gòu)進(jìn)行了優(yōu)化。結(jié)果表明,正八邊形多胞管的邊對(duì)邊布置形式的吸能效果優(yōu)于角對(duì)角布置形式,其S2S結(jié)構(gòu)吸能特性最優(yōu)。結(jié)構(gòu)參數(shù)對(duì)其吸能特性影響明顯,經(jīng)優(yōu)化后的S2S結(jié)構(gòu)比吸能提高了33。11%,最大壓潰力降低了3。78%,具有較好的吸能特性。

    關(guān)鍵詞:多胞薄壁管;軸向壓縮;吸能特性;優(yōu)化設(shè)計(jì)

    中圖分類號(hào):U463。8 文獻(xiàn)標(biāo)識(shí)碼:A

    金屬薄壁結(jié)構(gòu)以其成本低、吸能效率高等特點(diǎn),在航空和汽車工業(yè)上的應(yīng)用越來越廣泛。在過去數(shù)十年中,國(guó)內(nèi)外針對(duì)金屬薄壁結(jié)構(gòu)進(jìn)行了大量的研究,從試驗(yàn)、理論和數(shù)值分析方面探究了不同加載條件下的金屬薄壁結(jié)構(gòu)的吸能特性。然而,過去大部分研究工作都集中在圓形、方形或矩形截面金屬薄壁結(jié)構(gòu)上\[1-4\]。

    近年來,具有多個(gè)胞元的新型薄壁結(jié)構(gòu)受到關(guān)注\[5-10\]。Chen等人\[5\]于2001年率先在理論和數(shù)值分析方面研究了軸向加載下雙胞和三胞鋁制品結(jié)構(gòu)的沖擊優(yōu)勢(shì)。Zhang等人\[6\]比較研究了多胞結(jié)構(gòu)和泡沫填充結(jié)構(gòu)的吸能效率。Najafi等人\[7\]對(duì)不同截面方形多胞結(jié)構(gòu)的吸能特性進(jìn)行了研究。亓昌等人\[8\]研究了錐形多胞薄壁管斜向沖擊下的吸能特性。上述研究都表明:多胞薄壁結(jié)構(gòu)與單胞或填充材料結(jié)構(gòu)相比,具有更高的吸能效率。然而,到目前為止,關(guān)于非直角與板組成的多胞結(jié)構(gòu)的研究較少。

    湖南大學(xué)學(xué)報(bào)(自然科學(xué)版)2015年

    第10期白中浩等:正八邊形多胞薄壁管吸能特性仿真和優(yōu)化

    Rossi等人\[11\]研究發(fā)現(xiàn),在截面面積和厚度相同的情況下,八邊形薄壁結(jié)構(gòu)比四邊形薄壁結(jié)構(gòu)的平均碰撞力增加了25%,吸能變形量降低了20%。Fan等人\[12\]研究表明,八邊形薄壁管變形均勻,吸能特性較好。張宗華等人\[13\]研究了等周長(zhǎng)的不同截面形狀多邊形薄壁結(jié)構(gòu)的軸向沖擊性能,發(fā)現(xiàn)多邊形薄壁管的平均沖擊力和吸能特性隨著多邊形頂點(diǎn)的增加而提高,但八邊形以后基本趨于穩(wěn)定。作者之前設(shè)計(jì)了一種八邊形逐級(jí)吸能結(jié)構(gòu),有效地改進(jìn)了某越野車前縱梁的碰撞性能\[14-15\],現(xiàn)針對(duì)八邊形薄壁結(jié)構(gòu)以及多胞結(jié)構(gòu)具有結(jié)構(gòu)緊湊和高效吸能的特點(diǎn),研究了一種邊對(duì)邊(two side to side,S2S)布置的正八邊形多胞管S2S結(jié)構(gòu),該結(jié)構(gòu)同時(shí)結(jié)合了八邊形管的結(jié)構(gòu)緊湊和多胞管的高比吸能特性。通過一系列仿真探討了不同截面形狀和結(jié)構(gòu)參數(shù)對(duì)吸能特性的影響,并根據(jù)吸能特性指標(biāo)對(duì)正八邊形多胞管S2S結(jié)構(gòu)進(jìn)行了優(yōu)化,優(yōu)化后其吸能特性明顯提高,為此類結(jié)構(gòu)的優(yōu)化設(shè)計(jì)和工程應(yīng)用提供參考。

    1問題描述

    1。1物理模型

    本文研究的對(duì)象為正八邊形多胞管S2S結(jié)構(gòu),如圖1所示。長(zhǎng)度L為210 mm,正八邊形截面邊長(zhǎng)d為40 mm,壁厚t都為1 mm。底端采用固定支承,上端承受剛性墻以恒定速度軸向壓縮,速度大小為V。

    圖1正八邊形多胞管軸向壓縮示意圖

    Fig。1Schematic diagram of octagon multicell

    tubes under axial compression

    1。2結(jié)構(gòu)吸能特性指標(biāo)

    一般來說,分析結(jié)構(gòu)吸能特性的關(guān)鍵參數(shù)主要有:總吸能量(E)、比吸能(SEA)、平均壓潰力(Fav)、最大壓潰力(Fpeak)等。

    比吸能SEA是評(píng)價(jià)一個(gè)結(jié)構(gòu)的能量吸收效率的指標(biāo),用來描述結(jié)構(gòu)單位質(zhì)量吸收的沖擊能量,可表示為:

    SEA=EM。 (1)

    式中:M為結(jié)構(gòu)的質(zhì)量;E為結(jié)構(gòu)吸收的總能量,其值可以通過對(duì)力位移曲線積分得到。

    E=∫d0F(x)dx。 (2)

    式中:d為結(jié)構(gòu)的有效變形量。

    平均壓潰力為Fav,其表達(dá)式如式(3)所示。由此式可以計(jì)算出平均壓潰力隨時(shí)間的變化趨勢(shì)。

    Fav=E(d)d。 (3)

    最大壓潰力是指薄壁結(jié)構(gòu)在壓潰過程中的初始?jí)簼⑤d荷峰值,即薄壁結(jié)構(gòu)與剛性墻初始接觸時(shí)峰值載荷Fpeak。

    2仿真模型

    2。1有限元建模

    采用有限元軟件Abaqus/Explicit,建立了正八邊形多胞薄壁管模型。如圖2所示,頂部剛性平面以恒定的速度5 mm/min軸向壓縮多胞管,多胞管底部與另一剛性平面焊接連接,底部剛性平面固定。多胞管采用四節(jié)點(diǎn)殼單元S4R模擬,綜合考慮精度和求解效率,殼單元沿厚度方向取5個(gè)積分點(diǎn),單元特征長(zhǎng)度取2。5 mm。多胞結(jié)構(gòu)的材料為H300Lad鋼,密度ρ=7 850 kg/m3,彈性模量E=210 GPa,泊松比ν=0。3。H300Lad鋼材的真實(shí)應(yīng)力應(yīng)變曲線如圖3所示。

    圖2正八邊形多胞管有限元模型

    Fig。2Finite element model of octagon multicell tube

    應(yīng)變?chǔ)?/p>

    圖3真實(shí)應(yīng)力應(yīng)變曲線

    Fig。3Ture strssstrain of specimens

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

    為了驗(yàn)證有限元模型,根據(jù)文獻(xiàn)\[12\]的實(shí)驗(yàn)?zāi)P蛥?shù),對(duì)正八邊形單胞管進(jìn)行準(zhǔn)靜態(tài)軸向壓潰仿真。圖4為本文有限元仿真結(jié)果與文獻(xiàn)\[12\]中的實(shí)驗(yàn)結(jié)果對(duì)比,從圖中可以看出,力響應(yīng)曲線的峰值、峰值出現(xiàn)時(shí)刻以及波動(dòng)趨勢(shì)基本一致,變形模式方面存在一定的誤差,但誤差在可接受范圍內(nèi)。因此,本文所建模型有效,可用于后續(xù)的研究。

    圖4正八邊形單胞管仿真與試驗(yàn)結(jié)果對(duì)比

    Fig。4Simulation and experimental results

    comparison of octagon singlecell tubes

    3仿真結(jié)果與分析

    3。1截面形狀對(duì)吸能特性的影響分析

    圖5為幾種幾何參數(shù)相同而截面形狀不同的正八邊形薄壁管示意圖,單胞管命名為SO(singlecell octagon);2種不同截面形狀的四胞管分別命名為S2S(two side to side)和C2C(two corner to corner);2種不同截面形狀的八胞管分別命名為S4S(four side to side)和C4C(four corner to corner)。通過對(duì)不同截面形狀的多胞管進(jìn)行仿真分析,仿真結(jié)果如表1所示。

    圖5正八邊形薄壁管不同截面形狀示意圖

    Fig。5Different cross sections shape of octagon

    thinwalled tubes

    表1多胞管數(shù)值分析結(jié)果

    Tab。1Results of simulation for multicell tubes

    樣本

    E/kJ

    Fav/kN

    Fpeak/kN

    SEA/(kJ

    Wingdings^B@ kg-1)

    SO

    6。29

    42。05

    118。10

    11。94

    C2C

    11。90

    86。41

    190。21

    14。58

    S2S

    13。44

    95。31

    191。77

    16。58

    C4C

    18。14

    144。73

    268。74

    17。46

    S4S

    22。06

    175。12

    274。24

    22。21

    從表1可以看出,平均壓潰力Fav隨著胞管數(shù)增加而增大。S4S結(jié)構(gòu)的平均壓潰力是SO結(jié)構(gòu)的4倍多。不同截面形狀的同胞數(shù)結(jié)構(gòu)的平均壓潰力也不同,如同胞數(shù)的S2S結(jié)構(gòu)與C2C結(jié)構(gòu)相比,平均壓潰力增加了10。3%。這是由于不同截面形狀導(dǎo)致出現(xiàn)不同的變形模式,當(dāng)結(jié)構(gòu)中連通邊為3的角單元在夾角小于120°時(shí),隨著夾角的增大,平均壓潰力明顯增加\[16\]。

    最大壓潰力Fpeak隨著薄壁管胞數(shù)的增加而明顯增大,S4S結(jié)構(gòu)的最大壓潰力Fpeak最大,達(dá)到了274。24 kN,是SO結(jié)構(gòu)的2倍多,比S2S結(jié)構(gòu)增加了43%。同時(shí),S2S與C2C結(jié)構(gòu)以及S4S與C4C結(jié)構(gòu)相比,最大壓潰力都略大些。其中,C4C結(jié)構(gòu)的最大壓潰力比S2S結(jié)構(gòu)增加了40。14%。

    多胞管的比吸能明顯高于單胞管,這是由于薄壁管胞數(shù)增加后,管內(nèi)部增加了交叉板,這需要更多的能量來使之變形。此外,S2S與C2C結(jié)構(gòu)以及S4S與C4C結(jié)構(gòu)相比,比吸能分別增加了13。7%和27。2%,這與結(jié)構(gòu)的截面布置形式有關(guān),同胞數(shù)管邊對(duì)邊的布置形式(side to side,見圖5(a)和圖5(c))相比角對(duì)角布置形式(corner to corner,見圖5(b)和圖5(d)),增加了拐角的數(shù)目,隨著拐角數(shù)目的增大,產(chǎn)生塑性變形需要耗散的總能量就越多。然而,C4C結(jié)構(gòu)的比吸能比S2S結(jié)構(gòu)只增加了5。3%,說明截面形式不同對(duì)多胞結(jié)構(gòu)比吸能影響較大,邊對(duì)邊的布置形式比角對(duì)角布置形式的吸能效果更好。

    3。2結(jié)構(gòu)參數(shù)對(duì)吸能特性的影響分析

    在正八邊形多胞管S2S結(jié)構(gòu)基準(zhǔn)模型的基礎(chǔ)上,研究結(jié)構(gòu)幾何參數(shù)對(duì)多胞管吸能特性的影響。依據(jù)汽車車身結(jié)構(gòu)中典型薄壁管件尺寸,各參數(shù)取值如表2所示。

    表2正八邊形多胞管幾何參數(shù)取值

    Tab。2Geometry parameters of octagon multicell tubes

    數(shù)值

    壁厚t/mm

    邊長(zhǎng)d/mm

    下限值

    0。8

    25

    上限值

    2。4

    45

    步長(zhǎng)

    0。2

    5

    3。2。1壁厚對(duì)S2S結(jié)構(gòu)的影響

    保持截面尺寸不變,首先分析壁厚對(duì)S2S結(jié)構(gòu)的影響。圖6為準(zhǔn)靜態(tài)壓縮下壁厚t對(duì)多胞管S2S結(jié)構(gòu)比吸能的影響。可以看出,隨著壁厚的增加,比吸能逐漸增大,增大幅度逐漸減小。如壁厚t從1。2 mm增加到1。8 mm,比吸能增加了39。65%;t從1。8 mm增大到2。4 mm,比吸能只增加了13。77%。

    壁厚對(duì)最大壓潰力的影響如圖7所示,最大壓潰力隨著壁厚的增加而單調(diào)上升,近似呈線性正相關(guān)。

    壁厚t/mm

    圖6壁厚對(duì)比吸能的影響

    Fig。6Effect of wall thickness on SEA

    壁厚t/mm

    圖7壁厚對(duì)最大壓潰力的影響

    Fig。7Effect of wall thickness on peak force

    3。2。2截面尺寸對(duì)S2S結(jié)構(gòu)的影響

    圖8為截面邊長(zhǎng)對(duì)S2S結(jié)構(gòu)比吸能的影響。從圖中可知,隨著截面邊長(zhǎng)的增大,比吸能單調(diào)減小,近似呈線性負(fù)相關(guān)。

    截面尺寸對(duì)最大壓潰力的影響如圖9所示,隨著截面尺寸的增大,最大壓潰力略微增加。如邊長(zhǎng)d從35 mm增大到45 mm,最大壓潰力只增長(zhǎng)了6。87%。

    邊長(zhǎng)d/mm

    圖8截面尺寸對(duì)比吸能的影響

    Fig。8Effect of section size on SEA

    邊長(zhǎng)d/mm

    圖9截面尺寸對(duì)最大壓潰力的影響

    Fig。9Effect of section size on peak force

    總的來說,各結(jié)構(gòu)參數(shù)對(duì)正八邊形多胞管S2S結(jié)構(gòu)的吸能特性影響相互關(guān)聯(lián),彼此制約,針對(duì)具體問題,有必要對(duì)S2S結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì)得到滿足需求的最優(yōu)結(jié)構(gòu)。

    4結(jié)構(gòu)優(yōu)化設(shè)計(jì)

    4。1優(yōu)化問題定義

    目前,提高薄壁結(jié)構(gòu)吸能特性普遍采用2種方法:1)采用先進(jìn)的高強(qiáng)度特性的材料;2)對(duì)壁厚、截面尺寸等參數(shù)進(jìn)行優(yōu)化。在碰撞過程中,常以最大吸能量為設(shè)計(jì)準(zhǔn)則,但是,隨著不可再生能源的日益枯竭,輕量化設(shè)計(jì)已成為當(dāng)前汽車安全性設(shè)計(jì)主要考慮的因素之一,因此比吸能更加全面地表征了結(jié)構(gòu)的吸能性。此外,碰撞過程中薄壁結(jié)構(gòu)的最大碰撞力也直接影響到車內(nèi)乘員的安全性。

    上述仿真結(jié)果表明,壁厚t和截面尺寸d對(duì)多胞管S2S結(jié)果的吸能特性影響明顯。為了進(jìn)一步提高多胞管S2S結(jié)構(gòu)的吸能特性,以多胞管S2S結(jié)構(gòu)的壁厚t和截面尺寸d為優(yōu)化參數(shù);以比吸能SEA,平均壓潰力Fav以及最大壓潰力Fpeak為優(yōu)化目標(biāo),定義多胞管S2S結(jié)構(gòu)的多目標(biāo)優(yōu)化問題為:

    maxFavSEAminFpeaks。t。0。8 mm≤t≤2。4 mm25 mm≤d≤45 mm (4)

    為了提高求解效率,優(yōu)化過程采用二次多項(xiàng)式響應(yīng)面模型代替有限元仿真分析,并通過優(yōu)拉丁超立方試驗(yàn)設(shè)計(jì)(DOE)方法設(shè)計(jì)9組仿真試驗(yàn),建立響應(yīng)面模型。利用非支配遺傳算法(NSGAⅡ)\[17\]尋求多目標(biāo)優(yōu)化問題(式(4))的帕雷托最優(yōu)解,優(yōu)化流程如圖10所示。

    4。2構(gòu)建近似模型

    近似模型的方法已經(jīng)被證明是一種特別有效的方法,被廣泛應(yīng)用于薄壁管的優(yōu)化設(shè)計(jì)中,如響應(yīng)面模型(RSM)\[18\]。本文基于優(yōu)拉丁超立方試驗(yàn)設(shè)計(jì)的試驗(yàn)結(jié)果,利用逐步回歸方法得到優(yōu)化目標(biāo)關(guān)于兩個(gè)設(shè)計(jì)變量的二階多項(xiàng)式響應(yīng)面模型的多元回歸方程為:

    Fav=-121。046+73。687t+4。782d-

    3。579t2+3。786td-0。12d2; (5)

    SEA=23。04+23。466t-0。866d-

    4。235t2+0。011td+0。005d2;(6)

    Fpeak=44。371-18。104t-2。013d+

    23。895t2-0。048td+7。197d2。 (7)

    式中:Fav,SEA和Fpeak分別為正八邊形多胞管S2S結(jié)構(gòu)的平均壓潰力、比吸能和最大壓潰力;t和d分別為多胞管的壁厚和截面邊長(zhǎng)。

    圖10優(yōu)化流程圖

    Fig。10Optimization process sketch map

    平均壓潰力Fav的回歸方程的決定系數(shù)R2為99。97%,比吸能SEA的回歸方程的決定系數(shù)R2為99。93%,最大壓潰力Fpeak的回歸方程決定系數(shù)R2為99。99%。由R2值可以看出,這3個(gè)二階多項(xiàng)式響應(yīng)面模型的擬合精度較好,能夠較好地滿足預(yù)測(cè)精度的要求。

    表3給出了5個(gè)隨機(jī)樣本點(diǎn)的設(shè)計(jì)目標(biāo)響應(yīng)面模型計(jì)算值和有限元分析值。從表中可以看出,響應(yīng)面模型的計(jì)算值與有限元的分析值都相當(dāng)接近,相對(duì)誤差不超過5%,因此,設(shè)計(jì)目標(biāo)的響應(yīng)面模型精度都很高,可以替代有限元模型用于后續(xù)的多目標(biāo)優(yōu)化設(shè)計(jì)過程。

    4。3優(yōu)化結(jié)果分析

    多目標(biāo)優(yōu)化問題最明顯的特征就是不存在唯一的全局最優(yōu)解,而是存在一個(gè)最優(yōu)解集,稱為Pareto解集或Pareto前沿?;诖私饧?,設(shè)計(jì)者可以根據(jù)具體問題確定最終優(yōu)化方案。

    本文基于響應(yīng)面模型,采用非支配遺傳算法對(duì)回歸方程進(jìn)行多目標(biāo)優(yōu)化求解得到的Pareto前沿如圖11所示。圖中三角形點(diǎn)、

    圓形點(diǎn)和五角星點(diǎn)所示分別為最大壓潰力、比吸能和平均壓潰力的最優(yōu)解。由于3個(gè)目標(biāo)之間存在相互制約的矛盾,在此以比吸能為主要目標(biāo),兼顧最大壓潰力和平均壓潰力,選取權(quán)重分別為0。5,0。3,0。2。本多目標(biāo)優(yōu)化問題的最優(yōu)解如圖11中五邊形點(diǎn)所示,此時(shí)設(shè)計(jì)點(diǎn)取整后d=30 mm,t=1。0 mm。

    5結(jié)論

    研究了一種用于碰撞吸能的正八邊形多胞薄壁管S2S結(jié)構(gòu),采用有限元仿真分析了截面形狀(C2C,S2S,C4C和S4S)和結(jié)構(gòu)參數(shù)對(duì)吸能特性的影響,并對(duì)多胞管S2S結(jié)構(gòu)的結(jié)構(gòu)參數(shù)進(jìn)行了優(yōu)化。研究結(jié)果表明:

    1)與單胞管、多胞管C2C,C4C和S4S結(jié)構(gòu)相比,多胞管S2S結(jié)構(gòu)的吸能特性更好;邊對(duì)邊布置形式的吸能效果優(yōu)于角對(duì)角布置形式。

    2)減小截面尺寸、增加壁厚可以提高多胞管的比吸能,但壁厚的增加使得最大壓潰力也隨之增大;最大壓潰力與壁厚呈近似線性正相關(guān),可以通過減小壁厚和增大截面尺寸降低最大壓潰力。

    3)結(jié)合優(yōu)拉丁超立方試驗(yàn)設(shè)計(jì)方法和非支配遺傳算法,對(duì)多胞管S2S結(jié)構(gòu)的吸能特性進(jìn)行了優(yōu)化設(shè)計(jì),優(yōu)化結(jié)果比初始設(shè)計(jì)的吸能特性有了較大提高,比吸能比優(yōu)化前提高了33。11%,最大壓潰力降低了3。78%,吸能效果更好。

    參考文獻(xiàn)

    [1]ABRAMOWICZ W, JONES N。 Dynamic progressive buckling of circular and square tubes\[J\]。 International Journal of Impact Engineering, 1986, 4(4):243-270。

    \[2\]HUANG X, LU G, YU T, et al。 On the axial splitting and curling of circular metal tubes\[J\]。 International Journal of Mechanical Sciences,2002, 44(11):2369-2391。

    \[3\]TARIGOPULA V, LANGSETH M, HOPPERSTAD O S, et al。Axial crushing of thinwalled highstrength steel sections\[J\]。 International Journal of Impact Engineering, 2006, 32(5):847-882。

    \[4\] 顧紅軍,趙國(guó)志,陸廷金,等。軸向沖擊下薄壁圓柱殼的屈曲行為的實(shí)驗(yàn)研究\[J\]。振動(dòng)與沖擊,2004,23(4):58-59。

    GU Hongjun, ZHAO Guozhi, LU Tingjin, et al。 Buckling of thinwall cylindrical shell under axial impact\[J\]。 Journal of Vibration and Shock, 2004, 23(4):58-59。 (In Chinese)

    \[5\]CHEN W, WIERZBICKI T。 Relative merits of singlecell, multicell and foamfilled thinwalled structures in energy absorption\[J\]。ThinWalled Structures, 2001, 39(4): 287-306。

    \[6\]ZHANG X, CHENG G。 A comparative study of energy absorption characteristics of foamfilled and multicell square columns\[J\]。International Journal of Impact Engineering, 2007, 34(11):1739-1752。

    \[7\]NAJAFI A, RAISROHANI M。 Mechanics of axial plastic collapse in multicell, multicorner crush tubes\[J\]。 ThinWalled Structures, 2011, 49(1):1-12。

    \[8\] 亓昌,董方亮,楊姝,等。錐形多胞薄壁管斜向沖擊吸能特性仿真研究\[J\]。振動(dòng)與沖擊,2012,31(24):102-107。

    QI Chang, DONG Fangliang, YANG Shu。 Simulation and optimization for blastresistant performances of a graded aluminum foam sandwich structure\[J\]。 Journal of Vibration and Shock, 2012, 31(24):102-107。 (In Chinese)

    \[9\]BAI Z, GUO H, JIANG B, et al。 A study on the mean crushing strength of hexagonal multicell thinwalled structures\[J\]。 ThinWalled Structures, 2014,80: 38-45。

    \[10\]JIANG B, ZHU F, JIN X, et al。 Computational modeling of the crushing behavior of SKYDEX material using homogenized material laws\[J\]。Composite Structures, 2013,106:306-316。

    \[11\]ROSSI A, FAWAZ Z, BEHDINAN K。 Numerical simulation of the axial collapse of thinwalled polygonal section tubes\[J\]。 ThinWalled Structures, 2005, 43(10):1646-1661。

    \[12\]FAN Z, LU G, LIU K。 Quasistatic axial compression of thinwalled tubes with different crosssectional shapes\[J\]。 Engineering Structures, 2013, 55:80-89。

    \[13\]張宗華,劉樹田。多邊形薄壁管動(dòng)態(tài)軸向沖擊的耐撞性研究\[C\]//2007年中國(guó)汽車工程學(xué)會(huì)年會(huì)論文集。北京:機(jī)械工業(yè)出版社,2007:437-443。

    ZHANG Zonghua, LIU Shutian。 Crashworthiness of dynamical axial crushing of polygonal thinwalled tubes\[C\]//2007SAEChina Congress Proceedings。 Beijing: China Machine Press, 2007:437-443。 (In Chinese)

    \[14\]白中浩,陳天志,曹立波,等?;谡辉O(shè)計(jì)的汽車前縱梁吸能結(jié)構(gòu)的優(yōu)化\[J\]。汽車工程,2010,32(11):935-939。

    BAI Zhonghao, CHEN Tianzhi, CAO Libo, et al。 Optimization of energyabsorbing structure of vehicle front side rail based on orthogonal design\[J\]。 Automotive Engineering, 2010,32(11):935-939。(In Chinese)

    \[15\]曹立波,崔崇禎,白中浩,等。八邊形逐級(jí)吸能梁的設(shè)計(jì)與優(yōu)化\[J\]。湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2010,37(1):29-34。

    CAO Libo, CUI Chongzhen, BAI Zhonghao, et al。 Design and optimization of the octagonal step by step energy\[J\]。 Journal of Hunan University: Natural Sciences, 2010, 37(1):29-34。(In Chinese)

    \[16\]JONES N。 Energyabsorbing effectiveness factor\[J\]。 International Journal of Impact Engineering, 2010, 37(6):754-765。

    \[17\]白中浩,盧靜?;诙嗄繕?biāo)優(yōu)化的兒童座椅匹配方法\[J\]。湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2013,40(3):52-58。

    BAI Zhonghao, LU Jing。 Study of matching for a child seat based on multiobjective optimization\[J\]。 Journal of Hunan University: Natural Sciences, 2013, 40(3):52-58。(In Chinese)

    \[18\]YIN H, WEN G, HOU S, et al。 Crushing analysis and multiobjective crashworthiness optimization of honeycombfilled single and bitubular polygonal tubes\[J\]。 Materials and Design, 2011, 32(8/9):4449-4460。

    猜你喜歡
    優(yōu)化設(shè)計(jì)
    導(dǎo)彈舵面的復(fù)合材料設(shè)計(jì)與分析
    航空兵器(2016年4期)2016-11-28 21:47:29
    礦井主排水系統(tǒng)的優(yōu)化設(shè)計(jì)與改造
    科技資訊(2016年19期)2016-11-15 08:34:13
    數(shù)據(jù)挖掘?qū)虒W(xué)管理的優(yōu)化設(shè)計(jì)
    如何實(shí)現(xiàn)小學(xué)數(shù)學(xué)課堂練習(xí)設(shè)計(jì)優(yōu)化
    淺析人機(jī)工程學(xué)在家具創(chuàng)作中的作用
    試析機(jī)械結(jié)構(gòu)優(yōu)化設(shè)計(jì)的應(yīng)用及趨勢(shì)
    汽車行李箱蓋鉸鏈機(jī)構(gòu)的分析及優(yōu)化
    東林煤礦保護(hù)層開采卸壓瓦斯抽采優(yōu)化設(shè)計(jì)
    橋式起重機(jī)主梁結(jié)構(gòu)分析和優(yōu)化設(shè)計(jì)
    對(duì)無線傳感器網(wǎng)絡(luò)MAC層協(xié)議優(yōu)化的研究與設(shè)計(jì)
    科技視界(2016年22期)2016-10-18 15:25:08
    国产野战对白在线观看| 美国免费a级毛片| 麻豆乱淫一区二区| 男人添女人高潮全过程视频| 国产成人精品久久久久久| 国产国语露脸激情在线看| 一级毛片电影观看| 久久人人爽人人片av| 国产又爽黄色视频| 免费人成视频x8x8入口观看| 精品一区二区三区四区五区乱码| 亚洲国产欧美网| 每晚都被弄得嗷嗷叫到高潮| 国产黄a三级三级三级人| 国产精品久久视频播放| 精品久久久久久久久久久久久 | 身体一侧抽搐| 日本精品一区二区三区蜜桃| 欧美人与性动交α欧美精品济南到| 日韩欧美国产在线观看| 麻豆成人av在线观看| 欧美在线黄色| 母亲3免费完整高清在线观看| 成年女人毛片免费观看观看9| 91在线观看av| 精品人妻1区二区| 变态另类丝袜制服| 热99re8久久精品国产| 中文字幕久久专区| 国内少妇人妻偷人精品xxx网站 | 久久久久久久久久黄片| 男人操女人黄网站| 中文在线观看免费www的网站 | 国产激情欧美一区二区| 黄色毛片三级朝国网站| 精品久久久久久久人妻蜜臀av| 亚洲国产日韩欧美精品在线观看 | 麻豆成人av在线观看| 999久久久国产精品视频| 男女午夜视频在线观看| 精华霜和精华液先用哪个| av在线播放免费不卡| 黑丝袜美女国产一区| 在线观看日韩欧美| 看黄色毛片网站| 18禁观看日本| av在线天堂中文字幕| 精品无人区乱码1区二区| 久久国产亚洲av麻豆专区| 一进一出抽搐gif免费好疼| 亚洲男人天堂网一区| 18禁裸乳无遮挡免费网站照片 | av中文乱码字幕在线| 国产亚洲欧美在线一区二区| 亚洲精华国产精华精| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 嫩草影视91久久| 女性生殖器流出的白浆| 男女床上黄色一级片免费看| 久久午夜亚洲精品久久| 啦啦啦韩国在线观看视频| 欧美亚洲日本最大视频资源| 亚洲成av片中文字幕在线观看| 国产色视频综合| www.熟女人妻精品国产| 麻豆成人av在线观看| 国产91精品成人一区二区三区| 免费搜索国产男女视频| 99热6这里只有精品| 久99久视频精品免费| 欧美一级毛片孕妇| 一级毛片女人18水好多| 黑人操中国人逼视频| 精品国产国语对白av| 国产精品乱码一区二三区的特点| 国产精品国产高清国产av| 999久久久精品免费观看国产| 国产私拍福利视频在线观看| 国产精品美女特级片免费视频播放器 | 在线永久观看黄色视频| 久久久久九九精品影院| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看| 亚洲av第一区精品v没综合| 欧美另类亚洲清纯唯美| 亚洲精品在线美女| 色综合站精品国产| 色综合欧美亚洲国产小说| 亚洲国产精品合色在线| 黄片小视频在线播放| 国产亚洲av嫩草精品影院| 男女那种视频在线观看| 人人妻,人人澡人人爽秒播| av在线天堂中文字幕| 亚洲第一欧美日韩一区二区三区| 又紧又爽又黄一区二区| 99在线人妻在线中文字幕| 久99久视频精品免费| 国产亚洲精品综合一区在线观看 | 免费观看人在逋| 精品久久蜜臀av无| 亚洲国产高清在线一区二区三 | 国产精品永久免费网站| 精品电影一区二区在线| 久久人妻福利社区极品人妻图片| svipshipincom国产片| 中文亚洲av片在线观看爽| avwww免费| 久久久水蜜桃国产精品网| 欧美日韩乱码在线| 女性生殖器流出的白浆| 国产精品精品国产色婷婷| 久久久精品欧美日韩精品| 18禁美女被吸乳视频| 国产成人av教育| 午夜成年电影在线免费观看| 一区二区三区国产精品乱码| 午夜两性在线视频| 88av欧美| 国产精品亚洲av一区麻豆| 欧美不卡视频在线免费观看 | 1024视频免费在线观看| 国产激情欧美一区二区| av中文乱码字幕在线| 满18在线观看网站| 最新美女视频免费是黄的| 亚洲熟女毛片儿| 长腿黑丝高跟| 欧美黄色淫秽网站| 亚洲五月色婷婷综合| 长腿黑丝高跟| av片东京热男人的天堂| 欧美日韩瑟瑟在线播放| 国产精品一区二区三区四区久久 | 精品一区二区三区视频在线观看免费| 一个人免费在线观看的高清视频| 久久久久国产一级毛片高清牌| 亚洲精品一卡2卡三卡4卡5卡| 深夜精品福利| 狂野欧美激情性xxxx| 精品福利观看| 色综合站精品国产| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区视频了| 成人亚洲精品一区在线观看| 日本在线视频免费播放| 国产精品影院久久| 又大又爽又粗| 波多野结衣巨乳人妻| 美女大奶头视频| 最近最新免费中文字幕在线| 热re99久久国产66热| 国产成人系列免费观看| 久久香蕉激情| 国产精品亚洲美女久久久| 黄色女人牲交| 亚洲欧美激情综合另类| 欧美精品亚洲一区二区| 国产一区二区激情短视频| 欧美中文日本在线观看视频| 亚洲成人免费电影在线观看| 叶爱在线成人免费视频播放| 国产在线精品亚洲第一网站| 亚洲国产精品999在线| 久久精品成人免费网站| 女人爽到高潮嗷嗷叫在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 首页视频小说图片口味搜索| 欧美亚洲日本最大视频资源| 夜夜爽天天搞| 日日爽夜夜爽网站| 在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 欧美乱码精品一区二区三区| 亚洲成av片中文字幕在线观看| 99久久精品国产亚洲精品| 日本熟妇午夜| 女生性感内裤真人,穿戴方法视频| 青草久久国产| 久久久久久久久免费视频了| 国产一区在线观看成人免费| 夜夜爽天天搞| 观看免费一级毛片| 99re在线观看精品视频| 久久香蕉激情| 国语自产精品视频在线第100页| 亚洲精品美女久久av网站| 性欧美人与动物交配| 国产亚洲av高清不卡| 视频区欧美日本亚洲| 国产成人啪精品午夜网站| 妹子高潮喷水视频| 日韩欧美 国产精品| 国产又色又爽无遮挡免费看| 欧美日韩乱码在线| 国产区一区二久久| 久久精品亚洲精品国产色婷小说| 欧美一区二区精品小视频在线| 男人舔女人下体高潮全视频| 两性午夜刺激爽爽歪歪视频在线观看 | 母亲3免费完整高清在线观看| 久久香蕉激情| 老司机福利观看| 免费在线观看亚洲国产| 国产激情偷乱视频一区二区| 一级黄色大片毛片| 亚洲自拍偷在线| 午夜免费观看网址| 久久精品91无色码中文字幕| 女人被狂操c到高潮| 波多野结衣巨乳人妻| 成人av一区二区三区在线看| 国产高清激情床上av| 免费人成视频x8x8入口观看| 精品欧美国产一区二区三| 色播亚洲综合网| 亚洲avbb在线观看| 国产视频一区二区在线看| 女生性感内裤真人,穿戴方法视频| 一进一出抽搐动态| 在线观看舔阴道视频| 亚洲欧美日韩无卡精品| 90打野战视频偷拍视频| 人人妻人人澡人人看| 国产成人精品久久二区二区免费| 婷婷精品国产亚洲av在线| 成在线人永久免费视频| e午夜精品久久久久久久| 色婷婷久久久亚洲欧美| 欧美午夜高清在线| 久9热在线精品视频| 国产亚洲精品av在线| 久久狼人影院| 嫁个100分男人电影在线观看| 看黄色毛片网站| av欧美777| 18美女黄网站色大片免费观看| 伦理电影免费视频| 国产精品1区2区在线观看.| 国产成人系列免费观看| 男女视频在线观看网站免费 | 午夜福利视频1000在线观看| 免费一级毛片在线播放高清视频| 国产一区二区三区视频了| 长腿黑丝高跟| 后天国语完整版免费观看| 国内久久婷婷六月综合欲色啪| 女人爽到高潮嗷嗷叫在线视频| 亚洲av成人不卡在线观看播放网| 美女高潮到喷水免费观看| 最近最新免费中文字幕在线| 成人欧美大片| 免费看a级黄色片| 大型av网站在线播放| 黄色毛片三级朝国网站| 国产精品一区二区三区四区久久 | 男人舔女人的私密视频| 国产精品 国内视频| 精品卡一卡二卡四卡免费| 亚洲男人天堂网一区| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 在线看三级毛片| 老汉色∧v一级毛片| 国产精品久久久久久人妻精品电影| 国产精品野战在线观看| 国产一卡二卡三卡精品| 亚洲国产精品久久男人天堂| 国产真人三级小视频在线观看| 国产国语露脸激情在线看| 午夜两性在线视频| 久久精品成人免费网站| 亚洲欧美精品综合一区二区三区| 性色av乱码一区二区三区2| 日日摸夜夜添夜夜添小说| 久久天堂一区二区三区四区| 国产久久久一区二区三区| 村上凉子中文字幕在线| av片东京热男人的天堂| av天堂在线播放| 给我免费播放毛片高清在线观看| 一级黄色大片毛片| av视频在线观看入口| 91在线观看av| 19禁男女啪啪无遮挡网站| www.熟女人妻精品国产| 欧美黑人精品巨大| 校园春色视频在线观看| 精华霜和精华液先用哪个| 老司机午夜十八禁免费视频| 韩国精品一区二区三区| 一进一出抽搐动态| 大型黄色视频在线免费观看| 午夜精品久久久久久毛片777| 国产视频一区二区在线看| 大型av网站在线播放| 久久久久免费精品人妻一区二区 | 中文在线观看免费www的网站 | 亚洲精华国产精华精| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产精品麻豆| 国内精品久久久久精免费| 欧美午夜高清在线| 日韩欧美国产在线观看| 欧美午夜高清在线| 中文在线观看免费www的网站 | 久久精品91无色码中文字幕| 国产精品免费视频内射| 亚洲国产欧美一区二区综合| 日本熟妇午夜| 最近最新中文字幕大全免费视频| 色婷婷久久久亚洲欧美| 成人午夜高清在线视频 | 久久香蕉国产精品| 色综合欧美亚洲国产小说| 激情在线观看视频在线高清| 十八禁网站免费在线| 在线观看日韩欧美| 久久精品亚洲精品国产色婷小说| 亚洲美女黄片视频| 叶爱在线成人免费视频播放| 国产久久久一区二区三区| 国产99白浆流出| 国产一区二区在线av高清观看| 国产极品粉嫩免费观看在线| 一本大道久久a久久精品| 国产精品久久久人人做人人爽| 真人一进一出gif抽搐免费| 波多野结衣巨乳人妻| 久久国产精品男人的天堂亚洲| 大香蕉久久成人网| 国产99白浆流出| www.熟女人妻精品国产| 99久久久亚洲精品蜜臀av| 99久久无色码亚洲精品果冻| 国产黄片美女视频| 亚洲黑人精品在线| 亚洲精品一卡2卡三卡4卡5卡| 国产成人精品无人区| 老司机深夜福利视频在线观看| 黄色视频,在线免费观看| 一区二区三区国产精品乱码| 久久久精品欧美日韩精品| 一区二区三区国产精品乱码| 嫩草影院精品99| 国产成人精品久久二区二区免费| 亚洲av电影在线进入| 一级毛片高清免费大全| 免费高清在线观看日韩| 91九色精品人成在线观看| 国产精品爽爽va在线观看网站 | 亚洲成人国产一区在线观看| 久久国产精品人妻蜜桃| 亚洲男人天堂网一区| 观看免费一级毛片| 十八禁网站免费在线| 精品一区二区三区av网在线观看| 亚洲精品色激情综合| 中文在线观看免费www的网站 | 一区二区三区精品91| 精品少妇一区二区三区视频日本电影| 十八禁人妻一区二区| 国产黄片美女视频| 久久天躁狠狠躁夜夜2o2o| 白带黄色成豆腐渣| 一边摸一边抽搐一进一小说| 一区二区三区激情视频| 亚洲成av人片免费观看| 丝袜在线中文字幕| 国产成人av教育| 97碰自拍视频| 亚洲一区高清亚洲精品| 精华霜和精华液先用哪个| 亚洲欧美精品综合久久99| 国产在线观看jvid| 亚洲av电影不卡..在线观看| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 亚洲,欧美精品.| 国产精品一区二区免费欧美| АⅤ资源中文在线天堂| 人人妻人人看人人澡| videosex国产| 女人爽到高潮嗷嗷叫在线视频| 精品国产乱码久久久久久男人| 久久精品国产综合久久久| 久久久久久国产a免费观看| 成人精品一区二区免费| 极品教师在线免费播放| 欧美精品啪啪一区二区三区| 国产精品久久视频播放| 国产99白浆流出| 1024手机看黄色片| 日本免费一区二区三区高清不卡| x7x7x7水蜜桃| 亚洲欧美日韩无卡精品| 国产成人啪精品午夜网站| 久久精品国产亚洲av香蕉五月| 成人18禁在线播放| 99精品欧美一区二区三区四区| 怎么达到女性高潮| 久久久久久大精品| 久久精品亚洲精品国产色婷小说| 免费人成视频x8x8入口观看| 中文亚洲av片在线观看爽| 两性夫妻黄色片| 亚洲人成网站在线播放欧美日韩| 久久午夜亚洲精品久久| 日韩精品青青久久久久久| xxxwww97欧美| 巨乳人妻的诱惑在线观看| 欧美另类亚洲清纯唯美| 国产精品,欧美在线| 精品久久久久久久久久免费视频| 19禁男女啪啪无遮挡网站| 热re99久久国产66热| 久久中文字幕人妻熟女| 国产精品国产高清国产av| 一个人观看的视频www高清免费观看 | 又大又爽又粗| 丝袜在线中文字幕| 亚洲成a人片在线一区二区| 午夜福利18| 一卡2卡三卡四卡精品乱码亚洲| 一a级毛片在线观看| 成人18禁在线播放| 久久精品影院6| 中文字幕高清在线视频| 黄频高清免费视频| 婷婷精品国产亚洲av| 丝袜人妻中文字幕| 国产爱豆传媒在线观看 | 久久久久久久久中文| 久久人妻福利社区极品人妻图片| 国产黄片美女视频| 久久九九热精品免费| 久久精品国产清高在天天线| 国产片内射在线| 波多野结衣高清作品| 国产1区2区3区精品| 丰满的人妻完整版| 窝窝影院91人妻| 99国产综合亚洲精品| 亚洲一区二区三区不卡视频| 亚洲精华国产精华精| 动漫黄色视频在线观看| 国产成人av教育| 亚洲国产精品999在线| 亚洲av电影不卡..在线观看| 亚洲av美国av| 色综合婷婷激情| 欧美激情高清一区二区三区| 欧美精品亚洲一区二区| 麻豆成人午夜福利视频| 久久久国产欧美日韩av| 久久久久精品国产欧美久久久| 亚洲av成人av| 国产成+人综合+亚洲专区| 免费搜索国产男女视频| 一区二区三区激情视频| 中文字幕精品亚洲无线码一区 | 又大又爽又粗| 亚洲 国产 在线| 亚洲一区高清亚洲精品| 久久精品夜夜夜夜夜久久蜜豆 | 真人做人爱边吃奶动态| 两个人视频免费观看高清| 亚洲专区国产一区二区| 国产极品粉嫩免费观看在线| 波多野结衣av一区二区av| 亚洲成国产人片在线观看| 久久精品夜夜夜夜夜久久蜜豆 | av电影中文网址| 亚洲精品中文字幕一二三四区| 久热爱精品视频在线9| 1024手机看黄色片| 亚洲国产毛片av蜜桃av| 欧美大码av| 欧美性猛交黑人性爽| 欧美日本视频| 麻豆av在线久日| 亚洲一区二区三区色噜噜| 热99re8久久精品国产| 最好的美女福利视频网| 日韩欧美在线二视频| 免费高清视频大片| 欧美绝顶高潮抽搐喷水| 日韩av在线大香蕉| 国产99白浆流出| av片东京热男人的天堂| 精品国产美女av久久久久小说| 亚洲精品av麻豆狂野| 免费女性裸体啪啪无遮挡网站| 久热爱精品视频在线9| 欧美日韩福利视频一区二区| 亚洲国产精品999在线| 久久午夜综合久久蜜桃| 美女高潮到喷水免费观看| 中文字幕精品免费在线观看视频| 免费人成视频x8x8入口观看| 少妇粗大呻吟视频| 人妻久久中文字幕网| 美女高潮喷水抽搐中文字幕| 久久精品国产综合久久久| 最近最新免费中文字幕在线| 亚洲全国av大片| 亚洲欧美精品综合久久99| 久久精品国产清高在天天线| 精品久久久久久久人妻蜜臀av| 满18在线观看网站| 丝袜在线中文字幕| 特大巨黑吊av在线直播 | 丁香六月欧美| 精品免费久久久久久久清纯| 精品国产超薄肉色丝袜足j| 亚洲人成伊人成综合网2020| 在线永久观看黄色视频| 热re99久久国产66热| 欧美黑人巨大hd| 999久久久国产精品视频| 校园春色视频在线观看| 日韩三级视频一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲av嫩草精品影院| 老熟妇仑乱视频hdxx| 欧美一级毛片孕妇| 一二三四社区在线视频社区8| 女人爽到高潮嗷嗷叫在线视频| 日韩高清综合在线| 2021天堂中文幕一二区在线观 | 深夜精品福利| 日韩国内少妇激情av| 香蕉国产在线看| 国产在线观看jvid| 免费在线观看日本一区| 亚洲精品中文字幕一二三四区| 黄网站色视频无遮挡免费观看| 夜夜夜夜夜久久久久| 一进一出好大好爽视频| 搡老熟女国产l中国老女人| www.www免费av| 免费看a级黄色片| 最新美女视频免费是黄的| 嫩草影院精品99| 亚洲人成网站在线播放欧美日韩| e午夜精品久久久久久久| 免费电影在线观看免费观看| 久久精品亚洲精品国产色婷小说| 欧美成狂野欧美在线观看| 成人三级黄色视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美一区二区三区黑人| 少妇粗大呻吟视频| 日韩欧美在线二视频| 成人亚洲精品一区在线观看| 在线观看免费视频日本深夜| 国产日本99.免费观看| 国产主播在线观看一区二区| 亚洲av片天天在线观看| 嫁个100分男人电影在线观看| 97碰自拍视频| 最新美女视频免费是黄的| 日本免费a在线| 亚洲精品在线美女| 亚洲第一欧美日韩一区二区三区| 哪里可以看免费的av片| 午夜激情av网站| 很黄的视频免费| 99热只有精品国产| 超碰成人久久| 可以在线观看的亚洲视频| 99re在线观看精品视频| 亚洲 欧美一区二区三区| cao死你这个sao货| 国产精品影院久久| 日本一本二区三区精品| 狠狠狠狠99中文字幕| 天天躁夜夜躁狠狠躁躁| 午夜久久久在线观看| 久久久国产精品麻豆| 中文亚洲av片在线观看爽| 狂野欧美激情性xxxx| 亚洲精品粉嫩美女一区| 变态另类成人亚洲欧美熟女| 97超级碰碰碰精品色视频在线观看| 18禁国产床啪视频网站| 欧美 亚洲 国产 日韩一| 精品熟女少妇八av免费久了| 午夜福利一区二区在线看| 亚洲国产日韩欧美精品在线观看 | 久久国产乱子伦精品免费另类| 人人澡人人妻人| 国产精品久久久久久人妻精品电影| 国产精品久久久久久亚洲av鲁大| 久久国产精品人妻蜜桃| 大型av网站在线播放| 久久久久亚洲av毛片大全| 亚洲人成网站高清观看| 人成视频在线观看免费观看| 精品人妻1区二区| 精品国产乱子伦一区二区三区| 夜夜爽天天搞| 久久人人精品亚洲av| 国产激情欧美一区二区| 这个男人来自地球电影免费观看| 成人三级黄色视频| 精品久久久久久久人妻蜜臀av| 一区二区三区精品91| avwww免费| 中文资源天堂在线| 禁无遮挡网站| 午夜福利欧美成人| 免费在线观看日本一区| 久久精品成人免费网站|