Télécharger equiv1.eso

Retour à la liste

Numérotation des lignes :

  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