CONDICIONES DE CONTORNO

CONDICIONES DE CONTORNO

Valores de las variables que se prescriben en el contorno del mallado.

Pueden ser:

 Reales (p.ej., superficie del terreno, contorno de un talud)

 Artificiales: limitación del tamaño del modelo y del número de zonas

- Planos de simetría

- Líneas de contorno: truncamiento del modelo suficientemente alejado

 IDENTIFICACIÓN DE LA SIMETRÍA

Sin simetría 1/2 simetría 1/4 simetría 1/8 simetría

Planos de Simetría respecto de:

•Geometría
•Estado de tensión in situ
•Ley del material (grado de anisotropía)
•Materiales no homogéneos (interfases, inclusiones, etc.)
•Secuencia de modelización (excavación, relleno)
•Condiciones de contorno
•Medidas del sostenimiento
•Agua / Temperatura (acoplamientos)
 CLASIFICACIÓN
Las Condiciones de Contorno Mecánicas en un modelo numérico pueden ser:

 Condiciones de contorno de esfuerzos: fijar el estado tensional

con comando APPLY

 Condiciones de contorno de desplazamientos: fijar la velocidad de los

nodos con comando FIX
o comando APPLY

 Condiciones de contorno periódicas: asociar nodos opuestos

con comando ATTACH

 CONDICIONES DE CONTORNO DE TENSIONES

APPLY keyword value < gradient gx gy gz> <range …>

Aplica una condición a cualquier contorno del modelo, sea externo o interno. Las condiciones

de esfuerzo aplicadas a caras son:

keyword sxx: componente xx del tensor de esfuerzos aplicada en el borde del modelo

syy: componente yy del tensor de esfuerzos aplicada en el borde del modelo

szz: componente zz del tensor de esfuerzos aplicada en el borde del modelo

sxy: componente xy del tensor de esfuerzos aplicada en el borde del modelo

sxz: componente xz del tensor de esfuerzos aplicada en el borde del modelo

syz: componente yz del tensor de esfuerzos aplicada en el borde del modelo

nstress: componente de esfuerzos aplicada en dirección normal al plano de borde

dstress: componente de esfuerzos aplicada en dirección del dip direction del plano de borde

sstress: componente de esfuerzos aplicada en dirección del rumbo del plano de borde

value valor que toma la variable en el origen

gradient gx: incremento lineal de la variable en la dirección x Ecuación del gradiente:

Szz Szz 0  x gx  y gy  z gz
gy: incremento lineal de la variable en la dirección y

gz: incremento lineal de la variable en la dirección z

range Para planos paralelos a los ejes, se recomienda utilizar rangos tipo x, y o z. Para planos inclinados, se recomienda

emplear el rango plane, con sus correspondientes keywords.

 CONDICIONES DE CONTORNO DE FUERZAS

APPLY keyword value < gradient gx gy gz> <range …>

Las condiciones de fuerzas aplicadas a nodos son:

keyword xforce: componente x de la fuerza aplicada al nodo

yforce: componente y de la fuerza aplicada al nodo

zforce: componente z de la fuerza aplicada al nodo

value valor que toma la variable en el origen

gradient gx: incremento lineal de la variable en la dirección x

gy: incremento lineal de la variable en la dirección y

gz: incremento lineal de la variable en la dirección z

range Para planos paralelos a los ejes, se recomienda utilizar rangos tipo x, y o z. Para planos inclinados, se recomienda

emplear el rango plane, con sus correspondientes keywords.

CONDICIONES DE CONTORNO DE DESPLAZAMIENTOS

APPLY keyword value < gradient gx gy gz> <range …>

Aplica una condición a cualquier contorno del modelo, sea externo o interno:

keyword nvel: componente de la velocidad aplicada en la dirección normal

dvel: componente de la velocidad aplicada en la dirección de buzamiento

svel: componente de la velocidad aplicada en la dirección del rumbo

value valor que toma la variable en el origen

gradient gx: incremento lineal de la variable en la dirección x

gy: incremento lineal de la variable en la dirección y

gz: incremento lineal de la variable en la dirección z

 CONDICIONES DE CONTORNO DE DESPLAZAMIENTOS

FIX keyword <range ... >

Aplica una condición de contorno de desplazamientos, fijando el valor de la variable

especificada en el valor actual (dado que opera sobre velocidades, es necesario

inicializar previamente las velocidades a cero, para fijarlas posteriormente).

keywords x fija la componente X del vector velocidad en el rango seleccionado

y fija la componente Y del vector velocidad en el rango seleccionado

z fija la componente Z del vector velocidad en el rango seleccionado

Esta condición de contorno puede ser anulada mediante el comando FREE, que tiene la

misma sintaxis.
 CONDICIONES DE CONTORNO PERIÓDICAS

“ROLLERS” “ATTACH”
 CONDICIONES DE CONTORNO PERIÓDICAS
Ejemplo de Aplicación a un cruce de galerías

F L A C 3 D 2 .0 0
Step 2716 Model Projection
17:01:51 Mon Apr 17 2000

Center: Rotation:
X: 2.250e+001 X: 90.000
Y: 1.500e+001 Y: 0.000
Z: 3.000e+000 Z: 0.000
Dist: 1.394e+002 Size: 4.437e+001

Plane Origin: Plane Norm al:

X: 0.000e+000 X: 0.000e+000
Y: 0.000e+000 Y: 0.000e+000
Z: 2.000e+000 Z: 1.000e+000

C ontour of SZZ
Plane: on
Gradient Calculation
-1.0734e+007 to -9.0000e+006
-9.0000e+006 to -8.0000e+006
-8.0000e+006 to -7.0000e+006
-7.0000e+006 to -6.0000e+006
-6.0000e+006 to -5.0000e+006
-5.0000e+006 to -4.0000e+006
-4.0000e+006 to -3.0000e+006
-3.0000e+006 to -2.0000e+006
-2.0000e+006 to -1.0222e+006
Interval = 1.0e+006
D isplacement
Plane: on
Maxim um = 7.128e-003
Linestyle

Itasca Consulting Group, Inc.

Minneapolis, MN USA
 CONDICIONES DE CONTORNO PERIÓDICAS
Ejemplo de Aplicación a columnas de jet-grouting

3 .0 0 FL A C 3 D 3 .0 0
e Step 97764 M odel Projection
11:38:18 Wed Feb 22 2006

n: Center: Rotation:
000 X: -5.800e+000 X: 0.000
00 Y: 5.000e-001 Y: 0.000
00 Z: 5.140e+000 Z: 0.000
4 Dist: 2.930e+002 Size: 3.820e+001
2.500
B lock G roup
terreno
solera
solera2
solera3
jet
SE L G eom etry
M agfac = 0.000e+000

Itasca Consulting Group, Inc.

M inneapolis, M N USA
CONDICIONES DE CONTORNO PERIÓDICAS
Ejemplo de Aplicación a columnas de jet-grouting

1m
 CONSIDERACIONES

• Las condiciones de contorno de desplazamiento subestiman los valores

de tensión y deformación calculados, mientras que las de esfuerzos los
sobrestiman

• Las dos tipologías de condiciones de contorno acotan la solución real por

lo que, para independizarse de su influencia, es necesario alejar los bordes
una distancia conveniente.
INFLUENCIA DE LA DISTANCIA DE LOS CONTORNOS EN TÚNELES (1)
P
P a b 2  a 2 (1  2 ) 2
r a : ur a    b
b 2 G b2  a2    (1  2 )
ur a a
a   2
ura  b
P a   1
b  : ura  a
2 G

1.2




1.15 
Ur/Ur_inf

1.1

1.05

1
0 2 4 6 8 10 12 14 16 18 20
Radio contorno/Radio excavación
INFLUENCIA DE LA DISTANCIA DE LOS CONTORNOS EN TÚNELES (2)

P0 a 2 P0 a 3
Cilindro: ur  Esfera: ur 
2 G r 4 G r 2
1.6

1.5

1.4

1.3

1.2

1.1
Ur/Ur_inf

1.0

0.9

0.8

0.7

0.6

0.5

0.4
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Radio contorno / Radio excavacion
Cilindro/Presion Cilindro/Desplazamiento Esfera/Presion Esfera/Desplazamiento
ESTABLECIMIENTO DE LOS CONTORNOS EN TALUDES

Dimensiones del modelo

>w / 2 w >2w

>H / 2

Las dimensiones del modelo de un talud deben ser aproximadamente 2H x 3W

para evitar efectos de contorno.
 EJEMPLO DE APLICACIÓN

rest tunel_pro.sav

fix x ran x -.1 .1

fix y
Rangos por coordenadas
fix z ran z -40.1 -39.9

F L A C 3 D 2 .0 0 F L A C 3 D 2 .0 0
Settings: Model Projection Settings: Model Projection
18:24:23 Tue Dec 12 2000 10:11:34 W ed Jan 10 2001

Center: Rotation: Center: Rotation:

X: 2.000e+001 X: 0.000 X: 2.000e+001 X: 0.000
Y: 5.000e+000 Y: 0.000 Y: 5.000e+000 Y: 0.000
Z: 0.000e+000 Z: 0.000 Z: 0.000e+000 Z: 90.000
Dist: 2.262e+002 Size: 8.999e+001 Dist: 2.412e+002 Size: 9.596e+001

Sketch Sketch
Linestyle Linestyle
Fixity C onditions Fixity C onditions
Linestyle Linestyle

Itasca Consulting Group, Inc. Itasca Consulting Group, Inc.

Minneapolis, MN USA Minneapolis, MN USA

plo set rot 0 0 0 plo set rot 0 0 90

plo sk gpfix red plo sk gpfix red
 EJEMPLO DE APLICACIÓN

apply sxx -.8e6 grad 0. 0. 2.e4 ran x 39.9 40.1

set g 0
step 1 Comprobación
del apply

F L A C 3 D 2 .0 0
Step 1 Model Projection
18:29:14 Tue Dec 12 2000
z=40 sxx=0
Center: Rotation:
X: 2.000e+001 X: 0.000
Y: 5.000e+000 Y: 0.000
Z: 0.000e+000 Z: 0.000
Dist: 2.262e+002 Size: 8.999e+001
1
Sketch
Linestyle 2E4
FA P
Maxim um = 7.500e+006
Linestyle

z=0 sxx=0.8 MPa

z=-40 sxx=1.6 MPa

Itasca Consulting Group, Inc.
Minneapolis, MN USA

plo sk fap red

 EJEMPLO DE APLICACIÓN

;;;Excavacion del avance y simulación del efecto frente

mo null ran group avance

mo null ran group sostenimiento z 0. 10.
apply nstre -.64e6 ran z -.1 40. x .01 10. cyl end1 0. 0. 0. end2 0. 10. 0. rad 6.

F L A C 3 D 2 .0 0 F L A C 3 D 2 .0 0
Settings: Model Projection Step 1 Model Projection
18:42:53 Tue Dec 12 2000 18:43:45 Tue Dec 12 2000

Center: Rotation: Center: Rotation:

X: 6.149e+000 X: 0.000 X: 6.149e+000 X: 0.000
Y: 5.000e+000 Y: 0.000 Y: 5.000e+000 Y: 0.000
Z: 2.218e+000 Z: 0.000 Z: 2.218e+000 Z: 0.000
Dist: 2.262e+002 Size: 1.510e+001 Dist: 2.262e+002 Size: 1.510e+001

B lock G roup Sketch

destroza Linestyle
terreno
sostenim iento FA P
Maxim um = 4.400e+005
Linestyle

Itasca Consulting Group, Inc. Itasca Consulting Group, Inc.

Minneapolis, MN USA Minneapolis, MN USA

plo blo group plo sk fap red