Télécharger la plaquette Cast3M


fr en

Télécharger @LISPA16.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @LISPA16  PROCEDUR  ALANIER   95/06/13    21:33:48     1710           
  2. 'DEBP' @LISPA16 CHARG*TABLE DON*TABLE ;
  3.  
  4. 'MESS' '                                               ' ;
  5. 'MESS' 'PROCEDURE DE CALCUL EN FATIGUE/FLUAGE SELON A16' ;
  6. 'MESS' '                                               ' ;
  7.  
  8. 'SI' (('EXISTE' DON REPRISE) 'EGA' VRAI) ;
  9.   TABCON = DON.'EVCON' ;
  10.   TABCAR = DON.'EVCAR' ;
  11.   TABNCY = DON.'EVNCY' ;
  12. 'SINON' ;
  13.   TABCAR = 'TABLE' ;
  14.   TABNCY = 'TABLE' ;
  15.   TABCON = 'TABLE' ;
  16.   TOTAL    = MOT 'TOTAL' ;
  17.   FATIGUE  = MOT 'FATIGUE' ;
  18.   FLUAGE   = MOT 'FLUAGE' ;
  19.   TABCAR.'TOTAL'   = 'TABLE' ;
  20.   TABCAR.'FATIGUE' = 'TABLE' ;
  21.   TABCAR.'FLUAGE'  = 'TABLE' ;
  22.   DELTSIG  = MOT 'DELTSIG' ;
  23.   CHSIG_A  = MOT 'CHSIG_A' ;
  24.   CHSIG_B  = MOT 'CHSIG_B' ;
  25.   DELTSIG1 = MOT 'DELTSIG1' ;
  26.   CHSIG_A1 = MOT 'CHSIG_A1' ;
  27.   CHSIG_B1 = MOT 'CHSIG_B1' ;
  28.   SDEFRAS  = MOT 'SDEFRAS'  ;
  29.   SDEFRBS  = MOT 'SDEFRBS'  ;
  30.   SDEFR    = MOT 'SDEFR' ;
  31.   VALR     = MOT 'VALR' ;
  32.   VALQ     = MOT 'VALQ' ;
  33.   VALPSI   = MOT 'VALPSI' ;
  34.   VALKDJ   = MOT 'VALKDJ' ;
  35.   VALDJ    = MOT 'VALDJ' ;
  36.   VALDKEFF = MOT 'VALDKEFF' ;
  37.   VALDAFAT = MOT 'VALDAFAT' ;
  38.   EPSILON  = MOT 'EPSILON' ;
  39.   TABCON.'DELTSIG'  = 'TABLE' ;
  40.   TABCON.'CHSIG_A'  = 'TABLE' ;
  41.   TABCON.'CHSIG_B'  = 'TABLE' ;
  42.   TABCON.'DELTSIG1' = 'TABLE' ;
  43.   TABCON.'CHSIG_A1' = 'TABLE' ;
  44.   TABCON.'CHSIG_B1' = 'TABLE' ;
  45.   TABCON.'SDEFRAS'  = 'TABLE' ;
  46.   TABCON.'SDEFRBS'  = 'TABLE' ;
  47.   TABCON.'SDEFR'    = 'TABLE' ;
  48.   TABCON.'VALR'     = 'TABLE' ;
  49.   TABCON.'VALQ'     = 'TABLE' ;
  50.   TABCON.'VALPSI'   = 'TABLE' ;
  51.   TABCON.'VALKDJ'   = 'TABLE' ;
  52.   TABCON.'VALDJ'    = 'TABLE' ;
  53.   TABCON.'VALDKEFF' = 'TABLE' ;
  54.   TABCON.'VALDAFAT' = 'TABLE' ;
  55.   TABCON.'EPSILON'  = 'TABLE' ;
  56. 'FINSI' ;
  57.  
  58. 'SI' (('EXISTE' DON FATUNI) 'EGA' VRAI) ;
  59.   'SI' ((DON.FATUNI) 'EGA' VRAI)  ;
  60.     CALFAT = 1 ;
  61.     'MESS' '(OPTION CALCUL EN FATIGUE UNIQUEMENT)' ;
  62.   'SINON' ;
  63.     CALFAT = 0 ;
  64.   'FINSI' ;
  65. 'SINON' ;
  66.   CALFAT = 0 ;
  67. 'FINSI' ;
  68.  
  69.  
  70. 'SI' (('EXISTE' DON A316) 'EGA' VRAI) ;
  71.   'SI' ((DON.A316) 'EGA' VRAI)  ;
  72.     A316 = 1 ;
  73.   'SINON' ;
  74.     A316 = 0 ;
  75.   'FINSI' ;
  76. 'SINON' ;
  77.   A316 = 0 ;
  78. 'FINSI' ;
  79.  
  80. 'SI' (('EXISTE' DON RELAX) 'EGA' VRAI) ;
  81.   'SI' ((DON.RELAX) 'EGA' VRAI ) ;
  82.     RELAX = 1 ;
  83.   'SINON' ;
  84.     RELAX = 0 ;
  85.   'FINSI' ;
  86. 'SINON' ;
  87.   RELAX = 0 ;
  88. 'FINSI' ;
  89.  
  90. 'SI' (('EXISTE' DON NREA) 'EGA' VRAI) ;
  91.   KREA = DON.'NREA' ;
  92. 'SINON' ;
  93.   KREA = 0 ;
  94. 'FINSI' ;
  95.  
  96. REA = DON.'NBREA';
  97. 'REPETER' BOUREA REA ;
  98.  
  99.  
  100.   'SI' (('EXISTE' DON AMAX) 'EGA' VRAI) ;
  101.     'SI' ((DON.'AMAX') 'EGA' VRAI) ;
  102.       'MESS' 'ARRET DU CALCUL' ;
  103.       'QUIT' BOUREA ;
  104.     'FINSI' ;
  105.   'FINSI' ;
  106.  
  107.   KREA = KREA + 1 ;
  108.  
  109.  
  110.  
  111. ****************************************************************        
  112. *               RESTITUTION DES DONNEES                        *
  113. ****************************************************************
  114.  
  115. ** -- MAILLAGES
  116.  
  117.   FISS1   = DON.'MFIS' ;
  118.   NBFI = 'NBEL' FISS1 ;
  119.   PREMA = FAUX ;
  120.  
  121. ** -- DELTA-CHARGEMENT POUR LA FATIGUE
  122.  
  123.   'SI' ('EXISTE' CHARG DP) ;
  124.           DCHA1 = CHARG.'DP' ;
  125.   'SINON' ;
  126.           CHARGI = CHARG.'DC' ;
  127.           DCHA1 = CHARGI * 1.E-6 ;
  128.   'FINSI' ;
  129.  
  130.  
  131.   DCHA  = CHARG.'DC' ;
  132.  
  133.  
  134. ****** CHARGEMENT POUR LE FLUAGE
  135.  
  136.  
  137.   CHAA = CHARG.'CA' ;
  138.  
  139.   CHAB = CHARG.'CB' ;
  140.  
  141.   'SI' ('EXISTE' CHARG PA) ;
  142.           CHAA1 = CHARG.'PA' ;
  143.   'SINON' ;
  144.           CHARGI = CHARG.'CA' ;
  145.           CHAA1 = CHARGI * 1.E-6 ;
  146.   'FINSI' ;
  147.  
  148.   'SI' ('EXISTE' CHARG PB) ;
  149.           CHAB1 = CHARG.'PB' ;
  150.   'SINON' ;
  151.           CHARGI = CHARG.'CB' ;
  152.           CHAB1 = CHARGI * 1.E-6 ;
  153.   'FINSI' ;
  154.  
  155.  
  156. ****** DIMENSIONS DE LA FISSURE
  157. ** -- MATERIAU
  158.   OBJF = DON.'OBJ' ;
  159.   MATF = DON.'MAT' ;
  160.   CARF = DON.'CAR' ;
  161.   MATFC = MATF 'ET' CARF ;
  162.   OBJ = DON.'OBJT' ;
  163.   MAT = DON.'MATT' ;
  164.   CAR = DON.'CART' ;
  165.   MAT = MAT 'ET' CAR ;
  166.  
  167.   'SI' ('NON' ('EXISTE' DON REPRISE)) ;
  168.  
  169. *************************************************************
  170. * RE-ECRITURE DES CHAMPS MODELE MATERIAU ET CARACTERISTIQUE *
  171. *************************************************************
  172. 'MESS' '                                                       ';
  173. 'MESS' '--> RE-ECRITURE DES CHAMPS MODELE MATERIAU ET
  174.  CARACTERISTIQUE';
  175.     KR = 0 ;KRF = 0. ;
  176.     'REPETER' BOUCR NBFI ;
  177.       KR = KR + 1 ;KRF = KRF + 1. ;
  178.       ELC = FISS1 'ELEM' KR ;
  179.       MOE = 'REDU' OBJF ELC ;
  180.       MONE = 'EXTR' MOE 'FORM' MECANIQUE ;
  181.       CHFIC = 'MANU' 'CHML' MONE 'SCAL' KRF 'TYPE' 'SCALAIRE'
  182.               'STRESSES' ;
  183.       MAE = 'REDU' MATF ELC ;
  184.       EY = 'EXTR' MAE 'YOUN' 1 1 2 ;
  185.       PO = 'EXTR' MAE 'NU' 1 1 2 ;
  186.       MANE = 'MATE' MONE 'YOUN' EY 'NU' PO ;
  187.       CAE = 'REDU' CARF ELC ;
  188.       XV = 'EXTR' CAE 'VX' 1 1 2 ;
  189.       YV = 'EXTR' CAE 'VY' 1 1 2 ; 
  190.       ZV = 'EXTR' CAE 'VZ' 1 1 2 ;
  191.       EP = 'EXTR' CAE 'EPAI' 1 1 2 ;
  192.       PROF = 'EXTR' CAE 'FISS' 1 1 2 ;
  193.       CARNE = 'CARA' MONE 'EPAI' EP 'FISS' PROF 
  194.              'VX' XV  'VY' YV  'VZ' ZV ;
  195.       CAR0E = 'CARA' MONE 'EPAI' EP 'FISS' 0. 
  196.              'VX' XV  'VY' YV  'VZ' ZV ;
  197.       'SI' (KR 'EGA' 1) ;
  198.         CAR1 = CARNE ;
  199.         MOD1 = MONE ;
  200.         MAT1 = MANE ;
  201.         CAR0 = CAR0E ;
  202.         CHFIC1 = CHFIC ;
  203.       'SINON';
  204.         CAR1 = CAR1 'ET' CARNE ;
  205.         MOD1 = MOD1 'ET' MONE ;
  206.         MAT1 = MAT1 'ET' MANE ;
  207.         CAR0 = CAR0 'ET' CAR0E ;
  208.         CHFIC1 = CHFIC1 'ET' CHFIC ;
  209.       'FINSI';
  210.     'FIN' BOUCR ;
  211.  
  212.     SELEC = 'TABLE' ;
  213.     KJ = 0 ;KJF = 0. ;
  214.     'REPETER' BOUCJ NBFI ;
  215.       KJ = KJ + 1 ;KJF = KJF + 1. ;
  216.       SELECS = 'MASQ' CHFIC1 'SUPERIEUR' (KJF + 0.1) ;
  217.       SELECI = 'MASQ' CHFIC1 'INFERIEUR' (KJF - 0.1) ;
  218.       SELECC = SELECS + SELECI ;
  219.       SELEC.KJ = 'MASQ' SELECC 'INFERIEUR' 0.5 ;
  220.     'FIN' BOUCJ ;
  221.  
  222.     OBJF = MOD1 ;
  223.     CARF = CAR1 ;
  224.     'SI' (CALFAT 'EGA' 0) ;
  225.       CARF1 = CAR1 ;CARF2 = CAR1 ;
  226.       TABCAR.'FATIGUE'. 0 = CAR1 ;
  227.       TABCAR.'FLUAGE'. 0 = CAR1 ;
  228.     'FINSI' ;
  229.     TABCAR.'TOTAL'. 0 = CAR1 ;
  230.     MATF = MAT1 ;
  231.     MATFC = MATF 'ET' CARF ;
  232.     MATF0 = MATF 'ET' CAR0 ;
  233.  
  234.     DON.'SEL' = SELEC ;
  235.  
  236.   'SINON' ;
  237.  
  238.     SELEC = DON.'SEL' ;
  239.     'SI' (CALFAT 'EGA' 0) ;
  240.       CARF1 = DON.'CAR1' ;CARF2 = DON.'CAR2' ;
  241.     'FINSI' ;
  242.   'FINSI' ;
  243.  
  244. ****************************************************************
  245. *                   CREATION DES RIGIDITES                     *
  246. ****************************************************************
  247.  
  248.   BLOCA = DON.'BLO' ;
  249.   OBJTF = OBJ 'ET' OBJF ;
  250.   MATTF = MAT 'ET' MATFC ;
  251.   'SI'  ('NON' ('EXISTE' DON REPRISE)) ;
  252.     MATT0 = MAT 'ET' MATF0 ;
  253.     RIGIK0 = 'RIGI' OBJTF MATT0 ;
  254.     RIGK0 = RIGIK0 'ET' BLOCA ;
  255.   'FINSI' ;
  256.   RIGIK = 'RIGI' OBJTF MATTF ;
  257.   RIGK  = RIGIK 'ET' BLOCA ;
  258.  
  259. ****************************************************************
  260. *                         RESOLUTIONS                          *
  261. ****************************************************************
  262.  
  263.   'SI'  ('NON' ('EXISTE' DON REPRISE)) ;
  264.   'MESS' '--> RESOLUTIONS DU CHARGEMENT NOMINALE' ;
  265.   'MESS' '                                      ' ;
  266.     D_RE0 = 'RESO' RIGK0 DCHA ;
  267. *    RE_A0 = 'RESO' RIGK0 CHAA ;
  268.     RE_B0 = 'RESO' RIGK0 CHAB ;
  269.     'SI' (('MAXI' ('ABS' (DCHA - DCHA1))) < 1.E-10) ;
  270.       D_RE01 = D_RE0 ;
  271.     'SINON' ;
  272.       D_RE01 = 'RESO' RIGK0 DCHA1 ;
  273.     'FINSI' ;
  274.     'SI' (('MAXI' ('ABS' (CHAA - CHAA1))) < 1.E-10) ;
  275. *      RE_A01 = RE_A0 ;
  276.     'SINON' ;
  277. *      RE_A01 = 'RESO' RIGK0 CHAA1 ;
  278.     'FINSI' ;
  279.     'SI' (('MAXI' ('ABS' (CHAB - CHAB1))) < 1.E-10) ;
  280.       RE_B01 = RE_B0 ;
  281.     'SINON' ;
  282.       RE_B01 = 'RESO' RIGK0 CHAB1 ;
  283.     'FINSI' ;
  284.   'FINSI' ;
  285.  
  286. *  'MESS' 'RESOLUTIONS DU CHARGEMENT' ;
  287.   D_RE = 'RESO' RIGK DCHA ;
  288.   RE_A = 'RESO' RIGK CHAA ;
  289.   RE_B = 'RESO' RIGK CHAB ;
  290.   'SI' (('MAXI' ('ABS' (DCHA - DCHA1))) < 1.E-10) ;
  291.     D_RE1 = D_RE ;
  292.   'SINON' ;
  293.     D_RE1 = 'RESO' RIGK DCHA1 ;
  294.   'FINSI' ;
  295.   'SI' (('MAXI' ('ABS' (CHAA - CHAA1))) < 1.E-10) ;
  296.     RE_A1 = RE_A ;
  297.   'SINON' ;
  298.     RE_A1 = 'RESO' RIGK CHAA1 ;
  299.   'FINSI' ;
  300.   'SI' (('MAXI' ('ABS' (CHAB - CHAB1))) < 1.E-10) ;
  301.     RE_B1 = RE_B ;
  302.   'SINON' ;
  303.     RE_B1 = 'RESO' RIGK CHAB1 ;
  304.   'FINSI' ;
  305.  
  306.  
  307.  
  308. ****************************************************************
  309. *             CALCULS DES CHAMPS DE CONTRAINTES                *
  310. ****************************************************************
  311.  
  312.   'SI' ('NON' ('EXISTE' DON REPRISE)) ;
  313.  
  314.     'OUBLIER' RIGIK0 ;
  315.     'OUBLIER' RIGK0 ;
  316.  
  317.     SIGNO   = 'SIGMA' OBJF MATF0 D_RE0 ;
  318. *    SIGANO  = 'SIGMA' OBJF MATF0 RE_A0 ;
  319.     SIGBNO  = 'SIGMA' OBJF MATF0 RE_B0 ;
  320.     SIGNO1  = 'SIGMA' OBJF MATF0 D_RE01 ;
  321. *    SIGANO1 = 'SIGMA' OBJF MATF0 RE_A01 ;
  322.     SIGBNO1 = 'SIGMA' OBJF MATF0 RE_B01 ;
  323.  
  324.     'OUBLIER' D_RE0 ;
  325.     'OUBLIER' RE_B0 ;
  326.     'OUBLIER' D_RE01 ;
  327.     'OUBLIER' RE_B01 ;
  328.  
  329.     DON.'DNOP' = SIGNO1 ;
  330.     DON.'DNO'  = SIGNO ;
  331. *    DON.'NO1'  = SIGANO1 ;
  332. *    DON.'NO'   = SIGANO ;
  333.     DON.'NOB1' = SIGBNO1 ;
  334.     DON.'NOB'  = SIGBNO ;
  335.  
  336.   'SINON' ;
  337.  
  338.     SIGNO1  = DON.'DNOP' ;
  339.     SIGNO   = DON.'DNO' ;
  340. *    SIGANO1 = DON.'NO1' ;
  341. *    SIGANO  = DON.'NO' ;
  342.     SIGBNO1 = DON.'NOB1' ;
  343.     SIGBNO  = DON.'NOB' ;
  344.  
  345.   'FINSI' ;
  346.  
  347.  
  348.   'OUBLIER' RIGIK ;
  349.   'OUBLIER' RIGK ;
  350.  
  351.   DSIG  = 'SIGMA' OBJF MATFC D_RE ;
  352.   SIGA  = 'SIGMA' OBJF MATFC RE_A ;
  353.   SIGB  = 'SIGMA' OBJF MATFC RE_B ;
  354.   DSIG1 = 'SIGMA' OBJF MATFC D_RE1 ;
  355.   SIGA1 = 'SIGMA' OBJF MATFC RE_A1 ;
  356.   SIGB1 = 'SIGMA' OBJF MATFC RE_B1 ;
  357.  
  358.   'OUBLIER' D_RE ;
  359.   'OUBLIER' D_RE1 ;
  360.   'OUBLIER' RE_A ;
  361.   'OUBLIER' RE_A1 ;
  362.   'OUBLIER' RE_B ;
  363.   'OUBLIER' RE_B1 ;
  364.  
  365. ** -- FISSURE
  366.  
  367.   A0 = CHAN 'STRESSES' OBJF MATFC ;
  368.   A = 'EXCO' A0 'FISS' SCAL ;
  369.   W = 'EXCO' A0 'EPAI' SCAL ;
  370.   E = 'EXCO' A0 'YOUN' SCAL ; 
  371.   ESTAR = E ; EY = 'MAXI' E ;
  372.  
  373.   'SI' (('EXISTE' DON ALIMIT) 'EGA' VRAI) ;
  374.     COEFIMP = DON.ALIMIT ;
  375.   'SINON' ;
  376.     COEFIMP = 0.7 ;
  377.   'FINSI' ;
  378.  
  379.   'SI' (('EXISTE' DON REACTUA) 'EGA' VRAI) ;
  380.     COEFREA = 1. + ((DON.REACTUA) / 100.)  ;
  381.   'SINON' ;
  382.     COEFREA = 1.1 ;
  383.   'FINSI' ;
  384.  
  385.   'SI' ((COEFIMP > 1.) 'OU' (COEFIMP < 0.)) ;
  386.     'MESS' 'LE CONTENU DE DON.ALIMIT DOIT ETRE COMPRIS
  387.  ENTRE 0 ET 1' ;
  388.   'FINSI' ;
  389.  
  390.   'SI' ('EXISTE' DON ATOT) ;
  391.     ADEB = DON.ATOT ;
  392.   'SINON' ;
  393.     ADEB = A ;
  394.   'FINSI' ;
  395.  
  396.   'SI' (CALFAT 'EGA' 0) ;
  397.     A1 = CHAN 'STRESSES' OBJF CARF1 ;
  398.     A2 = CHAN 'STRESSES' OBJF CARF2 ;
  399.     'SI' (('EXISTE' DON ATOT1) 'ET' ('EXISTE' DON ATOT2)) ;
  400.       ADEB1 = DON.ATOT1 ;ADEB2 = DON.ATOT2 ;
  401.     'SINON' ;
  402.       ADEB1 = 'EXCO' A1 'FISS' SCAL ;
  403.       ADEB2 = 'EXCO' A2 'FISS' SCAL ;
  404.     'FINSI' ;
  405.   'FINSI' ;
  406.  
  407.   AFIN = ADEB * COEFREA ;
  408.   WIMP = W * COEFIMP ;
  409.  
  410.  
  411.  
  412.  
  413. ****** LOIS DE COMPORTEMENT
  414.  
  415.  
  416.  
  417. ****** LOIS DE PROPAGATION
  418.  
  419.   'SI' (A316 'EGA' 1) ;
  420. *   LOI DE PROPAGATION EN FATIGUE
  421.     CFAT = 1.4E-10 ; NFAT = 3. ;
  422. *   LOI DE PROPAGATION EN FLUAGE
  423.     CFLU = 2.17E-2 ; NFLU = 0.81 ;
  424. *   LOI DE COMPORTEMENT 
  425. *   COURBE DE TRACTION ( A 525 °C ) 
  426.     PROGTE = 'PROG' 0. 0.00066 0.002 0.01  0.05  0.31 ; 
  427.     PROGTS = 'PROG' 0. 100.    135.  164.5 283.5 619. ;
  428. *   COURBE CYCLIQUE 
  429.     PROGCS = 'PROG' 0. PAS 10. 800. ;
  430.     PROGCE =  (100. * (PROGCS * ((2. * (1. + 0.3)) / (3. * EY))))
  431.     + ((PROGCS / 718.) ** (1. / 0.319)) ;
  432.     PROGCE = PROGCE * (1. / 100.) ;
  433. *   MODULE D'YOUNG CYCLIQUE
  434.     ECYCL = (3. * EY) / (2. * ( 1. + 0.3)) ;
  435. *   LIMITE D'ELASTICITE
  436.     R002 = 114. ;
  437. *   COEFFICIENT DE FLUAGE
  438.     C1 = 1.8803E-12;  C2 = 0.40053 ;  C3 = 6.1073E23 ;
  439.     N1 = 4.125 ;      N3 = -8.2276 ;
  440.     KFP = 2.1218E-30 ;  XFP = -1.4967 ;   YFP = 10.2217 ;
  441.     CFS = 4.1482E-29 ;  NFS = 9.06 ;
  442.   'SINON' ;
  443. *   LOI DE PROPAGATION EN FATIGUE
  444.     CFAT = DON.FAT.'COEFFICIENT' ; NFAT = DON.FAT.'EXPOSANT' ;
  445. *   COURBE CYCLIQUE
  446.     PROGCE = DON.CC.'EPS'        ;PROGCS = DON.CC.'SIG' ;
  447. *   MODULE D'YOUNG CYCLIQUE
  448.     ECYCL = ('EXTR' PROGCS 1) / ('EXTR' PROGCE 1) ;
  449. *   COURBE DE TRACTION OU CYCLIQUE REDUITE
  450.     PROGTE = DON.CT.'EPS'      ; PROGTS = DON.CT.'SIG' ;
  451. *   LIMITE D'ELASTICITE
  452.     R002 = DON.'R002' ;
  453.     'SI' (CALFAT 'EGA' 0) ;
  454. *     LOI DE PROPAGATION EN FLUAGE
  455.       CFLU = DON.FLU.'COEFFICIENT' ; NFLU = DON.FLU.'EXPOSANT' ;
  456.       CEPSS = DON.EPS.'COEFFSEC' ; NEPSS = DON.EPS.'EXPSEC' ;
  457.     'FINSI' ;
  458.   'FINSI' ;
  459.  
  460. **** COURBES POUR LA LOI DE NEUBER
  461.  
  462. * CYCLIQUE
  463.   PROGNS = PROGCS ;
  464.   PROGNSE = PROGCS * PROGCE ;
  465. * TRACTION
  466.   PROGNST = PROGTS ;
  467.   PROGNSET = PROGTS * PROGTE ;
  468.  
  469. ****** PARAMETRES
  470.  
  471.   TM = DON.TPSM ;
  472.   'SI' ('EXISTE' DON TPSTOT ) ;
  473.     TTOTAL = DON.TPSTOT ;
  474.   'SINON' ;
  475.     TTOTAL = 0. ;
  476.   'FINSI' ;
  477.  
  478.   'SI' ('EXISTE' DON NBCYCLE ) ;
  479.     NBCYT = DON.NBCYCLE ;
  480.   'SINON' ;
  481.     NBCYT = 0 ;
  482.   'FINSI' ;
  483.  
  484.   'SI' (('EXISTE' DON AMAX) 'EGA' VRAI) ;
  485.     'SI' ((DON.AMAX) 'EGA' VRAI ) ;
  486.       'MESS' 'LA FISSURE EST TRAVERSANTE' ;
  487.       'QUIT' BOUREA ;
  488.     'FINSI' ; 
  489.   'FINSI' ;
  490.  
  491.  
  492.  
  493. ****** INITIALISATION DES CHAMPS CONSTANTS
  494.  
  495.   CHR02 = 'MANU' 'CHML' OBJF 'SCAL' R002 'TYPE' 'SCALAIRE'
  496.           'STRESSES' ;
  497.   CHAM1 = 'MANU' 'CHML' OBJF 'SCAL' 1. 'TYPE' 'SCALAIRE'
  498.           'STRESSES' ;
  499.   CHAM0 = 'MANU' 'CHML' OBJF 'SCAL' 0. 'TYPE' 'SCALAIRE'
  500.           'STRESSES' ;
  501.   NB_PGAUS = CHAM1 'MASQ' 'SUPERIEUR' 0.5 'SOMME' ;
  502.  
  503.  
  504. ****************************************************************
  505. *    EXTRACTION DES COMPOSANTES DES CHAMPS                     *
  506. ****************************************************************
  507.  
  508. *** FATIGUE ***
  509.  
  510.   'SI' ('NON' ('EXISTE' DON REPRISE)) ;
  511.  
  512. *---- DELTA-CHARGEMENT NOMINAL COMPLET
  513.  
  514.     DELT_NNO = 'EXCO' SIGNO 'NZZ' 'SCAL' ;
  515.     DELT_MNO = 'EXCO' SIGNO 'MXX' 'SCAL' ;
  516.  
  517.     'OUBLIER' SIGNO ;
  518.  
  519.  
  520. *---- DELTA-CHARGEMENT NOMINAL PRIMAIRE
  521.  
  522.     DELT_N1N = 'EXCO' SIGNO1 'NZZ' 'SCAL' ;
  523.     DELT_M1N = 'EXCO' SIGNO1 'MXX' 'SCAL' ;
  524.  
  525.     'OUBLIER' SIGNO1 ;
  526.  
  527.   'FINSI' ;
  528.  
  529.  
  530. *---- DELTA-CHARGEMENT PRIMAIRE
  531.  
  532.   DELT_N1 = 'EXCO' DSIG1 'NZZ' 'SCAL' ;
  533.   DELT_M1 = 'EXCO' DSIG1 'MXX' 'SCAL' ;
  534.   CH_NZZ = 'EXCO' DSIG1 'NZZ' ;CH_MXX = 'EXCO' DSIG1 'MXX' ;
  535.   CH_KI = 'EXCO' DSIG1 'KI' ;
  536.   DSIG1FIS = CH_NZZ 'ET' CH_MXX 'ET' CH_KI ;
  537.  
  538. *---- DELTA-CHARGEMENT 
  539.  
  540.   DELT_N = 'EXCO' DSIG 'NZZ' 'SCAL' ;
  541.   DELT_M = 'EXCO' DSIG 'MXX' 'SCAL' ;
  542.   CH_NZZ = 'EXCO' DSIG 'NZZ' ;CH_MXX = 'EXCO' DSIG 'MXX' ;
  543.   CHKI = 'EXCO' DSIG 'KI' 'SCAL' ;
  544.   CH_KI = 'EXCO' DSIG 'KI' ;
  545.   DSIGFIS = CH_NZZ 'ET' CH_MXX 'ET' CH_KI ;
  546.  
  547.   'OUBLIER' DSIG ;
  548.  
  549. *** FLUAGE ***
  550.  
  551.   'SI' ('NON' ('EXISTE' DON REPRISE)) ;
  552.  
  553. *---- CHARGEMENT NOMINAL B COMPLET
  554.  
  555.     NBNO = 'EXCO' SIGBNO 'NZZ' 'SCAL' ;
  556.     MBNO = 'EXCO' SIGBNO 'MXX' 'SCAL' ;
  557.  
  558.     'OUBLIER' SIGBNO ;
  559.  
  560. *---- CHARGEMENT NOMINAL B PRIMAIRE
  561.  
  562.     NBNO1 = 'EXCO' SIGBNO1 'NZZ' 'SCAL' ;
  563.     MBNO1 = 'EXCO' SIGBNO1 'MXX' 'SCAL' ;
  564.  
  565.     'OUBLIER' SIGBNO1 ;
  566.  
  567.   'FINSI' ;
  568.  
  569.  
  570. *---- CHARGEMENT A COMPLET
  571.  
  572.   NA = 'EXCO' SIGA 'NZZ' 'SCAL' ;
  573.   MA = 'EXCO' SIGA 'MXX' 'SCAL' ;
  574.   KIA = 'EXCO' SIGA 'KI' 'SCAL' ;
  575.   CH_NZZ = 'EXCO' SIGA 'NZZ' ;CH_MXX = 'EXCO' SIGA 'MXX' ;
  576.   CH_KI = 'EXCO' SIGA 'KI' ;
  577.   SIGAFIS = CH_NZZ 'ET' CH_MXX 'ET' CH_KI  ;
  578.  
  579.   'OUBLIER' SIGA ;
  580.  
  581. *---- CHARGEMENT B COMPLET
  582.  
  583.   NB = 'EXCO' SIGB 'NZZ' 'SCAL' ;
  584.   MB = 'EXCO' SIGB 'MXX' 'SCAL' ;
  585.   KIB = 'EXCO' SIGB 'KI' 'SCAL' ;
  586.   CH_NZZ = 'EXCO' SIGB 'NZZ' ;CH_MXX = 'EXCO' SIGB 'MXX' ;
  587.   CH_KI = 'EXCO' SIGB 'KI' ;
  588.   SIGBFIS = CH_NZZ 'ET' CH_MXX 'ET' CH_KI ;
  589.  
  590.  'OUBLIER' SIGB ;
  591.  
  592. *---- CHARGEMENT A PRIMAIRE
  593.  
  594.   NA1 = 'EXCO' SIGA1 'NZZ' 'SCAL' ;
  595.   MA1 = 'EXCO' SIGA1 'MXX' 'SCAL' ;
  596.   KIA1 = 'EXCO' SIGA1 'KI' 'SCAL' ;
  597.   CH_NZZ = 'EXCO' SIGA1 'NZZ' ;CH_MXX = 'EXCO' SIGA1 'MXX' ;
  598.   CH_KI = 'EXCO' SIGA1 'KI' ;
  599.   SIGA1FIS = CH_NZZ 'ET' CH_MXX 'ET' CH_KI ;
  600.  
  601.   'OUBLIER' SIGA1 ;
  602.  
  603. *---- CHARGEMENT B PRIMAIRE
  604.  
  605.   NB1 = 'EXCO' SIGB1 'NZZ' 'SCAL' ;
  606.   MB1 = 'EXCO' SIGB1 'MXX' 'SCAL' ;
  607.   KIB1 = 'EXCO' SIGB1 'KI' 'SCAL' ;
  608.   CH_NZZ = 'EXCO' SIGB1 'NZZ' ;CH_MXX = 'EXCO' SIGB1 'MXX' ;
  609.   CH_KI = 'EXCO' SIGB1 'KI' ;
  610.   SIGB1FIS = NB1 'ET' MB1 'ET' CH_KI ;
  611.  
  612.   'OUBLIER' SIGB1 ;
  613.  
  614.   'MENAGE' ;
  615.  
  616. ****************************************************************
  617. *          VALEURS DE FATIGUE (CALCULS PRELIMINAIRES)          *
  618. ****************************************************************
  619.  
  620. **** CALCUL DU SIGMA NOMINAL REEL
  621.  
  622.   'SI' ('NON' ('EXISTE' DON REPRISE)) ;
  623.  
  624.     DS_NOS = @SIGREF DELT_NNO DELT_MNO W CHAM0 OBJF ;
  625.     DS_NO  = 'ABS' DS_NOS ;
  626.  
  627.     DS_NO1S = @SIGREF DELT_N1N DELT_M1N W CHAM0 OBJF ;
  628.     DS_NO1  = 'ABS' DS_NO1S ;
  629.  
  630.     DE_NO1 = 'IPOL' DS_NO1 PROGCS PROGCE ;
  631.     DE_NOR = (DS_NO - DS_NO1 + (ECYCL * DE_NO1)) * (ECYCL ** -1);
  632.     DS_NOR = 'IPOL' DE_NOR PROGCE PROGCS ;
  633.  
  634.     DON.'DELSIGNO' = DS_NO ;
  635.     DON.'DELSIGNOR' = DS_NOR ;
  636.  
  637.   'SINON' ;
  638.  
  639.      DS_NO  = DON.'DELSIGNO' ;
  640.      DS_NOR = DON.'DELSIGNOR' ;
  641.  
  642.   'FINSI' ;
  643.  
  644. ** -- CALCUL DES CONTRAINTES PRIMAIRES PRES DU DEFAUT 
  645.  
  646.    DS_DEF1S = @SIGREF DELT_N1 DELT_M1 W A OBJF ;
  647.    DS_DEF1  = 'ABS' DS_DEF1S ;
  648.  
  649.    DE_DEF1 = 'IPOL' DS_DEF1 PROGCS PROGCE ;
  650.  
  651. ** -- CALCUL DES CONTRAINTES RELLES PRES DU DEFAUT 
  652.  
  653.    DS_DEFS = @SIGREF DELT_N DELT_M W A OBJF;
  654.    DS_DEF  = 'ABS' DS_DEFS ; 
  655.  
  656.    DE_DEFP = (DS_DEF - DS_DEF1 + (ECYCL * DE_DEF1)) * (ECYCL ** -1) ;
  657.  
  658.    SE = DS_DEF * DE_DEFP ;
  659.    DS_DEFR = 'IPOL' SE PROGNSE PROGNS ;
  660.    DE_DEFR = 'IPOL' DS_DEFR PROGCS PROGCE ;
  661.  
  662.  
  663. ****************************************************************
  664. *          VALEURS DE FLUAGE (CALCULS PRELIMINAIRES)           *
  665. **************************************************************** 
  666.  
  667.   'SI' ('NON' ('EXISTE' DON REPRISE)) ;
  668.  
  669. ** -- CALCUL DE LA CONTRAINTE NOMINALE PRIMAIRE EN B
  670.  
  671.     S_NOB1S = @SIGREF NBNO1 MBNO1 W CHAM0 OBJF ;
  672.     S_NOB1  = 'ABS' S_NOB1S ;
  673.  
  674.     E_NOB1 = 'IPOL' S_NOB1 PROGTS PROGTE ;
  675.  
  676. ** -- CALCUL DE LA CONTRAINTE NOMINALE REELLE
  677.  
  678.     S_NOBS = @SIGREF NBNO MBNO W CHAM0 OBJF ;
  679.     S_NOB  = 'ABS' S_NOBS ;
  680.  
  681.     E_NORB = (S_NOB - S_NOB1 + (E * E_NOB1)) * (E ** -1) ;
  682.     S_NORB = 'IPOL' E_NORB PROGTE PROGTS ;
  683.  
  684.     DON.'SIGNOB' = S_NOB ;
  685.     DON.'SIGNORB' = S_NORB ;
  686.  
  687.   'SINON' ;
  688.  
  689.     S_NOB  = DON.'SIGNOB' ;
  690.     S_NORB = DON.'SIGNORB' ;
  691.  
  692.   'FINSI' ;
  693.  
  694. ** -- CALCUL DES CONTRAINTES PRIMAIRES EN A
  695.  
  696.    SIGA1P = KIA1 'MASQ' 'SUPERIEUR' 0. ;
  697.    SIGA1N = (KIA1 'MASQ' 'INFERIEUR' 0.) * (-1.) ;
  698.    SIGA1S = SIGA1P + SIGA1N ;
  699.  
  700.    S_DEFA1S = @SIGREF NA1 MA1 W A OBJF;
  701.    TESTA11 = S_DEFA1S * SIGA1S ;
  702.    TESTA12 = TESTA11 'MASQ' 'INFERIEUR' 0. 'SOMME' ;
  703.    'SI' ((TESTA12 / NB_PGAUS) > 0.05) ;
  704.     'MESS' 'LES SIGNES ISSUS DES S_REFA1 ET DES KA1 NE SONT PAS
  705.  COHERENTS SUR ' TESTA12 'POINTS DE GAUSS' ;
  706.    'FINSI' ;
  707.  
  708.    S_DEFA1  = 'ABS' S_DEFA1S ;
  709.    E_DEFA1 = 'IPOL' S_DEFA1 PROGTS PROGTE ;
  710.  
  711. ** -- CALCUL DES CONTRAINTES DEFAUT REEL EN A
  712.  
  713.    SIGAP = KIA 'MASQ' 'SUPERIEUR' 0. ;
  714.    SIGAN = (KIA 'MASQ' 'INFERIEUR' 0.) * (-1.) ;
  715.    SIGAS = SIGAP + SIGAN ;
  716.  
  717.    S_DEFAS = @SIGREF NA MA W A OBJF ;
  718.    TESTA1 = S_DEFAS * SIGAS ;
  719.    TESTA2 = TESTA1 'MASQ' 'INFERIEUR' 0. 'SOMME' ;
  720.    'SI' ((TESTA2 / NB_PGAUS) > 0.05) ;
  721.     'MESS' 'LES SIGNES ISSUS DES S_REFA ET DES KA NE SONT PAS
  722.  COHERENTS SUR ' TESTA2 'POINTS DE GAUSS' ;
  723.    'FINSI' ;
  724.  
  725.    S_DEFA  = 'ABS' S_DEFAS ;
  726.    E_DEFA =  (S_DEFA - S_DEFA1 + (E * E_DEFA1)) * (E ** -1) ;
  727.    SE = S_DEFA * E_DEFA ;
  728.    S_DEFRA = 'IPOL' SE PROGNSET PROGNST ;
  729.  
  730. ** -- CALCUL DES CONTRAINTES PRIMAIRES EN B
  731.  
  732.    SIGB1P = KIB1 'MASQ' 'SUPERIEUR' 0. ;
  733.    SIGB1N = (KIB1 'MASQ' 'INFERIEUR' 0.) * (-1.) ;
  734.    SIGB1S = SIGB1P + SIGB1N ;
  735.  
  736.    S_DEFB1S = @SIGREF NB1 MB1 W A OBJF ;
  737.    TESTB11 = S_DEFB1S * SIGB1S ;
  738.    TESTB12 = TESTB11 'MASQ' 'INFERIEUR' 0. 'SOMME' ;
  739.    'SI' ((TESTB12 / NB_PGAUS) > 0.05) ;
  740.     'MESS' 'LES SIGNES ISSUS DES S_REFB1 ET DES KB1 NE SONT PAS
  741.  COHERENTS SUR ' TESTB12 'POINTS DE GAUSS' ;
  742.    'FINSI' ;
  743.  
  744.    S_DEFB1  = 'ABS' S_DEFB1S ;
  745.    E_DEFB1 = 'IPOL' S_DEFB1 PROGTS PROGTE ;
  746.  
  747. ** -- CALCUL DES CONTRAINTES DEFAUT REEL EN B
  748.  
  749.    SIGBP = KIB 'MASQ' 'SUPERIEUR' 0. ;
  750.    SIGBN = (KIB 'MASQ' 'INFERIEUR' 0.) * (-1.) ;
  751.    SIGBS = SIGBP + SIGBN ;
  752.  
  753.    S_DEFBS = @SIGREF NB MB W A OBJF ;
  754.    TESTB1 = S_DEFBS * SIGBS ;
  755.    TESTB2 = TESTB1 'MASQ' 'INFERIEUR' 0. 'SOMME' ;
  756.    'SI' ((TESTB2 / NB_PGAUS) > 0.05) ;
  757.     'MESS' 'LES SIGNES ISSUS DES S_REFB ET DES KB NE SONT PAS
  758.  COHERENTS SUR ' TESTB2 'POINTS DE GAUSS' ;
  759.    'FINSI' ;
  760.  
  761.    S_DEFB  = 'ABS' S_DEFBS ;
  762.    E_DEFB =  (S_DEFB - S_DEFB1 + (E * E_DEFB1)) * (E ** -1) ;
  763.    SE = S_DEFB * E_DEFB ;
  764.    S_DEFRB = 'IPOL' SE PROGNSET PROGNST ;
  765.  
  766.  
  767.  
  768. ****************************************************************
  769. *         CONTRIBUTION DE LA FATIGUE (SUR 1 CYCLE)             *
  770. ****************************************************************
  771.  
  772.  
  773. *- CORRECTION DE PLASTICITE :
  774.  
  775.   DS2DEFR = DS_DEFR ** 2 ;CHR022 = CHR02 ** 2 ;
  776.   CHR022 = CHR022 * 4. ;
  777.   PSI_DNO = (DS2DEFR + CHR022) ;PSI_DNO = PSI_DNO ** -1 ;
  778.   PSI = 0.5 * DS2DEFR * PSI_DNO ;
  779.  
  780. *- METHODE 'JS' :
  781.   DEDS = DS_DEFR ** -1 ;DEDS = DE_DEFR * DEDS ;
  782.   DEDSE = ECYCL * DEDS ;PSIDEDS = PSI + DEDSE ;
  783.   NORNO = DS_NO ** -1 ;NORNO = DS_NOR * NORNO ;
  784.   NORNO2 = NORNO ** 2 ;
  785.   K_DJ = PSIDEDS * NORNO2  ;
  786.  
  787. *- ON PREND LA VALEUR ABSOLUE DES 'DELTA-FACTEUR D INTENSITE DE CONTRAINTE'
  788.  
  789. *  DKI = CHKI 'MASQ' 'SUPERIEUR' 0. ;
  790. *  DKI = CHKI * DKI ;
  791.   DKI = 'ABS' CHKI ;
  792.  
  793.   DKI2 = DKI ** 2 ;
  794.   DJ_EL = DKI2 * (ECYCL ** -1);
  795.  
  796.   DJ = DJ_EL * K_DJ ;
  797.  
  798. **** CALCUL DE R ET CHOIX DE Q
  799. **** ---- CALCUL DE R
  800.  
  801.   S_DEFRAS = S_DEFRA * SIGAS ;S_DEFRBS = S_DEFRB * SIGBS ;
  802.   DSOMM = S_DEFRAS + S_DEFRBS ;DSOMM = DSOMM / 2. ;
  803.   DS2DEFR = 0.5 * DS_DEFR ;
  804.   RESU1 = DSOMM - DS2DEFR ;
  805.   RESUS1 = RESU1 'MASQ' 'INFERIEUR' 0. ;
  806.   RESU2 = DSOMM + DS2DEFR ;
  807.   RESUS2 = RESU2 'MASQ' 'INFERIEUR' 0. ;
  808.  
  809. *- LES ALTERNANCES NEGATIF/NEGATIF SONT MISES A ZERO
  810.  
  811.   RESU0 = RESUS1 * RESUS2 ;RESUF = RESU0 'MASQ' 'INFERIEUR' 0.99 ;
  812.  
  813.   R = RESU1 * (RESU2 ** (-1)) ;
  814.  
  815.  
  816.  
  817. **** ---- CHOIX DE Q
  818.  
  819.   CHRNEG = R 'MASQ' 'INFERIEUR' 0. ;
  820.   CHRPOS = R 'MASQ' 'SUPERIEUR' 0. ;
  821.   RS2 = 0.5 * R ;RS2 = CHAM1 - RS2 ;
  822.   RS1 = CHAM1 - R ;RS1 = RS1 ** -1 ;
  823.   QPOS = (RS2 ** -1) * CHRPOS ;
  824.   QNEG = (RS2 * RS1) * CHRNEG ;
  825.  
  826.   Q = QPOS + QNEG ;Q = Q * RESUF ;
  827.  
  828. **** CALCUL DE DELTA-K EFFICACE  
  829.  
  830.   MULTI = 1.E-3 ** 0.5 ;DJCY = ECYCL * DJ ;DJCY = DJCY ** 0.5 ;
  831.   DK_EFF = Q * DJCY * MULTI ;
  832.  
  833. **** CALCUL DE LA PROPAGATION EN FATIGUE PRIMAIRE SUR 1 CYCLE
  834.  
  835.   DA_FAT1 = DK_EFF ** NFAT ;
  836.   DA_FAT = CFAT * DA_FAT1 * 1000. ;
  837.  
  838.  
  839.  
  840. ****************************************************************
  841. *              CALCUL DE PROPAGATION SUR N CYCLE               *
  842. **************************************************************** 
  843.  
  844. * INITIALISATION D'ACTUALISATION
  845.  
  846.   ATOTAL = ADEB ;
  847.   ATOTAL1 = ADEB1 ;ATOTAL2 = ADEB2 ; 
  848.   NBCY = 0. ;
  849.   EF = CHAM0 ;EFP = CHAM0 ;EFS = CHAM0 ;EF1 = CHAM0 ;
  850.   KAC = 1 ; ACC = 1. ; KCRIT = 0 ;
  851.  
  852.   'SI' ('EXISTE' DON LR_NBCYT) ;
  853.     L_NBCYT = DON.LR_NBCYT ;
  854.   'SINON' ;
  855.     L_NBCYT = 'PROG' ;
  856.   'FINSI' ;
  857.  
  858.   'SI' ('EXISTE' DON LR_DAFAT) ;
  859.     L_DA_FAT = DON.LR_DAFAT ;
  860.   'SINON' ;
  861.     L_DA_FAT = 'PROG' ;
  862.   'FINSI' ;
  863.  
  864.   'SI' ('EXISTE' DON LR_ATOT) ;
  865.     L_ATOT = DON.LR_ATOT ;
  866.   'SINON' ;
  867.     L_ATOT = 'PROG' ;
  868.   'FINSI' ;
  869.  
  870.   'SI' ('EXISTE' DON LR_SEFB) ;
  871.     L_SEFB = DON.LR_SEFB ;
  872.   'SINON' ;
  873.     L_SEFB = 'PROG' ;
  874.   'FINSI' ;
  875.  
  876.   'SI' ('EGA' CALFAT 0) ;
  877.  
  878.     'SI' (('EXISTE' DON LR_ATOT1) 'ET' ('EXISTE' DON LR_ATOT2)) ;
  879.       L_ATOT1 = DON.LR_ATOT1 ;
  880.       L_ATOT2 = DON.LR_ATOT2 ;
  881.     'SINON' ;
  882.       L_ATOT1 = 'PROG' ;
  883.       L_ATOT2 = 'PROG' ;
  884.     'FINSI' ;
  885.  
  886.     'SI' ('EXISTE' DON LR_DAFLU) ;
  887.       L_DA_FLU = DON.LR_DAFLU ;
  888.     'SINON' ;
  889.       L_DA_FLU = 'PROG' ;
  890.     'FINSI' ;
  891.  
  892.     'SI' ('EXISTE' DON LR_EPFLU) ;
  893.       L_EP_FLU = DON.LR_EPFLU ;
  894.     'SINON' ;
  895.       L_EP_FLU = 'PROG' ;
  896.     'FINSI' ;
  897.  
  898. *--------------------------------------------------------------*
  899. *                 CALCUL EN FATIGUE/FLUAGE                     *
  900. *--------------------------------------------------------------*
  901.  
  902. *- MISE A ZERO DES FACTEUR D INTENSITE DE CONTRAINTE NULS DU AU CHARGT B
  903.  
  904.     RKI = KIB 'MASQ' 'SUPERIEUR' 0. ; 
  905.     KIB = KIB * RKI ;
  906.  
  907.     COEF = (KIB ** 2) * ((S_NORB * (S_NOB ** -1)) ** 2) ;
  908.  
  909. *- CALCUL DE SIG-EFFICACE-B ET MISE A 'EPSILON' DES VALEURS NEGATIVES
  910.  
  911.     S_EFFBC =  ((S_DEFRAS + S_DEFRBS) / 2.) + (0.5 * DS_DEFR) ;
  912.     SEFMAX = 'MAXI' S_EFFBC ;SEFLIM = 'ABS' (SEFMAX * 1.E-6) ;
  913.     CHEFFBP = 'MASQ' S_EFFBC 'SUPERIEUR' SEFLIM ;
  914.     CHEFFBN = 'MASQ' S_EFFBC 'INFERIEUR' SEFLIM ;
  915.     S_EFFBCP = S_EFFBC * CHEFFBP ;
  916.     S_EFFBCN = SEFLIM * CHEFFBN ;
  917.     S_EFFBC = S_EFFBCP + S_EFFBCN  ;
  918.  
  919.  
  920.     'SI' ('EXISTE' DON PROPA) ;
  921.       NBFISS = DON.PROPA ;
  922.     'SINON' ;
  923.       NBFISS = 500 ;
  924.     'FINSI' ;
  925.  
  926.     'SI' ('EXISTE' DON EPSF1) ;
  927.       EF1 = DON.EPSF1 ;
  928.     'SINON' ;
  929.       EF1 = CHAM0 ;
  930.     'FINSI' ;
  931.  
  932.  
  933.     WARN = AFIN 'MASQ' 'SUPERIEUR' 'SOMME' WIMP ;
  934.     'SI' (WARN > 0.5) ;
  935.       ALIM = WIMP ;TRAV = VRAI ;
  936.     'SINON' ;
  937.       ALIM = AFIN ;TRAV = FAUX ;
  938.     'FINSI' ;
  939.  
  940.     S_EFFB = S_EFFBC ; 
  941.     CPTG = 0 ;
  942.     'REPETER' CYCLAGE NBFISS ; 
  943.       CPTG = CPTG + 1 ;KCRIT = KCRIT + 1 ; TMV = ACC * TM ;
  944.  
  945.       'MENAGE' ;
  946.  
  947.  
  948. **** CALCUL DE LA PROPAGATION EN FLUAGE
  949.  
  950.       'SI' (A316 'EGA' 1) ;
  951.  
  952.         NEPSP = YFP ;
  953.         CEPSS = CFS ;
  954.         NEPSS = NFS ;
  955.         TFP = C3 * (S_EFFBC ** N3) ;
  956.         CHP = TFP 'MASQ' 'SUPERIEUR' TTOTAL ;
  957.         CHS = TFP 'MASQ' 'INFERIEUR' TTOTAL ;
  958.         NBP = TFP 'MASQ' 'SUPERIEUR' 'SOMME' TTOTAL ;
  959.         NBS = TFP 'MASQ' 'INFERIEUR' 'SOMME' TTOTAL ;
  960.  
  961. ****** CALCUL AVEC RELAXATION - A316 (FLUAGE PRIMAIRE ET SECONDAIRE)
  962.  
  963.         'SI' ('EGA' RELAX 1) ;
  964.           INDIC = 1 ;
  965.           NBR = 100 ;SUBD = TMV / NBR ;TIC = NBCY * TM ;
  966.           'SI' (NBCY < 0.1) ;
  967.             S_EFFBDP = S_EFFB ;
  968.             S_EFFBDS = S_EFFB ;
  969.           'FINSI' ;
  970.           'SI' (NBP 'NEG' 0) ;
  971.             CPT = 0 ;AC = 1. ;TC = 0. ;KA = 1 ;STA = 0 ;QT = 0 ;
  972.             S_EFFBD = S_EFFBDP ;
  973.             'REPETER' RELAX1 NBR ;
  974.               'SI' (QT 'EGA' 0) ;
  975.                 DT = SUBD * AC ;
  976.               'FINSI' ;
  977.               CPT = CPT + 1 ;TC = TC + DT ;STA = STA + 1 ;
  978.               EFEI = S_EFFBD ** N1 ;
  979.               TTC = TIC + TC ;C2P = C2 - 1. ;EFTT1 = TTC ** C2P ;
  980.               EFE = C1 * C2 * EFTT1 * EFEI * DT ;
  981.               EFP = EFP + EFE ;
  982.               EFPCUM = EFP ;
  983.               CEPSP1 = EFPCUM ** XFP ;
  984.               CEPSP = KFP * CEPSP1 * 0.01 ;YFPP = 1. - YFP ;
  985.               YFPPI = 1. / YFPP ;
  986.               S_EFFBF1 = S_EFFBD ** YFPP ;
  987.               S_EFFBF2 = YFPP * EY * CEPSP * DT ;
  988.               S_EFFBF3 = S_EFFBF1 - S_EFFBF2 ;
  989.               S_EFFBF = S_EFFBF3 ** YFPPI ;
  990.               S_EFFBI1 = S_EFFBD + S_EFFBF ;
  991.               S_EFFBI = S_EFFBI1 / 2. ;S_EFFBI2 = S_EFFBI ** NEPSP ;
  992.               EPS_PT = CEPSP * S_EFFBI2 ;S_EFFBI1 = S_EFFBI ** -1 ;
  993.               CSTAR1 = EPS_PT * S_EFFBI1 ;
  994.               CSTAR = COEF * CSTAR1 * 1.E-3 ;
  995.               DA_FLUEI = CSTAR ** NFLU ;
  996.               DA_FLUE = CFLU * DA_FLUEI * DT * 1000. ;
  997.               'SI' (CPT 'EGA' 1) ;
  998.                 DA_FLU = DA_FLUE ;
  999.                 DEFP = EFE ;
  1000.               'SINON';
  1001.                 DA_FLU = DA_FLU + DA_FLUE ;
  1002.                 DEFP = DEFP + EFE ;
  1003.               'FINSI';
  1004.               'SI' (QT 'EGA' 1) ;
  1005.                 'QUIT' RELAX1 ;
  1006.               'FINSI' ;
  1007.               RESID = TMV - TC ;
  1008.               'SI' (RESID < DT) ;
  1009.                 DT = RESID ;QT = 1 ;
  1010.               'FINSI' ;
  1011.               'SI' (STA > 1.1) ;
  1012.                 MDA_FLUE = 'MAXI' DA_FLUE ;MDAFLUE2 = 'MAXI' DA_FLUE2;
  1013.                 MDAFLUE2 = MDAFLUE2 ** -1 ;
  1014.                 CONTIMI = MDA_FLUE * MDAFLUE2 ;
  1015.                 M_EFE = 'MAXI' EFE ;M_EFE2 = 'MAXI' EFE2 ;
  1016.                 M_EFE2 = M_EFE2 ** -1 ;
  1017.                 EFETIMI = M_EFE * M_EFE2 ;
  1018.                 MS_EFFBI = 'MAXI' S_EFFBI ;MSEFFBI2 = 'MAXI' S_EFFB2 ;
  1019.                 MSEFFBI2 = MSEFFBI2 ** -1 ;
  1020.                 SIGTIMI = MS_EFFBI * MSEFFBI2 ;
  1021.                 'SI' ((CONTIMI > 0.99) 'ET' (SIGTIMI > 0.99)
  1022.                     'ET' (EFETIMI > 0.99))  ;
  1023.                   KA = KA + 1 ; AC = KA ** 2; STA = 0 ;
  1024.                 'FINSI' ;
  1025.               'FINSI' ;
  1026.               DA_FLUE2 = DA_FLUE ;
  1027.               S_EFFB2 = S_EFFBI ;
  1028.               EFE2 = EFE ;
  1029.               S_EFFBD = S_EFFBF ;
  1030.             'FIN' RELAX1 ;
  1031.             'SI' (CPT 'EGA' NBR) ;
  1032.               MESS 'RELAXATION SEVERE' ;
  1033.             'FINSI' ;
  1034.             S_EFFBDP = S_EFFBF ;
  1035.           'SINON' ;
  1036.             EPS_PT = CHAM1 ;
  1037.             DA_FLU = CHAM1 ;
  1038.             EFP    = CHAM1 ;
  1039.             DEFP   = CHAM1 ;
  1040.             S_EFFBDP = S_EFFB ;
  1041.           'FINSI' ;
  1042.  
  1043.           S_EFBDP = S_EFFBDP * CHP ;
  1044.           EPS_PTP = EPS_PT * CHP ;
  1045.           DA_FLUP = DA_FLU * CHP ;
  1046.           EF_P    = EFP    * CHP ;
  1047.           DEF_P   = DEFP   * CHP ;
  1048.  
  1049.           'SI' (NBS 'NEG' 0) ;  
  1050.             CPT = 0 ;AC = 1. ;TC = 0. ;KA = 1 ;STA = 0 ;QT = 0 ;
  1051.             S_EFFBD = S_EFFB ;NEPSSP = 1. - NEPSS ;
  1052.             NEPSSPP = NEPSSP ** -1 ;
  1053.             'REPETER' RELAX2 NBR ;
  1054.               'SI' (QT 'EGA' 0) ;
  1055.                 DT = SUBD * AC ;
  1056.               'FINSI' ;
  1057.               CPT = CPT + 1 ;TC = TC + DT ;
  1058.               STA = STA + 1 ;VALCH2 = EY * CEPSS * NEPSSP * DT ;
  1059.               CH2 = 'MANU' 'CHML' OBJF 'SCAL' VALCH2
  1060.                   'TYPE' 'SCALAIRE' 'STRESSES' ;
  1061.               EFEI = S_EFFBD ** NEPSS ;
  1062.               EFE = 100. * CEPSS * EFEI * DT ;
  1063.               EFS = EFS + EFE ;S_EFFBF1 = S_EFFBD ** NEPSSP ;
  1064.               S_EFFBF2 = S_EFFBF1 - CH2 ;
  1065.               S_EFFBF =  S_EFFBF2 ** NEPSSPP;
  1066.               S_EFFBI1 = S_EFFBD + S_EFFBF ;
  1067.               S_EFFBI = S_EFFBI1 / 2. ;S_EFFBI2 = S_EFFBI ** NEPSS ;
  1068.               EPS_PT = CEPSS * S_EFFBI2 ;S_EFFBI1 = S_EFFBI ** -1 ;
  1069.               CSTAR1 = EPS_PT * S_EFFBI1 ;
  1070.               CSTAR = COEF * CSTAR1 * 1.E-3 ;DA_FLUEI = CSTAR ** NFLU ;
  1071.               DA_FLUE = CFLU * DA_FLUEI * DT * 1000. ;
  1072.               'SI' (CPT 'EGA' 1) ;
  1073.                 DA_FLU = DA_FLUE ;
  1074.                 DEFS   = EFE ;
  1075.               'SINON';
  1076.                 DA_FLU = DA_FLU + DA_FLUE ;
  1077.                 DEFS   = DEFS + EFE ;
  1078.               'FINSI';
  1079.               'SI' (QT 'EGA' 1) ;
  1080.                 'QUIT' RELAX2 ;
  1081.               'FINSI' ;
  1082.               RESID = TMV - TC ;
  1083.               'SI' (RESID < DT) ;
  1084.                 DT = RESID ;QT = 1 ;
  1085.               'FINSI' ;
  1086.               'SI' (STA > 1.1) ;
  1087.                 MDA_FLUE = 'MAXI' DA_FLUE ;MDAFLUE2 = 'MAXI' DA_FLUE2;
  1088.                 MDAFLUE2 = MDAFLUE2 ** -1 ;
  1089.                 CONTIMI = MDA_FLUE * MDAFLUE2 ;
  1090.                 M_EFE = 'MAXI' EFE ;M_EFE2 = 'MAXI' EFE2 ;
  1091.                 M_EFE2 = M_EFE2 ** -1 ;
  1092.                 EFETIMI = M_EFE * M_EFE2 ;
  1093.                 M_EFFBI = 'MAXI' S_EFFBI ;M_EFFB2 = 'MAXI' S_EFFB2 ;
  1094.                 M_EFFB2 = M_EFFB2 ** -1 ;
  1095.                 SIGTIMI = M_EFFBI * M_EFFB2 ;
  1096.                 'SI' ((CONTIMI > 0.99) 'ET' (SIGTIMI > 0.99)
  1097.                     'ET' (EFETIMI > 0.99))  ;
  1098.                   KA = KA + 1 ; AC = KA ** 2; STA = 0 ;
  1099.                 'FINSI' ;
  1100.               'FINSI' ;
  1101.               DA_FLUE2 = DA_FLUE ;
  1102.               S_EFFB2 = S_EFFBI ;
  1103.               EFE2 = EFE ;
  1104.               S_EFFBD = S_EFFBF ;
  1105.             'FIN' RELAX2 ;
  1106.             'SI' (CPT 'EGA' NBR) ;
  1107.               MESS 'RELAXATION SEVERE' ;
  1108.             'FINSI' ;
  1109.             S_EFFBDS = S_EFFBF ;
  1110.           'SINON' ;
  1111.             EPS_PT = CHAM1 ;
  1112.             DA_FLU = CHAM1 ;
  1113.             EFS    = CHAM1 ;
  1114.             DEFS   = CHAM1 ;
  1115.             S_EFFBDS = S_EFFB ;
  1116.           'FINSI';
  1117.  
  1118.           S_EFBDS = S_EFFBDS * CHS ;
  1119.           EPS_PTS = EPS_PT * CHS ;
  1120.           DA_FLUS = DA_FLU * CHS ;
  1121.           EF_S    = EFS    * CHS ;
  1122.           DEF_S   = DEFS   * CHS ;
  1123.  
  1124.           S_EFFB = S_EFBDP + S_EFBDS ;
  1125.           EPS_PT = EPS_PTP + EPS_PTS ;
  1126.           DA_FLU = DA_FLUP + DA_FLUS ;
  1127.           EF     = EF_P + EF_S ;
  1128.           DEF = DEF_P + DEF_S ;
  1129.  
  1130.           EFS = EF ;
  1131.  
  1132.         'SINON' ;
  1133.  
  1134. ******** CALCUL SANS RELAXATION - FLUAGE PRIMAIRE ET SECONDAIRE - A316
  1135.  
  1136.           INDIC = 2 ;
  1137.           NBR = 100 ;S_EFFBM1 = S_EFFB ** -1 ;
  1138.           'SI' (NBP 'NEG' 0) ;
  1139.             TIC = NBCY * TM ;SUBD = TMV / NBR ;
  1140.             EFT2 = S_EFFB ** N1 ;EPS_PT1 = S_EFFB ** NEPSP ;
  1141.             CPT = 0 ;AC = 1. ;TC = 0. ;KA = 1 ;STA = 0 ;QT = 0 ;
  1142.             'REPETER' FLUA1 NBR ;
  1143.               'SI' (QT 'EGA' 0) ;
  1144.                 DT = SUBD * AC ;
  1145.               'FINSI' ;
  1146.               CPT = CPT + 1 ;STA = STA + 1 ;TC = TC + DT ;
  1147.               TTC = TIC + TC ;EFT1 = TTC ** C2 ;
  1148.               EF = C1 * EFT1 * EFT2 ;EFCUM = EF ;
  1149.               CEPSP1 = EFCUM ** XFP ;
  1150.               CEPSP =  KFP * CEPSP1 * 0.01 ;
  1151.               EPS_PT = CEPSP * EPS_PT1 ;
  1152.               CSTAR = COEF * EPS_PT * S_EFFBM1 * 1.E-3 ;
  1153.               DA_FLUEI = CSTAR ** NFLU ;
  1154.               DA_FLUE = CFLU * DA_FLUEI * DT * 1000. ;
  1155.               'SI' (CPT 'EGA' 1) ;   
  1156.                 DA_FLU = DA_FLUE ;
  1157.               'SINON';
  1158.                 DA_FLU = DA_FLU + DA_FLUE ;
  1159.               'FINSI';
  1160.               'SI' (QT 'EGA' 1) ;
  1161.                 'QUIT' FLUA1 ;
  1162.               'FINSI' ;
  1163.               RESID = TMV - TC ;
  1164.               'SI' (RESID < DT) ;
  1165.                 DT = RESID ;QT = 1 ;
  1166.               'FINSI' ;
  1167.               'SI' (STA > 1.1) ;
  1168.                 DA_FLMA1 = 'MAXI' DA_FLUE ;
  1169.                 DA_FLMA2 = 'MAXI' DA_FLUE2 ;DA_FLMA2 = DA_FLMA2 ** -1 ;
  1170.                 CONTIMI = DA_FLMA1 * DA_FLMA2 ;
  1171.                 'SI' (CONTIMI > 0.99) ;
  1172.                   KA = KA + 1 ; AC = KA ** 2; STA = 0 ;
  1173.                 'FINSI' ;
  1174.               'FINSI' ;
  1175.               DA_FLUE2 = DA_FLUE ;
  1176.             'FIN' FLUA1 ;
  1177.           'SINON';
  1178.             EPS_PT = S_EFFB ;
  1179.             DA_FLU = S_EFFB ;
  1180.             EF     = S_EFFB ;
  1181.           'FINSI';
  1182.  
  1183.           S_EFFBP = S_EFFB * CHP ;
  1184.           EPS_PTP = EPS_PT * CHP ;
  1185.           DA_FLUP = DA_FLU * CHP ;
  1186.           EF_P    = EF     * CHP ;
  1187.  
  1188.           'SI' (NBS 'NEG' 0) ;
  1189.             EPS_PT1 = S_EFFB ** NEPSS ;
  1190.             EPS_PT = CEPSS * EPS_PT1 ;
  1191.             CSTAR = COEF * EPS_PT * S_EFFBM1 * 1.E-3 ;
  1192.             DA_FLUI = CSTAR ** NFLU ;
  1193.             DA_FLU = CFLU * DA_FLUI * TMV * 1000. ;
  1194.             EFT1 = S_EFFB ** NEPSS ;EFT1 = 100. * CEPSS * EFT1 * TMV ;
  1195.             EF = EF1 + EFT1 ;
  1196.           'SINON' ;
  1197.             EPS_PT1 =  S_EFFB ;
  1198.             EPS_PT = S_EFFB ;
  1199.             DA_FLU = S_EFFB ;
  1200.             EF     = S_EFFB ;
  1201.           'FINSI';
  1202.  
  1203.           S_EFFBS = S_EFFB * CHS ;
  1204.           EPS_PTS = EPS_PT * CHS ;
  1205.           DA_FLUS = DA_FLU * CHS ;
  1206.           EF_S    = EF     * CHS ;
  1207.  
  1208.  
  1209.           S_EFFB = S_EFFBP + S_EFFBS ;
  1210.           EPS_PT = EPS_PTP + EPS_PTS ;
  1211.           DA_FLU = DA_FLUP + DA_FLUS ;
  1212.           EF = EF_P + EF_S ;
  1213.  
  1214.           EF1 = EF ;
  1215.  
  1216.         'FINSI' ;
  1217.  
  1218.       'SINON' ;
  1219.  
  1220.         'SI' (RELAX 'EGA' 1) ;
  1221.           INDIC = 3 ;
  1222. ******** CALCUL AVEC RELAXATION - FLUAGE SECONDAIRE 
  1223.  
  1224.           'SI' (NBCY < 0.1) ;
  1225.             S_EFFBD = S_EFFB ;
  1226.           'FINSI' ;
  1227.  
  1228.           NBR = 100 ; 
  1229.           SUBD = TMV / NBR ;NEPSSP = 1. - NEPSS ;
  1230.           NEPSSPP = 1. /NEPSSP ; 
  1231.  
  1232.           CPT = 0 ;AC = 1. ;TC = 0. ;KA = 1 ;STA = 0 ;QT = 0 ;
  1233.           'REPETER' RELAX2 NBR ;
  1234.             'SI' (QT 'EGA' 0) ;
  1235.               DT = SUBD * AC ;
  1236.             'FINSI' ;
  1237.             CPT = CPT + 1 ;TC = TC + DT ;STA = STA + 1 ;
  1238.             VALCH2 = EY * CEPSS * NEPSSP * DT ;
  1239.             CH2 = 'MANU' 'CHML' OBJF 'SCAL' VALCH2
  1240.                 'TYPE' 'SCALAIRE' 'STRESSES' ;
  1241.             EFEI = S_EFFBD ** NEPSS ;
  1242.             EFE = 100. * CEPSS * EFEI * DT ;
  1243.             EF = EF + EFE ;S_EFFBF1 = S_EFFBD ** NEPSSP ;
  1244.             S_EFFBF2 = S_EFFBF1 - CH2 ;
  1245.             S_EFFBF = S_EFFBF2 ** NEPSSPP ;
  1246.             S_EFFBI1 = S_EFFBD + S_EFFBF ;
  1247.             S_EFFB = S_EFFBI1 / 2. ;
  1248.             EPS_PT1 = S_EFFB ** NEPSS ;
  1249.             EPS_PT = CEPSS * EPS_PT1 ;S_EFFB1 = S_EFFB ** -1 ;
  1250.             CSTAR1 = EPS_PT * S_EFFB1 ;
  1251.             CSTAR = COEF * CSTAR1 * 1.E-3 ;
  1252.             DA_FLUE1 = CSTAR ** NFLU ;
  1253.             DA_FLUE = CFLU * DA_FLUE1 * DT * 1000. ;
  1254.             'SI' (CPT 'EGA' 1) ;
  1255.               DA_FLU = DA_FLUE ;
  1256.               DEF = EFE ;
  1257.             'SINON';
  1258.               DA_FLU = DA_FLU + DA_FLUE ;
  1259.               DEF = DEF + EFE ;
  1260.             'FINSI';
  1261.             'SI' (QT 'EGA' 1) ;
  1262.                'QUIT' RELAX2 ;
  1263.             'FINSI' ;
  1264.             RESID = TMV - TC ;
  1265.             'SI' (RESID < DT) ;
  1266.               DT = RESID ;QT = 1 ;
  1267.             'FINSI' ;
  1268.             'SI' (STA > 1.1) ;
  1269.               MDA_FLUE = 'MAXI' DA_FLUE ;MDAFLUE2 = 'MAXI' DA_FLUE2;
  1270.               MDAFLUE2 = MDAFLUE2 ** -1 ;
  1271.               CONTIMI = MDA_FLUE * MDAFLUE2 ;
  1272.               M_EFE = 'MAXI' EFE ;M_EFE2 = 'MAXI' EFE2 ;
  1273.               M_EFE2 = M_EFE2 ** -1 ;
  1274.               EFETIMI = M_EFE * M_EFE2 ;
  1275.               M_EFFB = 'MAXI' S_EFFB ;M_EFFB2 = 'MAXI' S_EFFB2 ;
  1276.               M_EFFB2 = M_EFFB2 ** -1 ;
  1277.               SIGTIMI = M_EFFB * M_EFFB2 ;
  1278.               'SI' ((CONTIMI > 0.99) 'ET' (SIGTIMI > 0.99)
  1279.                   'ET' (EFETIMI > 0.99))  ;
  1280.                 KA = KA + 1 ; AC = KA ** 2; STA = 0 ;
  1281.               'FINSI' ;
  1282.             'FINSI' ;
  1283.             DA_FLUE2 = DA_FLUE ;
  1284.             S_EFFB2 = S_EFFB ;
  1285.             EFE2 = EFE ;
  1286.             S_EFFBD = S_EFFBF ;
  1287.           'FIN' RELAX2 ;
  1288.           S_EFFBD = S_EFFBF ;
  1289.  
  1290.         'SINON' ;
  1291.  
  1292. ******** CALCUL SANS RELAXATION - FLUAGE SECONDAIRE
  1293.  
  1294.           INDIC = 4 ;
  1295.           EPS_PT1 = S_EFFB ** NEPSS ;
  1296.           EPS_PT = CEPSS * EPS_PT1 ;
  1297.           S_EFFB1 = S_EFFB ** -1 ;
  1298.           CSTAR1 = EPS_PT * S_EFFB1 ;
  1299.           CSTAR = COEF * CSTAR1 * 1.E-3 ;
  1300.           DA_FLU1 = CSTAR ** NFLU ;
  1301.           DA_FLU = CFLU * DA_FLU1 * TMV * 1000. ;
  1302.           EFT1 = S_EFFB ** NEPSS ;
  1303.           EFT2 = NBCY + ACC ;EFT2 = EFT2 * TM ;
  1304.           EF = 100. * CEPSS * EFT1 * EFT2 ;
  1305.  
  1306.         'FINSI' ;
  1307.       'FINSI' ;
  1308.  
  1309. *** INCREMENTATION DES VALEURS
  1310.       DA_FATC = ACC * DA_FAT ;
  1311.  
  1312.       ATOTAL = ATOTAL + DA_FATC + DA_FLU ;
  1313. *list-reel de la profondeur courante maxi
  1314.       ATOTALMA = 'MAXI' ATOTAL ;PATOTAL = 'PROG' ATOTALMA ;
  1315.       L_ATOT = L_ATOT 'ET' PATOTAL ;
  1316.  
  1317.       ATOTAL1 = ATOTAL1 + DA_FATC ;
  1318. *list-reel de la profondeur courante maxi
  1319.       ATOT1MA = 'MAXI' ATOTAL1 ;PATOT1 = 'PROG' ATOT1MA ;
  1320.       L_ATOT1 = L_ATOT1 'ET' PATOT1 ;
  1321.  
  1322.       ATOTAL2 = ATOTAL2 + DA_FLU ;
  1323. *list-reel de la profondeur courante maxi
  1324.       ATOT2MA = 'MAXI' ATOTAL2 ;PATOT2 = 'PROG' ATOT2MA ;
  1325.       L_ATOT2 = L_ATOT2 'ET' PATOT2 ;
  1326.  
  1327.       TTOTAL = TTOTAL + TMV ;
  1328.       NBCY = NBCY + ACC ;
  1329.       NBCYT = NBCYT + ACC ;
  1330. *list-reel des numero de cycle
  1331.       PNBCYT = 'PROG' NBCYT ;
  1332.       L_NBCYT = L_NBCYT 'ET' PNBCYT ;
  1333.  
  1334.       'SI' (KCRIT > 1.1) ;
  1335.         'SI' (A316 'EGA' 1) ;
  1336.            DA_FLUPM = 'MAXI' DA_FLUP ;
  1337.            DA_FLUOP = DA_FLUO * CHP ;DA_FLUOM = 'MAXI' DA_FLUOP ;
  1338.            CRIT = DA_FLUPM / DA_FLUOM ;
  1339.         'SINON' ;
  1340.           DA_FLUPM = 'MAXI' DA_FLU ;DA_FLUOM = 'MAXI' DA_FLUO ;
  1341.           CRIT = DA_FLUPM / DA_FLUOM ;
  1342.         'FINSI' ;
  1343.       'SINON';
  1344.         CRIT = 0 ;
  1345.       'FINSI';
  1346.       S_EFFBM = 'MAXI' S_EFFB ;
  1347. *list-reel des CONTRAINTES EFFICACES maxi
  1348.       PS_EFFBM = 'PROG' S_EFFBM ;
  1349.       L_SEFB = L_SEFB 'ET' PS_EFFBM ;
  1350.       DA_FATM = 'MAXI' DA_FATC ;
  1351. *list-reel des DA_FATIGUE maxi
  1352.       DA_FATMA = 'MAXI' DA_FAT ;PDA_FAT = 'PROG' DA_FATMA ;
  1353.       L_DA_FAT = L_DA_FAT 'ET' PDA_FAT ;
  1354.  
  1355. *list-reel des DA_FLUAGE maxi
  1356.       DA_FLUO = DA_FLU ;
  1357.       DA_FLUM = 'MAXI' DA_FLU ;DA_FLUMC = DA_FLUM / ACC ;
  1358.       PDA_FLU = 'PROG' DA_FLUMC ;
  1359.       L_DA_FLU = L_DA_FLU 'ET' PDA_FLU ;
  1360.  
  1361. *list-reel des EPSILON-FLUAGE maxi
  1362.       EFCML = EF ;
  1363.       EFCMLM = 'MAXI' EFCML ;PEFCMLM = 'PROG' EFCMLM ;
  1364.       L_EP_FLU = L_EP_FLU 'ET' PEFCMLM ;
  1365.  
  1366.       ATOTAM  = 'MAXI' ATOTAL ;
  1367.       'SI' (NBCY < 1.1) ;
  1368.         'MESS' '                                                ' ;
  1369.         'SI' (INDIC 'EGA' 1) ;
  1370.           'MESS' 'A316    -  RELAXATION (PRIMAIRE ET SECONDAIRE)' ;
  1371.           NBFI3 = NBFI * 3. ;RAT3 = NBP / NBFI3 ;RAT300 = RAT3 * 100. ;
  1372.           NBP% = ENTI RAT300 ;
  1373.           'MESS' 'POURCENTAGE DE POINTS DE GAUSS EN FLUAGE PRIMAIRE :'
  1374.  NBP% '%' ;
  1375.         'FINSI';
  1376.         'SI' (INDIC 'EGA' 2) ;
  1377.           'MESS' 'A316  -  PAS DE RELAXATION' ;
  1378.           NBFI3 = NBFI * 3. ;RAT3 = NBP / NBFI3 ;RAT300 = RAT3 * 100. ;
  1379.           NBP% = ENTI RAT300 ;
  1380.           'MESS' 'POURCENTAGE DE POINTS DE GAUSS EN FLUAGE PRIMAIRE :'
  1381.  NBP% '%' ;
  1382.         'FINSI';
  1383.         'SI' (INDIC 'EGA' 3) ;
  1384.           'MESS' 'MATERIAU QUELCONQUE  -  RELAXATION (SECONDAIRE)' ;
  1385.         'FINSI';
  1386.         'SI' (INDIC 'EGA' 4) ;
  1387.           'MESS' 'MATERIAU QUELCONQUE  -  PAS DE RELAXATION' ;
  1388.         'FINSI';
  1389.         'MESS' '                                              ' ;
  1390.         'MESS' 'VALEUR MAXIMUM DE D_A DONNEE POUR' ACC 'CYCLES' ;
  1391.         'MESS' '                                              ' ;
  1392.      'MESS' '   CYCLE   SS-CYCLE  SM-EFFB      DA_FAT       DA_FLU 
  1393.       ATOTAL' ;
  1394.         'MESS' '                                                ' ;
  1395.       'SINON' ;
  1396.         'SI' (KCRIT 'EGA' 1) ;
  1397.           'MESS' '                                                ' ;
  1398.      'MESS' '   CYCLE   SS-CYCLE  SM-EFFB      DA_FAT       DA_FLU 
  1399.       ATOTAL' ;
  1400.           'MESS' '                                                ' ;
  1401.         'FINSI' ;
  1402.         'SI' ((INDIC 'EGA' 1) 'OU' (INDIC 'EGA' 2)) ;
  1403.           'SI' (NBP 'NEG' NBPP) ;
  1404.             NBFI3 = NBFI * 3. ;RAT3 = NBP / NBFI3 ;RAT300 = RAT3 * 100.;
  1405.             NBP% = ENTI RAT300 ;
  1406.             'MESS' 'POURCENTAGE DE POINTS DE GAUSS EN FLUAGE PRIMAIRE :'
  1407.  NBP% '%' ;
  1408.             'MESS' '                                                ' ;
  1409.      'MESS' '   CYCLE   SS-CYCLE  SM-EFFB      DA_FAT       DA_FLU 
  1410.       ATOTAL' ;
  1411.             'MESS' '                                                ' ;
  1412.           'FINSI' ;
  1413.         'FINSI' ;
  1414.       'FINSI' ;
  1415.       NBPP = NBP ;
  1416.       ENBCYT = 'ENTI' NBCYT ;ENBCY = 'ENTI' NBCY ;
  1417.       BB = CHAIN ENBCYT*7 ENBCY*17 S_EFFBM*30
  1418.  DA_FATM*43 DA_FLUM*56 ATOTAM*69 ;
  1419.       'MESS' BB ;
  1420.       'MESS' '                                                     ' ;
  1421.  
  1422. *- TEST DE REACTUALISATION OU DE PERCEMENT
  1423.  
  1424.  
  1425.       RA = ATOTAL 'MASQ' 'SUPERIEUR' 'SOMME' ALIM ;
  1426.       'SI' (RA > 0.5) ;
  1427.         A_DEPAS = ATOTAL - ALIM ;A_DEPASM = 'MAXI' A_DEPAS ;
  1428.         DA_TOT = DA_FATC + DA_FLU ;DA_TOTM = 'MAXI' DA_TOT ;
  1429.         DA_TOTM1 = DA_TOTM ** -1 ;
  1430.         NBCYF = NBCY ;NBCYFT = NBCYT ;
  1431.         NBCYS = ACC * A_DEPASM * DA_TOTM1 ;
  1432.         NBCY = NBCY - NBCYS ;NBCYT = NBCYT - NBCYS ;
  1433.         DIMC = DIME L_NBCYT ;
  1434.         D_NBCYT = 'ENLE' L_NBCYT DIMC ;
  1435.         P_NBCYT = 'PROG' NBCYT ;
  1436.         L_NBCYT = D_NBCYT 'ET' P_NBCYT ;
  1437.         TTROP = NBCYS * TM ;
  1438.         TTOTAL = TTOTAL - TTROP ;
  1439.         DA_FATRO = NBCYS * DA_FAT ;NBAC = NBCYS / ACC ;
  1440.         DA_FLURO = NBAC * DA_FLU ;DA_TROP = DA_FATRO + DA_FLURO ;
  1441.         ATOTAL = ATOTAL - DA_TROP ;
  1442.         D_ATOT = 'ENLE' L_ATOT DIMC ;M_ATOT = 'MAXI' ATOTAL ;
  1443.         P_ATOT = 'PROG' M_ATOT ;
  1444.         L_ATOT = D_ATOT 'ET' P_ATOT ;
  1445.         ATOTAL1 = ATOTAL1 - DA_FATRO ;
  1446.         D_ATOT1 = 'ENLE' L_ATOT1 DIMC ;M_ATOT1 = 'MAXI' ATOTAL1 ;
  1447.         P_ATOT1 = 'PROG' M_ATOT1 ;
  1448.         L_ATOT1 = D_ATOT1 'ET' P_ATOT1 ;
  1449.         ATOTAL2 = ATOTAL2 - DA_FLURO ;
  1450.         D_ATOT2 = 'ENLE' L_ATOT2 DIMC ;M_ATOT2 = 'MAXI' ATOTAL2 ;
  1451.         P_ATOT2 = 'PROG' M_ATOT2 ;
  1452.         L_ATOT2 = D_ATOT2 'ET' P_ATOT2 ;
  1453.         ATOTAM = 'MAXI' ATOTAL ;
  1454.         'SI' (INDIC 'EGA' 1) ;
  1455.           'SI' (NBCY > 1.1) ;
  1456.             D_EFFBP = S_EFFBP0 - S_EFFBDP ;C_EFFBP = D_EFFBP * NBAC ;
  1457.             D_EFFBS = S_EFFBS0 - S_EFFBDS ;C_EFFBS = D_EFFBS * NBAC ;
  1458.             S_EFFBDP = S_EFFBDP + C_EFFBP ;
  1459.             S_EFFBDS = S_EFFBDS + C_EFFBS ;
  1460.           'FINSI' ;
  1461.           TTROP1 = NBCYF * TM ;TTROP2 = NBCY * TM ;
  1462.           TTROC1 = TTROP1 ** C2 ;TTROC2 = TTROP2 ** C2 ;
  1463.           DTROC = TTROC1 - TTROC2 ;S_EFBDP1 = S_EFFBDP ** N1 ;
  1464.           CORECP = C1 * DTROC * S_EFBDP1 ;CORECP = CORECP * CHP ;
  1465.           CORECS = DEFS * NBAC ;CORECS = CORECS * CHS ;
  1466.           COREC = CORECP + CORECS ;
  1467.           EF = EF - COREC ;
  1468.           S_EFFBDP = S_EFFBDP * CHP ;S_EFFBDS = S_EFFBDS * CHS ; 
  1469.           S_EFFB = S_EFFBDP + S_EFFBDS ;
  1470.         'FINSI' ;
  1471.         'SI' (INDIC 'EGA' 2) ;
  1472.           TMF = NBCYF * TM ;TMN = NBCY * TM ;
  1473.           TMF1 = TMF ** C2 ;TMN1 = TMN ** C2 ;D_TM = TMF1 - TMN1 ;
  1474.           D_NBCY = NBCYFT - NBCYT ;T_NBCY = D_NBCY * TM ;
  1475.           EFT2 = S_EFFB ** N1 ;EPS_PT1 = S_EFFB ** NEPSS ;
  1476.           CORECP = C1 * D_TM * EFT2 ;CORECP = CORECP * CHP ;
  1477.           CORECS = 100. * CEPSS * EPS_PT1 * T_NBCY ;
  1478.           CORECS = CORECS * CHS ;
  1479.           COREC = CORECP + CORECS ;
  1480.           EF = EF - COREC ;
  1481.         'FINSI' ;
  1482.         'SI' (INDIC 'EGA' 3) ;
  1483.           D_EFFBD = S_EFFB0 - S_EFFBD ;D_EFFBD = D_EFFBD * NBAC ;
  1484.           S_EFFB =  S_EFFBD + D_EFFBD ;
  1485.           COREC = DEF * NBAC ;
  1486.           EF = EF - COREC ;
  1487.         'FINSI' ;
  1488.         'SI' (INDIC 'EGA' 4) ;
  1489.           S_EFFB1 = S_EFFB ** NEPSS ;N_TROP = NBCYF - NBCY ;
  1490.           TTROP = N_TROP * TM
  1491.           COREC = 100. * CEPSS * S_EFFB1 * TTROP ;
  1492.           EF = EF - COREC ;
  1493.         'FINSI' ;
  1494.         D_EP_FLU = 'ENLE' L_EP_FLU DIMC ;
  1495.         EFTN = EF ;EFTNM = 'MAXI' EFTN ;P_EFTNM = 'PROG' EFTNM ;
  1496.         L_EP_FLU = D_EP_FLU 'ET' P_EFTNM ;
  1497.         D_L_SEFB = 'ENLE' L_SEFB DIMC ;M_S_EFFB = 'MAXI' S_EFFB ;
  1498.         P_S_EFFB = 'PROG' M_S_EFFB ;
  1499.         L_SEFB = D_L_SEFB 'ET' P_S_EFFB ;
  1500.         'QUIT' CYCLAGE ;
  1501.       'FINSI' ;
  1502.  
  1503.       'SI' (INDIC 'EGA' 1) ;
  1504.         S_EFFBP0 = S_EFFBDP ;S_EFFBS0 = S_EFFBDS ;
  1505.       'FINSI' ;
  1506.       'SI' (INDIC 'EGA' 3) ;
  1507.         S_EFFB0 = S_EFFBD ;
  1508.       'FINSI' ;
  1509.  
  1510.  
  1511. *- TEST POUR ACCELERATION DU CALCUL PAR ACCROISSEMENT DES ITERATIONS
  1512.  
  1513.       'SI' (CRIT > 0.99) ;
  1514.         KAC = KAC + 1 ;
  1515.         ACC = KAC ** 2 ;
  1516.         'MESS' '                                              ' ;
  1517.         'MESS' 'VALEUR MAXIMUM DE D_A DONNEE POUR' ACC 'CYCLES' ;
  1518.         'MESS' '                                              ' ;
  1519.         KCRIT = 0 ;
  1520.       'FINSI' ;
  1521.  
  1522.     'FIN' CYCLAGE ;
  1523.  
  1524.  
  1525.     DON.'EPSF1' = EF ;
  1526.  
  1527.     TABCON.'EPSILON' .  (KREA)     = DON.'EPSF1' ;
  1528.  
  1529.     DON.LR_ATOT1 = L_ATOT1 ;
  1530.     DON.LR_ATOT2 = L_ATOT2 ;
  1531.     DON.LR_DAFLU = L_DA_FLU ;
  1532.     DON.LR_EPFLU = L_EP_FLU ;
  1533.     DON.LR_SEFB  = L_SEFB ;
  1534.  
  1535.     DON.ATOT1 = ATOTAL1 ;
  1536.     DON.ATOT2 = ATOTAL2 ;
  1537.  
  1538.     'SI' (NBFISS 'EGA' CPTG) ;
  1539.      'MESS' '***************************************************';
  1540.      'MESS' 'LE NOMBRE D ITERATIONS DOIT ETRE AUGMENTE EN FLUAGE';
  1541.      'MESS' '> AUGMENTER LA VALEUR CONTENU SOUS l INDICE "PROPA"';
  1542.      'MESS' '***************************************************';
  1543.       'QUIT' BOUREA ;
  1544.       DON.'AMAX' = VRAI ;
  1545.     'FINSI';
  1546.  
  1547.     EPSM = 'MAXI' EF ;
  1548.     'SI' ('EGA' TRAV VRAI) ;
  1549.       'MESS'     '  CYCLE    SS-CYCLE    ATOTAL     EPSILON FLUAGE %';
  1550.       'MESS'     '  -----    --------    ------     --------------  ';
  1551.        BB = CHAIN ('ENTI' NBCYT)*7 ('ENTI' NBCY)*17 ATOTAM*31 EPSM*47 ;
  1552.       'MESS' '                                                     ';
  1553.       'MESS' BB ;
  1554.       'MESS' '                                          ' ;
  1555.       'MESS' '                                          ' ;
  1556.       'MESS' '                                          ' ;
  1557.       'MESS' '                                          ' ;
  1558.       'MESS' '**************************************' ;
  1559.       'MESS' COEFIMP '* EPAISSEUR ATTEINT'            ;
  1560.       'MESS' '**************************************' ;
  1561.       'MESS' '                          ' ;
  1562.       DON.'AMAX' = VRAI ;
  1563.     'SINON';
  1564.       'MESS' '  CYCLE    SS-CYCLE    ATOTAL     EPSILON FLUAGE %';
  1565.       'MESS' '  -----    --------    ------     --------------  ';
  1566.        BB = CHAIN ('ENTI' NBCYT)*7 ('ENTI' NBCY)*17 ATOTAM*31 EPSM*47 ;
  1567.       'MESS' '                                                     ';
  1568.       'MESS' BB ;
  1569.       'MESS' '                                          ' ;
  1570.       'MESS' '                                          ' ;
  1571.       'MESS' '                                          ' ;
  1572.       'MESS' '                                          ' ;
  1573.       'MESS' '                                          ' ;
  1574.       'MESS' '                                          ' ;
  1575.       'MESS' '                                          ' ;
  1576.       'MESS' '************************************************' ;
  1577.       'MESS' KREA 'EME ACTUALISATION DE FISSURE EFFECTIVE'   ;
  1578.       'MESS' '************************************************' ;
  1579.       'MESS' '                                             ' ;
  1580.       DON.'AMAX' = FAUX ;
  1581.     'FINSI';
  1582.  
  1583.     'SI' (('MAXI' EF) > 1.) ;
  1584.       'MESS' 'ATTENTION : LE EPSILON FLUAGE A DEPASSE 1%' ;
  1585.     'FINSI' ;
  1586.  
  1587.   'SINON';
  1588.  
  1589. *--------------------------------------------------------------*
  1590. *          CALCUL EN FATIGUE (EN L'ABSENCE DE FLUAGE)          *
  1591. *--------------------------------------------------------------*
  1592.  
  1593.     DA_STABL = DA_FAT ;
  1594.     DA_STAF = DA_STABL 'MASQ' 'INFERIEUR' 1.E-20 ;
  1595.     DA_STABL = ((CHAM1 * 1.E-20) * DA_STAF) + DA_STABL ;
  1596.     WARN = AFIN 'MASQ' 'SUPERIEUR' 'SOMME' WIMP ;
  1597.     'SI' (WARN > 0) ;
  1598.       ALIM = WIMP ;
  1599.       TRAV = VRAI ;
  1600.     'SINON';
  1601.       ALIM = AFIN ;
  1602.       TRAV = FAUX ;
  1603.     'FINSI';
  1604.     NBCYR = ('MINI' ((ALIM - ATOTAL) * (DA_STABL ** -1))) ;
  1605.     'SI' (('ENTI' NBCYR) 'EGA' 0) ;
  1606.       'MESS' 'PROPAGATION TRES IMPORTANTE' ;
  1607.       'MESS' 'IL CONVIENDRAIT D AUGMENTER LE COEFFICIENT DE 
  1608.  REACTUALISATION' ;
  1609.     'FINSI' ;
  1610.     ATOTAL = ATOTAL + (DA_STABL * NBCYR) ;
  1611. *list-reel des profondeurs maxi
  1612.     L_ATOT = L_ATOT 'ET' ('PROG' ('MAXI' ATOTAL)) ;
  1613.     NBCY = NBCYR ;
  1614.     NBCYT = NBCYT + NBCYR ;
  1615. *list-reel des numero de cycle
  1616.     L_NBCYT = L_NBCYT 'ET' ('PROG' NBCYT) ;
  1617. *list-reel des DA_FATIGUE
  1618.     L_DA_FAT = L_DA_FAT 'ET' ('PROG' ('MAXI' DA_FAT)) ;
  1619.     TTOTAL = TTOTAL + (NBCYR * TM) ;
  1620.     DA_STABM = 'MAXI' DA_STABL ;
  1621.     ATOTAM = 'MAXI' ATOTAL ;
  1622.     'SI' TRAV ;
  1623.       'MESS' '                          ' ;
  1624.       'MESS' '  CYCLE  SS-CYCLE     DA_FATIGUE     ATOTAL ' ;
  1625.       BB = CHAIN ('ENTI' NBCYT)*7 ('ENTI' NBCY)*17 DA_STABM*32
  1626.          ATOTAM*45 ;
  1627.       'MESS' BB ;
  1628.       'MESS' '                                      ' ;
  1629.       'MESS' '**************************************' ;
  1630.       'MESS' COEFIMP '* EPAISSEUR ATTEINT' ;
  1631.       'MESS' '**************************************' ;
  1632.       DON.'AMAX' = VRAI ;
  1633.     'SINON';
  1634.       'MESS' '                                             ' ;
  1635.      'MESS' '  CYCLE  SS-CYCLE     DA_FATIGUE     ATOTAL ' ;
  1636.       BB = CHAIN ('ENTI' NBCYT)*7 ('ENTI' NBCY)*17 DA_STABM*32
  1637.          ATOTAM*45 ;
  1638.       'MESS' BB ;
  1639.       'MESS' '                                           ' ;
  1640.       'MESS' '                                           ' ;
  1641.       'MESS' '                                           ' ;
  1642.       'MESS' '                                           ' ;
  1643.       'MESS' '                                           ' ;
  1644.       'MESS' '                                           ' ;
  1645.       'MESS' '************************************************' ;
  1646.       'MESS' KREA 'EME ACTUALISATION DE FISSURE EFFECTIVE'      ;
  1647.       'MESS' '************************************************' ;
  1648.       DON.'AMAX' = FAUX ;
  1649.     'FINSI';
  1650.   'FINSI';
  1651.  
  1652.   DON.LR_NBCYT = L_NBCYT ;
  1653.   DON.LR_DAFAT = L_DA_FAT ;
  1654.   DON.LR_ATOT  = L_ATOT ;
  1655.  
  1656.   DON.ATOT  = ATOTAL ;
  1657.  
  1658. *** --- NOUVEL OBJET CARACTERISTIQUE
  1659.  
  1660.   'MESS' '                                              ' ;
  1661.   'MESS' 'CONSTRUCTION D UN NOUVEL OBJET CARACTERISTIQUE' ;
  1662.   'MESS' '                                              ' ;
  1663.   DON.'NREA' = KREA ;
  1664.   KF = 0 ;ATOTALE = CHAM0 ;ATOTALE1 = CHAM0 ;ATOTALE2 = CHAM0 ;
  1665.   'REPETER' BOUCF NBFI ;
  1666.     KF = KF + 1 ;
  1667.     ELC = FISS1 'ELEM' KF ;
  1668.     CARFE = 'REDU' CARF ELC ;
  1669.     EP = 'EXTR' CARFE 'EPAI' 1 1 2 ;
  1670.     ATOTR = 'REDU' ATOTAL ELC ;
  1671.     PROF1 = 'EXTR' ATOTR 'SCAL' 1 1 1 ;
  1672.     PROF2 = 'EXTR' ATOTR 'SCAL' 1 1 2 ;
  1673.     PROF3 = 'EXTR' ATOTR 'SCAL' 1 1 3 ;
  1674.     PROFM = (PROF1 + PROF2 + PROF3) / 3. ;
  1675.     'SI' (PROFM > EP) ;
  1676.       PROFM = EP - 1.E-8 ;
  1677.     'FINSI';
  1678.     'SI' (CALFAT  'EGA' 0) ;
  1679.       ATOTR1 = 'REDU' ATOTAL1 ELC ;
  1680.       PROF11 = 'EXTR' ATOTR1 'SCAL' 1 1 1 ;
  1681.       PROF12 = 'EXTR' ATOTR1 'SCAL' 1 1 2 ;
  1682.       PROF13 = 'EXTR' ATOTR1 'SCAL' 1 1 3 ;
  1683.       PROFM1 = (PROF11 + PROF12 + PROF13) / 3. ;
  1684.       ATOTE1 = SELEC.KF * PROFM1 ;
  1685.       ATOTALE1 = ATOTALE1 + ATOTE1 ;
  1686.       ATOTR2 = 'REDU' ATOTAL2 ELC ;
  1687.       PROF21 = 'EXTR' ATOTR2 'SCAL' 1 1 1 ;
  1688.       PROF22 = 'EXTR' ATOTR2 'SCAL' 1 1 2 ;
  1689.       PROF23 = 'EXTR' ATOTR2 'SCAL' 1 1 3 ;
  1690.       PROFM2 = (PROF21 + PROF22 + PROF23) / 3. ;
  1691.       ATOTE2 = SELEC.KF * PROFM2 ;
  1692.       ATOTALE2 = ATOTALE2 + ATOTE2 ;
  1693.     'FINSI' ;
  1694.     ATOTE = SELEC.KF * PROFM ;
  1695.     ATOTALE = ATOTALE + ATOTE ;
  1696.   'FIN' BOUCF ;
  1697.   CHFIS0 = CARF - ('EXCO' CARF 'FISS') ;
  1698.   'SI' (CALFAT  'EGA' 0) ;
  1699.     CHFISS1 = 'MANU' 'CHML' OBJF 'FISS' ATOTALE1 'TYPE'
  1700.              'CARACTERISTIQUES' 'RIGIDITE' ;
  1701.     CAR11 = CHFIS0 + CHFISS1 ;
  1702.     DON.'CAR1' = CAR11 ;
  1703.     TABCAR.'FATIGUE'. KREA = CAR11 ;
  1704.     CHFISS2 = 'MANU' 'CHML' OBJF 'FISS' ATOTALE2 'TYPE'
  1705.              'CARACTERISTIQUES' 'RIGIDITE' ;
  1706.     CAR12 = CHFIS0 + CHFISS2 ;
  1707.     DON.'CAR2' = CAR12 ;
  1708.     TABCAR.'FLUAGE'. KREA = CAR12 ;
  1709.   'FINSI' ;
  1710.   CHFISS = 'MANU' 'CHML' OBJF 'FISS' ATOTALE 'TYPE'
  1711.            'CARACTERISTIQUES' 'RIGIDITE' ;
  1712.   CAR1 = CHFIS0 + CHFISS ;
  1713.   DON.'TPSTOT' = TTOTAL ;
  1714.   DON.'NBCYCLE' = NBCYT ;
  1715.   DON.'CAR' = CAR1 ; 
  1716.   DON.'MAT' = MATF ;
  1717.   DON.'OBJ' = OBJF ;
  1718.   DON.'REPRISE' = VRAI ;
  1719.  
  1720.   TABCAR.'TOTAL'. KREA = CAR1 ;
  1721.   TABNCY.KREA = NBCY ;
  1722.  
  1723.   TABCON.'DELTSIG'.   (KREA - 1) = DSIGFIS ;
  1724.   TABCON.'CHSIG_A'.   (KREA - 1) = SIGAFIS ;
  1725.   TABCON.'CHSIG_B'.   (KREA - 1) = SIGBFIS ;
  1726.   TABCON.'DELTSIG1'.  (KREA - 1) = DSIG1FIS ;
  1727.   TABCON.'CHSIG_A1'.  (KREA - 1) = SIGA1FIS ;
  1728.   TABCON.'CHSIG_B1'.  (KREA - 1) = SIGB1FIS ;
  1729.   TABCON.'SDEFRAS'.   (KREA - 1) = S_DEFRAS ;
  1730.   TABCON.'SDEFRBS'.   (KREA - 1) = S_DEFRBS ;
  1731.   TABCON.'SDEFR'.     (KREA - 1) = DS_DEFR ;
  1732.   TABCON.'VALR' .     (KREA - 1) = R ;
  1733.   TABCON.'VALQ' .     (KREA - 1) = Q ;
  1734.   TABCON.'VALPSI' .   (KREA - 1) = PSI ;
  1735.   TABCON.'VALKDJ' .   (KREA - 1) = K_DJ ;
  1736.   TABCON.'VALDJ' .    (KREA - 1) = DJ ;
  1737.   TABCON.'VALDKEFF' . (KREA - 1) = DK_EFF ;
  1738.   TABCON.'VALDAFAT' . (KREA - 1) = DA_FAT ;
  1739.  
  1740.   DON.'EVNCY' = TABNCY ;
  1741.   DON.'EVCAR' = TABCAR ;
  1742.   DON.'EVCON' = TABCON ;
  1743.  
  1744.  
  1745.  
  1746. 'FIN' BOUREA ;
  1747. *
  1748. *
  1749. *--------------------------------------------------------------*
  1750. *                           TRACES                             *
  1751. *--------------------------------------------------------------*
  1752. *
  1753. *--- Maxi des champs CHSIG-B
  1754. *
  1755.   nn = don.lr_nbcyt ;
  1756.   ff = done.lr_sefb ; 
  1757.   titr 'EVOLUTION DE LA CONTRAINTE EFFICACE (B)' ;
  1758.   EVCHSIGB = EVOL MANU 'CYCLES' nn 'SigB (MPa)' ff ;
  1759.   DESS EVCHSIGB ;
  1760. *
  1761. *--- da maxi de fluage
  1762. *
  1763.   fu = don.lr_daflu ;
  1764.   titr 'EVOLUTION DU DA_FLUAGE MAXI' ; 
  1765.   EVDAFLU = EVOL MANU 'CYCLES' nn 'dAflu (mm/cy)' fu ;
  1766.   DESS EVDAFLU ;
  1767. *
  1768. *--- profondeur maxi
  1769. *
  1770.   tatot = don.lr_atot ;
  1771.   titr 'EVOLUTION DE LA PROFONDEUR SOUS CONTRIBUTION FATIGUE/FLUAGE' ;  
  1772.   EVDAATOT = EVOL MANU 'CYCLES' nn 'dAmax (mm/cy)' tatot ;
  1773.   DESS EVDAATOT ;
  1774.  
  1775. *--- epsfluage cumulé maxi
  1776. *
  1777.   eps = don.lr_epflu ;
  1778.   titr 'EVOLUTION DU EPSILON FLUAGE MAXI' ;
  1779.   EVEPFLU = EVOL MANU 'CYCLES' nn 'EPSflu %' eps ;
  1780.   DESS EVEPFLU ;
  1781.  
  1782. *  MENAGE ;
  1783.  
  1784. 'FINP' DON ; 
  1785.  

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