在结构中求热应力

kidfly

我在这个长方体上（0.004*0.004*0.05m）EX7.2e10,，PRXY0.3，APLX2.4e-5。全是国际单位
施加了20摄氏度的参考温度以及520度的温度载荷，然后我限制了长方体Z轴方向两个面位移为0，限制了X轴方向两个面位移为1e-6，限制了Y轴方向两个面位移为1e-6！而实际上，长方体自由膨胀的话，Z轴方向会有6e-4的位移，在截面方向也会有4.8e-5的位移！所以我想这样的话长方体应该由于膨胀受阻会有很大的应力产生，但是模拟结果却只有0.277e8，这个值比我只限制Z轴方向算得的结果0.864e9还小很多！真是不知道错哪了，大家帮忙看看吧，谢谢了！
图片在附件里
第一张是加载完后的图
第二张是等效应力图，位移有些夸张，没用真实的比例！
第三张是位移图
/FILNAME,stress test,0  
/TITLE,test
/PREP7  
!*  
ET,1,SOLID45
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,7.2e10
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,1,REFT,,,20
MPDATA,ALPX,1,,2.4e-5   
FINISH  
/SOL
FINISH  
/PREP7  
BLOCK,,0.004,,0.004,,0.05,  
/USER,  1   
/VIEW,  1,  0.562065102477    , -0.101562903645E-01,  0.827030634465
/ANG,   1, -0.146174033205  
/REPLO  
/AUTO,1
/REP,FAST   
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ESIZE,0.001,0,  
MSHAPE,0,3D
MSHKEY,1
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
TUNIF,520,  
TREF,20,
SAVE
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
!*  
/GO
DA,P51X,UZ,
APLOT   
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
/USER,  1   
/VIEW,  1,  0.955138723028    , -0.259023894253E-01, -0.295023873602
/ANG,   1,  -4.72666555355  
/REPLO  
/VIEW,  1, -0.809364904384    ,  0.630827492704E-01, -0.583908398891
/ANG,   1,  -3.84128444134  
/REPLO  
/VIEW,  1, -0.668360624291    , -0.663734407485E-01, -0.740870192585
/ANG,   1,  -3.49732276154  
/REPLO  
/VIEW,  1, -0.772552880947    ,  0.256674967239    , -0.580758131525
/ANG,   1,  -2.66902155023  
/REPLO  
/VIEW,  1, -0.707759526500    ,  0.229998736482    , -0.667964844783
/ANG,   1,  -4.18490168691  
/REPLO  
FLST,2,1,5,ORDE,1   
FITEM,2,3   
!*  
/GO
DA,P51X,UY,1e-6
FLST,2,1,5,ORDE,1   
FITEM,2,4   
!*  
/GO
DA,P51X,UY,1e-006   
FLST,2,1,5,ORDE,1   
FITEM,2,4   
FLST,2,1,5,ORDE,1   
FITEM,2,3   
DADELE,P51X,UY  
FLST,2,1,5,ORDE,1   
FITEM,2,3   
!*  
/GO
DA,P51X,UY,-1e-006  
FLST,2,1,5,ORDE,1   
FITEM,2,5   
!*  
/GO
DA,P51X,UX,-1e-006  
FLST,2,1,5,ORDE,1   
FITEM,2,6   
!*  
/GO
DA,P51X,UX,1e-006   
ALLSEL,ALL  
/SHOW,JPEG,,0   
JPEG,QUAL,75,   
JPEG,ORIENT,HORIZ   
JPEG,COLOR,2
JPEG,TMOD,1
/GFILE,800,
!*  
/CMAP,_TEMPCMAP_,CMP,,SAVE  
/RGB,INDEX,100,100,100,0
/RGB,INDEX,0,0,0,15
/REPLOT
/CMAP,_TEMPCMAP_,CMP
/DELETE,_TEMPCMAP_,CMP  
/SHOW,CLOSE
/DEVICE,VECTOR,0
!*  
ALLSEL,ALL  
solve   
FINISH  
/POST1  
!*  
/EFACET,1   
PLNSOL, S,EQV, 2,1.0
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
/DIST,1,1.08222638492,1
/REP,FAST   
!*  
/DSCALE,ALL,1.0
/EFACET,1   
PLNSOL, S,EQV, 2,1.0
!*  
/DSCALE,ALL,AUTO
/EFACET,1   
PLNSOL, S,EQV, 2,1.0
/SHOW,JPEG,,0   
JPEG,QUAL,75,   
JPEG,ORIENT,HORIZ   
JPEG,COLOR,2
JPEG,TMOD,1
/GFILE,800,
!*  
/CMAP,_TEMPCMAP_,CMP,,SAVE  
/RGB,INDEX,100,100,100,0
/RGB,INDEX,0,0,0,15
/REPLOT
/CMAP,_TEMPCMAP_,CMP
/DELETE,_TEMPCMAP_,CMP  
/SHOW,CLOSE
/DEVICE,VECTOR,0
!*  
!*  
/EFACET,1   
PLNSOL, U,SUM, 2,1.0
/SHOW,JPEG,,0   
JPEG,QUAL,75,   
JPEG,ORIENT,HORIZ   
JPEG,COLOR,2
JPEG,TMOD,1
/GFILE,800,
!*  
/CMAP,_TEMPCMAP_,CMP,,SAVE  
/RGB,INDEX,100,100,100,0
/RGB,INDEX,0,0,0,15
/REPLOT
/CMAP,_TEMPCMAP_,CMP
/DELETE,_TEMPCMAP_,CMP  
/SHOW,CLOSE
/DEVICE,VECTOR,0
!*  
FINISH  
/PREP7  
!*  
/REPLOT,RESIZE  
FINISH  
! /EXIT,MODEL   
这个我还用梁单元beam189算过，位移约束越多，应力却变小了！真是搞不懂！