在结构中求热应力
我在这个长方体上（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算过，位移约束越多，应力却变小了！真是搞不懂！
页:
[1]