Télécharger la plaquette Cast3M


fr en

Télécharger ba1d.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : ba1d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************
  5. * DESCRIPTION
  6. ******************************************************
  7. * Local test - BA1D model
  8. *
  9. * Loading path considered:
  10. *
  11. * - : Cyclic bending
  12. *
  13. ******************************************************
  14. * AUTHOR
  15. ******************************************************
  16. * Developped by : Benjamin RICHARD 
  17. *                 CEA-DEN/DANS/DM2S/SEMT/EMSI
  18. *                 Benjamin.Richard@cea.fr
  19. *           
  20. *                 Laurent MOUTOUSSAMY
  21. *                 CEA-DEN/DANS/DM2S/SEMT/EMSI
  22. *                 Laurent.Moutoussamy@cea.fr
  23. ******************************************************
  24. graph = 'N';
  25. ******************************************************
  26.  
  27. ******************************************************
  28. * OPTIONS
  29. ******************************************************
  30. OPTI ELEM SEG2 DIME 2;
  31.  
  32. ******************************************************
  33. * GEOMETRICAL PARAMETERS
  34. ******************************************************
  35. DENSITE 1.0;
  36.  
  37. P1 = 0.0 0.0;
  38. P2 = 0.0 1.0;
  39.  
  40. L1 = D 1 P1 P2;
  41.  
  42. CL1 = BLOQ DEPL ROTA P1;
  43. CL2 = BLOQ UX        P2;
  44.  
  45. ******************************************************
  46. * MODEL
  47. ******************************************************
  48. MO1 = MODE L1 MECANIQUE ELASTIQUE PLASTIQUE BA1D POUT;
  49.  
  50. ******************************************************
  51. * MATERIAL
  52. ******************************************************
  53. MA1 = MATE MO1 YOUN 39.0E3 NU 0.0001 UELA 1.0E-3 FPLA 100 
  54. PFIS 0.3  HCIN 10.0 QFRA 0.5 APIH 1.0
  55. BPIH 5.0E10 SECT 1.0 INRZ 1.0 ;
  56.  
  57. ******************************************************
  58. * LOADING
  59. ******************************************************
  60. LI1 = PROG 0. 1. 2. 3.;
  61. LI2 = PROG 0. 2.5E-3 0. 2.5E-3;
  62. EV1 = EVOL MANU LI1 LI2;
  63.  
  64. DEP = DEPI CL2 1.0;
  65. CH1 = CHAR DIMP DEP EV1;
  66.  
  67. ******************************************************
  68. * COMPUTED PSEUDO TIMES
  69. ******************************************************
  70. LIC = PROG 0. PAS 0.01 3.0;
  71.  
  72. ******************************************************
  73. * PASAPAS PROCEDURE
  74. ******************************************************
  75. TAB1 = TABLE;
  76. TAB1.'BLOCAGES_MECANIQUES' = CL1 ET CL2;
  77. TAB1.'MODELE' = MO1;
  78. TAB1.'CHARGEMENT'=CH1;
  79. TAB1.'CARACTERISTIQUES' = MA1;
  80. TAB1.'TEMPS_CALCULES' = LIC;
  81. TAB1.'MOVA' = ((MOTS DP0) ET (MOTS DM0));
  82. TAB1.'HYPOTHESE_DEORMATIONS' = 'QUADRATIQUE';
  83. PASAPAS TAB1;
  84.  
  85. ******************************************************
  86. * POST TREATMENT
  87. ******************************************************
  88. NN = DIME TAB1.DEPLACEMENTS;
  89. NN = NN - 1;
  90.  
  91. ******************************************************
  92. * DEFORMED MESH
  93. ******************************************************
  94. DEF1 = DEFO L1 TAB1.DEPLACEMENTS.NN 0. BLEU;
  95. DEF2 = DEFO L1 TAB1.DEPLACEMENTS.NN    ROUG;
  96.  
  97. si (ega graph 'O');
  98. TRAC (DEF1 ET DEF2);
  99. finsi;
  100.  
  101. ******************************************************
  102. * LOAD/DISPLACEMENT CURVE
  103. ******************************************************
  104. EV2 = @GLOBAL TAB1 EV1 CL2 FX;
  105.  
  106. si (ega graph 'O');
  107. dess ev2;
  108. finsi;
  109.  
  110. ******************************************************
  111. * IMPLEMENTATION TEST
  112. ******************************************************
  113. aa = EXTR EV2 'ORDO';
  114. aaa = EXTR aa 290;
  115.  
  116. message 'erreur relative';
  117. err = ABS(aaa - 98.981) / 98.981;
  118. list err;
  119.  
  120. si (> err 1E-3);
  121. erreur(5);
  122. finsi;
  123.  
  124. FIN;
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  

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