ANSYS接觸分析(中文)

1、精品文檔接觸問題（參考ansys勺中文幫助文件）當兩個分離的表面互相碰觸并共切時， 就稱它們牌接觸狀態(tài)。 在一般的物理意義中，牌接觸狀態(tài)的表面有下列特點：1、 不互相滲透；2、 能夠互相傳遞法向壓力和切向摩擦力；3、 通常不傳遞法向拉力。接觸分類：剛性體柔性體、柔性體柔性體實際接觸體相互不穿透， 因此， 程序必須在這兩個面間建立一種關(guān)系， 防止它們在有限元分析中相互穿過。罰函數(shù)法。接觸剛度 lagrange 乘子法，增加一個附加自由度（接觸壓力） ，來滿足不穿透條件將罰函數(shù)法和 lagrange 乘子法結(jié)合起來，稱之為增廣 lagrange 法。三種接觸單元：節(jié)點對節(jié)點、節(jié)點對面、面對面。接觸

2、單元的實常數(shù)和單元選項設(shè)置：fkn法向接觸剛度。這個值應該足夠大，使接觸穿透量??；同時也應該足夠小，使問題沒有病態(tài)矩陣。fkn值通常在0.110之間，對于體積變形問題，用值1.0（默認） ，對彎曲問題，用值0.1 。ftoln最大穿透容差。穿透超過此值將嘗試新的迭代。這是一個與接觸單元下面的實體單元深度（h）相乘的比例系數(shù)，缺省為 0.1。此值太小，會引起收斂 困難。icont初始接觸調(diào)整帶。它能用于圍繞目標面給出一個“調(diào)整帶”，調(diào)整帶內(nèi)任 何接觸點都被移到目標面上；如果不給出icontfi, ansys艮據(jù)模型的大小提供 一個較小的默認值（ 0.03 =pinb：指定近區(qū)域接觸范圍（球形區(qū)）

3、。當目標單元進入pinball區(qū)時，認為它處于近區(qū)域接觸， pinball 區(qū)是圍繞接觸單元接觸檢測點的圓（二維）或球（三維）?？梢杂脤嵆?shù)pinb調(diào)整球形區(qū)（此方法用于初始穿透大的問題是必要的） pmin和pmax初始容許穿透容差。這兩個參數(shù)指定初始穿透范圍，ansysf巴整個目標面（連同變形體）移到到由pmin和pmax旨定的穿透范圍內(nèi)，而使其成為 閉合接觸的初始狀態(tài)。初始調(diào)整是一個迭代過程，ansysr多使用20個迭代步把目標面調(diào)整到pmin和pmam圍內(nèi)，如果無法完成，給出警告，可能需要修改 幾何模型。taumax接觸面的最大等效剪應力。給出這個參數(shù)在于，不管接觸壓力值多大， 只要等效

4、剪應力達到最大值 taumax就會發(fā)生滑動。該剪應力極限值通常用于 接觸壓力會變得非常大的情況。cnof指定接觸面偏移。+cnof1力口過盈、-cnofs少過盈或產(chǎn)生間隙、cno前 與幾何穿透組合應用。fkop接觸張開彈簧剛度。針對不分離或綁定接觸模型，需要設(shè)置實常數(shù) fkop 該常數(shù)為張開接觸提供了一個剛度值。fkofffi止接觸面的分離；fkopr認為1.0, 用于建立粘結(jié)模型，用一個較小值（ 1e-5 ）去建立軟彈簧模型。fkt:切向接觸剛度。作為初值，可以采用 fkt= 0.01*fkn,這是大多數(shù)ansys 接觸單元的缺省值。cohe粘滯力。即沒有法向壓力時開始滑動的摩擦應力值。fa

5、ct dc定義摩擦系數(shù)變化規(guī)律mu= muk*（1+ （fact-1） exp（ dc*vtfs/deltat ）式中：mui4動摩擦系數(shù)（用戶自己定義）fact= mus/muk用戶自己定義）mu容靜摩擦系數(shù)do衰減系數(shù)（用戶自己定義）vtfs/deltat =表面間的相對速度注意：動摩擦系數(shù)由被指定為材料屬性（mu,由m騎令或gui定義。缺省值:fact= 1, mu容 muk0, dc= 0keyopt 的介紹，以 target170,conta173 為例：首先介紹 170 的 keyoptkeyop（t 3） ：定義接觸行為keyop（t 1） ：單元階數(shù)（是否含有中節(jié)點）keyop

6、t1） =0:低階單元（不含中節(jié)點）keyopt1） =1：高階單元（含中節(jié)點）keyop（t 2） ：剛體目標面約束條件=0時，自動約束選項，每一個載荷步的末尾，程序內(nèi)部將是性面重新 設(shè)置約束。滿足以下條件，剛性面則缺省為自動約束沒有明確定義邊界條件；目標面與其它單元沒有聯(lián)系； 沒有定義耦合或約束方程。=1時，用戶定義選項。conta173 的 keyopt ：keyopt1）:自由度選項。=0時，結(jié)構(gòu)：ux, uy和uz; =1時，結(jié)構(gòu)和 熱；=2時，temp（用于純熱接觸問題）keyopt2）:選擇接觸算法。=0時，增廣的拉格朗日法（缺省選項）,推 薦于一般應用，它對罰剛度不太敏感，但是

7、也要求給出一個穿透容差。=1時，罰函數(shù)法。 它推薦應用于單元非常扭曲、 大摩擦系數(shù)和用增廣的拉格朗日法收斂 行為不好的問題。keyopt4）:選擇接觸檢查點。=0時，高斯點（缺省選項，推薦）；=1 時，節(jié)點；ansyso寸面單元默認用高斯積分點作為接觸檢查點。keyopt5）:自動cnof3整。允許ansyss于初始狀態(tài)自動給定 cnofs 導致“剛好接觸”配置。=0時，不進行自動調(diào)整；=1時，閉合間隙；=2時, 減小穿透；=3時，閉合間隙/減小穿透。keyop（t7） ：時間步控制選項。（只有在 solution control 中打開基于接觸狀 態(tài) 變化 的 時 間 步 預 測 ， 此 選

8、 項 才 起作 用 。 solution>unabridged menu>load step opts>solution ctrl=0時，不控制，不影響自動時間步長。對靜力問題自動時間步打開時此選項 一般是足夠的。=1時，自動二分，如果接觸狀態(tài)變化明顯，時間步長將二分，對于動力問 題自動二分通常是足夠的；=2時，合理值。比自動細分更耗時的算法；=3時，最小值。此選項為下一子步預測最小時間增量 （很耗機時，不推薦）。 keyopt8）:防止偽接觸選項。=0時，不防止；=1時，檢測并忽略偽接觸。keyop t9）:初始穿透間隙控制。=0時，包括幾何穿透/間隙和cnof =1時,

9、忽略幾何穿透/間隙和cnof =2時，包括幾何穿透/間隙和cnof且在第一個 載荷步中漸變；=3時，忽略幾何穿透/間隙，包括cnof =4時，忽略幾何穿 透/間隙，包括cnof且在第一個載荷步中漸變。keyopt10）:接觸剛度更新控制。=0,閉合狀態(tài)的接觸剛度不進行任何更新； =1,每一載荷步更新閉合狀態(tài)的接觸剛度（fknb fkt,由用戶指定）；=2,與 =1同，止匕外，在每一子步，程序自動更新接觸剛度（根據(jù)變形后下伏單元的剛 度） 。keyop（t 11） ：殼、梁單元厚度影響。如果已經(jīng)創(chuàng)建了一個梁或殼單元模型， 接觸表面能夠偏置， 用于考慮梁或殼的厚度。=0,在中面接觸（默認）；=1,

10、在指定表面的頂部或底部。注意：當用shell18仰元時，由于大應變變形引起的厚度改變也被考慮。keyop（t 12） ：創(chuàng)立不同的接觸表面相互作用模型。=0,標準的接觸行為，張開時法向壓力為 0;=1,粗糙接觸行為，不發(fā)生滑動（類似無限摩擦系數(shù)）；=2,不分離，允許滑動；=3,綁定接觸，目標面和接觸面一旦接觸就粘在一起；=4,不分離接觸（總是），初始位于pinball區(qū)域內(nèi)或已經(jīng)接觸的接觸檢查 點在法向不分離；=5,綁定接觸（總是），初始位于pinball區(qū)域內(nèi)或已經(jīng)接觸的接觸檢查點 總是與目標面綁定在一起；=6,綁定接觸（初始接觸），只在初始接觸的地方采用綁定，初始張開的地 方保持張開。3歡

11、迎下載 。5.1 概述接觸問題是一種高度非線性行為，需要較多的計算機資源。為了進行切實有 效的計算，理解問題的物理特性和建立合理的模型是很重要的。接觸問題存在兩個較大的難點：其一，在用戶求解問題之前，用戶通常不知道接觸區(qū)域。隨載荷、材料、邊界條件和其它因素的不同，表面之間可以接觸或者分開，這往往在很大程度上是難以預料的，并且還可能是突然變化的。其二，大多數(shù)的接觸問題需要考慮摩擦作用， 有幾種摩擦定律和模型可供挑選， 它們都是非線性的。摩擦效應可能是無序的，所以摩擦使問題的收斂性成為一個難點。注意 - 如果在模型中，不考慮摩擦，且物體之間的總是保持接觸，則可以應用約束方程或自由度藕合來代替接觸。

12、約束方程僅在小應變分析（ nlgeom，off）中可用。見ansy制odeling and meshing guide中的 12, coupling and constraint equations 。除了上面兩個難點外，許多接觸問題還必須涉及到多物理場影響，如接觸區(qū) 域的熱傳導、電流等。5.1.1 顯式動態(tài)接觸分析能力除了本章討論的隱式接觸分析外，ansy變在ansys/ls-dyna提供了顯式接觸分析功能。顯式接觸分析對于短時間接觸-碰撞問題比較理想。關(guān)于ansys/ls-dyna更多的信息參見 ansys/ls-dyna users guide。5.2 一般接觸分類接觸問題分為兩種基本類

13、型：剛體一柔體的接觸，柔體一柔體的接觸。在剛體一柔體的接觸問題中，接觸面的一個或多個被當作剛體，（與它接觸的變形體相比，有大得多的剛度） 。一般情況下，一種軟材料和一種硬材料接觸時，可以假定為剛體接柔體的接觸， 許多金屬成形問題歸為此類接觸。 柔體接柔體的接觸是一種更普遍的類型，在這種情況下，兩個接觸體都是變形體（ 有相似的剛度）柔體接柔體接觸的一個例子是栓接法蘭。精品文檔5.3 ansys接觸分析功能ansysfc持三種接觸方式：點一點，點一面，面一面接觸。每種接觸方式使 用不同的接觸單元集，并適用于某一特定類型的問題。為了給接觸問題建模，首先必須認識到模型中的哪些部分可能會相互接觸。 如果

14、相互作用的其中之一是一點，模型的對應組元是一個節(jié)點。如果相互作用的 其中之一是一個面，模型的對應組元是單元，如梁單元、殼單元或?qū)嶓w單元。有 限元模型通過指定的接觸單元來識別可能的接觸對，接觸單元是覆蓋在分析模型 接觸面之上的一層單元，至于ansy皺用的接觸單元和使用它們的過程，后面會分類詳述，然后論述ansysg觸單元和他們白功能。參見ansys elements reference 和 ansys theory reference 。點-百 八、八、點-面面-面conta 12conta52cont178aonta26conta48conta49contac171,1：target 169c

15、onta174 targe170點-百 八、八、yyy點-面yyy面-面yyyy2-dyyyyyy3-dyyyy滑動小小小大大大大大i曲面yy圓柱間隙yy純lagrange乘子y增加lagrange乘子yyyyy接觸剛度用戶定 義半自 動用戶定 義用戶定 義用戶定 義半自動半自動自動網(wǎng)格 工具eintf日ntfnonegcgengcgen esurfesur低階yyyyyyyy高階yyy表5-1 ansys接觸分析功能ct剛體-柔 體yyyyyyyy柔體-柔 體yyyyyyy熱接觸yyyy5.3.1 面一面的接觸單元ansysaf剛體一柔體和柔體一柔體的面一面的接觸單元。 這些單元應用“目 標

16、”面和“接觸”面來形成接觸對。分另i用targe16或targe17來模擬2d和3d目標面。用 conta1z1 conta1z2conta1z3conta17來模擬接觸面。為了建立一個“接觸對”，給目標單元和接觸單元指定相同的實常數(shù)號。參 見 5.4。這些面-面接觸單元非常適合于過盈裝配安裝接觸或嵌入接觸，鍛造，深拉 問題。與點一面接觸單元相比，面一面接觸單元有許多優(yōu)點：支持面上的低階和高階單元（即角節(jié)點或有中節(jié)點的單元）；支持有大滑動和摩擦的大變形。計算一致剛度陣，可用不對稱剛度陣 選項；提供為工程目的需要的更好的接觸結(jié)果，如法向壓力和摩擦應力；沒有剛體表面形狀的限制，剛體表面的光滑性不是

17、必須的，允許有自 然的或網(wǎng)格離散引起的表面不連續(xù)；與點一面接觸單元比，需要較少的接觸單元，因而只需較小的磁盤空 間和cpu時間，并具有高效的可視化；允許多種建?？刂?，例如：綁定接觸，不分離接觸，粗糙接觸；漸變初始穿透；目標面自動移動到初始接觸；平移接觸面（考慮梁和單元的厚度），用戶定義的接觸偏移；5歡在下載精品文檔支持熱 - 力耦合分析。使用這些單元來做為剛性目標面，能模擬2d和3d中的直線（面）和曲線（面）， 通常用簡單的幾何形狀例如圓、拋物線、球、圓錐、圓柱來模擬曲面。更復雜的剛體形狀或普通可變形體，可以應用特殊的前處理技巧來建模，參見 5.4 。面- 面接觸單元不能很好地應用于點 - 點

18、或點 - 面接觸問題，如管道或鉚頭裝配。 在這種情況下， 應當應用點 - 點或點 - 面接觸單元。 用戶也可以在大多數(shù)接觸區(qū)域應用面- 面接觸單元，而在少數(shù)接觸角點應用點 - 點接觸單元。面面接觸單元只支持一般的靜態(tài)或瞬態(tài)分析，屈曲、模態(tài)、譜分析或子結(jié)構(gòu)分析。 不支持諧響應分析、 縮減或模態(tài)疊加瞬態(tài)分析， 或縮減或模態(tài) 疊加諧響應分析。本章后面將分別討論ansy加同接觸分析類型的能力。5.3.2 點一面接觸單元接一面接觸單元主要用于給點一面接觸行為建模，例如兩根梁的相互接觸（ 梁端或尖角節(jié)點 ） ，鉚頭裝配部件的角點。如果通過一組節(jié)點來定義接觸面，生成多個單元，那么可以通過點一面接觸 單元來模

19、擬面一面的接觸問題。面既可以是剛性體也可以是柔性體。這類接觸問題的一個典型例子是插頭插到插座里。使用這類接觸單元，不需要預先知道確切的接觸位置，接觸面之間也不需要保持一致的網(wǎng)格。 并且允許有大的變形和大的相對滑動， 雖然這一功能也可以單 擬小的滑動。contact48 contact4單元是點一面的接觸單元。這2種單元支持大滑 動、大變形、以及接觸部件間不同的網(wǎng)格。用戶也可以用這2 種單元來進行熱-機械耦合分析，其中熱在接觸實體之間的傳導非常重要。應用contact26單元用來模擬柔性點一剛性面的接觸。對有不光滑剛性面 的問題，不推薦采用contact2科元，因為在這種環(huán)境下，可能導致接觸的丟

20、 失。在這種情況下，contact4逝過使用偽單元算法，能提供較好的建模能力（參 見 ansys theory reference ） ，但如果目標面嚴重不連續(xù)，依然可能失敗。5.3.3 點一點接觸單元接一點接觸單元主要用于模擬點一點的接觸行為。為了使用點一點接觸單元， 用戶需要預先知道接觸位置， 這類接觸問題只能適用于接觸面之間有較小相對滑動的情況（ 即使在幾何非線性情況下 ） 。其中一個例子是傳統(tǒng)的管道裝配單型，其中接觸點總是在管端和約束之間。接一點接觸單元也可以用于模擬面一面的接觸問題，如果兩個面上的節(jié)點一一對應， 相對滑動又可以忽略不計， 兩個面位移 （ 轉(zhuǎn)動 ） 保持小量， 那么可以

21、用點一點的接觸單元來求解面一面的接觸問題，過盈裝配問題是一個用點一點的接觸單元來模擬面一面接觸問題的典型例子。另一個點一點接觸單元的應用是表面應力的精確分析，如透平機葉片的分 析。ansys勺conta17四元是大多數(shù)點-點接觸問題的最好選擇。它比其他單元 提供了范圍更廣的選項和求解類型。contac1建口 contac52i元保留的理由，在很大程度上是為了與已有單型的向下兼容。5.4 面一面的接觸分析用戶可以應用面-面接觸單元來單擬剛體 -柔體或柔體之間的接觸。從菜單（preprocessorcreatecontact paircontact wizard） 進入接觸向?qū)?，為大多?shù)接觸問題建

22、立接觸對提供了簡單的方法。 接觸向?qū)⒅笇в脩艚⒔佑|對的整個過程。每個對話框中的hel限鈕對其應用及選項作了詳細說明。在用戶未對單型的任何區(qū)域分網(wǎng)之前，接觸向?qū)Р荒軕?。如果用戶希望建立剛體 - 柔體單型， 則在進入接觸向?qū)埃?僅對用作柔體接觸面的部分分網(wǎng) （ 不對 剛體目標面分網(wǎng) ） 。如用戶希望建立柔體- 柔體接觸單型，則應在進入接觸向?qū)?，對所有用作接觸面的部件進行分網(wǎng) （ 包括目標面 ） 。下面諸節(jié)將論述不用接觸向?qū)斫⒔佑|面和目標面的方法。5.4.1 應用面 -面接觸單元在涉及到兩個邊界的接觸問題中，很自然把一個邊界作為“目標”面，而把另一個作為“接觸”面。對剛體一柔體的接觸，

23、目標面總是剛性面，接觸面總是柔性面。對柔體一柔體的接觸，目標面和接觸面都與變形體關(guān)聯(lián)。這兩個面合起來叫作“接觸對”。使用 targe16的conta17貨conta1721l元來定義2-d 接觸對。使用targe17的conta17緘conta17令元來定義3-d接觸對。程序 通過相同的實常數(shù)號來識別每一個接觸對。5.4.2 接觸分析的步驟典型面一面接觸分析的基本步驟如下，后面將對每一步驟進行詳細解釋。1、建立幾何單型并劃分網(wǎng)格；2、識別接觸對； 3、指定接觸面和目標面；4、定義目標面； 5、定義接觸面；6、設(shè)置單元關(guān)鍵選項和實常數(shù)；7、定義控制剛性目標面的運動（ 僅適用于剛體- 柔體接觸 ）

24、 ；8、施加必須的邊界條件；9、定義求解選項和載荷步；10、求解接觸問題；11、查看結(jié)果。5.4.3 建立幾何模型并劃分網(wǎng)格在這一步，用戶需要建立代表接觸體的幾何實體模型。與其它分析一樣，需要設(shè)置單元類型、實常數(shù)、材料特性。用恰當?shù)膯卧愋徒o接觸體劃分網(wǎng)格。參見 ansys modeling and meshing guide 。命令： ameshvmeshgui： main menupreprocessormesh5.4.4 識別接觸對用戶必須判斷模型在變形期間哪些地方可能發(fā)生接觸。 一旦已經(jīng)判斷出潛在的接觸面， 就應該通過目標單元和接觸單元來定義它們， 目標和接觸單元將跟蹤變形階段的運動。

25、 構(gòu)成一個接觸對的目標單元和接觸單元通過共享的實常數(shù)號聯(lián)系起來。接觸區(qū)域可以任意定義，然而為了更有效地進行計算（主要指cpuh問），用戶可能想定義更小的局部化的接觸區(qū)域， 但要保證它足以描述所需要的所有接觸 行為。 不同的接觸對必須通過不同的實常數(shù)號來定義， 即使實常數(shù)沒有變化。 但 不限制允許的面的數(shù)目。9歡迎下載 。精品文檔圖5-1局部接觸區(qū)域由于幾何模型和潛在變形的多樣性，有時候一個接觸面的同一區(qū)域可能與多 個目標面產(chǎn)生接觸關(guān)系。在這種情況下，應該定義多個接觸對（使用多組覆蓋接觸單元）。每個接觸對有不同的實常數(shù)號。見圖5-1 o5.4.5 指定接觸面和目標面接觸單元被限制不得穿透目標面。

26、但是，目標單元可以穿透接觸面。對于剛 體-柔體接觸，目標面總是剛體表面，而接觸面總是柔體表面。對于柔體-柔體接 觸，選擇那一個面作為接觸面或目標面可能會引起穿透量的不同，從而影響求解結(jié)果。這可參照下面的論述：如凸面預期與一個平面或凹面接觸，則平面 /凹面應當指定為目標面；如一個面有較密的網(wǎng)格，而相比較之下，另一個面網(wǎng)格較粗，則較密網(wǎng)格的面應當是接觸面，而較粗網(wǎng)格的面則為目標面；如一個面比另一個面剛，則較柔的面應當指定為接觸面，而較剛的面 則為目標面；如果高階單元位于一個外表面，而低階單元位于另一個面，則前者應指定為接觸面，后者則為目標面；如果一個面明顯地比另一個面大（如一個面包圍其他面），則較

27、大的面 應指定為目標面。上面的論述對于不對稱接觸是正確的。但不對稱接觸可能不能滿足模型需要。 下面一小節(jié)祥細論述不對稱接觸和對稱接觸的差異， 并簡要說明需要對稱接 觸的一些場合。5.4.6 不對稱接觸與對稱接觸不對稱接觸定義為所有的接觸單元在一個面上， 而所有的目標單元在另一個面上的情況。有時候也稱為“單向接觸”。這在模擬面-面接觸時最為有效。但是，在某些環(huán)境下，不對稱接觸不能滿足要求。在這些情況下，可以把任一個面指定為目標面和接觸面。然后在接觸的面之間生成二組接觸對( 或僅是一個接觸對，如自接觸情況) 。這就稱為對稱接觸，有時也稱為“雙向接觸”。顯然，對稱接觸不如非對稱接觸效率高。但是，許多

28、分析要求應用對稱接觸( 典型地，是為了減少穿透) 。要求對稱接觸的情況如下：接觸面和目標面區(qū)分不十分清楚；二個面都有十分粗糙的網(wǎng)格。對稱接觸算法比非對稱接觸算法在更多的面上施加了接觸約束條件。如果二個面上的網(wǎng)格相同并且足夠密， 則對稱接觸算法可能不會顯著改變運行，而事實上可能更費cpuh問。在這種情況下，拾取一個面為目標面，而另一個面為接觸面。 在任何接觸模型中， 可以混合不同的接觸對： 剛體 -柔體或柔體 -柔體接觸對；對稱接觸或非線稱接觸。但在一個接觸對中只能有一種類型。5.4.7 定義目標面目標面可以是2d或3d的剛體或柔體的面。對于柔體目標面，一般應用esurf 命令來沿現(xiàn)有網(wǎng)格的邊界

29、生成目標單元。 也可以按相同的方法來生成柔體接觸面( 見 5.4.8) 。用戶不應當應用下列剛性目標面作為柔體接觸面： arc,carc,circ, cyl1, cone, sphere pilo。對于剛體目標面的情況論述如下。在2d情況下，剛性目標面的形狀可以通過一系列直線、圓弧和拋物線來描 述，所有這些都可以用targe169單元來表示。另外，可以使用它們的任意組合 來描述復雜的目標面。在3d情況下，目標面的形狀可以通過三角面、圓柱面、 圓錐面和球面來描述，所有這些都可以用tapge170單元來表示。對于一個復雜 的、任意形狀的目標面，可以使用低階/ 高階三角形和四邊形來給它建模。5.4.

30、7.1 控制節(jié)點剛性目標面可能會與“控制 (pilot)節(jié)點”聯(lián)系起來，它實際上是一個只有一個節(jié)點的單元， 其運動控制整個目標面的運動， 因此可以把控制節(jié)點作為剛性目標的控制器。整個目標面的力 /力矩和轉(zhuǎn)動/ 位移可以只通過控制節(jié)點來表示。控制節(jié)點可能是目標單元中的一個節(jié)點， 也可能是一個任意位置的節(jié)點。 只有當需要轉(zhuǎn)動或力矩載荷時， 控制節(jié)點的位置才是重要的。 如果用戶定義了控制節(jié)點，ansysi序只在控制節(jié)點上檢查邊界條件，而忽略其它節(jié)點上的任何約束。注意 - 當前的接觸向?qū)Р恢С稚煽刂乒?jié)點。用戶可以在接觸向?qū)舛x控制節(jié)點。5.4.7.2 基本圖元用戶可以使用基本幾何圖元，如圓、圓柱、

31、圓錐、球，來模擬目標面（ 它需要實常數(shù)來定義半徑） 。也可以組合圖元與一般的直線、拋物線、三角形和四邊形來定義目標面。5.4.7.3 單元類型和實常數(shù)在生成目標單元之前，首先必須定義單元類型（2維的targe16第元，或3維的targe17解元）。命令： etgui： main menupreprocessorelement type add/edit/delete隨后必須設(shè)置目標單元的實常數(shù)。命令： realgui： main mennpreprocessorreal constants對于targe169單元和targe170單元，僅需設(shè)置實常數(shù) r1和r2（如果需要 的話 ） 。關(guān)于目標

32、單元、單元形狀、實常數(shù)的完整描述，參見 ansys elements reference 中 targe169i元和 targe170i元的論述。注意 - 只有在使用直接生成法建立目標單元時， 才需要指定實常數(shù)r1、 r2。另外除了直接生成法，用戶也可以使用ansy制格劃分工具生成目標單元，下面 解釋這兩種方法。5.4.7.4 使用直接生成法建立剛性目標單元為了直接生成目標單元，使用下面的命令和菜單。命令： tshapgui： main menupreprocessormodeling-createelementselem attributes隨后指定單元形狀，可能的形狀有：直線 (2d)拋物

33、線 (2-d)順時針的圓弧(2-d)反時針的圓弧(2-d)圓(2-d)三角形 (3-d)圓柱 (3-d)圓錐 (3-d)球(3-d)控制節(jié)點 (2-d 和 3-d)一旦用戶指定目標單元形狀，所有以后生成的單元都將保持這個形狀，除非用戶指定另外一種形狀。注意-不能在同一個目標面上混合 2d和3d目標單元。注意 - 不能在同一個目標面上混合剛體目標單元和柔體目標單元。在求解期間，ansys寸具有下伏單元的目標單元指定為可變形狀態(tài)，而對沒有下伏單元 的目標單元指定為剛體狀態(tài)。 如果刪除柔性表面的下伏單元的一部分， 在求解時 會出現(xiàn)一個錯誤。用戶可以用標準的ansyst接生成技術(shù)生成節(jié)點和單元。參見a

34、nsys modeling and meshing guide 9。命令： negui： main menupreprocessor modeling- createnodesmain menupreprocessor modeling- createelements在建立單元之后，可以通過列表單元來驗證單元形狀。命令： elistgui： utility menulistelementsnodes+attributes5.4.7.5 使用ansy啊格劃分工具生成剛性目標單元13歡迎下載。精品文檔用戶也可以用標準的ansy制格劃分功能讓程序自動地生成目標單元。ansys 程序會基于體模型生成合

35、適的目標單元形狀，而忽略tsha騎令的選項。為了生成一個控制(pilot)節(jié)點，使用下面的命令或 gui路徑：命令：kmeshgui: main menupreprocessormeshing-meshkeypoints注意：kmesh總是生成控制節(jié)點。為了生成一個2d剛性目標單元，使用下面的命令和gui路徑。ansyse每條 線上生成一條單一的線，在b-樣條曲線上生成拋物線線段，在每條圓弧和倒角 線上生成圓弧線段，參見 圖5-2。如果所有的圓弧形成一個封閉的圓，ansys 成一個單一的圓，參見 圖5-3。但是，如果圍成封閉圓的弧是從外部輸入(如 iges)的幾何實體，則ansyst能無法生成

36、一個單一的圓。命令：lmeshgui: main menupreprocessormeshing-meshlinespifftnli 12wot sc wd tdi g1!圖5-2 ansys幾何實體和相應的剛性目標單元圖5-3從圓弧線段生成單一的圓為了生成3d的目標單元，使用下面的命令或 gui路徑。命令：ameshgui: main menupreprocessor-meshing-meshareas如果實體模型的表面部分形成了一個完整的球、圓柱或圓錐，那么ansy鄰序通過amesk令，自動生成一個基本3d目標單元。因為生成較少的單元，從 而使用戶分析計算更有效率。對任意形狀的表面，應該使

37、用amesk令來生成目標單元。在這種情況下，網(wǎng)格形狀的質(zhì)量不重要。而目標單元的形狀是否能較 好地模擬剛性面的表面幾何形狀顯得更重要。在所有可能的面上，推薦使用映射網(wǎng)格。如果在表面邊界上沒有曲率，則在 網(wǎng)格劃分時，指定那條邊界分為一份。剛體 tareg16單元總是在一根線上按一 個單元這樣來分網(wǎng)，而忽略lesize命令的設(shè)置。缺省的單元形狀是四邊形。如 果要用三角形目標單元，應用mshape。圖5-4示出任意目標面的網(wǎng)格布局。下面的命令或gui路徑，將盡可能地生成一個映射網(wǎng)格（如果不能進行映射， 它將生成自由網(wǎng)格）。命令：mshkfy, 2gui: main menupreprocessor-m

38、eshing-mesh-areas-target surf圖5-4任意目標面的網(wǎng)格布局如果目標面是平面（或接近平面），用戶可以選擇低階目標單元（3節(jié)點三角形 或4節(jié)點四邊形單元）。如果目標面是曲面，用戶應該選擇高階目標單元（6節(jié)點 三角形或8節(jié)點四邊形單元）。為此在目標單元定義中設(shè)置 keyopt（1）=1注意-低階單元致使獲取穿透和間隙時 cpu勺開銷較小；但是，分網(wǎng)后的 面可能不夠光滑。高階單元則在獲取穿透和間隙時 cpu勺開銷較大；但是需要較 少的單元就可以離散整個目標曲面。注意-如果通過程序分網(wǎng)（kmesh lmesh esurf命令）來建立目標單元, 則忽略tsha項令，而ansys

39、1動選擇合適的形狀。5.4.7.5.1 建模和網(wǎng)格劃分的一些訣竅一個目標面可能由兩個或多個不連續(xù)的區(qū)域組成。用戶應該盡可能地通過定義多個目標面，來使接觸區(qū)域限于局部（每個目標面有一個不同的實常數(shù)號）。剛 性面上的形狀不限制，不要求光滑。但是，要保證剛性目標面上曲面的離散足夠。 過粗的網(wǎng)格離散可能導致收斂問題。 如果剛性面有一個尖銳的凸角，求解大的滑 動問題時很難獲得收斂結(jié)果。為了避免這些建模問題，在實體模型上使用線或面 的倒角來使尖角光滑化，或者在曲率突變的區(qū)域使用更細的網(wǎng)格或使用高階單 元，見圖5-5 。圖5-5凸角的光滑化5.4.7.5.2 檢驗目標面的節(jié)點號順序（接觸方向）目標面的節(jié)點號

40、順序是重要的，因為它定義了接觸方向。對 2d接觸問題， 當沿著目標線從第一個節(jié)點移向第二個節(jié)點時， 變形體的接觸單元必須位于目標 面的右邊。見圖5-6。圖5-6正確的節(jié)點順序?qū)?d接觸問題，目標三角形單元號應使剛性面的外法線方向指向接觸面。 外法線通過右手法則來定義。為了檢查法線方向，顯示單元坐標系。15欠迎下載精品文檔命令： / psymb ， esys， 1gui： utility menuplotctrlssymbols如果單元法向不指向接觸面，選擇該單元，反轉(zhuǎn)表面法線的方向。命令： esurf,revegui： main menupreprocossorcreateelementssu

41、rf to surf或重新定向單元的法向：命令： enormgui： main menupreprocessorcreatemove/modifyshell normals注意 - 在目標元素（ 如完整的圓、圓柱、圓錐、球） 上的接觸，只能在這些目標元素的外表面上出現(xiàn)。5.4.8 定義柔體的接觸面為了建立柔體的接觸面，對于 2d接觸必須使用接觸單元confa171或 confa172t觸單元；對于3d接觸必須使用conta173或conta174接觸單元。程序通過組成柔體表面的接觸單元來定義接觸面。 接觸單元與下伏柔體單元有同樣的幾何特性。接觸單元與下伏柔體單元必須處于同一階次（ 低階或高階）

42、 ，以使在邊上的節(jié)點協(xié)調(diào)。 高階接觸單元可以通過消除中節(jié)點而與低階下層單元匹配。下伏單元可能是實體單元、殼單元、2d梁單元。接觸面可以在殼或梁單元任何一邊。 下伏單元也可以是超單元。 但是， 軸對稱調(diào)和單元不能用作下伏單元。與目標面單元一樣，用戶必須定義接觸面的單元類型，然后選擇正確的實常數(shù)號 （ 在每個接觸對中， 實常數(shù)號必須與它所對應的目標面的實常數(shù)號相同 ） ， 最 后生成接觸單元。5.4.8.1 單元類型下面簡單描述四種類型的接觸單元。參見 ansys elements reference 。conta171這是2d 2個節(jié)點的低階線單元，可位于 2d實體、殼或梁 單元（如beam3

43、plane4減shell51）的表面。conta172這是2d 3節(jié)點的高階拋物線形單元，可位于有中節(jié)點的2d實體或梁單元（如plane82或visco88）的表面。conta173這是3d 4節(jié)點的低階四邊形單元，可位于 3d實體或殼單 元（如solid45或shell181）的表面。可退化成3節(jié)點的三角形單元。conta174這是3d 8節(jié)點的高階四邊形單元，可位于有中節(jié)點的 3d 實體或殼單元（如solid92、solid95或shell93）的表面?？赏嘶? 節(jié)點的三角形單元。命令： etgui： main menupreprocessorelement typeadd/edit/d

44、elete5.4.8.2 實常數(shù)和材料特性在定義了單元類型之后，需要選擇正確的實常數(shù)集。一個接觸對中的接觸面和目標面必須有相同的實常數(shù)號。每個接觸對必須有不同的實常數(shù)號。ansys用下伏單元的材料特性來計算一個合適的接觸（或罰）剛度。在下層單元有用 tb 命令定義的塑性材料特性 （不論激活與否 ） 的情況下， 接觸的法向剛 度可能按照系數(shù)100降低。ansys自動為切向（滑動）剛度定義一個與mu和法向 剛度成正比的缺省值。 如果下伏單元是一個超單元， 接觸單元的材料必須與超單 元形成時的原始結(jié)構(gòu)單元相同。5.4.8.3 生成接觸單元既可以通過直接生成法生成接觸單元， 也可以在下層單元的外表面上

45、自動生成接觸單元。 推薦采用自動生成法， 這種方法更為簡單和可靠。 可以通過下面三 個步驟，來自動生成接觸單元。1、選擇節(jié)點選擇已分網(wǎng)的柔體表面的節(jié)點。對每一個面，檢查節(jié)點排列。如果用戶確定某一部分節(jié)點永遠不會接觸到目標面， 用戶可以忽略它以便減少計算時間。 然而 用戶必須保證設(shè)有漏掉可能會接觸到目標面的節(jié)點。命令： nselgui： utility menuselectentities2、生成接觸單元命令： esurfgui： main menupreprocessorcreateelementssurf to surf如果接觸單元是附在已用實體單元劃分網(wǎng)格的面或體上， 程序會自動決定接觸計

46、算所需的外法向。如果下層單元是梁或殼單元，則必須指明哪個表面（ 上表面或下表面） 是接觸面。命令： esurf,top （ 或 botiom）17歡迎下載。精品文檔gui: main menupreprocessorcreateelementssurf to surf使用top（缺?。┥山佑|單元，它們的外法向與梁或殼單元的法向相同；使 用botiom生成接觸單元，則它們的外法向與梁或殼單元的法向相反。必須確保梁上的單元或殼單元有一致的法向。如果下伏單元是實體單元，則 top或 botto憂項不起作用3、檢查接觸單元外法向。當程序進行是否接觸的檢查時，接觸面的外法線 方向至關(guān)重要。對于3d單元

47、，按節(jié)點順序號以右手法則來決定單元的外法向。 接觸面的外法向應該指向目標面。否則，在開始分析計算時，程序可能會認為是有過度穿透的面，而很難找到初始解。在這些情況下，程序一般會立即停止執(zhí)行。 圖5-7說明正確和不正確的外法向。命令：/ psymresysgui: utility menuplotctrlssymbols圖5-7定義接觸單元的外法向當發(fā)現(xiàn)單元的外法線方向不正確時，必須通過反轉(zhuǎn)所選擇的不正確單元的節(jié) 點號來改變它們：命令：esurf,revegui: main menupreprocessorcreateelementssurf to surf或重新定向單元法向：命令：enormgu

48、i: main menupreprocessorcreatemove/modifyshell normals5.4.9 設(shè)置實常數(shù)和單元關(guān)鍵選項程序使用20個實常數(shù)和數(shù)個單元關(guān)鍵選項，來控制面一面接觸單元的接觸， 參見ansys elements reference 中對接觸單元的描述。5.4.9.1 實常數(shù)在20個實常數(shù)中，兩個（r1和r2）用來定義目標面單元的幾何形狀。剩下的 用來控制接觸面單元。r1 和 r2定義目標單元幾何形狀。fkn定義法向接觸剛度因子。ftoln是基于單元厚度的一個系數(shù)，用于計算允許的穿透。icont定義初始閉合因子。pinb定義“pinball 區(qū)域。pmin和

49、pmax定義初始穿透的容許范圍。taumar 指定最大的接觸摩擦。cnof指定施加于接觸面的正或負的偏移值。fkop指定在接觸分開時施加的剛度系數(shù)。fkt指定切向接觸剛度。cohe制定滑動抗力粘聚力。tcc指定熱接觸傳導系數(shù)。fhtg指定摩擦耗散能量的熱轉(zhuǎn)換率。sbct指定 stefan-boltzman常數(shù)。rdvf指定輻射觀察系數(shù)。fwgt指定在接觸面和目標面之間熱分布的權(quán)重系數(shù)。fact靜摩擦系數(shù)和動摩擦系數(shù)的比率。dc靜、動摩擦衰減系數(shù)。命令： rgui： main menu preprocessorreal constant對實常數(shù) fkn, ftoln, icont, pinb,

50、pmax, pmin, fkop 和 fkt,用戶既 可以定義一個正值， 也可以定義一個負值。 程序?qū)⒄底鳛楸壤蜃樱?將負值作 為絕對值。程序?qū)⑾路鼏卧暮穸茸鳛?icon ftoln pinb, pmax和pmin的參考值。例如icon = 0.1表明初始閉合因子是“ 0.1*下層單元的厚度”。然而, icon = -0.1則表示真實調(diào)整帶是0.1單位。如果下伏單元是超單元，則將接 觸單元的最小長度作為厚度。參見圖5-8 。圖5-8下層單元的厚度2彼迎下載在模型中，如果單元尺寸變化很大，而且在實常數(shù)如icont, ftoln, pinb,pmax, pmin中應用比例系數(shù)，則可能會出現(xiàn)問

51、題。因為從比例系數(shù)得到的實際 結(jié)果，取決于下層單元的厚度，這就可能引起大、小單元之間的重大變化。如果 出現(xiàn)這一問題，請用絕對值代替比例系數(shù)。tcc, fhtg, sbct, rdvf和 fwgt僅用于熱接觸分析keyopt(1)=1。5.4.9.2 單元關(guān)鍵選項每種接觸單元都包括數(shù)個關(guān)鍵選項。對大多的接觸問題，缺省的關(guān)鍵選項是 合適的。而在某些情況下，可能需要改變?nèi)笔≈?。下面是可以控制接觸行為的一 些關(guān)鍵選項：keyopt(1)自由 度接觸算法(罰函數(shù)十拉格朗日乘子或罰函數(shù))keyopt(2)存在超單元時的應力狀態(tài)(僅2d)keyopt(3)接觸檢測點的位置(僅低階接觸單元) keyopt(4

52、)cnof1動調(diào)整keyopt(5)時間步控制keyopt偽接觸預防keyopt(8)初始穿透或間隙的影響keyopt(9)法向和切向接觸剛度修正方法控制 keyopt(10)殼的厚度影響 keyopt(11)接觸面行為 (粗糙、綁定等) keyopt(12)命令： keyoptetgui： main menupreprocessorelemant typeadd/edit/delete5.4.9.3 選擇接觸算法contact wizard basiccontact algorithm對面一面接觸單元，程序可以使用增進的拉格朗日方法或罰函數(shù)方法。通過 單元關(guān)鍵字keyopt(2沫指定。增進的

53、拉格朗日方法是為了找到精確的拉格朗日乘子(即接觸力)，而對罰函數(shù)進行一系列修正迭代。 與罰函數(shù)的方法相比， 拉格朗日方法容易得到良態(tài)條件，對接觸剛度的敏感性較小。然而，在有些分析中，增進的拉格朗日方法可能需要更多的迭代，特別是在變形后網(wǎng)格變得太扭曲時。使用拉格朗日方法的同時應使用實常數(shù)ftoln ftoln為拉格朗日方法指定容許的最大穿透。 如果程序發(fā)現(xiàn)穿透大于此值時， 即使不平衡力和位移增量已經(jīng)滿足了收斂準則，總的求解仍被當作不收斂處理。ftlon的缺省值為0.1。用戶可以改變這個值， 但要注意， 如果此值太小， 可能會造成太多的迭代次數(shù)或者不 收斂。5.4.9.4 確定接觸剛度contac

54、t wizardbasicnormal penaltycontact stiffness所有的接觸問題都需要定義接觸剛度， 兩個表面之間穿透量的大小取決于接觸剛度。 過大的接觸剛度可能會引起總剛矩陣的病態(tài)， 從而造成收斂困難。 一般來說， 應該選取足夠大的接觸剛度以保證接觸穿透小到可以接受， 但同時又應該讓接觸剛度足夠小以不致引起總剛矩陣的病態(tài)而保證收斂性。ansyss序根據(jù)下伏柔體單元的材料特性，來估計一個缺省的接觸剛度值。用戶可用實常數(shù)fkn 來為接觸剛度指定一個比例因子或指定一個絕對值。比例因子一般在0.01 和 10之間； 對于大變形問題， 選 1 是比較好的； 而對于彎曲為主的問題，通常為0.010.1。用戶應當總是檢驗以使穿透到達極小值，而又避 免過多的迭代次數(shù)。精品文檔注意-ftoln和fkn從一個荷載步到另一個荷載步中，都可以修改。也可 以在重啟動中修改。這時，必須定義 keyopt（10）=1,2為了確定一個較好的接觸剛度值，可能需要一些經(jīng)驗。用戶可以按下面的步 驟來進行嘗試：1、開始時取一個較低的值。低估值要比高估值好，因為由一個較低的接觸 剛度導致的穿透問題，比過高的接觸剛度導致的收斂性困難，要容易解 決。2、對前幾個子步進行計算分析，直到最終荷載的一個比例 （剛好