Télécharger cas_test_dp2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : cas_test_dp2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *** essai de traction simple
  5. *** Cas test
  6.  
  7. * Type & Elements de representation *
  8. opti dime 3 elem cub8 ;
  9.  
  10. * Caracteristiques geometriques *
  11. la=1.;
  12. mla=1. ;
  13. crit = 0.00002 ;
  14.  
  15. * Defintion des noeuds *
  16. P1 = 0. 0. 0. ;
  17. P2 = 0. la 0. ;
  18.  
  19. * Definition des contours *
  20. L1 = droit 1 P1 P2 ;
  21.  
  22. * Generation de contours *
  23. coteg = L1 TRAN (0. 0. la) DINI mla DFIN mla;
  24.  
  25. * Creation du volume *
  26. geom = VOLU coteg TRAN (la 0. 0.) DINI la DFIN la ;
  27.  
  28. * Defintion des surfaces supplementaires *
  29. coted=FACE geom 2 ;
  30.  
  31. * Epurer le volume *
  32. elim crit geom ;
  33. elim crit (geom et coteg);
  34. elim crit (geom et coted);
  35.  
  36. * Affichage *
  37. coted = coted coul rouge ;
  38. coteg = coteg coul bleu ;
  39.  
  40. * constantes materiaux
  41.  
  42. yg = 31.25e9 ;
  43. nub = 0.2 ;
  44. gf = 250 ;
  45. ltr = 3.5e6 ;
  46. lcs =35e6 ;
  47. lbi = 1.16 * lcs ;
  48. sigy = 20.e6 ;
  49. epm = -2.e-3;
  50. epu = -6.e-3 ;
  51.  
  52. mod1 = MODE geom MECANIQUE ELASTIQUE PLASTIQUE DRUCKER_PRAGER_2 ;
  53.  
  54. jaco1 = jaco mod1 ;
  55. jaco2 = chan 'RIGIDITE' mod1 jaco1 ;
  56. lcar = jaco2**(1./3.) ;
  57.  
  58. mat1 = MATE mod1 'YOUN' yg 'NU' nub 'GF' gf
  59. 'LTR' ltr 'LCS' lcs 'LBI' lbi 'SIGY' sigy 'EPM' epm
  60. 'EPU' epu 'LCAR' lcar ;
  61.  
  62. *** Conditions aux limites
  63.  
  64. cl1 = BLOQ P1 UZ UY;
  65. cl2 = BLOQ L1 UZ ;
  66. cl3 = BLOQ COTEG UX ;
  67. cldepl = BLOQ COTED UX ;
  68.  
  69. cltot = cl1 et cl2 et cl3 et cldepl ;
  70.  
  71. *** chargement
  72.  
  73. fo1 = DEPI cldepl 1. ;
  74. prog1 = PROG 0. 1. ;
  75. prog2 = PROG 0. 7.5e-4 ;
  76.  
  77. step1 = 1./20. ;
  78. list1 = PROG 0. PAS step1 1 ;
  79.  
  80. evol1 = EVOL MANU prog1 prog2 ;
  81.  
  82.  
  83. char1 = CHARGEMENT fo1 evol1 DIMP ;
  84.  
  85. tab1 = TABLE ;
  86. tab1 . PLASTIQUE = vrai ;
  87. tab1 . MODELE = mod1 ;
  88. tab1 . CARACTERISTIQUES=mat1 ;
  89. tab1 . CHARGEMENT = char1 ;
  90. tab1 . BLOCAGES_MECANIQUES = cltot ;
  91. tab1 . TEMPS_CALCULES = list1 ;
  92. tab1 . TEMPS_SAUVES = list1 ;
  93.  
  94. PASAPAS tab1 ;
  95.  
  96. ind = inde (tab1.contraintes);
  97. nbind1 = dime ind;
  98. nbind1 = nbind1 - 1 ;
  99. ************************************************************
  100.  
  101. *definition des fonctions *
  102.  
  103. sigxx1=prog;
  104. dex1=prog;
  105. dev = prog ;
  106.  
  107. i = 1;
  108.  
  109. repeter bou1 nbind1;
  110.  
  111. cont1 = (tab1).contraintes.i;
  112.  
  113. sigxx1=sigxx1 et (prog(extr cont1 SMXX 1 1 1));
  114.  
  115. dep1 = (tab1) . deplacements . i ;
  116.  
  117. eps0 = EPSI dep1 mod1;
  118. depx1 = (EXTR eps0 EPXX 1 1 1);
  119. depy1 = (EXTR eps0 EPYY 1 1 1);
  120. depz1 = (EXTR eps0 EPZZ 1 1 1) ;
  121.  
  122. * opti donn 5 ;
  123. ************************************************************
  124. * deformations suivant les 3 composantes *
  125.  
  126. devv = (depx1 + depy1 + depz1) ;
  127. dev = dev et (prog devv) ;
  128. dex1 = dex1 et (prog depx1);
  129. i = i+1;
  130.  
  131. fin bou1;
  132.  
  133. sixx = evol manu 'depl_xx' dex1 'sigxx' sigxx1;
  134. sivv = evol manu 'defo_vol' dev 'sigxx' sigxx1 ;
  135.  
  136. *dess sivv;
  137. *dess sixx;
  138. fin ;
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  

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