24小时热门版块排行榜

返回列表

proejhy

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2469.3
帖子: 97
在线: 69.4小时
虫号: 2467919
注册: 2013-05-16
专业: 传动机械学

[求助] 为什么三辊卷板的这个例子在我的机器里总也不收敛呢？已有1人参与

非线性分析综合应用实例—钢板卷制成圆筒
/CLEAR !清除数据库，新建文件
/FILNAM, EXAMPLE20 !指定工作名
/CONFIG, NRES,2000 !设置最大子步数

!前处理
/PREP7 !进入前处理器
/PNUM, VOLU,ON !打开体号
ET,1,SHELL43 !选择单元类型，壳单元用于划分上下辊
ET,2,SOLID186 !实体单元用于划分钢板
MP,EX,1,2E11 !定义材料模型1的弹性模量
MP,DENS,1,7800 !定义材料模型1的密度
MP,NUXY,1,0.3 !定义材料模型1的泊松比
MP,EX,2,2E11 !定义材料模型2的弹性模量
MP,DENS,2,7800 !定义材料模型2的密度
MP,NUXY,2,0.3 !定义材料模型2的泊松比
TB,BKIN,2,1
TBTEMP,0
TBDATA, ,240E6,2E8 !定义材料模型2的屈服极限、切向模量
R,1,0.02 !定义实常数1，壳单元的厚度
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
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    !指定线96（钢板厚度）被划分为1段
MAT,2 !为钢板划分单元指定属性
TYPE,2
MSHKEY, 1 !映射网格
MSHAPE,0 !指定单元形状为六面体
ESIZE,0.05 !指定单元边长度
VMESH,ALL !对块划分网格
ET,4,TARGE170    !指定单元类型，用于创建接触对
ET,5,CONTA174
KEYOPT,5,12,1 !设置单元5接触表面无滑动
R,10,,,0.1 !定义实常数，KFN=0.1
R,11,,,0.1
R,12,,,0.1
ASEL,S,,,40 !在钢板上表面和上辊表面建立接触对
NSLA,S,1
NSEL,U,LOC,Z,0.8,1.8
MAT,3
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,3
REAL,10
TYPE,5
ESURF
ALLS
ASEL,S,,,39 !在钢板下表面和左下辊表面建立接触对
NSLA,S,1
NSEL,U,LOC,Z,0.8,1.8
MAT,4
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,4
REAL,11
TYPE,5
ESURF
ALLS
ASEL,S,,,39 ! 在钢板下表面和右下辊表面建立接触对
NSLA,S,1
NSEL,U,LOC,Z,0.8,1.8
MAT,4
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,4
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 !施加约束，上辊下降0.08 m
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 !解答
!模拟卸载
ACEL, 0,9.8    !施加重力加速度
CSYS,0 !激活全球直角坐标系
NSEL,S,LOC,X !选择上辊轴线上所有节点
NSEL,R,LOC,Y
D,ALL,UY,-0.06 !上辊向上移动，卸载
D,ALL,UX
D,ALL,UZ
D,ALL,ROTX
D,ALL,ROTY
ALLS
KBC,0
NSUBST,20,300,10
TIME,2.1
SOLVE
FINI
!查看结果
/POST1
SET,LAST !读最后载荷子步计算结果
PLNS,U,SUM
ANTIME,50,0.5, ,1,2,0,2.1 !用动画查看卷制圆筒过程中的变形情况
FINI

回复此楼