Télécharger crplaw.eso

Retour à la liste

Numérotation des lignes :

crplaw
  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.  
  44.  

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