Accueil | Cast3M

* fichier : ricjoi_2d.dgibi 
************************************************************************
************************************************************************
******************************************************* 
*  Cas test de l'implantation numérique du modele 
*  RICJOI 2D LOCAL 
******************************************************* 
* Développé par :  
*     Benjamin Richard 
* Contact :  
* Benjamin.Richard@lmt.ens-cachan.fr 
******************************************************* 
* Les cas de charges sont entrés :  
* - 1 : Essai cyclique mono cadrant 
* - 2 : Essai cyclique bi cadrant (cisaillement) 
* - 3 : Essai effet unilatteral 
* - 4 : Essai frottement 
* - 5 : Essai variations stiffness 
* - 6 : Essai compression normale avec rouille monotone  
******************************************************* 
* Choix du cas de charge 
ncas  = 1; 
graph= MOT 'N';
 
* 
*-------------- Options de calcul --------------------- 
* 
OPTION DIME 2 ELEM SEG2 MODE PLAN DEFO ; 
 
* 
*-------------- Definition de la geometrie ------------ 
* 
A = 0.00 0.00 ; 
B = 1.00 0.00; 
C = 0.00 0.00 ; 
D = 1.00 0.00 ; 
 
AB =  A DROIT 1 B ; 
CD =  C DROIT 1 D ; 
 
OPTION ELEM RAC2 ; 
 
JOINT = RACC 0.001 CD AB ; 
 
* 
*------- Définition des conditions aux limites -------- 
*----------- et des déplacements imposés -------------- 
* 
CL1 = BLOQ AB UX ; 
CL2 = BLOQ AB UY ; 
 
SI (OU (EGA ncas 3) (EGA ncas 6)); 
CL3 = BLOQ CD UY ; 
SINON; 
CL3 = BLOQ CD UX ; 
FINSI; 
 
CLJOINT =  CL1 ET CL2 ET CL3 ; 
 
* 
*----- Définition du modèle --------------------------- 
* 
MODJOI = MODE JOINT MECANIQUE ELASTIQUE PLASTIQUE INTJOI JOI2 ; 
 
* 
*----- Paramètres matériaux --------------------------- 
*  
 
* Rigidite normale 
knn = (70000E+6); 
 
* Rigidite tangentielle 
ksn = (knn/2.4); 
 
* Coefficient de dilatation 
alpnn = 1.0E-10; 
 
* Fragilite 
adn = 1.1E-5; 
 
* Seuil initial en energie pour l'endommagement 
y0n = 50; 
 
* Coefficient de dissymétrie entre modes 
alpan = 6; 
 
* Coefficient d'écrouissage 1 
gainn = 2.0E9; 
 
* Coefficient d'écrouissage 2 
aainn = 5.0E-7; 
 
* Coefficient de Gurson modifié 1 
q1con = 3.5; 
 
* Coefficient de Gurson modifié 2 
q2con = 0.9; 
 
* Coefficient de Gurson modifié 3 
q3con = 0.1; 
 
* Contrainte de raidissement seuil 
sycon = -1.0E6; 
 
* Coefficient loi de contraction 1 
ncoen = 2.0; 
 
* Coefficient loi de contraction 2  
kcoen = 1.0E+10; 
 
* Degré de corrosion (en perte de section) 
SI (EGA ncas 6); 
tcn = 1.0E-10; 
SINON; 
tcn = 0.0; 
FINSI; 
 
* Indicateur de phase de calcul 
gonfn = 0; 
 
MAJOI = MATE MODJOI KS ksn KN knn  ALPN alpnn 
         AD    adn Y0 y0n ALPA alpan        
         GAIN  gainn AAIN aainn                
         Q1CO  q1con Q2CO q2con Q3CO q3con   
         SYCO sycon NCOE ncoen KCOE kcoen  
         TC tcn GONF gonfn; 
 
* 
*----- Définition des cas de charge ------------------- 
* 
 
* 
* ----------            DEFINITION DU CHARGEMENT       ---------- 
* 
* 
* ----------   DEFINITION DE LA PROGRESSION DU CHARGEMENT ---------- 
* 
 
* Essai cyclique mono cadrant 
SI (EGA ncas 1); 
  LI1  = PROG 0. 1.      2.     3.      ; 
  LI2  = PROG 0. 5.0E-3 0.0E0   10.0E-3 ; 
  LIS1 = PROG 0.0 PAS 0.01 3.0; 
FINSI; 
 
* Essai cyclique bi cadrant (cisaillement) 
SI (EGA ncas 2); 
  LI1  = PROG 0. 1. 2. 3. 4. 5. 6. 7. 8. 9; 
  LI2  = PROG 0. -3.0E-3 5.0E-3 -7.0E-3 9.0E-3 -11.0E-3 13.0E-3 
             -15.0E-3 17E-3 0.0; 
  LIS1 = PROG 0.0 PAS 0.01 9.0; 
FINSI; 
 
* Essai effet unilatteral 
SI (EGA ncas 3); 
  LI1  = PROG 0. 1. 2. 3.; 
  LI2  = PROG 0. 5.0E-3 -1.0E-3 10E-3; 
  LIS1 = PROG 0. PAS 0.01 3.0; 
FINSI; 
 
* Essai frottement 
SI (EGA ncas 4); 
  LI1  = PROG 0. pas 1. 7.; 
  LI2  = PROG 0. 5.0E-3 0.0E-3 6e-3 0.0e-3 8e-3 0.0e-3 10E-3; 
  LIS1 = PROG 0. PAS 0.01 7.0; 
FINSI; 
 
* Essai variations stiffness 
SI (EGA ncas 5); 
  LI1  = PROG 0. pas 1. 1.; 
  LI2  = PROG 0. 10.0E-3 ; 
  LIS1 = PROG 0. PAS 0.01 1.0; 
FINSI; 
 
* Essai compression normale avec rouille monotone 
SI (EGA ncas 6); 
  LI1  = PROG 0. pas 1. 1.; 
  LI2  = PROG 0. -1.0E-4 ; 
  LIS1 = PROG 0. PAS 0.01 1; 
FINSI; 
 
  RIG1 = DEPI CL3 1.0 ; 
  EV1  = EVOL MANU LI1 LI2; 
  CAR1 = CHAR 'DIMP' RIG1 EV1 ;                                                                              
 
* 
*----------- Calcul par l'operateur PASAPAS ------------ 
* 
   TAB                       = TABLE ; 
   TAB.'MOVA'                = 'DINT'; 
   TAB.'BLOCAGES_MECANIQUES' = CLJOINT; 
   TAB.'CARACTERISTIQUES'    = MAJOI; 
   TAB.'MODELE'              = MODJOI; 
   TAB.'TEMPS_CALCULES'      = LIS1; 
   TAB.'CHARGEMENT'          = CAR1; 
 
   PASAPAS TAB ; 
* 
* ---------- POST TRAITEMENT ---------- 
* 
 
* 
*----------- Courbe effort-deplacement ----------------- 
* 
 dep0 = TAB.'DEPLACEMENTS'; 
 tim0 = tab. TEMPS; 
 sig0 = TAB.'CONTRAINTES'; 
 var0 = TAB.'VARIABLES_INTERNES'; 
 ind0 = INDE dep0; 
 time = PROG (tim0 . 0); 
 
 uzp0  = PROG 0.; 
 fzp0  = PROG 0.; 
 
 ep1 = prog 0.; 
 si1 = prog 0.; 
 nn = (('DIME' ind0) - 1); 
 
 j = 1; 
 
 REPE lab1 (('DIME' ind0) - 1); 
   j = j + 1; 
   indj  = ind0.j; 
   depj  = dep0.indj; 
   sigj  = sig0.indj; 
   epsj  = epsi depj modjoi majoi; 
   vecj  = bsigma sigj modjoi majoi; 
   time  = time et (prog tim0. &lab1); 
 
   SI (EGA ncas 6); 
   uzp0  = uzp0 ET (PROG (EXTR (exco depj uy) scal (NOEUD 3))); 
   fzp0  = fzp0 ET (PROG (EXTR (exco vecj fy) scal (NOEUD 2))); 
   FINSI; 
 
   SI (EGA ncas 3); 
   uzp0  = uzp0 ET (PROG (EXTR depj 'UY' C)); 
   fzp0  = fzp0 ET (PROG (EXTR vecj 'FY' A)); 
   FINSI; 
 
   SI (OU (OU (OU (EGA ncas 1) (EGA ncas 2)) 
                  (EGA ncas 3)) (EGA ncas 4)); 
   uzp0  = uzp0 ET (PROG (EXTR depj 'UX' C)); 
   fzp0  = fzp0 ET (PROG (EXTR vecj 'FX' A)); 
   FINSI; 
 
 FIN lab1; 
  aa = extr fzp0 ( dime fzp0);
 
  err =  abs (aa +  9.32668E+06) / 9.32668E+06;
  message ' erreur relative ' err;
  si (err > 1.e-3);
    erreur (5);
     finsi;
 EV2 = EVOL MANU uzp0 (-1.0*fzp0); 
 si ( ega graph 'O'); 
 DESS EV2; 
 finsi;
 FIN ;