!aus dem Forum von jabba90

fini  $/clear
*abbr,jabba,/input,jabba,txt

*SET,r1,200 
*SET,r2,50  
*SET,r3,160 
*SET,d,40  
*SET,h,8
*SET,m2,1000
*SET,emodul,2e5 
*SET,querkon,0.3
/prep7  
/title,Kreisringplatte  
et,1,shell63
mp,ex,1,emodul  
mp,nuxy,1,querkon  
r,1,h  
k,1 
!circle,,r1  
circle,1,r1 
circle,1,r2 
circle,1,r3 
!kplot  
!/pnum,kp,1  
!kplot  
a,2,3,4,5  
! /kplot
!kplot  
a,6,7,8,9  
aplot  
!/pnum,area,1
!aplot  
asba,1,2
!CSLIST,ALL  
csys,1  
!kplot  
circle,10,d/2  
a,14,15,16,17  
!/pnum,area,1
!aplot  
!alist  
!aplot,1 
!aplot,3 
agen,3,1,,,,120 
!aplot  
!CSLIST,ALL  
asba,3,1
!aplot  
asba,5,2
!aplot  
asba,1,4
!lplot  
!/pnum,line,1
!lplot  
esize,20
lsel,s,,,1,4
lesize,all,,,15 
lsel,s,,,5,8
lesize,all,,20  
lsel,s,,,13,24  
lesize,all,,20  
amesh,all    
csys,1  
nrotat,all
  
!-----Lösung------
/SOL
nsel,s,loc,x,r1
d,all,ux  
d,all,uy
d,all,uz
nsel,s,loc,x,r2    
*get,nanz,node,,count  
*SET,mlast,m2*2*3.1416*r2/(4*(nanz-1))  
lsel,s,loc,x,r2    
f,all,my,mlast 
allsel
solve 
!-----Auswertung---
/post1
pldi
plns,s,eqv

!---Pfad mit Bohrung-----
path,mitBohrung,2,
ppath,1,,r1,0,0
ppath,2,,r2,0,0
pdef,vspa,s,eqv  !von Mises
plpath,vspa
!pdef,uz,u,z
!plpath,uz

!---Pfad ohne Bohrung-----
path,ohneBohrung,2,
ppath,1,,-r1,0,0
ppath,2,,-r2,0,0
pdef,vspa,s,eqv  !von Mises
plpath,vspa

!---Pfade betrachten---
/pbc,path,1
plns,s,eqv