袁鵬飛　江丙云　劉俊磊

摘要：? 為研究某交流充電樁IK10沖擊仿真結果與試驗結果的差異性，用ANSA創(chuàng)建交流充電樁IK10沖擊仿真模型。將該模型文件導入Abaqus中進行仿真計算，得到沖擊頭的加速度曲線、沖擊頭的侵入位移曲線、樁體與地面之間的反作用力曲線，再對交流充電樁進行IK10沖擊試驗，將得到的數(shù)據(jù)與仿真結果進行對比，發(fā)現(xiàn)二者的數(shù)據(jù)誤差約為10%，在可接受的誤差范圍內，可以為同類產品的研發(fā)提供參考。

關鍵詞：? 充電樁; IK10沖擊; 加速度; 侵入位移; 反作用力; 對標

中圖分類號：? U491.8;TM910.6文獻標志碼：? B

Research of charging pile IK10 impact simulation

and test comparison

YUAN Pengfei JIANG Bingyun LIU Junlei

（1. Guochuang Mobile Energy Innovation Center （Jiangsu） Co.， Ltd.， Changzhou 213000， Jiangsu， China;

2. Wanbang Digital Energy Co.， Ltd， Changzhou 213000， Jiangsu， China）

Abstract： In order to study the difference between the simulation results? and the test results of IK10 impact of an AC pile， ANSA is used to create the simulation model of IK10 impact? of a charging pile， and the model file is imported into Abaqus for simulation calculation. The acceleration curve， the intrusion displacement curve of impactor， and the reaction force curve between pile and ground are obtained. The IK10 impact test on the charging pile is carried out， and the simulation results and the test results is compared. There is a deviation of about 10%， which belongs to the acceptable range. It can provide reference for research and development of similar product.

Key words： charge pile; IK10 impact; acceleration; intrusion displacement; reaction force; comparison

基金項目：? 國家重點研發(fā)計劃（2018YFB0106300）

作者簡介： 袁鵬飛（1988—），男，河南舞陽人，工程師，碩士，研究方向為充電設備和汽車仿真，（E-mail）ypf.1059566@163.com0引言

汽車產業(yè)在人們的生產生活中扮演著重要角色。作為一種清潔能源，電能正在逐漸取代石油作為汽車的動力來源，這也是節(jié)能減排與改善地球環(huán)境的必然趨勢[1]。目前，很多國家都在積極致力于電動汽車的研發(fā)與生產，我國的電動汽車行業(yè)發(fā)展速度較快。

國家大力支持新能源汽車的發(fā)展，積極支持其配套設施建設[2]，預計會投入超過600億元人民幣進行充電樁建設，其中充電設備的投入占比至少為1/4[3]。如此大的市場份額，促使各個企業(yè)在充電樁產品研發(fā)過程中投入大量精力。只有研發(fā)出安全且性能優(yōu)異的充電樁，企業(yè)才能在日趨激烈的競爭中立于不敗之地。

1IK沖擊

IK即電器設備外殼對外界機械碰撞的防護等級，是充電樁的一項重要測試項目，在充電樁的安全評估方面至關重要。在充電樁的研發(fā)過程中必須進行IK沖擊仿真。

本文對充電樁進行IK10（即以5 kg的鋼球或半球體沖擊頭，在0.4 m的高度豎直下落到目標物表面，沖擊能量為20 J）沖擊仿真，并在充電樁制造完成后進行IK10沖擊試驗，將結果進行對標，校正仿真參數(shù)，以應用到其他同類產品的研發(fā)中。IK等級與沖擊能量對應見表1。

2顯式動力學仿真

仿真過程由前處理、計算、后處理等3個步驟組成。前處理包括幾何數(shù)據(jù)處理、網格劃分、創(chuàng)建連接和邊界條件、施加載荷，計算是對仿真模型運行求解得到結果，后處理是將計算結果以云圖、曲線、等值線等形式展示出來。

作為一款世界通用的工程仿真軟件，Abaqus常用的仿真計算方法有隱式算法和顯式動力學算法。其中，顯式動力學算法適于求解復雜的非線性動力學問題和準靜態(tài)問題，特別適用于模擬瞬時的動態(tài)過程，如沖擊和爆炸問題[4]。

顯式動力學應用中心差分法對方程進行顯式時間積分，應用一個增量步的動力學條件，計算下一個增量步的動力學條件。在顯式算法中，質量矩陣總是采用對角矩陣或者集中質量矩陣，所以求解加速度簡單，不需要求解聯(lián)立方程[5-6]。

顯式動力學算法最大的優(yōu)點是有較好的穩(wěn)定性，不用直接求解切線剛度，不需要進行平衡迭代，計算速度快，時間步長只要取得足夠小，一般不存在收斂性問題，因此需要的內存也比隱式算法少。

3有限元模型的建立

3.1網格劃分

充電樁由裝飾蓋、中蓋、底殼、燈條、PCB組件等組成，由于各部件結構比較復雜，使用ANSA軟件將這些部件劃分為四面體網格。為減少計算時間，保證一定的計算精度，本文的網格平均尺寸為5 mm，得到的仿真模型的節(jié)點總數(shù)約為77 500個，單元總數(shù)約為239 700個。

3.2材料與屬性

本文研究的是IK沖擊仿真，屬于瞬態(tài)過程，需要輸入材料在不同應變率下的拉伸應力-應變數(shù)據(jù)（見圖1）。對于缺失的材料參數(shù)，可由材料供應商提供或委外測試。

在ANSA中創(chuàng)建3個材料：PC+ASA、FR4、PC，材料信息包括密度、泊松比、彈性模量和拉伸應力-應變數(shù)據(jù)。裝飾蓋、中蓋、底殼的材質均為PC+ASA，PCB組件材質為FR4，燈條材質為PC。由于沖擊頭在整個沖擊過程中不發(fā)生變形，因此沖擊頭采用剛體建模，對其余各部件賦予材料和屬性，均采用SOLID SECTION屬性，材料信息見表2。

3.3模型建立

交流充電樁的塑膠部件之間，以及PCB板組件與底殼之間都通過螺栓連接，本文采用剛性Coupling單元模擬螺栓連接，采用剛性板模擬剛性地面。在沖擊試驗中不需要考慮沖擊頭的變形，因此用剛體球模擬半球形沖擊頭[7]，在剛體球的中心位置設置加速度傳感器單元，剛體球的質量為5 kg，剛體球與樁體之間創(chuàng)建面面接觸，其他部件之間創(chuàng)建全局自動接觸，全局接觸的滑動摩擦因數(shù)設為0.2，圖2為沖擊仿真模型和局部細節(jié)示意。

3.4邊界條件

將充電樁平放于固定的剛性板上，為減少仿真計算時間，在ANSA中將剛體球放置在樁體上表面2 mm的位置。根據(jù)自由落體運動知識和能量守恒定律，通過式（1）換算得到剛體球自由下落到樁體表面時的速度，即

mgh=mv²/2 （1）

式中：m為剛體球的質量，kg;g為重力加速度，m/s²;h為跌落高度，m;v為剛體球撞擊樁體表面時的速度，m/s。

給剛體球施加豎直向下2.8 m/s的初速度，模擬其從0.4 m的高度自由下落到充電樁表面的過程。

3.5計算及結果后處理

在ANSA中設置好仿真模型，導出以.inp結尾的計算文件，然后打開Abaqus求解器求解界面（DOS窗口），用Abaqus JOB=文件名.inp命令進入計算環(huán)節(jié)[8]。Abaqus求解器計算完成后，可在META中直接讀取.odb結果文件。仿真結果顯示，樁體的最大應變約為10.796%，位于裝飾蓋和中蓋與鋼體球接觸的區(qū)域，小于其材料PC+ASA的斷裂伸長率42.6%。因此，裝飾蓋和中蓋等塑膠部件的破損、開裂風險較小。通過META提取剛體球的加速度隨時間的變化曲線、剛體球的侵入位移隨時間的變化曲線、樁體與地面間的反作用力隨時間變化的曲線等信息，圖3為IK10沖擊仿真充電樁外殼的應變云圖。

4IK10沖擊試驗

對充電樁進行IK10沖擊試驗，沖擊試驗裝置示意見圖4。將充電樁自由放置在水平放置的鋼板基座上，在基座底部設置一個測力傳感器（見圖5），用于捕捉在沖擊作用下充電樁與剛板基座間的作用力隨時間的變化歷程。

將質量為5 kg的半球狀沖擊頭安裝在落錘臺架上，沖擊頭位于充電樁上部0.4 m處（考慮到豎直滑軌存在一定的阻尼，所以沖擊頭的高度會稍大于0.4 m）。沖擊頭沿著滑軌下落，保證沖擊方向豎直向下。在沖擊頭的中心位置設置2個加速度傳感器（見圖6），用于捕捉沖擊頭在沖擊過程的加速度隨時間的變化歷程。測力傳感器和加速度傳感器連接到數(shù)據(jù)采集設備（見圖7），采用高速攝像機記錄整個沖擊過程。試驗前在充電樁的裝飾蓋上噴一層白色自噴漆，以便標記沖擊的位置。

5IK10沖擊試驗結果

試驗得到沖擊頭的2個加速度傳感器的平均值隨時間變化的曲線見圖8，沖擊頭侵入位移隨時間變化的曲線見圖9，底部力傳感器捕捉的樁體與剛板基座之間的作用力隨時間變化的曲線見圖10。

6仿真與試驗結果對比

充電樁IK10沖擊試驗結果顯示樁體沒有出現(xiàn)破損、裂縫、脫落等現(xiàn)象，與仿真結果一致。沖擊試驗得到的原始加速度數(shù)據(jù)存在許多尖峰，對其進行濾波處理，可得到比較符合實際的加速度曲線。濾波后的沖擊頭加速度曲線見圖11。

沖擊頭加速度的試驗與仿真結果對比見圖12。由于仿真計算中沖擊頭只有1次沖擊充電樁的過程，沒有發(fā)生回落并與充電樁的二次接觸，因此仿真的加速度曲線只有1個峰值。試驗中沖擊頭反彈后回落，再次與充電樁接觸，因此出現(xiàn)2個較大的加速度峰值。本次只對比第一個加速度峰值的結果。試驗從沖擊頭開始下落的時刻開始計時，而仿真從沖擊頭即將接觸充電樁時開始計時，因此圖13中仿真的加速度曲線是將時間軸平移得到的。

由圖12可知，沖擊測試得到的沖擊頭的平均加速度最大值約為169.09 m/s²，沖擊仿真得到的沖擊頭最大加速度約為161.28 m/s²，且加速度曲線的趨勢也很接近，二者的誤差保持在合理誤差范圍（10%）內，說明仿真結果與測試結果有很好的一致性。

沖擊頭的最大侵入位移仿真和試驗示意見圖13，二者結果對比見圖14，其中仿真結果曲線通過時間軸平移得到。由圖14可知，沖擊測試中沖擊頭的最大侵入位移約為27.8 mm，沖擊仿真中沖擊頭的最大侵入位移約為25.9 mm，仿真誤差保持在合理誤差范圍（10%）內，說明仿真結果與測試結果有很好的一致性。需要注意，仿真的沖擊過程時間比試驗的沖擊時間偏短。

圖15為樁體與剛性地面間的反作用力的試驗與仿真結果對比。與加速度的曲線類似，試驗中由于沖擊頭的反彈回落，出現(xiàn)多個反作用力的峰值，本文只對比沖擊頭與樁體第一次接觸時的反作用力，仿真結果曲線通過時間軸平移得到。

由圖15可知，沖擊測試的樁體與剛性地面間的反作用力約為2 125.9 N，仿真得到的反作用力約為2 001.0 N，二者的誤差保持在合理誤差范圍（10%）內，說明仿真結果與測試結果有很好的一致性。

7結束語

對某款交流充電樁進行IK10沖擊仿真和試驗，發(fā)現(xiàn)樁體破壞和開裂的風險均較小。為更精細化地對比仿真與試驗結果，分別提取沖擊頭的加速度變化曲線、侵入位移變化曲線、充電樁與剛性地面間的作用力變化曲線。對比3種曲線發(fā)現(xiàn)，曲線趨勢基本一致，數(shù)值偏差保持在可接受的誤差范圍（10%）內，證明該仿真方法在評估充電樁IK10沖擊性能方面有很好的參考性。參考文獻：

[1]王旭， 齊向東. 電動汽車智能充電樁的設計與研究[J]. 機電工程， 2014， 31（3）： 393-396. DOI： 10.3969/j.issn.1001-4551.2014.03.027.

[2]吳春陽， 黎燦兵， 杜力， 等. 電動汽車充電設施規(guī)劃方法[J]. 電力系統(tǒng)自動化， 2010， 34（24）： 36-39.

[3]劉康華， 張會， 韓娟娟. 交直流一體充電樁的結構設計[J]. 科學技術創(chuàng)新， 2018（35）： 158-160.

[4]石亦平， 周玉蓉. Abaqus有限元分析實例詳解[M]. 北京： 機械工業(yè)出版社， 2008： 7-8.

[5]SHERWOOD J A， FROST C C. Constitutive modeling and simulation of energy absorbing polyurethane foam under impact loading[J]. Polymer Engineering and Science， 1992， 32（16）： 1138-1146.

[6]曾翔. 沖擊和快速加載作用下鋼筋混凝土梁柱構件性能試驗與數(shù)值模擬研究[D]. 長沙： 湖南大學， 2014.

[7]朱東俊， 葛亮， 劉瑩， 等. 基于Abaqus的復合材料板沖擊特性分析[J]. 計算機輔助工程， 2013， 22（S2）： 21-25.

[8]劉明卓. 基于Abaqus的汽車座椅塑料件有限元分析與結構優(yōu)化[J]. 塑料工程學報， 2011， 18（4）： 116-119. DOI：10.3969/j.issn.1007-2012.2011.04.023.（編輯陳鋒杰）