Télécharger camflo.eso

Retour à la liste

Numérotation des lignes :

  1. C CAMFLO SOURCE AM 13/12/16 21:15:23 7832
  2. SUBROUTINE CAMFLO(IFOUR,SIG,XINV,PI,XMAT,DFDS,DFDPI,H)
  3. *----------------------------------------------------------------
  4. * MODELE CAM-CLAY
  5. * QUANTITES NECESSAIRES A L'ECOULEMENT
  6. *----------------------------------------------------------------
  7. * ENTREE :
  8. * IFOUR CF CCOPTIO
  9. * SIG CONTRAINTES
  10. * PI VARIABLE INTERNE
  11. * XMAT PROPRIETES MATERIELLES
  12. *
  13. * SORTIE :
  14. * DFDS DF / DSIG
  15. * DFDPI DF / DPI
  16. * H LOI D'ECROUISSAGE
  17. *----------------------------------------------------------------
  18. *
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21. DIMENSION SIG(*),XINV(*),XMAT(*),DFDS(*)
  22. *
  23. E0 = XMAT(5)
  24. COHE = XMAT(7)
  25. GM = XMAT(8)
  26. XLAMD= XMAT(9)
  27. XKAPA= XMAT(10)
  28. *
  29. TETA=(1.D0+E0)/(XLAMD-XKAPA)
  30. UNGM2= 1.D0/(GM*GM)
  31. AUX=(COHE-PI)*0.5D0 -XINV(1)
  32. Y =AUX**2 + (XINV(2)/GM)**2
  33. Y = SQRT(Y)
  34. UNSY=1.D0/Y
  35. FAC =AUX/3.D0
  36. *
  37. DFDS(1)=UNSY*(UNGM2*(SIG(1)-0.5D0*(SIG(2)+SIG(3))) - FAC)
  38. DFDS(2)=UNSY*(UNGM2*(SIG(2)-0.5D0*(SIG(1)+SIG(3))) - FAC)
  39. DFDS(3)=UNSY*(UNGM2*(SIG(3)-0.5D0*(SIG(1)+SIG(2))) - FAC)
  40. DFDS(4)=UNSY*UNGM2*3.D0*SIG(4)
  41. IF(IFOUR.GE.1) THEN
  42. DFDS(5)=UNSY*UNGM2*3.D0*SIG(5)
  43. DFDS(6)=UNSY*UNGM2*3.D0*SIG(6)
  44. ENDIF
  45. *
  46. DFDPI=-0.5D0*(1.D0 + AUX*UNSY)
  47. *
  48. H=TETA*(PI+COHE)*AUX*UNSY
  49. *
  50. RETURN
  51. END
  52.  
  53.  
  54.  

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