1 : C TAKEFO SOURCE CHAT 05/01/13 03:31:49 5004 2 : SUBROUTINE TAKEFO(SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,TRAC, 3 : . NCOURB,SIGF,VARF,DEFP,KERRE) 4 : C----------------------------------------------------------------------- 5 : C 6 : C ROUTINE DE CALCUL DES CONTRAINTES EN CISAILLEMENT 7 : C CAS DU TAKEDA TRILINEAIRE SUR L'EFFORT TRANCHANT 8 : C 9 : C ENTREE 10 : C SIG0 CONTRAINTES INITIALES 11 : C DEPST INCREMENT OF GENERALIZED STRAIN 12 : C XMAT PROPRIETES DU MATERIAU 13 : C XCAR CARACTERISTIQUE DU MATERIAU 14 : C VAR0 VARFABLES INTERNES AU DEBUT 15 : C TRAC COURBE DE CHARGEMENT 16 : C SORTIE 17 : C SIGF CONTRAINTES A LA FIN 18 : C VARF VARFABLES INTERNES A LA FIN 19 : C DEFP DEFORMATION PLASTIQUE 20 : C 21 : C materiau takeda PP 17/9/92, AVRIL 92 A.P. et P.P. 22 : C PP renewed new ECOU 24/11/95 23 : C----------------------------------------------------------------------- 24 : IMPLICIT INTEGER(I-N) 25 : IMPLICIT REAL*8(A-H,O-Z) 26 : * 27 : DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),TRAC(*), 28 : . SIGF(*),VARF(*) ,DEFP(*) 29 : C 30 : C----------------------------------------------------------------------- 31 : C CARACTERISTIQUES 32 : C ================ 33 : C 34 : C INRX = XCAR(1) X-INERTIAL TORSION 35 : C INRY = XCAR(2) Y-INERTIAL BENDING 36 : C INRZ = XCAR(3) Z-INERTIAL BENDING 37 : C SECT = XCAR(4) SECTION 38 : C SRDY = XCAR(5) TRANSVERSE SHEAR'S REDUCED SECTION 39 : C SECZ = XCAR(6) Z-TRANSVERSE SHEAR SECTION <------ 40 : C----------------------------------------------------------------------- 41 : C VARIABLES INTERNES 42 : C ================== 43 : C 44 : C IFR1 = VARF(1) 45 : C GATA = VARF(2) 46 : C UEP1 = VARF(3) 47 : C PEP1 = VARF(4) 48 : C UEN1 = VARF(5) 49 : C PEN1 = VARF(6) 50 : C UCP1 = VARF(7) 51 : C UCN1 = VARF(8) 52 : C E1P1 = VARF(9) 53 : C E1N1 = VARF(10) 54 : C E2P1 = VARF(11) 55 : C E2N1 = VARF(12) 56 : C USP1 = VARF(13) 57 : C USN1 = VARF(14) 58 : C UTP1 = VARF(15) 59 : C PTP1 = VARF(16) 60 : C UTN1 = VARF(17) 61 : C PTN1 = VARF(18) 62 : C EP11 = VARF(19) 63 : C EN11 = VARF(20) 64 : C DES0 = VARF(21) 65 : C----------------------------------------------------------------------- 66 : C ON RECUPERE LES PARAMETRES DU MODELE 67 : C ==================================== 68 : C 69 : C E0 = XMAT(1) 70 : C NU = XMAT(2) 71 : C SFDP = XMAT(6) 72 : C SFDN = XMAT(7) 73 : C PINP = XMAT(8) 74 : C PINN = XMAT(9) 75 : C SRDP = XMAT(10) 76 : C SRDN = XMAT(11) 77 : C 78 : IF(NCOURB.EQ.4)THEN 79 : C 80 : C LA COURBE EST SYMETRIQUE 81 : C 82 : GAZ=TRAC(3)/TRAC(4) 83 : CMOP=TRAC(3) 84 : CMON=-CMOP 85 : YCUP=TRAC(6) 86 : YCUN=-YCUP 87 : YMOP=TRAC(5) 88 : YMON=-YMOP 89 : GAYP=(TRAC(7)-TRAC(5))/(TRAC(8)-TRAC(6)) 90 : GAYN=GAYP 91 : C 92 : XMAT(7)=XMAT(6) 93 : XMAT(9)=-XMAT(8) 94 : XMAT(11)=XMAT(10) 95 : C 96 : ELSEIF(NCOURB.EQ.7)THEN 97 : C 98 : C LA COURBE EST NON SYMETRIQUE 99 : C 100 : GAZ=(TRAC(9)-TRAC(5)) / (TRAC(10)-TRAC(6)) 101 : CMOP=TRAC(9) 102 : CMON=TRAC(5) 103 : YCUP=TRAC(12) 104 : YCUN=TRAC(4) 105 : YMOP=TRAC(11) 106 : YMON=TRAC(3) 107 : GAYP=(TRAC(13)-TRAC(11))/(TRAC(14)-TRAC(12)) 108 : GAYN=(TRAC(3) -TRAC(1) )/(TRAC(4) -TRAC(2)) 109 : C 110 : ELSE 111 : C 112 : C LA COURBE EST ERRONEE 113 : C 114 : KERRE=31 115 : RETURN 116 : C 117 : ENDIF 118 : C 119 : C----------------------------------------------------------------------- 120 : C 121 : C TEST ON ELASTIC STIFFNESS 122 : C 123 : XMODG=0.5D0 / (1.D0 + XMAT(2)) 124 : GAZZ0 = XMAT(1)*XCAR(6)*XMODG 125 : IF(ABS(GAZ-GAZZ0)/GAZ.GT.1.D-6)THEN 126 : KERRE=36 127 : RETURN 128 : ENDIF 129 : C 130 : XEFFZ = SIG0(3) 131 : C 132 : C APPELE AU PROGRAME QUI CALCULE L'EFFORT TRANCHANT ( XEFFZ ) 133 : C CORRESPONDENT A L' INCREMENT LINEAIRE DE MOMENT - DP1 - 134 : C >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 135 : DO 1 I=1,21 136 : VARF(I)=VAR0(I) 137 : 1 CONTINUE 138 : C 139 : C CONVERSION ENTIER 140 : C 141 : IFR1=INT(VARF(1)) 142 : C 143 : C INITIALISATION DE IFR1 ET DE VARF(.,2) 144 : C 145 : IFC1=0 146 : IF(IFR1.EQ.0) THEN 147 : IFR1 = 1 148 : VARF(2)=GAZZ0 149 : ENDIF 150 : C 151 : C DEFORMATIONS GENERALISEES 152 : C ========================= 153 : C 154 : C INCREMENT DE MOMENT AUTOUR DE YY 155 : C 156 : C DP1 = D_SIGMA / LINEAR_STIFFNESS * CURRENT STIFFNESS 157 : C 158 : C DP1 = DSIGT(3) / GAZZ0 * VARF(2) 159 : DP1 = DEPST(3) * VARF(2) 160 : C 161 : C 162 : C=================================================================== 162 : CALL DDNSTH(IFR1,IFC1,DP1,XEFFZ,VARF(21), 163 : * VARF(2),GAZZ0,GAYP,GAYN,CMOP,CMON,YMOP, 164 : * YMON,YCUP,YCUN,XMAT(6),XMAT(7),XMAT(8),XMAT(9), 165 : * XMAT(10),XMAT(11),VARF(3),VARF(4),VARF(5), 166 : * VARF(6),VARF(7),VARF(8),VARF(9), 167 : * VARF(10),VARF(11),VARF(12),VARF(13), 168 : * VARF(14),VARF(15),VARF(16),VARF(17), 169 : * VARF(18),VARF(19),VARF(20)) 170 : C=================================================================== 171 : C 172 : C CONVERSION ENTIER 173 : C 174 : VARF(1)=DBLE(IFR1) 175 : C 176 : C REMPLISSAGE DU TABLEAU CONT (PARTIE ELASTIQUE) 177 : C 178 : CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2)) 179 : C SIGF(1)=SIG0(1) + DSIGT(1) 180 : SIGF(1)=SIG0(1) + XMAT(1)*XCAR(4)*DEPST(1) 181 : C SIGF(2)=SIG0(2) + DSIGT(2) 182 : SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2) 183 : C SIGF(4)=SIG0(4) + DSIGT(4) 184 : SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4) 185 : C SIGF(5)=SIG0(5) + DSIGT(5) 186 : SIGF(5)=SIG0(5) + XMAT(1)*XCAR(2)*DEPST(5) 187 : C SIGF(6)=SIG0(6) + DSIGT(6) 188 : SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6) 189 : C 190 : SIGF(3)=XEFFZ 191 : C 192 : C LES NOUVAUX VALEURS DES VARIABLES INTERNES 193 : C 194 : C ----- NOTE: VARF(2) CONTAINS THE CURRENT (GA) THAT WILL BE USED 195 : C TO COMPUTE THE TANGENT STIFFNESS MATRIX 196 : C 197 : RETURN 198 : END 199 : 200 :
© Cast3M 2003 - All rights reserved.
Disclaimer