Télécharger dechin.eso

Retour à la liste

Numérotation des lignes :

  1. C DECHIN SOURCE CB215821 16/04/21 21:16:20 8920
  2. C DECHIN SOURCE INSL 24/10/96
  3. SUBROUTINE DECHIN(SIGRX,EPSRX,ETS,EDC,EPST,IFISU,JFISU,RTM,
  4. 1 EPSC,DEFR,TANG,IPLA,EQSTR,EPSEQ,EX,RB,ALPHA,EPO1,EMAX,ICAL,
  5. 2 IBB1,IGAU1,IDIR)
  6. C
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8(A-H,O-Z)
  9. C
  10. C CE SOUS-PROGRAMME CALCULE LA CONTRAINTE REELLE S1X ET INITIALISE
  11. C LES VARIABLES INTERNES LORSQU'IL Y A DECHARGE A PARTIR D'UN
  12. C POINT INTEGRE. ON SE BRANCHE ALORS SUR LA COURBE 1-D CORRESPONDANT
  13. C A ICU=3 (PARABOLE-SOFTENING LINEAIRE)
  14. C
  15. C PF : POINT FOCAL
  16. C EPST : DEFORMATION EN COMPRESSION CORRESPONDANT A -RBT
  17. C APRES DECHARGE
  18. C EPSC : DEFORMATION EN COMPRESSION A L'INSTANT DE LA DECHARGE
  19. C RTMC : RESISTANCE EN TRACTION MODIFIEE si decharge en compression
  20. C RTM : RESISTANCE EN TRACTION ACTUELLE
  21. C EDC : PENTE DECHARGE EN COMPRESSION
  22. C DEFR : DEFORMATION RESIDUELLE APRES DECHARGE EN COMPRESSION
  23. C ETS : PENTE EN TENSION-SOFTENING
  24. C TANG : PENTE ACTUELLE
  25. C IPLA : INDICATEUR DE PLASTICITE/ =0 comportement pre-pic en compression
  26. C =2 softening en compression
  27. C =4 rupture en compression (epsi>emax)
  28. C
  29. C------------------------------------------------------------------
  30. ZERO=0.D0
  31. UN=1.D0
  32. RBT=ALPHA*RB
  33. PF=RB/EX
  34. EPOC=-1.D0*EPO1
  35. IF(EPSRX.LT.(-1.D0*EMAX)) THEN
  36. C
  37. C ON SUPPOSE QUE LE POINT EST ROMPU EN COMPRESSION
  38. C
  39. IPLA=4
  40. SIGRX =ZERO
  41. RTM =ZERO
  42. TANG=ZERO
  43. ETS =ZERO
  44. EQSTR=ZERO
  45. RETURN
  46. ENDIF
  47. C
  48. IF(EPSRX .LE. EPOC) THEN
  49. C**********************************************************************
  50. C POINT EN REGIME POST-PIC (COMPRESSION)
  51. C**********************************************************************
  52. IPLA=2
  53. AA=-1.D0*RB/(EMAX-EPO1)
  54. BB=AA*EMAX
  55. DELTA=SIGRX-(AA*EPSRX+BB)
  56. IF(DELTA .GT. ZERO) THEN
  57. C
  58. C POINT EN DECHARGE A PARTIR DE LA COURBE DE COMPRESSION
  59. C
  60. CALL PARDEC(EPSRX,SIGRX,EPSC,EQSTR,EDC,IPLA,DEFR,RB,ALPHA,EPO1,
  61. & EX,EMAX,ICAL)
  62. EPSEQ=ABS(EPSC)
  63. GOTO 100
  64. ELSE
  65. C
  66. C POINT EN CHARGE (SOFTENING)
  67. C COURBE DE COMPRESSION TRANSLATEE
  68. C
  69. SEQ=ABS(SIGRX)
  70. IPLA=2
  71. CALL EPSEQU(EPEQTR,SEQ,IPLA,ICAL,1,EX,RB,EPO1,EPO1,EMAX)
  72. SIGRX=-1.D0*SEQ
  73. EPSRX=-1.D0*EPEQTR
  74. EQSTR=SEQ
  75. EPSEQ=EPEQTR
  76. ENDIF
  77. RETURN
  78. ENDIF
  79. C**********************************************************************
  80. C POINT EN REGIME PRE-PIC (COMPRESSION)
  81. C***********************************************************************
  82. DELTA=1.D0
  83. DEFRLI = EPSRX-SIGRX/EX
  84. IF(ABS(DEFRLI).LE.1.D-9) DEFRLI=0.D0
  85. IF(EPSRX.LE.ZERO) THEN
  86. EPEQ=ABS(EPSRX)
  87. CALL CDCY(EPEQ,SEQ2,IPL,TAN,EX,RB,EMAX,EPO1,ICAL)
  88. DELTA=-1.D0*SEQ2-SIGRX
  89. ENDIF
  90. IF(DELTA.LT.ZERO.OR.(DEFRLI.LT.ZERO.AND.EPSRX.GE.ZERO)) THEN
  91. C
  92. C POINT EN DECHARGE A PARTIR DE LA COURBE DE COMPRESSION
  93. C
  94. CALL PARDEC(EPSRX,SIGRX,EPSC,EQSTR,EDC,IPLA,DEFR,RB,ALPHA,EPO1,
  95. & EX,EMAX,ICAL)
  96. EPSEQ=ABS(EPSC)
  97. GOTO 100
  98. ENDIF
  99. C--------------------------------------------------------------
  100. C POINT CONSIDERE INITIALLEMENT TENDU
  101. C--------------------------------------------------------------
  102. EPSF0=RTM/EX
  103. IF(EPSRX.GE.ZERO .AND. SIGRX.GE.ZERO) THEN
  104. IF(SIGRX.GE.RTM) THEN
  105. SIGRX=RTM
  106. IFISU=1
  107. JFISU=1
  108. ENDIF
  109. IF(EPSRX.GE.EPSF0) THEN
  110. EPSRX=EPSF0+(RTM-SIGRX)/ETS
  111. IFISU=1
  112. JFISU=1
  113. ELSE
  114. EPSRX=SIGRX/EX
  115. ENDIF
  116. RETURN
  117. ENDIF
  118. C==============================================================
  119. C POINT INITIALLEMENT COMPRIME
  120. C COURBE DE COMPRESSION TRANSLATEE
  121. C
  122. SEQ=ABS(SIGRX)
  123. IPLA=0
  124. CALL EPSEQU(EPEQTR,SEQ,IPLA,ICAL,1,EX,RB,EPO1,EPO1,EMAX)
  125. SIGRX=-1.D0*SEQ
  126. EPSRX=-1.D0*EPEQTR
  127. EQSTR=SEQ
  128. EPSEQ=EPEQTR
  129. RETURN
  130. C==============================================================
  131. 100 CONTINUE
  132. C
  133. C RESISTANCE EN TRACTION MODIFIEE
  134. C
  135. EPOC=-1.D0*EPO1
  136. IF(ABS(DEFR).LT.1.D-10) DEFR=0.D0
  137. IF(DEFR.LT.EPOC) THEN
  138. RTMC=0.D0
  139. ELSE
  140. AS2 = RBT/(EPO1+RBT/EX)
  141. IF(ABS(EDC-AS2).GT.1.D-5) THEN
  142. EPS1 = (EDC*DEFR-AS2*EPOC)/(EDC-AS2)
  143. RTMC = EDC*(EPS1-DEFR)
  144. ELSE
  145. RTMC = 0.D0
  146. ENDIF
  147. C
  148. IF(RTMC.GT.RBT.AND.RTMC.LT.(1.00001D0*RBT)) RTMC = RBT
  149. IF(RTMC.GT.RBT.OR.RTMC.LT.0.D0) THEN
  150. WRITE(*,94) IBB1,IGAU1,IDIR
  151. WRITE(*,*) 'RTMC = ',RTMC, 'RBT = ',RBT
  152. WRITE(*,*) 'DEFR = ',DEFR,'EPO1 = ',EPO1
  153. WRITE(*,*) 'EPSRX= ',EPSRX,' SIGRX= ',SIGRX
  154. WRITE(*,*) 'EDC = ',EDC,' AS2 = ',AS2,' EPS1 = ',EPS1
  155. WRITE(*,*) 'ARRET ---> RTMC ?? !!! (SUBR DECHINT2) '
  156. RTMC=RBT
  157. STOP
  158. ENDIF
  159. ENDIF
  160. IF(RTMC .LT. ZERO) RTMC=ZERO
  161. IF(RTMC.LT.RTM) RTM=RTMC
  162. C
  163. C CALCUL DE LA DEFORMATION CORRESPONDANT A -RBT
  164. C
  165. IF(ABS(EDC).LT.1.D-5) EDC=EX
  166. EPST=DEFR-RBT/EDC
  167. C
  168. IF(EDC.LT.0.D0) THEN
  169. WRITE(*,94) IBB1,IGAU1,IDIR
  170. WRITE(*,*) 'RTMC = ',RTMC, 'RBT = ',RBT
  171. WRITE(*,*) 'DEFR = ',DEFR,'EPO1 = ',EPO1
  172. WRITE(*,*) 'EPSRX= ',EPSRX,' SIGRX= ',SIGRX
  173. WRITE(*,*) 'EDC = ',EDC,' AS2 = ',AS2,' EPS1 = ',EPS1
  174. WRITE(*,*) 'ARRET ---> EDC<0 ?? !!! (SUBR DECHINT2) '
  175. STOP
  176. ENDIF
  177. C
  178. C Dans le cas ou la decharge s'effectue dans la branche descendante
  179. C en-dessous du niveau |RBT| on m‚morise EPSC
  180. C
  181. IF(IPLA .EQ. 2 .AND. EQSTR .LE. RBT) EPST=EPSC
  182. C
  183. 94 FORMAT(' ELEME=',I3,'POINT=',I2,' IDIR=',I2)
  184. RETURN
  185. END
  186.  
  187.  
  188.  
  189.  
  190.  

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