Télécharger la plaquette Cast3M


fr en

Télécharger dynar.eso

Retour à la liste

Numérotation des lignes : 

dynar

  1. C DYNAR     SOURCE    OF166741  25/11/04    21:15:43     12349          
  2.       SUBROUTINE DYNAR(WRK52,WRK53,WRK54,IECOU)
  3.  
  4. C=======================================================================
  5. C
  6. C        MODELE VISCOPLASTIQUE VISCOENDOMMAGEABLE POUR LA DYNAMIQUE
  7. C                    - PROGRAMME PRINCIPAL -
  8. C                          VERSION 2.0
  9. C
  10. C=======================================================================
  11. C     CREATION :F.GATUINGT
  12. C=======================================================================
  13. C=========== Hypopthese: au temps t=0 tout est nul   ===================
  14. C=======================================================================
  15.         IMPLICIT REAL*8(A-H,O-Z)
  16.  
  17. -INC PPARAM
  18. -INC CCOPTIO
  19. -INC DECHE
  20.         DIMENSION EPSPTDT(6),EPSTE(6),EPSPOPT(6),EPSPO(6),SIGPO(6),
  21.      &     sigt(6),sigtdt(6),EPST(6),EPSTDT(6),EPSPT(6),EPSTDTE(6),
  22.      &   S(3,3),EPSIPP(3),ARRAY(3,3),SIGTDTE(6),DFIDSIG(6),EPSPOPTDT(6)
  23. C
  24.         PARAMETER (XUnDemi=0.5D0,XZERO=0.D0,UN=1.D0,DEUX=2.D0)
  25.         PARAMETER (XPETIT=1.D-12)
  26.         PARAMETER (TROIS=3.D0)
  27. C
  28. C - Recuperation des donnees materiau
  29. C
  30.         Eini=XMAT(1)
  31.         XNU=XMAT(2)
  32.         D0=XMAT(5)
  33.         Q1=XMAT(6)
  34.         Q2=XMAT(7)
  35.         Q3=XMAT(8)
  36.         SIGM0=XMAT(9)
  37.         XN=XMAT(10)
  38.         XKPLAS = XMAT(12)
  39.         XMPLAS = XMAT(11)
  40.         XK=XMAT(13)
  41.         XMDt=XMAT(14)
  42.         XNDt=XMAT(15)
  43.         At=XMAT(16)
  44.         Bt=XMAT(17)
  45.         XMDc=XMAT(18)
  46.         XNDc=XMAT(19)
  47.         Ac=XMAT(20)
  48.         Bc=XMAT(21)
  49.         ED0=XMAT(22)
  50. C - Pas de temps
  51.         DELTAT=TEMPF-TEMP0
  52.         deltat = max (1d-30,deltat)
  53. C - Contrainte au temps T
  54.         SIGT(1)=SIG0(1)
  55.         SIGT(2)=SIG0(2)
  56.         SIGT(3)=SIG0(3)
  57.         SIGT(4)=SIG0(4)
  58.         SIGT(5)=SIG0(5)
  59.         SIGT(6)=SIG0(6)
  60. C - Deformation au temps T
  61.         EPST(1)=EPST0(1)
  62.         EPST(2)=EPST0(2)
  63.         EPST(3)=EPST0(3)
  64.         EPST(4)=EPST0(4)*0.5D0
  65.         EPST(5)=EPST0(5)*0.5D0
  66.         EPST(6)=EPST0(6)*0.5D0
  67. C - Increment de deformation
  68.         DEPST(1)=DEPST(1)
  69.         DEPST(2)=DEPST(2)
  70.         DEPST(3)=DEPST(3)
  71.         DEPST(4)=DEPST(4)*0.5D0
  72.         DEPST(5)=DEPST(5)*0.5D0
  73.         DEPST(6)=DEPST(6)*0.5D0
  74. C======================================================================
  75. C - Subcycling si le pas de temps est trop grand=======================
  76. C======================================================================
  77.         IF (DELTAT.GT.1.D-3) THEN
  78.           coef=100
  79.         ELSE
  80.           coef=max (1.d0 , deltat/ 1.d-5)
  81.         ENDIF
  82.         nbiter=int(coef)
  83.         DO i=1,6
  84.           DEPST(i)=DEPST(i)/nbiter
  85.         ENDDO
  86.         DELTAT=DELTAT/nbiter
  87.         l=0
  88.         DO WHILE (l.LT.nbiter)
  89. C - Debut des boucles
  90. C=======================================================================
  91. C - Recuperation des variables historiques
  92. C=======================================================================
  93.         DT=VAR0(1)
  94.         fT=VAR0(2)+D0
  95.         EPSPT(1)=VAR0(3)
  96.         EPSPT(2)=VAR0(4)
  97.         EPSPT(3)=VAR0(5)
  98.         EPSPT(4)=VAR0(6)
  99.         EPSPT(5)=VAR0(7)
  100.         EPSPT(6)=VAR0(8)
  101.         SIGMT=VAR0(9)+SIGM0
  102.         dtt=VAR0(10)
  103.         dct=VAR0(11)
  104.         ES=VAR0(12)+ED0
  105. C======================================================================
  106. C - Deformation au temps t+dt
  107. C======================================================================
  108.         EPSTDT(1)=EPST(1)+DEPST(1)
  109.         EPSTDT(2)=EPST(2)+DEPST(2)
  110.         EPSTDT(3)=EPST(3)+DEPST(3)
  111.         EPSTDT(4)=EPST(4)+DEPST(4)
  112.         EPSTDT(5)=EPST(5)+DEPST(5)
  113.         EPSTDT(6)=EPST(6)+DEPST(6)
  114.         TRACEPSTDT=EPSTDT(1)+EPSTDT(2)+EPSTDT(3)
  115. C======================================================================
  116. C - Vitesse de Deformation
  117. C======================================================================
  118.         DO i=1,6
  119.           EPSPO(i)=DEPST(i)/DELTAT
  120.         ENDDO
  121.         TRACEPSPO=EPSPO(1)+EPSPO(2)+EPSPO(3)
  122. C=======================================================================
  123. C - Module d'elasticite tangent de la Matrice
  124. C=======================================================================
  125.         ETAN=(Eini)/XN*(SIGMt/SIGM0)**(1-XN)
  126. C=======================================================================
  127. C - Coefficients de compressibilite et cisaillement de la matrice sans les pores
  128. C=======================================================================
  129.         XKm=Eini/(TROIS*(UN-DEUX*XNU))
  130.         XGm=Eini/(DEUX*(UN+XNU))
  131. C=======================================================================
  132. C - Coeficients de compressibilite et cisaillement de la matrice avec les pores (Mori-Tanaka)
  133. C=======================================================================
  134.         XKporo=4.D0*XKm*XGm*(UN-fT)/(4.D0*XGm+3.D0*XKm*fT)
  135.         XGporo=XGm*(UN-fT)/(UN+fT*
  136.      &  (6.D0*XKm+12.D0*XGm)/(9.D0*XKm+8.D0*XGm))
  137. C ======================================================================
  138. C - Coefficients de compressibilite et cisaillement de la matrice endommagee
  139. C ======================================================================
  140.         XKendo=(UN-Dt)*XKporo
  141.         XGendo=(UN-Dt)*XGporo
  142. C ======================================================================
  143. C Calcul des invariants de la contrainte
  144. C ======================================================================
  145.         XI1t=SIGT(1)+SIGT(2)+SIGT(3)
  146.         XJ2t=0.5D0*((SIGT(1)-SIGT(2))**2
  147.      &  +(SIGT(2)-SIGT(3))**2+(SIGT(3)-SIGT(1))**2
  148.      &  +6.D0*(SIGT(4)**2+SIGT(5)**2+SIGT(6)**2))
  149. C======================================================================
  150. C - Deformation elastique au temps t+dt
  151. C======================================================================
  152. C Calcul du multiplicateur plastique au temps t
  153.         Fit=XJ2t/(SIGMT**2)+2*Q1*fT*COSH((Q2*XI1t)/(2*SIGMT))
  154.      &     -(1+(Q3*fT)**2)
  155.         IF (Fit.LT.0.) THEN
  156.           Fit=0.
  157.         ENDIF
  158.         xlambdapo=(fT/(1-fT))*(Fit/XKPLAS)**XMPLAS
  159.  
  160. C - Calcul des deformations plastique au temps t
  161. C Derivees partielles de fi
  162.         DFIDSIG(1)=0.5*(4*SIGT(1)-2*SIGT(2)-2*SIGT(3))/SIGMT**2
  163.      &  +Q1*Q2*ft*SINH(Q2*XI1t/(2*SIGMt))/SIGMt
  164.         DFIDSIG(2)=0.5*(-2*SIGT(1)+4*SIGT(2)-2*SIGT(3))/SIGMT**2
  165.      &  +Q1*Q2*ft*SINH(Q2*XI1t/(2*SIGMt))/SIGMt
  166.         DFIDSIG(3)=0.5*(-2*SIGT(1)-2*SIGT(2)+4*SIGT(3))/SIGMT**2
  167.      &  +Q1*Q2*ft*SINH(Q2*XI1t/(2*SIGMt))/SIGMt
  168.         DFIDSIG(4)=6/(SIGMt**2)*SIGT(4)
  169.         DFIDSIG(5)=6/(SIGMt**2)*SIGT(5)
  170.         DFIDSIG(6)=6/(SIGMt**2)*SIGT(6)
  171. C
  172.         DO i=1,6
  173.           EPSPOPT(i)=xlambdapo*DFIDSIG(i)
  174.         ENDDO
  175.         TRACEPSPOPT=EPSPOPT(1)+EPSPOPT(2)+EPSPOPT(3)
  176. C - Déformation viscoplastique
  177.         DO i=1,6
  178.           EPSPTDT(i)=EPSPT(i)+EPSPOPT(i)*DELTAT
  179.         ENDDO
  180.         TRACEPSPTDT=EPSPTDT(1)+EPSPTDT(2)+EPSPTDT(3)
  181. C - Deformation elastique
  182.         DO i=1,6
  183.           EPSTDTE(i)=EPSTDT(i)-EPSPTDT(i)
  184.         ENDDO
  185.         TRACEPSTDTE=EPSTDTE(1)+EPSTDTE(2)+EPSTDTE(3)
  186. C======================================================================
  187. C - Deformation equivalente elastique au sens de Mazars a t
  188.         ARRAY(1,1)=EPSTDTE(1)
  189.         ARRAY(2,2)=EPSTDTE(2)
  190.         ARRAY(3,3)=EPSTDTE(3)
  191.         ARRAY(1,2)=EPSTDTE(4)
  192.         ARRAY(2,3)=EPSTDTE(5)
  193.         ARRAY(1,3)=EPSTDTE(6)
  194.         ARRAY(3,1)=ARRAY(1,3)
  195.         ARRAY(3,2)=ARRAY(2,3)
  196.         ARRAY(2,1)=ARRAY(1,2)
  197. C
  198.         CALL JACOB3(ARRAY,3,EPSIPP,S)
  199.         EPSE=MAX( XZERO , EPSIPP(1) )**2 +
  200.      1       MAX( XZERO , EPSIPP(2) )**2 +
  201.      2       MAX( XZERO , EPSIPP(3) )**2
  202.         EPSE=SQRT(EPSE)
  203. C - Nouveau seuil de Mazars
  204.         ES=max (ES, EPSE)
  205. C ==============================================================
  206. C - Prediction élastique de la contrainte (pas d'évolution de d)
  207.         DO i=1,3
  208.           SIGPO(i)=XKendo*(TRACEPSPO-TRACEPSPOPT)+2*XGendo*
  209.      &    ((EPSPO(i)-1./3.*TRACEPSPO)-(EPSPOPT(i)-1./3.*TRACEPSPOPT))
  210.         ENDDO
  211.         DO i=4,6
  212.           SIGPO(i)=2*XGendo*((EPSPO(i))
  213.      &     -(EPSPOPT(i)))
  214.         ENDDO
  215. C
  216.         DO i=1,6
  217.           SIGTDTE(i)=SIGT(i)+SIGPO(i)*DELTAT
  218.         ENDDO
  219. C ============================================================
  220. C - Calcul des parametres avec la contrainte predite elastique
  221. C - Calcul des invariants à t
  222.         XI1tdte=SIGTDTE(1)+SIGTDTE(2)+SIGTDTE(3)
  223. C
  224.         XJ2tdte=0.5*((SIGTDTE(1)-SIGTDTE(2))**2
  225.      &  +(SIGTDTE(2)-SIGTDTE(3))**2+(SIGTDTE(3)-SIGTDTE(1))**2
  226.      &  +6*(SIGTDTE(4)**2+SIGTDTE(5)**2+SIGTDTE(6)**2))
  227. C
  228. C Derivees partielles de fi
  229.         DFIDSIG(1)=0.5*(4*SIGTDTE(1)
  230.      &  -2*SIGTDTE(2)-2*SIGTDTE(3))/SIGMT**2
  231.      &  +Q1*Q2*ft*SINH(Q2*XI1tdte/(2*SIGMt))/SIGMt
  232.         DFIDSIG(2)=0.5*(-2*SIGTDTE(1)+4*SIGTDTE(2)
  233.      &  -2*SIGTDTE(3))/SIGMT**2
  234.      &  +Q1*Q2*ft*SINH(Q2*XI1tdte/(2*SIGMt))/SIGMt
  235.         DFIDSIG(3)=0.5*(-2*SIGTDTE(1)-2*SIGTDTE(2)
  236.      &  +4*SIGTDTE(3))/SIGMT**2
  237.      &  +Q1*Q2*ft*SINH(Q2*XI1tdte/(2*SIGMt))/SIGMt
  238.         DFIDSIG(4)=6/(SIGMt**2)*SIGTDTE(4)
  239.         DFIDSIG(5)=6/(SIGMt**2)*SIGTDTE(5)
  240.         DFIDSIG(6)=6/(SIGMt**2)*SIGTDTE(6)
  241. C
  242.         TRACEDFIDSIG=DFIDSIG(1)+DFIDSIG(2)+DFIDSIG(3)
  243. C
  244.         DFIDSIGM=-2*XJ2tdte/(SIGMt**3)-Q1*Q2*ft*XI1tdte
  245.      &  *SINH(Q2*XI1tdte/(2*SIGMt))/(SIGMt**2)
  246. C
  247.         DFIDft=2*Q1*COSH(Q2*XI1tdte/(2*SIGMt))-2*ft*Q3*Q3
  248. C===================================================================
  249. C - Calcul de la surface seuil fi à l'aide de la contrainte prédite
  250.         Fi=XJ2tdte/(SIGMT**2)+2*Q1*fT*COSH((Q2*XI1tdte)/(2*SIGMT))
  251.      &     -(1+(Q3*fT)**2)
  252. C
  253. C======================================================================
  254. C========= Organigramme du calcul de la contrainte a t+Dt =============
  255. C======================================================================
  256. C
  257. C Test sur le signe de la deformation equivalente
  258.         IF(EPSE.GT.XZERO) THEN
  259.           IF(fi.LT.XZERO)THEN
  260.             IF(EPSE.LT.ES) THEN
  261. C =====================================================================
  262. C - Elasticite Seule
  263.               dtdt=dt
  264.               ftdt=ft
  265.               dttdt=dtt
  266.               dctdt=dct
  267.               sigmtdt=sigmt
  268.               DO i=1,6
  269.                 sigtdt(i)=sigtdte(i)
  270.               ENDDO
  271.             ELSE
  272. C =====================================================================
  273. C - Visco-Endommagement Seul
  274.               sigmtdt=sigmt
  275.               ftdt=ft
  276.               CALL VISCOE(XMAT,EPSIPP,EPSE,DTT,DCT,DT
  277.      &                    ,DELTAT,DTTDT,DCTDT,DTDT)
  278. C
  279.               Dpo=(Dtdt-Dt)/DELTAT
  280. C
  281. C - Raideur endommagee et ses derivees
  282.               XKendo=(1-Dtdt)*XKporo
  283.               XGendo=(1-Dtdt)*XGporo
  284.               XKpo=-Dpo*XKporo
  285.               XGpo=-Dpo*XGporo
  286. C
  287. C - Calcul de la contrainte
  288.               DO i=1,3
  289.                SIGPO(i)=XKendo*(TRACEPSPO-TRACEPSPOPT)+2*XGendo
  290.      &         *((EPSPO(i)-1./3.*TRACEPSPO)-(EPSPOPT(i)
  291.      &         -1./3.*TRACEPSPOPT))
  292.      &         +XKpo*(TRACEPSTDT-TRACEPSPTDT)+2*XGpo
  293.      &          *((EPSTDT(i)-1./3.*TRACEPSTDT)-(EPSPTDT(i)
  294.      &         -1./3.*TRACEPSPTDT))
  295.               ENDDO
  296.               DO i=4,6
  297.                 SIGPO(i)=2*XGendo*((EPSPO(i))
  298.      &          -(EPSPOPT(i)))+2*XGpo
  299.      &        *((EPSTDT(i))-(EPSPTDT(i)))
  300.               ENDDO
  301. C
  302.               DO i=1,6
  303.                 SIGTDT(i)=SIGT(i)+SIGPO(i)*DELTAT
  304.               ENDDO
  305.            ENDIF
  306.          ELSE
  307.            IF(EPSE.LT.ES) THEN
  308. C =====================================================================
  309. C - Plasticite sans endommagement
  310.              Dtdt=Dt
  311.              dttdt=dtt
  312.              dctdt=dct
  313. C - Evolution de la contrainte dans la matrice
  314.              EPSPOPM=-xlambdapo*DFIDSIGM
  315.              SIGPOM=(Etan*Eini/(Eini-Etan))*EPSPOPM
  316.              SIGMTDT=SIGMT+SIGPOM*DELTAT
  317. C - Evolution de la porosite
  318.              fpo=XK*((UN-ft)*ft*TRACEPSPOPT)
  319.              ftdt=ft+fpo*DELTAT
  320.              IF (ftdt.LT.XPETIT) THEN
  321.                 ftdt=XZERO
  322.              ENDIF
  323. C - Raideur endommagee et ses derives
  324.              XKporo=4*XKm*XGm*(1-fTDT)/(4*XGm+3*XKm*fTDT)
  325.              XGporo=XGm*(1-fTDT)/(1+fTDT*(6*XKm+12*XGm)/(9*XKm+8*XGm))
  326.              XKendo=(1-Dtdt)*XKporo
  327.              XGendo=(1-Dtdt)*XGporo
  328.              XKpo=(1-Dtdt)*(-4*XKm*XGm*(4*XGm+3*XKm)
  329.      &        /(4*XGm+3*XKm*fTDT)**2)*fpo
  330.              XGpo=(1-Dtdt)*(-5*XGm*(32*XGm**2+60*XGm*XKm+27*XKm**2)
  331.      &        /(4*(2+3*fTDT)*XGm+3*(3+2*fTDT)*XKm)**2)*fpo
  332. C =============================
  333. C - Correction de la contrainte
  334.              DO i=1,3
  335.                SIGPO(i)=XKendo*(TRACEPSPO-TRACEPSPOPT)+2*XGendo
  336.      &         *((EPSPO(i)-1./3.*TRACEPSPO)-(EPSPOPT(i)
  337.      &         -1./3.*TRACEPSPOPT))
  338.      &         +XKpo*(TRACEPSTDT-TRACEPSPTDT)+2*XGpo*
  339.      &         ((EPSTDT(i)-1./3.*TRACEPSTDT)-(EPSPTDT(i)
  340.      &          -1./3.*TRACEPSPTDT))
  341.              ENDDO
  342.              DO i=4,6
  343.                SIGPO(i)=2*XGendo*((EPSPO(i))
  344.      &          -(EPSPOPT(i)))+2*XGpo*
  345.      &         ((EPSTDT(i))-(EPSPTDT(i)))
  346.              ENDDO
  347. C
  348.              DO i=1,6
  349.                 SIGTDT(i)=SIGT(i)+SIGPO(i)*DELTAT
  350.              ENDDO
  351.            ELSE
  352. C =====================================================================
  353. C - Plasticite avec endommagement
  354. C - Evolution de l'endommagement
  355.               CALL VISCOE(XMAT,EPSIPP,EPSE,DTT,DCT,DT
  356.      &                    ,DELTAT,DTTDT,DCTDT,DTDT)
  357. C
  358.              Dpo=(Dtdt-Dt)/DELTAT
  359.              IF (Dpo.LE.0.) THEN
  360.                Dtdt=Dt
  361.                Dpo=0.
  362.              ENDIF
  363. C - Evolution de la contrainte dans la matrice
  364.              EPSPOPM=-xlambdapo*DFIDSIGM
  365.              SIGPOM=(Etan*Eini/(Eini-Etan))*EPSPOPM
  366.              SIGMTDT=SIGMT+SIGPOM*DELTAT
  367. C - Evolution de la porosite
  368.              fpo=XK*((1-ft)*ft*TRACEPSPOPT)
  369.              ftdt=ft+fpo*DELTAT
  370.              IF (ftdt.LT.0.) THEN
  371.                 ftdt=0.
  372.              ENDIF
  373. C - Raideur endommagee et ses derives
  374.              XKporo=4*XKm*XGm*(1-fTDT)/(4*XGm+3*XKm*fTDT)
  375.              XGporo=XGm*(1-fTDT)/(1+fTDT*(6*XKm+12*XGm)/(9*XKm+8*XGm))
  376.              XKendo=(1-dtdt)*XKporo
  377.              XGendo=(1-dtdt)*XGporo
  378.              XKpo=-Dpo*XKporo+(1-Dtdt)*
  379.      &        (-4*XKm*XGm*(4*XGm+3*XKm)/(4*XGm+3*XKm*fTDT)**2)*fpo
  380.              XGpo=-Dpo*XGporo
  381.      &        +(1-Dtdt)*(-5*XGm*(32*XGm**2+60*XGm*XKm+27*XKm**2)
  382.      &        /(4*(2+3*fTDT)*XGm+3*(3+2*fTDT)*XKm)**2)*fpo
  383. C =============================
  384. C - Correction de la contrainte
  385.              DO i=1,3
  386.                SIGPO(i)=XKendo*(TRACEPSPO-TRACEPSPOPT)+2*XGendo
  387.      &         *((EPSPO(i)-1./3.*TRACEPSPO)-(EPSPOPT(i)
  388.      &          -1./3.*TRACEPSPOPT))
  389.      &         +XKpo*(TRACEPSTDT-TRACEPSPTDT)+2*XGpo*
  390.      &        ((EPSTDT(i)-1./3.*TRACEPSTDT)-(EPSPTDT(i)
  391.      &         -1./3.*TRACEPSPTDT))
  392.              ENDDO
  393.              DO i=4,6
  394.                SIGPO(i)=2*XGendo*((EPSPO(i))
  395.      &           -(EPSPOPT(i)))+2*XGpo*
  396.      &        ((EPSTDT(i))-(EPSPTDT(i)))
  397.              ENDDO
  398. C
  399.              DO i=1,6
  400.                 SIGTDT(i)=SIGT(i)+SIGPO(i)*DELTAT
  401.              ENDDO
  402.            ENDIF
  403.          ENDIF
  404.         ELSE
  405.          IF(fi.LT.XZERO)THEN
  406. C =====================================================================
  407. C - Elasticite Seule
  408.            Dtdt=Dt
  409.            ftdt=ft
  410.            dttdt=dtt
  411.            dctdt=dct
  412.            sigmtdt=sigmt
  413.            DO i=1,6
  414.              sigtdt(i)=sigtdte(i)
  415.            ENDDO
  416.          ELSE
  417. C =====================================================================
  418. C - Plasticite sans endommagement
  419.            Dtdt=Dt
  420.            dttdt=dtt
  421.            dctdt=dct
  422. C - Evolution de la contrainte dans la matrice
  423.            EPSPOPM=-xlambdapo*DFIDSIGM
  424.            SIGPOM=(Etan*Eini/(Eini-Etan))*EPSPOPM
  425.            SIGMTDT=SIGMT+SIGPOM*DELTAT
  426. C - Evolution de la porosite
  427.            fpo=XK*((1-ft)*ft*TRACEPSPOPT)
  428.            ftdt=ft+fpo*DELTAT
  429.            IF (ftdt.LT.0.) THEN
  430.              ftdt=0.
  431.            ENDIF
  432. C - Raideur endommagee et ses derives
  433.            XKporo=4*XKm*XGm*(1-fTDT)/(4*XGm+3*XKm*fTDT)
  434.            XGporo=XGm*(1-fTDT)/(1+fTDT*(6*XKm+12*XGm)/(9*XKm+8*XGm))
  435.            XKendo=(1-Dtdt)*XKporo
  436.            XGendo=(1-Dtdt)*XGporo
  437.            XKpo=(1-Dtdt)*(-4*XKm*XGm*(4*XGm+3*XKm)
  438.      &      /(4*XGm+3*XKm*fTDT)**2)*fpo
  439.            XGpo=(1-Dtdt)*(-5*XGm*(32*XGm**2+60*XGm*XKm+27*XKm**2)
  440.      &      /(4*(2+3*fTDT)*XGm+3*(3+2*fTDT)*XKm)**2)*fpo
  441. C =============================
  442. C - Correction de la contrainte
  443.            DO i=1,3
  444.              SIGPO(i)=XKendo*(TRACEPSPO-TRACEPSPOPT)+2*XGendo
  445.      &       *((EPSPO(i)-1./3.*TRACEPSPO)-(EPSPOPT(i)
  446.      &       -1./3.*TRACEPSPOPT))
  447.      &       +XKpo*(TRACEPSTDT-TRACEPSPTDT)+2*XGpo*
  448.      &      ((EPSTDT(i)-1./3.*TRACEPSTDT)-(EPSPTDT(i)
  449.      &       -1./3.*TRACEPSPTDT))
  450.            ENDDO
  451.            DO i=4,6
  452.              SIGPO(i)=2*XGendo*((EPSPO(i))
  453.      &        -(EPSPOPT(i)))+2*XGpo*
  454.      &      ((EPSTDT(i))-(EPSPTDT(i)))
  455.            ENDDO
  456. C
  457.            DO i=1,6
  458.              SIGTDT(i)=SIGT(i)+SIGPO(i)*DELTAT
  459.            ENDDO
  460.          ENDIF
  461.         ENDIF
  462. C
  463. C======================================================================
  464. C========= Fin de l'organigramme du calcul de la contrainte a t+dt ====
  465. C======================================================================
  466. C
  467. C =====================================================================
  468. C - Evolution des variables historiques
  469.         VAR0(1)=DTDT
  470.         VAR0(2)=fTDT-D0
  471.         VAR0(3)=EPSPTDT(1)
  472.         VAR0(4)=EPSPTDT(2)
  473.         VAR0(5)=EPSPTDT(3)
  474.         VAR0(6)=EPSPTDT(4)
  475.         VAR0(7)=EPSPTDT(5)
  476.         VAR0(8)=EPSPTDT(6)
  477.         VAR0(9)=SIGMTDT-SIGM0
  478.         VAR0(10)=dttdt
  479.         VAR0(11)=dctdt
  480.         VAR0(12)=ES-ED0
  481. C =====================================================================
  482.         SIGT(1)=SIGTDT(1)
  483.         SIGT(2)=SIGTDT(2)
  484.         SIGT(3)=SIGTDT(3)
  485.         SIGT(4)=SIGTDT(4)
  486.         SIGT(5)=SIGTDT(5)
  487.         SIGT(6)=SIGTDT(6)
  488. C =====================================================================
  489.         EPST(1)=EPSTDT(1)
  490.         EPST(2)=EPSTDT(2)
  491.         EPST(3)=EPSTDT(3)
  492.         EPST(4)=EPSTDT(4)
  493.         EPST(5)=EPSTDT(5)
  494.         EPST(6)=EPSTDT(6)
  495. C =====================================================================
  496. C  - Compteur du Subcycling
  497.         l=l+1
  498.         ENDDO
  499. C =====================================================================
  500. C ========Fin du Subcycling=================
  501. C =====================================================================
  502. C - Evolution des variables historiques
  503.         VARF(1)=DTDT
  504.         VARF(2)=fTDT-D0
  505.         VARF(3)=EPSPTDT(1)
  506.         VARF(4)=EPSPTDT(2)
  507.         VARF(5)=EPSPTDT(3)
  508.         VARF(6)=EPSPTDT(4)
  509.         VARF(7)=EPSPTDT(5)
  510.         VARF(8)=EPSPTDT(6)
  511.         VARF(9)=SIGMTDT-SIGM0
  512.         VARF(10)=dttdt
  513.         VARF(11)=dctdt
  514.         VARF(12)=ES-ED0
  515. C======================================================================
  516. C - Mise a jour des contraintes
  517.         SIGF(1)=SIGTDT(1)
  518.         SIGF(2)=SIGTDT(2)
  519.         SIGF(3)=SIGTDT(3)
  520.         SIGF(4)=SIGTDT(4)
  521.         SIGF(5)=SIGTDT(5)
  522.         SIGF(6)=SIGTDT(6)
  523. C
  524.         END
  525.  
  526.  
  527.  

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