Télécharger la plaquette Cast3M


fr en

Télécharger mazars2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  mazars2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *   Christian La Borderie 20 janvier 2010
  5. *
  6. *
  7. *
  8. *TESTE LA VERSION MAZARS MODIFIÉE AVEC PRISE EN COMPTE DE GF
  9. *ON SOLLICITE UN ÉLÉMENT DE TAILLE VARIABLE EN TRACTION ET ON VÉRIFIE
  10. *QUE L'ÉNERGIE DE FISSURATION OBTENUE CORRESPOND À CELLE QUI ÉTAIT VOULUE
  11. *
  12. OPTI ECHO 0;
  13. COMPLET=FAUX;
  14. GRAPH=FAUX;
  15. *
  16. *
  17. *
  18. OPTI DIME 2 ELEM QUA4 MODE PLAN CONT ;
  19. SI COMPLET;N=3;SINON;N=1;FINSI;
  20. TAB2=TABLE;
  21. TAB2.BTRA=TABLE;
  22. TAB2.TAILLE=TABLE;
  23. TAB2.EVTRAC=TABLE;
  24. TAB2.GFCAL=TABLE;
  25. TAB2.VARI=TABLE;
  26. PROGG=PROG;
  27. PROGD=PROG;
  28. REPETER BOU1 N;
  29.   A=10.**(-1*&BOU1);
  30.   TAB2.TAILLE.I=A;
  31.   PROGD=INSE PROGD &BOU1 A;
  32.   DENSITE A;
  33.   P1=0. 0.;
  34.   P2=A 0.;
  35.   P3=A A;
  36.   P4=0. A;
  37.   D1=D P1 P2;D2=D P2 P3;D3=D P3 P4;D4=D P4 P1;
  38.   S1=DALLER D1 D2 D3 D4;
  39.   MOD1=MODE S1 MECANIQUE ELASTIQUE ISOTROPE ENDOMMAGEMENT MAZARS;
  40.   FT=3.D+6;
  41.   YG=30.D+9;
  42.   GF=100.;
  43.   XNU=0.2;
  44.   EPS0=FT/YG;
  45.   CHAM1=MANU CHML MOD1 BTRA 1. TYPE CARACTERISTIQUES RIGIDITE;
  46.   CHAM2=INTG CHAM1 MOD1 ELEM BTRA;
  47.   CHH=CHAM2**(1./2.); 
  48.   GF2=GF - (A*(EPS0**2)/(2.*YG));
  49.   CHAM31= CHH*FT/GF2 ;
  50.   MAT1=MATE MOD1 YOUN YG NU 0.2D0 ATRA -1. BTRA CHAM31 KTR0 EPS0 
  51.                             ACOM 2000. BCOM 20000. RHO 2.5E+3 ;
  52.   CL1=BLOQ P1 UY;
  53.   CL2=BLOQ D2 UX;
  54.   CL3=BLOQ D4 UX;
  55.   CLTOT=CL1 ET CL2 ET CL3;
  56.   F1=DEPI CL2 1.;
  57.   EVOL1=EVOL MANU (PROG 0. 1.)   (PROG 0.  1.);
  58.   CHAR1=CHARGEMENT F1 EVOL1 DIMP;
  59.   TAB1=TABLE;
  60.   TAB1.BLOCAGES_MECANIQUES=CLTOT;
  61.   TAB1.CARACTERISTIQUES=MAT1;
  62.   TAB1.GRAND_DEPLACEMENTS=VRAI;
  63.   TAB1.MODELE=MOD1;
  64.   TAB1.MOVA='D';
  65.   TAB1.CHARGEMENT=CHAR1;
  66.   TAB1.MAXITERATION=100;
  67.  
  68.   SI (&BOU1 EGA 1);
  69.   TAB1.TEMPS_CALCULES=PROG 0. (A*EPS0) PAS (A*0.1*EPS0) (A*10*EPS0)
  70.                                       PAS (A*EPS0)    (A*30*EPS0);
  71.  
  72.   FINSI;
  73.   SI (&BOU1 EGA 2);
  74.   TAB1.TEMPS_CALCULES=PROG 0. (A*EPS0) PAS (A*0.1*EPS0) (A*10*EPS0)
  75.                                       PAS (A*EPS0)    (A*100*EPS0)
  76.                                       PAS (A*10*EPS0)  (A*300*EPS0);
  77.  
  78.   FINSI;
  79.   SI (&BOU1 EGA 3);
  80.   TAB1.TEMPS_CALCULES=PROG 0. (A*EPS0) PAS (A*0.1*EPS0) (A*10*EPS0)
  81.                                       PAS (A*EPS0)    (A*100*EPS0)
  82.                                       PAS (A*10*EPS0)  (A*300*EPS0)
  83.                                       PAS (A*100*EPS0) (A*3000*EPS0);
  84.   FINSI;
  85.  
  86.  
  87.   PASAPAS TAB1;
  88.  
  89.   N1=DIME TAB1.TEMPS_CALCULES;
  90.   PROGT=TAB1.TEMPS_CALCULES;
  91.   TITRE 'COURBE CONTRAINTE - DEPLACEMENT';
  92.   EV1=(@GLOBAL TAB1 EVOL1 CL2 FX) / A;
  93.   SI (&BOU1 EGA 1);
  94.     EVT=EV1;
  95.     SINO;
  96.     EVT=EV1 ET EVT;
  97.   FINSI;
  98.   L1=ENLE (EXTR EV1 ABSC) 1;
  99.   L2=ENLE (EXTR EV1 ORDO) 1;
  100.   EV1=EVOL  MANU L1 L2;
  101.   G2=INTG EV1;
  102.   PROGG=INSE PROGG &BOU1 G2;
  103.   TAB2.GFCAL.I=G2;
  104. FIN BOU1;
  105.  
  106. *
  107. *
  108. *
  109. MESSAGE 'VALEURS DE GF2 OBTENUES';
  110. LIST PROGG;
  111. IERR=0;
  112. REPETER BOU3 N;
  113.  XERR=((EXTR PROGG &BOU3) - GF2)/GF2;
  114.  SI (XERR > 0.01); IERR=5;FINSI;
  115.  MESSAGE 'TAILLE  10-' &BOU3 '  ERREUR DE  CALCUL : ' (100*XERR) '%';
  116. FIN BOU3;
  117. ERREUR IERR;
  118. SI GRAPH;
  119.   DESS EVT;
  120. FINSI;
  121. FIN;
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  

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