Suscríbete a DeepL Pro para poder traducir archivos de mayor tamaño.

Más información disponible en www.DeepL.com/pro.

6. NorSand de calidad industrial (VBA)

Dra. Dawn Shuttle, Ing.

Mike Jefferies, Ing.
Enero, 2015
 La mecánica del suelo tiene que esperar 30 minutos...

Ahora está en el límite de lo que se debe hacer en hojas de

cálculo Excel. Ahora tenemos que hacer las cosas como un
ingeniero ...

◼ Necesidad de conocer algunas técnicas VBA

◼ Empaquetado para usted en "código abierto"

◼ Rápida visión general de VBA y luego te lanzas...

24 de julio de 2023 2
 Malas prácticas, frecuentes en geotecnia.

◼ Programar ecuaciones en una hoja de cálculo es como escribir código

informático
de los años 60... ¿qué significa una variable llamada "C3"?

◼ Los lenguajes informáticos modernos han evolucionado para utilizar

variables en "inglés sencillo", tanto para facilitar la lectura como el
mantenimiento. . . "C3" => "eta_max"

◼ Los lenguajes modernos también permiten comentarios adyacentes a

los nombres de las variables para mejorar la legibilidad.... eta_maxas
por eqn 6.73"

◼ Microsoft ha proporcionado una serie de lenguajes baratos y fáciles de

usar
◼ Quick Basic QB45 (en los viejos tiempos del DOS)
◼ Visual Basic VB3 (Windows 3.11)
 ◼ Ahora incluye un "entorno de desarrollo integrado" con Office: VBA

24 de julio de 2023 3
 Código VBA para OCC

24 de julio de 2023 4
 Referencia VBA

24 de julio de 2023 5
 Cargar programa

◼ Encuentra tu trabajo del Ejercicio 1

◼ Abrir NorTxl2(r10)_Nerlerk_270.xls
◼ ¡Activar macros al abrir !

24 de julio de 2023 6
 Interfaz de modelización

Utilizar para simular con NorSand

Utilizar para ver una p r u e b a

24 de julio de 2023 7
 Estructura del programa VBA

◼ Al pulsar el botón...

◼ Obtener valores de la hoja de cálculo Params & Plots

◼ Borrar hoja de salida

◼ Hacer todas las cosas inteligentes todo en VBA (subrutinas y

funciones)

◼ Colocar los resultados en una única matriz maestra

◼ Imagen de esa matriz de nuevo a la hoja de cálculo Excel Txl

SimResults

◼ Gráficos Excel trazados desde la hoja de salida en Parámetros y

gráficos
24 de julio de 2023 8
Acceso a VBA : Desarrollador ("Alt" +"F11")

"Activar macros"
 Código VBA

24 de julio de 2023 10
 Funciones frente a subrutinas

◼ Las funciones devuelven un único valor basado en sus

valores de entrada
◼ Función Mi (Mc, chi, lode, psi)
◼ Función e_crit(sigM,....)

◼ Las subrutinas no devuelven nada... Los resultados se

guardan en
variables
◼ Se utiliza para la mayoría de los modelos de estados críticos
24 de julio de 2023 11
 Mi mala práctica de programación...

◼ A los discípulos de ComSci les gusta "encapsulación"

◼ Todo lo que se utiliza en una función/subrutina se pasa a ella en la
"llamada"
◼ Práctica de programación muy "segura
◼ Pero, conduce a un código muy largo de leer
◼ StrainRates (ep1, ep2,ep3, sig1, sig2, sig3, NorProp(), Dp, z2, z3)

◼ Las variables globales facilitan el código en inglés

◼ GetPlasticStrainRates
◼ Entradas, propiedades, salidas deben mantenerse "globalmente"
◼ Accesible a cualquier rutina
◼ Pero, puede ser modificado por cualquier rutina
(accidentalmente...)
24 de julio de 2023 12
 EJERCICIO 5 - Calibración NorSand

◼ Aprenda a calibrar NorSand

◼ Utilice Nerlerk 270 como ejemplo
◼ Propiedades según sus valores en el Ejercicio 1
◼ Ahora piensa en la rigidez drenada y en todo el comportamiento no
drenado
◼ Comprender cómo afecta Gmax al comportamiento no drenado

◼ Aprenda a modificar NorSand

◼ NS es un marco con opciones
◼ Mmobilized es un tema abierto en Geotech
◼ La LCR puede adoptar diferentes formas

24 de julio de 2023 13
 Calibrado NorSand - Finalidad

◼ Finalidad de la calibración ?

◼ Posteriormente obtener ψ precisos de CPT.

◼ Cálculos "manuales" de ingeniería (por ejemplo, su /σv0 para análisis de
estabilidad)
◼ Análisis FE / FD

◼ ¿Cuál NO es el objetivo de la calibración?

◼ Conseguir una coincidencia fantástica entre las curvas

experimentales y teóricas para una sola prueba... si esta calibración
no es representativa de la mayoría de las pruebas...

24 de julio de 2023 14
 Calibrado NorSand - Finalidad

◼ Se requiere un ÚNICO CONJUNTO de propiedades

del material
Ya calibrados en el Ejercicio 1: Γ, λ, Μ, Ν, Χ

◼ Las mismas propiedades deberían ajustarse razonablemente bien a todas

las pruebas
◼ En caso necesario, adapte sus "mejores" ajustes a las condiciones in situ.
◼ Marcar las excepciones para un examen más detenido

◼ Abierto para calibración:

◼ Endurecimiento del plástico H
◼ Módulo de cizallamiento elástico: G
◼ Parámetro de estado real de la prueba (error de la prueba de laboratorio)

◼ Esto requiere un JUICIO DE INGENIERÍA.......

y aunque encontrarás que existen ajustes tontos, no hay un único ajuste
 "correcto".
24 de julio de 2023 15
 Datos Nerlerk "Resumen y propiedades"

Nerlerk 270/1 Arena de Golder (1989)

Propiedades del índice Propiedades del suelo
...del Ejercicio 1
D50 270 micras Gamma = 0,8 con estos valores indicados en el cuadro "Prope
multas = 1.90% lambda10 = 0,1 con estos valores indicados en el cuadro "Prope
emin = 0.536 Mtc = 1,20 líneas de tendencia con estos valores mostrados en el
"Suelo Prope
SG = 2.66 Ntc = 0,20 con estos valores indicados en el cuadro "Prope
chi_tc = 3,00 líneas de tendencia con estos valores mostrados en el
"Suelo Prope
véase el cuadro 3 del
informe
Como probado Inicial A dilatación máxima (=Dmin) al final de la prueba (no
CS)
Test No. ID de p0 e0 psi0 Dmin eta_max psi p e
1 prueba 200.1 0.714 0.144
CID_G151
CID-G152 no hay datos digitales
2 CID_G154 disponibles 0.127
49.5 0.758
3 CID_G155 501.3 0.639 0.109
4 CID-G156 200.9 0.660 0.090
5 CID_G157 200.0 0.592 0.022

Como probado Inicial estado crítico al final de la prueba (no

e CS)
n
p0 e0 psi0 p_c e_c p e
6 CIU_G101 500.0 0.804 0.274
7 CIU_G103 501.1 0.813 0.283
8 CIU_G105 501.0 0.777 0.247
9 CIU_G107 699.9 0.747 0.232
 10 CIU_G108 498.1 0.793 0.263
CIU-G104 700 0.761 0.246

24 de julio de 16
2023
Modelización de las pruebas de drenaje

3: Variar H0 para ajustar la prueba

2: Propiedades elásticas - establecidas a partir de
d a t o s geofísicos

1: variables de "estado" - elegir para la prueba las

condiciones iniciales
cambio menor en ψ del valor de laboratorio ok
24 de julio de 2023 17
 1: Elección de las variables de estado iniciales......
Propiedades del suelo....

CSL G= 0.7 ---

l= 0,030 en base e
para CSL curvo a= 0.97
ec = a - b (p/p_ref)^c b= 0.21
c= 0.28
Idealización CSL como semilog Cambiar el tipo
de CSL
Plasticidad
Mtc = 1.2
N =0 ,20 (típ. 0,2 - 0,3)
ctc =3 ,00 (a menudo se toma como 4,0)
H0 = 50
Hy = 200
(H = H0 - Hy .y ...) 80 típicamente ~ 67
Elasticidad
57
Δψ0 ± 0.02
@ p0 =
Gmáx 57 MPa
G_exp =1 exponente
elástico
n= 0.15
(Ir ...) 285 --- pequeños cambios
k= 0.005 ---
Estado inicial del suelo...
y 0= -0.15
a ψ0 permitido
0.609
=> e0 =
p0 200kPa η0 = 0,0 para
biash anisótropo =0
0
(sig1...)
.00 ---
200 estas pruebas
OCR ("R") =1 .00 ---
24 de julio de 2023 18 puede aumentar el
OCR
para CIU en caso necesario
 2: Módulo de cizallamiento elástico Nerlerk
Módulo de cizallamiento elástico: MPa
0 50 100 150 200
0
G de Seismic VSP, Tarsuit P-45

G de SBP descarga-recarga, Tarsuit P-45

5
G de SBP descarga-recarga, Amauligak I-
Profundidad por debajo de la

65 Nerlerk Gmax tendencia

superficie del núcleo: m

10

15

20 ⎛ ⎞0.8
p'
( )
Núcle
o
G enMPa = pref ⎟
⎠
66⎜
Berma

25
⎝
24 de julio de 2023 19
 3: Calibración - ajuste H...
Propiedades del suelo....

CSL G= 0.7 ---

l= 0.030 en base e
para CSL curvo a= 0.97
ec = a - b (p/p_ref)^c b= 0.21
c= 0.28
Idealización CSL como semilog Cambiar el tipo
de CSL
Plasticidad
Mtc = 1.2
Primera N =0 ,20 (típ. 0,2 - 0,3)
ctc =3 ,00 (a menudo se toma como 4,0)
pasada elige H0 = 50
Hy = 200

Hψ = 0
(H = H0 - Hy .y ...) 80 típicamente ~ 67
Elasticidad
Gmáx @ p0 = 57 MPa 57
G_exp =1 exponente elástico
n= 0.15
(Ir ...) 285 ---
k= 0.005 ---
Estado inicial del suelo...
y 0= -0.15
=> e0 = 0.609
p0 =200 kPa
sesgo anisotrópico h0 = 0.00 ---
(sig1...) 200
OCR ("R") =1 ,00 ---
24 de julio de 2023 20
 Ejemplo de ajuste de la compresión drenada
800
NorSand
800

700

600 600
tensión desviadora, q:
500

400 400

300
kPa

200 200
NorSand
100
CID_G156

0 0
0 5 10 15 20
deformación
axial: %.
1

-1

-2
volumétrica: %.
deformación

-3

-4

-5

24 de julio de 2023 21
 Calibración NorSand - finalmente a H......
Propiedades del suelo....

CSL G= 0.7 ---

l= 0,030 en base e
para CSL curvo a= 0.97
ec = a - b (p/p_ref)^c b= 0.21
c= 0.28
Idealización CSL como .semilog Cambiar el tipo de
CSL
Plasticidad
Mtc = 1.2
N= 0,20 (típ. 0,2 - 0,3)
ctc = 3,00 (a menudo se toma como 4,0)
H0 = 50
Hy = 200
(H = H0 - Hy .y ...) 80 normalmente ~ 67
Elasticidad
Gmax @ p0 57 MPa 57
= G_exp 1 exponente elástico
= 0.15
n= 285 ---
(Ir ...) 0.005 ---
Estado inicial del k=
suelo... -0.15
y 0= 0.609
=> e0 = 200 kPa
p0 =
biash anisotrópico 0 = 0.00 ---
(sig1...) OCR 200
("R") = 1.00 ---

24 de julio de 2023 22
 Calibración NorSand - H

Conjuntos de parámetros utilizados en las simulaciones (pegue los "valores" de la hoja "params & plots" para crear un registro de
simulaciones)
Conjunto de datos: 1 2 3 4 5 6 7 8 9
CID_G151 CID_G154 CID_G155 CID_G156 CID_G157 CIU_G101 CIU_G103 CIU_G105 CIU_G10
Mtc =
N=
c tc =
H0 =
Hy =
= H0 - Hy .y ...)
HCAP =
En la hoja "Resumen y propiedades
Gmax @ p0 =
G_exp =
n= Pegar como "sólo valores"
(Ir ...)
k=

y 0=
=> e0 =
p0 =

=
Observa que debajo de esta pasta hay
una pista sobre Δψ0 y un gráfico de Δψ0
K0

(sig1
OC

frente a p'

24 de julio de 2023 23
 Comprobación del sesgo en el error de prueba ψ

Resumen del "error" de ajuste que muestra las tolerancias de laboratorio

0.030
0 0.02
0 Δψ± 0,02 800 0.02
0.020
0 -0.02
800 -0.02
0.010
Error en psi0 (-

Precisión de la
medición del 0.000
contenido de agua 0 100 200 300 400 500 600 700 800 900
mediante
)

procedimientos
ASTM

-0.010
Sesgo en el error con p′ = CSL subóptimo
-0.020

-0.030
Tensión media efectiva inicial, p' (kPa)
24 de julio de 2023 24
 Procedimiento de calibrado: pruebas en vacío

◼ Empezar con pruebas drenadas

◼ Mejor ajuste cada uno a su vez

◼ Copie los parámetros de ajuste en "resumen" para ver las tendencias

emergentes

◼ Guarde una copia de su ajuste utilizando 'pegar especial' + 'metafile

mejorado'.

◼ Cuando esté satisfecho con el escurrido, pase al no escurrido.

◼ Debe conservar la tendencia H del drenado
◼ Gmax abierto a preguntas
24 de julio de 2023 25
 Ejemplo de resultado de calibración y tendencias
800 800

700

600 600

500 Línea de tendencia ajustada

400
NorSand
400

300
utilizando H0 , Hψ
250

CID_G156
200

endurecimiento, H
200 200

100
150
0 0

módulo de
0 5 10 15 20 0 100 200 300 400 500 600
deformación axial: %. tensión media efectiva, p': kPa 100
1 0.9

50
0
0.8

-1
0
-0.25-0.2-0.15-0.1-0.05 0 0.05 0.1 0.15
volumétrica: %.

índice de
deformación

parámetro de estado inicial

vacío

0.7
-2
CSL

-3
0.6

-4

0.5
-5 10 100 1000
tensión media efectiva, p': kPa

24 de julio de 2023 26
 Te toca calibrar
pruebas drenadas Pruebas 1 - 5

24 de julio de 2023 27
 Pruebas no drenadas...

Gmax y OCR

24 de julio de 2023 28
Modelización de ensayos no drenados

2: Utilizar la tendencia H de la calibración drenada para

comenzar
3: Variar Gmax para ajustar la trayectoria de
tensión y la fragilidad
1: variables de "estado" - elegir para la prueba las
condiciones iniciales
cambio menor en ψ del valor de laboratorio ok
24 de julio de 2023 29
 Efecto de Gmax
400 40

350
Utilización de Gmax a partir de 350

datos geofísicos
300 300
tensión desviadora, q:

250 250

200 200
NorSand
150 150
CIU_G103
kPa

400 400
100 100

50 50
350
Usando G < Gmax para ajustar la 350

prueba
300
tensión desviadora, q:
0 0

300
0 5 10 15 20 0 100 200 300 400 50
250 250
deformación axial: %. tensión media efectiva, p': kPa

200 200
NorSand
kPa

150
CIU_G10

100

50

0
0

24 de julio de 2023 30
 En realidad el problema no es Gmax...

◼ Condición no drenada en la descomposición elástico-plástica

. - de p
◼ dee = Þ dp = - K de p
v v v

K = 2(1+
◼ Elasticidad isótropa G
v)
(3 -
6v)
◼ Opciones
◼ Aceptar un comportamiento demasiado
frágil y utilizar Gmax
◼ Reducir κ / λ ~ 0.25
Gmax para ajustar los datos
◼ Reducir la relación de Poisson para
ajustar los datos
24 de julio de 2023 31
 Compresión isotrópica de la arena Erksak

K = p (1+e) / κ

κ / λ = 0,36 κ / λ = 0,25

κ ...es el
comportamiento
carga-descarga en
edómetro o triaxial

κ / λ = 0,27 κ / λ = 0,38

24 de julio de 2023 32
 Relaciones elásticas (isótropo)

G = K (3 - 6ν) /(2 +
4
2ν)
Módulo elas c ra o: K / G

Ticino
Arena
0
0.0 0.1 0.2 0.3 0.4 0.5

Ra de Poisson o , n

24 de julio de 2023 33
 EL PROBLEMA: Arenas
10000
Delta del Fraser, arenas y limos aluviales,

intercalados Arena Erksak, colocada

hidráulicamente (Tarsiut P-45) Arena Ticino,

pruebas de laboratorio (Bellotti et al, 1996)

/p

Playa de residuos (arena y limo)

1000
GmIrax
=
IR

G de κ
elasticida Ir_perfiles.xls
100
10 100 1000
d 10000
Tensión media efectiva, p' (kPa)
24 de julio de 2023 34
 EL PROBLEMA Limos
10000
Delta del Fraser, arenas y limos aluviales,

Coquitlam Silt intercalado, depositado

hidráulicamente

Estériles de limo (Faro)

1000
/p

Lodos de escombrera (limo arenoso)

GmIrax
=

100
IR

Elasticidad
κ
Ir_perfiles.xls
10
10 100 1000 10000
Tensión media efectiva, p' (kPa)
24 de julio de 2023 35
 Típicamente κ/λ circa 0.25 - ajustar usando Gmax

κ/λ

24 de julio de 2023 36
 Wroth tiene el mismo κ/λ...

Λ = (λ-κ)/λ

24 de julio de 2023 37
 ELASTICIDAD: CÓMO PROCEDER

◼ Calibración con pruebas de laboratorio

◼ Pruebas CID... utilizar G = Gmax

◼ Datos CIU para refinar Gmax con ν = 0,15 (0,2 ?)

◼ Ajustar... G ~ Gmax / 3
◼ Kludge de ingeniería, ya que los conocimientos actuales son
deficientes
◼ Necesidad de aprender qué fudge funciona para el suelo en
cuestión.
24 de julio de 2023 38
 Sobreconsolidación

◼ Geológico: la relación entre la mayor tensión efectiva vertical

en
la historia del suelo a la tensión efectiva vertical actual

◼ Mecánica: a qué distancia se encuentra de la superficie de

fluencia

La fluencia en descarga borra parcialmente la memoria de la tensión

máxima pasada Si está dentro de la superficie de fluencia, sólo tiene

tensiones elásticas = trayectoria de tensión vertical

24 de julio de 2023 39
 Definición de OCR

200
180
OCR = pimg / pimg
160
140
120
desviadora, q

100
tensión

80
60
40
20
0
0 50 100 150 200 250 300 350 400 450
tensión media efectiva, p'

24 de julio de 2023 40
 Ejemplo de OCR en datos

160 160

140

120 120
tensión desviadora, q:

100
CIU_L_200
80 80

60
kPa

40 40

20

0 0
0 5 10 15 20 0 50 100 150 200 250
deformación axial: %. tensión media efectiva, p': kPa

24 de julio de 2023 41
 Ejemplo de ajuste

24 de julio de 2023 42
 Calibrado: comprobación de la elasticidad
250.0

200.0

150.0
G a p0 (MPa)

100.0

50.0

0.0
0 100 200 300 400 500 600 700
p0 (kPa)

24 de julio de 2023 43
Te toca calibrar
ensayos no drenados
Ensayos 6 -
10
24 de julio de 2023 44
 Calibración NorSand - pase final

◼ ¡CALIBRACIÓN COMPLETADA !

24 de julio de 2023 45
 Modificar NorSand

24 de julio de 2023 46
 Relación de fricción movilizada

Mtc

Mi

24 de julio de 2023 47
 Ir a la función en VBA y modificar

24 de julio de 2023 48
 Vuelva a ejecutar el código después de cambiar M_psi

¿Cuál es el efecto del cambio de

idealización de M en el
comportamiento calculado del
suelo?

¿Necesita cambiar H y/o

G a los datos de mejor ajuste ?

24 de julio de 2023 49