Télécharger fefp_powder_bcn.dgibi
* fichier : fefp_powder_bcn.dgibi ************************************************************************ ************************************************************************ * COMPACTION OF A FLANGED COMPONENT * *------------- OPCIONES GENERALES -------------------------------- * * *------------- CREACION DE LA GEOMETRIA -------------------------- * r1 = 6.3D0; r2 = r1 + 4.6D0; r3 = r2 + 15.6D0; h1 = 0.D0; h2 = h1 + 13.7D0; h3 = h2 + 11.7D0; p1 = r1 h1;p2 = r2 h1;p3 = r2 h2; p4 = r3 h2;p5 = r3 h3;p6 = r2 h3; p7 = r1 h3;p8 = r1 h2;paux = r3 h3; den1 = 5; den2 = 10; den3 = 10; den4 = 8; l12 = d den1 p1 p2;l23 = d den2 p2 p3; l38 = d den1 p3 p8;l81 = d den2 p8 p1; l34 = d den3 p3 p4;l45 = d den4 p4 p5; l56 = d den3 p5 p6;l63 = d den4 p6 p3; l83 = d den1 p8 p3;l36 = d den4 p3 p6; l67 = d den1 p6 p7;l78 = d den4 p7 p8; malla1 = daller l12 l23 l38 l81; malla2 = daller l34 l45 l56 l63; malla3 = daller l83 l36 l67 l78; malla = malla1 et malla2 et malla3 ; llbase = l12 et l34; llsupe = l56 et l67; lllate = l23 et l45 et l78 et l81; * *------------- MODELO ------------------------------------------- * E_ELAS = 2000.D0; NUHOR = 0.37D0; ROHOR = 2.5D3; SIGINI = 90.D0; EXPN1 = 1.D0; EXPN2 = 2.7D0; ETAINI = 0.489D0; mod1 = 'MODELI' malla mecanique elastique plastique powder_fefp Q8RI; mat1 = mater mod1 youn E_ELAS nu NUHOR rho ROHOR sigy SIGINI nnn1 EXPN1 nnn2 EXPN2 eta0 ETAINI; * *-------------- CONDICIONES DE CONTORNO ------------------------- * * *-------------- MOVIMIENTO IMPUESTO ----------------------------- * valor = -6.06D0; * *--------------- RESOLUCION ------------------------------------- * maxiter = 50; t0 = table; t0.MODELE = mod1; t0.CARACTERISTIQUES = mat1; t0.BLOCAGES_MECANIQUES = rigcont et rigsupe; t0.CHARGEMENT = evfut; t0.ACCELERATION = maxiter; t0.MAXITERATION = maxiter; t0.PRECISION = 1.E-10; t0.CONVERGENCE_FORCEE = faux; t0.K_TANGENT = vrai; t0.NITER_KTANGENT = 1; t0.PRECISINTER = 1.E-9; t0.NITERINTER_MAX = 100; t0.FEFP_FORMULATION = vrai; t0.UPDATE_LAGRANGIAN = vrai; pas 0.001 0.01 pas 0.005 0.1 pas 0.010 0.2; t0.TEMPS_SAUVES = t0.TEMPS_CALCULES; t0.'PROCESSEURS' = 'MONO_PROCESSEUR'; pasapas t0; * *--------------- POSTPROCESO ------------------------------------ * i = 1; repeter blocdefi imax; i=i+1; fin blocdefi; des = valor * t0.TEMPS_CALCULES; err = (aux + 1.9003465E+04)/1.9003465E+04 ; list err; SI ((ABS err) < 0.5e-5) ; SINON ; FINSI ; FIN;
© Cast3M 2003 - Tous droits réservés.
Mentions légales