Télécharger la plaquette Cast3M


fr en

Télécharger biflex.eso

Retour à la liste

Numérotation des lignes : 

biflex

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

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