Télécharger mazars_traction_cyclique.dgibi

Retour à la liste

* fichier : mazars_traction_cyclique.dgibi
************************************************************************
*      Cas test elementaire pour la loi d'endommagement de Mazars      *
*                  Essai de traction simple cyclique                   *
*                                                                      *
* Le test est realise avec plusieurs modes de calcul :                 *
*   1) 3D volumique                                                    *
*   2) 3D poutre a fibre                                               *
*   3) 2D contraintes planes                                           *
*   4) 2D axisymetrique                                                *
*   5) 2D poutre a fibre                                               *
*                                                                      *
* ITRAC1 : VRAI => affichage resultats                                 *
* IPSC1  : VRAI => affichage resultats dans un PostScript Couleur      *
************************************************************************
 
** Options pour les traces
ITRAC1 = FAUX ;
IPSC1  = VRAI ;
SI IPSC1 ;
  OPTI 'TRAC' 'PSC' ;
FINSI ;
 
** Sortie sur Erreur 1152 "Solution calculee non conforme par rapport a la solution de reference"
IERR = FAUX ;
 
** Liste des modes de calcul a traiter
lcas = LECT 1 2 3 4 5 ;
*nh145313
*lcas = LECT 1 ; COMM '3D volumique' ;
*lcas = LECT 2 ; COMM '3D poutre a fibre' ;
*lcas = LECT 3 ; COMM '2D contraintes planes' ;
*lcas = LECT 4 ; COMM '2D axisymetrique' ;
*lcas = LECT 5 ; COMM '2D poutre a fibre' ;
 
** Parametres geometriques
long = 1. ;
epai = 0.1 ;
 
** Parametres materiau
you   = 20.E9 ;
nu    = 0.2 ;
epd0  = 8.E-5 ;
at    = 1. ;
ac    = 1.5 ;
bt    = 8000. ;
bc    = 1550. ;
beta  = 1. ;
 
** Paramètres chargement
umax  = 5.E-4 ;
 
*nh145313
** Legende courbes
  tdess = 'TABLE' ;
  tdess.'TITRE' = 'TABLE' ;
  tdess. 1 = 'MARQ LOSA' ;
  tdess.'TITRE'. 1 = MOT 'Calcul' ;
  tdess. 2 = 'MARQ CROI' ;
  tdess.'TITRE'. 2 = MOT 'Reference' ;
 
 
 
 
 
**************** C A S  #1  :  3 D   V O L U M I Q U E ****************
SI (EXIS lcas 1) ;
 
** Options generales
  OPTI 'DIME' 3 'MODE' 'TRID' 'ELEM' 'CUB8' ;
 
** Maillage
  p1      = 0. 0. 0. ;
  p2      = 0. long 0. ;
  l12     = DROI 1 p1 p2 ;
  sgauche = l12 TRAN 1 (0. 0. long) ;
  mail    = sgauche VOLU 'TRAN' 1 (long 0. 0.) ;
  sdroite = mail FACE 2 ;
  p3      = sdroite POIN 'PROC' (long 0. 0.) ;
  mes1    = MESU mail ;
*nh145313
cadr1 = COUL ((0. 1.05 1.) ET (1.15 1.05 1.)) 'BLAN' ;
rep1 = @REPERE (-0.55 -0.55 0.) 'TURQ' ;
SI ITRAC1 ;
  TRAC 'QUAL' (cadr1 ET rep1 ET mail ET (sgauche COUL 'VERT') ET (sdroite COUL 'ROUG')) 'TITR' '[3D V] Maillage' ;
FINSI ;
 
** Modele et caracteristiques materiau
  mo = MODE mail 'MECANIQUE' 'ELASTIQUE' 'ENDOMMAGEMENT' 'MAZARS' ;
  ma = MATE mo 'YOUN' you 'NU' nu 'KTR0' epd0 'ACOM' ac 'BCOM' bc 'ATRA' at 'BTRA' bt 'BETA' beta ;
 
** Blocages
  blgauche = BLOQ 'UX' sgauche ;
  bldroite = BLOQ 'UX' sdroite ;
  blrig    = (BLOQ 'UY' 'UZ' p1) ET (BLOQ 'UZ' p2) ;
  bl       = blgauche ET bldroite ET blrig ;
 
** Chargements
  didroite = DEPI bldroite umax ;
  ec1      = EVOL 'MANU' (PROG 0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.)
                         (PROG 0. 0.2 0.  0.4 0.  0.6 0.  0.8 0.  1.  0.) ;
  cha      = CHAR 'DIMP' didroite ec1 ;
 
** Tracé des blocages et du chargement
  re0 = @REPERE (-0.55 -0.55 0.) ;
  mor = MODE re0 'MECANIQUE' 'BARR' ;
  mar = MATE mor 'YOUN' 1. 'NU' 1. 'SECT' 1. ;
  rir = RIGI mor mar ;
  vf  = VECT didroite 'FLX' ' ' ' ' 'ROUGE' ;
  SI ITRAC1 ;
    TRAC ((RIGI mo ma) ET blgauche ET blrig ET rir) vf 'NOLE' 'TITR' '[3D V] Blocages et chargement' ;
  FINSI ;
 
** Resolution
  t                               = TABL ;
  t . 'MODELE'                    = mo ;
  t . 'CARACTERISTIQUES'          = ma ;
  t . 'BLOCAGES_MECANIQUES'       = bl ;
  t . 'CHARGEMENT'                = cha ;
  t . 'TEMPS_CALCULES'            = PROG 0. 'PAS' 0.01 0.1
                                            'PAS' 0.01 0.3
 
*                                           'PAS' 0.01 0.5
                                            'PAS' 0.01 0.46
                                            'PAS' 0.002 0.466
                                            'PAS' 0.00066 0.468
                                            'PAS' 0.002 0.47
 
*                                           'PAS' 0.01 0.7
                                            'PAS' 0.01 0.65
                                            'PAS' 0.005 0.75
 
                                            'PAS' 0.01 0.9
                                            'PAS' 0.01 1.0 ;
  t . 'MOVA'                      = MOT 'D' ;
  t . 'MES_SAUVEGARDES'           = TABL ;
  t . 'MES_SAUVEGARDES' . 'DEFTO' = VRAI ;
  PASAPAS t ;
 
** Post traitement
* courbes : endommagement moyen vs temps
*           contrainte moyenne vs deformation moyenne
*           force de reaction vs deplacement impose
  tt = t . 'TEMPS' ;
  tc = t . 'CONTRAINTES' ;
  te = t . 'DEFORMATIONS' ;
  tv = t . 'VARIABLES_INTERNES' ;
  tu = t . 'DEPLACEMENTS' ;
  tr = t . 'REACTIONS' ;
  def0 = DEFO (cadr1 ET rep1 ET mail) (tu . 0) 0. ;
  ltps = PROG ;
  lend = PROG ;
  lsig = PROG ;
  leps = PROG ;
  ldep = PROG ;
  lrea = PROG ;
  REPE b1 (DIME tt) ;
    tps1 = tt . (&b1 - 1) ;
    ltps = ltps ET tps1 ;
    end1 = (INTG mo (tv . (&b1 - 1)) 'D')    / mes1 ;
    sig1 = (INTG mo (tc . (&b1 - 1)) 'SMXX') / mes1 ;
    eps1 = (INTG mo (te . (&b1 - 1)) 'EPXX') / mes1 ;
    u1   = EXTR (tu . (&b1 - 1)) 'UX' p3 ;
    rea1 = 0. ;
    SI (NEG &b1 1) ;
      rea1 = @TOTAL (tr . (&b1 - 1)) sgauche 'FX' ;
    FINSI ;
    lend = lend ET end1 ;
    lsig = lsig ET sig1 ;
    leps = leps ET eps1 ;
    ldep = ldep ET u1 ;
    lrea = lrea ET rea1 ;
  FIN b1 ;
  ev0 = EVOL 'ROUG' 'MANU' 'Temps' ltps 'D'    lend ;
  ev1 = EVOL 'ROUG' 'MANU' 'EPXX'  leps 'SMXX' lsig ;
  ev2 = EVOL 'ROUG' 'MANU' 'UX'    ldep 'FX'   lrea ;
*nh145313
  ltps1 =  ENLE ltps 1 ;
  nb1 = DIME ltps1 ;
  un1 = PROG nb1 * 1.D0 ;
  lepd0 = PROG nb1 * epd0 ;
  lat1 = PROG nb1 * at ;
  lep1 = IPOL ltps1 ec1 ;
  lepe = lep1 * umax / long ;
*mess 'lepe avant seuillage 0'; list lepe;
  seuil0 = lepe MASQ 'EGAL' 0 ;
*mess 'test seuil0';list seuil0;
  lepe0 = seuil0 * 1.E-20 ;
*mess 'lepe0';list lepe0;
  lepe = lepe + lepe0 ;
*mess 'lepe apres seuillage 0'; list lepe;
  dt = un1 - ((epd0*(un1-lat1)/lepe) + (at*(EXP(bt*(lepd0-lepe))))) ;
  seuil1 = dt MASQ 'EGSUP' 0. ;
  dt = seuil1 * dt ; 
* pas de décroissance du dommage
  dtp = 0.;
  dt2 = PROG ;
  REPE bdt (DIME dt) ;
    dti = EXTR dt &bdt ;
    ddc = dti - dtp ;
    SI (ddc < 0.) ;
      dti = dtp ;
    FINS ;
    dt2 = dt2 ET (PROG dti) ;
    dtp = dti ;
  FIN bdt ;
  dt = dt2 ;
  lsig1 = (un1 - dt) * you * lepe ;
* contrainte nulle = 1.E-20
  seuilsc = lsig MASQ 'EGINF' 1.E-6 ;
  sc0 = seuilsc * 1.E-20 ;
  lsig = (lsig * (lsig MASQ 'EGSUP' 1.E-6)) + sc0 ;
  seuilsr = lsig1 MASQ 'EGINF' 1.E-6 ;
  sr0 = seuilsr * 1.E-20 ;
  lsig1 = (lsig1 * (lsig1 MASQ 'EGSUP' 1.E-6)) + sr0 ;
*mess 'Sig.calcul';list (ENLE lsig 1);
*mess 'Sig.Ref.';list lsig1;
  an0 = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'D'    dt    ;
  an1 = EVOL 'TURQ' 'MANU' 'EPXX'  lepe  'SMXX' lsig1 ;
  SI ITRAC1 ;
    DESS (ev0 ET an0) 'TITR' '[3D V] Endommagement vs Temps' 'LEGE' 'SE' tdess ;
    DESS (ev1 ET an1) 'XBOR' 0. 6.E-4 'TITR' '[3D V] Contrainte vs Deformation' 'LEGE' 'NE' tdess ;
    DESS ev2 'TITR' '[3D V] Reaction vs Deplacement' 'LEGE' 'SE' tdess ;
  FINSI ;
* evolutions temporelles a tous les points de Gauss
  ng = 8 ;
  tleg = TABL ;
  tleg . 1 = MOT 'MARQ CARR NOLI' ;
  tleg . 2 = MOT 'MARQ LOSA NOLI' ;
  tleg . 3 = MOT 'MARQ ROND NOLI' ;
  tleg . 4 = MOT 'MARQ ETOI NOLI' ;
  tleg . 5 = MOT 'MARQ TRID NOLI' ;
  tleg . 6 = MOT 'MARQ TRIU NOLI' ;
  tleg . 7 = MOT 'MARQ TRIL NOLI' ;
  tleg . 8 = MOT 'MARQ TRIR NOLI' ;
  tleg . 'TITRE' = TABL ;
  REPE b ng ;
    tleg . 'TITRE' . &b = CHAI 'Calc_Pt_Gauss' ' ' &b ;
  FIN b ;
* endommagement
  evend = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evend = evend ET (EVOL 'ROSE' 'TEMP' t 'VARIABLES_INTERNES' 'D' 1 1 &b) ;
  FIN b ;
*nh145313
  tleg . (ng+1) = MOT 'MARQ CROI' ;
  tleg . 'TITRE' . (ng+1) = MOT 'Reference' ;
  SI ITRAC1 ;
    DESS (evend ET an0) 'LEGE' 'SE' tleg 'TITR' '[3D V] D vs Temps' ;
  FINSI ;
* contrainte
  evsmxx = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evsmxx = evsmxx ET (EVOL 'BLEU' 'TEMP' t 'CONTRAINTES' 'SMXX' 1 1 &b) ;
  FIN b ;
*nh145313
  ansmxx = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'SMXX' lsig1 ;
  SI ITRAC1 ;
    DESS (evsmxx ET ansmxx) 'LEGE' 'NE' tleg 'TITR' '[3D V] SMXX vs Temps' ;
  FINSI ;
* deformation
  evepxx = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evepxx = evepxx ET (EVOL 'VERT' 'TEMP' t 'DEFORMATIONS' 'EPXX' 1 1 &b) ;
  FIN b ;
  SI ITRAC1 ;
    DESS evepxx 'LEGE' 'NO' tleg 'TITR' '[3D V] EPXX vs Temps' ;
  FINSI ;
* deformee
  REPE b2 (DIME tt) ;
    tps1 = tt . (&b2 - 1) ;
*   trace de la deformee (pour controle visuel)
    def1 = DEFO mail (tu . (&b2 - 1)) 333. 'ROUG' ;
    tit1 = CHAI '[3D V] Deformee au temps' ' ' tps1 ;
    SI ITRAC1 ;
*      TRAC (def0 ET def1) 'TITR' tit1 ;
    FINSI ;
  FIN b2 ;
  SI ITRAC1 ;
    TRAC (def0 ET def1) 'TITR' tit1 ;
  FINSI ;
* comparaison Calcul / Reference sur la contrainte moyenne
  ler1 = ABS ((ENLE lsig 1) - lsig1) / lsig1 ;
  mer1 = MAXI ler1 ;
  eer1 = EVOL 'ORAN' 'MANU' 'EPXX'  lepe  'Ec.rel.SMXX' ler1 ;
  tit1 = CHAI '[3D V] Ecart relatif Calcul/Reference Contrainte vs Deformation, max' ' ' mer1 ;
  SI ITRAC1 ;
    DESS eer1 'TITR' tit1 ;
  FINSI ;
  OPTI ECHO 0 ;
  SAUT 1 LIGN ;
  MESS ' **************** C A S  #1  :  3 D   V O L U M I Q U E ****************';
  SAUT 1 LIGN ;
  SI (mer1 < 1.E-10) ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' < 1.E-10 => OK' ;
    SAUT 1 LIGN ;
  SINO ;
    MESS ' Val.abs. Ecart relatif max entre Calcul et Reference' ' ' mer1 ' > 1.E-10 => ECHEC' ;
    IERR = VRAI ;
    SAUT 1 LIGN ;
  FINS ;
  SAUT 1 LIGN ;
  OPTI ECHO 1 ;
FINSI ;
 
 
 
 
 
*********** C A S  #2  :  3 D   P O U T R E   A   F I B R E ************
SI (EXIS lcas 2) ;
 
** Options generales
  OPTI 'DIME' 3 'MODE' 'TRID' 'ELEM' 'CUB8' ;
 
** Maillage
  p1s     = (-0.5 * epai) (-0.5 * epai) 0. ;
  p2s     = (-0.5 * epai) ( 0.5 * epai) 0. ;
  l12s    = DROI 1 p1s p2s ;
  mails   = l12s TRAN 1 (epai 0. 0.) ;
  p1      = 0. 0. 0. ;
  p2      = long 0. 0. ;
  mail    = DROI 1 p1 p2 ;
  mes1    = MESU mails ;
*nh145313
cadr1 = COUL ((-0.1 0. 0.) ET (1.2 0. 0.)) 'BLAN' ;
rep1 = @REPERE (-0.2 -0.2 0.) (PROG 0.2 0.2 0.2) 'TURQ' ;
SI ITRAC1 ;
  TRAC 'QUAL' (cadr1 ET rep1 ET (mails COUL 'VERT') ET (mail COUL 'ROUG')) 'TITR' '[3D PaF] Maillages de la Poutre et de sa Section' ;
FINSI ;
 
 
** Modele et caracteristiques materiau
  mos = MODE mails 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'MAZARS' 'QUAS' ;
  mas = MATE mos 'YOUN' you 'NU' nu 'KTR0' epd0 'ACOM' ac 'BCOM' bc 'ATRA' at 'BTRA' bt 'BETA' beta
                 'ALPY' 0.66 'ALPZ' 0.66 ;
  mo = MODE mail 'MECANIQUE' 'ELASTIQUE' 'SECTION' 'PLASTIQUE' 'SECTION' 'TIMO' ;
  ma = MATE mo 'MODS' mos 'MATS' mas 'VECT' (0. 1. 0.) ;
 
** Blocages
  blgauche = BLOQ 'DEPL' 'ROTA' p1 ;
  bldroite = BLOQ 'UX' p2 ;
  bl       = blgauche ET bldroite ;
*nh145313
  bldroit2 = BLOQ 'ROTA' p2 ;
  bl       = bl ET bldroit2 ;
 
** Chargements
  didroite = DEPI bldroite umax ;
  ec1      = EVOL 'MANU' (PROG 0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.)
                         (PROG 0. 0.2 0.  0.4 0.  0.6 0.  0.8 0.  1.  0.) ;
  cha      = CHAR 'DIMP' didroite ec1 ;
 
** Tracé des blocages et du chargement
  re0 = @REPERE (-0.2 -0.2 0.) (PROG 0.2 0.2 0.2) ;
  mor = MODE re0 'MECANIQUE' 'BARR' ;
  mar = MATE mor 'YOUN' 1. 'NU' 1. 'SECT' 1. ;
  rir = RIGI mor mar ;
  vf = VECT didroite 'FLX' ' ' ' ' 'ROUG' 200. ;
  SI ITRAC1 ;
    TRAC ((RIGI mo ma) ET blgauche ET rir) vf 'NOLE' 'TITR' '[3D PaF] Blocages et chargement' ;
  FINSI ;
 
** Resolution
  t                               = TABL ;
  t . 'MODELE'                    = mo ;
  t . 'CARACTERISTIQUES'          = ma ;
  t . 'BLOCAGES_MECANIQUES'       = bl ;
  t . 'CHARGEMENT'                = cha ;
  t . 'TEMPS_CALCULES'            = PROG 0. 'PAS' 0.01 0.1
                                            'PAS' 0.01 0.3
 
*                                           'PAS' 0.01 0.5
                                            'PAS' 0.01 0.46
                                            'PAS' 0.002 0.466
                                            'PAS' 0.00066 0.468
                                            'PAS' 0.002 0.47
 
*                                           'PAS' 0.01 0.7
                                            'PAS' 0.01 0.65
                                            'PAS' 0.005 0.75
 
                                            'PAS' 0.01 0.9
                                            'PAS' 0.01 1.0 ;
  t . 'MES_SAUVEGARDES'           = TABL ;
  t . 'MES_SAUVEGARDES' . 'DEFTO' = VRAI ;
  PASAPAS t ;
 
** Post traitement
* courbes : endommagement moyen vs temps
*           contrainte moyenne vs deformation moyenne
*           force de reaction vs deplacement impose
* evolutions temporelles a tous les points de Gauss
  ng = 4 ;
  tleg = TABL ;
  tleg . 1 = MOT 'MARQ CARR NOLI' ;
  tleg . 2 = MOT 'MARQ LOSA NOLI' ;
  tleg . 3 = MOT 'MARQ ROND NOLI' ;
  tleg . 4 = MOT 'MARQ ETOI NOLI' ;
  tleg . 'TITRE' = TABL ;
  REPE b ng ;
    tleg . 'TITRE' . &b = CHAI 'Calc_Pt_Gauss' ' ' &b ;
  FIN b ;
*nh145313
  tleg . (ng+1) = MOT 'MARQ CROI' ;
  tleg . 'TITRE' . (ng+1) = MOT 'Reference' ;
  tt = t . 'TEMPS' ;
  tv = t . 'VARIABLES_INTERNES' ;
  tu = t . 'DEPLACEMENTS' ;
  tr = t . 'REACTIONS' ;
  def0 = DEFO (cadr1 ET rep1 ET mail) (tu . 0) 0. ;
  ltps = PROG ;
  lend1 = PROG ;
  lend2 = PROG ;
  lend3 = PROG ;
  lend4 = PROG ;
  lsig1 = PROG ;
  lsig2 = PROG ;
  lsig3 = PROG ;
  lsig4 = PROG ;
  leps1 = PROG ;
  leps2 = PROG ;
  leps3 = PROG ;
  leps4 = PROG ;
  ldep = PROG ;
  lrea = PROG ;
  REPE b1 (DIME tt) ;
    tps1 = tt . (&b1 - 1) ;
    ltps = ltps ET tps1 ;
    SI (EGA &b1 1) ;
      lend1 = lend1 ET 0. ;
      lend2 = lend2 ET 0. ;
      lend3 = lend3 ET 0. ;
      lend4 = lend4 ET 0. ;
      lsig1 = lsig1 ET 0. ;
      lsig2 = lsig2 ET 0. ;
      lsig3 = lsig3 ET 0. ;
      lsig4 = lsig4 ET 0. ;
      leps1 = leps1 ET 0. ;
      leps2 = leps2 ET 0. ;
      leps3 = leps3 ET 0. ;
      leps4 = leps4 ET 0. ;
    SINON ;
      lend1 = lend1 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 1) ;
      lend2 = lend2 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 2) ;
      lend3 = lend3 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 3) ;
      lend4 = lend4 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 4) ;
      lsig1 = lsig1 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 1) ;
      lsig2 = lsig2 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 2) ;
      lsig3 = lsig3 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 3) ;
      lsig4 = lsig4 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 4) ;
      leps1 = leps1 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 1) ;
      leps2 = leps2 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 2) ;
      leps3 = leps3 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 3) ;
      leps4 = leps4 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 4) ;
    FINSI ;
    u1   = EXTR (tu . (&b1 - 1)) 'UX' p2 ;
    rea1 = EXTR (tr . (&b1 - 1)) 'FX' p1 ;
    ldep = ldep ET u1 ;
    lrea = lrea ET rea1 ;
  FIN b1 ;
  ev0 = (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend1) ET (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend2) ET
        (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend3) ET (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend4) ;
  ev1 = (EVOL 'ROUG' 'MANU' 'EPXX'  leps1 'SMXX' lsig1) ET (EVOL 'ROUG' 'MANU' 'EPXX' leps2 'SMXX' lsig2) ET
        (EVOL 'ROUG' 'MANU' 'EPXX'  leps3 'SMXX' lsig3) ET (EVOL 'ROUG' 'MANU' 'EPXX' leps4 'SMXX' lsig4) ;
  ev2 = EVOL 'ROUG' 'MANU' 'UX' ldep 'FX' lrea ;
*nh145313
  ltps1 =  ENLE ltps 1 ;
  nb1 = DIME ltps1 ;
  un1 = PROG nb1 * 1.D0 ;
  lepd0 = PROG nb1 * epd0 ;
  lat1 = PROG nb1 * at ;
  lep1 = IPOL ltps1 ec1 ;
  lepe = lep1 * umax / long ;
*mess 'lepe avant seuillage 0'; list lepe;
  seuil0 = lepe MASQ 'EGAL' 0 ;
*mess 'test seuil0';list seuil0;
  lepe0 = seuil0 * 1.E-20 ;
*mess 'lepe0';list lepe0;
  lepe = lepe + lepe0 ;
*mess 'lepe apres seuillage 0'; list lepe;
  dt = un1 - ((epd0*(un1-lat1)/lepe) + (at*(EXP(bt*(lepd0-lepe))))) ;
  seuil1 = dt MASQ 'EGSUP' 0. ;
  dt = seuil1 * dt ; 
* pas de décroissance du dommage
  dtp = 0.;
  dt2 = PROG ;
  REPE bdt (DIME dt) ;
    dti = EXTR dt &bdt ;
    ddc = dti - dtp ;
    SI (ddc < 0.) ;
      dti = dtp ;
    FINS ;
    dt2 = dt2 ET (PROG dti) ;
    dtp = dti ;
  FIN bdt ;
  dt = dt2 ;
  lsan1 = (un1 - dt) * you * lepe ;
* contrainte nulle = 1.E-20
  seuilsr = lsan1 MASQ 'EGINF' 1.E-6 ;
  sr0 = seuilsr * 1.E-20 ;
  lsan1 = (lsan1 * (lsan1 MASQ 'EGSUP' 1.E-6)) + sr0 ;
*mess 'Sig.Ref.';list lsan1;
  an0 = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'D'    dt    ;
  an1 = EVOL 'TURQ' 'MANU' 'EPXX'  lepe  'SMXX' lsan1 ;
  SI ITRAC1 ;
    DESS (ev0 ET an0) 'LEGE' 'SE' tleg 'TITR' '[3D PaF] Endommagement vs Temps' ;
    DESS (ev1 ET an1) 'LEGE' 'NE' tleg 'TITR' '[3D PaF] Contrainte vs Deformation' ;
    DESS ev2 'TITR' '[3D PaF] Reaction vs Deplacement' 'LEGE' 'SE' tdess ;
  FINSI ;
* deformee
  REPE b2 (DIME tt) ;
    tps1 = tt . (&b2 - 1) ;
*   trace de la deformee (pour controle visuel)
    def1 = DEFO mail (tu . (&b2 - 1)) 333. 'ROUG' ;
    tit1 = CHAI '[3D PaF] Deformee au temps' ' ' tps1 ;
    SI ITRAC1 ;
*      TRAC (def0 ET def1) 'TITR' tit1 ;
    FINSI ;
  FIN b2 ;
  SI ITRAC1 ;
    TRAC (def0 ET def1) 'TITR' tit1 ;
  FINSI ;
* comparaison Calcul / Reference sur la contrainte moyenne
  lsig = (lsig1 + lsig2 + lsig3 + lsig4) / 4. ;
* contrainte nulle = 1.E-20
  seuilsc = lsig MASQ 'EGINF' 1.E-6 ;
  sc0 = seuilsc * 1.E-20 ;
  lsig = (lsig * (lsig MASQ 'EGSUP' 1.E-6)) + sc0 ;
*mess 'Sig.calcul';list (ENLE lsig 1);
  ler1 = ABS ((ENLE lsig 1) - lsan1) / lsan1 ;
  mer1 = MAXI ler1 ;
  eer1 = EVOL 'ORAN' 'MANU' 'EPXX'  lepe  'Ec.rel.SMXX' ler1 ;
  tit1 = CHAI '[3D PaF] Ecart relatif Calcul/Reference Contrainte vs Deformation, max' ' ' mer1 ;
  SI ITRAC1 ;
    DESS eer1 'TITR' tit1 ;
  FINSI ;
  OPTI ECHO 0 ;
  SAUT 1 LIGN ;
  MESS ' *********** C A S  #2  :  3 D   P O U T R E   A   F I B R E ************';
  SAUT 1 LIGN ;
  SI (mer1 < 1.E-10) ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' < 1.E-10 => OK' ;
    SAUT 1 LIGN ;
  SINO ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' > 1.E-10 => ECHEC' ;
    IERR = VRAI ;
    SAUT 1 LIGN ;
  FINS ;
  SAUT 1 LIGN ;
  OPTI ECHO 1 ;
FINSI ;
 
 
 
 
 
******* C A S  #3  :  2 D   C O N T R A I N T E S   P L A N E S ********
SI (EXIS lcas 3) ;
 
** Options generales
  OPTI 'DIME' 2 'MODE' 'PLAN' 'CONT' 'ELEM' 'QUA4' ;
 
** Maillage
  p1      = 0. 0. ;
  p2      = 0. long ;
  sgauche = DROI 1 p1 p2 ;
  mail    = sgauche TRAN 1 (long 0.) ;
  sdroite = mail COTE 3 ;
  p3      = sdroite POIN 'PROC' (long 0.) ;
  mes1    = MESU mail ;
*nh145313
cadr1 = COUL ((0. 1.05) ET (1.2 1.05)) 'BLAN' ;
rep1 = @REPERE (-0.2 -0.2) (PROG 0.5 0.5) 'TURQ' ;
SI ITRAC1 ;
  TRAC 'QUAL' (cadr1 ET rep1 ET mail ET (sgauche COUL 'VERT') ET (sdroite COUL 'ROUG')) 'TITR' '[2D Plan cont.] Maillage' ;
FINSI ;
 
** Modele et caracteristiques materiau
  mo = MODE mail 'MECANIQUE' 'ELASTIQUE' 'ENDOMMAGEMENT' 'MAZARS' ;
  ma = MATE mo 'YOUN' you 'NU' nu 'KTR0' epd0 'ACOM' ac 'BCOM' bc 'ATRA' at 'BTRA' bt 'BETA' beta 'DIM3' epai ;
 
** Blocages
  blgauche = BLOQ 'UX' sgauche ;
  bldroite = BLOQ 'UX' sdroite ;
  blrig    = BLOQ 'UY' p1 ;
  bl       = blgauche ET bldroite ET blrig ;
 
** Chargements
  didroite = DEPI bldroite umax ;
  ec1      = EVOL 'MANU' (PROG 0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.)
                         (PROG 0. 0.2 0.  0.4 0.  0.6 0.  0.8 0.  1.  0.) ;
  cha      = CHAR 'DIMP' didroite ec1 ;
 
** Tracé des blocages et du chargement
  re0 = @REPERE (-0.2 -0.2) (PROG 0.5 0.5) ;
  mor = MODE re0 'MECANIQUE' 'BARR' ;
  mar = MATE mor 'YOUN' 1. 'NU' 1. 'SECT' 1. ;
  rir = RIGI mor mar ;
  vf = VECT didroite 'FLX' ' ' 'ROUG' ;
  SI ITRAC1 ;
    TRAC ((RIGI mo ma) ET blgauche ET blrig ET rir) vf 'NOLE' 'TITR' '[2D Plan cont.] Blocages et chargement' ;
  FINSI ;
 
** Resolution
  t                               = TABL ;
  t . 'MODELE'                    = mo ;
  t . 'CARACTERISTIQUES'          = ma ;
  t . 'BLOCAGES_MECANIQUES'       = bl ;
  t . 'CHARGEMENT'                = cha ;
  t . 'TEMPS_CALCULES'            = PROG 0. 'PAS' 0.01 0.1
                                            'PAS' 0.01 0.3
 
*                                           'PAS' 0.01 0.5
                                            'PAS' 0.01 0.46
                                            'PAS' 0.002 0.466
                                            'PAS' 0.00066 0.468
                                            'PAS' 0.002 0.47
 
*                                           'PAS' 0.01 0.7
                                            'PAS' 0.01 0.65
                                            'PAS' 0.005 0.75
 
                                            'PAS' 0.01 0.9
                                            'PAS' 0.01 1.0 ;
  t . 'MOVA'                      = MOT 'D' ;
  t . 'MES_SAUVEGARDES'           = TABL ;
  t . 'MES_SAUVEGARDES' . 'DEFTO' = VRAI ;
  PASAPAS t ;
 
** Post traitement
* courbes : endommagement moyen vs temps
*           contrainte moyenne vs deformation moyenne
*           force de reaction vs deplacement impose
  tt = t . 'TEMPS' ;
  tc = t . 'CONTRAINTES' ;
  te = t . 'DEFORMATIONS' ;
  tv = t . 'VARIABLES_INTERNES' ;
  tu = t . 'DEPLACEMENTS' ;
  tr = t . 'REACTIONS' ;
  def0 = DEFO (cadr1 ET rep1 ET mail) (tu . 0) 0. ;
  ltps = PROG ;
  lend = PROG ;
  lsig = PROG ;
  leps = PROG ;
  ldep = PROG ;
  lrea = PROG ;
  REPE b1 (DIME tt) ;
    tps1 = tt . (&b1 - 1) ;
    ltps = ltps ET tps1 ;
    end1 = (INTG mo (tv . (&b1 - 1)) 'D')    / mes1 ;
    sig1 = (INTG mo (tc . (&b1 - 1)) 'SMXX') / mes1 ;
    eps1 = (INTG mo (te . (&b1 - 1)) 'EPXX') / mes1 ;
    u1   = EXTR (tu . (&b1 - 1)) 'UX' p3 ;
    rea1 = 0. ;
    SI (NEG &b1 1) ;
      rea1 = @TOTAL (tr . (&b1 - 1)) sgauche 'FX' ;
    FINSI ;
    lend = lend ET end1 ;
    lsig = lsig ET sig1 ;
    leps = leps ET eps1 ;
    ldep = ldep ET u1 ;
    lrea = lrea ET rea1 ;
  FIN b1 ;
  ev0 = EVOL 'ROUG' 'MANU' 'Temps' ltps 'D'    lend ;
  ev1 = EVOL 'ROUG' 'MANU' 'EPXX'  leps 'SMXX' lsig ;
  ev2 = EVOL 'ROUG' 'MANU' 'UX'    ldep 'FX'   lrea ;
*nh145313
  ltps1 =  ENLE ltps 1 ;
  nb1 = DIME ltps1 ;
  un1 = PROG nb1 * 1.D0 ;
  lepd0 = PROG nb1 * epd0 ;
  lat1 = PROG nb1 * at ;
  lep1 = IPOL ltps1 ec1 ;
  lepe = lep1 * umax / long ;
*mess 'lepe avant seuillage 0'; list lepe;
  seuil0 = lepe MASQ 'EGAL' 0 ;
*mess 'test seuil0';list seuil0;
  lepe0 = seuil0 * 1.E-20 ;
*mess 'lepe0';list lepe0;
  lepe = lepe + lepe0 ;
*mess 'lepe apres seuillage 0'; list lepe;
  dt = un1 - ((epd0*(un1-lat1)/lepe) + (at*(EXP(bt*(lepd0-lepe))))) ;
  seuil1 = dt MASQ 'EGSUP' 0. ;
  dt = seuil1 * dt ; 
* pas de décroissance du dommage
  dtp = 0.;
  dt2 = PROG ;
  REPE bdt (DIME dt) ;
    dti = EXTR dt &bdt ;
    ddc = dti - dtp ;
    SI (ddc < 0.) ;
      dti = dtp ;
    FINS ;
    dt2 = dt2 ET (PROG dti) ;
    dtp = dti ;
  FIN bdt ;
  dt = dt2 ;
  lsig1 = (un1 - dt) * you * lepe ;
* contrainte nulle = 1.E-20
  seuilsc = lsig MASQ 'EGINF' 1.E-6 ;
  sc0 = seuilsc * 1.E-20 ;
  lsig = (lsig * (lsig MASQ 'EGSUP' 1.E-6)) + sc0 ;
  seuilsr = lsig1 MASQ 'EGINF' 1.E-6 ;
  sr0 = seuilsr * 1.E-20 ;
  lsig1 = (lsig1 * (lsig1 MASQ 'EGSUP' 1.E-6)) + sr0 ;
*mess 'Sig.calcul';list (ENLE lsig 1);
*mess 'Sig.Ref.';list lsig1;
  an0 = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'D'    dt    ;
  an1 = EVOL 'TURQ' 'MANU' 'EPXX'  lepe  'SMXX' lsig1 ;
  SI ITRAC1 ;
    DESS (ev0 ET an0) 'TITR' '[2D Plan cont.] Endommagement vs Temps' 'LEGE' 'SE' tdess ;
    DESS (ev1 ET an1) 'TITR' '[2D Plan cont.] Contrainte vs Deformation' 'LEGE' 'NE' tdess ;
    DESS ev2 'YBOR' -1.8E5 0. 'TITR' '[2D Plan cont.] Reaction vs Deplacement' 'LEGE' 'SE' tdess ;
  FINSI ;
* evolutions temporelles a tous les points de Gauss
  ng = 4 ;
  tleg = TABL ;
  tleg . 1 = MOT 'MARQ CARR NOLI' ;
  tleg . 2 = MOT 'MARQ LOSA NOLI' ;
  tleg . 3 = MOT 'MARQ ROND NOLI' ;
  tleg . 4 = MOT 'MARQ ETOI NOLI' ;
  tleg . 'TITRE' = TABL ;
  REPE b ng ;
    tleg . 'TITRE' . &b = CHAI 'Calc_Pt_Gauss' ' ' &b ;
  FIN b ;
* endommagement
  evend = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evend = evend ET (EVOL 'ROSE' 'TEMP' t 'VARIABLES_INTERNES' 'D' 1 1 &b) ;
  FIN b ;
*nh145313
  tleg . (ng+1) = MOT 'MARQ CROI' ;
  tleg . 'TITRE' . (ng+1) = MOT 'Reference' ;
  SI ITRAC1 ;
    DESS (evend ET an0) 'LEGE' 'SE' tleg 'TITR' '[2D Plan cont.] D vs Temps' ;
  FINSI ;
* contrainte
  evsmxx = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evsmxx = evsmxx ET (EVOL 'BLEU' 'TEMP' t 'CONTRAINTES' 'SMXX' 1 1 &b) ;
  FIN b ;
*nh145313
  ansmxx = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'SMXX' lsig1 ;
  SI ITRAC1 ;
    DESS (evsmxx ET ansmxx) 'LEGE' 'NE' tleg 'TITR' '[2D Plan cont.] SMXX vs Temps' ;
  FINSI ;
* deformation
  evepxx = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evepxx = evepxx ET (EVOL 'VERT' 'TEMP' t 'DEFORMATIONS' 'EPXX' 1 1 &b) ;
  FIN b ;
  SI ITRAC1 ;
    DESS evepxx 'LEGE' 'NO' tleg 'TITR' '[2D Plan cont.] EPXX vs Temps' ;
  FINSI ;
* deformee
  REPE b2 (DIME tt) ;
    tps1 = tt . (&b2 - 1) ;
*   trace de la deformee (pour controle visuel)
    def1 = DEFO mail (tu . (&b2 - 1)) 333. 'ROUG' ;
    tit1 = CHAI '[2D Plan cont.] Deformee au temps' ' ' tps1 ;
    SI ITRAC1 ;
*      TRAC (def0 ET def1) 'TITR' tit1 ;
    FINSI ;
  FIN b2 ;
  SI ITRAC1 ;
    TRAC (def0 ET def1) 'TITR' tit1 ;
  FINSI ;
* comparaison Calcul / Reference sur la contrainte moyenne
  ler1 = ABS ((ENLE lsig 1) - lsig1) / lsig1 ;
  mer1 = MAXI ler1 ;
  eer1 = EVOL 'ORAN' 'MANU' 'EPXX'  lepe  'Ec.rel.SMXX' ler1 ;
  tit1 = CHAI '[2D Plan cont.] Ecart relatif Calcul/Reference Contrainte vs Deformation, max' ' ' mer1 ;
  SI ITRAC1 ;
    DESS eer1 'TITR' tit1 ;
  FINSI ;
  OPTI ECHO 0 ;
  SAUT 1 LIGN ;
  MESS ' ******* C A S  #3  :  2 D   C O N T R A I N T E S   P L A N E S ********';
  SAUT 1 LIGN ;
  SI (mer1 < 1.E-10) ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' < 1.E-10 => OK' ;
    SAUT 1 LIGN ;
  SINO ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' > 1.E-10 => ECHEC' ;
    IERR = VRAI ;
    SAUT 1 LIGN ;
  FINS ;
  SAUT 1 LIGN ;
  OPTI ECHO 1 ;
FINSI ;
 
 
 
 
 
************ C A S  #4  :  2 D   A X I S Y M E T R I Q U E *************
SI (EXIS lcas 4) ;
 
** Options generales
  OPTI 'DIME' 2 'MODE' 'AXIS' 'ELEM' 'QUA4' ;
 
** Maillage
  p1      = 0. 0. ;
  p2      = long 0. ;
  sbas    = DROI 1 p1 p2 ;
  mail    = sbas TRAN 1 (0. long) ;
  shaut   = mail COTE 3 ;
  p3      = shaut POIN 'PROC' (0. long) ;
*nh145313
cadr1 = COUL ((0. 1.2) ET (1. 1.2)) 'BLAN' ;
rep1 = @REPERE (-0.2 -0.2) (PROG 0.5 0.5) 'TURQ' ;
SI ITRAC1 ;
  TRAC 'QUAL' (cadr1 ET rep1 ET mail ET (sbas COUL 'VERT') ET (shaut COUL 'ROUG')) 'TITR' '[2D Axi.] Maillage' ;
FINSI ;
 
** Modele et caracteristiques materiau
  mo = MODE mail 'MECANIQUE' 'ELASTIQUE' 'ENDOMMAGEMENT' 'MAZARS' ;
  ma = MATE mo 'YOUN' you 'NU' nu 'KTR0' epd0 'ACOM' ac 'BCOM' bc 'ATRA' at 'BTRA' bt 'BETA' beta ;
  zer = ZERO mo 'NOEUD' ;
  un  = zer + 1. ;
  mes1 = INTG mo un 'SCAL' ;
 
** Blocages
  blbas    = BLOQ 'UZ' sbas ;
  blhaut   = BLOQ 'UZ' shaut ;
  bl       = blbas ET blhaut ;
 
** Chargements
  dihaut   = DEPI blhaut umax ;
  ec1      = EVOL 'MANU' (PROG 0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.)
                         (PROG 0. 0.2 0.  0.4 0.  0.6 0.  0.8 0.  1.  0.) ;
  cha      = CHAR 'DIMP' dihaut ec1 ;
 
** Tracé des blocages et du chargement
  re0 = @REPERE (-0.2 -0.2) (PROG 0.5 0.5) ;
  mor = MODE re0 'MECANIQUE' 'BARR' ;
  mar = MATE mor 'YOUN' 1. 'NU' 1. 'SECT' 1. ;
  rir = RIGI mor mar ;
  vf = VECT dihaut ' ' 'FLX' 'ROUG' ;
  SI ITRAC1 ;
    TRAC ((RIGI mo ma) ET blbas ET rir) vf 'NOLE' 'TITR' '[2D Axi.] Blocages et chargement' ;
  FINSI ;
 
** Resolution
  t                               = TABL ;
  t . 'MODELE'                    = mo ;
  t . 'CARACTERISTIQUES'          = ma ;
  t . 'BLOCAGES_MECANIQUES'       = bl ;
  t . 'CHARGEMENT'                = cha ;
  t . 'TEMPS_CALCULES'            = PROG 0. 'PAS' 0.01 0.1
                                            'PAS' 0.01 0.3
 
*                                           'PAS' 0.01 0.5
                                            'PAS' 0.01 0.46
                                            'PAS' 0.002 0.466
                                            'PAS' 0.00066 0.468
                                            'PAS' 0.002 0.47
 
*                                           'PAS' 0.01 0.7
                                            'PAS' 0.01 0.65
                                            'PAS' 0.005 0.75
 
                                            'PAS' 0.01 0.9
                                            'PAS' 0.01 1.0 ;
  t . 'MOVA'                      = MOT 'D' ;
  t . 'MES_SAUVEGARDES'           = TABL ;
  t . 'MES_SAUVEGARDES' . 'DEFTO' = VRAI ;
  PASAPAS t ;
 
** Post traitement
* courbes : endommagement moyen vs temps
*           contrainte moyenne vs deformation moyenne
*           force de reaction vs deplacement impose
  tt = t . 'TEMPS' ;
  tc = t . 'CONTRAINTES' ;
  te = t . 'DEFORMATIONS' ;
  tv = t . 'VARIABLES_INTERNES' ;
  tu = t . 'DEPLACEMENTS' ;
  tr = t . 'REACTIONS' ;
  def0 = DEFO (cadr1 ET rep1 ET mail) (tu . 0) 0. ;
  ltps = PROG ;
  lend = PROG ;
  lsig = PROG ;
  leps = PROG ;
  ldep = PROG ;
  lrea = PROG ;
  REPE b1 (DIME tt) ;
    tps1 = tt . (&b1 - 1) ;
    ltps = ltps ET tps1 ;
    end1 = (INTG mo (tv . (&b1 - 1)) 'D')    / mes1 ;
    sig1 = (INTG mo (tc . (&b1 - 1)) 'SMZZ') / mes1 ;
    eps1 = (INTG mo (te . (&b1 - 1)) 'EPZZ') / mes1 ;
    u1   = EXTR (tu . (&b1 - 1)) 'UZ' p3 ;
    rea1 = 0. ;
    SI (NEG &b1 1) ;
      rea1 = @TOTAL (tr . (&b1 - 1)) sbas 'FZ' ;
    FINSI ;
    lend = lend ET end1 ;
    lsig = lsig ET sig1 ;
    leps = leps ET eps1 ;
    ldep = ldep ET u1 ;
    lrea = lrea ET rea1 ;
  FIN b1 ;
  ev0 = EVOL 'ROUG' 'MANU' 'Temps' ltps 'D'    lend ;
  ev1 = EVOL 'ROUG' 'MANU' 'EPZZ'  leps 'SMZZ' lsig ;
  ev2 = EVOL 'ROUG' 'MANU' 'UZ'    ldep 'FZ'   lrea ;
*nh145313
  ltps1 =  ENLE ltps 1 ;
  nb1 = DIME ltps1 ;
  un1 = PROG nb1 * 1.D0 ;
  lepd0 = PROG nb1 * epd0 ;
  lat1 = PROG nb1 * at ;
  lep1 = IPOL ltps1 ec1 ;
  lepe = lep1 * umax / long ;
*mess 'lepe avant seuillage 0'; list lepe;
  seuil0 = lepe MASQ 'EGAL' 0 ;
*mess 'test seuil0';list seuil0;
  lepe0 = seuil0 * 1.E-20 ;
*mess 'lepe0';list lepe0;
  lepe = lepe + lepe0 ;
*mess 'lepe apres seuillage 0'; list lepe;
  dt = un1 - ((epd0*(un1-lat1)/lepe) + (at*(EXP(bt*(lepd0-lepe))))) ;
  seuil1 = dt MASQ 'EGSUP' 0. ;
  dt = seuil1 * dt ; 
* pas de décroissance du dommage
  dtp = 0.;
  dt2 = PROG ;
  REPE bdt (DIME dt) ;
    dti = EXTR dt &bdt ;
    ddc = dti - dtp ;
    SI (ddc < 0.) ;
      dti = dtp ;
    FINS ;
    dt2 = dt2 ET (PROG dti) ;
    dtp = dti ;
  FIN bdt ;
  dt = dt2 ;
  lsig1 = (un1 - dt) * you * lepe ;
* contrainte nulle = 1.E-20
  seuilsc = lsig MASQ 'EGINF' 1.E-6 ;
  sc0 = seuilsc * 1.E-20 ;
  lsig = (lsig * (lsig MASQ 'EGSUP' 1.E-6)) + sc0 ;
  seuilsr = lsig1 MASQ 'EGINF' 1.E-6 ;
  sr0 = seuilsr * 1.E-20 ;
  lsig1 = (lsig1 * (lsig1 MASQ 'EGSUP' 1.E-6)) + sr0 ;
*mess 'Sig.calcul';list (ENLE lsig 1);
*mess 'Sig.Ref.';list lsig1;
  an0 = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'D'    dt    ;
  an1 = EVOL 'TURQ' 'MANU' 'EPZZ'  lepe  'SMZZ' lsig1 ;
  SI ITRAC1 ;
    DESS (ev0 ET an0) 'TITR' '[2D Axi.] Endommagement vs Temps' 'LEGE' 'SE' tdess ;
    DESS (ev1 ET an1) 'TITR' '[2D Axi.] Contrainte vs Deformation' 'LEGE' 'NE' tdess ;
    DESS ev2 'YBOR' -6.E6 0. 'TITR' '[2D Axi.] Reaction vs Deplacement' 'LEGE' 'SE' tdess ;
  FINSI ;
* evolutions temporelles a tous les points de Gauss
  ng = 4 ;
  tleg = TABL ;
  tleg . 1 = MOT 'MARQ CARR NOLI' ;
  tleg . 2 = MOT 'MARQ LOSA NOLI' ;
  tleg . 3 = MOT 'MARQ ROND NOLI' ;
  tleg . 4 = MOT 'MARQ ETOI NOLI' ;
  tleg . 'TITRE' = TABL ;
  REPE b ng ;
    tleg . 'TITRE' . &b = CHAI 'Calc_Pt_Gauss' ' ' &b ;
  FIN b ;
* endommagement
  evend = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evend = evend ET (EVOL 'ROSE' 'TEMP' t 'VARIABLES_INTERNES' 'D' 1 1 &b) ;
  FIN b ;
*nh145313
  tleg . (ng+1) = MOT 'MARQ CROI' ;
  tleg . 'TITRE' . (ng+1) = MOT 'Reference' ;
  SI ITRAC1 ;
    DESS (evend ET an0) 'LEGE' 'SE' tleg 'TITR' '[2D Axi.] D vs Temps' ;
  FINSI ;
* contrainte
  evsmzz = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evsmzz = evsmzz ET (EVOL 'BLEU' 'TEMP' t 'CONTRAINTES' 'SMZZ' 1 1 &b) ;
  FIN b ;
*nh145313
  ansmzz = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'SMZZ' lsig1 ;
  SI ITRAC1 ;
    DESS (evsmzz ET ansmzz) 'LEGE' 'NE' tleg 'TITR' '[2D Axi.] SMZZ vs Temps' ;
  FINSI ;
* deformation
  evepzz = VIDE 'EVOLUTIO' ;
  REPE b ng ;
    evepzz = evepzz ET (EVOL 'VERT' 'TEMP' t 'DEFORMATIONS' 'EPZZ' 1 1 &b) ;
  FIN b ;
  SI ITRAC1 ;
    DESS evepzz 'LEGE' 'NO' tleg 'TITR' '[2D Axi.] EPZZ vs Temps' ;
  FINSI ;
* deformee
  REPE b2 (DIME tt) ;
    tps1 = tt . (&b2 - 1) ;
*   trace de la deformee (pour controle visuel)
    def1 = DEFO mail (tu . (&b2 - 1)) 333. 'ROUG' ;
    tit1 = CHAI '[2D Axi.] Deformee au temps' ' ' tps1 ;
    SI ITRAC1 ;
*      TRAC (def0 ET def1) 'TITR' tit1 ;
    FINSI ;
  FIN b2 ;
  SI ITRAC1 ;
    TRAC (def0 ET def1) 'TITR' tit1 ;
  FINSI ;
* comparaison Calcul / Reference sur la contrainte moyenne
  ler1 = ABS ((ENLE lsig 1) - lsig1) / lsig1 ;
  mer1 = MAXI ler1 ;
  eer1 = EVOL 'ORAN' 'MANU' 'EPZZ'  lepe  'Ec.rel.SMZZ' ler1 ;
  tit1 = CHAI '[2D Axi.] Ecart relatif Calcul/Reference Contrainte vs Deformation, max' ' ' mer1 ;
  SI ITRAC1 ;
    DESS eer1 'TITR' tit1 ;
  FINSI ;
  OPTI ECHO 0 ;
  SAUT 1 LIGN ;
  MESS ' ************ C A S  #4  :  2 D   A X I S Y M E T R I Q U E *************';
  SAUT 1 LIGN ;
  SI (mer1 < 1.E-10) ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' < 1.E-10 => OK' ;
    SAUT 1 LIGN ;
  SINO ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' > 1.E-10 => ECHEC' ;
    IERR = VRAI ;
    SAUT 1 LIGN ;
  FINS ;
  SAUT 1 LIGN ;
  OPTI ECHO 1 ;
FINSI ;
 
 
 
 
 
*********** C A S  #5  :  2 D   P O U T R E   A   F I B R E ************
SI (EXIS lcas 5) ;
 
** Options generales
  OPTI 'DIME' 2 'MODE' 'PLAN' 'ELEM' 'QUA4' ;
 
** Maillage
  p1s     = (-0.5 * epai) (-0.5 * epai) ;
  p2s     = (-0.5 * epai) ( 0.5 * epai) ;
  l12s    = DROI 1 p1s p2s ;
  mails   = l12s TRAN 1 (epai 0.) ;
  p1      = 0. 0. ;
  p2      = long 0. ;
  mail    = DROI 1 p1 p2 ;
  mes1    = MESU mails ;
*nh145313
cadr1 = COUL ((-0.1 0.) ET (1.2 0.)) 'BLAN' ;
rep1 = @REPERE (-0.2 -0.2) (PROG 0.2 0.2) 'TURQ' ;
SI ITRAC1 ;
  TRAC 'QUAL' (cadr1 ET rep1 ET (mails COUL 'VERT') ET (mail COUL 'ROUG')) 'TITR' '[2D PaF] Maillages de la Poutre et de sa Section' ;
FINSI ;
 
** Modele et caracteristiques materiau
  mos = MODE mails 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'MAZARS' 'QUAS' ;
  mas = MATE mos 'YOUN' you 'NU' nu 'KTR0' epd0 'ACOM' ac 'BCOM' bc 'ATRA' at 'BTRA' bt 'BETA' beta
                 'ALPY' 1. ;
  mo = MODE mail 'MECANIQUE' 'ELASTIQUE' 'SECTION' 'PLASTIQUE' 'SECTION' 'TIMO' ;
  ma = MATE mo 'MODS' mos 'MATS' mas ;
 
** Blocages
  blgauche = BLOQ 'DEPL' 'ROTA' p1 ;
  bldroite = BLOQ 'UX' p2 ;
  bl       = blgauche ET bldroite ;
*nh145313
  bldroit2 = BLOQ 'ROTA' p2 ;
  bl       = bl ET bldroit2 ;
 
** Chargements
  didroite = DEPI bldroite umax ;
  ec1      = EVOL 'MANU' (PROG 0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.)
                         (PROG 0. 0.2 0.  0.4 0.  0.6 0.  0.8 0.  1.  0.) ;
  cha      = CHAR 'DIMP' didroite ec1 ;
 
** Tracé des blocages et du chargement
  re0 = @REPERE (-0.2 -0.2) (PROG 0.2 0.2) ;
  mor = MODE re0 'MECANIQUE' 'BARR' ;
  mar = MATE mor 'YOUN' 1. 'NU' 1. 'SECT' 1. ;
  rir = RIGI mor mar ;
  vf = VECT didroite 'FLX' ' ' 'ROUG' 200. ;
  SI ITRAC1 ;
    TRAC ((RIGI mo ma) ET blgauche ET rir) vf 'NOLE' 'TITR' '[2D PaF] Blocages et chargement' ;
  FINSI ;
 
** Resolution
  t                               = TABL ;
  t . 'MODELE'                    = mo ;
  t . 'CARACTERISTIQUES'          = ma ;
  t . 'BLOCAGES_MECANIQUES'       = bl ;
  t . 'CHARGEMENT'                = cha ;
  t . 'TEMPS_CALCULES'            = PROG 0. 'PAS' 0.01 0.1
                                            'PAS' 0.01 0.3
 
*                                           'PAS' 0.01 0.5
                                            'PAS' 0.01 0.46
                                            'PAS' 0.002 0.466
                                            'PAS' 0.00066 0.468
                                            'PAS' 0.002 0.47
 
*                                           'PAS' 0.01 0.7
                                            'PAS' 0.01 0.65
                                            'PAS' 0.005 0.75
 
                                            'PAS' 0.01 0.9
                                            'PAS' 0.01 1.0 ;
  t . 'MES_SAUVEGARDES'           = TABL ;
  t . 'MES_SAUVEGARDES' . 'DEFTO' = VRAI ;
  PASAPAS t ;
 
** Post traitement
* courbes : endommagement moyen vs temps
*           contrainte moyenne vs deformation moyenne
*           force de reaction vs deplacement impose
* evolutions temporelles a tous les points de Gauss
  ng = 4 ;
  tleg = TABL ;
  tleg . 1 = MOT 'MARQ CARR NOLI' ;
  tleg . 2 = MOT 'MARQ LOSA NOLI' ;
  tleg . 3 = MOT 'MARQ ROND NOLI' ;
  tleg . 4 = MOT 'MARQ ETOI NOLI' ;
  tleg . 'TITRE' = TABL ;
  REPE b ng ;
    tleg . 'TITRE' . &b = CHAI 'Calc_Pt_Gauss' ' ' &b ;
  FIN b ;
*nh145313
  tleg . (ng+1) = MOT 'MARQ CROI' ;
  tleg . 'TITRE' . (ng+1) = MOT 'Reference' ;
  tt = t . 'TEMPS' ;
  tv = t . 'VARIABLES_INTERNES' ;
  tu = t . 'DEPLACEMENTS' ;
  tr = t . 'REACTIONS' ;
  def0 = DEFO (cadr1 ET rep1 ET mail) (tu . 0) 0. ;
  ltps = PROG ;
  lend1 = PROG ;
  lend2 = PROG ;
  lend3 = PROG ;
  lend4 = PROG ;
  lsig1 = PROG ;
  lsig2 = PROG ;
  lsig3 = PROG ;
  lsig4 = PROG ;
  leps1 = PROG ;
  leps2 = PROG ;
  leps3 = PROG ;
  leps4 = PROG ;
  ldep = PROG ;
  lrea = PROG ;
  REPE b1 (DIME tt) ;
    tps1 = tt . (&b1 - 1) ;
    ltps = ltps ET tps1 ;
    SI (EGA &b1 1) ;
      lend1 = lend1 ET 0. ;
      lend2 = lend2 ET 0. ;
      lend3 = lend3 ET 0. ;
      lend4 = lend4 ET 0. ;
      lsig1 = lsig1 ET 0. ;
      lsig2 = lsig2 ET 0. ;
      lsig3 = lsig3 ET 0. ;
      lsig4 = lsig4 ET 0. ;
      leps1 = leps1 ET 0. ;
      leps2 = leps2 ET 0. ;
      leps3 = leps3 ET 0. ;
      leps4 = leps4 ET 0. ;
    SINON ;
      lend1 = lend1 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 1) ;
      lend2 = lend2 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 2) ;
      lend3 = lend3 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 3) ;
      lend4 = lend4 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VHNI' 1 1 4) ;
      lsig1 = lsig1 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 1) ;
      lsig2 = lsig2 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 2) ;
      lsig3 = lsig3 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 3) ;
      lsig4 = lsig4 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VONS' 1 1 1) 'SMXX' 1 1 4) ;
      leps1 = leps1 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 1) ;
      leps2 = leps2 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 2) ;
      leps3 = leps3 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 3) ;
      leps4 = leps4 ET (EXTR (EXTR (tv . (&b1 - 1)) 'VAIS' 1 1 1) 'VISX' 1 1 4) ;
    FINSI ;
    u1   = EXTR (tu . (&b1 - 1)) 'UX' p2 ;
    rea1 = EXTR (tr . (&b1 - 1)) 'FX' p1 ;
    ldep = ldep ET u1 ;
    lrea = lrea ET rea1 ;
  FIN b1 ;
  ev0 = (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend1) ET (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend2) ET
        (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend3) ET (EVOL 'ROUG' 'MANU' 'Temps' ltps 'D' lend4) ;
  ev1 = (EVOL 'ROUG' 'MANU' 'EPXX'  leps1 'SMXX' lsig1) ET (EVOL 'ROUG' 'MANU' 'EPXX' leps2 'SMXX' lsig2) ET
        (EVOL 'ROUG' 'MANU' 'EPXX'  leps3 'SMXX' lsig3) ET (EVOL 'ROUG' 'MANU' 'EPXX' leps4 'SMXX' lsig4) ;
  ev2 = EVOL 'ROUG' 'MANU' 'UX' ldep 'FX' lrea ;
*nh145313
  ltps1 =  ENLE ltps 1 ;
  nb1 = DIME ltps1 ;
  un1 = PROG nb1 * 1.D0 ;
  lepd0 = PROG nb1 * epd0 ;
  lat1 = PROG nb1 * at ;
  lep1 = IPOL ltps1 ec1 ;
  lepe = lep1 * umax / long ;
*mess 'lepe avant seuillage 0'; list lepe;
  seuil0 = lepe MASQ 'EGAL' 0. ;
*mess 'test seuil0';list seuil0;
  lepe0 = seuil0 * 1.E-20 ;
*mess 'lepe0';list lepe0;
  lepe = lepe + lepe0 ;
*mess 'lepe apres seuillage 0'; list lepe;
  dt = un1 - ((epd0*(un1-lat1)/lepe) + (at*(EXP(bt*(lepd0-lepe))))) ;
  seuil1 = dt MASQ 'EGSUP' 0 ;
  dt = seuil1 * dt ; 
* pas de décroissance du dommage
  dtp = 0.;
  dt2 = PROG ;
  REPE bdt (DIME dt) ;
    dti = EXTR dt &bdt ;
    ddc = dti - dtp ;
    SI (ddc < 0.) ;
      dti = dtp ;
    FINS ;
    dt2 = dt2 ET (PROG dti) ;
    dtp = dti ;
  FIN bdt ;
  dt = dt2 ;
  lsan1 = (un1 - dt) * you * lepe ;
* contrainte nulle = 1.E-20
  seuilsr = lsan1 MASQ 'EGINF' 1.E-6 ;
  sr0 = seuilsr * 1.E-20 ;
  lsan1 = (lsan1 * (lsan1 MASQ 'EGSUP' 1.E-6)) + sr0 ;
*mess 'Sig.Ref.';list lsan1;
  an0 = EVOL 'TURQ' 'MANU' 'Temps' ltps1 'D'    dt    ;
  an1 = EVOL 'TURQ' 'MANU' 'EPXX'  lepe  'SMXX' lsan1 ;
  SI ITRAC1 ;
    DESS (ev0 ET an0) 'LEGE' 'SE' tleg 'TITR' '[2D PaF] Endommagement vs Temps' ;
    DESS (ev1 ET an1) 'LEGE' 'NE' tleg 'TITR' '[2D PaF] Contrainte vs Deformation' ;
    DESS ev2 'TITR' '[2D PaF] Reaction vs Deplacement' 'LEGE' 'SE' tdess ;
  FINSI ;
* deformee
  REPE b2 (DIME tt) ;
    tps1 = tt . (&b2 - 1) ;
*   trace de la deformee (pour controle visuel)
    def1 = DEFO mail (tu . (&b2 - 1)) 333. 'ROUG' ;
    tit1 = CHAI '[2D PaF] Deformee au temps' ' ' tps1 ;
    SI ITRAC1 ;
*      TRAC (def0 ET def1) 'TITR' tit1 ;
    FINSI ;
  FIN b2 ;
  SI ITRAC1 ;
    TRAC (def0 ET def1) 'TITR' tit1 ;
  FINSI ;
* comparaison Calcul / Reference sur la contrainte moyenne
  lsig = (lsig1 + lsig2 + lsig3 + lsig4) / 4. ;
* contrainte nulle = 1.E-20
  seuilsc = lsig MASQ 'EGINF' 1.E-6 ;
  sc0 = seuilsc * 1.E-20 ;
  lsig = (lsig * (lsig MASQ 'EGSUP' 1.E-6)) + sc0 ;
*mess 'Sig.calcul';list (ENLE lsig 1);
  ler1 = ABS ((ENLE lsig 1) - lsan1) / lsan1 ;
  mer1 = MAXI ler1 ;
  eer1 = EVOL 'ORAN' 'MANU' 'EPXX'  lepe  'Ec.rel.SMXX' ler1 ;
  tit1 = CHAI '[2D PaF] Ecart relatif Calcul/Reference Contrainte vs Deformation, max' ' ' mer1 ;
  SI ITRAC1 ;
    DESS eer1 'TITR' tit1 ;
  FINSI ;
  OPTI ECHO 0 ;
  SAUT 1 LIGN ;
  MESS ' *********** C A S  #5  :  2 D   P O U T R E   A   F I B R E ************';
  SAUT 1 LIGN ;
  SI (mer1 < 1.E-10) ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' < 1.E-10 => OK' ;
    SAUT 1 LIGN ;
  SINO ;
    MESS ' Ecart relatif max entre Calcul et Reference' ' ' mer1 ' > 1.E-10 => ECHEC' ;
    IERR = VRAI ;
    SAUT 1 LIGN ;
  FINS ;
  SAUT 1 LIGN ;
  OPTI ECHO 1 ;
FINSI ;
 
** Sortie sur Erreur 1152 en cas de non conformite
SI IERR ;
  ERRE 1152 ;
FINSI ;
 
 
 
 
 
FIN ;