1 : C CCOTRA SOURCE PV 11/03/07 21:15:28 6885 2 : SUBROUTINE CCOTRA(WRK52,WRK2,NCOURB,WRK53) 3 : C COTRAC SOURCE AM 97/12/23 21:41:23 3015 4 : c SUBROUTINE COTRAC(WRK0,WRK2,NCOURB,KERRE) 5 : IMPLICIT INTEGER(I-N) 6 : IMPLICIT REAL*8(A-H,O-Z) 7 : -INC SMEVOLL 8 : -INC SMLREEL 9 : -INC DECHE 10 : SEGMENT/WRK2/(TRAC(LTRAC)*D) 11 : DATA NCOMAX/130/ 12 : KERRE=0 13 : MEVOLL=nint(XMAT(5)) 14 : IF(MEVOLL.EQ.0) THEN 15 : KERRE=37 16 : RETURN 17 : ENDIF 18 : SEGACT MEVOLL 19 : KEVOLL=IEVOLL(1) 20 : SEGACT KEVOLL 21 : MLREEL=IPROGX 22 : MLREE1=IPROGY 23 : * on ne desactive pas les segments pour reduire la contention sur esope en // 24 : *** SEGDES KEVOLL 25 : SEGACT MLREEL,MLREE1 26 : NBPOIX=PROG(/1) 27 : NBPOIY=MLREE1.PROG(/1) 28 : IF(NBPOIX.GT.NCOMAX+1) KERRE=31 29 : IF(NBPOIX.LT.3) KERRE=32 30 : IF(KERRE.EQ.0) GO TO 693 31 : 692 CONTINUE 32 : C 33 : C ERREUR 34 : C 35 : * on ne desactive pas les segments pour reduire la contention sur esope en // 36 : *** SEGDES MLREEL,MLREE1 37 : *** SEGDES MEVOLL 38 : GO TO 777 39 : 693 CONTINUE 40 : C 41 : C TESTS SUR LES POINTS DE LA COURBE ( A COMPLETER ) 42 : C 43 : YUNG=XMAT(1) 44 : IF(YUNG.EQ.0.D0) THEN 45 : KERRE=34 46 : GO TO 692 47 : ENDIF 48 : PSIG=MLREE1.PROG(1) 49 : PEPS=PROG(1) 50 : IF(PSIG.NE.0.D0.OR.PEPS.NE.0.D0) THEN 51 : KERRE=35 52 : GO TO 692 53 : ENDIF 54 : C 55 : C VERIF DE LA LIMITE ELASTIQUE NON NULLE 56 : C 57 : PSIG =MLREE1.PROG(2) 58 : IF(PSIG.EQ.0.D0) THEN 59 : KERRE=30 60 : GO TO 692 61 : ENDIF 62 : PEPS=PROG(2) 63 : IF(PEPS.EQ.0.D0) THEN 64 : KERRE=35 65 : GO TO 692 66 : ENDIF 67 : PENTE=PSIG/PEPS 68 : RA=ABS(PENTE-YUNG)/YUNG 69 : IF(RA.GT.1.D-2) THEN 70 : KERRE=36 71 : GO TO 692 72 : ENDIF 73 : C 74 : C VERIF DE LA PENTE 75 : C 76 : DO 100 I=3,NBPOIX 77 : I1=I-1 78 : DEPS=PROG(I)-PROG(I1) 79 : IF(DEPS.EQ.0.D0) THEN 80 : KERRE=33 81 : GO TO 692 82 : ENDIF 83 : PENTE=(MLREE1.PROG(I)-MLREE1.PROG(I1))/DEPS 84 : IF(PENTE.GE.YUNG) THEN 85 : KERRE=33 86 : GO TO 692 87 : ENDIF 88 : 100 CONTINUE 89 : NCOURB=1 90 : TRAC(1)=MLREE1.PROG(2) 91 : TRAC(2)=0.D0 92 : DO 695 I=3,NBPOIX 93 : PEPS=PROG(I) 94 : PSIG=MLREE1.PROG(I) 95 : NCOURB=NCOURB+1 96 : TRAC(2*NCOURB-1)=PSIG 97 : TRAC(2*NCOURB)=PEPS-PSIG/YUNG 98 : 695 CONTINUE 99 : C 100 : * on ne desactive pas les segments pour reduire la contention sur esope en // 101 : *** SEGDES MLREEL,MLREE1 102 : *** SEGDES MEVOLL 103 : RETURN 104 : C 105 : 777 CONTINUE 106 : RETURN 107 : END 108 : 109 : 110 : 111 : 112 : 113 : 114 : 115 : 116 :
© Cast3M 2003 - All rights reserved.
Disclaimer