Télécharger la plaquette Cast3M


fr en

Télécharger concyc1.eso

Retour à la liste

Numérotation des lignes : 

concyc1

  1. C CONCYC1   SOURCE    OF166741  25/11/04    21:15:37     12349          
  2.           SUBROUTINE CONCYC1(wrk52,wrk53,wrk54,nvari,iecou)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de mettre en oeuvre le
  6. C                        modele CONCYC (M. Vassaux & B. Richard) pour representer
  7. C                le comportement 2D/3D d'un beton ordinaire sous
  8. C                sollicitations cycliques
  9. C
  10. C Traits       : - Endommagement scalaire
  11. C                - Boucles hysteretiques
  12. C                - Deformations permanentes
  13. C                - Refermeture lineaire des fissures
  14. C                - Effet unilatteral complet
  15. C                - Non localite portant sur le seuil d'endommagement (Fd)
  16. C
  17. C Auteur       : B. Richard
  18. C
  19. C Co-auteur    : M. Vassaux
  20. C
  21. C Date         : 2010 - 2011
  22. C====&===1=========2=========3=========4=========5=========6=========7==
  23. C
  24. C-----DECLARATION GENERALE----------------------------------------------
  25. C
  26.       IMPLICIT REAL*8(A-H,O-Z)
  27.       REAL*8 DEPST2(4),DEPST3(6),VAR1(42)
  28. C
  29. C-----LISTE DES INCLUDES A CHARGER--------------------------------------
  30. C
  31.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. -INC DECHE
  35. C
  36. C-----DEFINITION DE L ISOTROPIE-----------------------------------------
  37. C
  38.       CMATE = 'ISOTROPE'
  39. C
  40. C-----PARAMETERES-------------------------------------------------------
  41. C
  42.  
  43. C     Dans le cas d'un calcul local      :
  44. C                --> NSUB = NSUB0
  45. C
  46. C     Dans le cas d'un calculs non local :
  47. C                --> calcul du seuil : NSUB = 1
  48. C                --> calcul moyenne  : NSUB = NSUB0
  49.  
  50.  
  51. C     Choix de NSUB0 de maniere adaptative
  52.       XTEMP = 0.0D0
  53.       DO I=1,(2*IDIM)
  54.          XTEMP = MAX(ABS(DEPST(I)),XTEMP)
  55.       ENDDO
  56.  
  57.       IF (XTEMP.LE.1.0D-5) THEN
  58.          NSUB0 = 10
  59.       ELSEIF ((XTEMP.GT.1.0D-5).AND.(XTEMP.LE.1.0D-4)) THEN
  60.          NSUB0 = 100
  61.       ELSEIF ((XTEMP.GT.1.0D-4).AND.(XTEMP.LE.1.0D-3)) THEN
  62.          NSUB0 = 1000
  63.       ELSEIF (XTEMP.GT.1.0D-3) THEN
  64.          NSUB0 = 1000
  65.       ENDIF
  66.  
  67.       IF (ISTEP.EQ.0) THEN
  68.          NSUB = NSUB0
  69.       ELSEIF (ISTEP.EQ.1) THEN
  70.          NSUB = 1
  71.       ELSEIF (ISTEP.EQ.2) THEN
  72.          NSUB = NSUB0
  73.       ENDIF
  74.  
  75. C     PAS DE SUBSTEPPING, tant que l'on ne sait pas gérer ISTEP depuis le script...
  76.       NSUB = 1
  77.  
  78. C     RICBET3D
  79. C-----------------------------------------------------------------------
  80. C
  81. C        On forme le vecteur des increments de deformations subdivisees
  82. C      DO K=1,6
  83. C          DEPST3(K)=DEPST(K)/NSUB
  84. C      ENDDO
  85.  
  86. C        On retient les varibles internes de depart
  87. C      DO K=1,35
  88. C          VAR1(K)=VAR0(K)
  89. C      ENDDO
  90.  
  91. C        On commence le substepping
  92. C      DO J=1,NSUB
  93.  
  94.           CALL CONCYC2(XMAT,SIG0,SIGF,VAR0,VARF,DEPST,ISTEP)
  95.  
  96. C          DO I=1,35
  97. C             VAR1(I)=VARF(I)
  98. C          ENDDO
  99. C
  100. C      ENDDO
  101.  
  102.  
  103.       RETURN
  104.       END
  105.  
  106.  
  107.  

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