Télécharger la plaquette Cast3M


fr en

Télécharger czeril.eso

Retour à la liste

Numérotation des lignes : 

czeril

  1. C CZERIL    SOURCE    OF166741  25/11/04    21:15:41     12349          
  2.  
  3.       SUBROUTINE CZERIL(wrk52,wrk53,wrk54,wrk2,wrk3,
  4.      &  IB,IGAU,NBPGAU,necou,ecou,iecou,xecou)
  5.  
  6. *_________________________________________________________________
  7. *
  8. *
  9. *    ENTREES :
  10. *    ---------
  11. *
  12. * CELLES DE ECOINC PLUS:
  13. * DT = PAS DE TEMPS DU CALCUL
  14. *
  15. *    SORTIE :
  16. *    --------
  17. *
  18. * CELLES DE ECOINC
  19. * TRAC(2*NCOURB) LE TABLEAU DE LA COURBE DE TRACTION
  20. *
  21. * ============================================================
  22. *       ICI IL FAUT PROGRAMMER EN FORTRAN PUR
  23. *=============================================================
  24.  
  25.       IMPLICIT INTEGER(I-N)
  26.       IMPLICIT REAL*8(A-H,O-Z)
  27.  
  28. -INC PPARAM
  29. -INC CCOPTIO
  30.  
  31. -INC DECHE
  32.  
  33. -INC TECOU
  34.  
  35.       SEGMENT WRK2
  36.         REAL*8 TRAC(LTRAC)
  37.       ENDSEGMENT
  38.  
  39.       SEGMENT WRK3
  40.         REAL*8 WORK(LW),WORK2(LW2)
  41.       ENDSEGMENT
  42.  
  43.       DIMENSION RDEPS(6)
  44.  
  45.       NSTRS1 = iecou.NSTRSS
  46.  
  47. *   Traitement des deplacements suivant le type de calcul
  48. *===========================================================
  49.       IF (iecou.MFR1.EQ.1 .OR. iecou.MFR1.EQ.31)  THEN
  50. *
  51. *---> 1 formulation massive
  52. *---> 2 formulation quasi incompressible
  53. *---> MASSIF 3D
  54. *
  55.         IF (NSTRS1 .EQ. 6) THEN
  56.           DO 110 I=1,NSTRS1
  57.             RDEPS(I)=DEPST(I)
  58.  110      CONTINUE
  59.         ELSE IF (NSTRS1.EQ.4.AND.((IFOUR.EQ. 0)
  60.      &           .OR.(IFOUR.EQ.-1).OR.(IFOUR.EQ.-2))) THEN
  61. *
  62. *---> Calcul en mode deformations planes ou axisymetrique
  63. *---> Calcul en mode contraintes planes
  64. *
  65.           DO 115 I=1,NSTRS1
  66.             RDEPS(I)=DEPST(I)
  67.  115      CONTINUE
  68.           RDEPS(5)=0.D0
  69.           RDEPS(6)=0.D0
  70.         ENDIF
  71.  
  72.       ELSE
  73.         KERRE = 99
  74.         RETURN
  75.       ENDIF
  76.  
  77. *   Passage des deformations de cisaillement exprimées
  78. *   en GAMA aux déformations de cisaillement exprimées
  79. *   en déformations
  80.       DO 117 I=4,6
  81.          RDEPS(I)=RDEPS(I)*0.5D0
  82.  117  CONTINUE
  83.  
  84. *   Données du materiau
  85. *===========================================================
  86.       YOUNG=XMAT(1)
  87.       XNU=XMAT(2)
  88.       DYG0=XMAT(5)
  89.       C1=XMAT(6)
  90.       C2=XMAT(7)
  91.       C3=XMAT(8)
  92.       C4=XMAT(9)
  93.       C5=XMAT(10)
  94.       XM0=XMAT(11)
  95.       XK0=XMAT(12)
  96.       XL0=XMAT(13)
  97.       TYPE0=XMAT(14)
  98.  
  99. ***************************************************
  100. *  Preparation au calcul de la courbe de traction
  101. ***************************************************
  102. *
  103. *---> Déformation plastique équivalente EPSP
  104. *
  105.       EPSP0=var0(1)
  106. *
  107. *---> Cas contraintes planes: déformation dans l'épaisseur
  108. *
  109.       iter00=0
  110.       IF (IFOUR.EQ.-2) THEN
  111.         tr00=-1.D0*XNU/(1.D0-XNU)
  112.         RDEPS(3)=tr00*(RDEPS(1)+RDEPS(2))
  113.       ENDIF
  114.  98   treps0=RDEPS(1)+RDEPS(2)+RDEPS(3)
  115.  
  116. *---> Vitesse de deformation equivalente
  117. *      DT=5.D-7
  118.       DEPS2=RDEPS(1)*RDEPS(1)+RDEPS(2)*RDEPS(2)+RDEPS(3)*RDEPS(3)
  119.       DEPS2=DEPS2-RDEPS(1)*RDEPS(2)-RDEPS(2)*RDEPS(3)
  120.       DEPS2=DEPS2-RDEPS(3)*RDEPS(1)
  121.       DEPS3=RDEPS(4)*RDEPS(4)+RDEPS(5)*RDEPS(5)+RDEPS(6)*RDEPS(6)
  122.       DEPS2=DEPS2+3.D0*DEPS3
  123.       DEPS0=((DEPS2)**(0.5D0))*2.D0/3.D0
  124.       DEPS0=DEPS0+(treps0/3.D0)
  125.       IF (ABS(DT).LT.1D-10) THEN
  126.          DEPS=0.D0
  127.       ELSE
  128.          DEPS=DEPS0/DT
  129.       END IF
  130.  
  131. *---> Deformation maximale
  132.       EPSM0=10.D0*DEPS0
  133.  
  134. *---> Nombre de points de la courbe de traction
  135.       necou.NCOURB=100
  136.  
  137. *---> Partie constante de la courbe de traction
  138.       Y0=(XK0*(XL0**(-0.5D0)))+DYG0
  139.  
  140. *---> Partie exponentielle de la courbe de traction
  141. *
  142.       IF (ABS(DEPS).LT.1D-10) THEN
  143.         COEFF0=0.D0
  144.       ELSE
  145.         COEFF0=EXP(C4*LOG(DEPS)-C3)
  146.       END IF
  147.  
  148. *---> Debut de la courbe de traction
  149.       TRAC(1)=0.D0
  150.       TRAC(2)=0.D0
  151.       TRAC(3)=Y0
  152.       TRAC(4)=0.D0
  153.  
  154. ***************************************************
  155. *  Courbe de traction dans le cas CFC
  156. ***************************************************
  157.  
  158.       IF (TYPE0.EQ.0.) THEN
  159. *
  160. *---> Courbe de traction avant le point courant
  161. *
  162.          EPSP=0.D0
  163.          DEPS=EPSP0/48.D0
  164.          DO 100 I=3,50
  165.             EPSP=EPSP+DEPS
  166. *---> Calcul de la limite d'elasticite dans le cas cfc
  167.             TRAC(2*I-1)=C2*COEFF0*SQRT(EPSP)+Y0
  168.             TRAC(2*I)=EPSP
  169.  100     CONTINUE
  170. *
  171. *---> Courbe de traction apres le point courrant
  172. *
  173.          DEPS=EPSM0/(NCOURB-50.D0)
  174.          EPSP=EPSP0
  175.          DO 200 I=51,NCOURB
  176.             EPSP=EPSP+DEPS
  177. *---> Calcul de la limite d'elasticite dans le cas cfc
  178.             TRAC(2*I-1)=C2*COEFF0*SQRT(EPSP)+Y0
  179.             TRAC(2*I)=EPSP
  180.  200     CONTINUE
  181.  
  182. *****************************************************
  183. *  Courbe de traction dans le cas CC
  184. *****************************************************
  185. *
  186.       ELSE
  187. *
  188. *---> Courbe de traction avant le point courant
  189. *
  190.          EPSP=0.D0
  191.          DEPS=EPSP0/48.D0
  192.          DO 300 I=3,50
  193.             EPSP=EPSP+DEPS
  194. *---> Calcul de la limite d'elasticite dans le cas cc
  195.             TRAC(2*I-1)=(C1*COEFF0)+(C5*(EPSP**XM0))+Y0
  196.             TRAC(2*I)=EPSP
  197.  300     CONTINUE
  198. *
  199. *---> Courbe de traction apres le point courrant
  200. *
  201.          DEPS=EPSM0/(NCOURB-50.D0)
  202.          EPSP=EPSP0
  203.          DO 400 I=51,NCOURB
  204.             EPSP=EPSP+DEPS
  205. *---> Calcul de la limite d'elasticite dans le cas cc
  206.             TRAC(2*I-1)=(C1*COEFF0)+(C5*(EPSP**XM0))+Y0
  207.             TRAC(2*I)=EPSP
  208.  400     CONTINUE
  209.  
  210.       END IF
  211.  
  212. *********************************************
  213. *  Fin du calcul de la courbe de traction
  214. *  CALCUL DE L'ECOULEMENT PLASTIQUE
  215. *   MODELE PLASTIQUE ISOTROPE
  216. *********************************************
  217.       inpla0 = INPLAS
  218.       if (inpla0.eq.50) INPLAS = 5
  219.  
  220.       CALL CCOINC(wrk52,wrk53,wrk54,wrk2,wrk3,
  221.      &            IB,IGAU,NBPGAU,ecou,necou,iecou)
  222.  
  223.       INPLAS = inpla0
  224. *
  225. *  Cas contraintes planes
  226. *
  227.       IF (IFOUR.EQ.-2) THEN
  228.         iter00=iter00+1
  229.         RDEPS3=DEFP(3)+tr00*(RDEPS(1)-DEFP(1)+RDEPS(2)-DEFP(2))
  230.         RDEPS4=ABS(RDEPS(3))
  231.         RDEPS5=ABS(RDEPS3)
  232.         IF (RDEPS4.LE.1.D-10) THEN
  233.           IF ((iter00.LE.200).AND.(RDEPS5.GT.1.E-10)) THEN
  234.             RDEPS(3)=RDEPS3
  235.             GOTO 98
  236.           ENDIF
  237.         ELSE
  238.           ERR0=ABS(RDEPS3-RDEPS(3))/RDEPS4
  239.           IF ((iter00.LE.200).AND.(ERR0.GT.1.D-5)) THEN
  240.             RDEPS(3)=RDEPS3
  241.             GOTO 98
  242.           ENDIF
  243.         ENDIF
  244.       ENDIF
  245.  
  246.       RETURN
  247.       END
  248.  
  249.  
  250.  

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