* fichier : mrsl_bcn.dgibi
* TRIAXIAL TEST WITH A NON-HOMOGENEOUS SAMPLE
* TEST: MRS-Lade model
*
fin; probleme // :
*
*------------- OPCIONES GENERALES --------------------------------
*
GRAPH = 'N' ;
OPTION DIME 2 TRAC x ELEM qua4 MODE axis;
*
precomp = 600.D0;
valdesp = -0.05 * 1.5;
*
npoi = 10;
npoi2 = 10;
*
*------------- CREACION DE LA GEOMETRIA -------------------------
*
p1 = 0. 0.;
p2 = 1. 0.;
p3 = 1. 1.5;
p4 = 0. 1.5;
l12 = d npoi p1 p2;
l23 = d npoi2 p2 p3;
l34 = d npoi p3 p4;
l41 = d npoi2 p4 p1;
malla = daller l12 l23 l34 l41;
mall1 = malla elem 1;
da = 1. / (flot npoi);
db = 1.5 / (flot npoi2);
aa = da 0.;
bb = da db;
cc = 0. db;
la2 = d 1 aa p2;
l4c = d 1 p4 cc;
lcb = d 1 cc bb;
lba = d 1 bb aa;
contaux = la2 et l23 et l34 et l4c et lcb et lba;
mall2 = malla inclus contaux 'BARY';
*
*-------------- MODELO ------------------------------------------
*
ROHOR = 2.5D3;
E_ELAS = 146.D3;
NUHOR = 0.2D0;
PCC = 0.D0;
PAA = 1.D0;
QAA = 1.D0;
EXM = 7.423D-2;
EWW = 0.7D0;
XK1 = 0.2D0;
XK2 = 0.7256D0;
ETB = 1.12D0;
EXV = 1.15D0;
XEP = 7.5D-1;
XNN = -0.5D0;
CCO = 4.D-3;
EXL = 1.0867654D0;
PCO = 5000.D0;
EXR = 1.592D0;
CCA = 0.000159D0;
PHIs = 0.D0;
ALPas = 0.8D0;
MOD1b = modl MALL2 mecanique elastique plastique mrs_lade;
MAT1b = mater mod1b youn E_ELAS nu NUHOR rho ROHOR pc PCC pa PAA
qa QAA expm EXM e EWW k1 XK1 k2 XK2 etab ETB
expv EXV epsi XEP n XNN ccon CCO expl EXL
pcap PCO expr EXR ccap CCA phi PHIs alp ALPas;
*
etb2 = etb*0.9D0;
pco2 = pco*0.9D0;
MOD1a = modl MALL1 mecanique elastique plastique mrs_lade;
MAT1a = mater mod1a youn E_ELAS nu NUHOR rho ROHOR pc PCC pa PAA
qa QAA expm EXM e EWW k1 XK1 k2 XK2 etab ETB2
expv EXV epsi XEP n XNN ccon CCO expl EXL
pcap PCO2 expr EXR ccap CCA phi PHIs alp ALPas;
modi = mod1a et mod1b;
mati = mat1a et mat1b;
*
*-------------- PROBLEMA: PRECOMPRESION -------------------------
*
RIGconti = (bloq uz l12) et (bloq ur l41);
valorp = pres mass modi precomp (l23 et l34);
evolu0 = EVOL MANU t (PROG 0. 1.) level (PROG 0. 1.);
evfut0 = CHAR DIMP evolu0 valorp;
maxiter = 10;
ti = TABLE ;
ti.MODELE = MODi;
ti.CARACTERISTIQUES = MATi;
ti.BLOCAGES_MECANIQUES = RIGconti;
ti.CHARGEMENT = evfut0;
*ti.ACCELERATION = maxiter;
*ti.MAXITERATION = maxiter;
ti.PRECISION = 1.E-8;
*ti.CONVERGENCE_FORCEE = faux;
*ti.K_TANGENT = VRAI;
ti.TEMPS_CALCULES = PROG 0. 1.;
ti.TEMPS_SAUVES = PROG 0. 1.;
PASAPAS ti ;
desppre = ti.deplacements.1;
form desppre;
*
*-------------- PROBLEMA: TRIAXIAL ------------------------------
*
mod1 = mod1a et mod1b;
mat1 = mat1a et mat1b;
RIGcont = (bloq uz l12) et (bloq ur l41);
*
*-------------- MOVIMIENTO IMPUESTO -----------------------------
*
rigsupe = (bloq uz l34);
movyy = (depi rigsupe valdesp);
evolu2 = EVOL MANU t (PROG 0. 1.) level (PROG 0. 1.);
evfut2 = CHAR DIMP evolu2 movyy;
valorp = pres mass mod1 precomp l23;
evolu1 = EVOL MANU t (PROG 0. 1.) level (PROG 1. 1.);
evfut1 = CHAR DIMP evolu1 valorp;
*
*-------------- RESOLUCION --------------------------------------
*
tiempos = prog 0. pas 0.04 0.48 pas 0.02 1.;
maxiter = 15;
t0 = TABLE ;
t0.CONTRAINTES = TABLE ;
t0.CONTRAINTES.0 = ti.CONTRAINTES.1;
t0.MODELE = MOD1;
t0.CARACTERISTIQUES = MAT1;
t0.BLOCAGES_MECANIQUES = RIGcont et rigsupe;
t0.CHARGEMENT = evfut2 et evfut1;
*t0.ACCELERATION = maxiter;
*t0.MAXITERATION = maxiter;
t0.PRECISION = 1.E-10;
*t0.CONVERGENCE_FORCEE = faux;
t0.K_TANGENT = VRAI;
t0.TEMPS_CALCULES = tiempos;
t0.TEMPS_SAUVES = tiempos;
PASAPAS t0 ;
*
*-------------- POSTPROCESO -------------------------------------
*
imax = 38;
fuer = prog 0.;
i = 1;
repeter blocdefi imax;
aux = reac rigsupe t0.deplacements.i;
aux = resu aux;
pbas = extr aux MAIL ;
pbas = POIN 1 pbas ;
aux = extr aux fz pbas;
aux = -1.D0 * aux;
fuer = fuer et (prog aux);
i=i+1;
fin blocdefi;
law = evol manu t (-1.0*(tiempos*valdesp)) f(t) fuer;
SI (NEG GRAPH 'N') ;
dessin law;
FINSI ;
err = (aux - 3.14205E+03)/3.14205E+03 ;
list err;
SI ((ABS err) < 0.5e-5) ;
ERRE 0 ;
SINON ;
ERRE 5 ;
FINSI ;
FIN;
© Cast3M 2003 - All rights reserved.
Disclaimer