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. -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.  
  44.  

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