24小时热门版块排行榜

返回列表

小雒1231

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 2
在线: 39分钟
虫号: 1781919
注册: 2012-04-27
专业: 交通工程

[求助] 我要做一个接触分析的，算不出结果，求高手帮我看看

我写了两段命令，一个算不出结果，另一个不收敛，求高手帮我看看。下面是命令：
h=6 !.....结构层厚度
h1=0.30 !.....面层厚度
h2=0.30 !.....级配碎石层厚度
h3=h-h1-h2  !.....土基厚度
b=0.189    !..........轮胎尺寸
d=0.32    !..........双轮间距
w=5    !..........模型宽度
FINISH
/prep7
*afun,deg
et,1,plane42
blc4,0,0,w,h3,0
blc4,0,h3,w,h2,0
blc4,0,h2+h3,w,h1,0

wprota,0,0,90
wpave,w/2-d/2-b/2,h,0
asbw,all

wpave,w/2-d/2+b/2,h,0
asbw,all

wpave,w/2+d/2-b/2,h,0
asbw,all

wpave,w/2+d/2+b/2,h,0
asbw,all

asel,s,loc,y,0,h3+h2-0.0001
aglue,all

asel,s,loc,y,h2+h3+0.0001,h
aglue,all

MP,ex,1,1200e6
mp,prxy,1,0.35

mp,ex,2,300e6
mp,prxy,2,0.3
tb,dp,2
TBDATA,1,54000
TBDATA,2,47

mp,ex,4,40e6
mp,prxy,4,0.35

asel,s,loc,y,0,h3
AATT,4,,,,0

asel,s,loc,y,h3,h3+h2
AATT,2,,,,0

asel,s,loc,y,h3+h2,h
AATT,1,,,,0

lsel,s,length,,h1
lesize,all,,,5

lsel,s,length,,h2
LESIZE,all,,,5

lsel,s,length,,h3
lesize,all,,,50

lsel,s,length,,b
lesize,all,,,3

lsel,s,length,,w/2-d/2-b/2
lesize,all,,,20

lsel,s,length,,d-b
lesize,all,,,2

allsel,all
amesh,all

/COM, CONTACT PAIR CREATION - START

MP,MU,1,0.5
MAT,1
!MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,2,169
ET,3,172
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
RMORE,10.0
KEYOPT,3,3,0
KEYOPT,3,4,2
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,6
KEYOPT,3,2,2
! 生成目标面

asel,s,mat,,2
lsla,r
lsel,r,loc,y,h2+h3
CM,_TARGET,LINE
TYPE,2
NSLL,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! 生成接触面
lsel,all
asel,s,mat,,1
lsla,r
lsel,r,loc,y,h2+h3
CM,_CONTACT,LINE
TYPE,3
NSLL,
ESLN,S,0
ESURF,all
ALLSEL,all
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/COM, CONTACT PAIR CREATION - END
/MREP,EPLOT

/solu

!ANTYPE,0

!NLGEOM,ON
!NROPT,FULL,,off

!pred,on
!AUTOTS,1
!TIME,1
!NSUBST,50,,,1
!TIME,0.1
!NSUBST,10,20,5

antype,0
nlgeom,on
autots,on

lsel,s,loc,y,0
d,all,all,0
allsel,all
lsel,s,loc,x,0
lsel,a,loc,x,w
d,all,,ux,0
allsel,all

asel,s,mat,,1
lsla,r
lsel,r,loc,y,h
lsel,r,length,,b
SFL,all,pres,0.7e6
allsel,all
solve

一上命令算不出结果，好像是荷载没加上一样。
第二个
h=6 !.....结构层厚度
h1=0.30 !.....面层厚度
h2=0.30 !.....级配碎石层厚度
h3= 0.30  !.....半刚性基层厚度
h4=h-h1-h2-h3 !.....土基厚度
b=0.189    !..........轮胎尺寸
d=0.32    !..........双轮间距
w=5    !..........模型宽度
FINISH
/prep7
*afun,deg
et,1,plane42
blc4,0,0,w,h4,0
blc4,0,h4,w,h3,0
blc4,0,h4+h3,w,h2,0
blc4,0,h4+h3+h2,w,h1,0

wprota,0,0,90
wpave,w/2-d/2-b/2,h,0
asbw,all

wpave,w/2-d/2+b/2,h,0
asbw,all

wpave,w/2+d/2-b/2,h,0
asbw,all

wpave,w/2+d/2+b/2,h,0
asbw,all

asel,s,loc,y,0,h4+h3+h2-0.0001
aglue,all

asel,s,loc,y,h2+h3+h4+0.0001,h
aglue,all

MP,ex,1,1200e6
mp,prxy,1,0.35

mp,ex,2,300e6
mp,prxy,2,0.3
tb,dp,2
TBDATA,1,5400
TBDATA,2,47

mp,ex,3,1500e6
mp,prxy,3,0.25

mp,ex,4,40e6
mp,prxy,4,0.35

asel,s,loc,y,0,h4
AATT,4,,,,0

asel,s,loc,y,h4,h4+h3
AATT,3,,,,0

asel,s,loc,y,h4+h3,h4+h3+h2
AATT,2,,,,0

asel,s,loc,y,h4+h3+h2,h
AATT,1,,,,0

lsel,s,length,,h1
lesize,all,,,5

lsel,s,length,,h2
LESIZE,all,,,5

lsel,s,length,,h3
lesize,all,,,5

lsel,s,length,,h4
lesize,all,,,50

lsel,s,length,,b
lesize,all,,,3

lsel,s,length,,w/2-d/2-b/2
lesize,all,,,20

lsel,s,length,,d-b
lesize,all,,,2

allsel,all
amesh,all

/COM, CONTACT PAIR CREATION - START

MP,MU,1,0.5
MAT,1
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,2,169
ET,3,172
R,3,,,1.0,0.1,,
NROPT,UNSYM
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
RMORE,10.0
KEYOPT,3,3,0
KEYOPT,3,4,2
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,6
KEYOPT,3,2,2
! 生成目标面

asel,s,mat,,2
lsla,r
lsel,r,loc,y,h2+h3+h4
CM,_TARGET,LINE
TYPE,2
NSLL,S,1
ESLN,S,0
ESURF,all
CMSEL,S,_ELEMCM
! 生成接触面
lsel,all
asel,s,mat,,1
lsla,r
lsel,r,loc,y,h2+h3+h4
CM,_CONTACT,LINE
TYPE,3
NSLL,
ESLN,S,0
ESURF,all
ALLSEL,all
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/COM, CONTACT PAIR CREATION - END
/MREP,EPLOT

/solu

!ANTYPE,0
!NLGEOM,1
!TIME,1
!AUTOTS,1
!NSUBST,20, , ,1
!TIME,0.1
!NSUBST,10,20,5

antype,0
nlgeom,on
autots,on

nsel,s,loc,y,0
d,all,all,0
allsel,all
nsel,s,loc,x,0
nsel,a,loc,x,w
d,all,ux,0
allsel,all

asel,s,mat,,1
lsla,r
lsel,r,loc,y,h
lsel,r,length,,b
SFL,all,pres,0.7e6
solve

这段命令好像不收敛，我对比两段命令，找不出错误，求高手帮我看看，谢谢了。
我是新手，没有太多金币，谢谢了。

回复此楼