1 : C GAUV1 SOURCE CHAT 05/01/13 00:16:17 5004 2 : SUBROUTINE GAUV1(DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL, 3 : 1 MFR,IFOU,IB,IGAU,EPAIST,MELE,NPINT,NBGMAT,NELMAT,SECT, 4 : 2 LHOOK,TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI, 5 : 3 SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,TRAC,NCOURB, 6 : 4 NBPGAU,DSIGT,SIGF,VARF,DEFP,KERRE) 7 : *_______________________________________________________________________ 8 : * 9 : IMPLICIT INTEGER(I-N) 10 : IMPLICIT REAL*8(A-H,O-Z) 11 : -INC CCOPTIO 12 : * 13 : DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),TRAC(*), 14 : . SIGF(*),VARF(*) ,DEFP(*) 15 : DIMENSION VALMAT(*),VALCAR(*),TXR(IDIM,*),CRIGI(*), 16 : . DDAUX(LHOOK,*),DDHOMU(LHOOK,*),XLOC(IDIM,*), 17 : . XGLOB(IDIM,*),D1HOOK(LHOOK,*),ROTHOO(LHOOK,*) 18 : DIMENSION E(4),CONT(4),DCONT(4),NUMETA(2),EPSPLA(4),EPS0(2), 19 : . DSIGT(*),SIGPA(2),SIGPB(2),SIGNA(2),SIGNB(2),SIGSUP(2),SIGINF(2), 20 : . EPSPA(2),EPSPB(2),EPSNA(2),EPSNB(2),EPSUP(2),EPINF(2), 21 : . SIGND(2),SIGNC(2),SIGPC(2),SIGPD(2),EPSND(2),EPSNC(2),EPSPC(2), 22 : . EPSPD(2),SIGBET(2,8),EPSBET(2,8) 23 : CHARACTER*8 CMATE 24 : * 25 : * INITIALISATIONS 26 : * 27 : KERRE=0 27 : CALL ZDANUL(SIGBET,16) 27 : CALL ZDANUL(EPSBET,16) 27 : CALL ZDANUL(EPSPLA, 4) 28 : IPRECI=0 29 : IF(NCOURB.NE.10.AND.NCOURB.NE.18) THEN 30 : KERRE=31 31 : RETURN 32 : ENDIF 33 : IF(NCOURB.EQ.10) IECRAS=0 34 : IF(NCOURB.EQ.18) IECRAS=1 35 : * 36 : * ON RECUPERE LES DONNEES 37 : * 38 : * D'ABORD LES CARACTERISTIQUES 39 : * 40 : TORS = XCAR(1) 41 : XINY = XCAR(2) 42 : XINZ = XCAR(3) 43 : SB = XCAR(4) 44 : SCY = XCAR(5) 45 : SCZ = XCAR(6) 46 : VX = XCAR(7) 47 : VY = XCAR(8) 48 : VZ = XCAR(9) 49 : DX = XCAR(10) 50 : DY = XCAR(11) 51 : DZ = XCAR(12) 52 : * 53 : * ENSUITE LES MATERIAUX 54 : * 55 : YOUNG=XMAT(1) 56 : XNU =XMAT(2) 57 : IF(IECRAS.EQ.1) THEN 58 : * SELON Y 59 : SIGPA(1)=TRAC(11) 60 : SIGPB(1)=TRAC(13) 61 : SIGPC(1)=TRAC(15) 62 : SIGPD(1)=TRAC(17) 63 : EPSPA(1)=TRAC(12) 64 : EPSPB(1)=TRAC(14) 65 : EPSPC(1)=TRAC(16) 66 : EPSPD(1)=TRAC(18) 67 : SIGNA(1)=TRAC( 7) 68 : SIGNB(1)=TRAC( 5) 69 : SIGNC(1)=TRAC( 3) 70 : SIGND(1)=TRAC( 1) 71 : EPSNA(1)=TRAC( 8) 72 : EPSNB(1)=TRAC( 6) 73 : EPSNC(1)=TRAC( 4) 74 : EPSND(1)=TRAC( 2) 75 : * SELON Z 76 : SIGPA(2)=TRAC(29) 77 : SIGPB(2)=TRAC(31) 78 : SIGPC(2)=TRAC(33) 79 : SIGPD(2)=TRAC(35) 80 : EPSPA(2)=TRAC(30) 81 : EPSPB(2)=TRAC(32) 82 : EPSPC(2)=TRAC(34) 83 : EPSPD(2)=TRAC(36) 84 : SIGNA(2)=TRAC(25) 85 : SIGNB(2)=TRAC(23) 86 : SIGNC(2)=TRAC(21) 87 : SIGND(2)=TRAC(19) 88 : EPSNA(2)=TRAC(26) 89 : EPSNB(2)=TRAC(24) 90 : EPSNC(2)=TRAC(22) 91 : EPSND(2)=TRAC(20) 92 : ELSE 93 : * SELON Y 94 : SIGPA(1)=TRAC( 7) 95 : SIGPB(1)=TRAC( 9) 96 : EPSPA(1)=TRAC( 8) 97 : EPSPB(1)=TRAC(10) 98 : SIGNA(1)=TRAC( 3) 99 : SIGNB(1)=TRAC( 1) 100 : EPSNA(1)=TRAC( 4) 101 : EPSNB(1)=TRAC( 2) 102 : * SELON Z 103 : SIGPA(2)=TRAC(17) 104 : SIGPB(2)=TRAC(19) 105 : EPSPA(2)=TRAC(18) 106 : EPSPB(2)=TRAC(20) 107 : SIGNA(2)=TRAC(13) 108 : SIGNB(2)=TRAC(11) 109 : EPSNA(2)=TRAC(14) 110 : EPSNB(2)=TRAC(12) 111 : ENDIF 112 : SIGTO =XMAT(6) 113 : SIGCO =XMAT(7) 114 : * 115 : * PUIS LES EFFORTS 116 : * 117 : * EFFX = SIG0(1) 118 : * EFFY = SIG0(2) 119 : * EFFY = SIG0(3) 120 : * XMOMX= SIG0(4) 121 : * XMOMY= SIG0(5) 122 : * XMOMZ= SIG0(6) 123 : * 124 : * ENFIN LES VARIABLES INTERNES 125 : * 126 : EPSEQ =VAR0(1) 127 : EPS0(1) =VAR0(2) 128 : SIGINF(1)=VAR0(3) 129 : EPINF(1) =VAR0(4) 130 : SIGSUP(1)=VAR0(5) 131 : EPSUP(1) =VAR0(6) 132 : NUMETA(1)=INT(VAR0(7)) 133 : EPS0(2) =VAR0(8) 134 : SIGINF(2)=VAR0(9) 135 : EPINF(2) =VAR0(10) 136 : SIGSUP(2)=VAR0(11) 137 : EPSUP(2) =VAR0(12) 138 : NUMETA(2)=INT(VAR0(13)) 139 : * 140 : * CALCUL DE L INCREMENT DE CONTRAINTE 141 : * 141 : CALL CALSIG(DEPST,DDAUX,NSTRS,CMATE,VALMAT,VALCAR, 142 : 1 N2EL,N2PTEL,MFR,IFOU,IB,IGAU,EPAIST,NBPGAU, 143 : 2 MELE,NPINT,NBGMAT,NELMAT,SECT,LHOOK,TXR,XLOC, 144 : 3 XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,IRTD) 145 : * 146 : IF(IRTD.NE.1) THEN 147 : KERRE=69 148 : GOTO 998 149 : ENDIF 150 : * 151 : * 152 : * 153 : * CALCUL DES CONTRAINTES INITIALES ET DES INCREMENTS 154 : * 155 : DIVX=XCAR(10)/XCAR(1) 156 : DIVY=XCAR(11)/XCAR(2) 157 : DIVZ=XCAR(12)/XCAR(3) 158 : IF(DIVX.EQ.0.) DIVX=1.D-10/SQRT(XCAR(1)*XCAR(4)) 159 : IF(DIVY.EQ.0.) DIVY=1.D-10/SQRT(XCAR(2)*XCAR(4)) 160 : IF(DIVZ.EQ.0.) DIVZ=1.D-10/SQRT(XCAR(3)*XCAR(4)) 161 : CONT(1) = SIG0(5) * DIVY 162 : CONT(2) = SIG0(6) * DIVZ 163 : CONT(3) = SIG0(4) * DIVX 164 : CONT(4) = SIG0(1) / XCAR(4) 165 : DCONT(1) = DSIGT(5) * DIVY 166 : DCONT(2) = DSIGT(6) * DIVZ 167 : DCONT(3) = DSIGT(4) * DIVX 168 : DCONT(4) = DSIGT(1) / XCAR(4) 169 : * 170 : E(1)= YOUNG 171 : E(2)= YOUNG 172 : E(3)= YOUNG/2.D0/(1.D0+XNU) 173 : E(4)= YOUNG 174 : * 174 : CALL CASTET(SIGPA,SIGPB,SIGNA,SIGNB,SIGTO,SIGCO,EPSPA,EPSPB, 175 : . EPSNA,EPSNB,SIGSUP,SIGINF,EPSUP,EPINF,EPS0,E,CONT,DCONT, 176 : . NUMETA,EPSEQ,IECRAS,SIGND,EPSND,SIGNC,EPSNC,SIGPC,EPSPC, 177 : . SIGPD,EPSPD,SIGBET,EPSBET,IPRECI,IIMPI,KERRE) 178 : * 179 : IF(KERRE.NE.0) RETURN 180 : * 181 : * ON REMET LES EFFORTS 182 : * 183 : SIGF(1) = DCONT(4) * XCAR(4) 184 : SIGF(2) = SIG0(2)+DSIGT(2) 185 : SIGF(3) = SIG0(3)+DSIGT(3) 186 : SIGF(4) = DCONT(3)/ DIVX 187 : SIGF(5) = DCONT(1)/ DIVY 188 : SIGF(6) = DCONT(2)/ DIVZ 189 : * 190 : * PUIS LES VARIABLES INTERNES 191 : * 192 : VARF(1) = EPSEQ 193 : VARF(2) = EPS0(1) 194 : VARF(3) = SIGINF(1) 195 : VARF(4) = EPINF(1) 196 : VARF(5) = SIGSUP(1) 197 : VARF(6) = EPSUP(1) 198 : VARF(7) = REAL(NUMETA(1)) 199 : VARF(8) = EPS0(2) 200 : VARF(9) = SIGINF(2) 201 : VARF(10)= EPINF(2) 202 : VARF(11)= SIGSUP(2) 203 : VARF(12)= EPSUP(2) 204 : VARF(13)= REAL(NUMETA(2)) 205 : * 206 : * ENFIN LES DEFORMATIONS PLASTIQUES 207 : * 208 : DEFP(1) = EPSPLA(4) 209 : DEFP(2) = 0.D0 210 : DEFP(3) = 0.D0 211 : DEFP(4) = EPSPLA(3) / DX 212 : DEFP(5) = EPSPLA(1) / DY 213 : DEFP(6) = EPSPLA(2) / DZ 214 : * 215 : 998 RETURN 216 : END 217 :
© Cast3M 2003 - All rights reserved.
Disclaimer