yan1984 发表于 2008-4-25 23:13

请教接触单元生死问题

面面接触,相同的接触区域,可不可以同时设置两组接触单元,然后在计算的时候利用杀死和激活接触单元单元来控制用哪一组接触单元参与该阶段的计算?

[ 本帖最后由 rodge 于 2008-4-27 16:39 编辑 ]

yan1984 发表于 2008-4-26 00:07

ET,3,169 ET,4,171 KEYOPT,4,2,0 !运算法则选择:扩展拉郎(默认) KEYOPT,4,9,0 KEYOPT,4,10,0 KEYOPT,4,12,0 MP,MU,6,0.364 r,3, RMODIF,3,1,,,kn1,-v, RMODIF,3,12,ks1,co1 MP,MU,7,1.1918 r,4, RMODIF,4,1,,,kn2,-v, RMODIF,4,12,ks2,co2 ! Generate the target surface !定义目标面 LSEL,S,line,,1,6 CM,_TARGET,LINE TYPE,3 mat,6 real,3 NSLL,S,1 !选择已选面上的所有节点 ESLN,s,1 !选择已选节点上的单元(选择外围梁单元) ESURF !覆盖上169单元 /psymb,esys,1 alls ! Generate the contact surface LSEL,S,,,7,12 CM,_CONTACT,LINE TYPE,4 mat,6 real,3 NSLL,S,1 !选择已选面上的所有节点 ESLN,s,1 !选择已选节点上的单元(选择里面梁单元) ESURF,ALL !覆盖上171单元 /psymb,esys,1 ALLSEL /psymb,esys,0 !!第二组接触 ! Generate the target surface !定义目标面 LSEL,S,line,,1,6 CM,_TARGET,LINE TYPE,3 mat,7 real,4 NSLL,S,1 !选择已选面上的所有节点 ESLN,s,1 !选择已选节点上的单元(选择外围梁单元) ESURF !覆盖上169单元 /psymb,esys,1 alls ! Generate the contact surface LSEL,S,,,7,12 CM,_CONTACT,LINE TYPE,4 mat,7 real,4 NSLL,S,1 !选择已选面上的所有节点 ESLN,s,1 !选择已选节点上的单元(选择里面梁单元) ESURF,ALL !覆盖上171单元!上面是我的两组接触单元的定义,希望有人帮我看看错在哪里了?可不可以进行接触单元生死来控制计算过程?在此谢过了!

htbbzzg 发表于 2008-4-26 13:22

除了点对点接触单元外,其它的接触单元多数可以进行生死单元的设置。
不过,你定义的好像不是接触单元,而是表面效应单元。

yan1984 发表于 2008-4-26 16:59

谢谢htbbzzg,接触单元不是这么定义的吗?我要想在两个面(平面模型中是两条曲线)上定义面面接触该怎么定义呢?能不能再提示一下?谁知道帮帮忙吧

yan1984 发表于 2008-4-27 16:27

怎么没人指点一下啊?是不是我上面的命令完全错的,朽木不可雕了啊?

rodge 发表于 2008-4-27 16:39

是这个方法定义的

yan1984 发表于 2008-4-27 23:01

谢谢rodge!
再问一下我上面定义的相互接触的梁单元在建模的时候需不需要考虑它们的厚度因素啊?

如果考虑它们厚度的话,就要建成不重合的线了。这样的话,计算就出错,说发生了刚体位移

rodge 发表于 2008-4-28 00:24

发生刚体位移是因为没有约束好,或者是接触部分发生渗透,与原来的接触点不再接触了,可以设置接触的接触行为来控制

16443 发表于 2008-4-28 09:21

回复 7楼 的帖子

beam和shell单元在接触计算时候都要考虑厚度,可以在几何模型中考虑,也可以在接触参数中设置偏移 。

yan1984 发表于 2008-4-29 00:35

真的是太感激楼上两位大哥的帮助了!
按照我对两位意思的理解,在我原来的基础上,采用了考虑单元厚度建模,并设置了KEYOPI(11)=1和KEYOPI(12)=2。
我的模型是这样的:
相互接触的两个圆环,里面是厚度为45厘米的梁单元,外面是10厘米的梁单元。材料都是混凝土
外面的梁单元与最外层的平面单元相联.
里面的梁单元没有约束。计算过程中梁单元会有指向圆心的均布力
计算后会出现下面的错误:
*** WARNING ***                         CP =      75.984   TIME= 00:11:00
There are 39 small equation solver pivot terms.                        

*** ERROR ***                           CP =      76.047   TIME= 00:11:00
One or more elements have become highly distorted.Excessive         
distortion of elements is usually a symptom indicating the need for   
corrective action elsewhere.Try incrementing the load more slowly   
(increase the number of substeps or decrease the time step size).You
may need to improve your mesh to obtain elements with better aspect   
ratios.Also consider the behavior of materials, contact pairs,      
and/or constraint equations.If this message appears in the first      
iteration of first substep, be sure to perform element shape checking.
两位大哥有时间的话,再帮我指导一下。先谢了!
页: [1]
查看完整版本: 请教接触单元生死问题