1 : C INITT1 SOURCE AM 08/12/19 21:16:46 6237 2 : SUBROUTINE INITT1(YOUNG,NYOUNG,XNU,NXNU,EN,NEN,EM,NEM,EK,NEK,ALF2, 3 : &NALF2,BET2,NBET2,ER,NER,EA,NEA,EQ,NEQ,ALPHAT,NALFAT,RHO,NRHO, 4 : &SIGY,NSIGY,XMAT,NCOMAT,T0,TINF,TSUP,ITEST) 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 20 : C ALF2(NALF2) | caracteristiques en fonction de T 21 : C BET2(NBET2) | intervenant dans les lois 22 : C ER(NER) | d'evolution 23 : C EA(NEA) | 24 : C EQ(NEQ) | 25 : C T0 temperature d'interpolation 26 : C ITEST = 1 listes de temperatures identiques pour coef. non lin. 27 : C = 0 listes de temperatures differentes 28 : C---------------------------------------------------------------------- 29 : C 30 : C-------------------------------------------------------------------- 31 : C Sortie: XMAT(NCOMAT) tableau des parametres a la temperature T0 32 : C TINF,TSUP bornes entre lesquelles est comprise T0 quand 33 : C les listes de temperatures definissant la variation des 34 : C coefficients non lineaires sont identiques 35 : C-------------------------------------------------------------------- 36 : IMPLICIT INTEGER(I-N) 37 : IMPLICIT REAL*8(A-H,O-Z) 38 : DIMENSION YOUNG(*),XNU(*),SIGY(*),RHO(*),ALPHAT(*),XMAT(*) 39 : DIMENSION EN(*),EM(*),EK(*),ALF2(*),BET2(*),ER(*),EA(*) 40 : DIMENSION EQ(*) 41 : C 42 : C write(6,*) 'Entrée dans INITT1' 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(NALF2,ALF2,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(NBET2,BET2,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(NEQ,EQ,T0,Y,YPRIM,TIN,TSU) 82 : XMAT(10)=Y 82 : CALL DERTRA(NALFAT,ALPHAT,T0,Y,YPRIM,TIN,TSU) 83 : XMAT(11)=Y 84 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 85 : TINF=TIN 86 : TSUP=TSU 87 : ENDIF 87 : CALL DERTRA(NRHO,RHO,T0,Y,YPRIM,TIN,TSU) 88 : XMAT(12)=Y 88 : CALL DERTRA(NSIGY,SIGY,T0,Y,YPRIM,TIN,TSU) 89 : XMAT(13)=Y 90 : IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN 91 : TINF=TIN 92 : TSUP=TSU 93 : ENDIF 94 : RETURN 95 : END 96 : 97 :
© Cast3M 2003 - All rights reserved.
Disclaimer