Télécharger colbec.eso

Retour à la liste

Numérotation des lignes :

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

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