Télécharger bareps.eso

Retour à la liste

Numérotation des lignes :

  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.  
  17. -INC PPARAM
  18. -INC CCOPTIO
  19. DIMENSION DEFO(*),COORD(3,*),V(3),W(3),DEPLA(*)
  20. C
  21. C
  22. KERRE=0
  23. V(3)=0.
  24. DO 1 I=1,IDIM
  25. 1 V(I)=COORD(I,2)-COORD(I,1)
  26. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  27. XL=SQRT(XL)
  28. IF(XL.EQ.0.) GO TO 999
  29. W(3)=0.
  30. DO 2 I=1,IDIM
  31. 2 W(I)=(DEPLA(I+IDIM)-DEPLA(I))/XL
  32. EPS=(W(1)*V(1)+W(2)*V(2)+W(3)*V(3))/XL
  33. IF(IDEUX.EQ.0) GO TO 100
  34. C
  35. C SECOND ORDRE
  36. C
  37. EPS=EPS+(W(1)*W(1)+W(2)*W(2)+W(3)*W(3))*0.5
  38. 100 CONTINUE
  39. DEFO(1)=EPS
  40. DEFO(2)=EPS
  41. RETURN
  42. 999 KERRE=1
  43. RETURN
  44. END
  45.  
  46.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales