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

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