基于流固雙向耦合的轎車氣動與流致振動特性

胡興軍 王澤偉 于果 劉飛 蘭巍 余天明 桑濤

摘要：傳統(tǒng)流體力學(xué)數(shù)值仿真不考慮車身彈性結(jié)構(gòu)與外流場之間的相互作用，致使得到的結(jié)果與真實情況不符.以某實車模型為研究對象，對其進(jìn)行考慮流固耦合效應(yīng)的CFD仿真，并在氣動阻力、氣動升力、流場結(jié)構(gòu)等方面與傳統(tǒng)流體數(shù)值仿真結(jié)果對比，結(jié)果表明：流固耦合效應(yīng)對氣動升力影響較大，隨車速增加兩種仿真方法差異率可達(dá)到38%，直接關(guān)系到車輛穩(wěn)定性及安全性.利用流固耦合CFD數(shù)值仿真探究整車風(fēng)激振特性，證明了實車振動幅值主要影響因素為風(fēng)激振頻率及作用力大小.進(jìn)一步通過對車輛彈性結(jié)構(gòu)的剛度優(yōu)化改善汽車風(fēng)激振現(xiàn)象，從而提高乘員的舒適度.

關(guān)鍵詞：汽車空氣動力學(xué)；雙向流固耦合；數(shù)值仿真；流場；安全性

中圖分類號：U461.1文獻(xiàn)標(biāo)志碼：A

Automobile Aerodynamic and Flow Induced Vibration Characteristics Based on Bidirectional Fluid-structure Interaction

HU Xingjun，WANG Zewei，YU Guo，LIU Fei，LAN Wei，YU Tianming，SANG Tao

（State Key Laboratory7 of Automotive Simulation and Control，Jilin University，Changchun 130022，China）

Abstract：The traditional numerical simulation of fluid mechanics does not consider the interaction between the elastic structure of the car body and the external flow field，which makes the results inconsistent with the real situation. Taking a real vehicle model as the research object，the CFD simulation considering the effect of fluid-structure interaction is carried out，and the results of aerodynamic drag，aerodynamic lift，flow field structure and other aspects are compared with the traditional fluid numerical simulation results. The results show that the fluid-solid coupling effect has a greater impact on the aerodynamic lift，and the difference between the two simulation methods can reach 38% with the increase of vehicle speed，which is directly related to vehicle stability and safety. Using fluidstructure coupling CFD numerical simulation to explore the characteristics of wind-induced vibration of the whole vehicle，it is proved that the main influencing factors of the actual vehicle vibration amplitude are the wind-induced vibration frequency and the magnitude of the force. Moreover，the stiffness of the vehicle elastic structure is optimized to improve the wind-induced vibration of the vehicle，so as to improve the comfort of passengers.

Key words：automobile aerodynamics；bidirectional fluid-solid interaction；numerical simulation；flow field；safety

汽車空氣動力學(xué)由于影響汽車的燃油經(jīng)濟性，操縱穩(wěn)定性及安全性等已經(jīng)成為汽車開發(fā)中不可缺少的一部分，傳統(tǒng)計算流體力學(xué)仿真只考慮流體對車輛的氣動特性影響，忽略了車輛振動與外流場的耦合作用，從而導(dǎo)致得到的數(shù)據(jù)與實際情況存在一定偏差.尤其在高速行駛工況下，流固耦合效應(yīng)加劇將更加影響車輛行駛的穩(wěn)定性.因此，在計算流體力學(xué)中考慮流固耦合效應(yīng)，采用流體與固體軟件協(xié)同仿真的方法[1]，使數(shù)值仿真結(jié)果更為精準(zhǔn)，可以更好地開展車輛氣動特性和振動特性研究.對提高車輛的穩(wěn)定性、舒適性設(shè)計提供參考依據(jù).

許多研究人員進(jìn)行了流固耦合相關(guān)研究，文獻(xiàn)［2］在時域內(nèi)對剪切三角翼的流固耦合效應(yīng)進(jìn)行了研究.文獻(xiàn)［3］對汽車外后視鏡風(fēng)致振機理與特性采用流固耦合手段進(jìn)行了探究.文獻(xiàn)［4］建立了風(fēng)-車（列車）-橋耦合系統(tǒng)模型，研究環(huán)境風(fēng)對耦合系統(tǒng)的振動影響，計算結(jié)果表明流致振動現(xiàn)象嚴(yán)重影響車輛穩(wěn)定性.李田[5]利用任意拉格朗日歐拉法實現(xiàn)列車的流固耦合仿真，指出了風(fēng)激振對高速列車運行安全性的影響.當(dāng)前的流固耦合數(shù)值仿真研究對象一般為列車、機翼或建筑等，對汽車實車研究多采用單純空氣動力學(xué)、多體動力學(xué)方法或者單向流固耦合數(shù)值仿真研究.

針對以上所述關(guān)于車輛流固耦合氣動及流致振動研究的不足，本文以某實車模型為研究對象，將考慮雙向流固耦合效應(yīng)的數(shù)值仿真與傳統(tǒng)計算流體力學(xué)仿真進(jìn)行對比，探究兩者流場及氣動力差異的產(chǎn)生機理，并進(jìn)行整車風(fēng)激振特性研究.

1流固耦合理論基礎(chǔ)

流固耦合數(shù)值仿真需要將流體域與固體域聯(lián)立協(xié)同求解，兩相交界面上的數(shù)據(jù)傳遞是關(guān)鍵問題所在.由于固體域與流體域計算采用的分別是拉格朗日和歐拉描述方式，在仿真過程中，當(dāng)固體產(chǎn)生較大位移時，兩者交界面上的節(jié)點不再重合，因此仿真出錯.為解決此問題，采用任意拉格朗日歐拉（ALE）法[6]即引入?yún)⒖甲鴺?biāo)系，可以連續(xù)追蹤邊界，交界面的網(wǎng)格節(jié)點跟隨固體結(jié)構(gòu)節(jié)點運動，防止網(wǎng)格畸變. 此外，流固耦合交接面上需要滿足位移與力平衡的動力學(xué)條件［7-8］.

d_f=d_s（1）

式中：d_f，d_s為流體節(jié)點與固體節(jié)點位移.

n·τ_f=n·τ_s（2）

式中：τ_（f/s）為應(yīng)力（流體/固體）；n為單位法向量分量.

節(jié)點坐標(biāo)可以在流體中通過已知的運動學(xué)條件求得，在流固交界面上將積分計算出的合力施加到對應(yīng)固體節(jié)點上[9]：

F_（t）=∫h^dτ_f·d_s（3）

式中：h^d為固體節(jié)點位移；d_s為微元面積.

流固耦合體系里解向量為：

x=（X_f，X_s）（4）

X_s/f是耦合交界面節(jié)點上的解，分別位于固體與流體上.

已知d_s=d_s（X_s），τ_f=τ_f（X_f），耦合系統(tǒng)的有限元方程式如下所示：

式中，F(xiàn)_f是流體方程；F_s是結(jié)構(gòu)方程.

進(jìn)行流固耦合瞬態(tài)仿真計算時，固體的變形會影響流場結(jié)構(gòu)，因此在流體計算時需要用動網(wǎng)格技術(shù)來完成網(wǎng)格的更新.相關(guān)動網(wǎng)格參數(shù)設(shè)置完成后，變形區(qū)域的網(wǎng)格即會進(jìn)行自適應(yīng)變形.本文中動網(wǎng)格方法采用迭代法，主要分為彈簧光順，動態(tài)分層和局部網(wǎng)格重劃三個模型，彈簧光順模型核心思想為在計算中，假定網(wǎng)格節(jié)點之間為平衡狀態(tài)下的彈簧系統(tǒng)，節(jié)點位移之后根據(jù)Hooke定律經(jīng)過迭代確定新平衡狀態(tài)下的網(wǎng)格節(jié)點坐標(biāo).動態(tài)分層模型多用結(jié)構(gòu)化網(wǎng)格中，在網(wǎng)格運動邊界根據(jù)其運動狀態(tài)增減網(wǎng)格層數(shù).局部網(wǎng)格重劃模型主要應(yīng)用于非結(jié)構(gòu)網(wǎng)格中，對運動區(qū)域的網(wǎng)格進(jìn)插值計算進(jìn)而重新劃分.本文采用的主要為局部網(wǎng)格重劃模型.

迭代法可以較好地處理結(jié)構(gòu)大位移運動，適合車身姿態(tài)變化的情況，但不足的是由于需要進(jìn)行大量插值計算，迭代動網(wǎng)格方法會延長計算時間，降低計算效率［10-14］.

2基于流固耦合效應(yīng)的車輛數(shù)值仿真

本章考慮流固耦合效應(yīng)，對實車在流場中的氣動力及風(fēng)激振現(xiàn)象展開研究.

2.1模型及網(wǎng)格參數(shù)

采用某轎車實車模型，長度×寬度×高度為4.4m×2.0m×1.5m，正投影面積：2.18m².對車輛的細(xì)節(jié)進(jìn)行簡化及平整化處理，并將實車的懸架及輪胎等效為彈性單元，使用貼體性好的5～20mm尺寸漸變?nèi)切尉W(wǎng)格對模型進(jìn)行網(wǎng)格劃分，如圖1所示.在輪胎和車體之間加入彈性阻尼單元以代替懸架和車輪等彈性結(jié)構(gòu)，車身及其部件的形變忽略不計，并將車輪底部設(shè)置為固定端以便研究車身在流固耦合仿真中的位移情況，彈性阻尼單元根據(jù)實車條件等效剛度設(shè)為20 000 N/m，阻尼150N/（m·s^-1），固體計算采用ABAQUS軟件中瞬態(tài)隱式方法求解.

流固耦合計算對計算資源有一定要求，在保證精度的前提下，也要考慮硬件條件.流體軟件使用STAR-CCM+，并對流體網(wǎng)格進(jìn)行網(wǎng)格無關(guān)性驗證，網(wǎng)格數(shù)量與阻力系數(shù)的關(guān)系見圖2.隨著流體網(wǎng)格數(shù)量逐漸增加，可以發(fā)現(xiàn)，當(dāng)網(wǎng)格數(shù)超過900萬時，阻力系數(shù)逐漸穩(wěn)定在0.332左右，另外邊界層厚度為面網(wǎng)格尺寸的1/3，層數(shù)設(shè)定為5層.最終確定流體域尺寸設(shè)為10倍車長，9倍車寬，5倍車高[15-17]，并設(shè)置兩個加密區(qū)，最終使用20～300 mm的流體網(wǎng)格[18-22].如圖3所示，網(wǎng)格數(shù)量大約為910萬.

本文采用瞬態(tài)流固耦合數(shù)值仿真，使用k-ωSST湍流模型[23]，時間步長設(shè)為0.001 s，流體與固體的數(shù)據(jù)傳遞通過車身作為耦合交界面實現(xiàn)[24].

2.2氣動性能分析

本節(jié)對整車氣動阻力和氣動升力展開研究，如圖4所示分別為是否考慮流固耦合效應(yīng)數(shù)值仿真中阻力與升力對比圖.圖5為升阻力差異率曲線，圖6 為兩種仿真方法波動情況.圖7為速度為120km/h時（常用高速工況）兩種仿真方法的氣動升力對比. （FSI為流固耦合，NO-FSI為非流固耦合）

在高速工況下，懸架受時變氣動載荷影響以及輪胎行駛過程中發(fā)生形變，進(jìn)而車輛姿態(tài)改變量變大，改變了車輛周圍外流場結(jié)構(gòu)，氣動力受此影響繼而改變.

由圖4圖5可以看出，是否考慮流固耦合效應(yīng)對車輛氣動阻力影響較小，不同車速下差異率均在4%以下并且流固耦合數(shù)值仿真阻力值普遍較傳統(tǒng)數(shù)值仿真偏小，相反地，流同耦合對氣動升力的影響較氣動阻力更加嚴(yán)重，數(shù)值仿真結(jié)果明顯小于傳統(tǒng)CFD數(shù)值仿真，隨車速增大，差異率增大并且都在5%以上，最大可達(dá)到38.9%.氣動升力直接影響車輛對地面的附著力，由圖5圖6可知，考慮流固耦合效應(yīng)數(shù)值仿真變化與傳統(tǒng)數(shù)值仿真相比更加劇烈，因此流固耦合CFD更能真實反映車輛行駛情況，為車輛穩(wěn)定性設(shè)計提供參考.

為了更加明顯體現(xiàn)氣動力改變的原因，使用車頭下移到最大時車底部速度矢量來對比分析兩種仿真方法結(jié)果，如圖7所示.

傳統(tǒng)CFD仿真結(jié)果車底部流速保持在一個穩(wěn)定的值，而流固耦合CFD仿真，車底部速度隨時間變化，由圖8可知，F(xiàn)SI車底部流速大于傳統(tǒng)仿真流速. 由伯努利原理，流速變大導(dǎo)致壓力變小，即為流固耦合仿真升力小于非耦合仿真的原因.

2.3風(fēng)激振特性分析

車身振動屬于流致振動中的風(fēng)激振，車體在流場中，雷諾數(shù)大于特定值時，會產(chǎn)生周期性的渦脫，作用在車體使其產(chǎn)生振動.主要分析車輛的氣動升力頻譜以及車身振動頻譜.圖9為車輛氣動升力頻譜分析，圖10為車身振動頻譜分析.

由圖9可知，氣動升力頻譜有兩階頻率，一階為氣動升力與車體的激勵頻率，隨流固耦合效應(yīng)而發(fā)生改變.車速在100km/h以上時振動頻率逐漸下降至5.4Hz，其他工況下一階振動頻率基本保持穩(wěn)定在5.6Hz，幅值則隨車速逐漸增加.二階為氣動升力自振頻率，隨車速增加，二階振動幅值先增后減而頻率一直呈現(xiàn)增加的趨勢.

由圖10可知，整車一階頻率為固有頻率，二階頻率為風(fēng)激振頻率.對整車進(jìn)行模態(tài)分析得到車輛振動頻率為4.86 Hz，對比仿真得到的一階頻率4.88 Hz，差別僅為0.4%.可以證明流固耦合仿真方法準(zhǔn)確性.車身振動幅值均隨車速增加而增大，一階固有頻率不隨車速變化而改變.對比圖4圖5、圖9圖10，當(dāng)車速較低小于60 km/h時，氣動升力較小導(dǎo)致流固耦合效應(yīng)較弱，此時氣動升力體現(xiàn)出自振特性，頻率距離車身一階固有頻率較遠(yuǎn)，因此振幅較小.隨著車速提高，氣動升力使車身姿態(tài)改變量增大，流固耦合效應(yīng)增強，車身振幅增大，該情況下幅值影響因素主要為激勵作用力大小.當(dāng)速度繼續(xù)增大到100 km/h以上時，固有頻率與氣動升力頻率接近，升力振幅與車身振幅都迅速增大，因為結(jié)構(gòu)阻尼影響，最終車身振幅達(dá)到一定閾值并趨于穩(wěn)定，在此工況下，外界激勵頻率成為決定振幅的主要因素.綜上，如果車輛高速工況下，長時間車身處在強烈振動下，會影響車輛氣動性能和造成部件的疲勞破壞.因此根據(jù)考慮流固耦合效應(yīng)的數(shù)值仿真方法可以準(zhǔn)確獲得車輛各種工況下流致振動情況，為進(jìn)一步優(yōu)化汽車結(jié)構(gòu)造型方面提供參考.

2.4汽車風(fēng)激振的改善措施

車輛風(fēng)致振動會在一定程度上影響車輛穩(wěn)定性和乘員的舒適性，因此從車輛內(nèi)部的彈性結(jié)構(gòu)入手，尋找最合適的懸架參數(shù)改善風(fēng)激振現(xiàn)象.

表1為人體舒適度與加權(quán)加速度均方根a_w的評價標(biāo)準(zhǔn)[25]，a_w為加權(quán)加速度均方根值.

如圖11所示，當(dāng)車速小于100 km/h，氣動力耦合效應(yīng)較弱，人體無不舒適反應(yīng)，當(dāng)車速逐漸增大至120 km/h，流固耦合效應(yīng)增強，車輛與流場產(chǎn)生共振，隨著振動強度的增加，人體不適感也隨之增加. 當(dāng)車速繼續(xù)增加，流固耦合效應(yīng)更加劇烈，振動強度增加導(dǎo)致人體不適感大幅增加.

隨著車輛高速化趨勢發(fā)展，流致振動（風(fēng)激振）對乘員舒適度和安全性影響越發(fā)增強，因此對此類振動尋找抑制的手段是有必要的.本節(jié)僅考慮懸架彈性剛度的對風(fēng)激振情況的影響，從而得到不同剛度下的a_w值，如圖12所示.

由圖12可知，通過改變彈性元件剛度可有效提高人體舒適度.由于僅考慮了風(fēng)激振的情況，如果要獲得最佳懸架參數(shù)，需對其他因素/如路面，發(fā)動機激勵）進(jìn)行綜合考慮.

3結(jié)論

1）流固耦合效應(yīng)對氣動升力影響嚴(yán)重，對氣動阻力影響較小，對二者的影響均隨車速增加而增大.

2）流固耦合數(shù)值仿真對氣動升力的預(yù)測更貼近真實情況，可以為車輛的穩(wěn)定性和安全性提供更準(zhǔn)確的參考.

3）模態(tài)分析和流固耦合仿真分別得到的車身振動頻率差異率僅為0.4%，證明流固耦合仿真的準(zhǔn)確性.影響車身振動的主要因素為外界激勵頻率及作用力大小，獲得各工況下的車輛風(fēng)激振特性，可以對車身結(jié)構(gòu)進(jìn)行進(jìn)一步優(yōu)化設(shè)計，依此對振動情況采取改善措施并提高乘員的舒適性.

參考文獻(xiàn)

[1]朱文峰，林佩劍，周輝.高速流固耦合效應(yīng)下車窗密封機理建

模與分析[J].汽車工程，2015，37（12）：1395-1399.

ZHU W F，LIN P J，ZHOU H. Modeling and analysis of window sealing mechanism under high-speed fluid-solid coupling effect [J]. Automotive Engineering，2015，37（12）：1395-1399. （In Chinese）

[2] KUMAR A A，MANOJ N，ONKAR A K，et al. Fluid-structure interaction analysis of a cropped delta wing [J]. Procedia Engineering，2016，144：1205-1212 .

[3]蘭巍，和生泰，胡興軍，等.汽車外后視鏡風(fēng)致振機理與特性[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2019，46（10）：10-18.

LAN W，HE S T，HU X J，et al. Wind-induced vibration mechanism and characteristics of automobile exterior mirrors [J]. Journal of Hunan University（Natural Sciences），2019，46（10）：10-18.（In Chinese）

[4]李永樂，向活躍，強士中.風(fēng)-列車-橋系統(tǒng)耦合振動研究綜述

[J].中國公路學(xué)報，2018，31（7）：24-37.

LI Y L，XIANG H Y，QIANG S Z. Review on coupling vibration of wind-vehicle-bridge systems[J]. China Journal of Highway and Transport，2018，31（7）：24-37. （In Chinese）

[5]李田.高速列車流固耦合計算方法及動力學(xué)性能研究[D].成都：西南交通大學(xué)，2012.

LI T. Approaches and dynamic performances of high-speed train fluid-structure [ D] . Chengdu：Southwest Jiaotong University，2012 （In Chinese）

[6] KUHL E，HULSHOFF S，DE BORST R. An arbitrary Lagrangian Eulerian finite-element approach for fluid-structure interactionphenomena [J]. International Journal for Numerical Methods in Engineering，2003，57（1）：117-142.

[7]錢若軍，董石麟，袁行飛.流固耦合理論研究進(jìn)展[J].空間結(jié)構(gòu)，2008，14（1）：3-15.

QIAN R J，DONG S L，YUAN X F. Advances in research on fluidstructure interaction theory[J]. Spatial Structures，2008，14 （1）：3-15.（In Chinese）

[8]蘇波，錢若軍，袁行飛.流固耦合界面信息傳遞理論和方法研究進(jìn)展[J].空間結(jié)構(gòu)，2010，16（1）：3-10.

SU B，QIAN R J，YUAN X F. Advances in research on theory and method of data exchange on coupling interface for FSI analysis [J]. Spatial Structures，2010，16（1）：3-10.（In Chinese）

[9]陶莉莉.基于流固耦合的高速客車氣動特性研究[D].濟南：山東大學(xué)，2014.

TAO L L. Research on aerodynamic characteristics of high speed bus based on fluid-structure interaction[D]. Jinan：Shandong Uni- versity，2014.（In Chinese）

[10] HELMIG J，KEY F，BEHR M，et al. Combining boundaryconforming finite element meshes on moving domains using a sliding mesh approach[J]. International Journal for Numerical Methods in Fluids，2021，93（4）：1053-1073.

[11] FOSSATI M，KHURRAM R A，HABASHI W G. An ALE mesh movement scheme for long-term in-flight ice accretion[J]. International Journal for Numerical Methods inFluids，2012，68（8）：958-976.

[12] SEVILLA R，GIL A J，WEBERSTADT M. A high-order stabilised ALE finite element formulation for the Euler equations on deformable domains[J]. Computers & Structures，2017，181：89-102.

[13] RAZZAGHI M M，MIRSAJEDI S M. A moving mesh method with defining deformable layers [J]. Progress in Computational Fluid Dynamics，an International Journal，2017，17（2）：63.

[14] PRUST L J. Moving and reactive boundary conditions in movingmesh hydrodynamics[J]. Monthly Notices of the Royal Astronomical Society，2020，494（4）：4616-4626.

[15]康順.計算域?qū)FD模擬結(jié)果的影響[J].工程熱物理學(xué)報，2005，26（S1）：57-60.

KANG S. Influence of computational domain on the cfd results[J]. Journal of Engineering Thermophysics，2005，26（S1）：57-60.（In Chinese）

[16] SONG K S，KANG S O，JUN S O，et al. Aerodynamic design optimization of rear body shapes of a Sedan for drag reduction [J]. International Journal of Automotive Technology，2012，13（6）：905-914.

[17] KANG S O，JUN S O，PARK H I，et al. Actively translating a rear diffuser device for the aerodynamic drag reduction of a passenger car[J]. International Journal of Automotive Technology，2012，13（4）：583-592.

[18]楊博，傅立敏.轎車外流場網(wǎng)格生成策略及數(shù)值模擬[J].農(nóng)業(yè)機械學(xué)報，2007，38（4）：8-11.

YANG B，F(xiàn)U L M. Mesh generation strategies of the external flow field around a Sedan and the numerical simulation research[J]. Transactions of the Chinese Society for Agricultural Machinery，2007，38（4）：8-11.（In Chinese）

[19]王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社，2004.

WANG F J. Computational fluid dynamics analysis-CFD software principle and application [M]. Beijing：Tsinghua University Press，2004.（In Chinese）

[20]傅立敏.汽車空氣動力學(xué)數(shù)值計算[M].北京：北京理工大學(xué)出版社，2000.

FU L M. Numerical calculation of automobile aerodynamics[M]. Beijing：Beijing Institute of Technology Press，2000.（In Chinese）

[21]陳矛章.粘性流體動力學(xué)基礎(chǔ)[M].北京：高等教育出版社，2004.

CHEN M Z. Fundamentals of viscous fluid dynamics[M]. Beijing：Higher Education Press，2004.（In Chinese）

[22]張兆順，崔桂香，許春曉.湍流理論與模擬[M].北京：清華大學(xué)出版社，2005.

ZHANG Z S，CUI G X，XU C X. Theory and modeling of turbulence [M]. Beijing：Tsinghua University Press，2005. （In Chinese）

[23]祝志文，袁濤，陳政清，等.基于CFD仿真的平板非線性氣動力系統(tǒng)特征研究[J].湖南大學(xué)學(xué)報（自然科學(xué)版），2017，44（1）：32-38.

ZHU Z W，YUAN T，CHEN Z Q，et al. Investigation on characteristics of nonlinear aerodynamic system of thin plate based on CFD simulations[J]. Journal of Hunan University（Natural Sciences），2017，44（1）：32-38.（In Chinese）

[24] SCHROECK D，KRANTZ W，WIDDECKE N，et al. Unsteady aerodynamic properties of a vehicle model and their effect on driver and vehicle under side wind conditions[J]. SAE International Journal of Passenger Cars - Mechanical Systems，2011，4（1）：108-119.

[25]國家質(zhì)量監(jiān)督檢驗檢疫總局，中國國家標(biāo)準(zhǔn)化管理委員會.汽車平順性試驗方法：GB/T4970—2009[S].北京：中國標(biāo)準(zhǔn)出版社，2010.

General Administration of Quality Supervision，Inspection and Quarantine of the People's Republic of China，Standardization Administration of the People's Republic of China. Method of running test—Automotive ride comfort：GB/T 4970—2009[S]. Beijing：Standards Press of China，2010. （In Chinese）