Télécharger la plaquette Cast3M


fr en

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

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