混合動力汽車行星耦合傳動系統(tǒng)的圖論建模及動力學分析*

楊亞聯(lián)，米 嬌，胡曉松，秦大同

(1.重慶大學，機械傳動國家重點實驗室，重慶 400030;2.北京理工大學機械與車輛學院，北京 100081)

前言

圖論是數(shù)學拓撲學的一個重要分支，它在社會科學、自然科學等各領域均有很多應用。文獻［1］中以數(shù)目合成理論為基礎進行自運動鏈合成行星齒輪系，所得結果雖不完整，但開啟了圖論理論應用于行星齒輪構造合成與分析的先河，有關行星齒輪系的研究方法甚多，其中以圖論方法最具有一般性與系統(tǒng)性。而行星齒輪系由于結構緊湊，易于實現(xiàn)多動力耦合的優(yōu)點，在混合動力汽車動力傳動系統(tǒng)中得到了廣泛的應用。

國內(nèi)學者應用圖論對行星齒輪系進行了大量研究，包括行星齒輪的構造合成、運動分析和傳動效率分析:其中，文獻［2］中將傳動方案結構的可行性判別轉化為圖的平面性判別，從而排除構件干涉的傳動簡圖，這種方法主要針對的是已知傳動比，如何確定行星齒輪變速器的傳動方案的問題;文獻［3］中建立了行星單元體和拓撲單元回路的概念，在圖論模型的基礎上提出了計算輪系自由度和傳動比的方法，主要用于分析行星輪系的傳動性能和功率流特性;文獻［4］中建立了行星變速機構的圖論模型并進行結構幾何矛盾判別，通過轉換得到構件分析圖和轉速分析圖等，但對于混合動力汽車行星耦合傳動方案的設計，該理論還須進一步完善。因此，隨著速度耦合傳動混合動力系統(tǒng)的應用，如何應用圖論理論對混合動力汽車中眾多具有行星耦合機構的混合動力傳動系統(tǒng)的方案進行建模和分析，成為一個值得繼續(xù)研究的課題。

本文中基于圖論理論，提出了具有行星耦合機構的混合動力系統(tǒng)的整體傳動方案的圖論模型，根據(jù)其鄰接矩陣輸出鄰接碼作為有效刪除同構方案的判定方法得到混合動力汽車傳動系統(tǒng)構型的圖論模型，并基于圖論理論和圖論模型進一步開展混合動力傳動系統(tǒng)的運動學和動力學分析，建立了相應的轉速和轉矩矩陣，為混合動力汽車傳動系統(tǒng)的方案設計和分析奠定了理論基礎。

1 混合動力汽車行星耦合系統(tǒng)的圖論建模

圖論的基本原理是研究能夠被抽象成為一組點和連接它們的表示特定關系的一組線所構成的圖形。所建立的圖論模型應具備構件信息、邏輯信息、層次信息和同構信息［3］?；旌蟿恿ο到y(tǒng)是指系統(tǒng)需要多種形式能量或多種動力源搭配工作的系統(tǒng)，在混合動力汽車傳動系統(tǒng)中混合動力對應為發(fā)動機和電機。為了更好地描述由發(fā)動機、電機、動力耦合機構、制動和輸出構件構成的復雜混合動力系統(tǒng)，將行星耦合混合動力系統(tǒng)模型的層次信息分為3層:首先，將行星齒輪機構的行星輪放在第Ⅰ層，稱為行星層;然后，將與太陽輪共軸回轉的其他行星齒輪機構的構件(行星架和齒圈)放在第Ⅱ層，稱為太陽層;最后，將輸入、輸出與制動構件放在第Ⅲ層，稱為輸入/輸出層，見圖1。

在建立行星耦合機構的圖論模型中，為了更有效地表示基本行星耦合機構的構造關系，以文獻［5］中發(fā)展的復接頭運動鏈圖論表示法為基礎，提出行星齒輪系各構件與圖論模型的符號轉化關系如下:

(1)行星耦合機構各個構件轉化為實心圓點;

(2)齒輪副連接轉化為虛線連接(若為齒輪內(nèi)嚙合用單箭頭虛線表示，若為齒輪外嚙合用雙箭頭虛線表示);

(3)同一回轉軸的兩個及以上的回轉副(即R、H、S共軸回轉副)轉化為通過點劃線連接;

(4)單一回轉副(即P與H的回轉副)轉化為通過雙實線連接。

基本行星耦合機構簡圖轉化為其所對應的圖論模型如圖2所示。在此圖論模型的基礎上加入輸入/輸出層，其符號轉化關系如下:

(1)輸入/輸出和制動構件轉化為空心圓點;

(2)行星排與輸入輸出構件的連接關系轉化為實線連接;

(3)通過離合器連接的轉化為雙箭頭實線。

根據(jù)上述分層結構和構件與符號的轉化關系，以豐田公司2004年的Prius混合動力汽車THS-Ⅰ傳動系統(tǒng)為例建模，其傳動系統(tǒng)構型及對應的圖論模型如圖3所示。

由于圖論模型中構件和線段的符號類型較多，不利于直觀快速分析圖論模型，提出如下簡化方法:(1)將行星排的構件和輸入/輸出層構件的符號統(tǒng)一為實心圓點;(2)將雙實線和點劃線統(tǒng)一用實線表示。簡化后的圖論模型如圖4所示。

2 混合動力汽車行星耦合系統(tǒng)圖論方案同構判定

在所有圖論表述的可能方案中，部分不同構型描述的傳動結構在運動學和動力學特性上可能實質是一樣的，因此必須對混合動力系統(tǒng)進行同構判定。本文中應用行星排鄰接矩陣和構件鄰接矩陣的組合來描述混合動力系統(tǒng)的圖論模型，通過鄰接矩陣輸出編碼來判斷傳動系統(tǒng)的同構方案。

2.1 鄰接矩陣表示法

鄰接矩陣是表示圖論中點和線的拓撲關系的有效工具，通過鄰接矩陣可以很方便地將圖論模型轉化為數(shù)學模型。首先，以圖論模型的行星層和太陽層為研究對象，建立行星耦合機構鄰接矩陣A。(1)點i和點j以單一回轉副邊連接，則aij=1;(2)點i和點j以齒輪副邊相連接，內(nèi)嚙合則aij=E，外嚙合則aij=F;(3)點i和點j以p個點共用同一回轉副連接，則 aij=p;(4)點 i和點 j不連接，則aij=0。

然后以圖論模型的太陽層和輸入/輸出層為研究對象，建立混合動力系統(tǒng)的構件鄰接矩陣B。(1)點m和點n以實線連接，則bmn=1;(2)點m和點n以雙箭頭實線連接即通過離合器連接，則bmn=2;(3)點m和點n不連接，則bmn=0。

以豐田公司Camry混合動力汽車傳動系統(tǒng)為例，如圖5所示。

根據(jù)上述圖論模型可得到行星耦合機構的鄰接矩陣A和構件鄰接矩陣B:

2.2 圖論模型的同構判定

對具有行星耦合機構的混合動力系統(tǒng)進行同構判定，首先應確定圖論模型各個層次的構件信息，本文中描述的混合動力系統(tǒng)的圖論模型的各層構件總數(shù)Ni如表1所示。

表1 不同構型圖論模型各層次的構件信息

根據(jù)各層次的構件數(shù)目就可分別確定行星耦合機構和混合動力傳動系統(tǒng)的可能方案的數(shù)目，再對圖論模型進行同構判定，首先根據(jù)行星排鄰接矩陣A對行星耦合機構進行同構判定，同構是指兩圖論模型的點與點的鄰接關系和點與邊的附隨關系是完全相同的［6］。是否同構可通過鄰接碼來判斷，由于鄰接矩陣A為上三角對稱矩陣，因此，取其矩陣上三角的所有數(shù)字依序排列的集合作為A的鄰接碼輸出，即由圖5的圖論模型得到鄰接矩陣A即式(1)，輸出的鄰接碼為［55510/555F0/55EE/50F/01/0］，然后按照下述步驟完成行星耦合機構的同構判定。

(1)列舉所有可能的位移圖形，即為除去行星齒輪圖論模型內(nèi)所有齒輪副對邊的不含任何回路的圖形;

(2)對每個位移圖形合成所有可能的含齒輪回轉對的齒輪圖形;

(3)刪除違反基本行星齒輪系的原則的圖形;

(4)由鄰接矩陣得到鄰接碼并進行同構判定;

(5)得到最終的行星齒輪系圖論方案并還原為行星齒輪機構圖。

以與Camry相同的3自由度雙排行星耦合機構為例，對表2所示位移圖形進行同構判定可知，2－2和2－3所示圖形與2－4的鄰接碼相同，因此為同構可以刪除。對此圖論模型得到的不同構最終模型為2－1和2－4，其中2－1即為Camry的行星耦合機構。同理，可以對不同自由度下的不同圖論模型的位移圖形進行同構判定。

表2 3自由度位移圖形的同構判定

然后根據(jù)構件鄰接矩陣B對混合動力傳動系統(tǒng)的可能方案進行同構判定?；旌蟿恿ο到y(tǒng)存在同構方案是由于系統(tǒng)可能存在具有相同性質的構件，例如存在雙電機或多個制動器的情況，這時須根據(jù)構件鄰接矩陣B輸出編碼來判定同構，取構件鄰接矩陣B的性質相同的構件所在行的數(shù)字構成集合為B矩陣鄰接碼，即圖5中對應的 B矩陣鄰接碼為［100000/000100］，若兩行數(shù)字對換后鄰接碼相同，則說明這兩個構型為同構，則可以刪除。

以如圖6所示的混合動力系統(tǒng)為例，取B矩陣的MG1和MG2所在行的數(shù)字集合構成鄰接碼B1、B2，B1中兩行元素互換后所得的鄰接碼若與B2相同，說明此兩種方案同構可以刪除其中之一。

3 混合動力傳動系統(tǒng)的運動學與動力學分析

3.1 基于圖論的運動學分析

有關行星耦合機構的運動學分析的研究方法有輪系值法［7］、圖解法［7］、建表法［8］和基本回路法［9］。基本回路法已成為一種系統(tǒng)的圖論理論推導方法。

圖論理論的基本回路法主要應用于行星耦合機構的圖論模型，因此對于具有單排、雙排甚至多排行星耦合機構的混合動力傳動系統(tǒng)，只須對傳動系統(tǒng)的行星耦合機構的圖論模型部分運用基本回路法即可開展運動學分析，由于基本回路法能更方便地得到行星輪的轉速與其他齒輪構件的轉速關系，因此運用此法也對加入包括行星輪在內(nèi)的各構件的轉動慣量用于動力學分析奠定了很好的基礎。

基本回路指的是圖論模型中一對齒輪副所對應的點與兩對回轉副所對應的點構成的三角形封閉回路，如圖7所示。結合本文中提出的圖論模型可知基本回路具備以下特點:(1)每個基本回路由3個實心圓點、一條實線、一條點劃線和一條虛線組成;(2)兩個相鄰回路都含有屬于同一行星耦合機構的行星架和行星輪;(3)每個基本行星耦合機構的基本回路數(shù)L必為偶數(shù)且等于帶箭頭的虛線數(shù);(4)行星耦合機構的構件數(shù)N與回轉副數(shù)M的關系為:N=M+1。

以圖2的基本行星耦合機構的圖論模型的第Ⅰ層和第Ⅱ層為研究對象，可得到單排行星耦合機構的基本回路圖，如圖7所示。

基本回路由3個基本構件1、2、3組成，其中1和2表示一對嚙合齒輪，2和3表示回轉副，則可得到該基本回路的轉速方程通式為

式中:ω1、ω2、ω3為基本構件 1、2、3 的轉速;i12為傳動比，i12=±Z2/Z1，±表示齒輪的內(nèi)嚙合或外嚙合。

以Prius的單排行星耦合混合動力系統(tǒng)為例，其圖論模型如圖3所示，結合圖7所示的單行星排基本回路的特征可知其行星耦合機構包括兩個基本回路(P，S，H)、(P，R，H)，則其基本回路方程為

雙排行星耦合混合動力系統(tǒng)以美國通用公司的AHS雙模傳動系統(tǒng)在低速模式下圖論模型為例，如圖8所示。其中發(fā)動機接第1排行星輪的行星架，電機MG1接第1排行星輪的太陽輪，第1排行星輪系的齒圈和第2排行星輪的太陽輪固聯(lián)，第2排行星輪的齒圈與機架固聯(lián)，第2排的太陽輪與電機MG2相連，第2排的行星架輸出動力到車輪。

根據(jù)上述圖論理論可知，雙行星耦合機構的基本回路分別為(P1，R1，H1)、(P1，S1，H1)、(P2，R2，H2)、(P2，S2，H2)，根據(jù)前述理論推導，可得到基本回路的運動方程組為

同時 ωr1= ωs2= ωmg2，ωs1= ωmg1，ωr2=0，ωh1=ωe，ωh2=ωout，ωs1=ωmg1，聯(lián)立并整理式(7)可得雙排行星耦合機構運動學矩陣:

由式(6)可知兩個方程解4個未知量，式(8)有4個方程，解6個未知量，行星排有兩個自由度，只要知道兩個確定的動力輸入轉速，則行星排其余所有的轉速均可確定，方程組均有解。

3.2 行星耦合混合動力系統(tǒng)的動力學分析

鑒于現(xiàn)有圖論理論行星耦合機構的動力學分析中忽略了行星耦合機構各構件的轉動慣量，也少有針對整體傳動系統(tǒng)的動力學分析，因此本文中在考慮到包含行星輪在內(nèi)的行星排各構件的轉動慣量的前提下，對行星耦合機構進行受力分析。首先，以Prius的THS單排行星耦合混合動力傳動系統(tǒng)為例，將圖4的圖論模型倒置并根據(jù)圖論的連接關系將各構件的符號一一對應的排列，則可得到如圖9所示的圖論模型。

將圖9的圖論模型分別加入行星排各構件和輸入輸出各構件的轉動慣量模塊以及受力分析即可得到如圖10所示的動力學模型。其中，設Ji為各構件的轉動慣量，ri為各齒輪半徑，F(xiàn)i為行星輪與其他齒輪之間的作用力，Ts、Tr、Th為作用在太陽輪、齒圈或行星架上的轉矩，設順時針為正。

取發(fā)動機為受力分析主體，有

取MG1電機為受力分析主體，有

取行星耦合機構后的動力傳動系統(tǒng)為受力分析主體，有

式中:m整車質量;r為車輪半徑;i為行星排至車輪之間的傳動比;Tf為車輪阻力矩折算到行星輪系輸出端的轉矩。Je、Jmg1、Jmg2分別為發(fā)動機、電機1和電機 2的慣量;ωe、ωmg1、ωmg2分別為發(fā)動機、電機 1和電機2的轉速;Te、Tmg1、Tmg2分別為發(fā)動機、電機1和電機2的轉矩;Ts、Th、Tr分別為單排行星輪系太陽輪、行星架和齒圈上的轉矩。

然后取行星輪系進行受力分析，可得THS系統(tǒng)的單排行星耦合機構的動力學方程:

又根據(jù)圖論連接關系可知:

結合上述各式，最后可得到Prius混合動力傳動系統(tǒng)的動力學矩陣:

同理，以圖8所示的GM-AHS傳動結構為例，將圖8轉換為圖11所示的動力學模型。

取發(fā)動機為受力分析主體，有

分別取電機1和電機2為受力分析主體，有

取行星耦合機構后的動力傳動系統(tǒng)為受力分析主體，有

取圖中第1排行星輪系為受力分析的主體，并進行受力分析，如圖12所示。

以行星架為參考坐標系對行星輪進行分析，由達朗伯原理，行星輪要受力平衡，須增加一虛擬的慣性力Fg1有

式中:mp1為行星輪質量;rh1為行星輪旋轉中心到行星排旋轉中心的距離。選順時針方向為正，在行星架旋轉坐標系上觀察行星輪運動:

結合式(19)～式(21)，有

分別以太陽輪、行星架和齒圈為分析對象，有

聯(lián)立式(23)～式(25)并帶入式(21)得

其中 rh1=rp1+rs1=rr1－ rp1，和式(19)、式(22)一起帶入式(26)并化簡，有

同理，取第2排行星輪系受力分析，有

因行星輪的慣量Jp相對其他行星部件的慣量小，如果忽略行星輪自轉的慣量Jp，只考慮折算到行星架上公轉的慣量，那么相當于行星架的轉動慣量增加了mpr2h，式(27)和式(28)可以進一步簡化為

聯(lián)立式(15)～式(18)、式(29)和式(30)并帶入式(31)，得

由式(32)可知，有6個方程解14個變量;結合式(8)，當2自由度行星排系統(tǒng)輸入兩個確定的輸入轉速，則其余所有的運動可求，去掉4個轉速變量;實際運行的系統(tǒng)，道路阻尼和動力源的功率均為已知，可再消去4個輸入輸出轉矩變量，因此6個方程求解6個轉矩變量，式(32)有解。同理，式(14)也有解。上述動力學模型可以完整地描述行星耦合機構和整個傳動系統(tǒng)的轉速轉矩關系。

3 結論

(1)基于圖論理論在對行星耦合機構建模的基礎上，增加了輸入/輸出構件層，建立了行星耦合混合動力系統(tǒng)的圖論模型。

(2)以Camry行星耦合混合動力系統(tǒng)為例，在行星排鄰接矩陣的基礎上增加了輸入輸出構件的構件鄰接矩陣，通過鄰接矩陣輸出鄰接碼判別傳動方案的同構模型，是基于圖論理論的混合動力系統(tǒng)方案設計的有效途徑。

(3)基于圖論理論基本回路法對行星排的運動學分析方法，分別以單排Prius-THS和雙排GM-AHS混合動力系統(tǒng)為例，建立了行星耦合混合動力系統(tǒng)的轉速分析模型;并考慮行星排轉動慣量，進行了行星排的受力分析，推導了行星耦合混合動力傳動的動力學轉速和轉矩模型，為進一步傳動方案的性能分析和篩選奠定了基礎。

［1］ Johnson R C，Towfigh K．Creat Design of Epicyclic Gear Trains U-sing Number Synthesis［J］．ASME Transactions．Journal of Engineering for Industry，1967:309 －314．

［2］ 鄧曉龍，朱大林，方子帆．三自由度串聯(lián)行星變速箱傳動方案設計的圖論法［J］．機械傳動，2002，26(2):11 －13．

［3］ 薛隆泉，雷亞榮，宋念龍．周轉輪系系統(tǒng)化的數(shù)學模型［J］．機械工程學報，2008，44(4)．

［4］ 李宏才，閆清東，李慎龍．行星傳動方案結構幾何矛盾圖論判別方法［J］．北京理工大學學報，2010，30(9)．

［5］ Yan H．S，Hsu C．H．Contracted Grahs of Kinematic Chains with Multile Joints［J］．Mathematical and Computer Modelling，1988，10(9):681－695．

［6］ Buchsbaum F，F(xiàn)reudenstein F．Synthesis of Kinematic Structure of Geared Kinematic Chain and Other Mechanisms［J］．Journal of Mechanisms，1970，5(3):357 －392．

［7］ 饒振剛．行星齒輪傳動設計［M］．北京:化學工業(yè)出版社，2003．

［8］ White G．Multiple-stage，Split-power Transmissions［J］．Journal of Mechanisms，1970，5(4):505 －520．

［9］ Freudenstein F，Yang A T．Kinematics and Statics of a Coupled Spur-gear Train［J］．Mechanism and Machine Theory，1972(7):263－275．

［10］ 秦大同，謝里陽．現(xiàn)代機械設計手冊［M］．北京:化學工業(yè)出版社，2011．

［11］ Liu Jinming．Modeling，Configuration and Control Optimization of Power-split HEVs［D］．Ann Arbor:Univ．of Michigan，2007．

［12］ Sherman Y T Lang．Graph-theoretic Modelling of Epicyclic Gear Systems［J］．Mechanism and Machine Theory，2005，40:511 －529．

［13］ Roland Mathis A，Yves Remond．Kinematic and Dynamic Simulation of Epicyclic Gear Trains［J］．Mechanism and Machine Theory，2009，44:412 －424．

［14］ 楊亞聯(lián)，馮荊松．基于圖論的行星輪系的運動性分析［J］．重慶理工大學學報，2010，24(7)．