曲柄滑塊機構基于連桿-滑塊間隙流體潤滑的動態(tài)分析雷戈瑞.丹尼爾，凱迪.凱文機械設計部，坎皮納斯大學，機械工程學院，郵件信箱 6122，13083-970 坎皮納斯，SP ，巴西摘要傳統(tǒng)的曲柄滑塊連桿機構廣泛的應用于機械系統(tǒng)中。軸承在機械中的使用特別作用在于減小摩擦，主要用于特殊條件的潤滑，如曲柄滑塊的鉸鏈連接。這類軸承屬于一類交替旋轉運動的軸承?？紤]到曲柄滑塊系統(tǒng)的動力學與軸承潤滑現象是相互作用的，本文為這種特殊問題提供了一種數學模型。兩種數學模型被使用與動態(tài)系統(tǒng)的分析。第一種模型（Eksergian 的運動方程），假設，連桿一端與軸承表面接觸時代表系統(tǒng)，在這種情況下，連桿的行為與剛性軸承相同（無間隙）。第二種模型（拉格朗日方法），當連桿一端在以在滑塊的間隙中為潤滑方式時，則代表系統(tǒng)，在這種情況下，在關系到連桿滑塊移動，則呈現出多自由度問題。流體動力潤滑的數學模型介紹了獲得系統(tǒng)動態(tài)行為中更為現實的結果。1.介紹在眾多功能中潤滑系統(tǒng)是評價機械最佳性能的一個至關重要的因素，例如，潤滑和保護部件，減少摩擦，清洗和冷卻的內部機制。本系統(tǒng)的操作需要適當的校準。因此，過多或者潤滑不足會影響機器的動態(tài)行為，會嚴重損害所設計的機構。許多機器的機構，包含了大量的流體動壓軸承，而軸承連桿-滑塊聯合尤為重要，因為這是一個新的流體動力軸承，被稱為流體動力軸承交替運動。不同于傳統(tǒng)的動壓軸承，這種軸承不是一個完整的旋轉。因此，隨著近來來日益增長的需要，研究這一特定類型的流體動力軸承和機器的動態(tài)行為的影響。目前為止，大多數關于滑塊連接的鉸鏈的流體動力軸承正在被研究所所研究著。事實上，該研究所已開發(fā)和建造的設備來研究摩擦潤滑滑塊銷軸承超過 15 年。比如，Takiguchi [ 1 ]等人研究一種浮式活塞在汽油動力汽車引擎上采用旋轉運動。三年后，他 [ 2 ]開發(fā)一個測量的設備，確定在連桿 –滑塊聯合的基礎上測得的摩擦力的流體動力軸承的潤滑狀態(tài)。在另一項研究中，在武藏技術研究所進行，有 Suhara 等人 [ 3 ]檢查潤滑條件下軸承的滑銷的汽油發(fā)動的汽車發(fā)動機，作為考慮的參數分析的長度，內徑和活塞銷的材料。最近，張等人 [ 4 ]開發(fā)的工具，探討在活塞銷軸承磨損。2005，利吉爾、拉戈分析了活塞銷軸承的行為。一年后，這些作者 [ 7 ]表達餓了連桿 –活塞接頭在四沖程發(fā)動機軸承運行的總的看法，強調在軸承供油。如前所述，流體動壓軸承的往復運動的滑塊曲柄機構一部分。因此，這一機制可靠的的數學模型必須考慮連桿 –滑塊聯合軸承的行為。由于曲柄滑塊機構的廣泛應用，許多研究都集中在設計的數學模型，分析了機械系統(tǒng)的動力學行為。施瓦布，梅嘉德，梅耶爾 [ 8 ]比較了滑塊曲柄連桿機構的動力學行為，考慮到不同的連桿 –滑塊聯合模型如耗散赫茲接觸模型，影響模型和動壓軸承模型在上述工作，在考慮連桿彈性與剛性構件的同時，作者還比較了動態(tài)行為的機制。那些作者所獲得的結果表明，連桿的假設為彈性元件，以及在關節(jié)的潤滑條件，往往顯著降低振動機構的動態(tài)響應。弗洛雷斯等人 [ 9 ]分析了曲柄滑塊連桿機構的動態(tài)行為建模，連桿 –滑塊聯合作為一個干接點無摩擦，干摩擦接觸和混合模型考慮了小偏心率和高偏心率的摩擦接觸的流體動力潤滑。他們的研究結果表明，由于在動態(tài)響應的振蕩幅度，干接觸摩擦模型比干接觸模型無摩擦更現實的。對于混合模型，他們獲得的結果顯示最小的振蕩的動態(tài)響應，但在文獻中沒有結果能支持這一觀察。此后不久，弗洛雷斯等人 [ 10 ]分析了滑塊曲柄機構考慮 –干接觸影響的動力學，在連桿 –滑塊摩擦和流體動力潤滑鉸鏈。得到的結果與干接觸模型無摩擦表現出較高的振蕩機構的動態(tài)響應，使該模型比摩擦干接觸模型不現實。當軸承模型被認為是，結果是那些與理想的滑塊曲柄機構鉸鏈獲得非相似。erkaya，su ,uzmay [ 11 ]進行了運動學和動力學的一個改進的滑塊曲柄機構和一個額外的 –偏心連桿之間的連桿和曲柄銷。他們將滑塊曲柄機構開發(fā)得到的結果進行了評價，同時與常規(guī)滑塊曲柄機構比較。這種比較表明，雖然傳統(tǒng)的和改進的滑塊曲柄機構的行程和缸內氣體壓力的改性機理相同，具有比常規(guī)的機制更高的扭矩輸出。khemili 和 Romdhane[ 12 ]分析了平面柔性滑塊曲柄含間隙機構的動力學行為，進行仿真和實驗測試，比較利用亞當斯軟件模擬的數值實驗結果。他們發(fā)現，間隙的存在影響了系統(tǒng)的動態(tài)響應，和耦合器的靈活性作為機制的懸掛系統(tǒng)。 estupi? 和 santos[ 13 ]建立一個線性往復壓縮機的數學型。他們檢查了機械系統(tǒng)的動態(tài)行為，考慮到基于多體動力學的機械部件的動力學（剛性部件）和有限元法（彈性元件）。他們還評估了流體動力軸承的影響，描述在鉸鏈處用雷諾茲方程的流體動力。他們的研究結果顯示，得到當活塞接近上死點，由于傾斜振蕩的曲柄的軌道的非線性行為的增加是最大力和最小油膜厚度，這是由曲柄銷的長度的影響。這項工作涉及的滑塊曲柄機構和軸承的流體動力潤滑條件在滑塊連桿鉸鏈處的動態(tài)行為分析。為此，分析是基于兩個不同的模型，其應用程序依賴于潤滑條件（偏心值）的軸承。第一個模型（1DOF）認為，連桿 –滑塊聯合作為一種理想的聯合應用時，連桿的一端與軸承表面接觸 [ 14 ]。第二個模型（3DOF）時，連桿一端的滑塊潤滑孔的間隙。在這種情況下，在關系到連桿滑塊，表征問題的多個自由度，通過流體的流體動力潤滑連桿和滑塊之間發(fā)生的相互作用。因此，解決的辦法是從一個由滑塊曲柄機構 –理想鉸鏈與連接桿 –滑塊聯合軸承滑塊曲柄機構的模型的混合模型得到的。因此，在這種分析的一個重要參數是考慮最小油膜厚度在流體動力潤滑。根據弗洛雷斯 [ 9 ]，高偏心率（低最小油膜厚度），壓力使表面的彈性變形，這可以是相同的順序作為潤滑油膜厚度。這些情況與那些在流體動力潤滑獲得一個更現實的分析，可以基于彈流潤滑理論。因為這個原因，這項工作被認為是最小油膜厚度超過 10%的徑向間隙的流體動力潤滑狀態(tài)，這代表 0.9 的偏心率。從混合模型得出的曲柄滑塊機構的動態(tài)響應與從常規(guī)的模型相比（理想滑塊曲柄機構鉸鏈）。此外，在軸承的壓力分布一個周期期間獲得的潤滑條件。它強調的是，通過與 Bannwart [ 12 ]合作以前開發(fā)相比，作者在目前的工作中使用的流體動力學模型更為重要。2.方法在這項工作中，一個平面曲柄滑塊機構模型，被用來確定其動態(tài)行為。然而，滑液壓動壓軸承被認為是在連桿滑塊的一個鉸鏈，代替?zhèn)鹘y(tǒng)的曲柄滑塊機構 –。這種假設認為這樣當軸承不考慮間隙時，銷不局限于只有一個方向的運動。本節(jié)描述的混合模型（自由度）的平面曲柄滑塊機構，考慮連桿 –滑塊聯合滑動軸承，這是用來分析在潤滑條件下的動態(tài)行為。該數學模型用于分析接觸狀況的動態(tài)行為是傳統(tǒng)的滑塊曲柄機構（1 自由度），是 Doughty [ 14 ]提出的。2.1 平面曲柄滑塊機構的運動學分析平面曲柄滑塊機構的運動學研發(fā)考慮了圖 1 中的方案。圖 1a 顯示了平面滑塊曲柄機構和圖 1b 描繪了自由體圖的連桿銷和滑塊孔。R 是曲柄的長度，L 是連接桿的長度，Q 是曲軸的角位移，一個是連桿的角位移，XP 和 YP 是連桿銷線性位移，XPT 和 YPT 是滑塊的直線位移， FXP 和 FYP 是分別在 x 和 y 方向上的流體動力。圖 1c 顯示的是連接的平面滑塊曲柄機構桿液壓動壓軸承軸承。OH 是中心的軸承，OP 是連桿銷中心，RH 是軸承的半徑，RP 是連桿銷半徑，E 是偏心的，hmix 是最小油膜厚度，hmax 是最大油膜厚度。此外，偏心率（ε）之間的偏心率（E）和徑向間隙（CR），在徑向游隙的軸承半徑和連桿銷半徑之間的差異根據這個數字，滑塊的銷的位置可以被表示為公式 1 所示： Xp = R:cos（q）+ L:cos（A） (1)Yp = R:sin（q）?L:sin （A） :公式 1 對連桿銷，速度和加速度的派生公式，得到： 圖 1.平面滑塊曲柄機構，（a）普通視圖，（b）擴展視圖，（C）的流體動力軸承在[ KP ]是的滑塊銷速度系數矩陣，和[LqP]和[LAP]是滑塊的銷的速度系數的偏導數矩陣，定義為：2.2. 對曲柄連桿系統(tǒng) –質心運動學分析這種分析的目的是確定的組件的質量中心的動態(tài)行為，使它找到一個滑塊曲柄機構 [ 14 ]的運動方程。圖 2a 和 b 顯示的分別是曲柄連桿與曲柄質心桿子系統(tǒng)和連桿。Upm 和 VPM 是曲軸的質量中心的坐標參考系統(tǒng)（PM）在位于曲柄（UM，VM），UPb 和 VPb 是質心坐標（PB）的參考系統(tǒng)位于連桿（UB，VB），XPM 和 YPM 的曲軸的質量中心直線位移（PM）在慣性參考系統(tǒng)（X，Y）和 XPB 和有排搏是連桿的質心線位移（PB）在慣性參考系統(tǒng)（X，Y）。根據圖 2a，曲柄的質量中心位置可以被描述為：速度可寫式的派生形式（4）為：考慮在圖 2b 的方案，得到同樣的連桿質心位移和速度。圖 2.曲柄連桿系統(tǒng)（a）對曲軸質量中心子系統(tǒng)，（b）在質量上的連桿中心子系統(tǒng)為了確定曲柄連桿系統(tǒng)的動態(tài)行為，線性和角的這個子系統(tǒng)的所有組件的速度必須分開 [ 14 ]。這些速度是需要確定的動能（旋轉和平移）的子系統(tǒng)，因此，采用拉格朗日方法得到的運動方程。因此，根據公式 7 和公式 9，它可能是寫：在[KC]代表曲柄連桿系統(tǒng)的廣義速度系數矩陣。導出[KC]相對獨立的變量 Q 和，有：在[LQ]和[LA]是廣義速度系數的偏導數矩陣。最后，分別考慮到 MM，MB，IM E IB 曲柄連桿質量，質量，轉動慣量和曲柄連桿轉動慣量，質量矩陣可以得到：2.3．在曲柄連桿系統(tǒng)的廣義力通過虛功的概念，得到了廣義的激振力，在這里提出的混合模型的情況下，由于流體動力潤滑在連桿銷 –滑塊間隙應用于耗散粘滯力。根據強 [ 14 ]，廣義力可確定使用廣義速度系數矩陣的連接點，在外力施加的條件下。因此，廣義力是由公式 14 給出了：由于液壓動力施加在滑動鉸鏈（FXP 和 FYP）中，導致廣義力發(fā)生。在這項工作中，液壓動力是由滑動軸承的潤滑模型確定的振蕩運動，通過邦瓦特 [ 13 ]的提議。2.4．的曲柄連桿系統(tǒng)的動能和勢能曲柄連桿系統(tǒng)動能的質量矩陣和廣義速度系數矩陣的確定，根據式 15：圖 3 顯示考慮到這些組件的質量中心，子系統(tǒng)由曲柄、連桿組成。根據圖 3A，潛在的能量可以寫成：3.結果本研究的目的是確定模型的適用性。計算機模擬在 FORTRAN?執(zhí)行，允許在流體動力軸承的潤滑狀態(tài)分析。圖 4 中的圖表明求解過程在動態(tài)分析中的應用。首先，在軸承的流體動力分析的初始條件。連桿銷和滑塊加速度的確定和液壓動力計算。因此，對初始條件求解運動方程，并與加速度的水動力進行了評價。因此，連桿的銷和滑塊的位移和速度計算到下一個時間步長（T +ΔT）。在確定了新的速度和位移，銷軸承偏心必須驗證。如果偏心率小于 0.9 時，該銷是在潤滑條件下。在這種情況下，新的液壓力計算，先后對銷和滑塊的加速度的重新確定。在流體動力潤滑條件下，考慮在第 2.7 節(jié)潤滑模型，對系統(tǒng)的動態(tài)行為進行評估，。然而，當滑塊銷超過閾值（偏心 ε= 0.9），它被認為是與軸承的表面接觸。在這種情況下，它是假定系統(tǒng)的動態(tài)行為是由傳統(tǒng)的滑塊曲柄機構的，在連桿 –滑塊聯合沒有間隙。因此，滑塊的曲柄連桿機構 –常規(guī)模型應用到一個反向運動發(fā)生在滑塊。在滑塊的反向運動后，滑塊銷葉片軸承表面，改變潤滑條件。至于運動微分方程的分辨率，傳統(tǒng)的模型并不復雜，因此可以解決比較快（在 2 GHz 處理器，3 GB RAM 的計算機系統(tǒng) 10 分鐘）。然而，由于高數值剛度的潤滑條件下的運動方程是高度復雜的。因此，一個特定的數值積分方法必須根據這些方程的特征選擇。在這項工作中，對 3 自由度的解決方案?；旌夏Ｐ褪峭ㄟ^使用一個僵硬的初始值問題的積分多步法的發(fā)現，rowmap 方法 [ 8 ]。然而，盡管這種方法的效率，對微分方程的解還需要計算成本高（約 48 小時，在 2 GHz 處理器，3 GB RAM 的計算機系統(tǒng)）。列于表 1 和表 2 參數分別表示軸承的幾何形狀和 –滑塊曲柄機構。從傳統(tǒng)的內燃發(fā)動機，得到了這些參數。表 3 給出的初始條件，在計算機模擬中考慮并與文獻一致，在操作過程中發(fā)現了一些實驗值。仿真我們開始考慮在表 1 和表 2 給出的物理參數，初始條件列于表 3 和 20μM 軸承徑向間隙?；瑝K的曲柄連桿機構的動力學與連桿軸承滑塊鉸鏈的分析在這項工作中的計算機模擬，考慮流體動力潤滑的閾值條件。因此，流體動力潤滑條件下高達 0.9 的偏心率是極低的膜的厚度，這是流體動力潤滑不充分的條件。因此，最大限度的模擬偏心率的值是 0.9，這表明銷在假設是接觸支承面上不是液壓動力條件。圖 5a 和 b 分別顯示在第一和第二圈獲得的軌道。如圖 5a 表示，銷的初始偏心率是 0.6。當運動開始，銷向左移動直到超過偏心極限，之后它可能仍然在與軸承表面接觸到滑塊的運動倒置。在這種情況下，即，當銷與軸承的接觸面，動力學是基于傳統(tǒng)的滑塊曲柄機構的 [ 14 ]的數學模型分析。當運動倒置，銷從軸承表面的軸承上的區(qū)域。表一流體動力軸承的參數。參數和單位 符號 數值動水軸承半徑[m]軸的半徑(連桿銷)[m]水動力軸承寬度[m]徑向間隙[μm]絕對粘度[Pa.s]質量密度[kg/m3] 表二曲柄滑塊機構的參數參數和單位 符號 數值曲柄的長度[m]連桿長度[m]曲柄慣性[kg.m]連桿慣性[kg.m]曲柄質量[kg]連桿質量[kg]滑塊質量[kg]曲軸質量中心在軸的位置[m]曲軸中心在軸虛擬機中的位置[m]該連接在軸 UB 桿質心位置[m]該連接在軸的 VB 桿質心位置[m]表三采用計算機模擬的初始條件。初始條件和裝置 符號 模擬 1 模擬 2 模擬 3曲柄的角位置[rad]曲柄角速度[rad/s]連桿的角位置[rad]連桿角速度[rad/s]滑塊的直線位置[m]滑塊的線性速度[m] 接近表面的軸承，銷和滑塊推向相反的方向上的流體動力，使它們相對于彼此移動。這使得銷回到軸承的下部區(qū)域，再次超過偏心極限。銷并保持與表面接觸到底圈。在第二圈，活塞銷的動態(tài)行為類似于第一圈。圖 6 說明了曲軸的動態(tài)行為。可以看出，曲柄的速度和加速度的兩個系統(tǒng)，即是相似的，滑動軸承滑塊曲柄機構（3 自由度）和常規(guī)滑塊曲柄機構（1 自由度）。然而，請注意，速度和加速度都傾向于作為一個在 3 自由度模型的時間函數的兩個應用模型之間的過渡的增加。圖 5.軌道滑塊的聯合模擬, (a)第一圈,(b)第二圈圖 6 曲柄分析模擬 (a)角速度,角加速度(b)圖 7.連桿分析模擬 (a)角位移、角速度(b)當銷超過偏心極限，最后重點分析了在水動力條件的 1 自由度的初始條件。傳統(tǒng)的機制（接觸）。因此，在這一點上的運動，是由上述自由度模型給定初始條件的動態(tài)行為描述。因此，在混合模型與傳統(tǒng)的過渡，輸入的速度不同，給出了兩種模型在 t = 0 的初始條件。因此，主要作用是在模型的每個過渡曲柄的角速度的增大，這是相對于 1 自由度位移的傳統(tǒng)的模型?？梢钥吹皆趫D 6 b,曲軸的角加速度小斷層與傳統(tǒng)的模型。雖然曲柄不能直接連接到水動力軸承,這些不再生產的動態(tài)響應是由銷在潤滑條件。相同的行為由施瓦布報道,Meijaard 和 Meijers[8]。根據 Schwab[8],與液壓動力軸承的動態(tài)行為曲柄滑塊機構在曲柄滑塊聯合常規(guī)機制類似,所代表的動態(tài)響應的平滑曲線。然而,反應略有不同,只有當銷進入軸承。這穿越涉及高速度和期刊中心,因此, 它增加了高峰。圖 7 描述了連桿的位移和速度。注意連桿的位移、速度和軸承滑塊曲柄機構（3 自由度）和常規(guī) –滑塊曲柄機構（ 1 自由度）是相似的。圖 8 所示，滑塊的動態(tài)行為的兩個系統(tǒng)是相似的（1DOF 和 3 自由度）。圖 9 顯示了連桿的角加速度和滑塊的直線加速。注意連桿和滑塊的加速度與軸承滑塊曲柄機構達到（3 自由度）類似于傳統(tǒng)的滑塊曲柄機構實現加速度（1自由度）。然而，在與軸承表明振蕩由于在流體動力潤滑條件的鉸鏈行為的滑塊曲柄機構達到加速度。圖 8.滑塊分析模擬 (a)線性位移、(b)線性速度圖 9.在模擬加速度分析 ,(a)連桿加速度, (b)滑塊加速度如圖 9 顯示，在潤滑時更重要的是在滑塊的加速度和連桿軸承的原因引起不連續(xù)，因為這些組件的交互。當銷是不是在潤滑條件下，它被認為是接觸支承面。的數學模型，用來分析在接觸條件的動態(tài)行為是傳統(tǒng)的滑塊曲柄機構模型。因此，機構的動態(tài)行為是傳統(tǒng)。4.結論這項工作涉及的滑塊曲柄機構和在潤滑條件下滑塊鉸鏈的動態(tài)行為分析。為此，考慮連桿滑塊鉸鏈的流體動力軸承，建立一個曲柄滑塊連桿機構的數學模型，。在這項工作中，數學模型的計算機模擬了連 –滑塊鉸鏈銷的動態(tài)行為。結果表明，在操作過程中進行的銷兩個條件。在某些時刻銷進行流體動力潤滑條件和別人是在與軸承的表面接觸。數值結果表明，不論初始條件軸承參數多少，銷斷裂的油膜，進入與軸承的表面接觸，如前所述的規(guī)定，達到一個高的偏心銷，使油膜壓力甚至分裂它。這種行為表明，彈流潤滑模型的混合模型接觸條件的使用（3 DOF。）可能更適合本文的分析，在這種情況下，直接接觸可以防止甚至當銷不動壓潤滑。然而，評估的動態(tài)響應和壓力分布潤滑時，所提出的混合模型可以應用。的彈流模型的實現需要從流體的混合模型彈流狀態(tài)過渡的進一步研究。對于軸承的徑向間隙，可以得出結論，在徑向間隙的增加會導致油膜厚度隨之增加，降低水動力和促進與軸承表面引腳的接觸。相反，在初始角速度的增加導致的流體動力作用于銷增加，因為這些力量的初始角速度成正比。這些較強的液壓動力增加銷的軸承表面的距離，防止觸碰這面。連桿滑塊鉸鏈是相對復雜的流體動力軸承由于使用的兩種流體動力潤滑和滑塊曲柄機構 –模型的非線性方程組 –滑塊曲柄機構問題的解決方案。這組方程提出了高數值的剛度，使數值積分方法在微分方程的成功解決問題的一個重要因素的選擇。致謝作者希望感謝 FAPESP，CAPES 和 CNPQ 提供的資金支持這項工作的。參考文獻[ 1 ] M. Takiguchi, M. Oguri, T. Someya,，旋轉在汽油發(fā)動機活塞銷的運動研究，SAE 論文 938142，底特律，美國，1993.[ 2 ]M. Takiguchi, K. Nagasawa, T. Suhara, M. Hiruma,摩擦小端連接一個汽車發(fā)動機連桿軸承潤滑特性，秋天技術會議 ASME 2 卷（1996）1 –6。T. Suhara, S. Ato, M. Takiguchi, S. Furuhama,，本文對古浜莊，活塞銷座軸承的汽車發(fā)動機的摩擦與潤滑特性，SAE 論文 970840，底特律，美國，1997。[ 4 ]C. Zhang, H.S. Cehng, L. Qiu, K.W. Knipstein, J. Bolyard,，活塞銷/行為劃傷孔軸承的混合潤滑，部分。1，實驗研究，賽寶反。46 卷（2004）193–199。。[ 5 ]C. Zhang, H.S. Cehng, J.O. Wang，活塞銷/行為劃傷孔軸承的混合潤滑，部分。2，實驗研究，賽寶。反。47 卷（2004）149 –156。丹尼爾的學生，,機械原理 46（2011）1434 –14521451[ 6 ] J.L. Ligier, P. Ragot，活塞 –銷：磨損和旋轉運動，SAE 紙 2005-01-1651，美國，底特律，2005。[ 7 ] J.L. Ligier, P. Ragot，小端連桿潤滑， SAE 紙 2006-01-1101，美國，底特律，2006。[ 8 ] A.L. Schwab, J.P. Meijaard, P. Meijers，在剛性和彈性力學系統(tǒng)的動態(tài)分析的旋轉關節(jié)間隙模型的比較，機械原理（9）37（2002）895 –913。[ 9 ] P. Flores, J. Ambrosio, J.P. Claro，潤滑關節(jié)的平面多體機械系統(tǒng)的動態(tài)分析，多體系統(tǒng)動力學 12 卷（2004）47 –74。[ 10 ] P. Flores, J. Ambrosio, J.C.P. Claro, H.M. Lankarani, C.S. Koshy，在包括間隙和潤滑關節(jié)的機械系統(tǒng)動力學研究，機械原理 41 卷（2006）247 –261。[ 11 ] S. Erkaya, S. Su, I. Uzmay 的滑塊曲柄機構， –偏心連接器和行星齒輪動態(tài)分析，機械原理 42 卷（2007）393 –408。[ 12 ] I. Khemili, L. Romdhane，一個靈活的曲柄滑塊 –含間隙機構的動力學分析，歐洲/固體力學學報 27 卷（2008）882 –898。[ 13 ] E.A. Estupi?an, I.F. Santos，建模的封閉式壓縮機使用不同的約束方程，以適應多體動力學和流體動力學潤滑，對機械科學與工程第三十一號 1巴西社會雜志（2009）35 –46。[ 14 ] S. Doughty，力學的機器，約翰威利父子，美國，1988。Analysis of the dynamics of a slider–crank mechanism with hydrodynamiclubrication in the connecting rod–slider joint clearanceGregory B. Daniel, Katia L. Cavalca Department of Mechanical Design, University of Campinas, Faculty of Mechanical Engineering, Postal Box 6122, 13083-970 Campinas, SP, BrazilabstractThe conventional slider-connecting rod–crank mechanism is widely applied in mechanical systems. The use of hydrodynamic bearings in the mechanism joints is of particular interest in reducing friction, mainly in special conditions of lubrication such as the connecting rod–slider joint. This bearing belongs to a class of bearings with alternating rotational motion. This paper proposes a mathematical model for this particular problem, considering the dynamics of the slider-connecting rod–crank system interacting with the lubrication phenomenon in bearings with alternating motion. Two models were used to analyze the dynamics of the system. The first model (by Eksergian Equation of Motion) represents the system when the connecting rod end is in contact with the bearing surface, assuming, in this condition, the same behavior as that of rigid bearings (without clearance). The second model (by Lagrange Method) represents the system when the connecting rod end is in the hydrodynamic lubrication mode in the slider bore clearance. In this condition, the slider moves in relation to the connecting rod, presenting a problem of multi-degrees-of-freedom. The mathematic model of hydrodynamic lubricationwas introduced to obtain more realistic results of the system's dynamic behavior.1. IntroductionThe lubrication system is a crucial element in the optimal performance of machines in general in view of its numerous functions, e.g., lubricating and protecting components, reducing friction, and cleaning and cooling internal mechanisms. The operation of this system requires proper calibration. Therefore, excessive or insufficient lubrication influences the dynamic behavior of the machine and can seriously damage the mechanisms involved. The mechanisms of many machines contain numerous hydrodynamic bearings, but the hydrodynamic bearing of the connecting rod–slider joint is particularly important because it is part of a new class of hydrodynamic bearing called a Hydrodynamic Bearing with Alternating Motion. Unlike traditional hydrodynamic bearings, this type of hydrodynamic bearing does not make a complete rotation. Therefore, recent years have seen a growing need to study the behavior of this specific type of hydrodynamic bearing and its influence on the dynamic behavior of the machine. Most of the research into the hydrodynamic bearing of the connecting rod–slider joint has so far been conducted by the Musashi Institute of Technology. In fact, this institute has developed and constructed devices to investigate lubrication and friction in the hydrodynamic bearing of slider pins for more than 15 years. Takiguchi et al. [1], for example, studied the rotating motion of a floating piston pin applied in gasoline-powered automotive engines. Three years later, Takiguchi et al. [2]developed a measuring device that determines the status of lubrication in the hydrodynamic bearing of the connecting rod–slider joint based on measured friction force. In another study conducted at the Musashi Institute of Technology, Suhara et al. [3]examined lubrication conditions in the hydrodynamic bearing of the slider pin of gasoline-powered automotive engines, considering as parameters of analysis the length, internal diameter and material of the piston pin. More recently, Zhang et al. [4,5]developed tools to investigate wear in hydrodynamic bearings of the piston pin. In 2005, Ligier and Ragot [6]analyzed the behavior of the hydrodynamic bearing of the piston pin. The following year, these authors [7] presented a general view of the operation of the hydrodynamic bearing of the connecting rod–piston joint in four-stroke engines, emphasizing the oil feed in the bearing.As mentioned earlier, hydrodynamic bearings with alternating motion are part of slider–crank mechanisms. Hence, a reliable mathematical model of this mechanism must take into account the behavior of the bearing in the connecting rod–slider joint. Owing to the wide application of the slider–crank mechanism, many studies have focused on devising a mathematical model and analyzing the dynamic behavior of this mechanical system. Schwab, Meijaard and Meijers [8]compared the dynamic behavior of the slider–crank mechanism, considering different connecting rod–slider joint models such as a Hertz contact model with dissipation, an impact modeland a hydrodynamic bearing model. In the aforementioned work, the authors also compared the dynamic behavior of the mechanism considering the connecting rod as an elastic and rigid component. The results obtained by those authors show that the assumption of a connecting rod as an elastic component, as well as the lubrication condition in the joint, tends to significantly decrease the vibrations in the dynamic response of the mechanism. Flores et al. [9] analyzed the dynamic behavior of the slider–crank mechanism, modeling the connecting rod–slider joint as a dry contact without friction, a dry contact with friction and as a hybrid model which considershydrodynamic lubrication for small eccentricity and dry contact with friction for high eccentricity. Their results show that the dry contact model with friction is more realistic than the dry contact model without friction, due to the magnitude of the oscillations in the dynamic response. As for the hybrid model, the results they obtained showed the smallest oscillation in the dynamic response, but there are no results in the literature that can support this observation. Soon thereafter, Flores et al. [10]analyzed the dynamics of the slider–crank mechanism considering the effects of dry contact, friction and hydrodynamic lubrication on the connecting rod–sliderjoint. The results obtained with the dry contact model without friction showed high oscillations in the dynamic response of the mechanism, making this model less realistic than the dry contact model with friction. When the model with hydrodynamic bearing was considered, the results were very similar to those obtained in the slider–crank mechanism with ideal joints. Erkaya, Su and Uzmay [11]analyzed the kinematics and dynamics of a modified slider–crank mechanism with an additional eccentric link between the connecting rod and crank pin. The results they obtained with the slider–crank mechanism they developed were evaluated and compared to that of a conventional slider–crank mechanism. This comparison indicated that although the conventional and modified slider–crank mechanisms have the same stroke and the same gas pressure in the cylinder, the modified mechanism has a higher torque output than the conventional mechanism. Khemili and Romdhane [12]analyzed the dynamic behavior of a planar flexible slider–crank mechanism with clearance, performing simulated and experimental tests and comparing the experimental results withthe numerical results of the simulations using ADAMS software. They found that the presence of clearance affects the dynamic response of the system, and that the coupler's flexibility acts as a suspension system in the mechanism. Estupi?an and Santos [13] developed a mathematical model for a linear reciprocating compressor. They checked the dynamic behavior of the mechanical system, considering the dynamics of the mechanical components based on multibody dynamics (rigid components) and the finite element method (flexible components). They also evaluated the influence of the hydrodynamic bearing, describing the hydrodynamic forces in the joints using the Reynolds equation. Their results showed that maximum forces and minimum film thickness are obtained when the piston is close to the top dead center and that the nonlinear behavior of the orbits increases due to the tilting oscillations of the crank, which are influenced by the length of the crank pin.This work involves an analysis of the dynamic behavior of the slider–crank mechanism and the hydrodynamic lubrication condition of the bearing in the connecting rod–slider joint. For this reason, the analysis is based on two different models, whose application depends on the lubrication conditions (eccentricity value) of the bearing. The first model (1dof) considers the connecting rod–slider joint as an ideal joint and is applied when the connecting rod end is in contact with the bearing surface [14]. The second model (3dof) is applied when the connecting rod end is in hydrodynamic lubrication in the slider bore clearance. In this condition, the slider moves in relation to the connecting rod, characterizing a problem of multiple degrees of freedom in which the interaction between the connecting rod and the slider occurs through the hydrodynamic forces of the lubricant fluid. Therefore, the solution is obtained from a hybrid model composed of a model of the slider–crank mechanism with ideal joints and a model of the slider–crank mechanism with hydrodynamic bearing in the connecting rod–slider joint.Thus, an important parameter in this analysis is the minimum oil film thickness considered in the hydrodynamic lubrication. According to Flores [9], for high eccentricity (low minimum oil film thickness), the pressure causes elastic deformation of the surfaces, which can be of the same order as the lubricant film thickness. These circumstances differ from those obtained in hydrodynamic lubrication and a more realistic analysis can be made based on the theory of elastohydrodynamic lubrication. For this reason, this work considered the hydrodynamic lubrication condition for a minimum oil film thickness exceeding 10% of the radial clearance, which represents an eccentricity ratio of 0.9.The dynamic responses of the slider–crank mechanism obtained from the hybrid model are compared with those obtained from the conventional model (slider–crank mechanism with ideal joints). Moreover, the pressure distribution in the bearing isobtained during one period of lubrication condition. It is important to emphasize that the hydrodynamic model used in the present work was developed previously by the authors in cooperation with Bannwart [12].2. MethodologyIn this work, a planar slider–crank mechanism was modeled in order to determine its dynamic behavior. However, instead of the conventional slider–crank mechanism, a hydrodynamic bearing was considered in the connecting rod–slider joint. Thisassumption considers that the pin is not restricted to only one direction of motion, which is the case when the bearing is considered without clearance.This section describes the hybrid mathematical model (3DOF) of the planar slider–crank mechanism considering a hydrodynamic bearing in the connecting rod–slider joint, which is used to analyze the dynamic behavior in the lubricationcondition. The mathematical model used to analyze the dynamic behavior in the contact condition was the conventional slider–crank mechanis