/PREP7   
  /TITLE, ELEMENT BIRTH/DEATH IN A FIXED BAR WITH THERMAL LOADING
  ET,1,LINK1          !二维杆单元
  MP,EX,1,30E6        !材料特性
  MP,ALPX,1,.00005
  MP,EX,2,30E6
R,1,1.0
  N,1  
  N,4,10   
  FILL
  E,1,2     
  EGEN,3,1,-1         !生成三个单元
  FINISH   
  /SOLU   
  ANTYPE,STATIC
  D,1,ALL,,,4,3       !固定杆的两端
  TREF,0              !参考温度0度
  TUNIF,100           !均匀温度载荷
  NROPT,FULL
  OUTPR,BASIC,ALL
  SOLVE   
/post1

etable,s1,ls,1

 

这个例子当中，结构没有发生变形却产生了15000的应力，而弹性模量又不是无穷大，这是因为热膨胀产生的变形和压应力产生的结构应变相抵消的结果

 δ＝a△TL＝PL/EA

 P=a△TEA

σ=P/A=a△TE

△T是同一位置点不同时刻的温度差