求助:关于混凝土分层浇筑应力场问题

springxiaozhu

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天,个人感觉温度场已经基本合格,但是应力场的计算老有问题,不知道问题出在哪,希望有这方面的专家给予指点

bytide

本来就不能计算变弹模应力累计的

springxiaozhu

我看好多人做过类似的～～对于一次性浇筑的

分层的时候，该如何考虑？

采用应力增量法计算考虑徐变的温度应力的有人做过～