Télécharger la plaquette Cast3M


fr en

Télécharger breva1.eso

Retour à la liste

Numérotation des lignes : 

breva1

  1. C BREVA1    SOURCE    AM        11/12/19    21:15:00     7227
  2.       SUBROUTINE BREVA1(SIG0,DEPST,VAR0,XMAT,XCAR,SIGF,VARF,DEFP)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    DRIVER OF THE BILINEAR PLASTIC BENDING MODEL
  6. C
  7. C  INPUT:
  8. C
  9. C    SIG0             INITIAL GENERALIZED STRESS
  10. C    DEPST            INCREMENT OF GENERALIZED STRAIN
  11. C    XMAT             MATERIAL CHARACTERISTICS
  12. C    XCAR             GEOMETRIC CHARACTERISTICS
  13. C    VAR0             INITIAL INTERNAL VARIABLE
  14. C
  15. C  OUTPUT:
  16. C
  17. C    SIGF             FINAL GENERALIZED STRESS
  18. C    VARF             FINAL INTERNAL VARIABLE
  19. C    DEFP             INCREMENT OF PLASTIC DEFORMATION
  20. C
  21. C  GENERALIZED STRESS (.=0 OR F):
  22. C
  23. C    MOMY = SIG.(5)
  24. C
  25. C  MATERIAL CHARACTERISTICS:
  26. C
  27. C    E    = XMAT(1)   ELASTIC YOUNG MODULUS
  28. C    NU   = XMAT(2)   POISSON RATIO
  29. C    ETAN = XMAT(5)   PLASTIC MODULUS
  30. C    YIEF = XMAT(6)   YIELDING FORCE
  31. C
  32. C  GEOMETRIC CHARACTERISTICS:
  33. C
  34. C    INRX = XCAR(1)   X-INERTIAL TORSION
  35. C    INRY = XCAR(2)   Y-INERTIAL BENDING        <-----
  36. C    INRZ = XCAR(3)   Z-INERTIAL BENDING
  37. C    SECT = XCAR(4)   SECTION
  38. C    SRDY = XCAR(5)   TRANSVERSE SHEAR'S REDUCED SECTION
  39. C    SRDZ = XCAR(6)   TRANSVERSE SHEAR'S REDUCED SECTION
  40. C
  41. C  INTERNAL VARIABLE (.=0 OR F):
  42. C
  43. C    EPSE = VAR.(1)   CUMULATED PLASTIC AXIAL STRAIN
  44. C    EITA = VAR.(2)   TANGENT X-AXIAL STIFFNESS
  45. C    MOEY = VAR.(3)   ELASTIC AXIAL FORCE
  46. C    EPSX = VAR.(4)   AXIAL STRAIN
  47. C
  48.       IMPLICIT INTEGER(I-N)
  49.       IMPLICIT REAL*8(A-H,O-Z)
  50.  
  51. -INC PPARAM
  52. -INC CCOPTIO
  53. C
  54.       DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),
  55.      .          SIGF(*),VARF(*) ,DEFP(*)
  56. C
  57. C-----------------------------------------------------------------------
  58. C
  59. C  ELASTIC STRESS INCREMENT FOR ALL STRESS COMPONENT EXCEPT EFFX
  60. C
  61.       CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2))
  62.       SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2)
  63.       SIGF(3)=SIG0(3) + CISA*XCAR(6)*DEPST(3)
  64.       SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4)
  65.       SIGF(5)=SIG0(5) + XMAT(1)*XCAR(2)*DEPST(5)
  66.       SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6)
  67. C
  68. C  ELASTIC AND PLASTIC AXIAL STIFFNESS
  69. C
  70.       EA     = XMAT(1)*XCAR(4)
  71.       ETA    = XMAT(5)*XCAR(4)
  72. C
  73. C
  74. C  AXIAL STRAIN INCREMENT
  75. C
  76.       DEPSX = DEPST(1)
  77. C
  78. C  CALL TO THE SCALAR BILINEAR PLASTIC MODEL COMPUTING MOMY
  79. C
  80.       FEP0 = SIG0(1) - VAR0(3)
  81.       FE0  = VAR0(3)
  82.       CALL BILINE(DEPSX,XMAT(6),EA, ETA, VARF(2),
  83.      >          FE0,FEP0,SIGF(1),DEPSE)
  84. C
  85. C  FINAL INTERNAL VARIABLE
  86. C
  87.       VARF(1)=VAR0(1)+DEPSE
  88.       VARF(3)=FE0
  89.       VARF(4)=VAR0(4)+DEPSX
  90. C
  91. C  PLASTIC DEFORMATION INCREMENT
  92. C
  93.       DEFP(1) = DEPSE
  94.       DEFP(2) = 0.D0
  95.       DEFP(3) = 0.D0
  96.       DEFP(4) = 0.D0
  97.       DEFP(5) = 0.D0
  98.       DEFP(6) = 0.D0
  99.       RETURN
  100.       END
  101.  
  102.  
  103.  
  104.  
  105.  

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