Télécharger cbifle.eso

Retour à la liste

Numérotation des lignes :

cbifle
  1. C CBIFLE SOURCE FD218221 19/12/04 21:15:07 10398
  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.  

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