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