Télécharger la plaquette Cast3M


fr en

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 1;
  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. YU = 1000. ;   XNU = 0.3 ;
  44. LTRAC = EVOL MANU 'EPS' ( PROG 0. 1.E-3 10.001 )
  45.                   'SIG' ( PROG 0. 1. 30.91026920);
  46. LEPS1 = LTRAC extr absc ;
  47. LSIG1 = LTRAC extr ordo ;
  48. lsig1 = lsig1 enle 1 ;
  49. leps1 = (leps1 enle 1) - (lsig1 / yu) ;
  50. LECRO = evol vert manu eps leps1 sig lsig1 ;
  51. si (neg graph 'N')
  52.   dess (ltrac et lecro) titr ' Courbes de traction et d ecrouissage (vert)' ;
  53. fins ;
  54. MA = MATE MO YOUN YU     NU XNU ECRO LECRO ;
  55. BL1 = BLOQ  DEPL P1 ;
  56. BL2 = BLOQ  DEPL P2 ;
  57. BL3 = BLOQ UX L3 ;
  58. BL4 = BLOQ UY L3 ;
  59. FF = DEPI BL3 1. ;
  60. *                                                                               
  61. *==================================================================*            
  62. *         Initialisation de la table pour appel a PASAPAS          *
  63. *==================================================================*            
  64. *                                                                               
  65. TAB1 = TABLE;
  66. tab1.grands_deplacements=vrai;
  67. TAB1 . BLOCAGES_MECANIQUES= BL1 ET BL2 ET BL3 ET BL4 ;
  68. TAB1 . MODELE = MO;
  69. TAB1 . CARACTERISTIQUES = MA ;
  70. EV1 = EVOL MANU X (PROG 0  5 10 15 20 )
  71.  Y  (PROG 0  5 10 5 0   );
  72. TAB1 . CHARGEMENT = CHAR MECA FF EV1 ;
  73. TAB1 . TEMPS_CALCULES = PROG 0.1   PAS 0.1 5.;
  74. TAB1.'HYPOTHESE_DEFORMATIONS' = MOT 'JAUMANN' ;
  75. PASAPAS TAB1 ;
  76. *                                                                               
  77. *==================================================================*            
  78. *                   affichage des resultats                        *            
  79. *==================================================================*            
  80. *                                                                               
  81. MESS ' RESULTATS TEST GDEF2 ';
  82. MESS ' -------------------- ';                                                  
  84. MESS '   On cherche le cisaillement ';
  85. MESS '   Il est compare a une solution analytique   ';
  87. *                                                                               
  88. *----------- r{sultats --------------------------------------------------       
  89. *                                                                               
  90. AA = TAB1 . DEPLACEMENTS ;
  91. BB = TAB1 . CONTRAINTES  ;
  92. N = DIME AA ;
  93. SP12= PROG ;
  94. SPTH12= PROG ;
  95. EP = PROG ;
  96. EP = TAB1 . TEMPS_CALCULES ;
  97. REPE BOU (N - 1 ) ;
  98. CSX = CHAN CHPO  MO ( BB . &BOU ) ;
  99. DX = EXTR ( AA . &BOU ) UX P3 ;
  100. S12 = EXTR CSX SMXY P1 ;
  101. SP12 = SP12 ET ( PROG S12 ) ;
  102. SDX = 0.577350269 +  DX  ;
  103. SPTH12= SPTH12 ET ( PROG  SDX ) ;
  104. FIN BOU ;
  105. EV12 = EVOL ROSE  MANU 'DEPLA' EP 'SIG12' SP12 ;
  106. FVTH12 = EVOL TURQ  MANU 'DEPLA' EP 'SIG12' SPTH12 ;
  107. SI (NEG GRAPH 'N');                                                             
  108.  DESS  ( FVTH12 ET EV12) ;
  109. FINSI;
  110. *
  111. *   COMPARAISON SOL ANALYTIQUE
  112. *   ecart
  113. *
  114. ER12 =(MAXI ABS ( SP12 - SPTH12 ))/(MAXI ABS SPTH12);
  115. *
  117. MESS ' Ecart  relatif en contraintes  sig12     :' ER12 ;
  118. *                                                                               
  119. *                                                                               
  120. *==================================================================*            
  121. *                  code fonctionnement                             *            
  122. *==================================================================*            
  123. *                                                                               
  124. * L'ecart maximum entre valeurs de reference et calculee doit etre
  125. * inferieur a 5.d-2
  126. *                                                                               
  127. SAUT 2 LIGNE;                                                                   
  128.     SI (ER12 <EG 0.05 ) ;
  129.         ERRE 0;
  130.     SINON;
  131.         ERRE 5;
  132.     FINSI;
  133. *                                                                               
  135. *                                                                               
  136. FIN;                                                                            
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  

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