Télécharger pougf2.eso

Retour à la liste

Numérotation des lignes :

  1. C POUGF2 SOURCE AM 10/09/20 21:15:22 6756
  2. SUBROUTINE POUGF2(COORD,DEPLA,CARAC,GRAD,WORK)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE CALCUL DES GRADIENTSDE FLEXION POUR LA POUTRE DROITE EN 2D
  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 DE FLEXION AUX 2 NOEUDS 1 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,*),P(2,2),
  21. . XE(2),YE(2),DEPLA(*)
  22. C
  23. DO 2 I=1,2
  24. XE(I)=COORD(1,I)
  25. 2 YE(I)=COORD(2,I)
  26. XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2
  27. XL=SQRT(XL)
  28. RL=1.D0/XL
  29. RL2=RL*RL
  30. C
  31. C CALCUL DE LA MATRICE DE PASSAGE
  32. C
  33. CALL POUPA2(XE,YE,P,KERRE)
  34. C
  35. C PASSAGE DES DEPLACEMENTS EN AXES LOCAUX
  36. C
  37. CALL POUVE2(DEPLA,WORK,P,1)
  38. C
  39. C CALCUL DES GRADIENTS DE FLEXION AUX 2 NOEUDS
  40. C
  41.  
  42. GRAD(1)=6.*(WORK(5)-WORK(2))*RL2-4.*WORK(3)*RL-2.*WORK(6)*RL
  43. GRAD(2)=6.*(WORK(2)-WORK(5))*RL2+2.*WORK(3)*RL+4.*WORK(6)*RL
  44. RETURN
  45. END
  46.  
  47.  
  48.  
  49.  
  50.  

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