ANSYS-高級(jí)接觸問(wèn)題.ppt

1、ANSYS 高級(jí)接觸問(wèn)題,接觸問(wèn)題概述 在工程中會(huì)遇到大量的接觸問(wèn)題，如齒輪的嚙合、法蘭聯(lián)接、機(jī)電軸承接觸、卡頭與卡座、密封、板成形、沖擊等等。接觸是典型的狀態(tài)非線性問(wèn)題，它是一種高度非線性行為。接觸例子如圖1：,分析中常常需要確定兩個(gè)或多個(gè)相互接觸物體的位移、接觸區(qū)域的大小和接觸面上的應(yīng)力分布。 接觸分析存在兩大難點(diǎn)： 在求解之前，你不知道接觸區(qū)域的范圍；表面之間是接觸還是分開(kāi)是未知的；表面之間突然接觸或突然不接觸會(huì)導(dǎo)致系統(tǒng)剛度的突然變化。 大多數(shù)接觸問(wèn)題需要計(jì)算摩擦。摩擦是與路徑有關(guān)的現(xiàn)象，摩擦響應(yīng)還可能是雜亂的，使問(wèn)題求解難以收斂。,ANSYS 高級(jí)接觸問(wèn)題,1 接觸分類,剛?cè)?一個(gè)表面

2、是完全剛性的除剛體運(yùn)動(dòng)外無(wú)應(yīng)變、應(yīng)力和變形，另一表面為軟材料構(gòu)成是可變形的。 只在一個(gè)表面特別剛硬并且不關(guān)心剛硬物體的應(yīng)力時(shí)有效。 柔柔 兩個(gè)接觸體都可以變形。,2 接觸單元,ANSYS采用接觸單元來(lái)模擬接觸問(wèn)題： 跟蹤接觸位置； 保證接觸協(xié)調(diào)性（防止接觸表面相互穿透）； 在接觸表面之間傳遞接觸應(yīng)力（正壓力和摩擦）。 接觸單元就是覆蓋在分析模型接觸面上的一層單元。 在 ANSYS 中可以采用三種不同的單元來(lái)模擬接觸： 面一面接觸單元； 點(diǎn)一面接觸單元； 點(diǎn)一點(diǎn)接觸單元。,不同的單元類型具有完全不同的單元特性和分析過(guò)程。 1. 面一面接觸單元用于任意形狀的兩個(gè)表面接觸 不必事先知道接觸的準(zhǔn)確位置

3、； 兩個(gè)面可以具有不同的網(wǎng)格； 支持大的相對(duì)滑動(dòng)； 支持大應(yīng)變和大轉(zhuǎn)動(dòng)。 例如: 面一面接觸可以模擬金屬成型，如軋制過(guò)程。,2 接觸單元,2 接觸單元,2. 點(diǎn)一面接觸單元用于某一點(diǎn)和任意形狀的面的接觸 可使用多個(gè)點(diǎn)面接觸單元模擬棱邊和面的接觸； 不必事先知道接觸的準(zhǔn)確位置； 兩個(gè)面可以具有不同的網(wǎng)格； 支持大的相對(duì)滑動(dòng)； 支持大應(yīng)變和大轉(zhuǎn)動(dòng)。 例：點(diǎn)面接觸可以模擬棱邊和面之間的接觸,2 接觸單元,2 接觸單元,3. 點(diǎn)點(diǎn)接觸單元用于模擬單點(diǎn)和另一個(gè)確定點(diǎn)之間的接觸。 建立模型時(shí)必須事先知道確切的接觸位置； 多個(gè)點(diǎn)點(diǎn)接觸單元可以模擬兩個(gè)具有多個(gè)單元表面間的接觸； 每個(gè)表面的網(wǎng)格必須是相同的；

4、相對(duì)滑動(dòng)必須很?。?只對(duì)小的轉(zhuǎn)動(dòng)響應(yīng)有效。 例如: 點(diǎn)一點(diǎn)接觸可以模擬一些面的接觸。如地基和土壤的接觸,2 接觸單元,2 接觸單元,3 關(guān)于耦合和約束方程的應(yīng)用,如果接觸模型沒(méi)有摩擦，接觸區(qū)域始終粘在一起，并且分析是小撓度、小轉(zhuǎn)動(dòng)問(wèn)題，那么可以用耦合或約束方程代替接觸。 使用耦合或約束方程的優(yōu)點(diǎn)是分析還是線性的,接觸問(wèn)題的一般特性,1 接觸剛度 1、所有的 ANSYS 接觸單元都采用罰剛度（接觸剛度）來(lái)保證接觸界面的協(xié)調(diào)性,在數(shù)學(xué)上為保持平衡，需要有穿透值 然而，物理接觸實(shí)體是沒(méi)有穿透的 分析者將面對(duì)困難的選擇： 小的穿透計(jì)算精度高，因此接觸剛度應(yīng)該大； 然而，太大的接觸剛度會(huì)產(chǎn)生收斂困難：模

5、型可能會(huì)振蕩，接觸表面互相跳開(kāi)。 接觸剛度是同時(shí)影響計(jì)算精度和收斂的最重要的參數(shù)。你必須選定一個(gè)合適的接觸剛度。 除了在表面間傳遞法向壓力外，接觸單元還傳遞切向運(yùn)動(dòng)（摩擦）。采用切向罰剛度保證切向的協(xié)調(diào)性。（圖12） 作為初值，可采用：Ktangent=0.01 Knormal 切向罰剛度與法向罰剛度以同樣的方式對(duì)收斂性和計(jì)算精度產(chǎn)生影響。,2、接觸剛度的選取 選定一個(gè)合適的接觸剛度值需要一些經(jīng)驗(yàn)。 對(duì)于面一面接觸單元，接觸剛度通常指定為基體單元?jiǎng)偠鹊囊粋€(gè)比例因子。 開(kāi)始估計(jì)時(shí)，選用 FKN = 1.0 大面積實(shí)體接觸 FKN = 0.01-0.1 較柔軟（彎曲占主導(dǎo)的部分） 另外，也可以指定

6、一個(gè)絕對(duì)剛度值，單位：（力/長(zhǎng)度）/ 面積。 點(diǎn)一點(diǎn)（除 CONTA178）和點(diǎn)面接觸單元需要為罰剛度 KN 輸入絕對(duì)值： 初始估計(jì)時(shí)： 對(duì)于大變形: 0.1*E KN 1.0*E 對(duì)于彎曲: 0.01*E KN 0.1*E E 為彈性模量,3、選取接觸剛度的指導(dǎo)： Step 1. 開(kāi)始采用較小的剛度值 Step 2. 對(duì)前幾個(gè)子步進(jìn)行計(jì)算 Step 3. 檢查穿透量和每一個(gè)子步中的平衡迭代次數(shù) 在粗略的檢查中，如以實(shí)際比例顯示整個(gè)模型時(shí)就能觀察到穿透，則穿透可能太大了，需要提高剛度重新分析。 如果收斂的迭代次數(shù)過(guò)多（或未收斂），降低剛度重新分析。 注意：罰剛度可以在載荷步間改變，并且可以在重

7、啟動(dòng)中調(diào)整。 牢記：接觸剛度是同時(shí)影響計(jì)算精度和收斂性的最重要的參數(shù)。如果收斂有問(wèn)題，減小剛度值，重新分析 在敏感的分析中，還應(yīng)該改變罰剛度來(lái)驗(yàn)證計(jì)算結(jié)果的有效性。 在分析中減小剛度范圍，直到結(jié)果（接觸壓力、最大SEQV 等）不再明顯改變。,2 摩擦,1、兩個(gè)接觸體的剪切或滑動(dòng)行為可以是無(wú)摩擦的或有摩擦的 無(wú)摩擦?xí)r允許物體沒(méi)有阻力地相互滑動(dòng)； 有摩擦?xí)r，物體之間會(huì)產(chǎn)生剪切力。 2、摩擦消耗能量，并且是路徑相關(guān)行為。 為獲得較高的精度，時(shí)間步長(zhǎng)必須?。▓D2-1）,式中： 摩擦系數(shù) 一旦所受剪力超過(guò) FT，兩物體將發(fā)生相對(duì)滑動(dòng)。 4、彈性庫(kù)侖摩擦模型：允許粘著和滑動(dòng)。,圖21 3、ANSYS 中，

8、摩擦采用庫(kù)侖模型，并有附加選項(xiàng)可處理復(fù)雜的粘著和剪切行為。 庫(kù)侖法則是宏觀模型，表述物體間的等效剪力 FT 不能超過(guò)正壓力 FN 的一部分： FT = FN,3 自動(dòng)時(shí)間步、控制,接觸單元的 Keyopt（7）選項(xiàng)控制時(shí)間步的預(yù)報(bào)。 0無(wú)控制：不影響時(shí)間步尺寸。當(dāng)自動(dòng)時(shí)間步開(kāi)關(guān)打開(kāi)時(shí)，對(duì)于靜態(tài)問(wèn)題通常選此項(xiàng)。 1自動(dòng)縮減：如果接觸狀態(tài)改變較大，時(shí)間步二分。對(duì)于動(dòng)態(tài)問(wèn)題，自動(dòng)縮減通常是充分的。 2合理的：比自動(dòng)縮減費(fèi)用更昂貴的算法。為保持一個(gè)合理的時(shí)間載荷增量，需要在接觸預(yù)測(cè)中選擇此項(xiàng)。適用于靜態(tài)分析和連續(xù)接觸時(shí)瞬態(tài)分析。 3最小值：該選項(xiàng)為下一子步、預(yù)報(bào)時(shí)間增量的最小值（計(jì)算費(fèi)用十分昂貴，建議

9、不用）。這個(gè)選項(xiàng)在碰撞和斷續(xù)接觸分析中是有用的。,接觸分析中自動(dòng)時(shí)間步的其它注意事項(xiàng)：,與所有其它非線性分析一樣，對(duì)接觸問(wèn)題，時(shí)間步長(zhǎng)是非常有力的提高收斂性的工具。 采用足夠小的時(shí)間步長(zhǎng)以獲得收斂。 對(duì)于沖擊瞬態(tài)分析，必須使用足夠數(shù)量的計(jì)算步以描述表面間的動(dòng)量轉(zhuǎn)移。 對(duì)于路徑相關(guān)現(xiàn)象（如接觸摩擦），相對(duì)較小的最大時(shí)間步長(zhǎng)對(duì)計(jì)算精度是必須的。,1 概述 面面接觸單元，是模擬任意兩個(gè)表面間接觸的方法。表面可以具有任意形狀。是 ANSYS 中最通用的接觸單元。精度高、特性豐富還可使用接觸向?qū)?，建模方便。（其它接觸單元目前尚不能用向?qū)В?面面接觸單元在面的高斯點(diǎn)處傳遞壓力，這種先進(jìn)技術(shù)使面面接觸單元

10、具有很多優(yōu)點(diǎn)： 與低階單元和高階單元都兼容 提供更好的接觸結(jié)果（于后處理接觸壓力和摩擦應(yīng)力） 可考慮殼和梁的厚度，以及殼的厚度變化 半自動(dòng)接觸剛度計(jì)算 剛性表面由“控制節(jié)點(diǎn) pilot node”控制 熱接觸特性 眾多的高級(jí)選項(xiàng)來(lái)處理復(fù)雜問(wèn)題。 具有眾多的高級(jí)選項(xiàng)（20 個(gè)可用的實(shí)常數(shù)、2 個(gè)材料屬性和 30 個(gè)可用的單元選項(xiàng)）提供了豐富的特征庫(kù)，能夠用于模擬特殊的效果和處理困難的收斂情況。 然而眾多的選項(xiàng)的智能缺省選項(xiàng)可以有效求解許多接觸問(wèn)題而不需要用戶介入太多。 通常的做法是：開(kāi)始使用高級(jí)選項(xiàng)之前，先試著采用缺省設(shè)置：只指定罰剛度，穿透容差和子步數(shù)，然后進(jìn)行分析。只在采用缺省設(shè)置遇到困難時(shí)

11、才采用高級(jí)選項(xiàng)。 所有的高級(jí)選項(xiàng)也可以通過(guò)接觸向?qū)?lái)控制。,2 面面接觸單元,2 面面接觸單元,使用面面接觸單元計(jì)算剛?cè)帷⑷崛峤佑|分析。 把一個(gè)面指定為目標(biāo)面（Target），另一個(gè)面指定為接觸面（contant），合起來(lái)叫接觸對(duì)。 接觸單元被約束不能侵入目標(biāo)面，然而目標(biāo)單元能侵入接觸面。 2D目標(biāo)單元 TARGE169 :,2D 面面接觸單元 CONTA171 2D、2 節(jié)點(diǎn)低階單元，可用于二維實(shí)體、殼、梁?jiǎn)卧谋砻?CONTA172 2D、3 節(jié)點(diǎn)高階單元，可用于帶中間節(jié)點(diǎn)的二維實(shí)體單元表面 3D 目標(biāo)單元 TARGE170,2 面面接觸單元,3 面一面接觸分析步驟、實(shí)例 （不通過(guò)接觸向?qū)?/p>

12、創(chuàng)建接觸對(duì)）,Step 1. 建立基體有限元模型 設(shè)置基體單元類型、實(shí)常數(shù)、材料特性 給基體分網(wǎng)： 命令：AMESH VMESH Step 2. 指定接觸面和目標(biāo)面 對(duì)于剛一柔接觸，目標(biāo)面總是剛性面； 對(duì)于柔柔接觸，目標(biāo)面和接觸面的不同選擇會(huì)產(chǎn)生不同的穿透（圖3-1），并且影響求解精度。,圖3-1 接觸面和目標(biāo)面確定準(zhǔn)則 如凸面和平面或凹面接觸，應(yīng)指定平面或凹面為目標(biāo)面； 如一個(gè)面上的網(wǎng)格較粗而另一個(gè)面上的網(wǎng)格較細(xì)，應(yīng)指定粗網(wǎng)格面為目標(biāo)面； 如一個(gè)面比另一個(gè)面的剛度大，應(yīng)指定剛度大的面為目標(biāo)面； 如一個(gè)面為高階單元而另一面為低階單元，應(yīng)指定低階單元面為目標(biāo)面； 如一個(gè)面比另一個(gè)面大，應(yīng)指定大的

13、面為目標(biāo)面。,例如: 超彈密封,Step 3. 設(shè)置單元選項(xiàng)和實(shí)常數(shù) 接觸對(duì)由實(shí)常數(shù)號(hào)來(lái)定義，接觸單元和目標(biāo)單元必須具有相同的實(shí)常數(shù)。 Step 4. 建立目標(biāo)單元（網(wǎng)格） 此步中所采用的方法依賴于目標(biāo)面是剛性的還是柔性的。 剛性目標(biāo)面采用： 直接生成 (E 命令) 自動(dòng)劃分 (LMESH, AMEAH) 可變形目標(biāo)面采用 Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf(ESURF) 對(duì)于直接生成剛性目標(biāo)面，在建立目標(biāo)單元之前需要要指定附加的單元屬性 TSHAP,剛性目標(biāo)面的自動(dòng)劃分不需要 TS

14、HAP。ANSYS 能根據(jù)實(shí)體模型確定合適的目標(biāo)單元形狀。 劃分線 (LMESH) 2D 剛性目標(biāo)面 劃分面 (AMESH) 3D 剛性目標(biāo)面 創(chuàng)建關(guān)鍵點(diǎn)(KMESH)控制節(jié)點(diǎn)（Pilot） 剛性目標(biāo)面能與控制點(diǎn)聯(lián)系起來(lái)。Pilot 實(shí)際上是只有一個(gè)節(jié)點(diǎn)的單元，通過(guò)這個(gè)節(jié)點(diǎn)的運(yùn)動(dòng)可以控制整個(gè)目標(biāo)面的運(yùn)動(dòng)。ANSYS 只在 Pilot 節(jié)點(diǎn)上檢查邊界條件而忽略其它節(jié)點(diǎn)的約束。 對(duì)可變形體目標(biāo)面建立目標(biāo)單元的步驟是： 1.先選擇可變形體表面上的節(jié)點(diǎn) 2.然后在可變形體上建立目標(biāo)單元 Main Menu Preprocessor Modeling Create Elements Surf/Conta

15、ct Surf to Surf,ANSYS 將根據(jù)基體的網(wǎng)格確定目標(biāo)單元形狀和外法線方向。 檢查外法線方向（這在自動(dòng)劃分剛性目標(biāo)面時(shí)非常重要）圖 3-3 打開(kāi)單元坐標(biāo)系標(biāo)志并重繪單元 /PSYMS,ESYS,1 目標(biāo)單元外法線方向應(yīng)該指向接觸面。如果單元法向不指向接觸面，用命令使之反轉(zhuǎn)： ESURF,REVE,例：Seal.dat (圖3-3) Step 5. 建立接觸面單元 設(shè)置接觸單元屬性、選擇可變形體表面節(jié)點(diǎn)，并在可變形體上建立接觸單元（過(guò)程與在可變形體上建立目標(biāo)單元相同） Main Menu Preprocessor Modeling Create Elements Surf/Cont

16、act Surf to Surf (ESURF) 這些接觸單元與基體有同樣的階數(shù)（低階或高階）。 注意，在殼或梁?jiǎn)卧辖⒛繕?biāo)單元或接觸單元時(shí)，可以選擇要在梁或殼單元的頂層還是底層建立單元。,圖3-3,在選擇柔體表面上的節(jié)點(diǎn)時(shí)，如果你確定某一部分節(jié)點(diǎn)永遠(yuǎn)不會(huì)接觸到目標(biāo)面，可以忽略它，以減少計(jì)算時(shí)間。 接觸面的外法向應(yīng)指向目標(biāo)面。如果發(fā)現(xiàn)外法線方向不正確，用下列命令修改之 ESURF,REVE Step 6. 在有限元模型上施加邊界條件 如果目標(biāo)面是剛性面，目標(biāo)面將會(huì)自動(dòng)固定。 定義了 Pilot 點(diǎn) ANSYS 只檢查該點(diǎn)的邊界條件，忽略目標(biāo)面上其它節(jié)點(diǎn)約束。控制點(diǎn)能控制目標(biāo)面的運(yùn)動(dòng)。 對(duì) S

17、eal.dat 施加的邊界條見(jiàn)圖 33。,Step 7. 定義求解選項(xiàng)和載荷步，以下是默認(rèn)設(shè)置 推薦使用N.L求解自動(dòng)控制 使用不帶自適應(yīng)下降的 full Newton-Raphson 法求解 時(shí)間步必須足夠小。使用自動(dòng)時(shí)間步。 子步數(shù)的最大值（NSBMX）應(yīng)較大，最小值（NSBMIN）應(yīng)較小 Step 8. 求解 Step 9. 后處理 結(jié)果包括位移、應(yīng)力、應(yīng)變和接觸等信息。 接觸壓力、摩擦應(yīng)力、總應(yīng)力、接觸侵入、接觸間隙距離、滑動(dòng)距離和接觸狀態(tài)都可以從 /POST1 或/POST26 中得到。,面一面接觸分析實(shí)例（建立接觸對(duì)不通過(guò)接觸向?qū)В?實(shí)例1：彈性環(huán)裝配,軸對(duì)稱 施加位移載荷：在 L

18、45 線上施加 0.4的 Y 向位移 打開(kāi)幾何 N.L 開(kāi)關(guān) (NLGEOM,on) 設(shè) Time = 0.4 并為自動(dòng)時(shí)間步給出子步數(shù)(20,500,10) 給出輸出控制（要求輸出每一子步結(jié)果） 求解并查看輸出和監(jiān)視文件 重啟動(dòng)分析亦將 Y 向偏移量增加到 0.55使咬接裝配的第 2 個(gè)齒咬合。 問(wèn)題以彎曲為主，設(shè) FKN = 0.1,GUI 方式： Step 1. 恢復(fù)數(shù)據(jù)庫(kù)文件 Snap.db（Snap.db 數(shù)據(jù)庫(kù)包含此有限元模型的完整幾何模型、材料、邊界條件。但未定義目標(biāo)單元與接觸單元）。 Utility Menu File Resume from 選 Snap.db 【OK】 基體

19、單元: Plane42,1 軸對(duì)稱 材料: MAT,1 EX = 0.175e6; NUXY = 0.35 ; MU = 0.0 MAT,2 EX = 90000; NUXY = 0.35 約束節(jié)點(diǎn) 163 UY = 0 182 UY = 0,Step 2. 創(chuàng)建目標(biāo)面單元類型 Main Menu Preprocessor Element Typer Add/Edit/Delete 【Add】 contact 2D targe169 Element typer reference number = 2 【OK】 或命令: /PREP 7 ET,2,Targe169 Step 3. 建立接觸面單

20、元類型 Main Menu Preprocessor Element Typer Add/Edit/Delete 【Add】 contact 2nd Surf 171 Element typer reference number = 3 【OK】 【close】 或命令: ET,3,conta171,Step 4.指定接觸法向剛度 Main Menu Preprocessor Real Constants Add/Edit/Delete 【Add】 Type 3 CONTA171 【OK】 Real Constant Set No. = 1 Normal Penalty Stiffness F

21、KN = 0.1(對(duì)彎曲問(wèn)題采用初值0.1) Penetration tolerance FTOLN = 0.1 (不用此,不收斂) 或命令：R,1,0.1,0.1,確定罰剛度 FKN 值通常在 0.01-10 之間，對(duì)于體積變形問(wèn)題用 1.0（默認(rèn)），對(duì)彎曲問(wèn)題用 0.1 確定侵入容差：侵入容差 (FTOLN) 是與接觸單元下面的實(shí)體單元深度（h） 相乘的比例因子。 若此值太小會(huì)引起收斂困難，絕對(duì)不要用太小的容差!增大罰剛度（FKN）將減少侵入。 將 FKN 增大 100 倍會(huì)相應(yīng)地減少侵入，但是接觸壓力只改變 5。 如不收斂可調(diào)整 FKN 或 FTOLN 值重新運(yùn)行。檢查侵入和每個(gè)子步的平

22、衡迭代數(shù)，如果收斂受侵入容差的驅(qū)使可能是FKN 值估計(jì)不足或 FTOLN 值太小。如果需要多次迭代才能使殘值收斂而不是侵入。FKN 值可能估計(jì)得太高。,Step 5. 創(chuàng)建目標(biāo)單元 (1) 為目標(biāo)面選擇線 Utility Menu Select Entities 拾取線(圖3-5) 【OK】 或命令： LSEL,S,2,4 LSEL,A,15,18 圖3-6 LSEL,A,63 (2).選擇附于線上的全部節(jié)點(diǎn)(圖3-6) Utility Menu Select Entities Nodes Attached to lines,all 【OK】 或命令： NSLL,S,1 NPlot,圖3-5,

23、圖3-6,(3).設(shè)置單元屬性 Main Menu Preprocessor Create Element Elem Attributes 或命令： Type,2 Mat,1 Real,1,(4).創(chuàng)建目標(biāo)單元 Main Menu Preprocessor Modeling Create Elements Surf / Contact Surf to Surf 注意：如果基體單元是 2D 平面或 3D 實(shí)體單元 “Tlab” 無(wú)效，如果基體單元是 2D 單元，“Shape” 無(wú)效。 選 【pick all】 或命令： ESURF(圖3-7),圖3-7,Step 6. 創(chuàng)建接觸單元 (1) 為接

24、觸面選線 Utility Menu Select Entities Lines 圖3-8 By Num/pick From Full 【OK】 選線(圖3-8) 【OK】 或命令: LSEL,s,33,34 LSEL,a,43,44 (2) 選擇附于選定線上的全部節(jié)點(diǎn)圖3-9 Utility Menu Select Entities Nodes Attached to Lines, all From Full 【OK】 或命令: NSLL, S, 1 Nplot(圖3-9),圖3-8,圖3-9,(3) 設(shè)置單元屬性 Main Menu Preprocessor Create Element E

25、lem Attributes Element type number = 3 ConTA 171 Material number = 1 Real constant set number = 1 【OK】 或命令： Type,3 Mat,1 Real,1 (4) 創(chuàng)建接觸面單元圖3-11 Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf to Surf 圖3-10 【OK】 【PICK ALL】(圖3-10) 或命令: ESURF,Step 7. 選擇所有選項(xiàng)并畫單元 Utility Menu Select Ev

26、erything Utility Menu Plot Elements 或命令: Allsel, all Eplot(圖3-11) Step 8. 求解 /solu antype,static nlgeom,on solc, on time,0.4 nsubst,20,500,10 outres,all,all dl,45,uy,-0.4 solve dl,45,uy,-0.55 solve fini Step 9. 后處理,4 應(yīng)用接觸向?qū)?chuàng)建接觸對(duì),所有面面接觸的單元選項(xiàng)和參數(shù)都可以通過(guò)接觸向?qū)?lái)控制。使用接觸向?qū)нM(jìn)行接觸分析方便快捷 ： 自動(dòng)定義單元類型和實(shí)常數(shù) 快速得到接觸選項(xiàng)和參數(shù)

27、快速顯示和反轉(zhuǎn)接觸法向 使用接觸向?qū)П仨毷紫葘?duì)基體進(jìn)行網(wǎng)格劃分，否則不能激活接觸向?qū)А?下面仍以彈性環(huán)裝配接觸問(wèn)題為例。采用接觸向?qū)瓿山佑|對(duì)的創(chuàng)建來(lái)說(shuō)明接觸向?qū)У氖褂梅椒ā?Step 1. 恢復(fù)數(shù)據(jù)文件 Utility Menu File Resume from 選 Snap.db 【OK】(圖4-1) 此數(shù)據(jù)庫(kù) Snap.db 包括此有限元模型的幾何、基體 單元、分網(wǎng)、材料和邊界條件。但不包括接觸單元與目標(biāo)單元及接觸有關(guān)的參數(shù)。,Step 2.啟動(dòng)接觸向?qū)?Main Menu Preprocessor Modeling Create Contact Pair,Step 3. 創(chuàng)建接觸對(duì),圖

28、4-3 圖4-4 【Pick Target】(圖43),拾取線(見(jiàn)圖35) 【OK】，再定義接觸面 返回圖44 【Next】 (圖45),圖45,【opening setting】圖46,圖46,Step 4. 設(shè)置接觸參數(shù)【OK】 返回圖45【create】,【Finish】返回圖42,退出 Step 5.求解 求解過(guò)程同前,圖42,5 實(shí)例 2 平面密封圈計(jì)算,驗(yàn)證為接觸剛度估計(jì)一個(gè)合適的初始值的重要性 驗(yàn)證接觸分析中摩擦的應(yīng)用 彈性庫(kù)侖摩擦 驗(yàn)證基于接觸單元臨界狀態(tài)變化的時(shí)間步預(yù)測(cè)如何會(huì)對(duì)計(jì)算效率有害 模型描述：2D 超彈平面密封圈 對(duì)稱，密封圈壓縮模擬,分析此模型采用兩種方法： 1.

29、不使用接觸向?qū)В⒚魑募?Seal.inp 2. 使用接觸向?qū)?chuàng)建接觸對(duì)，并用 GUI 方式創(chuàng)建 Pilot 節(jié)點(diǎn)。 本例具有： 幾何 N.L（大應(yīng)變與大變形） 材料 N.L（超彈） 接觸 N.L 數(shù)據(jù)文件: Seal.inp,Step 1. 建模：?jiǎn)卧愋?、?shí)常數(shù)、材料特性、基體分網(wǎng),/prep7 et,1,56 ! HYPER56 2D 4node U_P Hyperelastic Solid mp,nuxy,0.49 tb,mooney,1 tbdata,1,80 ! C10 = 80 tbdata,2,20 ! C01 = 20 k,1 $k,2,0.333,0 $k,3,0.86

30、7,0.867 k,4,1.1,0.867 k,5,1.1,1 $k,6,0.8,1 $k,7,0.267,0.133 k,8,0,0.133 l,1,2 *repeat,7,1,1 ! 將l,1,2命令重復(fù)7次 l,8,1 lfil,1,2,0.20 $lfil,2,3,0.15 $lfil,5,6,0.20 lfil,6,7,0.15 $lfil,7,8,0.05 $lfil,8,1,0.05 al,all ! 應(yīng)用所有選擇的線生成面 k,25,-0.6,0 $k,26,1.1,0 $k,27,-0.6,1.0 k,28,1.1,1.0 lstr,26,25 ! L9 lstr,27,28

31、 ! L10,圖5-2,! * 基體分網(wǎng) * lesize,8,2 $lesize,13,4 $lesize,14,4 esize,0.035 type,1 mat,1 amesh,all save,seal,db Step 2. 創(chuàng)建接觸對(duì)，先定義單元和接觸參數(shù) et,2,169 ! Target169 2D Target Segment et,3,171 ! Conta171 2D Surface to Surface Contact mp,mu,1,0 ! 用材料特性定義摩擦，本例先無(wú)摩擦計(jì)算然后再有摩擦計(jì)算，比較結(jié)果 r,1,0.1,0.1 ! 用實(shí)常數(shù)定義接觸高級(jí)選項(xiàng)對(duì)于彎曲為主的問(wèn)

32、題，采用接觸剛度FKN = 0.1作為初始估計(jì)值, FTOLN = 0.1 lagrange穿透比例系數(shù)（缺省值） r,2,0.1,0.1 ! 如果FKN = 1.0則不收斂,! * 創(chuàng)建接觸對(duì)1（定義實(shí)常數(shù)和分網(wǎng)） *,type,2 ! 目標(biāo)面 1 Target169 real,1 ! Target169的實(shí)常數(shù) mat,1 ! 接觸面 1 定義材料與超彈單元同 lesize,15,1 ! 目標(biāo)面1（L15）分網(wǎng)尺寸(圖53) lmesh,15 ! 目標(biāo)面1（L15）分網(wǎng),采用自動(dòng)分網(wǎng)，此剛性面自動(dòng)約束。檢查外法線方向 【OK】 ! * 接觸面1：實(shí)常數(shù)（與目標(biāo)面1一致）、分網(wǎng) lsel,s

33、,line,1,3,1 ! 選L1,L2,L3 lsel,a,line,9,10,1 ! 再加L9,L10 lsel,a,line,14 ! 再加L14 nsll,s,1 ! 選擇和所選線相連系的節(jié)點(diǎn) type,3 ! Conta171這里未發(fā)布Real,1，實(shí)常數(shù)未重新定義就是前面的Real,1，MAT也是前面的MAT, 1 esurf ! 生成接觸單元 lsel,all ! * 檢查外法線方向【OK】 *,圖53,! * 創(chuàng)建接觸對(duì)2 * type,2 ! 目標(biāo)面2 Target169 real,2 ! MAT未重新定義，就是前面的MAT, 1 lesize,16,1 ! 目標(biāo)面2 lme

34、sh,16 ! L16分網(wǎng)（目標(biāo)2分網(wǎng)） kmesh,27 ! 指定Pilot點(diǎn) ! * 檢查外法線方向【OK】 * type,3 ! 接觸面2 Conta171 Real,2； Mat,1 lsel,s,line,5,7, 1 ! 選L5, L6, L7 lsel,a,line,11,13,1 ! 再加L11,L12,L13 nsll,s,1 esurf ! 生成接觸單元,其Real,2； Mat,1 lsel,all nsel,all ! * 檢查外法線方向【OK】 *,! * 創(chuàng)建接觸對(duì)2 *,! 剛性面1 被約束 ! 剛性面2 隨Pilot點(diǎn)移動(dòng) ! L4對(duì)稱約束 dl,4,ux,0

35、n_load=node(kx(27),ky(27),0) ! 為定義剛性面2的位移做準(zhǔn)備 finish,Step 3. 施加邊界條件,/solu nlgeom,on solc,on ! N.L求解自動(dòng)控制打開(kāi)（缺省） time,0.85 ! 載荷步、結(jié)束時(shí)間 d,n_load,uy,-0.85 nsubst,25,500,10 outres,all,all monitor,var3,n_load,fy,Step 4.定義求解選項(xiàng)和載荷步,Step 5.求解 Solve Step 6.查看結(jié)果 /post1 pldsp,2 ! 變形圖 plnsol,s,eqv,0,1 ! Von Mises 應(yīng)

36、力云圖 plnsol,cont,pres,0, ! 接觸壓力圖 plnsol,epto,eqv ! 繪等效總應(yīng)變圖54 save,seal,db,定義反力變量、繪載荷變形圖,Utility Menu Plot Elements Main Menu Time Hist Postpro Define Variables 【Add】 Reaction force 【OK】 拾取節(jié)點(diǎn)263(控制點(diǎn)) 【OK】 User-Specified Load = FORCE: FY 【OK】 【close】 或命令： /POST 26 RFORCE,2,263,F,Y,FORCE Main Menu Time

37、Hist Postpro Graph Variables 1st Variable to gragh = 2 【OK】 或命令： PLVAR, 2 (圖54),圖54,2、用接觸向?qū)?chuàng)建接觸對(duì)，用 GUI 方式創(chuàng)建 Pilot 節(jié)點(diǎn)（剛性面控制點(diǎn)），然后求解。 Step 1. 恢復(fù)數(shù)據(jù)庫(kù)文件 Seal.db（包括基體的幾何、單元、分網(wǎng)；沒(méi)有選接觸單元與目標(biāo)單元，未定義接觸對(duì)） Utility Menu File Resume from 選 Seal.db 【OK】 Step 2. 啟動(dòng)接觸向?qū)?Main Menu Preprocessor Modeling Create Contact Pai

38、r 圖 42 【contact Wizard】,剛性目標(biāo)面 - 使用 Pilot 節(jié)點(diǎn), 創(chuàng)建接觸對(duì)1,圖 55,選線(1)( 接觸對(duì)1的目標(biāo)面)【OK】返回圖55【next】,圖56 選線(2) (接觸對(duì)1的接觸面)【OK】；返回圖56【next】,Step 3.設(shè)置接觸參數(shù),圖57,圖58 圖59,定義摩擦 (圖59) 【OK】返回圖 57 Coefficient friction 0.2 【Create】,511 圖512 接觸對(duì) 1 (圖511),圖510圖,Step 5.創(chuàng)建目標(biāo)單元控制點(diǎn),圖514,(1).設(shè)置單元屬性用命令流： Type,4 MAT,1 REAL,4 TSHAP,

39、PIL10 GUI: Main Menu Preprocessor Modeling Create Element Elem Attributes,(2).繪關(guān)鍵點(diǎn) Utility Menu Plot Keypoints Keypoints (3).建立目標(biāo)單元控制點(diǎn) Main Menu Preprocessor Meshing Mesh Keypoints 選 K27 【OK】 或命令流： KMESH,27 Step 6.施加位移約束 控制點(diǎn) 27 上施加 UY = -0.85 Main Menu Solution Define Loads Apply Structural Displace

40、ment On Keypoints 選 K27 【OK】 Lab2 DOFs to be constrained: UY Value: -0.5 【OK】,Step 5.創(chuàng)建目標(biāo)單元控制點(diǎn),施加對(duì)稱邊界條件： Main Menu Solution Define Loads Apply Structural Displacement Symmetry B.C. On Lines 選線L4 【OK】 Step 7.求解控制、求解 Main Menu Solution Analysis Type Soln Controls,Step 5.創(chuàng)建目標(biāo)單元控制點(diǎn),Main Menu Solution So

41、lve LS_Current Step 8.后處理 同1 無(wú)摩擦(MU = 0) Von Mises = 145.096 有摩擦(MU = 0.2) Von Mises = 142.038 接觸剛度取 1.0 不收斂 0.1 收斂 /POST26 Main Menu TimeHist Postpro 選控制點(diǎn)27， 繪圖(見(jiàn)圖5-4)。,Step 5.創(chuàng)建目標(biāo)單元控制點(diǎn),面一面接觸具有 20 個(gè)可用實(shí)常數(shù)，2 個(gè)材料屬性和 30 個(gè)可用單元選項(xiàng)設(shè)置。能夠模擬特殊的效果和處理困難的收斂情況。 通常分析開(kāi)始先用缺省值，只指定罰剛度和穿透容差及子步數(shù)。 只有在缺省設(shè)置遇到困難時(shí)才用高級(jí)選項(xiàng)。 所有的

42、單元選項(xiàng)和參數(shù)都可以通過(guò)接觸向?qū)?lái)控制；也可以通過(guò)實(shí)常數(shù)和單元選項(xiàng)來(lái)指定。 Conta178 接觸單元、實(shí)常數(shù)選項(xiàng)見(jiàn)圖 6-1 及 6-2。,圖6-1,6 面一面接觸單元高級(jí)接觸選項(xiàng),6 面一面接觸單元高級(jí)接觸選項(xiàng),圖 6-2,進(jìn)入接觸選項(xiàng)菜單用下列過(guò)程：,Main Menu Preprocessor Modeling Create Contact Pair 選 【property】 1、第 1 個(gè)重要選項(xiàng)：如何保證接觸協(xié)調(diào)性： 使用修正的拉朗日法（缺省） 純罰函數(shù)方法 在接觸向?qū)е?Basic Contact algorithm: Augmented Lagrange Method Pena

43、lty Method 對(duì)大多數(shù)模型，修正的拉格朗日方法能很好地工作。 罰函數(shù)法推薦用于具有變形很大的單元，很大的摩擦系數(shù)和/或用修正的拉格朗日方法時(shí)收斂性很差的情況。,1、最重要的選項(xiàng)是法向罰剛度或接觸剛度,對(duì)于大面積接觸采用起始值（因子）1.0，對(duì)于柔性接觸采用 0.1。 大值對(duì)應(yīng)較高精度；小值對(duì)應(yīng)較好的收斂性。 有時(shí)最好先以較小的接觸剛度進(jìn)行分析，然后在一系列載荷步中逐漸增大剛度“漸變”的接觸剛度，提高收斂性圖6-2 在最后的載荷步逐漸提高到一個(gè)剛硬的值將提高計(jì)算精度?？赏ㄟ^(guò)向?qū)У幕颈砘騿卧匦圆藛卧O(shè)置以允許程序更新接觸剛度 接觸向?qū)?Basic Normal Penalty Stiff

44、ness 1.0 factor constant 向?qū)?Basic contact stiffness update: Each load step (PAIR ID based) none Each load step (允許用戶指定剛度變化) Each substep (允許自動(dòng)和用戶指定變化),2、第 2 個(gè)重要選項(xiàng)穿透容差。,缺省情況下，穿透容差是一個(gè)因子乘以基體單元厚度。 對(duì)于變化很大的網(wǎng)格密度，采用因子會(huì)在接觸表面的某些部分產(chǎn)生太小的容差，這時(shí)采用絕對(duì)值可能更好。 不要使用太小的容差，因?yàn)樗偸菍?duì)收斂性有害。,3、對(duì)于臨界接觸狀態(tài)變化的自動(dòng)時(shí)間步控制,或在單元選項(xiàng)中控制： Cont

45、act time/load prediction K7: No prediction,1 不控制：不影響自動(dòng)時(shí)間步，對(duì)靜力問(wèn)題自動(dòng)時(shí)間步打開(kāi)時(shí)此選項(xiàng)一般是足夠的。2 自動(dòng)二分：如果接觸狀態(tài)變化明顯，時(shí)間步長(zhǎng)將二分。對(duì)于動(dòng)力問(wèn)題自動(dòng)二分通常是足夠的。3 合理值：比自動(dòng)細(xì)分更耗時(shí)的算法。4 最小值：此選項(xiàng)為下一步預(yù)測(cè)最小時(shí)間增量（很耗計(jì)算時(shí)間，不推薦）,4、Pinball 區(qū)域影響接觸狀態(tài)的確定和其它許多接觸特性,Pinball 區(qū)域是環(huán)繞接觸單元的園（2D）或球（3D），描述接觸單元周圍 “遠(yuǎn)” 和 “近” 區(qū)域的邊界 (圖6-3)。 在缺省情況下，Pinball 區(qū)域半徑是 4基體單元厚度（剛

46、柔）或 2基體單元厚度（柔一柔） 可以為 Pinball 半徑指定一個(gè)不同的值。,圖6-3,也可用實(shí)常數(shù)PINB調(diào)整球形區(qū)（對(duì)于初始值侵入大的問(wèn)題是必要的）。,5、幾種不同的接觸模式這些選項(xiàng)使你能夠模擬特殊的物理現(xiàn)象。,或單元選項(xiàng): Behavior contact surface K12: standard 這些選項(xiàng)包括： 標(biāo)準(zhǔn)：正常的接觸閉合和打開(kāi)行為，具有正常的粘著/滑動(dòng)摩擦行為。 粗糙：正常接觸閉合和打開(kāi)行為，但不發(fā)生滑動(dòng)（類似于具有無(wú)限摩擦系數(shù)） 不分離（滑動(dòng)）：目標(biāo)面和接觸面一旦建立接觸就不再分離（允許滑動(dòng)） 綁定：目標(biāo)面和接觸面一旦接觸就粘在一起 不分離（永遠(yuǎn)）：初始位于 Pin

47、ball 區(qū)域內(nèi)或已經(jīng)接觸的接觸檢查點(diǎn)在法向不分離（允許滑動(dòng)） 綁定接觸（永遠(yuǎn)）：初始位于 Pinball 區(qū)域內(nèi)或已經(jīng)接觸的接觸檢查點(diǎn)在剩余的分析過(guò)程中綁定在一起（Design Space 缺省值） 綁定接觸（初始接觸）：只在初始接觸的地方采用綁定，初始分開(kāi)的地方保持分開(kāi)。 計(jì)算實(shí)例：懸臂梁端部旋轉(zhuǎn)（不分離行為）,6、影響某些表面行為的選項(xiàng),Contact Opening stiffness（分開(kāi)時(shí)的間隙剛度）保證不分離和綁定行為，它通過(guò)使用當(dāng)存在間隙也具有非零剛度的彈簧來(lái)連接表面。 缺省情況下，此彈簧剛度等于法向罰剛度，其效果類似于法向罰剛度 剛度太小精度低；剛度太大會(huì)引起收斂問(wèn)題。,摩擦

48、系數(shù)影響基本摩擦行為：,Contact cohesion 表示當(dāng)沒(méi)有法向壓力時(shí)開(kāi)始滑動(dòng)的摩擦應(yīng)力值。 摩擦導(dǎo)致非對(duì)稱剛度陣。因?yàn)榉菍?duì)稱矩陣很難計(jì)算（因此導(dǎo)致求解變慢），程序自動(dòng)控制執(zhí)行對(duì)稱求解，利用此算法可以解決多數(shù)含摩擦接觸問(wèn)題。 有時(shí)，采用非對(duì)稱矩陣能獲用更好的收斂性。 如果遇到收斂緩慢問(wèn)題可以用不對(duì)稱求解選項(xiàng)。 記住：這種情況必須使用稀疏或波前求解器。 對(duì)于每個(gè)支持非對(duì)稱矩陣的單元，此選項(xiàng)也可以由下列菜單激活： Main Menu Solution Unabridged Menu Analysis Options 設(shè)置 Newton-Raphson 選項(xiàng)為 Full N_R unsymm

49、,7、初始穿透,有幾種技術(shù)可以模擬初始穿透接觸問(wèn)題（如過(guò)盈裝配）?？梢允褂贸跏紟缀未┩福蛑付ㄆ屏?，或二者皆有。(圖6-5) 指定偏移量（CNOF）,或在實(shí)常數(shù)中指定偏移量（CNOF）contact surface offset CNOF: 0.025,正的 CNOF 加大初始穿透 負(fù)的 CNOF 減小初始穿透或?qū)е麻g隙 CNOF 可與幾何穿透組合 自動(dòng) CNOF 調(diào)整 允許 ANSYS 基于初始穿透自動(dòng)給定 CNOF 值。 導(dǎo)致 “剛好接觸” 配置 ICONT 缺省為 0,或單元選項(xiàng): Auto CNOF/ICONT adjustment K5: No. Auto. Adjust,7、初始

50、穿透,初始穿透選項(xiàng)包括： Include everything：包括由幾何模型和指定偏移量 （如果有的話）引起的初始穿透 Exclude everything：忽略所有初始穿透效應(yīng)。 Include with ramped effects：漸變初始穿透，以提高收斂性。 Include offset only：只包括由偏移量指定的基本初始穿透。 Include offset only w/ramp：只包括由偏移量指定的基本初始穿透，且漸變初始穿透以提高收斂性,7、初始穿透,如果模型包含初始幾何穿透，接觸力將立即“階躍”到一個(gè)大值。 載荷突變經(jīng)常導(dǎo)致收斂困難，期望有一種機(jī)制能夠?qū)⒊跏即┩感?yīng)漸變到

51、零。 Include with ramped effects 和 Include off set only w/ramp選項(xiàng)通過(guò)在第一載荷步，將初始穿透漸變?yōu)榱憧朔諗坷щy。為求得好的結(jié)果，在第一載荷步不應(yīng)施加其它載荷 (圖6-6)。 計(jì)算實(shí)例：初始穿透,圖6-6,7、初始穿透,初始不接觸的兩個(gè)（或多個(gè)）物體的靜力分析中，在接觸建立前可能產(chǎn)生剛體運(yùn)動(dòng) (圖6-7)。 此例中圓柱體沒(méi)有施加位移約束，面由力控制。 圓柱體的約束由圓柱體和平板之間的接觸建立。 求解過(guò)程中兩個(gè)物體分離，剛度矩陣奇異。 ANSYS 將產(chǎn)生一個(gè)負(fù)主元警告。有幾個(gè)選項(xiàng)可以解決由于初始不相連物體引起的剛體模式：,圖6-7,8、

52、剛體模式,三個(gè)高級(jí)接觸特性允許調(diào)整初始接觸條件以防止剛體模式： (1).自動(dòng) CONF 調(diào)整 程序計(jì)算 CNOF 以清除間隙。 (2).初始接觸環(huán)（ICONT） 將調(diào)整帶內(nèi)接觸表面上的節(jié)點(diǎn)移到目標(biāo)面上 (3).初始允許穿透范圍 (PMIN Mat, 1 TYPE,2 MAT,2 REAL,2 LSEL,LOC,X,1.9 NSLL,1 ! 1.9線上的節(jié)點(diǎn) ESURF ! 創(chuàng)建目標(biāo)單元Target169 ! * TYPE,3 ! Mat, 2 ; Real, 2 LSEL,LOC,X,2 NSLL,S,1 ESURF ! 創(chuàng)建接觸面conta172 ALLSEL,ALL,/SOLUTION D

53、,NODE(1.5,0,0),UY,0 D,NODE(2.4,0,0),UY,0 NLGEOM,ON NSUBST,10,50,5 FINISH /PBC,U,1 /NUMBER,1 /PNUM,MAT,1 EPLOT SAVE,interference,db ! 包含有限元模型、接觸對(duì)、材料、邊界條件，求解控制,例3. 緊配合環(huán)（初始穿透）,分析步驟： Step 1.恢復(fù)數(shù)據(jù)庫(kù)文件 Utility Menu File Resume from 選 interference.db 【OK】 Step 2. 重定位輸出文件 Utility Menu File Switch Output to Fi

54、le Interference.output 【OK】 或命令： /Output,interference,Output Step 3.求解接觸分析 /Solu solve Step 4.重定位輸出到輸出窗口 Utility Menu File Switch Output to Output Window 或命令：/OUTPUT, TERM Step 5.繪制徑向應(yīng)力 Main Menu General Postproc Plot Results Contour Plot Nodal Solu Stress X-direction SX 【OK】,例3. 緊配合環(huán)（初始穿透）,或命令: /PO

55、ST1 PLNSOL,S,X ! 結(jié)果：徑向應(yīng)力0 Step 6.列表顯示接觸結(jié)果 Main Menu General Postproc List Results Nodal Solution 選 Nonlinear items contact CONT 【OK】 或命令： PRNSOL,CONT STAT = 0“張開(kāi)”幾乎不接觸 STAT = 1 “張開(kāi)”但幾乎接觸 STAT = 2 “閉合”并滑動(dòng) STAT = 3 “閉合”并粘著,例3. 緊配合環(huán)（初始穿透）,Step 7.檢查輸出文件內(nèi)容（interference.output） * NOTE * No contact was det

56、ected for a contact pair specified by real constant set 2 大初始穿透問(wèn)題可能導(dǎo)致無(wú)法檢測(cè)到接觸或檢測(cè)到接觸面求解存在收斂困難。在此情況下推薦在第一個(gè)載荷步中包含漸進(jìn)化效應(yīng)。Step 8.為初始穿透選項(xiàng)指定漸進(jìn)化效應(yīng) 可通過(guò)conta172單元選項(xiàng)設(shè)置 可通過(guò)接觸向?qū)е蠭nitial Adjustment的設(shè)置 (1) Main Menu Preprocessor Element Type Add/Edit/Delete 選 Type 3 conta172 【Option】 initial penetration/gap K9 = inc

57、lude_ramp 【OK】 【close】 或命令: /prep7 KEYOPT,3,9,2,例3. 緊配合環(huán)（初始穿透）,(2)或向?qū)?Initial Adjustment Initial penetration: Include everything with ramped effect 【OK】,例3. 緊配合環(huán)（初始穿透）,Step 9.重定位輸出到文件 Utility Menu File Switch Output to File Ramped.output 【OK】 Step 10. 分析 /Solu Solve Step 11.重定位輸出到輸出窗口 Utility Menu F

58、ile Switch Output to Output Window 命令：/output,TERM Step 12.繪徑向應(yīng)力 Main Menu General Postproc Plot Results Contour Plot Nodal Solu Stress X-direction SX 【OK】 或命令: /POST 1 PLNSOL,S,X,例3. 緊配合環(huán)（初始穿透）,Step 13.將軸對(duì)稱模型擴(kuò)展為整個(gè)模型顯示 Utility Menu PlotCtrls Style Symmetry Expansion 2D Axi-Symmetric Select expansion amount = Full expansion 【OK】 或命令： /EXPAND,36,Axis,10 /Repxzlot ISO 顯示 /View,1,1,1 /Replot Step 14. 檢查