【求助】《ANSYS14.0超级学习手册》三辊卷板机仿真错误
本帖最后由 liujj15 于 2016-3-7 15:25 编辑手头有这本书，这个算例是17章第一个，我是初学者，跟着提示做下来感觉问题很多。
命令流如下，我自己加了一些注释，大家忙的话可以帮我跑一下看看是不是我自己计算机的问题（做到第一个solve就行了），稍微有点时间的话可以详细看看里面哪里有问题，多谢多谢！
我自己做到 0.332858795s 处就判定不收敛了，再ANSYS14.0和ANSYS15.0都做过了
<p><p>/BATCH
! /COM,ANSYS RELEASE 14.0 UP20111024 11:37:55 09/10/2012
/input,start140,ans,'D:\Program Files\ANSYS Inc\v140\ANSYS\apdl\',,,,,,,,,,,,,,,,1
/PREP7
! /PNUM,VOLU,ON
ET,1,SHELL181
ET,2,SOLID186
MP,EX,1,2E11
MP,DENS,1,7800
MP,NUXY,1,0.3
MP,EX,2,2E11
MP,DENS,2,7800
MP,NUXY,2,0.3
! 材料2双线性随动强化
TB,BKIN,2,1
TBTEMP,0
TBDATA,,240E6,2E8
R,1,0.02 ! 疑问1：这个实数的定义似乎是后面定义shell181的厚度，不过似乎这种方法是没有用的？现在应该是都用section了？
! 下面是画出三个轧辊
CYLIND,0.38/2,0,0.3,2.3,0,360
CYLIND,0.38/4,0,0.3,0,0,360
CYLIND,0.38/4,0,2.6,2.3,0,360
VGLUE,ALL
VSEL,NONE
WPOFF,0.54/2,-0.4
CYLIND,0.3/2,0,0.3,2.3,0,360
CYLIND,0.3/4,0,0.3,0,0,360
CYLIND,0.3/4,0,2.6,2.3,0,360
VGLUE,ALL
VSEL,NONE
WPOFF,-0.54
CYLIND,0.3/2,0,0.3,2.3,0,360
CYLIND,0.3/4,0,0.3,0,0,360
CYLIND,0.3/4,0,2.6,2.3,0,360
VGLUE,ALL
! 因为轧辊需要转动，我看到的一个文献上解释说：solid没有办法转动，只能用shell，所以这里删掉了volume只保留了area
ALLS
VDELE,ALL
! 定义硬点
HPTCREATE,AREA,5,,COORD
HPTCREATE,AREA,10,,COORD,0,0,2.6
HPTCREATE,AREA,34,,COORD,-0.54/2,-0.4,2.6
HPTCREATE,AREA,29,,COORD,-0.54/2,-0.4,0
HPTCREATE,AREA,22,,COORD,0.54/2,-0.4,2.6
HPTCREATE,AREA,17,,COORD,0.54/2,-0.4,0
! 定义网格划分参数并划分网格，用了三角网格
MAT,1
TYPE,1
REAL,1
ESIZE,0.05
SMRTSIZE,8
MSHAPE,1
AMESH,ALL
! 画出被卷的钢板并进行网格划分
BLOCK,0.7,-1.0,0.3/2+.02,0.3/2+.02+0.03,0.5,2.1
LESIZE,96,,,1
MAT,2
TYPE,2
MSHKEY,1
MSHAPE,0
ESIZE,0.05
VMESH,ALL
! 定义4，5号单元用来接触，疑问2：漏了编号3，不过似乎没有什么影响
ET,4,TARGE170
ET,5,CONTA174
! 定义conta174的12号关键字为rough
KEYOPT,5,12,1
! 定义三个real，主要为了设置刚度系数为0.1
R,10,,,0.1
R,11,,,0.1
R,12,,,0.1
! 选择40号面，应该是钢板的上表面，并进一步选择了这个上表面除开z坐标0.8~1.8范围内的节点
! 相当于只是选择了钢板的两侧范围内的节点，下面还要选择轧辊上面两侧的节点，创建接触对
ASEL,S,,,40
NSLA,S,1
NSEL,U,LOC,Z,0.8,1.8
MAT,2
REAL,10
TYPE,4
ESURF
ALLS
ASEL,S,,,3,4,1
NSLA,S,1
NSEL,U,LOC,Z,0.9,1.7
NSEL,U,LOC,Z,0.3,0.4
NSEL,U,LOC,Z,2.2,2.3
MAT,1
REAL,10
TYPE,5
ESURF
ALLS
ASEL,S,,,39
NSLA,S,1
NSEL,U,LOC,Z,0.8,1.8
MAT,2
REAL,11
TYPE,4
ESURF
ALLS
ASEL,S,,,27,28,1
NSLA,S,1
NSEL,U,LOC,Z,0.9,1.7
NSEL,U,LOC,Z,0.3,0.4
NSEL,U,LOC,Z,2.2,2.3
MAT,1
REAL,11
TYPE,5
ESURF
ALLS
ASEL,S,,,39
NSLA,S,1
NSEL,U,LOC,Z,0.8,1.8
MAT,2
REAL,12
TYPE,4
ESURF
ALLS
ASEL,S,,,15,16,1
NSLA,S,1
NSEL,U,LOC,Z,0.9,1.7
NSEL,U,LOC,Z,0.3,0.4
NSEL,U,LOC,Z,2.2,2.3
MAT,1
REAL,12
TYPE,5
ESURF
ALLS
FINI
! 设置求解载荷及各项参数
/SOLU
D,ALL,UZ
D,ALL,ROTX
D,ALL,ROTY
NSEL,S,LOC,X
NSEL,R,LOC,Y
D,ALL,UY,-0.08
D,ALL,UX
D,ALL,UZ
ALLS
CSYS,4
ASEL,S,LOC,Y,-2,0.3/2+0.01
NSLA,S,1
D,ALL,ALL
ALLS
ANTYPE,TRANS
NEQIT,100
LNSRCH,ON
NLGEOM,ON
TIME,1
AUTOT,ON
NSUBST,30,80,20
KBC,0
OUTRES,ALL,ALL
SOLVE
!以下模拟下辊转动，驱动钢板前进，但是因为我第一步都做不过，所以以下就根本没有做了
CSYS,4
ASEL,S,LOC,Y,-2,0.3/2+0.01
NSLA,S,1
DDELE,ALL,ALL
ALLS
D,ALL,UZ
D,ALL,ROTX
D,ALL,ROTY
NSEL,S,LOC,X,0
NSEL,R,LOC,Y
D,ALL,ROTZ,-3.1415926
D,ALL,UX
D,ALL,UY
ALLS
WPOFF,0.54
NSEL,S,LOC,X,0
NSEL,R,LOC,Y
D,ALL,UX
D,ALL,UY
D,ALL,ROTZ,-3.1415926
TIME,2
ALLS
KBC,0
NSUBST,200,1000,100
SOLVE
! /COM,ANSYS RELEASE 15.0 UP20131014 08:33:41 03/04/2016
! /GRA,POWER
! /GST,ON
! /PLO,INFO,3
! /GRO,CURL,ON
! /CPLANE,1
! /REPLOT,RESIZE
WPSTYLE,,,,,,,,0
! LGWRITE,'lizi-17-1','lgw','C:\Users\CeShi\Desktop\',COMMENT
有什么其他信息需要的可以在下面留言~我继续添加和解释 更新：查出了一个问题，钢板和轧辊没有接触上，留有的空隙导致了刚体位移，把CONTA174的关键字5设为1（自动闭合），情况好了一点，但是运算到一多半仍然出现不收敛的情况
页:
[1]