Télécharger la plaquette Cast3M


fr en

Télécharger tifolo.eso

Retour à la liste

Numérotation des lignes : 

tifolo

  1. C TIFOLO    SOURCE    CHAT      05/01/13    03:39:23     5004
  2.       SUBROUTINE TIFOLO(XL,V,CONT,LHOOK,DDHOOK,IDEUX)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE
  6. C    DE TIMOSHENKO MULTI-FIBRE D APRES TEDEL
  7. C
  8. C  ENTREE
  9. C    XL      LONGUEUR DE LA POUTRE DROITE
  10. C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  11. C    DDHOOK(LHOOK,LHOOK) MATRICE DE HOOK
  12. C    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  13. C  SORTIE
  14. C    CONT    CONTRAINTES AU MILIEU 6 PAR NOEUD
  15. C
  16. C-----------------------------------------------------------------------
  17. C     Pierre Pegon (ISPRA) Juillet/Aout 1993
  18. C-----------------------------------------------------------------------
  19.       IMPLICIT INTEGER(I-N)
  20.       IMPLICIT REAL*8(A-H,O-Z)
  21.       DIMENSION CONT(*),V(*),DDHOOK(LHOOK,LHOOK),EPSI(6)
  22.       DATA P5/0.5D0/
  23. C
  24.       RL=1.D00/XL
  25. C
  26. C  DEFORMATIONS GENERALISEES
  27. C
  28.       EPSI(1) = (V(7)-V(1))*RL
  29.       EPSI(2) = (V(8)-V(2))*RL - (V(6)+V(12))*P5
  30.       EPSI(3) = (V(9)-V(3))*RL + (V(5)+V(11))*P5
  31.       EPSI(4) = (V(10)-V(4))*RL
  32.       EPSI(5) =-(V(5)-V(11))*RL
  33.       EPSI(6) =-(V(6)-V(12))*RL
  34. C
  35. C  TERMES DU SECOND ORDRE
  36. C
  37.       IF(IDEUX.NE.0) THEN
  38.       ENDIF
  39. C
  40. C  CONTRAINTES GENERALISEES
  41. C
  42.       DO I=1,6
  43.         CONT(I)=0.D0
  44.         DO J=1,6
  45.           CONT(I)=CONT(I)+DDHOOK(I,J)*EPSI(J)
  46.         END DO
  47.       END DO
  48. C
  49.       RETURN
  50.       END
  51.  
  52.  

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