求助:关于混凝土分层浇筑应力场问题
FINISH/CLEAR,all
/filname,suanli-1
/PREP7
!*
ET,1,PLANE55
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,173.9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,0.9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,2400
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,300.98
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,0.935
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2700
RECTNG,-30,30,,12,
RECTNG,-90,-30,-40,0,
RECTNG,-30,30,-40,0,
RECTNG,30,90,-40,0,
*do,i,1,8
k,100000,0,i*1.5,0
k,100001,1,i*1.5,0
k,100002,0,i*1.5,1
kwplane,-1,100000,100001,100002
alls
asbw,all
*enddo
aglue,all
NUMMRG,ALL !
NUMCMP, ALL !
CSYS,0 !
FLST,5,3,5,ORDE,3
FITEM,5,1
FITEM,5,-2
FITEM,5,6
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
FLST,5,4,5,ORDE,3
FITEM,5,3
FITEM,5,-5
FITEM,5,7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 1, , 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,3,5,ORDE,3
FITEM,5,1
FITEM,5,-2
FITEM,5,6
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 2, , 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,10,4,ORDE,5
FITEM,5,1
FITEM,5,3
FITEM,5,-9
FITEM,5,19
FITEM,5,-20
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,4, , , , ,0
!*
FLST,5,8,4,ORDE,8
FITEM,5,11
FITEM,5,-12
FITEM,5,14
FITEM,5,-15
FITEM,5,17
FITEM,5,-18
FITEM,5,21
FITEM,5,-22
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,3, , , , ,0
!*
MSHAPE,0,2D
MSHKEY,1
!*
FLST,5,7,5,ORDE,2
FITEM,5,1
FITEM,5,-7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,5,11,4,ORDE,10
FITEM,5,2
FITEM,5,6
FITEM,5,11
FITEM,5,-12
FITEM,5,14
FITEM,5,-15
FITEM,5,17
FITEM,5,-18
FITEM,5,20
FITEM,5,-22
LSEL,S, , ,P51X
nsll,s,1
sf,all,conv,1440,20
allsel
FLST,5,3,4,ORDE,3
FITEM,5,1
FITEM,5,7
FITEM,5,9
LSEL,S, , ,P51X
nsll,s,1
sf,all,hflux,0
FINISH
allsel
/SOL
NROPT,FULL
ESEL,S,MAT,,2 !
NSLE,S
IC,all,TEMP,20 !
allsel
esel,s,mat,,1
ekill,all
nsle,s
d,all,temp,14 !
ic,all,temp,14
*do,layer,1,8,1
!layer=1
nsel,s,loc,y,(layer*1.5-1.5),layer*1.5
esln,s,1
ddele,all,all
ealive,all
allsel
lsel,s,loc,y,layer*1.5
nsll,s,1
sf,all,conv,1440,20
*do,dday,layer*7-6,layer*7,1
allsel
nsel,s,loc,y,-0.1,0.1
nsel,r,loc,x,-300.1,30.1
*SET,he02,2280*32.4*(dDAY/(0.9+dDAY)-(dDAY-1)/(0.9+dDAY-1))
bf,all,hgen,he02
*do,loop,1,layer,1
*SET,Tday,dday-(loop-1)*7
*SET,sday,dday-loop*7
allsel
nsel,s,loc,y,loop*1.5-1.4,loop*1.5-0.1
*SET,he03,2280*32.4*(TDAY/(0.9+TDAY)-(TDAY-1)/(0.9+TDAY-1))
bf,all,hgen,he03
*if,dday,le,7*loop,then
allsel
nsel,s,loc,y,loop*1.5-0.1,loop*1.5+0.1
*SET,he00,2280*32.4*(TDAY/(0.9+TDAY)-(TDAY-1)/(0.9+TDAY-1))
bf,all,hgen,he00
*elseif,dday,gt,7*loop,and,dday,le,7*loop+3,then
allsel
nsel,s,loc,y,loop*1.5-0.1,loop*1.5+0.1
*SET,he02,2280*32.4*(TDAY/(0.9+TDAY)-(TDAY-1)/(0.9+TDAY-1))
*SET,he01,2280*32.4*(sDAY/(0.9+sDAY)-(sDAY-1)/(0.9+sDAY-1))
bf,all,hgen,0.1*he01+0.9he02
*else
allsel
nsel,s,loc,y,loop*1.5-0.1,loop*1.5+0.1
*SET,he02,2280*32.4*(TDAY/(0.9+TDAY)-(TDAY-1)/(0.9+TDAY-1))
*SET,he01,2280*32.4*(sDAY/(0.9+sDAY)-(sDAY-1)/(0.9+sDAY-1))
bf,all,hgen,0.5*(he02+he01)
*endif
*enddo
allsel
ANTYPE,4
LUMPM,0
TIME,dDAY
autots,0
DELTIM,1, , ,1
kbc,0
tsres,erase
outres,all,all,
solve
*enddo
allsel
lsel,s,loc,y,layer*1.5
nsll,s,1
sfdele,all,conv
sfdele,all,conv
!sfdele,all,conv
allsel
*enddo
allsel
FINISH
/PREP7
!*
ET,1,PLANE42
*do,i,1,56
mp,dens,i,2400
mp,ex,i,4.656e10*(1-exp(-4.02*(i**0.335)))
mp,prxy,i,0.167
mp,alpx,i,1.03e-5
mp,reft,i,25
!acel,0,9.8
*enddo
mp,dens,1000,2700.0
mp,ex,1000,5.0e10
mp,prxy,1000,0.167
mp,alpx,1000,0.8e-5
mp,reft,1000,25
!acel,0,9.8
allsel
FINISH
/SOL
antype,trans
nropt,full
nlgeom,on
lsel,s,loc,x,-90
lsel,a,loc,x,90
lsel,a,loc,y,-40
dl,all,,all,0
allsel
esel,s,mat,,2
mpchg,1000,all,
acel,0,9.8
allsel
esel,s,mat,,1
ekill,all
nsle,s,1
d,all,all,0
allsel
*do,dday,1,56
*do,layer,1,8,1
!*do,dday,layer*7-6,layer*7
*if,dday-(layer-1)*7,eq,1,then
nsel,s,loc,y,layer*1.5-1.5,layer*1.5
ddele,all,all
esln,s,1
ealive,all
acel,0,9.8
*endif
allsel
*do,loop,1,layer
*if,dday,gt,(loop-1)*7,then
*SET,tday,dday-(loop-1)*7
nsel,s,loc,y,loop*1.5-1.5,loop*1.5
esln,s,1
mpchg,tday,all,
allsel
*endif
*enddo
*enddo
ldread,temp,,,dday,,'suanli-1','rth',' '
convtol,f,,0.2,2,200
time,dday
autots,0
deltim,1,,,1
kbc,0
tsres,erase
outres,all,all,
solve
!*enddo
*enddo
上面是我做小例子的一些命令流,前面建模部分都是半命令流形式............
做的是分层浇筑,间隔6天,个人感觉温度场已经基本合格,但是应力场的计算老有问题,不知道问题出在哪,希望有这方面的专家给予指点 本来就不能计算变弹模应力累计的 我看好多人做过类似的～～对于一次性浇筑的
分层的时候，该如何考虑？
采用应力增量法计算考虑徐变的温度应力的有人做过～
页:
[1]