Télécharger la plaquette Cast3M


fr en

Télécharger rhmc_bcn.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  rhmc_bcn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * VERTICAL MOVEMENT OF A PILE
  5. * TEST: Rounded Hyperbolic Mohr-Coulomb model
  6. *
  7. *
  8. *-------------  OPCIONES GENERALES  --------------------------------
  9. *
  10. GRAPH = 'N' ;
  11. OPTION DIME 2 TRAC x ELEM qua4 MODE axis;
  12. *
  13. *-------------  CREACION DE LA GEOMETRIA  -------------------------
  14. *
  15. p1 = 7.5 0.;
  16. p2 = 345. 0.;
  17. p3 = 345. 5.;
  18. p4 = 7.5  5.;
  19. npoi1=1;
  20. l23 = d npoi1 p2 p3;
  21. l41 = d npoi1 p4 p1;
  22. xkk = 0.1*1.;
  23. ykk = 40.*1.;
  24. *xkk = 10.*1.;
  25. *ykk = 100.*1.;
  26. l12 = d p1 p2 dini xkk dfin ykk;
  27. l34 = d p3 p4 dini ykk dfin xkk;
  28. l43 = inve l34;
  29. malla = daller l12 l23 l34 l41;
  30. *
  31. *--------------  MODELO  ------------------------------------------
  32. *
  33.  E_ELAS = 1040.;
  34.  NUHOR  = 0.3;
  35.  COHES  = 1.;
  36. * 30 grados
  37.  ANGFR  = 0.52359878D0;
  38.  MOD1   = modeli malla mecanique elastique plastique rh_coulomb;
  39.  MAT1   =  mater mod1 youn E_ELAS nu NUHOR cohe COHES  phi ANGFR;
  40. *
  41. *--------------  MATRICES DE RIGIDEZ  -----------------------------
  42. *
  43.  rig1 = rigi mod1 mat1;
  44. *
  45. *--------------  CONDICIONES DE CONTORNO  -------------------------
  46. *
  47.  rigcont1 = (bloq uz l23) et (bloq ur (l23 et l41));
  48.  rig1b    = rig1 et rigcont1;
  49. * 
  50.  rigco2a  = rela ur l12 - ur l43;
  51.  rigco2b  = rela uz l12 - uz l43;
  52.  rigcont2 = rigco2a et rigco2b;
  53.  rigcont  = rigcont1 et rigcont2;
  54.  rig1     = rig1b et rigcont2; 
  55. *
  56. *--------------  MOVIMIENTO IMPUESTO  -----------------------------
  57. *
  58.  rigsupe = (bloq uz l41);
  59.  valor   = 0.2D0;
  60.  movyy   = (depi RIGsupe valor);
  61.  evol1   = EVOL MANU t (PROG 0. 1.) level (PROG 0. 1.);
  62.  evfut   = CHAR DIMP evol1 movyy;
  63. *
  64. *--------------  RESOLUCION  --------------------------------------
  65. *
  66.  maxiter = 15;
  67.  t0 = table;
  68.  t0.MODELE              = mod1;
  69.  t0.CARACTERISTIQUES    = mat1;
  70.  t0.BLOCAGES_MECANIQUES = rigcont et rigsupe;
  71.  t0.CHARGEMENT          = evfut;
  72.  t0.ACCELERATION        = maxiter;  
  73.  t0.MAXITERATION        = maxiter;
  74.  t0.PRECISION           = 1.E-10;
  75.  t0.CONVERGENCE_FORCEE  = faux;  
  76.  t0.K_TANGENT           = vrai; 
  77.  t0.TEMPS_CALCULES      = prog 0. pas 0.3 0.3 pas 0.05 1.;
  78.  t0.TEMPS_SAUVES        = t0.TEMPS_CALCULES;
  79.  t0.HYPOTHESE_DEFORMATIONS = 'LINEAIRE';
  80.  pasapas t0;
  81. *
  82. *--------------  POSTPROCESO  -------------------------------------
  83. *
  84.  imax   = (dime t0.TEMPS_SAUVES) - 1;
  85.  fuer   = prog 0.;
  86.  i      = 1;
  87.  repeter blocdefi imax;
  88.    aux  = reac rigsupe t0.deplacements.i;
  89.    aux  = resu aux;
  90.    pbas = extr aux MAIL ; 
  91.    pbas = POIN 1 pbas ;
  92.    aux  = extr aux fz pbas;
  93.    fuer = fuer et (prog aux);
  94.    i=i+1;
  95.  fin blocdefi;
  96.  law = evol manu t t0.TEMPS_SAUVES f(t) fuer;
  97.  SI (NEG GRAPH 'N') ;
  98.    dessin law;
  99.  FINSI ;
  100.  err = (aux - 4.7193694E+02)/4.7193694E+02 ;
  101.  list err;
  102.  SI ((ABS err) < 0.5e-5) ;
  103.    ERRE 0 ;
  104.  SINON ;
  105.    ERRE 5 ;
  106.  FINSI ;
  107.  FIN;
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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