TALK=f;RUN(1,1)
  DISPLAY
  Problem: Torsion of the long thin Beam (Plate).
   1. Comparison - a analytical solution (Timoshenko, Paragraph 108)
   2. Bottom of the beam is fixed (U1=V1=W1 = 0),
      Top - set tangential stress TauYZ for V1.
      
  ENDDIS

 
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
TEXT(3D TORSION of thin Beam; s705)
libref=705
title
  Declarations and settings
REAL(LX,LY,LZ,POISSON,YOUNG)
REAL(Momnt,TAUZ,GG,SHRYZ) 
LX=2.e-3
LZ=50.e-3
LY=10.e-3 

YOUNG   = 1/0.5E-11   ! Young's modulus
POISSON=0.3           ! Poisson's ratio
GG = YOUNG/(2.0*(1.0+POISSON))
  *********************************
Momnt = 10.0 * LX  
 
TAUZ= 3*Momnt/LY/(LX**3)/GG
TAUZ
SHRYZ = 2*GG*TAUZ

  ******* plus  SHRXZ (see Timoshenko, Paragraph 108 !)
SHRYZ = 2*SHRYZ 

SHRYZ 
  *********************************



INTEGER(NXBODY,NYBODY,NZBODY)

 ************************************************************
  Group 2. Time dependence
 STEADY  =    T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 NXBODY = 6
 GRDPWR(X,NXBODY,LX,1)  

 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NYBODY = 8
 GRDPWR(Y,NYBODY,LY,1)  
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZBODY = 12
 GRDPWR(Z,NZBODY,LZ,1)  
  ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(DISX,DISY,DISZ)
 SOLUTN(DISX  ,Y,Y,Y,N,N,Y)
 SOLUTN(DISY  ,Y,Y,Y,N,N,Y)
 SOLUTN(DISZ  ,Y,Y,Y,N,N,Y)
 
 STORE(PRPS,DRH1,VISL)     ! Necessary to save storage DRH1,VISL 

 STORE(STRX,STRY,STRZ,STXY,STXZ,STYZ)
 STORE(EPSY,EPSX,EPSZ)
 STORE(U1T,V1T,U1/T,V1/T)

 ************************************************************
  GROUP 8. ITERATION NUMBERS ETC

 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=1         
  160    7800.0    0.3       473.0   43.0      1.0e-5   0.5e-11 
 
 ************************************************************
  GROUP 11. INITIAL VALUES
fiinit(disx)=0.0
fiinit(disy)=0.0
fiinit(disz)=0.0
 
FIINIT(PRPS)=160

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
char(FormU,FormV)  
FormU=-:TAUZ:*ZG*(YG-0.5*:LY:)
FormV=:TAUZ:*ZG*(XG-0.5*:LX:)
 
PATCH(DOWNZ,LWALL,1,NX,1,NY,1,1,1,1)            ! z=0 - fixed ALL
COVAL(DOWNZ,DISZ,1,0.0)
PATCH(DOWNXY,CELL,1,NX,1,NY,1,1,1,1)
(SOURCE of DISX at DOWNXY is COVAL(FIXVAL,:FormU:)) 
(SOURCE of DISY at DOWNXY is COVAL(FIXVAL,:FormV:)) 


PATCH(UP,HIGH,1,NX,1,NY,NZ,NZ,1,1)            ! UP - Moment of two force
(SOURCE of DISY at UP is COVAL(FIXFLU,:SHRYZ:*(XG-0.5*:LX:))) 


 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =      4500
 ISG21=LSWEEP
 ************************************************************
  GROUP 17. RELAXATION
     #CONPROM
 ************************************************************
  GROUP 19. DATA TRANSMITTED TO GROUND
 STRA    =    T
 PARSOL  =    F
 ISG52   =    3   !   probe & res

 ************************************************************
  GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1   ! graphic-mode
NYPRIN = 1
NXPRIN = 1
NZPRIN = 1
IXMON = NX-2
IYMON = 2
IZMON = NZ-2
      #conprom 

  inform7begin
   **** CALCULATE analytical solution ***
(STORED VAR U1T IS :FormU:)
(STORED VAR V1T IS :FormV:)
(STORED VAR U1/T IS DISX/(U1T+1.e-20))
(STORED VAR V1/T IS DISY/(V1T+1.e-20))

output(drh1,n,n,n,n,n,n)  
output(prps,n,n,n,n,n,n)  
output(visl,n,n,n,n,n,n)  
output(u1t,n,n,n,n,n,n)  
output(v1t,n,n,n,n,n,n)  



  inform7end
    restrt(all)

STOP