fini
/clear
/DIRECTORY, , FileName, , 
/filname,rod1 
/title,rod1--Statisch-mechanisch (B5)
/units,MPA

!define parameter
!%%%%%%%%%%%%%%%%
radius=2
length=180
emodul=2.1e5
poisson=0.3
density=7.85e-9


/prep7
!create and mesh rod 
!%%%%%%%%%%%%%%%%%%%

csys,1
et,1,solid45
mp,ex,1,emodul 
mp,prxy,1,poisson 
mp,dens,1,density 

!create rod
cyl4,0,0,radius,,,,length

!mesh rod
smrtsize,2
vsweep,all

!create BC
!%%%%%%%%%

!fix rear face
da,2,all

!front face distributed coupling 
et,2,mpc184,1
type,2
nm=100000
n,nm
nsel,s,loc,z
nsel,u,node,,nm
*get,nn,node,,count
*get,na,node,,num,min
ni=na
*do,i,1,nn
e,nm,ni
ni=ndnext(ni)
*enddo
allsel
!front face moment
f,nm,mz,2000


/solu
nlgeom,on
solve

/post1
rsys,1
!/contour,,,-0.87676e-4,,-0.87621e-4
plnsol,u,z