1 : C INITT SOURCE CHAT 05/01/13 00:38:23 5004 2 : SUBROUTINE INITT(YOUNG,NYOUNG,XNU,NXNU,EN,NEN,EM,NEM,EK,NEK,ALPHA, 3 : &NALPHA,BETA,NBETA,ER,NER,EA,NEA,ALPHAT,NALFAT,RHO,NRHO,SIGY,NSIGY, 4 : &FKX,NKX,NPTEUR,NC,XMAT,NCOMAT,T0,TINF,TSUP,ITEST,TRUC,NCOURB) 5 : C 6 : C---------------------------------------------------------------- 7 : C Objet: Initialiser les parametres scalaires caracteristiques du 8 : C materiau a la temperature T0 9 : C---------------------------------------------------------------- 10 : C 11 : C--------------------------------------------------------------------- 12 : C Entree: YOUNG(NYOUNG) module d'Young en fonction de T 13 : C XNU(NXNU) coefficient de Poisson en fonction de T 14 : C SIGY(NSIGY) limite elastique en fonction de T 15 : C RHO(NRHO) masse volumique en fonction de T 16 : C ALPHAT(NALFAT) coeff de dilatation thermique en fonction de T 17 : C EN(NEN) 18 : C EM(NEM) 19 : C EK(NEK) courbes des autres coefficients caracteristiques 20 : C ALPHA(NALPHA) en fonction de T intervenant 21 : C BETA(NBETA) dans les lois d'evolution 22 : C ER(NER) 23 : C EA(NEA) 24 : C FKX(NKX) fonction k(X) tabulee en fonction de la temperature 25 : C NPTEUR(NC) tableau de pointeurs sur les courbes de k(X) 26 : C T0 temperature d'interpolation 27 : C ITEST = 1 listes de temperatures identiques pour coef. non lin. 28 : C = 0 listes de temperatures differentes 29 : C---------------------------------------------------------------------- 30 : C 31 : C-------------------------------------------------------------------- 32 : C Sortie: XMAT(NCOMAT) tableau des parametres a la temperature T0 33 : C TINF,TSUP bornes entre lesquelles est comprise T0 quand 34 : C les listes de temperatures definissant la variation des 35 : C coefficients non lineaires sont identiques 36 : C-------------------------------------------------------------------- 37 : IMPLICIT INTEGER(I-N) 38 : IMPLICIT REAL*8(A-H,O-Z) 39 : DIMENSION YOUNG(*),XNU(*),SIGY(*),RHO(*),ALPHAT(*),XMAT(*) 40 : DIMENSION EN(*),EM(*),EK(*),ALPHA(*),BETA(*),ER(*),EA(*) 41 : DIMENSION FKX(*),NPTEUR(*),TRUC(*) 42 : C 43 : TINF=T0 44 : TSUP=T0 44 : CALL DERTRA(NYOUNG,YOUNG,T0,Y,YPRIM,TIN,TSU) 45 : XMAT(1)=Y 45 : CALL DERTRA(NXNU,XNU,T0,Y,YPRIM,TIN,TSU) 46 : XMAT(2)=Y 46 : CALL DERTRA(NEN,EN,T0,Y,YPRIM,TIN,TSU) 47 : XMAT(3)=Y 48 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 49 : TINF=TIN 50 : TSUP=TSU 51 : ENDIF 51 : CALL DERTRA(NEM,EM,T0,Y,YPRIM,TIN,TSU) 52 : XMAT(4)=Y 53 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 54 : TINF=TIN 55 : TSUP=TSU 56 : ENDIF 56 : CALL DERTRA(NEK,EK,T0,Y,YPRIM,TIN,TSU) 57 : XMAT(5)=Y 58 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 59 : TINF=TIN 60 : TSUP=TSU 61 : ENDIF 61 : CALL DERTRA(NALPHA,ALPHA,T0,Y,YPRIM,TIN,TSU) 62 : XMAT(6)=Y 63 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 64 : TINF=TIN 65 : TSUP=TSU 66 : ENDIF 66 : CALL DERTRA(NBETA,BETA,T0,Y,YPRIM,TIN,TSU) 67 : XMAT(7)=Y 68 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 69 : TINF=TIN 70 : TSUP=TSU 71 : ENDIF 71 : CALL DERTRA(NER,ER,T0,Y,YPRIM,TIN,TSU) 72 : XMAT(8)=Y 73 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 74 : TINF=TIN 75 : TSUP=TSU 76 : ENDIF 76 : CALL DERTRA(NEA,EA,T0,Y,YPRIM,TIN,TSU) 77 : XMAT(9)=Y 78 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 79 : TINF=TIN 80 : TSUP=TSU 81 : ENDIF 81 : CALL DERTRA(NALFAT,ALPHAT,T0,Y,YPRIM,TIN,TSU) 82 : XMAT(10)=Y 82 : CALL DERTRA(NRHO,RHO,T0,Y,YPRIM,TIN,TSU) 83 : XMAT(12)=Y 83 : CALL DERTRA(NSIGY,SIGY,T0,Y,YPRIM,TIN,TSU) 84 : XMAT(13)=Y 85 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 86 : TINF=TIN 87 : TSUP=TSU 88 : ENDIF 89 : X=0. 89 : CALL VALPAR(FKX,NPTEUR,NKX,NC,X,Y,YPRIM,T0,TIN,TSU,TRUC,NCOURB) 90 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 91 : TINF=TIN 92 : TSUP=TSU 93 : ENDIF 94 : RETURN 95 : END 96 :
© Cast3M 2003 - All rights reserved.
Disclaimer