注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

htbbzzg的博客

与朋友们分享 CAE 资料和经验

 
 
 

日志

 
 

一个计算不收敛的接触问题命令流 - 请大家帮助分析一下问题在哪里,做个练习  

2011-11-09 09:21:37|  分类: 问与答 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

本例来自一个网友的接触分析课题,他在分析中遇到了问题,请我帮他查找分析失败的原因。这里将它作为一个实例,和大家一起来查找问题所在,希望能对有类似问题的网友起到一点帮助作用。< xmlnamespace prefix ="o" ns ="urn:schemas-microsoft-com:office:office" />

以下是该网友的部分说明:

现在我做了一个接触问题,简而言之,就是一个圆柱体(2)的底面和一个带有球冠顶面的圆柱体(1)的接触问题,试了好久 一直不成功。

带有球冠顶面的圆柱体(1)的下底面为固结  圆柱体(2)的上面施加均匀力,合力为 1000 N

为何接触一直计算不下去呢?每次都是不收敛。

以下是该网友所提供的命令流 (删除了一些无用的语句),感兴趣的朋友们不妨一试。

        我找到的问题会发布到下一篇日志中,请大家先做个练习吧。

/clear
/batch
/units,mks          !  单位系统采用m N
/NERR, 1000000, 9999999, ,  , 

d_1=119.75e-3
h_1=40.6e-3  

r_ding_1=663e-3
!r_ding_1=64.45e-3
!r_min=64.450438731527100e-3
!h_ding_1=2.71e-3

h_2=25.4e-3
d_2=d_1

h_ding_1=r_ding_1-sqrt(r_ding_1**2-(0.5*d_1)**2)


! ******************加刚体的平板***********************!

h_3=0.5*h_1
d_3=d_1

! ******************加刚体的平板完毕***********************!

/prep7
csys,1 !柱坐标
csys,0
k,1,0,0,0
k,2,d_1/2,0,0
k,3,d_1/2,h_1-h_ding_1,0
k,4,0,h_1,0
k,5,d_2/2,h_1,0
k,6,d_2/2,h_1+h_2,0
k,7,0,h_1+h_2,0

k,8,d_3/2,h_1+h_2+h_3,0
k,9,0,h_1+h_2+h_3,0

csys,0

l,1,4
l,2,3
l,1,2
larc,4,3,1,r_ding_1
l,4,5
l,5,6
l,7,6
l,4,7

al,1,2,3,4
al,5,6,7,8

l,7,9
l,6,8
l,9,8
al,7,10,11,9

vrotat,1,,,,,,1,4,90                         ! 旋转面 1,生成体和单元
vrotat,2,,,,,,1,4,90                         ! 旋转面 1,生成体和单元
vrotat,3,,,,,,1,4,90                         ! 旋转面 1,生成体和单元

et,1,SOLID186    
MP,EX,1,2.1e11       
MP,NUXY,1,0.3

et,2,plane42

et,7,SOLID186
mp,ex,7,2.11e30
mp,nuxy,7,0.3

esys,0
meshnum=10
div_factor=1

 

type,1
mat,1

vsel,s,VOLU,,2
LSEL, s,line ,  , 5,7,2
LSEL, a,line ,  , 17
LSEL, a,line ,  , 19
LESIZE, all,  ,  ,meshnum  , div_factor,  , ,  ,

LSEL, s,line ,  , 6,8,2
LSEL, a,line ,  , 18
LESIZE, all,  ,  ,6  , div_factor,  , ,  ,

LSEL, s,line ,  , 20
LSEL, a,line ,  , 21
LESIZE, all,  ,  ,1*meshnum  , 1,  , ,  ,
vmesh, 2

type,7
mat,7
allsel,all,all
vsel,s,VOLU,,3
LSEL, s,line ,  , 11
LSEL, a,line ,  , 22
LSEL, a,line ,  , 24
LESIZE, all,  ,  ,meshnum  , div_factor,  , ,  ,

LSEL, s,line ,  , 9,10
LSEL, a,line ,  , 23
LESIZE, all,  ,  ,4  , 1,  , ,  ,

LSEL, s,line ,  , 25
LSEL, a,line ,  , 26
LESIZE, all,  ,  ,1*meshnum  , 1,  , ,  ,
vmesh, 3

ALLSEL, all, all

NUMCMP,ALL 
NUMMRG,ALL, , , ,LOW

type,1
mat,1

vsel,s,VOLU,,1
LSEL, s,line ,  , 3,4
LSEL, a,line ,  , 12
LSEL, a,line ,  , 14
LESIZE, all,  ,  ,meshnum  , div_factor,  , ,  ,

LSEL, s,line ,  , 1,2
LSEL, a,line ,  , 13
LESIZE, all,  ,  ,6  , 1/div_factor,  , ,  ,


LSEL, s,line ,  , 15
LSEL, a,line ,  , 16
LESIZE, all,  ,  ,1*meshnum  , 1,  , ,  ,
vmesh, 1

!* 
! /COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
! /GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1  
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3 
ET,4,170
ET,5,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,5,4,0
KEYOPT,5,5,0
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,0
KEYOPT,5,10,2  
KEYOPT,5,11,0  
KEYOPT,5,12,0  
KEYOPT,5,2,1
KEYOPT,4,5,0

! Generate the target surface  
ASEL,S,,,8 
CM,_TARGET,AREA
TYPE,4 
NSLA,S,1
ESLN,S,0
ESLL,U 
ESEL,U,ENAME,,188,189  
NSLE,A,CT2 
ESURF  
CMSEL,S,_ELEMCM

! Generate the contact surface 
ASEL,S,,,4 
CM,_CONTACT,AREA
TYPE,5 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF  
ALLSEL 
ESEL,ALL
ESEL,S,TYPE,,4 
ESEL,A,TYPE,,5 
ESEL,R,REAL,,3 
 
ESEL,ALL
ESEL,S,TYPE,,4 
ESEL,A,TYPE,,5 
ESEL,R,REAL,,3 
CMSEL,A,_NODECM
CMDEL,_NODECM  
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM  
CMSEL,S,_KPCM  
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM  
CMSEL,S,_AREACM
CMDEL,_AREACM  
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM  
! /GRES,cwz,gsav 
CMDEL,_TARGET  
CMDEL,_CONTACT 
! /COM, CONTACT PAIR CREATION - END 

!   contact pair end
! /COM, CONTACT PAIR PROPERTIES - START  
RMODIF,3,5,0.2 
KEYOPT,5,2,0
KEYOPT,5,5,3
! /COM, CONTACT PAIR PROPERTIES - END

!* 
fini

/sol

!* 
allsel,all,all
asel,s,AREA,,1
asel,a,AREA,,2
asel,a,AREA,,3
asel,a,AREA,,7
asel,a,AREA,,11
asel,a,AREA,,15
Da, all, SYMM

!* 
allsel,all,all
asel,s,AREA,,6
NSLA,s,1
D, all, all

!* 
allsel
nsel,s,loc,x,0
nsel,r,loc,z,0
d,all,ux
d,all,uz

!* 
!allsel
!nsel,s,loc,z,0
!d,all,uz

allsel,all,all
asel,s,AREA,,14
NSLA,s,1 
!* 
!SFA, all,  , pres, -1*45*1000*9.8/(sqrt(2)*(3.14*(0.5*d_2)**2))
!sf,all,pres,-1*45*1000*9.8/(sqrt(2)*(3.14*(0.5*d_2)**2))
F, all, fy, -45*1000*9.8/(sqrt(2)*(3.14*(0.5*d_2)**2))
allsel,all,all

!* 
ANTYPE,4
!* 
TRNOPT,FULL
LUMPM,0
!* 
NSUBST,50,50,10
OUTRES,ERASE
OUTRES,ALL,5
AUTOTS,1
KBC,0  
PSTRES,1
TIME,1

save

solve


感兴趣的网友可以自己先试试看,下一次我会发布一个参考的解答。
  评论这张
 
阅读(1697)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016