1 : C CRPLAW SOURCE CHAT 05/01/12 22:31:35 5004 2 : SUBROUTINE CRPLAW(PENTE,EPSI,SIGMA,TET,TIME,DT,ICLE,NCOURB) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : C 6 : C CE SOUS PROGRAMME CALCULE LA VITESSE DE FLUAGE 7 : C 8 : C********** ENTREES 9 : C EPSI DEF INELASTIQUE EQUIVALENT 10 : C SIGMA CONTRAINTE EQUIVALENTE 11 : C TET(I) VALEUR DU CHAMP I A L'INSTANT OU ON FAIT LE CALCUL 12 : C LES CHAMPS SONT RANGES DS LE TABLEAU TET DS L'ORDRE 13 : C DS-LEQUEL ILS SONT DONNES : A CETTE EXCEPTION PRES 14 : C S'IL Y A 1 CHAMP DE TEMPERATURES C'EST LE PREMIER 15 : C PAR EX SI 2 CHAMPS TEMPERATURE ET DOSE 16 : C TET(1)=TEMPERATURE 17 : C TET(2)=DOSE 18 : C TIME TEMPS 19 : C DT PAS DE TEMPS 20 : C ICLE =0 CALCUL DE LA VITESSE 21 : C ICLE =1 CALCUL DE DELTA EPSILON INELASTIQUE 22 : C NCOURB NUMERO DU MATERIAU INELASTIQUE 23 : C 24 : C********* SORTIE PENTE 25 : C CAS ICLE=0 PENTE = VITESSE DE FLUAGE 26 : C CAS ICLE =1 PENTE=DELTA EPSILON INELASTIQUE 27 : C 28 : C EXEMPLE ON PROGRAMME ICI EPSILON POINT = A * SIGMA**N 29 : DATA BETA,DELTA,TZERO/3.34,.1,19.72/ 30 : DATA SI,GAMMA/1.44956E+5,1.3392939E+6/ 31 : C 32 : BID1 = EPSI 33 : BID2 = TET 34 : BID3=TIME 35 : BID4=DT 36 : IBID1=ICLE 37 : IBID2=NCOURB 38 : A=1.D-7 39 : XN=5.0D0 40 : PENTE=A*SIGMA**XN 41 : RETURN 42 : END 43 :
© Cast3M 2003 - All rights reserved.
Disclaimer