finish
/clear,start
/title, Torsionstest


/prep7

!Geometrie

cylind,0,10.2,-20,20,0,360
cylind,10,20,-10,10,0,360

/pnum,area,1
aplot


!Material

et,1,solid45
mp,ex, 1, 500   	!E-Modul für Material1
mp,nuxy, 1, 0.3  	!nuxy für Material 1

et,2,solid45
mp,ex, 2, 10000   	!E-Modul für Material2
mp,nuxy, 2, 0.3  	!nuxy für Material 2


et,3,MPC184 
KEYOPT,3,1,1		!rigid beam


!Vernetzung

mat,1
mshkey, 0			!free meshing
mshape,1
ESIZE,4
type,1				
vmesh,1

mat,2
mshkey, 0			!free meshing
mshape,1
ESIZE,4
type,1				
vmesh,2

!Kontaktpaarung

/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,1,0.8 
MAT,1   
R,3 
REAL,3  
ET,4,170
ET,5,174
KEYOPT,5,5,1
KEYOPT,5,9,0
KEYOPT,5,10,1   
R,3,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,3  
ASEL,A,,,4  
CM,_TARGET,AREA 
TYPE,4  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
ASEL,S,,,9  
ASEL,A,,,10 
CM,_CONTACT,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
/VIEW,  1,  0.565220912931    ,  0.518141951051    ,  0.641914509999
/ANG,   1,   12.4536590662  
/REPLO  
/VIEW,  1,  0.365525665943    ,  0.383593816112    ,  0.848084177295
/ANG,   1,   4.06813160951  
/REPLO  

	

!Randbedingungen

esel,s,mat,,1
nsle,s,1
nsel,r,loc,z,-15,-30
d,all,all
alls

esel,s,mat,,1
nsle,s,1
nsel,r,loc,z,15,30
d,all,all
alls

esel,s,mat,,2
nsle,s,1
d,all,uz
alls

!Solution

/solu

!Uebermass wegrechnen

nlgeom,on
kbc,0
nsub,5,10,2
time,1
neqit,100
nrop,unsymm
alls

solve

!Erzeugung Beamelemente

/prep7

n,100000,0,0,0	
csys,1
nsel,s,loc,x,20

!Beamelemente erzeugen

type,3
*get,nanz,node,,count   	!Anzahl der Knoten

narb=0   	        	!Arbeits-Knoten
*do,i,1,nanz		 	!Schleife mit Nanz durchläufen  

narb=ndnext(narb)        	
e,100000,narb
*enddo

alls
finish

!Momenteinpraegung

/solu
F,100000,MZ,100000
gplot
kbc,0
nsub,5,10,2
neqit,100
nrop,unsymm
alls

solve