Esope source files

   1 : C PODEF2    SOURCE    PV        08/09/11    21:16:16     6150           
   2 : C PODEF2    SOURCE    TC1       92/03/16    21:03:26     501                    
   3 :       SUBROUTINE PODEF2(V,DEFO,XL,IDEUX,xinz,sredy,you,xnu)                                        
   4 : C-----------------------------------------------------------------------        
   5 : C                                                                               
   6 : C    ROUTINE DE CALCUL DES DEFORMATIONS DE LA POUTRE DROITE EN 2D                     
   7 : C                                                                               
   8 : C  ENTREE                                                                       
   9 : C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE          
  10 : C    XL      LONGUEUR                                                           
  11 : C    IDEUX   INDICE POUR DEFORMATIONS DU 2-EME ORDRE                            
  12 : C                                                                               
  13 : C  SORTIE                                                                       
  14 : C    DEFO    DEFORMATIONS AUX 2 NOEUDS    6 PAR NOEUD                           
  15 : C            ( PAS D'EFFET D'EFFORT TRANCHANT || )                              
  16 : C                                                                               
  17 : C-----------------------------------------------------------------------        
  18 :       IMPLICIT INTEGER(I-N)
  19 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  20 :       DIMENSION DEFO(*),V(*)                                                    
  21 : C                                                                               
  22 :       RL=1.D00/XL                                                               
  23 :       RL2=RL*RL                                                                 
  24 :           X2=0.                                                                 
  25 :           X3=0.  
  26 :            YOUN=YOU
  27 :            CISAIL=YOU/(2.D0*(1.D0+XNU))
  28 : 
  29 :            IF( SREDY. NE.0.) X3= YOUN/CISAIL*RL2*XINZ/SREDY
  30 :              X13=1./(1.+12.*X3)                                                 
  31 :              X33=(1.+3.*X3)*X13                                                 
  32 :              X63=(1.-6.*X3)*X13                                                 
  33 :              X12=1./(1.+12.*X2)                                                 
  34 :              X32=(1.+3.*X2)*X12                                                 
  35 :              X62=(1.-6.*X2)*X12                                                 
  36 : C                                                                               
  37 : C  DEFORMATIONS GENERALISEES                                                    
  38 : C                                                                               
  39 :       EPS  =(V(4)-V(1))*RL                                                                                               
  40 :       XSIZ1=                                                                    
  41 :      .RL2*(-6.*V(2)*X13-4.*XL*X33*V(3)+6.*V(5)*X13-2.*XL*X63*V(6))             
  42 :       XSIZ2=                                                                    
  43 :      .RL2*( 6.*V(2)*X13+2.*XL*X63*V(3)-6.*V(5)*X13+4.*XL*X33*V(6))             
  44 : C                                                                               
  45 : C  REMPLISSAGE DU TABLEAU DEFO                                                  
  46 : C                                                                               
  47 :       DEFO(1)=EPS                                                               
  48 :       DEFO(2)=0.D0                                                                                                                    
  49 :       DEFO(3)=XSIZ1                                                             
  50 :       DEFO(4)=EPS                                                               
  51 :       DEFO(5)=0.D0                                                                                                                     
  52 :       DEFO(6)=XSIZ2                                                            
  53 :       IF(IDEUX.EQ.0) RETURN                                                     
  54 : C                                                                               
  55 : C  TERMES DU SECOND ORDRE                                                       
  56 : C                                                                               
  57 :       AUX1=0.5D0*(EPS**2+V(3)**2)                                       
  58 :       AUX2=0.5D0*(EPS**2+V(6)**2)                                     
  59 :       DEFO(1)=DEFO(1)+0.5D0*(AUX1+AUX2)                                         
  60 :       DEFO(4)=DEFO(4)+0.5D0*(AUX1+AUX2)                                         
  61 :       RETURN                                                                    
  62 :       END                                                                       
  63 :  
  64 :  
  65 :  
  66 :