Numérotation des lignes :

1. C COUPLA SOURCE CB215821 16/04/21 21:16:06 8920
2. C COUPLA SOURCE INSL 24/10/96
3. SUBROUTINE COUPLA(EPS1,COURBE,IPCOUR,YOUN,SIGX,PENTE)
4. C
5. IMPLICIT INTEGER(I-N)
6. IMPLICIT REAL*8(A-H,O-Z)
7. DIMENSION COURBE(IPCOUR)
8. C
9. C **************************************
10. C * CAS ELASTOPLASTIQUE ECROUISSABLE *
11. C * COURBE INTRODUITE POINT PAR POINT *
12. C * CALCUL DE LA PENTE DE LA COURBE *
13. C **************************************
14. C
15. A1=1.D0
16. IF(EPS1.LT.0.D0) A1=-1.D0
17. EPSX=ABS(EPS1)
18. SIGX1=ABS(SIGX)
19. C
20. IF(EPSX.LE.COURBE(4).AND.SIGX1.LT.COURBE(3)) THEN
21. C --------------------------------------------
22. C * PREMIER TRONCON DE LA COURBE DE TRACTION *
23. C --------------------------------------------
24. PENTE=YOUN
25. SIGX=A1*YOUN*EPSX
26. EPPX=0.D0
27. ELSE
28. C
29. DO 1 I=6,IPCOUR,2
30. C --------------------------------------------------------
31. C * RECHERCHE DU POINT DE LA COURBE CORRESPONDANT A SET *
32. C --------------------------------------------------------
33. IF(EPSX.LE.COURBE(I)) THEN
34. A=COURBE(I-3)
35. B=COURBE(I-2)
36. C=COURBE(I-1)
37. D=COURBE(I)
38. GOTO 3
39. ENDIF
40. 1 CONTINUE
41. C
42. PENTE=0.D0
43. SIGX=0.D0
44. EPPX=EPS1-SIGX/YOUN
45. GOTO 4
46. C ------------------------------------------------------------
47. C * PENTEC = PENTE DE LA COURBE DE TRACTION (BP= A-B*PENTEC) *
48. C ------------------------------------------------------------
49. 3 CONTINUE
50. PENTE=(C-A)/(D-B)
51. SIGX=A1*(A+(PENTE*(EPSX-B)))
52. EPPX=EPS1-SIGX/YOUN
53. ENDIF
54. 4 CONTINUE
55. C
56. RETURN
57. END
58.
59.
60.
61.
62.
63.

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