亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        材料應(yīng)變率效應(yīng)對2類吸能結(jié)構(gòu)碰撞性能的影響

        2018-11-13 06:51:14馮悅肖守訥朱濤楊冰楊皓杰
        關(guān)鍵詞:圓錐效應(yīng)誘導(dǎo)

        馮悅,肖守訥,朱濤,楊冰,楊皓杰

        ?

        材料應(yīng)變率效應(yīng)對2類吸能結(jié)構(gòu)碰撞性能的影響

        馮悅,肖守訥,朱濤,楊冰,楊皓杰

        (西南交通大學(xué) 牽引動力國家重點實驗室,四川 成都,610031)

        軌道車輛吸能結(jié)構(gòu)變形時伴隨著明顯的應(yīng)變率效應(yīng),為了研究材料應(yīng)變率效應(yīng)對碰撞仿真結(jié)果的影響,以應(yīng)變率敏感材料Q235為研究對象,首先,通過實驗數(shù)據(jù)對比和理論解對比2種方法,證明考慮應(yīng)變率效應(yīng)的Q235材料模型的可靠性;然后,研究材料應(yīng)變率效應(yīng)對擠壓式和壓潰式圓錐管吸能結(jié)構(gòu)碰撞載荷與吸能、變形模式等耐撞性能的影響;在此基礎(chǔ)上,提出2種受材料應(yīng)變率效應(yīng)的影響不明顯、變形有序可控的誘導(dǎo)式吸能結(jié)構(gòu)。研究結(jié)果表明:對于不同結(jié)構(gòu),材料應(yīng)變率效應(yīng)影響因子不同,考慮材料應(yīng)變率強化效應(yīng)的擠壓式吸能結(jié)構(gòu)碰撞力放大因子比壓潰式圓錐管結(jié)構(gòu)的大;對于不同結(jié)構(gòu),材料應(yīng)變率效應(yīng)影響因子隨速度變化的幅度不同,擠壓式吸能結(jié)構(gòu)的應(yīng)變率效應(yīng)放大因子隨速度變化幅度比壓潰式的??;材料應(yīng)變率效應(yīng)對第二類吸能結(jié)構(gòu)的變形模式影響顯著,應(yīng)變率強化效應(yīng)減少了應(yīng)變率較大的變形,增加或促進了應(yīng)變率較小的變形;誘導(dǎo)式吸能結(jié)構(gòu)的變形模式受材料應(yīng)變率效應(yīng)的影響不明顯,變形有序可控,在工程設(shè)計中可考慮采用誘導(dǎo)結(jié)構(gòu),以減少材料應(yīng)變率效應(yīng)對結(jié)構(gòu)變形模式的影響。

        軌道車輛;碰撞性能;吸能結(jié)構(gòu);應(yīng)變率效應(yīng)

        吸能結(jié)構(gòu)在軌道車輛碰撞過程中主要通過塑性變形吸收能量,其吸能性能與軌道車輛碰撞響應(yīng)關(guān)系密切,是軌道車輛被動安全中非常重要的部分[1?7]。由于碰撞過程中沖擊速度較大,在塑性變形過程中伴隨著明顯的應(yīng)變率效應(yīng)。然而,目前進行軌道車輛碰撞仿真時一般未考慮材料的應(yīng)變率效應(yīng),已有研究表明材料的應(yīng)變率效應(yīng)會對結(jié)構(gòu)碰撞性能產(chǎn)生明顯的影響。廖其紅等[8]針對閉口帽型薄壁梁在軸向碰撞載荷下的壓潰折曲過程進行了研究,對比有限元計算結(jié)果與實驗結(jié)果發(fā)現(xiàn)考慮材料應(yīng)變率效應(yīng)會使結(jié)構(gòu)吸能特性與實驗結(jié)果更接近,應(yīng)變率是造成碰撞實驗和準靜態(tài)壓縮實驗結(jié)果差異的最主要原因。孟廣成等[9]在對護欄的碰撞仿真研究中發(fā)現(xiàn)考慮應(yīng)變率影響比不考慮應(yīng)變率影響的仿真結(jié)果與實驗結(jié)果更加接近。王自力 等[10?13]研究了材料應(yīng)變率效應(yīng)對汽車和船舶碰撞中簡單結(jié)構(gòu)和復(fù)雜結(jié)構(gòu)的碰撞影響,結(jié)果表明材料的應(yīng)變率效應(yīng)對簡單結(jié)構(gòu)的影響比較明顯,會改變結(jié)構(gòu)變形模式、失效模式等,但對復(fù)雜結(jié)構(gòu)的影響不明顯,同時,材料的應(yīng)變率強化效應(yīng)會增大結(jié)構(gòu)的壓潰力,提高結(jié)構(gòu)的承載能力和吸能結(jié)構(gòu)的吸能能力。而楊超等[14]在對軌道車輛吸能結(jié)構(gòu)應(yīng)變率效應(yīng)影響的研究中得出5083H111鋁合金的負應(yīng)變率敏感性會減少結(jié)構(gòu)的吸能量。目前關(guān)于材料應(yīng)變率效應(yīng)影響的研究多集中于壓潰式結(jié)構(gòu),而且大多未對其影響進行深入分析,而軌道車輛吸能結(jié)構(gòu)種類繁多,變形模式多樣,僅研究壓潰式結(jié)構(gòu)受材料應(yīng)變率效應(yīng)的影響是不夠的。根據(jù)CALLADINE等[15]的研究指出,吸能結(jié)構(gòu)可以根據(jù)碰撞力位移曲線的特征分為2類:第一類吸能結(jié)構(gòu)的碰撞力位移曲線在沖擊發(fā)生開始階段迅速增加,隨著沖擊進行,碰撞力波動很小,基本保持平穩(wěn);第二類吸能結(jié)構(gòu)的碰撞力迅速增加之后有1個明顯峰值,下降后隨著沖擊進行明顯波動。這2類吸能結(jié)構(gòu)的變形模式具有不同的碰撞速度敏感性[16]。軌道車輛中常用的第一類吸能結(jié)構(gòu)主要有擠壓式、膨脹式、翻卷式等,第二類吸能結(jié)構(gòu)主要為壓潰式[17?19]。對于第一類吸能結(jié)構(gòu),翻卷式吸能結(jié)構(gòu)的翻卷變形對材料、結(jié)構(gòu)要求嚴格,變形模式不易發(fā)生;膨脹式吸能結(jié)構(gòu)軸向占用空間大;擠壓式吸能結(jié)構(gòu)的變形穩(wěn)定、有序,是理想的吸能結(jié)構(gòu)。第二類吸能結(jié)構(gòu)主要是各種截面的壓潰結(jié)構(gòu),不同截面壓潰式吸能結(jié)構(gòu)中圓錐管的壓潰可控性更強,綜合耐撞性能更優(yōu)。本文作者主要研究對比材料的應(yīng)變率強化效應(yīng)對擠壓式吸能結(jié)構(gòu)和壓潰式圓錐管吸能結(jié)構(gòu)軸向碰撞性能的影響,并分析產(chǎn)生差異的原因,進而指導(dǎo)工程中吸能結(jié)構(gòu)設(shè)計。

        1 材料模型

        為了更清晰地表現(xiàn)材料應(yīng)變率效應(yīng)對結(jié)構(gòu)碰撞性能的影響,使用典型的應(yīng)變率敏感材料Q235。材料模型為彈塑性本構(gòu)模型,通過廣義Cowper?Symonds模型考慮應(yīng)變率效應(yīng),公式為

        式中:d和s分別為高應(yīng)變率下的動態(tài)流動應(yīng)力和準靜態(tài)流動應(yīng)力;和為應(yīng)變率相關(guān)的材料參數(shù)。對于Q235低碳鋼,屈服極限為235 MPa,切線模量設(shè)為2 GPa。在考慮應(yīng)變率強化效應(yīng)時,=6 844 s?1,=3.91[13];未考慮應(yīng)變率強化效應(yīng)時,=0 s?1,=0。

        2 模型驗證

        利用顯式動力學(xué)軟件LS-DYNA進行有限元仿真計算。為保證計算結(jié)果的可靠性,首先從2個方面驗證本文有限元仿真計算方法的可靠性。

        2.1 與實驗結(jié)果比較

        采用上述材料模型,將同等條件下是否考慮該材料應(yīng)變率效應(yīng)的仿真結(jié)果與實驗結(jié)果進行對比。有限元仿真計算結(jié)構(gòu)與文獻[13]中沖擊實驗結(jié)構(gòu)一致,均是長度為500 mm、外徑為250 mm、壁厚為6 mm的圓筒,且在距底板30 mm處設(shè)置寬為10 mm、深為 3 mm的初始缺陷,沖擊質(zhì)量為1 030 kg,初始速度為17.8 m/s。有限元仿真計算模型如圖1所示,壓潰載荷計算結(jié)果對比如圖2所示。

        由圖2可知:在考慮Q235材料的應(yīng)變率強化效應(yīng)時,仿真結(jié)果與實驗結(jié)果的峰值壓潰載荷基本吻合,仿真結(jié)果的平均壓潰載荷大于實驗值,相對誤差為6.8%;而未考慮材料的應(yīng)變率強化效應(yīng)時,初始峰值載荷明顯小于實驗值,平均壓潰載荷偏小,壓潰位移偏大。經(jīng)對比發(fā)現(xiàn):在考慮材料的應(yīng)變率強化效應(yīng)后,仿真計算得到的峰值載荷以及載荷?位移曲線的變化趨勢與實驗結(jié)果更接近。由于有限元仿真計算中使用的材料本構(gòu)模型通過實驗結(jié)果擬合得到,與實際實驗中材料本身力學(xué)性能存在微小差異,故結(jié)構(gòu)變形與實際變形不完全相符。但由于相對誤差較小,本文使用的材料模型、有限元模型及計算方法均可以用于工程仿真研究。

        圖1 與實驗對比的有限元模型

        2.2 與理論解比較

        薄壁圓管軸向壓潰的塑性破損模式主要包括3種:圓環(huán)模式、鉆石模式和混合模式。國外對圓管壓潰變形的理論研究較多,故本文先以圓管壓潰結(jié)構(gòu)為研究對象,通過有限元仿真計算研究圓管壓潰結(jié)構(gòu)受材料應(yīng)變率強化效應(yīng)的影響,并與理論解進行比對。

        圓管外徑為150 mm,壁厚為3 mm,軸向長度為200 mm。使用的單元為Belytschko-Tsay薄殼單元,沿厚度方向取4個積分點,單元基本尺寸為2 mm,沖擊速度設(shè)定為恒速,依次為3,6,9,12和15 m/s。

        式中:為沖擊速度;m為圓管直徑。由此可以計算出在某一沖擊速度下的理論平均應(yīng)變率,并根據(jù)式(4)和仿真計算的未考慮材料應(yīng)變率效應(yīng)時的平均壓潰載荷計算考慮后結(jié)構(gòu)的理論動態(tài)平均壓潰載荷[20]。

        式中:d為動態(tài)平均壓潰載荷,即考慮了材料應(yīng)變率強化效應(yīng)后的平均壓潰載荷;s為靜態(tài)平均壓潰載荷,即不考慮材料應(yīng)變率效應(yīng)時的平均壓潰載荷。

        在計算平均壓潰載荷時,有效壓潰行程為吸能效率最大值所對應(yīng)的壓縮行程[21]。

        將動態(tài)平均壓潰載荷的理論計算值與仿真計算值進行對比,見表1。從表1可見:動態(tài)平均壓潰載荷理論計算結(jié)果與考慮材料應(yīng)變率效應(yīng)的仿真計算結(jié)果的最大相對誤差在10%以內(nèi),對工程研究來說精度 足夠。

        綜上可知:本文采用的材料模型、有限元模型及計算方法是可靠的,可以用于研究材料應(yīng)變率效應(yīng)對其他吸能結(jié)構(gòu)碰撞性能的影響。

        表1 平均壓潰載荷對比

        3 計算模型與邊界條件

        對于2種吸能結(jié)構(gòu),參照軌道車輛吸能結(jié)構(gòu)常用軸向長度,設(shè)計軸向長度為470 mm,發(fā)生塑性變形的薄壁厚度為3 mm,單元相關(guān)設(shè)置與模型驗證時保持一致。

        對于擠壓式吸能結(jié)構(gòu),小端外徑為135 mm,大端外徑為145 mm,通過圓弧光滑過渡;對于壓潰式圓錐管:小端外徑為105 mm,大端外徑為145 mm。2類吸能結(jié)構(gòu)模型如圖3所示。

        圖3 2類吸能結(jié)構(gòu)模型

        本文主要研究2類吸能結(jié)構(gòu)在軸向正面碰撞下的吸能特性。具體邊界條件設(shè)置為:法蘭座和安裝座固定不動,防爬齒部分以恒速沖擊擠壓管和壓潰管,沖擊速度分別取5,10和15 m/s,并且在每個模型中沿母線標記距離防爬齒板分別為75,250和370 mm的前、中、后部3個不同位置的單元用于碰撞性能分析。

        4 計算結(jié)果及分析

        4.1 碰撞載荷與吸能

        4.1.1 計算結(jié)果

        2種不同吸能結(jié)構(gòu)的碰撞力?位移曲線見圖4,其中擠壓式吸能結(jié)構(gòu)有效擠壓行程均按384 mm(防爬齒板接觸法蘭座)計算,壓潰式吸能結(jié)構(gòu)的有效壓潰行程取值方法與模型驗證時一致。

        從圖4可見:擠壓式吸能結(jié)構(gòu)的擠壓力?位移曲線比較平穩(wěn),僅存在微小波動;在材料應(yīng)變率強化效應(yīng)影響下,瞬時擠壓力與平均擠壓力整體提升,但碰撞載荷歷程未發(fā)生明顯改變。

        壓潰式圓錐管吸能結(jié)構(gòu)的壓潰力?位移曲線受材料應(yīng)變率強化效應(yīng)的影響比較顯著,在考慮材料應(yīng)變率強化效應(yīng)之后,結(jié)構(gòu)峰值壓潰力和平均壓潰力明顯增大,壓潰力歷程出現(xiàn)明顯差異,峰值和谷值的數(shù)量和位置發(fā)生改變,褶皺形成的時間和數(shù)量發(fā)生改變。

        (a) 擠壓式;(b) 壓潰式

        2種不同結(jié)構(gòu)在有效吸能過程中的總吸能見圖5。從圖5可知:在考慮材料的應(yīng)變率強化效應(yīng)后,擠壓式吸能結(jié)構(gòu)和壓潰式圓錐管結(jié)構(gòu)的總吸能均比未考慮時的大。

        將2種不同結(jié)構(gòu)在有效吸能過程中的平均碰撞力進行對比,發(fā)現(xiàn)在考慮材料的應(yīng)變率強化效應(yīng)之后,擠壓式吸能結(jié)構(gòu)的平均擠壓力和壓潰式圓錐管結(jié)構(gòu)的平均壓潰力也均比未考慮時的大。

        對于擠壓式吸能結(jié)構(gòu),不同沖擊速度下有效擠壓行程保持不變,而擠壓總吸能可表示為有效擠壓行程與平均擠壓力的乘積,故平均擠壓力與總吸能的放大因子相同。

        對于壓潰式圓錐管結(jié)構(gòu),有效壓潰行程因考慮材料應(yīng)變率強化效應(yīng)而增加,故總吸能的放大因子略大于平均壓潰力放大因子,壓潰式圓錐管結(jié)構(gòu)的平均壓潰力如表2所示。

        1—未考慮應(yīng)變率,擠壓式;2—考慮應(yīng)變率,擠壓式; 3—未考慮應(yīng)變率,壓潰式;4—考慮應(yīng)變率,壓潰式。

        表2 壓潰式圓錐管吸能結(jié)構(gòu)平均壓潰力統(tǒng)計

        4.1.2 討論與分析

        對不同速度下的總吸能放大因子和平均應(yīng)變率進行比較,其中平均應(yīng)變率取值為有效吸能過程中結(jié)構(gòu)母線上典型位置處單元應(yīng)變率的統(tǒng)計平均值,如圖6所示。從圖6可見:擠壓式吸能結(jié)構(gòu)的總吸能放大因子和平均應(yīng)變率隨沖擊速度的增加變化幅度很小,壓潰式吸能結(jié)構(gòu)比前者變化幅度更大;在速度以5 m/s的間隔變化時,壓潰式吸能結(jié)構(gòu)的應(yīng)變率以6 s?1的速度增加,擠壓式吸能結(jié)構(gòu)的增加速度只有1 s?1。為進一步確認擠壓式吸能結(jié)構(gòu)沖擊速度與總吸能放大因子的關(guān)系,對該結(jié)構(gòu)進行速度為30 m/s的碰撞仿真,此時,平均應(yīng)變率為6.65 s?1,放大因子為1.336,與其他速度相比明顯提高。仿真結(jié)果表明:隨速度同等增大,壓潰式圓錐管結(jié)構(gòu)比擠壓式結(jié)構(gòu)的平均應(yīng)變率和總吸能放大因子變化更明顯。

        1—放大因子,擠壓式;2—放大因子,壓潰式; 3—平均應(yīng)變率,擠壓式;4—平均應(yīng)變率,壓潰式。

        雖然擠壓式吸能結(jié)構(gòu)平均應(yīng)變率及其總吸能放大因子對速度變化不敏感,擠壓力歷程受材料應(yīng)變率強化效應(yīng)的影響也不明顯,但擠壓式吸能結(jié)構(gòu)因材料應(yīng)變率強化效應(yīng)帶來的平均擠壓力放大因子比壓潰式吸能結(jié)構(gòu)壓潰力放大因子更大。下面對2種結(jié)構(gòu)分別從碰撞力理論計算角度進行分析。

        對于壓潰式圓錐管結(jié)構(gòu),其在圓環(huán)變形、鉆石變形或者混合變形模式下,假設(shè)材料為理想剛塑性,平均壓潰力均可以表示成如下形式[22?23]:

        式中:為結(jié)構(gòu)平均壓潰力;為材料的靜態(tài)屈服應(yīng)力;為壁厚;和為與應(yīng)力無關(guān),與直徑、厚度和錐角相關(guān)的系數(shù),對于不同變形模式,該系數(shù)不同;為與結(jié)構(gòu)相關(guān)的相應(yīng)實數(shù)。式(5)與圓管壓潰平均壓潰力的表達形式一致,故在結(jié)構(gòu)有效壓潰行程不變情況下,結(jié)構(gòu)的平均壓潰力或總吸能的放大因子與材料流動應(yīng)力放大因子一致,動、靜態(tài)平均壓潰載荷之間均可通過式(4)計算平均壓潰力理論放大 因子。

        圓錐管的平均應(yīng)變率沒有理論計算公式,故取仿真結(jié)果的平均值,同時利用式(4)計算平均壓潰力放大因子并與仿真結(jié)果進行對比,見表3。從表3可見:圓錐管的平均壓潰力放大因子的理論計算結(jié)果與仿真計算結(jié)果基本吻合,最大相對誤差為1.0%,說明有限元仿真能準確地表征薄壁圓錐管在動態(tài)軸向壓縮時的碰撞性能。

        但對于擠壓式吸能結(jié)構(gòu),材料應(yīng)變率效應(yīng)使結(jié)構(gòu)平均擠壓力放大因子比材料自身流動應(yīng)力放大因子大得多,見表4。

        表3 壓潰式圓錐管理論與仿真平均壓潰力放大因子對比

        表4 結(jié)構(gòu)平均擠壓力與材料流動應(yīng)力放大因子對比

        擠壓式吸能結(jié)構(gòu)可視為向內(nèi)膨脹式吸能結(jié)構(gòu)[24],擠壓式和膨脹式吸能結(jié)構(gòu)的變形過程均包括彎曲、縮徑(擴徑)、摩擦3部分。緩沖力模型不能用類似于式(5)的形式表達[25],與壓潰式結(jié)構(gòu)由材料應(yīng)變率效應(yīng)引起的碰撞力放大因子存在差異。

        與壓潰式吸能結(jié)構(gòu)僅依靠塑性變形吸能不同,擠壓式吸能結(jié)構(gòu)能量吸收由塑性變形吸能與摩擦耗能2部分組成。不同吸能的放大因子對比見圖7。從圖7可見:考慮了材料的應(yīng)變率強化效應(yīng)時,摩擦耗能的放大因子大于塑性變形吸能的放大因子,總吸能放大因子提高。

        1—總吸能;2—塑性變形吸能;3—摩擦耗能。

        綜上所述,由于存在擠壓式吸能結(jié)構(gòu)特殊的變形模式,同時考慮擠壓管與法蘭座之間的摩擦耗能,故在考慮材料應(yīng)變率強化效應(yīng)后結(jié)構(gòu)的平均擠壓力放大因子大于材料流動應(yīng)力放大因子。但以上僅為定性分析結(jié)果,針對擠壓式吸能結(jié)構(gòu)平均擠壓力受材料應(yīng)變率強化效應(yīng)的影響有待深入研究。

        以上對比結(jié)果表明考慮材料應(yīng)變率強化效應(yīng)后結(jié)構(gòu)的平均壓潰力或擠壓力放大因子與結(jié)構(gòu)的變形模式密切相關(guān)。在變形過程中,應(yīng)變較大的結(jié)構(gòu)通常伴隨著較高的應(yīng)變率,如壓潰式結(jié)構(gòu)等第二類吸能結(jié)構(gòu)對速度的敏感性較強,變形應(yīng)變率及壓潰力或總吸能放大因子隨速度改變發(fā)生明顯變化;而對于應(yīng)變較小、變形穩(wěn)定的擠壓式等第一類吸能結(jié)構(gòu),應(yīng)變率通常較小,應(yīng)變率及相應(yīng)的放大因子對速度變化不敏感。但考慮材料應(yīng)變率強化效應(yīng)引起的結(jié)構(gòu)吸能或平均碰撞力的放大因子與結(jié)構(gòu)速度敏感性不成正相關(guān),而是取決于結(jié)構(gòu)變形機理及吸能方式。

        4.2 初期變形模式

        4.2.1 計算結(jié)果

        考慮材料的應(yīng)變率強化效應(yīng)會對結(jié)構(gòu)的初始變形產(chǎn)生影響。無論是擠壓式吸能結(jié)構(gòu)還是壓潰式圓錐管吸能結(jié)構(gòu),在考慮材料應(yīng)變率強化效應(yīng)之后,變形初期產(chǎn)生塑性應(yīng)變的單元數(shù)目均會增多,如圖8和圖9所示。

        4.2.2 討論與分析

        對于擠壓式吸能結(jié)構(gòu),在考慮材料的應(yīng)變率強化效應(yīng)時,擠壓后端壓潰鼓脹產(chǎn)生的塑性應(yīng)變比不考慮時的小,但產(chǎn)生塑性應(yīng)變的單元數(shù)目增多,見圖8所示的擠壓后端應(yīng)變云圖。由于壓潰鼓脹變形應(yīng)變率較大,材料的流動應(yīng)力增加明顯,增強了端部的抵抗能力,使壓潰鼓脹變形減小,軸向壓縮變形增加。同理,前端擠壓變形包括彎曲、軸向壓縮以及徑向收縮,彎曲變形應(yīng)變率最大,故與擠壓后端一致,彎曲變形減小,徑向收縮變形與軸向壓縮變形增大,總體表現(xiàn)為單元塑性應(yīng)變減小但產(chǎn)生塑性應(yīng)變的單元數(shù)目 增多。

        對于壓潰式圓錐管結(jié)構(gòu),在3種速度下,未考慮材料應(yīng)變率強化效應(yīng)時,初始時刻均已形成塑性鉸,如圖9所示。該部分應(yīng)變和應(yīng)變率較大,在考慮材料應(yīng)變率強化效應(yīng)后繞塑性鉸轉(zhuǎn)動變形所需要的力明顯增大,所以,在同樣的軸向力作用下,繞塑性鉸的轉(zhuǎn)動變形一部分轉(zhuǎn)變?yōu)閼?yīng)力增加較小的軸向壓縮變形,表現(xiàn)為未形成初始塑性鉸或形成滯后,應(yīng)變向后傳遞,在應(yīng)變云圖上表現(xiàn)為產(chǎn)生初始塑性應(yīng)變的單元數(shù)目增多。

        (a) 10 m/s,未考慮應(yīng)變率;(b) 10 m/s,考慮應(yīng)變率; (c) 15 m/s,未考慮應(yīng)變率;(d) 15 m/s,考慮應(yīng)變率

        (a) 5 m/s,未考慮應(yīng)變率;(b) 5 m/s,考慮應(yīng)變率; (c) 10 m/s,未考慮應(yīng)變率;(d) 10 m/s,考慮應(yīng)變率; (e) 15 m/s,未考慮應(yīng)變率;(f)15 m/s,考慮應(yīng)變率

        壓潰式圓錐管結(jié)構(gòu)前部、中部和后部3單元的應(yīng)變時間歷程見圖10。從圖10可見:每處單元產(chǎn)生塑性鉸之前即應(yīng)變明顯增大之前,考慮材料應(yīng)變率強化效應(yīng)時單元的應(yīng)變都比未考慮時的更大,增加的應(yīng)變即為結(jié)構(gòu)中軸向壓縮變形應(yīng)變的增加值。材料的應(yīng)變率強化效應(yīng)減少了繞塑性鉸轉(zhuǎn)動變形的應(yīng)變而增加了軸向壓縮變形的應(yīng)變,這與前面的分析結(jié)果一致。

        速度/(m·s?1):(a) 5;(b) 10;(c) 15

        4.3 整體變形模式

        4.3.1 計算結(jié)果

        從整體變形模式分析,擠壓式吸能結(jié)構(gòu)變形模式比較穩(wěn)定,是否考慮材料的應(yīng)變率效應(yīng)對結(jié)構(gòu)變形模式影響較小。但隨著沖擊速度增加,在吸能結(jié)構(gòu)遠離法蘭座的擠壓后端會出現(xiàn)圖11中類似于壓潰結(jié)構(gòu)的壓潰鼓脹變形,而且速度越大,壓潰鼓脹變形越明顯。

        對于壓潰式圓錐管結(jié)構(gòu),在不同的沖擊速度下,是否考慮材料應(yīng)變率強化效應(yīng)對變形模式影響顯著,如圖12所示。

        圖11 速度為15 m/s時后端鼓脹變形

        (a) 5 m/s,未考慮應(yīng)變率;(b) 5 m/s,考慮應(yīng)變率; (c) 10 m/s,未考慮應(yīng)變率;(d) 10 m/s,考慮應(yīng)變率; (e) 15 m/s,未考慮應(yīng)變率;(f) 15 m/s,考慮應(yīng)變率

        4.3.2 討論與分析

        對于擠壓式吸能結(jié)構(gòu),由于隨速度增加,單位時間內(nèi)變形增多,靠近法蘭座的擠壓前端的擠壓變形所需要的擠壓力大于擠壓后端壓潰鼓脹變形所需要的壓潰力,故表現(xiàn)為擠壓前端變?yōu)樗矔r固定端,結(jié)構(gòu)產(chǎn)生擠壓后端瞬時壓潰鼓脹變形。而對于繞塑性鉸轉(zhuǎn)動的壓潰鼓脹變形,在褶皺完全形成之前,隨著變形增大,壓潰力迅速增大;當(dāng)壓潰力大于擠壓前端的擠壓力時,開始發(fā)生擠壓前端的擠壓變形,并一直以擠壓變形到碰撞結(jié)束。

        對于壓潰式圓錐管結(jié)構(gòu),在5 m/s與10 m/s的沖擊速度下,變形模式因考慮材料應(yīng)變率強化效應(yīng)由有序變?yōu)闊o序,單個塑性鉸的形成所需要的軸向高度更大,結(jié)構(gòu)出現(xiàn)失穩(wěn)而未發(fā)生屈曲吸能的材料增加。由于考慮了材料的應(yīng)變率強化效應(yīng),材料變形過程中應(yīng)變率較大的部分材料的流動應(yīng)力明顯增大。薄壁圓錐管的軸向壓潰變形過程可以用折板的動力學(xué)行為進行解釋,主要包括應(yīng)變率較小的軸向壓縮與應(yīng)變率較大的繞塑性鉸的轉(zhuǎn)動。故在考慮材料的應(yīng)變率強化效應(yīng)后,繞塑性鉸的轉(zhuǎn)動變形由于應(yīng)變率較高,材料流動應(yīng)力增加顯著,所需要的轉(zhuǎn)動力矩增大,較大的力矩會繼續(xù)向結(jié)構(gòu)后部傳遞,在宏觀上表現(xiàn)為塑性鉸的內(nèi)折點更接近軸線,即由移動至,圓錐管的直徑內(nèi)縮現(xiàn)象更明顯,如圖13所示。對于15 m/s的沖擊速度,考慮材料的應(yīng)變率強化效應(yīng)時,褶皺形成的時間和位置發(fā)生明顯改變,也缺失了第一個階段的圓環(huán)變形模式,如圖12(e)和(f)所示。

        圖13 速度為10 m/s時考慮應(yīng)變率效應(yīng)的結(jié)構(gòu)變形

        對比2種不同類型的吸能結(jié)構(gòu)可知:圓錐管壓潰式等第二類吸能結(jié)構(gòu)的變形模式受材料應(yīng)變率強化效應(yīng)的影響更顯著;擠壓式吸能結(jié)構(gòu)等第一類吸能結(jié)構(gòu)的變形模式與是否考慮材料應(yīng)變率效應(yīng)關(guān)系不大。在不同吸能結(jié)構(gòu)的變形中,變形形式多樣,包括彎曲、壓縮、拉伸等,各種變形的應(yīng)變率不同。因考慮材料應(yīng)變率強化效應(yīng),材料流動應(yīng)力增大,使結(jié)構(gòu)應(yīng)變率大的變形減少,應(yīng)變率小的變形增多或者提前發(fā)生。

        5 誘導(dǎo)式吸能結(jié)構(gòu)

        在工程應(yīng)用中期望吸能結(jié)構(gòu)具有可控的變形模式。對于壓潰式吸能結(jié)構(gòu),實際碰撞中材料的應(yīng)變率效應(yīng)會使結(jié)構(gòu)變形不可控,需要采取一定的誘導(dǎo)措施使材料應(yīng)變率強化效應(yīng)增加的應(yīng)力弱于結(jié)構(gòu)自身的誘導(dǎo)能力,即結(jié)構(gòu)仍然以有序可控的模式變形。故針對前面所研究圓錐管壓潰變形,提出2種不同的誘導(dǎo)方式,并研究在誘導(dǎo)結(jié)構(gòu)下變形模式受材料應(yīng)變率強化效應(yīng)的影響。

        2種誘導(dǎo)式吸能結(jié)構(gòu)見圖14。

        誘導(dǎo)方式1:圓錐管周向間隔90°開設(shè)直徑為 20 mm的圓形誘導(dǎo)孔,軸向間距為50 mm。

        誘導(dǎo)方式2:錐管管壁預(yù)彎變形。

        由于速度越大,材料的應(yīng)變率效應(yīng)越明顯,變形模式的差異越明顯,故將15 m/s沖擊速度下的2種誘導(dǎo)方式及原有結(jié)構(gòu)的計算結(jié)果進行對比。

        (a) 誘導(dǎo)方式1 ;(b) 誘導(dǎo)方式2

        不同誘導(dǎo)方式下壓潰力?位移曲線見圖15。從圖15可見:在誘導(dǎo)方式下,壓潰力?位移曲線近似向上平移,明顯比無誘導(dǎo)結(jié)構(gòu)圓錐管的差異減小,褶皺的數(shù)目和出現(xiàn)的位置基本一致。

        1—誘導(dǎo)孔,未考慮應(yīng)變率;2—誘導(dǎo)孔,考慮應(yīng)變率; 3—預(yù)彎,未考慮應(yīng)變率;4—預(yù)彎,考慮應(yīng)變率。

        在沖擊速度為15 m/s時,不同結(jié)構(gòu)變形模式見圖16。從圖16可見:在同一時刻,誘導(dǎo)式吸能結(jié)構(gòu)的變形模式受材料應(yīng)變率強化效應(yīng)的影響比較小。開誘導(dǎo)孔與預(yù)彎變形均相當(dāng)于人為地約束了褶皺形成的數(shù)量和位置。雖然材料應(yīng)變率效應(yīng)會使流動應(yīng)力發(fā)生變化,但結(jié)構(gòu)自身的剛度差值大于應(yīng)變率效應(yīng)引起的應(yīng)力差值,結(jié)構(gòu)自身的誘導(dǎo)起主導(dǎo)作用,導(dǎo)致結(jié)構(gòu)變形。對于預(yù)彎變形的誘導(dǎo)方式,其變形模式受材料應(yīng)變率效應(yīng)的影響最小,因為其誘導(dǎo)方式迎合了結(jié)構(gòu)的變形模式,完成了結(jié)構(gòu)的前期變形,其壓潰力?位移曲線表明其已經(jīng)近似為第一類吸能結(jié)構(gòu)。

        (a) 無誘導(dǎo)吸能結(jié)構(gòu),15 m/s,未考慮應(yīng)變率; (b) 無誘導(dǎo)吸能結(jié)構(gòu),15 m/s,考慮應(yīng)變率; (c) 誘導(dǎo)孔式吸能結(jié)構(gòu),15 m/s,未考慮應(yīng)變率; (d) 誘導(dǎo)孔式吸能結(jié)構(gòu),15 m/s,考慮應(yīng)變率; (e) 預(yù)彎式吸能結(jié)構(gòu),15 m/s,未考慮應(yīng)變率; (f) 預(yù)彎式吸能結(jié)構(gòu),15 m/s,考慮應(yīng)變率

        誘導(dǎo)式吸能結(jié)構(gòu)可以有效減小材料應(yīng)變率強化效應(yīng)對結(jié)構(gòu)變形模式的影響,但誘導(dǎo)結(jié)構(gòu)一般是通過預(yù)變形或者開誘導(dǎo)孔的方式,預(yù)變形因為已有的塑性變形減小了吸能能力,誘導(dǎo)孔則因為材料減少削弱了吸能能力,在設(shè)計中需要綜合考慮,對結(jié)構(gòu)進行優(yōu)化 設(shè)計。

        6 結(jié)論

        1) 材料的應(yīng)變率強化效應(yīng)會提高結(jié)構(gòu)的承載能力與吸能能力,但不同結(jié)構(gòu)的影響因子不同,這取決于結(jié)構(gòu)的變形模式。材料應(yīng)變率強化效應(yīng)對擠壓式吸能結(jié)構(gòu)碰撞載荷和吸能的放大因子比壓潰式圓錐管吸能結(jié)構(gòu)的大,工程中不可忽略該類結(jié)構(gòu)受材料應(yīng)變率效應(yīng)的影響。

        2) 隨著速度增加,壓潰式圓錐管吸能結(jié)構(gòu)平均應(yīng)變率和應(yīng)變率效應(yīng)放大因子的變化速度比擠壓式吸能結(jié)構(gòu)的大,即壓潰式圓錐管吸能結(jié)構(gòu)平均應(yīng)變率和應(yīng)變率效應(yīng)放大因子對速度變化更敏感。

        3) 結(jié)構(gòu)的完整變形和初期變形受材料應(yīng)變率效應(yīng)的影響與結(jié)構(gòu)類型相關(guān),對變形中應(yīng)變率差異不明顯的第一類吸能結(jié)構(gòu)的變形模式影響較小,對應(yīng)變率差異比較大的第二類吸能結(jié)構(gòu)的變形影響顯著。材料的應(yīng)變率強化效應(yīng)會減小應(yīng)變率大的變形而增加或促進應(yīng)變率小的變形。

        4) 通過誘導(dǎo)方式可以有效減少材料應(yīng)變率效應(yīng)對結(jié)構(gòu)變形模式的影響,提高結(jié)構(gòu)變形的有序性和可控性,但在一定程度上減小了結(jié)構(gòu)的吸能能力,在設(shè)計時需綜合考慮。

        [1] 姚松, 田紅旗. 車輛吸能部件的薄壁結(jié)構(gòu)碰撞研究[J]. 中國鐵道科學(xué), 2001, 22(2): 55?60. YAO Song, TIAN Hongqi. Crash research on thin-shelled structure as vehicle energy-absorbing components[J]. China Railway Science, 2001, 22(2): 55?60.

        [2] WOLTER W. 鐵道車輛的防碰撞要求、設(shè)計原理和初步結(jié)果[J]. 閻鋒, 譯. 國外鐵道車輛, 2004, 41(2): 23?30. WOLTER W. Requirements, design principles and initial results of crashworthy rail vehicles[J]. YAN Feng, trans. Foreign Rolling Stock, 2004, 41(2): 23?30.

        [3] KOO J S, KWON T S, CHO H J. 韓國高速列車防撞設(shè)計與評估[J]. 周穎, 周艷芳, 譯. 中國鐵道科學(xué), 2004, 25(1): 1?7. KOO J S, KWON T S, CHO H J. Korean high-speed railway anti-impact structure design and appraisal[J]. ZHOU Yin, ZHOU Yanfang, trans. China Railway Science, 2004, 25(1): 1?7.

        [4] 謝素超, 田紅旗, 周輝. 耐沖擊地鐵車輛設(shè)計及整車碰撞研究[J]. 鐵道科學(xué)與工程學(xué)報, 2008, 5(5): 65?70. XIE Suchao, TIAN Hongqi, ZHOU Hui. The design of crashworthy subway vehicle and crash research of whole car-body[J]. Journal of Railway Science and Engineering, 2008, 5(5): 65?70.

        [5] 李健, 高廣軍, 董海鵬, 等. 帶隔板薄壁方管的耐撞性研究[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2014, 45(7): 2481?2488. LI Jian, GAO Guangjun, DONG Haipeng, et al. Research on crashworthiness of thin-walled square tubes with diaphragms[J]. Journal of Central South University(Science and Technology), 2014, 45(7): 2481?2488.

        [6] XIE Suchao, ZHOU Hui, LIANG Xifeng, et al. Contrastive analysis and crashworthiness optimization of two composite thin-walled structures[J]. Journal of Central South University, 2014, 21(11): 4386?4394.

        [7] 朱濤, 肖守訥, 楊超, 等. 機車車輛被動安全性研究綜述[J]. 鐵道學(xué)報, 2017, 39(5): 22?32. ZHU Tao, XIAO Shoune, YANG Chao, et al. State-of-the-art development of passive safety of rolling stocks[J]. Journal of the China Railway Society, 2017, 39(5): 22?32.

        [8] 廖其紅, 黃宏成, 薛量. 薄壁鋼結(jié)構(gòu)碰撞性能仿真中材料模型研究[J]. 機械設(shè)計與研究, 2001, 17(4): 9, 64?66. LIAO Qihong, HUANG Hongcheng, XUE Liang. On material model for thin-walled sheet metal crashworthiness simulation[J]. Machine Design and Research, 2001, 17(4): 9, 64?66.

        [9] 孟廣成, 邰永剛. 材料應(yīng)變率對護欄碰撞仿真結(jié)果影響研究[J]. 公路交通科技(應(yīng)用技術(shù)版), 2012(6): 401?403. MENG Guangcheng, TAI Yonggang. Study on the effect of material strain rate on the simulation result of guardrail collision[J]. Journal of Highway and Transportation Research and Development(Application Technology), 2012(6): 401?403.

        [10] 王自力, 顧永寧. 應(yīng)變率敏感性對船體結(jié)構(gòu)碰撞性能的影響[J]. 上海交通大學(xué)學(xué)報, 2000, 34(12): 1704?1707. WANG Zili, GU Yongning. Effect of strain-rate sensitivity on the behavior of ship structure in collision[J]. Journal of Shanghai Jiaotong University, 2000, 34(12): 1704?1707.

        [11] 曾必強, 胡遠志, 謝書港. 材料應(yīng)變率強化效應(yīng)對結(jié)構(gòu)碰撞響應(yīng)的影響[C]//第八屆國際汽車交通安全學(xué)術(shù)會議論文集. 蕪湖, 2010: 134?140.ZENG Biqiang, HU Yuanzhi, XIE Shugang. The influences of material strain-rate strengthen effect on the response of structure collision[C]//Proceedings of the 8th International Conference on Vehicle Traffic Safety. Wuhu, China, 2010: 134?140.

        [12] 張永康, 張紅兵. 不同構(gòu)型吸能元件的軸向抗撞擊性能分析[J]. 機械設(shè)計與制造, 2011(11): 213?215. ZHANG Yongkang, ZHANG Hongbing. Analysis of axial anti-impact performance for different configuration power-absorbing[J]. Machinery Design & Manufacture, 2011(11): 213?215.

        [13] 楊喆, 魏延鵬, 鄒金龍, 等. 薄壁圓筒沖擊屈曲吸能行為的實驗與數(shù)值模擬研究[J]. 兵工學(xué)報, 2014, 35(S2): 244?250. YANG Zhe, WEI Yanpeng, ZOU Jinlong, et al. Experimental and numerical research on the energy absorption behavior of thin-walled tube under impact load[J]. Acta Armamentarii, 2014, 35(S2): 244?250.

        [14] 楊超, 朱濤, 肖守訥. 列車車體鋁合金動態(tài)力學(xué)性能及其對吸能的影響[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2015, 46(7): 2744?2749. YANG Chao, ZHU Tao, XIAO Shoune. Dynamic mechanical properties of aluminum alloy used in car bodies of trains and effect on energy absorption[J]. Journal of Central South University(Science and Technology), 2015, 46(7): 2744?2749.

        [15] CALLADINE C R, ENGLISH R W. Strain?rate and inertia effects in the collapse of two types of energy?absorbing structure[J]. International Journal of Mechanical Sciences, 1984, 26(11): 689?701.

        [16] 余同希, 盧國興. 材料與結(jié)構(gòu)的能量吸收[M]. 華云龍, 譯. 北京: 化學(xué)工業(yè)出版社, 2006: 149?152. YU Tongxi, LU Guoxing. Energy absorption of structures and materials[M]. HUA Yunlong, trans. Beijing: Chemical Industry Press, 2006: 149?152.

        [17] 許平, 邵恒, 嚴佳麗. 基于自適應(yīng)響應(yīng)面法的膨脹式吸能結(jié)構(gòu)耐撞性優(yōu)化設(shè)計[J]. 振動與沖擊, 2017, 36(11): 118?123, 146. XU Ping, SHAO Heng, YAN Jiali. Crashworthiness optimization design of expanding type energy absorption devices based on adaptive response surface method[J]. Journal of Vibration and Shock, 2017, 36(11): 118?123, 146.

        [18] 郝林坡, 羅云華, 熊佳霖, 等. 翻卷管吸能啟動載荷的研究[J]. 鍛壓技術(shù), 2015, 40(12): 166?168, 175. HAO Linpo, LUO Yunhua, XIONG Jialin, et al. Research on starting load of inversion tube energy absorber[J]. Forging & Stamping Technology, 2015, 40(12): 166?168, 175.

        [19] 洪武, 徐迎, 金豐年, 等. 薄壁圓錐管軸向壓縮吸能特性研究[J]. 振動與沖擊, 2015, 34(5): 88?94. HONG Wu, XU Ying, JIN Fengnian, et al. Energy absorbing characteristics of tapered circular tubes under axial compression[J]. Journal of Vibration and Shock, 2015, 34(5): 88?94.

        [20] ABRAMOWICZ W, JONES N. Dynamic axial crushing of circular tubes[J]. International Journal of Impact Engineering, 1984, 2(3): 263?281.

        [21] 項燕飛, 楊黎明. 不同管狀和材料對結(jié)構(gòu)能量吸收性能的影響[J]. 寧波大學(xué)學(xué)報(理工版), 2015, 28(2): 85?90. XIANG Yanfei, YANG Liming. Effects of section and material of tubes on energy absorption[J]. Journal of Ningbo University(Natural Science & Engineering Edition), 2015, 28(2): 85?90.

        [22] MAMALIS A G, MANOLAKOS D E, VIEGELAHN G L, et al. The modeling of the progressive extensible plastic collapse of thin-wall shells[J]. International Journal of Mechanical Sciences, 1988, 30(3): 249?261.

        [23] MAMALIS A G, MANOLAKOS D E, SAIGAL S, et al. Extensible plastic collapse of thin-wall frusta as energy absorbers[J]. International Journal of Mechanical Sciences, 1986, 28(4): 219?229.

        [24] ALMEIDA B P P, ALVES M L, ROSA P A R, et al. Expansion and reduction of thin-walled tubes using a die: experimental and theoretical investigation[J]. International Journal of Machine Tools & Manufacture, 2006, 46(12/13): 1643?1652.

        [25] 吳鴻超, 梁增友, 馮陽, 等. 薄壁金屬管在中高速沖擊下的緩沖特性研究[J]. 應(yīng)用力學(xué)學(xué)報, 2016, 33(2): 325?331, 377?378. WU Hongchao, LIANG Zengyou, FENG Yang, et al. Study about the cushioning properties of the thin-walled metal tube which working in high-speed impact environment[J]. Chinese Journal of Applied Mechanics, 2016, 33(2): 325?331, 377?378.

        (編輯 陳燦華)

        Influence of material strain rate effect on crashworthiness of two kinds of energy-absorbing structure

        FENG Yue, XIAO Shoune, ZHU Tao, YANG Bing, YANG Haojie

        (State Key Laboratory of Traction Power, Southwest Jiaotong University, Chengdu 610031, China)

        The deformation of energy-absorbing structure on rail vehicles is accompanied by obvious strain rate effect. In order to study the influence of material strain rate effect on collision simulation results, the material Q235 which has obvious strain rate effect was studied. Firstly, the reliability of Q235 material model considering strain rate effect was proved by comparison with experimental results and theoretical results. Then the influence of material strain rate effect on the impact load and energy absorption, deformation modes of squeezed and crushed cone energy-absorbing structures were studied. On the basis of this, two kinds of induced energy-absorbing structures were proposed, which were almost not affected by material strain rate effect and the deformation was orderly and controllable. The results show that the influence factor of material strain rate effect varies with different structures. And the influence factor of squeezed structure is larger than crushed cone energy-absorbing structure considering the material strain rate enhancement effect. For different structures, the influence factor of material strain rate effect changes with the change of velocity. The influence factor of the squeezed structure does not increase as greatly as crushed cone energy-absorbing structure when velocity increases. The material strain rate effect has significant influence on the deformation mode of the second type of energy-absorbing structure, and the strain rate enhancement effect reduces the deformation with high strain rate and increases the deformation with low strain rate. The deformation mode of induced energy-absorbing structure is almost not affected by material strain rate effect, and the deformation is orderly and controllable. It can be considered in the engineering design to reduce the influence of material strain rate on the deformation mode.

        rail vehicles; crash worthiness; energy-absorbing structure; strain rate effect

        10.11817/j.issn.1672?7207.2018.10.032

        O313.4;U270.2

        A

        1672?7207(2018)10?2625?11

        2017?11?12;

        2018?01?18

        國家自然科學(xué)基金資助項目(51675446);國家重點研發(fā)計劃項目(2016YFB1200403)(Project(51675446) supported by the National Natural Science Foundation of China; Project(2016YFB1200403) supported by the National Key Research and Development Program of China)

        朱濤,博士,副研究員,從事機車車輛設(shè)計與理論研究;E-mail:zhutao034@swjtu.cn

        猜你喜歡
        圓錐效應(yīng)誘導(dǎo)
        圓錐擺模型的探究與拓展
        鈾對大型溞的急性毒性效應(yīng)
        齊次核誘導(dǎo)的p進制積分算子及其應(yīng)用
        圓錐截線與玫瑰線
        “圓柱與圓錐”復(fù)習(xí)指導(dǎo)
        計算法在圓錐保持架收縮模組合沖頭設(shè)計中的應(yīng)用
        哈爾濱軸承(2021年4期)2021-03-08 01:00:50
        同角三角函數(shù)關(guān)系及誘導(dǎo)公式
        懶馬效應(yīng)
        續(xù)斷水提液誘導(dǎo)HeLa細胞的凋亡
        中成藥(2017年12期)2018-01-19 02:06:52
        大型誘導(dǎo)標在隧道夜間照明中的應(yīng)用
        无遮高潮国产免费观看韩国| 色偷偷88888欧美精品久久久| 乱中年女人伦av三区| 日本一区二区三区在线视频观看 | 国产视频网站一区二区三区| 一区二区三区成人av| 无套无码孕妇啪啪| 国产97色在线 | 亚洲| 亚洲AV无码精品一区二区三区l| 亚洲一区二区三区精品久久| 一本色道久久hezyo无码| 亚洲一线二线三线写真| 精品国产品欧美日产在线| 国产洗浴会所三级av| 国内精品视频一区二区三区八戒| 成人精品综合免费视频| 红杏性无码免费专区| 国产性感午夜天堂av| 国产伦人人人人人人性| 亚洲欧美日韩精品高清| 正在播放东北夫妻内射| 国产亚洲欧美日韩综合综合二区| 亚洲处破女av一区二区| 国产亚洲精品综合一区| 国产自偷自偷免费一区| 女同性恋精品一区二区三区| av在线一区二区精品| 日本另类αv欧美另类aⅴ| 国产免费破外女真实出血视频| 久久国产高潮流白浆免费观看| 手机在线免费av资源网| 亚洲精品~无码抽插| 久久精品日韩av无码| 少妇被日到高潮的视频| 久久熟妇少妇亚洲精品| 成人片黄网站色大片免费观看cn | 在线看不卡的国产视频| 国产精品毛片无遮挡高清| 色偷偷久久一区二区三区| 人人爽亚洲aⅴ人人爽av人人片| 中文字幕高清不卡视频二区|