一種基于動態剛度系數和阻尼系數的轉子軸心軌跡求解方法，該方法包括下列步驟：根據滑動軸承的所受到的動載荷形式，并將固-液界面的速度滑移通過Navier速度滑移模型進行模擬，建立不定常狀態下的油膜潤滑雷諾方程，采用有限差分法并結合載荷增量法和擾動壓力法對不定常工況雷諾進行求解，得到軸承的四個動剛度系數和四個阻尼系數，將軸承油膜力通過四個動剛度系數和四個阻尼系數表示；分析轉子所受到的動載荷形式，建立轉子的運動平衡方程，并引入油膜力的油膜動態系數的表達形式，而運動平衡方程的解即為轉子軸心軌跡坐標，從而實現了轉子動態特征系數和軸心軌跡的聯立求解，同時分析軸承-轉子系統的動態性能及系統的穩定性。
【技術實現步驟摘要】
一種基于動態剛度系數和阻尼系數的轉子軸心軌跡求解方法
本專利技術涉及一種仿真方法，適用于求解滑動軸承的動態剛度系數、阻尼系數和轉子軸心軌跡，實現轉子系統的動態特性及穩定性的同步分析。
技術介紹
滑動軸承在工作過程中運行平穩、可靠并且無噪聲。在滑動軸承內，軸瓦表面和軸頸表面被潤滑油分開而不發生直接接觸，這大大減小了轉子運行過程中的摩擦損失和表面磨損，而且油膜還具有一定的吸振能力。實際情況下油膜的厚度很小，一般為um數量級，所以油液在軸承中的流動屬于微尺度的研究領域，使用微尺度流體流動的研究方法研究油膜的流動規律已成為一種必然的趨勢。目前的研究表明，微尺度流體流動在機理上迥異于常規尺度流動，微尺度流動所基于的物理因素與宏觀流動不同，一方面是尺寸的縮小引起的尺寸效應，表現為表面積與體積比增大，粘滯力、界面粗糙度、梯度參數效應對流體的流動影響增強，另一方面為一些對宏觀流體流動可以忽略的影響因素在微尺度流動中影響作用逐漸增強，某些因素將不能忽略。微尺度流體流動與宏觀尺度的差異性在實驗研究中也得到了證實。這些物理因素作用程度的改變或者新的物理因素的參與使得流體力學理論中流體流動的基本方程和邊界條件需要進行一定的修正，建立能夠表述微尺度流體流動特性的數學模型。在滑動軸承的研究中表現為對潤滑描述方程N-S方程進行修正。實際中油膜起著非線性彈簧和阻尼的作用；而在一些情況下也可以將油膜化簡為線性剛度的彈簧和阻尼，用線性化的油膜剛度和阻尼來衡量軸承-轉子系統內油膜的動力學特性，通過四個動態剛度系數和四個阻尼系數來衡量，這八個參數對轉子的動力學計算和系統穩定非常重要。軸心軌跡是一種評判轉子系統穩定性的方法，用一種更加直觀的方式反映了轉子瞬時運動情況及運動穩定性。轉子在工作的過程中，其運動會受到各種不平衡擾動的作用，比如轉子偏心質量造成不平衡動載荷、轉子部件的缺損造成給轉子帶來的不平衡動載荷，以及切削過程中切削力也會成為轉子不平衡動載荷的來源。轉子在單一不平衡動載荷或者復合動載荷的擾動下會產生不同軸心形狀的運動軌跡，對軸心軌跡進行分析能夠實現轉子系統在動載荷擾動下的動態性能、穩定狀態的評估。
技術實現思路
針對滑動軸承動態性能求解中問題，本專利技術提供一種基于動態剛度系數和阻尼系數的轉子軸心軌跡求解方法，本專利技術避免了在求解非線性油膜力時的大計算量，同時充分利用軸承的動態特征系數，將動態特征系數應用與轉子軸心軌跡的求解，實現軸承動態性能及系統穩定性的同步分析。一種基于動態剛度系數和阻尼系數的轉子軸心軌跡求解方法，該方法包括以下步驟：(1)建立考慮軸承內固-液界面速度滑移情況下的滑動軸承油膜潤滑雷諾方程，該方程沒有將轉子受的動載荷忽略，如下式式中，u為旋轉主軸表面的線速度；b為固-液界面的滑移系數，p為油膜壓力；η為潤滑油的動力粘度；h為油膜厚度；x為油膜流動周向方向坐標；z為油膜流動的軸向方向坐標；其中油膜厚度對時間的導數可以表示為式中，為位置角，記為從轉子軸心和軸承軸心的連心線為起點，按逆時針方法旋轉；(2)對雷諾方程(1)和油膜厚度隨時間的變化方程(2)采取無量綱化處理，取油膜特征壓強為p0，特征滑移長度和油膜特征厚度為軸承半徑間隙h0，油膜特征軸向長度為軸承長度的一半L/2，油膜特征周向長度為軸承半徑R，令p＝p0P，其中b＝b0h0，h＝Hh0，τ＝w0t，w0為轉子旋轉的角速度，不定常工況下的無量綱雷諾方程為：無量綱油膜厚度隨時間變化為：(3)采用有限差分方法對方程(3)線性化處理，得到如下線性化方程：其中，F＝D+E式中，R為軸承半徑，L為軸承長度，Δλ為軸承軸向方向網格長度，為軸承周向方向網格長度，(i，j)為油膜位置坐標，Hi,j表示(i，j)處無量綱油膜厚度；應用超松弛迭代法將線性化方程(5)在matlab軟件平臺上進行數值求解，即得擾動載荷下油膜壓力分布；(4)載荷增量法計算油膜動剛度。載荷增量法的計算思想就是將轉子軸心在平衡位置的基礎上沿不同的方向產生微小位移，求出在該移動后位置處的油膜力，應用油膜力的增量值和微小位移的比值用于計算油膜動剛度。1)動剛度系數Kxx和Kyx計算如圖1所示，O為軸承的中心位置，O'為轉子的平衡位置，轉子沿x方向分別產生微小的位移擾動Δx后的位置為O1和O2。轉子在平衡位置O'沿x方向向左產生微小擾動后，軸心位置移動到O1處，此時，x方向和y方向的油膜力分別為Fx1和Fy1；轉子在平衡位置O'沿x方向向右產生微小擾動后，軸心位置移動到O2處，此時，x方向和y方向的油膜力分別為Fx2和Fy2；根據油膜剛度的計算公式，油膜動剛度可以表示為：2)動剛度系數Kxy和Kyy計算如圖2，O為軸承的中心位置，O'為轉子的平衡位置，轉子沿y方向分別產生微小的位移擾動Δy后的位置為O3和O4。轉子在平衡位置O'沿y方向向下產生微小擾動后，軸心位置移動到O3處，此時，x方向和y方向的油膜力分別為Fx3和Fy3；轉子在平衡位置O'沿y方向向上產生微小擾動后，軸心位置移動到O4處，此時，x方向和y方向的油膜力分別為Fx4和Fy4；根據油膜剛度的計算公式，油膜動剛度可以表示為：(5)擾動壓力法計算油膜阻尼。如圖3所示，轉子在外載荷的作用下平衡于O'位置，在(Δe,Δθ)的擾動下，軸心位置移動到O1處，此時油膜力沿Δe的方向和垂直于Δe的方向分別為Fe和Fθ，油膜沿Δe的方向和垂直于Δe的方向的阻尼可以定義為：將方程(10)和(11)分別對和求導，并結合公式(12)、(13)、(14)和(15)可得：其中和是油膜壓力p對擾動量和求偏導數。將沿Δe的方向和垂直于Δe的方向的阻尼經過坐標變換，變換到(x,y)方向：根據無量綱雷諾方程(3)得擾動壓力和的微分形式為：式中：算子Rey表示(6)轉子運動方程建立。如圖4，O為軸承中心位置，O'為轉子中心位置，轉子在x，y方向所受到的動載荷分別為Qx和Qy，t時刻對應的油膜力在x，y方程的油膜力分別為轉子以w0的角速度做回轉運動，此時油膜力在x，y方向的分量分別為Fx(w0t)和Fy(w0t)，Mg為軸承承受的的轉子重量，則轉子軸心的運動方程為：(7)不平衡動載荷轉子軸心軌跡方程建立。轉子由于制造誤差、裝配誤差及材料缺陷原因難免會使得轉子的軸心位置和轉子的質心不重合，使得轉子存在偏心質量，偏心質量是轉子不平衡動載荷的主要來源之一。偏心質量會導致轉子軸心沿著一定的軌跡運動，通過分析轉子軸心軌跡能夠分析轉子的動力學行為；若轉子受到單一的偏心質量動載荷影響，eg為轉子的質量偏心距，所以作用在轉子上的不平衡載荷為：Qx＝Megw02sin(w0t)(25)Qy＝Megw02cos(w0t)(26)將公式(23)和(24)中的油膜力Fx(w0t)和Fy(w0t)用油膜的動剛度和阻尼系數，所以轉子的運動方程可以寫成：無量綱形式為：式中：當量質量無量綱偏心質量距轉子當量重量(8)轉子軸心軌跡坐標求解。無量綱軸心軌跡坐標(X,Y)就是方程式(29)和(30)的解，解的形式為X＝λ1cosτ+λ2sinτ(31)Y＝λ3cosτ+λ4sinτ(32)其中λ1，λ2，λ3，λ4可以通過下式求得：通過式(33)解出的四個值帶入方程(31)和(32)中，就可以得到轉子軸心運動的坐標(X,Y)本文檔來自技高網...
【技術保護點】
一種基于動態剛度系數和阻尼系數的轉子軸心軌跡求解方法，其特征在于：該方法包括以下步驟，(1)建立考慮軸承內固?液界面速度滑移情況下的滑動軸承油膜潤滑雷諾方程，該方程沒有將轉子受的動載荷忽略，如下式∂∂x={[-h312η-h2b2η]∂p∂x+hu2}+∂∂z[(-h312η-h2b2η)∂p2η]=-∂h∂t]]>式中，u為旋轉主軸表面的線速度；b為固?液界面的滑移系數，p為油膜壓力；η為潤滑油的動力粘度；h為油膜厚度；x為油膜流動周向方向坐標；z為油膜流動的軸向方向坐標；其中油膜厚度對時間的導數可以表示為式中，為位置角，記為從轉子軸心和軸承軸心的連心線為起點，按逆時針方法旋轉；(2)對雷諾方程(1)和油膜厚度隨時間的變化方程(2)采取無量綱化處理，取油膜特征壓強為p0，特征滑移長度和油膜特征厚度為軸承半徑間隙h0，油膜特征軸向長度為軸承長度的一半L/2，油膜特征周向長度為軸承半徑R，令p＝p0P，其中b＝b0h0，h＝Hh0，z=λL2(-1≤λ≤1),]]>(X,Y)=(x,y)h0,]]>(X.,Y.)=(x.,y.)h0w0,]]>τ＝w0t，w0為轉子旋轉的角速度，不定常工況下的無量綱雷諾方程為：無量綱油膜厚度隨時間變化為：(3)采用有限差分方法對方程(3)線性化處理，得到如下線性化方程：Pi,j=12Ai,j+2Ci,j[(Ai,j+Bi,j)Pi+1,j+(Ai,j-Bi,j)Pi-1,j+Ci,jPi,j+1+Ci,jPi,j-1-Fi,jPi,j]---(5)]]>其中，Ci,j=β2Hi,j3(1+6b0)Δλ2]]>F＝D+Eβ=2RL]]>式中，R為軸承半徑，L為軸承長度，Δλ為軸承軸向方向網格長度，為軸承周向方向網格長度，(i，j)為油膜位置坐標，Hi,j表示(i，j)處無量綱油膜厚度；應用超松弛迭代法將線性化方程(5)在matlab軟件平臺上進行數值求解，即得擾動載荷下油膜壓力分布；(4)載荷增量法計算油膜動剛度；載荷增量法的計算思想就是將轉子軸心在平衡位置的基礎上沿不同的方向產生微小位移，求出在該移動后位置處的油膜力，應用油膜力的增量值和微小位移的比值用于計算油膜動剛度；1)動剛度系數Kxx和Kyx計算O為軸承的中心位置，O'為轉子的平衡位置，轉子沿x方向分別產生微小的位移擾動Δx后的位置為O1和O2；轉子在平衡位置O'沿x方向向左產生微小擾動后，軸心位置移動到O1處，此時，x方向和y方向的油膜力分別為Fx1和Fy1；轉子在平衡位置O'沿x方向向右產生微小擾動后，軸心位置移動到O2處，此時，x方向和y方向的油膜力分別為Fx2和Fy2；根據油膜剛度的計算公式，油膜動剛度可以表示為：Kxx=Fx1-Fx22Δx---(6)]]>Kyx=Fy1-Fy22Δx---(7)]]>2)動剛度系數Kxy和Kyy計算O為軸承的中心位置，O'為轉子的平衡位置，轉子沿y方向分別產生微小的位移擾動Δy后的位置為O3和O4；轉子在平衡位置O'沿y方向向下產生微小擾動后，軸心位置移動到O3處，此時，x方向和y方向的油膜力分別為Fx3和Fy3；轉子在平衡位置O'沿y方向向上產生微小擾動后，軸心位置移動到O4處，此時，x方向和y方向的油膜力分別為Fx4和Fy4；根據油膜剛度的計算公式，油膜動剛度可以表示為：Kxy=Fx3-Fx42Δy---(8)]]>Kyy=Fy3-Fy42Δy---(9)]]>(5)擾動壓力法計算油膜阻尼；轉子在外載荷的作用下平衡于O'位置，在(Δe,Δθ)的擾動下，軸心位置移動到O1處，此時油膜力沿Δe的方向和垂直于Δe的方向分別為Fe和Fθ，油膜沿Δe的方向和垂直于Δe的方向的阻尼可以定義為：cee=(∂Fe∂e.)0---(12)]]>ceθ=(∂Fee∂θ.)0---(13)]]>cθe=(∂Fθ∂e.)0---(14)]]>cθθ=(∂Fθe∂θ.)0---(15)]]>將方程(10)和(11)分別對和求導，并結合公式(12)、(13)、(14)和(15)可得：其中和是油膜壓力p對擾動量和求偏導數；將沿Δe的方向和垂直于Δe的方向...
【技術特征摘要】
1.一種基于動態剛度系數和阻尼系數的轉子軸心軌跡求解方法，其特征在于：該方法包括以下步驟，(1)建立考慮軸承內固-液界面速度滑移情況下的滑動軸承油膜潤滑雷諾方程，該方程沒有將轉子受的動載荷忽略，如下式式中，u為旋轉主軸表面的線速度；b為固-液界面的滑移系數，p為油膜壓力；η為潤滑油的動力粘度；h為油膜厚度；x為油膜流動周向方向坐標；z為油膜流動的軸向方向坐標；其中油膜厚度對時間的導數可以表示為式中，為位置角，記為從轉子軸心和軸承軸心的連心線為起點，按逆時針方法旋轉；(2)對雷諾方程(1)和油膜厚度隨時間的變化方程(2)采取無量綱化處理，取油膜特征壓強為p0，將軸承半徑間隙h0參數作為油膜的特征滑移長度和特征厚度，油膜特征軸向長度為軸承長度的一半L/2，油膜特征周向長度為軸承半徑R，令p＝p0P，其中b＝b0h0，h＝Hh0，τ＝w0t，w0為轉子旋轉的角速度，不定常工況下的無量綱雷諾方程為：無量綱油膜厚度隨時間變化為：(3)采用有限差分方法對方程(3)線性化處理，得到如下線性化方程：其中，F＝D+E式中，R為軸承半徑，L為軸承長度，Δλ為軸承軸向方向網格長度，為軸承周向方向網格長度，(i，j)為油膜位置坐標，Hi,j表示(i，j)處無量綱油膜厚度；應用超松弛迭代法將線性化方程(5)在matlab軟件平臺上進行數值求解，即得擾動載荷下油膜壓力分布；(4)載荷增量法計算油膜動剛度；載荷增量法的計算思想就是將轉子軸心在平衡位置的基礎上沿不同的方向產生微小位移，求出在該移動后位置處的油膜力，應用油膜力的增量值和微小位移的比值用于計算油膜動剛度；1)計算動剛度系數Kxx和KyxO為軸承的中心位置，O'為轉子的平衡位置，轉子沿x方向分別產生微小的位移擾動Δx后的位置為O1和O2；轉子在平衡位置O'沿x方向向左產生微小擾動后，軸心位置移動到O1處，此時，x方向和y方向的油膜力分別為Fx1和Fy1；轉子在平衡位置O'沿x方向向右產生微小擾動后，軸心位置移動到O2處，此時，x方向和y方向的油膜力分別為Fx2和Fy2；根據油膜剛度的計算公式，油膜動剛度可以表示為：2)計算動剛度系數Kxy和KyyO為軸承的中心位置，O'為轉子的平衡位置，轉子沿y方向分別產生微小的位移擾動Δy后的位置為O3和O4；轉子在平衡位置O'沿y方向向下產生微小擾動后，軸心位置移動到O3處，此時，x方向和y方向的油膜力分別為Fx3和Fy3；轉子在平衡位置O'沿y方向向上產...
【專利技術屬性】
技術研發人員：陳東菊，邊艷華，周帥，范晉偉，
申請(專利權)人：北京工業大學，
類型：發明
國別省市：北京;11
還沒有人留言評論。發表了對其他瀏覽者有用的留言會獲得科技券。