FINISH
/CLEAR

/PREP7 
CYL4, , ,1.5, , , ,30 
ET,1,SOLID187 
vmesh,all

MPTEMP,,,,,,,, 
MPTEMP,1,0 
MPDATA,EX,1,,640000
MPDATA,PRXY,1,,0.3 
 
FINISH 
/SOL

LSCLEAR,ALL
ANTYPE,STATIC,NEW
CSYS,0 
ASELE,,loc,z,30
*GET,afest,area,,num,max
DA,afest,all,0
ALLSEL


moment=100000
ALLSEL
NSEL,,LOC,Z,0
*VGET,mvect,NODES,,NLIST
ALLSEL
*GET,manz,PARM,mvect,DIM,1
einzmom=moment/manz
*DO,mvakt,1,manz,1
 *GET,naktx,NODE,mvect(mvakt),LOC,X
 *GET,nakty,NODE,mvect(mvakt),LOC,Y
 betrag2=naktx**2+nakty**2
 fn=einzmom/betrag2
 fnx=fn*nakty
 fny=-fn*naktx
 f,mvect(mvakt),FX,fnx
 f,mvect(mvakt),FY,fny
*ENDDO 

SOLVE

/POST1
ALLSEL
NSEL,NONE
NSEL,A,NODE,,1
NSEL,A,NODE,,2
NSEL,A,NODE,,10
NSEL,A,NODE,,20
NSEL,A,NODE,,30
NSEL,A,NODE,,40

PRNSOL,U



 PRINT U    NODAL SOLUTION PER NODE
 
  ***** POST1 NODAL DEGREE OF FREEDOM LISTING *****                            
 
  LOAD STEP=     1  SUBSTEP=     1                                             
   TIME=    1.0000      LOAD CASE=   0                                         
 
  THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATE SYSTEM  
 
    NODE      UX          UY          UZ          USUM  
       1  -61.895     -39.110      2.9252      73.274    
       2  -59.606     -41.401     -4.4600      72.711    
      10  -60.748     -37.123      5.3210      71.391    
      20  -59.012     -41.322     -5.0646      72.219    
      30  -60.152     -40.043     -1.1132      72.270    
      40  -60.675     -37.830      3.8509      71.606    

 MAXIMUM ABSOLUTE VALUES
 NODE          1           2          10           1
 VALUE   -61.895     -41.401      5.3210      73.274