1 : C SAUFL1 SOURCE AM 09/12/08 21:15:23 6582 2 : SUBROUTINE SAUFL1(VARF,NVARI,CMAT,NMAT,SIGF, 3 : $ SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,KERRE) 4 : C 5 : C ROUTINE APPELE DANS CAMECA 6 : C_______________________________________________________________________ 7 : C 8 : C SOUS PROGRAMME QUI CALCUL LA DEFORMATION 9 : C DE FLUAGE PROPRE A LA FIN DU PAS DE TEMPS 10 : C_______________________________________________________________________ 11 : C 12 : IMPLICIT INTEGER(I-N) 13 : IMPLICIT REAL*8(A-H,O-Z) 14 : DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS) 15 : DIMENSION EPSFLU(NSTRS) 16 : DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS) 17 : DIMENSION EPSFL1(10),EPSFL2(10),EPSFL3(10) 18 : DIMENSION EPSFL4(10),EPSFL5(10),EPSFL6(10) 19 : DIMENSION EPSFLA(10),EPSFLD(10) 20 : DIMENSION SIGFI(10) 21 : C 22 : C TEST DE LA DIMENSION 23 : C 24 : NDI=10 25 : IF(NSTRS.GT.NDI) THEN 26 : PRINT *, ' SAUFL1 - DIMENSION INSUFFISANTE ' 27 : KERRE = 99 28 : RETURN 29 : ENDIF 30 : 31 : C 32 : C CALCUL DES DEFORMATIONS DE FLUAGE 33 : C POUR LES 6 CHAINES DE KELVIN VOIGT 34 : C 35 : C Chaine 1 36 : C 36 : CALL SAUFLD(VARF,NVARI,CMAT,NMAT,SIGF, 37 : $ SIGEF,NSTRS,IFOUR,1,EPSFL1,AN,BN,CN) 38 : C 39 : C Chaine 2 40 : C 40 : CALL SAUFLD(VARF,NVARI,CMAT,NMAT,SIGF, 41 : $ SIGEF,NSTRS,IFOUR,2,EPSFL2,AN,BN,CN) 42 : C 43 : C Chaine 3 44 : C 44 : CALL SAUFLD(VARF,NVARI,CMAT,NMAT,SIGF, 45 : $ SIGEF,NSTRS,IFOUR,3,EPSFL3,AN,BN,CN) 46 : C 47 : C Chaine 4 48 : C 48 : CALL SAUFLD(VARF,NVARI,CMAT,NMAT ,SIGF, 49 : $ SIGEF,NSTRS,IFOUR,4,EPSFL4,AN,BN,CN) 50 : C 51 : C Chaine 5 52 : C 52 : CALL SAUFLD(VARF,NVARI,CMAT,NMAT,SIGF, 53 : $ SIGEF,NSTRS,IFOUR,5,EPSFL5,AN,BN,CN) 54 : C 55 : C Chaine 6 56 : C 56 : CALL SAUFLD(VARF,NVARI,CMAT,NMAT,SIGF, 57 : $ SIGEF,NSTRS,IFOUR,6,EPSFL6,AN,BN,CN) 58 : C 59 : C Amortisseur 60 : C 60 : CALL SAUFLA(VARF,NVARI,CMAT,NMAT,SIGF, 61 : $ SIGEF,NSTRS,IFOUR,EPSFLA,AN,BN,CN,SIGFI) 62 : C 63 : C FLUAGE DESSICCATION 64 : C 64 : CALL SAUFL2(VARF,NVARI,CMAT,NMAT,SIGF, 65 : $ SIGEF,NSTRS,IFOUR,EPSFLD,AN,BN,CN,SIGFI) 66 : C 67 : C CALCUL DEFORMATION FLUAGE 68 : C 69 : DO 10 I=1,NSTRS 70 : EPSFLU(I) = EPSFL1(I) + EPSFL2(I) + EPSFL3(I) 71 : $ + EPSFL4(I) + EPSFL5(I) + EPSFL6(I) 72 : $ + EPSFLA(I) + EPSFLD(I) 73 : 10 CONTINUE 74 : C 75 : END 76 : 77 :
© Cast3M 2003 - All rights reserved.
Disclaimer