Télécharger biflex.eso

Retour à la liste

Numérotation des lignes :

  1. C BIFLEX SOURCE BP208322 15/06/22 21:15:22 8543
  2. SUBROUTINE BIFLEX(WRK0,WRK1,NSTRS,NVARI,NMATT,ICARA,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.  
  34. -INC PPARAM
  35. -INC CCOPTIO
  36. -INC SMCHAML
  37. -INC SMELEME
  38. -INC SMCOORD
  39. -INC SMMODEL
  40. -INC SMINTE
  41. -INC CCHAMP
  42. -INC SMLREEL
  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**2)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. CALL FCOUL1(DEPSI,IPMODL,IPSIG0,IPVAR0,IPCAR,
  110. $ SIGMA,IPSIGF,IPVARF,IRETO,NSTRS)
  111. IF(IRETO.EQ.0) THEN
  112. print*,'erreur dans fcoul1'
  113. KERRE=56
  114. END IF
  115. C
  116. C calcul de l'increment de deformation plastique
  117. C
  118. DO I=1,NSTRS
  119. SIGF (I)=SIGMA(I)
  120. DSIGT(I)=DSIGT(I) - (SIGF(I) - SIG0(I))
  121. END DO
  122. C
  123. CALL MATVE1 (DDHOOK,DSIGT,NSTRS,NSTRS,DEFP,1)
  124. DO I=1,NSTRS
  125. END DO
  126. C
  127. C sortie des variables internes
  128. C
  129. VARF(1)=DBLE(IPSIGF)
  130. VARF(2)=DBLE(IPVARF)
  131. ELSE
  132. print*,'erreur dans invalm'
  133. KERRE=56
  134. END IF
  135. ELSE
  136. print*,'erreur dans dohtif'
  137. KERRE=56
  138. END IF
  139. C
  140. 998 RETURN
  141. END
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

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