Télécharger concyc1.eso

Retour à la liste

Numérotation des lignes :

  1. C CONCYC1 SOURCE PV 17/12/08 21:16:57 9660
  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. -INC CCOPTIO
  32. -INC DECHE
  33. C
  34. C-----DEFINITION DE L ISOTROPIE-----------------------------------------
  35. C
  36. CMATE = 'ISOTROPE'
  37. C
  38. C-----PARAMETERES-------------------------------------------------------
  39. C
  40.  
  41. C Dans le cas d'un calcul local :
  42. C --> NSUB = NSUB0
  43. C
  44. C Dans le cas d'un calculs non local :
  45. C --> calcul du seuil : NSUB = 1
  46. C --> calcul moyenne : NSUB = NSUB0
  47.  
  48.  
  49. C Choix de NSUB0 de maniere adaptative
  50. XTEMP = 0.0D0
  51.  
  52. DO I=1,(2*IDIM)
  53. XTEMP = MAX(ABS(DEPST(I)),XTEMP)
  54. ENDDO
  55.  
  56. IF (XTEMP.LE.1.0D-5) THEN
  57. NSUB0 = 10
  58. ELSEIF ((XTEMP.GT.1.0D-5).AND.(XTEMP.LE.1.0D-4)) THEN
  59. NSUB0 = 100
  60. ELSEIF ((XTEMP.GT.1.0D-4).AND.(XTEMP.LE.1.0D-3)) THEN
  61. NSUB0 = 1000
  62. ELSEIF (XTEMP.GT.1.0D-3) THEN
  63. NSUB0 = 1000
  64. ENDIF
  65.  
  66. IF (ISTEP.EQ.0) THEN
  67. NSUB = NSUB0
  68. ELSEIF (ISTEP.EQ.1) THEN
  69. NSUB = 1
  70. ELSEIF (ISTEP.EQ.2) THEN
  71. NSUB = NSUB0
  72. ENDIF
  73.  
  74. C PAS DE SUBSTEPPING, tant que l'on ne sait pas gérer ISTEP depuis le script...
  75. NSUB = 1
  76.  
  77. C RICBET3D
  78. C-----------------------------------------------------------------------
  79. C
  80.  
  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.  
  97. C DO I=1,35
  98. C VAR1(I)=VARF(I)
  99. C ENDDO
  100. C
  101. C ENDDO
  102.  
  103.  
  104. RETURN
  105. END
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  

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