接触问题求教

YANGXHHURT

各位高手，小弟刚刚学习接触，想请教两个问题
（1）ansys瞬态动力中能模拟两杆的的高速撞击吗？
（2）接触面在受力后互相分离是什么原因造成的呢？
请各位高手指教啊

graininear

1）冲击问题可以考虑用ansys中的ls-dyna模块计算
2）受力与原来预想的变形不一样，原因有很多，可能是接触面参数设置的不合理。具体问题具体分析，建议上传命令

YANGXHHURT

graininear 发表于 2013-6-4 12:22
1）冲击问题可以考虑用ansys中的ls-dyna模块计算
2）受力与原来预想的变形不一样，原因有很多，可能是接触 ...

1）模型中要用到20节点的单元所以没有在dyna中模拟
2）接触的设置用的接触向导生成的，主要设置的参数就是FKN，

graininear

YANGXHHURT 发表于 2013-6-4 23:04
1）模型中要用到20节点的单元所以没有在dyna中模拟
2）接触的设置用的接触向导生成的，主要设置的参数就 ...

除了FKN，还有单元的各种关键字定义，都可能会产生影响的

YANGXHHURT

graininear 发表于 2013-6-4 23:59
除了FKN，还有单元的各种关键字定义，都可能会产生影响的

非常感谢你的回复，谢谢斑竹

YANGXHHURT

graininear 发表于 2013-6-4 23:59
除了FKN，还有单元的各种关键字定义，都可能会产生影响的

你好版主，我把命令流贴上，你能帮我看看哪里有问题吗
/PREP7  

ET,1,solid45      !3-D   
MP,EX,1,2.7e10  
MP,nuxy,1,0.2   
MP,dens,1,2500  

BLC4,,,0.25,0.25,5
BLOCK,0,0.25,0,0.25,5,10

esize,0.05        !将边长0.25分为5等分
vmesh,all   

ET,2,TARGE170    !3-D目标单元（指压杆）
ET,3,CONTA173    !3-D试件上节点表面接触单元
r,2              !建立一个空的实常数(硬接触)，以用于建立接触
r,2,,,1.0,0.1,0
rmore,,,1.0e20,0.0,1.0
keyopt,3,4,0
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,1
keyopt,3,10,1
keyopt,3,11,0
keyopt,3,12,4
keyopt,3,2,0

ASEL,S,,,7  !选择目标表面  A7   
NSLA,S,1    !选择面上的所有节点,必须是这样
ESLN,S,0
TYPE,2
real,2
ESURF,all

  

ESEL,S,TYPE,,2    !查看单元类型为2的单元，目的是确认是否为一个面
eplot

ASEL,S,,,2  !选择目标表面  A2   
NSLA,S,1    !选择面上的所有节点,必须是这样
ESLN,S,0
TYPE,3
real,2
ESURF,all


ESEL,S,TYPE,,3    !查看单元类型为3的单元，目的是确认是否为一个面
eplot   

allsel


*dim,force,table,5,1,1              !5行1列
*tread,force,'book1','txt','',1     !跳过第一行，因为文件force.txt的第一行为标识符号
*vplot,force(1,0),force(1,1)


/SOLU
ANTYPE,TRANS
TRNOPT,FULL
!NLGEOM,1     !打开“大变形”开关，此时是进行非线性分析
NLGEOM,0      !关闭“大变形”开关，此时是进行线性分析

OUTRES,ALL,ALL
AUTOTS,0                 !关闭自动时间步长   由于下面*do命令有时间增量
KBC,0                    !斜坡荷载    若采用stepping load结果也是一致的。

tm_start=1e-12     !开始时间必须大于0
tm_end=40e-3       !瞬态结束时间
tm_incr=2e-3      !时间增量

/PSF,PRES,NORM,2,1,1


NSEL,ALL
NSEL,S,LOC,Z,40-0.001,40+0.001  
NPLOT   
CM,nodes,NODE    !选择桩顶节点，建立节点组元nodes

*do,tm,tm_start,tm_end,tm_incr    ! tm循环从tm_start到tm_end,增量是tm_incr
   time,tm                        !时间值
   SF,nodes,PRES,force(tm)        !施加 时间～荷载  并查看施加的压力方向是否正确（压力是指向单元面为正）
   !最好通过节点组元来施加荷载，不能使用在area上施加的方法
   allsel
   solve                          !开始计算
*enddo