Télécharger la plaquette Cast3M


fr en

Télécharger plaqp.eso

Retour à la liste

Numérotation des lignes : 

plaqp

  1. C PLAQP     SOURCE    CB215821  16/04/21    21:18:02     8920
  2.       SUBROUTINE PLAQP(EC,QT,FNE,SNNP10,EPCO,XLAM3,SNNP11,EPCOF,XLAM3F)
  3. C-----------------------------------------------------------------------
  4. C
  5. C      DANS LE CADRE DU MODELE 3D MOHR COULOMB ELEMENTS JOINTS JOI1
  6. C       PLASTICITE EN COMPRESSION ECROUISSAGE ISOTROPE LINEAIRE
  7. C
  8. C-----------------------------------------------------------------------
  9. C ENTREES
  10. C     EC      = RAIDEUR A L'ORIGINE ELASTIQUE   °      °
  11. C     QT      = RAIDEUR TANGENTE RELIANT        F A EPS_TOTAL
  12. C     FNE     = LIMITE D'ELASTICITE
  13. C
  14. C     SNNP10  = ESTIMATION ELASTIQUE DE LA CONTRAINTE FINALE
  15. C     EPCO    = DEFORMATION PLASTIQUE EN COMPRESSION AU DEBUT DU PAS
  16. C     XLAM3   = MULTIPLICATEUR PLASTIQUE AU DEBUT DU PAS
  17. C
  18. C SORTIES
  19. C     SNNP11  = CONTRAINTE FINALE A LA FIN DU PAS
  20. C     EPCOF   = DEFORMATION PLASTIQUE EN COMPRESSION A LA FIN DU PAS
  21. C     XLAM3F  = MULTIPLICATEUR PLASTIQUE A LA FIN DU PAS
  22. C-----------------------------------------------------------------------
  23. C
  24.       IMPLICIT INTEGER(I-N)
  25.       IMPLICIT REAL*8(A-H,O-Z)
  26. C
  27.       PARAMETER(PRECIS=1D-6,XZER=0D0,IMAX=200)
  28. C
  29. C
  30.       H=EC*QT/(EC-QT)
  31. C
  32.       SNNP11=SNNP10
  33.       EPCOF=EPCO
  34.       XLAM3F=XLAM3
  35. C
  36.       F3=ABS(SNNP10)-H*XLAM3-FNE
  37.  
  38.       IF (F3.LE.PRECIS) GOTO 900
  39. C
  40. C
  41.       I=0
  42.    10 CONTINUE
  43.       I=I+1
  44.       IF (I.GT.IMAX) THEN
  45.          PRINT *,' LE CRITERE PLASTIQUE EN COMPRESSION EST VIOLE '
  46.          PRINT *,' DIVERGENCE APRES ',IMAX,' ITERATIONS INTERNES'
  47.          STOP
  48.       ENDIF
  49. C
  50.       DF3DF=SNNP11/(ABS(SNNP11))
  51.       DLAM3=F3/(EC+QT)
  52.       XLAM3F=XLAM3F+DLAM3
  53.       DEPCO=DLAM3*DF3DF
  54.       EPCOF=EPCOF+DEPCO
  55.       SNNP11=SNNP11-EC*DLAM3*DF3DF
  56.       F3=ABS(SNNP11)-H*XLAM3F-FNE
  57.       IF (ABS(F3).GT.PRECIS) GOTO 10
  58. C
  59.   900 CONTINUE
  60. C
  61.       RETURN
  62.       END
  63.  
  64.  
  65.  
  66.  

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