Télécharger flua.eso

Retour à la liste

Numérotation des lignes :

  1. C FLUA SOURCE AM 09/12/08 21:15:14 6582
  2. SUBROUTINE FLUA(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  3. $ AN,BN,CN,NDI)
  4. C
  5. C_______________________________________________________________________
  6. C
  7. C SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATION
  8. C DE FLUAGE PROPRE
  9. C D APRES LE MODELE AMORTISSEUR
  10. C_______________________________________________________________________
  11. C
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8(A-H,O-Z)
  14. DIMENSION VARF(NVARI),CMAT(NMAT),EPSFLU(6)
  15. DIMENSION AN(*),BN(NDI,*),CN(NDI,*)
  16. C
  17. CALL ZERO(AN,NSTRS,1)
  18. CALL ZERO(BN,NDI,NSTRS)
  19. CALL ZERO(CN,NDI,NSTRS)
  20. C
  21. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  22. C
  23. VIDV = CMAT(23)
  24. TDT = CMAT(NMAT)
  25. XNU = CMAT(24)
  26. C
  27. C RECUPERATION DES VARIABLES INTERNES INITIALES
  28. C
  29. C REMPLISSAGE DES VECTEURS DE DEFORMATION DE FLUAGE
  30. C
  31. EPSFLU(1) = VARF(45)
  32. EPSFLU(2) = VARF(46)
  33. EPSFLU(3) = VARF(47)
  34. EPSFLU(4) = VARF(48)
  35. EPSFLU(5) = VARF(49)
  36. EPSFLU(6) = VARF(50)
  37. C
  38. C
  39. C CONSTRUCTION DE LA MATRICE
  40. C
  41. ADR = 0.D0
  42. BDR = TDT/(2.D0*VIDV)
  43. CDR = TDT/(2.D0*VIDV)
  44. C
  45. DO 10 I=1,NSTRS
  46. AN(I) = ADR * EPSFLU(I)
  47. 10 CONTINUE
  48. C
  49. IF(IFOUR.EQ.-2) THEN
  50. C
  51. C - CONTRAINTES PLANES
  52. C
  53. BN(1,1)=1.D0
  54. BN(1,2)=-XNU
  55. BN(2,1)=BN(1,2)
  56. BN(2,2)=BN(1,1)
  57. BN(3,3)=1.D0+XNU
  58. ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  59. C
  60. C - DEFORMATION PLANE OU AXISYMETRIQUE
  61. C
  62. BN(1,1)=1.D0
  63. BN(1,2)=-XNU
  64. BN(1,3)=BN(1,2)
  65. BN(2,1)=BN(1,2)
  66. BN(2,2)=BN(1,1)
  67. BN(2,3)=BN(1,2)
  68. BN(3,1)=BN(1,2)
  69. BN(3,2)=BN(1,2)
  70. BN(3,3)=BN(1,1)
  71. BN(4,4)=1.D0+XNU
  72. ELSE
  73. C
  74. C - CALCUL TRIDIMENSIONEL
  75. C
  76. BN(1,1)=1.D0
  77. BN(1,2)=-XNU
  78. BN(1,3)=BN(1,2)
  79. BN(2,1)=BN(1,2)
  80. BN(2,2)=BN(1,1)
  81. BN(2,3)=BN(1,2)
  82. BN(3,1)=BN(1,2)
  83. BN(3,2)=BN(1,2)
  84. BN(3,3)=BN(1,1)
  85. BN(4,4)=1.D0+XNU
  86. BN(5,5)=BN(4,4)
  87. BN(6,6)=BN(4,4)
  88. ENDIF
  89. C
  90. DO 20 I=1,NSTRS
  91. DO 20 J=1,NSTRS
  92. CN(I,J) = CDR * BN(I,J)
  93. BN(I,J) = BDR * BN(I,J)
  94. 20 CONTINUE
  95. C
  96. END
  97.  
  98.  
  99.  

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