Esope source files

   1 : C BAREPS    SOURCE    CHAT      05/01/12    21:33:53     5004           
   2 :       SUBROUTINE BAREPS(COORD,DEPLA,DEFO,IDEUX)                                 
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE CALCUL DES DEFORMATIONS POUR LA BARRE                             
   5 : C                                                                               
   6 : C  ENTREE                                                                       
   7 : C          DEPLA   DEPLACEMENTS                                                 
   8 : C          COORD   COORDONNEES DES 2 NOEUDS                                     
   9 : C          IDEUX   INDICE POUR DEFORMATIONS DU 2-EME ORDRE                      
  10 : C                                                                               
  11 : C  SORTIE                                                                       
  12 : C          DEFO    DEFORMATION AUX 2 NOEUDS                                     
  13 : C-----------------------------------------------------------------------        
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  16 : -INC CCOPTIO                                                                    
  17 :       DIMENSION DEFO(*),COORD(3,*),V(3),W(3),DEPLA(*)                           
  18 : C                                                                               
  19 : C                                                                               
  20 :       KERRE=0                                                                   
  21 :       V(3)=0.                                                                   
  22 :       DO 1 I=1,IDIM                                                             
  23 :   1   V(I)=COORD(I,2)-COORD(I,1)                                                
  24 :       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))                                        
  25 :       XL=SQRT(XL)                                                               
  26 :       IF(XL.EQ.0.) GO TO 999                                                    
  27 :       W(3)=0.                                                                   
  28 :       DO 2 I=1,IDIM                                                             
  29 :   2   W(I)=(DEPLA(I+IDIM)-DEPLA(I))/XL                                          
  30 :       EPS=(W(1)*V(1)+W(2)*V(2)+W(3)*V(3))/XL                                    
  31 :       IF(IDEUX.EQ.0) GO TO 100                                                  
  32 : C                                                                               
  33 : C  SECOND ORDRE                                                                 
  34 : C                                                                               
  35 :       EPS=EPS+(W(1)*W(1)+W(2)*W(2)+W(3)*W(3))*0.5                               
  36 :  100  CONTINUE                                                                  
  37 :       DEFO(1)=EPS                                                               
  38 :       DEFO(2)=EPS                                                               
  39 :       RETURN                                                                    
  40 :  999  KERRE=1                                                                   
  41 :       RETURN                                                                    
  42 :       END                                                                       
  43 :