Esope source files

   1 : C PODEFO    SOURCE    CHAT      11/03/22    21:15:13     6910           
   2 :       SUBROUTINE PODEFO(V,DEFO,XL,IDEUX,carac,you,xnu)                                        
   3 : C-----------------------------------------------------------------------        
   4 : C                                                                               
   5 : C    ROUTINE DE CALCUL DES DEFORMATIONS DE LA POUTRE DROITE                     
   6 : C                                                                               
   7 : C  ENTREE                                                                       
   8 : C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE          
   9 : C    XL      LONGUEUR                                                           
  10 : C    IDEUX   INDICE POUR DEFORMATIONS DU 2-EME ORDRE                            
  11 : C                                                                               
  12 : C  SORTIE                                                                       
  13 : C    DEFO    DEFORMATIONS AUX 2 NOEUDS    6 PAR NOEUD                           
  14 : C            ( PAS D'EFFET D'EFFORT TRANCHANT || )                              
  15 : C                                                                               
  16 : C-----------------------------------------------------------------------        
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 :       DIMENSION DEFO(*),V(*),carac(*)                                                    
  20 : C                                                                               
  21 :       RL=1.D00/XL                                                               
  22 :       RL2=RL*RL                                                                 
  23 :           X2=0.                                                                 
  24 :           X3=0. 
  25 :           xiny=carac(2)
  26 :           xinz=carac(3) 
  27 :           sredy=carac(5) 
  28 :           sredz=carac(6)
  29 :           cisail= you/(2.*(1+xnu))
  30 : * of
  31 : *          if( sredy.ne.0.) X2=You/cisail*rl2*xiny/sredy 
  32 : *          if( sredy.ne.0.) X2=You/cisail*rl2*xinz/sredz                                                   
  33 :           if( sredy.ne.0.) X3=You/cisail*rl2*xinz/sredy 
  34 :           if( sredz.ne.0.) X2=You/cisail*rl2*xiny/sredz 
  35 : * of                                                  
  36 :              X13=1./(1.+12.*X3)                                                 
  37 :              X33=(1.+3.*X3)*X13                                                 
  38 :              X63=(1.-6.*X3)*X13                                                 
  39 :              X12=1./(1.+12.*X2)                                                 
  40 :              X32=(1.+3.*X2)*X12                                                 
  41 :              X62=(1.-6.*X2)*X12                                                 
  42 : C                                                                               
  43 : C  DEFORMATIONS GENERALISEES                                                    
  44 : C                                                                               
  45 :       EPS  =(V(7)-V(1))*RL                                                      
  46 :       XSIX =(V(10)-V(4))*RL                                                     
  47 :       XSIY1=                                                                    
  48 :      .RL2*( 6.*V(3)*X12-4.*XL*X32*V(5)-6.*V(9)*X12-2.*XL*X62*V(11))             
  49 :       XSIY2=                                                                    
  50 :      .RL2*(-6.*V(3)*X12+2.*XL*X62*V(5)+6.*V(9)*X12+4.*XL*X32*V(11))             
  51 :       XSIZ1=                                                                    
  52 :      .RL2*(-6.*V(2)*X13-4.*XL*X33*V(6)+6.*V(8)*X13-2.*XL*X63*V(12))             
  53 :       XSIZ2=                                                                    
  54 :      .RL2*( 6.*V(2)*X13+2.*XL*X63*V(6)-6.*V(8)*X13+4.*XL*X33*V(12))             
  55 : C                                                                               
  56 : C  REMPLISSAGE DU TABLEAU DEFO                                                  
  57 : C                                                                               
  58 :       DEFO(1)=EPS                                                               
  59 :       DEFO(2)=0.D0                                                              
  60 :       DEFO(3)=0.D0                                                              
  61 :       DEFO(4)=XSIX                                                              
  62 :       DEFO(5)=XSIY1                                                             
  63 :       DEFO(6)=XSIZ1                                                             
  64 :       DEFO(7)=EPS                                                               
  65 :       DEFO(8)=0.D0                                                              
  66 :       DEFO(9)=0.D0                                                              
  67 :       DEFO(10)=XSIX                                                             
  68 :       DEFO(11)=XSIY2                                                            
  69 :       DEFO(12)=XSIZ2                                                            
  70 :       IF(IDEUX.EQ.0) RETURN                                                     
  71 : C                                                                               
  72 : C  TERMES DU SECOND ORDRE                                                       
  73 : C                                                                               
  74 :       AUX1=0.5D0*(EPS**2+V(5)**2+V(6)**2)                                       
  75 :       AUX2=0.5D0*(EPS**2+V(11)**2+V(12)**2)                                     
  76 :       DEFO(1)=DEFO(1)+0.5D0*(AUX1+AUX2)                                         
  77 :       DEFO(7)=DEFO(7)+0.5D0*(AUX1+AUX2)                                         
  78 : C*OF?    DEFO(1)=DEFO(1)+AUX1
  79 : C*OF?    DEFO(7)=DEFO(7)+AUX2                                         
  80 :       RETURN                                                                    
  81 :       END                                                                       
  82 :  
  83 :  
  84 :