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*YG*(EPS0**2)/2.);
  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 BETA 1.06 ;
  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 ) coul &bou1;
  93. SI (&BOU1 EGA 1);
  94. EVT=EV1;
  95. SINO;
  96. EVT=EV1 ET EVT;
  97. FINSI;
  98. G2=INTG EV1;
  99. PROGG=INSE PROGG &BOU1 G2;
  100. TAB2.GFCAL.I=G2;
  101. FIN BOU1;
  102.  
  103. *
  104. *
  105. *
  106. MESSAGE 'VALEURS DE GF2 OBTENUES';
  107. LIST PROGG;
  108. IERR=0;
  109. REPETER BOU3 N;
  110. XERR=((EXTR PROGG &BOU3) - GF)/GF;
  111. SI (XERR > 0.01); IERR=5;FINSI;
  112. MESSAGE 'TAILLE 10-' &BOU3 ' ERREUR DE CALCUL : ' (100*XERR) '%';
  113. FIN BOU3;
  114. ERREUR IERR;
  115. SI GRAPH;
  116. DESS EVT;
  117. FINSI;
  118. FIN;
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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