Accueil | Cast3M

************************************************************************
************************************************************************
************************************************************************
*______________________________________________________________________*
*                                                                      *
*      Cas test de l'implantation Nonlocal Stress Based (NLSB)         *
*______________________________________________________________________*
*----------------------------------------------------------------------*
* Description du cas test : Chargement uniaxial sur un cube            *
* -----------------------                                              *
*                                                                      *
* L'objectif de l'essai est de verifier le bon fonctionnement de la    *
* regularisation nonlocale base sur l'etat de contrainte dans un calcul*
* nonlineaire (PASAPAS). Le modeles d'endommagement considere   pour   *
* le test de validation est le modele MAZARS.                          *
*                                                                      *
* On remarquera, que la regularisation n'est pas dependante du modele  *
* et qu'elle peut-etre applique a n'importe quel modele a partir du    *
* moment ou la variable a regulariser est fournie.                     * 
*                                                                      *
* C. Giry, F. Dufour, J. Mazars                                        *
* -----------------------------                                        *
* (fichier base sur ricrag_3D.dgibi)                                   *
************************************************************************
graph = 'N';
COMPLET = FAUX ;
 
* Choix du modele d'endommagement teste
* -------------------------------------
* nmod = 1; modele MAZARS
 
nmod  = 1;
*
*-------------- Options de calcul ---------------------
*
OPTION DIME 3 ELEM CUB8;
 
nbel  = 5;
tele  = 1./nbel; 
tdef  = 0.02;
nbelh = 2;
*
*-------------- Definition de la geometrie ------------
*
P1 = 0. 0. 0.;
P2 = 1. 0. 0.;
P3 = 1. 1. 0.;
P4 = 0. 1. 0.;
 
L1 = P1 DROIT nbel P2 ;
L2 = P2 DROIT nbel P3 ;
L3 = P3 DROIT nbel P4 ;
L4 = P4 DROIT nbel P1 ;
 
SURF1 = DALL L1 L2 L3 L4;
 
P11 = tdef tdef (0.5-(0.5*tele));
P21 = (1-tdef) tdef (0.5-(0.5*tele));
P31 = (1-tdef) (1-tdef) (0.5-(0.5*tele));
P41 = tdef (1-tdef) (0.5-(0.5*tele));
 
L11 = P11 DROIT nbel P21 ;
L21 = P21 DROIT nbel P31 ;
L31 = P31 DROIT nbel P41 ;
L41 = P41 DROIT nbel P11 ;
 
SURF11 = DALL L11 L21 L31 L41;
 
P12 = tdef tdef (0.5+(0.5*tele));
P22 = (1-tdef) tdef (0.5+(0.5*tele));
P32 = (1-tdef) (1-tdef) (0.5+(0.5*tele));
P42 = tdef (1-tdef) (0.5+(0.5*tele));
 
L12 = P12 DROIT nbel P22 ;
L22 = P22 DROIT nbel P32 ;
L32 = P32 DROIT nbel P42 ;
L42 = P42 DROIT nbel P12 ;
 
SURF12 = DALL L12 L22 L32 L42;
 
P13 = 0. 0. 1.;
P23 = 1. 0. 1.;
P33 = 1. 1. 1.;
P43 = 0. 1. 1.;
 
L13 = P13 DROIT nbel P23 ;
L23 = P23 DROIT nbel P33 ;
L33 = P33 DROIT nbel P43 ;
L43 = P43 DROIT nbel P13 ;
 
SURF2 = DALL L13 L23 L33 L43;
 
 
VOLTO = SURF1 VOLU nbelh SURF11 ;
VOLT1 = SURF11 VOLU 1 SURF12 ;
VOLT2 = SURF12 VOLU nbelh SURF2 ;
 
VOLTOT = VOLTO ET VOLT1 ET VOLT2;
*
*------- Definition des conditions aux limites --------
*----------- et des deplacements imposes --------------
*
CL  = BLOQ SURF1 UZ;
CLL = BLOQ P1 'DEPL';
CL1 = BLOQ SURF2 UZ;
CL2 = bloq P2 UY;
CL  = CL et CL2;
 
D1 = DEPI CL1 1;
*
*----- Definition du modele ---------------------------
*
MOD1 = MODE VOLTOT MECANIQUE ELASTIQUE ISOTROPE ENDOMMAGEMENT MAZARS 
            'NON_LOCAL' 'SB'  'V_MOYENNE' ('MOTS' 'EPTI') ;
* Parametres
* ----------
EBETON  = 34100E6; 
NUBETON = 0.16;  
 
* Longueur interne
lc      = 0.3;
 
* Module d'Young
youngn = EBETON;
 
* Coefficient de Poisson
nun    = NUBETON;
 
* Résistance en traction
ftn    = 3.6e6;
 
***************** Modele MAZARS **************
 
* Seuil de deformation 1ere fissuration 
 y0    = ftn/ebeton ; 
 
* Parametres At et Bt (comportement post-pic en traction)
 at = 0.95  ; 
 bt = 10000.;
 
* Parametres At et Bt (comportement non lineaire en compression) 
 ac = 1.25 ; 
 bc = 1000.;
 
* Parametre beta (correction pour le cisaillement) 
 bet   = 1.06   ;
 
*
* Material
*    
 
mat1 = MATE mod1 YOUN youngn NU   nun KTR0 y0     
                 ATRA at     BTRA bt  
                 ACOM ac     BCOM bc  BETA bet 
                'LCAR' LC   'SBFT' FTN ;
*
*----- Definition du chargement -------------------
*
LI1  = PROG 0. 1. ;
LI2  = PROG 0. 2.e-3;
'SI' ('NON' COMPLET) ;
  LIS1 = PROG 0. PAS 0.005 0.15 ;
'SINON' ;
  LIS1 = PROG 0. PAS 0.005 1.00 ;
'FINSI' ;
EV = EVOL MANU LI1 LI2 ;
CHA1 = CHAR 'DIMP' D1 EV ;
*
*----------- Calcul par l'operateur PASAPAS ------------
*
CO1    = CONNEC mod1 MAT1 'NORMAL';
CHAMFT = MANU 'CHML' mod1 'FT' ftn 'STRESSES';
 
TAB1                       = TABLE ;
TAB1.'BLOCAGES_MECANIQUES' = CL ET CLL ET CL1;
TAB1.'MODELE'              = MOD1;
TAB1.'MOVA'                = 'MOT' 'D   ';
TAB1.'CHARGEMENT'          = CHA1;
TAB1.'CARACTERISTIQUES'    = MAT1;
TAB1.'TEMPS_CALCULES'      = LIS1;
*TAB1.'CONN'                = CO1;
**TAB1.'NON_LOCAL'           = 'MOT' 'SB' ;
**TAB1.'LVIAM'               = 'MOTS' 'EPTI' ;
**TAB1.'LCAR'                = LC;
**TAB1.'NLOC_SB_MODL'        = MOD1;
**TAB1.'NLOC_SB_FT'          = CHAMFT;
**TAB1.'NLOC_SB_FT'          = ftn ;
 
PASAPAS TAB1 ;
 
dd = dime tab1.deplacements;
 
*
*----------- Courbe effort-deplacement -----------------
*
si (ega graph 'O');
 ev2=@global tab1 CL1 EV fz;
 DESS EV2;
 endf = EXCO (tab1.variables_internes.(dd-1)) 'D';
 trac endf mod1;
finsi;
 
endf        = EXCO (tab1.variables_internes.(dd-1)) 'D';
'SI' ('NON' COMPLET) ;
  dlim        = 0.99 ;
'SINO' ;
  dlim        = 0.999 ;
'FINSI' ;
mailend     = endf ELEM 'SUPERIEUR' dlim;
XX1 YY1 ZZ1 = COOR MAILEND;
Zmin        = MINI ZZ1;
Zmax        = MAXI ZZ1;
 
message 'Taille zone completement endommagee';
err = Zmax-Zmin;
list err;
 
si (err > (1.1*tele));
erreur(5);
finsi;
 
fin;