Télécharger find_f3d.eso

Retour à la liste

Numérotation des lignes :

  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