第25例非線性分析綜合應(yīng)用實(shí)例 鋼板卷制成圓筒本例介紹了綜合利用ANSYS非線性分析功能模擬將鋼板卷制成圓筒的方法和步驟。25.1 問題描述將鋼板卷制成圓筒一般要使用卷板機(jī)。圖25-1所示為對(duì)稱式三輥卷板機(jī)，該機(jī)器將鋼板卷制成圓筒時(shí)分為三個(gè)步驟：首先，上輥下降使鋼板發(fā)生撓曲，鋼板撓曲線服形成圓筒；最終，圓筒卷制完成，上輥上升卸下筒體。25-1對(duì)稱式三輥卷板機(jī)用ANSYS模擬將鋼板卷制成圓筒，相應(yīng)地也分為三個(gè)步驟。由于其次個(gè)步驟需要模擬上、下輥轉(zhuǎn)動(dòng)，而ANSYSSOLIDn單元不支持大轉(zhuǎn)動(dòng)，位移邊界條件不能施加大的轉(zhuǎn)動(dòng)角度，所以上、下輥需要用殼單元建立有限元模型。上、下輥與鋼板的作用需要用接觸模材料非線性和構(gòu)造非線性三種非線性。用ANSYS模擬將鋼板卷制成圓筒，計(jì)算結(jié)果可以得到圓筒直徑與上輥下壓量的關(guān)系，上、下輥受力大小，上、下輥的變形，下輥驅(qū)動(dòng)力矩及卸載回彈等重25.2/CLEAR/FILNAM,EXAMPLE25/CONFIG,NRES,2023/PREP7/PNUM,VOLU,ONET,1,SHELL181ET,2,SOLID186MP,EX,1,2E11MP,DENS,1,7800MP,NUXY,1,0.3MP,EX,2,2E11MP,DENS,2,7800MP,NUXY,2,0.3TB,BKIN,2,1TBTEMP,0TBDATA,,240E6,0SECTYPE,1,SHELLSECDATA,0.02CYLIND,0.38/2,0,0.2,1.7,0,360

要數(shù)據(jù)。由于分析過程簡潔，步驟較多，所以本例只承受命令流法執(zhí)行命令。!去除數(shù)據(jù)庫，建文件!指定任務(wù)名為?EXAMPLE25“!設(shè)置最大子步數(shù)!前處理!進(jìn)入前處理器!翻開體號(hào)!選擇單元類型，殼單元用于劃分上、下輥!實(shí)體單元用于劃分鋼板!定義材料模型1的彈性模量!定義材料模型1的密度!定義材料模型1的泊松比!定義材料模型2的彈性模量!定義材料模型2的密度!定義材料模型2的泊松比!定義材料模型2的屈服極限、切向模量!定義截面!殼厚度創(chuàng)立上輥，圓柱體CYLIND,0.38/4,0,0.2,0,0,360CYLIND,0.38/4,0,1.7,1.9,0,360VGLUE,ALLVSEL,NONEWPOFF,0.54/2,-0.37CYLIND,3600.3/2,0,0.2,1.7,0,CYLIND,0.3/4,0,0.2,0,0,360CYLIND,0.3/4,0,1.7,1.9,0,360WPOFF,-0.54VGLUE,ALLVSEL,NONECYLIND,3600.3/2,0,0.2,1.7,0,CYLIND,0.3/4,0,0.2,0,0,360CYLIND,0.3/4,0,1.7,1.9,0,360ALLS VGLUE,ALLVDELE,ALL !!HPTCREATE,AREA,5,,COORD

!黏結(jié)!偏移工作平面原點(diǎn)到右下輥中心!!偏移工作平面原點(diǎn)到左下輥中心!!在上輥端面中心處創(chuàng)立硬點(diǎn)八、、 HPTCREATE,AREA,10,,COORD,0,0,1.91.9

HPTCREATE,AREA,34,,COORD,-0.54/2,-0.37,

!在左下輥端面中心處創(chuàng)建硬點(diǎn)HPTCREATE,AREA,29,,COORD,-0.54/2,-0.37,0HPTCREATE,AREA,22,,COORD,0.54/2,-0.37,1.9 !在右下輥端面中心處創(chuàng)硬點(diǎn)HPTCREATE,AREA,17,,COORD,0.54/2,-0.37,0MAT,1TYPE,1SECN,1ESIZE,0.04SMRTSIZE,7MSHAPE,1AMESH,ALLBLOCK,0.7,-0.5,0.3/2,0.3/2+0.03,0.5,

!為上、下輥劃分單元指定屬性，材料模型!單元類型!殼橫截面!指定單元邊邊長度!指定智能尺寸級(jí)別!指定單元外形為三角形!對(duì)全部面劃分單元LESIZE,96,,,2MAT,2

1.4

!創(chuàng)立鋼板，塊!指定線96〔鋼板厚度〕被劃分為兩段TYPE,2MSHKEY,1MSHAPE,0ESIZE,0.05VMESH,ALLET,4,TARGE170ET,5,CONTA174,,,,,1KEYOPT,5,12,1

!為鋼板劃分單元指定屬性!映射網(wǎng)格!指定單元外形為六面體!指定單元邊長度!對(duì)塊劃分網(wǎng)格!指定單元類型，用于創(chuàng)立接觸對(duì)!設(shè)置單元5接觸外表無滑動(dòng)R,10,,,0.1!定義實(shí)常數(shù)，KFN==0.1R,11,,,0.1R,12,,,0.1ASEL,S,,,40!在鋼板上外表和上輥外表建立接觸對(duì)NSLA,S,1NSEL,U,LOC,MAT,3Z,0.8,1.1REAL,10TYPE,4ESURFALLSASEL,S,,,3,4,1NSLA,S,1NSEL,U,LOC,Z,0.2,0.5NSEL,U,LOC,Z,0.8,1.1NSEL,U,LOC,Z,1.4,1.7MAT,3REAL,10TYPE,5ESURFALLSASEL,S,,,39！在鋼板下外表和左下輥外表建立接觸對(duì)NSLA,S,1NSEL,U,LOC,Z,0.8,MAT,41.1REAL,11TYPE,4ESURFALLSASEL,S,,,27,28,1NSLA,S,1NSEL,U,LOC,Z,0.2,0.5NSEL,U,LOC,Z,0.8,1.1NSEL,U,LOC,Z,1.4,1.7MAT,4REAL,11TYPE,REAL,11TYPE,5ESURFALLSASEL,S,,,39NSLA,對(duì)S,1NSEL,U,LOC,Z,0.8,1.1MAT,4REAL,12TYPE,4ESURFALLSASEL,S,,,15,16,1NSLA,S,1NSEL,U,LOC,Z,0.2,0.5NSEL,U,LOC,Z,0.8,1.1NSEL,U,LOC,Z,1.4,1.7MAT,4REAL,12TYPE,5ESURFALLSFINI!首先模擬下輥不動(dòng)，上輥下降/SOLU

!進(jìn)入求解器D,ALL,UZD,ALL,UZD,ALL,ROTXD,ALL,ROTYNSEL,S,LOC,XNSEL,R,LOC,YD,ALL,UY,-0.03

!選擇上輥軸線上的節(jié)點(diǎn)!施加約束，上輥下降0.03mD,ALL,UZALLSCSYS,4ASEL,S,LOC,Y,-2,0.3/2ASEL,U,,,39NSLA,S,1D,ALL,ALLALLSLNSRCH,ONNLGEOM,ONTIME,1AUTOT,ONNSUBST,15,20,10KBC,0OUTRES,ALL,ALLSOLVE!以下模擬下輥轉(zhuǎn)動(dòng)，驅(qū)動(dòng)鋼板前進(jìn)CSYS,4ASEL,S,LOC,Y,-2,0.3/2ASEL,U,,,39

!選擇下輥上的全部面!選擇下輥上的全部節(jié)點(diǎn)!約束掉全部自由度!瞬態(tài)分析!翻開線性搜尋!翻開大變形選項(xiàng)!指定載荷步時(shí)問!翻開自動(dòng)載荷步長!指定子步數(shù)目!斜坡載荷!輸出全部子步全部項(xiàng)舊的結(jié)果!解答!選擇下輥上的全部面NSLA,NSLA,S,1!選擇下輥上的全部節(jié)點(diǎn)!刪除上一載荷步施加在下輥上的全部DDELE,ALL,ALL約束ALLSD,D,ALL,UZ!在全部節(jié)點(diǎn)上施加約束D,ALL,ROTXD,ALL,ROTYNSEL,S,LOC,X,0NSEL,R,LOC,Y

!選擇左下輥軸線上的全部節(jié)點(diǎn)D,D,ALL,ROTZ,-3.1415926/2!轉(zhuǎn)動(dòng)1/4周D,ALL,UX!約束掉移動(dòng)D,ALL,UYALLSWPOFF,0.54NSEL,S,LOC,X,0NSEL,R,LOC,Y

!偏移工作平面!選擇右下右下輥軸線上的全部節(jié)點(diǎn)D,D,ALL,UX!約束掉移動(dòng)D,ALL,UYD,ALL,ROTZ,-3.1415926/2TIME,2ALLSKBC,0NSUBST,800,1000,300

!轉(zhuǎn)動(dòng)1/4周SOLVE!解答!模擬卸載!施加重力加速度ACEL,0,9.8!激活全球直角坐標(biāo)系CSYS,0!選擇上輥軸線上的全部節(jié)點(diǎn)NSEL,S,LOC,XNSEL,R,LOC,Y!上輥向上移