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

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