1 : C TUFIPL SOURCE CHAT 05/01/13 03:56:07 5004 2 : SUBROUTINE TUFIPL(XM,XP,DELTAM,DELTAP,THETA,XJP,DLAM1,DLAM2,EP, 3 : &FIP,PRECIS,XM0,XP0,RAYOM,XLAM0,WORK,XJ1C,T,THETA0,EPAI,YOUN,KERRE, 4 : &MAPL) 5 : IMPLICIT INTEGER(I-N) 6 : IMPLICIT REAL*8 (A-H,O-Z) 7 : C*********************************************************************** 8 : C CE SOUS PROGRAMME EFFECTUE LA PARTIE DE L ECOULEMENT PLASTIQUE 9 : C PENDANT LAQUELLE IL N Y A PAS PROPAGATION 10 : C IL EST APPELE PAR TUFPLA 11 : C 12 : C ENTREES:XM VALEUR DU MOMENT AU DEBUT DU PAS 13 : C XP VALEUR DE L EFFORT NORMAL AU DEBUT DU PAS 14 : C DELTAM INCREMENT DE MOMENT 15 : C DELTAP INCREMENT D EFFORT NORMAL 16 : C THETA ANGLE DE LA FISSURE AU DEBUT DU PAS 17 : C XJP VALEUR DE LA PARTIE PLASTIQUE DE J 18 : C 19 : C SORTIES:XM VALEUR DU MOMENT A LA FIN DU PAS 20 : C XP VALEUR DE L EFFORT NORMAL A LA FIN DU PAS 21 : C DELTAM PARTIE DE L INCREMENT DE MOMENT RESTANT 22 : C A ECOULER EN PLASTICITE 23 : C DELTAP PARTIE DE L INCREMENT D EFFORT NORMAL 24 : C RESTANT A ECOULER EN PLASTICITE 25 : C======================================================================= 26 : -INC CCOPTIO 27 : SEGMENT/WORK/(SIG(NCOURB)*D,XLAM(NCOURB)*D) 28 : XMS=XM 29 : XPS=XP 30 : DMS=DELTAM 31 : DPS=DELTAP 32 : XJPS=XJP 33 : IF(IIMPI.EQ.999) WRITE (IOIMP,*) 'ENTREE DANS TUFIPL' 33 : CALL TUFIPA(XM,XP,DELTAM,DELTAP,THETA0,XJP,EP,FIP, 34 : $DLAM1,DLAM2,PRECIS,XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE, 35 : $THETA) 35 : CALL TUFIJE(XM,XP,THETA0,XJE,RAYOM,EPAI,YOUN) 36 : IF(MAPL.EQ.14) THEN 36 : CALL TUFIC2(XJP,XJE,THETA,IRE,PRECIS,XJ1C,T,RAYOM,THETA0) 37 : ELSE 38 : IRE=0 39 : ENDIF 40 : IF(IRE.EQ.0) THEN 41 : IF(IIMPI.EQ.999) WRITE (IOIMP,*) 'SORTIE DE TUFIPL' 42 : RETURN 43 : ENDIF 44 : XL1=0.75D0 45 : XM1=XMS 46 : XP1=XPS 47 : DM1=XL1*DMS 48 : DP1=XL1*DPS 48 : CALL TUFIPA(XM1,XP1,DM1,DP1,THETA,XJPS,EP,FIP,DLAM1,DLAM2,PRECIS, 49 : &XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,THETA) 50 : F1=CRIT2(XM1,XP1,THETA,XJPS,XJ1C,T,RAYOM,THETA0,EPAI,YOUN) 51 : XL2=0.5D0 52 : XM2=XMS 53 : XP2=XPS 54 : DM2=XL2*DMS 55 : DP2=XL2*DPS 55 : CALL TUFIPA(XM2,XP2,DM2,DP2,THETA,XJPS,EP,FIP,DLAM1,DLAM2,PRECIS, 56 : &XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,THETA) 57 : F2=CRIT2(XM2,XP2,THETA,XJPS,XJ1C,T,RAYOM,THETA0,EPAI,YOUN) 58 : 1 CONTINUE 59 : FP=(F2-F1)/(XL2-XL1) 60 : XL3=XL2-F2/FP 61 : XM3=XMS 62 : XP3=XPS 63 : DM3=XL3*DMS 64 : DP3=XL3*DPS 64 : CALL TUFIPA(XM3,XP3,DM3,DP3,THETA,XJPS,EP,FIP,DLAM1,DLAM2,PRECIS, 65 : &XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,THETA) 66 : F3=CRIT2(XM3,XP3,THETA,XJPS,XJ1C,T,RAYOM,THETA0,EPAI,YOUN) 67 : CRI=ABS((F3-F2)/F2) 68 : IF(CRI.GE.PRECIS) THEN 69 : XL1=XL2 70 : XL2=XL3 71 : F1=F2 72 : F2=F3 73 : GO TO 1 74 : ELSE 75 : XJP=XJPS 76 : XM=XM3 77 : XP=XP3 78 : DELTAM=(1.D0-XL3)*DMS 79 : DELTAP=(1.D0-XL3)*DPS 80 : IF(IIMPI.EQ.999) WRITE (IOIMP,*) 'SORTIE DE TUFIPL' 81 : RETURN 82 : ENDIF 83 : END 84 :
© Cast3M 2003 - All rights reserved.
Disclaimer