| 查看: 794 | 回复: 2 | ||
[求助]
为什么三辊卷板的这个例子在我的机器里总也不收敛呢？ 已有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 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有73人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
18255025093
金虫 (小有名气)
- 应助: 14 (小学生)
- 金币: 1372.8
- 红花: 2
- 帖子: 220
- 在线: 58小时
- 虫号: 1737801
- 注册: 2012-04-05
- 性别: GG
- 专业: 机械工程
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
proejhy: 金币+5, ★有帮助 2015-08-30 21:59:14
感谢参与，应助指数 +1
proejhy: 金币+5, ★有帮助 2015-08-30 21:59:14
|
这个例子我也做过，本来就不收敛，一直在收敛线上下震荡，不过还是能看到结果，准确度不敢恭维。书上例子很多都是这样的。 发自小木虫Android客户端 |
2楼2015-08-30 19:05:07
3楼2015-08-30 22:01:28
50