Télécharger la plaquette Cast3M


fr en

Télécharger equiv1.eso

Retour à la liste

Numérotation des lignes : 

equiv1

  1. C EQUIV1    SOURCE    CB215821  16/04/21    21:16:40     8920
  2. C EQUIV1     SOURCE    INSL       24/10/96
  3.       SUBROUTINE EQUIV1(DEFR0,EPSEQC,EQSTRC,EX,RB,EPO1,ICAL)
  4. C
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7. C
  8. C            ********************************************
  9. C            **  CALCUL  APPROCHE DE EPSEQ1 ET EQSTR1 ***
  10. C            ** COURBE DE COMPRESSION  INSA-HYPERBOLE ***
  11. C            ********************************************
  12. C
  13.       IF(ICAL.EQ.0) THEN
  14.         NBIT=200
  15.         PREC=1.D-8
  16.         EPEQ0=ABS(DEFR0)
  17.         EPEQ1=EPO1
  18.         DO 10 I=1,NBIT
  19.           EPEQ=EPEQ0+(EPEQ1-EPEQ0)/2.D0
  20.           SEQ1=EX*(EPEQ-ABS(DEFR0))
  21.           SEQ2=EX*EPEQ/(1.D0+(EX/RB-2.D0/EPO1)*EPEQ+(EPEQ/EPO1)**2)
  22.           CONT=ABS((SEQ1-SEQ2)/SEQ1)
  23.           IF(CONT.LT.PREC) GOTO 20
  24.           IF(SEQ1.GT.SEQ2) THEN
  25.             EPEQ1=EPEQ
  26.           ELSE
  27.             EPEQ0=EPEQ
  28.           ENDIF
  29.  10     CONTINUE
  30.  20     CONTINUE
  31.       ELSE
  32.         EPEQ=EPO1
  33.         SEQ2=RB
  34.       ENDIF
  35. C
  36.       IF(EPEQ.GT.EPO1) EPEQ=EPO1
  37.       IF(SEQ2.GT.RB) SEQ2=RB
  38.       EPSEQC=-1.D0*EPEQ
  39.       EQSTRC=-1.D0*SEQ2
  40.       DEFR0 = EPSEQC-EQSTRC/EX
  41. C
  42.       RETURN
  43.       END
  44.  
  45.  
  46.  
  47.  
  48.  

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