Télécharger pougf1.eso

Retour à la liste

Numérotation des lignes :

  1. C POUGF1 SOURCE AM 10/09/20 21:15:21 6756
  2. SUBROUTINE POUGF1(COORD,DEPLA,CARAC,GRAD,WORK)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE CALCUL DES GRADIENTS DE FLEXION 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 DE FLEXION 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. RL=1.D0/XL
  32. RL2=RL*RL
  33. C
  34. C CALCUL DE LA MATRICE DE PASSAGE
  35. C
  36. CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)
  37. C
  38. C PASSAGE DES DEPLACEMENTS EN AXES LOCAUX
  39. C
  40. CALL POUVEC(DEPLA,WORK,P,1)
  41. C
  42. C CALCUL DES GRADIENTS DE FLEXION AUX 2 NOEUDS
  43. C
  44. GRAD(1)=(WORK(10)-WORK(4))*RL
  45. GRAD(4)=GRAD(1)
  46. GRAD(2)=6.*(WORK(3)-WORK(9))*RL2-4.*WORK(5)*RL-2.*WORK(11)*RL
  47. GRAD(5)=6.*(WORK(9)-WORK(3))*RL2+2.*WORK(5)*RL+4.*WORK(11)*RL
  48. GRAD(3)=6.*(WORK(8)-WORK(2))*RL2-4.*WORK(6)*RL-2.*WORK(12)*RL
  49. GRAD(6)=6.*(WORK(2)-WORK(8))*RL2+2.*WORK(6)*RL+4.*WORK(12)*RL
  50. RETURN
  51. END
  52.  
  53.  
  54.  
  55.  

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