Télécharger la plaquette Cast3M


fr en

Télécharger timolo.eso

Retour à la liste

Numérotation des lignes : 

timolo

  1. C TIMOLO    SOURCE    CHAT      05/01/13    03:40:40     5004
  2.       SUBROUTINE TIMOLO(V,CONT,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
  3.      .                  XINX,XINY,XINZ,IDEUX)
  4. C-----------------------------------------------------------------------
  5. C
  6. C    ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE
  7. C    D APRES TEDEL
  8. C
  9. C  ENTREE
  10. C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  11. C    XL      LONGUEUR
  12. C    SD      SECTION DROITE
  13. C    SREDY   SECTION REDUITE A L EFFORT TRANCHANT TY
  14. C    SREDZ   SECTION REDUITE A L EFFORT TRANCHANT TZ
  15. C    XINX    INERTIE IX
  16. C    XINY    INERTIE IY
  17. C    XINZ    INERTIE IZ
  18. C    YOUN    MODULE D YOUNG
  19. C    CISAIL  MODULE DE CISAILLEMENT
  20. C    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  21. C  SORTIE
  22. C    CONT    CONTRAINTES AU MILIEU 6 PAR NOEUD
  23. C
  24. C    PP 4/9/92
  25. C-----------------------------------------------------------------------
  26.       IMPLICIT INTEGER(I-N)
  27.       IMPLICIT REAL*8(A-H,O-Z)
  28.       DIMENSION CONT(*),V(*)
  29.       DATA P5/0.5D0/
  30. C
  31.       RL=1.D00/XL
  32.       RL2=RL*RL
  33.       RL3= RL2*RL
  34. C
  35. C  SECTION "REDUITE"
  36. C
  37.       IF ( SREDZ. NE.0. ) THEN
  38.            SRZ=SREDZ
  39.       ELSE
  40.            SRZ=SD
  41.       ENDIF
  42.       IF ( SREDY. NE.0. ) THEN
  43.            SRY=SREDY
  44.       ELSE
  45.            SRY=SD
  46.       ENDIF
  47. C
  48. C  DEFORMATIONS GENERALISEES
  49. C
  50.       EPS  = (V(7)-V(1))*RL
  51.       EPTY = (V(8)-V(2))*RL - (V(6)+V(12))*P5
  52.       EPTZ = (V(9)-V(3))*RL + (V(5)+V(11))*P5
  53.       XSIX = (V(10)-V(4))*RL
  54.       XSIY =-(V(5)-V(11))*RL
  55.       XSIZ =-(V(6)-V(12))*RL
  56. C
  57. C  TERMES DU SECOND ORDRE pour la membrane
  58. C
  59.       IF(IDEUX.NE.0) THEN
  60.       AUX1=0.5D0*(EPS**2+((V(8)-V(2))*RL)**2 +((V(9)-V(3))*RL)**2)
  61.       EPS = EPS + AUX1
  62.       ENDIF
  63. C
  64. C  CONTRAINTES GENERALISEES
  65. C
  66.       SIG  =YOUN*SD*    EPS
  67.       TRANY=CISAIL*SRY* EPTY
  68.       TRANZ=CISAIL*SRZ* EPTZ
  69.       TOR  =CISAIL*XINX*XSIX
  70.       XMOMY=YOUN*XINY*  XSIY
  71.       XMOMZ=YOUN*XINZ*  XSIZ
  72. C
  73. C  REMPLISSAGE DU TABLEAU CONT
  74. C
  75.       CONT(1)=SIG
  76.       CONT(2)=TRANY
  77.       CONT(3)=TRANZ
  78.       CONT(4)=TOR
  79.       CONT(5)=XMOMY
  80.       CONT(6)=XMOMZ
  81.       RETURN
  82.       END
  83.  
  84.  
  85.  

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