Télécharger fefp_powcap_bcn.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fefp_powcap_bcn.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * COMPACTION OF A FLANGED COMPONENT
  5. *
  6. *------------- OPCIONES GENERALES --------------------------------
  7. *
  8. OPTION DIME 2 TRAC x ELEM qua8 MODE axis;
  9. *
  10. *------------- CREACION DE LA GEOMETRIA --------------------------
  11. *
  12. r1 = 6.3D0; r2 = r1 + 4.6D0; r3 = r2 + 15.6D0;
  13. h1 = 0.D0; h2 = h1 + 13.7D0; h3 = h2 + 11.7D0;
  14. p1 = r1 h1;p2 = r2 h1;p3 = r2 h2;
  15. p4 = r3 h2;p5 = r3 h3;p6 = r2 h3;
  16. p7 = r1 h3;p8 = r1 h2;paux = r3 h3;
  17. den1 = 5; den2 = 10; den3 = 10; den4 = 8;
  18. l12 = d den1 p1 p2;l23 = d den2 p2 p3;
  19. l38 = d den1 p3 p8;l81 = d den2 p8 p1;
  20. l34 = d den3 p3 p4;l45 = d den4 p4 p5;
  21. l56 = d den3 p5 p6;l63 = d den4 p6 p3;
  22. l83 = d den1 p8 p3;l36 = d den4 p3 p6;
  23. l67 = d den1 p6 p7;l78 = d den4 p7 p8;
  24. malla1 = daller l12 l23 l38 l81;
  25. malla2 = daller l34 l45 l56 l63;
  26. malla3 = daller l83 l36 l67 l78;
  27. malla = malla1 et malla2 et malla3 ;
  28. elim malla 1.D-4;
  29. llbase = l12 et l34; llsupe = l56 et l67;
  30. lllate = l23 et l45 et l78 et l81;
  31. *
  32. *------------- MODELO -------------------------------------------
  33. *
  34. E_ELAS = 2000.D0;
  35. NUHOR = 0.37D0;
  36. ROHOR = 2.5D3;
  37. SIGINI = 90.D0;
  38. EXPN1 = 1.D0;
  39. EXPN2 = 2.7D0;
  40. ETAINI = 0.489D0;
  41. COHES = 15.D0;
  42. ANGFR = 0.52359878D0 * 2.;
  43. ANGFR0 = 0.52359878D0 * 1.5;
  44. NNC = 2.D0;
  45. mod1 = 'MODELI' malla mecanique elastique plastique
  46. powdercap_fefp Q8RI;
  47. mat1 = mater mod1 youn E_ELAS nu NUHOR rho ROHOR
  48. sigy SIGINI nnn1 EXPN1 nnn2 EXPN2
  49. eta0 ETAINI cohe COHES
  50. phi ANGFR phi0 ANGFR0 nnnc NNC ;
  51. *
  52. *-------------- CONDICIONES DE CONTORNO -------------------------
  53. *
  54. rigcont = (bloq uz llbase) et (bloq ur lllate);
  55. *
  56. *-------------- MOVIMIENTO IMPUESTO -----------------------------
  57. *
  58. rigsupe = bloq uz llsupe;
  59. valor = -6.06D0;
  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 = 50;
  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.NITER_KTANGENT = 1;
  78. t0.PRECISINTER = 1.E-9;
  79. t0.NITERINTER_MAX = 100;
  80. t0.FEFP_FORMULATION = vrai;
  81. t0.UPDATE_LAGRANGIAN = vrai;
  82. t0.TEMPS_CALCULES = prog 0. pas 0.0001 0.001
  83. pas 0.0005 0.01
  84. pas 0.005 0.1
  85. pas 0.010 0.2;
  86. t0.TEMPS_SAUVES = t0.TEMPS_CALCULES;
  87. t0.'PROCESSEURS'='MONO_PROCESSEUR';
  88. pasapas t0;
  89. *
  90. *--------------- POSTPROCESO ------------------------------------
  91. *
  92. imax = (dime t0.TEMPS_CALCULES) - 1;
  93. fuer = prog 0.;
  94. i = 1;
  95. geo0 = form;
  96. repeter blocdefi imax;
  97. aux = reac rigsupe t0.deplacements.i;
  98. aux = resu aux;
  99. pbas = extr aux MAIL ;
  100. pbas = POIN 1 pbas ;
  101. aux = extr aux fz pbas;
  102. fuer = fuer et (prog aux);
  103. i=i+1;
  104. fin blocdefi;
  105. des = valor * t0.TEMPS_CALCULES;
  106. law = evol manu t des f(t) fuer;
  107. err = (aux + 2.1063258E+04)/2.1063258E+04 ;
  108. list err;
  109. SI ((ABS err) < 0.5e-5) ;
  110. ERRE 0 ;
  111. SINON ;
  112. ERRE 5 ;
  113. FINSI ;
  114. FIN;
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  

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