abaqus第七講：abaqus接觸分析解讀課件

第七講接觸非線性

王慎平北京怡格明思軟件技術(shù)有限公司

第七講接觸非線性王慎平北京怡格明思軟1概述簡介主/從方法和隱含假設(shè)定義雙面接觸概述定義面的規(guī)則局部表面行為接觸問題中點的相對滑動接觸問題中調(diào)整初始節(jié)點位置接觸輸出完全約束的接觸行為剛體概述簡介2簡介什么是接觸?當(dāng)兩個實體接觸時，力通過它們的接觸表面?zhèn)鬟f。在某些情況下，只傳遞垂直接觸表面的力。如果存在摩擦，沿接觸表面的切向傳遞切向力。一般目標(biāo)：確定接觸面積和傳遞的應(yīng)力。接觸是嚴(yán)重不連續(xù)形式的非線性行為,是一類特殊的不連續(xù)約束。或者施加約束（表面不可以互相穿透），或者忽略約束。簡介什么是接觸?3為什么要定義接觸？除非用戶指定可能會接觸的表面和/或節(jié)點，ABAQUS現(xiàn)在還不能檢測接觸。為什么要定義接觸？4“DetroitEdisonpipewhipexperiment,”ABAQUS例子手冊2.1.2節(jié)接觸例子-間隙接觸用節(jié)點與節(jié)點接觸的方式為點接觸建模?！癉etroitEdisonpipewhipexpe5 赫茲接觸接觸面之間的相對位移很小。分布表面的接觸。例子選自“Coolantmanifoldcovergasketedjoint,”ABAQUS例子手冊5.1.3節(jié) 赫茲接觸例子選自“Coolantmanifoldco6 變形體之間的有限滑動接觸這是最一般類型的接觸。例子：螺紋連接。這些問題一般包含初始的過盈配合由于過盈精度引起的接觸壓力的分布 變形體之間的有限滑動接觸由于過盈精度引起的接觸壓力的分布7

自接觸單個表面與它自身接觸稱為自接觸。在ABAQUS的二維和三維模型中可用。在分析過中，當(dāng)表面嚴(yán)重變形時，使用自接觸將非常方便。對于某些問題，在分析之前不可能確定單個的接觸區(qū)域，或者確定接觸區(qū)域是非常困難的。把單個接觸表面作為接觸對定義自接觸，而不是通過兩個不同的表面定義。最小主應(yīng)力云圖SURF1(剛體)SURF2例子：橡膠墊片的壓縮(選自“Self-contactinrubber/foamcomponents:rubbergasket,”ABAQUS例子手冊1.1.17節(jié)) 自接觸最小主應(yīng)力云圖SURF1SURF2例子：橡膠墊片的壓8例子：金屬成型模擬變形體與剛體的接觸:表面之間的有限滑動（大位移）。變形組件之間的有限應(yīng)變。典型例子：成型模擬(剛體砧/模具、可變形組件。)例子：金屬成型模擬變形體與剛體的接觸:典型例子：9使用剛體表面減少計算時間

對于兩個相互接觸的物體，如果其中一個物體比另外的物體剛硬許多，可以將較為剛硬的物體指定為剛體；如，金屬成型過程中的砧。ABAQUS不需要求解剛體的變形。在參考點，最多利用六個自由度就可以模擬剛體的運動。有三種方法，可以以幾何的方式定義剛體表面：定義解析剛體表面。使用單元類型組合定義剛體（包括剛體單元），并將物體聲明為剛體（一般剛體功能）。寫用戶子程序(RSURFU)定義表面(只有ABAQUS/Standard可用)。使用剛體表面減少計算時間10主/從方法和隱含假設(shè)主/從方法和隱含假設(shè)11因為接觸表面下面的實體被離散化，接觸表面也必須被離散化。接觸表面分為:主面從面因為接觸表面下面的實體被離散化，接觸表面也必須被離散化。12ABAQUS/Standard中使用的嚴(yán)格的主/從接觸算法。從屬表面的節(jié)點不能穿透到主控表面。主控表面的節(jié)點可以穿透到從屬表面。ABAQUS/Standard中使用的嚴(yán)格的主/從接觸算法。13在多數(shù)情況下，ABAQUS/Explicit使用平衡的主/從算法。兩次應(yīng)用單純的主-從接觸算法，并平均化。在第二次施加接觸約束時，交換主/從表面。使接觸物體之間的侵徹達(dá)到最小化。在多數(shù)情況下，ABAQUS/Explicit使用平衡的主/從14接觸方向總是與主控表面垂直。沿主控表面的法向檢查接觸條件。沿主控表面的法向傳遞接觸力。沿接觸表面的切向傳遞摩擦力。接觸方向總是與主控表面垂直。15定義雙面接觸概述定義雙面接觸概述16定義接觸的三個步驟：

定義接觸表面。定義接觸對。定義接觸屬性。123定義接觸的三個步驟：12317輸入文件中，完整定義接觸句法的例子：*SURFACE,NAME=ASURFSLIDER,S1*SURFACE,NAME=BSURFBLOCK,S3*CONTACTPAIR,INTERACTION=FRIC1ASURF,BSURF*SURFACEINTERACTION,NAME=FRIC11.0,*FRICTION0.4,輸入文件中，完整定義接觸句法的例子：18定義表面利用*SURFACE選項定義表面。利用每個單元集的表面標(biāo)識符指定面。可以用單元集的名字或單元號指定表面。*SURFACE,NAME=ASURFSLIDER,S1*SURFACE,NAME=BSURFBLOCK,S3接觸發(fā)生在單元集SLIDER的底部(S1)接觸發(fā)生在單元集的BLOCK頂部(S3)定義表面接觸發(fā)生在單元集SLIDER的底部(S1)接觸發(fā)生19定義表面接觸屬性

*SURFACEINTERACTION選項定義表面接觸屬性。定義表面行為屬性，比如摩擦。對于二維問題，定義接觸面在垂直平面方向的厚度。定義表面接觸屬性20*CONTACTPAIR,INTERACTION=FRIC1ASURF,BSURF*SURFACEINTERACTION,NAME=FRIC11.0,*FRICTION0.4,*SURFACEINTERACTION的子選項，列出表面的基本屬性接觸面在垂直平面方向的厚度*SURFACEINTERACTION的子選項，列出表面的21定義接觸對定義表面和接觸屬性之后，就可以定義“接觸對”。在分析過程中，接觸對指定了兩個相互接觸的表面。在ABAQUS/Standard中，第一個表面為從屬表面，第二個表面為主控表面。在ABAQUS/Explicit中，表面的順序一般不影響接觸計算。*CONTACTPAIR,INTERACTION=FRIC1ASURF,BSURF定義接觸對*CONTACTPAIR,INTERACTI22定義面的規(guī)則定義面的規(guī)則23定義面的規(guī)則在變形體或劃分網(wǎng)格的剛體上，使用*SURFACE,TYPE=ELEMENT選項定義表面。通過指定單元表面標(biāo)識符定義表面。讓ABAQUS自動確定連續(xù)體單元的“自由表面”。使用*SURFACE,TYPE=[SEGMENTS|CYLINDER|REVOLUTION]選項和*RIGIDBODY選項定義解析剛體表面。使用*SURFACE,TYPE=NODE選項，定義可能接觸的節(jié)點表面。定義面的規(guī)則在變形體或劃分網(wǎng)格的剛體上，使用*SURFAC24在實體單元上定義表面使用表面標(biāo)識符例子：4-節(jié)點四面體單元(CPE4,CAX4,等等)*SURFACE,NAME=EXAMPLE11,S41,S12,S12,S2...在實體單元上定義表面25使用自動的表面定義*SURFACE,NAME=EXAMPLE2ELSET1,沒有面標(biāo)識符使用自動的表面定義沒有面標(biāo)識符26在結(jié)構(gòu)單元上定義表面(殼、薄膜、剛體、梁)結(jié)構(gòu)單元的法向表示預(yù)期的接觸方向。法向基于單元局部的節(jié)點編號。正的法線方向在=SPOS表面。負(fù)的法線方向在=SNEG表面。殼和薄膜(S4R,S8R,M3D4,等)2-D桁架和梁(B21,T2D2,等)在結(jié)構(gòu)單元上定義表面(殼、薄膜、剛體、梁)殼和薄膜(S27在所定義的表面集之內(nèi)，表面法向的定義必需一致。*ELEMENT,TYPE=B21,ELSET=BOTTOM10,1,211,2,312,3,6*ELEMENT,TYPE=B21,ELSET=TOP20,4,521,5,6*ELSET,ELSET=BEAMSBOTTOM,TOP*SURFACE,NAME=SURF1BOTTOM,SPOSTOP,SNEG*SURFACE,NAME=SURF1BEAMS,SPOS在所定義的表面集之內(nèi)，表面法向的定義必需一致。*SURFAC28解析剛體表面所有的剛體表面將被作為主控表面；剛體表面的自由度與剛體的參考節(jié)點相關(guān)聯(lián)，而不是表面上的點。使用*SURFACE選項，有三種可用的解析表面：使用TYPE=SEGMENTS定義二維剛體表面。使用TYPE=CYLINDER定義三維剛體表面，在垂直于平面方向無限延伸。使用TYPE=REVOLUTION定義旋轉(zhuǎn)的三維表面。解析剛體表面29基于節(jié)點的表面與指定單元表面定義接觸表面的方法不同，基于節(jié)點的表面只包含節(jié)點?；诠?jié)點的表面一般被用作從屬表面?；诠?jié)點的表面30基于節(jié)點表面的例子：定義包含接觸節(jié)點的表面前面定義的表面線：基于節(jié)點的表面球：基于單元的表面*SURFACE,TYPE=NODE,NAME=STRINGSSTRINGS,*CONTACTPAIR,INTERACTION=SMOOTHSTRINGS,BALL基于節(jié)點表面的例子：定義包含接觸節(jié)點的表面前面定義的表面線：31表面定義的限制接觸表面下所有的單元必須協(xié)調(diào)。它們必須：維數(shù)相同（二維或三維）。對于二維表面：所有單元為平面的或軸對稱的（不可以兩種都有）。插值階數(shù)相同（一階或二階）。所有單元為可變形的或所有單元為剛體（不可以兩種都有）。對于ABAQUS/Standard中的主控表面和ABAQUS/Explicit中的所有表面，還有附加的限制：在不離開表面、穿過表面或通過單個點穿越表面的前提下，必須可以跨過表面中的兩個任意的點。表面定義的限制32對于三維情況，限制的例子。需要穿過或離開表面才能跨過。合法的、封閉的3-D主表面合法的、簡單連接的3-D主表面不合法的3-D主表面不合法的3-D主表面不可以通過一個點跨過對于三維情況，限制的例子。需要穿過或離開表面才能跨過。合法的33局部表面行為局部表面行為34一般的接觸建模包括不同方向的接觸相互作用：主控表面的法向表面的切向一般的接觸建模包括不同方向的接觸相互作用：35接觸法向行為硬接觸在所有接觸問題中，“硬”接觸是默認(rèn)的局部行為。硬接觸可以通過經(jīng)典的拉格朗日乘子方法實施，或者通過增廣拉格朗日方法實施：

壓力－間距關(guān)系*surfaceinteraction,name=...*surfacebehavior,augmentedlagrange接觸法向行為壓力－間距關(guān)系*surfaceinteract36除硬接觸外，其它的方法作為*SURFACEINTERACTION

選項的子選項*SURFACE

BEHAVIOR選項用于指定：軟接觸（指數(shù)或表格方式

表達(dá)的壓力-間距關(guān)系）沒有分離的接觸其它選項：間距相關(guān)的粘性阻尼(*CONTACTDAMPING)。帶有過盈或拉伸接觸力的接觸(*CONTACTCONTROLS;只有ABAQUS/Standard可用)。除硬接觸外，其它的方法37接觸的切向行為兩個實體之間的接觸面可能存在摩擦剪切應(yīng)力。如果剪切應(yīng)力達(dá)到某個臨界值，實體之間可能會發(fā)生相對滑動；否則它們粘在一起。接觸的切向行為38摩擦可以引起高度的非線性效應(yīng)。很難獲得解。除非物理上必要，盡量不要使用。摩擦是非保守的。在ABAQUS/Standard中，摩擦將導(dǎo)致系統(tǒng)方程的不對稱。對于較大的摩擦系數(shù)(

>0.2)，*STEP,UNSYMM=YES選項將自動啟用。使用UNSYMM=NO將減緩收斂速度，但是解將是正確的（如果得到）。它可能會使用更少的磁盤空間。在ABAQUS/Explicit中，過大的摩擦系數(shù)將不會引起類似的問題，因為不需要求解系統(tǒng)方程組。摩擦可以引起高度的非線性效應(yīng)。39默認(rèn)情況下，ABAQUS使用庫倫摩擦模型。臨界摩擦應(yīng)力取決于接觸壓力：cr=p?；揪浞ǎ?FRICTION摩擦系數(shù)可以是相對滑動速度、壓力、溫度和場變量的函數(shù)。出于計算的考慮，ABAQUS/Standard中默認(rèn)的摩擦模型對理想行為作一定的近似：在達(dá)到不可恢復(fù)滑動發(fā)生之前，允許小量的彈性滑動。默認(rèn)情況下，ABAQUS使用庫倫摩擦模型。40可以定義混合的靜動摩擦模型。作為滑動速度函數(shù)的摩擦系數(shù)m，從ms（靜摩擦系數(shù)）到mk（動摩擦系數(shù)）以指數(shù)方式衰退。還有其它可用的摩擦行為，包括用戶自定義的摩擦模型。：

(在ABAQUS/Standard中為FRIC在ABAQUS/Explicit中為VFRIC)?？梢远x混合的靜動摩擦模型。41接觸問題中點的相對滑動接觸問題中點的相對滑動42兩種滑動距離選項：有限滑動：默認(rèn)為有限滑動，有限滑動是最通用的滑動方式。 在接觸表面之間，允許任意大的滑動和旋轉(zhuǎn)。小滑動： 在接觸表面之間，允許小的相對滑動。 只要接觸表面之間沒有大的相對移動，允許接觸表面

之間有大的轉(zhuǎn)動。

相比有限滑動，具有較小的計算費用。兩種滑動距離選項：43有限滑動在有限滑動問題中的主控表面，最好使用延伸到角點之外的定義。它們將防止從屬節(jié)點“跌落”到主表面的后面或被主表面的后面捕獲。對于ABAQUS/Standard中的接觸模擬，這是尤為重要的。不合適的定義(修整過的)合適的(未修整的)主控表面定義有限滑動不合適的定義(修整過的)合適的(未修整的)主控表44小滑動*CONTACTPAIR,INTERACTION=FRIC1,SMALLSLIDING

小滑動問題不需要通用的有限滑動算法。在整個分析過程中，從屬節(jié)點與主控表面上預(yù)先確定的、固定數(shù)量的節(jié)點相互作用。只是在分析開始時進(jìn)行接觸搜尋，所以不需要復(fù)雜的有限滑動搜尋算法。優(yōu)點：相對有限滑動，有更低的計算費用，尤其是三維問題。小滑動45小滑動接觸約束在分析開始時，ABAQUS將確定哪些節(jié)點是相互接觸的。右圖為小滑動選項基礎(chǔ)之上，節(jié)點101和它的主控表面BSURF之間可能的接觸演化方式。小滑動接觸約束46主控表面上，以內(nèi)插值方式得到節(jié)點201和202之間部分的法向量指向從屬節(jié)點101。在分析過程中，節(jié)點101不可以穿過垂直于法向量的平面。因為接觸線沒有端點，所以從屬節(jié)點不可能“滑出”接觸線。在分析過程中，如果使用NLGEOM選項，平面的方向?qū)⒒诠?jié)點201和202的變形而更新。主控表面上，以內(nèi)插值方式得到節(jié)點201和202之間部分的法向47接觸問題中調(diào)整初始節(jié)點位置接觸問題中調(diào)整初始節(jié)點位置48在分析之前，可以調(diào)整接觸面中節(jié)點的初始位置，并且不產(chǎn)生任何的應(yīng)力和應(yīng)變。ABAQUS/Standard允許用戶通過指定絕對距離或節(jié)點集的方式調(diào)整節(jié)點位置。*CONTACTPAIR,ADJUST=[distance|nodesetname]在分析之前，可以調(diào)整接觸面中節(jié)點的初始位置，并且不產(chǎn)生任何的49指定需要調(diào)整的絕對距離：*CONTACTPAIR,INTERACTION=FRIC1,ADJUST=a在調(diào)整范圍內(nèi)，初始張開的從屬節(jié)點將被移到主控表面上。調(diào)整距離沿主控表面的法向度量。所有初始過盈的從屬節(jié)點被重新放在主控表面上。指定需要調(diào)整的絕對距離：50指定需要調(diào)整的從屬節(jié)點集：在分析過程中，不在節(jié)點集

中并帶有初始過盈的從屬節(jié)點

將保持過盈，并產(chǎn)生應(yīng)變。

這些應(yīng)變將被施加在接觸面上。

指定

ADJUSTnodeset的例子：*NSET,NSET=CONNODE,GENERATE1,8,1*CONTACTPAIR,INTERACTION=RIG,ADJUST=CONNODE指定需要調(diào)整的從屬節(jié)點集：51ABAQUS/Explicit不允許接觸表面的初始過盈。在分析之前，接觸面上的節(jié)點將被自動調(diào)整，以刪除任何的初始過盈。在隨后的分析步中，這樣的調(diào)整將引起應(yīng)變。大量的調(diào)整將引起初始單元形狀的嚴(yán)重扭曲。如果發(fā)現(xiàn)這樣的錯誤信息，最好需要運行datacheck分析，并利用ABAQUS/Viewer查找問題。ABAQUS/Explicit不允許接觸表面的初始過盈。52接觸輸出接觸輸出53在使用*CONTACTOUTPUT、*CONTACTFILE或*CONTACTPRINT選項之后，所有接觸表面的接觸信息都是可用的。對于ABAQUS/Standard和ABAQUS/Explicit：接觸應(yīng)力：CSTRESS(接觸壓力CPRESS和摩擦剪切應(yīng)力CSHEAR1和CSHEAR2)對于ABAQUS/Standard，還可以得到：接觸位移：CDISP(接觸張開距離COPEN和相對切向運動CSLIP1和CSLIP2)對于ABAQUS/Explicit，還可以得到：滑動速度：FSLIPR累計的滑動距離：FSLIP在使用*CONTACTOUTPUT、*CONTACTF54附加的可用輸出變量包括：輸出變量描述CAREA接觸中從屬表面的全部面積CFN

CFS由從屬表面上的接觸壓力和摩擦剪應(yīng)力各自引起的全部力向量CMN

CMS由從屬表面上的接觸壓力和摩擦剪應(yīng)力各自引起的、關(guān)于原點的全部彎矩向量CFTCFN和CFS的向量和CMTCMN和CFS的向量和XN點的坐標(biāo)，從屬表面上的接觸壓力對該點的全部彎矩為零XS點的坐標(biāo)，從屬表面上的摩擦力對該點的全部彎矩為零XT點的坐標(biāo)，從屬表面上的接觸壓力和摩擦力對該點的全部彎矩為零附加的可用輸出變量包括：輸出變量描述CAREA接觸中從屬表面55與ABAQUS/Standard中接觸分析相關(guān)的，用于生成與打印輸出相關(guān)的兩個其它選項為：*PREPRINT,CONTACT=YES:在前處理過程中，控制打印到輸出文件的輸出給出內(nèi)部生成的接觸單元的細(xì)節(jié)推薦用于小滑動接觸問題，用于驗證主-從節(jié)點的相互作用用于檢測表面定義和相互接觸的正確性*PRINT,CONTACT=YES:在分析過程中，用于控制信息文件的輸出給出迭代過程的細(xì)節(jié)用于分析在接觸計算過程中出現(xiàn)的數(shù)值問題與ABAQUS/Standard中接觸分析相關(guān)的，用于生成與56主/從方法中網(wǎng)格密度的考慮主/從方法中網(wǎng)格密度的考慮57從屬表面的網(wǎng)格應(yīng)該比主控表面的網(wǎng)格更密。如果網(wǎng)格密度相同，從屬表面所占的部分應(yīng)該選擇較軟的材料。不正確主面的網(wǎng)格較密穿透到從屬表面的部分較多。正確主面的網(wǎng)格較稀穿透到從屬表面的部分較少。從屬表面的網(wǎng)格應(yīng)該比主控表面的網(wǎng)格更密。不正確正確58ABAQUS/Standard中的接觸邏輯ABAQUS/Standard中的接觸邏輯59接觸算法檢查平衡6進(jìn)行迭代3確定接觸狀態(tài)1不施加約束2檢查接觸狀態(tài)的變化4點張開；嚴(yán)重不連續(xù)迭代5施加約束2點接觸；嚴(yán)重不連續(xù)迭代5開始增量步張開接觸p<0h>0接觸增量步?jīng)]有變化不收斂7收斂8接觸算法檢查平衡6進(jìn)行迭代3確定接觸狀態(tài)1不施加約束2檢查接60在ABAQUS/Explicit中定義接觸在ABAQUS/Explicit中定義接觸61通用接觸通用接觸用于多組件，并具有復(fù)雜拓?fù)潢P(guān)系的模型建模容易使用“自動接觸”節(jié)省生成模型的時間通用接觸算法一般比雙面接觸算法快在“ABAQUS/Explicit:AdvancedTopics”課程中將詳細(xì)討論通用接觸。通用接觸通用接觸用于多組件，并具有復(fù)雜拓?fù)潢P(guān)系的模型建模在“62常用的、最簡單的接觸定義：*CONTACT*CONTACTINCLUSIONS,ALLELEMENTBASED整個模型的“自動接觸”分析中的接觸定義逐漸變得更加詳細(xì)全局/局部摩擦系數(shù)和其它接觸屬性用戶控制的接觸厚度（尤其是殼）以接觸對的方式制定接觸域（而不是ALLELEMENTBASED）常用的、最簡單的接觸定義：63可以為指定的接觸區(qū)域獨立的分配接觸屬性*CONTACT*CONTACTINCLUSIONS,ALLELEMENTBASED*CONTACTPROPERTYASSIGNMENT,,prop_1 (以全局的方式重新制定屬性)alum_surf,steel_surf,prop_2 (局部修改)alum_surf,alum_surf,prop_3 (局部修改)可以為指定的接觸區(qū)域獨立的分配接觸屬性*CONTACT64將接觸對作為部分的歷程數(shù)據(jù)將接觸對作為部分的歷程數(shù)據(jù)65在輸入文件中，對于ABAQUS/Explicit，將接觸對定義作為歷程數(shù)據(jù)的一部分。*HEADING..*STEP*DYNAMIC,EXPLICIT,200E-3*CONTACTPAIRASURF,BSURF...*ENDSTEP可以根據(jù)需要，以分析步為基礎(chǔ)，定義接觸對、或刪除接觸對。在輸入文件中，對于ABAQUS/Explicit，將接觸對定66

*CONTACTPAIR選項可以包含OP參數(shù)，OP參數(shù)可以的值為ADD或DELETE。例子：*STEP*DYNAMIC,EXPLICIT..*CONTACTPAIRASURF,BSURF*ENDSTEP*STEP..*CONTACTPAIR,OP=DELETEASURF,BSURF*CONTACTPAIR,OP=ADDBSURF,CSURF*ENDSTEP刪除包含表面ASURF和BSURF的接觸對。添加包含表面BSURF和CSURF的接觸對。*CONTACTPAIR選項可以包含OP參數(shù)，OP參數(shù)可67實施接觸約束實施接觸約束68可以用下面的算法之一實施接觸約束：運動依從（只有接觸對算法可用）罰函數(shù)多數(shù)情況下，運動依從算法和罰函數(shù)算法將得到幾乎一致的結(jié)果。然而，在某些情況下，一種方法可能比另一種方法更可取?？梢杂孟旅娴乃惴ㄖ粚嵤┙佑|約束：69運動依從接觸默認(rèn)的運動接觸公式達(dá)到的計算精度與接觸條件相一致。在多數(shù)情況下，它工作得很好。但是在某些情況下，如抖動接觸，使用罰函數(shù)接觸會更容易得到收斂的解。不能為剛體-剛體接觸建模。只對接觸對算法可用。運動依從接觸70罰函數(shù)接觸罰函數(shù)接觸算法中接觸約束的嚴(yán)格性低于運動學(xué)算法。罰函數(shù)算法可以處理更一般類型的接觸；比如，剛體之間的接觸。因為罰函數(shù)算法在模型中引入附加剛度，該附加剛度將影響穩(wěn)定時間增量。通用接觸算法中使用罰函數(shù)接觸約束。如果需要在接觸對算法中使用罰函數(shù)接觸約束：在*CONTACTPAIR選項中，設(shè)置MECHANICALCONSTRAINT參數(shù)為PENALTY。罰函數(shù)接觸71ABAQUS/Explicit自動計算由侵徹距離引入的“彈簧”剛度或“罰”剛度。但是必須考慮相關(guān)的影響：應(yīng)該盡量減少對穩(wěn)定時間增量的影響。在所有的分析中，允許的侵徹不可以太大。對于接觸對算法：通過在*CONTACTCONTROLS選項指定SCALEPENALTY參數(shù)。對于通用接觸算法：可以使用*CONTACTCONTROLSASSIGNMENT,TYPE=SCALEPENALTY選項縮放罰剛度。ABAQUS/Explicit自動計算由侵徹距離引入的“彈簧72第七講接觸非線性

王慎平北京怡格明思軟件技術(shù)有限公司

第七講接觸非線性王慎平北京怡格明思軟73概述簡介主/從方法和隱含假設(shè)定義雙面接觸概述定義面的規(guī)則局部表面行為接觸問題中點的相對滑動接觸問題中調(diào)整初始節(jié)點位置接觸輸出完全約束的接觸行為剛體概述簡介74簡介什么是接觸?當(dāng)兩個實體接觸時，力通過它們的接觸表面?zhèn)鬟f。在某些情況下，只傳遞垂直接觸表面的力。如果存在摩擦，沿接觸表面的切向傳遞切向力。一般目標(biāo)：確定接觸面積和傳遞的應(yīng)力。接觸是嚴(yán)重不連續(xù)形式的非線性行為,是一類特殊的不連續(xù)約束。或者施加約束（表面不可以互相穿透），或者忽略約束。簡介什么是接觸?75為什么要定義接觸？除非用戶指定可能會接觸的表面和/或節(jié)點，ABAQUS現(xiàn)在還不能檢測接觸。為什么要定義接觸？76“DetroitEdisonpipewhipexperiment,”ABAQUS例子手冊2.1.2節(jié)接觸例子-間隙接觸用節(jié)點與節(jié)點接觸的方式為點接觸建模?！癉etroitEdisonpipewhipexpe77 赫茲接觸接觸面之間的相對位移很小。分布表面的接觸。例子選自“Coolantmanifoldcovergasketedjoint,”ABAQUS例子手冊5.1.3節(jié) 赫茲接觸例子選自“Coolantmanifoldco78 變形體之間的有限滑動接觸這是最一般類型的接觸。例子：螺紋連接。這些問題一般包含初始的過盈配合由于過盈精度引起的接觸壓力的分布 變形體之間的有限滑動接觸由于過盈精度引起的接觸壓力的分布79

自接觸單個表面與它自身接觸稱為自接觸。在ABAQUS的二維和三維模型中可用。在分析過中，當(dāng)表面嚴(yán)重變形時，使用自接觸將非常方便。對于某些問題，在分析之前不可能確定單個的接觸區(qū)域，或者確定接觸區(qū)域是非常困難的。把單個接觸表面作為接觸對定義自接觸，而不是通過兩個不同的表面定義。最小主應(yīng)力云圖SURF1(剛體)SURF2例子：橡膠墊片的壓縮(選自“Self-contactinrubber/foamcomponents:rubbergasket,”ABAQUS例子手冊1.1.17節(jié)) 自接觸最小主應(yīng)力云圖SURF1SURF2例子：橡膠墊片的壓80例子：金屬成型模擬變形體與剛體的接觸:表面之間的有限滑動（大位移）。變形組件之間的有限應(yīng)變。典型例子：成型模擬(剛體砧/模具、可變形組件。)例子：金屬成型模擬變形體與剛體的接觸:典型例子：81使用剛體表面減少計算時間

對于兩個相互接觸的物體，如果其中一個物體比另外的物體剛硬許多，可以將較為剛硬的物體指定為剛體；如，金屬成型過程中的砧。ABAQUS不需要求解剛體的變形。在參考點，最多利用六個自由度就可以模擬剛體的運動。有三種方法，可以以幾何的方式定義剛體表面：定義解析剛體表面。使用單元類型組合定義剛體（包括剛體單元），并將物體聲明為剛體（一般剛體功能）。寫用戶子程序(RSURFU)定義表面(只有ABAQUS/Standard可用)。使用剛體表面減少計算時間82主/從方法和隱含假設(shè)主/從方法和隱含假設(shè)83因為接觸表面下面的實體被離散化，接觸表面也必須被離散化。接觸表面分為:主面從面因為接觸表面下面的實體被離散化，接觸表面也必須被離散化。84ABAQUS/Standard中使用的嚴(yán)格的主/從接觸算法。從屬表面的節(jié)點不能穿透到主控表面。主控表面的節(jié)點可以穿透到從屬表面。ABAQUS/Standard中使用的嚴(yán)格的主/從接觸算法。85在多數(shù)情況下，ABAQUS/Explicit使用平衡的主/從算法。兩次應(yīng)用單純的主-從接觸算法，并平均化。在第二次施加接觸約束時，交換主/從表面。使接觸物體之間的侵徹達(dá)到最小化。在多數(shù)情況下，ABAQUS/Explicit使用平衡的主/從86接觸方向總是與主控表面垂直。沿主控表面的法向檢查接觸條件。沿主控表面的法向傳遞接觸力。沿接觸表面的切向傳遞摩擦力。接觸方向總是與主控表面垂直。87定義雙面接觸概述定義雙面接觸概述88定義接觸的三個步驟：

定義接觸表面。定義接觸對。定義接觸屬性。123定義接觸的三個步驟：12389輸入文件中，完整定義接觸句法的例子：*SURFACE,NAME=ASURFSLIDER,S1*SURFACE,NAME=BSURFBLOCK,S3*CONTACTPAIR,INTERACTION=FRIC1ASURF,BSURF*SURFACEINTERACTION,NAME=FRIC11.0,*FRICTION0.4,輸入文件中，完整定義接觸句法的例子：90定義表面利用*SURFACE選項定義表面。利用每個單元集的表面標(biāo)識符指定面?？梢杂脝卧拿只騿卧栔付ū砻妗?SURFACE,NAME=ASURFSLIDER,S1*SURFACE,NAME=BSURFBLOCK,S3接觸發(fā)生在單元集SLIDER的底部(S1)接觸發(fā)生在單元集的BLOCK頂部(S3)定義表面接觸發(fā)生在單元集SLIDER的底部(S1)接觸發(fā)生91定義表面接觸屬性

*SURFACEINTERACTION選項定義表面接觸屬性。定義表面行為屬性，比如摩擦。對于二維問題，定義接觸面在垂直平面方向的厚度。定義表面接觸屬性92*CONTACTPAIR,INTERACTION=FRIC1ASURF,BSURF*SURFACEINTERACTION,NAME=FRIC11.0,*FRICTION0.4,*SURFACEINTERACTION的子選項，列出表面的基本屬性接觸面在垂直平面方向的厚度*SURFACEINTERACTION的子選項，列出表面的93定義接觸對定義表面和接觸屬性之后，就可以定義“接觸對”。在分析過程中，接觸對指定了兩個相互接觸的表面。在ABAQUS/Standard中，第一個表面為從屬表面，第二個表面為主控表面。在ABAQUS/Explicit中，表面的順序一般不影響接觸計算。*CONTACTPAIR,INTERACTION=FRIC1ASURF,BSURF定義接觸對*CONTACTPAIR,INTERACTI94定義面的規(guī)則定義面的規(guī)則95定義面的規(guī)則在變形體或劃分網(wǎng)格的剛體上，使用*SURFACE,TYPE=ELEMENT選項定義表面。通過指定單元表面標(biāo)識符定義表面。讓ABAQUS自動確定連續(xù)體單元的“自由表面”。使用*SURFACE,TYPE=[SEGMENTS|CYLINDER|REVOLUTION]選項和*RIGIDBODY選項定義解析剛體表面。使用*SURFACE,TYPE=NODE選項，定義可能接觸的節(jié)點表面。定義面的規(guī)則在變形體或劃分網(wǎng)格的剛體上，使用*SURFAC96在實體單元上定義表面使用表面標(biāo)識符例子：4-節(jié)點四面體單元(CPE4,CAX4,等等)*SURFACE,NAME=EXAMPLE11,S41,S12,S12,S2...在實體單元上定義表面97使用自動的表面定義*SURFACE,NAME=EXAMPLE2ELSET1,沒有面標(biāo)識符使用自動的表面定義沒有面標(biāo)識符98在結(jié)構(gòu)單元上定義表面(殼、薄膜、剛體、梁)結(jié)構(gòu)單元的法向表示預(yù)期的接觸方向。法向基于單元局部的節(jié)點編號。正的法線方向在=SPOS表面。負(fù)的法線方向在=SNEG表面。殼和薄膜(S4R,S8R,M3D4,等)2-D桁架和梁(B21,T2D2,等)在結(jié)構(gòu)單元上定義表面(殼、薄膜、剛體、梁)殼和薄膜(S99在所定義的表面集之內(nèi)，表面法向的定義必需一致。*ELEMENT,TYPE=B21,ELSET=BOTTOM10,1,211,2,312,3,6*ELEMENT,TYPE=B21,ELSET=TOP20,4,521,5,6*ELSET,ELSET=BEAMSBOTTOM,TOP*SURFACE,NAME=SURF1BOTTOM,SPOSTOP,SNEG*SURFACE,NAME=SURF1BEAMS,SPOS在所定義的表面集之內(nèi)，表面法向的定義必需一致。*SURFAC100解析剛體表面所有的剛體表面將被作為主控表面；剛體表面的自由度與剛體的參考節(jié)點相關(guān)聯(lián)，而不是表面上的點。使用*SURFACE選項，有三種可用的解析表面：使用TYPE=SEGMENTS定義二維剛體表面。使用TYPE=CYLINDER定義三維剛體表面，在垂直于平面方向無限延伸。使用TYPE=REVOLUTION定義旋轉(zhuǎn)的三維表面。解析剛體表面101基于節(jié)點的表面與指定單元表面定義接觸表面的方法不同，基于節(jié)點的表面只包含節(jié)點?；诠?jié)點的表面一般被用作從屬表面。基于節(jié)點的表面102基于節(jié)點表面的例子：定義包含接觸節(jié)點的表面前面定義的表面線：基于節(jié)點的表面球：基于單元的表面*SURFACE,TYPE=NODE,NAME=STRINGSSTRINGS,*CONTACTPAIR,INTERACTION=SMOOTHSTRINGS,BALL基于節(jié)點表面的例子：定義包含接觸節(jié)點的表面前面定義的表面線：103表面定義的限制接觸表面下所有的單元必須協(xié)調(diào)。它們必須：維數(shù)相同（二維或三維）。對于二維表面：所有單元為平面的或軸對稱的（不可以兩種都有）。插值階數(shù)相同（一階或二階）。所有單元為可變形的或所有單元為剛體（不可以兩種都有）。對于ABAQUS/Standard中的主控表面和ABAQUS/Explicit中的所有表面，還有附加的限制：在不離開表面、穿過表面或通過單個點穿越表面的前提下，必須可以跨過表面中的兩個任意的點。表面定義的限制104對于三維情況，限制的例子。需要穿過或離開表面才能跨過。合法的、封閉的3-D主表面合法的、簡單連接的3-D主表面不合法的3-D主表面不合法的3-D主表面不可以通過一個點跨過對于三維情況，限制的例子。需要穿過或離開表面才能跨過。合法的105局部表面行為局部表面行為106一般的接觸建模包括不同方向的接觸相互作用：主控表面的法向表面的切向一般的接觸建模包括不同方向的接觸相互作用：107接觸法向行為硬接觸在所有接觸問題中，“硬”接觸是默認(rèn)的局部行為。硬接觸可以通過經(jīng)典的拉格朗日乘子方法實施，或者通過增廣拉格朗日方法實施：

壓力－間距關(guān)系*surfaceinteraction,name=...*surfacebehavior,augmentedlagrange接觸法向行為壓力－間距關(guān)系*surfaceinteract108除硬接觸外，其它的方法作為*SURFACEINTERACTION

選項的子選項*SURFACE

BEHAVIOR選項用于指定：軟接觸（指數(shù)或表格方式

表達(dá)的壓力-間距關(guān)系）沒有分離的接觸其它選項：間距相關(guān)的粘性阻尼(*CONTACTDAMPING)。帶有過盈或拉伸接觸力的接觸(*CONTACTCONTROLS;只有ABAQUS/Standard可用)。除硬接觸外，其它的方法109接觸的切向行為兩個實體之間的接觸面可能存在摩擦剪切應(yīng)力。如果剪切應(yīng)力達(dá)到某個臨界值，實體之間可能會發(fā)生相對滑動；否則它們粘在一起。接觸的切向行為110摩擦可以引起高度的非線性效應(yīng)。很難獲得解。除非物理上必要，盡量不要使用。摩擦是非保守的。在ABAQUS/Standard中，摩擦將導(dǎo)致系統(tǒng)方程的不對稱。對于較大的摩擦系數(shù)(

>0.2)，*STEP,UNSYMM=YES選項將自動啟用。使用UNSYMM=NO將減緩收斂速度，但是解將是正確的（如果得到）。它可能會使用更少的磁盤空間。在ABAQUS/Explicit中，過大的摩擦系數(shù)將不會引起類似的問題，因為不需要求解系統(tǒng)方程組。摩擦可以引起高度的非線性效應(yīng)。111默認(rèn)情況下，ABAQUS使用庫倫摩擦模型。臨界摩擦應(yīng)力取決于接觸壓力：cr=p?；揪浞ǎ?FRICTION摩擦系數(shù)可以是相對滑動速度、壓力、溫度和場變量的函數(shù)。出于計算的考慮，ABAQUS/Standard中默認(rèn)的摩擦模型對理想行為作一定的近似：在達(dá)到不可恢復(fù)滑動發(fā)生之前，允許小量的彈性滑動。默認(rèn)情況下，ABAQUS使用庫倫摩擦模型。112可以定義混合的靜動摩擦模型。作為滑動速度函數(shù)的摩擦系數(shù)m，從ms（靜摩擦系數(shù)）到mk（動摩擦系數(shù)）以指數(shù)方式衰退。還有其它可用的摩擦行為，包括用戶自定義的摩擦模型。：

(在ABAQUS/Standard中為FRIC在ABAQUS/Explicit中為VFRIC)?？梢远x混合的靜動摩擦模型。113接觸問題中點的相對滑動接觸問題中點的相對滑動114兩種滑動距離選項：有限滑動：默認(rèn)為有限滑動，有限滑動是最通用的滑動方式。 在接觸表面之間，允許任意大的滑動和旋轉(zhuǎn)。小滑動： 在接觸表面之間，允許小的相對滑動。 只要接觸表面之間沒有大的相對移動，允許接觸表面

之間有大的轉(zhuǎn)動。

相比有限滑動，具有較小的計算費用。兩種滑動距離選項：115有限滑動在有限滑動問題中的主控表面，最好使用延伸到角點之外的定義。它們將防止從屬節(jié)點“跌落”到主表面的后面或被主表面的后面捕獲。對于ABAQUS/Standard中的接觸模擬，這是尤為重要的。不合適的定義(修整過的)合適的(未修整的)主控表面定義有限滑動不合適的定義(修整過的)合適的(未修整的)主控表116小滑動*CONTACTPAIR,INTERACTION=FRIC1,SMALLSLIDING

小滑動問題不需要通用的有限滑動算法。在整個分析過程中，從屬節(jié)點與主控表面上預(yù)先確定的、固定數(shù)量的節(jié)點相互作用。只是在分析開始時進(jìn)行接觸搜尋，所以不需要復(fù)雜的有限滑動搜尋算法。優(yōu)點：相對有限滑動，有更低的計算費用，尤其是三維問題。小滑動117小滑動接觸約束在分析開始時，ABAQUS將確定哪些節(jié)點是相互接觸的。右圖為小滑動選項基礎(chǔ)之上，節(jié)點101和它的主控表面BSURF之間可能的接觸演化方式。小滑動接觸約束118主控表面上，以內(nèi)插值方式得到節(jié)點201和202之間部分的法向量指向從屬節(jié)點101。在分析過程中，節(jié)點101不可以穿過垂直于法向量的平面。因為接觸線沒有端點，所以從屬節(jié)點不可能“滑出”接觸線。在分析過程中，如果使用NLGEOM選項，平面的方向?qū)⒒诠?jié)點201和202的變形而更新。主控表面上，以內(nèi)插值方式得到節(jié)點201和202之間部分的法向119接觸問題中調(diào)整初始節(jié)點位置接觸問題中調(diào)整初始節(jié)點位置120在分析之前，可以調(diào)整接觸面中節(jié)點的初始位置，并且不產(chǎn)生任何的應(yīng)力和應(yīng)變。ABAQUS/Standard允許用戶通過指定絕對距離或節(jié)點集的方式調(diào)整節(jié)點位置。*CONTACTPAIR,ADJUST=[distance|nodesetname]在分析之前，可以調(diào)整接觸面中節(jié)點的初始位置，并且不產(chǎn)生任何的121指定需要調(diào)整的絕對距離：*CONTACTPAIR,INTERACTION=FRIC1,ADJUST=a在調(diào)整范圍內(nèi)，初始張開的從屬節(jié)點將被移到主控表面上。調(diào)整距離沿主控表面的法向度量。所有初始過盈的從屬節(jié)點被重新放在主控表面上。指定需要調(diào)整的絕對距離：122指定需要調(diào)整的從屬節(jié)點集：在分析過程中，不在節(jié)點集

中并帶有初始過盈的從屬節(jié)點

將保持過盈，并產(chǎn)生應(yīng)變。

這些應(yīng)變將被施加在接觸面上。

指定

ADJUSTnodeset的例子：*NSET,NSET=CONNODE,GENERATE1,8,1*CONTACTPAIR,INTERACTION=RIG,ADJUST=CONNODE指定需要調(diào)整的從屬節(jié)點集：123ABAQUS/Explicit不允許接觸表面的初始過盈。在分析之前，接觸面上的節(jié)點將被自動調(diào)整，以刪除任何的初始過盈。在隨后的分析步中，這樣的調(diào)整將引起應(yīng)變。大量的調(diào)整將引起初始單元形狀的嚴(yán)重扭曲。如果發(fā)現(xiàn)這樣的錯誤信息，最好需要運行datacheck分析，并利用ABAQUS/Viewer查找問題。ABAQUS/Explicit不允許接觸表面的初始過盈。124接觸輸出接觸輸出125在使用*CONTACTOUTPUT、*CONTACTFILE或*CONTACTPRINT選項之后，所有接觸表面的接觸信息都是可用的。對于ABAQUS/Standard和ABAQUS/Explicit：接觸應(yīng)力：CSTRESS(接觸壓力CPRESS和摩擦剪切應(yīng)力CSHEAR1和CSHEAR2)對于ABAQUS/Standard，還可以得到：接觸位移：CDISP(接觸張開距離COPEN和相對切向運動CSLIP1和CSLIP2)對于ABAQUS/Explicit，還可以得到：滑動速度：FSLIPR累計的滑動距離：FSLIP在使用*CONTACTOUTPUT、*CONTACTF126附加的可用輸出變量包括：輸出變量描述CAREA接觸中從屬表面的全部面積CFN

CFS由從屬表面上的接觸壓力和摩擦剪應(yīng)力各自引起的全部力向量CMN

CMS由從屬表面上的接觸壓力和摩擦剪應(yīng)力各自引起的、關(guān)于原點的全部彎矩向量CFTCFN和CFS的向量和CMTCMN和CFS的向量和XN點的坐標(biāo)，從屬表面上的接觸壓力對該點的全部彎矩為零XS點的坐標(biāo)，從屬表面上的摩擦力對該點的全部彎矩為零XT點的坐標(biāo)，從屬表面上的接觸壓力和摩擦力對該點的全部彎矩為零附加的可用輸出變量包括：輸出變量描述CAREA接觸中從屬表面127與ABAQUS/Standard中接觸分析相關(guān)的，用于生成與打印輸出相關(guān)的兩個其它選項為：*PREPRINT,CONTACT=YES:在前處理過程中，控制打印到輸出文件的輸出給出內(nèi)部生成的接觸單元的細(xì)節(jié)推薦用于小滑動接觸問題，用于驗證主-從節(jié)點的相互作用用于檢測表面定義和相互接觸的正確性*PRINT,CONTACT=YES:在分析過程中，用于控制信息文件的輸出給出迭代過程的細(xì)節(jié)用于分析在接觸計算過程中出現(xiàn)的數(shù)值問題與ABAQUS/Standard中接觸分析相關(guān)的，用于生成與128主/從方法中網(wǎng)格密度的考慮主/從方法中網(wǎng)格密度的考慮129從屬表面的網(wǎng)格應(yīng)該比主控表面的網(wǎng)格更密。如果網(wǎng)格密度相同，從屬表面所占的部分應(yīng)該選擇較軟的材料。不正確主面的網(wǎng)格較密穿透到從屬表面的部分較多。正確主面的網(wǎng)格較稀穿透到從屬表面的部分較少。從屬表面的網(wǎng)格應(yīng)該比主控表面的網(wǎng)格更密。不正