Télécharger la plaquette Cast3M


fr en

Télécharger pougr1.eso

Retour à la liste

Numérotation des lignes : 

pougr1

  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