Télécharger cdcy.eso

Retour à la liste

Numérotation des lignes :

cdcy
  1. C CDCY SOURCE CHAT 05/01/12 21:51:18 5004
  2. C CDCY SOURCE INSL 24/10/96
  3. SUBROUTINE CDCY(EPEQ,SEQ,IPL,PENT,EX,TU,EMAX,EPO,ICAL)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. C
  8. C *******************************
  9. C *** COURBE DE COMPRESSION ***
  10. C *** UNIAXIALE ***
  11. C *******************************
  12. C-------------------------------------------------------------------
  13. IF(EPEQ .GE.EMAX) THEN
  14. SEQ =0.D0
  15. IPL = 4
  16. PENT=0.D0
  17. GOTO 10
  18. ENDIF
  19. C-------------------------------------------------------------------
  20. IF(ICAL.EQ.1) THEN
  21. C ******************************
  22. C ** PRE-PIC LINEAIRE **
  23. C ******************************
  24. IF(EPEQ.LT.EPO) THEN
  25. SEQ =EX*EPEQ
  26. PENT=EX
  27. IPL = 0
  28. ELSE
  29. SEQ =TU*(EMAX-EPEQ)/(EMAX-EPO)
  30. PENT=-TU/(EMAX-EPO)
  31. IPL =2
  32. ENDIF
  33. GOTO 10
  34. ENDIF
  35. C-------------------------------------------------------------------
  36. C ******************************
  37. C ** COURBE DE COMPRESSION ***
  38. C ** INSA-LBS1 ***
  39. C ******************************
  40. IF(EPEQ.LT.EPO) THEN
  41. B3=EX/TU-2.D0/EPO
  42. C3=1.D0/EPO/EPO
  43. DIN3=1.D0+B3*EPEQ+C3*EPEQ*EPEQ
  44. SEQ=EX*EPEQ/DIN3
  45. PENT=EX*(1.D0-C3*EPEQ*EPEQ)/DIN3/DIN3
  46. IPL = 0
  47. ELSE
  48. SEQ=TU*(EMAX-EPEQ)/(EMAX-EPO)
  49. PENT=-TU/(EMAX-EPO)
  50. IPL=2
  51. ENDIF
  52. C-------------------------------------------------------------------
  53. 10 CONTINUE
  54. RETURN
  55. END
  56.  
  57.  
  58.  
  59.  

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