1 : C ETAGE SOURCE CHAT 05/01/12 23:45:54 5004 2 : C====================================================================== 3 : C ETAGE - D. COMBESCURE et P. PEGON - ELSA- 1996 4 : C====================================================================== 5 : C 6 : C MODELE GLOBAL D'ETAGE 7 : C (Sur des elements de poutre TIMO - Effort tranchant/Cisail.) 8 : C 9 : SUBROUTINE ETAGE(WRK0,WRK1,WRK2,NCURVP,NCURVN,KERRE) 10 : IMPLICIT INTEGER(I-N) 11 : IMPLICIT REAL*8(A-H,O-Z) 12 : C 13 : C======================================================================= 14 : C CETTE ROUTINE EST APPELE DANS ECO100 15 : C 16 : C 17 : C WRK0 = Segment materiaux 18 : C WRK1 = Segment contraintes 19 : C WRK2 = Segment courbe 20 : C NCURVP = Longueur courbe sens positif 21 : C NCURVN = Longueur courbe sens negatif 22 : C 23 : C====================================================================== 24 : C XCAR = Caracteristique de la section 25 : C DEPST = Increment de deformation axiale 26 : C SIG0 = Contrainte initiale 27 : C VAR0 = Variables internes initiales 28 : C SIGF = Contrainte finale 29 : C VARF = Variables internes finales 30 : C DEFP = Deformation plastique 31 : C 32 : C======================================================================= 33 : C 34 : REAL*8 INRX,INRY,INRZ 35 : PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16) 36 : C 37 : C======================================================================= 38 : C VARIABLES ET SEGMENTS NECESSAIRES 39 : C========================================================================= 40 : -INC CCOPTIO 41 : C Segment des materiaux 42 : SEGMENT WRK0 43 : REAL*8 XMAT(NCXMAT) 44 : ENDSEGMENT 45 : C Segment des contraintes 46 : SEGMENT WRK1 47 : REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS) 48 : REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI) 49 : REAL*8 DEFP(NSTRS),XCAR(ICARA) 50 : ENDSEGMENT 51 : C Segment de la courbe 52 : SEGMENT WRK2 53 : REAL*8 TRAC(LTRAC) 54 : ENDSEGMENT 55 : C 56 : C Lecture materiau 57 : C 58 : YOUN = XMAT(1) 59 : XNU = XMAT(2) 60 : XDELAP = XMAT(5) 61 : XDELAN = (-1.D0)*XMAT(6) 62 : XDMAXP = XMAT(7) 63 : XDMAXN = XMAT(8) 64 : XBETA = XMAT(9) 65 : XALPH = XMAT(10) 66 : XTETA = XMAT(11) 67 : C 68 : C 69 : IF (IDIM.EQ.3) THEN 70 : INRX = YOUN*XCAR(1) 71 : INRY = YOUN*XCAR(2) 72 : INRZ = YOUN*XCAR(3) 73 : SECT = YOUN*XCAR(4) 74 : SRDY = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(5) 75 : SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(6) 76 : C 77 : XDDEP = DEPST(3) 78 : XFOR0 = SIG0(3) 79 : C 80 : ELSE 81 : C 82 : SECT = YOUN*XCAR(1) 83 : INRZ = YOUN*XCAR(2) 84 : SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(3) 85 : C 86 : XDDEP = DEPST(2) 87 : XFOR0 = SIG0(2) 88 : ENDIF 89 : C 90 : C Appel de la loi 91 : C 91 : CALL ETAGEL(XDDEP,XFOR0,XFORF,VAR0,VARF,SECZ, 92 : & XDELAP,XDELAN,XDMAXP,XDMAXN,XBETA,XALPH,XTETA, 93 : & WRK2,NCURVP,NCURVN,KERRE) 94 : C 95 : C================================================================== 96 : C 97 : C On remplit les tableaux avant de sortir 98 : C 99 : C================================================================== 100 : IF (IDIM.EQ.3) THEN 101 : SIGF(1) = SIG0(1) + SECT*DEPST(1) 102 : SIGF(2) = SIG0(2) + SRDY*DEPST(2) 103 : SIGF(3) = XFORF 104 : SIGF(4) = SIG0(4) + INRX*DEPST(4) 105 : SIGF(5) = SIG0(5) + INRY*DEPST(5) 106 : SIGF(6) = SIG0(6) + INRZ*DEPST(6) 107 : ELSE 108 : SIGF(1) = SIG0(1) + SECT*DEPST(1) 109 : SIGF(2) = XFORF 110 : SIGF(3) = SIG0(3) + INRZ*DEPST(3) 111 : ENDIF 112 : C 113 : C===================================================== 114 : C FIN DE LA ROUTINE DU MODELE 115 : C====================================================== 116 : RETURN 117 : 118 : END 119 : 120 : 121 :
© Cast3M 2003 - All rights reserved.
Disclaimer