Télécharger la plaquette Cast3M


fr en

Télécharger find_f3d.eso

Retour à la liste

Numérotation des lignes : 

find_f3d

  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.
  41.       call sigf3D_implicit(esigmae,x1,x2,sigmaf,H66,P2,pi2,ag,i6,fc,fb,
  42. c          ***************
  43.      .            parahot3,idimpara3,lerror,lcp,esigmae1,U)
  44.  
  45.       if (lerror) then
  46.         write(2,*) '>Error in sigf3d'
  47.         return
  48.       endif
  49.  
  50.       rkappat = rkappait+x1
  51.       rkappac = rkappaic+x2
  52.  
  53.       tau_t = fhardct(rkappat,ft)
  54.       tau_c = fhardcc(rkappac,fc,fc0,rkappa1,ac,bc,lerror)
  55.       if (lerror) then
  56.         write(2,*) '>Error in compressive hardening'
  57.         return
  58.       endif
  59.  
  60.       f(1) = Ranfunc(sigmaf,i3,i6,parahot3,idimpara3,lcp)
  61.       f(1) = f(1)-tau_t
  62.  
  63.       f(2) = DPfunc(sigmaf,alpha)
  64.       f(2) = f(2)-(1d0-alpha)*tau_c
  65.  
  66.       return
  67.       END
  68.  
  69.  
  70.  
  71.  
  72.  

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