ansys接触代码整理.docx
《ansys接触代码整理.docx》由会员分享,可在线阅读,更多相关《ansys接触代码整理.docx(54页珍藏版)》请在冰豆网上搜索。
ansys接触代码整理
接触分析代码整理
目录
1./TITLE,AnalysisofaAxisContactingaholeinaDisc!
定义标题1
2.非线性分析综合应用实例—冲击3
3.非线性分析综合应用实例—钢板卷制成圆筒5
4.!
本实例的轴为一等直径空心轴,盘为等厚度圆盘,其结构及尺寸如图20.1所示。
9
5:
ansys中contact为什么没有keyopt312
6:
求助一个用Ansys做出来的跟土木建筑结构有关的,13
7:
桩17
8:
在做接触式逆止器的ansys分析,怎么定义接触对?
22
9:
ANSYS不收敛22
10.利用MPC技术对SOLID-SHELL单元进行连接实例—简支梁27
11.SHSD用于壳-实体装配实例An31
1./TITLE,AnalysisofaAxisContactingaholeinaDisc!
定义标题
/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,CONTACTPAIRCREATION-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
!
*Generatethetargetsurface下面创建目标面
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
!
*Generatethecontactsurface下面创建接触面
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
2.非线性分析综合应用实例—冲击
/CLEAR
/FILNAM,EXAMPLE19
/PREP7
ET,1,SOLID45
ET,2,SHELL43
ET,3,COMBIN14
ET,4,CONTA173
ET,5,TARGE170
MP,EX,1,2E11
MP,NUXY,1,0.3
MP,DENS,1,7800
TB,BKIN,1,1
TBTEMP,0
TBDATA,,240E6,2E8
R,1,0.06
R,2,100000
R,3,,,0.5!
法向接触刚度因子
/VIEW,1,1,1,1
K,1,-0.8,0,-1.5
K,2,-0.8,0,1.5
K,3,0.8,0,1.5
K,4,0.8,0,-1.5
A,1,2,3,4
TYPE,2
REAL,1
ESIZE,0.08
AMESH,1
BLOCK,-0.1,0.1,0.001,0.201,-0.1,0.1
TYPE,1
ESIZE,0.05
VMESH,1
N,950,-0.8,-0.5,-1.5
N,951,0.8,-0.5,-1.5
N,952,0.8,-0.5,1.5
N,953,-0.8,-0.5,1.5
TYPE,3
REAL,2
E,1,950
E,60,951
E,40,952
E,2,953
NSEL,S,,,950,953,1
D,ALL,ALL
ALLSEL
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,-0.3,0.3
NSEL,R,LOC,Z,-0.3,0.3
REAL,3
TYPE,5
ESURF
ALLS!
alls,lclear,all,来清理网格,再重新划分单元
ASEL,S,loc,y,0.001
NSLA,S,1!
nsla,type,nkey:
选择与选中面相关的节点
REAL,3
TYPE,4
ESURF
ALLS
FINI
/SOLU
ANTYPE,TRANS
ACEL,0,9.8
VSEL,S,,,1
NSLV,S,1
IC,ALL,UY,,-8!
IC应该是定义初始条件的一个命令
D,ALL,UX
D,ALL,UZ
ALLS
KBC,0
LNSRCH,ON
AUTOT,ON
OUTRES,ALL,ALL
TIME,0.0013
DELTIM,2.5E-5,1E-6,4E-5
SOLVE
FINI
/POST1
PLNSOL,S,EQV,0,1.0
FINI
/POST26
NSOL,2,880,U,Y
DERIV,3,2,1
DERIV,4,3,1
PLVAR,2
FINI
3.非线性分析综合应用实例—钢板卷制成圆筒
/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.08m
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
4.!
本实例的轴为一等直径空心轴,盘为等厚度圆盘,其结构及尺寸如图20.1所示。
由于模型和载荷都是轴对称的,可以用轴对称方法进行分析。
这里为了后处理时观察结果更直观,我们采用整个模型的四分之一进行建模分析,最后将其进行扩展,来观察整个结构的变形及应力分布、变化情况。
盘和轴用同一种材料,其性质如下:
弹性模量:
EX=2.1E5
泊松比:
NUXY=0.3
接触摩擦系数:
MU=0.2
/TITLE,AnalysisofaAxisContactingaholeinaDisc!
定义标题
/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,CONTACTPAIRCREATION-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,UNSYMUNSYM:
!
完全牛顿拉夫逊法,且允许非对称刚阵指定牛顿拉夫逊法求解的选项
!
*Generatethetargetsurface下面创建目标面
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
!
*Generatethecontactsurface下面创建接触面
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
5:
ansys中contact为什么没有keyopt3
/COM,CONTACTPAIRCREATION-START
/GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,FKN,0.1,0,
!
RMORE,,,1.0E20,0.0,0.5,
!
RMORE,0.0,0,1.0,,1.0,0.5
!
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,2
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
!
Generatethetargetsurface
ASEL,S,,,18
asel,a,,,10
CM,area_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF
ESEL,ALL
!
Generatethecontactsurface
ASEL,S,,,2
ASEL,A,,,6
CM,area_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ALLSEL,ALL
/GRES,cwz,gsav
/COM,CONTACTPAIRCREATION-END
6:
求助一个用Ansys做出来的跟土木建筑结构有关的,
比如梁柱等,在一定荷载作用下产生变形的模型演示的命令流
/prep7
et,1,plane42
et,2,solid45
mp,ex,1,2.5e10!
桩的弹性模量
mp,nuxy,1,0.2!
桩的泊松比
mp,dens,1,2500!
桩的密度
mp,ex,2,2.5e8!
土的弹性模量
mp,nuxy,2,0.4!
土的泊松比
mp,dens,2,2000!
土的密度
tb,dp,2
tbdata,1,19,32,30!
粘聚力c为19,摩擦角为32度,膨胀角为30
RECTNG,0,1,0,8,!
面1
RECTNG,1,5,0,6,!
面2
RECTNG,0,5,0,-16,!
面3
/pnum,area,1
/pnum,line,1
asel,s,,,2,3,1
aglue,all
numcmp,all
allsel
a