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