Télécharger saufla.eso

Retour à la liste

Numérotation des lignes :

saufla
  1. C SAUFLA SOURCE AM 09/12/08 21:15:26 6582
  2. SUBROUTINE SAUFLA(VARF,NVARI,CMAT,NMAT,SIGF,
  3. $ SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,SIGFI)
  4. C
  5. C_______________________________________________________________________
  6. C
  7. C SOUS PROGRAMME QUI CALCUL LA DEFORMATION
  8. C DE FLUAGE PROPRE A LA FIN DU PAS DE TEMPS
  9. C_______________________________________________________________________
  10. C
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8(A-H,O-Z)
  13. DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS)
  14. DIMENSION EPSFLU(6),SIGFI(NSTRS)
  15. DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  16. C
  17. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  18. C
  19. TDT = CMAT(NMAT)
  20. C
  21. C RECUPERATION DES VARIABLES INTERNES INITIALES
  22. C
  23. C FLUAGE PROPRE
  24. C
  25. EPSFLU(1) = VARF(45)
  26. EPSFLU(2) = VARF(46)
  27. EPSFLU(3) = VARF(47)
  28. EPSFLU(4) = VARF(48)
  29. EPSFLU(5) = VARF(49)
  30. EPSFLU(6) = VARF(50)
  31. C
  32. C TEST SI LE TEMPS EST NUL
  33. C
  34. IF (TDT.EQ.0.D0) THEN
  35. GOTO 30
  36. ENDIF
  37. C
  38. C INITIALISATION DES MATRICES
  39. C
  40. CALL ZERO(AN,NSTRS,1)
  41. CALL ZERO(BN,NSTRS,NSTRS)
  42. CALL ZERO(CN,NSTRS,NSTRS)
  43. CALL ZERO(SIGFI,NSTRS,1)
  44. C
  45. DO 10 I=1,NSTRS
  46. SIGFI(I) = SIGF(I)
  47. 10 CONTINUE
  48. C
  49. C_______________________________________________________________________
  50. C
  51. C FLUAGE PROPRE AMORTISSEUR
  52. C
  53. NDI=NSTRS
  54. CALL FLUA(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  55. $ AN,BN,CN,NDI)
  56. C
  57. DO 20 I=1,NSTRS
  58. EPSFLU(I) = EPSFLU(I) + AN(I)
  59. DO 20 J=1,NSTRS
  60. EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  61. 20 CONTINUE
  62. C
  63. 30 CONTINUE
  64. C
  65. C SAUVEGARDE DES DEFORMATIONS
  66. C
  67. C FLUAGE PROPRE
  68. C
  69. VARF(45) = EPSFLU(1)
  70. VARF(46) = EPSFLU(2)
  71. VARF(47) = EPSFLU(3)
  72. VARF(48) = EPSFLU(4)
  73. VARF(49) = EPSFLU(5)
  74. VARF(50) = EPSFLU(6)
  75. C
  76. END
  77.  
  78.  

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