Télécharger tufadv.eso

Retour à la liste

Numérotation des lignes :

  1. C TUFADV SOURCE CHAT 05/01/13 03:54:58 5004
  2. C TUFADV SOUCE BROC 94/04/09 21:19:18
  3. C----------------------------------------------------------------
  4. SUBROUTINE TUFADV(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,
  5. & SIGPT,EPSVPT,VARPT,XMAT,XCAR,
  6. & NSTRS,NVARI,IFOUR,INPLAS,NCOMAT )
  7. C -----------------------------------------------------------------
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. C Include contenant quelques constantes dont XPI :
  11. -INC CCREEL
  12. DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*)
  13. DIMENSION SIGPT(*),EPSVPT(*),VARPT(*),XMAT(*),XCAR(*)
  14. DIMENSION DTUFI(2,2),CTUFI(2,2)
  15. C----------------------------------------------------------------
  16. C CETTE ROUTINE CALCULE LES VALEURS FINALES DE SIG EPSV VAR A
  17. C PARTIR DES VALEURS INITIALES DE L'INCREMENT DE TEMPS ET DES
  18. C DERIVEES SIGPT,EPSVPT,VARPT
  19. C----------------------------------------------------------------
  20.  
  21. THETA = XCAR(9)/2
  22. EPAIS = XCAR(2)
  23. RAYOM = XCAR(1)-(EPAIS/2)
  24. YOUN = XMAT(1)
  25. ANG=XPI*THETA/180.0D0
  26.  
  27. C EVALUATION DE L INCREMENT D ELONGATION DE FLUAGE
  28.  
  29. CALL TUFIHO(ANG,DTUFI,CTUFI,YOUN,RAYOM,EPAIS)
  30. DF = DTUFI(1,1)
  31. DEF= DTUFI(1,2)
  32. DE = DTUFI(2,2)
  33. CM = CTUFI(1,1)
  34. CMP= CTUFI(1,2)
  35. CP = CTUFI(2,2)
  36. PHIF = CM*SIG(6) + CMP*SIG(1)
  37. ELF = CMP*SIG(6) + CP*SIG(1)
  38. RAPF=ELF/PHIF
  39.  
  40. DO 10 I=1,6
  41. SIG1(I)=SIG(I)
  42. 10 CONTINUE
  43.  
  44. C________________________________________________________________
  45. C CALCUL DES DEFORMATIONS VISQUEUSES ELONGATION & ROTATION
  46. C DUES AU FLUAGE
  47. C________________________________________________________________
  48. DO 20 I =2,6
  49. EPSV1(I) =EPSV(I)+TAU*EPSVPT(I)
  50. 20 CONTINUE
  51. EPSV1(1)=RAPF*EPSV1(6)
  52. EPSVPT(1)=RAPF*EPSVPT(6)
  53. C________________________________________________________________
  54. C MISE A JOUR DES VARIABLES INTERNES
  55. C________________________________________________________________
  56. DO 30 I=1,5
  57. VAR1(I)=VARPT(I)
  58. 30 CONTINUE
  59. VAR1(6)=VAR(6)+TAU*VARPT(6)
  60. VAR1(7)=TAU*VARPT(7)
  61. SIG1(1) = SIG(1) + TAU*(SIGPT(1) - DEF*EPSVPT(6)-DE*EPSVPT(1))
  62. SIG1(6) = SIG(6) + TAU*(SIGPT(6) - DF*EPSVPT(6)-DEF*EPSVPT(1))
  63. 99 CONTINUE
  64. RETURN
  65. END
  66.  
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  
  73.  

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