Numérotation des lignes :

defela
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