1 : C TAKEPO SOURCE CHAT 05/01/13 03:31:59 5004 2 : SUBROUTINE TAKEPO(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 DE LA ROTULE 7 : C CAS DU TAKEDA TRILINEAIRE SUR LE MOMENT 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 SRDZ = XCAR(6) TRANSVERSE SHEAR'S REDUCED SECTION 40 : C----------------------------------------------------------------------- 41 : C VARIABLES INTERNES 42 : C ================== 43 : C 44 : C IFR1 = VARF(1) 45 : C EITA = 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 : 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(5) 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(5) / EIYY0 * VARF(2) 158 : DP1 = DEPST(5) * 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(4)*DEPST(1) 180 : C SIGF(2)=SIG0(2) + DSIGT(2) 181 : SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2) 182 : C SIGF(3)=SIG0(3) + DSIGT(3) 183 : SIGF(3)=SIG0(3) + CISA*XCAR(6)*DEPST(3) 184 : C SIGF(4)=SIG0(4) + DSIGT(4) 185 : SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4) 186 : C SIGF(6)=SIG0(6) + DSIGT(6) 187 : SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6) 188 : C 189 : SIGF(5)=XMOMY 190 : C 191 : C LES NOUVAUX VALEURS DES VARIABLES INTERNES 192 : C 193 : C ----- NOTE: VARF(2) CONTAINS THE CURRENT (EI) THAT WILL BE USED 194 : C TO COMPUTE THE TANGENT STIFFNESS MATRIX 195 : C 196 : RETURN 197 : END 198 : 199 :
© Cast3M 2003 - All rights reserved.
Disclaimer