Esope source files

   1 : C KRITER    SOURCE    CHAT      05/01/13    01:07:35     5004           
   2 :       SUBROUTINE KRITER(IC,ALFADV,ALFAD1,ALFAD2,VMELAS,DPELAS,
   3 :      .DPELA1,DPELA2,PORELA,PENTE1,PENTE2,PENTE3,DPLAM,S,FS,FC,KERRE)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : -INC CCOPTIO
   8 : C
   9 :       DIMENSION S(*)
  10 : C
  11 :       TRS=S(1)+S(2)+S(3)
  12 :       SEQ2=AVM(S,S)
  13 :       SEQ=SQRT(SEQ2)
  14 :       GO TO(100,200,300,400,500),IC
  15 :       WRITE(IOIMP,10) IC
  16 :       KERRE=640
  17 :       RETURN
  18 : C
  19 : C CRITERE DE LA POROSITE (1)
  20 : C
  21 :   100 FS=-TRS/3.D0
  22 :       PORELA=PORELA+PENTE3*DPLAM
  23 :       FC=FS-PORELA
  24 :       RETURN
  25 : C
  26 : C CRITERE DE DRUCKER DUCTILE (2)
  27 : C
  28 :   200 FS=ALFADV*TRS+SEQ
  29 :       FC=FS-DPELAS
  30 :       RETURN
  31 : C
  32 : C CRITERE DE VON MISES (3)
  33 : C
  34 :   300 FS=SEQ
  35 :       VMELAS=VMELAS+PENTE1*DPLAM
  36 :       FC=FS-VMELAS
  37 :       RETURN
  38 : C
  39 : C CRITERE DE DRUCKER FRAGILE (4)
  40 : C
  41 :   400 FS=ALFAD1*TRS+SEQ
  42 :       FC=FS-DPELA1
  43 :       RETURN
  44 : C
  45 : C CRITERE DE DRUCKER FRAGILE ECROUI. (5)
  46 : C
  47 :   500 FS=ALFAD2*TRS+SEQ
  48 :       DPELA2=DPELA2+PENTE2*DPLAM
  49 :       FC=FS-DPELA2
  50 : C
  51 :    10 FORMAT(1X,'ERREUR DANS KRITER IL FAUT DETERMINER (IC) IC=',I4)
  52 : C
  53 :       RETURN
  54 :       END
  55 :