GRLWEAP™

Principios Fundamentales de la Ecuación de la Onda

Frank Rausche, Jorge Beim

2007-2017, GRL Engineers, Inc.
CONTENIDO
• Antecedentes.
• El modelo de Ecuación de Onda de la “varilla”.
• El análisis numérico de la Ecuación de Onda.
• El modelo del suelo en la Ecuación de Onda.
• El flujo del programa.
 Aplicación del GRLWEAP

• ¿Cuándo?
• Antes de comenzar la hinca de los pilotes.
• O después de haber realizado las pruebas iniciales de pilotes
(análisis refinado).
• ¿Para qué?
• Formular el criterio de hinca.
• Seleccionar el equipo de hinca adecuado.
• Determinar los esfuerzos en el pilote.
• Calcular el descenso final o el conteo de golpes.
• Determinar la relación entre conteo de golpes y capacidad.
 HISTORIA DEL ANÁLISIS DE ECUACIÓN DE LA ONDA
• Concebido en los años 1940 por E.A.L. Smith
• Junto con IBM, Smith desarrolló un programa para hacer los
cálculos: “primera aplicación no militar del computador”.
• Premisa básica: Generar el gráfico de capacidad (“Bearing
Graph”) utilizando modelos mejorados del pilote (llevando
en cuenta la elasticidad) y del suelo (comportamiento
elasto-plástico con amortiguación).
• Permitir cálculo de esfuerzos.
• Mejoramientos posteriores:
• Modelo realístico de martillos Diésel.
• Comparación con medidas reales de campo.
• Desarrollo de constantes de suelo más confiables.
• Análisis Driveability (“hincabilidad”).
• etc.
 Formula Dinámica
Energía disipada en el suelo =
Energía suministrada por el martillo

Ru (s + sl) = ηWr h
sl … desplazamiento “perdido” (empírico o
medido),
η … eficiencia(?) del sistema de hinca
 Formulas comunes en los EUA
1. Formula Engineering News (EN)
Ru (s + sl) = Wr h …... Capacidad última

R = Ru/(FS) …………...……. Carga segura

R = Wr h/[6 (s + sl)] ………. con FS = 6

Usando h en pies y sl = 0.1 pulgadas:

R = 2 Wr h/(s + 0.1) …….. Formula EN

Se puede esperar entonces que la capacidad última sea

Ru = 4 Wr h/(s + 0.1)
Formulas comunes en los EUA
2. Formula de Gates

Ru = 7 (Wrh)½ log(10Golpes/pulgada) - 550

Ru … Capacidad (kN)

Wr … peso del mazo (kN)

h … altura de caída real (m)

log … logaritmo base 10

“Bearing Graphs” a partir de Formulas Dinámicas
Vulcan 506

Bearing Graphs from Formulas

4500
4000
3500
Capacity in kN

3000
2500 Ru-Gates
2000 ENR - inferred
1500
1000
500
0
0 5 10 15 20
Blows /25 mm
 EL MODELO DE LA ECUACIÓN DE ONDA

• El Análisis de Ecuación de Onda calcula el

desplazamiento de cualquier punto de una barra
elástica delgada en cualquier momento.
• El cálculo se basa en las características de la varilla:
• Longitud
• Área de sección
• Módulo elástico
• Densidad de masa
 Fundamentos del GRLWEAP

• Para el análisis de hinca de pilote, la “varilla”

es el sistema Martillo + Sistema de Hinca +
Pilote.

• Se supone que la varilla es elástica(?) y

esbelta(?)

• El suelo está representado por fuerzas de

resistencia que actúan en la interface del
suelo con el pilote.
 Longitud del segmento de pilote

Para posibilitar cálculos realistas, la varilla se

divide en segmentos
• De longitud aproximada ΔL = 1 m (3.3 ft)
• Con masa m=ρAΔL
• Y rigidez k=EA/ΔL

Hay N = L / Δ L segmentos de pilote

Nota: Los segmentos no necesitan tener la misma longitud

 El modelo del pilote y del suelo
Densidad de masa, 
Módulo, E
ΔL= L/N  1m
Área, A

Resorte (resistencia estática)

Masa mi Rigidez ki
Amortiguador (r. dinámica)
Resistencia del suelo

• La resistencia del suelo ralentiza el movimiento

del pilote y causa su rebote.
• Un pilote que se mueve muy lentamente sólo
encuentra resistencia estática.
• Un pilote que se mueve rápidamente también
encuentra resistencia dinámica.
• La resistencia estática a la hinca, sin embargo,
difiere de la resistencia bajo cargas estáticas.
 El modelo del suelo
Ki-1,Rui-1
Segmento
i-1 Ji-1

ki,Rui
SUELO RÍGIDO
ENVOLVIENDO Segmento
LA INTERFACE Ji
SUELO/PILOTE i

ki+1,Rui+1
Segmento
Ji+1
i+1
El modelo del suelo

Resistencia total del suelo

Rtotal = Rsi +Rdi

Segmento
i

ui Fijo
vi
Rsi Modelo de resistencia estática de fuste
Rui Parámetros Rui, qi

Resorte elástico con

ksi = Rui /qi deformación máxima q
(quake)
1 ui
Deslizador plástico
quake, qi rígido con resistencia Rui

Macizo de
suelo
inamovible
La resistencia estática de fuste durante el rebote

Rsi
-Rui Rui
qi

qi

Quake de fuste recomendado:

2.5 mm; 0.1 pulgadas

ui
La resistencia estática de punta durante el rebote
Rsi
Rui
qi
Para martillos
normales de
impacto

qi

Quake de punta recomendado:

1 mm: roca dura
2.5 mm para pilotes no taponados
ui D/120 a D/60 para pilotes taponados
Quakes de punta recomendados

Pilotes no Pilotes taponados

taponados (displacement)
(non- D/120: suelos muy qi Rui Rsi
displacement) densos/duros
2.5 mm D/60: suelos más
blandos/sueltos qi

ui
D
Modelo de resistencia dinámica de fuste
Smith y Smith-viscoso; Parámetro Jsi

Segmento
de pilote Referencia Rdi = RsiJsi vi
i fija (suelo)
Jsi: Factor de
amortiguamiento de
Smith

velocidad v

amortiguador
 Factores de amortiguamiento
recomendados, según Smith

Fuste
Arcilla: 0.65 s/m o 0.20 s/ft
Arena: 0.16 s/m o 0.05 s/ft
Limo: use un valor intermedio
Suelos estratificados: use una media ponderada

Punta
Todos suelos: 0.50 s/m o 0.15 s/ft
 Cálculo del descenso y conteo de golpes

• Una vez la punta del pilote rebote, la máxima penetración

de la punta es conocida, digamos 7.0mm o .27”
• Descenso final = Max desplazamiento punta – Quake
= 7.0 - 2.5 = 4.5 mm
= 0.27 - 0.1 = .17 pulgadas
• El conteo de golpes es el inverso del descenso final:
BCT = 1000/4.5 = 222 golpes/m
BCT = 12/.17 = 71 golpes/ft
Cálculo del descenso – método 1

R Desplazamiento máximo

Calculado
Ru
Extrapolado

Despl.
Descenso final Quake
 Cálculo alternativo del descenso a través de RSA

• Análisis de Esfuerzos Residuales (RSA en la sigla en inglés) a

veces también se llama Análisis de Múltiples Golpes.
• Analiza varios golpes consecutivos con esfuerzos y
desplazamientos iniciales calculados del estado final del
golpe precedente.
• Provee los esfuerzos residuales en el pilote al final del
golpe.
• Su uso en general resulta en mayores descensos (menores
conteos de golpes).
Opción de Esfuerzos Residuales (RSA)
El pilote y el suelo almacenan energía entre los golpes del martillo

Descenso para 2 golpes

Convergencia:
Golpes consecutivos
tienen la misma
compresión/descenso del pilote
 Flujo del programa – Bearing Graph
Distribuya Ru
Establezca constantes
Entrada del suelo Aumente Ru
Incremento de tiempo

Modele martillo y
sistema de hinca Análisis estático
Velocidad del mazo
Análisis dinámico S
¿Aumentar
R u?
Modele Pilote ➢ Esfuerzos
➢ Energías
➢ Velocidades N
Seleccione 1º
Ru Calcule descenso y Salida
conteo de golpes
Bearing Graph: V 506; HP 12x53; Arena/Arcilla
 Comparación de Bearing Graphs

Bearing Graphs from Formulas and Wave Equation

4500
4000
3500
Capacity in kN

Ru-Gates
3000
2500 ENR - inferred
2000 GW-Sand
1500
GW-Clay
1000
500
0
0 5 10 15 20
Blows/25 mm
 Comparación de capacidades

Conteo de golpes 5 8 golpes/25 mm:

Gates 1950 2250 kN
ENR 1930 2570 kN
GW-Arena 1750 1980 kN
GW-Arcilla 1390 1580 kN
 Ecuación de onda:
Diferencia de las formulas
• La ecuación de la onda modela los sistemas componentes
individuales y, por lo tanto, necesita de datos de entrada
más detallados.
• El análisis de ecuación de onda trabaja con capacidades
últimas. Los factores de seguridad tienen que ser aplicados
para determinar las cargas de trabajo.
• Las capacidades determinadas a partir del descenso final se
refieren al instante de la hinca y, por lo tanto, no llevan en
cuenta las perdidas o ganancias con el tiempo.
Ejemplo modalidad “Bearing Graph”
• Vamos a elegir un martillo
• Un pilote típico
• Un perfil de suelo común
 Flujo del programa – Driveability

Entrada Calcule Ru para

1º “gain/loss”

Modele el martillo
y el sistema
de hinca Análisis Próximo G/L

Primera
profundidad de Aumentar
análisis G/L?

Longitud y
N
modelo del pilote
Aumentar Salida
Prof.?
Aumente
profundidad N
Driveability
La modalidad Driveability no trabaja con la
resistencia estática de largo plazo

• ¿CUÁL ES EL RU DURANTE LA HINCA?

• ¿Vamos a tener recuperación (setup) del suelo a lo
largo del fuste con el tiempo? Podemos ganar
25% en arena, ¡o hasta 300% en arcilla!
• Tal vez sea necesario analizar con máxima pérdida
de fricción, o con pérdida parcial, o sin ninguna
pérdida.
 Análisis modalidad Driveability con
tiempo de recuperación (setup) variable

Tiempo de setup
• Factor de setup, SF Ru

• Tiempo de setup
Ru/SF

Tiempo
Ru
• Distancia límite
Ru/SF
Distancia
límite
Desplazamiento
Factores de recuperación (setup)
Soil Type Setup Factor
Clay 2
Silt – Clay 1
Silt 1.5
Sand – Clay 1.2
Fine Sand 1
Sand - Gravel 1
 Análisis modalidad Inspector Chart
JB Instrumentacao Eletronica Ltda 08-May-2014
Inspector chart example GRLWEAP Version 2010

20 20
BANUT S8000
Capacity 2000.0 kN
Ram Weight 78.50 kN
Efficiency 0.950
16 16

Compressive Stress (MPa)

Helm et Weight 1.96 kN

Tension Stress (MPa)

Ham m er Cushion 706 kN/m m

• Capacidad fija
Pile Cushion 383 kN/m m
COR of P.C. 0.500
12 12 Skin Quake 2.540 mm
Toe Quake 6.666 mm

• Altura de caída 8 8
Skin Dam ping
Toe Dam ping
Pile Length
0.494
0.490
30.00 m
sec/m
sec/m

Pile Penetration 30.00 m

(stroke) variable Pile Top Area 2736.33 cm 2

Skin Friction

• Para cada stroke,

4 4 Pile Model Distribution

el programa calcula 0 0

2.00

el descenso (o
conteo de golpes) y 1.60

esfuerzos.
Stroke (m)

1.20

0.80

0.40

0.00
0 40 80 120 160 200 240
Res. Shaft = 10 %
Set (m m /10 bl) (Proportional)
 RESUMEN
• GRLWEAP se basa en el modelo de Smith con extensiones
importantes tales como:
• Modelos realistas del martillo (ECH, OED, CED, VIB)
• Modelos de resortes no lineales para interfaces y slacks (holguras)
• Modelos alternativos del suelo
• Análisis de esfuerzos residuales
• El análisis de ecuación de onda trabaja con la “resistencia
estática a la hinca” (SRD en la sigla en inglés), más un
amortiguamiento (damping) o resistencia dinámica.
• Dos importantes modalidades adicionales de análisis son el
“Driveability” y el “Inspector’s Chart”.