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

    軸向壓縮下泡沫鋁填充薄壁圓管吸能特性研究

    2016-04-15 09:52:28劉志芳路國(guó)運(yùn)閆慶榮
    關(guān)鍵詞:數(shù)值模擬

    黃 睿,劉志芳,路國(guó)運(yùn),閆慶榮

    (太原理工大學(xué) 應(yīng)用力學(xué)與生物醫(yī)學(xué)工程研究所,太原 030024)

    ?

    軸向壓縮下泡沫鋁填充薄壁圓管吸能特性研究

    黃睿,劉志芳,路國(guó)運(yùn),閆慶榮

    (太原理工大學(xué) 應(yīng)用力學(xué)與生物醫(yī)學(xué)工程研究所,太原 030024)

    摘要:利用有限元軟件ABAQUS對(duì)7種不同幾何尺寸的泡沫鋁填充管進(jìn)行準(zhǔn)靜態(tài)軸向壓縮的數(shù)值仿真分析,系統(tǒng)地研究了管的高度、壁厚、直徑以及泡沫鋁的填充對(duì)圓管吸能性能的影響;與實(shí)驗(yàn)對(duì)比,分析了7種泡沫鋁填充試件的平均載荷、初始峰值載荷、比吸能(ESA)和壓縮力效率(ECF)等吸能評(píng)價(jià)指標(biāo)。研究結(jié)果表明,泡沫鋁填充管在準(zhǔn)靜態(tài)軸向壓縮時(shí),管的壁厚與直徑對(duì)管的吸能性能有顯著影響;管的高度對(duì)其吸能性能影響較小,但高度的增加可以增加管的總吸能;泡沫鋁的填充提高了管的承載力、總吸能、比吸能和壓縮力效率。本研究成果可為設(shè)計(jì)理想的緩沖吸能裝置提供一定的技術(shù)依據(jù)。

    關(guān)鍵詞:泡沫鋁填充管;軸向壓縮;能量吸收;數(shù)值模擬

    薄壁金屬管是工程結(jié)構(gòu)中常用的一種構(gòu)件,也是一種理想的緩沖吸能元件。因其結(jié)構(gòu)簡(jiǎn)單、易于加工、具有良好的能量吸收特性,而廣泛應(yīng)用于航空航天、交通運(yùn)輸、石油化工等領(lǐng)域。但是,薄壁金屬單管結(jié)構(gòu)的承載力較弱,且變形不穩(wěn)定,容易發(fā)生歐拉失穩(wěn)。因此,為了提高管件的吸能效率,泡沫填充管作為一種新型的緩沖吸能元件被廣泛地應(yīng)用于飛行器設(shè)計(jì)、汽車(chē)耐撞性等領(lǐng)域。泡沫金屬是一種輕質(zhì)高效的材料,它有著幾乎不變的平應(yīng)力,而且壓縮行程較長(zhǎng),具有較高的吸能效率,將其填充到薄壁金屬管中,可構(gòu)成一種性能良好的抗沖擊吸能裝置,具有廣闊的應(yīng)用前景。因而對(duì)泡沫金屬填充薄壁管狀結(jié)構(gòu)的研究具有十分重要的意義。

    近年來(lái),國(guó)內(nèi)外學(xué)者對(duì)于薄壁金屬管及其泡沫填充結(jié)構(gòu)的吸能特性做了大量的研究,并取得了一些研究成果。早期ANDREWS,et al[1],GUILLOW,et al[2]通過(guò)大量的實(shí)驗(yàn),觀察薄壁金屬圓管的軸向壓潰變形,發(fā)現(xiàn)變形模式的不同主要與管件的長(zhǎng)度與直徑比值L/D以及直徑與厚度的比值D/d有關(guān);CHEN,et al[3]通過(guò)有限元模擬法預(yù)測(cè)了軸向沖擊下圓管的初始峰值載荷,討論了影響初始峰值載荷大小的因素,如管的幾何尺寸、材料的屬性和沖擊速度等;KARAGIOZOVA,et al[4]研究了圓管結(jié)構(gòu)的幾何尺寸、材料參數(shù)以及邊界條件對(duì)結(jié)構(gòu)的屈曲變形模式和吸能的影響;SEITZBERGER,et al[5]對(duì)泡沫填充管的軸向壓縮特性分別進(jìn)行了實(shí)驗(yàn)和數(shù)值模擬研究,發(fā)現(xiàn)填充泡沫對(duì)圓管吸能性能的提高有很大的作用;桂良進(jìn)等[6-7]對(duì)軸向壓縮及軸向沖擊下泡沫填充管的吸能特性進(jìn)行了研究,系統(tǒng)地分析了結(jié)構(gòu)的幾何尺寸、材料參數(shù)、加載速率對(duì)吸能的影響,并對(duì)結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì);SANTOSA,et al[8]通過(guò)實(shí)驗(yàn)和數(shù)值模擬的方法對(duì)軸向壓縮下的泡沫鋁填充管的力學(xué)響應(yīng)進(jìn)行了研究,在模擬中用動(dòng)態(tài)顯式分析法模擬準(zhǔn)靜態(tài)實(shí)驗(yàn),所預(yù)測(cè)出的變形模式以及壓潰力,與實(shí)驗(yàn)得到的結(jié)果吻合較好?;谠诟咚倭熊?chē)上的應(yīng)用,張志新[9]設(shè)計(jì)了管長(zhǎng)度不同的泡沫鋁填充裝置,并將它帶入到列車(chē)模型中,通過(guò)數(shù)值模擬分析發(fā)現(xiàn),該裝置對(duì)保護(hù)乘客安全起到很大的作用。

    綜上所述,前人的研究都只單獨(dú)考慮幾何尺寸或是填充泡沫對(duì)薄壁金屬管吸能性能的影響。在此基礎(chǔ)上,筆者同時(shí)考慮管的幾何尺寸的改變和泡沫鋁填充效果對(duì)薄壁金屬管吸能性能的影響。首先通過(guò)軸向壓縮實(shí)驗(yàn)給出了兩組試件的變形模態(tài)圖與載荷-位移曲線,利用ABAQUS軟件建立了與實(shí)驗(yàn)對(duì)應(yīng)的有限元模型,并與實(shí)驗(yàn)進(jìn)行對(duì)比,驗(yàn)證了有限元模型的可靠性與有效性;之后利用初始峰值載荷、平均載荷、比吸能、壓縮率等吸能評(píng)價(jià)指標(biāo)系統(tǒng)地分析了圓管的壁厚、高度、直徑和泡沫鋁的填充對(duì)其能量吸收的影響。

    1數(shù)值模擬

    為研究準(zhǔn)靜態(tài)軸向載荷作用下泡沫鋁填充薄壁金屬管的吸能性能,制備了7種不同幾何尺寸的泡沫鋁填充薄壁金屬圓管與空?qǐng)A管。其中,金屬外管的高度為h=30,45,60 mm;壁厚為d=0.9,1.2,1.5 mm;外徑為D=30,35,40 mm。針對(duì)不同的薄壁金屬外管,內(nèi)部填充泡沫鋁的高度與金屬圓管相同,直徑為薄壁金屬圓管的內(nèi)徑。

    1.1有限元模型

    利用有限元軟件ABAQUS/Explicit,對(duì)泡沫鋁填充薄壁金屬圓管在準(zhǔn)靜態(tài)軸向載荷作用下的吸能性能進(jìn)行了數(shù)值模擬,并與空管進(jìn)行了對(duì)比。有限元模型主要由上端剛性板、中間泡沫鋁填充圓管結(jié)構(gòu)或圓管、下端剛性板3部分組成。圓管或者泡沫鋁填充管放置在下端剛性板的中間,上端剛性板向下移動(dòng)實(shí)現(xiàn)對(duì)管結(jié)構(gòu)的準(zhǔn)靜態(tài)軸向加載。下端的剛性板約束了全部6個(gè)自由度,而上端剛性板約束了除Z軸方向的平動(dòng)自由度以外的其他自由度。

    圓管采用4節(jié)點(diǎn)殼單元(S4R),沿厚度方向設(shè)置5個(gè)積分點(diǎn)。泡沫鋁采用8節(jié)點(diǎn)六面體實(shí)體單元(C3D8R)。上、下的剛性板為離散剛體模型,采用六面體網(wǎng)格劃分。圓管與泡沫鋁之間,圓管與上、下剛板之間,以及圓管自身在壓縮過(guò)程中的接觸采用通用接觸,接觸面間的摩擦系數(shù)為0.2。

    1.2有限元模型的可靠性驗(yàn)證

    為了驗(yàn)證建立的有限元模型的有效性與計(jì)算的可靠性,與實(shí)驗(yàn)結(jié)果進(jìn)行了比較。實(shí)驗(yàn)使用的薄壁金屬圓管材料為鋁合金AA6061-T6,彈性模量E=68.2 GPa,泊松比μ=0.33,密度ρ=2.7 mg/mm3,材料的應(yīng)力-應(yīng)變曲線如圖1所示。作為填充材料的泡沫鋁,密度ρ=0.33 mg/mm3,彈性模量E=80 MPa,泊松比μ=0.01,泡沫鋁材料的應(yīng)力-應(yīng)變曲線如圖2所示,其平臺(tái)應(yīng)力σp=3.6 MPa。實(shí)驗(yàn)中泡沫鋁填充物自由放入圓管內(nèi),其柱面與管內(nèi)壁不做任何粘結(jié)處理。實(shí)驗(yàn)在萬(wàn)能材料試驗(yàn)機(jī)上完成,試件中心放在底板上,上板以3 mm/min的速度向下移動(dòng),通過(guò)計(jì)算機(jī)自動(dòng)記錄載荷-位移曲線。空管與填充管的變形通過(guò)照相機(jī)以間隔40~60 s記錄。

    圖1 鋁合金的拉伸應(yīng)力-應(yīng)變關(guān)系曲線Fig.1 Tensile stress-strain curve of aluminum alloy

    圖2 泡沫鋁的應(yīng)力-應(yīng)變關(guān)系曲線Fig.2 Stress-strain curve of aluminum foam

    圖3 圓管的實(shí)驗(yàn)與模擬變形模態(tài)對(duì)比Fig.3 Comparison of experimental and simulativedeformation modes of tube

    圖4 圓管的實(shí)驗(yàn)與模擬載荷-位移曲線對(duì)比Fig.4 Comparison of experimental and simulativeforce-displacement curves of tube

    圖3和圖4分別給出了直徑為20 mm,高度為30 mm,壁厚為1.7 mm的鋁合金圓管在準(zhǔn)靜態(tài)壓縮實(shí)驗(yàn)與有限元模擬中的變形模態(tài)與載荷-位移曲線圖。由圖可知,圓管的準(zhǔn)靜態(tài)實(shí)驗(yàn)和數(shù)值模擬中變形模態(tài)相似,都產(chǎn)生了2個(gè)軸對(duì)稱(chēng)的褶皺,并且都是向外折疊,圓管的下端都是向外翻出,沒(méi)有形成第三個(gè)完整的褶皺。數(shù)值模擬的載荷-位移曲線與準(zhǔn)靜態(tài)實(shí)驗(yàn)吻合較好。

    圖5和圖6分別給出了直徑為30 mm,高度為30 mm,壁厚為1.2 mm,泡沫鋁密度為0.33 mg/mm3的泡沫鋁填充管在準(zhǔn)靜態(tài)壓縮實(shí)驗(yàn)和數(shù)值模擬中的變形模態(tài)與載荷-位移曲線。從圖中可以看出,數(shù)值模擬與準(zhǔn)靜態(tài)實(shí)驗(yàn)中的變形模態(tài)基本相同,均產(chǎn)生了2個(gè)軸對(duì)稱(chēng)的褶皺,載荷-位移曲線吻合較好。

    圖5 泡沫鋁填充管的實(shí)驗(yàn)與模擬變形模態(tài)對(duì)比Fig.5 Comparison of experimental and simulativedeformation modes of foam-filled tube

    圖6 泡沫鋁填充管實(shí)驗(yàn)與模擬載荷-位移曲線對(duì)比Fig.6 Comparison of experimental and simulativeforce-displacement curves of foam-filled tube

    2數(shù)值模擬結(jié)果與討論

    2.1圓管的承載力和吸能

    分別對(duì)不同高度h、不同壁厚d、不同直徑D的鋁合金圓管進(jìn)行準(zhǔn)靜態(tài)軸向壓縮載荷下的數(shù)值模擬,分析了幾何尺寸對(duì)其承載力和吸能性能的影響。圓管的材料為鋁合金,與實(shí)驗(yàn)中的材料參數(shù)相同。

    圖7給出了直徑30 mm,高度45 mm,壁厚d分別為0.9,1.2,1.5 mm的3種不同圓管的載荷-位移曲線、吸能-位移曲線,以及平均載荷-位移曲線。從載荷-位移曲線a可以看出,圓管的軸向壓縮過(guò)程分為明顯的3個(gè)階段:初始彈性階段,穩(wěn)定漸進(jìn)屈曲階段,密實(shí)化階段。隨著壁厚的增大,圓管的承載力也增大,其初始峰值載荷、平均載荷以及總能量吸收都有提高。

    圖7 不同壁厚d的圓管準(zhǔn)靜態(tài)軸向壓縮曲線Fig.7 Curves of tubes with different thicknesses d under quasi-static axial compression

    圖8 不同高度h的圓管準(zhǔn)靜態(tài)軸向壓縮曲線Fig.8 Curves of tubes with different heights h under quasi-static axial compression

    圖8為直徑30 mm,壁厚0.9 mm,高度h分別為30,45,60 mm的3種不同圓管的載荷-位移曲線、吸能-位移曲線以及平均載荷-位移曲線。由圖可知,圓管的軸向壓縮位移隨著高度的增加而增大。在進(jìn)入密實(shí)階段之前,不同高度的圓管的初始峰值載荷、平均載荷以及總能量吸收基本相同。因此在選擇圓管的高度時(shí)主要考慮所需要的壓縮行程。

    圖9為壁厚1.5 mm,高度45 mm,直徑D分別為30,35,40 mm的3種不同圓管的載荷-位移曲線、吸能-位移曲線以及平均載荷-位移曲線。從圖中可以看出,圓管的直徑越大,壓縮到同一位移時(shí)的壓潰力、平均載荷以及能量吸收都越大,管的初始峰值載荷也隨著直徑的增大而增大。

    圖9 不同直徑D的圓管準(zhǔn)靜態(tài)軸向壓縮曲線Fig.9 Curves of tubes with different diameters D under quasi-static axial compression

    2.2泡沫鋁填充圓管的承載力和吸能

    泡沫鋁填充管中的外部金屬圓管與上一節(jié)中的圓管幾何尺寸相同,泡沫鋁的直徑為管的內(nèi)徑,高度與外管的高度相同,其外部的圓管與泡沫鋁的材料參數(shù)與實(shí)驗(yàn)中試件的材料參數(shù)相同。

    圖10 不同壁厚d的泡沫填充管準(zhǔn)靜態(tài)軸向壓縮曲線Fig.10 Curves of foam-filled tubes with different thicknesses d under quasi-static axial compression

    圖11 不同高度h的泡沫填充管準(zhǔn)靜態(tài)軸向壓縮曲線Fig.11 Curves of foam-filled tubes with different heights h under quasi-static axial compression

    圖12 不同直徑D的泡沫填充管準(zhǔn)靜態(tài)軸向壓縮曲線Fig.12 Curves of foam-filled tubes with different diameters D under quasi-static axial compression

    圖10—圖12分別給出了不同的壁厚d、不同高度h、不同直徑D的泡沫鋁填充圓管的載荷-位移曲線、吸能-位移曲線以及平均載荷-位移曲線。可以看出,與圓管相比,泡沫鋁填充管的軸向載荷波動(dòng)較小,在達(dá)到初始峰值載荷之后,呈逐漸上升趨勢(shì)。其原因可能是在軸向壓縮過(guò)程中,圓管和填充物泡沫鋁之間存在相互作用,這種相互作用效應(yīng)隨著壓縮位移的增大而增大,并且泡沫鋁本身也隨著壓縮位移的增大而逐漸被壓實(shí)。泡沫鋁填充管的平均載荷隨著位移的增大而逐漸增大,不同于空管的平均載荷在漸進(jìn)屈曲階段幾乎保持不變的情況。由圖10可知,隨著填充管壁厚的增加,其初始峰值載荷、平均載荷和總吸能增加。由圖11可以看出,填充管的高度對(duì)其初始峰值載荷、平均載荷和總吸能的影響較小,在密實(shí)化之前,不同高度填充管的壓縮載荷、總吸能與平均載荷幾乎重合,表明填充管的吸能與管的高度成線性關(guān)系。由圖12可知,隨著填充管直徑的增加,其初始峰值載荷、總吸能,平均載荷都增大。

    2.3吸能指標(biāo)

    為評(píng)價(jià)泡沫鋁填充管與圓管的耐撞性,引入了以下幾種參數(shù)[10]:

    1) 初始峰值載荷Fmax;

    2) 壓縮到某一位移時(shí)的總吸收能量E;

    3) 比吸能ESA(Specific Energy Absorption),即單位質(zhì)量的試件所吸收的能量:

    (1)

    式中:m為對(duì)應(yīng)模型的質(zhì)量;F為管在某時(shí)刻的軸向壓縮載荷;l為該時(shí)刻載荷所對(duì)應(yīng)的位移值。

    4) 壓縮到某一位移時(shí)的平均載荷Fm:

    (2)

    5) 壓縮力效率ECF(Crush Force Efficiency),

    定義為初始峰值載荷Fmax與平均載荷Fm之比:

    (3)

    表1給出了不同尺寸的泡沫鋁填充管和圓管在壓縮距離為初始管長(zhǎng)的80%時(shí)的各項(xiàng)吸能指標(biāo)。

    表1 圓管及泡沫填充管的吸能評(píng)價(jià)指標(biāo)

    2.3.1不同壁厚的泡沫填充管與圓管的吸能性能

    考慮直徑為30 mm,高度為45 mm,壁厚d分別為0.9,1.2,1.5 mm的泡沫填充管與空?qǐng)A管的吸能指標(biāo)。表1給出了壓縮距離為原管長(zhǎng)的80%時(shí)的比吸能、平均載荷、初始峰值載荷及壓縮力效率,可以看出,空管的比吸能、平均載荷、初始峰值載荷及壓縮力效率都隨著管壁厚的增大而增大。每一種規(guī)格的空管填充了泡沫鋁之后,其各項(xiàng)吸能指標(biāo)都較原來(lái)有所提高;但不同壁厚的管填充泡沫之后吸能指標(biāo)提高的程度有所差別,壁厚越小的管件,對(duì)應(yīng)的泡沫鋁填充管的各項(xiàng)吸能指標(biāo)提升的比例越大,即泡沫填充效果越明顯。如,壁厚為0.9 mm的圓管填充泡沫鋁之后其比吸能、平均載荷、初始峰值載荷及壓縮力效率較原來(lái)分別提高了17.7%,125.8%,4.4%,116.5%;壁厚為1.2 mm的填充管相應(yīng)指標(biāo)分別提高了5.3%,76.3%,3.2%,70.8%;壁厚為1.5 mm的填充管對(duì)應(yīng)的數(shù)值為1.9%,55.7%,2.2%,53%。

    2.3.2不同高度的泡沫填充管與圓管的吸能性能

    為對(duì)比直徑為30 mm,壁厚為0.9 mm,管的高度h分別為30,45,60 mm的泡沫填充管與空?qǐng)A管的吸能性能,表1給出了其壓縮距離為初始管長(zhǎng)的80%時(shí)的比吸能、平均載荷、初始峰值載荷及壓縮力效率值。從表中可以看出,泡沫鋁填充管的比吸能、平均載荷、初始峰值載荷與壓縮力效率均高于空管。隨著高度的增加,這些指標(biāo)的變化較小,表明高度對(duì)泡沫鋁填充管與空管的吸能性影響較小。因而增加泡沫鋁填充管或空管的高度,可以增加其可利用的壓縮距離,從而使整體結(jié)構(gòu)的總能量吸收得到提高,對(duì)其他吸能指標(biāo)的影響較小。該組對(duì)比中泡沫鋁填充管的壓縮力效率都大于1,這是由于該組模型壁厚都較薄;管的壁厚越薄,填充結(jié)構(gòu)受軸向壓縮時(shí)泡沫鋁的承載力與吸能所占比例越大;隨著壓縮位移的增大,管和泡沫之間的相互作用也隨之增大,泡沫鋁也慢慢被壓實(shí),使得壓縮載荷在達(dá)到峰值之后呈現(xiàn)逐漸上升趨勢(shì),因而壓縮力效率較高。

    2.3.3不同直徑的泡沫填充管與圓管的吸能性能

    考慮高度為45 mm,壁厚為1.5 mm,直徑D分別為30,35,40 mm的泡沫鋁填充管與空?qǐng)A管的吸能指標(biāo)。從表1給出的泡沫填充管與空?qǐng)A管的比吸能、平均載荷、初始峰值載荷及壓縮力效率對(duì)比中可以看出,泡沫鋁填充管的各項(xiàng)吸能指標(biāo)均高于空管。隨著直徑的增大,泡沫鋁填充管的比吸能、平均載荷與壓縮力效率較空管增加的比例增大。如,直徑為30 mm的管件填充泡沫鋁之后其比吸能、平均載荷及壓縮力效率較空管分別提高了1.9%,55.7%,53%;直徑為35 mm的填充管相應(yīng)指標(biāo)分別提高了10.4%,78.3%,70.8%;直徑為40 mm的填充管對(duì)應(yīng)的數(shù)值為20.8%,109.2%,105.9%??展艿某跏挤逯递d荷以及平均載荷都隨著管直徑的增大而增大,相反的比吸能以及壓縮力效率都隨著管的直徑的增大而減小。

    3結(jié)論

    利用有限元軟件ABAQUS對(duì)泡沫鋁填充薄壁金屬圓管和圓管在軸向準(zhǔn)靜態(tài)壓縮時(shí)的吸能性能進(jìn)行了數(shù)值模擬。首先,通過(guò)與實(shí)驗(yàn)結(jié)果的對(duì)比驗(yàn)證了有限元模型的有效性與計(jì)算的可靠性;之后,根據(jù)建立的模型系統(tǒng)地研究了管的高度、壁厚、直徑以及泡沫鋁的填充對(duì)圓管吸能性能的影響;利用初始峰值載荷、平均載荷、總吸能與比吸能等吸能評(píng)價(jià)指標(biāo)對(duì)結(jié)構(gòu)的吸能特性進(jìn)行了分析。結(jié)果表明,鋁合金圓管受到軸向壓縮時(shí)的載荷-位移曲線分為明顯的三個(gè)階段:初始彈性階段,塑性變形階段,密實(shí)化階段。圓管在準(zhǔn)靜態(tài)軸向壓縮時(shí)的初始峰值載荷、平均載荷、總吸能、比吸能隨著管壁厚的增加而增大;初始峰值載荷、平均載荷以及總吸能隨著管直徑的增大而增大,但是比吸能和壓縮力效率隨管直徑的增大而減小;當(dāng)管的高度增大時(shí),其比吸能、平均載荷、初始峰值載荷及壓縮力效率都沒(méi)有發(fā)生明顯變化,即高度對(duì)空管的吸能效率影響不是很大,與空管相比,泡沫鋁填充管的比吸能、平均載荷、初始峰值載荷及壓縮力效率均高于對(duì)應(yīng)的空管,但增大的程度隨著管的幾何參數(shù)的不同而變化。填充管的壁厚越小,其各項(xiàng)吸能指標(biāo)較空管提升的比例越大;填充管的直徑越大,其比吸能、平均載荷、壓縮力效率增大的比例越大;填充管的高度對(duì)泡沫填充管的吸能性能影響較小。研究表明,泡沫鋁填充薄壁金屬圓管是一種理想能量吸收元件。

    參考文獻(xiàn):

    [1]ANDREWS K R F,ENGLAND G L,GHANI E.Classification of the axial collapse of cylindrical tubes under quasi-static loading[J].International Journal of Mechanical Sciences,1983,25(9/10):687-696.

    [2]GUILLOW S R,LU G,GRZEBIETA R H.Quasi-static axial compression of thin-walled circular aluminium tubes[J].International Journal of Mechanical Sciences,2001,43(9):2103-2123.

    [3]CHEN D H,VSHIJIMA K.Estimation of the initial peak load for circular tubes subjected to axial impact[J].Thin Walled Structures,2011,49(7):889-898.

    [4]KARAGIOZOVA D,NORMAN J.Dynamic effects on buckling and energy absorption of cylindrical shells under axial impact[J].Thin-Walled Structures,2001,39(7):583-610.

    [5]SEITZBERGER M,RAMMERSTORFER F G,DEGISCGER H P,et al.Crushing of axially compressed steel tubes filled with aluminium foam[J].Acta Mechanica,1997,125(1/2/3/4):93-105.

    [6]桂良進(jìn),范子杰,王青春.泡沫填充圓管的軸向壓縮能量吸收特性[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,43(11):1526-1529.

    [7]桂良進(jìn),范子杰,王青春.泡沫填充圓管的動(dòng)態(tài)軸向壓縮吸能特性[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,44(5):709-712.

    [8]SANTOSA S,WIERZBICKI T.Experimental and numerical studies of foam-filled sections[J].International Journal of Impact Engineering,2000,24(5):509-534.

    [9]張志新,田愛(ài)琴,車(chē)全偉,等.高速列車(chē)車(chē)體端部吸能結(jié)構(gòu)研究[J].機(jī)車(chē)電傳動(dòng),2013(1):43-47.

    [10]李志斌,虞吉林,鄭志軍,等.薄壁管及其泡沫金屬填充結(jié)構(gòu)耐撞性的實(shí)驗(yàn)研究[J].實(shí)驗(yàn)力學(xué),2012,27(1):77-86.

    (編輯:龐富祥)

    Study on Energy Absorption Properties of Thin-walled Circular Tubes Filled with Aluminum Foam under Axial Compression

    HUANG Rui,LIU Zhifang,LU Guoyun,YAN Qingrong

    (InstituteofAppliedMechanicsandBiomedicalEngineering,TaiyuanUniversityofTechnology,Taiyuan030024,China)

    Abstract:Numerical simulation analysis of tubes filled with aluminum foam of seven different geometric sizes under quasi-static axial compression was conducted using the finite element software ABAQUS. The influence of the height,thickness and diameter of tubes on the energy absorption properties of the tubes and the foam-filled tubes was studied systematically.Comparative analysis was conducted on the evaluation indexes of energy absorption such as average load,initial peak load, specific energy absorption (ESA) and compression force efficiency (ECF) for seven types of foam-filled tubes. The results show that the thickness and diameter of the tubes have a significant effect on energy absorption properties of the structures when the tubes with aluminum foam filler are under quasi-static axial compression.The effect of height is not evident but the increase of the height can improve the total energy absorption of the tube. Aluminum foam filler improves the bearing capacity,total energy absorption, specific energy absorption and compression force efficiency of the tubes. The results of the research can provide some technical basis for the design of ideal energy absorption devices.

    Key words:tubes filled with aluminum foam;axial compression;energy absorption;numerical simulation

    中圖分類(lèi)號(hào):O347.3

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

    DOI:10.16355/j.cnki.issn1007-9432tyut.2016.01.020

    作者簡(jiǎn)介:黃睿(1991-),女,山西運(yùn)城人,碩士生,主要從事泡沫鋁夾芯沖擊動(dòng)力學(xué)研究,(E-mail)rui152@126.com通訊作者:閆慶榮(1960-),女,高級(jí)工程師,主要從事沖擊動(dòng)力學(xué)研究,(E-mail)yanqingrong@tyut.edu.cn,(Tel)0351-6010560

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目:漿砌多級(jí)構(gòu)建(層合)殼體的沖擊動(dòng)力學(xué)行為研究(11372209);山西省自然科學(xué)基金資助項(xiàng)目:柱殼結(jié)構(gòu)中的非線性彌散波與波致失效的研究(2013011005-2)

    收稿日期:2015-04-30

    文章編號(hào):1007-9432(2016)01-0101-07

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無(wú)支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開(kāi)切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    欧美日韩国产亚洲二区| 乱人视频在线观看| 亚洲av免费高清在线观看| 91麻豆av在线| 国产免费男女视频| 久久国产乱子免费精品| 俄罗斯特黄特色一大片| 美女 人体艺术 gogo| 精品久久久久久久久久免费视频| 亚洲成人免费电影在线观看| 高清日韩中文字幕在线| 少妇猛男粗大的猛烈进出视频 | 十八禁网站免费在线| 最近最新免费中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 18+在线观看网站| 两个人的视频大全免费| 免费高清视频大片| 午夜精品在线福利| 啦啦啦观看免费观看视频高清| 欧美最新免费一区二区三区| 俄罗斯特黄特色一大片| 少妇人妻精品综合一区二区 | 国产欧美日韩一区二区精品| 成人一区二区视频在线观看| 久久热精品热| 变态另类丝袜制服| 日韩 亚洲 欧美在线| 日韩大尺度精品在线看网址| 在线天堂最新版资源| 婷婷精品国产亚洲av在线| 美女大奶头视频| 亚洲一级一片aⅴ在线观看| 久久久久久久久久黄片| 国产欧美日韩一区二区精品| 国内精品久久久久精免费| 亚洲av五月六月丁香网| 久久精品夜夜夜夜夜久久蜜豆| 欧美性感艳星| 成年免费大片在线观看| 日韩国内少妇激情av| 人人妻,人人澡人人爽秒播| 婷婷六月久久综合丁香| 久久精品综合一区二区三区| 亚洲五月天丁香| 精品99又大又爽又粗少妇毛片 | 国产精华一区二区三区| 国产av一区在线观看免费| 国产黄片美女视频| 国产91精品成人一区二区三区| 天天躁日日操中文字幕| 国产又黄又爽又无遮挡在线| 我的老师免费观看完整版| 99热只有精品国产| 人妻少妇偷人精品九色| netflix在线观看网站| 嫩草影视91久久| 国产精品三级大全| 国产成人aa在线观看| 成人欧美大片| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 国产精品永久免费网站| 永久网站在线| 国产精品日韩av在线免费观看| 国产人妻一区二区三区在| 国产av麻豆久久久久久久| 亚洲成人久久性| 国产av不卡久久| 真人一进一出gif抽搐免费| 精品久久久噜噜| 九九在线视频观看精品| 亚洲国产欧洲综合997久久,| 午夜精品久久久久久毛片777| 国产色婷婷99| 日韩欧美国产一区二区入口| 午夜亚洲福利在线播放| 亚洲av不卡在线观看| 国语自产精品视频在线第100页| 久久亚洲精品不卡| a在线观看视频网站| 观看免费一级毛片| 亚洲第一区二区三区不卡| 在线免费观看的www视频| 国内精品宾馆在线| 美女被艹到高潮喷水动态| 成人毛片a级毛片在线播放| 男女边吃奶边做爰视频| 桃色一区二区三区在线观看| 国内精品一区二区在线观看| 久久精品国产亚洲av香蕉五月| 欧美黑人巨大hd| 久久精品国产亚洲av天美| av福利片在线观看| 日韩欧美三级三区| 国产欧美日韩一区二区精品| 人人妻人人看人人澡| 午夜免费成人在线视频| 久久久久久久精品吃奶| 亚洲在线自拍视频| 波多野结衣高清无吗| 国产精品自产拍在线观看55亚洲| 亚洲一区高清亚洲精品| 三级男女做爰猛烈吃奶摸视频| 久久久久久伊人网av| 国产亚洲av嫩草精品影院| 亚洲在线观看片| 国产女主播在线喷水免费视频网站 | 亚洲真实伦在线观看| 18禁在线播放成人免费| 精品午夜福利视频在线观看一区| 国产白丝娇喘喷水9色精品| 国产免费av片在线观看野外av| www.www免费av| 精品久久久噜噜| 国产乱人视频| 一级毛片久久久久久久久女| 校园春色视频在线观看| 欧美日韩黄片免| 亚洲国产欧美人成| 亚洲一区高清亚洲精品| 在线观看66精品国产| 免费电影在线观看免费观看| 亚洲人成伊人成综合网2020| 又黄又爽又免费观看的视频| 91麻豆精品激情在线观看国产| 禁无遮挡网站| 免费av观看视频| 精品久久久久久,| 一个人观看的视频www高清免费观看| 高清在线国产一区| 男女啪啪激烈高潮av片| 亚洲七黄色美女视频| 别揉我奶头 嗯啊视频| 精品一区二区三区视频在线观看免费| 久久精品综合一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区二区三区色噜噜| 欧美日韩中文字幕国产精品一区二区三区| 人妻久久中文字幕网| 亚洲av二区三区四区| 亚洲av美国av| 最新在线观看一区二区三区| av天堂中文字幕网| 欧美一区二区国产精品久久精品| 蜜桃久久精品国产亚洲av| 亚洲成人免费电影在线观看| 99九九线精品视频在线观看视频| 国产av一区在线观看免费| av在线蜜桃| a级毛片a级免费在线| 免费看日本二区| .国产精品久久| 精品午夜福利在线看| av女优亚洲男人天堂| 国产精华一区二区三区| 国产亚洲精品av在线| 久久精品久久久久久噜噜老黄 | 此物有八面人人有两片| 日本 欧美在线| 99热精品在线国产| 精品人妻视频免费看| av视频在线观看入口| 午夜久久久久精精品| 日韩欧美 国产精品| 中文字幕熟女人妻在线| 欧洲精品卡2卡3卡4卡5卡区| 日本免费一区二区三区高清不卡| 亚洲av熟女| 免费黄网站久久成人精品| 日本-黄色视频高清免费观看| 色综合站精品国产| 欧美极品一区二区三区四区| 人妻夜夜爽99麻豆av| 国产精品女同一区二区软件 | 精品久久国产蜜桃| 亚洲人成网站在线播| 日本黄大片高清| 欧美性感艳星| 日韩大尺度精品在线看网址| 精品国产三级普通话版| 精品人妻熟女av久视频| 在线观看66精品国产| 国国产精品蜜臀av免费| 日本在线视频免费播放| 午夜福利18| 国模一区二区三区四区视频| 尾随美女入室| 日本精品一区二区三区蜜桃| 一级av片app| av.在线天堂| 亚洲狠狠婷婷综合久久图片| 日本在线视频免费播放| 亚洲精品一区av在线观看| 天天躁日日操中文字幕| 日韩欧美在线乱码| 精品一区二区免费观看| 欧美丝袜亚洲另类 | 简卡轻食公司| 日本成人三级电影网站| 国产精品自产拍在线观看55亚洲| 国产高清三级在线| 亚洲精品日韩av片在线观看| 亚洲四区av| 熟妇人妻久久中文字幕3abv| 长腿黑丝高跟| 九九爱精品视频在线观看| 国产免费一级a男人的天堂| 观看美女的网站| 国产 一区 欧美 日韩| 男人舔奶头视频| 午夜福利在线观看吧| 看片在线看免费视频| 国内揄拍国产精品人妻在线| 中文字幕精品亚洲无线码一区| 黄色日韩在线| 免费观看人在逋| 真实男女啪啪啪动态图| 3wmmmm亚洲av在线观看| 欧美xxxx性猛交bbbb| 99久久九九国产精品国产免费| 女生性感内裤真人,穿戴方法视频| 欧美另类亚洲清纯唯美| 日韩人妻高清精品专区| 久久九九热精品免费| 成人永久免费在线观看视频| 床上黄色一级片| 欧美日韩亚洲国产一区二区在线观看| 亚洲真实伦在线观看| 日本一二三区视频观看| 欧美精品啪啪一区二区三区| 国产精品日韩av在线免费观看| 亚洲精品乱码久久久v下载方式| 国产精品综合久久久久久久免费| videossex国产| 国产精品人妻久久久久久| 国产蜜桃级精品一区二区三区| 无人区码免费观看不卡| 九九在线视频观看精品| av视频在线观看入口| 国产伦精品一区二区三区四那| 国产不卡一卡二| 可以在线观看的亚洲视频| 在线观看66精品国产| 偷拍熟女少妇极品色| 国内精品久久久久精免费| 欧美bdsm另类| 国产69精品久久久久777片| 国内精品美女久久久久久| 欧美成人性av电影在线观看| 欧美潮喷喷水| 日本 av在线| 三级毛片av免费| 熟女电影av网| 欧美激情在线99| 99久久精品一区二区三区| 国产精品av视频在线免费观看| 热99re8久久精品国产| 国产熟女欧美一区二区| 久久6这里有精品| 中文字幕高清在线视频| 天堂动漫精品| 十八禁国产超污无遮挡网站| 久久99热这里只有精品18| 最新在线观看一区二区三区| 久久久午夜欧美精品| 久久国内精品自在自线图片| 免费av不卡在线播放| 欧美性猛交╳xxx乱大交人| 国产高清三级在线| 久久久久久伊人网av| 91精品国产九色| 又粗又爽又猛毛片免费看| 桃红色精品国产亚洲av| 黄色日韩在线| 色播亚洲综合网| 免费av观看视频| 亚洲av成人av| 最近最新中文字幕大全电影3| 美女大奶头视频| 亚洲人与动物交配视频| 精品一区二区三区人妻视频| 亚洲欧美清纯卡通| av福利片在线观看| 国产精品久久久久久av不卡| 一个人看视频在线观看www免费| 国产精品日韩av在线免费观看| 观看免费一级毛片| 91麻豆精品激情在线观看国产| 女的被弄到高潮叫床怎么办 | 别揉我奶头 嗯啊视频| 少妇人妻一区二区三区视频| 欧美色欧美亚洲另类二区| 两个人的视频大全免费| 国内精品宾馆在线| 亚洲电影在线观看av| 少妇猛男粗大的猛烈进出视频 | 成人欧美大片| 免费高清视频大片| 观看免费一级毛片| 色吧在线观看| 日本在线视频免费播放| 草草在线视频免费看| 成年免费大片在线观看| 自拍偷自拍亚洲精品老妇| 熟女人妻精品中文字幕| 1000部很黄的大片| 久久午夜亚洲精品久久| 国产熟女欧美一区二区| 俺也久久电影网| 欧美不卡视频在线免费观看| 狠狠狠狠99中文字幕| 国产精品不卡视频一区二区| 乱人视频在线观看| 中国美白少妇内射xxxbb| 国产一区二区激情短视频| 天堂网av新在线| 亚洲综合色惰| 又爽又黄a免费视频| 亚洲美女黄片视频| 国产精品久久视频播放| 午夜视频国产福利| 内射极品少妇av片p| 国产成人一区二区在线| 午夜福利18| 日韩高清综合在线| 联通29元200g的流量卡| 国产精品无大码| www.www免费av| 国产在线精品亚洲第一网站| 97碰自拍视频| 国产高清有码在线观看视频| 久久人妻av系列| 日日啪夜夜撸| 99riav亚洲国产免费| 免费av观看视频| 亚洲av.av天堂| 欧美色欧美亚洲另类二区| 麻豆成人av在线观看| 欧美一区二区国产精品久久精品| 欧美极品一区二区三区四区| 12—13女人毛片做爰片一| 最近中文字幕高清免费大全6 | 亚洲最大成人中文| 97超视频在线观看视频| 麻豆精品久久久久久蜜桃| 亚洲专区国产一区二区| a级毛片a级免费在线| 久久久久久久久久久丰满 | 超碰av人人做人人爽久久| 真实男女啪啪啪动态图| 久久亚洲精品不卡| 女人被狂操c到高潮| 一区二区三区高清视频在线| 最近中文字幕高清免费大全6 | 精品午夜福利在线看| 日本一本二区三区精品| 国产男人的电影天堂91| 欧美潮喷喷水| a级毛片a级免费在线| 少妇的逼水好多| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 久久精品综合一区二区三区| 亚洲国产日韩欧美精品在线观看| 少妇裸体淫交视频免费看高清| 在线观看舔阴道视频| 最近视频中文字幕2019在线8| 波多野结衣巨乳人妻| 国国产精品蜜臀av免费| 亚洲成av人片在线播放无| 亚洲无线在线观看| 亚洲性久久影院| 超碰av人人做人人爽久久| 伦精品一区二区三区| 超碰av人人做人人爽久久| 亚洲,欧美,日韩| 男女之事视频高清在线观看| 看十八女毛片水多多多| 俺也久久电影网| 亚洲电影在线观看av| 亚洲最大成人手机在线| 老熟妇仑乱视频hdxx| 99精品久久久久人妻精品| 亚洲成人久久性| 熟女人妻精品中文字幕| 最好的美女福利视频网| 老司机午夜福利在线观看视频| 精华霜和精华液先用哪个| 中出人妻视频一区二区| 九九久久精品国产亚洲av麻豆| 蜜桃久久精品国产亚洲av| 99久久精品热视频| 亚洲在线观看片| 少妇的逼好多水| 成人一区二区视频在线观看| 最近视频中文字幕2019在线8| 久久久久免费精品人妻一区二区| 亚洲av熟女| 又黄又爽又免费观看的视频| 日韩欧美精品免费久久| 免费在线观看影片大全网站| 午夜精品在线福利| 又紧又爽又黄一区二区| 中亚洲国语对白在线视频| 三级男女做爰猛烈吃奶摸视频| 一级毛片久久久久久久久女| 真人做人爱边吃奶动态| 欧美性猛交黑人性爽| 中国美白少妇内射xxxbb| 观看美女的网站| 黄色一级大片看看| 俺也久久电影网| 色哟哟哟哟哟哟| 一进一出抽搐动态| 少妇人妻一区二区三区视频| 日本在线视频免费播放| 亚洲国产欧洲综合997久久,| 午夜日韩欧美国产| 亚洲七黄色美女视频| 国产白丝娇喘喷水9色精品| 久久久久国内视频| 在线天堂最新版资源| 欧美激情在线99| 搡老熟女国产l中国老女人| 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 黄片wwwwww| x7x7x7水蜜桃| 日本免费一区二区三区高清不卡| 欧美xxxx黑人xx丫x性爽| 免费看光身美女| 最近视频中文字幕2019在线8| 国产av麻豆久久久久久久| 亚洲精品色激情综合| 直男gayav资源| 国产一区二区三区av在线 | 999久久久精品免费观看国产| 国产高清有码在线观看视频| 午夜日韩欧美国产| 免费高清视频大片| 老女人水多毛片| 男女做爰动态图高潮gif福利片| 免费看日本二区| 神马国产精品三级电影在线观看| 给我免费播放毛片高清在线观看| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久| 黄色丝袜av网址大全| www日本黄色视频网| 热99re8久久精品国产| 一区二区三区高清视频在线| 国产亚洲av嫩草精品影院| 国语自产精品视频在线第100页| 一夜夜www| 中国美白少妇内射xxxbb| 男女边吃奶边做爰视频| 久久亚洲精品不卡| 久久久久久久久久久丰满 | 亚洲图色成人| 亚洲熟妇熟女久久| 欧美日韩国产亚洲二区| 久久天躁狠狠躁夜夜2o2o| 国产免费一级a男人的天堂| 在线看三级毛片| 男人舔女人下体高潮全视频| 精品久久久噜噜| or卡值多少钱| 国产乱人伦免费视频| 99久久精品国产国产毛片| 老熟妇仑乱视频hdxx| 性欧美人与动物交配| 国产蜜桃级精品一区二区三区| 日韩大尺度精品在线看网址| 免费看美女性在线毛片视频| 变态另类成人亚洲欧美熟女| 国产成人aa在线观看| 麻豆久久精品国产亚洲av| 欧美性感艳星| 日本欧美国产在线视频| 麻豆一二三区av精品| 黄色女人牲交| 91精品国产九色| 2021天堂中文幕一二区在线观| 18禁黄网站禁片午夜丰满| 麻豆成人av在线观看| 99久久无色码亚洲精品果冻| 国产探花极品一区二区| 欧美性感艳星| 国产精品亚洲美女久久久| 俄罗斯特黄特色一大片| 欧美日本视频| 亚洲狠狠婷婷综合久久图片| 午夜免费男女啪啪视频观看 | 香蕉av资源在线| 成人二区视频| 91av网一区二区| 精品午夜福利视频在线观看一区| 国产精品国产三级国产av玫瑰| 人人妻人人看人人澡| 91久久精品国产一区二区成人| 日本黄大片高清| 国产69精品久久久久777片| 色综合亚洲欧美另类图片| 免费在线观看日本一区| а√天堂www在线а√下载| 免费观看在线日韩| 99热这里只有是精品在线观看| 男女做爰动态图高潮gif福利片| 窝窝影院91人妻| 日本免费a在线| 亚洲中文字幕日韩| 国产视频内射| 高清日韩中文字幕在线| 在线观看舔阴道视频| 中国美女看黄片| 国内少妇人妻偷人精品xxx网站| 91麻豆av在线| 99在线人妻在线中文字幕| 亚洲欧美日韩东京热| 午夜精品在线福利| 91av网一区二区| 国产成人a区在线观看| 欧美一区二区亚洲| 99视频精品全部免费 在线| 亚洲人成伊人成综合网2020| 他把我摸到了高潮在线观看| 亚洲精华国产精华液的使用体验 | 观看美女的网站| 国产精品久久久久久久电影| 欧美丝袜亚洲另类 | 少妇的逼水好多| 日本 欧美在线| 男女啪啪激烈高潮av片| 国产精品永久免费网站| or卡值多少钱| 中文亚洲av片在线观看爽| 桃色一区二区三区在线观看| 99热只有精品国产| 一个人看视频在线观看www免费| 亚洲最大成人手机在线| 美女免费视频网站| 久久久精品大字幕| 国产伦人伦偷精品视频| 999久久久精品免费观看国产| 一级av片app| 十八禁网站免费在线| 少妇熟女aⅴ在线视频| 亚洲国产色片| 黄片wwwwww| 日韩一区二区视频免费看| 欧美zozozo另类| 久久热精品热| 久久精品国产鲁丝片午夜精品 | 亚洲精品一卡2卡三卡4卡5卡| 熟女人妻精品中文字幕| 99久久精品热视频| 韩国av一区二区三区四区| 亚洲精品久久国产高清桃花| 午夜免费成人在线视频| 亚洲人成网站高清观看| 两个人视频免费观看高清| 久久久久久九九精品二区国产| 人妻久久中文字幕网| 国国产精品蜜臀av免费| 欧美日韩综合久久久久久 | 1000部很黄的大片| 精品久久久久久成人av| a级毛片a级免费在线| 99视频精品全部免费 在线| 亚洲一区高清亚洲精品| 长腿黑丝高跟| 网址你懂的国产日韩在线| 亚洲人成伊人成综合网2020| 精品乱码久久久久久99久播| 国产精品乱码一区二三区的特点| 精品人妻熟女av久视频| 一进一出好大好爽视频| 国产免费男女视频| 成人性生交大片免费视频hd| 最新中文字幕久久久久| 午夜免费男女啪啪视频观看 | 国产在线精品亚洲第一网站| 岛国在线免费视频观看| 国产三级中文精品| 国产高清三级在线| 亚洲aⅴ乱码一区二区在线播放| 欧美+亚洲+日韩+国产| 国产综合懂色| 欧美成人性av电影在线观看| 亚洲精品亚洲一区二区| 九色国产91popny在线| 狠狠狠狠99中文字幕| 成年女人看的毛片在线观看| 麻豆成人午夜福利视频| 男女啪啪激烈高潮av片| 如何舔出高潮| 国产精品人妻久久久影院| 国产精品乱码一区二三区的特点| 国产成人aa在线观看| 我的女老师完整版在线观看| 一个人看视频在线观看www免费| 熟女人妻精品中文字幕| 国产91精品成人一区二区三区| 欧美人与善性xxx| 能在线免费观看的黄片| 免费电影在线观看免费观看| 亚洲国产精品合色在线| 国产中年淑女户外野战色| 亚洲一级一片aⅴ在线观看| 国产毛片a区久久久久| 少妇人妻精品综合一区二区 | 国产免费一级a男人的天堂| 亚洲第一电影网av| 超碰av人人做人人爽久久|