Télécharger la plaquette Cast3M


fr en

Télécharger cbifle.eso

Retour à la liste

Numérotation des lignes : 

cbifle

  1. C CBIFLE    SOURCE    CB215821  24/04/12    21:15:11     11897          
  2.       SUBROUTINE CBIFLE(wrk52,wrk53,wrk54,NSTRS1,NVARI,ICARA)
  3. C BIFLEX    SOURCE    CHAT      00/09/29    21:15:55     3957
  4. c      SUBROUTINE BIFLEX(WRK0,WRK1,NSTRS,NVARI,NMATT,ICARA,CMATE,KERRE)
  5. *--------------------------------------------------------------------*
  6. *     ecoulement du modele a fibre, interface avec ecoul2
  7. *--------------------------------------------------------------------*
  8. C
  9. C
  10. C     variables en entree
  11. C
  12. C     WRK0,KRK1  pointeurs sur des segments de travail
  13. C
  14. C     NSTRS1      nombre de composantes dans les vecteurs des contraintes
  15. C                                        et les vecteurs des deformations
  16. C
  17. C     NVARI      nombre de variables internes (doit etre egal a 2)
  18. C
  19. C     NMATT      nombre de constantes du materiau
  20. C
  21. C
  22. C     variables en sortie
  23. C
  24. C     VARF      variables internes finales dans WRK1
  25. C
  26. C     SIGF      contraintes finales dans WRK1
  27. C
  28. *--------------------------------------------------------------------*
  29. *     Pierre Pegon (ISPRA) Juillet/Aout 1993                         *
  30. *                          Fevrier 1996 nouvel ecou                  *
  31. *--------------------------------------------------------------------*
  32.       IMPLICIT INTEGER(I-N)
  33.       IMPLICIT REAL*8(A-H,O-Z)
  34. *
  35.  
  36. -INC PPARAM
  37. -INC CCOPTIO
  38. -INC SMCHAML
  39. -INC SMELEME
  40. -INC SMCOORD
  41. -INC SMMODEL
  42. -INC SMINTE
  43. -INC CCHAMP
  44. -INC SMLREEL
  45. -INC DECHE
  46. *
  47.       DIMENSION DEPSI(6),SIGMA(6)
  48.       DIMENSION  CRIGI1(12),CMASS1(12)
  49. *
  50.       INTEGER NSTRS1,NVARI
  51.       INTEGER IRTD,ISTRS
  52.       KERRE=0
  53. C
  54. C     recuperation des instants initial et final
  55.       TIME0 = TEMP0
  56.       TIMEF = TEMPF
  57.  
  58. C     calcul ou recuperation de la matrice elastique
  59. C
  60.       IPMODL = NINT(XMAT(1))
  61.       IPCAR  = NINT(XMAT(2))
  62.       MLREEL = NINT(XMAT(3))
  63.       IF(MLREEL.EQ.0)THEN
  64.         CALL FRIGIE(IPMODL,IPCAR,CRIGI1,CMASS1)
  65.         CALL DOHTIF(CRIGI1,CMATE,IFOUR,NSTRS1,DDHOOK,IRTD)
  66.       ELSE
  67.         SEGACT, MLREEL
  68.         IF (PROG(/1).NE.NSTRS1**2)THEN
  69.           IRTD = 0
  70.         ELSE
  71.           IRTD = 1
  72.           CALL DOHOOO(PROG,NSTRS1,DDHOOK)
  73.         ENDIF
  74.         SEGDES, MLREEL
  75.       ENDIF
  76. C
  77.       IF ( IRTD .EQ. 1) THEN
  78. C
  79. C     calcul de l incrément de contrainte (pour les defo plastique)
  80. C
  81.          CALL MATVE1 (DDHOOK,DEPST,NSTRS1,NSTRS1,DSIGT,2)
  82. C
  83. C     inversion de la matrice DDHOOK
  84. C
  85.          PREC=1.D-08
  86.          CALL INVALM(DDHOOK,NSTRS1,NSTRS1,IRTD,PREC)
  87.          IF (IRTD.EQ.0)THEN
  88. C
  89. C      chargement de l'increment de deformation
  90. C
  91.             DO I=1,NSTRS1
  92.               DEPSI(I)=DEPST(I)
  93.             ENDDO
  94. C
  95. C      recuperation du modele, des materiaux, des contraintes
  96. C      et des varint de la section
  97. C
  98.             IPSIG0=NINT(VAR0(1))
  99.             IPVAR0=NINT(VAR0(2))
  100.             IF(IPSIG0.EQ.0) CALL ZEROP(IPMODL,'CONTRAIN',IPSIG0)
  101.             IF(IPVAR0.EQ.0) CALL ZEROP(IPMODL,'VARINTER',IPVAR0)
  102. C
  103. C      ecoulement de la section
  104. C
  105.             CALL FCOUL1(DEPSI,IPMODL,IPSIG0,IPVAR0,IPCAR,TIME0,TIMEF,
  106.      $                  SIGMA,IPSIGF,IPVARF,IRETO,NSTRS1)
  107.             IF(IRETO.EQ.0) THEN
  108.               print*,'erreur dans fcoul1'
  109.               KERRE=56
  110.             END IF
  111. C
  112. C      calcul de l'increment de deformation plastique
  113. C
  114.             DO I=1,NSTRS1
  115.               SIGF (I)=SIGMA(I)
  116.               DSIGT(I)=DSIGT(I) - (SIGF(I) - SIG0(I))
  117.             END DO
  118. C
  119.             CALL MATVE1 (DDHOOK,DSIGT,NSTRS1,NSTRS1,DEFP,1)
  120.             DO I=1,NSTRS1
  121.             END DO
  122. C
  123. C      sortie des variables internes
  124. C
  125.             VARF(1)=DBLE(IPSIGF)
  126.             VARF(2)=DBLE(IPVARF)
  127.          ELSE
  128.              print*,'erreur dans invalm'
  129.              KERRE=56
  130.          END IF
  131.       ELSE
  132.          print*,'erreur dans dohtif'
  133.          KERRE=56
  134.       END IF
  135. C
  136.  998  RETURN
  137.       END
  138.  
  139.  
  140.  
  141.  

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