/prep7  
/title, Sample application of PSMESH

et,1,92
mp,ex,1,1e7 
mp,alpx,1,1.3e-5
mp,prxy,1,0.30
mp,ex,2,3e7 
mp,alpx,2,8.4e-6
mp,prxy,2,0.30
tref,70 

/foc,,-.09,.34,.42 
/dist,,.99
/ang,,-55.8 
/view,,.39,-.87,.31
/pnum,volu,1 
/num,1
cylind,0.5,,  -0.25,0,  0,180
cylind,0.5,,  1,1.25,   0,180
cylind,0.25,, 0,1,      0,180
wpoff,.05 
cylind,0.35,1, 0,0.75, 0,180
wpoff,-.1 
cylind,0.35,1, 0.75,1, 0,180
wpstyle,,,,,,,,0
vglue,all 
numc,all 
vplot
mat,1 
smrt,off 
vmesh,4,5
mat,2 
vmesh,1,3
/pnum,mat,1 
eplot
psmesh,,example,,volu,1,0,z,0.5,,,,elems   
CM,lines,LINE 
/dist,,1.1 
cmplot
/solu
eqslve,pcg,1e-8
asel,s,loc,y 
da,all,symm 
asel,all
dk,1,ux 
dk,12,ux 
dk,1,uz
sload,1,9,,force,100,1,2
/title,Sample application of PSMESH - preload only  
solve   
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Finally, we construct the actual solution of interest. We want to
!know what happens to the preload in the bolt, and the stress field around
!it, when the assembly temperature rises to 150° F.
!Both the preload and the stresses increase because, for a uniform
!temperature rise, there is greater thermal expansion in the aluminum plates
!than in the steel bolt. Any method for applying preload that did not
!allow the load to change would be unable to predict this result.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/post1
plnsol,s,z
/solu
antype,,restart 
tunif,150    
/title,Sample application of PSMESH - uniform 150°    
solve     

/post1   
plnsol,s,z