Télécharger la plaquette Cast3M


fr en

Télécharger j2_bcn.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  j2_bcn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * PERFORATED STRIP UNDER TRACTION
  5. * TEST: J2 perfect plasticity model
  6. *
  7. *
  8. *-------------  OPCIONES GENERALES  --------------------------------
  9. *
  10. GRAPH = 'N' ;
  11. OPTION DIME 2 TRAC x ELEM qua8 MODE plan defo;
  12. *
  13. *-------------  CREACION DE LA GEOMETRIA  --------------------------
  14. *
  15. p1 =  5.  0.;
  16. p2 = 10.  0.;
  17. p3 = 10. 18.;
  18. p4 =  0. 18.;
  19. p5 =  0.  5.;
  20. pce = 0.  0.;
  21. aux = 0.7071067 * 5.;
  22. p6 =  aux aux;
  23. den = 12;
  24. l12 = d den p1 p2;
  25. l23 = d den p2 p3;
  26. l36 = d den p3 p6;
  27. l61 = cerc den p6 pce p1;
  28. l34 = d den p3 p4;
  29. l45 = d den p4 p5;
  30. l56 = cerc den p5 pce p6;
  31. l63 = d den p6 p3;
  32. mall1 = daller l12 l23 l36 l61;
  33. mall2 = daller l34 l45 l56 l63;
  34. malla = mall1 et mall2;
  35. elim 1.E-3 malla;
  36. *
  37. *-------------  MODELO  -------------------------------------------
  38. *
  39.  NUHOR  = 0.2;
  40.  ROHOR  = 2.5e3;
  41.  E_ELAS = 7.e7;
  42.  E_PLAS = 0.e0;
  43.  K_PLAS = (E_PLAS*E_ELAS)/(E_ELAS-E_PLAS);
  44.  SININI = 2.43e5;
  45.  SINFIN = 0.;
  46.  VELOCI = 0.;
  47.  mod1   = model malla mecanique elastique plastique j2;
  48.  mat1   = mater mod1 youn E_ELAS  nu   NUHOR    rho  ROHOR
  49.                      sig0 SININI  sigi SINFIN   kiso K_PLAS
  50.                      velo VELOCI;
  51. * mod1   = model malla mecanique elastique plastique parfait;
  52. * mat1   = mater mod1 youn E_ELAS  nu   NUHOR    
  53. *                     sigy SININI  ;
  54. *
  55. *--------------  MATRICES DE RIGIDEZ  -----------------------------
  56. *
  57.  rig1 = rigi mod1 mat1;
  58. *
  59. *--------------  CONDICIONES DE CONTORNO  -------------------------
  60. *
  61.  rigcont1 = bloq uy l12;
  62.  rigcont2 = bloq ux l45;
  63.  rigcont  = rigcont1 et rigcont2;
  64.  rig2     = rig1 et rigcont;
  65. *
  66. *--------------  MOVIMIENTO IMPUESTO  -----------------------------
  67. *
  68.  rigsupe = bloq uy l34;
  69.  valor   = .1;
  70.  movyy   = depi rigsupe valor;
  71.  evol1   = evol manu t (prog 0. 1.) level (prog 0. 1.);
  72.  evfut   = char dimp evol1 movyy;
  73. *
  74. *---------------  RESOLUCION  -------------------------------------
  75. *
  76.  maxiter = 15;
  77.  t0 = table;
  78.  t0.MODELE              = mod1;
  79.  t0.CARACTERISTIQUES    = mat1;
  80.  t0.BLOCAGES_MECANIQUES = rigcont et rigsupe;
  81.  t0.CHARGEMENT          = evfut;
  82.  t0.ACCELERATION        = maxiter;  
  83.  t0.MAXITERATION        = maxiter;
  84.  t0.PRECISION           = 1.E-10;
  85.  t0.CONVERGENCE_FORCEE  = faux;  
  86.  t0.K_TANGENT           = vrai;
  87.  t0.TEMPS_CALCULES      = prog 0. pas 0.1  1.;
  88.  t0.TEMPS_SAUVES        = t0.TEMPS_CALCULES;
  89.  t0.HYPOTHESE_DEFORMATIONS = 'LINEAIRE' ;
  90.  pasapas t0;
  91. *
  92. *---------------  POSTPROCESO  ------------------------------------
  93. *
  94.  imax   = 10;
  95.  fuer   = prog 0.;
  96.  i      = 1;
  97.  repeter blocdefi imax;
  98.    aux  = reac rigsupe t0.deplacements.i;
  99.    aux  = resu aux;
  100.    pbas = extr aux MAIL ; 
  101.    pbas = POIN 1 pbas ;
  102.    aux  = extr aux fy pbas;
  103.    fuer = fuer et (prog aux);
  104.    i=i+1;
  105.  fin blocdefi;
  106.  law = evol manu t (prog 0. pas 0.02 0.2) f(t) fuer;
  107.  SI (NEG GRAPH 'N') ;
  108.    dessin law;
  109.   @cartoon t0 malla 1.;
  110.  FINSI ;
  111.  err = (aux - 1.40382E+06)/1.40382E+06 ;
  112.  list err;
  113.  SI ((ABS err) < 0.5e-5) ;
  114.    ERRE 0 ;
  115.  SINON ;
  116.    ERRE 5 ;
  117.  FINSI ;
  118.  FIN;
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  

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