Télécharger la plaquette Cast3M


fr en

Télécharger fibeta.eso

Retour à la liste

Numérotation des lignes : 

fibeta

  1. C FIBETA    SOURCE    CB215821  17/11/30    21:16:12     9639           
  2. C FIBETA      SOURCE    AM        96/12/24    21:19:32     2448
  3. C======================================================================
  4. C       FIBETA  - D. COMBESCURE et P. PEGON - ELSA- 1996
  5. C======================================================================
  6. C
  7. C                    MODELE GLOBAL D'ETAGE
  8. C           (Sur des elements de poutre TIMO - Fibre
  9. C                 (Contrainte de cisaillement/Cisail.)
  10. C
  11.       SUBROUTINE FIBETA(XMAT,XCAR,SIG0,VAR0,SIGF,VARF,DEPS,
  12.      . WRK2,NCURVP,NCURVN,KERRE)
  13.       IMPLICIT INTEGER(I-N)
  14.       IMPLICIT REAL*8(A-H,O-Z)
  15. C
  16. C=======================================================================
  17. C             CETTE ROUTINE EST APPELE DANS ECO100
  18. C
  19. C
  20. C WRK0      = Segment materiaux
  21. C WRK1      = Segment contraintes
  22. C WRK2      = Segment courbe
  23. C NCURVP    = Longueur courbe sens positif
  24. C NCURVN    = Longueur courbe sens negatif
  25. C
  26. C======================================================================
  27. C XCAR      = Caracteristique de la section
  28. C DEPS     = Increment de deformation axiale
  29. C SIG0      = Contrainte initiale
  30. C VAR0      = Variables internes initiales
  31. C SIGF      = Contrainte finale
  32. C VARF      = Variables internes finales
  33. C DEFP      = Deformation plastique
  34. C
  35. C=======================================================================
  36. C
  37.       PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
  38.        DIMENSION XMAT(*),XCAR(*),SIG0(*),VAR0(*),SIGF(*),VARF(*),DEPS(*)
  39. C
  40. C=======================================================================
  41. C      VARIABLES ET SEGMENTS NECESSAIRES
  42. C=========================================================================
  43.  
  44. -INC PPARAM
  45. -INC CCOPTIO
  46. C  Segment de la courbe
  47.        SEGMENT WRK2
  48.          REAL*8 TRAC(LTRAC)
  49.        ENDSEGMENT
  50. C
  51. C
  52. C  Lecture materiau
  53. C
  54.        YOUN = XMAT(1)
  55.        XNU = XMAT(2)
  56.        XDELAP = XMAT(5)
  57.        XDELAN = (-1.D0)*XMAT(6)
  58.        XDMAXP = XMAT(7)
  59.        XDMAXN = XMAT(8)
  60.        XBETA  = XMAT(9)
  61.        XALPH  = XMAT(10)
  62.        XTETA  = XMAT(11)
  63. C
  64. C
  65.        IF (IDIM.EQ.3) THEN
  66.         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(2)
  67. C
  68.         XDDEP = DEPS(3)
  69.         XFOR0 = SIG0(3)
  70. C
  71.        ELSE
  72. C
  73.         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(1)
  74. C
  75.         XDDEP = DEPS(2)
  76.         XFOR0 = SIG0(2)
  77.        ENDIF
  78. C
  79. C Appel de la loi
  80. C
  81.         CALL ETAGEL(XDDEP,XFOR0,XFORF,VAR0,VARF,SECZ,
  82.      & XDELAP,XDELAN,XDMAXP,XDMAXN,XBETA,XALPH,XTETA,
  83.      & WRK2,NCURVP,NCURVN,KERRE)
  84. C
  85. C==================================================================
  86. C
  87. C  On remplit les tableaux avant de sortir
  88. C
  89. C==================================================================
  90.        IF (IDIM.EQ.3) THEN
  91.         SIGF(1) = XZER
  92.         SIGF(2) = XZER
  93.         SIGF(3) = XFORF
  94.         SIGF(4) = XZER
  95.         SIGF(5) = XZER
  96.         SIGF(6) = XZER
  97.        ELSE
  98.         SIGF(1) = XZER
  99.         SIGF(2) = XFORF
  100.         SIGF(3) = XZER
  101.        ENDIF
  102. C
  103. C=====================================================
  104. C            FIN DE LA ROUTINE DU MODELE
  105. C======================================================
  106.       RETURN
  107.  
  108.       END
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  

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