Ensayos Presiometricos
Ensayos Presiometricos
Ensayos Presiometricos
INGENIEROS DE MINAS
- -
RAMON RODR~GUEZPONS-ESPARVER
Ingeniero de Minas
DIRECTOR
INGENIEROS D E MINAS
--
--- -
(Firma)
PRESIDENTE:
. - 5 p,A L P; T P R ~ H
VOCALES:
. An'll%~-1\510 YAGbp5
./&,,S ~u'c+zet-l;
VOCAL SECRETARIO:
. FiQiVUFL ?fi%ak
&&E?
CALIFICACI~N:
~ P T O C d N (-AL, ()e ( P ~ RL)IJA~~)YI.DAD).
EL PRESIDENTE LOS VOCALES
Este trabajo de tesis con1enz. gracias a una beca de Formacin de Personal Inkestigador
del Ministerio de Educacin y Ciencia. Con ella me fue posible realizar varias estancias
breves en las instalaciones que el Ismes, Spa tiene en Brgamo (1). De la colaboracion
entre el lsmes y el Departamento de Matemtica Aplicada y Mtodos lnt'ormticos de la
E.T.S.I. de Minas de Madrid y de la participacin conjunta en el proyecto europeo
Interclay II surgi la posibilidad de desarrollar el trabajo que a continuacin se presenta.
A lo largo de estos aos han sido muchas las personas que de una u otra forma me han
apoyado y acompaado. A todas ellas mi ms sincero agradecimiento.
Quiero agradecer a Dfia. Rta Pellegrini. ingeniero del Ismes, su disponibilidad y apoyo
durante mi permanencia en Brgamo.
1-0s datos que se emplean para llevar a cabo y validar esta metodologa proceden de un
conjunto de ensayos presiomtricos (SBP) drenados realizados en celda de calibracin
con arenas del Ticino [Belloti et al. 19871.
Con todo ello se ha desarrollado un cdigo de ordenador que incluye diferentes tcnicas
de estimacin autonitica de parmetros geotcnicos a partir del ensayo presiometrico.
Este cdigo puede ser fcilmente ampliado para la determinacin de otros parmetros
correspondientesa otros ensayos geotcnicos o geomecnicos.
ABSTRACT
The work presented in this research proposes a methodology for the automatic
interpretation of data provided by the pressuremeter test that can be used in an eficient
way in the geotechnical design. First this methodology requires mathematical models
capable of reproducing the response observed in the experiment in the considered soil -
sandy soils in our case-. and second, identification methods of material characteristic
parameters from experimental data.
Different parameter identification methods have been applied to two material behaviour
models of the pressuremeter response in sandy soils. One OSthem is an analytic non
linear elastic model [Femeira, 19921 and the other, a classical elastoplastic model
[Drucker-Prager, 19521.
The data employed to develop and validate the methodology carne from a set of self-
boring pressuremeter drained tests (SBP). c d e d out in a calibration chamber with
Xcino sand, reported by Belloti et al, 1987.
Hay que sefalar que los ensayos de laboratorio presentan ciertas limitaciones. Uno de
sus inconvenientes es el tamao de las muesuas. Normalmente son pequefas y por tanto
incapaces de reflejar :a respuesta de una gran niasa de material, en la que suele haber
fsuras. cementaciones. heterogeneidades. Es lo que suele generar el llamado efecto de
escala. Por este motivo se tiende a estudiar muestras del mayor tamaiio posible. Por otro
lado es muy dificil que la muestra llegue inalterada al laboratorio ya que suele sufrir los
procesos de extraccin. transporte. almacenamiento ... que ocasionan la prdida del
estado de tensiones, ndice de poros, humedad ...originales.
Los ensayos de campo o ensayos sobre el terreno presentan otro tipo de problemas. En
ellos es ms dificil definir y conocer las condiciones de contorno, las tensiones y las
deformaciones no suelen variar espacialmente de forma uniforme y las condiciones de
drenaje suelen no ser conocidas y varian espacial y temporalmente. Adems el terreno
facilmente sufre alteraciones al aplicar el insuumento de medida. Por otro lado es cierto
que en general los ensayos sobre el terreno recogen y reflejan el efecto de la gran escala,
el efecto Quetiene en cuenta las heterweneidades o defecto de uniformidad, que suelen
ser los que en ltima instancia defmen el comportamiento del medio a escala de obra
civil o minera.
Desde que L. Menard present en 1957 el primer presimetro para realizar ensayos en
campo, se han desarrollado diferentes prototipos con el fin de poder estudiar
adecuadamente el comportamiento de distintos medios -suelos arenosos o arcillosos y
rocas blandas- sometidos a condiciones distintas. En su conjunto los ensayos
presiomtricos se clasifican en cuatro grupos: presimetro de Menar o presimetro con
pre-perforacin (Pre-Bored Pressuremeter. PBP). presimetro autoperforante (Self-
Boring I'ressuremeter. SBP). presimetro de desplazamiento total (Full-Displacement
I'ressurenieter. FDP) y presiomeiro de empuje (Push-In Pressuremetcr. PIP). La
diferencia principal entre unos y otros radica en el modo como se introduce el
presimetro en el medio que se desea estudiar. Inicialmente esta clasificacin se reducia
a dos tipos de presimetro: el presimetro con pre-perforacin y el presinietro
autoperforante. El primero se lleva a cabo al introducir la sonda en una perforacin
previa. mientras que el segundo realiza l mismo la perforacin a medida que se va
introduciendo en el terreno con lo que se reducen las alteraciones del medio en los
alrededores de la pared interior del sondeo. Esta clasificacin se ampli al desarrollarse
el presimetro de desplazamientototal y el presimetro de empuje, en los que el ensayo
se lleva a cabo introduciendo el presimetro por empuje en el terreno.
I'or otro lado. con la llegada de las celdas de calibracin se ha abierto la posibilidad de
realizar ensayos presiomtricos en laboratorio. Este hecho ha permitido estudiar y
iiicjorar las tcnicas de insercin de la sonda y de toma de datos. concebir nuevos
ensayos y calibrar los distintos mtodos de interpretacin existentes que permiten
deducir los parmetros de comportamiento del medio.
Una vez definido el modelo de comportamiento, el siguiente paso ser detemiinar los
parimetros del mismo mediante tecnicas manuales o automticas. A este proceso se le
denomina retroanlisis o problema inverso. La denominacin de problema inverso surge
en contraposicin al problema directo. En ste ultimo, a partir del conocimiento de un
niodelo de comportamiento, de unos parmetros y de las acciones que actan sobre el
sistema se busca simular la respuesta del ensayo. Por el contrario, el problema inverso
consiste en estimar los parmetros del modelo de comportamiento a partir de medidas
sobre la respuesta del sistema y de la informacin previa disponible sobre los distintos
parmetros.
Con todo ello se desarroll un cdigo de ordenador que incluye diferentes tcnicas de
estimacin automtica no lineal de parmetros en suelos a partir del ensayo
presiom6trico que pueden ser aplicadas para la determinacin de parmetros
correspondientes a otros ensayos geotcnicos o geomecnicos, con slo cambiu dentro
del programa, el mdulo que corresponde al modelo del material utilizado. La
metodologa que se propone es una metodologa muy verstil qr incorpora tcnicas de
identificacin automtica de parmetros novedosas en geotecnia. Estas tcnicas han sido
verificadas por un conjunto de ensayos sintticos. haciendo tambin un anlisis de
incertidumbre de los pariimetros obtenidos. Posteriormente esta metodologa se ha
aplicado con xito a la interpretacin de ensayos reales.
Los datos que se emplean para llevar a cabo y validar esta metodologa proceden de un
conjunto de ensayos presiomtricos (SBP) drenados realizados en celda de calibracin
con arenas del Ticino [Belloii et al. 19871.
Esta memoria de tesis se desarrolla en seis captulos y tres anexos. En el primer capitulo
se describe sintticamente en qu consiste un ensayo presiomtrico, se hace una breve
clasificacin histrica de los diferentes mtodos de interpretacin y por ltimo. se hace
referencia de un modo ms exhaustivo y detallado a la metodologa de identificacin de
parmetros desarrollada, asi como a las etapas y objetivos de la tesis. En el Anexo 1 se
puede encontrar una breve descripcin de los diferentes tipos de presimetros que
existen junto a una revisin histrica de los mtodos de interpretacin ms utilizados
internacionalmente con el fin de exponer cul es el estado del arte en la realizacin e
interpretacin de ensayos presiomtricos.
Por ltimo. en el captulo sexto se presentan las principales conclusiones del trabajo. as
como las lneas futuras de investigacin derivadas de esta tesis.
PLANTEAMIENTOS BASICOS PARA EL
DESARROLLO DE UNA METODOLOGA DE
IDENTIFICACIN DE PARMETROS
Cuando a finales de los aos 50 apareci el primer presimetro ste haba sido diseiado
como una herramienta para analizar formaciones rocosas compactos y suelos in-situ.
Uno de los primeros presimetros fue el presimetro de Menard. Este presimetro
forma parte de la clase de presibmetros con pre-perforacin (pre-bored pressuremeter.
PBP) que producen una distorsin considerable en el medio. A pesar del cambio
significativo del estado de tensiones en el medio provocado por la perforacin. la
tcnica del presimetro con pre-perforacin puede resultar adecuada en suelos rigidos o
compactos y en rocas dbiles. Los resultados del ensayo se presentan en diagrama
presin-deformacin y muestran una curva caracterstica de carga denominada curva-s
(Figura Al-3. Anexo 1). La distorsin causada por la perforacin es, de alguna manera.
similar a la que genera la preparacin de Izs muestras de los ensayos de laboratorio. En
este sentido no se aportaba ninguna mejora. Precisamente por este motivo se dio poca
importancia a esta ttcnica durante los d o s 60. salvo en Francia donde s se continu
aplicando y desarrollando.
En esta breve pero intensa historia del desarrollo de los ensayos presiom6uicos. se
pueden identificar tres reas principales de investigacin:
- La segunda tiene que ver con el diseo del ensayo: ensayos con control de
deformaciones, control de tensiones, holding test (Pymh et al, 1985; Belloti
et al, 1986; Fioravante, 19881.
Este captulo comienza haciendo una breve descripcin del ensayo presiomtrico
indicando sus principales caractersticas y la informacin sobre el comportamiento del
terreno que se puede obtener con dl -poniendo especial atencin en los suelos arenosos
pues sobre ellos se han llevado a cabo los ensayos con los que se ha trabajado en esta
tesis-. A continuacin se hace referencia de un modo sinttico a un3 clasificacion de los
mtodos de interpretacin del ensayo presiomtnco presentes en la literatura tecnolgica
internacional. Por ultimo. se presenta la metodologa de identificaci6n de parmetros
desarrollada en esta tesis indicando las causas que la motivan y los objetivos que
persigue.
Los factores que ejercen una influencia sobre la resistencia al cone en medios arenosos.
para un mismo tipo de arena son los siguientes: densidad relativa (D,), tensin de
confinamiento en la celda (q). anisotropia del medio. tipo de ensayo e historia de carga
del suelo. De estos, los dos factores que ms influyen en el valor del Ungulo de
rozamiento interno de pico -parmetro de resistencia ms empleado en ingeniera- son 13
densidad relativa y la tensin de confinamiento.
DEFORMACION CORTANTE
Figura 1.2. Vnriacin del ngulo %ante el avance del proceso cortante eti
un elemenro de arena densa.
CAPITULO 1: PLANEAMIENTOS BASICOS PARA EL DESARROLLO DE UNA
METODOLOG~ADE IDENTIF~CACI~N DE PARAMETROS
.
. de defo-
1.2.2 C~ractcristieas .. de las arenas
Al igual que ocurre con la resistencia al wrte de las arenas. su deformabilidad depende
tambin de un conjunto de parmetros. Para un mismo tipo de arena las caractersticas
de su deformabilidad dependen bitsicamente de los siguientes factores: ndice de poros
inicial (e,,) n densidad relativa inicial (DR),
tensin de consolidacin de referencia (a,).
Iiisioria del estado tensio-deformativo, ngulo formado por la direccin de la tension
principal mxima apjicada sobre la muestra en el pasado y la direccin de la tension
principal mxima durante la fase desviadora, y tiempo transcurrido antes de la
apliccii>nde las tensiones desviadoras.
Cuando se trabaja con deformaciones cortantes (y) inferiores al 10' % se puede decir
[Ualdi et 1. 1986) que el mdulo cortante inicial depende nicamente de la densidad
relativa y de la tensin de consolidacin de referencia (G,=f(DR.o,)).Este mdulo se
puede considerar realmente como un nldulo elstico de las arenas ya que estas para
'
dek>rmacionesy 10" % y a paridad de DRy o,presentan un comportamiento elstico-
lineal. La influencia del resto de los pahnetros en el mdulo elitstico cortante para
deformaciones y > 1 0 ' ~% puede encontrarse en [Jambu. 19631 y en [Konder y Zelasko.
19633.
Por el tipo de ley constitutiva empleada se pueden distinguir entre: modelos elsticos
lineales [Baguelin et al, 1978; Mair y Wood.19871. rgido plsticos (Hughes. 19771
modelos elstico lineal perfectamente plsticos [Gibson y Anderson. 1961; Vesic. 1972;
Jefferies, 1988; Houlsby, 1986; Houlsby y Withen. 19881, modelos elsticos
hiperblicos [Ferreira, 19921, modelos no lineales perfectamente plsticos [Denby.
1978, Bahar. 19921, modelos elastoplsticos con endurecimiento o reblandecimiento
[Prevost y Hoeg, 1975; Juran y Mahmoodzadegan, 19891.
Por el tipo de material sobre el que se realiza el ensayo se pueden distinguir entre:
materiales de comportamiento cohesivo puro -arcilla-. materiales de comportamiento de
friccin pura -arena- y materiales de comportamiento cohesivo-fnccional -suelo natural-
Por el tipo de informacin del ensayo presiomtrico que se utiliza se pueden distinguir
aquellos mtodos que requieren nicamente los datos de la fase de carga. aquellos que
requieren nicamente los datos de la fase de descarga y aquellos otros que precisan de la
informacibn tanto de la curva de carga como de la de descarga. En los ensayos SBP la
curva de descarga puede considerarse que tiene la misma importancia que la
correspondiente a la fase de carga.
Desde un punto de vista matemtico se pueden distinguir entre: mtodos que emplean
soluciones en forma analitica [Baguelin. 1972; Palmer. 1972; Ladanyi. 19721 realizando
o no un ajuste con la cuna experimental. mtodos que a partir de los datos
experimentales y mediante algn tipo de construccin geomerrica obtienen los
parimetros [Wroth y Windle. 19751 y mtodos que se basan en soluciones numricas
aproximadas empleando tcnicas de diferencias finitas [Ladanyi. 1963; Manassero;
19891. elementos finitos.. .
Algunos de los mdtodos de interpretacin del ensayo presiomtrico que han surgido en
las iiltimas dcadas se presentan en al Anexo 1. Cada uno de ellos tiene un campo de
aplicacibn concreto. emplea diferentes modelos de comportamiento dependiendo de las
hiptesis de panida. presenta sus limitaciones y tienen sus ventajas e inconvenientes
sobre los dems.
Juego de
parametros iniciales
Nuevo juego de
parhetros (Algontmo
de minimizacin) Informacin
V adtc~onalsobre
A Clculo de la funcin * los parmetros
objetivo
4
Datos
experimentales
C. Convergencia
Fin
Un ensayo presiomtrico es una de las pruebas in-situ que mejor permiten ohtener
informacin acerca de los parmetros de deformabilidad. de resistencia y del estado
tensional de los suelos sometidos a estudio. A pesar de la enorme popularidad de este
tipo de ensayo es preciso indicar cules son sus aspectos menos favorables:
Con el lin de superar algunos de estos inconvenientes y sobre todo para minimizar la
alteracin del terreno ocasionada por la perforacin previa al comienzo del ensayo, se
desarroll en los aos setenta [Baguelin et al., 1972. 1973, 1974; Wroth and Hughes.
19731 el presimetro autoperforante (SBP). Dicho procediiniento pretende situar en la
posicin y profundidad deseada la sonda presiomtrica sin vaciar el estado tensional y
deformr?tivo del medio circundante. Se puede decir que existen dos tipos de
presimetros autoperforantes: el francs (denominado Pafsor) y el ingls (denominado
Camkometer). Ambos se basan en el mismo principio; en la parte inferior de la sonda.
por debajo de ia zona de expansin, se s i t b un cilindro hueco de paredes rigidas y con
el borde inferior afilado. A medida que se va introduciendo la sonda se tritura el terreno
que penetra en el interior del cilindro hueco mediante un dispositivo rotatorio y
mediante un flujo de agua canalizado que atraviesa la sonda se arrastra a la superficie.
De este modo el medio permanece prcticamente inalterado durante la introduccin del
presimetro y el ensayo puede comenzar con un medio no alterado alrededor de la
sonda. El presimetro empleado en los ensayos que se analizan en esta tesis es la sonda
Camkometer Mark VIII.
Los ensayos presiomtricos con los que se trabaja en esta tesis -como fuente de datos-
estn realizados en celda de calibracin por lo que la alteracin del medio alrededor de
presimetro es minima. Con ello se consiguen condiciones mucho m h controladns que
si el ensayo se realizara in-situ. Los materiales einpleados en los ensayos son distintas
muestras de un suelo arenoso de la regin del Ticino (Suiza) cuyas caracteristicas
fundamentales se citan en cl Anexo 111. Adems los ensayos se llevaron a cabo n
condiciones de drenaje por lo que se produjo una deformacin volumtrica que debe
tenerse en cuenta a la hora de modelizar el material. Aqui se considera nicamente la
informacin procedente de la curva de carga pues la curva de descarga. en ensayos
drenados, es dificil de interpretar.
1.4.2.1 introduccin
A la hora de afrontar y resolver cualquier problema de tensin-deformacin en el
continuo. en un medio homogneo es preciso conocer el equilibrio del sistema de
fuerzas actuantes. la relacin entre las deformaciones que tienen lugar, las condiciones
de contorno del problema concreto y adoptar una ley de comport-.?iento o ley
constitutiva que relacione las tensiones con las deformaciones.
1.4.2.2 Hiptesis
Un ensayo presiomtrico (SBP) consiste en definitiva en el estudio de la expansin de la
cavidad preciomtrica teniendo en cuenta una serie de hip6tesis que simplifican el
problema. Dichas hiptesis son normdmente las siguientes:
clt)
yoe el rC. ;-r
co~isi~lerand<r resulta,
Siniplificando &.
Diferenciando la expresin de r.
CAPJWLO 1: PLANTEAMIENTOS BASICOS PARA EL DESARROLLO DE UNA
METODOLOG~ADE IDENTIFICACINDE PARMETROS
dr-dr, dr-dr+&-J;
E, c-= --
dro 4 0 J.,,
y de la deformacion tangencial,
Ee =
-
urco inicio1 orco zenerico
orco inicial
TE~=E ; d e , +r&, =e
y entrando en la expresibn de la deformacin radial (1.2) con este valor de 4,
de donde.
Antes de comenzar el ensayo las condiciones iniciales de presin y deformacin son las
si yuieiites:
:
P~sicin(p)
sobre la
membrana
Entre las leyes de comportamiento que se han empleado o suelen emplearse para la
modelizaciin de ensayos presiomtricos figuran: ley elstica lineal. ley elstica no
lineal, ley rigido plstica, ley elstica lineal perfectamente plstica y ley elasto-plstica
con endurecimiento y10 reblandecimiento.
En esta tesis con las hiptesis mencionadas se ha modelizado el suelo arenoso del
Ticino mediante el modelo hiperblico de Ferreira (19923. Dicho modelo describe bien
el componamiento del medio bajo una gran variedad de condiciones de carga y tiene en
cuenta la informacin de la variacin de la deformacin volumetrica aue tiene lugar en
los ensayos drenados. Adems es un modelo sencillo definido por parmetros con claro
-
significado mecnico que pueden ser empleados posteriorniente en el d~seopeoiccnico.
Por otro lado. se ha considerado interesante realizar una simulacin numrica del ensayo
presiomtrico utilizando el modelo elastoplsiico de Drucker-Prager (19521. La
modelizacin e Drucker-Prager se ha llevado a cabo con el tin de poder comparar los
resultados obtenidos mediante la ley analitica de Ferreira con un modelo no analtico
que precisa el empleo de un mtodo numrico. De este modo se pueden llevar 3 cabo
anlisis en dos y tres dimensiones, considerar condiciones de contorno diversas, trabajar
con materiales no homogneos... Hay que 9 , alar que algunas de las hiptesis citadas
anteriormente (apartado 1.4.2.2) no se tienen por qu mantener cuando el material se
niodeliza con el modelo de Dnicker-Prager. Al realizarse en esta tesis un anlisis
bidimensional axisimtrico del problema mediante un cdigo de elementos finitos la
longitud de la cavidad no se considera infinita. el campo inicial de tensiones puede no
ser constante con la profundidad y el suelo puede suponerse rstratificado
horizontaimente.
m
-toa
1.4.3 Tcnicas de
. . .. lemas inv
ajuste automtico por el contrario permite obtener los parmetros del modelo empleando
algn algoritmo iterativo de optimizacibn de la funcin objetivo (definida en esta tesis a
partir de la teota de la mxima verosimilitud; [Edwards. 19721).
Los inetodos de optimizacin llama .os "clsicos" utilizan. bien la derivadas primeras
(mtodos de gradiente), bien las derivadas segundas (niiodos de tipo Nelnon). Ciertos
mtodos combinan ambos mtodos, como los mbdos de tipo cuasi-Newon o los
metodos de mtrica variable (OFP. mtodo de Davidon. Fletcher y Powell; BFGS, n~todo
de Broyden. Fletcher. Goldfarb y Shmo: IF11:rcher. 19811; [Kirsch. 19961). En general,
tales miodos de optimizacin dan buenos i-esultados para problema que estn bien
condicionados. Sin enibargo, un buen niunero de problemas de optimiracin que sc
encuentran en niodelos reales presentan dificultades. como los no diferenciables. los de
funciones no convexas o aquellos con gran sensibilidad a ciertas variables (problemas mal
condicionados).
.4simismo. cuando el modelo no est definido por una simple ecuacin analitca, como por
ejemplo en el caso de la ley elastoplstica de Drucker-Prager. aparecen dificultades
suplementarias. El valor de la funcin objetivo para cada valor de los parmetros no puede
obtenerse salvo si se emplea la simulacin numrica. Ademis puede tratarse de la
simulacin de un proceso descrito por un conjunto de ecuacianes diferenciales o de la
sin~ulacionde un sistema discreto. As, en muchos casos, ser muy dificil y pesado
calcular una buena estimacion de las derivadas de la funcion.
Es, en estas situaciones, donde los mtodos clsicos de optimizacin zncuentran lmites de
aplicacin y los mtodos directos [mtodo Simplex (Nelder y Mead, 1965); mtodo
Complex (Box, 1965); mtodo de los Algoritmos Genticos (MichUlwicz 1996)] resultan
interesantes. Los mtodos directos utilizan nicamente el valor de la huicin objetivo en
ciertos puntos para calcular un nuevo punto, o vector de parmetros que mejore el valor de
la funcin objetivo.
Por otra parte, una dificultad que aparece frecuentemente cuando se lleva a cabo una
minimizacin de una fwicin objetivo es que dicha funcin presente varios mnimos
locales que incluyan al minimo global o absoluto. El problema de encontrar un extremo
global, mnimo global en nuestro caso, suele ser un problema complicado [Horst. 19901.
Para lograrlo se suelen emplear tradicionalmente dos tcnicas: la primera consiste en
encontrar el mayor nmero de mnimos posibles, para lo cual se comienza con una
variedad muy grande de juegos de parhmetros iniciales, en muchos casos aleatoria. y se
termina seleccionando el extremo ms bajo obtenido -siempre que se haya encontrado ms
de uno-; la segunda t6cnica consiste en, una vez alcanzado el extremo local, perturbar el
juego de parhmetros que lo definen y ver si con los nuevos parjsnetros la funcin objetivo
vuelve al extremo inicial o bien alcanza uno nuevo. Adems, existen tcnicas automticas
como el mtodo de Recocido Simulado [Kirkpauick, 1984).
Tambin hay que tener eii cuenta que a la hora de buscar el mnimo global de una funcin
objetivo los valores de los parmetros que hacen que la huicin alcance dicho valor
extremo no pueden ser cualquiera sino que existen una serie de restricciones tisicas que se
deben verificar. Estas restricciones proceden en general de la informacin previa obtenida
en ciertos experimentos adicionales o bien por limitaciones de los parmetros en la
delinicin de la ley de coniportamiento.
As. para elegir un mtodo de optimizacin para problemas inversos se debe considerar en
primer lugar si se elige un mtodo que precise evaluar nicamente la funcin objetivo para
ciertos valores de los parmetros o bien, si se elige un mtodo que ademis de requerir
dicha evaluacin requiera la evaluacin de las derivadas parciales de dicha funcin
objetivo respecto de cada uno de los parrimetros del modelo. En general, los algoritmos
que emplean las derivadas de la tuncin objetivo son computacionalmente ms rpidos
que los que emplean informacin slo de la funcin si bien. en cienas r~:asiones. el
resultado obtenido no compensa el esferm de progmmacin requerido.
En esta tesis se han empleado bsicamente dos algoritmos de minimizacin. Uno de cada
tipo. Para el modelo hiperblico -analtico-, cuyas derivadas resultan computables, se ha
empleado un algoritmo clsico -algoritmo de Marquardt, [1963]- y otro directo -algoritmo
tipo complex, [Box,1965]- y para el modelo de Dmcker-Prager. implementado en un
nitodo de elemenws finitos. se emplea nicamente el algoritmo directo tipo complex.
Se analiza tambin el problema mencionado de encontrar los extremos globales. Para ello
se presentan los denominados mtodos de recocido simulado (simulated annealing). Estos
mtodos han demostrado su validez en la resolucin de ciertos problemas definidos tanto
en espacios discretos como continuos, en particular en el caso del modelo hiperblico.
1.4.4 ~ c s i o m 6 t r i c o s
Cada uno dc estos ejemplos sintticos se ha resuelto con los diferentes algoritmos
presentados, y con los dos modelos de comportamiento elegidos. Sin embargo se ha
Iiccho un anlisis ms exhaustivo con el modelo de comportamiento dc Ferreira al ser
este analtico. En el caso de la modelizacin de Drucker-Prager se han tenido en cuenta
tambin diferentes tipos de mallado y de elementos finiros.
1.5 CONCLUSIONES
Por ltimo decir que en este captulo se han presentado las bases del desarrollo que se ha
seguido en este trabajo de tesis con el fin de que queden claros los objetivos de los
captulos sucesivos. su motivacin y el lector los pueda abordar teniendo una visin
global del problema que se quiere afrontar.
O MODELIZACI~NMECANICA DEL ENSAYO PRESIOMTRICO
C A P ~ L 11:
CAPITULO 2
2: MODELIZACI~NMECANICA DEL ENSAYO PRESIOMTRICO
CAP~TULO
Sin embargo. ante ensayos presiomtricos drenados aparecen variables adicionales que
deben ser tomadas en consideracin, tales como la deformacin volumtrica (E,.) y el
cambio en la tensin principal normal. En este caso ya no se tiene nicamente la
condicion de cortante puro.
Se puede decir por tanto, que la respuesta que proporciona un ensayo presiomtrico
puede ser no drenada ylo drenada. La primera se produce cuando no se permite el
MECNICA DEL ENSAYO PRESIOMT~ICO
2: MODELIZACI~N
CAP~TULO
De entre los modelos citados se ha elegido, por las ventajas que se indican en el
prbximo apartado. el modelo hiperblico propuesto por Ferreira. Dicho modelo emplea
en su desanollo propiedades de la aproximacin llevada a cabo por Hughes et al (1977)
y por Houlsby et a1 (1986).
2.3.1 Modelo h i n e r w o de F e r r e i r m
Estas suposiciones son en esencia las formuladas por Gibson y Anderson (1961) y por
Jefferies (1988). exceptuando el tipo de niodelo y la restriccin sobre la magnitud de las
defoi;?iaciones.
En ensayos no drenados, la representacin del modelo hiperblico viene dada por una
nica curva cortante que describe la relacin tensin-deformacin en el plano tensin
cortante/deformacin de la cavidad. (Figura 2.4.a).
En los ensayos drenados no ocurre lo mismo, sino que aparece una familia de curvas
cortantes (Figura 2.4.b). cada una de las cuales es funcin de la tensin principal normal
efectiva (S,) que cambia a lo largo del ensayo. Durante la fases de expansin y
contraccin hay una migracin de una curva a otra en funcin del nivel alcanzado de
tensin principal normal efectiva. Para evitar estos saltos de una curva cortante a otra se
puede normalizar la tensin cortante en cada momento (1) con la tensin principal
normal efectiva correspondiente (S). De este modo se consigue que una nica curva
represente la fase de carga del ensayo presiomtrico.
A 4
,,E : E"mT
0 5,
1
y--;
'lensibn Tensibn
conante (1)
Ijj
/
conanie (1)
+
Defomaci6n de la cavidad (E) Defornacin de la cavidad (E)
Figuras 2.4. Modelo hiperblico en ensuyos no drenados (a) y drenados con crrurro
niveles de rensin principal normal e/ccrivu (b).
Convenio de s i ,
Las deformaciones y las tensiones normales de compresin son positivas. Las
deformaciones volumtricas son positivas si hay reduccin de volumen, en caso
contrario se& negativas.
CAP/TULO2: MODELIZACINMECANICADEL EPiSAYO PRESIOM~TRICO
. ., , .
p
En este desarrollo se considera una relacin lineal entre la deformacin volumerrica y la
delbrmacion de ia cavidad. Teniendo en cuenta el convenio de signos se tiene:
AV
E,. = - - = - S E
V" '
Expansin
Contraccin
scn V= -&,J&
Arena densa
Expansion
/--
- - - . A
Pmcesode carga
l\"
_I___)
1 1
ti' Deformacin cortante (7) *
contmcci'in "
Arena suelta
.,
La relacin tensindeformacin usada para representar el modelo hiperblico se define
a partir de (2.1.) mediante el ratio de tensiones (tsj wmo funcin de la deformacin de
la cavidad ( e ):
2: MODELIZACINM
CAP~TULQ .CANICA DEL ENSAYO P R E S I O M ~ I C O
$=-
fs',+ o',
*
Deformaci6n de la cavidad ( E )
de la c a v W
Se considera que cuando la deformacin de la cavidad presiom6trica es nula (E = O ) la
tensin eficaz radial coincide con la tensin efectiva horizontal inicial (o',=o',,).
tro en la fase de
De la definicin dada anteriormente del ratio cortante movilizado (ILS) se puede
determinar la siguiente relacin entre la tensin efectiva radial y la tensin efectiva
circunfcrencial.
dc donde:
luego.
Dicb i relacin se puede expresar en funcin del ngulo de fnccin de pico ($ =a)si se
tiene en cuenta la expresin que liga las tensiones con las deformaciones en una
I
expansin plhstica (- = constante = =n$' ) del modo siguiente.
S
'
S",
do', = ,A4,-'d r = - -dr
- --
r
(1 - A'(&))
0 ',
Suponiendo que la deformacion volumtnca que tiene lugar puede definirse mediante la
relacion lineal con la deformacin de la cavidad.
\Vorth and Windle (1975) presentaron la siguiente ecuacin que liga la deformacin de
la cavidad ( c ) con la coordenada radial ( r ):
MECANICA DEL ENSAYO PRESIOMTRJCO
2: MODELIZACI~N
CAP~TULO
S
2 sen v
=-=l-n
' I+senv
donde,
1 - senv
n =-
1+ senv
A partir de (2.4.). teniendo en cuenta estas dos ltimas expresiones y despreciando los
trminos de mayor orden de la deformacin de la cavidad -trminos c' y superiores- se
llega a la ecuacian diferencial:
como adems.
finalmente se tiene,
donde.
I'or otro lado el valor de ( t r ) puede determinarse tambin (Ferreira. 1992) en funcin del
riio final de ieiisiones ( ( r /.S),, ) mediante la expresin:
donde.
1+
K," (1 f ( 1 1 S)"~,
donde e es la base del logaritmo nepenano. Normalmente, para las arenas que se
estudian en esta tesis. se considera que valor del ngulo de friccin a volumen constante
(4 '& est prximo a 34'.
A partir de estos p,umeiros del modelo de Ferreira se pueden obtener los siguientes
parametros del suelo:
..
Mdulo elastico corta m i n h i d ( G , )
Este mdulo se determina a partir de la pendiente de la tangente inicial a la curva de la
relacin hiperblica (S,).Al comparar las curvas hiperblicas [ ( I 1S ) - E ] y [ ( I 1s) - y ]
Ferreira [1992] obtiene la relacin siguiente:
. .,
(4)
lo de f r i c c u u h .p s ~
Teniendo en cuenta el criterio de rotura de Mohr-Coulombse tiene:
2: MODELIZACINMECANICADEL ENSAYO P R E S I O M ~ ~ C O
CAP~TULO
Aneulo dev 1
Conocido el ngulo de friccin a volumen constante y por tanto K Y . el ngulo de
dilatancia se obtiene a partir de la ley de la dilatancia de Rosve [1962 y 19721 del
siguiente modo,
Historicarnenie se puede decir que fue Coulomb (1773) el que propuso el primer criterio
de fluencia para suelos. Dicho criterio se comporta como un criterio de rotura y viene
dado por la expresin,
MECANICA DEL ENSAYO PRESIOMTRICO
2: MODELIZACI~N
CAP~TULO
es decir,
o bien.
2: MODELIZACINMECANICADEL ENSAYO PRESIOMETRICO
CAP~TULO
donde o, y o> son respectivamente las tensiones principales mayor y menor. Como
puede apreciarse el criterio de Mohr-Coulomb no tiene en cuenta los efectos de la
tensin principal intermedia op
Cuando se cumple que o, =o,= o ,se obtiene de (2.8) la tensin hidrosttica principal
c m= -ccotg 4. indicando que el vrtice de la pirmide hexagonal se sita sobre el
espacio diagonal (Figura 2.9.) en el punto o, = o 2 = a, = -c cotg 4 .
El comportamiento del material antes de alcanzar el nivel de fluencia est definido por
la ley de Hook generalizada. que relaciona el tensor de tensiones (o,,)con el tensor de
deformaciones (E,,):
Las constantes de Lam estiin relacionadas con el mdulo elstico de Young (E) y con
el coeficiente de Poisson (v) por las siguientesexpresiones,
los que se cumpla f=K habrii deformaciones plsticas mientras que no las habr, sern
iinicamente elsticas, en aquellos estados de tensibn eii los que se cumpiaf<K.
se pueden definir los invariantes de tensin (J,. J , J,) > los invarimtes de iension
desviadora (-1, '. J? : .l. 1 referidos respectivamente al tensor de tensiones o,,y a! tensor
de tensiones desviadt .iu',,. Las espresiones de cada uno de ellas son las siguientes:
Eii los ~iiaterialesgranulares no iucede esto y el criterio de tluencia debe tener en cuenta
la componente Iiidrosttica de las tensiones. Modificando el criterio de fluencia de Von
Mises e iniroduciendo la influencia J e la componente hidrosttica de las tensiones en la
Ilucncia se llega :a la expresion de1 criterio de fluencia de Drucker-Pnger (1952):
donde:
MECANICA DEL ENSAYO PRESIOM&TRICO
CAPTULO2: MODELIZACI~N
2 sen 41 6c cos4
a= f'=
&(3-sen4) &(3-5en+)
Cuando se produce por primera vez la fluencia el comportamiento del material ser
parcialmente elstico y parcialmente plktico. Por tanto. se considera que los cambios en
las d~.lbrmacionespueden dividirse en dos componentes. una elstica (d~fll) y otra
plstica (c/~;,fl):
(2.1 O.)
Ley de flujo
La ley de flujo es la relacin entre la componente de la deformacin plstica y el
incremento de tensiones. Dicha ley gobierna el flujo plstico despues de la fluencia.
Esta relacin considera que el incremento de la deformacin plstica es proporcional al
gradiente con respecto a las tensiones de una cantidad llamada potencial plstico L). es
decir,
(2.1 l.)
donde d?, es el multiplicador plstico. En general se puede suponer que este potencial
va a ser tambien funcin de los invariantes J'? y J'3.
af
-
Superficie de fiuencia
proceso de carga
En algunos materiales. entre ellos en los suelos. piiede ocurrir que la superficie de
fluencia disn~inuyaal aumentar la deformacin plstica. Cuando sucede este fenmeno
-
ya no se Iiabla de endurecimiento sino de ablandamiento. En el modelo istropo que es
el modelo que se considera en este trabajo - el ablandamiento supone que la superficie
de fluencia se contrae a medida que aumenta la deformacin plstica. En este caso la
fluencia conlleva una rotura local y la superficie de fluencia se convierte en un criterio
de rotura.
Cuando tl/<Il se produce una descarga elstica y el punto que se est estudiando vuelve
al interior de la superticie de fluencia, donde el comportamiento es eiktico. Cuando
dfk0 indica que no hay aumento de carga por lo que el punto permanece sobre la
superficie de fluencia y se comporta plsticamente -comportamiento perfectamente
plistico-. Por ltimo cuando df>U tiene lugar una carga plstica y el punto que se est
estudiando permanece sobre la superficie de fluencia mientras esta se expande -
comportamiento plstico para un material con endurecimiento pliistico-.
- La tensin de fluencia K.
..
2.4.6 F o r m u l a c i o n l modelo e l a s t o m
En este apartado se desarrolla una formulacin matricial [Owven et al, 1986) de las
expresiones vistas hasta ahora referidas al modelo elastoplstico de Drucker-Pnger con
el iiii de poder llevar a cabo una resolucin numrica del problema elastoplstico que
nos ocupa, la expansin de la c a v i d J del presimetro en medios gmnulares.
diferenciando,
Esta expresin puede tambin escribirse como:
fi'd~- A A = O
u'do = AW
llamando.
A partir dc aqu la expresin de la relacin incremental completa que liga tensiones con
detbrniaciones elasioplsticas (2.12.) puede escribirse:
Esta ltima expresin permite escribir la relacin incremental completa (2.13.) anterior
como:
donde.
o bien.
y comparando trmino a trmino la ecuacin (2.17.) con esta ltima ecuacin se obtiene:
Una altcmativa dei tercer invariante de tensin desviadora ( J ; ) es !a priniera raz en (8)
de la expresin anterior considerando que 38 puede variar en un rango de A
panir de la naturaleza cclica de sen(38 +2m)se obtienen los tres posibles valores del
seno que definen las tres tensiones pnncipalris. Las tres tensiones principales
C A P ~ ~ U LMODELIZACIN
O~; MECNICA DEL ENSAYO PRESIOM&T~(ICO
se aprecia que no depende del tercer invariante. Sin embargo, cuando se quiere calcular
la matriz D,, es preciso expresar el vector de flujo uT en funcin de J, J , y O para que
sea coherente la nueva fomulacion. Es decir.
donde o
' = [o,.oy.o~.r,z.r~.rT 1 .
o bien,
c ~ ~ i n r2:t MODELIZACIN
o MECNICA DEL ENSAYO PRESIOMBTRLCO
cal
y sustituyendrr la ltima expresin de 7en la expresin del vector de flujo (2.19.) se
CYT
obtiene.
dc dondc,
> operando.
donde,
al'
c, = -
a,
2.4.8.3 Axisimetrla
Esta ha sido la simplificacin empleada en la modelizacin del ensayo presiomtrico
con el modelo elastoplstico y por unto se le presta una atencin mayor.
CAP~
2:U O MECANICA DEL ENSAYO PRESIOMTRICO
MODELIZACI~N
1.::L[", . % ..
Cargas y c. contorno
axisim6iricas
Las componentes del vector de flujo segun la formulacin alternativa empleada para la
programacin de las diferentes superficies de fluencia es la siguiente.
donde.
L'AP~TULO2: MODELIZACI~N
MECANICA DEL ENSAYO PRESIOM~TRICO
u,' = (1,l.O.l)
con los valores de las constantes C,, C2 y C,. para el caso de Dmcker-Prager. citados
anteriormente en (2.20.).
d,, = D.
de donde.
2: MODELIZACI~NMECNICADEL ENSAYO PRESIOM~TRICO
CAP~TULO
Ev(a, + u2+ )
con MI =
(l+v)(l-2v)
2.4.9 Ecuaciones cl .. .
corOperOMiento e m s t i c o de un solido .
Para ello. a partir de las ecuaciones r.onstiiuiivas y del principio de conservacin, que
establece que dentro de cualquier subdominio o del dominio C(Figura 2.13) de nuestro
problema el flujo neto que atraviesa el contorno do del subdominio debe ser igual a la
suma producida por las h e n a s interiores o mAsicas. se puede decir que localmente se
verifica:
donde o' = (o, ,o,.T, .a,) son las componentes mdial, axial. desviadora y tangencial
. ..
2.4.11 Simulacin numfrica dcl c m o n r c s ~ t r i c out ilizando el rnodel~
2.4.11.1 Introduccin
A continuacin se va a describir como se simula el comportamiento del material durante
el ensayo resolviendo numericarnente el problema en ecuaciones derivadas parciales
(2.24.). La metodologia numrica empleada en este trabajo ha sido el mtodo de
elementos finitos.
[.as dos primeras etapas se han descrito ya en los apartados 2.4.9 y 2.4.10
respecti\ariientr. El resto se estudian a ccntinuacion.
. . ., , . . ..
2.4.1 1.7.1 Prohlema de contorno v g a a c i o n a l . g s
trabriios virtuale~
Para establecer la formulacin variacioiial del problema de contorno se considerar en
principio un problema elstico ( 0 = D E ) y teniendo en cuenta (2.24.) se define el
residuo tp mediante:
CAPITULO~: MODELIZACI~NMECANICA DEL ENSAYO PRESIOMTRICO
dcnde r/r.z) son las fuerzas que acian por unidad de superficie sobre el contorno ?a,.
Ilc cstc modo las condiciones de contorno n~turalesquedan incluidas en la misma
definiciOii del problema.
I'or tanto el problenio de contorno \rariacional se esiiblece como: buscar una funciOn rr.
para cualquier funcin de prueba admisible i*.que \,erifique (2.25.) y que cumpla las
coiidicioiies de contorno esenciales sobre a,.
1.a ecuacion (3.25.) ser5 la formulacin dbil o variacional del problema de contorno en
dos diiiiensiones. La condicin que deben cumplir tanto 11 como t. para que ellas y sus
primeras derivadas sean suficientemente regulares es que pertenezcan al espacio H 1 Q).
(
Es decir que diclias funciones tengan sus derivadas prinieras de cuadrado integrable en
el dominio R.
Siempre que las fuiicioiies de pmeba. los datos y la solucin del problema variacional
sean sulicientemente suaves cualquier solucin del problema variacional sera solucion
del problema clisico o fuerte. Por el cantmio cualquier solucin del problema tuerte es
siempre so1uc;n del problema variacional.
Por otro lado se puede llegar a la formulacin variacional del problema de contorno
mediante el principio de los trabajos virtuales. Dicho principio establece que iin sistrma
esta en equilibrio ante fuerzas exteriores. si al imponer sobre el mismo unos
desplazamientos arbitrarios infinitesimales (virtuales. r), compatibles con la
2: MODELIZACI~XMECANICA DEL ENSAYO PRESIOMTRICO
CAP~TVLO
sustentacin del sistema y a partir de una posicin de equilibrio, el trabajo realizado por
las fuerzas exteriores es igual al realizado por las fuerzas y tensiones interiores.
. .
7.4.1 1.2.2 medianre l mtodo de elementos m s . A ~ r b c i nde
Galerkin.
El problema variacional es el punto de partida idneo para construir soluciones
aproximadas mediante el mtodo de eleme:itos finitos. El mtodo de Galerkin consiste
en buscar una solucin aproximada del problema de contorno variacional, en vez de en
el espacio H'(R) de funciones de pnieba admisibles, en un subespacio de ste de
dimensin finitaH1'(Qh). El dominio R,, se obtiene al dividir Q en un conjunto E de
elementos finitos de n nodos cada uno. Cada elemento, que a partir de ahora llamaremos
c. estar formado por los correspondientes n nodos y por las rectas que unen dichos
nodos y que constituyen el contorno del elemento. En los codos van a estar definidas las
incgnitas de nuestro problema. concretamente en nuestro caso. los desplazamientos u.
Ahora bien, interesa ademas conocer de un modo nico cul va a ser el campo de
despla7amie;iros en cada punto del elemento e en cuestion. Para ello se toman un
conjunto de funciones (base del subespacio Hh ) que nos definen el desplazamiento cn
el interior de cada elemento a partir de los desplazamientos nodales del mismo. Estas
funciones de desplazamiento o funciones de forma, permiten definir tambin de un
iiiado nico el estado de deformacin en cada punto del elemento. que junto con el
estado de deformacin inicial y las propiedades constitutivas del material permiten
definir el estado de tensiones en todo el elemento. Por ltimo, se establece un sistema de
fuerzas concentradas en los nodos que equilibra las tensiones en el contorno y las cargas
repartidas en cada elemento.
El vcctor 11. que define el desplazamiento en cualquier punto del elemento e del
dominio R y que constitu).e en este caso la incgnita del problema, puede ser
aproximado mediante un vector u'. Este vector u' define los desplazam;entos en
cualquier punto del elemento mediante unas funciones de posicin o funciones de forma
(.Y,) .ri partir de los desplazamientosnodales (u,). Es decir,
radial y axial del nodo i, y el vector a' los desplazamientos de los nodos del elemento c.
es decir.
4
ir' = 1A','u:
lil
Una vez definidos los desplazamientos para todos los puntos del eleniento pueden
determinarse las deformaciones en cualquier punto del mismo. Las deforn~aciones
pueden expresarse .segun lo visto anteriormente en (2.23.). como el producto de un
operador lineal S por el vector desplazamientos u. es decir:
considerando la Iiiptesis de pequeos desplazamieiitos. Donde E ' = (E, .e,. y, ,e, ) son
I d componentes de la deformacin radial, axial. cortiuittt y circunferencial
respectivamente.
donde el \.ector de tensiones a' = (u, ,G:.T, ,o,,) est fomiado por las componentes de
la tensi011radial, axial, cortante y circunferencial respectivamente, y U es la matriz
constitutiva que contiene la; propiedades (elsticas) del material.
- Las fuerzas msicas :1>) que son por definicin las que actiran por unidad de
volumen en las direccionesde los desplnzamientos ir en dicho punto.
- Las fuenas ( r ) que actan por unidad de superficie sobre el contorno del
elemento.
- Las herzas externas puntuales (p) sobre los nodos del elemento.
- Las fuerzas nodales equivalentes (y) que son aquellas que actuando sobre
los nodos son estticamente equivalentes a las tensiones en el contorno y las
fuerzas msicas que actan sobre el elemento.
MECANICA
2: MODELIZACI~N
CAP~TULO DEL ENSAYO PRESIOMETRICO
Una vez descritas las caractersticas de un e!emento finito cualquiera el siguiente paso
consiste en establecer un equilibrio entre las fuenxs nodales equivalentes y el resto de
f u e m s y tensiones actuantes sobre el elemento. Para ello una posibilidad es recurrir al
principio ya enunciado de los trabajos virtuales. Sea 6 u' el vector de desplazamiento
virtual de los nodos del elemento e. Dicho desplazamiento va a originar unos
desplazamientos (S t ) y unas deformaciones (6 E ) en el interior del elemento que
vendrh dados por las expresiones siguientes,
El trabajo externo realizado por las fuerzas nodales equivalentes (y) es igual a la suma
de los productos de las cowponentes de cada una de las fuerzas por sus correspondientes
desplazamientos, es decir.
El trabajo interno por unidad de volumen efectuado por las tensiones y fiiems masicas
motivado por el desplazamiento virtual 6 a'. viene dado por,
Por ltimo. el trabajo externo realizado por las f u e m i r I que actuan sobre el contorno
del elemento est definido por:
y el trabajo externo realizado por las fuerzas puntuales @) sobre los nodos por la
expresin,
C A ~ ~ T U L2:OM O D E L I Z A C I ~MECANICA
N DEL ENSAYO PRESIOMETRICO
Como esta espresi0n. por el principio de los trabajos virtuales, debe verificarse para
cualquier desplazamiento viriual arbitrario se puede escribir:
donde.
CAPITULO2: MODELIZAC16N MECANICA DEL ENSAYO PRESJOMETRICO
El principio de los trabajos virtuales aplicado a todo el continuo (Q) supone que el
trabajo total es igual a la suma de las contribuciones de todos los elementos.
donde las integrales implicadas en cada trmino podrin resolverse siempre que sus
integrandos sean finitos. Es decir, se considera que no se producen aiscontinuidades
entre elementos adyacentes. Si apareciesen esas discontinuidades habra que aadir una
contriburion igual al trabajo efectuado por las tensiones en las separaciones o
sol~pamientosde los elementos.
Para asegurar que las inc6gnitas (u,&) presenten continuidad. o como mucho saltos
finitos entre los elementos. es preciso que las funciones de forma sean de clase c". Con
ello se asegura la continuidad en los desplazamientos y que las deformaciones. que se
obtienen a partir de las derivadas primeras de funciones de forma continuas, presenten
en el caso mas desfavorable saltos. que al ser finitos, no impedirn que se pueda llevar a
cabo la integracion. En el caso de que la deformacin estuviese definida por derivadas
segundas de las funciones de forma entonces las funciones de fonna deberan ser al
menos de clase C.!.
A la Iiora de realizar los clculos para determinar las contribuciones de cada elemento r
a cada uno de los trminos del problema (2.28.) es evidente que dichos clculos s e r h
diferentes por ser diferentes los elementos en que se ha dividido el dominio R. Con el
fin de evitar la tediosidad en los clculos y homogeneizarlos se disea un elemento tipo,
sencillo y se define una transformacin invertible entre el elemento tipo y cualquier
elemento e, de tal modo que todos los clculos se lleven a cabo sobre dicho elemento
tipo que se define en un sistema (cq)de coordenadas locales (Figura 2.1 5).
A la hora de definir las funciones de forma se emplean las coordenadas locales (5.q).
iiormalizadas. Dichas coordenadas se eligen en el elemento rectangular tipo, de manera
que tomen los valores 1+ 1 en los lados del rectngulo. Es decir,
cAP/T(;LO 2: MoDELIZACIN MECANICA DEL ENSAYO PRESIOMETRICO
donde a y h son magnitudes que representan la mitad de la longitud dr los lados del
elemcnto c. de lados paralelos a los ejes (5.~1).y (&y,) las coordenadas del origen del
sistema de coordenadas local (Figura 2.15).
donde los valorcs de las coordenadas locales para cada uno de los nodos son los
expres~dosen la Figura 2.1 5 :
Mediante las reglas de derivacin en cadena se pueden obtener para el nodo i del
elemento e las derivadas de las funciones de forma respecto a (5.q).
o bien.
Por tanto las derivadas de l a funciones de forma respecto de las coordenadas globales.
que son las que interesan, se obtienen mediante:
.i\l aplicar esta misma tnnslormacin al volumen elemental de integracin ,!L este
pasari de expresarse en coordc.nadas globales como d., = 2n r'drclz a expresarse en
coordenadas locales como d , , = 2n r' de: Jdt 4 .
donde los ndices i,j vqran desde 1 al nmem n de nodo- del elemento, cuairo en el
elemento tipo con el que estamos trabajando.
Si se llama T y a l integrando 7;; = (B," )' D'B: 2n r' det .l entonces.
El siguiente p iso consiste en calcular las integrales. Dicha iniegracijn se lleva a cabo en
coordciiadas locales y se resuelve numricamente calculando la integral en uria serie de
puntos de coordenadas (6.q) conocidas llamados puntos de integracion de Gauss. Si se
considera como Iiasta ahora un elemento lineal de cuatro lados con nxn (n=2) puntos de
integracin en su interior se tendr:
donde N> y W, son las funciones de peso de la integracin y (<,q,) las posiciones de
los puntos de Gauss dentro del elemento tipo.
2- MODELIZACIN MECANICA DEL ENSAYO PRESIOM~TRICO
CAP~TULO
En In Figura 3.16. se muestra la situacin dentro del elemento e de los cuatro puntos de
integracin de Gauss de una regla de integracin nunrica 2.~2.Los anlisis que se
presentan en el prximo capitulo se han llevado a cabq tanto con elementos de cuatro
nodos y orden de integracin 2x2. como con elementos de ocho nodos y orden de
integracibn 3x3.
llna vez que se tiene la informacin correspondientc de cada elemento por separado se
pasa a ensanlblar toda la informacin. La informacin de la matriz de rigidez de todo el
doniiniu 0. correspondiente a un nodo comn a varios elenientos. se obtiene sumando
las contribuciones f:t cada elemento a ese nodo. El vector de fuerzas globalesj que
acta sobre cada nodo se ohtiene de igual terma. sun~andolas contribuciones de cada
eleniento a ese nodo. El sistema de ecuaciones se completara teniendo en cuenta las
coiidicioiies de contorno y las condiciones iniciales. Dicho sistema tiene tantas
ecuaciones conio nodos '-9 tiene el recinto R. la matriz de rigidez es una matriz
cuadrada (.A'.K\') 1. siniplificadamente se puede expresar como:
de donde.
,,
donde LJ it~cluyela ley de comportamiento elastoplstico del material cn cuestin. D,,
viene da,ia pnr (2.15.). A partir de la informacin correspondiente a cada uno de los
elementos del dominio R y siguiendo los pasos descritos para el caso lineal. se forma el
sistema de ecuaciones, que en este caso ser no lineal,
El sistema de ecuaciones planteado que se tiene que resolver viene dado por (2.30.).
Dicho sistema tendr que resolverse de modo iterativo. Para ello es preciso definir el
residuo en cada iteracin. Se denomina residuo de la iteracin i-sima a 1 t diferencia
entre la fuemi aplicada en la iteracion i v las fuerzas nodales equivalentes calculadas en
dicha iteracin. es decir.
Con el fin de no tener que invertir en cada iteracin la matriz de rigidez, el mtodo de
Newtoii-Raphson puede escribirse tambin en forma incremental:
De esta forma el sistema no lineal (2.30.) se reduce a la resolucin en crl.. ikracin del
sistema lineal (2.31.) que permitir obtener el vector intermedio y'. Con este vector y
con el desplazamientoen la iteracin i se obtiene el desplazamiento U"'.
Los sistemas lineales (2.3 1 .) y (2.32.) se resuelven por el mtodo de eliminacin directa
de Gauss. Este mtodo permite transformar un sistema de N eiuaciones en un sistema
triangular equivalente tras A'-1 transfomaciones. Este nuevo sistema se resuelve de
forma inmediata mediante un proceso de remonte. Se podran haber aplicado mtodos
MODELIZACI~NMECANICA DEL ENSAYO PRESIOMETRICO
CAP/TULO~:
Al comparar las fuerzas nodales equivalentes con las fuerzas aplicadas se obtienen las
fuerzas residuales y los residuos, que son necesarios para poder analizar la convergencia
del proceso iterativo no lineal. El criterio de convergencia empleado es el siguiente:
INCREM
taniris ecagaacucslam~l~la
cm
DII~esb--dc1~1. ni
-*@ans
4
ALGOR
cC.wdalarOlCmm-.-an~sia
sitpomapmmo*apmsiimads~
I
SnFFP
E~imea'rdBneOQMparamesear
mu m p a Y n s M . e r t C O c m > < > - ~ ~ 3
FRONT
Erll~dsicrnmLOn~~oPmsa~(onn~irmm
sinrapmamcspa"nrush.Isstaadr
~ r o S U M w ~ m ~ m a i n x m m ~ d e G a s s
r
RESIDU
I INVAR
E-h(yrmpdma
YIELDF y FOWPL
Odlrni"a*uQI*
Los datos obsenrados, es decir los datos procedentes del experimento en cuestin se
agrupan en vectores columna y, con k=l. 2 ..,,N. Para ordenar adecuadamente los datos
experimentales el subndice k puede considerarse como un parinietro de discretizacin
temporal. Cada vector columna yA=(yl. ,v2....yJ7' contendr los desplazamientos.
presiones. presiones de poro ... que se han medido en ese "instante" k. fin general no es
demasiado importante el sigiiificado fsico preciso de las magnitudes observadas sino
unicamente el hecho de que en el experimento .a habido una serie de magnitudes que se
han podido medir en diferentes instantes (k).
donde:
Cada columna de datos y* puede tener una depeiidencia funcional de los parimetros x
diferente, que se seaala con el subindice k de la funcin h. Es decir. si los datos
CAP~TULO3: TCNICAS
NUMERICAS DE OPTIMIZAC16N PARA LA IDENTIFICACIONDE
PARMETROS GEOTCNICOS
disponibles procediesen de experimentos distintos -no es este el caso que nos ocupa- se
tendran que emplear diferentes modelos de material. que veiidrian dados por el
subndice correspondiente k.
Cuando la dependencia del modelo sea lineal con los parmetros x se podr expresar el
modelo matemtico znterior del siguiente mudo.
El otro problema es determinar lo cerca que se encuentran los parmetros estimados .Y,.
de los parmetros verdaderos x,. Al ser normalmente desconocido el error numerico
entre los parmctros verdaderos y los estimados (x, -x,) el problema es desarrollar un
modelo de inccrtidumbre para dicho error.
donde el estimador de niinimos cuadrados x,,, ser una funcion lineal de y,. 1.a hncin
S, se dcnoniiiia tambin tuiicion objetivo y como puede apreciarse en (3.1.) est
brmada por una sutira de tCrminos al cuadrado denominados residuos r, es decir. r =
yk-If$.
Llamando Pi a:
se tiene la expresin,
donde Pi existe si IV, es definida positiva y Hkes de rango n. dimensin del vector de
parametros x.
En este apartado se describe una metodologia que define una funcin objetivo que
permite. al mininiizar su valor. estimar los parametros de comportamiento ms realistas
del modelo de coniportamiento que se est utilizando.
1-U \,erosiiiiiliiud de una hiptesis (valores fijos para los parmetros) se define como una
iiiagniiud proporcional a la probabilidad de observar las medidas obtenidas en el ensayo,
dados los partmetros [Ed\vards. 19721. En pocas palabras se puede decir que el mtodo
de mxima verosimilitud consiste en maximizar la verosimilitud de una cierta hiptesis
-con,juiito de parametros- a la vista de las observaciones medidas. Por tanto el mtodo de
iiijxirri~verosimilitud es un modo concreto de afrontar el problema inverso, la
estimacibn de los parmetros del inodelo de comportamiento a partir de las respuestas
del sistema y de la informacibn previa ponderada adecuadamente.
P(Y) = -
-~')'(c,)-'(yy')
1
y con ello, a partir de (3.2.). se puede escribir:
-
donde el residuo (y yj es el vecto~de diferencias entre las variables calculadas con el
modelo y y la variables medidas y , ni es el nmero de medidas y C,, es la matriz de
covarianzas de los errores de las medidas. Cuando estos errores son independientes entre
si, la matriz de covarianzas C, es diagonal.
donde los dos ltimos sumandos son fijos y no influyen por tanto en el proceso de
identificacin. Si ademis la estructura del error de las medidas (C,)se considera lija
entonces la funcin que habr que minimizar vendr dada nicamente por el primer
termino de la expresin anterior. es decir:
Si se supone adems que los errores de medida (C,) son independientes entre si, se
puede escribir la matriz de covarianzas del siguiente modo.
De tal modo que (3.4.). despreciando los dos ltimos trminos fijos. resulta ser:
Como a prior es muy dificil conocer el valor de las varianzas de las medidas (S,' ) y con
el fin de realizar el proceso de identificacin sin necesidad de aiiadir hiptesis sobre su
valor, normalmente se toman las varianzas como variables a identificar dentro del
proceso de identificacin [Carrera, 19861.
C~fin/tO
3: T ~ I C A NUMRICAS
S DE OPI~MIZACI~N
PARA LA lDENTlFlCACI6N DE
PARAMETROS GEOTCNICOS
Suponiendo que no hay relacin entre ia informacin previa y las medidas realizadas se
puede definir la verosimilitud de una hiptesis como una funcin proporcional al
producto de la distribucin de probabilidad de las medidas PO) por la distribucin de
probabilidad de los parinetros P(x), es decir.
CAP~TULO3: TCNICAS NUMRICAS DE OPTIMlZACldN PARA LA IDENTlFlCAClON DE
PARAMETROS
GEOTECNICOS
Lis) = k
3
1
exp[--(y -.t~')'((',)-'().-y')
\liqq
1 [t
exp - 7 ( x - x ' ) 1 ( C ' ~ ) - t ( x - x ' ) ]
4-1
o bien de minimizar la funcin soporte S(x).
Cuando se consideren fijas las estructuras de los errores de las medidas y de los
parinietros. al ser fijos tambien los ltimos tres sumandos de (3.6.). la expresin
anterior se reduce a la forma siguicnte,
Finalmente, si los errores de medida son independientes entre si y lo mismo sucede con
los errores de los parmetros, las matrices de covarianzas C .y C i sern diagonales y se
podr obtener la expresin de S(x) en funcin de las desviacrones tpicas de las medidas
y de los parimetros.
3.4 MINIMIZACINDE LA FUNCION OBJETIVO
3.4.2.1 Introduccin
Un mtodo de bsqueda directa es aquel que utiliza slo valores de la funcin objetivo
en distintos puntos para localizar el mnimo de la funcin. Hay varios mtodos clsicos
de este tipo, entre los que destacan: el de Hooke y Jeeves [1961], el mtodo del Simplex
[Nelder y Mead, 19651y el mtodo del Comples [Box 19651.
Los mtodos del Sir;.plex o del Complex pueden considerarse como mtodos volumtricos.
Por ejemplo. en d.el mtodo Simplex construye un tetraedn, que se va desplazando y
deformando (aumento o reduccin de volumen) hasta encerrar un punto ptimo en un
tetraedm de volumen ian pequeo como se quiera.
El mtodo de Simplex de Nelder y Mead es, sin duda el ms popular por su simplicidad.
su robustez y su eficacia prctica en problemas de tamao pequefo. Dicho mtodo ha
dado lugar a diversas modificaciones con el objetivo de mejorarlo [Parkinson et ai 1972;
Rikov. 1983; Torcmn. 1991; G m i a et al, 19921. Su principio es muy simple. En cada
iteracin, el \rrtiie de mayor valor asociado de la funcin objetivo se transforma, bien por
reflexin (simplex simetrico de volumen constante), bien por expansin (volumen
aumentado). bien por contraccin (reduccin de volumen). Si alguna de estas
CAP~TULO3: TCNICAS NUMRICAS DE OPTIMIZACI~NPARA LA IDENTIFICACION
DE
PARMETROS GEOTCNICOS
Otro mtodo que tambin emplea el Simplex, aunque basado en un principio diferente es
el mtodo baridntnco. En este caso es el baricentro del simplex el que se pretende que
mejore en cada iteracin. En cada etapa varios vrtices se transforman simultneamente
para desplazar el simplex y su baricentro [Enjalben et al, 19931. Se aplican etapas de
expansin, cor=tracciny reduccin pero siguiendo un esquema ligeramente diferente al
seguido en el algoritmo de Nelder y Mead. El objetivo perseguido por este nuevo mtodo
es niejorar la convergencia para los pnhlemas de gran dimensin (n > 10) en relacin con
el algoritmo de Nelder y Mead, para el cual las ventajas de clculo se degradan muy
rpidamente con la dimensin (n > 6).
- B) Se deterniina entre [x,. x,. x,. ....x,,,] el mximo valor de la funcin (/A). el
siguiente valor ms alto de la funcin (/A). el menor valor de la funcin (./;) y
los puntos donde estos valores se alcanzan. x,,. xF XI.
luego,
ii) Si j, 7.1; pero no se cumple quej, >fa. quiere decir que el valor de la
.
funcin en x, (.f;) es menor que los valores de la funcin en los dos puntos
mas desfavorables (/, <& <f;S, y por tanto se ha avanzado en el camino para
alcanzar el mnimo de la funcin. Se sustituye xn por x, se realiza el test de
convergencia (apartado 1) y si esta no se produce se vuelve al apartado (B).
es decir, x, + -1 (x, + x, ) .
2
Una vez obtenidos los nuevos puntos x, se calculan los valores def; para i = l .
2. ...(n+I) y se realiza el test de convergencia del siguiente apartado. Si no
converge se continua con el apartado (B).
- n+l
donde: 1 = (f i ( n + 1)
rrl I
T
Cijlculo de los n + l v&rticesdel
simplex iniaal
(xl. 12.x3. .
rn+l)
T
CBlculo del valor de la
func16n0bjeti~0(0 en cada
v&bce (fl. r;! f3. . fn+l)
T
b Determinar xl. xg. xh. fl. fg, fi
T
Calcular el cenlroide xo y fo
T No SI No
'frcn? , ~fr>fg? 'fr>fi? , Sustituir xh por
xr
SI Si
No
T h
Expansibny obtenci6n ContracUbny
de xe y fe obtencibn de xc y fc
-
T NO si
6 fe fl ? , T
Sustituir xh por xr Sustituir xh por xc 4
T
L fc < m ?
-
Si No
T T
Reduc66ndel sirnolex
Sustituir xh por xe
asignando a cada i,el
valor IxI+xIIR
v Evaluar
(fl. 12. B. . . fn+l)
Test de
Fin
Este mtodo es un mtodo iterativo que parte del conocimiento del nmero de
parmetros (n). del intervalo de variacin de cada uno de ellos -restricciones explcitas-.
del nmero de restricciones implicitas (m) y de un punto o vector de parametros de
partida (x,=x,) admisible. Por tanto el problema planteado ser:
min f ( x )
R, 0 ) 2 bl i = 1.2, ...m
1 5 x
1 1
Iii
1
j = 1.2. ...n
donde ir, y 1, son los valores superior e inferior respectivamente de cada una de las n
restriccionesexplicitas y los h, son las cotas superioresde las m condiciones implcitas.
v
Lectura de las restricciones
explicitas de cada padmetro
(0, uj,j= 1.. .n)
Inicializo el centroide xc=xl
- - ~
v
Conjunto de k=2n puntos
xij=~+j+mndom(ujJj)
LES
admisible xr? b Xr'(Xr+XCfl
- v
Actualizar xc.
v
Evaluar la funcion objetivo (0
en cada punto del complex. fi
-
v
Ordenar los puntos x i en
funcion de los valores fi
--
--Y- -
Comienzo proceso
iterativo
-
Cuando se dispone del complex y de los valores de la funcin objetivo en cada uno de
sus puntos se procede a ordenarlos en funcin del correspondiente valor de la funcin
objetivo.
l:na vez construido el primer complex comienza el proceso iterativo (Figura 3.8.) que
pretende dirigirse I~aciael punto donde la funcin objetivo olcanza un mnimo dentro de
la regin admiii:)le. El proceso que tiene lugar es el siguiente:
x, = (x, + x,;, f 2
iii) Se analiza si es vlido el nuevo punto. S i lo es se sigue c o n el paso D y si
no lo es se repite el paso C hasta encontrar un x, vlido.
v
Obtener el centroide xco de
los k-1 mejores puntos
v
Reflexin de xh en xco y
obtencion de xr y fr
Es NO Restriccin expliccta
admisible Si se violan las restncccones 11 o uj. A
xr? x@)=O+ lOexp(6J O x@)=uj-lOexp(-6)
v
Resificcin
impliuta
xr=(xc0+xry2
? No
fr fh?
v
Sustituir xh porxr
- v
Actualizar y ordenar los
d o r e s 11.12. .... Ik
v
No
Test de
convergenwa
(dm. a)
Fin
Eii la I'iyurli 3.9 puede observarse u11 organigrama del algoritmo de ajuste para el
modelo dc comporiamienio elastopliistico. En l se puede ver que el clculo numrico
de las trayectorias para cada uno de los nuevos vertices admitidos del complex -nuevos
iuegob de pariiietros- se realiza mediante el mtodo de elementos finitos adoptando
para el iiiatcrial el modelo dr. Drucker-13mger.
TCNICAS NUMBRICAS DE O~IMIZACION
CAP~TULOJ: PARA LA IDENTIFICACION DE
PARAMETROSGEOTCNICOS
IPaMmetros iniciales
C Generacibn o
modificacibn de los
"k" v6rtices del
complex admisible
Clculo de las
trayectorias descritas
por cada vertice del
complex mediante el 1
modelo de D.-P.
(F.E.M.) j
j
objetivo para cada v6itice
convergencia
1
Figttru 3. Y. Algoritmo de ajuste pura el modelo de Drucker-Prager (D.-P.).
3.4.2.5.1
. .
La tctica de este tipo de mtodos queda muy bien indicada -porla fbula de los tres
excursionistas [Charon et al, 19961. Dicha fbula comienza as. Erase una vez cuatro
excursionisras con poco sentido de la orientacin que se perdieron en la montaa en
CAP/WO 3: T~CNICAS
NUMRICAS DE OPTIMZACI&N PARA LA IDENTIFIcAcI~N DE
PARAMETROS GEOTCNICOS
cuatro lugares distintos. Siendo precavidos se han informado antes de salir y saben que
un equipo de rescate pasa con regularidad por el punto de menor altitud de la zona. Cada
uno dispona de un altmetro y una brjula aunque ninguno conocia la altitud del punto
ms bajo. A causa de la niebla su vi~ines limitada y no pueden ver a su alrededor ms
que el entorno inmediato. Por ello al llegar a un cruce slo pueden saber si un camino
sube o baja caminando sobre el. lo que les impide examinar a la vez los distintos
caminos que confluyen en el cruce. Cul fue la estrategia que sigui cada uno para
encontrar el lugar adecuado?
El primero. no sabiendo que direccin elegir. adopta el siguiente mtodo: en cada cruce
conocer la altitud a la que se encuentra y probar todos los caminos posibles.
examinando sucesi\mnente todas sus ramificaciones, y pudiendo volver hacia atrs para
explorar aquellas ramificacionesque todava no haba recorrido. Cuando haya terminado
de explorar todo. podr volver a la altitud ms baja. sabiendo con certeza que se
encuentra en el lugar adecuado para esperar el socorro. Sin embargo. esic excursionista
no Iiabia niedido bien la magnitud de su progecio pues si no ha muerto de cansancio.
comer continuamente por la montaa recomendo una de las numerosas ramificaciones
descubiertas despus de su salida.
El segundo prefiere otro mtodo: al llegar a cada cruce seguir el primer camino que se
presenie y que se dirija hacia abajo. Qu hay ms normal que descender para ir hacia
aba.jo!. Este excursionista tampoco tiene suene y el camino que ha tomado desciende
Iiacia una hondonada que consider el punto mas bajo. Estuvo esperando toda la noche
al equipo de socorro. Al amanecer pudo comprobar al desaparecer la niebla que se haba
equivocado. ponindose a caminar por fin haca el lugar adecuado que eslaba ms abajo.
..
3.4.2.5.2 Introduccion
Los denominados mtodos de recocido simulado son mtodos relativamente nuevos
[Kirkpatrick. 1984; Otten et al, 1989 Press et al, 1995; Charon et a1.19961 que se
emplean para resolver problemas donde es importante determinar los extremos globales
de la funcin objetivo. De hecho son mtodos que han demostrado su validez al haber
sido capaces de encontrar extremos globales en presencia de un gran nmero de
extremos locales no deseados [Press et al. 19951.
Desde un punto de vista prctico con esta merodologia se han resuelto problemas tpicos
como el problema del vendedor ambulante que consiste en encontrar el itinerario ms
corto que debe recorrer dicho vendedor pasando slo una vez por todas y cada una de
las vr ciudades donde quiere vender sus productos y retomar a la ciudad de origen. Este
es en realidad un problema de mi~mizacincombinatoria ya que el espacio en el que se
dcfine la funcin objetivo que se quiere minimizar no es realmente un espacio
n-diinensional de n variables o parhetros continuos sino que se trata mis bien de un
espacio discreto, el conjunto de las diferentes ordenaciones posibles de las n ciudades.
El nmero de elementos de este espacio es muy grande e indudablemente no es posible
considerar todos. Adems al tratarse de un espacio discreto deja de tener significado el
concepto de direccin de minimizacin que en un espacio continuo se emplea para
indicar la direccin que ms rpidamente conduce hacia la regin donde se encuentra el
mnimo.
3:
CAP~TULO TCNICAS
NUMRICAS DE OPTIMIZACIONPARA LA IDENTIFICACIN DE
PARAMETROS GEOT~CNICOS
El origen de la metodologa del recocido simulado hay qiie buscarlo en la analoga que
nrrsenta con fenmenos termodinmicos tales conio el vroceso de enfriamiento Y
posierior cristalizacin de los liquidas. o el camino que siguen los metales al enfriarse y
"recocerse".
I<ii general. aunque la analogia no es perfecta, se puede decir que la mayor parte de los
algoritnio, de minimizacin avanzan segn Ia direccin de enfriamiento mximo. Es
decir. a partir de un punto inicial buscan el camino ms corto y rpido en la direccin de
mxima minimizacin -en la direccin del mnimo-. lo que no siempre implica que se
alcance el mnimo global. sino que en muchos casos se obtienen minimos locales. En la
I'igura 3.10 se muestra cul es la esencia y el objetivo del recocido simulado.
3.4.2.5.3 5
. .. . .
A la hora de aplicar la metodologa del recocido simulado a la minimizacin de la
funciii objetivo resultante de la comparacin entre dos curvas de un ensayo
presiomtrico -una experimental y otra calculada numricamente- se considera que se
trabaja en un espacio continuo ndimensional. donde n representa el nmero de
parhetms del que depende la ley de wmportamiento del material.
En este caso el objetivo es buscar el mnimo global de una funcin objetivo V(x)) del
tipo (3.6.) que se ha definido al w m i e m del presente capitulo. Este mnimo se
alcanzar con un determinado juego de paremetros (x) que representa el estado del
sistema citado por Metmpolis. El partimetro de contml(7') sigue siendo algo parecido a
la temperatura que disminuye gradualmente siguiendo las indicaciones del protocolo de
recocido. El generador aleatorio de estados del sistema o de juegos de parmetros
CAP~WL3
O: TkCNICAS NUMRICAS DE 0PTIMIZACI~NPARA LA IDENTIFICACI~NDE
PARAMETROS GEOT~~CNICOS
El mtodo de recocido simclado al que aqu se alude [Press et al. 19951 esta basado en
una variante del mtodo de bsqueda directa Simplex [Nelder y Mead, 19651. ya
.I
estudiado en el apartado 3.4.2.3. Es decir. cada estado del sistema o conjunto de
partimetros 0)da lugar un simplex que est formado por n+l vrtices. Este simplex
mediante movimientos de reflexin. contraccin y expansin va avanzando hacia el
estado de mnima energia -conjunto de parmetros x que minimizan la funcin objetivo-
Pan un valor cualquiera de T,el simplex se expande hasta alcanzar una escala prxima
al tamaiio de la regin del espacio que puede ser "abarcada por el valor del parmetro
de control. Cuanto mayor sea T mayor ser la regin abarcada por el simplex en cada
caso. Dentro de esa regin se produce un movimiento aleatorio Brosniano que permite
ensayar nuevos puntos elegidos aleatoriamente dentro de ella. De este modo la
eficiencia con la que se explora unz regin es independiente de su estrechez -para un
valle elipsoidal. del mtio de sus ejes principales- y de su orientacin. Si el enfriamiento
( T ) se lleva a cabo lentamente la eficiencia ser tan buena que el simplex se ir
contrayendo hacia la regin que contiene el mnimo global buscado.
CAP~TULO3: TCNICAS NUMRICAS DE GPTIMIZACI~N PARA LA IDENTIFICACI~N DE
PARAMETROSGEOT~CNICOS
-v -
Cdlailo de los n+l vbltices del
simplex incaal
xl. x2, x3. .
m+/)
T
C&iailo del valor de la funa6n objetivo (0
en cada vebce (fl. f2. 13. fn+l)
v
Fluduao6n Irm~ca aleatona
Al valor de f en cada Verlice se le afiade una vanable aleatona
poslbva de distrcbua6n logaritmccaproporaonal a T
f = f + T log(random)
T
h Deternilnar xl. xg. x h n, fg. fh
v
Reduca6nde T
v
Reflexi6nde xh por un f a a ~
-1 y obtenabn de xr y fr
Fluciuaa6ntermcca aleatona
fr = fr - T log(random)
V No S, - - - -
Contracci6ncon facior
Lfrcn7 b fr>fg7 b0.5yobbnahdexcvfc
Si
-
fe = fc T bg(randM)
V 140
Expansioncon facior 2 y
v
obtena6n de xe y fe Suatulr xh por xc 4 fc c fh 7
fe = fe - Tbg(random) Si
-
T
No
T v Redwrr6n del s~mplex
fe fJ7 b Sustituir xh por u alrededor de xl
SI NO M
T
Suslituir xh por xe
M v -
V Evaluar
fl. R.f3, fn+l)
Test de
No convergema
3.4.2..Algoritmos genticos
La idea que subyace en la metodologia de los algoritmos genticos es la de la evolucin de
los seres vivos: los mejor adaptados tienen mas posibilidades de sobrevivir y por tanto de
reproducirse. Si se pone el ejemplo [Michalewicz 19961 de una poblacin de conejos, se
puede decir que dentro de ella habr conejos que sean mas rpidos y astutos. Otros en
cambio ser& ms lentos. Estos ltimos sern una presa mas fcil para sus depredadores,
mientras que los primeros tendrn nias posibilidades de vivir y de reproducirse. La
poblacin de conejos que ha sobrevivido durante un periodo de tiempo comenzar a
aparearse y --producirse. producind~seuna mezcla del material genetico: algunos conejos
lentos se aparearn con conejos rpidos, ot:os con conejos astutos y otros entre si; lo
mismo suceder con los conejos rpidos y con los astutos. Puede que incluso tenga lugar
una mutacin gentica provocada por la naturaleza. Se puede decir que las cras de conejo
surgidas de este proceso sern mas astutas y rpidas que las de la generacin anterior por
Iiaber sufrido menor depredacin los conejos dpidos y astutos.
Las entidades responsab!es de la produccin de los caracteres hereditarios son los genes y
el conjunto de genes de un individuo deline su gcnotipo. IJn gen es un segmento de
cromosoma -filaniento de ADN- que es el material genetico de todas las clulas. Los
cromosomas se situan en el nuclco de las clulas. Cada clula de cualquier organismo de
una especie tiene un nmero determinado de cromosonias pxsentados por pares. En cada
par uno de los cromosomas proviene del padre y el otro de la madre. Los cromosomas de
un mismo par se denominan homlogos.
)os mecanismos permiten fabricar clulas nuevas [Charon et al, 19961. El primero tiene
lugar por la ditisin celular: un3 clula duplica su material gentico antes de dividirse en
dos copias iguales a casi iguales -si se pmduce un error que afecta a un gen- en cuyo caso
se tratara de uiia mutacin -de ese gen-. Este mecanismo tiene lugar durante la
reproduccin asexual y se produce rara y aleatoriamente. El segundo es el mecanismo de la
reproduccion sexual o procreacin cuyo fruto es un nuevo individuo a partir de sus padres.
I'.n este mecanismo las clulas sexuales transmitidas por los padres aportan. cada una. la
mitad de los cromosomas del nuevo individuo lo que le proporciona e1 numero total de
cromosomas.
Los algoritmos genticos se inspiran mas o menos fielmente en estos mecanismos. de tal
modo que la solucin al problema planteado de minimizacin de una funcin objetivo
vendr dada por una determinada cadena de cromosomas con su ordenacin interna
camcteristica. Para alcanzar dicha solucin se lleva a cabo un proceso en el que se
mantiene una poblacin de soluciones potenciales. Es importante que este mtodo de
bsqueda tenga en cuenta por un lado explorar todo el espacio de bsqueda y por otro
aprovechar y explotar las regiones donde las soluciones parecen mas prometedoras.
Por ultimo decir que los algoritmos genticos se han usado can resultados satisfactorios
en problemas de control adaptativo. en el problemn del vendedor ambulante. en
problemas de control ptimo. elc [Grefenstette, 1985; Davis, 1987; Michalewict, 1990;
Forrest, 1993; Eshelman, 1995; Michalewicz, 1996 J.
CAP~TULO3: TECNICAS NUMBRICAS DE OPTIMIZACI~N
PARA LA IDENTIFICACION DE
PARAMETROS GEOTECNICOS
3.4.3 Mtodos de -b
. .
3.4.3.1 Introduccin
La clasificacin de los llamados mtodos "clisicos" o de bsqueda indirecta depende del tipo
de informacin que precisen de la funcin objetivo. As. se tienen por orden decreciente de
informacin que se precisa los mtodos de Newton, los mtodos de cuasi-Newton y los
mtodos del bmdiente. El mtodo de Newton requiere informacin del gradiente y de la
matriz hessiana de la funcin objetivo al realizar en cxda iteracin una aproximacin de
segundo orden de dicha funcin. En general la matriz hessiana es muy costosa de evaluar y
adems debe ser definida positiva para asegurar la convergencia. Los mtodos de cuasi-
Newton [Zijlsua, 19961 precisan informacin del gradiente de la tiuicin objetivo y
aproximan de forma muy simple la matriz hessiana mantenindola simtrica definida
positiva.
A continuacin se hace una breve referencia i estos mtodos, pues uno de ellos. el mtodo de
Marquardt (1 9631. es el que se ha empleado en la minimizacin de la funcin objetivo para el
niodelo de material dado por la ecuacin analtica hiperblica de Ferreira.
minf (x)
con,
m
f (x) = r'r = C ri2
i=l
3: TCNICAS NUMBRJCAS DE OPTIMIZACIN PARA LA IDENTIFICACIN DF
CAP~TULO
PARAMETROSGEOTCNICOS
Cada residuo es igual a la diferencia entre el valor de la variable observada (y' ) y el valor
de la variable calculado (y). La estimacin por mnimos cuadrados de los parmetros x
debe verificar que el gradiente d e n ) sea nulo Cg(x)=O) [Norton, 19881.
donde J , es la derivada parcial del residuo r, con respecto al parmetro x,. J se denomina
matnz jacobiana de la funcin,f(x) respecto a ios parmetros x y est formada por ni filas -
numero de datos u observaciones- y n coliimnas -nmero de parmetros x,.
Para localizar las races de las ecuaciones g(x)=O se scirle aplicar el mtodo de Nerrzon.
Con dicho mtodo el nuevo vector de parmetros x 'se obtiene del siguiente modo.
donde g '(x) es la ma'riz hessiana H(xj cuyos elementos vienen dados por la expresin,
Como se puede apreciar la matriz hessiana viene dada por la suma de dos matrices. El
primer sumando se obtiene a partir de la informacin de las primeras derivadas y el
segundo a partir de la informacin de los residuos y de las derivadas segundas de los
mismos.
Como se puede apreciar los mtodos de Gauss-Newton son muy parecidos a los mtodos
de Ne~vtondiferencindose nicam :nte en el clculo de la matriz hessiana. En el mtodo
de Newton sta se calcula exactamente mientras que en los mtodos de Gauss-Newton se
emplean aproximacionesde la misma.
m
t > i i n f ( x )= 1r,2 ( x ) = r 7'r
i=t
G.iO.)
Una condicin necesaria y suficiente para encontrar el mnimo de esta funcin es que el
gndiente se anule k(x)=U) y el hessiano sea definido positivo wash, 19871. Aplicando
en (3.10.) el mtodo de Newton para determinar las mices del sistema de ecuaciones
(3.7.).
J'JI = - ~ ' r
(3.1 l.)
Marquardt [1963: sugera que K fuera la matriz identidad o bien la matnz formada por
los elementos de la diagonal de H'. En este trabajo, siguiendo las indicaciones de Nash
11977 y 19791 se ha considerado que Kes una matriz diagonal cuyos elementos son los
elementos de la diagonal de H' ms la matriz identidad l. es decir.
Esta aproximacin presenta las ventajas de tener las propiedades de escalamiento de los
elementos de la diagonal de H' cuando estos son grandes, y de evitar las posibles
singularidadesque puedan aparecer en la resolucin de las ecuaciones de Gauss-Newton
~AP/TULO3 . TCNICAS NUMERICAS DE OPTIMIZACI~N
PARA LA IDENTIFICACI~NDE
I*ARAMETROS
GEOT~CNICOS
cuando algunos elementos de la diagonal de H' se hacen cero como consecuencia del
clculo computacional.
En esta tesis las ecuaciones de Gauss-Newton en el contexto de Marquardt se han
resuelto mediante un algoritmo de Cholesky que permite fiiiar cienos parmetros y
establecer intenralos vlidos para cada uno de ellos. Al igual que se apunt
anteriormente, estas ecuaciones permiten determinar la direccin de mxima
minimizacin (1). que a su vez proporciona el nuevo juego de parmetros x ',. Con stos
se calcula finalmente el nuevo valor de la funcin objetivoflx 1.
3.4.3.4.1 L l e f i n i c i d n e t r o de -
M
La estrategia seguida para definir lamda (A) mash. 19871 depende si el valor de la
funcin con los ndevos parmetros flx] es menor o mayor que su valor con los
parnietros viejos f(x). El proceso comienzz asignando a h un valor arbitrario en la
primera iteracin. en particular en este trabajo se ha asignado el valor de h = 0.00004.
Si fi.r '/>f(x) o la matriz hessiana resulta ser singular entonces se aumenta el valor d e l
multiplicndolo por 10. Por el contrario. si f(x ]<flx) entonces . se multiplica por U. 1.
Para evitar que se viole alguna restriccin (N(Nd987). es decir, que aparezcan valores
de los parmetros fuera de sus intervalos de variacin, y asumiendo que la direccin de
bsqueda (11 es aquella en la que el tamao de paso o incremento (SS) es positivo, se
puede detenninar cul es el incremento mximo que a partir del parmetro x,. a lo largo
de la direccin I,, permite encontrar un nuevo parmetro que no viole las resiricciones.
Es decir. si:
CAPITULO 3: T ~ C N I C A NUMRICAS
S DE OPTIMIZACI~NPARA LA IDENTlFlCACldN DE
PARAMETROS GEOT~CNICOS
Por otro lado se sabe qi:e la direccin de bsqueda que es igual a la contrara del
gradiente de una funcin, indica la direccin en la que la funcin se reduce ms
rpidamente (I=-~(x)). Por tanto. cuando est activa la condicin superior (x,,,) la
correspor~dientecomponente del gradiente es negativa (g, < O). y cuando est activa la
condicin inferior (x,,,,~)la correspondiente componente del gradiente es ,positiva
(g, > O). En la Figura 3.12 se puede apreciar el signo del gradiente dependiendo de qui.
condicin este activa.
l I . ,I
\ ; - 1
\; /, i
<. '
I I
, - \ \<.
3.4.3.4.4 p
Como se acaba de mencionar, cuando el parmetro, tiende a cero la direccin de
biisqueda del mnimo tiende a la de Gauss-Newton. Sin embargo aunque el algoritmo de
Maquardt es de convergencia m& lenta que el de Gauss-Newton proporciona una
seguridad mayor de convergencia. Es decir. cuando el vector inicial o puntb Iie partida
se encuentra lejos del mnimo. el algoritmo de Gauss-Newton puede dar lugar U
incrementos demasiado grandes que conduzcan a regiones con valores superiores de la
funcn objetivo. En esos casos la limitacin automdtica del tamao de paso (SS) y la
~nodificacinen la direccin de bsqueda ( 1 ) son propiedades muy tiles del algoritmo
de Marquardt. En la Figura 3.13 se presenta un algoritmo completo del mtodo de
minimizacin de Marquardt.
Por otro lado el algoritmo de Marquardt no requiere el clculo de las derivadas segundas
que suele ser siempre muy costoso. y la aproximacin de la matriz hessiana est mejor
condicionada que la empleada por el mtodo de Gauss-Newion.
Nmero de pa&netros (n)
Y nmero de restricciones
impllcilas (m)
7
Punto o vector inicial x
Calcular el valor de la
funcin objetivo f(x)
7
Valor iniual de ).
'v
Calcular la matnz J ~ J
Y el segundo miembro del
slslema de ecuauones JF
7
Completar la diagonal de J,J
(JTJ) + AI(JTJ)B + q
7
Resolverel sistema
de Gauss-Newton
por el m. de Cholesky
v
Obtener direccin de
minimizaciont.
t=-g
T
S%= 1
'
SS SS0 b Reductr por restncciones SS
Si No
7
SS = S% , 7
Nuevos PaIametros
d=x+sst
SI 7
Aumentar
Test de
3. 4 'f*) 'f(x) 4 Catcutarf(x> 4 conwQenoa
No No
7 Si
Disminuir v
?.
Fin
Cuando se trabaja con un modelo lineal. el residuo viene dado por la expresiii.
donde la variania de los errores de medida (S') se representa por la estimacin sigtiiente.
donde m es el nmero de datos u observaciones y n el nmero de par.imetros del
modelo.
Suponiendo que el residuo fuera nulo. la estimacion de los parmetros x viene dada por:
donde J es la matriz jawbiana y S' viene dada por $3.12.) con residuos r no lineales. En
el caso de existir parmetros fijos la expresin de S- es la siguiente,
donde la mariz Q(x) es una matriz dagonal que tiene como eleiiientos el recproco de la
raz cuadrada de los elementos de la diagonal de I/(x).
3: TCNICAS NUM~RICASDE OPTIMIZACINPARA LA IDENTIFICACINDE
CAP~TULO
PARAMETROS GEOTCNICOS
Estas son algunas de las aproximaciones que existen para estimar cuhles son los
intervalos de confianza de cada uno de los parmetros, es decir, regiones del espacio de
los parmetros donde estos se pueden encontrar con una determinada probabilidad. La
linealizacin de problemas no lineales proporciona medidas de la dispersin de los
parmetros de las que se pueden deducir los intervalos de confiaw correspondientes.
3.4.4.2.1 Radio de c u r v u
Una de las posibilidades para conocer una estimacin de la dispersin de cada
parmetro. cuando no se conoce la estructura de la funcin objetivo empleada en la
estimacin. es la busqueda a lo largo del eje de cada parmetro para determinar el radio
de curvatura y el ngulo 0 (Figura 3.14). Esta mtodo proporciona una medida de la
dispersin bastante tosca.
en la que ic(i)es un incremento para cada parhmetro x, dado por una expresin del
tipo.
El modo concreto de llevar a cabo esta bsqueda axial es definir en el punto de soliicin.
mediante un modelo cuadrtico -dos veces diferenciahle-. una superficie que da lugar a
una parbola a lo largo del eje de cada parhetro. Llamandof(U), f(-). f(+) a los valores
de la hincibn en Los tres puntos considerados para el eje de cada parametro se tiene,
como.
Para calcular el ngulo de inclinacin O en grados basta con aplicar la relacin siguiente
O = -arrtung(,f') = - u r c ~ u n ~ ( h )de
. tal modo que en las proximidades de las
restricciones superiores e inferiores O estar respectivamente prximo a -90 y 90 grados.
Por otro lado cuando la funcin objetivo es simtrica respecto a la solucin estimada se
tendra que h=O.
Por ltimo, decir que la raiz cuadrada de cada una de las componentes de V(x) podr ser
comparada de algn modo con la desviacin tpica correspondiente a cada parametro
obtenida con un anlisis de incertidumbre de un mtodo indirecto.
La determinacin del mnimo de una funcin plantea una serie de problemas que se
deben tener en cuenta, y que van ms all de la aplicacin simple de un algoritmo de
minimizacin. As por ejemplo desde un punto de vista terico pueden aparecer
problemas de unicidad, de existencia y de estabilidad.
La unicidad tiene que ver con el hecho de que puede haber ms de un vector de
parametros que permita obtener el mnimo global de la funcin objetivo. La existencia
se refiere a que no tiene por qu haber necesariamente un conjunto de parmetros que
minimicen la funcin. y el problema de la estabilidad se plantea cuando, suponiendo que
se cumplan los postulados de unicidad y existencia. pequeas variacones en los datos
observados dieran lugar a parmetros muy diferentes.
Una vez obtenidas dichas derivadas se han programado e implementado dentro del
programa de Marquardt, cuyo algoritmo se encuentra en la Figura 3.13. El test de
convergencia que emplea este programa es un test de igualdad en funcin de la precisin
de la mquina (l..!?). es decir. cuando el vector de parhetros nuevos (a') es "igual" al
vector ptimo anterior (a) el programa se detiene al no poder encontrar un valor ms
pequefio de la funcin objetivo.
,
CAPITULO 4
Todos ejemplos se ha tratado con los diferentes algoritmos presentados, y con los dos
modelos de comportamiento elegidos, el modelos elstico no lineal de Ferrein y el
modelo elastoplstico de Drucker-Prager. En el caso de la modelizacin de Drucker-
Prager se ha analizado el efecto de los diferentes tipos de mallado y de elementos finitos
utilizados.
Adems. con el fin de estudiar la influencia de los errores en las medidas. a los datos
iniciales de la presin del ensayo ficticio se ha sumado una cantidad aleatona. Dicha
cantidad resulta del producto de un nmero aleatorio comprendido entre O y 1 por la
precisin del presimetro respecto la presin. que es de t 10 Kpa. La variabiliu.id de los
parmetros obtenidos tras el proceso de minimizacin han permitido analizar la
influencia de dichos errores de medida en la estimacin paramtrica. Por otro lado, se ha
estudiado la influencia de los valores iniciales de los parmetros y del intervalo de
variacin -restricciones explcitas- de cada uno de ellos en el proceso itemtivo de
solucin del problema inverso planteado.
Hay que seaalar que se ha podido realkm un esnidio ms completo del comportamiento
de los diferentes a l-m i m o s de minimizacin de la funcin obietivo con el modelo de
comportamiento de Ferreira del ensayo presiombico. debido a su carcter de solucin
CAP~TULO
4: APLICACI~N
A WEMPLOS SINTTICOS
analtica y a la facilidad y rapidez con la que se obtienen nuevos partimetros del medio.
En el caso de modelizar el comportamiento de las arenas con una ley como la de
Druker-Prager, los clculos son m b tediosos y el tiempo de ejecucin requerido mucho
mayor. Este hecho no afecta tanto a la metodologa de caracterizacin desarrollada
(Figura 4.1). que funciona con independencia del modelo considerado, como a la
modelizacin del comportamiento del material.
Juego de
parametros iniciales
Informacion
adicional sobre
Calculo de la funcion los padmetros
objetivo
4
Datos
experimentales
NO T
4 C. Convergencia
Si
T
Fin
Por liimo comentar que todos los clculos que se presentan en este capitulo y en el
siguiente se han realizado en un PC con un procesador Pentium a 166 MHz.
o', = 21 KPa
(I / S ) , = 058
S, = lo0
CAP/WLO 4: A?LICACI~NA EJEMPLOS SINTETICOS
Ejemplo sint6tico
P
= la,
P,
o .-
Fiprru 4.7. Duros iniciales del enw+w con arenujicticia (modelo de Ferreira}.
4.3 MODELO DE COMPORTAMIENTO ELASTICO NO LINEAL
A continuacin se describen los distintos ensayos que se han realizado en funcin del
tipo de algoritmo empleado en el proceso de minimizacin de la funcin obietivo. Para
un mismo algoritmo se han llevado a cabo experimentos diferentes. escalando o no los
parmetros. introduciendo los posibles errores en las medidas y considerando diferentes
valores iniciales de los parametros en intervalos de variacin suficientemente amplios.
En algunos de los algoritmos que se citan a continuacin se han podido delimitar lo3
intervalos de variacin fsica posible para cada parmetro. En concreto se han impuesto
restricciones explcitas de los parmetros en los algoritmos de Marquardt y Complex.
niientras que se ha trabajado sin restricciones con el algoritmo de Recocido Simulado.
Los intervalos dc variacin considerados para cada parmetro se presentan en la Tabla
4.1.
("=)ult o I
Si SO 200
Tuhlu J. l . Inrrr~~ulos
de ~ririuci(inde los purcn~e~ro-Y
del nlodrlo ile Ferreiru
Con este algoritmo se han llevado a cabo mltiples anlisis inversos del ejemplo
sinttico tratando de estudiar si se puede reproducir la curva inicial de carga con
independencia del punto de partida. De todos esos anlisis nicamente se va a hacer
referencia en esta memoria a los que desde el punto de vista de los resultados obtenidos
se creen m& interesantes.
CAP~TULO
4: APLICACIN
A EJEMPLOS SMTlfTlCoS
Las ciiwas de carga correspondientes a los valores ptimos de los parmetros de los
cinco experimentos se representan en la Figura 4.3. En ella se puede apreciar que
existen dificultades en los tres ultimos casos de la Tabla 4.2.
- .
- - -- - -- - - - - -
En la Tabla 4.4 se encuentran los parmetros iniciales de los tres ltimos experimentos
(E.F.M.3, E.F.M.4 y E.F.M.5). Al trabajar con parmetros escalados se aade al nombre
de la simulacin las letras ES (E.F.M.ES.3, E.F.M.ES.4 y E.F.M.ES.5).
C A P ~ J L O 4: APLICACIN A EJEMPLOS SMTTICOS
MX)
En la Tabla 4.5 se puede ver como el valor de la funcin objetivo para esta simulacin
(E.F.M.ES.4) se reduce muchisimo y como los parmetros obtenidos se aproximan
bastante ms a los parmetros del ensayo ficticio que los obtenidos sin escalar los
parmetros (E.F.M 4).
Tabla 4.5. Resulrados del experimento E.F. M.4 con el algoritmo de Marqirard~con
parmelros escaladosy no escalados.
Tambin se aprecia que el tiempo y el coste para obtener esta solucin ha sido muy
elevado pues se han realizado 13465 iteraciones y 18833 evaluaciones de la funcin
objetivo. Por este motivo se han aplicado otros mtodos de optimizacin con el fin de
estudiar si es mas ventajoso su comportamiento.
Se han llevado a cabo los mismos experimentos que antes, cuyos parmetros iniciales
quedan reflejados en la Tabla 4.2, es decir en todas ellas se han considerado parmetros
sin escaiar. Los resultados obtenidos se muestran en la Figura 4.6 y en la Tabla 4.6. En
ellas se ha aiiadido una 'R' a la denominacin de los ensayos para hacer notar que los
datos iniciales de partida llevaban incorporado ruido. En la Figura 4.6 se observa que
nicamente E.F.M.R.1 y E.F.M.R.2 son capaces de reproducir la curva experimental.
Con el resto de los ensayos sucede lo que ya ocurria cuando no haba error en los datos
iniciales (Figura 4.3).
En la Tabla 4.6 figuran los resultados de cada uno de los experimentos. Se puede
apreciar que el valor de la funcin objetivo en E.F.M.R.1 y E.F.M.R.2, an siendo
considerablemente menor que el correspondiente a E.F.M.R.3, E.F.M.R.4 y E.F.M.R.5,
es mayor que el alcanzado al partir de datos sin ruido (Tabla 4.3). Tambin son mayores
el iimero de iteraciones necesario para la convergencia del proceso, el nmero de
evaluaciones de la funcin objetivo y el tiempo empleado. ~ o s ~ ~ a h n e tresultantes
ros
del proceso de minimizacin partiendo con datos con ruido se encuentran, en todos los
casos, mas alejados de los originales que cuando se parte de datos sin mido (Tabla 4.6).
+E.FmCo +
Rudo
4E.FMR.1
4 E 5 M R . 2
+E.FMR.3
-8-E.FMR.4
a EFMR.5
-. - - -
Con el fin de conocer c6mo influye el emor en las medidas sobre los parmetros optimos
finales que se obtienen tras el proceso de miriimizacion con el algoritmo de Marquardt
se han realizado tambin los cinco experimentos de la Tabla 4.2 introduciendo un error
aleatorio en las medidas de presion de I 20 KPa y I 30 KPa . Los resultados obtenidos
se encuentran en las Tablas 4.7 y 4.8 y en las Figuras 4.7 y 4.8.
Figura 1.7. Optimizacin a partir del algoritmo de Marquardt con ruido (& 20 W a ) .
Tabla 4.7. Resuliados de la optimizacin con el algori~mode Marquardi
con ruido (% 20 KPa).
-- -- --
- 600
+CFMR 1
EFMR 2
u - -EFMR 3
-EFMR 4
t E E U R 6
--
0.00 200 4.00 6.00 8.00 10 00
Deformaoldn de la Cavidad 1%)
Figura 4.8. Opfimzacinu pariir del algoritmo de Marqunrdr con ruido (i30 KPa).
CAP~TULO4: APLICACI~NA EJEMPLOS SINTTICOS
Como puede apreciarse tanto en las Figuras 4.6,4.7 y 4.8 como en las Tablas 4.6,4.7 y
4.8 el aumento de ruido en los datos experimentales apenas influye en los parmetros
ptimos que se obtienen en cada uno de los cinco experimentos analizados. Esto
demuestra que el algo-itmo de Marquardt es bastante estable. La tendencia en todos los
experimentos es similar y se mantienen los problemas con los experimentos E.F.M.R.3.
E.F.M.R.4 y E.F.M.R.S.
-0-E Frno
-eE F r n o +
Rudo
-+EFCR 1
+EFCR 2
D
i? +EFc ;
200
&EFCR 4
e
+EFCR 5
A cantinu~cion.en las Figuras 4.1 1 y 4.12 y en las Tablas 4.12 y 4.1 3 se presentan los
resultados obtenidos al considerar errores aleatonos en las presiones experimentales de
f 20KPayf 30KPa.
CAPITULO4: APLlCACldN A EJEMPLOS SINTETICOS
-m- E. FTDSD +
-e E.F.C.R. 1
E F.C.R. 2
-E.F.CR. 3
+EF.CR 4
EFCR 5
. - -- .-- - -
1
~
+E.F&zD +
Rudo lI30Kpal
+E.F.CR. 1
E.F.CR. 2
-s- EF.C.R. 3
+E I . C R 4
E.F.C.R. 5
Fig~ru4. I2. Op~imizacionu partir del al~oritmoComplex cun ruido (& 30 KPa).
Obsewando las Tablas 4.1 1.4.12 y 4.13 junto con las Figuras 4.10.4.1 1 y 4.12 se puede
afirmar que el ruido creciente no impide que en todos los casos se puedan reproducir las
curvas experimentales de la Tabla 4.9. Por otro lado, al variar la magnitud del error los
parheiros ptimos que se obtienen en cada experimento no guardan relacidn aparente.
Es decir, en cada experimento, con distintos niveles de mido, se consiguen buenos
ajustes aunque con parimetros diferentes.
4.33.1 Pariimetros sin escalar
Con este mtodo directo de optimizacion se han repetido todos los experimentos
llevados a cabo con los otros dos metodos vistos hasta ahora. En este apartado se
presentan iinicamente los resultados correspondientes a las casos que Iian dado ms
problemas (E.F.3, E.F.4 y E.F.5). Los parmetros de partida de cada uno de estos
experimentosse encuentran en la Tabla 4.14.
datos experimentales sin ruido (Tabla 4.15). El valor de la funcin objetivo es, al igual
que suceda en los casos anteriores w n ruido, superior al obtenido con los datos
experimentales limpios. Lo mismo sucede en general, con el nmero de iteraciones
(Iter) necesario para que se reduzca el parmetro de control T. con el nmero de
evaluaciones de la funcin objetivo y. por tanto, con el tiempo consumido en la
resolucin del problema inverso.
- -
- -
600 -
CC
Y +E Frncm +
!m
- -e E F A R 3
2 m3
B
a
-2 200
C
+E FAR S
e 100
P
O.00 2.00 4.00 6.00 8.00 10.00
Oeformaci6n de la Cavidad 1%)
-- ---- - - -.- - - -
OhsendolasTablas4.l6.4.17y4.l8junioconlasF~4.I4.4.15y4.16spmdc
afirmar que el ido craaaitc ao mpide que eo kis ors casos andiadm se obvagim
ajustes de las o u r v a ~ de La Tabla 4.14. Tambin se aprrcia que al
aumcntarelnivcldeanrrd~KgCI.da,bdxnrfiJdo.9i?akja&uczml
del valor i n i d G m 2827). Par ~ mM, o se quc aunque mrho br
parSmermaBptimoairaksalvahrlamsgniadddamr.oa~a~bimammniiarr
que las que se han pmucidoooa d dgormao Cooqpkn
4.4 MODELO DE COMPORTAMIENTO ELASTO-PLASTICO
- Tensin de fluencia K
En la Figura 4.17 puede observarse como para el mismo juego de parmerros inicial del
modelo de Drucker-Prager la modelizacin de 465 elementos tiene un comportamiento
ms rgido que la que emplea unicamente una fila de 17 elementos. es decir, esta ultima
presenta mayores deformaciones para un inismo nivel de cxga. Tambin puede
apreciarse que este fenbmeno resulta ms significativo a medida que crece la presin
aplicada.
- -
- O - --- - -- -
Figura 4.1 7. Dalos iniciales del ensayo con arenaficticia (modelo de Drucker-Prager).
A EJEMPLOS SWITTICOS
4: APLICACI~N
CAP~TULO
Las restricciones explcitas para cada uno de los parmetros o intervalos de variacin
posible se indican en la Tabla 4.20.
Las curvas resultantes de los cinco experimentos llevados a cabo con el algoritmo
Complei: se representan en la Figura 4.18. En elC puede apreciarse el buen ajuste
alcanzado en todos los casos (Tabla 4.21). Sin embargo se producen ciertas variaciones
en los parmetros finales del modelo lo que refleja que la solucibn del problema inverso
puede no ser nica.
- - - . -
m - -- -
*EFC 4
5
g
.- 100 - - - - - -
-- - -
Por otro lado, la incertidumbre de los parrnetros obtenidos en las simulaciones que han
empleado un m6todo directo de minimizacin de la funcin objetivo -mtodo Complex
y Recocido Simulado- viene dada por la varianza que se ha determinado por el mtodo
de jachknife. Este in6todo (apartado 3.4.4.3) lleva a cabo tantos anlisis de la misma
simulacin como datos iniciales mas uno tenga sta. Por este motivo unicamente se ha
aplicado este mtodo para calcular la incertidumbre de los parmetros del modelo de
Ferreira, ya que aplicarlo a los parmetros del modelo de Drucker-Prager supondra un
gai J computacional muy elevado.
En la Tabla 4.22 quedan reflejados los resultados del aniisis de incertidumbre llevado a
cabo con los experimentos en los que se ha empleado el mtodo de minimizacin de
Marquardt (Tabla 4.3). Puede observarse (Tabla 4.22) que los valores mas bajos de la
desviacin tipica se corresponden w n las simulaciones E.F.M.1 y E.F.M.2que son las
que obtenan los parmetros ms prximos, prcticamente identicos. a los del ensayo
ficticio (Tabla 4.3). En estos dos casos se dispone tambin de informacin del radio de
curvatura p y del ngulo de inclinacin 8. Cuanto menor es el radio de curvatura el valle
ser mas estrecho y viceversa. La desviacin tpica de K,eY es nula pues se ha supuesto
que este parmetro permanece fijo.
CAP/TULO 4: APLIcACIN A EJEMPLOS S I N T ~ C O S
En la Tabla 4.23 se muestran los resultados del anlisis de incertidumbre del ensayo
E.F.M.ES.4(Tabla 4.5). Se puede coniprobar tambin aqui que el parmetro (I/s),,, es el
que menor desviacibn tpica tiene. La desviacin tpica de K,ct' es nula pues se ha
supuesto que este parmetro permanece fijo en e s e caso.
En l a Tablas 4.24 y 4.25 se recogen los resultados del anlisis de incertidumbre llevado
a cabo con el mdtodo de jackknife sobre simulaciones en las que se ha minimizado la
funcin objetivo con un mdtodo directo.
experimentos sucede que las desviaciones tpicas de cada uno de los parmetros sean
menores que en el resto de los casos.
4.6 CONCLUSIONES
los datos originales. Al introducir un error en los datos de presin de partida los
resultados obtenidos -sin escalar los parmetros- son semejantes a los que se acaban de
citar: nicamente se han conseguido ajustes buenos con las simulaciones E.F.M.R.1 y
E.F.M.R.2 (Tabla 4.6) aunque con un coste computacional y un valor final de la funcin
objetivo es mayor. Al aumentar c'! nivel de error se comprueba que el algoritmo de
Marquardt es un algoritmo bastante estable.
desviacin tipica de cada parmetro es muy bajo. El paremetro del modelo de Femeira
que ha dado un mayor nivel de confianza con el algoritmo de Marquardt ha sido (th),,,..
-
APLICACION A EJEMPLOS REALES
CAP~TULO5: APLICACI~NA EJEMPLOS REALES
CAPITULO
5: APLICACINA EJEMPLOS REALES
Las muestras de arena sobre las que se ha realizado el ensayo presiomtrico se han
formado bajo distintas condiciones y presentan diferentes propiedades. Son muestras
normalmente consolidadas o sobreconsolidadas valiendo K, un valor particular en cada
caso, en las que una vez alcanzado el nivel de carga prefijado se imponen las
condiciones de c,.ntomo apropiadas.
Las tcnicas empleadas para la minimizacin de la funcin objetivo son las que se
describen en el tercer capitulo de esta tesis, es decir. el mtodo de Marquardt. el mtodo
Complex y el mtodo de Recocido Simulado. Los resultados obtenidos con cada una de
las tcnicas de identificacin de parmetros se comparan entre s y, en los casos que es
posible. con los resultados obtenidos por otros autores referidos a los mismos ensayos.
Tambin se lleva a cabo un estudio de la incertidumbre de los parhetros obtenidos en
cada caso con el fin de poder saber el grado de fiabilidad de cada uno de ellos.
Se presentan a continuacin los resultados obtenidos para cada uno de los tres ensayos
citados.
Este ensayo, al igual que los otros dos que se estudian en este captulo, es un ensayo
presiomtrico drenado llevado a cabo en una celda de calibracin. La historia de carga
del ensayo no 222 comienza con 150 KPa y llega hasta 734 KPa est descrita por 93
escalones de carga y la deformacin tangencial unitaria mxima de la membrana del
presimetro es del 10.19 %. En la Figura 5.1 se puede apreciar dicha historia de carga en
una representacin de eje vertical. presin radial sobre la pared de la cavidad
presiomhica, y de eje horizontal, deformacin tangencial en la pared de la cavidad. Las
condiciones de contorno en la celda de calibracin durante el ensayo se mantuvieron
constantes. es decir, oFcte y oh=cte.
Las principales carac'-risticas del ensayo no 222 quedan reflejadas en la Tabla 5.1. En
ella se observa que la muestra de arena tiene una densidad relativa (DR) del 46.2 (%)
despus de la fase de consolidaci6n. Los valotes de las tensiones eficaces vertical y
horizontal tras la fase de consolidacin evaluadas a la profundidad conespsndiente al
centro de la membrana del presimetro vienen dadas respectivamente por cr',=I 14 KPa
y o',,,=98 KPa. La presin de separacin de la membrana del presirt~etrodel soporte
rgido viene dada por p0=144 KPa. En condiciones ideales p, debera coincidir con a',,;
el hecho que esta coincidencia no tenga lugar puede deberse a varios motivos, entre
ellos se cita [Manassero, 198711a falta de precisin de los transductores que miden la
deformacin de la cavidad. El ndice de sobreconsolidacin tiene un valor de
OCR=5.497 y K,=0.85. cociente de a', y a',. Por ltimo el peso especifico de la
niuestra de arena en seco viene dado por y,=I .S19 t/m3.
. -
o -- - - - . ~-.
O 2 4 6 8 10 12
Deformacn de la Cavidad 1%)
(V=)ult
Kacv 2827 - -
Tabla 5.2. intervalos de variacin y valores iniciales de los parmetros del modelo de
Ferreira para el ensayo n0222.
Se ha trabajado tanto con parmetros del modelo sin escalar wmo escalados. Las curvas
de carga correspondientes a los valores ptimos de los parmetros en ambos casos se
representan en la Figura 5.2. En ella se aprecia el mal ajuste obtenido con los
parmetros sin escalar (Pres. Cav. Marq.) y un ajuste muy bueno al emplear parmetros
escalados (Pres. Cav. Marq. esc.)
.-.- - . - -- - - -~ -
Los parmetros ptimos correspondii 'ites a los dos aniilisis inversos quedan reflejados
en la Tabla 5.3. En ella puede observarse la diferencia en el valor de la funcin objetivo
(Fun. Obj.), que es mucho menor en el caso de escalar los parmetros. Este hecho se ha
conseguido con un coste mayor (ITN y no Evals.). En lo que respecta a los parmetros
del modelo la diferencia ms notable aparece en el parhetros (S,).
+Pres.Cav.(Exp)
+Pres.Cav. Complex
o
O 2 4 6 8 10 12
Deformacin de la Cavidad (%I
En la Tabla 5.4 se indican los parmetros del modelo de Ferreira que han conseguido el
ajuste ptimo de la curva de carga del ensayo no 222. En ella puede observarse que el
ninimo valor de la funcin objetivo es similar al obtenido con el algoritmo de
Marquardt escalando los parmetros (Tabla 5.3) aunque el coste temporal es mucho
menor con el algoritmo Complex. En lo que respecta a los parmetros ptimos son muy
similares a los obtenidos escalando los parmetros en el caso de Marquardt.
- - -
PI
13
-a
Sa-u
s::
~j
:o
6
a
o
O 2 4 6 8 10 12
Deformacibn da Is Cavidad 1%)
En la Tabla 5.6 se detallan los parmetros ptimos del modelo de Ferreira obtenidos en
cada uno de los anlisis y las caractersticas fiindamentales de cada uno de los anlisis
desarrolladospara el ensayo no222.
Tabla 5.6. Resrtlfadosde la oplimizu~rncon lo.?algoritmos de Recocido Simrrludo,
Cbmplex y Recorido Simulado. Ensayo no222.
El analisis de incertidumbre que se presenta para los tres ensayos sigue los pasos del que
se ha desarrollado en el capitulo anterior para el conjunto de ensayos ficticios. A
continuacion se presentan los resultados del anlisis de incenidumbre del ensayo no 722.
En la Tabla 5.7 se refleja la incenidumbre de los parmetros del modelo de Femira en
cada uno de los anlisis llevados a cabo.
Deformacin da la Cavidad %l
(VS)uit I O 5
Si 400 10 700
KaC' .2827 - -
A
Las curvas de carga correspondientesa los valores ptimos de los parmetros al realizar
el anlisis inverso con parmetros sin escalar y escalados se encuentran representados en
la Figura 5.8. En ella puede apreciarse que s u d e algo parecido a lo que suceda en la
Figura 5.2. Tambin aqui el ajuste obtenido sin escalar los parmetros es bastante malo
CApfTULO5: APLICACI~NA EJEMPLOS REALES
(Pres. Cav. Marq.). Al escalarlos, la curva de carga resultante (Pres. Cav. Marq. esc)
ajusta muy bien la curva de carga experimental (Pres. Cav. Exp).
O - - -
O 2 4 6 8 10 12
Deformacin de la Cavidad 1%)
Los parmetros ptimos que wrresponden con las curvas (Pres. Cav. Marq.) y (Pres.
Cav. Marq. esc) se detallan en la Tabla 5.10. En elle puede apreciarse como se consigue
un valor bastante mejor de la funci6n objetivo escalando los parmetros. Igual que en el
caso del ensayo no 222 el parmetro del modelo que presenta mayor diferencia en los
dos anlisis es Si.
En la Figura 5.9 se observa lo mismo que sucedia con el residuo en valor absoluto en el
ensayo no 222 (Figura 5.3). es decir ste resulta menor y ms homogneo cuando hay un
escalamiento de los parmetros (Residuo-Esc-Abs) que cuando no lo hay (Residuo-
Abs).
+RESIDUO ABS
+R E S I D U M S C
1-0s resultados obtenidos al resolver el anlisis inverso del ensayo no 225 medimte el
algoritmo Complex se presentan en la Figura 5.10. El ajuste obtenido es bastante bueno,
sobre todo al comienzo de la fase de carga. Sin embargo, sta es la zona donde los
resultados experimentalesestn ms influenciados por la insercin del presidmetro. Por
ello en ocasiones lo que se puede hacer es considerar nicamente los datos
experinientales a partir de un porcentaje de deformacin -por ejemplo 3%-. o bien
ponderarlos crecientemente de modo que la influencia de la instalacin del presimetro
en los parmetros finales sea mnima. Los resultados que se presentan se han obtenido
sin escalar los parmetros del modelo.
- -- . - --
Ensayo preslomtrico no 226 ICornpiex)
lo00
900 -
800
700
600
603
+Pres. Cav. I E q I
400
m +Pres. Cav. Complex
200
100
En la Tabla 5.11 se indican los parmetros del modelo de Ferreira que han conseguido el
ajuste ptimo de la curva de carga del ensayo no 225. El mnimo valor de la funcin
objetivo es prcticamente igual al obtenido con el algoritmo de Marquardt escalando los
parmetros (Tabla 5.9) aunque el coste es mucho menor con el algoritmo Complex. En
lo que respecta a los parinetros ptimos son tambien muy semejantes a los obtenidos
escalando los partimetros en el caso de Marquardt.
-
(5)
Normal 63.2 .S735 34.7 .2827 35 1119.4 1.19 27 255 13610.8 10
modo que se lleva a cabo un anfilisis inverso con tres parmetros sin restricciones.
Tambien en este caso el ajuste es mejor al comienzo de la fase de carga por lo que los
comentarios hechos en el apartado 5.3.2 encaminados a evitar la influencia del proceso
de insercin del presimetro son tambin aqul vlidos
-
+Pres. Cav. (Erq>l
Pres. Cav. Rec. Sim
o . -
O 2 4 6 8 10 12
Deformacin de la Cavidad 1%)
198
A A M P L O S REALES
5: APLICACI~N
CAP/TU~O
Los resultados obtenidos en el anlisis inverso del ensayo no 225 con los algoritmos de
Marquardt, Complex y Recocido Simulado se presentan conjuntamente en la Figura
5.12. En ella pueden dpreciarse cada uno de los ajustes y como estos son mejores al
comienzo de la fase de carga. Los resultados peores son los que se obtienen con el
algoritmo de Marquardt sin escalar los parmetros.
o. - --- . -
O 2 4 6 8 10 12
Deformacin d e Is Cavidad (96)
199
CAP~TULO
5: APLICAC~~N
A UEMPLOS REALES
n9228
5.4 ENSAYO PRESIOMTRICO
La historia de carga del ensayo presiomtrico no 228 se representa er! 'a Figura 5.13.
Dicho proceso de carga comienza con 208 KPa y llega hasta 1660 KPa, est descrita por
111 escalones de carga y la deformacin tangencia1 unitaria mxima que tiene lugar es
del 10.25 %.
CAP/TU&O5: APLICACINA EJEMPLOS REALES
---- - - - - -
5.4.1 con el
Los pasos seguidos con el ensayo no 228 han sido los mismos que ya se han descrito
para los ensayos no 222 y no 225. Se ha comenzado realizando un analisis inverso a
panir de la curva de carga del ensayo no 228 mediante el algoritmo de Marquardt. En la
Tabla 5.16 be presentan los valores iniciales de los parmetros del modelo de Ferreira,
junto con los intervalos de variacin de cada uno de ellos.
(VS)ult 1 1 0 5
KaCY 2827
1 -
I'n la Figura 5.14 se representan las curvas de carga correspondientes a los valores
bpiiinos de los parnietros al realizar el anilis. inverso con parmetros sin escalar
(iJres. Cat . Xlarq.) y escalados (Pws. Cav. M q .esc). En ambos casos los ajustes que
se consiguen son aparenremente buenos no habiendo apenas diferencia entre uno y otro.
o
O 2 4 6 8 1C 12
Defotmacidn de la Cavidad (%I
Las diferencias entre los dos an6lisis quedan reflejadas es la Tabla 5.17. En ella puede
apreciarse que en esta ocasibn no se consiguen mejores resultados escalando los
parmetros. Pues aunque los parmetros finales del modelo son semejantes, el valor
minimo de la funcibn objetivo es menor al no escalar los parmetros. Adems result~
tambin ms barato en tiempo.
-.
Los parmetros dptimos del modelo de Ferreira que han conseguido el mejor ajsste de la
curva de carga del ensayo no 228 se muestran en la Tabla 5.18. Como puede apreciarse
se obtiene un valor de la funcin objetivo inferior al conseguido con el algoritmo de
Marquardt (Tabla 5.17) a menor coste. En cuanto a los parmetros ptimos del modelo
como se puede apreciar en las Tablas 5.17 y 5.18 son muy semejantes.
. - -
o -.
O 2 4 6 8 10 12
Deformacin de la Cavidad 1%)
Ensayo
"'228
+' Ci u o*
ITN Evab. Fun. t
) S ) ~S
( ~ p ~( U i~ c
~ v o ( ~ p ~ (O)
)
Obj. (S)
Normal 174.4 643 99.7 -2827 40 59245 7.5 30 430 11939.9 13
5.3.3 -0
. .
dc Recocido Simulado
o . -
O 2 4 6 8 10 12
Deformacin de la Cavidad I%)
Ensayo a'bo + Cg U n*
Evais. Fua.
t
lter T
n.228 ( ~ p (Us),,lt
~ ) Si hcv (o) ((KP~) (O) (S)
Obj.
Nonnal 152.06 .63 869.2 .2827 39.05 73307. 627 1705 10921 87 30 1000
Los resultados obtenidos en el anhlisis inverso del ensayo no 228 con los algoritmos de
Marquardt. Complex y Recocido Simulado se presentan conjuniamente en la Figura
5.18. En ella pueden apreciarse cada uno de los njustes obtenidos y que todos ellos son
muy semejantes.
-t P ~ C~V.
S EWI
+Pms Cav. Maq (ex)
4Pms Cav Maq
-ct Pms Cev. Com>iex
. -
O 2 4 6 8 10
Deformacin de la Cavidad t%l
Las diferencias entre unos y otros pueden observarse con ma;or detenimiento en la
Tabla 5.20. En ella se han incluido los resultados para algunos parjmetros obtenidos por
otros autores que estn disponibles para este ensayo [Manassero, 1987; Ferreim, 19921.
Se puede apreciar que todos se corresponden bastante bien. El parmetro del modelo de
Ferrera que permanece ms estable en todos los anlisis es (I/s)~,.
CAPITULO5: APLICACI~NA EJEMPLOS REALES
Los resultados del anAlisis de incertidumbre de los parmetros ptimos del modelo de
Ferreira para el ensayo no228 quedan reflejados en la Tabla 5.21.
En ella se indica la fiabilidad de cada uno de los parmetros de la misma forma que se
ha citado para los dos ensayos anteriores. Puede apreciarse que de los anlisis llevados a
cabo con el mtodo de Marquardt presentan menores valores de o los parbetros
obtenidos sin escalar los parmetros. En cuanto a los resultados del anilisis de
incertidumbre para los paidmetros obtenidos a partir de los dos metodos directos
estudiados se observa la misma tendencia que en los ensayos no 222 y no 225, es decir.
los valores de 0 son menores para los parmetros obtenidos mediante el mtodo de
Recocido Simulado.
En la Tabla 5.22 se presentan los valores iniciales de los parmetros del modelo de
Drucker-lJrager. junto con los iniervalos de variacin de cada uno de ellos. Puede
observarse que en este caso no se ha fijado ningn parmetro, sino que se ha querido ver
como se comportaba la metodologia de optimizacin con todos los parmetros libres.
Tuhla 5.22. btten*ulosde variaciny valores iniciales de los parmetros del modelo de
Dritcker-Pruger en la modelizacin del ensayo no228.
En la Tabla 5.23 se muestran los parmetros 6ptimos obtenidos con el metodo Complex.
Llama la atencin el tiempo empleado, que muestm perfectamente el coste que implica
el no considerar una ley analtica de comportamiento. El hecho de que el valor final de
la funcin objetivo sea mucho menor que el correspondiente de la Tabla 5.1 8 se debe a
que -como ya se coment al final del capitulo anterior- con el modelo de Drucker-Prager
se ha calculado la funcin objetivo como la diferencia entre los valores experimentales y
calculados de las deformaciones y no, como ocurre con el modelo de Ferreira, con los
valores de las presiones sobre la pared de la cavidad.
5.6 CONCLUSIONES
Los resultados obienidos con cada uno de los mtodos de optimizacin (Marquardt,
Complex y Recocido Simulado) para cada uno de los ensayos han sido coherentes entre
s. En las Tablas 5.6, 5.13 y 5.20 pueJen verse respectivamente -de un modo sinttico-
los resultados de los ensayos no 222. no 225 y no 228. En general. puede decirse que los
resultados obtenidos con los mtodos directos (Complex y Recocido Simulado) sin
escalar los parmetros del modelo de comportamiento de Ferreira son mejores que los
obtenidos con el mtodo de Marquardt. Con ste nicamente se obtienen buenos
resultados escalando los parmetros del modelo. aunque el coste computacional es
siempre mayor que el empleado por los mtodos directos.
Comparando entre si los resultados obtenidos con los metodos Complex y Recocido
Simulado se puede afirmar que se obtienen menores valores de la Funcin objetivo con
el nitodo de Kecocido Simulado aunque con un coste mayor que el empleado por el
mtodo Complex. Tambin son mejores los resultados del anlisis de incertidumbre de
los parmetros ptimos obtenidos con el metodo de Recocido Simulado.
Por ultimo decir que se ha comprobado que la metodologia desarrollada funciona bien
con el modelo no analtico de Dmcker-Pmger. aunque nicamente se ha verificado con
el mtodo directo de optimizacin Complex. El coste computacional resulta muy
elevado. lo que ha motivado que slo se aplicase a la optimizacin de parmetros del
modelo de comportamientoen el ensayo no 228 y no se pudiese realizar un anlisis de
incertidumbre sobre los mismos.
6: CONCLUSIONESY L~JEAsFUTURAS DE INVESTIGAC16N
CAP~TULO
Bhsicamente se puede decir u-- la metodologia desarrollada tiene dos etapas: en primer
lugar se han planteado en el Captulo 11 dos modelos matemiticos de comportamiento
del material arenoso (Ferreira y Drucker-Prager) que tienen en cuenta los fenmenos
fisicos que se producen durante la realizacin del experimento. y en segundo lugar, en el
Captulo 111, se han analizado distintos mtodos para la identificacin de los parmetros
caractersticosdel material a partir de los datos experimentales disponibles.
Por tanto. se ha podido comprobar que tos mtodos directos estudiados constituyen una
alternativa sencilla a los metodos clsicos basados en el cAlculo del jacobano y pueden
resultar ms eficaces y prcticos que stos.
y Torczon, 19911 para los cuales se puede demostrar la convergencia global [Torczon.
19971 y, algoritmos que permitan abordar problemas geomtricamente no lineales con
grandes deformaciones [Mahnken y Stein. 19971.
El principal inconveniente que presenta este ensayo es la alteracin que sufre el medio
que circunda a la perforacin. En medios arcillosos rgidos o en rocas blandas la
TIPOS DE PRESI~METROS
Y ALGUNOS
INTERPRETACI~N
A1.2.3 M e t r o de- e
En 10s ltimos atios y en relacin sobre todo con los trabajos mar adentro (off-shore) se
Iia desarrollado el presimetro de desplazamiento total -FDP-[Hughes y Robertson.
1985: Withers et al. 19861, tambin llamado prcsiopenetrmetro [Jzquel et al., 19821
presimetro de cono -CP- [Campanella et al. 1990; Houlsby y Schnaid, 19941.
Este procedimiento produce una alteracion grande del volumen de suelo sobre el que se
va a realizar el ensayo presiomtrico -zona tensionalmente cargada- por lo que la
expansin de la membrana deber ser suficienteniente grande para recoger informacin
de la zona no alterada. La interpretacin de estos ensayos se realiza normalmente sobre
la curva de la fase de descarga o de contraccin de la membrana. Adems las medidas
de deformacin obtenidas con el presimetro se pueden relacionar con las medidas de
resistencia del penetrmetro para una misma profundidad. Se trata de un ensayo que
fcilmente puede repetirse y que no depende del personal que en cada caso realice el
HISTRICA DE ALGUNOS MTODOS DE
ANEXO P TIPOS DE PRESIOMETROS Y R E V I S I ~ N
INTERPRETACI~N
ensayo. En la Figura Al-6 [Femeira, 1992) puede verse una representacin ideal dr este
tipo de ensdyo.
v
ANEXO I: TIPOS DE PRESI~METROSY REVISIN HIST~RICADE ALGUNOS M&TODOS DE
INTERP~ACICIN
LL
defonnacibn de la cavidad (S)
Todos los metodos de interpretacin que han surgido en las ltimas dcadas. algunos de
los cuales se van a presentar brevemente a continuacin, tienen sus ventajas e
inconvenientes. parten de ciertas hiptesis, presentan sus limitaciones y tienen campos
de aplicacin concretos.
Dependiendo del tipo de material qus se estudie los mtodos de interpretacin debern
ser capaces de suministrar informacin sobre diferentes parrimetros. As, si el medio
donde se realiza el ensayo cs arcilloso ;e prestar atencin especial al mdulo cortante. a
la tensibn horizontal in-situ, a I z resistencia cortante no drenada. al coeficiente
horizontal de consolidacin, a las tensiones efwtivas y a la presin de poro. Si el medio
es arenoso interesar el mduio cortante. la tensin horizontal inicial in-situ, la presin v
refieren a la interpretacin tanto de ensayos sobre medios arenosos como sobre niedios
arcillosos, prestando especial atencin a los ensayos drenados por ser el tipo de ensayos
sobre los que se valida el mktodo de interpretacin propuesto es esta tesis. En principio
estos mdtodos se ordenan segn sea la ley de comportamiento utilizada una ley tensin-
deformacin clsica o bien, una ley que se obtiene a partir de los datos experimeniales.
Con todo ello el estado de tensiones y deformaciones en los alrededores del presimetro
puede ser determinado con tan slo dos parmetros: el r.idulo de Young y el
coeficiente de Poisson. Los resultados de este mtodo se han mostrado nicamente
aceptables para los primeros estadios del experimento presiomtrico. A medida que la
deformacin de la cavidad aumenta -hasta llegar a ser mayor que el 10%- los resultados
son menos fiables.
Por otro lado Gibson y Anderson desarrollaron a partir de las hiptesis anteriores una
ecuacin analtica para ensayos arenosos. En ella se usaha el criterio de Mohr-Coulornb
para sealar el comienzo de la fase plstica. La ecuacin analitica permite llevar a cabo
una representacin bilogartmica de la presin sobre la cavidad bersus la deformacin de
la misma resultando una lnea recta cuya pendienre permite determinar el ngulo de
friccin. Los valores obtenidos del gngulo de friccin suelen resultar aceptables. Sin
embargo al no considerar el cambio de volumen del suelo alrededor del presimetro
dwante el transcurso del ensayo drenado -dilatancia-, ha sido un mtodo poco popular
en ensayos arenosos.
M-li
Es un mtodo complejo y de difcil aplicacin. Los puntos de partida de esta
metodologa son bsicamente los de Gibson y Anderson (1961) con la diferencia de que
ste si tiene en cuenta los cambios de volumen del suelo alrededor del presimetro
durante el proceso cortante. El principal inconveniente reside en que precisa de
informacin sobre el promedio de deformacin volumtrica en la rotura para diferentes
presiones de confinamiento, siendo necesarios por tanto ensayos de laboratorio
especiales que deben ser realizados con muestras de densidad relativa similar a la del
suelo ensayado in-situ. aplicando historias de carga semejantes a las del ensayo
presiomtrico y dentro de un rango de tensiones de confinamiento que depender de la
profundidad a la que se realicen los ensayos. El bgulo de rozamiento interno se obtiene
a partir de los datos del ensayo mediante un proceso iterativo. El ngulo que se obtiene
al final del proceso es hincin de dos factores. El primero se denomina factor de
expansin de la cavidad y est relacionado con la tensin normal principal en el suelo
alrededor del presimetro y con el ngulo de iccin. Este factor comienza calculndose
a partir de un ngulo de friccin inicial y de una presin final o ltima definida a partir
de los datos del ensayo. El segundo factor se denomina indice de rigidez reducida y se
calcula a partir del promedio de deformacin volumtnca a rotura previsto. del ngulo
de ficcin y del mdulo cortante que se determina a partir del primer tramo de la curva
presiomtrica. Correlacionando ambos factores se obtienen los sucesivos valores del
ngulo de rozamiento interno. E1 proceso iterativo fmlVa cqmndola diferencia entre el
ltimo valor del ngulo de friccin y el nnterior es suficientemente pequea.
do de C-
Jefferies desarroll la metodologia de Gibson y Anderson y la extendio a la fase de
descarga de ensayos SBP en arcillas, considerando el ensayo ideal. es decir. que no se
produce ningn tipo de alteracin en el medio. Las hiptesis de partida son las misnias
de Gibson y inderson. Jefferies dedujo cuatro ecuaciones analticas para simular un
ensayo presiomtrico. Dos de ellas para la fase de carga -fase elstica y fase plstica- y
las otras dos para la fase de descarga -f*.se elstica y fase plstica-. Permite determinar
tres parmetros del suelo: presin horizontal inicial, resistencia cortante no drenada y
mdulo cortante elstico. Para ello emple un procedimiento matematico de ajuste
\ isual en la pantalla de un ordenador -modelizacin asistida por ordenador- de la curva
experimental con las curvas de las ecuaciones analticas para diferentes juegos de
parametros. El mejor juego de pdmetros ser el que produzca un mejor ajuste. Esta fue
una de las primeras metodologias que incluy la fase de dcscarga de un ensayo
prcsiomcitrico. Al suponer que no existe alteracin del medio esta metodologia es
aplicable slo a ensayos SRP en arcillas. El inconveniente principal que tiene este
mtodo es la excesiva extensin de la respuesta elstico lineal que no resulta realista con
el coniponamiento de las arcillas naturales.
e d - UWiW ( 1988)
Esta nietodologia asume todas las hiptesis de pariida de Gibson y Anderson y fue
desarrollada para los ensayos presiomtricos de desplazamiento total ("Full-
Displacrrnent Pressuremeter" (FDP)). Trata las deformaciones elsticas como
deformaciones pequefias -deformacin de Cauchy-. las deformaciones plsticas como
grandes deformaciones -deformacin de Hencky- y considera nicamente la fase de
descarga del ensayo presiomtrico ya que en este tipo de ensayos la fase de carga no
representa una respuesta verdadera de un suelo no alterado. Permite determinar tres
pararnetros: resistencia cortante no drenada. mdulo elstico cortante y presin
horizontal inicial. Este ltimo parmetro es el que se considera menos fiable dentro de
esta metodologa al resultar normalmente su valor excesivamente alto. El modulo
elstico inicial se obtiene del primer tramo de la curva de descarga -fase elstica,
pequeiias deformaciones- y la resistencia cortante no drenada se determina mediante un
procedimiento grfico a partir de los datos de la fase de descarga -fase plstica. grandes
deforniaciones-. Esta metodologia tambin presenta el inconveniente de extender
demasiado el comportamiento elstico lineal del material. Adems el mdulo elstico
cortante que determina no suele ser vlido en aplicaciones prcticas.
Mtodo de H o u l s b m
Este metodo parte de las mismas hiptesis que Hughes (1977) slo que emplea un
modelo elastico perfectamente plstico y adems considera la curva de descarga del l
ensayo para la determinacin de parmetros del suelo. En los ensayos SBP la curva de
descarga puede considerarse que tiene la misma importancia que la correspondiente a la
fase de carga. La curva completa resultante de un ensayo SBP consta de cuatro tramos:
elhtico de carga. plstico de carga, elstico de descarga y plstiw de descarga. Los
pbntos de comienzo de la deformacin plstica en cada caso se determinan con el
criterio de Mohr-Coulomb de manera semejante a wmo lo hacian Gibson y Anderson
(1961) en medios arenosos. El mdulo elstico cortante se obtiene direciamente para 1
XII
ANEXO 1: TIWS DE PRESI~METROSY REVISIN HISTRICA DE ALGUNOS M ~ O D O S
DE
INTERPRETACIN
cada tramo elbtico a partir de la curva experimental. Para determinar los panhetros
plsticos se aplica una tcnica numrica de ajuste. De este modo se obtienen dos
mdulos cortantes y dos ngulos de friccin. uno para la fase de carga y otro para la fase
de descarga. La ley de la dilatancia de Rowe (1962) se puede emplear para determinar
los ngulos de dilatancia, y la presin horizontal inicial se obtiene a partir de la presion
de separacin de la membrana. En los ensayos drenados es dificil reproducir
analticamente la curva de descarga debido a la presin de poro estatica que empuja la
membrana. Del mismo modo que en el mtodo de Hughes este metodo contempla una
variacin de volumen lineal con la deformacin cortante. Los valores de los parametros
que se obtienen son aceptables nicamente cuando se considera que el hgulo de
friccin en la descarga es igual al ngulo de friccin a volumen constante. Esta
suposicin indica. segun la ley de Rowe, que durante la fase de descarga no se produce
variacin de volumen, lo cual es bastante poco frecuente y restrictivo [Jewell et al.
19801.
Xlll
TIPOS DE PRESI~METROS ALGUNOS
-
INTERPRETACI~N
informacin de la curva de carga. los valores de los parmetros que se obtienen con esta
nietodologia corresponden al suelo alterado ms que al medio natural.
Este mtodo trata el ensayo presiomtrico a partir del problema dc expansin de una
cavidad cilndrica. bajo la hiptesis de deformacin plana. Considera tanto la fase de
carga del ensayo presiomtrico SBP como la de descarga. Emplea una rehcin tensin-
deformacin hiperblica mediante el modelo de Duncan ( 1970) para la modelizacin de
la Ese elstica. y el modelo elastoplhtico CJS desarrollado por Cambou. Jarafi y
Sidoro (1989) para la fase plstica. Este ltimo modelo describe bien el
comportamiento de los materiales granulares sometidos a historias de carga complejas.
Mediante un proceso de identificacin manual logra determinar algunos de los
parmetros del modelo CJS. Por otro lado. dispone tambien de la posibilidad de
considerar el "creeping" que se presenta en cienos ensayos de larga duracin mediante
un modelo reholopico no lineal visco-elstico.
variacin de volumen considera una funcin potencial plstica que permite las
deformaciones volumtricas durante el proceso cortante. De las dos componentes que
tiene la deformacin total -elstica y plstica- considera que sta est formada
nicamente por ladeformacin plstica. Adems adopta dos superficies de fluencia -tipo
Mohr-Coulomb- y dos respuestas de la deformacin volumtrica segun se trate de
arenas densas (dilatacin) o de arenas sueltas (contraccin). A panir de una tension
horizontal inicial conocida y de la pendiente de la curva de carga para cierto grado de
deformacin. establece una ecuacin en la que se relaciona la tensin cortante con la
deformacin circunferencial. Esta ecuacin junto con las ecuaciones que gobiernan la
respuesta plstica en trminos de tension -tensin de fluencia- y de deformacin
volumtrica -potencial plstico- permiten obtener de modo incremental la curva del
esfuerzo cortante y la historia de tensiones efectivas en la pared dc la cavidad. De este
modo se pueden determinar el hgulo de friccin de pico. el mbdulo cortante y el hngulo
de dilatancia. La validacin del mtodo se hizo en laboritorio mediante un mini-
presimetro dentro de una celda tnaxial -ensayos triaxiales con cilindro hueco-. No se
han interpretado ensayos in-situ con este mtodo. Sin embargo se puede decir que al
emplear informacin directa de la curva de carga los valores de los parmetros
obtenidos estarn influenciados por la alteracin del medio a causa del proceso de
insercin d t la sonda.
Una alternativa a este mtodo. debida tambin a Ladanyi, considera que la deformacin
volumtnca total antes de la rotura es aquel!a que en la representacin grfica
log(AV/C' + ed versus log@/pd da una lnea recta. Segn la notacin de Ladanyi AV/V
es la deformacin volumtrica del presimetro, e, es la deformacin volumtnca total
del suelo antes de la r o a p es Ic prcsn en la pared de la cavidad y p , la tensin
horizontal inicial en el interior del medio.
Esta metodologa pretende interpretar los resultados de ensayos SBP sobre arenas
directamente a partir de 10s datos experimentales. Los principios de partida empleados
en este mtodo son bsicamente los empleados por Hughes (1977). La innovacin reside
en el modo de afrontar el problema y en las tcnicas matemhticas empleadas. Basandose
en el mtodo de diferencias finitas y dando el punto inicial de la curva de carga (hsn).
dondep, es la presin de separacin de la membrana y E, la deformacin circunferencia1
inicial -que se considera nula inicialmente cn =O- ,logra seguir paso a paso toda la curva
de la fase de carga obteniendo de este proceso la ley constitutiva. De esta manera
consigue captar la no linealidad de las tensiones y los cambios de volumen durante el
proceso cortante. Con este mtodo se pueden determinar el mdulo cortante, el hgulo
de friccin y la dilatancia mediante la aplicacin de la ley de Rowe. La ley de
comportamiento se obtiene de la curva de carga de un ensayo SBP, es decir,
directamente de los datos experimentales aunque en ocasiones, para evitar dispersin en
los datos. ajusta la curva experimental con una funcin polinmica.
De los cuatro prrsimetros que se han comentado es este anexo nicamente se hace
mencibn a lo largo del desarrollo de la tesis a los dos primeros. es decir al presimetro
con preperforacin y al presimetro autoperforante. Los otros dos, el presimetro de
empuje y el presimetro de desplazamiento total se han citado para que se conozcan otro
tipo de dispositivos que se aplican a medios distintos.
Por otro lado esta breve descripcin histrica de los mtodos de interpretacin de los
ensayos presiomtricos drenados tambin pone de manifiesto que muchos mtodos
precisan adems de la informacin proveniente de los datos experimentalesdel ensayo.
de alguna informacin adicional para determinar los parmetros que en cada caso
caracterizan el comportamiento de las arenas ensayadas. Esta informacin adicional se
obtiene de experimentos en laboratorio con muestras inalteradas, muy dificiles de
conseguir especialmenteen terrenos arenosos.
Los ensayos con los que se ha trabajado son ensayos SBP drenados realizados en celda
de calibracin por lo que se ha puesto atencin especial en los procedimientos de
DE ALGUNOS MTODOS DE
ANEXO 1: TIPOS DE PRESIMETROSY REVISI6N HHIST~RICA
INTERPRETACI~N
'
interpretacidn que tengan en cuenta los cambios de volumen del medio en los
alrededores de la cavidad presiom6trica.
Etapa 1.- En el proceso iterativo que se describe las cargas aplicadas en la iteracin i-
esinia son en realidad las fuenas residuales de la iteracin anterior ( i -1). Estas fuenas
junto con la matriz de rigidez ( K ) g el vector de desplazamientos (u) forman el sistema
de ecuaciones que se resuelve en cada iteracin. La solucin del sistema da lugar a un
incremento de desplazamientos du' que a su vez permite calcular en cada punto de
integracin el incremento de deformacin &'.
Etapa 4.- Este paso depende de si se produce o no fluencia en los puntos de Gauss
durante la iteracin i-l sima. Por tanto. es preciso comprobar si
.
E'" > o, = o>+ Ifs;-' donde o'-'y o , son respectivamente la tensin efectiva y la
tensin de fluencia uniaxial del incremento considerado. En nuestro caso. al considerar
el criterio de fluencia de Drucker-Prager se tiene que E = d , +(J;)"' y a , = K .
donde J , es el primer invariante de tensin. J* es el segundo invatiante de tensiones
desviadoras y a y K son parmetros del modelo funcin de la cohesin p del angulo
de friccin interna ( apartado 2.4.3; criterio de fluencia). Por otro lado K es el modulo
de endurecimiento lineal y E;-' es la deformacin plstica efectiva al tinalizar la
iteracin i-l sima.
4.1.2.- Si 5:>o'-' en el punto de Gauss, que previamente habia sufrido fluencia, las
tensiones siguen creciendo. Por tanto todo el exceso de tensin dado por o: -o'-' debe
ser reconducido hasta la supeficie de fluencia, tal y como se refleja en la Figura AIl. 1.
La porcin de tensin que se debe modificar para que se satisfaga el criterio de fluencia
viene dado por el factor de reduccin R que en este caso vale 1.
Etapa 6.- La porcin de tensiones que resta. es decir, Rda: se debe transformar de
alguna forma en deformacibn plhstica. Es decir el punto A (Figura A11.2) al permitir que
el material se deforme plsticamente debe pasar a situarse sobre la superficie de
fluencia. En sentido estricto se puede decir que la carga aplicada en la nueva iteracin
hace que el punto C se mueva elsticamente hasta que se encuentra con la superficie de
tluencia B. Si no existiera la superficie de fluencia el material continuara su
comportamiento elstico hasta llegar a A. Sin embargo, es estado de tensin A no puede
moverse por fuera de la superficie de fluencia, que s61o podrh atravesar mientras se
satisfacen tanto las ecuaciones de equilibrio como la relacin constitutiva. Por tanto a
partir de las ecuaciones 2.16.2.17 y 2.1 8 se tiene.
o bien.
que indica la tensin total a' que satisface las condiciones elastoplsticas despus de
incrementar la tensin desde a'-'. En la Figura A11.2 se observa como el punto de
tension correspondiente a o' (D) puede no estar situado sobre la superficie de fluencia.
Esta situacin se puede solventar haciendo que los incrementos de carga aplicada sean
suficientementepequeos.
Otra posibilidad que tambin contempla esta subrutina es la de hacer que el punto D se
site sobre la superficie de fluencia al escalar o ' . Es decir. como la tensin efectiva 5'
debida a la tensin o' debe coincidir con el valor o , = o :+ fZp cuando el punto D
se encuentra sobre la superficiede fluencia. se toma como factor de escala el siguiente.
ANEXO 11: D&TERMMACI~N
DE LA DEFORMACI~NELASTOPLASTICA
Esto representa un escalamiento del vector o' que implica una reduccin porporcionill
de las componentes individuales de la tensidn. La condicibn de normalidad citada en el
apartado 2.4.4 se muestra evidente en la Figura A11.2 ya que DdA a = Drlc F .
En la Figura A11.3 puede apreciarse cul sera la posicin del punto final de tensisn (D)
y el punto D' que se obtiene sobre la superficie de fluencia al escalar o'.En este caso
si se divide el exceso de tensin elsiica Rdci (AB) en un nmero ni de intervalos. y
despuds de cada incremento se reduce el punto correspondiente a la superficie de
fluencia se obtendra el punto E (Figura A11.3, con m=3). Dicho punto podra corregirse
mediante un escalamiento obteniendo el punto E', que como puede apreciarse dista
bastante del D'.
(y) Otl
donde 5:.-a , da una idea del exceso de tensin AB y a '; es la tensin uniaxial
inicial.
Etapa 8.- Por ltimo se calculan las fuenus nodales equivalentes mediante la expresin
MVI
( Dimensionamientode variables y matrices (
-
I
Ciclo sobre cada uno de los elementos del dominio
1
Geometra y propiedades del material del elemento en cuestin
1
Inicializacina cero de la matriz de
1
Ciclo de integracin numrica sobre los puntos de integracin
1
Posicin (54) del punto de integracin correspondiente
1
Evaluacin de la matriz de elasticidad D
1
Para cada punto de integracin en coordenadas (C.?) proporciona
las funciones de forma N, y las derivadas a N, la5 y d N, I d .q
I
11 ( A partir de N, .a N, /aC y a N, l a q proporcona las derivadas parciales de las 1
l
funciones de forma con respecto a las coordenadas cartesianas @.y) o cilindricas
(r.z).la matriz jacobiana J. su inversa J'. el valor del determinante de J y las
coordenadas cartesianas o cilndricas del punto de gauss correspondiente
I
1
Calcula el volumen elemental de integracin dQ como(AW,W, multiplicado por
el espesor comspondiente o por 21 r (problemas axisimtricos). W: y I$ son
las funciones de peso para cada punto de integracin.
1
El ensayo presiomtrico es una de las pruebas in-situ que mejor permiten obtener
informacin acerca de los partimetros de deformacin. de resistencia y del estado
tensional de los suelos sometidos a estudio. En medios granulares se suele emplear el
presimetro autoperforante (Self-Boring Pressuremeter, SBP) con el fin de minimizar la
alteracin del medio al realizar la perforacin necesaria para introducir la sonda. Los
datos con los que se ha trabajado en este trabajo de tesis provienen de una serie de
ensayos presiometricos drenados en celda de calibracin realizados en los laboratorios
del ENEL-CRIS en MilAn [Belloti et al, 19871. A continuacin se hace una descripcin
de los dispositivos empleados en la realizacin de los ensayos y de cmo se han llevado
a cabo.
Entre la muestra de arena y la celda hay una membrana flexible que rodea lateralmente
y por la base a la muestra. La membrana que la rodea lateralmente esta sellada a una
placa de aluminio que constituye el limite superior de la muestra de arena. Por encima
de dicha placa se encuentra una tapa muy rgida que impide que se produzcan
deformaciones en sentido vertical. Tanto en el centro de la placa como de la tzp., igida
existe una cavidad de 12 cm. de dimetro para posibilitar el acceso fl; la sonda que
corresponda en cada caso.
Las presiones sobre la muestra se pueden aplicar tanto en sentido vertical como
horizontal. Las presiones en sentido vertical se aplican mediante un pistn ngido que se
encuentra en la base de la celda y que esta rodeado de una membrana flexible. La
cavidad que queda entre el pistn (Figura AIII.1) y la membrana se llena de agua
dcsgasiicada lo que permite aplicar una presin uniforme sobre la cara inferior de la
muestra. El desplazamiento mximo del pistn es de unos 10 cm. Las presiones en
sentido horizontal se consiguen al introducir agua a presin en la cavidad que rodea
lateralmente la muestra. La presin se consigue mediante dos depsitos de agua-aire
regulando la presion de aire.
Con el fin de evitar en lo posible la alteracin del medio que ocasiona la insercin del
presimetru, ste se introduce en la celda de calibracin antes que la arena. A
cuntinuacin se forma la muestra de arena mediante la tcnica de deposicin pluvial. De
este modo se consigue una muestra muy uniforme. la posibilidad de que el ensayo se
pueda repetir en las mismas condiciones y la obtencin de muestras de diferentes
densidades relativas (DR). Por consiguiente se pueden estudiar muestras dentro de un
amplio rango de densidades.
Esta tcnica de deposicin pluvial (Figura A111.2) sigue los pasos de Jacobsen (1976) y
Battaglio et al (1979) y para ello emplea un distribuidor fijo de arena que se coloca
sobre la celda de calibracin. Dicho distribuidor est formado por:
- Un difusor. Una vez que la arena atraviesa las placas perforadas experimenta
una cada de un metro aproximadamente por el interior de una cavidad
cilndrica del mismo dimetro que el depbsito hasta llegar a la parte superior
de la celda de calibracin donde est situado un difusor de acero. Este
difusor tiene la funcin de hacer que la deposicion de la a t n a sea uniforme
y as obtener una muestra homognea.
Distribuidor
de
arena
11 Arena - n, arena
Placas
wrfomda
alre
comprimido
Cavidad
/ cilindrica
-- Difusor
W H
I Chapa
celda
de Muestra r r/ cilindrica
calibracin de -Membrana
-
arena lateral
Membrana
inferior
Adems de la medida de la deformacin radial es preciso conocer la presin del gas que
se est inyectando. Dicha medida se realza en el interior de la cavidad de expansin
mediante una celda de presin. De este modo la medida de la presin no se ve afectada
por las prdidas de presin que tienen lugar entre la bomba de presin y la cavidad de
expansin. lo que oeurriria en el caso de considerar la presin del gas como la presin
del gas de inyeccin en la bomba.
De lo visto se deduce que van a ser cinco los parmetros de salida que esta sonda
presiom4trica proporciona en un ensayo de arena en celda de calibracin. Dichos
parmetros corresponden a los tres desplazamientos radiales medidos en la seccin
media de la membrana y que distan 120, a la media de los desplazamientos y a la
presibn de gas de expansin.
Este presimetro tambin tiene la capacidad de medir la presin de poro del medio. En
el presente trabajo al tratarse de ensayos sobre arenas secas los transductores de presin
de poro no se emplean.
A continuacin se hace una descripcin breve de los pasos seguidos para la realizacin
de los ensayos presiomtricos de arenas en celda de calibracin, es decir, de la insercin
de la sonda. de la formacin de la muestra de arena. de las condiciones de rontorno
aplicadas sobre ella. de la prueba de expansin propiamente dicha y de la determinacin
por ltimo. de la densidad de la muestra (Figura A111.3).
presiomtrica
Formacin
muestra de arena
contorno en la
expansin
celda de calibracin
y determinacin
de la densidad relativa
Los ensayos realizados han sido de dos tipos con el fin de estudiar si la posicion del
presimetro influye en los resultados del experimento. Ensayos con e! presimetro
colocado en su posicin natural. es decir. el sistema autoperforante en la base, y
ensayos con el presimetro invertido. es decir. el sistema de perforacin en la zona
superior. En ambos casos. para favorecer la deposicin uniforme de la arena alrededor
de la sonda, se ha colocado sobre la misma un cilindro metlico de su mismo dimetro
que llega hasta el difusor. Este cilindro hueco tambien sirve para alo;ar el cable que
sostiene a la sonda en su posicin natural, el conducto por donde se introduce el gas de
expansin y el resto de los cables de medida. De los ensayos llevados a cabo se deduce
que la posicin (natural - invertida) no afecta prcticamente a los resultados del
experimento.
Una vez alcanzado el Nvel de carga deseado se puede imponer una de las cuatro
condiciones de wntomo siguientes:
En algunas pruebas se han realizado durante la fase principal de carga varios ciclos de
descarga-recarga, y durante !a fase principal de descarga algn ciclo de carga-descarga.
En un ensayo SBP se admite que la insercin de la ronda se realiza con una alteracin
minima del terreno que circunda a la sonda. Por este motivo se puede medir la presion
geostatica horizontal total (o*)que existe en el medio a la profundidad donde se sita el
ANEXO 111: DEscRIPCIN TkCNICA DE LOS ENSAYOS Y CARACTERISTICAS GENERALES
DE LAS ARENAS
En los ensayos en celda de calibracin lo dicho hasta este momento acerca de la presin
total horizontal y la presin de separacin de la membrana del cilindro metlico debera
verificarse w n mayor motivo debido a la instalacin "ideal" llevada a cabo. Sin
embargo se ha obsemdo que los valores de o,,,y po no siempre coinciden. sobre todo
en niuestras de arena densa. Las posibles causas de esta no coincidencia son tres:
De las tres la que se ha considerado mas decisiva ha sido la ultima que se refiere a la
sonda empleada. Las dos primeras tienen que ver con la formacin de la muestra que en
principio debiera ser uniforme al haberse realizado en condiciones "ideales".
Las medidas resultaron menos dispares al modificar los transductores. es decir los
valores de o*,,y p, se aproximaron. Las diferencias permanecan en los ensayos donde
la pendiente inicial de la curva (a,,-%) era mayor, que corresponde a materiales poco
deformables. es decir, a arenas densas.
Por lo t a t o con el procedimiento descrito se logra que la alrencibn del medio durante la
realizacin del ensayo sea minima. o lo que es lo mismo que el ensayo se realice en
condiciones ideales. Ademb. gracias a la deposicin pluvial se consigue una muestra
hornogenea e istropa. Un ensayo presiomtrico, bajo unas condiciones de contorno
determinadas y sometido a una historia de carga pariicular, proporcionar como
respuesta una curva representada en los ejes presin sobre la cavidad versus
deformacin de la pared de la cavidad.
sonda. En todos los ensayos que se analizan en esta tesis se considera nicamente el
tramo de curva correspondiente a la fase de carga del ensayo.
AIII.5 CARACTER~STICAS
GENERALES DE LAS ARENAS
Desde un punto de vista ingenieril el suelo est formado por todos aquellos materiales
de la superficie de la corteza terrestre susceptibles de ser desplazados por una azada.
Estos materiales suelen estar formados par granos, agua y aire, es decir por
componentes de las tres fases: slida. liquida y gaseosa. El contenido 2c dgua y aire eii
suelo vana en funcin de las condiciones ambientales por lo que cuando se habla de
suelo se suele hablar propiamente de su fase slida. La fase liquida y gaseosa sc
engloban en el indice de poros o en la porosidad.
Los suelos granulares que tienen un camaiio de grano comprendido entre 2 inin. y 0.05 b
0.03 mm dependiendode la clasificacin que se use, se llaman arenas. Sobre este tipo de
material. las arciias. es sobre el que se ha de.sarrollado el trabajo que aqu se describe. En
el comportami-nto dc las arenas influyen distintos factores: carncreristicas del suelo.
estado de tensiones, grado de compactacin. historia de carga eic.
Indica el estado tensional al que se encuentra sometido la arena que se est estudiando.
Cuando la atena est sometida a una tensin horizontal considerable existe menor
tendencia al cambio de volumen y al asiento o subsidencia del suelo. Si en un suelo se
retira una sobrecarga vertical permanece un cierto porcentaje de tensin horiwntal (K,
coefiicient of lateral ea& pressure at rest). El coeficiente K, se define wmo el ratio
tension horiwntatensibn venical en un estado de equilibrio (deformacin nula).
ANEYO 111: DESCRIPCI~N T~~CNICADE LOS ENSAYOS Y CARACTER~STICASGENERALES
DE LAS ARENAS
Tambidn se pueden incluir aqu factores como la situacin y variacin del nivel fretico
o estados pretensionales o predefonnativos.
El nivel de compactacin tiene que ver con la disposicin, forma y orientacin de los
granos y con la densidad del medio. Es muy importante conocer cul es el nivel de
compactacin inicial y su variabilidad en el medio que se quiere estudiar. Si el medio
esta compuesto por materiales sueltos de gran vwiabilidad se producirn asentamientos
grandes. Por el contrario en depsitos de arenas densas apenas habr asentaniientos y la
arena se expander ante esfuems cortantes.
- Mica W).
DE LOS ENSAYOS Y CARICTERISTICAS G M E ~ L E S
ANEXO 111: DESCRJPCI~NT ~ N I C A
DE LAS ARENAS 1
Se puede apreciar que se asignan valores menores del ngulo de friccin a volumen
constante a las arenas de grano ms redondeado y valores mayores paro las arenas con
granos ms angulosos. En el caso de las arenas del Ticino (TS-4). arenas con buena
gradacion. se ha considerado un ngulo de friccin a volumen constante de 34' o 35"
dependiendo de quien las haya estudiado.
tomo fruto de las estancias beves llevadas a cabo en el lsmes se ha estudiado cul es el
comportamiento de este tipo de arena a partir de un wnjunto de ensayos triaxiales y
presiometricos. Con ellos se ha realizado un trabajo de comparacin tratando siempre
yue los ensayos comparados difirieran tan slo en-una o factor con el fin-de
apreciar la influencia de dicha variable en el wmponamiento de la arena.
Del estudio de comparacin llevado a cabo a partir de los grificos mencionados han
obtenido las siguientes conclusion~srelativas al comportamiento de las arenas TC-4:
- 1Jna muestra con una historia dc carga de expansin (EPC). a igualdad del
resto de los factores. precisa de mayor tensin para alcanzar la rotura que
cuando se carga mediante compresion (CPC).
- tina muesm consolidada con Cf&,D, a igualdad del resto de los factores,
rumpe a mayor tensin que una consolidada mediante CID. y adems la
rotura se produce w n menor defomiacin axiul.
- A igualdad del resto de lus factores. una muestra con una historia de carga
CPC rompe frgilmente, mientras que la cargada con P-const. o S-const.
presenta endurecimiento.
f'ara una misma tension de confinamiento la evolucin de una arena densa y de una
arena suelta al aumentar el ratio desviador (al/03)aparece reflejada en la Figura Ali1.6.
ANGXO 111: DESCRIPCI~NTTNICA DE LOS ENSAYOS Y CARACTER~STICAS GENERALES
DE LAS ARENAS
XLV
A LOS ENSAYOS Y CARACTERISTICAS GENERALES
ANWO 111: DESCRIPCI~NT ~ ~ N I c DE
DE LAS ARENAS
Por otro lado tambin puede afirmarse que una arena de alta densidad relativa puede
tener un comportamiento con endurecimiento cuando las presiones de confinamiento
son muy altas y. de la misma forma una arena de baja densidad relativa puede presentar
una rotura frgil si las presiones de confinamiento son suficientementepequefias.
1.a densidad relativa y la presin de confinamiento son los factores que m& influyen en
el ngulo de friccin de pico. El ndice de sobreconsolidacin (OCR) influye
principalmente en la deformabilidad. En lo que se refiere a la rotura. a paridad de
densidad relativa, una arena normaimente consolidada (NC) y una sobreconsolidada
(OC)presentan prcticamente la misma wna de rotura.
Los datos que se emplean en este apartado para estudiar el comportamiento de la arena
TC-1 corresponden a la respuesta se una serie de ensayos presiomiricos realizados en
celda de calibracin. A la hora de comparar las respuestas de los diferentes ensayos se
ha puesto especial atencin en la presion interna de la cavidad presiomkirira y eti las
tensicines vertical (SlGvo) y horizontal (SIGho) aplicadas en la celda de calibracin. Los
casos estudiados corresponden con aquellos en los que las condiciones de contorno son
presin venical (o,)y presin horizontal (a,,)constantes. La representacin que se ha
realizado para la wmparacibn de ensayos ha sido la ms frecuente, es decir, aquella en
la que se representa la deformacin tangencia1 o circunferencia1 (q,)con la presin del
iiiterii~rde la cavidad (p).