ANSYS接觸問題的計算方法及參數(shù)設置

1、 ANSYS接觸問題的計算方法及參數(shù)設置 接觸問題的關鍵在于接觸體間的相互關系，此關系又可分為在接觸前后的法向關系與切向關系。 法向關系： 在法向，必須實現(xiàn)兩點：1）接觸力的傳遞。2）兩接觸面間沒有穿透。 ANSYS通過兩種算法來實現(xiàn)此法向接觸關系：罰函數(shù)法和拉格朗日乘子法。 1罰函數(shù)法 是通過接觸剛度在接觸力與接觸面間的穿透值（接觸位移）間建立力與位移的線性關系： 接觸剛度*接觸位移=法向接觸力 對面面接觸單元17*，接觸剛度由實常數(shù)FKN來定義。 穿透值在程序中通過分離的接觸體上節(jié)點間的距離來計算。接觸剛度越大，則穿透就越小，理論上在接觸剛度為無窮大時，可以實現(xiàn)完全的接觸狀態(tài)，使穿透值等于

2、零。但是顯而易見，在程序計算中，接觸剛度不可能為無窮大（否則病態(tài)），穿透也就不可能真實達到零，而只能是個接近于零的有限值。 以上力與位移的接觸關系可以很容易地合并入整個結構的平衡方程組K*X=F中去。并不改變總剛K的大小。這種罰函數(shù)法有以下幾個問題必須解決： 1）接觸剛度FKN應該取多大？ 2）接觸剛度FKN取大些可以減少虛假穿透，但是會使剛度矩陣成為病態(tài)。 3）既然與實際情況不符合的虛假穿透既然是不可避免的，那么可以允許有多大為合適？ 因此，在ANSYS程序里，通常輸入FKN實常數(shù)不是直接定義接觸剛度的數(shù)值，而是接觸體下單元剛度的一個因子，這使得用戶可以方便地定義接觸剛度了，一般FKN取0.

3、1到1中間的值。當然，在需要時，也可以把接觸剛度直接定義，F(xiàn)KN輸入為負數(shù)，則程序將其值理解為直接輸入的接觸剛度值。 對于接近病態(tài)的剛度陣，不要使用迭代求解器，例如PCG等。它們會需要更多的迭代次數(shù)，并有可能不收斂?？梢允褂弥苯臃ㄇ蠼馄鳎缦∈枨蠼馄鞯?。這些求解器可以有效求解病態(tài)問題。 穿透的大小影響結果的精度。用戶可以用PLESOL,CONT,PENE來在后處理中查看穿透的數(shù)值大小。如果使用的是罰函數(shù)法求解接觸問題，用戶一般需要試用多個FKN值進行計算，可以先用一個較小的FKN值開始計算，例如0.1。因為較小的FKN有助于收斂，然后再逐步增加FKN值進行一系列計算，最后得到一個滿意的穿透值

4、。 FKN的收斂性要求和穿透太大產生的計算誤差總會是一對矛盾。解決此矛盾的辦法是在接觸算法中采用擴展拉格朗日乘子法。此方法在接觸問題的求解控制中可以有更多更靈活的控制?？梢愿斓膶崿F(xiàn)一個需要的穿透極限。 2拉格朗日乘子法與擴展拉格朗日乘子法 拉格朗日乘子法與罰函數(shù)法不同，不是采用力與位移的關系來求接觸力，而是把接觸力作為一個獨立自由度。因此這里不需要進行迭代，而是在方程里直接求出接觸力（接觸壓力）來。 Kx=F+Fcontact 從而，拉格朗日乘子法不需要定義人為的接觸剛度去滿足接觸面間不可穿透的條件，可以直接實現(xiàn)穿透為零的真實接觸條件，這是罰函數(shù)法所不可能實現(xiàn)的。使用拉格朗日乘子法有下列注意

5、事項： 1）剛度矩陣中將有零對角元，使有些求解器不克使用。只能使用直接法求解器，例如波前法或系數(shù)求解器。而PCG之類迭代求解器是不能用于有零主元問題的。 2）由于增加了額外的自由度，剛度陣變大了。 3）一個可能發(fā)生的嚴重問題，就是在接觸狀態(tài)發(fā)生變化時，例如從接觸到分離，從分離到接觸，此時接觸力有個突變，產生chattering（接觸狀態(tài)的振動式交替改變）。如何控制這種chattering，是純粹拉格朗日法所難以解決的。 因此，為控制chattering，ANSYS采用的是罰函數(shù)法與拉格朗日法混合的擴展拉格朗日乘子法。在擴展拉格朗日法中，可以采用實常數(shù)TOLN來控制最大允許穿透值。還有最大允許拉

6、力FTOL。這兩個參數(shù)只對擴展拉格朗日乘子法有效。 在擴展拉格朗日乘子法里，程序按照罰函數(shù)法開始，與純粹拉格朗日法類似，用TOLN來控制最大允許穿透值。如果迭代中發(fā)現(xiàn)穿透大于允許的TOLN值，（對178單元是TOLN，而對面面接觸單元171-174則是FTOLN）則將各個接觸單元的接觸剛度加上接觸力乘以拉格朗日乘子的數(shù)值。因此，這種擴展拉格朗日法是不停更新接觸剛度的罰函數(shù)法，這種更新不斷重復，直到計算的穿透值小于允許值為止。 盡管與拉格朗日法相比，擴展拉格朗日法的穿透并不是零，與罰函數(shù)法相比，可能迭帶次數(shù)會更多。擴展拉格朗日法有下列優(yōu)點： 1）較少病態(tài)，個接觸單元的接觸剛度取值可能更合理。 2

7、）與罰函數(shù)法相比較少病態(tài)，與單純的拉格朗日法相比，沒有剛度陣零對角元。因此在選擇求解器上沒有限制，PCG等迭代求解器都可以應用。 3）用戶可以自由控制允許的穿透值TOLN。（如果輸入了TOLN，而使用罰函數(shù)法，則程序忽略它） 切向關系： 摩擦的處理與法向接觸力類似。由于摩擦是非對稱的，使問題變的更為復雜。ANSYS缺省是做對稱求解，即使用對稱求解器作近似求解。但是可以改變幾個選項強迫做非對稱求解。非對稱求解更精確，但是計算量大許多。 參數(shù)設置 實常數(shù)FKN，F(xiàn)TOLN，ICONT，PINB，PMAX，和PMIN，你既可以定義一個正值也可以定義一個負值。既然程序將正值作為比例因子，將負值作為真實

8、值，程序將下面覆蓋原單元的厚度作為ICON，F(xiàn)TOLN，PINB，PMAX和PMIN的參考值。 那么：單位改變后，你只要保證你的FKN，F(xiàn)TOL的物理值大小不變就可以了。 作為非線性的接觸問題，F(xiàn)KN與FTOLN取值沒有公式可套，只能慢慢試出個合適的值來。 由于FKN越小越容易收斂，你可以先用一個比較小的FKN開始計算，收斂后再改大些，由于接觸剛度FKN越大則穿透越小，結果越合理，因此在適當時應該把FKN改大些重新計算直到一個滿意的結果。（可以用PLESOL,CONT,PENE來在后處理中查看穿透的數(shù)值大?。?注意：FKN與穿透值大小總是一對矛盾，在矛盾中達到合理的平衡是做接觸分析的關鍵。 理論上接觸問題的穿透應該是零，所以穿透越小則越精確，也即FKN越大則越精確，但是實際上不可能使穿透變?yōu)榱?，F(xiàn)KN太大了計算就會不收斂。 如果只改FKN得到的穿透值總是太大，不能滿意，再加大FKN就不收斂了。那么可以把求解方法改為擴展拉格朗日乘子法。（用KEYOPT改），再使用F