Télécharger la plaquette Cast3M


fr en

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.       'NON_LOCAL' 'MOYE' 'V_MOYENNE' (MOTS 'EPTI') ;
  41.  
  42.  
  43. mat1 = MATE mod1 'YOUN' 20.74e9 'NU' 0.2 'AC' 2.75
  44. 'AT' 1. 'BC' 2470. 'BT' 15000. 'EPD0' 10.e-5
  45. 'RC' 120.e6 'RT' 10.e6 'P' .4 'GAMA' 0.99 'ALFA' 0.5
  46. 'A' -1 'AH' 7.e-5 'BH' 2.e-2 'CH' 3.e-6 'BETA' 1. 'K0' 0.1
  47. 'LCAR' 0.05 ;
  48.  
  49. *** Conditions aux limites
  50.  
  51. cl1 = BLOQ P1 UZ UY;
  52. cl2 = BLOQ L1 UZ ;
  53. cl3 = BLOQ COTEG UX ;
  54. cldepl = BLOQ COTED UX ;
  55.  
  56. cltot = cl1 et cl2 et cl3 et cldepl ;
  57.  
  58.  
  59. *** Chargement
  60.  
  61. fo1 = DEPI cldepl 1. ;
  62. prog1 = PROG 0. 1. ;
  63. prog2 = PROG 0. -2.e-3 ;
  64.  
  65. step1 = 1./20. ;
  66. list1 = PROG 0. PAS step1 1. ;
  67.  
  68. evol1 = EVOL MANU prog1 prog2 ;
  69.  
  70. char1 = CHARGEMENT fo1 evol1 DIMP ;
  71.  
  72.  
  73. tab1 = TABLE ;
  74. tab1 . PLASTIQUE = vrai ;
  75. tab1 . MODELE = mod1 ;
  76. tab1 . CARACTERISTIQUES = mat1 ;
  77. tab1 . CHARGEMENT = char1 ;
  78. tab1 . BLOCAGES_MECANIQUES = cltot ;
  79. tab1 . TEMPS_CALCULES = list1 ;
  80. tab1 . TEMPS_SAUVES = list1 ;
  81. tab1 . MOVA = 'D' ;
  82. tab1 . PRECISION = 1.e-8;
  83. **tab1 . 'NON_LOCAL' = 'MOT' 'MOYE' ;
  84. **tab1 . 'LCAR' = 0.05 ;
  85. tab1 . 'CONN' = 'CONN' mod1 mat1 'NORMAL' 'INTERIEUR' ;
  86. *Tab1 . 'LVIAM' = mots EPTI ;
  87.  
  88. PASAPAS tab1 ;
  89.  
  90. ind = inde (tab1.contraintes);
  91. nbind1 = dime ind;
  92. nbind1 = nbind1 - 1 ;
  93. ************************************************************
  94.  
  95. *definition des fonctions *
  96.  
  97. sigxx1=prog;
  98. dex1=prog;
  99. dev = prog ;
  100.  
  101. i = 1;
  102.  
  103. repeter bou1 nbind1;
  104.  
  105.   cont1 = (tab1).contraintes.i;
  106.  
  107.   sigxx1=sigxx1 et (prog(extr cont1 SMXX 1 1 1));
  108.  
  109.   dep1 = (tab1) . deplacements . i ;
  110.  
  111.   eps0 = EPSI dep1 mod1;
  112.   depx1 = (EXTR eps0 EPXX 1 1 1);
  113.   depy1 = (EXTR eps0 EPYY 1 1 1);
  114.   depz1 = (EXTR eps0 EPZZ 1 1 1) ;
  115.  
  116. *  opti donn 5 ;
  117. ************************************************************
  118. * deformations suivant les 3 composantes *
  119.  
  120.   devv = (depx1 + depy1 + depz1) ;
  121.   dev = dev et (prog devv) ;
  122.   dex1 = dex1 et (prog depx1);
  123.  i = i+1;
  124.  
  125. fin bou1;
  126.  
  127. sixx = evol manu 'depl_xx' dex1 'sigxx' sigxx1;
  128. sivv = evol manu 'defo_vol' dev 'sigxx' sigxx1 ;
  129.  
  130. *dess sivv;
  131. *dess sixx;
  132.  
  133. fin ;
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  

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