前面幾節(jié)介紹了接觸分析中一些需要注意的重要問題,本節(jié)將給出兩個接觸分析模型的實例,請讀者自己判斷模型中是否有錯誤,如果有錯誤,應如何改正。隨書光盤中給出了錯誤的和正確的模型文件。
16.3.1 接觸分析綜合實例1
【綜合實例 16-1】
請打開隨書光盤以下文件夾中的模型文件:contact-example-1 wrongContact-Exam-ple-1-Wrong. cae
這是一個圓柱形剛體下壓坯料的平面應變模型,根據對稱性,只取右半部分建模(如圖16-19所示),分析的目的是獲得接觸區(qū)域的應力分布。請讀者檢查此模型是否存在錯誤。
此模型中已經定義了接觸分析模型的所有必備參數,看上去似乎一切正常,提交分析作業(yè)后,在ABAQUS/CAE中沒有出現任何錯誤信息或警告信息,分析可以順利完成,在后處理時也可以看到很漂亮的云紋圖。但是正如本書導言中所強調過的那樣,能夠建立模型、完成分析和得到計算結果,并不是就意味著萬事大吉了。如果建模方法或模型參數不正確,得到的結果就會是不準確的,有時甚至可能是完全錯誤的。
仔細檢查此模型,可以發(fā)現存在以下問題:
1)接觸定義中的位置誤差限度(ADJUST參數)過大。在Visualization功能模塊中顯示未變形圖,可以看到在初始狀態(tài)下,壞料頂部有一個節(jié)點的坐標發(fā)生了變化(如圖16-20所示)。這是因為在定義接觸時,將位置誤差限度設為了0.01,這個值太大(如圖16-21所示),其含義為,如果某個從面節(jié)點和主面的初始距離小于0.01就改變此節(jié)點的坐標,使其恰好位于主面上。
由于模型中的節(jié)點坐標可能存在數值誤差,使從面節(jié)點和主面之間出現一個本不應有的微小距離,因此一般情況下應設置一定的位置誤差限度,以保證從面節(jié)點和主面在初始狀態(tài)下就能夠相互接觸。但本實例中的位置誤差限度值設置得太大,造成了不應出現的節(jié)點坐標變化。應該將圖16-21中的位置誤差限度值改為小一些的值,例如0.00001。
圖16-20未變形圖:在初始狀態(tài)下,坯料頂部有一個節(jié)點的坐標發(fā)生了變化
圖16-21 接觸定義中的位置誤差限度設為0.01
2)網格密度。初學者最容易犯的錯誤就是,不假思索地為整個模型劃分均勻的網格。本實例的分析目的是得到接觸區(qū)域的應力分布,如果在接觸區(qū)域劃分如圖16-19所示的粗糙網格,得到的應力結果精度會非常差。正確的做法是,首先為接觸部位分割(partition)出形狀規(guī)則的區(qū)域,再分別為不同區(qū)域定義不同密度的網格種子(如圖16-22所示),從而生成局部細化的高質量網格(如圖16-23所示)。
圖16-22為接觸部位分割出形狀規(guī)則的區(qū)域,再分別為不同區(qū)域定義不同密度的網格種子
圖16-23局部細化的高質量網格
3)單元類型。在Visualization功能模塊中顯示變形圖,可以看到單元變?yōu)榻惶娉霈F的梯形形狀(如圖16-24所示),這說明模型中的減縮積分單元CPE4K出現了沙漏模式的數值問題,本書第9.2.1節(jié)“沙漏模式”詳細介紹了相關知識。
選擇菜單Result→HistoryOuiput,繪制偽應變能ALLAE和內能ALLIE的曲線,可以看到在分析結束時,偽應變能為60921.6,內能為1.24224c7,偽應變能為內能的0.5%。盡管這個比例小于1%,但單元已經出現了交替出現的梯形形狀,這仍然是一種不正常的位移結果正確的做法是,將單元類型改為非協(xié)調模式(incompatiblemodes)單元CPE4I。使用非協(xié)調模式單元時,應盡量保證關鍵區(qū)域的單元形狀是矩形的,否則會影響結果的計算精度。
4)對稱邊界條件。坯料左側沒有定義邊界條件,既然是根據對稱性只取一半建模,就必須在對稱面上定義對稱邊界條件XSYMM。
那么如果此模型不是1/2坯料,圖16-19中的左側邊就是坯料本身的邊界,是否就沒有問題了呢?這時模型就存在另外一個問題,坯料左上角的節(jié)點是自由的,在剛體的擠壓下會移動到剛體表面的上方,出現穿透現象。產生這種問題的原因是,模型違反了一個接觸分析的重要原則:從面節(jié)點在分析過程中應該始終處于主面法線的覆蓋范圍之內,而不能落到主面之外。
因此,如果不是對稱模型,就必須把剛體圓弧的左端畫得更長一些,以保證坯料左上角的節(jié)點始終處于主面的覆蓋范圍之內。
5)剛體參考點的位置。模型中的剛體參考點沒有位于剛體圓弧的圓心處。此模型中的剛體只發(fā)生平動,因此剛體參考點的位置不會影響分析結果。假如剛體會發(fā)生轉動,剛體參考點就必須位于轉動的軸心處。因此建議讀者在建模時養(yǎng)成習慣,讓剛體參考點位于剛體形心的位置,或載荷合力的作用點上。
題:
隨書光盤的以下文件夾中提供了正確的模型文件,其中更正了上述5個方面的問
contact-example-1correctContact-Example-l-Correct. cae
(內容、圖片來源:《ABAQUS有限元分析常見問題解答》,侵刪)
版權與免責聲明:
凡未注明作者、來源的內容均為轉載稿,如出現版權問題,請及時聯(lián)系我們處理。我們對頁面中展示內容的真實性、準確性和合法性均不承擔任何法律責任。如內容信息對您產生影響,請及時聯(lián)系我們修改或刪除。
[行業(yè)資訊] 達索系統(tǒng) SIMULIA Fe-safe技術論壇即將開講
2025-08-01
[ABAQUS] 哪些行業(yè)會用到碰撞安全仿真?Abaqus軟件在其中有何應用
2025-08-01
[ABAQUS] Abaqus在疲勞耐久分析中的優(yōu)勢有哪些
2025-08-01
[ABAQUS] Abaqus結構優(yōu)化仿真如何實現產品性能與經濟性的平衡?
2025-08-01
2025-07-30
[CST] SIMULIA電磁仿真線上技術匯|專題三:CST低頻仿真高
2025-07-30
[行業(yè)資訊] 基于組合優(yōu)化策略的白車身輕量化設計
2025-07-29
[行業(yè)資訊] 復合材料單螺栓連接結構預緊力的影響因素研究
2025-07-29
[行業(yè)資訊] 車輪損傷狀態(tài)下高速列車平穩(wěn)性和舒適度分析
2025-07-29
[有限元知識] SIMULIA 一文讀懂系列:高精度仿真技術破解風電噪聲難
2025-07-28
2023-08-29
2023-08-24
[ABAQUS] ABAQUS如何建模?ABAQUS有限元分析教程
2023-07-07
[ABAQUS] 有限元分析軟件abaqus單位在哪設置?【操作教程】
2023-09-05
[ABAQUS] ABAQUS軟件中Interaction功能模塊中的綁定接
2023-07-19
[ABAQUS] Abaqus怎么撤回上一步操作?Abauqs教程
2024-05-01
[ABAQUS] abaqus里面s11、s12和u1、u2是什么意思?s和
2023-08-30
[ABAQUS] Abaqus單位對應關系及參數介紹-Abaqus軟件
2023-11-20
[ABAQUS] ABAQUS軟件教程|場變量輸出歷史變量輸出
2023-07-18
[ABAQUS] Abaqus中的S、U、V、E、CF分別是什么意思?
2024-05-11
[行業(yè)資訊] 達索系統(tǒng) SIMULIA Fe-safe技術論壇即將開講
2025-08-01
[行業(yè)資訊] 基于組合優(yōu)化策略的白車身輕量化設計
2025-07-29
[行業(yè)資訊] 復合材料單螺栓連接結構預緊力的影響因素研究
2025-07-29
[行業(yè)資訊] 車輪損傷狀態(tài)下高速列車平穩(wěn)性和舒適度分析
2025-07-29
[有限元知識] SIMULIA 一文讀懂系列:高精度仿真技術破解風電噪聲難
2025-07-28
2025-07-28
[行業(yè)資訊] 制造企業(yè)是如何發(fā)掘仿真軟件價值的
2025-07-25
[行業(yè)資訊] 達索有限元仿真軟件Abaqus一套多少錢?
2025-07-24
[行業(yè)資訊] 通過XFlow和Simpack動力集中動車組單線隧道內列尾
2025-07-23
[行業(yè)資訊] SIMULIA的可靠航空通信與導航系統(tǒng)解決方案
2025-07-23