Télécharger colbec.eso

Retour à la liste

Numérotation des lignes :

  1. C COLBEC SOURCE BP208322 17/03/01 21:16:01 9325
  2. SUBROUTINE COLBEC(WRK52,WRK53,WRK54,WRKK2,NSTRS1,NVARI,ICARA)
  3. C CLBECO SOURCE AM 98/12/23 21:20:31 3409
  4. C
  5. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  6. C % %
  7. C % subroutine integrant le comportement du beton %
  8. C % %
  9. C % - determine les nouvelles valeurs des variables %
  10. C % d'histoire YLIM1 YLIM2 D1 et D2 %
  11. C % %
  12. C % %
  13. C % en fonction du vecteur des deformations %
  14. C % %
  15. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  16. C
  17. C
  18. C
  19. c SUBROUTINE CLBECO(WRK0,WRK1,WRKK2,WRK5,NSTRS,NVARI,NMATT,
  20. c 1 ISTEP,ICARA,KERRE)
  21. C
  22. C
  23. C
  24. C
  25. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  26. C % %
  27. C % %
  28. C % declaration des variables %
  29. C % %
  30. C % %
  31. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  32. C
  33. C
  34. C
  35. IMPLICIT INTEGER(I-N)
  36. IMPLICIT REAL*8(A-H,O-Z)
  37. -INC CCOPTIO
  38. -INC DECHE
  39. SEGMENT WRKK2
  40. REAL*8 EPSILI(NSTRS1)
  41. END SEGMENT
  42.  
  43. SEGMENT WRK3
  44. REAL*8 EPSILO(NSTRS1)
  45. END SEGMENT
  46. *
  47. REAL*8 YLIM1,YLIM2,DOM1,DOM2,Y1,Y2,Z1,Z2
  48. *
  49. REAL*8 EPS33(3,3),XROTA(3,3),EPSA(3),EPSN(3)
  50. REAL*8 SIGNN(3),SIGA(3),SIGMAN(3,3)
  51. INTEGER JSIGN(3)
  52. *
  53. *
  54. LOGICAL CONV
  55. REAL*8 UN,DEUX,ZERO
  56. PARAMETER (ZERO=0.D0)
  57. PARAMETER (UN=1.D0)
  58. PARAMETER (DEUX=2.D0)
  59. *
  60. * -----> RECUPERATION DES VARIABLES INTERNES
  61. *
  62. YLIMA1 = VAR0(1)
  63. YLIMA2 = VAR0(2)
  64. DOM1 = VAR0(3)
  65. DOM2 = VAR0(4)
  66. E = XMAT(1)
  67. *
  68. *
  69. *
  70. Z1=UN/(UN-DOM1)
  71. Z2=UN/(UN-DOM2)
  72. *
  73. * -----> ON DIAGONALISE LES ANCIENNES DEFORMATIONS
  74. *
  75. CALL ENDOCA(EPSILI,EPS33,2)
  76. CALL JACOB3(EPS33,2,EPSA,XROTA)
  77. * PRINT*,'ANCIENNES DEFORMATIONS PPALES',EPSA(1),EPSA(2),EPSA(3)
  78. *
  79. * -----> ON DIAGONALISE LES ANCIENNES CONTRAINTES
  80. *
  81. CALL ENDOCA(SIG0,EPS33,1)
  82. CALL JACOB3(EPS33,2,SIGA,XROTA)
  83. *
  84. * -----> ON CALCULE LES DEFORMATIONS TOTALES
  85. *
  86. SEGINI WRK3
  87. DO ISTRS=1,NSTRS1
  88. EPSILO(ISTRS)=EPSILI(ISTRS)+DEPST(ISTRS)
  89. * print*,'ISTRS',ISTRS,'EPSILO',EPSILO(ISTRS)
  90. END DO
  91. *
  92. * -----> ON LES DIAGONALISE
  93. *
  94. CALL ENDOCA(EPSILO,EPS33,2)
  95. CALL JACOB3(EPS33,2,EPSN,XROTA)
  96. * PRINT*,'DEFORMATIONS PPALES',EPSN(1),EPSN(2),EPSN(3)
  97. *
  98. * -----> ON INITIALISE Y1 ET Y2
  99. *
  100. IF (SIGA(1) .GE. ZERO) THEN
  101. JSIGN(1)=1
  102. ELSE
  103. JSIGN(1)=0
  104. END IF
  105. IF (SIGA(2) .GE. ZERO) THEN
  106. JSIGN(2)=1
  107. ELSE
  108. JSIGN(2)=0
  109. END IF
  110. TRSIG=SIGA(1)+SIGA(2)
  111. IF (TRSIG .GE. ZERO) THEN
  112. FSIG=TRSIG
  113. ELSE IF (TRSIG.GE.-SIGF) THEN
  114. FSIG=TRSIG*(UN+TRSIG/(DEUX*SIGF))
  115. ELSE
  116. FSIG=-SIGF*TRSIG/DEUX
  117. END IF
  118. YA1=( JSIGN(1)*SIGA(1)*SIGA(1)+JSIGN(2)*SIGA(2)*SIGA(2) )
  119. & *Z1*Z1/(DEUX*E)
  120. YA2=( (1-JSIGN(1))*SIGA(1)*SIGA(1)+
  121. & (1-JSIGN(2))*SIGA(2)*SIGA(2) )
  122. & *Z2*Z2/(DEUX*E)
  123. C
  124. C
  125. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  126. C % %
  127. C % DECOUPAGE DE L'INCREMENT DE DEFORMATION %
  128. C % %
  129. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  130. c
  131. NEPS=1
  132. CONV=.FALSE.
  133. DO WHILE ((.NOT. CONV).AND.(NEPS.LE.1000))
  134. SIG1 = SIGA(1)
  135. SIG2 = SIGA(2)
  136. YLIM1=YLIMA1
  137. YLIM2=YLIMA2
  138. IF (SIG1.GE.0.D0)THEN
  139. JSIGN(1)=1
  140. ELSE
  141. JSIGN(1)=0
  142. END IF
  143. IF (SIG2.GE.0.D0)THEN
  144. JSIGN(2)=1
  145. ELSE
  146. JSIGN(2)=0
  147. END IF
  148. Z1=1.D0/(1.D0-DOM1)
  149. Z2=1.D0/(1.D0-DOM2)
  150. Y1 = YA1
  151. Y2 = YA2
  152. DEPS1=(EPSN(1)-EPSA(1))/NEPS
  153. DEPS2=(EPSN(2)-EPSA(2))/NEPS
  154. IEPS=0
  155. CONV=.TRUE.
  156. DO WHILE (CONV .AND. (IEPS.LT.NEPS))
  157. IEPS=IEPS+1
  158. EPSP1=EPSA(1)+IEPS*DEPS1
  159. EPSP2=EPSA(2)+IEPS*DEPS2
  160. c
  161. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  162. C % %
  163. C % CALCUL DES VALEURS DES VARIABLES INTERNES %
  164. C % %
  165. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  166. C
  167. C
  168. c
  169. nmatt1 = nmatt
  170. CALL CMODCL(WRK52,nmatt1,EPSP1,EPSP2,YLIM1,YLIM2,Z1,Z2,Y1,Y2,
  171. 1 JSIGN,SIG1,SIG2,CONV)
  172. C
  173. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  174. C % %
  175. C % REACTUALISATION DES SEUILS %
  176. C % %
  177. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  178. c
  179. IF (CONV) THEN
  180. IF (Y1 .GT.YLIM1)THEN
  181. VARF(1)=Y1
  182. VARF(3)=1.D0-1.D0/Z1
  183. ELSE
  184. VARF(1)=VAR0(1)
  185. VARF(3)=VAR0(3)
  186. END IF
  187. IF (Y2 .GT. YLIM2)THEN
  188. VARF(2)=Y2
  189. VARF(4)=1.D0-1.D0/Z2
  190. ELSE
  191. VARF(2)=VAR0(2)
  192. VARF(4)=VAR0(4)
  193. END IF
  194. END IF
  195. END DO
  196. NEPS=NEPS*10
  197. END DO
  198. IF (.NOT. CONV)THEN
  199. print*,'valeur de neps et ieps',neps,ieps
  200. print*, 'NON CONVERGE DANS EULER'
  201. nmatt1 = nmatt
  202. CALL MODCLB(wrk52,nmatt1,EPSP1,EPSP2,YLIM1,YLIM2,Z1,Z2,Y1,Y2,
  203. 1 JSIGN,SIG1,SIG2,CONV)
  204. SEGSUP WRK3
  205. KERRE=56
  206. return
  207. END IF
  208. C
  209. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  210. C % %
  211. C % CALCUL DES CONTRAINTES DANS LE REPERE DE L'ELEMENT %
  212. C % %
  213. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  214. C
  215. SIGNN(1)=SIG1
  216. SIGNN(2)=SIG2
  217. SIGNN(3)=0.D0
  218. CALL PRODT2 (SIGMAN,SIGNN,XROTA,3)
  219. SIGF(1)=SIGMAN(1,1)
  220. SIGF(2)=SIGMAN(2,2)
  221. SIGF(3)=SIGMAN(3,3)
  222. SIGF(4)=SIGMAN(1,2)
  223. C
  224. C
  225. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  226. C % %
  227. C % CALCUL DES DEFORMATIONS ANELASTIQUES %
  228. C % %
  229. C %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  230. C
  231. CALL MATVE1 (DDHOOK,SIGF,NSTRS1,NSTRS1,EPSILI,1)
  232. DO ISTRS=1,NSTRS1
  233. EPINF(ISTRS)=EPSILO(ISTRS)-EPSILI(ISTRS)
  234. END DO
  235. SEGSUP WRK3
  236. RETURN
  237. END
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  

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