Télécharger compression_nloc.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : compression_nloc.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. tab1 . 'NON_LOCAL' = 'MOT' 'MOYE' ;
  82. tab1 . 'LCAR' = 0.05 ;
  83. tab1 . 'CONN' = 'CONN' mod1 0.05 'NORMAL' 'INTERIEUR' ;
  84. tab1 . 'LVIAM' = mots EPTI ;
  85.  
  86. PASAPAS tab1 ;
  87.  
  88. ind = inde (tab1.contraintes);
  89. nbind1 = dime ind;
  90. nbind1 = nbind1 - 1 ;
  91. ************************************************************
  92.  
  93. *definition des fonctions *
  94.  
  95. sigxx1=prog;
  96. dex1=prog;
  97. dev = prog ;
  98.  
  99. i = 1;
  100.  
  101. repeter bou1 nbind1;
  102.  
  103. cont1 = (tab1).contraintes.i;
  104.  
  105. sigxx1=sigxx1 et (prog(extr cont1 SMXX 1 1 1));
  106.  
  107. dep1 = (tab1) . deplacements . i ;
  108.  
  109. eps0 = EPSI dep1 mod1;
  110. depx1 = (EXTR eps0 EPXX 1 1 1);
  111. depy1 = (EXTR eps0 EPYY 1 1 1);
  112. depz1 = (EXTR eps0 EPZZ 1 1 1) ;
  113.  
  114. * opti donn 5 ;
  115. ************************************************************
  116. * deformations suivant les 3 composantes *
  117.  
  118. devv = (depx1 + depy1 + depz1) ;
  119. dev = dev et (prog devv) ;
  120. dex1 = dex1 et (prog depx1);
  121. i = i+1;
  122.  
  123. fin bou1;
  124.  
  125. sixx = evol manu 'depl_xx' dex1 'sigxx' sigxx1;
  126. sivv = evol manu 'defo_vol' dev 'sigxx' sigxx1 ;
  127.  
  128. *dess sivv;
  129. *dess sixx;
  130.  
  131. fin ;
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  

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