Télécharger saufl2.eso

Retour à la liste

Numérotation des lignes :

saufl2
  1. C SAUFL2 SOURCE AM 09/12/08 21:15:25 6582
  2. SUBROUTINE SAUFL2(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 VARIABLES INTERNES INITIALES
  18. C
  19. C FLUAGE PROPRE
  20. C
  21. EPSFLU(1) = VARF(51)
  22. EPSFLU(2) = VARF(52)
  23. EPSFLU(3) = VARF(53)
  24. EPSFLU(4) = VARF(54)
  25. EPSFLU(5) = VARF(55)
  26. EPSFLU(6) = VARF(56)
  27. C
  28. C TEST SI LE TEMPS EST NUL
  29. C
  30. TDT = CMAT(NMAT)
  31. IF (TDT.EQ.0.D0) THEN
  32. GOTO 30
  33. ENDIF
  34. C
  35. C INITIALISATION DES MATRICES
  36. C
  37. CALL ZERO(AN,NSTRS,1)
  38. CALL ZERO(BN,NSTRS,NSTRS)
  39. CALL ZERO(CN,NSTRS,NSTRS)
  40. CALL ZERO(SIGFI,NSTRS,1)
  41. C
  42. DO 10 I=1,NSTRS
  43. SIGFI(I) = SIGF(I)
  44. 10 CONTINUE
  45. C
  46. C_______________________________________________________________________
  47. C
  48. C FLUAGE PROPRE AMORTISSEUR
  49. C
  50. NDI=NSTRS
  51. CALL FLUDE1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  52. $ AN,BN,CN,NDI)
  53. C
  54. DO 20 I=1,NSTRS
  55. EPSFLU(I) = EPSFLU(I) + AN(I)
  56. DO 20 J=1,NSTRS
  57. EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  58. 20 CONTINUE
  59. C
  60. 30 CONTINUE
  61. C
  62. C SAUVEGARDE DES DEFORMATIONS
  63. C
  64. C FLUAGE DESSICCATION
  65. C
  66. VARF(45) = EPSFLU(1)
  67. VARF(46) = EPSFLU(2)
  68. VARF(47) = EPSFLU(3)
  69. VARF(48) = EPSFLU(4)
  70. VARF(49) = EPSFLU(5)
  71. VARF(50) = EPSFLU(6)
  72. C
  73. RETURN
  74. END
  75.  
  76.  

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