IMPLEMENTACIN DEL ENSAYO DE COLUMNA RESONANTE

DERLY ANDREA MURCIA LPEZ

LEIDY MAYERLY PERDOMO CRUZ

UNIVERSIDAD MILITAR NUEVA GRANADA

FACULTAD DE INGENIERA
PROGRAMA DE INGENIERA CIVIL
BOGOTA, D.C., 2012
IMPLEMENTACIN DEL ENSAYO DE COLUMNA RESONANTE

DERLY ANDREA MURCIA LPEZ

LEIDY MAYERLY PERDOMO CRUZ

Trabajo de grado presentado como requisito

parcial para optar al ttulo de Ingeniero Civil

Director: Ing. JAVIER FERNANDO CAMACHO TAUTA Ph.D.

UNIVERSIDAD MILITAR NUEVA GRANADA

FACULTAD DE INGENIERA
PROGRAMA DE INGENIERA CIVIL
BOGOTA, D.C., 2012
AUTORIDADES DE LA UNIVERSIDAD MILITAR NUEVA GRANDA

MAYOR GENERAL (R) EDUARDO ANTONIO HERRERA BERBEL

Rector

BRIGADIER GENERAL (R) ALBERTO BRAVO SILVA

Vicerrector General

Dra. SONIA OSPINO GMEZ

Vicerrectora Administrativa

Dra. MARTHA LUCA BAHAMN JARA

Vicerrectora Acadmica

Dra. JACQUELINE BLANCO BLANCO

Vicerrectora de Investigaciones

Ing. ERNESTO VILLARREAL SILVA Ph.D.

Decano de la Facultad de Ingeniera

Ing. JOSE GONZALO ROS MARN

Director Programa de Ingeniera Civil
 NOTA DE ACEPTACION

El trabajo de grado titulado IMPLEMENTACION DEL ENSAYO DE

COLUMNA RESONANTE presentado por las estudiantes Derly Andrea
Murcia Lopez y Leidy Mayerly Perdomo Cruz, como requisito parcial para
optar por el titulo de Ingeniero Civil fue revisado y aprobado por:

_________________________________________
Ing. JAVIER FERNANDO CAMACHO TAUTA, PhD
Director

_________________________________________
Ing. _______________________________
Jurado

Bogot, septiembre de 2012

 A mis padres
A mi hermana
A mis abuelos, tos y primos

Andrea
 A mis padres
A mis hermanos
A mi abuela, tos y primos
En memoria de mi abuelo Jorge y mi to Roberto

Mayerly
 AGRADECIMIENTOS

En primera instancia agradecemos a Dios por la oportunidad que nos dio

inicialmente de terminar con buen resultado nuestra carrera y en segundo
trmino por darnos la fortaleza de llevar a cabo este trabajo de grado.

A la familia Murcia Lpez, en cabeza de Mauricio Murcia, mi padre y mi

madre Karen Lpez, que siempre han estado presente en todo instante y
ms aun cuando he sentido necesitar su ayuda.

A la Familia Perdomo Cruz, a mi padre Eduardo Perdomo y mi madre

Margarita Cruz que con su apoyo incondicional lograron dar las palabras
adecuadas para no desistir en el alcance de un nuevo logro en mi vida
personal y profesional.

A nuestros hermanos: Maura Murcia, Daniel Perdomo y Luis Perdomo que

han sido el motor fundamental en los momentos en que sentamos no poder
ms. A nuestras familias en general, abuelos, tos y primos.

Un agradecimiento especial al Ingeniero Javier Camacho, que siempre

estuvo presente para resolver las dudas que se generaron con el paso de los
das durante desarrollo del trabajo, gracias por su paciencia, tiempo dedicado
y conocimiento transmitido.

VIII
Al Ingeniero John Bermdez quien comparti su tiempo contribuyendo en el
desarrollo de este trabajo de grado, estando dispuesto a aprender y
ensearnos temas que tal vez no son afines a nuestra carrera.

A Eduardo Dulcey, encargado del laboratorio de suelos de la Universidad

que siempre estuvo dispuesto a brindar su colaboracin en los momentos en
que se haca uso de las instalaciones del laboratorio.

Al grupo de Investigacin de Geotecnia y a todo el grupo de estudiantes de

Ingeniera Civil que hacen parte el semillero de investigacin que de alguna
forma hicieron posible este trabajo.

A nuestros amigos y compaeros, a aquellos que iniciaron y tal vez no

terminaron, a los que se quedaron atrasados y a quienes con mucho
esfuerzo lograron dar fin a sus materias junto a nosotras, gracias por las
palabras de apoyo y al tiempo que juntos dedicamos a finalizar uno de
nuestros primeros grandes logros.

IX
 Advertencia

La Universidad Militar Nueva Granada no se

hace responsable de las opiniones y
conceptos expresados por los autores en sus
respectivos trabajos de grado; slo vela
porque no se publique nada contrario al
dogma ni a la moral catlica y porque el
trabajo no contenga ataques personales y
nicamente se vea en el anhelo de buscar la
verdad cientfica.

X
 TABLA DE CONTENIDO

AGRADECIMIENTOS .................................................................................. VIII

TABLA DE CONTENIDO ............................................................................... XI
LISTA DE FIGURAS .................................................................................... XIII
LISTA DE TABLAS .................................................................................... XVII
RESUMEN ................................................................................................ XVIII
ABSTRACT ................................................................................................. XIX
1 INTRODUCCION ..................................................................................... 1
1.1 JUSTIFICACION ................................................................................ 3
1.2 OBJETIVOS....................................................................................... 4
2 PROPIEDADES DINMICAS DEL SUELO.............................................. 7
2.1 Introduccin ....................................................................................... 7
2.2 Tipos de ondas .................................................................................. 7
2.2.1 Ondas superficiales .................................................................... 8
2.2.2 Ondas de cuerpo ......................................................................... 9
2.3 Parmetros dinmicos del suelo ...................................................... 10
2.3.1 Mdulo de corte......................................................................... 10
2.3.2 Relacin de amortiguamiento .................................................... 14
2.3.3 Variacin del comportamiento del suelo a lo largo de las
deformaciones ........................................................................................ 17
2.3.4 Esfuerzo efectivo principal y relacin de vacos ........................ 20
2.4 Medicin de las propiedades dinmicas del suelo ........................... 24
2.4.1 Ensayos de campo .................................................................... 25
2.4.2 Ensayos de laboratorio.............................................................. 29
3 ENSAYO DE COLUMNA RESONANTE ................................................ 34
3.1 Introduccin ..................................................................................... 34
3.2 Equipo de columna resonante ......................................................... 35

XI
 3.3 Ensayo convencional de columna resonante (RCCT) ..................... 40
3.4 Ensayo en el dominio de la frecuencia (RCFD) ............................... 44
3.5 Ensayo de no resonancia (NR) ....................................................... 47
4 ADQUISICION Y PROCESAMIENTO DE SEALES ............................ 49
4.1 Introduccin ..................................................................................... 49
4.2 Digitalizacin .................................................................................... 50
4.3 Seales de tiempo continuo............................................................. 53
4.4 Seales de tiempo discreto .............................................................. 53
4.5 Representacin de seales en trminos de sus componentes de
frecuencia .................................................................................................. 55
4.6 Reduccin de ruido .......................................................................... 57
4.7 Transformada de Fourier ................................................................. 58
4.8 Funcin de transferencia y coherencia ............................................ 59
4.9 Sistemas .......................................................................................... 59
5 IMPLEMENTACION DEL ENSAYO DE COLUMNA RESONANTE ....... 61
5.1 Introduccin ..................................................................................... 61
5.2 Descripcin programa en LabVIEW ................................................. 61
5.3 Calibracin ....................................................................................... 63
5.3.1 Calibracin del equipo de columna resonante .......................... 63
5.3.2 Calibracin de los sensores para la adquisicin de datos ......... 71
5.4 Ejecucin del ensayo de columna resonante................................... 80
5.4.1 Datos de calibracin .................................................................. 80
5.4.2 Datos iniciales del espcimen ................................................... 81
5.4.3 Datos de estado ........................................................................ 85
5.4.4 Otras propiedades ..................................................................... 87
5.4.5 Resultados obtenidos ................................................................ 89
6 CONCLUSIONES ................................................................................... 98
7 RECOMENDACIONES ........................................................................ 100
BIBLIOGRAFA ........................................................................................... 102

XII
 LISTA DE FIGURAS

Figura 1. Tipos de ondas (Kramer, 1996) ....................................................... 9

Figura 2. Diagrama tensin-deformacin en un material lineal
(Santibaez, 2006)..................................................................... 11
Figura 3. Diagrama tensin-deformacin en un material no lineal con
sus mdulos de corte (Santibaez, 2006) .................................. 11
Figura 4. Curva de reduccin de mdulo de corte y variacin de la
razn de amortiguacin con el nivel de deformacin de
corte para una arcilla blanda (Universidad Politcnica de
Catalunya, 2004)........................................................................ 13
Figura 5. Curvas de reduccin del mdulo de corte para arenas (Bala &
Balan, 2006)............................................................................... 13
Figura 6. Variacin de la relacin de amortiguamiento de suelos
granulares finos con amplitud cclica de corte-deformacin
y plasticidad ndice. Efecto de la plasticidad de suelos en
respuesta cclica (Kramer, 1996). .............................................. 16
Figura 7. Modelo del comportamiento del suelo en conformidad con las
caractersticas de deformacin (Kramer, 1996) ......................... 18
Figura 8. Ciclo de histresis debido a cargas armnicas en un cuerpo
viscoelstico (Camacho-Tauta, 2011)........................................ 19
Figura 9. Variacin del mdulo de corte con relacin de vacos y
presin de confinamiento para formas redondeadas y
angulares secas (lneas continuas) y arenas de grano
(lneas discontinuas) (Camacho-Tauta, 2011). .......................... 21
Figura 10. Ensayo de reflexin ssmica (Universidad Politcnica de
Catalunya, 2004)........................................................................ 25
Figura 11. Ensayo de refraccin ssmica (Rinaldi) ........................................ 26

XIII
Figura 12. Mtodo de cross hole (Cedex, 2010) ........................................... 27
Figura 13 Ensayo down hole (Cedex, 2010) ................................................. 28
Figura 14 Ensayo de Columna Resonante (Moreno & Francisco, 1998) ...... 30
Figura 15 Ensayo de Pulso Ultrasnico (Lombillo & Villegas, 2010)............. 31
Figura 16 Ensayo de Elementos Bender (Sandox, 2010) ............................. 32
Figura 17 Ensayo Triaxial Cclico (Concha, Poblete, & Mauro, 2004) .......... 33
Figura 18. Sistema del equipo de columna resonante .................................. 36
Figura 19. Montaje fijo-libre de columna resonante (Camacho-Tauta,
2011) .......................................................................................... 37
Figura 20. Montaje extremo activo con platina en la parte superior
(Camacho-Tauta, 2011) ............................................................. 38
Figura 21. Modificacin del sistema fijo-libre con platina en la parte
superior de la muestra (Camacho-Tauta, 2011) ........................ 39
Figura 22. Sistema cuasi libre de columna resonante (Camacho-Tauta,
2011) .......................................................................................... 40
Figura 23. Ensayo convencional de columna resonante (Camacho-
Tauta, 2011); a) Torque aplicado; b) Rotacin generada; c)
Figura de Lissajous .................................................................... 41
Figura 24. Determinacin del mdulo de amortiguamiento (Camacho-
Tauta, 2011); a) Interrupcin de la seal de entrada; b)
Ciclo ........................................................................................... 43
Figura 25. Ensayo en el dominio de la frecuencia (Camacho - Tauta,
2008); a) Torque aplicado; b) Rotacin generada...................... 46
Figura 26. Funcin de transferencia (Camacho - Tauta, 2008); a)
Magnitud; b) Angulo de fase ...................................................... 46
Figura 27. Ensayo de no resonancia (Camacho-Tauta, 2011) a) Torque;
b) Rotacin................................................................................. 48
Figura 28.Propiedades dinmicas obtenidas con un ensayo de no
resonancia (Camacho-Tauta, 2011) : a) Mdulo de corte; b)
Relacin de amortiguamiento .................................................... 48

XIV
Figura 29. Pasos para el procesamiento de datos en una
Transformada de Fourier Rpida (FFT) por sus siglas en
ingles. (Avitabile, 2001) .............................................................. 49
Figura 30. Digitalizacin de una seal analgica, adaptado de
(Orfanidis, 1996), (Camacho-Tauta, 2011) ................................ 52
Figura 31. Representacin de una seal de tiempo discreto ........................ 54
Figura 32. Ejemplos de espectros e amplitud de seal (Kamen & Heck,
2008) .......................................................................................... 56
Figura 33. Control de ruido por apilado de seal en el tiempo
dominante. La relacion de seal para ruido aumenta con
relacion l numero de seales incrementadas. (Santamaria &
Fratta, 2005) .............................................................................. 57
Figura 34. Panel frontal (National Instruments Corporation,, 2012) .............. 62
Figura 35. Diagrama de bloques (National Instruments Corporation,,
2012) .......................................................................................... 63
Figura 36. Barra de Calibracin (J1) ............................................................. 64
Figura 37. Masa adicional (J2) ...................................................................... 64
Figura 38. Calibracin Torque-Corriente ....................................................... 68
Figura 39 Determinacin de la relacin de amortiguamiento del aparato.
a) Curva de decaimiento de vibraciones libres. b) Amplitud -
nmero de ciclos (Camacho-Tauta, 2011) ................................. 69
Figura 40. Determinacin de la relacin de amortiguamiento del aparato
mediante el mtodo half-power. Amplitud de
desplazamiento y corriente para frecuencias cercanas a la
resonancia. a) Frecuencia; b) frecuencia normalizada.
(Camacho-Tauta, 2011) ............................................................. 70
Figura 41. Funcin de aceleracin ................................................................ 71
Figura 42. Funcin de velocidad ................................................................... 72
Figura 43. Funcin de desplazamiento ......................................................... 72
Figura 44. Relacin desplazamiento rotacin ............................................... 73

XV
Figura 45. Calibracin Celda de Carga ......................................................... 75
Figura 46. Calibracin Presin de Cmara ................................................... 76
Figura 47. Calibracin Contrapresin ........................................................... 77
Figura 48. Calibracin Presin de Poros ...................................................... 78
Figura 49. Calibracin Deformacin Vertical ................................................. 79
Figura 50. Calibracin Cambio Volumtrico .................................................. 80
Figura 51. Factores de calibracin ................................................................ 81
Figura 52. Datos iniciales del espcimen. ..................................................... 84
Figura 53. Datos de estado. .......................................................................... 85
Figura 54. Panel de navegacin. Ensayo en el dominio de la frecuencia. .... 92
Figura 55. Panel de navegacin. Ensayo de No Resonancia. ...................... 94

XVI
 LISTA DE TABLAS

Tabla 1. Constantes en ecuaciones empricas propuestas en mdulos

de deformacin pequea (Camacho-Tauta, 2011)......................... 22
Tabla 2. Frecuencia de resonancia del sistema + Barra fn1 ......................... 65
Tabla 3. Frecuencia de resonancia del sistema + Barra + Masa (fn2) .......... 65
Tabla 4. Datos de calibracin para la obtencin de la inercia rotacional
de la masa activa ........................................................................... 66
Tabla 5. Frecuencia de resonancia del sistema con el cabezal libre ............ 67
Tabla 6. Calibracin Torque-Corriente .......................................................... 67
Tabla 7. Calibracin Celda de Carga ............................................................ 74
Tabla 8. Calibracin Presin de Cmara ...................................................... 75
Tabla 9. Calibracin Contrapresin ............................................................... 76
Tabla 10. Calibracin Presin de Poros ........................................................ 77
Tabla 11. Calibracin Deformacin Vertical .................................................. 78
Tabla 12. Calibracin Cambio Volumtrico ................................................... 79
Tabla 13. Parmetros de muestreo para el ensayo en el dominio de la
frecuencia. (Camacho-Tauta, 2011)............................................... 93
Tabla 14. Parmetros de muestreo para el ensayo de no resonancia.
(Camacho-Tauta, 2011) ................................................................. 97

XVII
 RESUMEN

El presente trabajo de grado desarrolla la teora correspondiente al ensayo

de columna resonante, el cual est regulado por la norma ASTM D-4015,
cuyo propsito es determinar el mdulo de corte, relacin de
amortiguamiento y deformacin de corte. Cabe resaltar que el equipo cuenta
con la posibilidad de ejecutar adicionalmente el ensayo de bender elements y
torsin cclica, aunque no sern desarrollados dentro del presente trabajo ya
que no se encuentran delimitados en el alcance de esta investigacin.

El propsito del trabajo fue lograr la implementacin del ensayo en las

instalaciones del Laboratorio de Suelos de la Universidad Militar Nueva
Granada; el deseo de llevar a cabo esta investigacin fue transmitido por el
Ingeniero Javier Fernando Camacho Tauta, tutor del trabajo de grado, luego
de haber desarrollado su tesis doctoral en la Universidad Tcnica de Lisboa
acerca del tema mencionado.

Palabras claves: Mdulo de corte, relacin de amortiguamiento, columna

resonante, deformaciones.

XVIII
 ABSTRACT

The present degree work develops the theory corresponding to the resonant
column test, which is regulated by the standard ASTM D-4015. The purpose
is to determine both the shear modulus, damping ratio and shear strain. It is
worth noting that the equipment has the possibility to execute further Bender
elements test and cyclic torsion test, however it is out of scope of this
investigation.

The purpose of this work was to achieve the test implementation in the
facilities of the Soil Laboratory of New Granada Military University; and the
desire to carry out this research was transmitter by engineer Javier Fernando
Camacho, tutor of this degree work, after having developed his doctoral
thesis at the Technical University of Lisbon on the subject mentioned.

Key words: shear modulus, damping ratio, resonant column, deformations.

XIX
1 INTRODUCCION

La dinmica de suelos, siendo una rama de la ingeniera civil, requiere de un

profundo anlisis para lograr comprender la incidencia que su estudio puede
generar en la actualidad; vibraciones a causa de maquinaria o el simple paso
de vehculos, sismos y la ejecucin de fundaciones logran generar
modificaciones en el estado actual de un terreno, stas son algunas de las
razones por lo cual se hace necesario implementar un ensayo cuyo propsito
sea el estudio de las propiedades dinmicas del suelo en el laboratorio de la
Universidad Militar Nueva Granada.

El estudio de la dinmica de suelos necesita del empleo de diferentes

herramientas de orden conceptual, numrico y experimental, siendo as la
contraprestacin al campo investigativo de gran utilidad ya que mediante
esto se logra entender el comportamiento del suelo ante solicitaciones
dinmicas. Este trabajo de grado desarrolla temticas que logran abarcar las
diferentes herramientas que permiten tener una nocin mas clara acerca del
tema mencionado anteriormente.

El ensayo de columna resonante permite medir las propiedades dinmicas

del suelo a pequeas deformaciones (mdulo cortante y relacin de
amortiguamiento), no sobra mencionar, que el equipo implementado en las
instalaciones del laboratorio de suelos de la Universidad permite ejecutar
otros ensayos los cuales tienen la misma finalidad del ensayo inicialmente
mencionado, estos son: bender elements y torsin cclica, aunque no tendrn
un desarrollo en este trabajo ya que no se encuentran contemplado dentro
del alcance planteado.

1
La tesis doctoral del Ingeniero Javier Fernando Camacho Tauta, tutor de este
trabajo de grado, Evaluation of the Small-StrainStiffness of Soilby Non-
Conventional Dynamic Testing Methods fue uno de los trabajos que motiv
el desarrollo del presente trabajo. Finalmente, este proyecto se encuentra
enmarcado dentro de un proyecto de investigacin desarrollado por el grupo
de Grupo de Investigacin en Geotecnia, adscrito al Programa de Ingeniera
Civil de la Universidad Militar Nueva Granada cuyo nombre es: Rigidez del
suelo a pequeas deformaciones. Evaluacin del efecto de la frecuencia

2
1.1 JUSTIFICACION

A travs de los aos los inconvenientes relacionados con la propagacin de

las ondas, sismos y vibraciones en general, a raz de la aplicacin de cargas
dinmicas han logrado dar como resultado hechos catastrficos en diferentes
lugares del mundo; un ejemplo de la problemtica mencionada son los
terremotos ocasionados en los ltimos aos en Japn, Chile, Mxico, entre
otros. Por sta y otras razones se ha visto la necesidad de estudiar lo
relacionado con la dinmica de suelos, para tener un amplio conocimiento y
de esta forma plantear ptimas alternativas de solucin.

No es fcil dar soluciones concretas a este tipo de problemas planteados, por

tal razn el enfoque trazado para este trabajo de grado se basa nicamente
en la implementacin del ensayo de columna resonante. Como ya se ha
mencionado anteriormente este proyecto est respaldado por el proyecto de
investigacin: ING-723 Rigidez del Suelo a Pequeas Deformaciones.
Evaluacin del Efecto de la Frecuencia., ejecutado por el Grupo de
Investigacin de Geotecnia de la Universidad.

Una vez logrado el objetivo de la implementacin del equipo se estar en la

capacidad de ejecutar mltiples tipos de ensayo para fines acadmicos,
investigativos o de venta de servicios si as se desea.

3
1.2 OBJETIVOS

Implementar el ensayo de Columna Resonante para determinar las

propiedades dinmicas del suelo.

Identificar los conceptos ms importantes en dinmica de suelos y que

intervienen en la ejecucin, procesamiento y anlisis de los resultados
producidos por el ensayo de columna resonante.

Contribuir al diseo del programa que contrala el ensayo de columna

resonante haciendo uso del software LabView.

Desarrollar un manual para el usuario que permita manejar de forma

adecuada el equipo y por ende dar confianza en la ejecucin de los ensayos
llevados a cabo por la comunidad estudiantil de la Universidad Militar Nueva
Granada.

4
1.3 ANTECEDENTES

La dinmica de suelos es una ciencia que da a da se ha ido enriqueciendo y

ha suministrado distintas herramientas en el campo cientfico e investigativo.
Esta ciencia hace uso de diferentes herramientas de carcter investigativo
llegando a incrementar el nivel de conocimiento sobre el comportamiento del
suelo ante cargas dinmicas a travs del desarrollo de ensayos tanto en
laboratorio como en campo.

(Richart & Woods, 1970) en su libro de Vibraciones de Suelos y Fundaciones

afirman el proceso de obtener valores representativos para las propiedades
criticas del suelo es probablemente la parte ms difcil del estudio de diseo.
Esta dificultad surge porque las propiedades crticas del suelo son
dependientes de distintos parmetros. Por ejemplo para un problema dado
dichas propiedades podran variar en factor de 10 en un depsito de suelo
que parece relativamente homogneo por su dependencia de parmetros
como la amplitud de deformacin., el estado de esfuerzo efectivo y el nmero
de ciclos de carga.

Se puede decir que un ejemplo de aplicacin de lo mencionado

anteriormente son los ensayos de laboratorio, se destaca entre estos, el
ensayo de columna resonante en el cual se obtienen los parmetros
dinmicos del suelo (mdulo cortante y relacin de amortiguamiento). Se han
realizado trabajos enfocados en esta temtica, uno de ellos y que ha servido
como gua en la elaboracin de este trabajo es la tesis doctoral Evaluation of
the Small-StrainStiffness of Soilby Non-Conventional DynamicTesting
Methods realizada por el Ingeniero (Camacho-Tauta, 2011), quien es el tutor
de la presente propuesta de trabajo de grado.

5
Por otro lado, uno de los registros ms destacados de medicin de
propiedades dinmicas del suelo se han hallado a travs de Bender
Elements, ensayo implementado a partir de 1978 con el fin de obtener el
mdulo de corte a pequeas deformaciones (Shirley & Hampton, 1978) con
el que inicialmente se determinaban las ondas de corte en el suelo, su
respectiva frecuencia de resonancia y de igual modo la velocidad de llegada
de cada una de ellas. Aspectos como este impulsaron el estudio y la
profundizacin del tema, motivo por el cual el Ingeniero (Guzmn A. , 1986)
interviene en el desarrollo del equipo de Columna Resonante y afirma que
este ensayo permite la determinacin de las propiedades dinmicas del
suelo en un espcimen de cilindro macizo o hueco en toda la gama de las
amplitudes de deformacin (Guzmn A. , 1986)

Existen otros investigadores que han realizado grandes aportes sobre

Dinmica de Suelos, pero no se debe dejar de lado el uso de elementos
fundamentales para la ejecucin de cualquier ensayo, en este caso
el Standard Test Methods for Modulus and Damping of Soils by the
Resonant-Column Method, ASTM International, (2002), aquella norma que
expone de manera especfica el procedimiento a seguir para la ejecucin del
ensayo a desarrollar en muestras de suelos inalteradas o modificada. Los
posibles sometimientos a esfuerzos en especmenes de suelo fueron
estructurados en el ASTM D 1415 donde de manera resumida se concluye
que el espcimen es sometido a carga axial dentro de una cmara que
posteriormente genera presiones alrededor de la muestra. La posibilidad de
controlar otras condiciones como la presin de poros, el grado de saturacin
y la temperatura se manifiestan de forma explcita en esta norma.

6
2 PROPIEDADES DINMICAS DEL SUELO

2.1 Introduccin

La dinmica de suelos es una extensa disciplina que estudia distintas

temticas tales como: comportamiento y respuesta del suelo durante la
aplicacin de carga, principales propiedades dinmicas (mdulo de corte,
relacin de amortiguamiento, deformaciones), comportamiento de ondas,
principalmente. Es importante saber que mediante este estudio se obtiene
informacin que contribuye a la ejecucin de proyectos y resolucin de
problemas como el caso de cimentacin de maquinaria o la compactacin
por vibracin, en general alteraciones que sufre el suelo por distintas causas.

Igualmente la dinmica de suelos al evaluar y analizar los campos frecuentes

de Ingeniera Ssmica, relaciona las propiedades dinmicas del suelo
anteriormente mencionadas con procesos y fenmenos indispensables como
lo son las ondas, elemento bsico que se mide y se valora especialmente su
velocidad de corte, por medio de las propiedades.

2.2 Tipos de ondas

Las ondas ssmicas son un tipo de onda elstica consistentes en la

propagacin de perturbaciones temporales del campo de tensiones que
generan pequeos movimientos en un medio.

Las ondas ssmicas pueden ser generadas por movimientos telricos

naturales, los ms grandes de los cuales pueden causar daos en zonas
donde hay asentamientos urbanos. Existe toda una rama de la sismologa

7
que se encarga del estudio de este tipo de fenmenos fsicos. Las ondas
ssmicas pueden ser generadas tambin artificialmente mediante el empleo
de explosivos o camiones vibradores. (Diaz, 2005). Segn su forma de
propagacin las ondas se pueden clasificar as:

2.2.1 Ondas superficiales

2.2.1.1 Ondas Rayleigh

(Rayleigh, 1885) Predijo la presencia de ondas superficiales diseando

matemticamente el movimiento de ondas planas en un espacio semi-infinito
elstico. Las ondas de Rayleigh ilustradas en la Figura 1D causan un
movimiento rodante parecido a las ondas del mar y sus partculas se mueven
en forma elipsoidal en el plano vertical, que pasa por la direccin de
propagacin. En la superficie el movimiento de las partculas es retrgrado
con respecto al avance de las ondas. Segn (Dobrin & Savit, 1988) la
velocidad de las ondas Rayleigh (VR) es menor que la velocidad de las ondas
s (Vs) y es aproximadamente VR = 0,9 x Vs (Klee, Griem, & Griem, 1999).

2.2.1.2 Ondas de Love

(Love, 1911) Descubri la onda superficial, que lleva su nombre estudiando

el efecto de vibraciones elsticas a una capa superficial. Las ondas de Love
(ver Figura 1C.) requieren la existencia de una capa superficial de menor
velocidad en comparacin a las formaciones subyacentes o es decir un
gradiente de velocidad positivo (velocidad se incrementa) con la profundidad.
Las ondas de Love son ondas de corte, que oscilan solo en el plano
horizontal, es decir las ondas de Love son ondas de cizalla horizontalmente
polarizadas.

8
2.2.2 Ondas de cuerpo

2.2.2.1 Ondas de compresin

Llamadas tambin ondas longitudinales u ondas "p. Bajo este tipo de onda,
las partculas de una onda p, longitudinal o de compresin oscilan en la
direccin de propagacin de la onda. Las ondas p, Figura 1A, son parecidas
a las ondas sonoras ordinarias. Las ondas p son ms rpidas que las ondas
s o es decir despus un temblor en un observatorio primeramente llegan las
ondas p, secundariamente las ondas s.

2.2.2.2 Ondas transversales

Conocidas tambin como ondas de cizalla u ondas "s". Bajo este tipo de
onda, las partculas oscilan perpendicularmente a la direccin de
propagacin. Se distingue las ondas sh, cuyas partculas oscilan en el plano
horizontal y perpendicular a la direccin de propagacin, y las ondas sv,
cuyas partculas oscilan en el plano vertical y perpendicular a la direccin de
propagacin. En las ondas s polarizadas sus partculas oscilan en un nico
plano perpendicular a su direccin de propagacin. (Kramer, 1996). Figura
1.B

Figura 1. Tipos de ondas (Kramer, 1996)

9
2.3 Parmetros dinmicos del suelo

El comportamiento de los suelos bajo cargas dinmicas ha llamado la

atencin notablemente en los ltimos aos, muestra de esto son los anlisis
de situaciones en donde distintas muestras de suelo son estudiadas. Prueba
de estos fenmenos es la alteracin en tres tipos de casos: alteracin de
esfuerzo cclico en los depsitos de suelo, alteracin de esfuerzo en la
subrrasante y la alteracin de esfuerzos cclicos causados por ondas de mar.

Estos ambientes que dependiendo de sus caractersticas se puede decir que

son naturales o artificiales y al presentar cargas cclicas, se relacionan con
parmetros importantes que facilitan y contribuyen a obtener un mejor
entendimiento del tema. A continuacin se darn a conocer algunas de las
caractersticas de dichos parmetros:

2.3.1 Mdulo de corte

El mdulo de corte (G) es un parmetro ingenieril importante en un amplio

rango de aplicaciones geotcnicas como es el caso de: diseo de
fundaciones, predicciones de asentamiento, mejorar el control del suelo y la
evaluacin de licuefaccin. Este mdulo describe el comportamiento del
suelo bajo cargas de trabajo y en movimientos pequeos y medios causados
por sismos.

El mdulo de corte se puede definir como la relacin entre la tensin de corte

() y la deformacin unitaria de corte (). La respuesta que presenta un suelo
frente a las deformaciones se puede obtener en un diagrama de tensin-
deformacin, en un material lineal la grfica que se presenta es una lnea
recta en tanto que un material no lineal genera una curva y su mdulo de

10
corte vara a lo largo de esta curva. Por consiguiente se definen varios
mdulos de corte: mdulo de corte mximo, mdulo de corte secante y
mdulo de corte tangente, representados en la Figura 2 y Figura 3
(Santibaez, 2006):

Figura 2. Diagrama tensin-deformacin en un material lineal (Santibaez, 2006).

Figura 3. Diagrama tensin-deformacin en un material no lineal con sus mdulos de corte (Santibaez, 2006)

Se han desarrollado numerosas tcnicas para investigar el mdulo tanto en

campo como en laboratorio, en general esto se puede dividir en diversas
categoras principales como los son los ensayos destructivos y no
destructivos y los ensayos de carga dinmica o monotnicos (Ibrahim,
Orense, Pender, Tai, & Kikkawa, 2011).

11
2.3.1.1 Mdulo de corte mximo

Clculos de deformaciones, diseo de fundaciones sujetas a cargas

dinmicas y evaluacin del potencial de licuacin son algunos de los
fenmenos de los cuales se encarga este parmetro, se halla mediante
mediciones de velocidad de propagacin de ondas de corte Vs aunque este
mtodo presenta limitaciones y se relaciona con la densidad del suelo s
(Narsilio & Santamarina, 2001).

Ecuacin 1

Donde:

Mdulo de corte mximo

Densidad del suelo
Velocidad de ondas de corte

2.3.1.2 Curvas de reduccin del mdulo

Es comn encontrar la representacin del mdulo de corte versus el mdulo

de corte mximo, conocida como curva de reduccin de mdulo donde se
observa que el mdulo de corte disminuye a medida que aumenta el nivel de
deformacin. La relacin de amortiguamiento tambin depende del nivel de
deformacin aumentando a medida que la deformacin de corte incrementa
(Universidad Politcnica de Catalunya, 2004). Figura 4.

12
 Figura 4. Curva de reduccin de mdulo de corte y variacin de la razn de amortiguacin con el nivel de
deformacin de corte para una arcilla blanda (Universidad Politcnica de Catalunya, 2004).

Las propiedades dinmicas del suelo se definen mediante la relacin de

amortiguamiento y las curvas de reduccin del mdulo principalmente, por
esta razn (Seed & Idriss, 1970) crearon este tipo de curvas para el caso de
arenas, arcillas y gravas. Se muestran en la Figura 5 tres tipos de curvas en
cuanto a lmite inferior, lmite superior y valores promedio para arenas.

Figura 5. Curvas de reduccin del mdulo de corte para arenas (Bala & Balan, 2006)

13
2.3.2 Relacin de amortiguamiento

Las propiedades de amortiguamiento juegan un papel importante en el

estudio de la dinmica de suelos y especficamente en la respuesta que dan
los distintos tipos de suelo frente a un sismo o algn otro fenmeno. Esta
relacin es bsicamente el resultado de dos factores principales:

Las propiedades de viscosidad y plasticidad del suelo

Amortiguamiento por radiacin en el que la energa de vibracin es
liberada con propagacin de ondas.

La razn de amortiguamiento expresa la capacidad del material para disipar

la energa y se calcula mediante la siguiente expresin (Universidad
Politcnica de Catalunya, 2004):

Ecuacin 2

Donde:

Relacin de amortiguamiento
Energa disipada
Energa de deformacin mxima
rea del lazo del ciclo de histresis
Mdulo de corte secante
Deformacin

Esta relacin se encuentra influenciada por la presin de confinamiento

efectiva particularmente en suelos de baja plasticidad. (Ishibashi & Zhang,
1993) Desarrollaron una expresin emprica para esta relacin en suelos

14
plsticos y no plsticos, con esta se puede computar el factor de reduccin
del mdulo, G/Gmax y por consiguiente la relacin de amortiguamiento est
dada por (Kramer, 1996):

[ ( ) ] Ecuacin 3

Donde:

Relacin de amortiguamiento
Mdulo de corte
Mdulo de corte mximo

Por otro lado el modelo visco-elstico es usado para representar las

propiedades de esta relacin en un suelo determinado. En el modelo de
Voigt la relacin de amortiguamiento incrementa con la frecuencia de
vibracin, mientras que en el modelo de Maxwell esta se reduce. El valor de
esta relacin para los diferentes tipos de suelo vara en el rango de 0.3 a 10.
(Bharat, 2009).

En cuanto a la disipacin de energa no histertica que toma lugar en las

deformaciones, se encuentran debajo del umbral de corte-deformacin
cclico lineal. Experimentalmente se ha observado que algo de esta energa
es disipada en niveles de deformacin muy baja, por consiguiente la relacin
de amortiguamiento es cero. Por encima del umbral la amplitud de los bucles
de histresis generada por suelos cargados cclicamente indica que la
relacin de amortiguamiento incrementa cuando la amplitud de deformacin
aumenta.

15
Las relaciones de amortiguamiento de suelos altamente plsticos son ms
bajas que las de suelos con plasticidad baja en la misma amplitud de
deformacin cclica, esto se puede evidenciar en la siguiente figura:

Figura 6. Variacin de la relacin de amortiguamiento de suelos granulares finos con amplitud cclica de corte-
deformacin y plasticidad ndice. Efecto de la plasticidad de suelos en respuesta cclica (Kramer, 1996).

Se observa que la curva de amortiguamiento en PI=0 (PI: ndice de

plasticidad) es idntica a la de amortiguamiento promedio que es usada para
los suelos de granos gruesos cuando son trabajos separadamente de los
suelos con granos finos. Este hecho sugiere que las curvas de
amortiguamiento de la figura anterior se pueden aplicar a ambos tipos de
suelo: de granos finos y de granos gruesos. El comportamiento de la relacin
de amortiguamiento en las gravas es similar al de las arenas.

16
2.3.3 Variacin del comportamiento del suelo a lo largo de las
deformaciones

La necesidad de evaluar las caractersticas de esfuerzo-deformacin sobre

un amplio rango de amplitud de corte-deformacin ha simulado el desarrollo
de una variedad de tcnicas de medicin tanto en laboratorio como en
campo. Sin embargo se ha logrado un progreso en cuanto a equipos de
ensayo, en la interpretacin de lmites, efectos de membranas y en el
entendimiento de capacidades y limitaciones de cada tcnica; la
comparacin entre predicciones basadas en resultados de laboratorio y la
realizacin en campo podra estar ser decepcionante. Esto se puede atribuir
primero al gran nmero de factores que controlan el comportamiento
mecnico de materiales granulares que hacen difcil la estabilizacin de las
relaciones de esfuerzo-deformacin.

Se afirma que la densidad y el estado de esfuerzo no son suficientes para

caracterizar el comportamiento de los materiales granulares. Otros factores
como estructura y su asociada anisotropa, historia de la relacin esfuerzo-
deformacin, envejecimiento y la trayectoria de esfuerzo, entre otros se
deben tener en cuenta, sin embargo se ha observado que ni la estructura ni
la presin de confinamiento parecen afectar el comportamiento cclico de
presin de poros de arenas y ensayos de deformacin controlados. (Alarcn,
Chameau, & Leonards, 1986).

Cuando el comportamiento del suelo permanece dentro del rango de

pequeas deformaciones, el uso de un modelo elstico es necesario y el
mdulo de corte es un parmetro clave para obtener correctamente un
modelo del comportamiento del suelo, esto se puede evidenciar en la
siguiente figura:

17
 Figura 7. Modelo del comportamiento del suelo en conformidad con las caractersticas de deformacin
(Kramer, 1996)

Cuando un problema dado est asociado con el rango medio de

deformacin, aproximadamente por debajo del nivel de 10 -3 el
comportamiento del suelo llega a ser elastoplstico y el mdulo de corte
tiende a disminuir en tanto aumenta la relacin de corte-deformacin, al
mismo tiempo la energa de disipacin en suelos es principalmente
independiente de la velocidad y de la naturaleza histertica y la relacin de
amortiguamiento puede ser usada para representar las propiedades
absorbentes de energa de los suelos.

Por lo general la relacin esfuerzo - deformacin es lineal aunque es

inelstica en el caso de deformaciones muy pequeas, esta regin de
comportamiento lineal es caracterizada por valores de deformacin dentro
del rango 0<<t donde el lmite superior t ha sido catalogado como el ciclo
lineal del umbral de deformacin-corte por Vucetic (1994). Dentro de esta
regin la no degradacin de rigidez es observada en la respuesta del suelo,

18
sin embargo el ciclo de histresis, en el plano esfuerzo-deformacin, muestra
un rea no nula.

El fenmeno de disipacin de energa en niveles de deformacin pequea es

causado por la existencia de un ensayo controlado de deformacin, de un
retraso de tiempo entre un ciclo de deformacin dado y un ciclo de esfuerzo,
la palabra histresis proviene de Grecia y significa retraso o demora. Este
tiempo es el responsable de la prdida de energa sobre un periodo finito de
tiempo que es tpico desde un punto de vista fenomenolgico de un
comportamiento viscoelstico, como se puede observar en la figura.

Figura 8. Ciclo de histresis debido a cargas armnicas en un cuerpo viscoelstico (Camacho-Tauta, 2011)

Los datos comunes de experimentacin no apoyan la existencia del

fenmeno de disipacin de energa instantneo para <t, este sera el
comportamiento tpico de un fenmeno elastoplstico.

De acuerdo con varios autores t en la prctica est sobre 10-3 % y

depende de ambos factores externos e intrnsecos hallados
experimentalmente. Otro factor importante del comportamiento del suelo en

19
pequeas deformaciones observado experimentalmente es la invariancia de
importantes propiedades del suelo como la rigidez, el coeficiente de
amortiguamiento y el nmero de avances de ciclos. Este hecho se evidencia
en la forma del ciclo de histresis porque no se usa para cambiar el
incremento de ciclos. (Presti & Turco).

2.3.4 Esfuerzo efectivo principal y relacin de vacos

Los efectos de diferentes factores sobre la velocidad de onda de corte en

arenas mediante ensayos de columna resonante fueron analizados por
(Hardin & Richard, 1963) que ms adelante computaron dichos datos con
curvas de regresin de la siguiente forma:

Ecuacin 4

Esta investigacin revel que dos factores son los que ms influyen en la
velocidad de onda de corte, son: el esfuerzo normal efectivo principal y la
relacin de vacos. Al remplazar la densidad de masa del suelo se obtiene:

( ) Ecuacin 5

El mdulo de corte inicial se puede expresar como:

Ecuacin 6

20
Al simplificar se obtiene:

Ecuacin 7

En donde:

Constantes empricas determinadas experimentalmente

= Relacin de vacos
Esfuerzo efectivo principal

En la Figura 9 se observa un ejemplo de los resultados hallados por (Hardin

& Richard, 1963).

Figura 9. Variacin del mdulo de corte con relacin de vacos y presin de confinamiento para formas
redondeadas y angulares secas (lneas continuas) y arenas de grano (lneas discontinuas) (Camacho-Tauta,
2011).

21
Al comparar dos tipos de suelos se evidencia que las constantes A y B
detienen la influencia del tipo del grano dentro del mdulo de corte. El valor
del exponente n es igual a 0.5 para los dos tipos de arena por consiguiente la
ecuacin anterior se puede expresar como:

Ecuacin 8

Donde es una funcin de la relacin de vacos. Al

analizar la funcin de la figura 11 se puede afirmar que el mdulo de corte
inicial aumenta con el esfuerzo efectivo principal y disminuye con la relacin
de vacos.

Ms tarde (Ishihara, 1996) reuni resultados experimentales que fueron

complementados por (Kokusho, 1981) con datos de distintos autores de
diferentes mtodos de ensayo para arenas, arcillas y gravas. La siguiente
tabla representa los datos disponibles de referencias citadas y otras
publicaciones recientes.

Tabla 1. Constantes en ecuaciones empricas propuestas en mdulos de deformacin pequea (Camacho-

Tauta, 2011)

Metho
Reference A F(e) n Soil Material
d*
Hardin and Round grained Ottawa sand
7000 0.5 RC
Richart (1963) [ ]
Angular grained crushed
Hardin and quartz
3300 0.5 RC
Richart (1963) [ ]
[ ]
Iwasaki and Eleven kinds of clean sand
1410
Tatsuoka 0.4 [ ] RC
0
(1977) [ ] kPa
Toyoura sand
Kokusho
8400 0.5 [ ] CTX
(1980)
[ ] kPa

22
 Three kinds of clean sand
Yu and Richart
7000 0.5 [ ] RC
(1984)
[ ] kPa
Silty sand. Fines = 41%: PI = 9
Gomes Correia
5000 0.5 [ ] RC
et al. (1993)
[ ]
Sand D50 = 0.65 mm: Ut =2.4
Gomes Correia
9000 0.5 [ ] RC
et al. (1994)
[ ]
Altun and Toyoura sand D50 = 0.17 mm;
Goktepe 5937 0.5 Uc =1.43; [ ] CTT
(2006) [ ] kPa
Low plasticity clays
Hardin and
3300 0.5 [ ] RC
Black (1968)
[ ] kPa
Kaolinite. PI = 35
Marcuson and
4500 0.5 [ ] RC
Wahls (1972)
[ ] kPa
Bentonite. PI = 60
Marcuson and
450 0.5 [ ] RC
Wahls (1972)
[ ] kPa
Undisturbed clays, PI = 40-85
Kokusho et
141 0.6 [ ] CTX
al.(1982)
[ ] kPa
Heelis and
Peat [ ]
Gomes Correia 350 0.5 RC
[ ] kPa
(1998)
Silty clay. PI = 40-42
Santos et al.
520 0.5 [ ] RC
(1999)
[ ] kPa
Ballast, D50 = 40 mm;
0.3
Prange (1981) 7230 Uc =3.0; [ ] RC
8
[ ] kPa
Kokusho and 1300 0.5 Crushed rock, D50 = 30 mm;
CTX
Esashi (1981) 0 5 Uc =10; [ ]

Kokusho and Round gravel, D50 = 10 mm;

8400 0.6 CTX
Esashi (1981) Uc =20; [ ]
* RC: Resonant-column; CTX: Cyclic triaxial; TX: Triaxial, CTT: Cyclic torsional

23
2.4 Medicin de las propiedades dinmicas del suelo

La medicin de propiedades dinmicas del suelo es una tarea importante en

la solucin de problemas relacionados con Geotecnia y con sismos
especficamente, para lograr dichas soluciones la Ingeniera ha desarrollado
distintas tcnicas de campo y de laboratorio. Muchas de estas tcnicas estn
orientadas hacia la medicin de propiedades con bajas deformaciones y
otras hacia propiedades con deformaciones ms largas.

Las propiedades del suelo que intervienen en la propagacin de las ondas y

otros fenmenos de bajas deformaciones incluyen rigidez, amortiguamiento,
relacin de Poisson y densidad; la rigidez y el amortiguamiento son las ms
importantes, ya que son evaluadas crticamente en muchos problemas de
ingeniera que estudian especficamente sismos, no solamente se estudian
en bajas deformaciones sino que tambin se evalan en deformaciones
intermedias y altas, esto se debe a que los suelos son materiales no lineales,
las otras propiedades tienen menos influencia y tienden a poseer rangos
estrechos relativamente.

En niveles de deformacin alta la influencia de la tasa y el nmero de ciclos

de carga en resistencia al corte son aspectos importantes al igual que las
caractersticas de cambio de volumen. La gran mayora de los ensayos de
campo y de laboratorio desarrollados se han encargado principalmente de la
medicin de las propiedades del suelo, tambin se han creado otros que son
versiones modificadas de ensayos usados para medir el comportamiento del
suelo bajo condiciones de carga monotnicas.

24
2.4.1 Ensayos de campo

2.4.1.1 Ensayo de reflexin ssmica

Permiten determinar la velocidad de propagacin de la onda y el espesor de

las capas superficiales a partir de la superficie del suelo o de construcciones
realizadas en el ocano, tanto el ensayo como su interpretacin son
conceptualmente muy simples. El ensayo de reflexin Figura 10, es til para
investigacin a grandes escalas y para estratigrafas de ondas muy
profundas. El ensayo es raramente usado para delineacin de capas de
suelo superficiales.

Figura 10. Ensayo de reflexin ssmica (Universidad Politcnica de Catalunya, 2004)

25
2.4.1.2 Ensayo de refraccin ssmica

Elimina la limitacin ms importante del ensayo de reflexin ssmica

mediante el uso de tiempos de llegada de las primeras ondas,
independientemente de la ruta y as alcanzar un receptor determinado. El
ensayo contiene medicin de los tiempos de viaje de las ondas p y/o s de
una fuente de impulso a una matriz lineal de puntos a lo largo de la superficie
de la tierra a distancias diferentes de la fuente. Este ensayo, Figura 11, es
ms usado que el ensayo de reflexin debido a que su aplicacin es ms
grande en ingeniera de sismos por su delineacin en la mayora de unidades
estratigrficas.

Gefono 1 Gefono 2 Gefono 3 Gefono 4

Trigger
Estrato 1

Estrato 2

Figura 11. Ensayo de refraccin ssmica (Rinaldi)

2.4.1.3 Ensayo ssmico de cross hole

Usado en dos o ms perforaciones para medir velocidades de propagacin

de onda a lo largo de rutas horizontales. El ensayo ms simple de Cross
Hole, Figura 12, consiste en dos perforaciones, una contiene una fuente de
impulso de energa y la otra un receptor, fijando ambos (fuente y receptor) a
la misma profundidad en cada perforacin, se mide la velocidad de
propagacin de onda del material entre las perforaciones a esa profundidad.

26
Al ensayar a varias profundidades se puede obtener un perfil de velocidad,
se desea por lo general el uso de ms de dos perforaciones para evitar
posibles imprecisiones generadas por mediciones de tiempo desde el
momento en que se inicia la perturbacin producida por el disparo que da la
orden de inicio, efectos y el sitio de anisotropa.

Este tipo de ensayo suministra datos de velocidad confiables usando fuentes

de impulso mecnicas y datos de profundidades ms grandes con fuentes
explosivas. (Kramer, 1996).

Figura 12. Mtodo de cross hole (Cedex, 2010)

2.4.1.4 Ensayo down hole

Est enfocado fundamentalmente a detectar la velocidad de propagacin de

las ondas tangenciales SH (vibracin de las partculas en la direccin
horizontal) en su trayectoria vertical a travs de una masa de suelo. El
ensayo consiste en generar ondas tangenciales mediante una fuente de

27
energa reversible que produzca un efecto de cizalla horizontal en la
superficie del terreno, y en registrar la llegada del impacto a lo largo de un
sondeo situado a una distancia comprendida entre 2 y 5 m del foco emisor
(Figura 13).

Al estar la fuente de energa cerca del sondeo receptor, las ondas que
interesa detectar se propagarn en una direccin muy prxima a la vertical
por lo que es lcito suponer una trayectoria directa desde el foco emisor a los
receptores para calcular la velocidad de propagacin de las ondas. Frente al
ensayo cross-hole presenta la ventaja de que tan slo exige la perforacin de
un sondeo para llevar a cabo los ensayos, y el inconveniente de que la
energa ssmica que se puede generar mediante el golpe de un martillo se
amortigua rpidamente en el terreno alcanzando una profundidad mxima de
15 a 20 m. (Cedex, 2010)

Figura 13 Ensayo down hole (Cedex, 2010)

28
2.4.2 Ensayos de laboratorio

2.4.2.1 Ensayo de columna resonante

El ensayo de la Columna Resonante, Figura 14, es un mtodo de laboratorio

para conocer las propiedades dinmicas de los depsitos de suelo, el cual
cubre la determinacin del mdulo de corte, las caractersticas de
amortiguamiento, la velocidad de onda de corte y el mdulo de Young, para
muestras cilndricas de suelo en condiciones alteradas o inalteradas. Este
mtodo est normalizado bajo la designacin D4015-92 de la Sociedad
Americana de Ensayos y Materiales (ASTM por sus iniciales en ingls).

El ensayo de la columna resonante consiste fundamentalmente en someter

muestras cilndricas macizas o huecas de suelo a cargas torsionales o
longitudinales armnicas por medio de un sistema electromagntico. (Alva,
2010).

Es posible desarrollar diferentes versiones de ensayos de columna resonante

utilizando diferentes condiciones de fijacin de los extremos y de
confinamiento de la muestra. Para este caso especifico el sistema de fijacin
de la muestra es fijo - libre, es decir que el sistema est fijo en un extremo y
libre en el otro y el confinamiento es isotrpico. Los ensayos se hicieron en
muestras cilndricas macizas de arcilla y la excitacin se realiz de tipo
torsional.

Equipos: Para la realizacin del ensayo se requiere de un mecanismo que

consta de:

- Los platos para fijacin y transmisin del movimiento a la muestra

- La celda de presin
- El sistema excitador compuesto por bobinas e imanes

29
 - El transductor para medir la respuesta
- La instrumentacin (Moreno & Francisco, 1998)

Figura 14 Ensayo de Columna Resonante (Moreno & Francisco, 1998)

2.4.2.2 Ensayo de pulso ultrasnico

Las velocidades de propagacin de onda pueden ser medidas en el

laboratorio por mtodos de ensayos de pulso ultrasnico. Los transmisores y
receptores se renen en platinas que pueden ser puestas al final de un
espcimen con una distancia que las separa. Los transmisores y receptores
son hechos de materiales piezoelctricos que presentan cambios en las
dimensiones cuando se someten a voltaje a travs de sus caras, un pulso
elctrico de alta frecuencia es aplicado al transmisor y causa rpidamente

30
una deformacin produciendo una onda de esfuerzo que viaja a travs del
espcimen hacia el receptor.

Cuando la onda de esfuerzo alcanza el receptor, este genera un pulso de

voltaje que es posteriormente medido, la distancia entre el transmisor y el
receptor es dividida por el tiempo de diferencia entre los pulsos de voltaje, de
esta forma se obtiene la velocidad de propagacin de la onda. El ensayo de
pulso ultrasnico ilustrado en la Figura 15, es usado por lo general para
materiales muy suaves como sedimentos de suelo de mar. (Kramer, 1996).

Figura 15 Ensayo de Pulso Ultrasnico (Lombillo & Villegas, 2010)

2.4.2.3 Ensayo de elementos bender piezoelctricos

Es un ensayo que sirve para medir la velocidad de onda de corte en

especmenes de laboratorio, estos elementos son construidos por la unin de
dos materiales piezoelctricos, all se aplica un voltaje especficamente en
las caras causando expansin en unas y contraccin en otras, generando
finalmente que el elemento se doble. Igualmente una perturbacin lateral
del elemento piezoelctrico producir un voltaje, entonces dicho elemento

31
puede ser usado como transmisor y receptor de onda. (Kramer, 1996). La
siguiente figura indica parte del montaje del ensayo:

Figura 16 Ensayo de Elementos Bender (Sandox, 2010)

2.4.2.4 Ensayo triaxial cclico

El ensayo cclico se ejecuta en una celda similar a la del ensayo triaxial con
excepcin de que la parte alta de la celda tiene un pistn ms largo que
puede transmitir cargas cclicas sobre la parte superior de la muestra a
ensayar. La Figura 17 muestra el equipo del ensayo que tiene que aplicar
una carga uniforme sinusoidal a una frecuencia entre 0.1 y 2 ciclos por
segundo (Hz). La frecuencia ms comn utilizada es de un ciclo por
segundo. Las mediciones de carga axial se realizan mediante una celda
electrnica y las presiones de poro se monitorean con un equipo de alta
precisin capaz de medir variaciones de +0.25 psi.

El objetivo principal de un programa de ensayo de triaxial cclico es el

determinar el potencial de licuacin de los suelos a carga o deformacin
controlada, en probetas de 4 6 de dimetro y 6 7 de altura,
respectivamente. Consta, bsicamente, de cuatro elementos:

32
1. Sistema hidrulico (aplicacin de carga).
2. Panel de control.
3. Sistema de aplicacin de presin (confinante y contrapresin).
4. Sistema automatizado del equipo, captura y procesamiento de datos.
(Castrellon, Gmez , Romo , & Flores)

Figura 17 Ensayo Triaxial Cclico (Concha, Poblete, & Mauro, 2004)

33
3 ENSAYO DE COLUMNA RESONANTE

3.1 Introduccin

Una de las pruebas de laboratorio ms conocidas en la actualidad encargada

de evaluar las propiedades dinmicas de los suelo es el ensayo de columna
resonante; la determinacin del mdulo de Young, del mdulo de corte y sus
respectivas relaciones de amortiguamiento pueden ser determinadas segn
el esquema de aplicacin de cargas que haya sido establecido, ya sea
aplicacin de fuerza axial o torque.

El ensayo consiste en el sometimiento de una probeta de suelo a una

excitacin torsional de tipo sinusoidal a frecuencia de resonancia. (Camacho
- Tauta, 2008). Este mtodo est normalizado bajo la designacin D4015-92
(2000) Standard Test Methods for Modulus and Damping of Soils by the
Resonant-Column Method de la Sociedad Americana de Ensayos y
Materiales (ASTM por sus iniciales en ingls).

Con el equipo de columna resonante es posible desarrollar diferentes

mtodos de ensayo una vez se haya modificado y establecido las nuevas
caractersticas de excitacin a las cuales ser sometida la muestra de suelo,
esto con el fin de analizar las propiedades dinmicas a diferentes frecuencias
de resonancia. (Camacho-Tauta, 2011).

Esta seccin inicia con una breve descripcin del equipo de columna
resonante, los diferentes montajes que se pueden llevar a cabo segn las
caractersticas de excitacin o el tipo de vibracin generadas a las muestras,

34
adicionalmente se hace una breve mencin de las ecuaciones fundamentales
necesarias para el clculo de las propiedades dinmicas del suelo; y
finalmente, una explicacin de los posibles mtodos de ejecucin del ensayo
de columna resonante.

3.2 Equipo de columna resonante

Como se ha mencionado anteriormente, el equipo de columna resonante

trabaja bajo dos condiciones especiales, aplicacin de fuerza axial o
induccin de torque, esto, segn el tipo de propiedades dinmicas de inters
a analizar. La Figura 18 presenta un esquema grafico del sistema que
conforma el equipo mencionado.

Para determinar el mdulo de Young y su relacin de amortiguamiento

asociada es necesario que el sistema actu bajo carga axial, de otro modo,
el mdulo de cortante y su respectiva relacin de amortiguamiento se
obtienen cuando el sistema le genera torque al espcimen de estudio.

Se conoce variedad de configuraciones para el equipo de columna resonante

fundamentadas en las condiciones de frontera de los extremos de las
muestras que sern ampliamente explicadas en el siguiente tem.

En 1937 (Woods, 1978) atribuye la primera aplicacin del mtodo a los

ingenieros japoneses Ishimoto e Ida. (Camacho-Tauta, 2011). Veinte aos
despus resurge la tcnica en manos de (Shannon, 1959); ao desde el cual
se ha venido generando grandes mejoras a los procesos desarrollados
desde tiempos atrs y que actualmente continan hacindolo.

35
Problemas relacionados con la no uniforme distribucin de la deformacin a
los largo de la longitud de la muestra y de la tensin de corte a travs del
radio de la muestra, fueron resueltos por (Drevinch, 1967) mediante la
adicin de una masa en la parte superior de la muestra y recurriendo a
especmenes huecos, respectivamente. (Camacho - Tauta, 2008).

Importantes trabajos realizados han permitido desarrollar mejoras

significativas al mtodo analizado, como lo es la incursin del anlisis en el
dominio de la frecuencia en el ensayo de columna resonante.

Otro aspecto motivo de investigacin, tiene que ver con la medicin

adecuada de la excitacin de entrada, ya que durante el ensayo, se produce
una fuerza electromagntica opuesta a la fuerza aplicada por las bobinas, lo
que resulta en un amortiguamiento adicional generado por el equipo que
causa errores en el clculo de la relacin de amortiguamiento del suelo.
(Camacho - Tauta, 2008).

Figura 18. Sistema del equipo de columna resonante

36
Tipos de montajes de columna resonante

Existe la posibilidad de llevar a cabo el ensayo de columna resonante

haciendo uso de diferentes montajes fundamentados en las condiciones
fsicas a las que pueden ser sometidas las muestras en sus extremos.
Caractersticas basadas en la aplicacin de cargas o excitaciones,
confinamiento del espcimen, presin de poros y cambios volumtricos se
pueden ver modificados a razn de la implementacin de alguno de los
posibles tipos de montajes que se pueden ejecutar

- Sistema fijo-libre

El primero de los casos a describir cosiste en un montaje fijo-libre en donde

el especimen no es sometido a rotacin debido a que la muestra est unida
al pedestal tal como se puede apreciar en la Figura 19.

Figura 19. Montaje fijo-libre de columna resonante (Camacho-Tauta, 2011)

37
En este sistema el pedestal se encuentra rgidamente anclado y tiene una
inercia rotacional de masa (JP) 100 veces mayor a la inercia rotacional de
masa de la muestra de suelo, as mismo la distribucin de ngulo de rotacin
a lo largo de la altura z obedece a una funcin sinusoidal. Para el primer
modo de vibracin, la rotacin que se corresponde a de la onda sinusoidal;
la mxima rotacin alcanzada se presenta en la parte superior de la muestra
con una inercia rotacional de masa (J). (Camacho-Tauta, 2011)

- Sistema con extremo activo

El segundo tipo corresponde a un montaje con un extremo activo, en la parte

superior del espcimen se localiza una platina con inercia rotacional de masa
(JA), tal como se indica en la Figura 20, (Camacho-Tauta, 2011), generando
una distribucin de la rotacin casi lineal en la totalidad de la altura de la
muestra de suelo. De esta forma se generan pruebas en donde se trata de
mantener el estado natural de las muestras, es decir, se logra que la tensin
de corte calculada concuerde con las tensiones de corte locales.

Figura 20. Montaje extremo activo con platina en la parte superior (Camacho-Tauta, 2011)

38
La parte superior de la masa cuenta con un motor de accionamiento
empleado para generar la torsin al espcimen. (Camacho-Tauta, 2011).

- Sistema modificado

La Figura 21 hace referencia a una modificacin del sistema fijo-libre de

extremo activo y con una platina en la parte superior de la muestra, explicado
anteriormente.

Figura 21. Modificacin del sistema fijo-libre con platina en la parte superior de la muestra (Camacho-Tauta,
2011)

(Camacho-Tauta, 2011) Expone que este sistema cuenta con un resorte y

amortiguador cuyas constantes estn dadas por KA y A respectivamente
generando de esta forma rigidez torsional y amortiguamiento al sistema. Esta
modificacin genera una mejor distribucin del ngulo de rotacin mientras el
sistema aplica hasta niveles medios tensin de corte; adicionalmente los
esfuerzos anisotrpico e isotrpicos estn en funcin del actuador que aplica
las cargas axiales.

39
 - Sistema cuasi libre

Propone un sistema ideal en el cual la muestra se encuentra aislada en los

dos extremos y a los lados, convirtindose en un estado difcil de lograr en la
prctica debido a las condiciones generadas por los equipos de
accionamiento como lo es el motor y el confinamiento generado por la
platina en la parte superior de la muestra. Figura 22.

Figura 22. Sistema cuasi libre de columna resonante (Camacho-Tauta, 2011)

(Kalinski & Thummaluru, 2005) Desarrollaron un sistema libre de columna

resonante para llevar a cabo la medicin de pequeas deformaciones con la
aplicacin de altas presiones de confinamiento. (Camacho-Tauta, 2011)

3.3 Ensayo convencional de columna resonante (RCCT)

El ensayo convencional de columna resonante (RCCT) es el resultado

experimental en el cual se logra la frecuencia de resonancia de un sistema
de forma manual, es decir que la aplicacin de la excitacin sinusoidal es con
una frecuencia constante. Es decisin del operador aumentar o no la

40
frecuencia obteniendo un nuevo valor constante con base a la observacin y
medicin del ngulo de fase entre el torque de entrada y la vibracin de
salida.

El torque sinusoidal aplicado al sistema produce una vibracin rotacional, la

cual en general se encuentra desfasada con relacin al torque. Cuando el
sistema est vibrando a la frecuencia de resonancia, las seales de entrada
y salida estn desfasadas 90. La figura de Lissajous, que se puede obtener
en un osciloscopio convencional, sirve para identificar esta condicin.
(Camacho - Tauta, 2008).

En la Figura 23 se puede observar el torque sinusoidal que es aplicado al

sistema (a), la rotacin generada (b) y la grfica del torque de entrada con la
rotacin generada, la cual tiene forma de elipse, figura de Lissajous (c).

Figura 23. Ensayo convencional de columna resonante (Camacho-Tauta, 2011); a) Torque aplicado; b) Rotacin
generada; c) Figura de Lissajous

41
Cuando se logra la condicin de resonancia en el sistema, es posible calcular
la incgnita que hasta el momento es desconocida, el mdulo cortante; para
desarrollar este proceso se emplea la Ecuacin 9 donde

Ecuacin 9

[ ( ) ]
( )

En este caso particular no es necesario el torque ni la rotacin aunque estas

resultan ser indispensables para el clculo de la relacin de amortiguamiento
del suelo llevndose a cabo con la Ecuacin 10.

Ecuacin 10
( ) ( )

Donde:

: Amplitud del torque aplicado sobre la muestra.

: Amplitud de la rotacin resultante.
: Angulo de fase entre rotacin y torque.
: Frecuencia angular de vibracin.
: Inercia rotacional de masa del espcimen.
: Altura del espcimen.
: Mdulo cortante del espcimen.
: Relacin de amortiguamiento del espcimen.
: Densidad de masa del suelo.
: Inercia rotacional de la masa activa del aparato.
: Relacin de amortiguamiento de la masa activa del aparato.
: Frecuencia angular de resonancia de la masa activa del aparato.
:

42
El mtodo de decaimiento libre es una alternativa ms para calcular la
relacin de amortiguamiento; en este sistema, (Camacho-Tauta, 2011), la
seal de entrada se interrumpe bruscamente registrando las vibraciones tal
como se muestra en la Figura 24(a). En la Figura 24 (b) se observa el
decaimiento de la amplitud la cual cumple con una funcin exponencial,
donde cada pico representa la amplitud de cada ciclo.

Figura 24. Determinacin del mdulo de amortiguamiento (Camacho-Tauta, 2011); a) Interrupcin de la seal
de entrada; b) Ciclo

De esta forma y por regresin es posible calcular Ecuacin 11 y

posteriormente la relacin de amortiguamiento Ecuacin 12.

Ecuacin 11

Ecuacin 12

El clculo obtenido por la Ecuacin 12, hace referencia a un sistema total

que incluye la muestra de suelo y la masa activa. El RCCT solo toma los
datos de resonancia para calcular las propiedades dinmicas.

43
De este modo, el montaje se puede tornar un tanto complejo en el momento
en el que se deba establecer la frecuencia de resonancia por parte del
operador, ya que se describe como un proceso manual y de tal forma
demanda mayor disposicin de tiempo.

3.4 Ensayo en el dominio de la frecuencia (RCFD)

El ensayo convencional de columna resonante, como ya se mencion en el

literal anterior, adems de demandar gran cantidad de tiempo es muy
propensa a generar errores en la medicin; esto debido a la cantidad de
ciclos que el operador debe aplicar para obtener la resonancia y que
posteriormente afectan las propiedades dinmicas de la muestra.

Si los datos prcticos generan una funcin de transferencia (el lado izquierdo
de la Ecuacin 10), por medio de un procedimiento numrico efectuado,
como el ajuste por mnimos cuadrados de la informacin en variable
compleja se puede encontrar el mdulo cortante (G) y la relacin de
amortiguamiento ( ), siendo estos los parmetros desconocidos del lado
derecho de la Ecuacin 10. (Camacho - Tauta, 2008)

El modelo matemtico a seguir para el desarrollo del ensayo en el dominio

de la frecuencia consiste en la aplicacin de la Ecuacin 13, en la cual se
toma el promedio de todos los ensayos realizados. En el instante en que se
obtiene la funcin de transferencia se puede proceder al clculo de las
propiedades dinmicas haciendo uso de la Ecuacin 10.

Ecuacin 13

44
Donde:

{ } { } Ecuacin 14

{ } { } Ecuacin 15

: Espectro cruzado entre (t) y T(t)

: Espectro de potencia deT(t)
{} : Transformada de Fourier
: Funcin Conjugada compleja

En este tipo de ensayo se implementa una prueba en la cual se utilizan

ondas de barrido sinusoidal como forma de excitacin; esta prueba de
barrido sinusoidal se caracteriza por tener una frecuencia con variacin lineal
respecto al tiempo; el ensayo se realiza con un valor inicial de frecuencia f start
y un valor final fstop, la sea sinusoidal se expresa de la siguiente forma:

[ ( ) ] Ecuacin 16

Donde:
Amplitud
Periodo de la seal

Este procedimiento se lleva a cabo en dos etapas, en donde inicialmente se

toma un intervalo de onda bastante grande con el propsito de poder
identificar el rango en el cual se encuentra el modo fundamental, una vez
hecho esto, se realiza un nuevo ensayo, en esta ocasin con frecuencias que
contienen el modo fundamental, logrando de esta forma un desarrollo ms
detallado de la resonancia del sistema y de la funcin de transferencia

45
(Camacho-Tauta, 2011); esta funcin de transferencia se obtiene a travs de
la Ecuacin 13.

A continuacin se presenta un ejemplo del ensayo de columna resonante en

el dominio de la frecuencia; en la Figura 25(a) se observa la variacin del
torque en funcin de la frecuencia en la Figura 25(b) la rotacin generada.

Figura 25. Ensayo en el dominio de la frecuencia (Camacho - Tauta, 2008); a) Torque aplicado; b) Rotacin
generada

En la Figura 26(a) y Figura 26(b) se observa la magnitud y el ngulo de fase

respectivamente.

Figura 26. Funcin de transferencia (Camacho - Tauta, 2008); a) Magnitud; b) Angulo de fase

46
3.5 Ensayo de no resonancia (NR)

En este caso en particular, para el ensayo de no resonancia, se puede

remplazar el mdulo cortante y amortiguamiento por el mdulo complejo, el
cual est representado como se muestra a continuacin:

Ecuacin 17

De esta forma la Ecuacin 10 podr ser remplazada por la Ecuacin 18,

aplicando una excitacin armnica con frecuencia constante y amplitud
conocida con el fin de medir la rotacin y ngulo de fase del sistema,
finalmente podr determinarse G* siendo la nica variable desconocida, con
esto se establece que no es necesario determinar las condiciones de
resonancia para conocer las propiedades dinmicas de suelo. (Camacho-
Tauta, 2011).

Ecuacin 18
( ) ( )

Este ensayo consiste en obtener diferentes frecuencias para deformaciones

iguales; de tal modo, se procede a la aplicacin de una excitacin sinusoidal
de baja amplitud, la cual se ir variando para obtener la rotacin o
deformacin deseada, es importante aclarar que dicha deformacin cuenta
con un rango de error el cual debe ser en cuanto se pueda menor.

Cuando se obtiene la rotacin deseada el ensayo se repite, obteniendo de

esta forma un grupo de puntos los cuales se encuentran inmersos en un
intervalo de frecuencias previamente establecido.

47
 Figura 27. Ensayo de no resonancia (Camacho-Tauta, 2011) a) Torque; b) Rotacin

En la Figura 27 se observan los valores obtenidos para el rango de

frecuencias indicado, la Figura 27(b) muestra claramente el intervalo de
frecuencias para el mismo valor de deformacin.

Figura 28.Propiedades dinmicas obtenidas con un ensayo de no resonancia (Camacho-Tauta, 2011) : a)

Mdulo de corte; b) Relacin de amortiguamiento

La Figura 28 muestra las propiedades dinmicas de un ensayo no

resonancia con el cual se puede concluir que tanto el mdulo cortante como
la relacin de amortiguamiento varan con relacin a la frecuencia.

48
4 ADQUISICION Y PROCESAMIENTO DE SEALES

4.1 Introduccin

Las modificaciones en el equipo original logran ser de mayor utilidad en el

momento en se implementan modernas tcnicas de adquisicin y
procesamiento de seales. Por medio de la Figura 29, (Avitabile, 2001)
describe de forma grafica el procedimiento de cmo un espectro es analizado
para la obtencin de la funcin de transferencia.

Figura 29. Pasos para el procesamiento de datos en una Transformada de Fourier Rpida (FFT) por sus siglas en
ingles. (Avitabile, 2001)

49
La adquisicin y procesamiento de seales analgicas se compone de las
siguientes etapas: (Orfanidis, 1996)

1) Digitalizacin de la seal analgica: La seal analgica continua se

mide peridicamente cada segundo ( ) durante un tiempo total
, donde es el numero de muestras. Seguido de esto, cada
muestra se cuantifica e un numero finito de bits. (Conversin A/D lo
cual corresponde a una seal Anlogo/Digital segn sus iniciales)

2) Las muestras digitalizadas son procesadas por un procesador de

seal, el cual desarrolla funciones de filtrado y estimacin del
espectro, entre otras.

3) Las muestras de salida resultantes pueden ser convertidas a forma

analgica nuevamente por analoga de reconstruccin. (Conversin
D/A)

A continuacin se desarrolla una descripcin ms detalladas de los procesos

que intervienen en la digitalizacin y procesamiento de seales.

4.2 Digitalizacin

Esta etapa de digitalizacin esta compuesta por dos operaciones, las cuales
estn definidas como: muestreo y cuantizacin. Inicialmente se define como
muestreo la operacin de reduccin de una seal continua a una seal
discreta donde . Debe cumplirse el teorema de
muestreo de Nyquist-Shannon, el cual establece dos condiciones que se
deben cumplir para una representacin exacta de la seal analgica por un
conjunto de muestras discretas, (Camacho-Tauta, 2011):

50
 a) El espectro de la seal analgica debe ser limitado a cierta frecuencia
mxima, fmax.
b) La velocidad de muestreo debe ser elegida de tal forma que sea al
menos el doble de la frecuencia mxima de la seal analgica:

Ecuacin 19

El cociente recibe por nombre: frecuencia de Nyquist o frecuencia de

plegado y define el punto de corte de frecuencias en las operaciones de
procesamiento de seales.

(Camacho-Tauta, 2011) Expone en su tesis doctoral, que como la mayora

de las seales no estn limitadas en banda, antes de los pre-filtros para la
digitalizacin se utilizan frecuencias por encima de la frecuencia de Nyquist.
Es necesario que la frecuencia de muestreo sea lo suficientemente alta para
que en el proceso de filtrado la seal ya filtrada contenga los componentes
mas importantes de la seal de entrada. Se recomienda adoptar un diez
veces mayor .

La cuantificacin es un proceso de aproximacin por redondeo de un rango

de valores continuos par un conjunto finito de valores enteros .
Cada valor entero esta representado por un nmero finito de bits B. El
conversor A/D esta caracterizado por un rango de escala R, el cual se divide
en partes iguales en los niveles de cuantificacin , (Camacho-Tauta,
2011), el Q de resolucin esta dado por la Ecuacin 20:

Ecuacin 20

51
A continuacin se menciona un ejemplo:

Si el conversor A/D es de 9 bits en la resolucin y el rango de escala

completo es de 20 Voltios, se tiene entonces una resolucin de:

Esta resolucin se debe ajustar de acuerdo con la amplitud que se espera de

la seal analgica, de tal forma la escala de medicin debe ser del mismo
orden de magnitud de la seal analgica, de lo contrarios la seal podra
verse distorsionada. (Camacho-Tauta, 2011)

La Figura 30 presenta la seal obtenida a travs del proceso de digitalizacin

la cual tiende a seguir el contorno de la seal analgica. Un proceso exitoso
de digitalizacin depende de la seleccin apropiada de la frecuencia de
muestreo y la amplitud de entrada.

Figura 30. Digitalizacin de una seal analgica, adaptado de (Orfanidis, 1996), (Camacho-Tauta, 2011)

52
4.3 Seales de tiempo continuo

Una seal de tiempo continuo es aquella que toma cualquier valor para
cualquier instante de tiempo; (Kamen & Heck, 2008) hacen referencia a una
seal como una funcin con valor real o escalar de la variable, que en
este caso es el tiempo (t).

Cuando se habla de valor real, se entiende que para un valor fijo cualquiera
de la variable del tiempo t, el valor de la seal en ese tiempo t, es un nmero
real; en el momento en el cual esta variable toma su valor del conjunto de los
nmeros reales, se dice que t es una variable de tipo continuo y que es
una seal de tiempo continuo o una seal analgica.

4.4 Seales de tiempo discreto

Una variable de tiempo es una variable de tiempo discreto, si slo toma

los valores discretos para algn intervalo de numero enteros de .
podra tomar valores enteros ; de tal forma para
. Una seal de tiempo discreto es aquella que esta en funcin del
tiempo, siendo esta una variable discreta . De forma ms sencilla, se
puede decir que una seal de tiempo discreto solo tiene valores en los
puntos de tiempo discreto , en la cual solo toma valores enteros. (Kamen
& Heck, 2008)

Para explicar de forma mas detallada, se denota una seal de tiempo

discreto como [ ], la cual se encuentra definida en los puntos de tiempo
. Cabe resaltar que la variable entera corresponde a instantes . De

53
igual manera los corchetes son utilizados para hacer referencia a una seal
de tiempo discreto [ ], mientras que los parntesis definen una seal de
tiempo continuo (Kamen & Heck, 2008). En cuanto a la forma de
visualizacin, una seal de tiempo discreto est en trmino de los valores
[ ] versus la variable de tiempo . Los valores graficados de [ ] se
representan con crculos rellenos y lneas verticales que los unen con el eje
del tiempo tal como lo indica la Figura 31, siendo esta la grfica que
corresponde a la seal de tiempo discreto [ ] dada por:

[ ] [ ] [ ] [ ] [ ]

Con [ ] para cualquier otro .

Figura 31. Representacin de una seal de tiempo discreto

Seales de tiempo discreto peridicas

Una seal de tiempo discreto [ ] es peridica si existe un entero positivo ,

tal que:

[ ] [ ] Para todos los enteros

54
De tal forma, [ ] es peridica si y solo si existe un entero positivo , tal que
[ ] se repita a si misma a cada instante , donde es conocido como
periodo. El periodo fundamental es aquel valor mas pequeo de , en el cual
la seal se repite. (Kamen & Heck, 2008)

4.5 Representacin de seales en trminos de sus componentes de

frecuencia

El contenido de frecuencia de una seal es un concepto importante en el

estudio de seales, el contenido de frecuencia se puede generar para una
gran cantidad de seales si se divide la seal en sus componentes de
frecuencia que estn dadas por sinusoides. En el siguiente ejemplo se
considera la seal de tiempo continuo x(t) definida por la suma finita de
sinusoides.

Ecuacin 21

Donde:
Entero positivo
Amplitudes de las funciones sinusoidales
Frecuencias de las sinusoides
Fases de las sinusoides

La seal anteriormente nombrada est formada por las frecuencias

y por las componentes de frecuencia de la seal
. Cabe notar que esta seal se caracteriza por completo por las
frecuencias , las amplitudes , (que especifican los
pesos relativos de las componentes de frecuencia y determinan la forma de
la seal) y las fases .

55
Considerando esta seal con igual a la frecuencia variable, las amplitudes
se pueden graficar contra . Dado que existe solamente un nmero finito
de frecuencias presentes en x(t), la grfica de contra constar de un
nmero finito de puntos graficados con las frecuencias presentes en x(t).
Por lo general se dibujan lneas verticales conectando los valores de con
los puntos , la grfica resultante es un ejemplo de un espectro de lnea y
se conoce como espectro de amplitud de la seal x(t), dicho espectro
muestra las magnitudes relativas de las distintas componentes de frecuencia
que forman la seal.

Adems de contar la seal con un espectro de amplitud, la seal cuenta con

un espectro de fase en donde se obtiene una grfica de la fase en grados
o radianes contra la variable de frecuencia . De igual manera se dibujan
lneas verticales conectando los valores de con los puntos de frecuencia
, de modo que el espectro de fase tambin resulta un espectro de lnea,
Figura 32.

Figura 32. Ejemplos de espectros e amplitud de seal (Kamen & Heck, 2008)

56
4.6 Reduccin de ruido

Apilado: Mejora y Resolucin de la relacin de seal para ruido.

El apilado de seal es una alternativa efectiva que rene seales claras por
encima del nivel de ruido de fondo. Este proceso consiste en medir la seal
diferentes veces y promediar todo el conjunto, el elemento i-th en la seal
principal significa el promedio de todos los elementos i-th en las seales
medidas. La suposicin principal es que el ruido tiene cero significado por lo
tanto el promedio reduce el nivel del ruido en la seal principal e incrementa
la relacin de seal para ruido del componente de correlacin.

La siguiente figura muestra una seal de ruido que es simulada por la adicin
de ruido aleatorio para un periodo sinusoidal en descomposicin. La
fluctuacin del ruido de fondo es la misma que la seal peridica por
consiguiente la relacin de seal para ruido es igual a uno, se puede
observar el efecto de apilado para un nmero creciente de seales.
(Santamaria & Fratta, 2005)

Figura 33. Control de ruido por apilado de seal en el tiempo dominante. La relacion de seal para ruido
aumenta con relacion l numero de seales incrementadas. (Santamaria & Fratta, 2005)

57
4.7 Transformada de Fourier

La funcin se conoce como la integral de Fourier o Transformada de

Fourier para la funcin que no es necesariamente peridica, y est
definida de la siguiente forma:

[ ] Ecuacin 22

Esta integral se usa para calcular cuando se conoce la . En general

el resultado de esta integral es un nmero complejo que tiene una parte Real
ms una parte Imaginaria y puede escribirse as:

; En su forma cartesiana Ecuacin 23

; En su forma polar Ecuacin 24

Donde:

: Espectro de magnitud de la

: Espectro de fase de la

La integral o transformada de Fourier permite cambiar de dominio una

funcin conocida que est en el dominio del tiempo, transformndola en
una en el dominio de la frecuencia angular.

La condicin para que exista es que la integral del valor absoluto de

debe ser finita:

Ecuacin 25

Si esto se cumple entonces existe . (Guzmn, 2008)

58
4.8 Funcin de transferencia y coherencia

La funcin de transferencia entre la entrada y la salida es

computada a travs del promedio de datos de diferentes ensayos mostrados
en la siguiente ecuacin, la barra sobre las funciones significa ese promedio.

Ecuacin 26

La funcin de transferencia es evaluada por la funcin de coherencia

definida por la ecuacin nombrada a continuacin y da una medida de la
dependencia lineal entre dos seales como una funcin de frecuencia.

Las frecuencias con valores de funcin de coherencia concluyen a la unidad

y es una relacin lineal perfecta entre las dos seales adems existe una alta
probabilidad de no tener contaminacin generada por ruido u otras no
linealidades.

| |
Ecuacin 27

4.9 Sistemas

(Kamen & Heck, 2008) Definen un sistema como una coleccin de uno o ms
dispositivos, procesos o algoritmos implementados por computadora, que
operan sobre una seal de entrada x para producir una seal de salida y.
Cuando la entrada y la salida son seales continuas en el tiempo, x(t) y y(t),

59
se dice que el sistema es un sistema de tiempo continuo o sistema
analgico. Por otro lado cuando las seales de entrada y de salida son
discretas, x [ ] y y[ ] se dice que es un sistema de tiempo discreto.

Existen diversos ejemplos de sistemas que procesan seales (procesadores

de seales) que operan sobre una seal de entrada para producir una seal
de salida deseada, un ejemplo comn es el filtro que elimina componentes
sinusoidales cuyas frecuencias se encuentran en un cierto rango o remueven
ruido presente en una seal. Igualmente existen sistemas que consisten en
dispositivos fsicos o procesos, algunos de ellos se nombran a continuacin:

Un circuito elctrico con entradas iguales a voltajes y/o corrientes y con

salidas iguales a voltajes y/o corrientes en diversos puntos del circuito.

Una refinera petrolera, con una entrada igual a la velocidad de flujo de

petrleo y con salida igual a la velocidad de flujo de gasolina.

Un sistema de manufactura, con entradas iguales a la velocidad de flujo de

materias primas, y con salidas iguales a la velocidad de produccin de
productos terminados.

60
5 IMPLEMENTACION DEL ENSAYO DE COLUMNA RESONANTE

5.1 Introduccin

El ensayo de columna resonante como herramienta base de estudio y

posterior anlisis de las propiedades dinmicas del suelo requiere para su
implementacin y ejecucin de una serie de factores caractersticos. Dichos
factores como diseo de un programa de medicin, calibracin del equipo y
procedimientos de anlisis se describirn de forma detallada en este
captulo, de igual forma se nombraran las principales caractersticas en
cuanto a los ensayos de inters: ensayo en el dominio de la frecuencia y el
ensayo de no resonancia.

La seleccin de informacin relevante y su posterior interpretacin en cuanto

a programas de diseo que relacionan la interfaz humana con el sistema de
operacin, contribuyen al desarrollo de este proyecto, un ejemplo claro es el
programa LabView, con el cual se elabor el programa que calcula e indica
los resultados para su continuo anlisis. En este captulo se describen todos
los componentes necesarios para lograr la obtencin de informacin til.

5.2 Descripcin programa en LabVIEW

Para el desarrollo del proyecto especficamente el registro de mediciones se

hizo uso del software para diseo de sistemas denominado LabVIEW,
mediante esta herramienta de trabajo se crearon dos programas: uno
enfocado hacia la parte dinmica y otro dirigido hacia la parte esttica.

61
Cada uno de los programas est compuesto por dos pneles, uno frontal y
un segundo conocido como diagrama de bloques, a continuacin se realizar
una breve descripcin de los pneles implementados en diseo de un
software:

Frontal: Es la interfaz con el usuario, es el espacio donde se definen los

controles e indicadores para llegar a la obtencin de resultados, la Figura 34
presenta una descripcin visual de lo que corresponde al panel frontal de un
programa previamente diseado; este panel tiene con funcin principal
responder a las ordenes que el diseador a introducido con anterioridad, por
lo cual requiere de un proceso que garantice la calidad de los resultados que
puede y debe entregar.

Figura 34. Panel frontal (National Instruments Corporation, 2012)

Diagrama de bloques: Es el cdigo que controla el programa (Figura 35), se

compone de una serie de conos interconectados para llevar a cabo las
funciones que el usuario desea; este panel cuenta con herramientas graficas
lo cual permite hacer mas fcil el proceso de programacin de igual manera
permite organizar los comandos y operaciones indicadas con el fin de facilitar
la comprensin del usuario.

62
LabView logra convertirse en un lenguaje de programacin un tanto ms fcil
de entender e implementar a comparacin de C++ o Visual Basic.

Figura 35. Diagrama de bloques (National Instruments Corporation, 2012)

5.3 Calibracin

La calibracin del equipo tiene un papel de gran importancia debido a que las
propiedades dinmicas calculadas a partir de las ecuaciones son
dependientes de los factores de calibracin; como consecuencia a esto se
tiene la calidad y confiabilidad de los resultados generados a partir de los
ensayos.

Este captulo se divide en dos secciones las cuales pretenden explicar los
factores principales de calibracin; en primer lugar se desarrollar lo
relacionado a la calibracin del equipo de columna resonante, en especial al
sistema fijo-libre mejorado y finalmente se pretende desarrollar lo referente a
la calibracin de los sensores para la adquisicin de datos.

5.3.1 Calibracin del equipo de columna resonante

Para la calibracin del equipo es necesario tener presente la inercia

rotacional de la masa activa (JA), la frecuencia de resonancia del aparato

63
(A) y posteriormente el factor torque-corriente (TFC); a continuacin se
realiza una breve descripcin de cada uno de los factores a tener en cuenta.

5.3.1.1 Inercia rotacional de la masa superior activa (JA)

Este procedimiento desarrollado por el fabricante est basado en lo

estipulado por (Tatsuoka & Silver, 1980); en donde se hace uso de una barra
metlica de calibracin y una masa de material idntico con el fin de obtener
la inercia rotacional de la masa activa, (Camacho-Tauta, 2011); tal como se
explica a continuacin:

Se debe determinar de forma geomtrica la inercia rotacional de la barra de

calibracin y de la masa adicional, (Figura 36) y(Figura 37); este
procedimiento se llev a cabo de forma manual y mecnica, para el segundo
mtodo de clculo empleado se utiliz el SoftwareSolid Works, en donde
finalmente fue el valor tenido en cuenta a causa de la mayor precisin que
presentaba su resultado; los valores corresponden a J1= 29,27 Kg cm2 y J2=
185.93 Kg cm2respectivamente.

Figura 36. Barra de Calibracin (J1) Figura 37. Masa adicional (J2)

64
A continuacin se fija la barra de calibracin al pedestal y al cabezal usando
tornillos; de esta forma se determina la frecuencia de resonancia (fn 1) y la
frecuencia de resonancia (fn2) para la barra de calibracin ms la masa
adicional; proceso que se lleva a cabo con diferentes amplitudes (u0), tal
como lo indican laTabla 2 y la Tabla 3.

Tabla 2. Frecuencia de resonancia del sistema + Barra fn1

FRECUENCIA DE
AMPLITUD (mV)
RESONANCIA (Hz)
20 164,62
50 163,26
100 162,18
200 162,09
500 161,22
1000 160,36
PROMEDIO 162,29

Tabla 3. Frecuencia de resonancia del sistema + Barra + Masa (fn2)

FRECUENCIA DE
AMPLITUD (mV)
RESONANCIA (Hz)
20 94,94
50 94,68
100 94,44
200 94,03
500 93,34
1000 92,62
PROMEDIO 94,01

La Ecuacin 28 permite calcular la inercia rotacional de la masa activa

superior, los valores se encuentran resumidos en la Tabla 4 para las
diferentes amplitudes adoptadas; finalmente se tom el valor promedio de
49.85 kg cm2, a continuacin se presenta la correspondiente solucin:

65
 Ecuacin 28
[ ]

[ ]

Tabla 4. Datos de calibracin para la obtencin de la inercia rotacional de la masa activa

AMPLITUD Fn1 Fn2 JA

(mV) Hz Hz Kg cm2
20 164,62 94,94 48,8
50 163,26 94,68 50,1
100 162,18 94,44 51,1
200 162,09 94,03 50,2
500 161,22 93,34 49,7
1000 160,36 92,62 49,2
PROMEDIO 162,29 94,01 49,85

5.3.1.2 Frecuencia de resonancia de aparato (A)

A partir del procedimiento descrito en la norma ASTM D-4015, sin masas o

especmenes dentro de la cmara triaxial, se debe medir la frecuencia de
resonancia del sistema para diferentes niveles de rotacin (Camacho-Tauta,
2011). El promedio obtenido fue de 6.38 Hz (40.09 Rad/s) para un nivel de
vibracin. Tabla 5

66
 Tabla 5. Frecuencia de resonancia del sistema con el cabezal libre

FRECUENCIA DE
AMPLITUD (mV)
RESONANCIA (Hz)
20 6,1
50 6,5
100 6,4
200 6,6
500 6,7
1000 6,0
PROMEDIO 6,38

Ecuacin 29

Para el clculo de la frecuencia de resonancia del aparato se utiliz la

Ecuacin 29 de la cual se obtuvo una a= 40.1 Rad/s.

5.3.1.3 Factor torque-corriente (TCF)

Para la calibracin del factor torque corriente en las bobinas se hace uso
del sistema del servomotor, este elemento tiene la capacidad de entregar un
torque proporcional a una corriente la cual corresponder a la corriente
aplicada a las bobinas. La Tabla 6 y la Figura 38 indican el proceso de
calibracin tenido en cuenta; el factor de correlacin por ser cercano a la
unidad nos indica un excelente ajuste.

Tabla 6. Calibracin Torque-Corriente

VOLTAJE (V) TORQUE (N-m)

0,5 0,22
1 0,43
1,5 0,65
2 0,87
2,5 1,09

67
 1,200

1,00

,800

Torque (N_m) ,600 Torque = 0,4349(Voltaje)

R = 0,9999
,400

,200

,00
,00 ,500 1,00 1,500 2,00 2,500 3,00
Voltaje (V)

Figura 38. Calibracin Torque-Corriente

5.3.1.4 Relacin de amortiguamiento

La relacin de amortiguamiento del aparato se obtiene a travs de las

vibraciones libres de la masa activa en el momento en que la excitacin se
interrumpe. HalfPower fue el mtodo utilizado para la verificacin de los
resultados.

Decaimiento de vibraciones libres

Este mtodo corresponde al procedimiento adoptado de la Norma ASTM-

4050 empleada para obtener las caractersticas de amortiguamiento del
aparato. Mientras el sistema vibra a frecuencia de resonancia se debe
interrumpir la fuente de excitacin e ir registrando las oscilaciones hasta el
reposo, una muestra de esto es la Figura 39a). Posterior a esto se identifican
los picos para terminar por calcular el decremento logartmico () por medio
de la Ecuacin 30. (Camacho-Tauta, 2011)

68
 0
1.0 10

0.5

u0 [mm]
u [mm]

u0=1.1119e-0.285Cycle
-0.5

-1
-1.0 10
0 0.5 1 1.5 0 2 4 6 8
t [s] Cycle

a) b)
Figura 39 Determinacin de la relacin de amortiguamiento del aparato. a) Curva de decaimiento de
vibraciones libres. b) Amplitud - nmero de ciclos (Camacho-Tauta, 2011)

( ) ( ) Ecuacin 30

Dnde:

u1 : Amplitud de vibracin del primer ciclo despus de que la alimentacin

se interrumpe

un+1 : Amplitud para el ciclo (n+1)

Un resultado similar se puede obtener calculando la pendiente de la amplitud

logartmica como funcin del numero de ciclos para los primeros siete ciclos
Figura 39. La relacin de amortiguamiento del aparto se calcula utilizando la
Ecuacin 31. (Camacho-Tauta, 2011)

Ecuacin 31

69
Mtodo HalfPower

El mtodo HalfPower, es un mtodo de la teora de la dinmica de

estructuras el cual esta fundamentado en la medicin de las amplitudes de
desplazamiento en vibracin forzada sobre un rango de frecuencias bajo
amplitud de torque constante (Clough & Penzien, 1995). Este mtodo toma
dos frecuencias f1 y f2 en donde la amplitud es igual que la amplitud en
resonancia sobre 2. La Ecuacin 32 perite calcular la relacin de
amortiguamiento del sistema.

Ecuacin 32

La Figura 40a) presenta dos curvas, la lnea azul indica que la corriente de
salida no es constante (voltaje constante), esto debido a que la variable no
es controlada por el usuario. La lnea negra es la amplitud de la respuesta en
el caso en que la corriente en las bobinas se mantenga en la misma
intensidad que la aplicada durante la resonancia. Posteriormente se pueden
identificar las frecuencias f1 y f2 la relacin de amortiguamiento haciendo uso
de la Ecuacin 32.

0.25 1.00 Constat Voltage

Constat Voltage Constant Torque
0.20 Constant Torque f2
0.707 f1
0.15
u0 [mm]

u0/u0 max

0.50
0.10

0.05 0.25

0 0
5 10 15 20
5 10 15 20
f = /2 [Hz]
f = /2 [Hz]

Figura 40. Determinacin de la relacin de amortiguamiento del aparato mediante el mtodo half-power.
Amplitud de desplazamiento y corriente para frecuencias cercanas a la resonancia. a) Frecuencia; b)
frecuencia normalizada. (Camacho-Tauta, 2011)

70
5.3.2 Calibracin de los sensores para la adquisicin de datos

5.3.2.1 Acelermetro

La Figura 41 indica una funcin sinusoidal de aceleracin en funcin del

tiempo, el propsito es poder visualizar la conversin de lecturas de
aceleracin a velocidad y desplazamiento por medio de sus derivadas para
finalmente obtener la rotacin en trminos de voltaje.

(mm/s2)
1,5
1
0,5
0
0 5 10 15 20
-0,5
-1
-1,5 t (s)
Figura 41. Funcin de aceleracin

Ecuacin 33

Donde:
V = Voltaje
G= Ganancia
S= Sensibilidad

Cabe aclarar la ganancia es el valor en el que puede ser graduado el equipo

(1, 10 100) para que el instrumento sea lo suficientemente sensible.
Adicional a esto y teniendo en cuenta que la velocidad es el resultado de la
derivada de la aceleracin, se puede presentar Figura 42, en la cual se

71
indican el desplazamiento resultante de la Ecuacin 34, relacionndola con la
frecuencia angular.

(mm/s)
1,5
1
0,5
0
-0,5 0 5 10 15 20
-1
-1,5
t (s)

Figura 42. Funcin de velocidad

Ecuacin 34

Con la derivad de la velocidad se obtiene la funcin de desplazamiento la

cual se puede observar en la Figura 43; la siguiente ecuacin corresponde a
la funcin mencionada:

Ecuacin 35

u (m)
1,5
1
0,5
0
-0,5 0 5 10 15 20

-1
-1,5
t (s)

Figura 43. Funcin de desplazamiento

72
Remplazando y en la Ecuacin 35 se puede obtener:

( ) Ecuacin 36

Finalmente, para determinar la rotacin en funcin del voltaje y la frecuencia

angular se puede apreciar la siguiente imagen:

Figura 44. Relacin desplazamiento rotacin

De esta forma se puede obtener la rotacin en funcin del voltaje que es lo

que finalmente se est buscando, as pues, se tiene:

( ) Ecuacin 37

Para realizar la calibracin del sensor de aceleracin se tiene en cuenta la

siguiente relacin:

Dada la sensibilidad del sensor:

73
Debido a que la tarjeta de adquisicin de datos lee valores en unidad de
voltaje se hace necesario sacar el inverso de la sensibilidad con el propsito
de determinar a cuantos corresponde un , de esta forma se obtiene:

La tarjeta de adquisicin cuenta con un rango de lectura comprendido entre -

a+ , valores que expresados en voltaje corresponden a:

5.3.2.2 Celda de carga

Para llevar a cabo el proceso de calibracin de este sensor se hace uso de

un manmetro el cual indica la presin a la cual est sometida la muestra,
esta presin est representada por un voltaje; la Tabla 7 y la Figura 45
indican lo valores de calibracin y su respectiva representacin grfica.

Tabla 7. Calibracin Celda de Carga

ANILLO FUERZA CELDA PRESIN

DIVISIONES N mV kPa
205 1863,372 -0,59188 130
402 3654,027 -1,03964 250
604 5490,130 -1,47237 330

74
 6000
5000
4000

Fuerza (N)
3000
Fuerza = -4118.4 (Voltaje) - 591.8
2000
R = 0.9997
1000
0
-2,00 -1,500 -1,00 -,500 ,00
Voltaje (V)

Figura 45. Calibracin Celda de Carga

5.3.2.3 Presin de cmara

Es necesario tener presente que antes de ejecutar la calibracin de la

presin de cmara es absolutamente necesario verificar que la campana con
la cual cuenta el equipo est en la posicin indicada con el fin de evitar algn
tipo de fuga; una vez comprobada la instalacin de este elemento se procede
a aplicar una presin controlada por un manmetro dentro de la cmara; esta
presin es leda por el sensor de presin entregando como resultado un valor
voltaje proporcional. La Tabla 8 y la Figura 46 indican los registros de la
calibracin realizada.

Tabla 8. Calibracin Presin de Cmara

VOLTAJE PRESIN DE CMARA

(mV) (kPa)
2,37 0,00
2,50 50,00
2,79 100,00
3,15 150,00
3,40 200,00
3,78 250,00
4,08 300,00
4,74 400,00

75
 450
400
Presin de cmara = 162.71 (Volyaje) - 363.93
350

Presin de cmara (kPa)

R = 0.9939
300
250
200
150
100
50
0
2,00 2,500 3,00 3,500 4,00 4,500 5,00
Voltaje (mV)

Figura 46. Calibracin Presin de Cmara

5.3.2.4 Contrapresin

La Tabla 9 y la Figura 47 indican los resultados obtenidos en el proceso de

calibracin; para llevar a cabo este procedimiento se utiliza una bomba de
presin la cual permite variar los valores de presin aplicados, dichos valores
son ledos por el sensor y posteriormente entregados en voltaje; cabe aclarar
que los procedimiento de calibracin tanto del sensor de presin de poros,
presin de cmara y contrapresin son similares.

Tabla 9. Calibracin Contrapresin

VOLTAJE PRESIN (kPa)

-0,13 0
0,35 100,00
0,93 200,00
1,39 300,00
1,90 400,00
2,40 500,00

76
 600

500

Contrapresin (kPa)
400

300

200

100
Contrapresin = 196,94 (Voltaje) + 25,46
0
R = 0,9995
-100
,00 ,500 1,00 1,500 2,00 2,500
Voltaje (V)

Figura 47. Calibracin Contrapresin

5.3.2.5 Presin de poros

Como ya se ha mencionado en el tem anterior, la calibracin del sensor se

efecta aplicando una presin controlada que posteriormente es leda por el
sensor y finalmente convertida a voltaje; la Tabla 10 y la Figura 48 hacen
referencia a la calibracin del sensor de presin de poros.

Tabla 10. Calibracin Presin de Poros

VOLTAJE PRESIN (kPa)

-0,07 0
0,21 100,00
0,42 200,00
0,68 300,00
0,92 400,00
1,18 500,00

77
 600

500

Presin de poros (kPa)

400

300

200

100 Presin de poros = 404,48 (Voltaje) + 24,606

R = 0,9995
0

-100
,00 ,200 ,400 ,600 ,800 1,00 1,200 1,400
Voltaje (V)

Figura 48. Calibracin Presin de Poros

5.3.2.6 Deformacin vertical

Es necesario utilizar un calibrador de potencimetro el cual permite generar

una deformacin contralada al sensor; este valor es ledo por el sensor en
voltaje, a continuacin se presentan los resultados de la calibracin, Tabla 11
y Figura 49.

Tabla 11. Calibracin Deformacin Vertical

DESPLAZAMIENTO
VOLTAJE (V)
VERTICAL (mm)
7,00 0,00
5,83 2,00
3,19 4,00
1,14 6,00
3,25 8,00

78
 7

Deformacion vertical (mm)

5 Deformacin vertical= -0,9709 (Voltaje) + 7,166
R = 0,9803
4

0
,00 2,00 4,00 6,00 8,00
Voltaje (V)

Figura 49. Calibracin Deformacin Vertical

5.3.2.7 Cambio volumtrico

La calibracin se desarrollada haciendo uso de una bureta con el fin de

controlar la cantidad de agua aplicada al cilindro de cambio volumtrico, al
introducir progresivamente el lquido se genera un desplazamiento el cual es
medido por el sensor, finalmente se realiza la relacin de voltaje y volumen
tal como se indica en la Tabla 12 y Figura 50.

Tabla 12. Calibracin Cambio Volumtrico

CAMBIO VOLUMTRICO
VOLTAJE (V)
(ml)

10,06 0,00
8,43 18,00
6,68 39,00
5,00 58,00
2,82 80,50
1,12 100,38

79
 120

100

Cambio Volumtrico (ml)

80 Cambio volumtrico = -11,199 (Voltaje) + 112,98
R = 0,9996
60

40

20

0
,00 2,00 4,00 6,00 8,00 10,00
Voltaje (V)

Figura 50. Calibracin Cambio Volumtrico

5.4 Ejecucin del ensayo de columna resonante

5.4.1 Datos de calibracin

Los datos de calibracin, indicados en la Figura 51 usando dgitos que

permiten distinguirlos como pasos, son datos constantes y son los primeros
que el usuario ingresa, estos datos son importantes porque en cierta medida
de ellos depende el grado de confiabilidad y la calidad de los ensayos. En la
Figura 51 (paso 1) se indican grficamente los factores de calibracin que
previamente ya fueron explicados:

- Inercia rotacional de la masa activa (JA)

- Frecuencia resonante torsional (f A)
- Localizacin del transductor de movimiento (B)
- Relacin de amortiguamiento del aparato ( A)

80
Una vez ingresados los datos el usuario define el nmero de intervalos,
tiempo de ejecucin y nmero de pruebas deseadas (paso 3 y 4).

Figura 51. Factores de calibracin

5.4.2 Datos iniciales del espcimen

Son los datos caractersticos de las muestras que van a ser evaluadas, (paso
1), se cuenta con una carpeta donde se guardar la informacin, el usuario
asignar un nombre para el ensayo, se registrara la fecha y hora de
ejecucin y una pequea descripcin, (paso 2). De igual forma el usuario
seleccionar el tipo de consolidacin que necesite. A continuacin se har
una breve descripcin de cada uno de ellos, la Figura 52 presenta las
pestaas, que sern explicadas, (paso 3):

81
5.4.2.1 Medidas iniciales

Dimetro externo: Es el ancho de cualquier cara de la muestra cilndrica.

Dimetro interno: En este tipo de ensayos se trabaja con muestras donde

este dimetro siempre va a ser cero debido a que no presentan ningn tipo
de orificio.

Longitud: Magnitud que expresa la distancia entre dos puntos o cada una de
las dimensiones de un cuerpo (Envl & KG, 2007), para este caso es la
distancia entre los extremos inferior y superior de las muestras.

Masa total: Cantidad de materia que contiene la muestra.

Contenido de humedad: relacin existente entre el peso del agua y el peso

de slidos.

Ecuacin 38

Dnde:
: Contenido de humedad
: Peso del agua
: Peso de slidos

Gravedad especifica de slidos: Es el peso unitario del material de la

muestra (arena) dividido por el peso unitario del agua.

Ecuacin 39

82
Dnde:
: Contenido de humedad
: Masa de slidos
: Volumen de slidos
: Densidad del agua

Gravedad especifica del agua: unidad basada en la densidad, el agua tiene

un valor de 1 g/cm3.

5.4.2.2 Lecturas iniciales

Desplazamiento axial: Es la distancia a la cual la muestra llega a

expandirse, deformacin que se presenta cuando se aplica una carga.

Ecuacin 40

Dnde:

: Deformacin axial

: Deformacin vertical

: Longitud inicial

Cambio volumtrico: Es aquel cambio que sufre la muestra con relacin al

volumen antes y despus de colocar una carga, (paso 4).

Ecuacin 41

83
Dnde:

: Cambio volumtrico

: Longitud inicial

: Deformacin volumtrica

De igual manera se cuenta con botones como save para guardar la

informacin, about para presentar la informacin relacionada con los autores
y el software y finalmente exit para salir, (paso 5).

Figura 52. Datos iniciales del espcimen.

84
5.4.3 Datos de estado

Este tipo de datos son variables y se van actualizando a medida que el

ensayo se est ejecutando, hacen referencia a las presiones utilizadas en
cada ensayo, a aquellas presiones a las que est sometida la muestra y las
deformaciones resultantes, as como los valores actuantes del especimen,
propiedades ndice y estados de esfuerzo. A continuacin se explicar cada
uno de estos datos, la Figura 53 indica grficamente la ubicacin de la
pestaa (paso 1):

Figura 53. Datos de estado.

85
5.4.3.1 Presiones

Presin de poros: Considerada como una presin hidrosttica. (Fine-Civil

Engineering Software, 2007). Es aquella presin que acta sobre los fluidos
en los espacios porosos de un suelo. (Pemex, 2009). (paso 2).

Ecuacin 42

Dnde:
: Presin de poros
Esfuerzo
: Esfuerzo efectivo

Presin lateral: Presin que se ejerce alrededor de la muestra.

Contrapresin: Es la presin de agua aplicada a la muestra durante el

proceso de saturacin y consolidacin.

5.4.3.2 Esfuerzos

Esfuerzo radial efectivo: Diferencia entre el esfuerzo total en dicha

direccin y la presin del agua que existe en los vacos del suelo (Alva,
2010).

Esfuerzo axial efectivo: Diferencia entre el esfuerzo total en torno al eje

vertical y la presin del agua que existe en los vacos del suelo.

Esfuerzo isotrpico efectivo: Esfuerzo en el cual las presiones que se

efectan deben ser iguales.

86
 Ecuacin 43

: Esfuerzo isotrpico
: Esfuerzos normales efectivos en las tres direcciones

Esfuerzo desviador: Es el valor de la diferencia entre los esfuerzos

principales. (Instituto Nacional de Vas), (paso 3).

Ecuacin 44

Dnde:
: Esfuerzo desviador
: Carga
: rea

5.4.4 Otras propiedades

Cambio volumtrico: Es aquel cambio que sufre la muestra con relacin al

volumen antes y despus de colocar una carga, (paso 4).

Carga vertical: Fuerza que se aplica en la longitud de la muestra.

Longitud: Magnitud que expresa la distancia entre dos puntos o cada una de
las dimensiones de un cuerpo (Envl & KG, 2007) , para este caso es la
distancia entre los extremos inferior y superior de las muestras.

87
rea: Superficie de la muestra comprendida dentro de un permetro
expresada en unidades de medidas conocidas como superficiales.
(Definicion.De, 2008)

( ) Ecuacin 45

Volumen: Magnitud fsica que mide la cantidad de espacio ocupado por la

muestra.

Ecuacin 46

Densidad de masa del suelo: Relacin existente entre la masa y el

volumen.

Ecuacin 47

Relacin de vacos: Relacin entre el volumen de vacos y el volumen de

slidos de una muestra determinada.

( ) Ecuacin 48

Dnde:
: Gravedad especifica
: Densidad total de masa
: Contenido de humedad

Exceso de presin de poros: Aumento de la presin del agua que hay en

los vacos del suelo respecto a un valor de referencia cuando es aplicado un
valor externo.

88
Finalmente se cuenta con una serie botones para guardar la informacin,
hacer que los sensores vuelvan a estar en cero, y uno para salir de esa
seccin.

5.4.5 Resultados obtenidos

5.4.5.1 Ensayo en el dominio de la frecuencia

Inicialmente se debe enviar una seal de tipo sinusoidal haciendo uso de un

generador de funciones (HAMEG), esta seal es enviada en forma de voltaje
y transmitida hacia un amplificador (BRIDGABLE), para finalmente ser
enviadas a las bobinas, en el panel de navegacin se puede seleccionar el
tipo de seal de entrada que el usuario desee analizar; para este
procedimiento debe cambiarse la opcin de la pestaa Test Type.

En el paso 1, los factores como amplitud, frecuencia de inicio y frecuencia

final pueden ser controlados y/o modificados por parte del operador o ser el
generador de funciones el que vare los datos. A continuacin se describen
cada uno de estos elementos:

Tipo de ensayo: El tipo de seal que ser enviada a las bobinas para
efectuar el ensayo, por defecto est predeterminado un ensayo tipo
SineSweep.

Amplitud: Este parmetro el sistema lo selecciona para llegar a una

deformacin de corte deseada.

Frecuencia inicial: Nmero de repeticiones por unidad de tiempo que

ingresa el usuario.

89
Frecuencia final: Nmero de repeticiones que el usuario desea para que
termine los ensayos.

Nmero de ensayos: Nmero de pruebas que el usuario quiere que se

realicen por cada ensayo para obtener un promedio que ser analizado al
final.

Ventaneo: Proceso en el cual se eliminan las discontinuidades generadas

entre el final de una seal y el comienzo de una nueva, para lograr esto es
necesario aplicar una frecuencia ms alta y multiplicar los datos por una
funcin en este caso rectangular para lograr la continuidad que se requiere y
reducir a cero el perodo de las seales.

Frecuencia de muestreo: Nmero de datos que recibe la tarjeta de

adquisicin de datos por segundo.

Los datos de las casillas del paso 2, en el cuales se calculan las propiedades
dinmicas de la muestra se irn actualizando segn la cantidad de ensayos
que el usuario previamente ha establecido, los resultados que se obtienen
son:

df: Distancia entre cada uno de los puntos que componen una seal en
trminos de frecuencia. Se cuenta con la opcin de seleccionarla o de
digitarla.

Tiempo por ensayo: Duracin de cada ensayo.

Muestras por ensayo: Nmero de ejemplares que se analizan por cada

ensayo.

90
Step: Nmero secuencial que marca cada cuanto se est obteniendo un
valor de frecuencia.

Frecuencia de resonancia: Es el valor ms alto de frecuencia obtenido en el

ensayo y se halla con la funcin de transferencia. Este valor se puede
visualizar en la casilla correspondiente o se puede obtener de manera
manual con ayuda del botn inferior llamado ajuste automtico en el cual el
usuario con ayuda del posicionamiento de dos cursores encuentra el valor
adecuado de resonancia.

Deformacin de corte: Deslizamiento relativo o rodadura entre partculas.

Mdulo de corte: Relacin entre la tensin de corte () y la deformacin

unitaria de corte ().

Velocidad de la onda s: Velocidad de propagacin de las partculas en el

medio; se desplazan perpendicularmente a la direccin de propagacin por lo
que estn asociadas con deformaciones de corte en el del terreno.

Relacin de amortiguamiento: Se obtiene a travs de las vibraciones libres

de la masa activa en el momento en que la excitacin se interrumpe.
HalfPower fue el mtodo utilizado para la verificacin de los resultados.

Las pestaas 3 a 7 de la Figura 54 presentan los resultados grficos

generados por la ejecucin del ensayo:

91
 Figura 54. Panel de navegacin. Ensayo en el dominio de la frecuencia.

3. Seales de tiempo: Corresponde a las seales en tiempo real del torque

y de aceleracin que se generan en la muestra. Para la obtencin de la seal
del torque se utilizan bobinas y para la aceleracin se hace uso de un
acelermetro.

4. Espectro: Representa la frecuencia e intensidad contenida en una seal.

Se puede observar en la figura que para cada frecuencia existe un valor de
torque y de aceleracin.

5. Espectro cruzado y coherencia: En esta ventana se observa la

evaluacin de correlacin que se genera entre la entrada y la salida de una
seal (coherencia), dicha correlacin debe ser constante y vara
generalmente entre 0 y 1.

92
6. Relacin entrada salida y ngulo de fase: Es aquella relacin entre la
seal de entrada (torque) y la seal de salida (aceleracin), esta relacin se
caracteriza por tener un ngulo de fase que se puede observar en la figura.

7. Funcin de transferencia: Permite el procesamiento y comparacin entre

la rotacin y el torque generado en el sistema para encontrar una frecuencia
de resonancia.

El usuario debe seleccionar la carpeta en la cual desea almacenar los

resultados obtenidos (8). Finalmente se da inicio al ensayo (9).

A continuacin se presenta la tabla que indica los parmetros de inicio para

la ejecucin del ensayo en el dominio de la frecuencia, (Camacho-Tauta,
2011):

Tabla 13. Parmetros de muestreo para el ensayo en el dominio de la frecuencia. (Camacho-Tauta, 2011)

VALOR VALOR VALOR

PARAMETRO SIMBOLO UNIDAD
MINIMO MAXIMO RECOMENDADO
Amplitud A V 0,005 2 (a)
Frecuencia
fstart Hz 5 fstop fr - 25
Inicial
Frecuencia
fstop Hz fstart 1000 fr + 25
Final
Numero de
Nt - 1 10
Ensayos
Ventana w - - - Rectangular
Frecuencia de
fs Hz >4fstop 50M 24k
muestreo
Rango del
ChA V 0,050 20 (b)
canal A
Rango del
ChB V 0,050 20 (b)
canal B
Delta de
df Hz fs/65536 fs/128 0,37
Frecuencia
(a) Depende del mdulo de corte deseado
(b) Debe ajustarse de acuerdo al nivel de entrada durante el ensayo

93
5.4.5.2 Ensayo de no resonancia

Figura 55. Panel de navegacin. Ensayo de No Resonancia.

El paso 1 indicado en la Figura 55 permite ingresar valores tales como lo que

a continuacin sern mencionados:

Deformacin de corte: Deslizamiento relativo o rodadura entre partculas.

Error: Es el valor que indica un grado de confiabilidad y exactitud del

ensayo, en este caso hace referencia a ese porcentaje que el usuario define
para obtener un valor apropiado de deformacin de corte. Se representa
mediante la siguiente ecuacin:

Ecuacin 49

94
Amplitud: Definida como la variacin del desplazamiento en forma peridica,
para este tipo de ensayos se observa una amplitud de pico a pico.

Frecuencia inicial: Nmero de repeticiones por unidad de tiempo que

ingresa el usuario.

Frecuencia final: Nmero de repeticiones que el usuario desea para que

termine los ensayos.

Frecuencia de paso: Es una iteracin que comienza con la frecuencia de

inicio hasta llegar a la frecuencia final.

Ciclos: Repeticin de la parte ms pequea de una onda.

Corte Inicial: Fragmento de la seal que se debe omitir para evaluar

solamente la fraccin definida y obtener datos ptimos.

Frecuencia de muestreo: Nmero de datos que recibe la tarjeta de

adquisicin de datos por segundo.

Nivel de vibracin: Es un indicador de las frecuencias tanto real como

deseada.

3. Diagrama de tiempo: Esta ventana se encarga de mostrar la funcin

generada por el torque generado en un transcurso de tiempo determinado,
de igual forma se puede observar la funcin generada por la relacin entre la
vibracin y el tiempo.

95
4. Figura de Lissajous: Es una representacin de los valores de la seal de
entrada (torque) comparada con los valores de la seal de salida
(aceleracin), por lo general tiene forma de elipse.

5. Amplitud: Parmetro necesario para hallar una deformacin de corte

deseada.

6. Relacin entrada salida y ngulo de fase: Es aquella relacin entre la

seal de entrada (torque) y la seal de salida (aceleracin), esta relacin se
caracteriza por tener un ngulo de fase que se puede observar en la figura.

7. Funcin de transferencia: Permite el procesamiento y comparacin entre

la rotacin y el torque generado en el sistema para encontrar una frecuencia
de resonancia.

8. Resultados: Finalmente despus de haber realizado una serie de clculos

y adquirir diferentes funciones, el usuario cuenta con una ventana donde se
muestran los datos fundamentales y de inters para su posterior anlisis e
interpretacin, el mdulo de corte y la relacin de amortiguamiento
relacionados con un intervalo de frecuencia.

La Tabla 14 hace referencia al rango promedio de los valores que pueden

tomar los parmetros que son motivo de estudio en el ensayo de no
resonancia.

96
 Tabla 14. Parmetros de muestreo para el ensayo de no resonancia. (Camacho-Tauta, 2011)

VALOR VALOR VALOR

PARAMETRO SIMBOLO UNIDAD
MINIMO MAXIMO RECOMENDADO
mdulo de
- 10-6 10-5 (a)
corte
Error e Hz 0 10 2
Amplitud A V 0,005 2 Automatica
Frecuencia
fstart Hz 5 fstop - fstep fr - 25
Inicial
Frecuencia
fstop Hz fstart + fstep 1000 fr + 25
Final
Frecuencia de
fstep Hz 0,1 fstop - fstart 0,5
Paso
Ciclos Nc - 1 50 100
Corte inicial ic % 0 100 30
Frecuencia de
fs Hz >4fstop 50M 24k
muestreo
Rango del
ChA V 0,050 20 (b)
canal A
Rango del
ChB V 0,050 20 (b)
canal B
(a) Depende del mdulo de corte deseado
(b) Debe ajustarse de acuerdo al nivel de entrada durante el ensayo

97
6 CONCLUSIONES

Se adquirieron resultados notables y satisfactorios en la realizacin del

proyecto, actividades importantes como: clasificacin e interpretacin de
informacin preliminar, identificacin de las principales propiedades
dinmicas del suelo, diseo del programa de adquisicin de datos,
calibracin del equipo y elaboracin del manual de funciones, son
algunos de los objetivos que se alcanzaron exitosamente.

El software implementado permite generar una variedad de seales en

un rango bastante amplio de frecuencias, controla y actualiza datos
usando como herramienta los sensores de presin y elementos que
conforman el equipo tales como el actuador neumtico y la celda de
carga conformando de esta forma un equipo de columna resonante de
gran eficiencia para la determinacin de las propiedades dinmicas del
suelo.

La implementacin del lenguaje de programacin LabView logr

simplificar y facilitar en gran medida la elaboracin del software que
permite poner en marcha el equipo; en primera instancia por ser un
programa bastante grfico logra hacer el desarrollo menos complejo; otro
punto a favor es que la implementacin del software requiere de una gran
cantidad de operaciones internas que finalmente pudieron ser
simplificadas con el uso de subrutinas que permiten integrar un grupo de
procedimientos previamente definidos.

98
Los parmetros dinmicos del suelo, tales como: el mdulo de corte, la
relacin de amortiguamiento y la deformacin de corte representan una
parte fundamental del comportamiento suelo, estas propiedades
determinadas por el ensayo de columna resonante forman parte de un
gran conjunto de elementos que al relacionarse y convertirse en
funciones llegan a generar informacin valiosa en el campo de la
ingeniera ssmica.

99
7 RECOMENDACIONES

Aquellos estudiantes que continen con el desarrollo del proyecto se les

recomienda llevar a cabo la calibracin de todos los sensores instalados
en el equipo (celda de carga, presin de poros, presin de cmara,
acelermetro, torque, cambio volumtrico, contrapresin y
desplazamiento vertical) debido a la sensibilidad que presenta cada uno
de los elementos mencionados; se recomienda igualmente que esta
calibracin se realice de forma peridica.

Para lograr una ptima realizacin de ensayos es importante tener un

manejo adecuado de los elementos que intervienen en la ejecucin de la
prueba; por consiguiente se propone estudiar de manera anticipada el
manual del usuario en el cual se describe detalladamente cada uno de
los comandos y resultados que entrega el programa.

Se recomienda como futuro proyecto de grado el estudio de materiales

de diferentes granulometras con el propsito de analizar sus resultados
y comportamiento ante condiciones idnticas de estudio.

Para la preparacin de las muestras motivo de estudio se sugiere tener

en cuenta las recomendaciones estipuladas por las normas
correspondientes para tal fin.

Es necesario contar con condiciones ptimas tanto de conexiones de

tierra como fsicas para la realizacin de las muestras, para garantizar
resultados con un margen de error mnimo.

100
Evitar cualquier tipo de daos en los cables al desconectarlos; existen
cables que no pueden ser desconectados como el caso de los Bender
Elements, por lo tanto impedir cualquier tipo de contacto en ellos. El
nico cable que se debe desconectar para el montaje de las muestras es
el del acelermetro.

Los orings, elementos que permiten el hermetismo en la cmara deben

estar siempre en su debido lugar.

El amplificador de seales (generador de funciones) debe cumplir con las

condiciones iniciales implementadas por el tcnico.

Cerciorarse de despejar el lugar de trabajo en el momento de levantar la

campana del equipo. Adems en el momento de su colocacin es
necesario observar que ningn cable sufra algn dao y que todo est en
su lugar.

Hacer uso de los elementos de proteccin en la ejecucin de los

ensayos.

Garantizar que cada vez que se realice un ensayo, los sensores estn en
cero

101
 BIBLIOGRAFA

Fine-Civil Engineering Software. (2007). Retrieved 08 15, 2012, from

http://www.finesoftware.es/quienes-somos/

Definicion.De. (2008). Retrieved 08 15, 2012, from http://definicion.de/area/

Ensayos Ingeniera Civil. (2011, 05 06). Retrieved 06 20, 2012, from

http://www.ensayosingenieriacivil.com/2011/05/el-ensayo-de-
penetracion-estandar-spt.html

National Instruments Corporation. (2012). Retrieved 07 2012, 23, from

http://sine.ni.com/cs/app/doc/p/id/cs-12742

Alarcn, A., Chameau, J., & Leonards, G. (1986). A New Apparatus for
Investigating the Stress-Strain Characteristics of Sands. Geotechnical
Testing Journal.

Alva, J. (2010, Octubre 08). Universidad Nacional Facultad de Ingeniera Civil

Postgrados. Retrieved Noviembre 16, 2011, from
http://es.scribd.com/doc/35684287/Dinamica-de-Suelos

Assimaki & Kausel. (2000). Model for Dynamic Shear Modulus and Damping
for Granular Soils. J. Geotech. Geoenviron.Eng.

Avitabile, P. (2001). Experimental Modal Analysis.

Bala, A., & Balan, S. (2006). Seismic Site Effects Based on In Situ Borehole
Measurements in Bucharest, Romania. Bucarest.

Berry, L. P. (1993). Mecnica de Suelos. Santa Fe de Bogot: McGraw Hill.

Bharat, B. P. (2009). Fundamentals of Soil Dynamic and Earthquake

Engineering.

Camacho - Tauta, J. F. (2008). Rigidez del suelo a pequeas deformaciones.

Evaluacin del efecto de la frecuencia.

Camacho-Tauta, J. F. (2011). Evaluation of the Small-Strain Stiffness of Soil

by Non-Conventional Dynamic Testing Methods. Lisboa.

102
Castrellon, O., Gmez , E., Romo , M., & Flores, R. (n.d.). Instrumentacin y
Automatizacin del Equipo Triaxial. Retrieved from
http://www.sismica.com.mx/ACERVO-
ART/Instrumentacion_equipo_triaxial_c.pdf

Cedex. (2010). Cedex. Retrieved 06 20, 2012, from Tcnicas Geofsicas

Cross-Hole y Down Hole: http://www.cedex.es/lg/geofis/cd_hole.htm

Clough, R., & Penzien, J. (1995). Dynamic of Structures, Computers &

Structures. Berkeley.

Concha, L., Poblete, F., & Mauro, E. (2004, 11). Revistas Electrnicas UACh.
Retrieved 06 22, 2012, from Revistas Electrnicas UACh:
http://www.google.com.mx/imgres?q=ensayo+triaxial+ciclico&um=1&hl
=es&sa=N&qscrl=1&nord=1&rlz=1T4TSNG_esCO425CO426&biw=13
66&bih=614&tbm=isch&tbnid=m-
xTxTFEVj6uUM:&imgrefurl=http://mingaonline.uach.cl/scielo.php%3Fpi
d%3DS0718-025X2004000100005%26script%3Dsc

Diaz, J. (2005). Dinmica de Suelos. Mexico: Limusa.

Dobrin, M., & Savit, C. (1988). Introduction to Geophysical Prospecting.

Singapore: McGraw Hill Book Co.

Drevinch, V. (1967). Effect of Strain History on the Dynamic Properties of

Sand, PhD Dissertation.

Envl, & KG. (2007). Diccionario Manual de la Lengua Espaola. Larousse

Editorial SL.

Ghent. (1986). Influence of Rotation of Principal Stress Directions on the

Cyclic Bahaviour of Sand 2nd International Symposium on Numerical
Models in Geomechanics.

Guzmn, A. (1986). A New Apparatus for Investigating the Stess-Strain

Characteristic of Sands.

Guzmn, M. C. (2008). Matemticas Especiales Para Ingeniera. Medelln:

Textos Acadmicos.

Hardin & Richard. (1963). Elastic Wave Velocities in Granular Soils. J. Soil
Mech. Found.

103
Hurtado, D. I. (2010, Octubre 08). Universidad Nacional, Facultad de
Ingeniera Civil, Postgrados. Retrieved Noviembre 16, 2011, from
http://es.scribd.com/doc/35684287/Dinamica-de-Suelos

Ibrahim, A., Orense, P., Pender, M., Tai, A., & Kikkawa, N. (2011). Laboratory
Determination of Small Strain Shear Modulus of Auckland Residual
Soil.

Instituto Nacional de Vas. (n.d.).

ftp://ftp.unicauca.edu.co/...07/.../Norma%20INV%20E-153-07.pdf.
Retrieved 09 07, 2012, from
ftp://ftp.unicauca.edu.co/...07/.../Norma%20INV%20E-153-07.pdf:
ftp://ftp.unicauca.edu.co/Facultades/FIC/IngCivil/Especificaciones_Nor
mas_INV-07/Normas/Norma%20INV%20E-153-07.pdf

Ishibashi, I., & Zhang, X. (1993). Unified Dynamic Shear Moduli and Damping
Ratios of Sand and Clay. Soils and Foundations.

Ishihara, K. (1996). Soil Behaviour in Earthquake Engineering. Tokyo.

Iwasaki, & Tatsuoka. (1977). Effects of Grain Size and Grading on Dynamic
Shear Moduli of Sands. Soils Found.

Kalinski, M., & Thummaluru, M. (2005). A New Free-Free Resonant Column

Device for Measurement of Gmax and Dmin at Higher Confining
Stresses. Geotechnical Testing Journal.

Kamen, E. W., & Heck, B. S. (2008). Fundamentos de Seales y Sistemas

Usando la WEB y MATLAB (Tercera Edicion ed.). Mexico: PEARSON.

Kausel, E. (2002). Seismic Simulation of Inelatic Soils via Frequency-

Dependent Moduli and Damping. J. Eng. Mech.

Klee, Griem, W., & Griem, S. (1999). Universidad de Atacama, Copiap-

Chile. Retrieved 10 09, 2011, from
http://www.geovirtual2.cl/geologiageneral/ggcap01c.htm

Kokusho. (1981). In Situ Dynamic Properties and Their Evaluation.

Preoceedings of the 10th International Conference on Soil Mechanics
and foundation Engineering.

Kramer. (1996). Geotechnical Eathquake Engineering.

104
Lombillo, I., & Villegas, L. (2010). Metodologas no Destructivas Aplicadas a
la Rehabilitacin Estructural del Patrimonio. Retrieved 06 22, 2012,
from
http://grupos.unican.es/gted/archivos/otros/ponencias/rehabend06_il.p
df

Love, A. E. (1911). Some Problems of Geodynamics. Cambridge: Cambridge

University Press.

Moreno & Francisco. (1998). GiGun. Retrieved Septiembre 15, 2011, from
Montaje del Ensayo de Columna Resonante:
http://www.ing.unal.edu.co/gigun/info/proy_columna.html

Narsilio, G., & Santamarina, C. (2001). Clasifficacin de Suelos: Fundamento

Fsico, Prcticas Actuales y Recomendaciones.

Orfanidis, S. (1996). Introduction to Signal Processing. Prentice-Hall.

Pemex. (2009). Presiones de Formacin. Retrieved 08 15, 2012, from

http://es.scribd.com/doc/20227315/Presiones-de-Formacion

Pestana & Whittle. (1995). Compression Model for Cohesionless Soils.

Geotechnique.

Pestana. (2006). Small-Strain Behaviour of Granular Soils I: Model for

Cemented and Uncemented sands and Gravels.

Presti, L., & Turco. (n.d.). Low Strain Stifness and Material Damping Ratio
Coupling in Soils.

Rayleigh, L. (1885). On Waves Propagated Along the Plane Surface of an

Elastic Solid. Londres.

Richart, & Woods. (1970). Vibrations of Soil and Foundations.

Richart, F. E. (1963). Elastic Wave Velocities in Granular Soils. J. Soil Mech.

Found.

Rinaldi, V. A. (n.d.). Ingeniera Geotcnica Aplicada a la Minera. Retrieved

06 20, 2012, from
icc.ucv.cl/geotecnia/18_ciclo.../clases/01ensayosgeofisicos.ppt

Sandox. (2010). Sandox. Retrieved 06 22, 2012, from Sandox:

http://www.sandoxcientifica.com/geotecnia.html

105
Santamaria, & Fratta. (2005). Discrete Signals and Inverse Problems.

Santibaez, D. H. (2006). Determinacin del Potencial de Licuefaccin de

Suelos no Cohesivos Saturados Bajo Cargas Ssmicas Usando el
Ensayo de Penetracin Estndar. Retrieved 08 21, 2012, from
http://es.scribd.com/doc/22621695/52/Modulo-de-corte-y-razon-de-
amortiguamiento

Sanz, J. (1975). Mecnica de Suelos, Reunin de Ingenieros. Barcelona:

Editions Eyrolles.

Seed, & Idriss. (1970). Soil Moduli and Damping Factors for Dynamic
Responce Analyses.

Shannon, W. (1959). Dynamic Triaxial Test n Sand. First Panamerican

Conference on Soil Mechanics and Foundation Engineering. Mexico
City.

Suarez, J. (n.d.). Ensayos In Situ. Retrieved 06 20, 2012, from

ocw.uis.edu.co/ingenieria.../ENSAYOS_DE_SUELOS_INSITU.pdf

Suriol, J., Lloret, A., & Josa, A. (2007). Reconocimiento Geotcnico del
Terreno.

Tatsuoka, F., & Silver, L. (1980, March). New Method for the Calibration of
the Inertia of Resonant Column Devices. Geotechnical Testing Journal,
3(1), 30-34.

Universidad Politcnica de Catalunya. (2004). UPCommons Portal dAccs

Obert al Coneixement de la UPC. Retrieved 08 08, 2012, from
UPCommons Portal dAccs Obert al Coneixement de la UPC:
http://upcommons.upc.edu/pfc/bitstream/2099.1/3363/6/36067-6.pdf

Whittle, A. (1995). Compression Model for Cohesionless Soils. Geotechnique.

Woods, R. (1978). Measurement of Dynamic Soil Properties. Earthquake

Engineering and Soil Dynamics. Pasadena, CA: ASCE.

Zhang, X. (1993). Unified Dynamic Shear Moduli and Damping Ratios of

Sand and Clay. Soils and Foundations.

106