Télécharger pougr1.eso

Retour à la liste

Numérotation des lignes :

  1. C POUGR1 SOURCE AM 10/09/20 21:15:23 6756
  2. SUBROUTINE POUGR1(COORD,DEPLA,CARAC,GRAD,WORK)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE CALCUL DES GRADIENTS POUR LA POUTRE DROITE
  5. C
  6. C ENTREE
  7. C DEPLA DEPLACEMENTS
  8. C CARAC TABLEAU DE CARACTERISTIQUES GEOMETRIQUES ET
  9. C MATERIELLES DE LA POUTRE ( 11 VALEURS )
  10. C COORD COORDONNEES DES 2 NOEUDS
  11. C
  12. C SORTIE
  13. C GRAD GRADIENTS AUX 2 NOEUDS 6 PAR NOEUD
  14. C TRAVAIL
  15. C WORK(300)
  16. C-----------------------------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. DIMENSION WORK(*)
  20. DIMENSION GRAD(*),CARAC(*),COORD(3,*),VECT(3),P(3,3),
  21. . XE(2),YE(2),ZE(2),DEPLA(*)
  22. C
  23. DO 1 I=1,3
  24. 1 VECT(I)=CARAC(6+I)
  25. DO 2 I=1,2
  26. XE(I)=COORD(1,I)
  27. YE(I)=COORD(2,I)
  28. 2 ZE(I)=COORD(3,I)
  29. XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2
  30. XL=SQRT(XL)
  31. C
  32. C CALCUL DE LA MATRICE DE PASSAGE
  33. C
  34. CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)
  35. C
  36. C PASSAGE DES DEPLACEMENTS EN AXES LOCAUX
  37. C
  38. CALL POUVEC(DEPLA,WORK,P,1)
  39. C
  40. C CALCUL DES GRADIENTS AUX 2 NOEUDS
  41. C
  42. GRAD(1)=(WORK(7)-WORK(1))/XL
  43. GRAD(4)=GRAD(1)
  44. GRAD(2)=WORK(6)
  45. GRAD(5)=WORK(12)
  46. GRAD(3)=-WORK(5)
  47. GRAD(6)=-WORK(11)
  48. RETURN
  49. END
  50.  
  51.  
  52.  
  53.  

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