劉瑩，翟斌，于成龍，趙前程

(華晨汽車工程研究院CAE工程室，遼寧沈陽 110141)

?

四連桿無骨雨刮力學性能仿真

劉瑩，翟斌，于成龍，趙前程

(華晨汽車工程研究院CAE工程室，遼寧沈陽 110141)

針對正向開發(fā)中的雨刮系統(tǒng)，應用有限元方法，校核四連桿無骨雨刮工作一個循環(huán)的機構運動，得到風擋的接觸壓力分布及軌跡，為刮片的疲勞計算提供基礎。

有限元；四連桿無骨雨刮；力學性能；機構運動

0　引言

生活中人們發(fā)現(xiàn)下雨天有些時候，雨刷在玻璃上留下擦拭不均的痕跡。刮刷效果設計中要求初期到50萬次不允許出現(xiàn)任何刮刷痕跡，100萬次可以出現(xiàn)局部模糊刮痕(即接觸壓力為0)，但往往這個指標很難實現(xiàn)。排除雨刷硬化的原因，設計初始的壓力分配效果對膠條功能和壽命有著重要的影響[1]。作者對電動無骨雨刮進行仿真分析，旨在模擬設計初始狀態(tài)刮刷效果及壓力分布，為設計提供參考。

1　分析思路

由于雨刮結構本身的力學原理是基于預應力的基礎之上，所以如何恰當?shù)貙嶋H結構轉化成力學模型是分析成功的關鍵，也是分析精度的保證。作者首先模擬雨刮由自由狀態(tài)到安裝狀態(tài)下簧片的變形導致的應力，再將應力以場的形式引入裝配體。同時通過彈簧的力學分析得到彈簧的剛度并考慮初始伸長量，將曲線偏移實現(xiàn)雨刮預壓。

2　力學模型的細節(jié)處理

2.1力學模型建立

2.1.1網(wǎng)格的劃分、單元的選取及接觸的設定

眾所周知：刮片的本體由于其長寬比的原因，通常采用平面應變單元進行二維分析。但由于其斷面的力隨時間和位置變化，三維分析更能精準地捕捉歷程信息。對于單元的劃分，作者采用SWEEP的形式將二維單元延伸得到三維實體，如圖1所示。這樣除了減少計算花費外，對于模型的收斂性也提供一定的保障。

圖1　刮片網(wǎng)格示意圖

膠條本體與簧片作者用節(jié)點耦合處理，同時將膠條的單元選取為雜交單元。接觸行為一共有3個，均以*CONTACT PAIR來實現(xiàn)，如圖2所示。其中膠條兩側的接觸在模型設置中設定，不同分析步中激活和關閉。

圖2　接觸位置

2.1.2膠條特性處理

由于橡膠本身的特性，在模擬仿真中作者采用超彈+黏彈本構，利用關鍵字來實現(xiàn)[2]。

*HYPERELASTIC, NEO HOOKE, TEST DATA INPUT, MODULI=LONG TERM

*VISCOELASTIC, TIME=RELAXATION TEST DATA, NMAX

而由于引用橡膠的本構關系，ABAQUS軟件會自動選取單元的控制語句，引入適當?shù)纳陈┛刂啤?/p>

2.2初始變形及初始應力

2.2.1獲取彈簧特性，模擬預張緊力

通常彈簧的剛度由供應商提供給OEM，也可以通過有限元對已有的CAD模型進行剛度的模擬，以此來匹配合適的剛度。如圖3所示，作者對彈簧施加軸向力得到軸向變形，進而得到剛度。彈簧剛度近似為線性行為，但為精確模擬預拉，作者偏移曲線，引入非線性材料實現(xiàn)這一效果，如圖4所示。

圖3　彈簧變形分析

圖4　彈簧力學性能

通過關鍵字*CONNECTOR ELASTICITY, COMPONENT=1, NONLINEAR 賦予彈簧單元CONNECTOR AXIAL非線性屬性，實現(xiàn)彈簧預壓的力學模型。值得一提的是：兩個刮臂由于其幾何特性，彈簧的性能和預伸長不同。

2.2.2簧片預壓，引入預應力場

抵抗簧片曲度回復的力，是影響刮擦效果的重要因素。通常用兩種方式實現(xiàn)簧片的預壓：一種是直接設定分析步，結果在后續(xù)的分析中繼承，如圖5所示。這樣做的好處是分析具備連貫性，缺點是增加分析難度。另一種方法是將本體的結果通過IMPORT的方式將應力場引入，這樣做的好處是可以迅速找到合適的載荷來對應后續(xù)分析的初始形態(tài)，如圖6所示。

圖5　刮片本體變形圖

圖6　初始應力場

2.3四連桿及刮臂連接的處理

電機的輸出軸與四連桿之間的搖臂，作者采用ABAQUS強大的CONNECTOE單元HINGE+BEAM來實現(xiàn)，而搖臂與拉桿之間用JOIN+ROTATION單元模擬，如圖7所示。由于雨刮系統(tǒng)是三維四連桿結構，這允許拉桿兩端有相對的Z向運動，而由于電機輸出端已經(jīng)采用JOIN+ROTATION單元，所以另一端為防止軸向旋轉采用UJOINT 單元來模擬(見圖8)，這樣可以達到模擬實際的效果。

圖7　曲柄球鉸連接副單元

圖8　拉桿萬向節(jié)連接副單元

刮臂的結構除彈簧以外，均為鉸鏈HINGE單元，包括刮片、簧片、雨刮臂之間，如圖9所示。其中彈簧單元的屬性，引用第2.2節(jié)中得到的彈簧特性。

圖9　刮臂連接副單元

2.4極限位置膠條的反轉

當雨刮系統(tǒng)工作半個循環(huán)，膠條單元面臨突然反轉的可能，如圖 10所示。這個時候可以采用隱式動力學處理嚴重接觸不連續(xù)的現(xiàn)象，使分析收斂性大大提高。同時前半個循環(huán)，膠條一側有接觸，而另一側只有后半循環(huán)才會接觸，這樣在前半分析步，可以利用關鍵字激活和關閉必要的接觸。

*MODEL CHANGE, ADD, TYPE=CONTACT PAIR

S3601;NSET4, S3602;FACET4

*MODEL CHANGE, REMOVE, TYPE=CONTACT PAIR

S3603;NSET5, S3604;FACET5

圖10　膠條反轉示意圖

3　后處理及腳本的應用

考察雨刮設計是否合理的一個重要評判是接觸壓力是否足夠?？梢栽谲浖暮筇幚砜吹侥硶r刻的接觸壓力，但如何將這些歷程累加成為云圖是作者需要考慮的問題。基于此作者對ABAQUS進行二次開發(fā)，利用腳本Sweep Pressure.py[3]實現(xiàn)壓

力軌跡，如圖11為不同限值的壓力云圖顯示，滿足設計要求。

圖11　接觸壓力累積變形圖

原理是在后處理中提取每個點的壓力，累積后通過場的形式施加給玻璃。如圖11所示，無骨雨刮的刮擦效果合格，但其沿刮片方向壓力梯度較大，這容易造成雨刮中部壓力較大易疲勞。建議工程調整簧片預彎，使壓力更為均勻。

通過圖12可知，利用ABAQUS OVERLAP功能顯示雨刮系統(tǒng)工作一個循環(huán)下的玻璃壓力分布穩(wěn)定無突變，從而判定系統(tǒng)分析的穩(wěn)定性。

圖12　一個工作循環(huán)玻璃應力分布

4　結論及展望

(1)實現(xiàn)整個雨刮系統(tǒng)的仿真模擬，成功得到接觸壓力軌跡和玻璃壓力分布，為前期設計提供指導。

(2)局限于常溫分析，并沒有考慮不同溫度下的力學性能。溫度場的模擬僅僅對材料提出更全面的要求即可，不存在技術問題。

(3)下一步工作應以此為基礎結合優(yōu)化軟件進行優(yōu)化設計，同時開展疲勞仿真工作。

【1】邵剛.汽車雨刷片設計[J].安徽電子信息職業(yè)技術學院學報,2012,11(5):43-47.

SHAO G.The Design of Vehicle Wiper Blade[J].Journal of Anhui Vocational College of Electronics & Information Technology,2012,11(5):43-47.

【2】ABAQUS幫助文檔[M].

【3】曹金鳳.Python語言在ABAQUS中的應用[M].北京：機械工業(yè)出版社,2011.

Simulation on Mechanical Property of Four-bar Linkage Boneless Wipers

LIU Ying, ZHAI Bin, YU Chenglong, ZHAO Qiancheng

(CAE Section, Brilliance Automobile Engineering Research Institute,Shenyang Liaoning 110141,China)

Aiming at the wipers system in forward development process, finite element method was used to check a circular motion of the four-bar linkage boneless wipers. The pressure distribution and trajectory of the windshield were gotten. It provides basis for the calculation of the blade fatigue.

FEA; Four-bar linkage boneless wipers; Mechanical properties; Mechanism motion

2016-05-18

劉瑩(1982—)，男，碩士，工程師，研究方向為結構耐久。E-mail:ying.l@brilliance-auto.com。

U463.85+5

A

1674-1986(2016)08-015-04