Download efcont.eso

Back to the list

   1 : C EFCONT    SOURCE    CHAT      05/01/12    23:28:22     5004           
   2 :       SUBROUTINE EFCONT(EPAIST,EXCEN,NSTRS,XSTRS)                               
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C     TRANSFORME LES EFFORTS EN CONTRAINTES POUR LES COQUES MINCES              
   6 : C     EXCENTREES.                                                               
   7 : C                                                                               
   8 : C  ENTREES                                                                      
   9 : C    XSTRS(NSTRS)  = EFFORTS                                                    
  10 : C    NSTRS       = NOMBRE D'EFFORTS                                             
  11 : C    EPAIST      = EPAISSEUR DE L'ELEMENT                                       
  12 : C    EXCEN       = EXCENTREMENT DE L'ELEMENT                                    
  13 : C  SORTIES                                                                      
  14 : C    XSTRS(NSTRS)  = CONTRAINTES                                                
  15 : C                                                                               
  16 : C=======================================================================        
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 :       DIMENSION XSTRS(*)                                                        
  20 :       DIMENSION WORK(8)                                                         
  21 : C                                                                               
  22 :       EE=6.D0/(EPAIST*EPAIST)                                                   
  23 :       NAUX=NSTRS/2                                                              
  24 : C                                                                               
  25 : C   CISAILLEMENT TRANSVERSE (COQ4)                                              
  26 :       IF (NSTRS.EQ.8) THEN                                                      
  27 :       WORK(7)=XSTRS(7)/EPAIST                                                   
  28 :       WORK(8)=XSTRS(8)/EPAIST                                                   
  29 :       NAUX=NAUX-1                                                               
  30 :       ENDIF                                                                     
  31 : C                                                                               
  32 : C   AUTRES CONTRAINTES                                                          
  33 :       DO 100 I=1,NAUX                                                           
  34 :          WORK(I)= XSTRS(I)/EPAIST                                               
  35 :          WORK(I+NAUX)= EE*(-EXCEN*XSTRS(I)+XSTRS(I+NAUX))                       
  36 :  100  CONTINUE                                                                  
  37 :       DO 200 I=1,NSTRS                                                          
  38 :          XSTRS(I)=WORK(I)                                                       
  39 :  200  CONTINUE                                                                  
  40 :       RETURN                                                                    
  41 :       END                                                                       
  42 :  

© Cast3M 2003 - All rights reserved.
Disclaimer