Télécharger gdef2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gdef2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************
  5. * *
  6. * TEST GDEF2 *
  7. * *
  8. * CISAILLEMENT PUR EN GRANDES DEFORMATIONS PLASTIQUES *
  9. * ( MODELE DE PLASTICITE ISOTROPE ) *
  10. * On compare avec les valeurs obtenues a une solution *
  11. * analytique *
  12. * *
  13. ********************************************************************
  14. *
  15. OPTI ECHO 0;
  16. TITRE ' TEST GDEF2 ' ;
  17. GRAPH = 'N';
  18. OPTI DIME 2 ELEM QUA4 MODE PLAN CONT;
  19. *
  20. *==================================================================*
  21. * geometrie : maillage *
  22. *==================================================================*
  23. *
  24. P1 = 0 0 ;
  25. P2 = 1. 0. ;
  26. P3 = 1. 1. ;
  27. P4 = 0. 1. ;
  28. L1 = P1 D 1 P2 ;
  29. L2 = P2 D 1 P3;
  30. L3 = P3 D 1 P4 ;
  31. L4 = P4 D 1 P1 ;
  32. SU = DALLER L1 L2 L3 L4 ;
  33. SI (NEG GRAPH 'N');
  34. TRAC SU ;
  35. FINSI;
  36. *
  37. *==================================================================*
  38. * modele - materiau - caracteristiques en Pa *
  39. * conditions aux limites *
  40. *==================================================================*
  41. *
  42. MO = MODE SU MECANIQUE ELASTIQUE PLASTIQUE ISOTROPE
  43. epsilon jaumann ;
  44. YU = 1000. ; XNU = 0.3 ;
  45. LTRAC = EVOL MANU 'EPS' ( PROG 0. 1.E-3 10.001 )
  46. 'SIG' ( PROG 0. 1. 30.91026920);
  47. MA = MATE MO YOUN YU NU XNU TRAC LTRAC ;
  48. BL1 = BLOQ DEPL P1 ;
  49. BL2 = BLOQ DEPL P2 ;
  50. BL3 = BLOQ UX L3 ;
  51. BL4 = BLOQ UY L3 ;
  52. FF = DEPI BL3 1. ;
  53. *
  54. *==================================================================*
  55. * Initialisation de la table pour appel a PASAPAS *
  56. *==================================================================*
  57. *
  58. TAB1 = TABLE;
  59. tab1.grands_deplacements=vrai;
  60. TAB1 . BLOCAGES_MECANIQUES= BL1 ET BL2 ET BL3 ET BL4 ;
  61. TAB1 . MODELE = MO;
  62. TAB1 . CARACTERISTIQUES = MA ;
  63. EV1 = EVOL MANU X (PROG 0 5 10 15 20 )
  64. Y (PROG 0 5 10 5 0 );
  65. TAB1 . CHARGEMENT = CHAR MECA FF EV1 ;
  66. TAB1 . TEMPS_CALCULES = PROG 0.1 PAS 0.1 5.;
  67. PASAPAS TAB1 ;
  68. *
  69. *==================================================================*
  70. * affichage des resultats *
  71. *==================================================================*
  72. *
  73. MESS ' RESULTATS TEST GDEF2 ';
  74. MESS ' -------------------- ';
  75. MESS ' On cherche le cisaillement ';
  76. MESS ' Il est compare a une solution analytique ';
  77. *
  78. *----------- r{sultats --------------------------------------------------
  79. *
  80. AA = TAB1 . DEPLACEMENTS ;
  81. BB = TAB1 . CONTRAINTES ;
  82. N = DIME AA ;
  83. SP12= PROG ;
  84. SPTH12= PROG ;
  85. EP = PROG ;
  86. EP = TAB1 . TEMPS_CALCULES ;
  87. REPE BOU (N - 1 ) ;
  88. CSX = CHAN CHPO MO ( BB . &BOU ) ;
  89. DX = EXTR ( AA . &BOU ) UX P3 ;
  90. S12 = EXTR CSX SMXY P1 ;
  91. SP12 = SP12 ET ( PROG S12 ) ;
  92. SDX = 0.577350269 + DX ;
  93. SPTH12= SPTH12 ET ( PROG SDX ) ;
  94. FIN BOU ;
  95. EV12 = EVOL ROSE MANU 'DEPLA' EP 'SIG12' SP12 ;
  96. FVTH12 = EVOL TURQ MANU 'DEPLA' EP 'SIG12' SPTH12 ;
  97. SI (NEG GRAPH 'N');
  98. DESS ( FVTH12 ET EV12) ;
  99. FINSI;
  100. *
  101. * COMPARAISON SOL ANALYTIQUE
  102. * ecart
  103. *
  104. ER12 =(MAXI ABS ( SP12 - SPTH12 ))/(MAXI ABS SPTH12);
  105. *
  106. MESS ' Ecart relatif en contraintes sig12 :' ER12 ;
  107. *
  108. *
  109. *==================================================================*
  110. * code fonctionnement *
  111. *==================================================================*
  112. *
  113. * L'ecart maximum entre valeurs de reference et calculee doit etre
  114. * inferieur a 5.d-2
  115. *
  116. SAUT 2 LIGNE;
  117. SI (ER12 <EG 0.05 ) ;
  118. ERRE 0;
  119. SINON;
  120. ERRE 5;
  121. FINSI;
  122. *
  123. *
  124. FIN;
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  

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