Numérotation des lignes :

find_f3d
C FIND_F3D  SOURCE    AM        12/02/28    00:28:19     7289      Subroutine find_f3d(f,esigmae,rkappait,rkappaic,x1,x2,     .      alpha,pi2,fc,ft,ag,fc0,rkappa1,ac,bc,H66,P2,lg,ntot,     .       i2,i6,fb,parahot3,idimpara3,lerror,lcp,esigmae1,U) c     Finds the value(s) of the function f(esigmae,Kappat,Kappac,dKappat,dKappac) c     f         *   function to be calculatedc     esigmae   *   vector of the trial effective stress = elastic predictorc     rkappait  *   initial hardening in tensionc     rkappaic  *   initial hardening in compressionc     x1        *   increment of hardening in tensionc     x2        *   increment of hardening in compression       IMPLICIT REAL*8 (A-B,D-H,O-Z)      implicit integer (I-K,M,N)      implicit logical (L)      implicit character*10 (C)       dimension f(i2)      dimension esigmae(i6),esigmae1(i6)****  dimension sigmaf(i6)      dimension sigmaf(6)c     updated stress vector      dimension H66(i6,i6)c     matrix of elasticity****  dimension P2(i6,i6),pi2(i6)      dimension P2(i6,i6),pi2(6)      dimension parahot3(idimpara3),U(6,6)  ******     MESSAGES D'ERREUR ( SUPPRESSION DES AUTOMATIC OBJECTS)      IF(I6.GT.6) PRINT *, ' FIND_F3D - ERREUR  I6 = ', I6, ' > 6 '*****        i3 = 3       lerror = .false.      call sigf3D_implicit(esigmae,x1,x2,sigmaf,H66,P2,pi2,ag,i6,fc,fb,c          ***************     .            parahot3,idimpara3,lerror,lcp,esigmae1,U)       if (lerror) then        write(2,*) '>Error in sigf3d'        return      endif       rkappat = rkappait+x1      rkappac = rkappaic+x2       tau_t = fhardct(rkappat,ft)      tau_c = fhardcc(rkappac,fc,fc0,rkappa1,ac,bc,lerror)      if (lerror) then        write(2,*) '>Error in compressive hardening'        return      endif       f(1) = Ranfunc(sigmaf,i3,i6,parahot3,idimpara3,lcp)      f(1) = f(1)-tau_t       f(2) = DPfunc(sigmaf,alpha)      f(2) = f(2)-(1d0-alpha)*tau_c       return      END     

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