發布時間：2016-05-12所屬分類：科技論文瀏覽：1次
摘 要： 對于多軸重型車輛，行駛路況復雜惡劣，載荷變化較大，油氣懸架以車高可調、連通平衡、剛度非線性等優越特性在多軸重型車輛領域越來越多地得到應用。本文是一篇 車輛工程論文 范文，主要論述了基于分數階的多支路油氣懸架模型研究。 摘 要：基于流體力學理論
對于多軸重型車輛，行駛路況復雜惡劣，載荷變化較大，油氣懸架以車高可調、連通平衡、剛度非線性等優越特性在多軸重型車輛領域越來越多地得到應用。本文是一篇車輛工程論文范文，主要論述了基于分數階的多支路油氣懸架模型研究。
摘 要：基于流體力學理論和理想氣體狀態方程，建立油氣懸架主要部件的模型，包括節流閥模型、溢流閥模型、蓄能器模型、油缸模型。根據系統中各部件的連通關系，建立多支路油氣懸架整數階模型。引入分數階理論，建立多支路油氣懸架分數階模型，并將分數階模型、整數階模型與臺架試驗對比，驗證了分數階模型的正確性。通過改變分數階階次，觀察系統阻尼力的變化規律，證明了分數階模型在研究剛度及阻尼非一致性問題時的優勢，為后續的多支路油氣懸架阻尼非一致性問題的研究奠定了基礎。
關鍵詞：分數階,多支路,油氣懸架,模型,階次影響
Abstract：The models had been built for major components of a hydropneumatic suspension system based on knowledge of hydromechanics and the equation of state for an ideal gas. According to the connection of components， the integer order model of multi-branch hydropneumatic suspension was built and then the fractional order model was built based on the fractional order theory. The feasibility and advantages of the fractional order theory were verified through the comparison of the bench tests on both models. The influence of change in fractional operators on the system characteristics was analyzed for the subsequent research on damping non-coincidence in multi-branch hydropneumatic suspension.
Key words：fractional order; multi-branch; hydropneumatic suspension; model; fractional influence
由于各種建模和分析方法都有其自身的特點和局限性，必須針對具體問題采用適當的方法對系統進行理論分析和參數識別。下面通過查閱文獻整理出一些比較有代表性的研究。Koenraad根據油氣獨立懸架油缸的內部結構，提出利用有關物理定律建立油氣懸架油缸初始模型 [1]。Joo在綜合考慮氣體多變過程、孔口紊流出流方程、庫侖摩擦力、油液的可壓縮性和緩沖限位塊等非線性因素后，建立了油氣懸架的非線性數學模型 [2]。Felez等人提出了一種建立多軸汽車起重機互連式油氣懸架模型的方法――功率鍵合圖法，指出在多軸起重機上油氣懸架相比傳統懸架具有諸多優勢 [3]。Yousefi 等人提出一種通過低階近似對油氣懸架非線性模型進行線性化處理的方法，通過設計一個包括模型特征參數的狀態空間，使狀態變量及其導數的近似值均達到最小[4]。
分數階微積分很早就已經被提出[5-7]，分數階微積分理論是研究任意階數的微分積分特性及應用的數學理論，是整數階微積分理論的延伸和推廣。
油氣懸架是一個多相介質的組合體，氣體隨溫度的變化、油液粘度、油液可壓縮性、外力的作用時間及加載歷史等都會影響到油氣懸架特性，這些特性十分符合粘彈性物質的特點。粘彈性材料的力學性質介于理想彈性體和牛頓流體之間，既表現出彈性，同時又具有粘性特性，其應力-應變響應依賴于時間和應變率，又與載荷和變形歷史有關，即應力-應變具有記憶性。而分數階粘彈性微積分模型可以用較少的參數來準確地描述大量復雜的粘彈性材料在廣泛頻率范圍內的動力學特性[8-10]。
1 多支路油氣懸架結構簡介
用于四軸重型車輛的互連式油氣懸架，其前兩橋與后兩橋結構相同，以一二橋的結構為例來說明。系統采用雙腔油缸，上腔為環形腔，下腔為無桿腔。一橋與二橋的同側油缸上下腔對應連通，同橋左右側油缸交叉連通。壓縮行程時，油缸下腔的油液經多條支路流向蓄能器和其它油缸。
為了實現多軸車輛的軸荷平衡及抗側傾功能，油氣懸架的各油缸之間是按照一定的規則連通的，這就不可避免存在單個油缸連通多條支路的結構。以單油缸三支路結構的油氣懸架為研究對象，對多支路油氣懸架的建模方法進行研究，單油缸三支路油氣懸架原理結構簡圖如圖1所示。
一個單作用油缸分別通過三條支路與三個蓄能器連通，每條支路中包括節流閥、壓縮溢流閥、伸張溢流閥。多支路油氣懸架系統中三條支路壓縮伸張行程油液所流過的閥見表1。
2 多支路油氣懸架模型
采用傳統的節流閥和溢流閥的數學模型，在以下建模中，以式(1)的形式來表述。
基于節流閥模型、溢流閥模型、蓄能器模型，可推導支路1數學模型如式(3)所示[11-12]。
式(1)～(4)中，QT1為支路1中節流閥T1的流量，m3/s;Q1為支路1總的流量，m3/s;dR1c為支路1中溢流閥R1-c的開口直徑，m;為支路1中溢流閥R1-c的閥芯最大位移，m;dR1r為支路1中溢流閥R1-r的開口直徑，m;為支路1中溢流閥R1-r的閥芯最大位移，m;pc為油缸內的壓力，Pa;pA1為蓄能器A1內氣體的壓力，Pa。 同理得支路2和支路3的數學模型[13]。
由油缸內油液壓力微分方程可得油缸模型如式(5)所示[14]。
式中，CHC為油缸的容積模量，m3/Pa;Vb為靜平衡位置的油缸容積，m3;Koil為油液體積模量，Pa。
綜合三條支路的模型和油缸模型，可得到單油缸多支路多蓄能油氣懸架系統模型�；谝呀⒌墓澚鏖y模型、溢流閥模型、蓄能器模型、油缸模型，采用模塊化建模方法，通過分析各部件模型之間的輸入輸出關系，建立多支路油氣懸架系統建模流程圖，如圖2所示。各支路閥組模型輸出的流量作為本支路蓄能器模型的輸入，各支路蓄能器模型輸出的壓力和油缸模型輸出的壓力作為本支路閥組模型的輸入，各閥組輸出的流量作為油缸模型的輸入。
3 多支路油氣懸架分數階模型
早先的學者對油氣懸架彈性力和阻尼力的研究，是將系統看做理想的彈性體和流體，彈性力與位移之間滿足一定的關系，阻尼力與速度之間滿足一定的關系。但是，油氣懸架是一個多相介質的組合體，既不是理想的彈性體也不是理想的流體，而是介于彈性體和流體之間的類似粘彈性體。因此，用分數階微積分來描述油氣懸架在廣泛頻率范圍內的系統特性是恰當的。
3.1 分數階微積分的定義及計算方法
分數階微積分常見的有三種定義方式：Riemann-Liouville定義、Grunwald-Letnikov定義和Caputo定義，在工程應用中最常用到的是Riemann-Liouville定義。
Riemann-Liouville分數階微分定義為
目前對分數階微積分常用的計算方法中，濾波逼近的方法效果較好，精度理想，并且可以對位置的信號進行分數階微積分，最具代表性的是由Oustaloup提出的濾波逼近分數階的方法[15-16]。研究表明，當逼近階次時，對不同的分數階次α，無論是幅頻特性還是相頻特性，都能很好地逼近實際分數階濾波器的頻率特性[17]。
3.2 多支路油氣懸架分數階模型
鑒于油氣懸架是一個多相介質的組合體，其特性不僅和外力大小有關，還和溫度的改變、力的作用時間及加載歷史都有關系，這些十分符合粘彈性物質的特點�；谝陨戏治觯梅謹惦A粘彈性微積分模型可以用較少的參數來準確地描述大量復雜粘彈性材料在廣泛頻率范圍內動力學特性的優勢，將油氣懸架系統整體等效為粘彈性系統進行研究，引入分數階模型描述彈性力和阻尼力項，建立油氣懸架分數階模型，對油氣懸架系統輸出力特性進行分析。
基于前面建立的多支路油氣懸架整數階模型，用分數階算子替換已建立的整數階模型中的位移項和速度項。觀察整數階模型的各方程，只有油缸模型的表達式中含有位移項和速度項，用分數階算子替換表達式中的位移項和速度項之后，如式(8)所示，其它部分的模型保持不變。
4 多支路油氣懸架分數階模型仿真分析
4.1 分數階階次和傳遞函數的確定
對于Oustaloup濾波逼近法，當階次時已經能滿足精度和計算量的要求，因此采用五階有理函數來逼近分數階函數�？紤]到車輛工作的有效頻率范圍，選取逼近頻率區間為[0.01 rad/s，100 rad/s]，即，，采用試錯法找出階次α和β的最優取值。階次α和β取值與對應的五階傳遞函數見表2。
4.2 分數階模型與試驗對比分析
仿真參數見表3。在1.0 Hz/0.03 m正弦激勵下，經仿真對比，當分數階模型中的階次α-0.1、β0.9時，仿真曲線與試驗曲線更接近，阻尼力仿真曲線與試驗曲線對比如圖3所示。
圖3中黑色實線為試驗測得的阻尼力曲線，紅色虛線為整數階模型仿真得到的阻尼力曲線，藍色點畫線為分數階模型仿真得到的阻尼力曲線，由兩條曲線對比可以看出，分數階模型的仿真曲線，其原點左右兩側的曲線不對稱性更明顯，阻尼力最大值不在原點處取得，分數階模型更能體現多支路油氣懸架阻尼力的實際特性。
由分數階模型仿真曲線和試驗曲線的對比，證明了油氣懸架分數階模型的可行性和準確性，體現了分數階的階次可以任意連續可調的優勢。
改變節流閥T1的節流孔直徑，觀察節流孔直徑變化對分數階模型的阻尼力-位移曲線、阻尼力-速度曲線的影響，驗證油氣懸架分數階模型的參數靈敏度，如圖4和圖5所示。
圖4為不同節流孔直徑下的分數階模型阻尼力-位移曲線，橫軸的上半部分為壓縮行程阻尼力曲線，從最低點-30 mm運動到最高點30 mm;下半部分為伸張行程阻尼力曲線，從最高點30 mm運動到最低點-30 mm。由圖4可以看出，隨著節流孔直徑的減小，分數階模型的阻尼力增大，且阻尼力曲線的左右非對稱性變得更加明顯，即壓縮行程下半段(-30～0 mm)的阻尼力小于壓縮行程上半段
(0～30 mm)的阻尼力，伸張行程上半段(30～0 mm)
的阻尼力小于伸張行程下半段(0～-30 mm)的阻尼力，即多支路油氣懸架存在阻尼非一致性問題。
圖5為不同節流孔直徑下的分數階模型阻尼力-速度曲線，第一象限為壓縮行程阻尼力曲線，第三象限為伸張行程阻尼力曲線。由圖5可知，隨著節流孔直徑的減小，分數階模型的阻尼力增大，并且在一個完整壓縮伸張行程中，阻尼力曲線不重合，中間有一定的包絡面積，隨著節流孔直徑的減小，阻尼力曲線之間的包絡面積逐漸增大，即多支路油氣懸架的阻尼力非一致性問題變得更加顯著。
4.3 分數階節次對阻尼力的影響及物理意義
經前面的分數階模型仿真與試驗對比，驗證了分數階模型的可行性與正確性。分數階的階次對模型的幅值和相位會產生影響，給系統加入1.0 Hz/30 mm的正弦激勵信號，分析分數階的階次對油缸活塞速度、阻尼力-時間曲線、阻尼力-位移曲線的影響，為分數階模型的階次選擇提供理論指導，并結合油氣懸架的實際情況，對分數階的物理意義進行闡述。 圖6為不同分數階階次下的阻尼力-時間曲線。由圖6可知，當阻尼分數階算子的階次β=1.0時，阻尼力在相位上和油缸速度是同步的，幅值在七條曲線中居中，此種情況表征了油氣懸架理想動態特性，即油氣懸架的系統輸出特性與油缸運動狀態同步。
當阻尼分數階算子的階次β<1.0時，阻尼力在相位上相對油缸速度存在延遲，隨著階次的減小，相位延遲增大，阻尼力幅值逐漸減小，此種情況表征了油氣懸架實際的動態特性，即油氣懸架的系統輸出特性與油缸運動狀態存在相位延遲。在實際中由于油液可壓縮性、油液受溫度的影響、氣體受溫度的影響、管道壓力損失、油缸與活塞的摩擦力等因素，阻尼力在相位上和油缸速度存在一定的延遲，且阻尼力幅值也會較理想情況有所減小，在[0， 1.0]區間內通過適當選取階次，可以達到相位及幅值與實際情況較好吻合的效果。分數階階次的選取體現了油液可壓縮性、油液受溫度的影響、氣體受溫度的影響、管道壓力損失等因素的綜合影響效果。當阻尼分數階算子的階次β趨于0時，相位延遲達到90°，經仿真研究，阻尼分數階算子的階次在區間[0.8， 1.0]內選取能達到與實際情況較好吻合的效果。
當階次β>1.0時，阻尼力在相位上相對油缸速度存在超前，隨著階次的增大，相位超前增大，阻尼力幅值逐漸增大。在實際中，油氣懸架輸出特性在相位上不會產生比油缸運動狀態超前的情況，因此，此種情況沒有實際的物理意義。
圖7為不同分數階階次下的阻尼力-位移曲線，零點上方為壓縮行程阻尼力曲線，零點下方為伸張行程阻尼力曲線。由圖7可知，當阻尼分數階算子的階次β=1.0時，以位移零點為分界線的阻尼力曲線左右部分基本是對稱的;當階次β<1.0時，在位移絕對值相同位置處(例如，看圖中位移為0.02 m和-0.02 m處)，壓縮行程左半部分曲線要比右半部分曲線幅值小，伸張行程左半部分曲線要比右半部分曲線幅值大，即左右部分的曲線是不對稱的，且隨著階次的減小，這種不對稱性變得更加明顯;當階次β>1.0時，沒有明確的物理意義，這里就不再討論。
圖8為不同分數階階次下的阻尼力-速度曲線，零點上方為壓縮行程阻尼力曲線，零點下方為伸張行程阻尼力曲線。當階次β=1.0時，在一個完整的壓縮和伸張行程，阻尼力-速度曲線重合為一條曲線;當階次β<1.0時，在一個完整的壓縮和伸張行程中，阻尼力-速度曲線不再重合為一條曲線，而是有一定的包絡面積，且隨著階次的減小，包絡面積逐漸增大，這就是多支路油氣懸架阻尼非一致性問題，在后續的工作中將對其進行深入研究。
5 結論
根據系統結構，采用模塊化的建模方法，建立了多支路油氣懸架整數階模型。針對油氣懸架是多相介質組合體的特點，將油氣懸架等效為粘彈性體，引入分數階微積分理論，建立了多支路油氣懸架分數階模型。將模型仿真與臺架試驗對比，并進行了后續的分析，可得到以下結論：
(1) 通過分數階模型仿真、整數階模型仿真與臺架試驗的對比，證明了多支路油氣懸架分數階模型準確性更高，更能反映多支路油氣懸架的實際特性。
(2) 改變節流孔直徑的大小，觀察多支路油氣懸架阻尼力的變化，結果表明，系統阻尼力隨節流孔直徑變化的靈敏度是滿足要求的。
(3) 改變分數階階次，當阻尼分數階算子的階次等于1時，表征了油氣懸架理想動態特性，即油氣懸架的系統輸出特性與油缸運動狀態同步;當阻尼分數階算子的階次小于1時，阻尼力在相位上相對油缸速度存在延遲，此種情況表征了油氣懸架實際的動態特性，即油氣懸架的系統輸出特性與油缸運動狀態存在相位延遲;當阻尼分數階算子的階次大于1時，阻尼力在相位上相對油缸速度存在超前，此種情況沒有實際的物理意義。
參考文獻(References)：
KOENRAAD R. A Nonlinear Parametric Model of an Automotive Shock Absorber [C]// SAE Paper 940869.
JOO F R. Dynamic Analysis of a Hydro-pneumatic Suspension System [D]. Montreal：Concordia University，1991.
FELEZ J，VERA C. Bond Graph Assisted Models for Hydro-pneumatic Suspensions in Crane Vehicles [J]. Vehicle System Dynamics，1987，16(5-6)：313-332.
YOUSEFI A，BORIS S. Order Reduction of Nonlinear Hydro-pneumatic Vehicle Suspension [C]//Proceedings of the 2006 IEEE International Conference on Control Applications，Munich，2006，1404-1408.
LEIBNIZ G W. Mathematical Schiften [M]. Hildeshcim：Olms Verlagsbuchhandlung，1962.
OLDHAM K B，SPANIER J. The Fractional Calculus [M]. New York：Academic Press，1974.
MANDELBROT B B. The Fractal Geometry of Nature [M]. New York：W. H. Freeman and Company，1982.
孫會來，金純，張文明，等. 工程車輛油氣懸架分數階建模與特性分析 [J]. 農業機械學報，2014，45(5)：16-21. Sun Huilai，Jin Chun，Zhang Wenming. Fractional Modeling and Characteristic Analysis of Hydro-Pneumatic Suspension for Construction Vehicle [J]. Transactions of the Chinese Society for Agricultural Machinery，2014， 45(5)：16-21.(in Chinese)
朱呈祥，鄒云. 分數階控制研究綜述 [J]. 控制與決策， 2009，24(2)：161-169.
Zhu Chengxiang，Zou Yun. Summary of Research on Fractional-Order Control [J]. Control and Decision， 2009，24(2)：161-169. (in Chinese)
車輛工程論文發表期刊推薦《鐵道車輛》（月刊）創刊于1963年，是由四方車輛研究所主辦的鐵道車輛專業的綜合性技術期刊。本刊辦刊宗旨是：認真貫徹執行“百花齊放、百家爭鳴”的方針，力求迅速、準確、全面反映和交流我國廣大鐵路車輛工作者在各個方面的研究成果和工作經驗，努力為鐵路運輸和生產的發展服務，為提高我國鐵路車輛的技術水平、實現鐵路運輸現代化做出貢獻。主要讀者對象：我國鐵路和路內外從事鐵道車輛試驗研究、設計制造和運用檢修的領導干部、科技人員、技術工人、管理干部和有關院校師生。
SCISSCIAHCI
