Télécharger sste1_bcn.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : sste1_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.SUBSTEPPING = vrai;
  80. t0.NMAXSUBSTEPS = 1000;
  81. t0.PRECISINTER = 1.E-12;
  82. t0.NITERINTER_MAX = 5;
  83. t0.TEMPS_CALCULES = prog 0. pas 0.3 0.3 pas 0.05 1.;
  84. t0.TEMPS_SAUVES = t0.TEMPS_CALCULES;
  85. pasapas t0;
  86. *
  87. *-------------- POSTPROCESO -------------------------------------
  88. *
  89. imax = (dime t0.TEMPS_SAUVES) - 1;
  90. fuer = prog 0.;
  91. i = 1;
  92. repeter blocdefi imax;
  93. aux = reac rigsupe t0.deplacements.i;
  94. aux = resu aux;
  95. pbas = extr aux MAIL ;
  96. pbas = POIN 1 pbas ;
  97. aux = extr aux fz pbas;
  98. fuer = fuer et (prog aux);
  99. i=i+1;
  100. fin blocdefi;
  101. law = evol manu t t0.TEMPS_SAUVES f(t) fuer;
  102. SI (NEG GRAPH 'N') ;
  103. dessin law;
  104. FINSI ;
  105. err = (aux - 4.71961E+02)/4.71961E+02 ;
  106. list (aux - 4.71961E+02);
  107. list err;
  108. SI ((ABS err) < 0.5e-5) ;
  109. ERRE 0 ;
  110. SINON ;
  111. ERRE 5 ;
  112. FINSI ;
  113. FIN;
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  

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