Télécharger la plaquette Cast3M


fr en

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