Télécharger biline.eso

Retour à la liste

Numérotation des lignes :

  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