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

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