Télécharger la plaquette Cast3M


fr en

Télécharger biline.eso

Retour à la liste

Numérotation des lignes : 

biline

  1. C BILINE    SOURCE    CHAT      05/01/12    21:40:08     5004
  2.       SUBROUTINE BILINE(DELTAD,FY,STIFE,STIFP,RIG,FE0,FEP0,F,DELDPL)
  3. C-----------------------------------------------------------------------
  4. C
  5. C SUBROUTINE FOR ELASTIC PLASTIC (WITH HARDENING) RELATIONSHIP.
  6. C
  7. C
  8. C A PERFECT E.P. SCALAR MODEL WITH INITIAL STIFNESS STIFEQ AND
  9. C YIELDING FYEQ IS ADDED IN PARALEL TO AN ELASTIC SCALAR MODEL
  10. C OF STIFNESS STIFP.
  11. C
  12. C INPUT
  13. C
  14. C    DELTAD  DISPLACEMENT INCREMENT
  15. C    FY      YIELDING FORCE (FY=STIFE*FYEQ/STIFEQ)
  16. C    STIFE   ELASTIC STIFNESS (STIFE=STIFEQ+STIFP)
  17. C    STIFP   PLASTIC STIFNESS
  18. C    FE0     INITIAL FORCE IN THE ELASTIC MODEL
  19. C    FEP0    INITIAL FORCE IN THE E.P. MODEL
  20. C
  21. C OUTPUT
  22. C
  23. C    RIG     FINAL STIFFNESS
  24. C    F       FINAL FORCE
  25. C    FE0     FINAL FORCE IN THE ELASTIC MODEL
  26. C    FEP0    FINAL FORCE IN THE ELASTIC MODEL
  27. C    DELDPL  PLASTIC DEFORMATION INCREMENT
  28. C
  29. C WARNING FE0 AND FEP0 SHOULD BE INITIALIZED !!
  30. C
  31. C-----------------------------------------------------------------------
  32. C VERY OLD (CRAZY ?) STUFF SLIGHTLY MODIFIED FOR DELDPL
  33. C AND TAKEN AS AN EXEMPLE OF "EXISTING" MODEL
  34. C-----------------------------------------------------------------------
  35.       IMPLICIT INTEGER(I-N)
  36.       IMPLICIT REAL*8(A-H,O-Z)
  37. C
  38. C  COMPUTATION OF THE STIFNESS AND YIELDING OF THE E.P. MODEL
  39. C
  40.       STIFEQ=STIFE-STIFP
  41.       FYEQ=FY*STIFEQ/STIFE
  42. C
  43. C  ELASTIC PREDICTION OF THE FORCE IN THE E.P. MODEL ...
  44. C
  45.       FEP=FEP0+STIFEQ*DELTAD
  46. C
  47. C  ... AND POSSIBLE RETURN
  48. C
  49.       IF (ABS(FEP).GT.FYEQ)THEN
  50.         FEP=SIGN(FYEQ,FEP)
  51.         RIG=STIFP
  52.         DELDPL=DELTAD-(FEP-FEP0)/STIFEQ
  53.       ELSE
  54.         RIG=STIFE
  55.         DELDPL=0.D0
  56.       ENDIF
  57.       FEP0=FEP
  58. C
  59. C  FORCE IN THE ELASTIC MODEL
  60. C
  61.       FE=FE0+STIFP*DELTAD
  62.       FE0=FE
  63. C
  64. C  TOTAL FORCE
  65. C
  66.       F=FE+FEP
  67.       RETURN
  68.       END
  69.  
  70.  
  71.  

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