Télécharger bilifo.eso

Retour à la liste

Numérotation des lignes :

bilifo
  1. C BILIFO SOURCE CHAT 05/01/12 21:40:04 5004
  2. SUBROUTINE BILIFO(SIG0,DEPST,VAR0,XMAT,XCAR,SIGF,VARF,DEFP)
  3. C-----------------------------------------------------------------------
  4. C
  5. C DRIVER OF THE BILINEAR PLASTIC TRANSVERSE SHEAR MODEL
  6. C
  7. C INPUT:
  8. C
  9. C SIG0 INITIAL GENERALIZED STRESS
  10. C DEPST INCREMENT OF GENERALIZED STRAIN
  11. C XMAT MATERIAL CHARACTERISTICS
  12. C XCAR GEOMETRIC CHARACTERISTICS
  13. C VAR0 INITIAL INTERNAL VARIABLE
  14. C
  15. C OUTPUT:
  16. C
  17. C SIGF FINAL GENERALIZED STRESS
  18. C VARF FINAL INTERNAL VARIABLE
  19. C DEFP INCREMENT OF PLASTIC DEFORMATION
  20. C
  21. C GENERALIZED STRESS (.=0 OR F):
  22. C
  23. C EFFZ = SIG.(3)
  24. C
  25. C MATERIAL CHARACTERISTICS:
  26. C
  27. C E = XMAT(1) ELASTIC YOUNG MODULUS
  28. C NU = XMAT(2) POISSON'S RATIO
  29. C EAYI = XMAT(5) PLASTIC MODULUS
  30. C YEFF = XMAT(6) YIELDING TRANSVERSE SHEAR FORCE
  31. C
  32. C GEOMETRIC CHARACTERISTICS:
  33. C
  34. C INRX = XCAR(1) X-INERTIAL TORSION
  35. C INRY = XCAR(2) Y-INERTIAL BENDING
  36. C INRZ = XCAR(3) Z-INERTIAL BENDING
  37. C SECT = XCAR(4) SECTION
  38. C SRDY = XCAR(5) TRANSVERSE SHEAR'S REDUCED SECTION
  39. C SECZ = XCAR(6) Z-TRANSVERSE SHEAR SECTION <------
  40. C
  41. C INTERNAL VARIABLE (.=0 OR F):
  42. C
  43. C EPSE = VAR.(1) CUMULATED PLASTIC TRANSVERSE SHEAR STRAIN
  44. C GATA = VAR.(2) TANGENT Z-SHEAR STIFFNESS
  45. C EFEZ = VAR.(3) ELASTIC SHEAR FORCE
  46. C CISZ = VAR.(4) TRANSVERSE SHEAR STRAIN
  47. C
  48. C-----------------------------------------------------------------------
  49. C P. PEGON 26/11/92 renewed new ECOU 24/11/95
  50. C-----------------------------------------------------------------------
  51. IMPLICIT INTEGER(I-N)
  52. IMPLICIT REAL*8(A-H,O-Z)
  53.  
  54. -INC PPARAM
  55. -INC CCOPTIO
  56. C
  57. DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),
  58. . SIGF(*),VARF(*) ,DEFP(*)
  59. C
  60. C-----------------------------------------------------------------------
  61. C
  62. C ELASTIC STRESS INCREMENT FOR ALL STRESS COMPONENT EXCEPT MOMY
  63. C
  64. CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2))
  65. C SIGF(1)=SIG0(1) + DSIGT(1)
  66. SIGF(1)=SIG0(1) + XMAT(1)*XCAR(4)*DEPST(1)
  67. C SIGF(2)=SIG0(2) + DSIGT(2)
  68. SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2)
  69. C SIGF(4)=SIG0(4) + DSIGT(4)
  70. SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4)
  71. C SIGF(5)=SIG0(5) + DSIGT(5)
  72. SIGF(5)=SIG0(5) + XMAT(1)*XCAR(2)*DEPST(5)
  73. C SIGF(6)=SIG0(6) + DSIGT(6)
  74. SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6)
  75. C
  76. C ELASTIC AND PLASTIC SHEAR STIFNESS
  77. C
  78. XMODG=0.5D0 / (1.D0 + XMAT(2))
  79. GAZ = XMAT(1)*XCAR(6)*XMODG
  80. GYAZ = XMAT(5)*XCAR(6)*XMODG
  81. C
  82. C Z-TRANSVERSE SHEAR INCREMENT
  83. C
  84. C DCISAZ = DSIGT(3) / GAZ
  85. DCISAZ = DEPST(3)
  86. C
  87. C CALL TO THE SCALAR BILINEAR PLASTIC MODEL COMPUTING MOMY
  88. C
  89. FEP0 = SIG0(3) - VAR0(3)
  90. FE0 = VAR0(3)
  91. CALL BILINE(DCISAZ,XMAT(6),GAZ, GYAZ, VARF(2),
  92. > FE0,FEP0,SIGF(3),DEPSE)
  93. C
  94. C FINAL INTERNAL VARIABLE
  95. C
  96. VARF(1)=VAR0(1)+DEPSE
  97. VARF(3)=FE0
  98. VARF(4)=VAR0(4)+DCISAZ
  99. C
  100. C PLASTIC DEFORMATION INCREMENT
  101. C
  102. DEFP(1) = 0.D0
  103. DEFP(2) = 0.D0
  104. DEFP(3) = DEPSE
  105. DEFP(4) = 0.D0
  106. DEFP(5) = 0.D0
  107. DEFP(6) = 0.D0
  108. RETURN
  109. END
  110.  
  111.  
  112.  

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