高發(fā)華,羅明軍
(奇瑞商用車(安徽)有限公司工程研究院,安徽蕪湖 241006)
減振器彈簧閥片的力學建模及CAE仿真分析
高發(fā)華,羅明軍
(奇瑞商用車(安徽)有限公司工程研究院,安徽蕪湖 241006)
針對用傳統(tǒng)材料力學分析難以精確獲得減振器工作狀態(tài)下伸張閥和壓縮閥彈簧閥片組的應力和應變值的難題,采用有限元法對其工作狀態(tài)下的壓力分布進行研究。根據(jù)減振器工作狀態(tài)下彈簧閥片組的內部受力關系,基于有限元軟件ANSYS,通過創(chuàng)建彈簧閥片之間的接觸對來模擬彈簧閥片之間的約束關系。為準確分析閥片組的受力關系,建立了伸張閥和壓縮閥的力學模型,并結合減振器工作特性曲線定義伸張閥和壓縮閥的載荷。通過對伸張閥和壓縮閥的強度和剛度進行分析,得到了工作狀態(tài)下減振器伸張閥和壓縮閥的應力、應變分布云圖,為后續(xù)優(yōu)化設計提供依據(jù)。
減振器;伸張閥;壓縮閥;力學模型;有限元分析
對于汽車減振器的研究,以前對減振器閥片的分析都是利用材料力學的公式推算,把減振器的閥片簡化,從而推算出減振器閥片的撓度公式[1]。這種方法因為簡化了閥片,對于形狀稍微復雜的減振器閥片就不再適用[2]。而隨著有限元分析的興起,利用計算機對減振器的閥片進行CAE分析可以得到更精確的減振器閥片的應力和應變情況[3]。在某汽車前減振器中,伸張閥由3片閥片和1個墊片組成的,而壓縮閥由3片閥片組成,所以只是利用材料力學的公式推導其撓度公式是難以實現(xiàn)的,而利用有限元軟件分析其應力和應變分布情況則可以實現(xiàn)[4-6]。作者擬運用商業(yè)通用有限元軟件ANSYS對某汽車前減振器的伸張閥和壓縮閥的應力和應變進行分析,考察閥片強度和剛度性能。
該汽車前減振器產品的活塞連桿組件中包括的閥門有流通閥和伸張閥,各個閥片的三維實體都是較規(guī)則的圓形片體,其中伸張閥由3個閥片和1個墊片組成[7-9]。該前減振器產品的底閥組件包括:底閥座、補償閥限位、補償閥彈簧、補償閥閥片、壓縮閥底座、壓縮閥片(3片)、銷。其中伸張閥和壓縮閥是減振器組件中關鍵的部件,是影響減振器性能的主要因素。
伸張閥和壓縮閥在使用過程中屬于大撓度變形問題,主要是徑向方向的彎曲變形。為了更加準確地研究伸張閥和壓縮閥在工作過程中的變形問題,基于UG軟件建立了汽車減振器彈簧閥片的實體模型,分別如圖1和圖2所示。
圖1 伸張閥
圖2 壓縮閥
2.1 模型的導入
首先將實體模型以*.stp格式保存,然后通過UG和有限元分析軟件ANSYS的接口將模型導入有限元分析環(huán)境中。模型在接口轉換過程中不會出現(xiàn)數(shù)據(jù)丟失導致模型失真現(xiàn)象,這是有限元分析高精度得到保障的前提。
2.2 定義屬性
在建立減振器的伸張閥和壓縮閥有限元模型之前,先對分析軟件ANSYS參數(shù)進行設置,包括分析類型、材料屬性等。分析類型為結構分析。材料均為線性材料[11],具體如表1所示。
表1 彈簧閥片材料參數(shù)表
2.3 網(wǎng)格劃分
對于減振器的彈簧閥片而言,它屬于薄板元件。在采用ANSYS軟件實現(xiàn)網(wǎng)格劃分時將閥片定義為殼單元,故采用SOLID45單元對彈簧閥片進行網(wǎng)格劃分。最終減振器彈簧閥片的伸張閥和壓縮閥有限元模型分別離散為12 500和12 722個
節(jié)點,51 531和52 838個單元。
2.4 邊界條件的施加
(1)接觸對的創(chuàng)建
由于壓縮閥和伸張閥都不是由單個閥片組成,是由多個閥片組成的,所以在分析過程中需考慮它們之間的約束關系。為能夠準確模擬閥片之間的關系,需要妥善考慮其內部受力情況,否則分析結果誤差過大。但接觸問題是一種高度非線性行為,需要較大的計算資源,為了進行實際有效的計算,理解問題的特性和建立合理的模型是很重要的。接觸問題存在兩個較大的難點:①在求解問題之前,不知道接觸區(qū)域,表面之間是接觸或分開是未知的、突變的,這由載荷、材料、邊界條件和其他因素而定;②大多的接觸問題需要計算摩擦,有多種摩擦和模型可以選擇,但摩擦使問題的收斂性變得困難。作者采用面-面接觸模擬彈簧閥片之間的連接關系,不但能較好地表征閥片之間的內部關系,而且減少了閥片組之間由于摩擦副帶來的大量計算工作量。
(2)添加邊界條件
圖3和圖4分別為該汽車前減振器伸張閥彈簧閥片和壓縮閥彈簧閥片的力學模型。伸張閥閥片簡化為受均布載荷q作用、內邊緣為固定夾緊的圓環(huán)薄板;壓縮閥閥片簡化為受均布載荷q作用、外邊緣為簡單鉸鏈支承的圓環(huán)薄板。
圖3 伸張閥彈簧閥片的力學模型
圖4 壓縮閥彈簧閥片的力學模型
根據(jù)減振器的實際使用工況下減振器的示功圖[12](如圖5所示),對閥片組的載荷進行定義。減振器在伸張至工作狀態(tài)下,伸張閥片受到油液通過6個均布的活塞孔而作用在伸張閥上的力,大小為50 N。減振器在壓縮工作狀態(tài)下,壓縮閥受到4個均布的力,大小為100 N。
圖5 示功圖
對于壓縮閥,減振器在壓縮運動的整個過程中,由于壓縮閥底座的作用,最內圈的邊緣都是完全固定的。而伸張閥則由于彈簧的作用,先會向下運動小段距離后,彈簧壓縮到一定程度就可以完全抵消伸張閥受到的油壓力。這里考慮的是彈簧壓縮之后的,所以伸張閥的最內圈也是完全固定的,得到伸張閥和壓縮閥有限元模型,分別如圖6和圖7所示。
圖6 伸張閥有限元模型
圖7 壓縮閥有限元模型
通過ANSYS有限元軟件,根據(jù)減振器在實際工況下工作原理,對減振器彈簧閥片的伸張閥和壓縮閥進行強度和剛度分析,得到了壓縮閥和伸張閥的應力和應變圖,分別如圖8—11所示。
圖8 壓縮閥應力分布
圖9 壓縮閥應變分布
圖10 伸張閥應力分布
圖11 伸張閥應變分布
由圖8—11可以明顯看出壓縮閥和伸張閥的最危險受力截面均在壓縮閥內圈的邊緣處,而最大應變均發(fā)生在壓縮閥的外圈,最大應力和應變值見表2,這些均與實際情況相符。
表2 彈簧閥片的最大應力和應變值
通過對汽車前減振器的強度和剛度進行分析,獲得了前減振器在工作狀態(tài)下伸張閥和壓縮閥彈簧閥片的應力應變分布云圖,得到以下結論:
(1)建立減振器彈簧閥片的實體模型,以*.stp格式通過接口技術導入ANSYS有限元軟件環(huán)境中,具有建模快、數(shù)據(jù)不丟失、模型不失真優(yōu)勢。
(2)通過對減振器彈簧閥片的受力分析,建立了彈簧閥片伸張閥和壓縮閥的物理力學模型,并結合減振器工作特性曲線對伸張閥和壓縮閥彈簧閥片有限元模型進行邊界定義。
(3)在減振器彈簧閥片有限元模型基礎上對它們進行了強度和剛度分析,計算結果與工作工況下的狀態(tài)較為一致,解決了傳統(tǒng)材料力學難以直接計算得到復雜彈簧閥片結構應力應變的難題,為汽車減振器的性能改善提供了可靠的依據(jù)。
【1】李世民,呂振化.汽車減震器阻力—速度特性的預測法[J].機械設計與制造,2004(6):76-78. LI S M,LV Z H.Prediction Technique of Automotive Shock Absorber
Force versus Velocity Characteristics[J].Machinery Design & Manufacture,2004(6):76-78.
【2】程小彪.基于UG的JL6370汽車減震器虛擬設計[D].武漢:武漢汽車工業(yè)大學,2000.
【3】張德慶.汽車減震器的建模與仿真設計[D].鎮(zhèn)江:江蘇大學,2007.
【4】SAMANTARAY A K.Modeling and Analysis of Preloaded Liquid Spring/Damper Shock Absorbers[J].Simulation Modeling Practice and Theory,2009,17(1):309-325.
【5】GALLARDO J M,SORIA L,HERRERA E J.Investigation of Service Failures in Automobile Shock Absorbers[J].Engineering Failure Analysis,2007,14(2):355-363.
【6】RIVAS A,BIERA J,SACRAMENTO G,et al.Development of a Thermal Model for Automotive Twin-tube Shock Absorbers[J].Applied Thermal Engineering,2005,25(11):1836-1853.
【7】劉延慶,張建武,張利,等.車用雙筒液壓減震器動態(tài)響應的試驗與仿真[J].上海交通大學學報,2002,36(8):1095-1099. LIU Y Q,ZHANG J W,ZHANG L,et al.Experiment and Simulation of Dynamic Response for Vehicle Twin Tube Hydraulic Shock Absorbers[J].Journal of Shanghai Jiaotong University,2002,36(8):1095-1099.
【8】李幼德,李靜,宋大鳳,等.汽車減震器彈簧閥片變形模型的研究[J].汽車工程,2003,25(3):287-290. LI Y D,LI J,SONG D F,et al.A Research on the Deformation of Spring Valve Plate in Automotive Shock Absorbers[J].Automotive Engineering,2003,25(3):287-290.
【9】袁長頌.有限元法在減振器建模中的應用[J].客車技術與研究,2007,29(2):20-22. YUAN C S.Application of Finite Element Method to Shock Absorber Modeling[J].Bus Technology and Research,2007,29(2):20-22.
【10】史廣奎,李檳,孟憲民.汽車設計中減震器相對阻尼系數(shù)的確定[J].汽車工程,1995(6):367-374.
【11】黃恒,程廣偉,鄧楚男.車用減震器的外特性建模與仿真[J].機車技術,2005(11):4-7. HUANG H,CHENG G W,DENG C N.Modeling and Simulation of External Characteristic of Automotive Shock Absorbers[J].Automobile Technology,2005(11):4-7.
【12】馮雪梅,劉佐民.汽車液力減震器技術的發(fā)展與現(xiàn)狀[J].武漢理工大學學報(交通科學與工程版),2003,27(3):340-343. FENG X M,LIU Z M.Development and Current Situation of Automotive Hydraulic Shock Absorber Technology[J].Journal of Wuhan University of Technology(Transportation Science & Engineering),2003,27(3):340-343.
【13】張旭.雙向減震器的結構及特性研究[D].大慶:大慶石油學院,2003.
Mechanics Modeling and CAE Simulation Analysis of Spring Valve Plate on Shock Absorber
GAO Fahua,LUO Mingjun
(Automotive Engineering Research Institute,Wuhu Chery Commercial Vehicle Company, Wuhu Anhui 241006,China)
According to difficult to accurately obtain stress and strain value of rebound valve and compression valve on shock absorber in working status by using traditional materials mechanics, finite element analysis was applied in the research. Constraint relation among spring valve plates on shock absorber was simulated by creating contact pair based on ANSYS according to internal force relations of spring valve plates group on shock absorber in working status. The mechanical models of rebound valve and compression valve were established in order to accurately analyze force relation of the valve plates group. The loads of rebound valve and compression valve were defined by combining performance curve of shock absorber. The stress and strain distribution of rebound valve and compression valve were obtained by analysis of strength and rigidity when shock absorber was in working status. It provides basis for subsequent optimization design.
Shock absorber; Rebound valve; Compression valve; Mechanics model;Finite element analysis
2017-01-13
國家自然科學基金資助項目(51405123);江西省教育部項目(贛教技字[2007]20號)
高發(fā)華(1979—),男,主要研究方向為汽車整車技術研究與工程開發(fā)。
羅明軍,E-mail:lmjlmh2008@163.com。
10.19466/j.cnki.1674-1986.2017.05.008
U461.1
A
1674-1986(2017)05-038-04