Télécharger la plaquette Cast3M


fr en

Télécharger desmorat.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : desmorat.dgibi 
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. OPTI dime 3 elem CUB8  echo 0;
  6. *-------------------définition de la géométrie--------------------
  7. P1 = 0. 0. 0.;
  8. P2 = 0. 0.1 0.;
  9. P3 = 0.1 0.1 0.;
  10. P4 = 0.1 0. 0.;
  11. V1 = 0. 0. 0.1;
  12. P5 = 0.1 0. 0.1;
  13. Li1 = D 1 P1 P2 ; 
  14. Li2 = D 1 P2 P3;
  15. Li3 = D 1 P3 P4;
  16. Li4 = D 1 P4 P1;
  17. Ltot = Li1 ET Li2 ET Li3 ET li4;
  18. Surf1 = SURF Ltot PLANE;
  19. Vovo = VOLU Surf1 1  TRANS V1;
  20. Surf2 = FACE 2 Vovo;
  21. Surf3 = Vovo POIN PLAN P1 P4 P5 0.001;
  22. Surf4 = Vovo POIN PLAN P3 P4 P5 0.001;
  23. ELIM 0.01 Vovo ;
  24. *-----------MODELE------------------    
  25.     Mod1= MODE Vovo MECANIQUE ELASTIQUE ENDOMMAGEMENT DESMORAT CUB8;
  26.     kappa0 = 5E-05 ;
  27.     Mat1 = MATE Mod1 RHO 2400 YOUN 42.E09 NU 0.2 'K0' kappa0 'A' 5E+3 
  28.     'a' 2.93E-4 
  29.     'etaC' 0. 'etaT' 3. 'Dc' 0.99; 
  30.  
  31.  
  32.     list mod1;
  33. *-------------CONDITIONS LIMITES-------------
  34. CL1 = BLOQ Li1 Ux;
  35. CL2 = BLOQ P2 Uy;
  36. CL3 = BLOQ Surf1 Uz ;
  37. Cldep = BLOQ Surf2 Uz;
  38. *-------------ASSEMBLAGE CONDITIONS LIMITES
  39. CLtot = CL1 ET CL2 ET CL3  ET  CLdep;
  40. *CLtot = CL3  ET  CLdep;
  41. *--------------CHARGEMENT---------------
  42. FO1 = DEPI Cldep -1;
  43. LI1 = PROG 0. 1e-2;
  44. LI2 = PROG 0  1e-3 ;
  45. EV = EVOL MANU LI1 LI2 ;
  46. CHA1 =  CHAR 'DIMP' FO1 EV;
  47. LIS1 = PROG 0. PAS 1e-5 0.3e-2 ;
  48. TAB1 = TABLE ;
  49. TAB1.'BLOCAGES_MECANIQUES' = CLtot; 
  50. TAB1.'MODELE' = Mod1;
  51. TAB1.'CHARGEMENT' = CHA1;
  52. TAB1.'CARACTERISTIQUES' = Mat1;
  53. TAB1.'TEMPS_CALCULES' = LIS1;
  54. TAB1.'MOVA' = 'D11';
  55. TAB1.DYNAMIQUE=faux ;
  56. *--------------RESOLUTION--------------
  57. PASAPAS TAB1;
  58.  
  59. n= dime (tab1.deplacements) - 1;
  60. var = tab1.variables_internes.n ;
  61. tra1= extr var 'D11 ' 1 1 1 ;
  62. tra2= extr var 'D22 ' 1 1 1 ;
  63. tra3= extr var 'D33 ' 1 1 1 ;
  64.  
  65. D11 = 0.90687;
  66. D22 = 0.90589;
  67. D33 = 4.18284E-22;
  68.  
  69. err1 = abs ((TRA1 - D11)/D11 * 100) ;
  70. err2 = abs ((TRA2 - D22)/D22 * 100) ;
  71. *err3 = abs ((TRA3 - D33)/D33 * 100) ;
  72. err3 = (abs ( d33 - tra3)) * 1.e9;
  73.  
  74. mess ' err1 ' err1  ' err2 ' err2 ' err3 ' err3;
  75. si ((err1 > .06) ou (err2 > .06) ou (err3 > .1));
  76.  erre 5;
  77. sinon ;
  78.  erre 0 ;
  79. finsi ;
  80.  
  81. fin;
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  

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