車輛蛇行運動是軌道車輛獨有的現(xiàn)象,由于車輪踏面近似為錐形,輪對運行時會發(fā)生橫向擺動,同時伴隨有繞其質(zhì)心的垂直軸來回搖頭轉(zhuǎn)動,運動軌跡類似“蛇”爬行過程的波形。蛇行運動存在臨界失穩(wěn)點,正常情況下輪對系統(tǒng)會在蠕滑力作用下自動對中,但一旦越過臨界失穩(wěn)點蛇行運動就會失控,輪軌蠕滑力不足以使輪對對中,輪對發(fā)生蛇行運動,橫向運動加劇,車輛振動平穩(wěn)性和運行安全性將嚴重變差,甚至有發(fā)生脫軌和翻車的可能,可以說蛇行運動是決定車輛能否高速運行的關鍵因素。

車輛出現(xiàn)蛇行失穩(wěn)的臨界狀態(tài)對應的速度稱為蛇行失穩(wěn)臨界速度,為了獲取準確的臨界速度,學者們將穩(wěn)定性和分叉理論引入蛇行運動問題中,建立了精細化的動力學模型,提出了許多定性以及定量的研究方法。曾京[1]和呂可維[2]分別運用打靶法和延續(xù)算法對非線性車輛蛇行運動Hopf分叉現(xiàn)象進行了數(shù)值求解,對輪對蛇行運動中的極限環(huán)性質(zhì)進行準確判別,可精準求解出非穩(wěn)定的虛環(huán),解決了龍格-庫塔法求解不穩(wěn)定極限環(huán)的困難。丹麥True[3]教授借助延續(xù)算法軟件包PATH 對非線性鐵道車輛系統(tǒng)的穩(wěn)定性與分岔問題進行了深入討論,取得了較為典型的理論性成果。Polach[4]結合鐵道車輛蛇行運動理論與試驗研究,詳細對比分析了幾種標準常用的構架失穩(wěn)評判標準的差異,指出車輛穩(wěn)定性的計算結果除了依賴于車輛結構參數(shù)之外,與輪軌接觸關系和數(shù)值計算方法也有關聯(lián)。高學軍研究了客車不對稱分岔行為[5],創(chuàng)新性地提出“合成分岔圖”的概念和構造方法。

但是上述這些研究都是在確定域中進行的,均沒有考慮隨機激擾的影響,這僅僅是理想情況下的理論分析,判定車輛系統(tǒng)是否失穩(wěn)通常都是基于系統(tǒng)部件橫向振動的衰減特性,但實際車輛線路運行時,由于持續(xù)存在隨機激擾,車輛振動響應曲線不會出現(xiàn)明顯的收斂或發(fā)散特征,此時確定性框架下的理論失穩(wěn)判定方法不再適用,已經(jīng)不能滿足高速車輛的穩(wěn)定性分析需求。我國對于線路運行過程中的車輛失穩(wěn)判定方法:對軸箱彈簧正上方對應的構架位置橫向加速度進行10 Hz低通濾波,若該橫向加速度峰值出現(xiàn)連續(xù)6次達到或超過極限值0.8g時,則判定轉(zhuǎn)向架橫向蛇行失穩(wěn)。但這一方法并無理論根據(jù),且由于我國列車運行環(huán)境、運行速度與國外的巨大差異,因此并不具備廣泛適用性。盡管也有許多學者提出了新的判定方法[6-8],但都有各自的局限性。

大量研究發(fā)現(xiàn),隨機激擾和非線性會出現(xiàn)相互耦合關系,兩者的耦合效應激發(fā)了激擾在非線性系統(tǒng)中的作用,這種作用在系統(tǒng)響應演化過程中具有決定性且不可控?；谶@一前提,我們需要跳出穩(wěn)定域框架,在隨機域中尋求新的突破,聚焦于對自然規(guī)律的描述更為本質(zhì)和真實的隨機系統(tǒng),探討隨機激擾對車輛動力學動態(tài)行為的影響,從機理上分析車輛隨機振動響應,研究隨機振動響應的內(nèi)在非線性特性及變化規(guī)律,為軌道車輛系統(tǒng)安全性能設計提供理論指導價值,在高鐵全面飛速發(fā)展的大背景下,這也是大勢所趨。

目前關于車輛隨機穩(wěn)定性的研究較少,比較典型的研究為劉偉渭等[9-10]和張波等[11-13]的研究工作,但這些研究工作都是基于超臨界輪對系統(tǒng),對亞臨界系統(tǒng)的隨機穩(wěn)定性并未做研究。本文在此基礎上,進一步深入研究隨機參激下亞臨界和超臨界分叉2種輪對系統(tǒng)的隨機穩(wěn)定性及隨機分叉,分析隨機激擾對穩(wěn)定性的影響,對比2種類型輪對系統(tǒng)的隨機分叉的差異,探討車輛隨機系統(tǒng)的失穩(wěn)判定方法,研究隨機臨界速度的計算方法。此外,本文關注的重點在于隨機激擾的影響,關于其他車輛參數(shù)的影響,已有大量研究,本文不做討論。

1 懸掛輪對系統(tǒng)隨機動力學建模與求解

1.1 模型建立

軌道車輛結構龐大,自由度較多,其非線性問題非常復雜,如果再考慮隨機影響,問題的復雜度將劇增。在車輛系統(tǒng)研究過程中,通常根據(jù)研究需要對車輛模型進行適當簡化,對于重點關注項點精細建模,對于非重點關注項點可做適當簡化。為了簡化理論研究的難度,本文選取自由度較少的懸掛輪對橫向動力學模型來作為分析對象。

建模過程中可作如下基本簡化[14]:

(1)構架與輪對之間的耦合效應可忽略,構架的搖頭以及側(cè)滾運動對輪對的影響可忽略。

(2)由于輪對側(cè)滾與橫移的耦合效應,僅考慮輪對搖頭和橫移2個自由度。

雖然近些年開始實施一系列措施，以促進黨員干部在村級事務管理、發(fā)家致富、農(nóng)村環(huán)境治理方面起到好的作用，但在具體實施過程中也出現(xiàn)了部分黨員不但起不到模范作用反而起到了反作用。部分黨組織書記安于現(xiàn)狀，對村內(nèi)事務管理存在“不求有功但求無過”思想，嚴重制約了村集體發(fā)展。

(3)不考慮輪對柔性,不考慮車輪與鋼軌脫離情形。

(4)暫不考慮一系懸掛的非線性特性,僅考慮輪軌間非線性特性。

(5)車輪踏面外形采用經(jīng)典的錐形踏面。

輪對系統(tǒng)建模時采用Kalker經(jīng)典線性輪軌蠕滑理論,輪軌非線性接觸關系可近似處理為一個非線性五次多項式,對應的2個非線性項系數(shù)α1、α2的選取依據(jù)來源于文獻[14]。車輛系統(tǒng)的隨機激擾主要包括隨機外激和隨機參激,前者主要包括氣動隨機激擾以及輪軌隨機激擾,后者主要是由于自身結構的隨機特性引起的激擾。目前商用動力學仿真軟件很容易實現(xiàn)外部激擾對車輛穩(wěn)定性的影響分析,但是關于隨機參激激擾的研究則鳳毛麟角,商用軟件也無法考慮參激激擾的影響,因此本文側(cè)重于研究隨機參激對穩(wěn)定性的影響,且此方法仍然適用于隨機外激。

聯(lián)合監(jiān)管機構雖抽調(diào)了相關部門的執(zhí)法人員參與現(xiàn)場監(jiān)管，并明確了各單位的職責，但具體的責任落實體系尚未建立起來。因管理涉及多部門，造成內(nèi)部管理綜合協(xié)調(diào)難度較大；再加上部分單位受利益驅(qū)動或成本投入影響，在某些問題上雖能達成一致意見，但很難具體落實到位。對因履行監(jiān)管職責不到位、工作不完全落實，導致偷采現(xiàn)象嚴重、采砂秩序混亂、影響防洪和通航安全的責任人，沒有相關的預警和追究措施。

實際車輛中,隨機參激也是無處不在的,車輛實際參數(shù)值和理論設計值通常會有一定差別,這一差別往往是隨機不可控的,而車輛系統(tǒng)一系懸掛參數(shù)對運動穩(wěn)定性的影響較大,因此研究一系懸掛參數(shù)的隨機變化對車輛穩(wěn)定性的影響意義重大。

實際上由于材料及結構的制造誤差,橡膠元件的老化以及懸掛元件的頻變特性都會引起參數(shù)的隨機變化,由于這些因素具有隨機性,無法準確描述其特性。為了描述參數(shù)的隨機性,通常的做法是結合大批量相應部件試驗結果統(tǒng)計分析參數(shù)分布規(guī)律,這種分布規(guī)律通常是參數(shù)理論值與一個隨機過程的疊加效應,這一隨機過程模型就是本文所說的隨機參激。鐵道車輛輪軌隨機激勵屬于強隨機激擾,主要呈現(xiàn)非高斯分步,而車輛一系懸掛參數(shù)的隨機激擾相對較弱且接近于高斯分布,因此工程分析時可近似處理為高斯白噪聲隨機激擾。這種方法有兩大優(yōu)勢:(1)可直接應用成熟的高斯白噪聲數(shù)學理論;(2)高斯白噪聲數(shù)學模型有利于理論分析,簡化理論分析難度。因此,本文選取高斯白噪聲經(jīng)典隨機過程來模擬輪對系統(tǒng)的一系隨機參激,高斯白噪聲的最鮮明特點是功率譜密度為常數(shù)且功率是無窮大。軌道車輛輪對系統(tǒng)搖頭和橫移耦合的隨機幾何模型示意圖見圖1。

圖1 輪對幾何模型示意

相關函數(shù)為

式中:v為車輛運行速度;ξ(t)為高斯白噪聲隨機參激激擾,該隨機過程均值為0;2D為強度。

首先，人際情緒管理受到社會標準和期待等文化因素的影響，不同文化背景下的個體在人際情緒管理方面存在顯著的差異.相比個體主義文化，集體主義文化中的人們更多地關注他人的情緒，但是自己遇到情緒問題時，尋求外部支持的意愿卻很小[23].因此，未來需要進一步開展跨文化研究，以增進對人際情緒管理的理解.

由圖1關系,可建立輪對橫移和搖頭運動耦合方程。

(1)輪對橫移運動

2.食品流通中控制關鍵控制點的對策。接受訂單。在該環(huán)節(jié)中，需要配備專門人員，使其負責這項工作。該環(huán)節(jié)控制對策為：針對儲存商品，確定相關內(nèi)容，主要包括商品名稱、儲存環(huán)境以及數(shù)量等，同時還需要判定商品是否易燃易爆以及有毒有害等。

(2)輪對搖頭運動

軌道車輛及輪對系統(tǒng)中的Hopf分叉有2種常見類型:超臨界和亞臨界分叉,為了研究這兩類系統(tǒng)的隨機動力學行為差異,分別選取2組輪對系統(tǒng)參數(shù)作為代表,各參數(shù)取值及物理意義詳見表1。關于隨機項控制系數(shù)β1和β2的選取,實際應用時可以通過部件大批量抽樣試驗統(tǒng)計激擾的能量來確定控制參數(shù),本文的參數(shù)取值僅僅只是一個示例,并不精確。

表1 兩類輪對系統(tǒng)參數(shù)

參數(shù)參數(shù)取值 亞臨界 超臨界輪對質(zhì)量m/kg 1 612 1 627輪對搖頭轉(zhuǎn)動慣量Iw/(kg·m3) 881 830輪對自旋轉(zhuǎn)動慣量Iy/(kg·m3) 100 100等效錐度λe 0.25 0.13車輪半徑r0/m 0.46 0.46軸重W/t 11.22 11.22一系懸掛橫向跨距之半l/m 1.02 1.02車輪滾動圓橫向跨距之半a/m 0.746 5 0.746 5縱向蠕滑系數(shù)f11/MN 13.7 1.5橫向蠕滑系數(shù)f22/MN 11.1 1.5橫向自旋蠕滑系數(shù)f23/MN 0 0自旋蠕滑系數(shù)f33/MN 0 0一系縱向定位剛度kx/(MN·m-1) 4 0.98一系橫向定位剛度ky/(MN·m-1) 1 0.98一系縱向阻尼cx/(k N·m·s-1) 0 0一系橫向阻尼cy/(k N·m·s-1) 0 0輪軌非線性控制參數(shù)α1 -1.6×1011 -1.6×1011輪軌非線性控制參數(shù)α2 1.6×1015 1.6×1015隨機參激項控制參數(shù)β1 1×102 1 060隨機參激項控制參數(shù)β2 1×102 1 060

上述輪對系統(tǒng)的運動方程表示成矩陣形式,即

式中:

1.2 輪對系統(tǒng)的隨機平均求解

方程(4)是隨機微分方程,其求解比較復雜,許多學者針對具體問題提出了相應的求解方法[15]并取得了很好的效果,但都不具備很好的推廣性。相對來說,隨機平均法應用最為廣泛,隨機平均法實際上是隨機意義上的Krylov-Bogoliubov-Mitropolsky(以下簡稱KBM)漸近法,許多學者進行了大量關于隨機平均法的研究,最為典型的就是浙江大學朱位秋院士的成果,朱院士在Hamilton框架下建立了一套完備的擬Hamilton系統(tǒng)的隨機平均法理論體系[16],可快速求解非線性系統(tǒng)的隨機響應,進而分析隨機穩(wěn)定性與隨機分叉等問題。

（1）從本次問題的處理來看，調(diào)節(jié)變壓器的分接頭后必須要測量各分接頭的直流電阻，如發(fā)現(xiàn)三相電阻不平衡，可以將分接開關進行旋轉(zhuǎn)處理去除其表面的氧化膜。由于本次#6B高廠變吊罩大修時，分接開關操作桿進行了拆裝，所以這項工作尤其必要，可確保分接開關觸頭接觸良好。

相比于能量和位移幅值,位移和速度隨時間的變化速率明顯更快,也就是說能量和位移幅值屬于慢變過程,而速度和位移屬于快變過程。由于慢變過程體現(xiàn)了系統(tǒng)長時間歷程內(nèi)的性態(tài),在研究中應重點關注。隨機平均求解的基本原理就是將快變過程隨時間變化的特性在時間域進行平均,從而將系統(tǒng)簡化為關于慢變過程的隨機平均方程。隨機平均法主要有以下兩大優(yōu)勢:(1)可對系統(tǒng)降維簡化問題的復雜度;(2)系統(tǒng)降維簡化信息不失真,可保留原系統(tǒng)的主要非線性特性。本文通過隨機平均法將輪對橫向位移和速度等快變過程組成的多維空間域映射成關于系統(tǒng)能量等慢變過程的一維域。

如果說養(yǎng)老地產(chǎn)硬件的投入比較容易到位，那么投資企業(yè)長久經(jīng)營，提供養(yǎng)老服務有時可能有心無力。由于開發(fā)養(yǎng)老地產(chǎn)資金回收期較長，投資企業(yè)一般不會考慮長期持有經(jīng)營，這又勢必影響?zhàn)B老地產(chǎn)本身的含金量。在實際運行中，投資養(yǎng)老地產(chǎn)的企業(yè)要量力而行，可考慮開發(fā)和經(jīng)營分開的做法：

取廣義位移q1=y,q2=φ,p1=,p2=,將輪對運動方程轉(zhuǎn)化為如下形式

前文提及,唐代的隋堤柳吟詠基本都集中在中晚唐時期,安史之亂前,除開元時期的王泠然外,鮮有詩歌涉及相關話題。為何中晚唐詩人對隋堤柳吟詠如此熱衷?這就不能不提到通濟渠、邗溝運道與唐代交通和中晚唐政治的關系了。與唐代詩歌中諸多的柳意象相比,隋堤柳可以說是直接由運河交通所催生的文學意象。中晚唐時期通濟渠、邗溝運道地位的巨大提升及由此帶來的運河懷古熱潮,與隋堤柳吟詠的勃興有直接關系。

式中

該輪對系統(tǒng)中,由于陀螺力[17]的存在,導致輪對系統(tǒng)變?yōu)榉潜Ｊ叵到y(tǒng),鑒于陀螺項的反對稱性,可取系統(tǒng)的Hamilton函數(shù)為

式中:H為對應輪對系統(tǒng)的廣義能量。

二十世紀初的中國，社會矛盾重重，東方文化與西方文化交匯，傳統(tǒng)意識與現(xiàn)代觀念發(fā)生碰撞。在這樣特殊的社會背景下，木蘭除了從父母那里接受傳統(tǒng)儒道文化的浸染，還接受了現(xiàn)代教育，她從十幾歲時就開始思考“男人與女人的分別這件事”，也經(jīng)常和父親談論“新時代的女子”的話題，是一個“新式的女孩子”。“木蘭有幾種女人所沒有的本領：第一，她會吹口哨兒；第二，她會唱京戲；第三，她收集古董，而且能鑒賞?！盵1]60可見，木蘭稱得上是東西方文化交融下，集傳統(tǒng)美德與現(xiàn)代性于一身的奇女子。

我國也出臺了很大畜牧獸醫(yī)方面的法律，在實際生活中也似乎起到了一些作用，但收效甚微。尤其是我國經(jīng)濟社會發(fā)展到今天，讓很多人和機構鉆了法律的空子，給畜牧獸醫(yī)的發(fā)展造成了很大的困擾。

依據(jù)隨機平均法,式(6)可轉(zhuǎn)化為It?(伊藤)隨機微分方程

式中:od B(t)表示隨機積分,且

根據(jù)Hamilton系統(tǒng)理論分類[18],該輪對隨機系統(tǒng)屬于擬不可積Hamilton系統(tǒng),利用擬不可積Hamilton系統(tǒng)隨即平均法對位移和速度等快變過程進行隨機平均,可得關于慢變過程能量H(t)的It?隨機平均微分方程

式中:

其中,

且R 滿足

2 輪對系統(tǒng)的隨機穩(wěn)定性與隨機分叉

2.1 輪對系統(tǒng)Hopf分叉

當不考慮隨機激擾時,即輪對系統(tǒng)(4)中隨機項滿足R=0。本文根據(jù)表1中的參數(shù)采用打靶法[1]分別求解兩類輪對系統(tǒng)橫向位移的Hopf分叉及極限環(huán)。為研究系統(tǒng)長時間歷程內(nèi)的性態(tài),通常做法是關注慢變過程——輪對橫向振動位移幅值隨速度增大過程中的變化情況,畫出Hopf分叉圖,見圖2。

圖2 輪對系統(tǒng)Hopf分叉

輪對系統(tǒng)發(fā)生了亞臨界Hopf分叉見圖2(a),系統(tǒng)既出現(xiàn)了Hopf分叉(Hopf分叉點A)也出現(xiàn)了鞍結分岔(鞍結分岔點B)。Hopf分叉點A對應的速度90.1 m/s為線性臨界速度,鞍結分岔點B對應的速度80.7 m/s為非線性臨界速度。當速度小于非線性臨界速度時,輪對對任何擾動均具有抗干擾作用,其穩(wěn)態(tài)運動都能收斂于平衡位置從而保持穩(wěn)定;當速度位于非線性臨界速度和線性臨界速度之間時,當輪對的初始擾動幅值較小(虛線以下)時,輪對系統(tǒng)穩(wěn)態(tài)運動逐漸收斂于平衡位置而保持穩(wěn)定,但當出現(xiàn)大幅初始擾動(虛線以上)時則逐漸趨于實線從而失穩(wěn),此時系統(tǒng)小激擾穩(wěn)定而大激擾失穩(wěn);當速度高于線性臨界速度時,系統(tǒng)完全失穩(wěn),任何擾動都能導致穩(wěn)態(tài)運動發(fā)散。

輪對超臨界Hopf分叉見圖2(b),Hopf分叉點A對應的速度為84.3 m/s,為臨界速度。當速度低于84.3 m/s時,輪對系統(tǒng)穩(wěn)定,對任何擾動均具有抗干擾效應,穩(wěn)態(tài)運動均收斂于平衡位置;當速度高于84.3 m/s時,輪對系統(tǒng)失穩(wěn),任何擾動都能導致穩(wěn)態(tài)運動發(fā)散直至出現(xiàn)極限環(huán)。

2.2 輪對系統(tǒng)的隨機穩(wěn)定性

隨機穩(wěn)定性分析可采用奇異邊界方法[19],這種方法的缺點是只適用于一維,對于高維問題則無能為力,但隨機平均法的降維功能恰好彌補了這一缺陷,為奇異邊界方法的實施奠定了基礎。

對于式(11)的一維能量擴散解過程H(t),其概率為1漸近穩(wěn)定性與其平穩(wěn)概率密度的存在性反映系統(tǒng)長期行為的性態(tài),而系統(tǒng)長期性態(tài)可通過邊界性態(tài)間接獲取。接下來將重點分析能量擴散過程H(t)在左右2個邊界上的性態(tài)。

對于右邊界H=∞,有|m(∞)|=∞,根據(jù)一維擴散過程邊界分類[20]可知,右邊界邊界類型為無窮遠處第二類奇異邊界。對式(14)關于θ 求期望,得

當H →∞時,R的低次項可忽略,H 應對應R的高次項,則有,代入式(12)可得

下標R表示右邊界,此時,有k24)H ＜0,且βR=1,根據(jù)邊界分類可知,右邊界H=∞為排斥自然邊界。

對于左邊界H=0,當H=0時,有σ2(0)=0,m(0)=0,則左邊界類型屬于第一類奇異邊界(套點)。此時對應的漂移和擴散系數(shù)滿足如下條件

根據(jù)第一類奇異邊界理論[20],邊界擴散指數(shù)αL、漂移指數(shù)βL 以及特征標值cL分別為(下標L表示左邊界)

當特征標值cL＜1時,則有左邊界H=0是吸引自然邊界,由于右邊界不受特征標值影響,這樣解曲線H(t)會由于排斥效應遠離左邊界而逐漸趨向于左邊界,說明能量逐漸穩(wěn)態(tài)收斂于0,此時系統(tǒng)解曲線H(t)是全局穩(wěn)定的。

籍此次活動，天王表也再次踐行公益之心，與百大名品中心溫暖攜手，一同向蚌埠特殊教育中心進行愛心捐贈，再次彰顯了品牌的責任擔當。

相反,當特征標值cL＞1時,左邊界也是排斥自然邊界,由于排斥效應解曲線H(t)會同時遠離左右邊界,解曲線H(t)最終會收斂于某個穩(wěn)態(tài)能量處(見圖3)。此時系統(tǒng)解曲線H(t)是隨機不穩(wěn)定的。

[4]中國科技信息研究所.2016年中國英文科技期刊引證報告[R].貴陽：2016年第8屆中國英文科技期刊研討會，2016

圖3 輪對隨機系統(tǒng)穩(wěn)定性邊界分類示意

穩(wěn)態(tài)能量的具體位置對應能量穩(wěn)態(tài)概率密度函數(shù)的峰值,而能量穩(wěn)態(tài)概率密度函數(shù)需要通過求解系統(tǒng)的FPK(Fokker Planck Kolmogorov equation)方程獲得,這項內(nèi)容涉及隨機P 分叉的研究內(nèi)容,本文暫不作討論,后續(xù)研究中會有涉及。

當cL=1時,左邊界為嚴格自然,此時邊界狀態(tài)處于吸引和排斥臨界狀態(tài),對應系統(tǒng)隨機穩(wěn)定性臨界狀態(tài),可導出系統(tǒng)的隨機穩(wěn)定域邊界。由式(18)可得

將表1中相關參數(shù)代入方程(19),可得兩類輪對系統(tǒng)的左邊界特征標值cL 與激擾強度D、速度v的關系(見圖4),2種情形下規(guī)律基本一致,特征標值cL 與速度v、激擾強度D 正相關,速度v和激擾強度D越大,cL越大。當速度v和激擾強度D 逐漸增大使特征標值大于1時,系統(tǒng)跨越絕對穩(wěn)定區(qū)域邊界,出現(xiàn)概率意義失穩(wěn)。

車輛穩(wěn)定性研究的首要任務是將車輛運行速度控制在穩(wěn)定安全域內(nèi),接下來將分別針對兩類輪對系統(tǒng),確定隨機穩(wěn)定域及穩(wěn)定區(qū)域邊界,為車輛安全性設計提供理論指導。

兩類輪對系統(tǒng)的隨機穩(wěn)定性臨界曲線見圖5,臨界曲線正好是絕對穩(wěn)定區(qū)域與隨機不穩(wěn)定區(qū)域的分界線,曲線上方為隨機不穩(wěn)定域,曲線下方為絕對穩(wěn)定域。車輛運行速度須控制在絕對穩(wěn)定區(qū)域內(nèi),若速度跨越穩(wěn)定域邊界,輪對發(fā)生隨機分叉,可能出現(xiàn)失穩(wěn),造成安全隱患。在穩(wěn)定性分析時需要根據(jù)車輛參數(shù)確定系統(tǒng)穩(wěn)定域,隨機激擾的激擾強度D 可以通過試驗模擬,只要確定了最大可能的激擾強度Dmax,然后將車速控制在Dmax對應的絕對穩(wěn)定區(qū)域內(nèi)即可保證車輛的安全性。

圖4 輪對隨機系統(tǒng)左邊界特征標值

圖5 輪對系統(tǒng)隨機穩(wěn)定臨界曲線

2.3 輪對系統(tǒng)的隨機分叉

隨機分叉與一般的確定性分叉不同,隨機分叉研究的重點主要集中于隨機激擾引起的非線性突變或躍遷現(xiàn)象。一般來說,隨機分叉主要有兩類:一類是隨機動態(tài)分叉(Dynamical bifurcation),又叫D 分叉;另一類是隨機維象分叉(Phenomenological bifurcation),也叫P 分叉。一般來說,動態(tài)分叉先于維象分叉,為了求解全局隨機臨界速度,本文只研究前者。

一直以來，我國的體育教育可能都不太被人們所重視，新課改中就要求中職院校在進行體育教育中必須要把學生的身心健康放在第一位，讓學生在學習體育的過程中感受到快樂，明白鍛煉身體的重要性，從而主動的去鍛煉身體，重視起體育課來。但是由于各種原因，學生沒有鍛煉身體，好好上體育課的意識，在一些中職院校中，更是沒有按規(guī)定完成體育課的教學任務。另一方面，教師也要在體育課程中重視學生的發(fā)展。畢竟體育教學對于學生的身體素質(zhì)的提高十分的重要，學院應該高度重視體育課的進度，密切關注學生的身體健康。

隨機動態(tài)分叉主要關注系統(tǒng)不變測度的動態(tài)特性,比較普遍的做法是通過研究Lyapunov指數(shù)這一不變測度的動態(tài)突變或躍遷變化,實現(xiàn)輪對的隨機動態(tài)分叉分析。因此,研究動態(tài)分叉的關鍵在于最大Lyapunov指數(shù)的計算。

實際上,由于隨機平均法保留原輪對系統(tǒng)的主要非線性特性,因此隨機平均后的平均It?隨機微分方程的最大Lyapunov指數(shù)可近似等價于原系統(tǒng)方程,則問題轉(zhuǎn)化為研究平均It?隨機微分方程(11)的最大Lyapunov指數(shù)。將平均It? 隨機微分方程(11)在H=0處線性化,得

由上式可求解得

則系統(tǒng)的最大Lyapunov指數(shù)為

當H →0時,式(14)中關于R的高次項可忽略不計,對該式R的低次項部分關于θ 求期望,可得

即R2=2H,則有

則可得輪對系統(tǒng)最大Lyapunov指數(shù)為

當λ＜0,輪對系統(tǒng)概率1漸近穩(wěn)定,此時系統(tǒng)的能量增長速率慢于能量消耗速率,系統(tǒng)的總能量會逐漸耗散并趨于0,系統(tǒng)絕對穩(wěn)定。最大Lyapunov指數(shù)λ的大小對應系統(tǒng)總能量的耗散速度,λ越小說明能量減小的越快,系統(tǒng)穩(wěn)定性能越好。

當λ＞0,輪對系統(tǒng)為概率1 漸近不穩(wěn)定,此時系統(tǒng)的能量消耗速率小于能量增長,總能量將逐漸累積,系統(tǒng)動能會逐漸增加,可能出現(xiàn)失穩(wěn)及分叉現(xiàn)象甚至發(fā)生脫軌。

當λ=0,輪對系統(tǒng)處于隨機穩(wěn)定臨界狀態(tài),λ 對應的零點為系統(tǒng)隨機動態(tài)分叉點。

輪對隨機系統(tǒng)的最大Lyapunov指數(shù)λ 這一特性有點類似于確定系統(tǒng)中的特征根實部,λ＜0是系統(tǒng)概率1漸近穩(wěn)定的充要條件。根據(jù)表1中的輪對參數(shù),分別畫出了不同激擾強度D 下兩類輪對系統(tǒng)的最大Lyapunov指數(shù)隨速度變化圖,見圖6。

由圖6可知,最大Lyapunov指數(shù)λ 與激擾強度D、運行速度v 具有正效應,當激擾強度D和運行速度v 增大時,最大Lyapunov指數(shù)會逐漸增大并由負變正,此時輪對系統(tǒng)由概率為1穩(wěn)定變?yōu)椴环€(wěn)定,系統(tǒng)發(fā)生隨機動態(tài)分叉。

圖6 最大Lyapounov指數(shù)隨速度變化

當λ=0時,為隨機動態(tài)分叉臨界點,此時對應的速度為動態(tài)分叉臨界速度,記為vD。

此外,分別針對兩類輪對系統(tǒng),對比分析奇異邊界法與最大Lyapounov指數(shù)法的分析結果,兩類輪對系統(tǒng)對比分析見圖7(圖中激擾強度D=0.05),對比結果顯示2種方法的分析結果完美契合,動態(tài)分叉臨界速度vD 與隨機穩(wěn)定域邊界對應的速度值正好相同,且兩類輪對系統(tǒng)結論一致,側(cè)面應證了奇異邊界法和最大Lyapunov指數(shù)法的有效性。

綜上所述,隨機參激對輪對系統(tǒng)臨界速度影響較大,當隨機激擾強度逐漸增大時,臨界速度顯著減小。隨機動態(tài)分叉前輪對系統(tǒng)動態(tài)特征為概率1意義下的隨機穩(wěn)定運動,此時系統(tǒng)絕對穩(wěn)定;隨機動態(tài)分叉后系統(tǒng)動態(tài)行為比較復雜,其具體特征需要通過P分叉進一步研究。

2.4 隨機動態(tài)分叉與Hopf分叉的對比

前文通過隨機動態(tài)分叉研究得到了隨機動態(tài)分叉臨界速度vD,但是隨機動態(tài)分叉概念不夠直觀,為了詳細說明隨機動態(tài)分叉以便于理解,接下來將隨機動態(tài)分叉與確定性Hopf分叉進行對比,探究隨機動態(tài)分叉和傳統(tǒng)分叉之間的聯(lián)系,建立隨機動態(tài)分叉和傳統(tǒng)Hopf分叉溝通的橋梁。

圖7 隨機穩(wěn)定性與動態(tài)分叉對比

隨機動態(tài)分叉臨界速度隨激擾強度的變化曲線見圖8。圖8中vDs和vDp分別為超臨界系統(tǒng)和亞臨界隨機動態(tài)分叉臨界速度,vH為超臨界系統(tǒng)Hopf分叉臨界速度,vA為亞臨界輪對系統(tǒng)鞍結分叉臨界速度。從圖8可以看出,當D→0時,超臨界輪對系統(tǒng)中動態(tài)分叉點逐漸趨近于Hopf分叉點,而亞臨界輪對系統(tǒng)中動態(tài)分叉點逐漸趨近于鞍結分叉點。當激擾強度時,隨機動態(tài)分叉逐漸退化為確定性分叉。需要說明的是超臨界輪對系統(tǒng)Hopf分叉點和亞臨界輪對系統(tǒng)鞍結分叉點都是各自輪對系統(tǒng)的全局(最小)臨界速度,當速度小于全局臨界速度時,系統(tǒng)絕對穩(wěn)定,也就是說隨機動態(tài)分叉臨界速度可定義為輪對系統(tǒng)的全局隨機臨界速度。在隨機激擾下,臨界速度出現(xiàn)向前漂移現(xiàn)象,且隨著激擾強度增大,動態(tài)分叉臨界速度顯著變小。

圖8 隨機動態(tài)分叉臨界速度隨激擾強度變化曲線

在傳統(tǒng)車輛穩(wěn)定性理論中,認為車輛臨界速度跟外部激擾無關,激擾只是改變了系統(tǒng)的初始狀態(tài),一旦激擾撤去,系統(tǒng)穩(wěn)態(tài)解仍會呈現(xiàn)Hopf分叉相應特征。這一觀點本身并沒有問題,只是考慮問題的角度不同而已,實際上外部激擾存在時,車輛系統(tǒng)不能再單純考慮車輛本身,應該將它們作為一個耦合的大系統(tǒng)來分析。這也打破了確定性車輛系統(tǒng)中的臨界速度的思維定勢,車輛的臨界速度不再是單一的固定值,而是隨激擾強度變化而改變的,這也可以從理論上解釋車輛在不同線路下的失穩(wěn)臨界速度不同的原因。

在仔豬白痢發(fā)病的后期，仔豬會出現(xiàn)不良狀況，上述方法可用于治療仔豬。此外，還可在靜脈注射或腹腔注射的協(xié)助下，利用添加一定量抗生素的葡萄糖鹽水治療患病仔豬，為防止機體酸中毒的情況發(fā)生，可以利用碳酸氫鈉（3%～5%）；如果患病仔豬心臟出現(xiàn)衰竭，則需要配合利用強心藥物進行及時治療。

車輛在實際運行中隨機激擾無處不在,常規(guī)情況下車輛受到的隨機激擾源主要是軌道不平順激擾,但相對來說軌道不平順激擾屬于弱隨機因素,一般情況下其激擾強度D 較小,根據(jù)前述理論分析結果可知,系統(tǒng)的臨界速度漂移量較小,符合實際車輛運行情況。但是當遇到大風、強震等極端激擾條件下,此時激擾強度D 較大,車輛系統(tǒng)的臨界速度會顯著降低,穩(wěn)定性能急劇惡化,而用Hopf分叉理論顯然無法解釋這一點。因此,常規(guī)運行時分析車輛穩(wěn)定性時,隨機激擾影響不大,但遇到大風等強激擾環(huán)境時,隨機激擾的影響將不可忽視。

此外,隨機動態(tài)分叉點對應的速度為系統(tǒng)的全局隨機臨界速度,最大Lyapunov指數(shù)法對隨機激擾具有魯棒性,這為新的線路失穩(wěn)判定方法提供了理論依據(jù)。最大Lyapunov指數(shù)可以通過時間響應序列直接計算,憑借這一優(yōu)勢,該方法將可以直接應用到線路失穩(wěn)判定中。

3 結論

本文運用隨機平均法研究了兩類輪對系統(tǒng)的隨機穩(wěn)定性以及隨機動態(tài)分叉,定義了全局隨機臨界速度,提出了以最大Lyapunov指數(shù)法為基礎的失穩(wěn)判定方法,獲得了以下主要結論:

(1)隨機動態(tài)分叉點對應的速度為輪對系統(tǒng)的全局隨機臨界速度,超臨界Hopf分叉系統(tǒng)中動態(tài)分叉點對應于Hopf分叉點,亞臨界Hopf分叉系統(tǒng)中動態(tài)分叉點對應于鞍結分叉點。

(2)動態(tài)分叉前輪對系統(tǒng)動態(tài)特征為概率1 意義下的穩(wěn)定運動,分叉后輪對系統(tǒng)具體特征需要通過P分叉進一步研究。

(3)隨機因素對輪對系統(tǒng)平衡點的穩(wěn)定性有決定性影響,隨機激擾使輪對系統(tǒng)的臨界速度出現(xiàn)向前漂移,隨著激擾強度增大,臨界速度漂移量顯著增大。

(4)輪對系統(tǒng)的臨界速度不再是單一的固定值,而是隨激擾強度變化而改變的,這也從理論上解釋了不同線路條件下失穩(wěn)臨界速度不相同的原因。

(5)常規(guī)運行時分析車輛穩(wěn)定性時,隨機激擾影響不大,但遇到大風等強激擾環(huán)境時,隨機激擾的影響將不可忽視。

(6)最大Lyapunov指數(shù)法可作為線路失穩(wěn)判定的新方法。

參考文獻：

[1]曾京.車輛系統(tǒng)的蛇行運動分叉及極限環(huán)的數(shù)值計算[J].鐵道學報,1996,18(3):13-19.

ZENG Jing.Numerical Computations of the Hunting Bifurcation and Limit Cycles for Railway Vehicle System[J].Journal of the China Railway Society,1996,18(3):13-19.

[2]呂可維.車輛系統(tǒng)非線性動力學問題研究[D].成都:西南交通大學,2004:87-91.

[3]TRUE H.On the Theory of Nonlinear Dynamics and Its Applications in Vehicle Systems Dynamics[J].Vehicle System Dynamics,1999,31(5/6):393-421.

[4]POLACH O.On Non-linear Methods of Bogie Stability Assessment Using Computer Simulations[J].Proceedings of the Institution of Mechanical Engineers,Part F:Journal of Rail and Rapid Transit,2006,220(1):13-27.

[5]高學軍.鐵道客車系統(tǒng)橫向運動對稱不對稱分岔行為與混沌研究[D].成都:西南交通大學,2010:25-37.

[6]賈璐,曾京,池茂儒.車輛系統(tǒng)橫向運動穩(wěn)定性評判的數(shù)值仿真研究[J].鐵道車輛,2011,49(9):1-6.

JIA Lu,ZENG Jing,CHI Maoru.The Numerical Simulation Research on Judgment of Lateral Motion Stability of the Vehicle[J].Rolling Stock,2011,49(9):1-6.

[7]朱慶龍.轉(zhuǎn)向架橫向穩(wěn)定性實時監(jiān)測評估方法研究[D].成都:西南交通大學,2013:60-85.

[8]孫麗霞.高速列車橫向運動穩(wěn)定性和動態(tài)脫軌理論分析及評價方法研究[D].北京:中國鐵道科學研究院,2014:111-125.

[9]劉偉渭,戴煥云,曾京.軌道隨機激勵對彈性約束輪對隨機穩(wěn)定性與隨機分岔的影響[J].振動與沖擊,2013,32(21):170-177.

LIU Weiwei,DAI Huanyun,ZENG Jing.Influence of Rail Stochastic Excitation on Stochastic Stability and Stochastic Bifurcation of an Elastic Constraint Wheelset System[J].Journal of Vibration and Shock,2013,32(21):170-177.

[10]劉偉渭,戴煥云,劉轉(zhuǎn)華,等.彈性約束輪對系統(tǒng)的隨機Hopf分岔研究[J].鐵道學報,2013,35(10):38-45.

LIU Weiwei,DAI Huanyun,LIU Zhuanhua,et al.Research on Stochastic Hopf Bifurcation of Elastic Constraint Wheelset System[J].Journal of the China Railway Society,2013,35(10):38-45.

[11]張波,曾京,劉偉渭.Gauss白噪聲參激下懸掛輪對系統(tǒng)的隨機穩(wěn)定性研究[J].振動與沖擊,2015,34(19):49-56.

ZHANG Bo,ZENG Jing,LIU Weiwei.Research on Stochastic Stability of Suspension Wheelset System under Gauss White Noise[J].Journal of Vibration and Shock,2015,34(19):49-56.

[12]ZHANG B,ZENG J,LIU W W.Research on Stochastic Stability and Stochastic Bifurcation of Suspended Wheelset[J].Journal of Mechanical Science and Technology,2015,29(8):3097-3107.

[13]張波.鐵道車輛系統(tǒng)隨機穩(wěn)定性及隨機分叉研究[D].成都:西南交通大學,2016:26-53.

[14]WAGNERU V.Nonlinear Dynamic Behaviour of a Railway Wheelset[J].Vehicle System Dynamics,2009,47(5):627-640.

[15]劉先斌,陳虬,陳大鵬.非線性隨機動力系統(tǒng)的穩(wěn)定性和分岔研究[J].力學進展,1996,26(4),437-452.

LIU Xianbin,CHEN Qiu,CHEN Dapeng.The Researchs on the Stability and Bifurcation of Nonlinear Stochastic Dynamical Systems[J].Advances in Mechanics,1996,26(4),437-452.

[16]朱位秋.非線性隨機動力學與控制——Hamilton 理論體系框架[M].北京:科學出版社,2003.

[17]張波,曾京,董浩.非線性輪對陀螺系統(tǒng)的穩(wěn)定性及分叉研究[J].振動.測試與診斷,2015,35(5):955-960.

ZHANG Bo,ZENG Jing,DONG Hao.Research on Stability and Bifurcation of Nonlinear Gyroscopic Wheelset System[J].Journal of Vibration,Measurement& Diagnosis,2015,35(5):955-960.

[18]ZHU W Q,YANG Y Q.Stochastic Averaging of Quasi-Nonintegrable-Hamiltonian Systems[J].Journal of Applied Mechanics,1997,64(1):157-164.

[19]劉先斌.隨機力學系統(tǒng)的分叉行為與變分方法研究[D].成都:西南交通大學,1995:52-57.

[20]葛根.矩形薄板振動的隨機分岔和可靠性研究[D].天津:天津大學,2009:18-23.

[21]齊德芬,田也壯,曹澤輝.隨機擾動下激勵系統(tǒng)的分岔與控制研究[C]//中國管理科學學術年會論文集.北京:中國優(yōu)選法統(tǒng)籌法與經(jīng)濟數(shù)學研究會,2008:23-28.