【文章內(nèi)容簡介】 常數(shù) ? mat,1 ! 接觸面 1定義材料與超彈 ? 單元同 ? lesize,15,1 ! 目標(biāo)面 1（ L15）分網(wǎng)尺寸 ? (圖 5－ 3) ? lmesh,15 ! 目標(biāo)面 1（ L15）分網(wǎng) ,采用自動分網(wǎng)，此剛性面自動約束。檢查外法線方向 【 OK】 ? ! *** 接觸面 1：實常數(shù)（與目標(biāo)面 1一致）、 分網(wǎng) lsel,s,line,1,3,1 ! 選 L1,L2,L3 ? lsel,a,line,9,10,1 ! 再加 L9,L10 ? lsel,a,line,14 ! 再加 L14 ? nsll,s,1 ! 選擇和所選線相連系的節(jié)點 ? type,3 ! Conta171這里未發(fā)布 Real,1，實常數(shù)未重新定義就是前面的 Real,1， MAT也是前面的 MAT, 1 ? esurf ! 生成接觸單元 ? lsel,all ? ! *** 檢查外法線方向 【 OK】 *** 圖 5－ 3 ? ! *** 創(chuàng)建接觸對 2 *** ? type,2 ! 目標(biāo)面 2 Target169 ? real,2 ! MAT未重新定義，就是前面的MAT, 1 ? lesize,16,1 ! 目標(biāo)面 2 ? lmesh,16 ! L16分網(wǎng)（目標(biāo) 2分網(wǎng)） ? kmesh,27 ! 指定 Pilot點 ? ! *** 檢查外法線方向 【 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】 *** ? Step ? ! 剛性面 1 被約束 ? ! 剛性面 2 隨 Pilot點移動 ? ! L4對稱約束 ? dl,4,ux,0 ? n_load=node(kx(27),ky(27),0) ! 為定義剛性面 2的位移做準備 ? finish ? Step ? /solu ? nlgeom,on ? solc,on ! （缺?。? ? time, ! 載荷步、結(jié)束時間 ? d,n_load,uy, ? nsubst,25,500,10 ? outres,all,all ? monitor,var3,n_load,fy ? Step ? solve ? Step ? /post1 ? pldsp,2 ! 變形圖 ? plnsol,s,eqv,0,1 ! Von Mises 應(yīng)力云圖 ? plnsol,cont,pres,0, ! 接觸壓力圖 ? plnsol,epto,eqv ! 繪等效總應(yīng)變圖 5－ 4 ? save,seal,db 定義反力變量 ,繪載荷變形圖 ? Utility Menu Plot Elements ? Main Menu Time Hist Postpro Define Variables ? 【 Add】 ? Reaction force ? 【 OK】 ? 拾取節(jié)點 263(控制點 ) ? 【 OK】 ? UserSpecified Load = FORCE: FY ? 【 OK】 ? 【 close】 ? 或命令： ? /POST 26 ? RFORCE,2,263,F,Y,FORCE ? Main Menu Time Hist Postpro Graph Variables ? 1st Variable to gragh = 2 ? 【 OK】 ? 或命令： PLVAR, 2 (圖 5－ 4) 圖 5－ 4 ? 用接觸向?qū)?chuàng)建接觸對，用 GUI方式創(chuàng)建Pilot節(jié)點（剛性面控制點），然后求解。 ? Step （包括基體的幾何、單元、分網(wǎng)；沒有選接接觸單元與目標(biāo)單元，未定義接觸對） ? Utility Menu File Resume from… ? 選 【 OK】 ? Step ? Main Menu Preprocessor Modeling Create Contact Pair ? 圖 4－ 2 【 contact Wizard】 1 創(chuàng)建接觸對 1 ? 圖 5－ 5 選線 (1)( 接觸對 1的目標(biāo)面 )【 OK】 返回圖 5－ 5【 next】 ? 圖 5－ 6 ? 選線 (2)( 接觸對 1的接觸面 )【 OK】 ；返回圖 5－ 6【 next】 Step ? 圖 5－ 7 ? 圖 5－ 8 圖 5－ 9 ? 定義摩擦 (圖 5－ 9) 【 OK】 返回圖 5－ 7 Coefficient friction 【 Create】 ? 5－ 11 圖 5－ 12 ? 接觸對 1(圖 5－ 11) 圖 5－ 10圖 Step (1).設(shè)置單元屬性 用命令流： Type,4 MAT,1 REAL,4 TSHAP,PIL10 GUI: Main Menu Preprocessor Modeling Create Element Elem Attributes ? 圖 5－ 14 ? (2).繪關(guān)鍵點 ? Utility Menu Plot Keypoints Keypoints ? (3).建立目標(biāo)單元控制點 ? Main Menu Preprocessor Meshing Mesh Keypoints ? 選 K27 【 OK】 ? 或命令流： KMESH,27 ? Step ? 控制點 27上施加 UY = ? Main Menu Solution Define Loads Apply Structural Displacement On Keypoints ? 選 K27 【 OK】 ? Lab2 DOFs to be constrained: UY ? Value: ? 【 OK】 ? 施加對稱邊界條件： ? Main Menu Solution Define Loads Apply Structural Displacement ? Symmetry . On Lines ? 選線 L4 【 OK】 ? Step 、求解 ? Main Menu Solution Analysis Type Sol′n Controls ? Main Menu Solution Solve LS_Current ? Step ? 同 1 ? 無摩擦 (MU = 0) Von Mises = ? 有摩擦 (MU = ) Von Mises = ? 接觸剛度取 不收斂 ? 收斂 ? /POST26 ? Main Menu TimeHist Postpro ? 選控制點 27， 繪圖 (見圖 54)。 面一面接觸具有 20個可用實常數(shù)， 2個材料屬性和 30個可用單元選項設(shè)置。能夠模擬特殊的效果和處理困難的收斂情況。 通常分析開始先用缺省值，只指定罰剛度和穿透容差及子步數(shù)。 只能在缺省設(shè)置遇到困難時才用高級選項。 所有的單元選項和參數(shù)都可以通過接觸向?qū)砜刂?；也可以通過實常數(shù)和單元選項來 指定 。 Conta178接觸單元、實常數(shù)選項見圖 61及 62。 ? 圖 61 圖 62 167。 6 面一面接觸單元高級接觸選項 進入接觸選項菜單用下列過程： ? Main Menu Preprocessor Modeling Create Contact Pair ? 選 【 property】 ? 第 1個重要選項：如何保證接觸協(xié)調(diào)性： ? 使用修正的拉朗日法（缺?。? ? 純罰函數(shù)方法 ? 在接觸向?qū)е?Basic Contact algorithm: Augmented Lagrange MethodPenalty Method ? －對大多數(shù)模型，修正的拉格朗日方法能很好地工作。 ? －罰函數(shù)法推薦用于具有變形很大的單元，很大的摩擦系數(shù)和 /或用修正的拉格朗日方法時收斂性很差的情況。 最重要的選項是法向罰剛度或接觸剛度。 ? －對于大面積接觸采用起始值（因子） ，對于柔性接觸采用 。 ? —大值對應(yīng)較高精度；小值對應(yīng)較好的收斂性。 ? 有時最好先以較小的接觸剛度進行分析，然后在一系列載荷步中逐漸增大剛度－ “ 漸變 ” 的接觸剛度，提高收斂性 ? 圖 62 ? 在最后的載荷步逐漸提高到一個剛硬的值將提高計算精度?？赏ㄟ^向?qū)У幕颈砘騿卧匦圆藛卧O(shè)置以允許程序更新接觸剛度 ? 接觸向?qū)?Basic Normal Penalty Stiffness ⊙ factor ○ constant ? 向?qū)?Basic contact stiffness update: Each load step (PAIR ID based)noneEach load step(允許用戶指定剛度變化 )Each substep(允許自動和用戶指定變化 ) 第 2個重要選項－穿透容差。 ? 缺省情況下，穿透容差是一個因子乘以基體單元厚度 。 ? －對于變化很大的網(wǎng)格密度，采用因子會在接觸表面的某些部分產(chǎn)生太小的容差，這時采用絕對值可能更好。 ? －不要使用太小的容差，因為它總是對收斂性有害。 對于臨界接觸狀態(tài)變化的自動時間步控制 －不控制：不影響自動時間步，對靜力問題自動時間步打開時此選項一般是足夠的。 －自動二分：如果接觸狀態(tài)變化明顯，時間步長將二分。對于動力問題自動二分通常是足夠的。 －合理值：比自動細分更耗時的算法。 －最小值：此選項為下一步預(yù)測最小時間增量（很耗計算時間，不推薦） ? 或在單元選項中控制： ? Contact time/load prediction K7: No prediction Pinball區(qū)域影響接觸狀態(tài)的確定和其它許多接觸特性。 ? Pinball區(qū)域是環(huán)繞接觸單元的園（ 2D）或球（ 3D），描述接觸單元周圍 “ 遠 ” 和 “ 近 ” 區(qū)域的邊界 (圖 63)。 ? 在缺省情況下， Pinball區(qū)域半徑是 4 基體單元厚度（剛－柔）或 2 基體單元厚度（柔一柔） ? 可以為 Pinball半徑指定一個不同的值。 圖 63 ? 也可用實常數(shù) PINB調(diào)整球形區(qū)（對于初始值侵入大的問題是必要的）。 幾種不同的接觸模式 這些選項使你能夠模擬特殊的物理現(xiàn)象。 ? 或單元選項 : Behavior contact surface K12: standard ? 這些選項包括： ? －標(biāo)準：正常的接觸閉合和打開行為，具有正常的粘著 /滑動摩擦行為。 ? －粗糙：正常接觸閉合和打開行為，但不發(fā)生滑動（類似于具有無限摩擦系數(shù)） ? －不分離（滑動）：目標(biāo)面和接觸面一旦建立接觸就不再分離（允許滑動） ? －綁定：目標(biāo)面和接觸面一旦接觸就粘在一起 ? －不分離（永遠）：初始位于 Pinball區(qū)域內(nèi)或已經(jīng)接觸的接觸檢查點在法向不分離（允許滑動） ? －綁定接觸（永遠）：初始位于 Pinball區(qū)域內(nèi)或已經(jīng)接觸的接觸檢查點在剩余的分析過程中綁定在一起（ Design Space缺省值） ? －綁定接觸（初始接觸）：只在初始