Télécharger reelle.eso

Retour à la liste

Numérotation des lignes :

  1. C REELLE SOURCE CB215821 16/04/21 21:18:10 8920
  2. C REELLE SOURCE NECHNECH 99/06/08
  3. SUBROUTINE REELLE(S1,DEP,DEPSITF1,DEPSITF2,DEPE,S2,ENDM,
  4. & NSTRS,IFIS,IPLA,BETJEF,NECH0,NECH1)
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. DIMENSION S1(4), S2(4), V1(4)
  9. DIMENSION DEP(4,4),DEPE(4,4)
  10. * COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,
  11. * & TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  12. * & ITER,ISIM,IBB,IGAU1,IZON
  13. * COMMON /DBETNECH0/DT,DC,ALFG,S0
  14. * COMMON /DBETNECH1/ENDL
  15. SEGMENT BETJEF
  16. REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  17. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  18. INTEGER ICT,ICC,IMOD,IVIS,ITER,
  19. & ISIM,IBB,IGAU,IZON
  20. ENDSEGMENT
  21. SEGMENT NECH0
  22. REAL*8 DT,DC,ALFG,S0,ENDO
  23. ENDSEGMENT
  24. SEGMENT NECH1
  25. REAL*8 ENDL
  26. ENDSEGMENT
  27. C
  28. *--------------------------------------------------------------------*
  29. * CALCUL LA MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE ET *
  30. * DE LA CONTRAINTE REELLE *
  31. * *
  32. * *
  33. * ENTREES: *
  34. * *
  35. * S1 CONTRAINTE EFFECTIVE *
  36. * NSTRS NOMBRE DE CONTRAINTES *
  37. * DPSTF1 DEFORMATION EQUIVALENTE EN TRACTION *
  38. * DPSTF2 DEFORMATION EQUIVALENTE EN COMPRESSION *
  39. * DEP MATRICE ELSATO-PLASTIQUE *
  40. * *
  41. * SORTIES: *
  42. * *
  43. * S2 CONTRAINTE REELLE *
  44. * DEPE MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE *
  45. * ENDO ENDOMMAGEMENT *
  46. *--------------------------------------------------------------------*
  47. C INITIALISATION
  48. C
  49. C S0 = 0.25
  50. ENDT = 0.D0
  51. ENDC = 0.D0
  52. ENDM = 0.D0
  53. NN=NSTRS-1
  54. CALL ZERO(DEPE,4,4)
  55. CALL ZERO(S2,4,1)
  56. SOM1 = 0.D0
  57. SOM2 = 0.D0
  58. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  59. C CALCUL DE L'ENDOMMAGEMENT
  60. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  61. C WRITE(*,*)'K1=',DEPSITF1
  62. C WRITE(*,*)'K2=',DEPSITF2
  63. C
  64. CALL DAMAG1(DEPSITF1,DEPSITF2,ENDC,
  65. & ENDT,IFIS,IPLA,3,BETJEF,NECH0,NECH1)
  66. C
  67. C
  68. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  69. C CALCUL DES VALEURS PROPRES DE LA MATRICE DE VECTEUR CONTRAINTE
  70. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  71. C
  72. CALL PRINC(S1,V1,NSTRS)
  73. C
  74. DO 10 I=1, NN
  75. SOM1=SOM1+(V1(I)+ABS(V1(I)))/2.D0
  76. SOM2=SOM2+ABS(V1(I))
  77. 10 CONTINUE
  78. IF(SOM2.EQ.0.D0) THEN
  79. SOM2=SOM2+1.D-10
  80. ENDIF
  81. RSIG=SOM1/SOM2
  82. RSIG=ABS(RSIG)
  83. S=S0+(1-S0)*RSIG
  84. C S = 1.D0
  85. C
  86. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  87. C CALCUL DE L'ENDOMMAGEMENT TOTAL
  88. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  89. C
  90. ENDM=1.D0-((1.D0-ENDC)*(1.D0-S*ENDT))
  91. C
  92. IF((ENDM.LT.0.D0).OR.(ENDM.GT.1D0))THEN
  93. WRITE(*,*)'ENDO=',ENDM
  94. C STOP
  95. ENDIF
  96. C
  97. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  98. C
  99. DO 20 I=1, NSTRS
  100. S2(I)=(1.D0-ENDM)*S1(I)
  101. 20 CONTINUE
  102. C
  103. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  104. C CALCUL DE LA MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE
  105. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  106. C
  107. DO 30 I=1, NSTRS
  108. DO 40 J=1, NSTRS
  109. DEPE(I,J)=(1.D0-ENDM)*DEP(I,J)
  110. 40 CONTINUE
  111. 30 CONTINUE
  112. C
  113. RETURN
  114. END
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  

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