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. CHAM31= CHH*FT/GF ;
  49. MAT1=MATE MOD1 YOUN YG NU 0.2D0 ATRA -1. BTRA CHAM31 KTR0 EPS0
  50. ACOM 2000. BCOM 20000. RHO 2.5E+3 ;
  51. CL1=BLOQ P1 UY;
  52. CL2=BLOQ D2 UX;
  53. CL3=BLOQ D4 UX;
  54. CLTOT=CL1 ET CL2 ET CL3;
  55. F1=DEPI CL2 1.;
  56. EVOL1=EVOL MANU (PROG 0. 1.) (PROG 0. 1.);
  57. CHAR1=CHARGEMENT F1 EVOL1 DIMP;
  58. TAB1=TABLE;
  59. TAB1.BLOCAGES_MECANIQUES=CLTOT;
  60. TAB1.CARACTERISTIQUES=MAT1;
  61. TAB1.GRAND_DEPLACEMENTS=VRAI;
  62. TAB1.MODELE=MOD1;
  63. TAB1.MOVA='D';
  64. TAB1.CHARGEMENT=CHAR1;
  65. TAB1.MAXITERATION=100;
  66.  
  67. SI (&BOU1 EGA 1);
  68. TAB1.TEMPS_CALCULES=PROG 0. (A*EPS0) PAS (A*0.1*EPS0) (A*10*EPS0)
  69. PAS (A*EPS0) (A*30*EPS0);
  70.  
  71. FINSI;
  72. SI (&BOU1 EGA 2);
  73. TAB1.TEMPS_CALCULES=PROG 0. (A*EPS0) PAS (A*0.1*EPS0) (A*10*EPS0)
  74. PAS (A*EPS0) (A*100*EPS0)
  75. PAS (A*10*EPS0) (A*300*EPS0);
  76.  
  77. FINSI;
  78. SI (&BOU1 EGA 3);
  79. TAB1.TEMPS_CALCULES=PROG 0. (A*EPS0) PAS (A*0.1*EPS0) (A*10*EPS0)
  80. PAS (A*EPS0) (A*100*EPS0)
  81. PAS (A*10*EPS0) (A*300*EPS0)
  82. PAS (A*100*EPS0) (A*3000*EPS0);
  83. FINSI;
  84.  
  85.  
  86. PASAPAS TAB1;
  87.  
  88. N1=DIME TAB1.TEMPS_CALCULES;
  89. PROGT=TAB1.TEMPS_CALCULES;
  90. TITRE 'COURBE CONTRAINTE - DEPLACEMENT';
  91. EV1=(@GLOBAL TAB1 EVOL1 CL2 FX) / A;
  92. SI (&BOU1 EGA 1);
  93. EVT=EV1;
  94. SINO;
  95. EVT=EV1 ET EVT;
  96. FINSI;
  97. L1=ENLE (EXTR EV1 ABSC) 1;
  98. L2=ENLE (EXTR EV1 ORDO) 1;
  99. EV1=EVOL MANU L1 L2;
  100. S1=INTG EV1;
  101. G2=EXTR S1 1;
  102. PROGG=INSE PROGG &BOU1 G2;
  103. TAB2.GFCAL.I=G2;
  104. FIN BOU1;
  105.  
  106. *
  107. *
  108. *
  109. MESSAGE 'VALEURS DE GF OBTENUES';
  110. LIST PROGG;
  111. IERR=0;
  112. REPETER BOU3 N;
  113. XERR=((EXTR PROGG &BOU3) - GF)/GF;
  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.  

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