Télécharger la plaquette Cast3M


fr en

Télécharger tufipl.eso

Retour à la liste

Numérotation des lignes : 

tufipl

  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.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29.       SEGMENT/WORK/(SIG(NCOURB)*D,XLAM(NCOURB)*D)
  30.       XMS=XM
  31.       XPS=XP
  32.       DMS=DELTAM
  33.       DPS=DELTAP
  34.       XJPS=XJP
  35.       IF(IIMPI.EQ.999) WRITE (IOIMP,*) 'ENTREE DANS TUFIPL'
  36.       CALL TUFIPA(XM,XP,DELTAM,DELTAP,THETA0,XJP,EP,FIP,
  37.      $DLAM1,DLAM2,PRECIS,XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,
  38.      $THETA)
  39.       CALL TUFIJE(XM,XP,THETA0,XJE,RAYOM,EPAI,YOUN)
  40.       IF(MAPL.EQ.14) THEN
  41.         CALL TUFIC2(XJP,XJE,THETA,IRE,PRECIS,XJ1C,T,RAYOM,THETA0)
  42.       ELSE
  43.         IRE=0
  44.       ENDIF
  45.       IF(IRE.EQ.0) THEN
  46.           IF(IIMPI.EQ.999) WRITE (IOIMP,*) 'SORTIE DE TUFIPL'
  47.           RETURN
  48.       ENDIF
  49.       XL1=0.75D0
  50.       XM1=XMS
  51.       XP1=XPS
  52.       DM1=XL1*DMS
  53.       DP1=XL1*DPS
  54.       CALL TUFIPA(XM1,XP1,DM1,DP1,THETA,XJPS,EP,FIP,DLAM1,DLAM2,PRECIS,
  55.      &XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,THETA)
  56.       F1=CRIT2(XM1,XP1,THETA,XJPS,XJ1C,T,RAYOM,THETA0,EPAI,YOUN)
  57.       XL2=0.5D0
  58.       XM2=XMS
  59.       XP2=XPS
  60.       DM2=XL2*DMS
  61.       DP2=XL2*DPS
  62.       CALL TUFIPA(XM2,XP2,DM2,DP2,THETA,XJPS,EP,FIP,DLAM1,DLAM2,PRECIS,
  63.      &XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,THETA)
  64.       F2=CRIT2(XM2,XP2,THETA,XJPS,XJ1C,T,RAYOM,THETA0,EPAI,YOUN)
  65. 1     CONTINUE
  66.       FP=(F2-F1)/(XL2-XL1)
  67.       XL3=XL2-F2/FP
  68.       XM3=XMS
  69.       XP3=XPS
  70.       DM3=XL3*DMS
  71.       DP3=XL3*DPS
  72.       CALL TUFIPA(XM3,XP3,DM3,DP3,THETA,XJPS,EP,FIP,DLAM1,DLAM2,PRECIS,
  73.      &XM0,XP0,RAYOM,XLAM0,WORK,EPAI,YOUN,KERRE,THETA)
  74.       F3=CRIT2(XM3,XP3,THETA,XJPS,XJ1C,T,RAYOM,THETA0,EPAI,YOUN)
  75.       CRI=ABS((F3-F2)/F2)
  76.       IF(CRI.GE.PRECIS) THEN
  77.            XL1=XL2
  78.            XL2=XL3
  79.            F1=F2
  80.            F2=F3
  81.            GO TO 1
  82.       ELSE
  83.            XJP=XJPS
  84.            XM=XM3
  85.            XP=XP3
  86.            DELTAM=(1.D0-XL3)*DMS
  87.            DELTAP=(1.D0-XL3)*DPS
  88.            IF(IIMPI.EQ.999) WRITE (IOIMP,*) 'SORTIE DE TUFIPL'
  89.            RETURN
  90.       ENDIF
  91.       END
  92.  
  93.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales