Télécharger compression.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : compression.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Type & Elements de representation *
  5. opti dime 3 elem cub8 ;
  6.  
  7. * Caracteristiques geometriques *
  8. la=1.;
  9. mla=1. ;
  10. crit = 0.00002 ;
  11.  
  12. * Defintion des noeuds *
  13. P1 = 0. 0. 0. ;
  14. P2 = 0. la 0. ;
  15.  
  16. * Definition des contours *
  17. L1 = droit 1 P1 P2 ;
  18.  
  19. * Generation de contours *
  20. coteg = L1 TRAN (0. 0. la) DINI mla DFIN mla;
  21.  
  22. * Creation du volume *
  23. geom = VOLU coteg TRAN (la 0. 0.) DINI la DFIN la ;
  24.  
  25. * Defintion des surfaces supplementaires *
  26. coted=FACE geom 2 ;
  27.  
  28. * Epurer le volume *
  29. elim crit geom ;
  30. elim crit (geom et coteg);
  31. elim crit (geom et coted);
  32.  
  33. * Affichage *
  34. coted = coted coul rouge ;
  35. coteg = coteg coul bleu ;
  36.  
  37. *** Modele
  38.  
  39. mod1 = MODE geom MECANIQUE ELASTIQUE PLASTIQUE_ENDOM 'ENDO_PLAS' ;
  40.  
  41. mat1 = MATE mod1 'YOUN' 20.74e9 'NU' 0.2 'AC' 2.75
  42. 'AT' 1. 'BC' 2470. 'BT' 15000. 'EPD0' 10.e-5
  43. 'RC' 120.e6 'RT' 10.e6 'P' .4 'GAMA' 0.99 'ALFA' 0.5
  44. 'A' -1 'AH' 7.e-5 'BH' 2.e-2 'CH' 3.e-6 'BETA' 1. 'K0' 0.1;
  45.  
  46.  
  47. *** Conditions aux limites
  48.  
  49. cl1 = BLOQ P1 UZ UY;
  50. cl2 = BLOQ L1 UZ ;
  51. cl3 = BLOQ COTEG UX ;
  52. cldepl = BLOQ COTED UX ;
  53.  
  54. cltot = cl1 et cl2 et cl3 et cldepl ;
  55.  
  56.  
  57. *** Chargement
  58.  
  59. fo1 = DEPI cldepl 1. ;
  60. prog1 = PROG 0. 1. ;
  61. prog2 = PROG 0. -2.e-3 ;
  62.  
  63. step1 = 1./20. ;
  64. list1 = PROG 0. PAS step1 1. ;
  65.  
  66. evol1 = EVOL MANU prog1 prog2 ;
  67.  
  68. char1 = CHARGEMENT fo1 evol1 DIMP ;
  69.  
  70.  
  71. tab1 = TABLE ;
  72. tab1 . PLASTIQUE = vrai ;
  73. tab1 . MODELE = mod1 ;
  74. tab1 . CARACTERISTIQUES = mat1 ;
  75. tab1 . CHARGEMENT = char1 ;
  76. tab1 . BLOCAGES_MECANIQUES = cltot ;
  77. tab1 . TEMPS_CALCULES = list1 ;
  78. tab1 . TEMPS_SAUVES = list1 ;
  79. tab1 . MOVA = 'D' ;
  80. tab1 . PRECISION = 1.e-8;
  81.  
  82. PASAPAS tab1 ;
  83.  
  84. ind = inde (tab1.contraintes);
  85. nbind1 = dime ind;
  86. nbind1 = nbind1 - 1 ;
  87. ************************************************************
  88.  
  89. *definition des fonctions *
  90.  
  91. sigxx1=prog;
  92. dex1=prog;
  93. dev = prog ;
  94.  
  95. i = 1;
  96.  
  97. repeter bou1 nbind1;
  98.  
  99. cont1 = (tab1).contraintes.i;
  100.  
  101. sigxx1=sigxx1 et (prog(extr cont1 SMXX 1 1 1));
  102.  
  103. dep1 = (tab1) . . i ;
  104.  
  105. eps0 = EPSI dep1 mod1;
  106. depx1 = (EXTR eps0 EPXX 1 1 1);
  107. depy1 = (EXTR eps0 EPYY 1 1 1);
  108. depz1 = (EXTR eps0 EPZZ 1 1 1) ;
  109.  
  110. * opti donn 5 ;
  111. ************************************************************
  112. * deformations suivant les 3 composantes *
  113.  
  114. devv = (depx1 + depy1 + depz1) ;
  115. dev = dev et (prog devv) ;
  116. dex1 = dex1 et (prog depx1);
  117. i = i+1;
  118.  
  119. fin bou1;
  120.  
  121. sixx = evol manu 'depl_xx' dex1 'sigxx' sigxx1;
  122. sivv = evol manu 'defo_vol' dev 'sigxx' sigxx1 ;
  123.  
  124. *dess sivv;
  125. *dess sixx;
  126.  
  127. fin ;
  128.  
  129.  
  130.  
  131.  
  132.  

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