劉樹成　邢慶坤　李為薇　張　靜

1.中國(guó)北方車輛研究所車輛傳動(dòng)重點(diǎn)實(shí)驗(yàn)室，北京，1000722.武漢大學(xué)，武漢，430072

?

車用液力變矩器混合流道CFD仿真方法

劉樹成1邢慶坤1李為薇2張靜1

1.中國(guó)北方車輛研究所車輛傳動(dòng)重點(diǎn)實(shí)驗(yàn)室，北京，1000722.武漢大學(xué)，武漢，430072

針對(duì)車用液力變矩器復(fù)雜動(dòng)態(tài)過程中工作相位隨時(shí)發(fā)生轉(zhuǎn)換，不能及時(shí)判斷相應(yīng)流場(chǎng)結(jié)構(gòu)的改變，難以對(duì)瞬時(shí)流場(chǎng)特性進(jìn)行準(zhǔn)確仿真的問題，基于傳統(tǒng)變矩器CFD流道模型和導(dǎo)輪空轉(zhuǎn)無葉片模型，建立了液力變矩器混合流道CFD仿真模型。該仿真模型可以自動(dòng)識(shí)別變矩器變矩、偶合和功率反傳等工作相位及其相位轉(zhuǎn)換過程，并根據(jù)導(dǎo)輪是否空轉(zhuǎn)自動(dòng)選擇相應(yīng)流道模型。對(duì)某變矩器進(jìn)行了一系列穩(wěn)態(tài)通用特性和動(dòng)態(tài)特性的仿真與試驗(yàn)研究，對(duì)比結(jié)果表明，液力變矩器混合流道CFD仿真方法對(duì)變矩器穩(wěn)態(tài)和動(dòng)態(tài)特性仿真精度較高，有效解決了變矩器復(fù)雜動(dòng)態(tài)過程難以快速實(shí)時(shí)仿真的問題，具有一定的工程實(shí)際意義。

流體力學(xué)；液力變矩器；CFD；混合流道

0　引言

液力變矩器是車輛液力傳動(dòng)系統(tǒng)的重要組成部分，其工作過程一般分為三個(gè)階段[1]：速比小于偶合工況轉(zhuǎn)換點(diǎn)，為變矩工況；速比大于偶合工況轉(zhuǎn)換點(diǎn)但小于1，為偶合工況；速比大于1，為功率反傳工況。變矩器工況的轉(zhuǎn)換稱為換相，是變矩器內(nèi)部常見的過渡工況。目前，液力變矩器內(nèi)流場(chǎng)CFD數(shù)值仿真方法在研究其穩(wěn)態(tài)和動(dòng)態(tài)特性研究中應(yīng)用廣泛[2]。但在復(fù)雜動(dòng)態(tài)過程中，變矩器內(nèi)部相位轉(zhuǎn)換頻繁，難以及時(shí)做出判斷，因而不能快速根據(jù)流場(chǎng)結(jié)構(gòu)的改變，對(duì)瞬時(shí)流場(chǎng)特性進(jìn)行準(zhǔn)確仿真。

由于在變矩工況，變矩器各工作輪轉(zhuǎn)速均為已知，所以進(jìn)行內(nèi)流場(chǎng)CFD仿真較為方便，而對(duì)于偶合工況和功率反傳工況，導(dǎo)輪處于空轉(zhuǎn)狀態(tài)，且轉(zhuǎn)速未知，進(jìn)行內(nèi)流場(chǎng)CFD仿真較為困難。對(duì)于導(dǎo)輪空轉(zhuǎn)工況，文獻(xiàn)[3-4]采用零轉(zhuǎn)矩試算法(簡(jiǎn)稱試算法)進(jìn)行仿真，即設(shè)定初值-計(jì)算-修改初值-計(jì)算……多次反復(fù)直到導(dǎo)輪轉(zhuǎn)矩接近于零，但流體網(wǎng)格數(shù)目較多或者瞬態(tài)求解時(shí)，每次試算耗時(shí)很長(zhǎng)，多次反復(fù)試算工作冗繁；文獻(xiàn)[5]在分析了導(dǎo)輪軸面和流面導(dǎo)流作用的基礎(chǔ)上，基于導(dǎo)輪空轉(zhuǎn)時(shí)流面導(dǎo)流作用近似為零，提出了導(dǎo)輪無葉片仿真方法，有效解決了導(dǎo)輪空轉(zhuǎn)工況的快速仿真問題。

本文針對(duì)變矩器復(fù)雜動(dòng)態(tài)過程的多工況轉(zhuǎn)換問題，將傳統(tǒng)變矩器流道模型和導(dǎo)輪無葉片模型有機(jī)結(jié)合，并根據(jù)仿真和實(shí)驗(yàn)數(shù)據(jù)制定模型選擇方法，建立變矩器混合流道模型進(jìn)行多工況動(dòng)態(tài)過程仿真，并進(jìn)行試驗(yàn)驗(yàn)證。

1　基本原理

(1)

(2)

2　偶合工況轉(zhuǎn)換點(diǎn)

由于變矩器的通用特性，變矩器導(dǎo)輪開始空轉(zhuǎn)的速比(偶合工況轉(zhuǎn)換點(diǎn))io隨泵輪轉(zhuǎn)速的變化而產(chǎn)生變化，因此需要首先得到io隨泵輪轉(zhuǎn)速nB變化的規(guī)律。

導(dǎo)輪無葉片和有葉片時(shí)的兩種流道模型如圖1所示。本文首先用常規(guī)流道模型，即導(dǎo)輪有葉片流道模型(圖1b)進(jìn)行仿真，以變矩器變矩比K開始變?yōu)?、導(dǎo)輪扭矩接近于0時(shí)的速比作為導(dǎo)輪開始空轉(zhuǎn)的速比點(diǎn)，并與從液力變矩器試驗(yàn)數(shù)據(jù)中確定的偶合器工況轉(zhuǎn)換點(diǎn)進(jìn)行對(duì)比驗(yàn)證，如圖2所示。

(a)導(dǎo)輪無葉片(b)導(dǎo)輪有葉片圖1　兩種液力變矩器流場(chǎng)仿真模型[5]

圖2　導(dǎo)輪開始空轉(zhuǎn)的偶合器工況轉(zhuǎn)換點(diǎn)與泵輪轉(zhuǎn)速的關(guān)系

圖2中給出了泵輪轉(zhuǎn)速范圍為500～3000 r/min內(nèi)仿真得到的偶合工況轉(zhuǎn)換點(diǎn)，還給出了泵輪轉(zhuǎn)速范圍為1000～2100 r/min內(nèi)由試驗(yàn)數(shù)據(jù)確定的偶合器工況轉(zhuǎn)換點(diǎn)，從圖中可見，隨著泵輪轉(zhuǎn)速的增大，偶合器工況轉(zhuǎn)換點(diǎn)io有所增加，雖然試驗(yàn)數(shù)據(jù)有限且由于試驗(yàn)中受溫度變化、油液黏度變化、測(cè)量誤差等因素影響使偶合器工況轉(zhuǎn)換點(diǎn)io產(chǎn)生一定波動(dòng)，但仿真曲線仍能與之基本吻合，因此可以認(rèn)為仿真模型計(jì)算得到的曲線函數(shù)就是所要求的io=io(nB)，其線性擬合公式為

(3)

3　混合流道CFD仿真模型

在研究變矩器動(dòng)態(tài)特性時(shí)，在導(dǎo)輪開始空轉(zhuǎn)前后，用單一傳統(tǒng)流道模型難以有效仿真變矩器全部運(yùn)轉(zhuǎn)特性。因此，基于式(3)確定的偶合工況轉(zhuǎn)換規(guī)律，將圖1所示傳統(tǒng)變矩器流道模型和導(dǎo)輪無葉片模型有機(jī)結(jié)合，建立車用液力變矩器混合流道CFD模型，模型結(jié)構(gòu)如圖3所示。其運(yùn)行過程：在用CFD模型進(jìn)行變矩器內(nèi)流場(chǎng)動(dòng)態(tài)特性仿真時(shí)，在每一個(gè)時(shí)間步內(nèi)，按式(3)和瞬時(shí)nB確定出偶合工況轉(zhuǎn)換點(diǎn)io，按照io進(jìn)行模型劃分，在i≤io時(shí)采用常規(guī)流道模型，即導(dǎo)輪有葉片流道模型(圖1b)，在i>io時(shí)采用導(dǎo)輪無葉片流道模型(圖1a)，并用CFD仿真模型計(jì)算得到液力扭矩與負(fù)載扭矩、機(jī)械損失扭矩和圓盤摩擦扭矩一起計(jì)算下一時(shí)間步開始時(shí)的泵輪和渦輪角加速度和轉(zhuǎn)速，直到設(shè)定的全部時(shí)間步走完為止。

圖3　液力變矩器混合流道模型動(dòng)態(tài)特性仿真模型架構(gòu)

基于Isight集成軟件平臺(tái)將上述混合流道的模型架構(gòu)集成為變矩器混合流道集成仿真平臺(tái)(圖4)，該平臺(tái)可以根據(jù)設(shè)定的動(dòng)態(tài)工況進(jìn)行變矩器動(dòng)態(tài)特性實(shí)時(shí)仿真。

圖4　基于Isight的變矩器混合流道集成仿真平臺(tái)

4　仿真實(shí)例與試驗(yàn)驗(yàn)證

以某三元件綜合式液力變矩器為研究對(duì)象，其有效直徑400 mm，采用的工作油密度為860 kg/m3，動(dòng)力黏度為0.0125 Pa·s。在泵輪恒定轉(zhuǎn)速(2000 r/min)下，渦輪初始載荷為0，對(duì)渦輪輸出軸按照設(shè)定的程序依次進(jìn)行平穩(wěn)加載，平穩(wěn)卸載，急劇加載、急劇卸載工況下的性能仿真和試驗(yàn)。

參照液力變矩器性能試驗(yàn)標(biāo)準(zhǔn)GB/T 7680-2005，設(shè)計(jì)圖5所示的變矩器性能試驗(yàn)臺(tái)，依次進(jìn)行變矩器通用穩(wěn)態(tài)特性試驗(yàn)和設(shè)定的動(dòng)態(tài)特性試驗(yàn)。在通用穩(wěn)態(tài)特性試驗(yàn)中，設(shè)定一系列泵輪轉(zhuǎn)速1000, 1200，…，2100 r/min，在每組恒定的泵輪轉(zhuǎn)速下，渦輪轉(zhuǎn)速按照速比i分別為0，0.1,0.2，…，1設(shè)定，依次測(cè)量變矩器輸入輸出的轉(zhuǎn)速和扭矩，并整理出變矩器原始特性。在動(dòng)態(tài)特性試驗(yàn)中，在恒定的泵輪轉(zhuǎn)速2000 r/min下，渦輪轉(zhuǎn)速?gòu)乃俦?.96開始，按照斜坡時(shí)間依次為30 s、10 s進(jìn)行減速到0，或者從速比0開始，按照斜坡時(shí)間依次為30 s、10 s進(jìn)行加速到速比0.96，采樣頻率80 Hz，記錄該動(dòng)態(tài)過程中輸入輸出的轉(zhuǎn)速和扭矩，并整理出動(dòng)態(tài)原始特性。

(a)試驗(yàn)臺(tái)布局圖

(b)性能試驗(yàn)臺(tái)架圖5　液力變矩器動(dòng)態(tài)性能試驗(yàn)臺(tái)

利用圖4所示的液力變矩器混合流道集成仿真平臺(tái)進(jìn)行相同過程下的穩(wěn)態(tài)與動(dòng)態(tài)性能仿真。

現(xiàn)將CFD仿真結(jié)果與相應(yīng)的試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證。其中穩(wěn)態(tài)通用特性仿真結(jié)果與試驗(yàn)數(shù)據(jù)的對(duì)比如圖6所示，部分動(dòng)態(tài)特性CFD仿真結(jié)果與試驗(yàn)數(shù)據(jù)對(duì)比結(jié)果如圖7所示，圖中nB、nT分別為泵輪和渦輪轉(zhuǎn)速；TB、TT分別為泵輪和渦輪扭矩；i、K、η分別為無因次的速比、變矩比和效率；下標(biāo)中，sim表示仿真結(jié)果，test表示試驗(yàn)結(jié)果。

(a)變矩器效率

(b)變矩器變矩比圖6　液力變矩器通用特性仿真與試驗(yàn)對(duì)比

從圖6中可見，仿真數(shù)據(jù)與試驗(yàn)數(shù)據(jù)絕大多數(shù)吻合良好，僅在速比接近于1的狹窄過渡工況仿真相對(duì)誤差較大。其中效率預(yù)測(cè)的最大相對(duì)誤差為11.9%，變矩比預(yù)測(cè)的最大相對(duì)誤差為12.7%，且在較常用的牽引工況，仿真預(yù)測(cè)的效率和變矩比相對(duì)誤差均小于10%，可見變矩器特性混合流道仿真模型對(duì)穩(wěn)態(tài)通用特性的預(yù)測(cè)精度較高。

由圖7所示的加減速動(dòng)態(tài)特性對(duì)比結(jié)果可見，本文所建立的混合流道仿真模型能夠很好地預(yù)測(cè)液力變矩器在動(dòng)態(tài)過程中的性能。取相應(yīng)時(shí)間點(diǎn)的仿真值與試驗(yàn)值比較，工作輪轉(zhuǎn)速預(yù)測(cè)的最大相對(duì)誤差為13.6%，液力扭矩預(yù)測(cè)的最大相對(duì)誤差為11.1%，可見模型對(duì)動(dòng)態(tài)特性預(yù)測(cè)的精度較高。

從上述動(dòng)態(tài)性能仿真與試驗(yàn)數(shù)據(jù)中，各自提取變矩器的原始特性曲線，對(duì)比結(jié)果如圖8所示，圖中仿真模型預(yù)測(cè)的變矩器動(dòng)態(tài)原始特性與試驗(yàn)吻合良好，這表明本文所建立的變矩器混合流道仿真平臺(tái)具有較好的可信度。

(a)轉(zhuǎn)速

(b)扭矩

(c)速比、變矩比和效率圖7　設(shè)定動(dòng)態(tài)工況的仿真與試驗(yàn)時(shí)間歷程對(duì)比

圖8　變矩器動(dòng)態(tài)原始特性仿真與試驗(yàn)結(jié)果

5　結(jié)束語(yǔ)

液力變矩器混合流道模型適用于變矩器和導(dǎo)輪空轉(zhuǎn)兩種工況仿真時(shí)，可按照相位切換原則，靈活選用具體的工況模型，即可簡(jiǎn)化變矩器動(dòng)態(tài)過程的CFD仿真過程，提高運(yùn)行效率。

無論是液力變矩器穩(wěn)態(tài)特性還是動(dòng)態(tài)特性，用液力變矩器混合流道模型的仿真結(jié)果與試驗(yàn)數(shù)據(jù)均能吻合良好。這證明了本文所建立的液力變矩器混合流道仿真模型具有較好的性能預(yù)測(cè)精度，有效解決了復(fù)雜動(dòng)態(tài)過程中變矩器內(nèi)流場(chǎng)結(jié)構(gòu)改變時(shí)的性能仿真問題。

[1]Pritchard E G D. Torque Converter Interactions in a Parallel Post Transmission Hybrid Drivline[D]. Raleigh:North Carolina State University,2010.

[2]劉春寶，馬文星，朱喜林. 液力變矩器三維瞬態(tài)流場(chǎng)計(jì)算[J]. 機(jī)械工程學(xué)報(bào)，2010,46(14)：161-166.

Liu Chunbao, Ma Wenxing, Zhu Xilin.3D Transient Calculation of Internal Flow Field for Hydrodynamic Torque Converter[J].Journal of Mechanical Engineering,2010,46(14): 161-166.

[3]才委. 雙渦輪液力變矩器轉(zhuǎn)矩分配特性研究[D]. 長(zhǎng)春：吉林大學(xué), 2009.[4]魏巍，閆清東，馬越. 液力變矩-減速裝置制動(dòng)性能流場(chǎng)數(shù)值模擬研究[C]//中國(guó)汽車工程學(xué)會(huì)越野車技術(shù)分會(huì)2008年學(xué)術(shù)年會(huì)論文集.西安, 2008:212-215.[5]劉樹成，魏巍，閆清東，等. 液力變矩器導(dǎo)輪空轉(zhuǎn)特性無葉片數(shù)值仿真[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版),2013, 43(1):22-27.

Liu Shucheng,Wei Wei,Yan Qingdong,et al. Simulation Method of Idling Characteristic of Hydrodynamic Torque Converter Stator Wheel[J].Journal of Jilin University(Engineering and Technology Edition), 2013, 43(1) :22-27.

[6]Versteeg H K, Malalasekera W. An Introduction to Computational Fluid Dynamics: the Finite Volume Method[M]. Harlow:Pearson Education, 2007.

[7]劉樹成，潘鑫，魏巍，等. 基于復(fù)雜性測(cè)度的變矩器流場(chǎng)仿真湍流模型穩(wěn)健性分析[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版),2013, 43(3):613-618.

Liu Shucheng, Pan Xin, Wei Wei, et al. Complexity-based Robustness Analysis of Turbulence Model in Torque Converter Flow Field Simulation[J]. Journal of Jilin University(Engineering and Technology Edition), 2013, 43(3):613-618.

[8]Menter F R. Two-equation Eddy-viscosity Turbulence Models for Engineering Applications[J]. AIAA Journal,1994,32(8):1598-1605.

[9]Menter F R. Review of the Shear-stress Transport Turbulence Model Experience from an Industrial Perspective[J]. International Journal of Computational Fluid Dynamics, 2009, 23(4): 305-316.

(編輯郭偉)

A Mixed Flow Channel CFD Simulation Method for Vehicle Torque Converter

Liu Shucheng1Xing Qingkun1Li Weiwei2Zhang Jing1

1.Science and Technology on Vehicle Transmission Laboratory， China North Vehicle Research Institute，Beijing,100072 2.Wuhan University,Wuhan,430072

In a complex dynamic process of hydrodynamic converter for vehicle, the work phase shift might occur at any time. Without identifying the changes of flow field in time, it could become difficult to simulate the transient flow performance. Based on traditional converter CFD channel model and bladeless channel model, a new hybrid-flow-channel CFD simulation model of hydrodynamic torque converter was established. This model might automatically identify the work phase and phase transition of torque changing, coupling and power reverse transfer. According to whether the stator wheel was idling or not, it automatically selected the appropriate model channel. A series of simulations and experimental studies of steady state and dynamic characteristics of a torque converter were conducted to validate the credibility of hybrid-flow-channal simulation model. The results show that the CFD simulation method of hybrid channels of the hydrodynamic torque converter is of high accuracy for the steady and dynamic simulation of the converter.

fluid mechanics; hydrodynamic torque converter; computational fluid dynamics(CFD);mixed flow channel

2015-06-01

產(chǎn)品創(chuàng)新國(guó)防基礎(chǔ)研究項(xiàng)目(VTDP-3101)

TH137.332DOI：10.3969/j.issn.1004-132X.2016.07.023

劉樹成，男，1986年生。中國(guó)北方車輛研究所車輛傳動(dòng)重點(diǎn)實(shí)驗(yàn)室工程師、博士。主要研究方向?yàn)檐囕v傳動(dòng)系統(tǒng)理論與技術(shù)。邢慶坤，男，1983年生。中國(guó)北方車輛研究所傳動(dòng)技術(shù)部工程師。李為薇，女，1994年生。武漢大學(xué)動(dòng)力與機(jī)械學(xué)院本科生。張靜，女，1988年生。中國(guó)北方車輛研究所傳動(dòng)技術(shù)部工程師。