



下載本文檔
版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、接觸問題之經驗談關鍵字:法向關系、切向關系、罰函數、拉格朗日乘子法、擴展拉格朗日乘子法。接觸問題的關鍵在于接觸體間的相互關系,也就是求解對象之間的相互聯系,此關系又可分為在接觸前后的法向關系與切向關系。法向關系(垂直于接觸面的接觸關系):在法向關系中,必須實現兩點:1)接觸力的傳遞;2)兩接觸面間沒有穿透。也就是說即有力的傳遞,同時又不能發生穿透現象。ANSYS通過兩種算法來實現此法向接觸關系:罰函數法和拉格朗日乘子法。1、 罰函數法是通過接觸剛度在接觸力與接觸面間的穿透值(接觸位移)間建立力與位移的線性關系:接觸剛度×接觸位移=法向接觸力對面面接觸單元(CONTA17*系列),接觸
2、剛度由實常數FKN來定義,穿透值在程序中通過分離的接觸體上節點間的距離來計算。接觸剛度越大,則穿透就越小,理論上在接觸剛度為無窮大時,可以實現完全的接觸狀態,即使穿透值等于零。但是顯而易見,在程序計算中,接觸剛度不可能為無窮大(否則會導致病態矩陣),穿透也就不可能真實達到零,而只能是個接近于零的有限值。以上力與位移的接觸關系可以很容易地合并入整個結構的平衡方程組【K】×【X】=【F】中去,并不改變總體剛度【K】的大小(【X】位移矩陣,【F】力矩陣);同時,采用罰函數法有以下幾個問題必須解決:1)接觸剛度FKN應該取多大?2)接觸剛度FKN取大些可以減少虛假穿透,但是也可能會使剛度矩陣
3、成為病態。3)既然與實際情況不符合的虛假穿透既然是不可避免的,那么可以允許有多大為合適?因此,在ANSYS程序里,通常輸入FKN實常數不是直接定義接觸剛度的數值,而是定義接觸面覆蓋之下的單元剛度的一個因子,這使得用戶可以方便地定義接觸剛度,一般FKN取0.1到1中間的值。當然,在需要時,也可以把接觸剛度直接定義,例如:FKN輸入為負數,則程序將其值理解為直接輸入的接觸剛度值。對于接近病態的剛度陣,不要使用迭代求解器,例如PCG等,它們會需要更多的迭代次數,同時有可能導致不收斂。這時可以使用直接法求解器,例如稀疏矩陣求解器等,這些求解器可以有效求解病態問題。穿透的大小影響結果的精度:用戶可以用P
4、LESOL,CONT,PENE來在后處理中查看穿透的數值大小。如果使用的是罰函數法求解接觸問題,用戶一般需要試用多個FKN值進行計算,可以先用一個較小的FKN值開始計算,例如0.1。因為較小的FKN有助于收斂,然后再逐步增加FKN值進行一系列計算,最后得到一個滿意的穿透值。FKN的收斂性要求和穿透太大產生的計算誤差總會是一對矛盾。解決此矛盾的辦法是在接觸算法中采用擴展拉格朗日乘子法。此方法在接觸問題的求解控制中可以有更多更靈活的控制,可以更快的實現一個需要的穿透極限。2、拉格朗日乘子法與擴展拉格朗日乘子法拉格朗日乘子法與罰函數法不同,不是采用力與位移的關系來求接觸力,而是把接觸力作為一個獨立自
5、由度。因此這里不需要進行迭代,而是在方程里直接求出接觸力(接觸壓力)來。Kx=F+Fcontact從而,拉格朗日乘子法不需要定義人為的接觸剛度去滿足接觸面間不可穿透的條件,可以直接實現穿透為零的真實接觸條件,這是罰函數法所不可能實現的,也是采用拉格朗日乘子法的優點所在,但是在使用拉格朗日乘子法有下列注意事項:1)若剛度矩陣中有零對角元,則使得某些求解器不可使用,而只能使用直接法求解器,例如波前法或系數求解器,而PCG之類迭代求解器是不能用于有零主元問題的。2)由于增加了額外的自由度,剛度矩陣變大了。3)對于一個可能發生的接觸突變的解除類問題,就是在接觸狀態發生變化時,例如從接觸到分離,從分離到
6、接觸,此時接觸力有個突變,產生chattering(接觸狀態的振動式交替改變)。如何控制這種chattering,是純粹拉格朗日法所難以解決的。因此,為控制chattering,ANSYS采用的是罰函數法與拉格朗日法相結合的擴展拉格朗日乘子法。在擴展拉格朗日法中,可以采用實常數TOLN來控制最大允許穿透值,還有最大允許拉力FTOL;這兩個參數只對擴展拉格朗日乘子法有效。在擴展拉格朗日乘子法里,程序按照罰函數法開始,與純粹拉格朗日法類似,用TOLN來控制最大允許穿透值。如果迭代中發現穿透大于允許的TOLN值(對CONTA178單元是TOLN,而對面面接觸單元CONTA171-174則是FTOLN),則將各個接觸單元的接觸剛度加上接觸力乘以拉格朗日乘子的數值。因此,這種擴展拉格朗日法是不停更新接觸剛度的罰函數法,這種更新不斷重復,直到計算的穿透值小于允許值為止。盡管與拉格朗日法相比,擴展拉格朗日法的穿透并不是零,與罰函數法相比,可能迭帶次數會更多;但是擴展拉格朗日法有下列優點:1)較少病態矩陣,對于接觸單元的接觸剛度取值可能更合理。2)與罰函數法相比較少病態,與單純的拉格朗日法相比,又沒有剛度陣零對角元,因此在選擇求解器上沒有限制,PCG等迭代求解器都可以應用。3)用戶可以自由控制允許的穿透值TOLN(如果輸入了TOLN,而使用罰函數法,則程序
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 信陽師范大學《液壓與氣壓傳動1》2023-2024學年第二學期期末試卷
- 煙臺汽車工程職業學院《波斯語報刊選讀》2023-2024學年第二學期期末試卷
- 江西工業貿易職業技術學院《中醫眼科學》2023-2024學年第一學期期末試卷
- 四川外國語大學成都學院《ERP供應鏈管理》2023-2024學年第二學期期末試卷
- 江蘇省海安市2025屆高三下第一次階段性檢測試題生物試題含解析
- 江西應用科技學院《PROE三維機械設計》2023-2024學年第二學期期末試卷
- 二零二五版美團會員服務協議
- 二零二五經營場地租賃協議書范例
- 二零二五版投資理財協議
- 二零二五版投資人入股協議書
- 項目建筑智能化工程施工招標文件模板
- 110kv線路施工方案
- 大東鞋業合同協議書
- 港口與船舶協同運營優化
- 蘇州工業園區應急管理系統招聘筆試真題2023
- 用所給詞的適當形式填空(專項訓練)人教PEP版英語六年級上冊
- NB-T42090-2016電化學儲能電站監控系統技術規范
- 勉縣房地產市場調研報告
- NBT 47013.2-2015 承壓設備無損檢測 第2部分:射線檢測
- 《大學生美育》 課件 第七章 藝術美
- 2024年巴西獸醫快速診斷市場機會及渠道調研報告
評論
0/150
提交評論