Télécharger la plaquette Cast3M


fr en

Télécharger beton.eso

Retour à la liste

Numérotation des lignes : 

beton

  1. C BETON     SOURCE    OF166741  25/11/04    21:15:06     12349          
  2.       SUBROUTINE BETON(SIG0,DEPST,VAR0,XMAT,IVAL,NCOMAT,XCAR,
  3.      1   DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL,IB,IGAU,EPAIST,
  4.      2   MELE,NPINT,SECT,LHOOK,TXR,XLOC,XGLOB,D1HOOK,
  5.      3   ROTHOO,DDHOMU,CRIGI,DSIGT,SIGF,VARF,DEFP,NBPGAU,KERRE
  6.      4, ecou,necou,iecou)
  7. C
  8. C----------------------------------------------------------------------
  9. C             PLASTICITE MODELE BETON
  10. C
  11. C ENTREES
  12. C     SIG0(NSTRS)   = CONTRAINTES INITIALES
  13. C     NSTRS         = NOMBRE DE CONTRAINTES
  14. C     DEPST(NSTRS)  = INCREMENT DE CONTRAINTES TOTALES
  15. C     VAR0(NVARI)   =  VARIABLES INTERNES DEBUT
  16. C     VAR0(   1 )   =  EPSE
  17. C     VAR0(   2 )   =  DAM1
  18. C     VAR0(   3 )   =  DAM2
  19. C     VAR0(   4 )   =  DAM3
  20. C     VAR0(   5 )   =  OUV1
  21. C     VAR0(   6 )   =  OUV2
  22. C     VAR0(   7 )   =  OUV3
  23. C     VAR0(   8 )   =  LAMD
  24. C     VAR0(   9 )   =  VF1(1)
  25. C     VAR0(  10 )   =  VF1(2)
  26. C     VAR0(  11 )   =  VF1(3)
  27. C     VAR0(  12 )   =  VF2(1)
  28. C     VAR0(  13 )   =  VF2(2)
  29. C     VAR0(  14 )   =  VF2(3)
  30. C     VAR0(  15 )   =  VF3(1)
  31. C     VAR0(  16 )   =  VF3(2)
  32. C     VAR0(  17 )   =  VF3(3)
  33. C     XMAT(NCOMAT)  =  COMPOSANTES DE MATERIAU
  34. C     IVAL(NCOMAT)  =  INDICE DES COMPOSANTES DE MATERIAU
  35. C     NCOMAT        = NOMBRE DE COMPOSANTES DE MATERIAU
  36. C     XCAR(ICARA)   =  CARACTERISTIQUES
  37. C     MFR           =  NUMERO DE LA FORMULATION DE L'ELEMENT FINI
  38. C                   = 1       MASSIF
  39. C                   = 3       COQUE MINCE ( COQ2 , COQ3  ET DKT )
  40. C                   = 5       COQUE EPAISSE ( COQ6 , COQ8 )
  41. C                   = 7       POUTRE
  42. C                   = 9       COQUE MINCE AVEC CISAILLEMENT TRANSVERSE
  43. C                                    ( COQ4 ET DST )
  44. C                   = 11      LIQUIDE
  45. C                   = 13      TUYAU
  46. C                   = 15      LINESPRING
  47. C                   = 17      TUYAU FISSURE
  48. C                   = 19      RACCORD MASSIF
  49. C                   = 21      RACCORD COQUE
  50. C                   = 23      SURFACE LIBRE
  51. C                   = 25      MEMBRANE
  52. C                   = 27      UNIAXIALE
  53. C                   = 29      THERMIQUE
  54. C                   = 31      INCOMPRESSIBLES
  55. C                   = 33      POREUX
  56. C                   = 35      JOINT
  57. C                   = 37      HOMOGENEISE
  58. C                   = 39      TUYO
  59. C                   = 41      TUYAU ACOUSTIQUE PURE
  60. C                   = 43      RACCORD TUYAU FLUIDE
  61. *    DDAUX = MATRICE DE HOOKE ELASTIQUE
  62. *    NSTRS = NBRE DE COMPOSANTES DES DEFORMATIONS
  63. *    CMATE = NOM DU MATERIAU
  64. *    VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU
  65. *    VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES
  66. *    N2EL  = NBRE D ELEMENTS DANS SEGMENT DE HOOKE
  67. *    N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE
  68. *    MFR   = NUMERO DE LA FORMULATION
  69. *    IFOU  =
  70. *    IB    = NUMERO DE L ELEMENT COURANT
  71. *    IGAU  = NUMERO DU POINT COURANT
  72. *    EPAIST= EPAISSEUR
  73. *    NBPGAU= NBRE DE POINTS DE GAUSS
  74. *    MELE  = NUMERO DE L ELEMENT FINI
  75. *    NPINT = NBRE DE POINTS D INTEGRATION
  76. *    NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES
  77. *    NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES
  78. *    SECT  = SECTION
  79. *    LHOOK = TAILLE DE LA MATRICE DE HOOKE
  80. *    TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI = TABLEAUX UTILISES
  81. *    POUR LE CALCUL DE LA MATRICE DE HOOKE
  82. C
  83. C  SORTIES
  84. C     SIGF(NSTRS)   = CONTRAINTES FINALES
  85. C     VARF(NVARI)   = VARIABLES INTERNES FINALES
  86. C     DEFP(NSTRS)   = DEFORMATIONS PLASTIQUES
  87. C     KERRE         = 0    TOUT OK
  88. C        1   SI DLAMBDA NEGATIF
  89. C        2   NOMBRE MAX D ITERATIONS INTERNES DEPASSE
  90. C       21   ON NE TROUVE PAS L INTERSECTION AVEC LA SURFACE DE CHARGE
  91. C       22   SIG0 A L EXTERIEUR DE LA SURFACE DE CHARGE
  92. C
  93. C-----------------------------------------------------------------------
  94. C      VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU  ET NECOU
  95. C
  96. C  IFOUR    INDICE DU TYPE DE PROBLEME
  97. C            -3  DEFORMATIONS PLANES GENERALISEES
  98. C            -2  CONTRAINTES PLANES
  99. C            -1  DEFORMATIONS PLANES
  100. C             0  AXISYMETRIQUE
  101. C             1  SERIE DE FOURIER
  102. C             2  TRIDIMENSIONNEL
  103. C  ITYP      TYPE DE FORMULATION MECANIQUE
  104. C  --------------- ATTENTION ---------------
  105. C     IL EST ACTIF APRES L APPEL DE VISAVI
  106. C  -----------------------------------------
  107. C       ITYP=1   CAS DES ELEMENTS MASSIFS
  108. C       ITYP=2   CAS DES COQUES
  109. C       ITYP=3   CAS DES MEMBRANES
  110. C       ITYP=4   CAS DES CABLES ET DES BARRES
  111. C       ITYP=5   CAS QUELCONQUE
  112. C       ITYP=6   CAS DES CONTRAINTES PLANES
  113. C       ITYP=7   CAS DES COQUES A NU=0. OU CONTRAINTES PLANES
  114. C       ITYP=8   CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES
  115. C       ITYP=9   CAS DES COQUES EPAISSES
  116. C       ITYP=10  CAS DES JOINTS
  117. C       ITYP=11  CAS DES POUTRES
  118. C       ITYP=12  CAS DES TUYAUX
  119. C       ITYP=13  CAS DES COQUES AVEC CISAILLEMENT TRANSVERSE
  120. C
  121. C-----------------------------------------------------------------------
  122. C
  123.       IMPLICIT INTEGER(I-N)
  124.       IMPLICIT REAL*8(A-H,O-Z)
  125.  
  126. -INC PPARAM
  127. -INC CCOPTIO
  128.  
  129. -INC TECOU
  130.  
  131.       DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),IVAL(*)
  132.       DIMENSION SIGF(*),VARF(*),DEFP(*)
  133.       DIMENSION SIGMA(6),DSIGMA(6),STOT(6),SIGFIN(6),ORMAT(1)
  134.       DIMENSION XLAMBD(6),DEFPLA(6),DEPS(6),DSIGT(*)
  135.       DIMENSION SPHER(6),AUXIL(6),DSIGZE(6)
  136.       DIMENSION VF1(3),VF2(3),VF3(3),XLTR(3),XLTT(3),EPTR(3),EPTT(3),
  137.      .          EPRS(3),OUVER(3),RT(3)
  138.       DIMENSION WF1(3),WF2(3),WF3(3),YLTR(3),YLTT(3),DETR(3),DETT(3),
  139.      .          DERS(3),FISSU(3),ST(3)
  140.       DIMENSION VALMAT(*),VALCAR(*)
  141.       DIMENSION TXR(IDIM,*),CRIGI(12)
  142.       DIMENSION DDAUX(LHOOK,*),DDHOMU(LHOOK,*)
  143.       DIMENSION XLOC(3,3),XGLOB(3,3)
  144.       DIMENSION D1HOOK(LHOOK,*),ROTHOO(LHOOK,*)
  145.       CHARACTER*8 CMATE
  146.       INTEGER IFIS
  147.  
  148.       IFIS=0
  149.       NSTRS = iecou.NSTRSS
  150.       MFR   = iecou.MFR1
  151. C
  152. C TEST DE CONSISTANCE DES DONNEES
  153. C
  154.       DO 100 I=1,NCOMAT
  155.         IF(XMAT(I).LT.0.D0) KERRE=437
  156.   100 CONTINUE
  157.       IF(KERRE.EQ.437) RETURN
  158. C
  159.       YOUN= 0.D0
  160.       XNU = 0.D0
  161.       RHO = 0.D0
  162.       ALPH= 0.D0
  163.       CTR1= 0.D0
  164.       ETR1= 0.D0
  165.       CTT1= 0.D0
  166.       ETT1= 0.D0
  167.       ERS1= 0.D0
  168.       CTR2= 0.D0
  169.       ETR2= 0.D0
  170.       CTT2= 0.D0
  171.       ETT2= 0.D0
  172.       ERS2= 0.D0
  173.       CTR3= 0.D0
  174.       ETR3= 0.D0
  175.       CTT3= 0.D0
  176.       ETT3= 0.D0
  177.       ERS3= 0.D0
  178.       BETR= 0.D0
  179.       CLCS= 0.D0
  180.       EPCS= 0.D0
  181.       CLBI= 0.D0
  182.       VF1(1)= 0.D0
  183.       VF1(2)= 0.D0
  184.       VF1(3)= 0.D0
  185.       VF2(1)= 0.D0
  186.       VF2(2)= 0.D0
  187.       VF2(3)= 0.D0
  188.       VF3(1)= 0.D0
  189.       VF3(2)= 0.D0
  190.       VF3(3)= 0.D0
  191. *
  192. *    VARIABLE IDECA POUR SE DECALER EN CAS D'EPAISSEUR EN CP
  193. *
  194.       IDECA=1
  195.       IF(IFOUR.EQ.-2) IDECA=2
  196. C
  197.       YOUN= XMAT(1)
  198.       XNU = XMAT(2)
  199.       RHO = XMAT(3)
  200.       ALPH= XMAT(4)
  201.       CTR1= XMAT(5+IDECA)
  202.       ETR1= XMAT(6+IDECA)
  203.       CTT1= XMAT(7+IDECA)
  204.       ETT1= XMAT(8+IDECA)
  205.       ERS1= XMAT(9+IDECA)
  206.       CTR2= XMAT(10+IDECA)
  207.       ETR2= XMAT(11+IDECA)
  208.       CTT2= XMAT(12+IDECA)
  209.       ETT2= XMAT(13+IDECA)
  210.       ERS2= XMAT(14+IDECA)
  211.       CTR3= XMAT(15+IDECA)
  212.       ETR3= XMAT(16+IDECA)
  213.       CTT3= XMAT(17+IDECA)
  214.       ETT3= XMAT(18+IDECA)
  215.       ERS3= XMAT(19+IDECA)
  216.       BETR= XMAT(20+IDECA)
  217.       CLCS= XMAT(21+IDECA)
  218.       EPCS= XMAT(22+IDECA)
  219.       CLBI= XMAT(23+IDECA)
  220.       VF1(1)= XMAT(24+IDECA)
  221.       VF1(2)= XMAT(25+IDECA)
  222.       VF1(3)= XMAT(26+IDECA)
  223.       VF2(1)= XMAT(27+IDECA)
  224.       VF2(2)= XMAT(28+IDECA)
  225.       VF2(3)= XMAT(29+IDECA)
  226.       VF3(1)= XMAT(30+IDECA)
  227.       VF3(2)= XMAT(31+IDECA)
  228.       VF3(3)= XMAT(32+IDECA)
  229. C
  230.       IF(IIMPI.EQ.9) WRITE(IOIMP,2007) (XMAT(I),I=1,32+IDECA)
  231.       IF(IIMPI.EQ.9) WRITE(IOIMP,2008) (IVAL(I),I=1,32+IDECA)
  232. C
  233. C  PETIT TEST SUR NU POUR CERTAINS CAS
  234. C
  235.       IF(MFR.EQ.3.AND.IFOUR.EQ.-2.AND.XNU.NE.0.D0) THEN
  236.         KERRE=38
  237.         RETURN
  238.       ENDIF
  239. C
  240. C VALEURS PAR DEFAUT
  241. C
  242. C                  POUR LA TRACTION
  243. C
  244.       IF(IVAL( 5+IDECA).EQ.0) CTR1=YOUN*1.2E-4
  245.       IF(IVAL( 6+IDECA).EQ.0) ETR1=3.D0*CTR1/YOUN
  246.       IF(IVAL( 7+IDECA).EQ.0) CTT1=0.D0
  247.       IF(IVAL( 8+IDECA).EQ.0) ETT1=ETR1
  248.       IF(IVAL( 9+IDECA).EQ.0) ERS1=0.D0
  249.       IF(IVAL(10+IDECA).EQ.0) THEN
  250.         CTR2=CTR1
  251.         ETR2=ETR1
  252.         CTT2=CTT1
  253.         ETT2=ETT1
  254.         ERS2=ERS1
  255.       ELSE
  256.         IF(IVAL(11+IDECA).EQ.0) ETR2=3.D0*CTR2/YOUN
  257.         IF(IVAL(12+IDECA).EQ.0) CTT2=0.D0
  258.         IF(IVAL(13+IDECA).EQ.0) ETT2=ETR2
  259.         IF(IVAL(14+IDECA).EQ.0) ERS2=0.D0
  260.       ENDIF
  261.       IF(IVAL(15+IDECA).EQ.0) THEN
  262.         CTR3=CTR1
  263.         ETR3=ETR1
  264.         CTT3=CTT1
  265.         ETT3=ETT1
  266.         ERS3=ERS1
  267.       ELSE
  268.       IF(IVAL(16+IDECA).EQ.0) ETR3=3.D0*CTR3/YOUN
  269.       IF(IVAL(17+IDECA).EQ.0) CTT3=0.D0
  270.       IF(IVAL(18+IDECA).EQ.0) ETT3=ETR3
  271.       IF(IVAL(19+IDECA).EQ.0) ERS3=0.D0
  272.       ENDIF
  273. *
  274.       IF(MFR.EQ.1.AND.(IFOUR.GE.-1.OR.IFOUR.EQ.-3)) THEN
  275. C
  276.         IF(IVAL(24+IDECA).EQ.0.AND.IVAL(27+IDECA).EQ.0
  277.      &                        .AND.IVAL(30+IDECA).EQ.0) NVEC=0
  278.         IF(IVAL(24+IDECA).NE.0.OR.IVAL(27+IDECA).NE.0
  279.      &                        .OR.IVAL(30+IDECA).NE.0)  NVEC=1
  280.         IF(IVAL(24+IDECA).NE.0.AND.IVAL(27+IDECA).NE.0) NVEC=2
  281. C
  282.         KAS=1
  283.         IF(IVAL(27+IDECA).NE.0.AND.NVEC.EQ.1) KAS=2
  284.         IF(IVAL(30+IDECA).NE.0.AND.NVEC.EQ.1) KAS=3
  285.         IF(IFOUR.EQ.-3.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) THEN
  286.           NVEC=1
  287.           KAS=3
  288.           VAR0(15)=0.D0
  289.           VAR0(16)=0.D0
  290.           VAR0(17)=1.D0
  291.           VF3(1)=0.D0
  292.           VF3(2)=0.D0
  293.           VF3(3)=1.D0
  294.         ENDIF
  295. C
  296.       ENDIF
  297. C
  298.       IF(IVAL(20+IDECA).EQ.0)  BETR=0.1D0
  299. C
  300. C VALEURS PAR DEFAUT POUR LE MODELE BETON CONTRAINTES PLANES
  301. C
  302. C                  POUR LA TRACTION
  303. C
  304.       IF(MFR.EQ.3.OR.MFR.EQ.9.OR.IFOUR.EQ.-2) THEN
  305.         IFIS=0
  306.         IF(IVAL(5+IDECA).NE.0.AND.IVAL(10+IDECA).NE.0) IFIS=1
  307.         ANGL=0.D0
  308.         IF(IVAL(24+IDECA).NE.0.AND.IVAL(25+IDECA).NE.0.AND.IFIS.EQ.1)
  309.      .     ANGL=ATAN2(VF1(2),VF1(1))
  310. C
  311. C                  POUR LA COMPRESSION
  312. C
  313.         IF(IVAL(21+IDECA).EQ.0) CLCS=YOUN*1.2D-3
  314.         IF(IVAL(22+IDECA).EQ.0) EPCS=10.D0*CLCS/YOUN
  315. C
  316. C                  BI-COMPRESSION
  317. C
  318.         IF(IVAL(23+IDECA).EQ.0) THEN
  319.           REP=EPCS*YOUN/CLCS
  320.           UMN=1.D0-XNU
  321.           ALF=((REP-UMN-(SQRT(0.5D0*REP*UMN)))/(2.D0*REP-UMN))*0.95D0
  322.           DDD=(1.D0-ALF)/(1.D0-2.D0*ALF)
  323.           CLBI=DDD*CLCS
  324.         ENDIF
  325.       ENDIF
  326. C
  327. C INITIALISATION
  328. C
  329.       JEPRIN = 0
  330.       UNIT=0.0174532925199432957692D0
  331.       ANGLE=0.D0
  332.       IF(IFIS.EQ.1) ANGLE=ANGL
  333. C
  334.       XLTR(1)=CTR1
  335.       XLTR(2)=CTR2
  336.       XLTR(3)=CTR3
  337.       XLTT(1)=CTT1
  338.       XLTT(2)=CTT2
  339.       XLTT(3)=CTT3
  340.       EPTR(1)=ETR1
  341.       EPTR(2)=ETR2
  342.       EPTR(3)=ETR3
  343.       EPTT(1)=ETT1
  344.       EPTT(2)=ETT2
  345.       EPTT(3)=ETT3
  346.       EPRS(1)=ERS1
  347.       EPRS(2)=ERS2
  348.       EPRS(3)=ERS3
  349. C
  350. C CALCUL ET TEST DES PARAMETRES DU BETON
  351. C
  352. C     CALCUL DES BRANCHES DESCENDANTES DANS LE DOMAINE DE LA TRACTION
  353. C     POUR LE MODELE BETON CONTRAINTES PLANES
  354. C
  355.       IF(MFR.EQ.3.OR.MFR.EQ.9.OR.IFOUR.EQ.-2) THEN
  356. C
  357.         IF(IFIS.EQ.0) THEN
  358.           R=ABS(CTR1-CTR2)
  359.           IF(R.GT.1.D-7) THEN
  360.             KERRE=441
  361.             RETURN
  362.           ENDIF
  363.         ENDIF
  364. C
  365.         EP1=ETR1-VAR0(5)
  366.         EP2=ETR2-VAR0(6)
  367. C
  368.         IF(EP1.LE.1.D-10) THEN
  369.           CTR1=0.D0
  370.           PENTT1=0.D0
  371.         ENDIF
  372. C
  373.         IF(EP2.LE.1.D-10) THEN
  374.           CTR2=0.D0
  375.           PENTT2=0.D0
  376.         ENDIF
  377. C
  378.         EP1=ETR1-CTR1/YOUN/(1-XNU)
  379.         EP2=ETR2-CTR2/YOUN/(1-XNU)
  380. C
  381.         IF(EP1.LT.1.D-10) THEN
  382.           KERRE=438
  383.           RETURN
  384.         ENDIF
  385. C
  386.         IF(EP2.LT.1.D-10) THEN
  387.           KERRE=439
  388.           RETURN
  389.         ENDIF
  390. C
  391.         IF(EP1.GT.1.D-10) PENTT1=-CTR1/ETR1
  392.         IF(EP2.GT.1.D-10) PENTT2=-CTR2/ETR2
  393. C
  394. C     CALCUL DES PARAMETRES DU DRUCKER POUR LE MODELE
  395. C     BETON CONTRAINTES PLANES
  396. C
  397.         DI1=CLCS-CLBI
  398.         IF(DI1.GE.0.D0) THEN
  399.           KERRE=443
  400.           RETURN
  401.         ENDIF
  402. C
  403.         DI2=CLCS-2.D0*CLBI
  404.         ALFAD2=DI1/DI2
  405.         DPELA2=(1.D0-ALFAD2)*CLCS
  406. C
  407.         R=DPELA2/(1.D0+ALFAD2)
  408.         IF(R.LT.CTR1.OR.R.LT.CTR2.AND.IIMPI.EQ.9) WRITE(IOIMP,900)
  409. C
  410.         R=ABS(PENTT1-PENTT2)
  411.         IF(R.GT.1.D-7.AND.CTR1.NE.0.D0.AND.CTR2.NE.0.D0.
  412.      .     AND.IIMPI.EQ.9) WRITE(IOIMP,901)
  413. C
  414.         EPSA=CLCS/YOUN
  415.         PENTE2=-CLCS/(EPCS-EPSA)
  416.         PENTE2=(YOUN*PENTE2)/(YOUN-PENTE2)
  417.         PENTE2=PENTE2*(ALFAD2-1.D0)*(ALFAD2-1.D0)
  418. C
  419.       ENDIF
  420. C
  421.       IF(BETR.GT.1.D0) THEN
  422.         KERRE=442
  423.         RETURN
  424.        ENDIF
  425. C
  426. C     CALCUL DES CONSTANTES
  427. C
  428.       EPS1=CTR1/YOUN
  429.       EPS2=CTR2/YOUN
  430.       EPS3=CTR3/YOUN
  431. C
  432.       TRTRAC=CTR1+CTR2+CTR3
  433. C
  434.       SIGTRA=TRTRAC*0.333333333333333333333333333333333D0
  435. C
  436. C     TEST SUR LES VALEURS DES PARAMETRES
  437. C
  438.       IF(MFR.EQ.3.OR.MFR.EQ.9.OR.IFOUR.EQ.-2) THEN
  439. C
  440.         IF(ETR1.LT.EPS1) THEN
  441.           KERRE=449
  442.           REAERR(1)=ETR1
  443.           REAERR(2)=EPS1
  444.           RETURN
  445.         ENDIF
  446. C
  447.         IF(ETR2.LT.EPS2) THEN
  448.           KERRE=450
  449.           REAERR(1)=ETR2
  450.           REAERR(2)=EPS2
  451.           RETURN
  452.         ENDIF
  453. C
  454.         IF(EPCS.LT.EPSA) THEN
  455.           KERRE=444
  456.           REAERR(1)=EPCS
  457.           REAERR(2)=EPSA
  458.           RETURN
  459.         ENDIF
  460. C
  461.         IF(CLCS.LT.SIGTRA) THEN
  462.           KERRE=445
  463.           REAERR(1)=CLCS
  464.           REAERR(2)=SIGTRA
  465.           RETURN
  466.         ENDIF
  467. C
  468.         ECR2=(2.D0*YOUN*(0.5D0-ALFAD2)*(0.5D0-ALFAD2))/(1.D0-XNU)
  469.         IF(ABS(PENTE2).GE.ECR2) THEN
  470.           KERRE=458
  471.           REAERR(1)=PENTE2
  472.           REAERR(2)=ECR2
  473.           RETURN
  474.         ENDIF
  475.       ENDIF
  476. C
  477.       IF(MFR.EQ.1.AND.(IFOUR.GE.-1.OR.IFOUR.EQ.-3)) THEN
  478. C
  479.         IF(CTR1.LT.CTT1) THEN
  480.           KERRE=446
  481.           REAERR(1)=CTR1
  482.           REAERR(2)=CTT1
  483.           RETURN
  484.         ENDIF
  485. C
  486.         IF(CTR2.LT.CTT2) THEN
  487.           KERRE=447
  488.           REAERR(1)=CTR2
  489.           REAERR(2)=CTT2
  490.           RETURN
  491.         ENDIF
  492. C
  493.         IF(CTR3.LT.CTT3) THEN
  494.           KERRE=448
  495.           REAERR(1)=CTR3
  496.           REAERR(2)=CTT3
  497.           RETURN
  498.         ENDIF
  499. C
  500.         IF(ETR1.LT.EPS1) THEN
  501.           KERRE=449
  502.           REAERR(1)=ETR1
  503.           REAERR(2)=EPS1
  504.           RETURN
  505.         ENDIF
  506. C
  507.         IF(ETR2.LT.EPS2) THEN
  508.           KERRE=450
  509.           REAERR(1)=ETR2
  510.           REAERR(2)=EPS2
  511.           RETURN
  512.         ENDIF
  513. C
  514.         IF(ETR3.LT.EPS3) THEN
  515.           KERRE=451
  516.           REAERR(1)=ETR3
  517.           REAERR(2)=EPS3
  518.           RETURN
  519.         ENDIF
  520. C
  521.         IF(ETR1.LT.ETT1.OR.ETT1.LT.EPS1) THEN
  522.           KERRE=452
  523.           REAERR(1)=ETR1
  524.           REAERR(2)=ETT1
  525.           REAERR(3)=EPS1
  526.           RETURN
  527.         ENDIF
  528. C
  529.         IF(ETR2.LT.ETT2.OR.ETT2.LT.EPS2) THEN
  530.           KERRE=453
  531.           REAERR(1)=ETR2
  532.           REAERR(2)=ETT2
  533.           REAERR(3)=EPS2
  534.           RETURN
  535.         ENDIF
  536. C
  537.         IF(ETR3.LT.ETT3.OR.ETT3.LT.EPS3) THEN
  538.           KERRE=454
  539.           REAERR(1)=ETR3
  540.           REAERR(2)=ETT3
  541.           REAERR(3)=EPS3
  542.           RETURN
  543.         ENDIF
  544. C
  545.         IF(ETR1.LT.ERS1) THEN
  546.           KERRE=455
  547.           REAERR(1)=ETR1
  548.           REAERR(2)=ERS1
  549.           RETURN
  550.         ENDIF
  551. C
  552.         IF(ETR2.LT.ERS2) THEN
  553.           KERRE=456
  554.           REAERR(1)=ETR2
  555.           REAERR(2)=ERS2
  556.           RETURN
  557.         ENDIF
  558. C
  559.         IF(ETR3.LT.ERS3) THEN
  560.           KERRE=457
  561.           REAERR(1)=ETR3
  562.           REAERR(2)=ERS3
  563.           RETURN
  564.         ENDIF
  565. C
  566.       ENDIF
  567. C
  568. C  CALCUL DE L INCREMENT DE CONTRAINTES
  569. C
  570.       CALL CALSIG(DEPST,DDAUX,NSTRS,CMATE,VALMAT,
  571.      1 VALCAR,N2EL,N2PTEL,MFR,necou.IFOURB,IB,IGAU,EPAIST,
  572.      2 NBPGAU,MELE,NPINT,iecou.NBGMAT,iecou.NELMAT,SECT,LHOOK,TXR,
  573.      3 XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,IRTD)
  574. *
  575.        IF(IRTD.NE.1) THEN
  576.           KERRE=69
  577.           GOTO 998
  578.        ENDIF
  579. *
  580. C
  581. C ECOULEMENT PLASTIQUE
  582. C
  583.   333 MCOD=1
  584.       CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  585.      .        SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  586.      .        DSIGZE,ICENT2,MCOD,IBOU,MFR,NSTRS,XCAR,CMATE,ecou,necou)
  587.       IF(ITYP.EQ.0) THEN
  588.         KERRE=269
  589.         RETURN
  590.       ENDIF
  591. *
  592. C
  593. C IMPRESSION DES PARAMETRES DU MODELE
  594. C
  595.       IF(IIMPI.EQ.9) THEN
  596.         WRITE(IOIMP,1000)
  597.         WRITE(IOIMP,1001) YOUN
  598.         WRITE(IOIMP,1002) XNU
  599.         WRITE(IOIMP,1003)
  600.         IF(ITYP.EQ.2.OR.ITYP.EQ.6.OR.ITYP.EQ.7) THEN
  601.           WRITE(IOIMP,1004) CTR1,CTR2
  602.           WRITE(IOIMP,1005) ETR1,ETR2
  603.           WRITE(IOIMP,1006) VAR0(5),VAR0(6)
  604.           WRITE(IOIMP,1007) (ANGLE/UNIT)
  605.         ENDIF
  606.         IF(MFR.EQ.1.AND.(IFOUR.GE.-1.OR.IFOUR.EQ.-3)) THEN
  607.           WRITE(IOIMP,1008) CTR1,CTR2,CTR3
  608.           WRITE(IOIMP,1009) CTT1,CTT2,CTT3
  609.           WRITE(IOIMP,1010) ETT1,ETT2,ETT3
  610.           WRITE(IOIMP,1011) ETR1,ETR2,ETR3
  611.           WRITE(IOIMP,1012) ERS1,ERS2,ERS3
  612.           WRITE(IOIMP,1013) VAR0(5),VAR0(6),VAR0(7)
  613.           IF((KAS.EQ.1.AND.NVEC.EQ.1).OR.NVEC.EQ.2)
  614.      .      WRITE(IOIMP,1014) (VF1(I),I=1,3)
  615.           IF((KAS.EQ.2.AND.NVEC.EQ.1).OR.NVEC.EQ.2)
  616.      .      WRITE(IOIMP,1015) (VF2(I),I=1,3)
  617.           IF((KAS.EQ.3.AND.NVEC.EQ.1).OR.NVEC.EQ.2)
  618.      .      WRITE(IOIMP,1016) (VF3(I),I=1,3)
  619.         ENDIF
  620.         WRITE(IOIMP,1017) BETR
  621.         IF(ITYP.EQ.2.OR.ITYP.EQ.6.OR.ITYP.EQ.7) WRITE(IOIMP,1018)
  622.         IF(MFR.EQ.1.AND.(IFOUR.GE.-1.OR.IFOUR.EQ.-3)) WRITE(IOIMP,1019)
  623.         IF(ITYP.EQ.2.OR.ITYP.EQ.6.OR.ITYP.EQ.7) THEN
  624.           WRITE(IOIMP,1020)
  625.           WRITE(IOIMP,1021) CLCS,EPCS
  626.           WRITE(IOIMP,1022) CLBI
  627.         ENDIF
  628.       ENDIF
  629. C
  630. C CAS DES CONTRAINTES PLANES :
  631. C             IFOUR = -2
  632. C AVEC DEUX TYPES DE FORMULATION MECANIQUE :
  633. C          CAS DES COQUES :             ITYP = 2 ==> ALFAH = 0
  634. C                                   OU  ITYP = 7 ==> ALFAH = 0
  635. C          CAS DES CONTRAINTES PLANES : ITYP = 6
  636. C
  637.       IF(ITYP.EQ.2.OR.ITYP.EQ.6.OR.ITYP.EQ.7) THEN
  638. C
  639. C DANS LE CAS DES COQUES ITYP = 2  , NOUS METTONS ALFAH = 0
  640. C
  641.         IF(ITYP.EQ.2.OR.ITYP.EQ.7) ALFAH=0.D0
  642. C
  643. C INITIALISATIONS DES PARAMETRES
  644. C
  645.         IF(CTR1.GT.0.D0) RT1=(1.D0-VAR0(2))*CTR1
  646.         IF(CTR1.LE.0.D0) RT1=0.D0
  647.         IF(CTR2.GT.0.D0) RT2=(1.D0-VAR0(3))*CTR2
  648.         IF(CTR2.LE.0.D0) RT2=0.D0
  649.         XLAMBD(1)=VAR0(5)
  650.         XLAMBD(2)=VAR0(6)
  651.         XLAMBD(3)=VAR0(8)
  652.         IF(IVAL(24+IDECA).EQ.0) VF1(1)=VAR0(9)
  653.         IF(IVAL(25+IDECA).EQ.0) VF1(2)=VAR0(10)
  654. C
  655.         IF(IFIS.EQ.0.AND.(XLAMBD(1).GT.0.D0.OR.XLAMBD(2).GT.0D0)) THEN
  656.           IFIS=1
  657.           ANGLE=ATAN2(VF1(2),VF1(1))
  658.         ENDIF
  659. C
  660.         DO 200 IBA=1,IBOU
  661.         STOT(IBA)=SIGMA(IBA)+DSIGMA(IBA)
  662.   200   CONTINUE
  663. C
  664. C ECOULEMENT PLASTIQUE
  665. C
  666.         CALL ECBECP(SIGMA,DSIGMA,XLAMBD,DEFPLA,YOUN,XNU,
  667.      .              IFIS,ANGLE,BETR,RT1,PENTT1,RT2,PENTT2,
  668.      .              DPELA2,ALFAD2,PENTE2,IDAM,SIGFIN,ITYP,KERRE)
  669. C
  670.         IF(IIMPI.EQ.9) THEN
  671.           WRITE(IOIMP,2000)
  672.           WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  673.           WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  674.           WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  675.           WRITE(IOIMP,2004) (SIGFIN(I),I=1,IBOU)
  676.           WRITE(IOIMP,2005) (DEFPLA(I),I=1,IBOU)
  677.           WRITE(IOIMP,2006) (XLAMBD(I),I=1,6)
  678.         ENDIF
  679. C
  680.         IF(KERRE.NE.0) THEN
  681.            IF(JEPRIN.EQ.0) THEN
  682.              KERRE=0
  683.              JEPRIN=1
  684.              IIMPI=9
  685.              GO TO 333
  686.            ELSE
  687.              IIMPI=0
  688.              RETURN
  689.            ENDIF
  690.         ENDIF
  691. C
  692. C PAS DE PLASTICITE COMPORTEMENT ELASTIQUE
  693. C
  694.         IF(IDAM.EQ.0) THEN
  695.           DO 300 IBA=1,IBOU
  696.           SIGFIN(IBA)=STOT(IBA)
  697.           DEFPLA(IBA)=0.D0
  698.   300     CONTINUE
  699.         ENDIF
  700. C
  701. C NOUS AVONS ENDOMMAGE LE MATERIAU
  702. C
  703. C MISE A JOUR DES VARIABLES INTERNES
  704. C
  705.         VARF(1)=SQRT(XLAMBD(1)*XLAMBD(1)+XLAMBD(2)*XLAMBD(2)+
  706.      .                XLAMBD(3)*XLAMBD(3))
  707.         IF(CTR1.GT.0.D0) VARF(2)=1.D0-RT1/CTR1
  708.         IF(CTR1.LE.0.D0) VARF(2)=1.D0
  709.         IF(CTR2.GT.0.D0) VARF(3)=1.D0-RT2/CTR2
  710.         IF(CTR2.LE.0.D0) VARF(3)=1.D0
  711.         VARF(5)=XLAMBD(1)
  712.         VARF(6)=XLAMBD(2)
  713.         VARF(8)=XLAMBD(3)
  714.         VARF(9)=COS(ANGLE)
  715.         VARF(10)=SIN(ANGLE)
  716. C
  717. C     CAS DES DEFORMATIONS PLANES , AXISYMETRIQUE ET TRIDIME
  718. C
  719.       ELSEIF (ITYP.EQ.1) THEN
  720. C
  721. C       INITIALISATIONS DES PARAMETRES
  722. C
  723.         IF(CTR1.GT.0.D0) RT(1)=(1.D0-VAR0(2))*CTR1
  724.         IF(CTR1.LE.0.D0) RT(1)=0.D0
  725.         IF(CTR2.GT.0.D0) RT(2)=(1.D0-VAR0(3))*CTR2
  726.         IF(CTR2.LE.0.D0) RT(2)=0.D0
  727.         IF(CTR3.GT.0.D0) RT(3)=(1.D0-VAR0(4))*CTR3
  728.         IF(CTR3.LE.0.D0) RT(3)=0.D0
  729.         OUVER(1)=VAR0(5)
  730.         OUVER(2)=VAR0(6)
  731.         OUVER(3)=VAR0(7)
  732.         IF(IVAL(24+IDECA).EQ.0) VF1(1)=VAR0(9)
  733.         IF(IVAL(25+IDECA).EQ.0) VF1(2)=VAR0(10)
  734.         IF(IVAL(26+IDECA).EQ.0) VF1(3)=VAR0(11)
  735.         IF(IVAL(27+IDECA).EQ.0) VF2(1)=VAR0(12)
  736.         IF(IVAL(28+IDECA).EQ.0) VF2(2)=VAR0(13)
  737.         IF(IVAL(29+IDECA).EQ.0) VF2(3)=VAR0(14)
  738.         IF(IVAL(30+IDECA).EQ.0) VF3(1)=VAR0(15)
  739.         IF(IVAL(31+IDECA).EQ.0) VF3(2)=VAR0(16)
  740.         IF(IVAL(32+IDECA).EQ.0) VF3(3)=VAR0(17)
  741. C
  742.         IF(NVEC.EQ.2) THEN
  743.           NBVECD=NVEC
  744.         ELSE
  745.           NBVECD=0
  746.           CALL NORME(VF1,P1)
  747.           CALL NORME(VF2,P2)
  748.           CALL NORME(VF3,P3)
  749.           IF(P1.GT.0.D0) NBVECD=NBVECD+1
  750.           IF(P2.GT.0.D0) NBVECD=NBVECD+1
  751.           IF(P3.GT.0.D0) NBVECD=NBVECD+1
  752.           IF(NBVECD.GT.2) NBVECD=2
  753.         ENDIF
  754. C
  755. C       CALCUL DE DEFORMATION PLASTIQUE
  756. C
  757.         CALL EPSIG(DSIGMA,DEPS,IFOUR,YOUN,XNU,ITYP,ORMAT,XCAR)
  758. C
  759.         IF(IIMPI.EQ.9) THEN
  760.           WRITE(IOIMP,*) 'AVANT L APPEL A TENSI3 AVANT L APPEL A AJUSTE'
  761.           WRITE(IOIMP,*) 'NBVECD',NBVECD
  762.           WRITE(IOIMP,*) 'ANGLE ',(ANGLE/UNIT)
  763.           WRITE(IOIMP,*) 'VF1   ',(VF1(I),I=1,3)
  764.           WRITE(IOIMP,*) 'VF2   ',(VF2(I),I=1,3)
  765.           WRITE(IOIMP,*) 'VF3   ',(VF3(I),I=1,3)
  766.         ENDIF
  767. C
  768.         MODE=1
  769.         CALL AJUSTE (VF1,VF2,VF3,XLTR,EPTR,XLTT,EPTT,EPRS,RT,OUVER,
  770.      .               WF1,WF2,WF3,YLTR,DETR,YLTT,DETT,DERS,ST,FISSU,
  771.      .               ANGLE,NBVECD,KAS,MODE)
  772. C
  773.         CALL TENSI3 (SIGMA,DSIGMA,DEPS,YOUN,XNU,BETR,
  774.      .               WF1,WF2,WF3,ANGLE,NBVECD,
  775.      .               YLTR,DETR,YLTT,DETT,DERS,ST,FISSU,
  776.      .               SIGFIN,DEFPLA,KERRE)
  777. C
  778.         IF(IIMPI.EQ.9) THEN
  779.           WRITE(IOIMP,*) 'APRES L APPEL A TENSI3 AVANT L APPEL A AJUSTE'
  780.           WRITE(IOIMP,*) 'NBVECD',NBVECD
  781.           WRITE(IOIMP,*) 'ANGLE ',(ANGLE/UNIT)
  782.           WRITE(IOIMP,*) 'W1    ',(WF1(I),I=1,3)
  783.           WRITE(IOIMP,*) 'W2    ',(WF2(I),I=1,3)
  784.           WRITE(IOIMP,*) 'W3    ',(WF3(I),I=1,3)
  785.         ENDIF
  786. C
  787.         MODE=2
  788.         CALL AJUSTE (VF1,VF2,VF3,XLTR,EPTR,XLTT,EPTT,EPRS,RT,OUVER,
  789.      .               WF1,WF2,WF3,YLTR,DETR,YLTT,DETT,DERS,ST,FISSU,
  790.      .               ANGLE,NBVECD,KAS,MODE)
  791. C
  792.         IF(IIMPI.EQ.9) THEN
  793.           WRITE(IOIMP,*) 'KAS   ',KAS
  794.           WRITE(IOIMP,*) 'SIGMA ',(SIGMA(I),I=1,IBOU)
  795.           WRITE(IOIMP,*) 'DSIGMA',(DSIGMA(I),I=1,IBOU)
  796.           WRITE(IOIMP,*) 'DEPS ',(DEPS(I),I=1,IBOU)
  797.           WRITE(IOIMP,*) 'SIGFIN',(SIGFIN(I),I=1,IBOU)
  798.           WRITE(IOIMP,*) 'DEFPLA',(DEFPLA(I),I=1,IBOU)
  799.           WRITE(IOIMP,*) 'OUVER ',(OUVER(I),I=1,3)
  800.           WRITE(IOIMP,*) 'RT    ',(RT(I),I=1,3)
  801.           WRITE(IOIMP,*) 'NBVECD',NBVECD
  802.           WRITE(IOIMP,*) 'ANGLE ',(ANGLE/UNIT)
  803.           WRITE(IOIMP,*) 'VF1   ',(VF1(I),I=1,3)
  804.           WRITE(IOIMP,*) 'VF2   ',(VF2(I),I=1,3)
  805.           WRITE(IOIMP,*) 'VF3   ',(VF3(I),I=1,3)
  806.         ENDIF
  807. C
  808.         IF(KERRE.NE.0) THEN
  809.            IF(JEPRIN.EQ.0) THEN
  810.              KERRE=0
  811.              JEPRIN=1
  812.              IIMPI=9
  813.              GO TO 333
  814.            ELSE
  815.              IIMPI=0
  816.              RETURN
  817.            ENDIF
  818.         ENDIF
  819. C
  820. C MISE A JOUR DES VARIABLES INTERNES
  821. C
  822.         VARF(1)=SQRT(OUVER(1)*OUVER(1)+OUVER(2)*OUVER(2)+
  823.      .                OUVER(3)*OUVER(3))
  824.         IF (CTR1.GT.0.D0) THEN
  825.           VARF(2)=1.D0-RT(1)/CTR1
  826.         ELSE
  827.           VARF(2)=1.D0
  828.         ENDIF
  829.         IF (CTR2.GT.0.D0) THEN
  830.           VARF(3)=1.D0-RT(2)/CTR2
  831.         ELSE
  832.           VARF(3)=1.D0
  833.         ENDIF
  834.         IF (CTR3.GT.0.D0) THEN
  835.           VARF(4)=1.D0-RT(3)/CTR3
  836.         ELSE
  837.           VARF(4)=1.D0
  838.         ENDIF
  839.         VARF(5)=OUVER(1)
  840.         VARF(6)=OUVER(2)
  841.         VARF(7)=OUVER(3)
  842.         VARF( 9)=VF1(1)
  843.         VARF(10)=VF1(2)
  844.         VARF(11)=VF1(3)
  845.         VARF(12)=VF2(1)
  846.         VARF(13)=VF2(2)
  847.         VARF(14)=VF2(3)
  848.         VARF(15)=VF3(1)
  849.         VARF(16)=VF3(2)
  850.         VARF(17)=VF3(3)
  851.       ELSE
  852.         KERRE=463
  853.         RETURN
  854.       ENDIF
  855. C
  856.       MCOD=2
  857.       CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  858.      .        SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  859.      .        DSIGZE,ICENT2,MCOD,IBOU,MFR,NSTRS,XCAR,CMATE,ecou,necou)
  860.  
  861. C
  862. C LES FORMATS D IMPRESSION
  863. C
  864.   900 FORMAT(1X,'ATTENTION LES CRITERES DE TRACTION SERONT INACTIFS',/,
  865.      * 1X,'LIMITE EN TRACTION TROP ELEVE SEUL LE DRUCKER EST ACTIF',/)
  866.   901 FORMAT(1X,'ATTENTION DEFORMATIONS A RUPTURE ET LIMITES EN',/,
  867.      * 1X,'TRACTION SONT INCOMPATIBLES POUR AVOIR UN MATERIAU',/,
  868.      * 1X,'ISOTROPE INITIALEMENT',/)
  869. C
  870.  1000 FORMAT(////,31X,'MATERIAU BETON',//)
  871.  1001 FORMAT(6X,'MODULE D''YOUNG',T51,'= ',1PE12.5)
  872.  1002 FORMAT(6X,'COEF. DE POISSON',T51,'= ',1PE12.5)
  873.  1003 FORMAT(/,1X,'DOMMAGE PAR TRACTION :',/)
  874.  1004 FORMAT(6X,'CONTRAINTES LIMITES :',/
  875.      *     11X,'DIRECTION(1) = ',1PE12.5,
  876.      *      5X,'DIRECTION(2) = ',1PE12.5/)
  877.  1005 FORMAT(6X,'DEFORMATIONS A RUPTURE EN TRACTION :',/
  878.      *     11X,'DIRECTION(1) = ',1PE12.5,
  879.      *      5X,'DIRECTION(2) = ',1PE12.5/)
  880.  1006 FORMAT(6X,'OUVERTURES DES FISSURES :',/
  881.      *     11X,'DIRECTION(1) = ',1PE12.5,
  882.      *      5X,'DIRECTION(2) = ',1PE12.5/)
  883.  1007 FORMAT(6X,'ANGLE DE FISSURATION',T51,'= ',1PE12.5/)
  884.  1008 FORMAT(6X,'CONTRAINTES LIMITES :',/
  885.      *     11X,'DIRECTION(1) = ',1PE12.5,
  886.      *      5X,'DIRECTION(2) = ',1PE12.5,5X,'DIRECTION(3) = ',1PE12.5/)
  887.  1009 FORMAT(6X,'CONTRAINTES DE TRANSITION :',/
  888.      *     11X,'DIRECTION(1) = ',1PE12.5,
  889.      *      5X,'DIRECTION(2) = ',1PE12.5,5X,'DIRECTION(3) = ',1PE12.5/)
  890.  1010 FORMAT(6X,'DEFORMATIONS DE TRANSITION :',/
  891.      *     11X,'DIRECTION(1) = ',1PE12.5,
  892.      *      5X,'DIRECTION(2) = ',1PE12.5,5X,'DIRECTION(3) = ',1PE12.5/)
  893.  1011 FORMAT(6X,'DEFORMATIONS A RUPTURE EN TRACTION :',/
  894.      *     11X,'DIRECTION(1) = ',1PE12.5,
  895.      *      5X,'DIRECTION(2) = ',1PE12.5,5X,'DIRECTION(3) = ',1PE12.5/)
  896.  1012 FORMAT(6X,'DEFORMATIONS RESIDUELLES :',/
  897.      *     11X,'DIRECTION(1) = ',1PE12.5,
  898.      *      5X,'DIRECTION(2) = ',1PE12.5,5X,'DIRECTION(3) = ',1PE12.5/)
  899.  1013 FORMAT(6X,'OUVERTURES DES FISSURES :',/
  900.      *     11X,'DIRECTION(1) = ',1PE12.5,
  901.      *      5X,'DIRECTION(2) = ',1PE12.5,5X,'DIRECTION(3) = ',1PE12.5/)
  902.  1014 FORMAT(6X,'COMPOSANTES DU VECTEUR VF1 :',/
  903.      *     11X,'DIRECTION(X) = ',1PE12.5,
  904.      *      5X,'DIRECTION(Y) = ',1PE12.5,5X,'DIRECTION(Z) = ',1PE12.5/)
  905.  1015 FORMAT(6X,'COMPOSANTES DU VECTEUR VF2 :',/
  906.      *     11X,'DIRECTION(X) = ',1PE12.5,
  907.      *      5X,'DIRECTION(Y) = ',1PE12.5,5X,'DIRECTION(Z) = ',1PE12.5/)
  908.  1016 FORMAT(6X,'COMPOSANTES DU VECTEUR VF3 :',/
  909.      *     11X,'DIRECTION(X) = ',1PE12.5,
  910.      *      5X,'DIRECTION(Y) = ',1PE12.5,5X,'DIRECTION(Z) = ',1PE12.5/)
  911.  1017 FORMAT(6X,'CISAILLEMENT RESIDUEL APRES FISSURATION  ',
  912.      *     T51,'= ',1PE12.5/)
  913.  1018 FORMAT(1X,'LA DECHARGE SERA EFFECTUEE AVEC LE MODULE ',
  914.      *'ELASTIQUE DANS LE DOMAINE DE LA TRACTION',/)
  915.  1019 FORMAT(1X,'LA DECHARGE SERA EFFECTUEE AVEC UN MODULE ',
  916.      *'D ENDOMMAGEMENT DANS LE DOMAINE DE LA TRACTION',/)
  917.  1020 FORMAT(/,1X,'DOMMAGE PAR COMPRESSION :',/)
  918.  1021 FORMAT(6X,'COMPRESSION SIMPLE :',/
  919.      *     11X,'CONTRAINTE LIMITE = ',1PE12.5,
  920.      *      5X,'DEFORMATION A RUPTURE   = ',1PE12.5/)
  921.  1022 FORMAT(6X,'COMPRESSION BIAXIALE :',/
  922.      *     11X,'CONTRAINTE LIMITE = ',1PE12.5/)
  923. C
  924.  2000 FORMAT(1X,'RESULTATS DE L ECOULEMENT PLASTIQUE')
  925.  2001 FORMAT(1X,'SIGMA  =',6(1X,1PE12.5))
  926.  2002 FORMAT(1X,'DSIGMA =',6(1X,1PE12.5))
  927.  2003 FORMAT(1X,'STOT   =',6(1X,1PE12.5))
  928.  2004 FORMAT(1X,'SIGFIN =',6(1X,1PE12.5))
  929.  2005 FORMAT(1X,'DEFPLA =',6(1X,1PE12.5))
  930.  2006 FORMAT(1X,'XLAMBD =',6(1X,1PE12.5))
  931.  2007 FORMAT(7(6(1X,E12.5),/))
  932.  2008 FORMAT(7(6(1X,I12),/))
  933. C
  934.  998  RETURN
  935.       END
  936.  
  937.  
  938.  

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