Télécharger la plaquette Cast3M


fr en

Télécharger curgcs.eso

Retour à la liste

Numérotation des lignes : 

curgcs

  1. C CURGCS    SOURCE    PV        22/04/22    21:15:08     11344          
  2.       SUBROUTINE CURGCS(WRK52,WRK53,WRK54,WRK4,NSTRSS,
  3.      1  IFOUR,DDT,IB,IGAU,HCAR,IVIS,T0,TF)
  4. C---------------------------------------------------------------------
  5. C             PLASTICITE MODELE BETON
  6. C
  7. C ENTREES
  8. C     SIG0(NSTRS)   = CONTRAINTES INITIALES
  9. C     NSTRS         = NOMBRE DE CONTRAINTES
  10. C     DEPST(NSTRS)  = INCREMENT DE DEFORMATIONS TOTALES
  11. C     VAR0(NVARI)   =  VARIABLES INTERNES DEBUT
  12. C     VAR0(  1  )   = IFIS    :Indicateur de fissuration (0 1 2)
  13. C     VAR0(  2  )   = ANGL    :Angle de fissuration
  14. C     VAR0(  3  )   = IPLA    :Indicateur d'etat en bicompres. (0 1 2 3)
  15. C     VAR0(  4  )   = SIG1    :Force d écrouissage de traction
  16. C     VAR0(  5  )   = SIG2    :Force d écrouissage de compression
  17. C     VAR0(  6  )   = EPS1    :Variable d écrouissage de traction
  18. C     VAR0(  7  )   = EPS2    :Variable d écrouissage de compression
  19. C     VAR0(  8  )   = TDEF      :Taux de déformation
  20. C     VAR0(  9  )   = TCON      :Taux de contrainte
  21. C     VAR0( 10  )   = SIGP(1):Contrainte plastique en mode
  22. C                              viscoplastique
  23. C     VAR0( 11  )   = SIGP(2):
  24. C     VAR0( 12  )   = SIGP(3):
  25. C     VAR0( 13  )   = SIGP(4):
  26. C     VAR0( 14  )   = DPSTV1 :Variable d'écrouissage plastique
  27. C                              en mode viscoplastique
  28. C     VAR0( 15  )   = DPSTV2 :
  29. C     VAR0( 16  )   = SIGV1  :Force d ecrouissage de traction
  30. C                             en mode viscoplastique
  31. C     VAR0( 17  )   = SIGV2  :Force d ecrouissage de compression
  32. C                             en mode viscoplastique
  33. C
  34. C     XMAT(NCOMAT)  =  COMPOSANTES DE MATERIAU
  35. C     NCOMAT        = NOMBRE DE COMPOSANTES DE MATERIAU
  36. C  SORTIES
  37. C     SIGF(NSTRS)   = CONTRAINTES FINALES
  38. C     VARF(NVARI)   = VARIABLES INTERNES FINALES
  39. C     KERRE         = 0    TOUT OK
  40. C---------------------------------------------------------------------
  41. C
  42. C  IFOUR    INDICE DU TYPE DE PROBLEME
  43. C            -2  CONTRAINTES PLANES
  44. C            -1  DEFORMATIONS PLANES
  45. C             0  AXISYMETRIQUE
  46. C             1  SERIE DE FOURIER
  47. C             2  TRIDIMENSIONNEL
  48. C---------------------------------------------------------------------
  49. C        COMPOSANTES DE MATERIAU
  50. C=====================================================================
  51. C  YOUN : Module d'Young
  52. C  XNU  : Coeficient de Poisson
  53. C  RHO  : Masse volumique (Facultatif)
  54. C  ALPH : Coeficient de dilation thermique (Facultatif)
  55. C  ALFA : Resis. tract. simple / resis. compr. simple
  56. C  BETA : Resis. compr. biax. / resis. compr. simple (Modele de NADAI)
  57. C  RB   : Resis. compr. simple
  58. C  GFC : Energie de rupt en compression
  59. C  GFT : Energie de rupt en traction
  60. C  ETA : parametre de viscosite
  61. C  ICT  : Choix de la courbe de traction (Calibrage)
  62. C  ICC : Choix de la courbe de compression (Calibrage)
  63. C
  64. C  IVIS : Modele visqueux ou non
  65. C            0 : non visqueux
  66. C            1 : viscoplastique
  67. C            2 : viscoelastoplastique
  68. C            3 : elastoplastique endommageable
  69. C            4 : elasto viscoplastique endommageable
  70. C
  71. C  IMOD : Choix du modele Beton
  72. C          1 = Modele BETON_INSA ( Von MISES + Sigma Max en trac.)
  73. C              Rheologie 2D
  74. C          2 = Modele BETON_INSA ( Von MISES + Sigma Max en trac.)
  75. C              Rheologie 3D
  76. C          3 = Modele BETON_INSA ( Druck Prager + Sigma Max en trac.)
  77. C              Rheologie 2D
  78. C          4 = Modele BETON_INSA ( Druck Prager + Sigma Max en trac.)
  79. C              Rheologie 3D
  80. C  ITER : nombre d'iterations internes admissibles
  81. C Pour modeles fluage
  82. C TZER = DATE MINI POUR APPLICATION DU PREMIER CHARGEMENT
  83. C PGDUR = CHOIX DE PROGRESSION DES PAS DE TEMPS DE DISCRETISATION DES COURBES DE RELAXATION
  84. C PGTZO = CHOIX DE PROGRESSION DES PAS DE TEMPS DE DISCRETISATION DES COURBES DE RELAXATION
  85. C ITYPE = CHOIX DU CODE DE CALCUL DE FLUAGE CEB-FIP, B3, EUROCODE2, ETC
  86. C TAU2 = SECOND TEMPS DE RELAXATION PERMETTANT LE CALCUL DE LA SUITE DES TEMPS DE RELAXATION
  87. C TAU1 = PREMIER TEMPS DE RELAXATION PERMETTANT LA REPONSE ELASTIQUE
  88. C DATSEC = DATE DE DEBUT DE SECHAGE
  89. C DATCUR = DATE DE FIN DE CURE DU BETON
  90. C NBRC = NOMBRE DE BRANCHE DE MAXWELL
  91. C DATCOU = DATE D APPLICATION DU CHARGEMENT
  92. C IFOR = 1 formulation solide
  93. C      = 2 liquide
  94. C---------------------------------------------------------------------
  95. C
  96.       IMPLICIT INTEGER(I-N)
  97.       IMPLICIT REAL*8(A-H,O-Z)
  98. -INC DECHE
  99.       REAL*8 HCAR
  100.       CHARACTER*40 FMT,TITRE
  101.       DIMENSION SIR(9,4),EPST(4),SIGP(4)
  102.       DIMENSION D(6,6),D1(6,6),STRN(4)
  103.       DIMENSION SIGR(4),VART(200),VV(36),SIGMF(4)
  104.       DIMENSION SIRL(8,4)
  105. C
  106.       SEGMENT WRK4
  107.         REAL*8 XE(3,NBNN)
  108.       ENDSEGMENT
  109.  
  110.       SEGMENT BETJEF
  111.        REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  112.      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  113.        INTEGER ICT,ICC,IMOD,IVISS,ITER,
  114.      &                ISIM,IBB,IGAU1,IZON
  115.       ENDSEGMENT
  116.       SEGMENT VISCO
  117.        REAL*8 DPSTV1,DPSTV2,SIGV1,SIGV2,ENDV
  118.       ENDSEGMENT
  119.       SEGMENT BETFLU
  120.        REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  121.      &         TP0,TZER
  122.        INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  123.       ENDSEGMENT
  124.       SEGMENT NECH0
  125.        REAL*8 DT,DC,ALFG,S0,ENDO
  126.       ENDSEGMENT
  127.       SEGMENT NECH1
  128.        REAL*8 DLMT
  129.       ENDSEGMENT
  130.       SEGMENT NECH2
  131.        REAL*8 ATR,GTR,ALPH0
  132.       ENDSEGMENT
  133.       SEGMENT NECH3
  134.        REAL*8 RBT,ALFAT,YOUNT,GFCT,GFTT,ALPH
  135.       ENDSEGMENT
  136. C
  137. C      NCOMAT = nmatt
  138. C      NBNN=XE(/2)
  139.       NVARI=VAR0(/1)
  140.       SEGINI BETJEF
  141.       SEGINI VISCO
  142.       SEGINI BETFLU
  143.       SEGINI NECH0
  144.       SEGINI NECH1
  145.       SEGINI NECH2
  146.       SEGINI NECH3
  147.       PDT = DDT
  148.       TP0=TEMP0
  149.       TP00=TEMP0
  150.       CAR = HCAR
  151.       CALL ZERO(SIGP,4,1)
  152.       CALL ZERO(SIR,9,4)
  153.       CALL ZERO(SIRL,8,4)
  154. C
  155. C TEST DE CONSISTANCE DES DONNEES
  156. C
  157.       YOUNT = 0.D0
  158.       YOUN = 0.D0
  159.       XNU  = 0.D0
  160.       RHO  = 0.D0
  161.       ALPH = 0.D0
  162.       ALFAT = 0.D0
  163.       ALFA = 0.D0
  164.       BETA = 0.D0
  165.       RBT   = 0.D0
  166.       RB   = 0.D0
  167.       GFCT  = 0.D0
  168.       GFTT  = 0.D0
  169.       GFC  = 0.D0
  170.       GFT  = 0.D0
  171.       ETA  = 0.D0
  172.       ITER = 0
  173.       IMOD = 0
  174.       ICT  = 0
  175.       ICC  = 0
  176.       DC   =0.D0
  177.       DT   =0.D0
  178.       IASO =0
  179.       ICOU =0
  180.       ENDO = 0.D0
  181.       ENDV =0.D0
  182. C
  183.       AA=0.D0
  184.       BB=0.D0
  185.       DK1=0.D0
  186.       DK2=0.D0
  187.       CNC1=0.D0
  188.       CNC1=0.D0
  189.       CNT1=0.D0
  190.       CNT2=0.D0
  191. C
  192.       YOUNT = XMAT( 1)
  193.       XNU  = XMAT( 2)
  194.       RHO  = XMAT( 3)
  195.       ALPH = XMAT( 4)
  196.       ALFAT = XMAT( 5)
  197.       BETA = XMAT( 6)
  198.       RBT   = XMAT( 7)
  199.       GFCT  = XMAT( 8)
  200.       GFTT  = XMAT( 9)
  201.       ITER = INT(REAL(XMAT(10)))
  202.       IMOD = INT(REAL(XMAT(11)))
  203.        IF (IVIS.EQ.0) THEN
  204.         ICT  = INT(REAL(XMAT(12)))
  205.         ICC  = INT(REAL(XMAT(13)))
  206.        ENDIF
  207.        IF (IVIS.EQ.1) THEN
  208.         ICT  = INT(REAL(XMAT(12)))
  209.         ICC  = INT(REAL(XMAT(13)))
  210.         ETA  = XMAT(14)
  211.         ISIM = 1
  212.        ENDIF
  213.        IF (IVIS.EQ.2) THEN
  214.         ICT  = INT(REAL(XMAT(12)))
  215.         ICC  = INT(REAL(XMAT(13)))
  216.         DATCOU  = XMAT(14)
  217.         NBRC = INT(REAL(XMAT(15)))
  218.         IFOR = INT(REAL(XMAT(16)))
  219.         ITYPE = INT(REAL(XMAT(17)))
  220.         DATCUR = XMAT(18)
  221.         DATSEC = XMAT(19)
  222.         E28= YOUN
  223.        ENDIF
  224.        IF (IVIS.EQ.3) THEN
  225.         DC   = XMAT(12)
  226.         DT   = XMAT(13)
  227.         IASO = INT(REAL(XMAT(14)))
  228.         ICOU = INT(REAL(XMAT(15)))
  229.        ENDIF
  230.        IF (IVIS.EQ.4) THEN
  231.         DC   = XMAT(12)
  232.         DT   = XMAT(13)
  233.         IASO = INT(REAL(XMAT(14)))
  234.         ICOU = INT(REAL(XMAT(15)))
  235.         ETA  = XMAT(16)
  236.        ENDIF
  237. C
  238. C---------------------------------------------------------------------
  239. C        INITIALISATION DES VARIABLES DU MODELE D'ENDOMMAGEMENT
  240. C---------------------------------------------------------------------
  241. C      A: CALCUL ASSOCIE OU NON
  242. C
  243.       IF(IASO.EQ.0) THEN
  244.        ALFG = (BETA-1.D0)/(2*BETA-1.D0)
  245.       ELSE
  246.        ALFG = 0.2
  247.       ENDIF
  248. C---------------------------------------------------------------------
  249. C      B: CALCUL COUPLE OU NON
  250. C
  251.       IF(ICOU.EQ.0) THEN
  252.        ATR=0.D0
  253.        GTR=0.D0
  254.       ELSE
  255.        ATR=2.075D0
  256.        GTR=0.235D0
  257.       ENDIF
  258. C---------------------------------------------------------------------
  259. C
  260.       DLMT=0.9999D0
  261. C       S0  = 0.25D0
  262.       IF(IVIS.EQ.4)THEN
  263. C attention en comportement visqueux refermeture pose probleme
  264.         S0  = 1.D0
  265.       ELSE
  266.          S0  = 0.25D0
  267. C        S0  = 1.D0
  268.       ENDIF
  269. C
  270. C---------------------------------------------------------------------
  271. C                    CAS DU COUPLAGE THERMIQUE-MECANIQUE
  272. C---------------------------------------------------------------------
  273. C       IF IASO=0 Utilisation des lois f(T) locales
  274. C       IF IASO=1 Utilisation des lois du fichier de donnée
  275. C
  276.        ITHHER = 1
  277.        IF(ITHHER.NE.0)THEN
  278.          CALL EVOLTH(T0,TF,BETJEF,NECH0,NECH3)
  279.        ELSE
  280. CC        write(*,*)'ITHHER=',ITHHER
  281.         RB   = RBT
  282.         ALFA = ALFAT
  283.         YOUN = YOUNT
  284.         GFC  = GFCT
  285.         GFT  = GFTT
  286.        ENDIF
  287. C
  288. C
  289. C---------------------------------------------------------------------
  290. C
  291. C     CORRESPONDANCE DES VARIABLES POUR URGCST
  292. C
  293.       CALL ZERO(D,6,6)
  294.       CALL ZERO(D1,6,6)
  295.       CALL ZERO(VART,NVARI,1)
  296.       CALL ZERO(SIGR,4,1)
  297. C
  298.       IVISS = IVIS
  299.       IBB = IB
  300.       IGAU1 = IGAU
  301.       EPAIST=1.D0
  302.       NSTRS=NSTRSS
  303.       IFOU=IFOUR
  304. *
  305.       DO 1 I=1,NSTRS
  306.         SIGR(I)=SIG0(I)/EPAIST
  307. *       WRITE(*,*)'SIG0',I,'=',SIG0(I)
  308.         STRN(I)=DEPST(I)
  309. *       WRITE(*,*)'DEPST',I,'=',DEPST(I)
  310.     1 CONTINUE
  311. *
  312.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
  313.         NSTRS=3
  314.         IFOU=-2
  315.         SIGR(3)=SIG0(4)
  316.         STRN(3)=DEPST(4)
  317.         STRN(4)=DEPST(3)
  318.       ENDIF
  319. *
  320.       IF((NSTRS.EQ.4).AND.IFOUR.EQ.-2.AND.
  321.      *         (IMOD.NE.2.AND.IMOD.NE.4)) THEN
  322.          NSTRS=3
  323.          SIGR(3)=SIG0(4)
  324.          STRN(3)=DEPST(4)
  325.          STRN(4)=DEPST(3)
  326.       ENDIF
  327. *
  328.       DO 22 II=1,NVARI
  329.         VART(II)=VAR0(II)
  330.    22 CONTINUE
  331. C------------------------------------------------------
  332. C   INITIALISATION DES VARIABLES INTERNES
  333. C------------------------------------------------------
  334. C       A: MODELE PLASTIQUE
  335. C
  336.         IFISU  = INT(REAL(VART(1)))
  337.         ANGL   = VART(2)
  338.         IPLA   = INT(REAL(VART(3)))
  339.         SIG1= VART(4)
  340.         SIG2= VART(5)
  341.         DPSTF1 = VART(6)
  342.         DPSTF2 = VART(7)
  343. C
  344. C------------------------------------------------------
  345. C       B: MODELE VISCOPLASTIQUE CLASSIQUE
  346. C
  347.         IF (IVIS.EQ.1) THEN
  348.         TDEF    = VART( 8)
  349.         TCON    = VART( 9)
  350.         SIGP(1) = VART(10)
  351.         SIGP(2) = VART(11)
  352.         SIGP(3) = VART(12)
  353.         SIGP(4) = VART(13)
  354.         DPSTV1 = VART(14)
  355.         DPSTV2 = VART(15)
  356.         SIGV1 = VART(16)
  357.         SIGV2 = VART(17)
  358.         ENDIF
  359. C
  360. C------------------------------------------------------
  361. C       C: MODELE VISCO-ELASTO PLASTIQUE
  362. C
  363.         IF (IVIS.EQ.2.AND.IFOR.EQ.1) THEN
  364. C modele visco elasto plastique solide de MAXWELL
  365. C-------------------------------------------------
  366.         MC = NBRC + 1
  367.       DO 101 I1 = 1,MC
  368.       DO 102 J1 = 1,NSTRS
  369.          NV = NSTRS * (I1 - 1) + J1
  370.          SIR(I1,J1) = VART(7 + NV)
  371.   102 CONTINUE
  372.   101 CONTINUE
  373. C
  374.       ELSE IF (IVIS.EQ.2.AND.IFOR.EQ.2) THEN
  375. C modele visco elasto plastique liquide de MAXWELL
  376. C-------------------------------------------------
  377.         MC = NBRC + 1
  378.       DO 131 I2 = 1,MC
  379.       DO 132 J2 = 1,NSTRS
  380.          NV = NSTRS * (I2 - 1) + J2
  381.          SIRL(I2,J2) = VART(7 + NV)
  382.   132 CONTINUE
  383.   131 CONTINUE
  384. C
  385.       ENDIF
  386. C
  387. C----------------------------------------------
  388. C       D: MODELE PLASTIQUE ENDOMMAGEABLE
  389. C
  390.         IF (IVIS.EQ.3) THEN
  391.          ENDO    = VART( 8)
  392.          DECO    = VART( 9)
  393.         ENDIF
  394. C
  395. C----------------------------------------------
  396. C       D: MODELE VISCOPLASTIQUE ENDOMMAGEABLE
  397. C
  398.         IF (IVIS.EQ.4) THEN
  399.         TDEF    = VART( 8)
  400.         TCON    = VART( 9)
  401.         SIGP(1) = VART(10)
  402.         SIGP(2) = VART(11)
  403.         SIGP(3) = VART(12)
  404.         SIGP(4) = VART(13)
  405.         DPSTV1 = VART(14)
  406.         DPSTV2 = VART(15)
  407.         SIGV1 = VART(16)
  408.         SIGV2 = VART(17)
  409.              ENDO  = VART(18)
  410.         ENDV  = VART(19)
  411.         ENDIF
  412. C  -------------------------------------
  413.       TETA=ANGL
  414. C---------------------------------------------------------------------
  415. C
  416.       GO TO (10,10,30,40),NSTRS
  417.    10 CONTINUE
  418.       KERRE=437
  419.       WRITE(*,*) '!! ATTENTION DANS BETON  NSTRS=',NSTRS
  420.       STOP
  421.       GO TO 1000
  422. C
  423.    40 CONTINUE
  424. C
  425.    30 CONTINUE
  426. C
  427. C---------------------------------------------------------------------
  428. C
  429.       IF(IMOD.GE.1.AND.IMOD.LE.4) THEN
  430. *     *-----------------------*
  431. *     *   MODELE URGC ST      *
  432. *     *-----------------------*
  433. *
  434.          CALL BONE(SIGR,SIGMF,STRN,IPLA,IFISU,SIG1,SIG2
  435.      A        ,NSTRS,D,D1,IFOU,SIGP,EPST,SIR,SIRL,
  436.      B         ENDO,ITHHER,T0,TF,BETJEF,VISCO,BETFLU,
  437.      C         NECH0,NECH1,NECH2,NECH3)
  438. C
  439.       ELSE
  440.          WRITE(*,*) '!! ATTENTION CE MODELE N EXISTE PAS (URGCST)'
  441.          STOP
  442.          RETURN
  443.       ENDIF
  444. C---------------------------------------------------------------------
  445. C
  446. C
  447.       IF((NSTRS.EQ.3).AND.IFOU.EQ.-2.AND.
  448.      *                                (IMOD.NE.2.AND.IMOD.NE.4)) THEN
  449.          SIGMF(4)=SIGMF(3)
  450.          SIGMF(3)=0.D0
  451.       ENDIF
  452. C
  453.       IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.3) THEN
  454.         SIGMF(3)=0.D0
  455.         SIGMF(4)=SIGMF(3)
  456.       ENDIF
  457. C
  458.       DO 2 I=1,NSTRSS
  459.         SIGF(I)=SIGMF(I)
  460. C        WRITE(*,*)'SIGMF(',I,')=',SIGMF(I)
  461.     2 CONTINUE
  462. C
  463. C---------------------------------------------------------------------
  464. C
  465. C---------------------------------------------------------------------
  466. C              ACTUALISATION DES VARIABLES INTERNES
  467. C---------------------------------------------------------------------
  468. C       A: MODELE PLASTIQUE
  469. C
  470.       VART( 1)=FLOAT(IFISU)
  471.       VART( 2)=TETA
  472.       VART( 3)=FLOAT(IPLA)
  473.       VART( 4)=SIG1
  474.       VART( 5)=SIG2
  475.       VART( 6)=DPSTF1
  476.       VART( 7)=DPSTF2
  477.       K1 = 7
  478. C
  479. C---------------------------------------------------------------------
  480. C       B: MODELE VISCOPLASTIQUE CLASSIQUE
  481. C
  482.       IF (IVIS.EQ.1) THEN
  483.       VART(8)=TDEF
  484.       VART(9)=TCON
  485.       VART(10)=SIGP(1)
  486.       VART(11)=SIGP(2)
  487.       VART(12)=SIGP(3)
  488.       VART(13)=SIGP(4)
  489.       VART(14)=DPSTV1
  490.       VART(15)=DPSTV2
  491.       VART(16)=SIGV1
  492.       VART(17)=SIGV2
  493.       K1 = 17
  494.       ENDIF
  495. C---------------------------------------------------------------------
  496. C       C: MODELE VISCO-ELASTO PLASTIQUE
  497. C
  498.       IF (IVIS.EQ.2.AND.IFOR.EQ.1) THEN
  499.         MC = NBRC + 1
  500.       DO 103 I1 = 1,MC
  501.       DO 104 J1 = 1,NSTRS
  502.          NV = NSTRS * (I1 - 1) + J1
  503.          VART(7 + NV) = SIR(I1,J1)
  504.   104 CONTINUE
  505.   103 CONTINUE
  506. C      DO 123 I1 = 1,MC
  507. C      DO 124 J1 = 1,MC
  508. C         K = MC * (I1 - 1) + J1
  509. C         VART(43 + K) = CODU(I1,J1)
  510. C  124 CONTINUE
  511. C  123 CONTINUE
  512.       K1=43
  513.       ENDIF
  514. C----------------------------------------
  515.       IF (IVIS.EQ.2.AND.IFOR.EQ.2) THEN
  516.       DO 108 I1 = 1,NBRC
  517.       DO 109 J1 = 1,NSTRS
  518.          NV = NSTRS * (I1 - 1) + J1
  519.               VART(7 + NV) = SIRL(I1,J1)
  520.   109 CONTINUE
  521.   108 CONTINUE
  522. C      DO 143 I1 = 1,NBRC
  523. C      DO 144 J1 = 1,NBRC
  524. C         K = NBRC * (I1-1) + J1
  525. C             VART(43 + K) = CODL(I1,J1)
  526. C  144 CONTINUE
  527. C  143 CONTINUE
  528.       K1=43
  529.         ENDIF
  530. C
  531. C---------------------------------------------------------------------
  532. C       D: MODELE PLASTIQUE ENDOMMAGEABLE
  533. C
  534.         IF (IVIS.EQ.3) THEN
  535.               VART(8) = ENDO
  536.          VART(9) = MAX(VAR0(9),ENDO)
  537.               K1= 9
  538.              ENDIF
  539. C
  540. C---------------------------------------------------------------------
  541. C       B: MODELE VISCOPLASTIQUE ENDOMMAGEABLE
  542. C
  543.       IF (IVIS.EQ.4) THEN
  544.       VART(8)=TDEF
  545.       VART(9)=TCON
  546.       VART(10)=SIGP(1)
  547.       VART(11)=SIGP(2)
  548.       VART(12)=SIGP(3)
  549.       VART(13)=SIGP(4)
  550.       VART(14)=DPSTV1
  551.       VART(15)=DPSTV2
  552.       VART(16)=SIGV1
  553.       VART(17)=SIGV2
  554.       VART(18)=ENDO
  555.       VART(19)=ENDV
  556.       K1 = 19
  557.       ENDIF
  558. C
  559.       GO TO (11,11,31,41),NSTRS
  560.    11 CONTINUE
  561.       KERRE=437
  562.       WRITE(*,*) '!! ATTENTION DANS BETON  NSTRS=',NSTRS
  563.       STOP
  564.       GO TO 1000
  565. C
  566.    41 CONTINUE
  567. C
  568.    31 CONTINUE
  569. C
  570. C     CONTRAINTES OU DEFORMATIONS PLANES
  571. C       SANS CISAILLEMENTS TRANSVERSAL
  572. C
  573.       L1=0
  574.       DO 66 I1 = 1,NSTRS
  575.       DO 67 J1 = 1,NSTRS
  576.          L1=L1+1
  577.          VV(L1) = D(J1,I1)
  578.    67 CONTINUE
  579.    66 CONTINUE
  580. *
  581.       L1=0
  582.       DO 62 I1 = 1,NSTRS
  583.       DO 63 J1 = 1,NSTRS
  584.         K=K1+(I1-1)*NSTRS+J1
  585.         L1=L1+1
  586.         VART(K)=VV(L1)
  587.    63 CONTINUE
  588.    62 CONTINUE
  589. C
  590.       DO 23 II=1,NVARI
  591.         VARF(II)=VART(II)
  592.    23 CONTINUE
  593. C
  594. C---------------------------------------------------------------------
  595. C
  596.  1000 CONTINUE
  597.       SEGSUP BETJEF
  598.       SEGSUP VISCO
  599.       SEGSUP BETFLU
  600.       SEGSUP NECH0
  601.       SEGSUP NECH1
  602.       SEGSUP NECH2
  603.       SEGSUP NECH3
  604. C
  605.       RETURN
  606.       END
  607.  
  608.  
  609.  
  610.  
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  
  619.  
  620.  
  621.  
  622.  

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