在ANSYS中正確地模擬過盈配合 (1) 轉貼
2017-03-02 by:CAE仿真在線 來源:互聯(lián)網
過盈配合在機械產品的裝配中使用的相當普遍。比如軸與軸承、軸與軸瓦、汽車的制動盤等,都是通過一定的過盈量來使兩個裝配部件緊密地連接起來。
下面討論如何在 ANSYS 中正確地模擬過盈配合。
過盈配合在有限元分析中是一種典型的非線性接觸行為。在有限元分析中設定了接觸,從本質上來講就是對相互接觸的兩個部件施加了某種約束,不同的接觸算法對于接觸約束的處理方法有所不同。接觸約束的理論算法的選擇,在 ANSYS 中是通過設置 contact 單元的 KEOPT(2) 選項來實現的。
在 ANSYS 中目前主要有 5 種接觸約束算法:
KEYOPT(2)=0 Augmented Lagrangian - 加強的拉格朗日算法, 這是 ANSYS 的缺省選擇;
KEYOPT(2)=1 Penalty function - 罰函算法;
KEYOPT(2)=2 Multipoint constraint (MPC) - 多點約束算法;
KEYOPT(2)=3 Lagrange multiplier on contact normal and penalty on tangent -
接觸法向采用拉格朗日乘子,接觸切向采用罰函數的綜合算法。
KEYOPT(2)=4 Pure Lagrange multiplier on contact normal and tangent - 法向和切向均采用拉格朗日乘子算法。
各種不同的約束算法各有其優(yōu)缺點,各有各自最適用的場合,具體情況需要具體對待。大部分情況下,默認選擇KEYOPT(2)=0 就夠用了。
過盈配合所致的接觸分析的難點在于如何確定初始接觸狀態(tài)。初始接觸狀態(tài)設置得不對,會導致錯誤的計算結果或者不準確的計算結果,下面舉兩個例子來說明。
例1.兩個圓柱體在幾何上是剛好接觸,劃分網格后有限元模型有間隙。如圖 1 所示。
這兩個圓柱體,在幾何上是剛好相切的,即處于幾何上剛好接觸的初始狀態(tài)。劃分網格后,由于在圓周上用小段直線代替了弧線,兩個圓柱體之間產生了一定的間隙,兩個圓柱體的有限元模型的初始狀態(tài)不再是接觸的。此時,如果接觸參數設置不當,就會因為初始約束不足,圓柱體出現剛體位移,得到錯誤的結果。
(說明:例 1 本來與設置過盈量是無關的,為了說明初始接觸狀態(tài)的重要性順帶說說。)
例2.有的人把兩個接觸部件的幾何位置設定一定的過盈量,想用這個過盈量來模擬過盈配合,這種做法是錯誤的,幾何上的過盈量不等于劃分網格后有限元模型的實際過盈量。
下面的圖 2 中,是一個孔類零件和一個軸類零件的截面圖,軸和孔在幾何位置上預設了過盈量。(內圈的紅色圓是孔邊界,外圈的藍色圓是軸邊界,軸和孔在幾何上是相互侵入的)。
在幾何上,圖 2 的軸和孔有一定的過盈配合量,其大小等于兩個圓的半徑之差,我們的本意是想用這個幾何位置上的過盈量來模擬過盈配合。
不幸的是,兩個部件劃分網格之后,實際的過盈量應該為單元之間的距離,即圖中靠得比較近的兩條線段之間的距離,顯然,這個距離不再等于我們預先設置的過盈量了。更何況,上面這個圖還是兩個部件的網格對得比較整齊的情況,如果網格對得不整齊,過盈量就和我們預設的差得更遠了。對于過盈配合來講,過盈量的數值變化對于過盈產生的應力的影響是很大的。
在 ANSYS 中,要正確的設置過盈配合,主要分 3 步:
第一步:設置 KEYOPT(9) = 4
KEYOPT(9) 的默認值為 0,意思是既考慮兩個接觸部件由于初始幾何位置造成的初始侵入量 (或者間隙),同時也考慮 CNOF 參數設置的偏移量。意即接觸部件的初始接觸狀態(tài)是由 CNOF 和初始侵入量 (或間隙) 共同決定的。在這種情況下,兩個接觸部件的初始幾何位置對初始接觸狀態(tài)是有影響的,這對于準確設置過盈量是很不利的。(前面例 2 已經說明了通過幾何位置設置初始過盈量是不準確的)。
在設置了 KEOPT(9) = 4 之后,程序在計算初始接觸狀態(tài)的時候就只考慮 CNOF 的設置值,不考慮接觸部件的幾何位置造成的侵入或間隙,而且過盈量是以 ramp 方式施加的。(ramp 施加方式即逐步施加)。
第二步:通過設置 Icont 實常數
劃分網格后,通常情況下,Target surface 和 contact surface 上的單元之間會有間隙或者過盈量,如果間隙或者過盈量在 Icont 設定誤差范圍內,間隙或者過盈量會被消除掉,程序會使 contact surface 和 target surface 上的單元處于剛好接觸的狀態(tài)。這個值的的具體設置可以參見幫助文檔,本文中設置為 0.2。
第三步:通過設置實常數 CNOF 來設置過盈量
在第二步中,通過 Icont 的設置,已經使得 Contact surface 上的單元和 Target surface 上面的單元處于剛好接觸的位置了,此時再設置 CNOF,CNOF 的值就是過盈量。(CNOF 的本意并不是過盈量,只是在有了前面的設定后,它的值就是我們所要的過盈量,其具體含義請參考 ANSYS 的幫助文檔)。
下面這個例子實際上是一本 ANSYS 書上的一個例子,這個例子的 PDF 版本在網上流傳甚廣,但是原書上的分析結果是錯誤的,具體錯誤之處,將在后面提及。
例3.一個簡單的軸和帶孔圓盤的過盈配合的實例。
圓盤的基本尺寸為:
內徑 Rpin = 35 mm (原書中此值為 34 mm),外徑 Rpout = 100 mm,盤高 Hp = 25mm;
軸的基本尺寸為:
內徑 Rain = 25 mm,外徑 Raout = 35 mm,軸長 La = 150 mm。
(原書中圓盤孔內徑為 Rpin = 34 mm,和軸在幾何上形成 1 mm 的過盈量。
由于結構是完全軸對稱的,故可只取四分之一模型分析之。
本例分析中,取過盈量 f = 0.01 mm,而且本例僅僅計算由于過盈配合所產生的應力。
按照本例各個物理量所取的單位,最終的計算結果中,應力單位應該為 MPa;
完整的命令流如下:
Finish
/clear,start
/TITLE,Contact analysis with initial interference
/PREP7
! 帶孔圓盤的基本尺寸;
Rpin=35
Rpout=100
Hp=25
! 軸的基本尺寸;
Rain=25
Raout=35
La=150
! 過盈量 f;
f=0.01
! 實體的單元類型為帶中間節(jié)點的2階六面體單元;
ET,1,solid186
MP,EX,1,2.1E5! 彈性模量;
MP,PRXY,1,0.3! poisson系數;
! 生成帶孔圓盤的1/4實體模型;
CYL4,0,0,Raout,0,Rpout,90,Hp
! 軸的 1/4 實體模型;
CYL4,0,0,Rain,0,Raout,90,La
! 把軸的位置沿著軸向移動一段距離;
VGEN, ,2, , , , ,-10, , ,1
!
!
!
!************對實體劃分網格 ****************************
!*
LESIZE,17, , ,15, , , , ,1
LESIZE,19, , ,15, , , , ,1
!*
LESIZE,18, , ,2, , , , ,1
LESIZE,20, , ,2, , , , ,1
!*
LESIZE,22, , ,20, , , , ,1
!*
LESIZE,5, , ,10, , , , ,1
LESIZE,7, , ,10, , , , ,1
!*
LESIZE,6, , ,8, , , , ,1
LESIZE,8, , ,8, , , , ,1
!*
LESIZE,10, , ,3, , , , ,1
!*
VSWEEP,ALL
!*************網格劃分完畢**********************************
!
!
/COM, CONTACT PAIR CREATION - START
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,0.2,0.2,0.9,0! Icont 系數設置為 0.2;
! 設置過盈量為 f;
RMORE,,,1.0E20,f,1.0,0
RMORE,0.0,0,1.0,,1.0,0
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,0
NROPT,UNSYM
KEYOPT,3,7,0
KEYOPT,3,8,0
!KEYOPT(9)=4;
! 不考慮初始幾何位置造成的過盈或者間隙,
! 只考慮 CNOF 參數設置的值,即過盈量。
KEYOPT,3,9,4
KEYOPT,3,10,2
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,9
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
/COM, CONTACT PAIR CREATION - END
FINISH
/SOL
FLST,2,4,5,ORDE,4
FITEM,2,5
FITEM,2,-6
FITEM,2,11
FITEM,2,-12
DA,P51X,SYMM! 設置對稱約束;
FLST,2,1,5,ORDE,1
FITEM,2,3
/GO
DA,P51X,ALL,
!*
ANTYPE,0
NLGEOM,1! 必須打開大變形效應;
NSUBST,1,0,0
AUTOTS,0
TIME,100
/STATUS,SOLU
SOLVE
FINISH
/POST1
SET,1,LAST,1,
PLNSOL,s,eqv! 查看單元的Von Mises stress;
esel,s,ename,,174
PLNSOL, CONT,PRES, 0,1.0! 在contact 單元上查看contact pressure;
SAVE
計算結果:
說明:修改 f 值即可修改過盈量。如果過盈量設置的過大,接觸部件可能會產生塑性變形,如果要考慮塑性變形,則應該修改材料屬性定義為塑性材料。
原書的錯誤在于:
附:原書上的命令流;
錯誤 1:試圖用幾何初始過盈量來模擬過盈配合,接觸部件在幾何位置上的過盈量不等于接觸部件劃分網格后實際的過盈量。
錯誤 2:他的所謂過盈量設置為 1 mm,這個值很不符合工程常理。對于配合尺寸只有35mm的軸和孔的過盈配合,這個值太大了,大得遠遠不符合工程常理,他的計算結果也可以證明這一點:按照那本書上的命令流運行計算,得到的最大接觸應力有 2000 多 MPa,最大 Von Mises 應力大約有 6000-7000 MPa,這個時候鋼早就屈服并產生塑形變形了,根本不可能產生這么大應力值,更何況,他的分析是在線彈性材料的基礎上進行分析的。
錯誤 3:他的 Contact 單元選用的是 Conta174,而實體單元卻選擇的是不帶中間節(jié)點的一階單元 SOLID185,這兩種單元的匹配是不對的。實體單元應該選擇帶中間節(jié)點的 SOLID186 單元。
下面是幫助文檔中關于 Conta174 單元的一段話:
The element is defined by eight nodes (the underlying solid or shell element has midside nodes)……
……If the underlying solid or shell elements do not have midside nodes, use CONTA173 (you may still use CONTA174 but you must drop all midside nodes).
大意為:如果接觸單元所依賴的實體單元 (或者 Shell 單元) 是帶中間節(jié)點的單元,接觸單元應該用 Conta174。如果實體單元 (或者 Shell 單元) 不帶中間節(jié)點,應該用 Conta173,如果你堅持使用 Conta174,則你必須 drop 所有的中間節(jié)點。
原書使用 SOLID185 計算出來的應力云圖分布規(guī)律是不對的 (且不說數值對否)。理由:這個結構是完全軸對稱的,在同一個高度上,沿著圓周方向各點的接觸應力值應該是一樣的,接觸應力值只能沿著軸向有變化。而他的接觸應力云圖分布在圓周方向是大,小,大,小…循環(huán)變化的,這顯然是錯的。
附:原書上的命令流;
/TITLE,Analysis of a Axis Contacting a hole in a Disc ! 定義標題
/PREP7
!*
ET,1,SOLID185 ! 定義單元類型
!*
MP,EX,1,2.1E5 ! 定義材料屬性
MP,PRXY,1,0.3
!*
CYL4,0,0,34,0,100,90,25 ! 創(chuàng)建四分之一圓環(huán)
CYL4,0,0,25,0,35,90,150
VGEN, ,2, , , , ,-10, , ,1 ! 移動軸
!*
LESIZE,17, , ,15, , , , ,1 ! 定義線的分網尺寸
LESIZE,19, , ,15, , , , ,1
!*
LESIZE,18, , ,2, , , , ,1
LESIZE,20, , ,2, , , , ,1
!*
LESIZE,22, , ,20, , , , ,1
!*
LESIZE,5, , ,10, , , , ,1
LESIZE,7, , ,10, , , , ,1
!*
LESIZE,6, , ,8, , , , ,1
LESIZE,8, , ,8, , , , ,1
!*
LESIZE,10, , ,3, , , , ,1
!*
VSWEEP,ALL ! 用掃掠方式對創(chuàng)建的體進行網格劃分
!*
/COM, CONTACT PAIR CREATION - START
MP,MU,1,0.2 ! 定義接觸摩擦系數
MAT,1
R,3 ! 定義接觸實常數
REAL,3
ET,2,170 ! 定義接觸單元類型
ET,3,174
R,3,,,0.1,0.1,,
NROPT,UNSYM
!* Generate the target surface 下面創(chuàng)建目標面
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
!* Generate the contact surface 下面創(chuàng)建接觸面
ASEL,S,,,9
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL
CMDEL,_TARGET
CMDEL,_CONTACT
ALLSEL,ALL
EPLOT
FINISH
!*
/SOLU ! 進入求解器
DA,5,SYMM ! 定義面的對稱位移邊條
DA,6,SYMM
DA,11,SYMM
DA,12,SYMM
DA,3,ALL, ! 定義面的位移約束條件
!*
ANTYPE,0 ! 指定分析類型為靜力分析
NLGEOM,1 ! 考慮大變形影響
AUTOTS,0
TIME,100
SOLVE !求解第一載荷步
!*
NSUBST,150,10000,10
OUTRES,ALL,ALL
AUTOTS,1
TIME,250
NSEL,S,LOC,Z,140 ! 選定軸向坐標為140的所有節(jié)點
D,ALL,UZ,40
ALLSEL,ALL
SOLVE ! 求解第二載荷步
!*
/EXPAND,4,POLAR,HALF,,90 !進行模型擴展
/REPLOT
!*
/POST1 ! 進入通用后處理器
SET,1,LAST,1, ! 指定查看的載荷步
PLNSOL,S,EQV,0,1 ! 查看等效應力的云圖
!*
SET, , ,1, ,120, ,
ESEL,S,ENAME,,174
EPLOT
PLNSOL,CONT,PRES,0,1
!*
PLNS,S,EQV
ANDATA,0.5, ,1,0,0,1,1,1 ! 查看動畫顯示
!*
/POST26
!*
RFORCE,2,925,F,Z, FZ_2 ! 定義約束反力變量
PLVAR,2, ! 繪制變量-時間曲線
FINISH
開放分享:優(yōu)質有限元技術文章,助你自學成才
相關標簽搜索:在ANSYS中正確地模擬過盈配合 (1) 轉貼 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓
編輯