这么早的帖子了……
过盈的话,没有试过【CNOF设置过盈】的方法,倒是直接画出几何体的尺寸这个方法屡试不爽,最早是在史亚杰和叶先磊的书《ANSYS工程分析软件应用实例》第八章上看到的,APDL我贴在后面吧,感觉接触问题挺难的,有些参数的设置基本上只有靠试……- /TITLE,Analysis of a Axis Contacting a hole in a Disc !定义标题
- /PREP7
- !*
- ET,1,SOLID185 !定义单元类型
- !*
- MP,EX,1,2.1E5 !定义材料属性
- MP,PRXY,1,0.3
- !*
- CYL4,0,0,34,0,100,90,25 !创建四分之一圆环
- CYL4,0,0,25,0,35,90,150
- VGEN, ,2, , , , ,-10, , ,1 !移动轴
- !*
- LESIZE,17, , ,15, , , , ,1 !定义线的分网尺寸
- LESIZE,19, , ,15, , , , ,1
- !*
- LESIZE,18, , ,2, , , , ,1
- LESIZE,20, , ,2, , , , ,1
- !*
- LESIZE,22, , ,20, , , , ,1
- !*
- LESIZE,5, , ,10, , , , ,1
- LESIZE,7, , ,10, , , , ,1
- !*
- LESIZE,6, , ,8, , , , ,1
- LESIZE,8, , ,8, , , , ,1
- !*
- LESIZE,10, , ,3, , , , ,1
- !*
- VSWEEP,ALL !用扫掠方式对创建的体进行网格划分
- !*
- /COM, CONTACT PAIR CREATION - START
- MP,MU,1,0.2 !定义接触摩擦系数
- MAT,1
- R,3 !定义接触实常数
- REAL,3
- ET,2,170 !定义接触单元类型
- ET,3,174
- R,3,,,0.1,0.1,,
- NROPT,UNSYM
- !* Generate the target surface 下面创建目标面
- ASEL,S,,,4
- CM,_TARGET,AREA
- TYPE,2
- NSLA,S,1
- ESLN,S,0
- ESURF,ALL
- !* Generate the contact surface 下面创建接触面
- ASEL,S,,,9
- CM,_CONTACT,AREA
- TYPE,3
- NSLA,S,1
- ESLN,S,0
- ESURF,ALL
- CMDEL,_TARGET
- CMDEL,_CONTACT
- ALLSEL,ALL
- EPLOT
- FINISH
- !*
- /SOLU 进入求解器
- DA,5,SYMM !定义面的对称位移边条
- DA,6,SYMM
- DA,11,SYMM
- DA,12,SYMM
- DA,3,ALL, !定义面的位移约束条件
- !*
- ANTYPE,0 !指定分析类型为静力分析
- NLGEOM,1 !考虑大变形影响
- AUTOTS,0
- TIME,100
- SOLVE !求解第一载荷步
- !*
- NSUBST,150,10000,10
- OUTRES,ALL,ALL
- AUTOTS,1
- TIME,250
- NSEL,S,LOC,Z,140 !选定轴向坐标为140的所有节点
- D,ALL,UZ,40
- ALLSEL,ALL
- SOLVE !求解第二载荷步
- !*
- /EXPAND,4,POLAR,HALF,,90 !进行模型扩展
- /REPLOT
- !*
- /POST1 !进入同样后处理器
- SET,1,LAST,1, !指定查看的载荷步
- PLNSOL,S,EQV,0,1 !查看等效应力的云图
- !*
- SET, , ,1, ,120, ,
- ESEL,S,ENAME,,174
- EPLOT
- PLNSOL,CONT,PRES,0,1
- !*
- PLNS,S,EQV
- ANDATA,0.5, ,1,0,0,1,1,1 !查看动画显示
- !*
- /POST26
- !*
- RFORCE,2,925,F,Z, FZ_2 !定义约束反力变量
- PLVAR,2, !绘制变量-时间曲线
- FINISH
复制代码
|