Télécharger defela.eso

Retour à la liste

Numérotation des lignes :

  1. C DEFELA SOURCE AM 09/12/08 21:15:10 6582
  2. SUBROUTINE DEFELA(DEP,EPSILO,NSTRS,SIGF,IFOUR,KERRE)
  3. C_______________________________________________________________________
  4. C
  5. C CALCUL DE LA DEFORMATION ELASTIQUE EN FIN DE PAS
  6. C_______________________________________________________________________
  7. C C
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. DIMENSION DEP(NSTRS,NSTRS)
  11. DIMENSION EPSILO(NSTRS)
  12. DIMENSION SIGF(NSTRS)
  13. DIMENSION TEMP1(3,3),TEMP2(10,10)
  14. C
  15. C MATRICE DEP A INVERSER
  16. C
  17. NDI =10
  18. CALL ZERO(TEMP2,NDI,NSTRS)
  19. C
  20. C
  21. IF (IFOUR.EQ.-2) THEN
  22. CALL ZERO(TEMP1,3,3)
  23. DO 10 I=1,3
  24. DO 10 J=1,3
  25. TEMP1(I,J)=DEP(I,J)
  26. 10 CONTINUE
  27. C
  28. CALL INVMA1(TEMP1,3,3,ISING)
  29. IF (ISING.EQ.1) THEN
  30. WRITE(*,*) 'Matrice singuliere dans defela'
  31. KERRE = 99
  32. RETURN
  33. ENDIF
  34. C
  35. TEMP2(1,1)=TEMP1(1,1)
  36. TEMP2(1,2)=TEMP1(1,2)
  37. TEMP2(2,1)=TEMP1(2,1)
  38. TEMP2(2,2)=TEMP1(2,2)
  39. TEMP2(3,3)=TEMP1(3,3)
  40. ELSE
  41. DO 30 I=1,NSTRS
  42. DO 30 J=1,NSTRS
  43. TEMP2(I,J)=DEP(I,J)
  44. 30 CONTINUE
  45. CALL INVMA1(TEMP2,NDI,NSTRS,ISING)
  46. IF (ISING.EQ.1) THEN
  47. WRITE(*,*) 'Matrice singuliere dans defela'
  48. KERRE = 99
  49. RETURN
  50. ENDIF
  51. ENDIF
  52. C
  53. C CALCUL DU VECTEUR DEFORMATION ELASTIQUE EN FIN DE PAS
  54. C
  55. DO 50 I=1,NSTRS
  56. EPSILO(I)=0.D0
  57. DO 50 J=1,NSTRS
  58. EPSILO(I)=EPSILO(I)+TEMP2(I,J)*SIGF(J)
  59. 50 CONTINUE
  60. C
  61. END
  62.  
  63.  
  64.  

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