Télécharger la plaquette Cast3M


fr en

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