汽車碰撞試驗中碰撞過程的運動及力學性能

王家順, 安 琦

(華東理工大學機械與動力工程學院,上海 200237)

汽車碰撞試驗中碰撞過程的運動及力學性能

王家順, 安 琦

(華東理工大學機械與動力工程學院,上海 200237)

以汽車碰撞試驗為研究對象,分析歸納了摩擦力做功及阻擋桿變形能的計算方法,根據能量守恒的原則,給出了碰撞過程的能量平衡模型,構建了對碰撞過程中滑車的運動以及阻擋桿的變形進行計算的數學模型,并提出了一種利用計算機對該模型進行計算求解的計算方法。以Matlab為計算工具,編制了計算程序,結合一個具體的算例,研究滑車質量、初始碰撞速度、內側柱銷中心距、阻擋桿截面尺寸等因素對滑車位移、沖撞力、速度、以及減速度的影響規律,得到了相應的變化曲線,并對其進行了分析,得出了有益的結論。

碰撞試驗; 減速裝置; 力學建模; 運動性能; 力學性能

在汽車安全裝置的開發測試過程中需要進行大量的汽車碰撞試驗,主要包括實車碰撞試驗和滑車碰撞試驗[1]?；嚺鲎苍囼灳哂性囼炠M用低、準備周期短等優點,故應用較為廣泛。

碰撞減速裝置是該類試驗的核心部件,其原理是通過滑車的碰撞頭沖撞阻擋桿,使其變形,從而形成迅速停車效應,并且滑車的減速度波形需滿足規范要求。實際應用中,由于缺乏理論指導,設計人員往往需要進行大量的試驗,才能獲得所需要的減速度波形,不僅效率低,而且難以把握。

對于汽車碰撞試驗中的碰撞問題已有很多研究。胡玉梅等[1]對減速裝置中鋼板的排列與減速度波形的關系進行了研究,并總結出該類裝置鋼板的布置方法;曹立波等[2-3]利用LS-DYNA研究碰撞過程中鋼板變形,得出截面尺寸等參數對滑車減速度的影響; Massenzio等[4]通過建立減速裝置的力學模型來求解碰撞過程滑車減速度的大小,但是該模型僅考慮了鋼板塑性變形的影響,模型過于簡化;Giavotto等[5]利用試驗測試和有限元仿真的方法驗證了利用鋼條塑性變形來吸收滑車動能的可行性,并建立數學模型,利用遺傳算法研究如何利用鋼條布置來獲得期望的減速度波形;葛如海等[6]通過仿真研究,分析了鋼筋類型等參數對減速過程的影響;張慧云等[7]根據試驗經驗,總結出復現不同減速度波形時鋼筋的擺放方式;Xu等[8-9]對薄壁管結構的碰撞性能進行了仿真研究,并得出不同參數對該結構碰撞性能的影響;Jimenez等[10]仿真研究蜂窩鋁機械塊的吸能效果,結果顯示該裝置對改善液壓型減速裝置的初始減速度波形有較好效果。這些研究主要是通過有限元軟件對碰撞過程中的減速裝置進行研究,對碰撞過程中滑車的減速度變化等缺少定量分析。

本文以碰撞試驗臺為研究對象,對碰撞過程的運動和阻擋桿的變形進行力學分析,建立了相關力學模型,通過數值計算,研究了滑車質量、初始碰撞速度、內側柱銷中心距、阻擋桿的截面尺寸等因素對碰撞過程中滑車減速度、停車位移等參數的影響。

1 碰撞過程的力學建模

1.1 汽車碰撞試驗原理

利用彈性繩作為彈射介質的汽車碰撞試驗臺模型如圖1所示,其工作原理為:卷揚機拖動滑車沿軌道向后運動,同時滑車拉伸彈性繩進行儲能。當滑車運動至釋放位置時,由釋放裝置釋放滑車,滑車在彈性繩的拉力下加速運動,直至彈性繩恢復原長,隨后滑車與減速裝置進行碰撞,通過檢測試驗假人在減速過程中各部位的傷害值來評價汽車安全裝置對汽車乘員的保護性能。

碰撞過程中,滑車的碰撞桿頭部與阻擋桿接觸,并且推動阻擋桿向前運動,此時,阻擋桿在柱銷處以及碰撞桿頭部倒角處被迫發生塑性彎曲變形。碰撞過程中不同時刻的碰撞狀態如圖2所示。

圖1 碰撞試驗臺模型Fig.1 Model of collision test bench

圖2 不同時刻的碰撞狀態Fig.2 Collision state at different time

1.2 碰撞過程的力學分析

碰撞過程中,滑車動能迅速減小,轉變成其他形式的能量而被逐漸耗散。為了便于分析,作如下假設:

(1) 忽略發熱等因素的影響;

(2) 滑車車輪與軌道之間的滾動摩擦阻力因子、阻擋桿與柱銷之間的滑動摩擦因子為定值;

(3) 碰撞過程中,碰撞桿頭部與阻擋桿始終保持接觸;

(4) 碰撞結束時,滑車無回彈運動;

(5) 碰撞過程阻擋桿只發生塑性變形,忽略其在變形過程中的伸長量;

(6) 忽略應變速率對材料性能的影響,塑性變形過程中應力-應變之間的本構關系簡化為線性強化剛塑性模型。

碰撞過程中,滑車水平方向受力分析如圖3所示。m為滑車質量,F′為阻擋桿對滑車上碰撞桿施加的力,在F′作用下滑車迅速減速,直至停止運動;Ff為滑車車輪與軌道之間的滾動摩擦力,運動過程中,車輪與軌道之間摩擦力做功為Wf1。

圖3 滑車受力分析Fig.3 Mechanics analysis of sled

碰撞過程中,阻擋桿受力分析如圖4所示。其中,F為碰撞桿對阻擋桿施加的力,FN1、FN2、FN3、FN4為阻擋桿在柱銷處的支反力,Ff1、Ff2、Ff3、Ff4為阻擋桿與柱銷之間的摩擦力。

將減速裝置看作近似完全對稱,計算得

(1)

(2)

運動過程中,摩擦力Ff1、Ff2、Ff3、Ff4做功用Wf2表示。

圖4 阻擋桿受力分析Fig.4 Mechanics analysis of block bar

對于滑車和阻擋桿組成的系統,F和F′是一對作用力和反作用力。根據假設分析可知,碰撞過程中,滑車的動能主要通過以下3種方式耗散:

(1) 滑車與軌道之間的摩擦力做功Wf1;

(2) 阻擋桿與柱銷之間的摩擦力做功Wf2;

(3) 阻擋桿發生塑性變形所需要的變形能Wp。故碰撞過程中,某一個時刻滑車速度由速度V1減至速度V2時,根據能量守恒定律,可得能量平衡式

(3)

1.3 摩擦力做功及阻擋桿變形能

1.3.1 滑車與軌道之間摩擦力做功 在運動過程中,滑車與軌道通過4個圓柱形鋼制車輪與軌道接觸,假設滑車車輪與軌道之間的滾動摩擦阻力因子為μ1,那么運動過程中,滑車單個車輪所受摩擦力為

(4)

如圖5所示,在滑車由速度V1減速至V2,滑車移動距離為L的過程中,滑車與軌道之間摩擦力做功

Wf1=4FfL=μ1mgL

(5)

1.3.2 變形協調關系 為計算碰撞過程中各因素所消耗滑車動能的量,需要明確滑車位移、柱銷處阻擋桿彎曲變形量、彎曲角度之間的變形協調關系。為簡化計算模型,假設在運動過程中,柱銷之間的阻擋桿不發生偏轉,其變形示意圖如圖5所示。

圖5 阻擋桿變形示意圖Fig.5 Diagram of block bar deformation

由圖5中的幾何關系,可得

(6)

(7)

(8)

故單側柱銷處阻擋桿彎曲變形量

(9)

(10)

(11)

1.3.3 阻擋桿與柱銷之間摩擦力做功計算 根據圖3滑車受力分析圖,由牛頓第二定律可得

Ma=F′+4Ff

(12)

其中,a為滑車的加速度。由式(4)、 (12)可得

F′=Ma-4Ff=Ma-μ1Mg

(13)

假設阻擋桿與柱銷之間的滑動摩擦因子為μ2,故當阻擋桿移動一個微量距離dl時,柱銷與阻擋桿之間的摩擦力做功為

dwf2=μ2(FN1+FN2+FN3+FN4)dl=

(14)

在滑車由速度V1減至V2、阻擋桿在柱銷中移動距離為l的過程中,阻擋桿與柱銷之間的摩擦力做功為

(15)

1.3.4 阻擋桿彎曲變形能的計算 為計算碰撞過程中阻擋桿變形所需要的能量,首先要明確塑性變形時,材料內部應變能密度的大小。

寬度為b的阻擋桿發生如圖6所示的彎曲變形時,由力學知識可知,與其中性層距離為y(-b/2≤y≤b/2)處的應變為

(16)

式中,ρ為曲率半徑。

圖6 彎曲變形Fig.6 Bending deformation

由文獻[11]可知,對于線性強化剛塑性模型,材料在塑性變形過程中,其應力與應變之間的本構關系為

σ=σs+Hε

(17)

式中:σs為材料的屈服強度;H為材料的切線模量。

由文獻[12]可知,阻擋桿塑性變形時的應變能密度為

(18)

碰撞過程中,阻擋桿會在柱銷處以及碰撞桿頭部倒角處發生塑性變形。將整個碰撞過程滑車位移離散為若干個長度為dL的區間,由圖5可知,該位移區間內,有長度為dl的阻擋桿通過O1A截面在柱銷處發生彎曲變形,所需變形能為dWp1;同時有長度為dl-ρ1dθ的阻擋桿脫離柱銷,隨臺車向前運動的過程中又恢復至原狀,該段長度的阻擋桿發生變形所需變形能量為dWp2;在碰撞桿頭部的單側倒角處有長度為ρ2dθ的阻擋桿發生彎曲變形,所需變形能為dWp3。假設阻擋桿寬度為b,高度為h,那么對于長度為dl的阻擋桿彎曲變形所需要的能量為

(19)

將式(16)、(18)帶入式(19)并積分可得

(20)

同理可得

(21)

(22)

滑車由速度V1運動至速度V2的過程中,滑車移動距離為L,單側柱銷處阻擋桿彎曲變形量為l、彎曲角度為θ。則該過程中,單側柱銷處阻擋桿彎曲變形所需要的變形能為

(23)

同理可得阻擋桿脫離柱銷后恢復變形所需要的變形能為

(24)

碰撞桿頭部單側倒角處阻擋桿變形所需要的變形能為

(25)

減速裝置完全對稱,所以碰撞過程中阻擋桿彎曲變形所消耗的總能量為

Wp=2(Wp1+Wp2+Wp3)

(26)

1.4 力學模型求解方法

為分析碰撞過程中滑車速度、減速度等參數的變化情況,本文采用如下方法對模型進行求解:

(1) 給定滑車質量M、柱銷半徑R1、碰撞桿頭部倒角半徑R2、阻擋桿截面尺寸b、h,柱銷圓心與碰撞桿頭部倒角圓心的距離l0以及柱銷之間的中心距l1、l2;

(2) 給定初值l=0、L=0、θ=0,設定滑車微小位移dL的大小,以及滑車初始碰撞速度V1的值;

(3) 計算滑車在當前狀態下滑車移動dL時摩擦力做功Wf1、Wf2及阻擋桿彎曲變形所需要的變形能Wp,根據式(3)計算滑車移動dL后的速度V2;

(4) 根據式(27)、(28)求解該位移區間內滑車的減速度an、時間tn,并計算本次運動結束時L、l、θ的大小,作為下次運動的初始狀態;

(27)

(28)

(5) 判斷V2≥10-8是否成立。若成立,將V2作為下一位移區間的初始碰撞速度V1,重復步驟(3)～(5),若不成立,繪制相應的變化曲線,結束

計算。

按照上述的計算方法,利用Matlab編程計算,其過程如圖7所示。

圖7 計算流程圖Fig.7 Calculation flow chart

2 算例研究

根據以上計算方法,以某一具體的碰撞試驗臺為算例進行計算研究,并分析各參數對滑車碰撞過程的影響?；嚺鲎苍囼灥闹饕獏狄姳?,阻擋桿材料為Q235,切線模量為6 100 MPa,截面為寬度b=20 mm、高度h=60 mm的矩形截面。

表1 滑車碰撞試驗參數

圖8示出了初始碰撞速度為60 km/h,滑車質量為m時對碰撞過程的影響。從圖中曲線可以看出:質量越大,碰撞位移越大,不同質量的滑車在碰撞過程中,其位移曲線變化大致相同,呈非線性變化;滑車質量對沖撞力影響很小,碰撞過程中,沖撞力先增大,然后保持基本不變,不同質量的滑車,其碰撞力變化曲線幾乎重合,從碰撞力變化曲線局部放大圖可以看出,隨著滑車質量的增加,碰撞過程中最大碰撞力稍有增加;滑車質量越大,滑車速度變化越慢,在碰撞過程中達到的最大減速度越小。

圖9示出了滑車質量M=900 kg時,滑車初始碰撞速度對碰撞過程的影響。從圖中曲線可以看出:初始碰撞速度越大,碰撞位移越大,位移曲線變化越快,達到最大沖撞力所需要的時間越短,但其對最大沖撞力影響較小;不同的初始碰撞速度,滑車的速度變化曲線大致相同;初始碰撞速度越大,滑車的減速度變化越快;初始碰撞速度對滑車所能達到的最大減速度影響不大。

圖10示出了初始碰撞速度為60 km/h時,內側柱銷中心距l2對碰撞過程的影響。從圖中的曲線可以看出:中心距越大,碰撞位移越大,沖撞力變化越慢,但其對最大沖撞力影響較小;中心距越小,滑車速度變化越快,滑車的減速度越大,但是中心距的變化對滑車所能達到的最大減速度影響不大。

圖8 滑車質量對碰撞過程運動參數的影響Fig.8 Influences of the sled quality on the motion parameters during collision

圖9 初始碰撞速度對碰撞過程運動參數的影響Fig.9 Influences of the initial collision speed on the motion parameters during collision

圖11示出了初始碰撞速度為50 km/h,阻擋桿矩形截面寬度b不變時,截面高度h對碰撞過程的影響。由圖中的曲線可以看出:阻擋桿截面高度越大,碰撞位移越小,沖撞力變化越快,且達到的最大沖撞力越大,碰撞過程滑車速度變化越快,得到的最大減速度越大。

圖11 截面高度對碰撞過程運動參數的影響Fig.11 Influences of the cross section height on the motion parameters during collision

3 結 論

(1) 以汽車碰撞試驗為研究對象,對碰撞過程中滑車的運動以及阻擋桿的變形進行了分析,根據能量守恒定律得到了碰撞過程的力學模型,并對模型中摩擦力做功及阻擋桿變形能建立了計算模型;構建了一種能夠對碰撞過程的位移、速度、減速度以及沖撞力進行計算的算法,以Matlab為計算工具進行編程,實現了對碰撞過程的運動學參數和沖撞力進行定量計算。

(2) 利用所建立的計算模型,以某一具體的汽車碰撞試驗臺為對象進行了算例研究,研究了滑車質量、初始碰撞速度、內側柱銷中心距、阻擋桿截面尺寸等參數對滑車碰撞過程的影響,得出了一系列曲線。研究結果表明:在其他參數保持不變的情況下,滑車的碰撞位移曲線變化大致相同,呈非線性變化,且滑車位移隨滑車質量的增加而增大,隨初始碰撞速度的增加而增大,隨內側柱銷中心距的增加而增大,隨阻擋桿截面高度的增加而減小;阻擋桿截面高度對滑車的最大沖撞力影響較大,其余參數對最大沖撞力影響較小;碰撞過程中,滑車質量越大,速度變化越慢,阻擋桿截面高度越大,速度變化越快,初始碰撞速度以及內側柱銷中心距對滑車速度變化影響較小;滑車的最大減速度隨滑車質量的增大而減小,隨初始碰撞速度的增加而增大,隨柱銷中心距的增大而減小,隨阻擋桿截面高度的增加而增大。

[1] 胡玉梅,曾繁林,張科峰,等.臺車試驗中鋼板布置的研究[J].汽車工程,2012,34(6):496-501.

[2] 曹立波,鐘志華,白中浩,等.臺車碰撞試驗用機械緩沖吸能裝置研究[J].汽車工程,2002,24(3):228-230.

[3] ZHANG Xu,JIANG Qi,LI Liwei,etal.Study on the crash energy absorption structure of a vehicle simulator crash testing equipment[J].Key Engineering Materials,2011,486:266-269.

[4] MASSENZIO M,MAUPAS A,BENNANI A,etal.Study of a device for controlling the pulses of sled testing[J].International Journal of Crashworthiness,2007,12(3):311-318.

[5] GIAVOTTO V,ANGHILERI M,CASTELLETTI L M L,etal.A braking system for tests with a prescribed deceleration pulse[J].Experimental Mechanics,2010,50(7):915-930.

[6] 葛如海,穆青,陳曉東,等.汽車模擬碰撞吸能器的仿真分析及試驗[J].江蘇大學學報(自然科學版),2006,27(4):324-327.

[7] 張慧云,柳立志,楊斌,等.減速臺車波形發生器的曲線復現及應用[J].汽車科技,2015(1):57-62.

[8] XU Fengxiang,WANG Chao.Dynamic axial crashing of tailor-welded blanks (TWBs) thin-walled structures with top-hat shaped section[J].Advances in Engineering Software,2016,96:70-82.

[9] MEHMET A G,MUHAMMED E C,BERTAN B,etal.The effect of geometrical parameters on the energy absorption characteristics of thin-walled structures under axial impact loading[J].International Journal of Crashworthiness,2010,15(4):377-390.

[10] JIMENEZ M,MARTINEZ J,FIGUEROA U,etal.Finite element simulation of mechanical bump shock absorber for sled tests[J].International Journal of Automotive Technology,2015,16(1):167-172.

[11] 尚福林,王子昆.塑性力學基礎[M].西安:西安交通大學出版社,2011:9-13.

[12] 沈超.關于塑性應變能密度斷裂準則研究[D].西安:西安理工大學,2013.

MotionandMechanicalPropertiesoftheCollisionProcessinVehicleCollisionTest

WANGJia-shun,ANQi

(SchoolofMechanicalandPowerEngineering,EastChinaUniversityofScienceandTechnology,Shanghai200237,China)

The paper took the vehicle collision test as the research object,analyzed the calculation method of the friction work and deformation energy of the block bar.Based on the principle of energy conservation,the energy balance model of the collision process is described,and the mathematical model is built to calculate the movement of the sled and the deformation of the block bar during the collision process,and a calculation method is proposed to solve this model by Matlab.With a specific example,the influences of the quality of sled,the initial collision speed,the center distance of the inner pin,and the cross section size of the block and so on,the displacement,collision force,speed and deceleration are studied,and some curves are obtained and analyzed.

collision test; deceleration device; mechanical modeling; motion properties; mechanical properties

1006-3080(2017)06-0863-08

10.14135/j.cnki.1006-3080.2017.06.017

2017-01-18

上海市設計學Ⅳ類高峰學科項目(DA17014)

王家順(1992-),男,山東人,碩士生,研究方向為機械設計及理論。

安 琦,E-mail:anqi@ecust.edu.cn

TH113.2+2

A