Esope source files

   1 : C FIND_F3D  SOURCE    AM        12/02/28    00:28:19     7289           
   2 :       Subroutine find_f3d(f,esigmae,rkappait,rkappaic,x1,x2,
   3 :      .      alpha,pi2,fc,ft,ag,fc0,rkappa1,ac,bc,H66,P2,lg,ntot,
   4 :      .       i2,i6,fb,parahot3,idimpara3,lerror,lcp,esigmae1,U)
   5 : 
   6 : c     Finds the value(s) of the function f(esigmae,Kappat,Kappac,dKappat,dKappac)
   7 : 
   8 : c     f         *   function to be calculated
   9 : c     esigmae   *   vector of the trial effective stress = elastic predictor
  10 : c     rkappait  *   initial hardening in tension
  11 : c     rkappaic  *   initial hardening in compression
  12 : c     x1        *   increment of hardening in tension
  13 : c     x2        *   increment of hardening in compression
  14 : 
  15 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
  16 :       implicit integer (I-K,M,N)
  17 :       implicit logical (L)
  18 :       implicit character*10 (C)     
  19 : 
  20 :       dimension f(i2)
  21 :       dimension esigmae(i6),esigmae1(i6)  
  22 : ****  dimension sigmaf(i6)
  23 :       dimension sigmaf(6)
  24 : c     updated stress vector
  25 :       dimension H66(i6,i6)      
  26 : c     matrix of elasticity
  27 : ****  dimension P2(i6,i6),pi2(i6)
  28 :       dimension P2(i6,i6),pi2(6)
  29 :       dimension parahot3(idimpara3),U(6,6)
  30 : 
  31 : 
  32 : *****
  33 : *     MESSAGES D'ERREUR ( SUPPRESSION DES AUTOMATIC OBJECTS)
  34 :       IF(I6.GT.6) PRINT *, ' FIND_F3D - ERREUR  I6 = ', I6, ' > 6 '
  35 : *****
  36 : 
  37 : 
  38 :       i3 = 3
  39 : 
  40 :       lerror = .false.
  40 :       call sigf3D_implicit(esigmae,x1,x2,sigmaf,H66,P2,pi2,ag,i6,fc,fb,
  41 : c          ***************
  42 :      .            parahot3,idimpara3,lerror,lcp,esigmae1,U)    
  43 : 
  44 :       if (lerror) then
  45 :         write(2,*) '>Error in sigf3d'
  46 :         return
  47 :       endif
  48 : 
  49 :       rkappat = rkappait+x1               
  50 :       rkappac = rkappaic+x2                     
  51 : 
  52 :       tau_t = fhardct(rkappat,ft)
  53 :       tau_c = fhardcc(rkappac,fc,fc0,rkappa1,ac,bc,lerror)
  54 :       if (lerror) then
  55 :         write(2,*) '>Error in compressive hardening'
  56 :         return
  57 :       endif
  58 : 
  59 :       f(1) = Ranfunc(sigmaf,i3,i6,parahot3,idimpara3,lcp)
  60 :       f(1) = f(1)-tau_t
  61 : 
  62 :       f(2) = DPfunc(sigmaf,alpha)
  63 :       f(2) = f(2)-(1d0-alpha)*tau_c
  64 :   
  65 :       return
  66 :       END
  67 : 
  68 : 
  69 : 
  70 :