找回密码
 注册
Simdroid-非首页
查看: 221|回复: 10

[接触分析] 求助:ANSYS中关于cohesive 单元是怎么定义怎么设置的?

[复制链接]
发表于 2011-8-5 18:59:55 | 显示全部楼层 |阅读模式 来自 陕西西安
ANSYS中好像也有cohesive单元,但是如何定义,如何设置的,请教一下,尽量详细具体一点,谢谢!!
发表于 2011-8-7 21:32:23 | 显示全部楼层 来自 日本
Simdroid开发平台
cohesive有两种,一种用inter单元做,一种用contact单元做。不知到楼主想要哪一种。
ansys的帮助里有说明,两者基于的理论不一样。路径如下:
// Theory Reference // 4. Structures with Material Nonlinearities // 4.13. Cohesive Zone Material Model
inter单元我当时试过,发现两个实体单元不会完全分离,回位后就变回原来的状态。而且貌似拉压都按cohesive准则来算,所以只能模拟分离的情况,实体单元如果受压会发生单元重叠。
contact单元则不会出现上述的情况。
另,inter单元要czmesh命令分网,而且这两种单元在定义材料参数是用tb和tbdata来定义的。所以貌似只能用命令操作,界面是无法操作的。
回复 1 不支持 0

使用道具 举报

发表于 2011-11-15 06:09:10 | 显示全部楼层 来自 加拿大
c_shn 发表于 2011-8-7 21:32
cohesive有两种,一种用inter单元做,一种用contact单元做。不知到楼主想要哪一种。
ansys的帮助里有说明, ...

我用inter单元做了,仔细阅读了TBDATA, 但是材料属性一直附加不到单元上去,你是怎么做到的?能否指教?
回复 不支持

使用道具 举报

发表于 2011-11-16 13:33:53 | 显示全部楼层 来自 日本
wangmingdi 发表于 2011-11-15 06:09
我用inter单元做了,仔细阅读了TBDATA, 但是材料属性一直附加不到单元上去,你是怎么做到的?能否指教? ...

楼主是用命令做的吗?我把我做的一个简单的例子贴上来,看看对你有没有帮助。

FINISH
/CLEAR,START
/TITLE,`PRACTICE`
/PREP7
LOCAL,45,,0,0,0,45,
LOCAL,315,,0,0,0,-45,
LOCAL,90,,0,0,0,90


ET,50,MESH200,7
ET,51,MESH200,11

MP,EX,1,156300
MP,EY,1,8900
MP,EZ,1,8900
MP,PRXY,1,0.34
MP,PRXZ,1,0.34
MP,PRYZ,1,0.45
MP,GXY,1,5600
MP,GXZ,1,5600
MP,GYZ,1,3070
TB,CZM,2
TBDATA,1,20,0.025,0.0086

!make the 2Dmodel

TYPE,50
RECTNG,0,1,0,1
ESIZE,1        !control the size of the elements
MSHKEY,1
AMESH,ALL  

!make the 3Dmodel

ASEL,S,,,ALL,,,0
CM,A1,AREA
TYPE,51
ESIZE,1                   !control the thickness
VEXT,A1,,,0,0,2,        !control the thickness
VSWEEP,ALL,,,
CMSEL,S,A1
ACLEAR,ALL

!add the element type


ET,1,SOLID186
ET,2,SOLID186
ET,102,INTER204


!change the element type

*DO,i,1,2
   ALLSEL,ALL
   NSEL,S,LOC,Z,i-1
   NSEL,A,LOC,Z,i
   NSEL,A,LOC,Z,i-0.5
   ESLN,R,1
   EMODIF,ALL,TYPE,i
   EMODIF,ALL,MAT,1
*ENDDO


!create and mesh an interface area by cohesive zone elements.

ALLSEL,ALL
NSEL,S,LOC,Z,1
ESEL,S,TYPE,,1
ESLN,R
CM,TARG,ELEM
ALLSEL,ALL
NSEL,S,LOC,Z,1
ESEL,S,TYPE,,2
ESLN,R
CM,CONT,ELEM
ALLSEL
MAT,2
TYPE,102
CZMESH,TARG,CONT,,,

NUMCMP,ELEM
NUMCMP,NODE

/SOLU
ANTYPE,STATIC
OUTRES,ALL,LAST
ALLSEL,ALL
CSYS,0
NSEL,S,LOC,Z,0
D,ALL,UZ,0
NSEL,S,LOC,X,0
D,ALL,UX,0
NSEL,S,LOC,Y,0
D,ALL,UY,0

ALLSEL,ALL

tim=0


*DO,i,1,200
tim=tim+0.005
TIME,tim
NSEL,S,LOC,Z,2
D,ALL,UZ,0.005*i
OUTRES,ALL,ALL
ALLSEL,ALL
SOLVE

*ENDDO


*DO,i,1,199
tim=tim+0.005
TIME,tim
NSEL,S,LOC,Z,2
D,ALL,UZ,0.95-0.005*i
OUTRES,ALL,ALL
ALLSEL,ALL
SOLVE

*ENDDO



FINISH
/POST1
回复 不支持

使用道具 举报

发表于 2011-11-22 02:32:50 | 显示全部楼层 来自 加拿大
c_shn 发表于 2011-11-16 13:33
楼主是用命令做的吗?我把我做的一个简单的例子贴上来,看看对你有没有帮助。

FINISH

多谢啦,你的例子怎么没有应力变化?
回复 不支持

使用道具 举报

发表于 2011-11-22 03:13:34 | 显示全部楼层 来自 加拿大
c_shn 发表于 2011-11-16 13:33
楼主是用命令做的吗?我把我做的一个简单的例子贴上来,看看对你有没有帮助。

FINISH

施加的载荷是如何实现随时间变化的?我刚刚接触,不太懂,请多指教,qq:25624092
回复 不支持

使用道具 举报

发表于 2011-11-22 16:39:09 | 显示全部楼层 来自 日本
wangmingdi 发表于 2011-11-22 03:13
施加的载荷是如何实现随时间变化的?我刚刚接触,不太懂,请多指教,qq:25624092 ...

楼主你是问这个例子的应力不随时间变化还是指应力分布没有变化?这个例子只有两个单元,而且只是加载了一个随时间变化的单向唯一载荷,也就是只有拉伸和压缩。所以不存在应力分布的变化。应力是随位移载荷变化而变化的,只是到最后两单元距离很大时,相当于分离状态。应力约为0,也就相当于没有变化。你可以通过main menu-〉general postproc-〉result viewer查看。载荷随时间变化有几种方法做,但我只用循环命令做,这个你可以问问其它人。就是这段
*DO,i,1,200
tim=tim+0.005
TIME,tim             !设时间点
NSEL,S,LOC,Z,2
D,ALL,UZ,0.005*i !设该时间点的载荷大小
OUTRES,ALL,ALL
ALLSEL,ALL
SOLVE
*ENDDO
回复 不支持

使用道具 举报

发表于 2011-11-24 00:39:19 | 显示全部楼层 来自 加拿大
c_shn 发表于 2011-11-22 16:39
楼主你是问这个例子的应力不随时间变化还是指应力分布没有变化?这个例子只有两个单元,而且只是加载了一 ...

哦,因为我算了一下,你的例子里面应力不会变化,不知道是不是我理解错误?
回复 不支持

使用道具 举报

发表于 2011-11-24 00:46:41 | 显示全部楼层 来自 加拿大
c_shn 发表于 2011-11-22 16:39
楼主你是问这个例子的应力不随时间变化还是指应力分布没有变化?这个例子只有两个单元,而且只是加载了一 ...

或者您能否给我一个比较详细的例子,多谢!
回复 不支持

使用道具 举报

发表于 2011-11-26 00:42:38 | 显示全部楼层 来自 日本
wangmingdi 发表于 2011-11-24 00:46
或者您能否给我一个比较详细的例子,多谢!

那试试这个吧,也是以前自己做练习留下的。不过很久之前了,是不是这个已经忘了,可能会出错。而且计算量挺大,要算很久。

FINISH
/CLEAR,START
/TITLE,`PRACTICE`
/PREP7
LOCAL,45,,0,0,0,45,
LOCAL,315,,0,0,0,-45,
LOCAL,90,,0,0,0,90


ET,50,MESH200,7
ET,51,MESH200,11

MP,EX,1,156300
MP,EY,1,8900
MP,EZ,1,8900
MP,PRXY,1,0.34
MP,PRXZ,1,0.34
MP,PRYZ,1,0.45
MP,GXY,1,5600
MP,GXZ,1,5600
MP,GYZ,1,3070
TB,CZM,2
TBDATA,1,20,0.025,0.0086

!make the 2Dmodel

TYPE,50
RECTNG,0,1,0,1
ESIZE,1        !control the size of the elements
MSHKEY,1
AMESH,ALL  

!make the 3Dmodel

ASEL,S,,,ALL,,,0
CM,A1,AREA
TYPE,51
ESIZE,1                   !control the thickness
VEXT,A1,,,0,0,2,        !control the thickness
VSWEEP,ALL,,,
CMSEL,S,A1
ACLEAR,ALL

!add the element type


ET,1,SOLID186
ET,2,SOLID186
ET,102,INTER204


!change the element type

*DO,i,1,2
   ALLSEL,ALL
   NSEL,S,LOC,Z,i-1
   NSEL,A,LOC,Z,i
   NSEL,A,LOC,Z,i-0.5
   ESLN,R,1
   EMODIF,ALL,TYPE,i
   EMODIF,ALL,MAT,1
*ENDDO


!create and mesh an interface area composed of cohesive zone elements.

ALLSEL,ALL
NSEL,S,LOC,Z,1
ESEL,S,TYPE,,1
ESLN,R
CM,TARG,ELEM
ALLSEL,ALL
NSEL,S,LOC,Z,1
ESEL,S,TYPE,,2
ESLN,R
CM,CONT,ELEM
ALLSEL
MAT,2
TYPE,102
CZMESH,TARG,CONT,,,

NUMCMP,ELEM
NUMCMP,NODE

/SOLU
ANTYPE,STATIC
OUTRES,ALL,LAST
ALLSEL,ALL
CSYS,0
NSEL,S,LOC,Z,0
D,ALL,UZ,0
NSEL,S,LOC,X,0
D,ALL,UX,0
NSEL,S,LOC,Y,0
D,ALL,UY,0

ALLSEL,ALL

tim=0

*DO,i,1,50
tim=tim+0.005
TIME,tim
NSEL,S,LOC,Z,2
D,ALL,UZ,i*0.005                   
ALLSEL,ALL


DELTIM,0.00005
OUTRES,ALL,ALL

SOLVE

*ENDDO

FINISH
/POST1
回复 不支持

使用道具 举报

发表于 2011-11-26 04:52:46 | 显示全部楼层 来自 加拿大
c_shn 发表于 2011-11-26 00:42
那试试这个吧,也是以前自己做练习留下的。不过很久之前了,是不是这个已经忘了,可能会出错。而且计算量 ...

多谢,you are very nice
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-29 17:35 , Processed in 0.032890 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表