Télécharger elasf1.eso

Retour à la liste

Numérotation des lignes :

  1. C ELASF1 SOURCE AM 09/12/08 21:15:12 6582
  2. SUBROUTINE ELASF1(DEP,AN,BN,CN,DEPST,SIGEF,
  3. $ NSTRS,DSIGT,KERRE)
  4. C
  5. C ROUTINE APPELE DANS CAMECA
  6. C_______________________________________________________________________
  7. C
  8. C SOUS PROGRAMME QUI CALCUL LES CONTRAINTES DE TEST FINALES
  9. C ELASTIQUES CORRIGE PAR LE FLUAGE TOTAL
  10. C_______________________________________________________________________
  11. C
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8(A-H,O-Z)
  14. DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  15. DIMENSION DEP(NSTRS,NSTRS),SIGEF(NSTRS),DSIGT(NSTRS)
  16. DIMENSION DEPST(NSTRS)
  17. *
  18. DIMENSION EFLU(10,10),TEMP(10,10)
  19. DIMENSION DEFLUN(10),DEFELA(10)
  20.  
  21. C
  22. C INITIALISATION DES VARIABLES
  23. C
  24. NDI = 10
  25. CALL ZERO(DSIGT,NSTRS,1)
  26. CALL ZERO(EFLU,NDI,NSTRS)
  27. CALL ZERO(TEMP,NDI,NSTRS)
  28. CALL ZERO(DEFELA,NSTRS,1)
  29. CALL ZERO(DEFLUN,NSTRS,1)
  30. C
  31. DO 5 I=1,NSTRS
  32. DO 5 J=1,NSTRS
  33. DO 5 K=1,NSTRS
  34. TEMP(I,J) = TEMP(I,J) + DEP(I,K) * CN(K,J)
  35. 5 CONTINUE
  36. C
  37. C CONSTRUCTION DE LA MATRICE D ELASTICITE CORRIGE PAR LE FLUAGE
  38. C
  39. C EQUATION : (2.4-1)
  40. C
  41. DO 10 I=1,NSTRS
  42. TEMP(I,I) = 1.D0 + TEMP(I,I)
  43. 10 CONTINUE
  44. C
  45. CALL INVMA1(TEMP,NDI,NSTRS,ISING)
  46. IF (ISING.EQ.1) THEN
  47. WRITE(*,*) ' '
  48. WRITE(*,*) 'ERREUR DANS LE CALCUL DE LA MATRICE DE HOOKE'
  49. WRITE(*,*) 'CORRIGE PAR LE FLUAGE DANS ELASF1'
  50. WRITE(*,*) 'MATRICE NON INVERSIBLE'
  51. WRITE(*,*) ' '
  52. KERRE=99
  53. RETURN
  54. ENDIF
  55. C
  56. C CALCUL DE (1 + CN)-1 * E0
  57. C
  58. DO 20 I=1,NSTRS
  59. DO 20 J=1,NSTRS
  60. DO 20 K=1,NSTRS
  61. EFLU(I,J) = EFLU(I,J) + TEMP(I,K)*DEP(K,J)
  62. 20 CONTINUE
  63. C
  64. C DO 30 I=1,NSTRS
  65. C DO 30 J=1,NSTRS
  66. C DEP(I,J) = EFLU(I,J)
  67. C 30 CONTINUE
  68. C
  69. C CALCUL DES DEFORMATIONS DE FLUAGE PUIS ELASTIQUES
  70. C
  71. DO 40 I=1,NSTRS
  72. DEFLUN(I) = DEFLUN(I) + AN(I)
  73. DO 40 J=1,NSTRS
  74. DEFLUN(I) = DEFLUN(I) + (BN(I,J)+CN(I,J)) * SIGEF(J)
  75. 40 CONTINUE
  76. C
  77. DO 50 I=1,NSTRS
  78. DEFELA(I) = DEPST(I) - DEFLUN(I)
  79. 50 CONTINUE
  80. C
  81. C CALCUL DE L INCREMENT DE CONTRAINTES ELASTIQUES
  82. C
  83. C EQUATION (2.4-2)
  84. C
  85. DO 60 I=1,NSTRS
  86. DO 60 J=1,NSTRS
  87. DSIGT(I) = DSIGT(I) + EFLU(I,J) * DEFELA(J)
  88. 60 CONTINUE
  89. C
  90. END
  91.  
  92.  
  93.  

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