Télécharger la plaquette Cast3M


fr en

Télécharger betaxi.eso

Retour à la liste

Numérotation des lignes : 

betaxi

  1. C BETAXI    SOURCE    OF166741  25/11/04    21:15:05     12349          
  2.       SUBROUTINE BETAXI(SIG0,NSTRS,DSIGT,VAR0,XMAT,IVAL,NCOMAT,XCAR,
  3.      1  SIGF,VARF,DEFP,MFR,KERRE,ECOU,NECOU)
  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     DSIGT(NSTRS)  = INCREMENT DE CONTRAINTES TOTALES
  11. C     VAR0(NVARI)   =  VARIABLES INTERNES DEBUT
  12. C     VAR0(   1 )   =  EPSE
  13. C     VAR0(   2 )   =  DAM1
  14. C     VAR0(   3 )   =  DAM2
  15. C     VAR0(   4 )   =  DAM3
  16. C     VAR0(   5 )   =  ANGL
  17. C     VAR0(   6 )   =  LAM1
  18. C     VAR0(   7 )   =  LAM2
  19. C     VAR0(   8 )   =  LAM3
  20. C     VAR0(   9 )   =  LAM4
  21. C     VAR0(  10 )   =  LAM5
  22. C     XMAT(NCOMAT)  =  COMPOSANTES DE MATERIAU
  23. C     IVAL(NCOMAT)  =  INDICE DES COMPOSANTES DE MATERIAU
  24. C     NCOMAT        = NOMBRE DE COMPOSANTES DE MATERIAU
  25. C     XCAR(ICARA)   =  CARACTERISTIQUES
  26. C     MFR           =  NUMERO DE LA FORMULATION DE L'ELEMENT FINI
  27. C                   = 1       MASSIF
  28. C                   = 3       COQUE MINCE ( COQ2 , COQ3  ET DKT )
  29. C                   = 5       COQUE EPAISSE ( COQ6 , COQ8 )
  30. C                   = 7       POUTRE
  31. C                   = 9       COQUE MINCE AVEC CISAILLEMENT TRANSVERSE
  32. C                                    ( COQ4 ET DST )
  33. C                   = 11      LIQUIDE
  34. C                   = 13      TUYAU
  35. C                   = 15      LINESPRING
  36. C                   = 17      TUYAU FISSURE
  37. C                   = 19      RACCORD MASSIF
  38. C                   = 21      RACCORD COQUE
  39. C                   = 23      SURFACE LIBRE
  40. C                   = 25      MEMBRANE
  41. C                   = 27      UNIAXIALE
  42. C                   = 29      THERMIQUE
  43. C                   = 31      INCOMPRESSIBLES
  44. C                   = 33      POREUX
  45. C                   = 35      JOINT
  46. C                   = 37      HOMOGENEISE
  47. C                   = 39      TUYO
  48. C                   = 41      TUYAU ACOUSTIQUE PURE
  49. C                   = 43      RACCORD TUYAU FLUIDE
  50. C  SORTIES
  51. C     SIGF(NSTRS)   = CONTRAINTES FINALES
  52. C     VARF(NVARI)   = VARIABLES INTERNES FINALES
  53. C     DEFP(NSTRS)   = DEFORMATIONS PLASTIQUES
  54. C     KERRE         = 0    TOUT OK
  55. C        1   SI DLAMBDA NEGATIF
  56. C        2   NOMBRE MAX D ITERATIONS INTERNES DEPASSE
  57. C       21   ON NE TROUVE PAS L INTERSECTION AVEC LA SURFACE DE CHARGE
  58. C       22   SIG0 A L EXTERIEUR DE LA SURFACE DE CHARGE
  59. C
  60. C-----------------------------------------------------------------------
  61. C      VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU  ET NECOU
  62. C
  63. C  IFOUR    INDICE DU TYPE DE PROBLEME
  64. C            -3  DEFORMATIONS PLANES GENERALISEES
  65. C            -2  CONTRAINTES PLANES
  66. C            -1  DEFORMATIONS PLANES
  67. C             0  AXISYMETRIQUE
  68. C             1  SERIE DE FOURIER
  69. C             2  TRIDIMENSIONNEL
  70. C  ITYP      TYPE DE FORMULATION MECANIQUE
  71. C  --------------- ATTENTION ---------------
  72. C     IL EST ACTIF APRES L APPEL DE VISAVI
  73. C  -----------------------------------------
  74. C       ITYP=1   CAS DES ELEMENTS MASSIFS
  75. C       ITYP=2   CAS DES COQUES
  76. C       ITYP=3   CAS DES MEMBRANES
  77. C       ITYP=4   CAS DES CABLES ET DES BARRES
  78. C       ITYP=5   CAS QUELCONQUE
  79. C       ITYP=6   CAS DES CONTRAINTES PLANES
  80. C       ITYP=7   CAS DES COQUES A NU=0. OU CONTRAINTES PLANES
  81. C       ITYP=8   CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES
  82. C       ITYP=9   CAS DES COQUES EPAISSES
  83. C       ITYP=10  CAS DES JOINTS
  84. C       ITYP=11  CAS DES POUTRES
  85. C       ITYP=12  CAS DES TUYAUX
  86. C       ITYP=13  CAS DES COQUES AVEC CISAILLEMENT TRANSVERSE
  87. C
  88. C-----------------------------------------------------------------------
  89. C
  90.       IMPLICIT INTEGER(I-N)
  91.       IMPLICIT REAL*8(A-H,O-Z)
  92.  
  93. -INC PPARAM
  94. -INC CCOPTIO
  95.  
  96. -INC TECOU
  97.  
  98.       DIMENSION SIG0(*),DSIGT(*),VAR0(*),XMAT(*),XCAR(*),IVAL(*)
  99.       DIMENSION SIGF(*),VARF(*),DEFP(*)
  100.  
  101.       DIMENSION SIGMA(6),DSIGMA(6),SIGFIN(6)
  102.       DIMENSION XLAMBD(5),DEFPLA(6)
  103.       DIMENSION SPHER(6),AUXIL(6),DSIGZE(6),ORMAT(1)
  104.       CHARACTER*8 CMATE
  105.  
  106.       CMATE='ISOTROPE'
  107.       KERRE=0
  108. C
  109. C TEST DE CONSISTANCE DES DONNEES
  110. C
  111.       DO 100 I=1,NCOMAT
  112.         IF (XMAT(I).LT.0.D0) KERRE=437
  113.   100 CONTINUE
  114.       IF (KERRE.EQ.437) RETURN
  115. C
  116.       YOUN= 0.D0
  117.       XNU = 0.D0
  118.       RHO = 0.D0
  119.       ALPH= 0.D0
  120.       CTR1= 0.D0
  121.       CTR2= 0.D0
  122.       CTR3= 0.D0
  123.       ANGL= 0.D0
  124.       CLCS= 0.D0
  125.       EPCS= 0.D0
  126.       CLCT= 0.D0
  127.       PCT = 0.D0
  128.       EPCT= 0.D0
  129.       CLCD= 0.D0
  130.       PCD = 0.D0
  131.       EPCD= 0.D0
  132.       CLPH= 0.D0
  133.       PENT= 0.D0
  134. C
  135.       YOUN= XMAT(1)
  136.       XNU = XMAT(2)
  137.       RHO = XMAT(3)
  138.       ALPH= XMAT(4)
  139.       CTR1= XMAT(5)
  140.       CTR2= XMAT(6)
  141.       CTR3= XMAT(7)
  142.       ANGL= XMAT(8)
  143.       CLCS= XMAT(9)
  144.       EPCS= XMAT(10)
  145.       CLCT= XMAT(11)
  146.       PCT = XMAT(12)
  147.       CLCD= XMAT(13)
  148.       PCD = XMAT(14)
  149.       EPCD= XMAT(15)
  150.       CLPH= XMAT(16)
  151.       PENT= XMAT(17)
  152. C
  153.       IF(IIMPI.EQ.9) WRITE(IOIMP,2007) (XMAT(I),I=1,17)
  154.       IF(IIMPI.EQ.9) WRITE(IOIMP,2008) (IVAL(I),I=1,17)
  155. C
  156. C  PETIT TEST SUR NU POUR CERTAINS CAS
  157. C
  158.       IF (MFR.EQ.1.AND.(IFOUR.GE.-1.OR.IFOUR.EQ.-3)) THEN
  159. C
  160. C VALEURS PAR DEFAUT
  161. C
  162. C                  DONNEES GENERALES
  163. C
  164.         IF(IVAL( 3).EQ.0) RHO =2.4D+3
  165.         IF(IVAL( 4).EQ.0) ALPH=1.2D-5
  166. C
  167. C                  POUR LA TRACTION
  168. C
  169.         IF(IVAL( 5).EQ.0) CTR1=YOUN*1.2E-4
  170.         IF(IVAL( 6).EQ.0) CTR2=CTR1
  171.         IF(IVAL( 7).EQ.0) CTR3=CTR1
  172. C
  173. C                  POUR LA COMPRESSION
  174. C
  175.         IF(IVAL( 9).EQ.0) CLCS=YOUN*1.2D-3
  176.         IF(IVAL(10).EQ.0) EPCS=10.D0*CLCS/YOUN
  177.         IF(IVAL(12).EQ.0) PCT=1.6D0*CLCS
  178. C
  179. C                  BI-COMPRESSION
  180. C
  181.         REP=EPCS*YOUN/CLCS
  182.         AAA=CLCT/CLCS
  183.         BBB=PCT/CLCS
  184.         CCC=CLCT/PCT
  185.         UMN=1.D0-XNU
  186.         ALF=((REP-UMN-(SQRT(0.5D0*REP*UMN)))/(2.D0*REP-UMN))*0.95D0
  187.         DDD=(1.D0-ALF)/(1.D0-2.D0*ALF)
  188.         CLBI=DDD*CLCS
  189. C
  190. C                CISAILLEMENT FRAGILE
  191. C
  192.         IF(IVAL(11).EQ.0) THEN
  193.           CCC=((1.D0+BBB)/BBB)+((3.D0*ALF)/(1.D0-ALF))
  194.           CLCT=CCC*PCT
  195.         ELSE IF(IVAL(12).EQ.0) THEN
  196.           BBB=((AAA-1.D0)*(1.D0-ALF))/(1.D0+2.D0*ALF)
  197.           PCT=BBB*CLCS
  198.         ENDIF
  199. C
  200. C                CISAILLEMENT DUCTILE
  201. C
  202.         IF(IVAL(14).EQ.0) PCD=1.125D0*PCT
  203.         CCC=CLCT/PCT
  204.         EEE=PCD/PCT
  205.         FFF=(((1.D0+2.D0*ALF)*(EEE-1.D0))/(EEE*(1.D0-ALF)))+(CCC/EEE)
  206.         FFF=FFF/1.05D0
  207.         IF(IVAL(13).EQ.0) CLCD=FFF*PCD
  208.         IF(IVAL(15).EQ.0) EPCD=5.D0*EPCS
  209. C
  210. C                PRESSION HYDROSTATIQUE
  211. C
  212.         BBB=PCT/CLCS
  213.         IF(IVAL(16).EQ.0) CLPH=(((CCC+2.D0)*BBB)/3.D0)*1.05D0*CLCS
  214.         IF(IVAL(17).EQ.0) PENT=YOUN/(9.D0*(1.D0-2.D0*XNU))
  215. C
  216. C INITIALISATION
  217. C
  218.         JEPRIN = 0
  219.         PREC=1.D-2
  220.         RFPR=1.D-3
  221.         RFEP=1.D-7
  222.         RFSG=YOUN*RFEP
  223. C
  224. C CALCUL ET TEST DES PARAMETRES DU BETON
  225. C
  226. C     CALCUL DES PENTES
  227. C
  228.         XK1=YOUN/(3.D0*(1.D0-2.D0*XNU))
  229.         XK2=PENT
  230.         EPSA=CLCS/YOUN
  231.         EPSB=(CLCT-PCT)/YOUN
  232.         EPSC=(CLCD-PCD)/YOUN
  233.         PENTE1=(CLCD-PCD-CLCT+PCT)/(EPCD-EPSB)
  234.         PENTE1=(YOUN*PENTE1)/(YOUN-PENTE1)
  235.         PENTE2=-CLCS/(EPCS-EPSA)
  236.         PENTE2=(YOUN*PENTE2)/(YOUN-PENTE2)
  237.         PENTE3=XK1*XK2/(XK1-XK2)
  238. C
  239.         CLCS=-CLCS
  240.         CLCSA=CLCS*RFEP
  241.         CLCT=-CLCT
  242.         CLCD=-CLCD
  243.         PCT=-PCT
  244.         PCD=-PCD
  245.         EPSA=-EPSA
  246.         EPSB=-EPSB
  247.         EPSC=-EPSC
  248.         EPCD=-EPCD
  249.         EPCS=-EPCS
  250. C
  251. C     CALCUL DES CONSTANTES
  252. C
  253.         XI1A=CLCS
  254.         XI1AA=CLCSA
  255.         XI1B=CLCT+2.D0*PCT
  256.         XI1C=CLCD+2.D0*PCD
  257.         SEQA=ABS(CLCS)
  258.         SEQAA=ABS(CLCSA)
  259.         SEQB=ABS(CLCT-PCT)
  260.         SEQC=ABS(CLCD-PCD)
  261. C
  262.         VMELAS=SEQB
  263. C
  264.         ALFADV=(SEQC-SEQB)/(XI1B-XI1C)
  265.         DPELAS=ALFADV*XI1B+SEQB
  266. C
  267.         ALFAD1=(SEQB-SEQAA)/(XI1AA-XI1B)
  268.         DPELA1=ALFAD1*XI1B+SEQB
  269. C
  270.         ALFAD2=(SEQB-SEQA)/(XI1A-XI1B)
  271.         DPELA2=ALFAD2*XI1B+SEQB
  272.         PENTE2=PENTE2*(ALFAD2-1.D0)*(ALFAD2-1.D0)
  273. C
  274.         PORELA=CLPH
  275. C
  276.         UMU=YOUN/(1.D0+XNU)/2.D0
  277.         ALAM=2.D0*UMU*XNU/(1.-2.D0*XNU)
  278.         CPLA=YOUN/(1.D0-(XNU*XNU))
  279. C
  280.         TRTRAC=CTR1+CTR2+CTR3
  281.         XLAMER=(ALFAD2*TRTRAC-DPELA2)/PENTE2
  282.         XLAMAX=(ALFAD2*DPELA1-ALFAD1*DPELA2)/(ALFAD1*PENTE2)
  283.         DP2MIN=DPELA1*ALFAD2/ALFAD1
  284. C
  285.         SIGTRA=TRTRAC*0.333333333333333333333333333333333D0
  286.         PMIN=-XI1B*0.333333333333333333333333333333333D0
  287.         SIG1=(DPELA2-(2.D0*ALFAD2+1.D0)*PCD)/(ALFAD2-1.D0)
  288.         SIG3=(DPELA2-(ALFAD2-1.D0)*CLCD)/(2.D0*ALFAD2+1.D0)
  289. C
  290. C     TEST SUR LES VALEURS DES PARAMETRES
  291. C
  292.         IF(CLCT.GT.CLCS.OR.CLCT.LT.CLCD) THEN
  293.           KERRE=640
  294.           WRITE(IOIMP,5010) ABS(CLCT),ABS(CLCS),ABS(CLCD)
  295.         ENDIF
  296. C
  297.         IF(PCT.LT.PCD) THEN
  298.           KERRE=640
  299.           WRITE(IOIMP,5011) ABS(PCT),ABS(PCD)
  300.         ENDIF
  301. C
  302.         IF(EPCS.GE.EPSA) THEN
  303.           KERRE=640
  304.           WRITE(IOIMP,5012) ABS(EPCS),ABS(EPSA)
  305.         ENDIF
  306. C
  307.         IF(EPCD.GE.EPSC) THEN
  308.           KERRE=640
  309.         WRITE(IOIMP,5013)ABS(EPCD),ABS(EPSC)
  310.         ENDIF
  311. C
  312.         IF(ALFADV.LT.0.D0) THEN
  313.           KERRE=640
  314.           WRITE(IOIMP,5017)
  315.         ENDIF
  316. C
  317.         IF(ALFAD1.LT.0.D0.OR.ALFAD2.LT.0.D0) THEN
  318.           KERRE=640
  319.           WRITE(IOIMP,5018)
  320.         ENDIF
  321. C
  322.         IF(ABS(CLCS).LT.SIGTRA) THEN
  323.           KERRE=640
  324.           WRITE(IOIMP,5014) ABS(CLCS),SIGTRA
  325.         ENDIF
  326. C
  327.         IF(PMIN.GE.PORELA) THEN
  328.           KERRE=640
  329.           WRITE(IOIMP,5015) PMIN
  330.         ENDIF
  331. C
  332.         IF(ALFAD1.GE.0.5D0) THEN
  333.           KERRE=640
  334.           WRITE(IOIMP,5019) ABS(CLCT),(4*ABS(PCT))
  335.          ENDIF
  336. C
  337.         IF(DPELA2.GE.DPELAS) THEN
  338.           KERRE=640
  339.           WRITE(IOIMP,5020) ABS(CLCD),ABS(SIG1),ABS(PCD),ABS(SIG3)
  340.         ENDIF
  341. C
  342.         IF(XK2.GE.XK1) THEN
  343.           KERRE=640
  344.           WRITE(IOIMP,5021) XK2,XK1
  345.         ENDIF
  346. C
  347.         ECR2=(2.D0*YOUN*(0.5D0-ALFAD2)*(0.5D0-ALFAD2))/(1.D0-XNU)
  348.         IF(ABS(PENTE2).GT.ECR2) THEN
  349.           KERRE=640
  350.           WRITE(IOIMP,5016) ABS(PENTE2),ECR2
  351.         ENDIF
  352. C
  353.         ECR3=(1.5D0/(1.D0+XNU))+(3.D0*ALFAD2*ALFAD2/(1.D0-2.D0*XNU))
  354.         ECR3=ECR3*YOUN
  355.         IF(ABS(PENTE2).GT.ECR3) THEN
  356.           KERRE=640
  357.           WRITE(IOIMP,5022) ABS(PENTE2),ECR3
  358.         ENDIF
  359. C
  360.         IF(KERRE.NE.0) RETURN
  361. C
  362. C ECOULEMENT PLASTIQUE
  363. C
  364.   333   MCOD=1
  365.         CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  366.      .          SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  367.      .        DSIGZE,ICENT2,MCOD,IBOU,MFR,NSTRS,XCAR,CMATE,ecou,necou)
  368.         IF(necou.ITYP.EQ.0) THEN
  369.           KERRE=269
  370.           RETURN
  371.         ENDIF
  372. C
  373. C IMPRESSION DES PARAMETRES DU MODELE
  374. C
  375.         IF(IIMPI.EQ.9) THEN
  376.           WRITE(IOIMP,1000)
  377.           WRITE(IOIMP,1001) YOUN
  378.           WRITE(IOIMP,1002) XNU
  379.           WRITE(IOIMP,1003) RHO
  380.           WRITE(IOIMP,1004) ALPH
  381.           WRITE(IOIMP,1005)
  382.           WRITE(IOIMP,1006) CTR1,CTR2,CTR3
  383.           WRITE(IOIMP,1007) ANGL
  384.           WRITE(IOIMP,1008)
  385.           WRITE(IOIMP,1009) ABS(CLCS),ABS(EPCS)
  386.           WRITE(IOIMP,1010) ABS(CLBI)
  387.           WRITE(IOIMP,1011)
  388.           WRITE(IOIMP,1012) ABS(CLCT),ABS(PCT)
  389.           WRITE(IOIMP,1013) ABS(CLCD),ABS(PCD),ABS(EPCD)
  390.           WRITE(IOIMP,1014)
  391.           WRITE(IOIMP,1015) CLPH,PENT
  392.         ENDIF
  393. C
  394. C INITIALISATIONS DES PARAMETRES
  395. C
  396.         EPSE=VAR0(1)
  397.         IF(CTR1.GT.0.D0) RT1=(1.D0-VAR0(2))*CTR1
  398.         IF(CTR1.LE.0.D0) RT1=0.D0
  399.         IF(CTR2.GT.0.D0) RT2=(1.D0-VAR0(3))*CTR2
  400.         IF(CTR2.LE.0.D0) RT2=0.D0
  401.         IF(CTR3.GT.0.D0) RT3=(1.D0-VAR0(4))*CTR2
  402.         IF(CTR3.LE.0.D0) RT3=0.D0
  403.         ANGLE=VAR0(5)
  404.         IF(IVAL( 8).NE.0) ANGLE=ANGL
  405.         XLAMBD(1)=VAR0(6)
  406.         XLAMBD(2)=VAR0(7)
  407.         XLAMBD(3)=VAR0(8)
  408.         XLAMBD(4)=VAR0(9)
  409.         XLAMBD(5)=VAR0(10)
  410. C
  411.         IFIS=0
  412.         IF((RT1.EQ.0.D0.OR.RT2.EQ.0.D0).AND.ANGLE.NE.0.D0) IFIS=1
  413.         IF(IVAL( 8).NE.0) IFIS=1
  414. C
  415.         IDED=0
  416.         IF(XLAMBD(3).GT.0.D0) IDED=IDED+1
  417.         IF(XLAMBD(5).GT.0.D0) IDED=IDED+2
  418. C
  419.         DO 200 IB=1,IBOU
  420.         STOT(IB)=SIGMA(IB)+DSIGMA(IB)
  421.   200   CONTINUE
  422. C
  423. C ECOULEMENT PLASTIQUE
  424. C
  425.         IF(IIMPI.EQ.9) THEN
  426.           WRITE(IOIMP,2000)
  427.           WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  428.           WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  429.           WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  430.           WRITE(IOIMP,2006) (XLAMBD(I),I=1,5)
  431.           WRITE(IOIMP,2010) RT1,RT2,RT3
  432.         ENDIF
  433. C
  434.         CALL DAMAGE(ALFADV,ALFAD1,ALFAD2,VMELAS,DPELAS,
  435.      .              DPELA1,DPELA2,PORELA,RT1,RT2,RT3,
  436.      .              PENTE1,PENTE2,PENTE3,XLAMBD,SIGMA,DSIGMA,STOT,
  437.      .              SIGFIN,ANGLE,YOUN,XNU,XLAMER,XLAMAX,DP2MIN,
  438.      .              PREC,RFSG,RFEP,RFPR,IFIS,IDED,IDAM,KERRE)
  439.  
  440. C
  441.         IF(IIMPI.EQ.9) WRITE(IOIMP,2011) IDAM
  442. C
  443.         IF(KERRE.NE.0) THEN
  444.            IF(JEPRIN.EQ.0) THEN
  445.              KERRE=0
  446.              JEPRIN=1
  447.              IIMPI=9
  448.              GO TO 333
  449.            ELSE
  450.              IIMPI=0
  451.              RETURN
  452.            ENDIF
  453.         ENDIF
  454. C
  455. C PAS DE PLASTICITE COMPORTEMENT ELASTIQUE
  456. C
  457.         IF(IDAM.EQ.1) THEN
  458.           DO 300 IB=1,IBOU
  459.           SIGFIN(IB)=STOT(IB)
  460.           DEFPLA(IB)=0.D0
  461.   300     CONTINUE
  462.         DEPSEQ=0.D0
  463. C
  464.         ELSE
  465. C
  466. C NOUS AVONS ENDOMMAGE LE MATERIAU
  467. C
  468. C       CALCUL DE DEFORMATION PLASTIQUE
  469. C
  470.         DO 400 IB=1,IBOU
  471.         DSIGP(IB)=STOT(IB)-SIGFIN(IB)
  472.   400   CONTINUE
  473.         CALL EPSIG(DSIGP,DEFPLA,IFOUR,YOUN,XNU,necou.ITYP,ORMAT,XCAR)
  474.         DEPSEQ=SQRT(AVM(DEFPLA,DEFPLA))
  475. C
  476.         ENDIF
  477. C
  478.         IF(IIMPI.EQ.9) THEN
  479.           WRITE(IOIMP,2000)
  480.           WRITE(IOIMP,2001) (SIGMA(I),I=1,IBOU)
  481.           WRITE(IOIMP,2002) (DSIGMA(I),I=1,IBOU)
  482.           WRITE(IOIMP,2003) (STOT(I),I=1,IBOU)
  483.           WRITE(IOIMP,2004) (SIGFIN(I),I=1,IBOU)
  484.           WRITE(IOIMP,2005) (DEFPLA(I),I=1,IBOU)
  485.           WRITE(IOIMP,2006) (XLAMBD(I),I=1,5)
  486.         ENDIF
  487. C
  488. C MISE A JOUR DES VARIABLES INTERNES
  489. C
  490.         VARF(1)=EPSE+DEPSEQ
  491.         IF(CTR1.GT.0.D0) VARF(2)=1.D0-RT1/CTR1
  492.         IF(CTR1.LE.0.D0) VARF(2)=1.D0
  493.         IF(CTR2.GT.0.D0) VARF(3)=1.D0-RT2/CTR2
  494.         IF(CTR2.LE.0.D0) VARF(3)=1.D0
  495.         IF(CTR3.GT.0.D0) VARF(4)=1.D0-RT3/CTR3
  496.         IF(CTR3.LE.0.D0) VARF(4)=1.D0
  497.         VARF(5)=ANGLE
  498.         VARF(6)=XLAMBD(1)
  499.         VARF(7)=XLAMBD(2)
  500.         VARF(8)=XLAMBD(3)
  501.         VARF(9)=XLAMBD(4)
  502.         VARF(10)=XLAMBD(5)
  503. C
  504.         IF(IIMPI.EQ.9) WRITE(IOIMP,2009) EPSE,DEPSEQ,VARF(1)
  505. C
  506.          MCOD=2
  507.         CALL VISAVI(SIG0,DSIGT,VAR0,SIGMA,DSIGMA,SPHER,AUXIL,
  508.      .          SIGF,DEFP,VARF,SIGFIN,DEFPLA,
  509.      .        DSIGZE,ICENT2,MCOD,IBOU,MFR,NSTRS,XCAR,CMATE,ecou,necou)
  510. C
  511.       ELSE
  512.         KERRE=639
  513.         RETURN
  514.       ENDIF
  515. C
  516. C LES FORMATS D IMPRESSION
  517. C
  518.  1000 FORMAT(//,31X,'MATERIAU BETON',//)
  519.  1001 FORMAT(6X,'MODULE D YOUNG  ',T51,'= ',1PD12.5)
  520.  1002 FORMAT(6X,'COEF. DE POISSON',T51,'= ',1PD12.5)
  521.  1003 FORMAT(6X,'MASSE VOLUMIQUE ',T51,'= ',1PD12.5)
  522.  1004 FORMAT(6X,'COEF. DE DILATATION THERMIQUE',T51,'= ',1PD12.5)
  523.  1005 FORMAT(/,1X,'DOMMAGE PAR TRACTION :',/)
  524.  1006 FORMAT(6X,'CONTRAINTES LIMITES :',/
  525.      *     11X,'DIRECTION(1) = ',1PD12.5,
  526.      *      5X,'DIRECTION(2) = ',1PD12.5,5X,'DIRECTION(3) = ',1PD12.5/)
  527.  1007 FORMAT(6X,'ANGLE DE FISSURATION',T51,'= ',1PD12.5/)
  528.  1008 FORMAT(/,1X,'DOMMAGE PAR COMPRESSION :',/)
  529.  1009 FORMAT(6X,'COMPRESSION SIMPLE :',/
  530.      *     11X,'CONTRAINTE LIMITE = ',1PD12.5,
  531.      *      5X,'DEFORMATION A RUPTURE   = ',1PD12.5/)
  532.  1010 FORMAT(6X,'COMPRESSION BIAXIALE :',/
  533.      *     11X,'CONTRAINTE LIMITE = ',1PD12.5/)
  534.  1011 FORMAT(/,1X,'DOMMAGE PAR CISAILLEMENT :',/)
  535.  1012 FORMAT(6X,'SOUS CONFINEMENT DE TRANSITION :',/
  536.      *     11X,'CONTRAINTE LIMITE = ',1PD12.5,
  537.      *      5X,'PRESSION DE CONFINEMENT = ',1PD12.5/)
  538.  1013 FORMAT(6X,'SOUS CONFINEMENT DUCTILE :',/
  539.      *     11X,'CONTRAINTE LIMITE = ',1PD12.5,
  540.      *      5X,'PRESSION DE CONFINEMENT = ',1PD12.5,/
  541.      *     11X,'DEFORMATION AU DEBUT DU PALIER PLASTIQUE = ',1PD12.5/)
  542.  1014 FORMAT(/,1X,'DOMMAGE PAR PRESSION HYDROSTATIQUE :',/)
  543.  1015 FORMAT(6X,'PRESSION LIMITE = ',1PD12.5,
  544.      *     12X,'MODULE TANGENT = ',1PD12.5,//)
  545. C
  546.  2000 FORMAT(1X,'RESULTATS DE L ECOULEMENT PLASTIQUE')
  547.  2001 FORMAT(1X,'SIGMA  =',6(1X,1PD12.5))
  548.  2002 FORMAT(1X,'DSIGMA =',6(1X,1PD12.5))
  549.  2003 FORMAT(1X,'STOT   =',6(1X,1PD12.5))
  550.  2004 FORMAT(1X,'SIGFIN =',6(1X,1PD12.5))
  551.  2005 FORMAT(1X,'DEFPLA =',6(1X,1PD12.5))
  552.  2006 FORMAT(1X,'XLAMBD =',5(1X,1PD12.5))
  553.  2007 FORMAT(7(6(1X,1PD12.5),/))
  554.  2008 FORMAT(7(6(1X,I12),/))
  555.  2009 FORMAT(1X,'EPSE0  =',1X,1PD12.5,1X,'DEPSEQ=',1X,1PD12.5,
  556.      .       1X,'EPSE   =',1X,1PD12.5)
  557.  2010 FORMAT(1X,'RT1    =',1X,1PD12.5,1X,'RT2   =',1X,1PD12.5,
  558.      .       1X,'RT3    =',1X,1PD12.5)
  559.  2011 FORMAT(1X,'INDICE D ENDOMMAGEMENT IDAM  =',I4)
  560. C
  561.  5010 FORMAT(1X,'ERREUR LA LIMITE LCT =',1PD12.5,' DOIT ETRE COMPRISE',
  562.      .     /,1X,'ENTRE LCS =',1PD12.5,' ET LCD =',1PD12.5)
  563.  5011 FORMAT(1X,'ERREUR LE CONFINEMENT PCT =',1PD12.5,' DOIT ETRE ',/,
  564.      .       1X,'INFERIEUR A PCD =',1PD12.5)
  565.  5012 FORMAT(1X,'ERREUR LA DEFORMATION A RUPTURE EN COMPRESSION SIMPLE'
  566.      .         ,' EST MAUVAISE',/,
  567.      .       1X,'EPCS =',1PD12.5,' DOIT ETRE SUPERIEURE A ',1PD12.5)
  568.  5013 FORMAT(1X,'ERREUR LA DEFORMATION A RUPTURE EN COMPRESSION DUCTILE'
  569.      .         ,' EST MAUVAISE',/,
  570.      .       1X,'EPCD =',1PD12.5,' DOIT ETRE SUPERIEURE A ',1PD12.5)
  571.  5014 FORMAT(1X,'ERREUR LA LIMITE EN COMPRESSION SIMPLE LCS =',1PD12.5,
  572.      .     /,1X,'DOIT ETRE SUPERIREURE A LA LIMITE EN TRACTION LCT =',
  573.      .           D12.5)
  574.  5015 FORMAT(1X,'ERREUR LA LIMITE EN PRESSION HYDROSTATIQUE SE TROUVE',
  575.      .          ' DANS LE DOMAINE FRAGILE',/,
  576.      .       1X,'IL FAUT AUGMENTER LPH POUR QU''ELLE SOIT SUPERIEURE ',
  577.      .          'A PMIN =',1PD12.5)
  578.  5016 FORMAT(1X,'ERREUR MODULE D ECROUISSAGE DU DRUCKER EST TROP GRAND',
  579.      .          ' =',1PD12.5,/,
  580.      .       1X,'PENTE2 < 2 * E * (0.5 - ALFAD2 ) ** 2 / (1-NU) ',
  581.      .          'DONC PENTE2 MAX VAUT =',1PD12.5)
  582.  5017 FORMAT(1X,'ERREUR LA VALEUR DE ( LCD - PCD ) DOIT ETRE ',
  583.      .          'SUPERIEURE A ( LCT - PCT )')
  584.  5018 FORMAT(1X,'ERREUR LA VALEUR DE ( LCT - PCT ) DOIT ETRE ',
  585.      .          'SUPERIEURE A LCS')
  586.  5019 FORMAT(1X,'ERREUR LA LIMITE LCT =',1PD12.5,'DOIT ETRE INFERIEURE',
  587.      .          ' A',1PD12.5)
  588.  5020 FORMAT(1X,'ERREUR LA LIMITE LCD =',1PD12.5,'DOIT ETRE INFERIEURE',
  589.      .          ' A',1PD12.5,'   OU BIEN  :',/,
  590.      .       1X,'LE CONFINEMENT PCD =',1PD12.5,' DOIT ETRE SUPERIEURE',
  591.      .          ' A',1PD12.5)
  592.  5021 FORMAT(1X,'ERREUR LA PENTE DE LA PRESSION HYDROSTATIQUE ',
  593.      .          'EST MAUVAISE',/,
  594.      .       1X,'PENT =',1PD12.5,' DOIT ETRE INFERIEURE A ',1PD12.5)
  595.  5022 FORMAT(1X,'ERREUR MODULE D ECROUISSAGE DU DRUCKER EST TROP GRAND',
  596.      .           ' =',1PD12.5,/,
  597.      .       1X,'PENTE2 < 1.5*E/(1+NU) + 3*ALFAD2**2*E/(1-2*NU) ',
  598.      .          'DONC PENTE2 MAX VAUT ',1PD12.5)
  599. C
  600.       RETURN
  601.       END
  602.  
  603.  
  604.  

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