接觸有限元分析及應(yīng)用

王月宏

（陜西能源職業(yè)技術(shù)學院 機電與信息工程學院，咸陽 712000）

1 接觸有限元分析

彈性接觸問題隸屬于邊界非線性問題，它所涉及的問題既有因接觸區(qū)內(nèi)的變化而引起的非線性問題，又有因接觸壓力的分布變化而引起的非線性問題和接觸間由于摩擦存在而產(chǎn)生的非線性問題，求解過程可以認為是反復迭代尋求精準接觸狀態(tài)的過程。因此，在求解之前先設(shè)定一個可能的潛在接觸狀態(tài)，然后把設(shè)定條件帶入求解方程，可以求解得到具體接觸位置的接觸內(nèi)力和位移，并由此結(jié)果判斷是否滿足條件。當求解值不滿足所需的接觸條件時，修改接觸點起初的接觸狀態(tài)進行重新求解，直至所有接觸點的接觸狀態(tài)都滿足接觸條件。

物理特征接觸狀態(tài)過程是依賴于時間，并且伴隨著材料的非線性以及幾何形狀的非線性變化的一個不確定過程。在一些特殊的接觸過程中，物體接觸界面的求解參數(shù)是未知的，如接觸區(qū)域、形狀和它的動力學和運動學的狀態(tài)等。這些求解的特點，決定了接觸問題一般采用增量法求解。

運用增量求解方法時，首先將施加載荷分為若干個小單元f0，f1，…fm，那么相應(yīng)的變形位移也可以分為若干小單元a0，a1，……am。它的一般步驟是假設(shè)第m步的載荷fm和相對應(yīng)的位移am已知，而下一步施加的載荷增加為fm+1=(fm+Δfm)，再求解am+1=(am+Δam)。假設(shè)每一步的載荷增量足夠小，那么解函數(shù)的收斂性可以保證。同時,可以求解得到在加載過程中每一階段的中間數(shù)值結(jié)果，以便研究結(jié)構(gòu)的位移和應(yīng)力等參數(shù)隨著載荷變化的具體情況。

接觸條件的約束關(guān)系均是不等式的，所以其被稱為單邊約束關(guān)系。此外，接觸區(qū)域的范圍和接觸狀態(tài)都是未知的，怎樣將接觸面條件巧妙地引入求解過程是求解接觸問題的關(guān)鍵。針對非線性接觸問題的特殊性，一般的求解過程都要經(jīng)過試探-校核的嘗試迭代方法進行逐步求解，每一個增量步的求解迭代過程可以歸納為：

（1）根據(jù)上一增量步的求解結(jié)果和當前增量步所給定的載荷條件，經(jīng)過接觸狀態(tài)條件的檢查和搜索，假定此增量步的第一次迭代的接觸面和接觸狀態(tài)已知。這里的接觸狀態(tài)一般有“滑動”狀態(tài)或“粘著”狀態(tài)。

（2）根據(jù)第一步的假設(shè)，對于接觸區(qū)域上的所有點，把運動學或者動力學的理論關(guān)系上的不等式約束條件改為等式約束條件，作為定解條件帶入方程并進行求解。

（3）利用接觸區(qū)域上的初次計算結(jié)果和步驟（2）中所提及的運動學或者動力學的不等式約束關(guān)系，對初步假設(shè)的某一種接觸狀態(tài)進行檢查校驗。如果接觸區(qū)域內(nèi)的所有點都滿足校核條件，即完成了本增量步的求解過程并進入下一個增量步的計算，否則就要審視并修改相應(yīng)的接觸狀態(tài)，返回到步驟（1）進行再次搜索和迭代求解，直至所有接觸區(qū)域上點的解狀態(tài)都滿足良好的校核條件，進而再進入下一個增量步進行求解[1]。

2 有限元軟件的接觸分析

2.1 ANSYS的接觸類型

ANSYS Workbench分析軟件共存在兩種接觸類型。它們分別是剛體-柔體類型和柔體-柔體接觸類型。剛-柔接觸，指是適用于兩個剛體之間相差較大的物體之間的接觸，一般認為剛度較大的物體為剛體。通常，一種軟（塑性）材料和一種硬（脆性）材料的接觸，即可被定義為剛-柔接觸類型。柔體-柔體之間的接觸是普遍應(yīng)用的一種接觸類型，一般需要假設(shè)兩個接觸物體都是變形體，它是指適用于彈性模量和結(jié)構(gòu)剛性都相等或者比較接近的物體之間的接觸。

2.2 ANSYS的接觸方式

ANSYS Workbench分析軟件存在三種接觸方式。它們分別是點與點的接觸方式、點與面的接觸方式以及面與面的接觸方式。點與點接觸方式主要用于模擬點和點之間的接觸行為，分析者需預先知道準確的接觸點位置，其僅僅適用于模擬分析接觸區(qū)域之間有微小的相對滑動的情況。點與面的接觸方式是指允許在接觸區(qū)域內(nèi)某一個節(jié)點和被接觸物體的接觸區(qū)域的某一單元體相接觸，接觸位置可以不用十分精確，而且接觸區(qū)域之間不需要有一致的網(wǎng)格，同時此種接觸方式允許有物體之間的熱傳導、有較大的變形和比較明顯的相對滑動等非線性行為。面與面的接觸方式，相對計算量較少，可適用于復雜的表面、大變形以及含有摩擦力的復雜接觸問題的求解。齒輪嚙合傳動過程中，接觸部位的剛度變化，會導致齒面的實際接觸部位位于接觸線附近有限的齒面上，齒面發(fā)生相互作用時，其接觸區(qū)域內(nèi)的接觸節(jié)點位置是不固定的，并且齒面的相對滑動也是存在的，因此點-點接觸方式是不適合斜齒輪嚙合接觸分析的。對于點-面接觸方式，可以把接觸面設(shè)定為一組有規(guī)律的節(jié)點，從而代替面-面的接觸方式。面-面接觸單元相比點-面接觸單元就有很大的優(yōu)勢，比如，它的接觸單元不受剛體接觸面形狀的限制，面-面接觸單元允許有接觸表面不連續(xù)性的存在；允許有大的位移滑動和摩擦的存在，協(xié)調(diào)剛度陣的計算，面-面接觸單元有不對稱的剛度選項。

2.3 ANSYS的接觸算法

ANSYS Workbench分析軟件有三種求解非線性接觸問題的算法：罰函數(shù)法（Penalty method）、拉格朗日乘子算法（Lagrange method）和增廣拉格朗日乘子法（Augmented Lagrange method）。這里主要介紹拉格朗日乘子算法。

為了滿足無穿透的良好接觸條件，拉格朗日乘子算法增加了一個單獨變量，即接觸壓力，也就是說它增加了一個自由度。此算法不需要人為設(shè)定初始接觸剛度大小，它的穿透量為零，也就是說，它是可以直接實現(xiàn)真實良好接觸條件的，計算結(jié)果比較精確。但是，該方法由于新增加的接觸壓力，就增加了整個求解系統(tǒng)的變量數(shù)目，并且對角線元素全部為零的子矩陣也存在于剛度陣中，為了保持計算精度，必須施加附加的操作，這就給整個計算過程帶來困難。在嚙合的過程中，一對圓柱直齒齒輪的接觸狀態(tài)是不斷發(fā)生變化的，并且接觸力也有突變，接觸狀態(tài)會出現(xiàn)周期性振動式的交替變化，單純的拉格朗日算法無法進行求解來滿足接觸狀態(tài)周期性變化和突變條件。此算法主要適用于一些特殊單元界面的接觸分析問題。該算法也有不足之處，接觸物體之間存在相對運動，但是它必須限制相對運動量，并且發(fā)生的接觸點準確位置需要預先知道，人們才能夠施加界面單元。

3 結(jié)論

非線性接觸問題具有復雜性，筆者研究了接觸有限元方法，得出求解非線性接觸問題的迭代步驟，并分析了有限元軟件的接觸類型、方式以及算法，然后以一對嚙合齒輪為例，利用ANSYS Workbench軟件進行接觸強度分析，

驗證了該方法的可靠性。