Télécharger mazars.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : mazars.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************************
  5. * *
  6. * TEST; MATERIAU ENDOMMAGEMENT MAZARS *
  7. * *
  8. * *
  9. * EPROUVETTE EN TRACTION AVEC DEPLACEMENTS IMPOSES *
  10. * *
  11. * CONTRAINTES PLANES. *
  12. * *
  13. * LES RESULTATS OBTENUS SONT COMPARES A LA SOLUTION ANALYTIQUE *
  14. * *
  15. ******************************************************************* *
  16.  
  17.  
  18.  
  19. opti echo 0 dime 2 elem qua4 mode plan cont;
  20. * ajout de option epsilon lineaire pour la precision des test!
  21. OPTION epsilon lineaire;
  22.  
  23. densite 1.;
  24. p1=0. 0.;
  25. p2=1. 0.;
  26. p3=1. 1.;
  27. p4=0. 1.;
  28. d1=d p1 p2;d2=d p2 p3;d3=d p3 p4;d4=d p4 p1;
  29. surf1=daller d1 d2 d3 d4 plane;
  30. at=1.e0;
  31. ac=1.2e0;
  32. bc=1.5e+03;
  33. bt=10000.;
  34. epsd0=1.e-04;
  35. YG=31.E+9;
  36. mod1=mode surf1 mecanique elastique isotrope endommagement mazars;
  37. mat1=MATE mod1 YOUN YG NU 0.2 RHO 2.3E+03
  38. btra bt ktr0 epsd0 atra at acom ac bcom bc beta 1.06e0;
  39. blo1=bloq ux p1;
  40. blo2=bloq uy d1;
  41. blo3=bloq uy d3;
  42. cltot=blo1 et blo2 et blo3 ;
  43. fo1 = depi blo3 1. ;
  44. evol1=evol manu temps (prog 0. 1. )
  45. depl (prog 0. 1 );
  46. char1=charg evol1 fo1 DIMP;
  47. lit1=prog 1.e-4 pas 25.e-6 5.e-4 ;
  48. tab1=table;
  49. tab1.mova='D ';
  50. tab1.blocages_mecaniques=cltot;
  51. tab1.caracteristiques=mat1;
  52. tab1.modele=mod1;
  53. tab1.chargement=char1;
  54. tab1.temps_calcules=lit1;
  55. pasapas tab1;
  56. *opti donn 5;
  57. tab2=table;
  58. tab2.1='MARQ CROI';
  59. ev1=@global tab1 blo3 evol1 fy ;
  60. *
  61. *construction du resultat theorique
  62. *
  63. nb= dime lit1;
  64. UN=prog nb*1.d0;
  65. LAT = prog nb * AT;
  66. LEPS=prog nb*epsd0;
  67. DT=UN - ( ( EPSD0* (UN-LAT) /lit1 ) + (AT * ( EXP (BT*(LEPS-lit1)))));
  68. lsig0=(un - dt) * YG * lit1;
  69. lsig1 = enle (extr ev1 ordo) 1;
  70. *
  71. *comparaison des deux resultats
  72. *
  73. err1=maxi (( lsig0 - lsig1) / lsig0);
  74. temps;
  75. si (err1 < 1.e-10);
  76. ERRE 0;
  77. sinon ;
  78. ERRE 5;
  79. finsi;
  80. fin;;
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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