ADAMS中的接觸和接觸摩擦作用機制實例詳解

高一佳

（陜西保利特種車制造有限公司，陜西 西安 710200）

ADAMS中的接觸和接觸摩擦作用機制實例詳解

高一佳

（陜西保利特種車制造有限公司，陜西 西安 710200）

文章首先從理論上詳細說明了 ADAMS軟件中的接觸和接觸摩擦在模型中的作用機制。然后結合拉臂式垃圾車動力學分析實例，研究了各參數變動對接觸造成的影響。最后給出了解決接觸求解無法收斂問題的一般原則。

ADAMS；接觸；接觸摩擦

CLC NO.:U461Document Code:AArticle ID:1671-7988 (2017)06-64-04

前言

即使是經驗豐富的分析工程師，在處理包含接觸和接觸摩擦的ADAMS仿真任務時，也時常遭遇頻繁的求解失敗。是什么原因讓日常生活中無處不在的接觸和摩擦在 CAE分析中變成了一個麻煩制造者。本文將試圖從理論與實際相結合的角度闡述ADAMS中接觸和接觸摩擦的作用機制及造成求解困難的原因，并給出改善接觸求解困難的一些建議。

1、什么是接觸

在ADAMS中，當兩個分離的表面互相碰觸并互切時，就稱它們處于接觸狀態。在一般的物理意義中，處于接觸狀態的表面有下列特點：

·不互相穿透；

·能夠傳遞法向壓力和切向摩擦力；

·通常不傳遞法向拉力。（因此，它們相互間可以自由地分開并遠離）。

圖1 接觸的三種狀態

接觸由如圖1所示三種狀態。在ADAMS中分別表述為Open，Closed and Stick，Closed and Slide。接觸具有強烈的非線性、非保守特性。隨著接觸狀態的改變，接觸表面的法向和切向剛度都有顯著的變化。剛度方面大的突變通常會導致嚴重的收斂困難，特別是當存在滑動時。

大多數的接觸問題還需要同時考慮摩擦，摩擦計算同樣

2、ADAMS中接觸力的計算

Adams/Solver有兩個幾何引擎用來檢測三維接觸，分別是Parasolid和RAPID。Parasolid是來自EDS/Unigraphics的一個幾何套件，支持到 Parasolid14.0。但是目前 RAPID是Adams的默認幾何引擎，當前支持到2.01版本。幾何引擎負責檢測兩個幾何體間的接觸，定位接觸點，并且計算接觸點的公法線。

在Adams中，有兩種計算接觸力的方法：基于回歸的接觸算法(Restitution-base contact)和基于碰撞函數的接觸算法(IMPACT-Function-based contact)。回復系數被用在一些特殊場合，就是有時候無法從材料手冊或物理實驗中獲得碰撞數值，但是可以獲得回復數值。一般來說，碰撞方法在數值上更平順，仿真速度也更快。只有在回復數值可以獲得但碰撞數值無法獲得的情況下才使用回復系數方法。下面只對應用較廣的IMPACT型接觸作出說明，基于IMPACT的接觸模型如下圖2所示。

圖2 接觸碰撞模型及參數設定

圖3 接觸力計算參數示例

Adams根據Hertz接觸理論，采用impact函數提供的非線性等效彈簧阻尼模型計算接觸力，接觸力計算參數如圖3所示。接觸力如下式 1-1所示，各參數的物理意義如圖3所示。

式中：

x-表示兩接觸物體的距離；

x˙-表示接觸的相對速度；

x1-表示x的自由距離。當x＜x1時，兩物體接觸并產生接觸壓力，反之不接觸；

k-表示接觸表面剛度，為一非負實數。

n-表示力的指數，為一非負實數。對于剛度大的接觸面，n＞1。反之n＜1；

cmax- 表示接觸過程中的最大阻尼系數，為一非負實數。

d-表示接觸過程中的最大阻尼系數，為一非負實數。

由式1-1可看出接觸力由兩部分組成。 k(x1- x)n為彈性力分量， -STEP(x, x1-d,cmax,x1,0)·x˙ 為粘性阻尼力分量。

在ADAMS中定義接觸時需輸入圖2右側所示的4個參數，下面將對這幾個參數的作用和物理意義作出詳細說明。

·k (Stiffness)：表示接觸表面的剛度。一般來說，剛度值越大，兩物體滲透的量越小，越接近實際情況。但是過大的剛度會造成病態矩陣和積分求解困難。如果剛度值過小，雖然能夠便于計算，但不能模擬兩個構件之間的真實接觸情況。剛度值的合理范圍可通過Hertz接觸理論進行估算。

·n (Force Exponent): 力指數，用來計算瞬時法向力中材料剛度項貢獻值的指數，表征了材料的非線性特性。對于剛度比較大的接觸，取n＞1，否則取n＜1。對于金屬常用1.3～1.5，對于橡膠可取2甚至3。一般用1.5。

·cmax(Damping): 定義接觸材料的阻尼屬性，用于表征接觸能量損失。通常取剛度值的 0.1～1﹪，一般取 10-100 N*sec/mm。C越大，滲透量曲線越平滑，碰撞力曲線越平滑。

·d (Penetration Depth): 定義全阻尼(full damping)時的最大穿透深度。理論上該值越小越逼近真實，但該值過小會影響到ADAMS的數值收斂性，一般可采用ADAMS中推薦的取值0.01 mm。注意，該值并不是最大穿透深度（阻尼達到最大值后，零件之間的相互穿透還可以繼續）。

3、ADAMS接觸中摩擦力的計算

在定義接觸時，摩擦力的定義對接觸力的仿真結果有很大影響，定義摩擦力可以使接觸力曲線更平滑，更趨近真實。接觸力中的摩擦力是接觸正壓力和摩擦系數的乘積。

圖4 摩擦系數與滑移速度之間的關系

Adams為接觸使用了一種基于速度的摩擦模型，即庫倫摩擦（Coulomb）。根據兩接觸物體的相對滑移速度的不同，摩擦在動摩擦與靜摩擦之間相互轉換。摩擦系數隨相對滑移速度的變化曲線如下圖4所示。非線性摩擦系數模型公式如下式1-2所示。

式中：是摩擦系數，是最大靜摩擦系數，是動摩擦系數，v是相對滑移速度，是靜摩擦轉變速度，是動摩擦轉變速度。在ADAMS中使用庫倫摩擦需指定如下圖5中所示的參數。

圖5 摩擦系數與滑移速度之間的關系

·Coulomb Friction: 指定CONTACT中的摩擦模型為動摩擦，而不是靜摩擦。

·Static Coefficient (MU_STATIC)，靜摩擦系數，取值范圍為≥0。當相對滑移速度小于時，為接觸點指定摩擦系數。取值過大會引起積分困難。

· Dynamic Coefficient (MU_DYNAMIC)，動摩擦系數，取值范圍為≥≥0。當相對滑移速度大于時，為接觸點指定摩擦系數。取值過大會引起積分困難。

·Stiction Transition Vel. (STICTION_ TRANSITION _VELOCITY)，靜摩擦轉換速度，取值范圍為≥≥0，單位為mm/sec。用來在庫倫摩擦模型中在接觸點計算摩擦力。當接觸點滑移速度減小時，ADAMS/SOLVER逐漸將摩擦系數從過渡到。當滑移速度等于指定的值時，有效摩擦系數是。

注意：過小的會引起數值積分困難。指定該值的一般經驗法則是＞5*ERROR。其中ERROR是積分器要求的精度，默認值是1E-3。

·——Friction Transition Vel. (FRICTION_TRANSITION _VELOCITY)，動摩擦轉換速度，取值范圍為≥≥0，單位為mm/sec。用來在庫倫摩擦模型中在接觸點計算摩擦力。當接觸點滑移速度增加時，ADAMS/SOLVER逐漸將摩擦系數從過渡到。當滑移速度等于指定的值時，有效摩擦系數是。

注意：過小的會引起數值積分困難。指定該值的一般經驗法則是＞5*ERROR。其中ERROR是積分器要求的精度，默認值是1E-3。對于大部分相對運動的實際零件，兩者之間的摩擦在絕大多數情況下都是動摩擦。所以一定要保證動摩擦系數的準確性，靜摩擦系數的準確性相對來說就不是那么重要了。靜摩擦轉變速度和動摩擦轉變速度的值一般都很小，具體值的準確性也不是太重要。

4、接觸摩擦參數設定實例

本節將以拉臂式垃圾車動力學分析為例來說明如何設置接觸和摩擦參數。仿真該作業過程需要定義三個接觸對。如圖6所示，接觸1為箱體滾輪與地面之間的接觸（鋼-混凝土），接觸2為箱體縱梁與車架導向輪之間的接觸（鋼-尼龍），接觸3為箱體吊環與拉臂鉤之間的接觸（鋼-鋼）。

圖6 拉臂車上的接觸對

滾輪的材質為鋼，路面為混凝土，掛鉤為鋼，縱梁導輪為尼龍，因此接觸對1和接觸對2為剛-柔接觸，接觸對3為剛-剛接觸。根據上述分析，設定如下接觸參數：

表1 各接觸對的接觸和摩擦參數值

為進一步說明各參數對接觸效果的影響，下面分別對各參數進行調整，檢查接觸力的變化情況，并與原始結果進行比較。結果如下圖7所示。

圖7 各接觸參數變動后接觸力變化

圖7-1為剛度k值增大3倍，最大接觸力增加了約40%。圖7-2為力指數n值增大2倍，最大接觸力增加了約40%，并在初始階段出現了振動。可看出，增大k和增大n有類似的效果。

圖7-3為最大阻尼系數Cmax值增大2倍，接觸力振幅減小了約64%。可看出增大cmax可以使接觸過程更加柔和。圖7-4為全阻尼穿透深度值增大100倍，可看出，d值的變化并不會顯著影響接觸力。在不造成求解速度降低和收斂困難的前提下，該值越小越好。

5、處理求解困難的一般原則

盡管已經充分了解了ADAMS接觸和摩擦的算法，但包含接觸和接觸摩擦的實際模型在求解時仍然時常失敗，下面給出處理此類問題的一般原則。

? 總是運用crawl-walk-run（爬-走-跑）方法。首先運行沒有接觸摩擦的模型，然后添加摩擦。

? 選擇適當的積分器。三種常用積分器GSTIFF、WSTIFF和 BDF計算的效率和穩定性各有不同，計算穩定性 BDF＞WSTIFF ＞GSTIFF，而計算效率GSTIFF ＞WSTIFF ＞ BDF。GSTIFF是系統默認的積分器。

? 選擇適當的積分格式。ADAMS 提供了三種積分格式：I3、SI2 和 SI1。I3 格式求解速度快但精度差；SI2 格式可以避免Jacobian矩陣的病態，而且考慮了約束方程。雖然求解時間變長但求解精度高；SI1 求解精度比SI2 還高，但它計算量太大，一般不采用。推薦使用使用S12公式。該公式的修正器比標準的 GSTIFF更穩定，并且也許能求解標準GSTIFF失敗的問題。

? 選擇適當的積分誤差。積分誤差決定了在求解動力學方程的過程中，某一步的預測值與校正值之間所能接受的差值。積分誤差過大，計算容易進行，但最終結果會產生過大誤差；積分誤差過小，求解時間太長。一般來講，積分誤差設為0.001 較為適宜。

? 選擇適當的接觸剛度。過大的接觸剛度會引起Jacobian矩陣病態，過小的接觸剛度又會造成相互穿透和結果誤差過大。應當在開始時設定適當的接觸剛度，然后逐漸增大該值直到結果誤差滿足要求或達到指定剛度值。

? 如果增加摩擦力引起了數字困難或者仿真減速，逐漸增大靜摩擦轉變速度和動摩擦轉變速度。減小摩擦系數通常也會對積分器有所幫助。減小阻尼系數也能夠減少接觸延遲并能夠幫助仿真。

6、結論

本文對 ADAMS中接觸和接觸摩擦作用機制在理論上進行了詳細說明和解釋，并結合具體工程實際應用，給出了設置合理的接觸和接觸參數的方法。并在最后給出了提高帶接觸摩擦仿真計算收斂性的原則和方法。

[1] ADAMS Optimization Guide [M]. MSC. Software Corporation, 1994(11).

[2] H.M. Lankarani and P.E. Nikravesh, “Continuous Contact Force Models for Impact Analysis in Multibody Systems”, Nonlinear Dynamics, 5: 193-207, 1994.

[3] Johnson, K. L.: Contact mechanics. Cambridge University Press, 6. Nachdruck der 1. Auflage, 2001.

[4] Popov, Valentin L.: Contact Mechanics and Friction. Physical Principles and Applications, Springer-Verlag, 2010, 362 p., ISBN 978-3-642-10802-0.

Research of Mechanism of ADAMS Contact and Contact Friction

Gao Yijia
( Shaanxi Ploy Special Vehicle Manufacturing Co., LTD., Shaanxi Xi'an 710200 )

Firstly This paper describes mechanism of action of ADAMS contact and contact friction in detail. Then a dynamics analysis of a hooklift garbage truck is used to illustrate how to set the contact parameters. Each parameter is changed individually, and variation of contact force is investigated. Finally tips for working with contacts is given.

ADAMS; contact; contact friction

U461

A

1671-7988 (2017)06-64-04

高一佳（1982-），男，工程師，就職于陜西保利特種車制造有限公司，主要從事軍用裝甲車輛結構強度及多體動力學仿真分析工作。也是高度不連續的。摩擦是路徑相關的(能量耗散)現象，需要用較小的時間步和精確的加載歷程。摩擦響應可能是混沌的，使求解收斂困難，特別是在低速時。摩擦困難經常會和已經存在的接觸求值困難相復合，使摩擦問題的收斂性變得更困難。

10.16638/j.cnki.1671-7988.2017.06.020