Télécharger la plaquette Cast3M


fr en

Télécharger beinsa.eso

Retour à la liste

Numérotation des lignes : 

beinsa

  1. C BEINSA    SOURCE    PV        05/04/25    21:15:01     5077
  2. C BEINSA     SOURCE    INSL       24/10/96
  3.       SUBROUTINE BEINSA(SIG0,NSTRSS,DEPST,VAR0,XMAT,IVAL,NMATT,
  4.      1 SIGF,VARF,KERRE,MELE,IFOUR,NVARI,XCAR,NCARR,MFR,EPIN0,EPINF,
  5.      2 DTEMP,XE,NBNN,CMATE,IBB,IGAU,wrk12)
  6. C---------------------------------------------------------------------
  7. C             PLASTICITE MODELE BETON
  8. C
  9. C ENTREES
  10. C     SIG0(NSTRS)   = CONTRAINTES INITIALES
  11. C     NSTRS         = NOMBRE DE CONTRAINTES
  12. C     DEPST(NSTRS)  = INCREMENT DE DEFORMATION TOTALES
  13. C     VAR0(NVARI)   =  VARIABLES INTERNES DEBUT
  14. C     VAR0(  1  )   = EPPE    :Deformation plastique equivalente
  15. C     VAR0(  2  )   = ALPT1   :Coef. multi. Defo. elas. au pic dir. 1
  16. C     VAR0(  3  )   = ALPT2   :Coef. multi. Defo. elas. au pic dir. 2
  17. C     VAR0(  4  )   = IPLA1   :Indicateur plas dir. 1
  18. C     VAR0(  5  )   = IFISU2  :Indicateur de fissure actuelle  (0 ou 1)
  19. C     VAR0(  6  )   = DTR1    :Resistance en traction dir. 1
  20. C     VAR0(  7  )   = DTR2    :Resistance en traction dir. 2
  21. C     VAR0(  8  )   = IFIS0   :Indicateur de fissuration etat precedent
  22. C     VAR0(  9  )   = JFISU   :Indicateur de fissure  (0 ou 1)
  23. C     VAR0( 10  )   = ANGL    :Angle de fissuration
  24. C     VAR0( 11  )   = IPLA    :Indicateur d'etat en bicompres. (0 2 4)
  25. C     VAR0( 12  )   = EQSTRC1 :Contrainte  equiv. de compres. dir. 1
  26. C     VAR0( 13  )   = EPSEQC1 :Deformation equiv. de compres. dir. 1
  27. C     VAR0( 14  )   = EQSTRT1 :Contrainte  equiv. de traction dir. 1
  28. C     VAR0( 15  )   = EPSEQT1 :Deformation equiv. de traction dir. 1
  29. C     VAR0( 16  )   = EQSTRC2 :Contrainte  equiv. de compres. dir. 2
  30. C     VAR0( 17  )   = EPSEQC2 :Deformation equiv. de compres. dir. 2
  31. C     VAR0( 18  )   = EQSTRT2 :Contrainte  equiv. de traction dir. 2
  32. C     VAR0( 19  )   = EPSEQT2 :Deformation equiv. de traction dir. 2
  33. C     VAR0( 20  )   = IDECH1  :Indicateur de decharge beton intergre
  34. C     VAR0( 21  )   = OUV1    :Indicateur ouverture fissure 1
  35. C     VAR0( 22  )   = EDC1    :Module (d'elasti.) initiale equiva. dir.1
  36. C     VAR0( 23  )   = EDC2    :Module (d'elasti.) initiale equiva. dir.2
  37. C     VAR0( 24  )   = ETS1    :Pente courbe de traction (Post-pic)   dir.1
  38. C     VAR0( 25  )   = ETS2    :Pente courbe de traction (Post-pic)   dir.2
  39. C     VAR0( 26  )   = EPF01   :Deformation au pic courbe de traction dir.1
  40. C     VAR0( 27  )   = EPF02   :Deformation au pic courbe de traction dir.2
  41. C     VAR0( 28  )   = OUV2    :Indicateur ouverture fissure 2
  42. C     VAR0( 29  )   = DEFR1   :defo residuelle en compres.1 (beton fissure)
  43. C     VAR0( 30  )   = DEFR2   :defo residuelle en compres.2 (beton fissure)
  44. C     VAR0( 31  )   = IPAS0   :Initial. variables internes 0=oui/1=non
  45. C     VAR0( 32  )   = EPSR(1) :Deform. initi. EpsX (repere local elem.)
  46. C     VAR0( 33  )   = EPSR(2) :Deform. initi. EpsY      ...
  47. C     VAR0( 34  )   = EPSR(3) :Deform. initi. GamXY     ...
  48. C     VAR0( 35  )   = EPSR(4) :Deform. initi. GamXZ     ...
  49. C     VAR0( 36  )   = EPSR(5) :Deform. initi. GamYZ     ...
  50. C     VAR0( 37  )   = DEP(1,1):Terme de matrice elastoplastique tangente
  51. C     VAR0( 38  )   = DEP(1,2):
  52. C     VAR0( 39  )   = DEP(1,3):
  53. C     VAR0( 40  )   = DEP(1,4):
  54. C     VAR0( 41  )   = DEP(1,5):
  55. C     VAR0( 42  )   = DEP(2,1):
  56. C     VAR0( 43  )   = DEP(2,2):
  57. C     VAR0( 44  )   = DEP(2,3):
  58. C     VAR0( 45  )   = DEP(2,4):
  59. C     VAR0( 46  )   = DEP(2,5):
  60. C     VAR0( 47  )   = DEP(3,1):
  61. C     VAR0( 48  )   = DEP(3,2):
  62. C     VAR0( 49  )   = DEP(3,3):
  63. C     VAR0( 50  )   = DEP(3,4):
  64. C     VAR0( 51  )   = DEP(3,5):
  65. C     VAR0( 52  )   = DEP(4,1):
  66. C     VAR0( 53  )   = DEP(4,2):
  67. C     VAR0( 54  )   = DEP(4,3):
  68. C     VAR0( 55  )   = DEP(4,4):
  69. C     VAR0( 56  )   = DEP(4,5):
  70. C     VAR0( 57  )   = DEP(5,1):
  71. C     VAR0( 58  )   = DEP(5,2):
  72. C     VAR0( 59  )   = DEP(5,3):
  73. C     VAR0( 60  )   = DEP(5,4):
  74. C     VAR0( 61  )   = DEP(5,5):
  75. C     XMAT(NMATT)  =  COMPOSANTES DE MATERIAU
  76. C     IVAL(NMATT)  =  INDICE DES COMPOSANTES DE MATERIAU
  77. C     NMATT        = NOMBRE DE COMPOSANTES DE MATERIAU
  78. C  SORTIES
  79. C     SIGF(NSTRS)   = CONTRAINTES FINALES
  80. C     VARF(NVARI)   = VARIABLES INTERNES FINALES
  81. C     DEFP(NSTRS)   = DEFORMATIONS PLASTIQUES
  82. C     KERRE         = 0    TOUT OK
  83. C---------------------------------------------------------------------
  84. C      VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU  ET NECOU
  85. C
  86. C  IFOUR    INDICE DU TYPE DE PROBLEME
  87. C            -2  CONTRAINTES PLANES
  88. C            -1  DEFORMATIONS PLANES
  89. C             0  AXISYMETRIQUE
  90. C             1  SERIE DE FOURIER
  91. C             2  TRIDIMENSIONNEL
  92. C---------------------------------------------------------------------
  93. C        COMPOSANTES DE MATERIAU
  94. C=====================================================================
  95. C  EX   : Module d'Young
  96. C  XNU  : Coeficient de Poisson
  97. C  RHO  : Masse volumique (Facultatif)
  98. C  ALPH : Coeficient de dilation thermique (Facultatif)
  99. C  ALPHA: Resis. tract. simple / resis. compr. simple
  100. C  RB   : Resis. compr. simple
  101. C  EMAX : Deformation ultime en compression (sigma = 0)
  102. C  EPUT : Deformation ultime en traction (sigma = 0)
  103. C  FTC  : Facteur de transfert de cisaillement
  104. C
  105. C  ICAL : Choix du comportement Pre-pic
  106. C  ICAL    0 = Comportement elastoplastique ecrouissable
  107. C          1 = Comportement elastique
  108. C
  109. C  ICU  : Choix de la courbe de compression (Calibrage)
  110. C          1 =  Parabole pre-pic et post-pic
  111. C          2 =  Parabole pre-pic et lineaire post-pic
  112. C          3 =  courbe INSA pre-pic et lineaire post-pic
  113. C          5 =  courbe de Krishnan pre-pic et  post-pic
  114. C          4 =  courbe de Krishnan pre-pic et lineaire post-pic
  115. C
  116. C  ILOI : Choix de la loi d'ecoulement
  117. C          0 = Plasticite Associee
  118. C          1 = Plasticite Non Associee
  119. C
  120. C  IOUV : Nombre d'iteration interne
  121. C
  122. C  IMOD : Choix du modele Beton
  123. C          0 = Modele BETON_INSA ( critere d'OTTOSEN)
  124. C          1 = Modele BETON_INSA ( critere de NADAI partout)
  125. C          2 = Modele BETON_INSA ( NADAI + Sigma Max en trac.)
  126. C
  127. C---------------------------------------------------------------------
  128. C
  129.       IMPLICIT INTEGER(I-N)
  130.       IMPLICIT REAL*8(A-H,O-Z)
  131.       CHARACTER*8 CMATE
  132.       DIMENSION SIG0(NSTRSS),DEPST(NSTRSS),VAR0(NVARI),XMAT(NMATT)
  133.       DIMENSION IVAL(*),SIGF(NSTRSS),VARF(NVARI),XE(3,NBNN)
  134.       DIMENSION EPIN0(NSTRSS),EPINF(NSTRSS),DEPIN(6),EPFLU(6)
  135.       DIMENSION STRN(6),SIGM(6),D(6,6),D1(6,6),D2(6,6),XCAR(NCARR)
  136.       DIMENSION EPSR(6),SIGR(6),VART(100),VV(36)
  137.  
  138.  
  139. C
  140.        segment wrk12
  141.        real*8 AA,BB,DK1,DK2,RB,ALPHA,EX,XNU,EMAX
  142.        real*8 EPUT,FTC,EPO,EPO1,ENGF,RMOY,PHIF,TEMP0
  143.        real*8 DTEMP1,TEMP1,POAR,SCT,TETA,DTR1,DTR2,EDC1
  144.        real*8 EDC2,ETS1,ETS2,EDT1,EDT2,OUV1,OUV2,TANG1
  145.        real*8 TANG2,DEFR1,DEFR2,EPSC1,EPSC2,EPST1,EPST2,EQSTR1
  146.        real*8 EQSTR2,EPSEQ1,EPSEQ2,EQSTR3,EPSEQ3,EPST3,EPSC3,DEFR3
  147.        real*8 RTM3,EDC3,ETS3,EDT3,OUV3,TANG3
  148.        integer ICU,ILOI,IOUV,ICAL,IFLU,IPLA2,IPLA1,IFISU2
  149.        integer IFISU1,JFISU,JFISU2,IPLA3,IFISU3,JFISU3,IBB1,IGAU1
  150.        endsegment
  151. c                     1  2  3   4   5   6    7  8   9    10   11  12
  152. c      COMMON /CINSA/ AA,BB,DK1,DK2,RB,ALPHA,EX,XNU,EMAX,EPUT,FTC,EPO,
  153. c             13    14  15   16    17    18    19    20   21   22
  154. c     1       EPO1,ENGF,RMOY,PHIF,TEMP0,DTEMP1,TEMP1,POAR,SCT,TETA,
  155. c             23   24   25   26   27   28   29   30   31   32   33
  156. c     2       DTR1,DTR2,EDC1,EDC2,ETS1,ETS2,EDT1,EDT2,OUV1,OUV2,TANG1,
  157. c             34     35    36   37    38    39     40    41     42
  158. c     3       TANG2,DEFR1,DEFR2,EPSC1,EPSC2,EPST1,EPST2,EQSTR1,EQSTR2,
  159. C              43     44     45    46     47    48     49    50   51
  160. c     4       EPSEQ1,EPSEQ2,EQSTR3,EPSEQ3,EPST3,EPSC3,DEFR3,RTM3,EDC3,
  161. c             52   53   54    55
  162. c     5       ETS3,EDT3,OUV3,TANG3,
  163. C             1    2   3     4    5    6     7     8      9
  164. c     6       ICU,ILOI,IOUV,ICAL,IFLU,IPLA2,IPLA1,IFISU2,IFISU1,
  165. C             10     11     12     13    14    15    16
  166. c    7       JFISU,JFISU2,IPLA3,IFISU3,JFISU3,IBB1,IGAU1
  167. *
  168.       IRTD =1
  169.       IBB1 =IBB
  170.       IGAU1=IGAU
  171. C
  172. C TEST DE CONSISTANCE DES DONNEES
  173. C
  174.       EX   = XMAT( 1)
  175.       XNU  = XMAT( 2)
  176.       RHO  = XMAT( 3)
  177.       ALPH = XMAT( 4)
  178.       IKA=4
  179.       IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND.IFOUR.EQ.-2) IKA=5
  180. *
  181.       IFOU=IFOUR
  182.       IF(MELE.EQ.28.OR.MFR.EQ.3) IFOU=-2
  183. C
  184.       ALPHA=0.1D0
  185.       RB   =EX*1.D-3
  186.       EMAX =10.D0*RB/EX
  187.       EPUT =3.D0*RB*ALPHA/EX
  188.       FTC  =0.1D0
  189.       ICU  =3
  190.       ILOI =0
  191.       IOUV =2
  192.       IFLU =0
  193.       PHIF =2.D0
  194.       TEMP1=0.D0
  195.       POAR =0.D0
  196.       ENGF =0.D0
  197. *
  198.       ICAL  = 0
  199.       DTEMP1=DTEMP
  200. C
  201.       AA=0.D0
  202.       BB=0.D0
  203.       DK1=0.D0
  204.       DK2=0.D0
  205. *
  206.       ALPHA= XMAT(IKA+1)
  207.       RB   = XMAT(IKA+2)
  208.       EMAX = XMAT(IKA+3)
  209.       EPUT = XMAT(IKA+4)
  210.       FTC  = XMAT(IKA+5)
  211.       ICAL = INT(REAL(XMAT(IKA+6)))
  212.       ENGF = XMAT(IKA+7)
  213. *
  214.       IFLU = INT(REAL(XMAT(IKA+8)))
  215.       RMOY = XMAT(IKA+9)
  216.       PHIF = XMAT(IKA+10)
  217.       TEMP1= XMAT(IKA+11)
  218.       POAR = XMAT(IKA+12)
  219.       EPO  = 2.4D0*RB/EX
  220. C
  221. C       VALEURS PAR DEFAUT
  222. C
  223.       IF(MELE.EQ.4.OR.MELE.EQ.8.OR.MELE.EQ.6.OR.MELE.EQ.10) THEN
  224.          IF(IVAL(IKA+7).EQ.0.AND.IVAL(IKA+4).EQ.0) ENGF =.15D0
  225.       ELSE
  226.          ENGF =.0D0
  227.       ENDIF
  228. *
  229.       IF(IVAL(IKA+8).EQ.0) IFLU=0
  230.       IF(IVAL(IKA+10).EQ.0) PHIF=2.D0
  231.       IF(IVAL(IKA+11).EQ.0) TEMP1=0.D0
  232.       IF(IVAL(IKA+12).EQ.0) POAR =.0D0
  233. C---------------------------------------------------------------------
  234.       IECR=0
  235.       IF(IECR.EQ.10) THEN
  236. *      WRITE(*,*) ' ** XMAT='
  237. *      WRITE(*,1991) (XMAT(I),I=1,NMATT)
  238. *      WRITE(*,*) ' ** IVAL='
  239. *      WRITE(*,1992) (IVAL(I),I=1,NMATT)
  240.       WRITE(*,404) EX,XNU,RHO,ALPH,ALPHA,RB,EMAX,EPUT,FTC,ICAL,ENGF
  241. 404   FORMAT('EX=',E9.3,' XNU=',E9.3,' RHO=',E9.3,' ALPH=',E9.3,
  242.      *' ALPHA=',E9.3,/,' RB=',E9.3,' EMAX=',E9.3,' EPUT=',E9.3,
  243.      &' FTC=',E9.3,' ICAL=',I1,' ENGF=',E9.3)
  244.       WRITE(*,405) NSTRSS,NMATT,MELE,NVARI,NCARR,MFR,NBNN
  245. 405   FORMAT(' NSTRSS=',I1,' NMATT=',I2,' MELE=',I2,/,' NVARI=',I2,
  246.      *' NCARR=',I2,' MFR=',I2,' NBNN=',I2)
  247.       ENDIF
  248. C---------------------------------------------------------------------
  249. C
  250.       IF(IFLU.EQ.20.OR.IFLU.EQ.21) THEN
  251.           IF(IVAL(IKA+10).EQ.0) PHIF= 60.D0
  252.       ENDIF
  253. *
  254.       IF(IFLU.EQ.30.OR.IFLU.EQ.31) THEN
  255.           IF(IVAL(IKA+9).EQ.0) RMOY= 0.159D0
  256.           IF(IVAL(IKA+10).EQ.0) PHIF=-0.88D0
  257.       ENDIF
  258. *
  259.       IF((IFLU.EQ.10.OR.IFLU.EQ.11.OR.IFLU.EQ.20.OR.IFLU.EQ.21).AND.
  260.      &   IVAL(IKA+9).EQ.0) THEN
  261.        WRITE(*,*)'!! ATTENTION POUR FAIRE DU FLUAGE BETON (LOI BPEL91)'
  262.        WRITE(*,*)' IL FAUT DONNER LA VALEUR DU RAYON MOYEN (RMOY) '
  263.        STOP
  264.       ENDIF
  265. C---------------------------------------------------------------------
  266. C
  267. C     CORRESPONDANCE DES VARIABLES POUR BETDJE
  268. C
  269.       CALL ZERO(D,6,6)
  270.       CALL ZERO(D1,6,6)
  271.       CALL ZERO(EPSR,6,1)
  272.       CALL ZERO(DEPIN,6,1)
  273.       CALL ZERO(EPFLU,6,1)
  274.       CALL ZERO(VART,100,1)
  275.       CALL ZERO(SIGR,6,1)
  276.       CALL ZERO(SIGM,6,1)
  277.       CALL ZERO(STRN,6,1)
  278. C
  279.       EPAIST=1.D0
  280.       NSTRS=NSTRSS
  281. *
  282.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.6) THEN
  283.         NSTRS=3
  284.         EPAIST=XCAR(1)
  285.       ENDIF
  286. *
  287.       DO 1 I=1,NSTRS
  288.         SIGR(I)=SIG0(I)/EPAIST
  289.         STRN(I)=DEPST(I)
  290.         EPFLU(I)=EPIN0(I)
  291.     1 CONTINUE
  292. *
  293.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
  294.         NSTRS=3
  295.         SIGR(3) =SIG0(4)
  296.         STRN(3) =DEPST(4)
  297.       ENDIF
  298. C---------------------------------------------------------------------
  299.       IF(IECR.EQ.10) THEN
  300.         WRITE(*,*) ' ** SIGR='
  301.         WRITE(*,1991) (SIGR(I),I=1,NSTRS)
  302.         WRITE(*,*) ' ** STRN='
  303.         WRITE(*,1991) (STRN(I),I=1,NSTRS)
  304.       ENDIF
  305. C------------------------------------------------------
  306. C   INITIALISATION DES VARIABLES INTERNES
  307. C------------------------------------------------------
  308.         DO  II=1,NVARI
  309.           VART(II)=VAR0(II)
  310.         END DO
  311. C
  312.         EPSPL  = VART( 1)
  313.         TANG1  = VART( 2)
  314.         TANG2  = VART( 3)
  315.         IPLA1  = INT(REAL(VART( 4)))
  316.         IFISU2 = INT(REAL(VART( 5)))
  317.         DTR1   = VART( 6)
  318.         DTR2   = VART( 7)
  319.         IFISU1 = INT(REAL(VART( 8)))
  320.         JFISU  = INT(REAL(VART( 9)))
  321.         ANGL   = VART(10)
  322.         JFISU2 = INT(REAL(VART(11)))
  323.         EQSTR1 = VART(12)
  324.         EPSEQ1 = VART(13)
  325.         EPST1  = VART(14)
  326.         EPST2  = VART(15)
  327.         EQSTR2 = VART(16)
  328.         EPSEQ2 = VART(17)
  329.         EPSC1  = VART(18)
  330.         EPSC2  = VART(19)
  331.         IPLA2  = INT(REAL(VART(20)))
  332.         OUV1   = VART(21)
  333.         EDC1   = VART(22)
  334.         EDC2   = VART(23)
  335.         ETS1   = VART(24)
  336.         ETS2   = VART(25)
  337.         EDT1   = VART(26)
  338.         EDT2   = VART(27)
  339.         OUV2   = VART(28)
  340.         DEFR1  = VART(29)
  341.         DEFR2  = VART(30)
  342.         IPAS0  = INT(REAL(VART(31)))
  343.         TEMP0  = VART(32)
  344. C
  345.         KD = 32
  346.         IF(NSTRS.EQ.4.AND.(IFOU.EQ.-1.OR.IFOU.EQ.0)) THEN
  347.           EQSTR3 = VART(KD+1)
  348.           EPSEQ3 = VART(KD+2)
  349.           EPST3  = VART(KD+3)
  350.           EPSC3  = VART(KD+4)
  351.           DEFR3  = VART(KD+5)
  352.           RTM3   = VART(KD+6)
  353.           EDC3   = VART(KD+7)
  354.           ETS3   = VART(KD+8)
  355.           EDT3   = VART(KD+9)
  356.           TANG3  = VART(KD+10)
  357.           OUV3   = VART(KD+11)
  358.           IPLA3  = INT(REAL(VART(KD+12)))
  359.           IFISU3 = INT(REAL(VART(KD+13)))
  360.           JFISU3 = INT(REAL(VART(KD+14)))
  361.           KD = 46
  362.         ENDIF
  363. C
  364.         EPSR(1)= VART(KD+1)
  365.         EPSR(2)= VART(KD+2)
  366.         EPSR(3)= VART(KD+3)
  367. C  -------------------------------------
  368.       IF(IFLU.NE.0.AND.(TEMP0+DTEMP-TEMP1).LE.0.D0) THEN
  369.         WRITE(*,*)'!! ATTENTION DANS LE MODELE DE FLUAGE BETON '
  370.         TEMPF = TEMP0+DTEMP-TEMP1
  371.         WRITE(*,*)' (LOI BPEL91) TEMPS NEGATIF TEMPF=',TEMPF
  372.         WRITE(*,*)' TEMP0=',TEMP0,' TEMP1=',TEMP1,' DTEMP=',DTEMP
  373.         STOP
  374.       ENDIF
  375. C  -------------------------------------
  376.       TETA=ANGL
  377.       IF((NSTRS.EQ.4.OR.NSTRS.EQ.6).AND.IFOU.EQ.-2) THEN
  378.          NSTRS=3
  379.          SIGR(3)=SIG0(4)
  380.          STRN(3)=DEPST(4)
  381.          STRN(4)=DEPST(3)
  382.          SIGR(4)=0.D0
  383.          EPFLU(3)=EPIN0(4)
  384.          EPFLU(4)=EPIN0(3)
  385.       ENDIF
  386. C
  387.       GO TO (10,10,30,40,50,60),NSTRS
  388.    10 CONTINUE
  389.       KERRE=437
  390.       WRITE(*,*) '!! ATTENTION DANS BETON  NSTRS=',NSTRS
  391.       STOP
  392. c+mdj
  393. C      GO TO 1000
  394. c+mdj
  395. C
  396.    60 CONTINUE
  397. C
  398. C     COMPORTEMENT TRIDIMENSIONNEL
  399. C
  400.         EPSR(4)= VART(KD+4)
  401.         EPSR(5)= VART(KD+5)
  402.         EPSR(6)= VART(KD+6)
  403.         GO TO 30
  404. C
  405.    50 CONTINUE
  406. C
  407. C     CONTRAINTES OU DEFORMATIONS PLANES
  408. C       AVEC TOUS LES CISAILLEMENTS
  409. C
  410.         EPSR(4)= VART(KD+4)
  411.         EPSR(5)= VART(KD+5)
  412.         GO TO 30
  413. C
  414.    40 CONTINUE
  415. C
  416. C     COMPORTEMENT AXISYMETRIQUE
  417. C
  418.         EPSR(4)= VART(KD+4)
  419. C
  420.    30 CONTINUE
  421. C
  422. C---------------------------------------------------------------------
  423.        IF(IFLU.EQ.0)  CALL ZERO(EPFLU,6,1)
  424.        IF(ICAL.EQ.1) EPO=RB/EX
  425.        EPO1=EPO
  426.        IF(ICAL.NE.1.AND.ICU.EQ.3) EPO1=0.85D0*EPO
  427. *     *-----------------------*
  428. *     *   MODELE D' OTTOSEN   *
  429. *     *-----------------------*
  430. *
  431.          CALL PAMAOT(ALPHA,AA,BB,DK1,DK2)
  432.          CALL BETDJE(EPSR,SIGR,SIGM,NSTRS,D,D1,IFOU,STRN,EPSPL,
  433.      A   DEPIN,EPFLU,XE,NBNN,MELE,wrk12)
  434. C---------------------------------------------------------------------
  435.       NSTRS1=NSTRSS
  436.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.GT.3) NSTRS1=3
  437. C
  438.       IF((NSTRS1.EQ.4.OR.NSTRS1.EQ.6).AND.IFOU.EQ.-2) THEN
  439.          SIGR(4)=SIGR(3)
  440.          SIGR(3)=0.D0
  441. *
  442.          DEPFZ   =DEPIN(4)
  443.          DEPIN(4)=DEPIN(3)
  444.          DEPIN(3)=DEPFZ
  445.       ENDIF
  446. C
  447.       DO 2 I=1,NSTRS1
  448.         SIGF(I) =SIGR(I)*EPAIST
  449. *       EPINF(I)=EPIN0(I)+DEPIN(I)
  450.         EPINF(I)=DEPIN(I)
  451.     2 CONTINUE
  452. *
  453.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.6) THEN
  454.         CALL ZERO(D1,6,6)
  455.         CALL ZERO(D2,6,6)
  456.         CALL ZERO(SIGM,6,1)
  457. *
  458.         CALL DOHCOM(XMAT,NMATT,CMATE,IFOUR,NSTRSS,D1,IRTD)
  459.         CALL HOOKMU(EPAIST,0.D0,NSTRSS,D1,D2)
  460.         CALL BST(D2,DEPST,NSTRSS,NSTRSS,SIGM)
  461. *
  462.         DO I = 4,NSTRSS
  463.            SIGF(I)=SIG0(I)+SIGM(I)
  464. *          EPINF(I)=EPIN0(I)+DEPIN(I)
  465.            EPINF(I)=DEPIN(I)
  466.         END DO
  467.       ENDIF
  468. *
  469.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
  470.          SIGF(3) =0.D0
  471.          SIGF(4) =SIGR(3)*EPAIST
  472.          DEPFZ   =DEPIN(4)
  473.          DEPIN(4)=DEPIN(3)
  474.          DEPIN(3)=DEPFZ
  475.       ENDIF
  476. C---------------------------------------------------------------------
  477.       IF(IECR.EQ.20) THEN
  478.       WRITE(*,*) ' ** IFLU =',IFLU,' RMOY=',RMOY,' PHIF=',PHIF,
  479.      & 'TEMP0= ',TEMP0,' DTEMP =',DTEMP
  480.       WRITE(*,*) ' **  DEPIN  *** '
  481.       WRITE(*,1991) (DEPIN(IC),IC=1,NSTRSS)
  482.       WRITE(*,*) ' ** EPIN0  / EPINF '
  483.       WRITE(*,1991) (EPIN0(IC),IC=1,NSTRSS),(EPINF(IC),IC=1,NSTRSS)
  484.       ENDIF
  485. C-----------------------------------------------------------------------
  486. C
  487.       VART( 1)=EPSPL
  488.       VART( 2)=TANG1
  489.       VART( 3)=TANG2
  490. c+mdj
  491.       VART( 4)=IPLA1
  492.       VART( 5)=IFISU2
  493. c+mdj
  494.       VART( 6)=DTR1
  495.       VART( 7)=DTR2
  496. c+mdj
  497.       VART( 8)=IFISU1
  498.       VART( 9)=JFISU
  499. c+mdj
  500.       VART(10)=TETA
  501. c+mdj
  502.       VART(11)=JFISU2
  503. c+mdj
  504.       VART(12)=EQSTR1
  505.       VART(13)=EPSEQ1
  506.       VART(14)=EPST1
  507.       VART(15)=EPST2
  508.       VART(16)=EQSTR2
  509.       VART(17)=EPSEQ2
  510.       VART(18)=EPSC1
  511.       VART(19)=EPSC2
  512. c+mdj
  513.       VART(20)=IPLA2
  514. c+mdj
  515.       VART(21)=OUV1
  516.       VART(22)=EDC1
  517.       VART(23)=EDC2
  518.       VART(24)=ETS1
  519.       VART(25)=ETS2
  520.       VART(26)=EDT1
  521.       VART(27)=EDT2
  522.       VART(28)=OUV2
  523.       VART(29)=DEFR1
  524.       VART(30)=DEFR2
  525.       VART(31)=1.D0
  526.       VART(32)=TEMP0
  527. C
  528.       KD = 32
  529.       IF(NSTRSS.EQ.4.AND.(IFOU.EQ.-1.OR.IFOU.EQ.0)) THEN
  530.         VART(KD+1) =EQSTR3
  531.         VART(KD+2) =EPSEQ3
  532.         VART(KD+3) =EPST3
  533.         VART(KD+4) =EPSC3
  534.         VART(KD+5) =DEFR3
  535.         VART(KD+6) =RTM3
  536.         VART(KD+7) =EDC3
  537.         VART(KD+8) =ETS3
  538.         VART(KD+9) =EDT3
  539.         VART(KD+10)=TANG3
  540.         VART(KD+11)=OUV3
  541.         VART(KD+12)=IPLA3
  542.         VART(KD+13)=IFISU3
  543.         VART(KD+14)=JFISU3
  544.         KD = 46
  545.       ENDIF
  546. C
  547.       VART(KD+1)=EPSR(1)
  548.       VART(KD+2)=EPSR(2)
  549.       VART(KD+3)=EPSR(3)
  550.       K1=KD+3
  551. C
  552.       GO TO (11,11,31,41,51,61),NSTRS
  553.    11 CONTINUE
  554.       KERRE=437
  555.       WRITE(*,*) '!! ATTENTION DANS BETON  NSTRS=',NSTRS
  556.       STOP
  557. c+mdj
  558. C      GO TO 1000
  559. c+mdj
  560. C
  561.    61 CONTINUE
  562. C
  563. C     COMPORTEMENT TRIDIMENSIONNEL
  564. C
  565.         VART(KD+4)=EPSR(4)
  566.         VART(KD+5)=EPSR(5)
  567.         VART(KD+6)=EPSR(6)
  568.         K1=KD+6
  569.         GO TO 31
  570. C
  571.    51 CONTINUE
  572. C
  573. C     CONTRAINTES OU DEFORMATIONS PLANES
  574. C       AVEC TOUS LES CISAILLEMENTS
  575. C
  576.         VART(KD+4)=EPSR(4)
  577.         VART(KD+5)=EPSR(5)
  578.         K1=KD+5
  579.         GO TO 31
  580. C
  581.    41 CONTINUE
  582. C
  583. C     COMPORTEMENT AXISYMETRIQUE
  584. C
  585.         VART(KD+4)=EPSR(4)
  586.         K1=KD+4
  587.         GO TO 31
  588.    31 CONTINUE
  589. C
  590. C     CONTRAINTES OU DEFORMATIONS PLANES
  591. C       SANS CISAILLEMENTS TRANSVERSAL
  592. C
  593.       L1=0
  594.       DO 66 I1 = 1,6
  595.       DO 67 J1 = 1,6
  596.          L1=L1+1
  597.          VV(L1) = D(J1,I1)
  598.    67 CONTINUE
  599.    66 CONTINUE
  600. *
  601.       L1=0
  602.       DO 62 I1 = 1,NSTRS
  603.       DO 63 J1 = 1,NSTRS
  604.         K=K1+(I1-1)*NSTRS+J1
  605.         L1=L1+1
  606.         VART(K)=VV(L1)
  607.    63 CONTINUE
  608.    62 CONTINUE
  609. C
  610.       DO 23 II=1,NVARI
  611.         VARF(II)=VART(II)
  612.    23 CONTINUE
  613. C---------------------------------------------------------------------
  614.       IF(IECR.EQ.10) THEN
  615.          WRITE(*,66770) IBB,IGAU
  616. 66770    format(///'SORTIE DE BEINSA : element  ',i6,2x,'point  ',i3//)
  617.          WRITE(*,66771) MELE,IFOU,MFR,NSTRS
  618. 66771    format('0  MELE=',i4,2x,'ifour=',i4,'mfr=',i2,'nstrs=',i2/)
  619.          WRITE(*,66774) (DEPST(I),I=1,NSTRSS)
  620. 66774    format(2x,'  depst '/(6(1x,1pe12.5)))
  621.          WRITE(*,66772) (SIGF(I),I=1,NSTRSS)
  622. 66772    format(2x,'  sigf '/(6(1x,1pe12.5)))
  623.          WRITE(*,66773) (VARF(I),I=1,NVARI)
  624. 66773    format(2x,'  varf '/(6(1x,1pe12.5)))
  625.       ENDIF
  626. C
  627. 1991  FORMAT(18(1X,E12.5))
  628. 1992  FORMAT(18(1X,I6))
  629. C---------------------------------------------------------------------
  630. C
  631.  1000 CONTINUE
  632. C
  633.       RETURN
  634.       END
  635.  
  636.  
  637.  
  638.  
  639.  
  640.  
  641.  
  642.  

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