Sources Esope | Cast3M

bareps
C BAREPS    SOURCE    CHAT      05/01/12    21:33:53     5004
      SUBROUTINE BAREPS(COORD,DEPLA,DEFO,IDEUX)
C-----------------------------------------------------------------------
C  ROUTINE DE CALCUL DES DEFORMATIONS POUR LA BARRE
C
C  ENTREE
C          DEPLA   DEPLACEMENTS
C          COORD   COORDONNEES DES 2 NOEUDS
C          IDEUX   INDICE POUR DEFORMATIONS DU 2-EME ORDRE
C
C  SORTIE
C          DEFO    DEFORMATION AUX 2 NOEUDS
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
      DIMENSION DEFO(*),COORD(3,*),V(3),W(3),DEPLA(*)
C
C
      KERRE=0
      V(3)=0.
      DO 1 I=1,IDIM
  1   V(I)=COORD(I,2)-COORD(I,1)
      XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
      XL=SQRT(XL)
      IF(XL.EQ.0.) GO TO 999
      W(3)=0.
      DO 2 I=1,IDIM
  2   W(I)=(DEPLA(I+IDIM)-DEPLA(I))/XL
      EPS=(W(1)*V(1)+W(2)*V(2)+W(3)*V(3))/XL
      IF(IDEUX.EQ.0) GO TO 100
C
C  SECOND ORDRE
C
      EPS=EPS+(W(1)*W(1)+W(2)*W(2)+W(3)*W(3))*0.5
 100  CONTINUE
      DEFO(1)=EPS
      DEFO(2)=EPS
      RETURN
 999  KERRE=1
      RETURN
      END