Télécharger plaqp.eso

Retour à la liste

Numérotation des lignes :

  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