XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Anlisis no lineal de estructuras de cables y aplicaciones en el
SAP 2000 versin 8.
Bachiller: Liliana Lourdes Snchez Fernndez.
PhD : Vctor Snchez Moya.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Introduccin
El anlisis estructural de estructuras de cables, sobre todo para puentes colgantes de gran luz,
requiere tomar en cuenta los efectos no lineales introducidos por los grandes desplazamientos
que se producen en este tipo de estructuras. El calculo de la lecha de montaje de un puente
colgante as como el efecto de cargas concentradas sobre puentes colgantes sin viga de
rigidez, son ejemplos que muestran la necesidad de considerar la no linealidad de la estructura.
El presente trabajo es parte de un trabajo mas completo que tiene por objeto el desarrollo de
puentes colgantes carrozables (vehculos menores) de gran luz basados en el uso de
estructuras de cables principalmente del sistema cuadricable utilizados en acueductos,
oleoductos, mineroductos, etc.
Esta primera parte esta orientada a mostrar mediante ejemplos sencillos los problemas de no
linealidad y la manera de resolverlos por medio de la ultima versin del SAP 2000, que
presenta opciones de clculo para anlisis no lineal geomtrico.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
1.- Ejemplo sencillo de no linealidad geomtrica por grandes desplazamientos:
Primero se describir el comportamiento no lineal mediante una solucin analtica (exacta),
posteriormente se realizara un calculo de la misma de manera iterativa, y finalmente se
realizara el calculo matricial, que es la manera como trabajan los programas de calculo como
el SAP 2000.
1.1.- Calculo analtico.
Para ilustrar lo escrito anteriormente se trabajara con la estructura que se muestra en la figura
1.
Datos:
L/2 =
P=
1 =
E=
A=
z=
Valores
5
14
5
1.61E+07
9.62E-04
0.43744332
Unidades
m
tn.
tn/m
m
m
Figura 1: Barras elsticas sometidas a una fuerza vertical
Donde:
Li = 5.01909919 m.
Li Sen1 = 0.43744332 m.
Del equilibrio en la posicin deformada es:
2TSen 2 = P
T=
PL f
P
=
2 Sen 2 2( Li Sen 1 + )
.....(1)
Se obtiene T:
Pero Lf segn las relaciones constitutivas se puede expresar como:
.....(2)
L f = Li +
TLi
EA
Y segn las relaciones de compatibilidad Lf se puede expresar como:
.....(3)
L f = ( Li Sen1 + ) 2 + ( Li Cos1 ) 2
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
De las ecuaciones (1),(2) y (3) se obtiene la siguiente expresin:
2 EA
2 EA x 2 + ( LiCos1 ) 2
Li
1 / 2
x=P
Donde
x = ( Li Sen1 + )
Calculando se obtienen los resultados de la tabla:
Lf =
5.037687 m
T=
=
57.35911 m
0.177609 m
Tabla 1: resultados de la solucin analtica
Si verificamos el equilibrio en la posicin inicial obtendremos un valor para T = 80.32 tn. Valor
que difiere del obtenido en la tabla 1.
1.2.- Calculo por iteraciones.
Otra manera de calcular estructuras con no linealidad, es mediante procesos iterativos, para
este ejemplo sencillo se ha realizado dos formas de calculo iterativo, la primera de ellas
consiste en una secuencia de pasos lineales, colocando al sistema un desplazamiento
pequeo u con el cual se pueden establecer relaciones lineales con la geometra del sistema
tal como se muestra en la Figura 2.
Los desplazamientos se incrementan hasta que la carga P que equilibra el sistema sea del
valor que deseamos, es decir P = 14 tn.
Figura 2: Desplazamientos incrementales
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Resultados para este proceso iterativo
u
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.01
0.002
0.002
0.002
0.002
0.0002
0.0002
0.0002
0.0002
0.000025
iteraciones
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
i
0.0873
0.0893
0.0912
0.0932
0.0952
0.0972
0.0992
0.1011
0.1031
0.1051
0.1071
0.1091
0.1110
0.1130
0.1150
0.1170
0.1189
0.1209
0.1213
0.1217
0.1221
0.1225
0.1225
0.1226
0.1226
0.1226
ui
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.1
0.11
0.12
0.13
0.14
0.15
0.16
0.17
0.172
0.174
0.176
0.178
0.1782
0.1784
0.1786
0.1788
0.178825
Ti
2.7101
5.4606
8.2716
11.1429
14.0746
17.0666
20.1188
23.2311
26.4035
29.6359
32.9282
36.2803
39.6921
43.1637
46.6948
50.2855
53.9355
54.6774
55.4215
56.1680
56.9169
56.9920
57.0671
57.1423
57.2175
57.2269
Pi
0.4724
0.9734
1.5072
2.0744
2.6758
3.3119
3.9835
4.6913
5.4359
6.2180
7.0383
7.8974
8.7960
9.7347
10.7142
11.7352
12.7983
13.1884
13.4113
13.6358
13.8622
13.9250
13.9479
13.9707
13.9936
14.0003
Tabla 2: Resultado del proceso iterativo por desplazamientos incrementales.
0.2
0.18
Desplazamientos u
0.16
0.14
0.12
0.1
0.08
0.06
0.04
0.02
0
0
10
15
20
25
30
Numero de iteraciones
Grafico 1: Desplazamientos vs. Numero de iteraciones
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Otra alternativa de calculo es realizar la secuencia de pasos lineales calculando las tensiones
en los elementos considerando la configuracin inicial y posteriormente volver a calcularlas con
la posicin corregida hasta que los valores converjan y se verifique el equilibrio.
Esta alternativa es mas recomendable, puesto que se emplea menor numero de iteraciones,
adems que nos da un resultado mas preciso.
Figura 3: Secuencia de calculo lineal
Los resultados obtenidos para este proceso son, los mostrados en la tabla 4 y en el grafico 2 se
muestra el proceso de convergencia.
Iteraciones
1
2
3
4
5
6
7
8
9
10
11
12
13
i
0.087266
0.133850
0.119802
0.123027
0.122230
0.122423
0.122376
0.122388
0.122385
0.122386
0.122385
0.122385
0.122385
Ti
80.315993
52.453838
58.569570
57.041866
57.412095
57.321700
57.343731
57.338359
57.339669
57.339350
57.339428
57.339409
57.339413
lf
5.04512635
5.03609736
5.03807922
5.03758415
5.03770413
5.03767483
5.03768197
5.03768023
5.03768066
5.03768055
5.03768058
5.03768057
5.03768057
ui
0.23583219
0.16445083
0.18081406
0.17676679
0.17775003
0.17751011
0.17756859
0.17755433
0.17755781
0.17755696
0.17755717
0.17755712
0.17755713
(i+1)
0.13384999
0.11980236
0.12302700
0.12222966
0.12242337
0.12237611
0.12238763
0.12238482
0.12238550
0.12238534
0.12238538
0.12238537
0.12238537
Tabla 3: Resultados obtenidos del proceso iterativo
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Desplazamient o (u
0.25
0.2
0.15
0.1
0.05
0
0
10
12
14
N de i teraci ones
Grafico 2: Desplazamientos vs. Numero de iteraciones
1.3.- Calculo matricial
Se describe el anlisis matricial de estructuras de barra considerando los efectos de no
linealidad para grandes desplazamientos, con el fin de ilustrar el procedimiento con el que
trabaja el SAP 2000.
En la figura 4 se muestra un elemento antes y despus de que haya experimentado
deformaciones, como se aprecia las deformaciones introducen cambios en la longitud del
miembro, producen una rotacin de cuerpo rgido del elemento y generan una curvatura en el
mismo.
El cambio de longitud y la rotacin de cuerpo rgido del elemento pueden introducirse
fcilmente en la matriz de rigidez en funcin de los desplazamientos en los extremos del
elemento.
Figura 4: Elongacin y rotacin axial de un cuerpo rgido
Los trminos que se introducen a la matriz de rigidez viene ha estar incluidos en la matriz de
rigidez geomtrica, la cual es:
u1 v1 u2
v2
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
0 1 0
1
F 0 1
0 1
[Kg ] =
0
L 1 0 1
0 1 0 1
Entonces la expresin matricial se transforma en:
s1
1
s
2 = EA 0
s3 L 1
0
s4
0 1 0 u1
0 0 0 v1 F
+
0 1 0 u2 L
0 0 0 v2
1 0 1 0 u1
0 1 0 1 v
1
1 0 1 0 u 2
0 1 0 1 v2
Donde S es la matriz de fuerzas que actan en el elemento, la primera matriz es la matriz de
rigidez elstica Ke, La segunda matriz es la matriz de rigidez geomtrica Kg. Se aplicara el
anlisis matricial al ejemplo sencillo, que se vio anteriormente.
Para ello se deben tener en cuenta la orientacin de los elementos, puesto que las matrices de
rigidez mostradas anteriormente han sido obtenidas para un elemento cuya posicin inicial es
la horizontal (ver Figura4).
0 2T
0
0 2 EA Cos ( + ) 2
+
P = L
2
0
(
)
Sen
u 2 L
Ensamblando y simplificando obtenemos:
1 0 0
0 1 u
Figura 5
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Efectuamos el calculo de acuerdo a la expresin de arriba y teniendo en cuenta la notacin de
la figura 5:
Z
Sen(1)
-0.437443
-0.736072
-0.635942
-0.615797
-0.615002
-0.615000
5.019099
5.053890
5.040280
5.037778
5.037681
5.037681
0.087156
0.145645
0.126172
0.122236
0.122080
0.122080
0.000000
0.006932
0.004220
0.003722
0.003702
0.003702
0.000
107.358
65.361
57.640
57.340
57.339
46.881
172.501
123.772
114.756
114.406
114.405
-14.000
17.272
2.493
0.091
0.000
0.000
-0.298628
0.100129
0.020145
0.000795
0.000001
0.000000
-0.298628
-0.198499
-0.178354
-0.177558
-0.177557
-0.177557
Tabla5: Resultados obtenidos del calculo matricial
Donde:
Lf L
L
R = P + 2TSen 1
k=
2 EA
2T
Sen 2 1 +
L
L
= R/k
1.4.- Uso del SAP 2000 versin 8 en el anlisis no lineal de estructuras sometidas a
grandes desplazamientos.
Las versiones anteriores de SAP 2000 permitan realizar anlisis no lineal geomtrico de
estructuras por pandeo, a lo que se le conoce como el efecto P- delta, esta alternativa de
anlisis no contempla el caso en el que la no linealidad es originada por los grandes
desplazamientos.
En Diciembre del ao 2002 ingreso al mercado la versin 8 del programa SAP 2000, la cual
presenta la alternativa de calculo no lineal para estructuras sometidas a grandes
desplazamientos, a si como el anlisis por etapas o procesos constructivos.
Para apreciar las ventajas que presenta la ultima versin de SAP 2000 frente al calculo no
lineal y anlisis de estructuras de cables, se presentara el anlisis del ejemplo anterior que es
una clara ilustracin de no linealidad geomtrica debida a grandes desplazamientos, as como
se vera la aplicacin practica de este programa en el calculo de cables colgantes.
Se definieron en el programa el tipo de material y la seccin con la que se va trabajar, de este
modo se ingresan los valores de E =1.61e7 tn / m2 y A = 9.62E-04 m2, y se define el tipo de
anlisis que se va emplear , en este caso, es anlisis no lineal con largos desplazamiento, los
resultados que se presentan en las figuras siguientes son el resultado del anlisis no lineal con
largos desplazamientos (figuras 6 y 7).
Anlisis no lineal empleando el SAP 2000
En vista que el anlisis lineal no dio los resultados esperados, puesto que no cumple con la
condicin de equilibrio, se usara un anlisis no lineal para el calculo de este ejemplo sencillo.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Se ha definido en el Entorno del SAP 2000 versin 8, un tipo de anlisis esttico de naturaleza
no lineal, escogindose la opcin de grandes desplazamientos.
Los resultados de este anlisis son los que se muestran en las figuras 6 y 7:
Figura 6: Anlisis no lineal (grandes desplazamientos) en el SAP 2000. Desplazamientos
en el nudo central.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
10
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Figura 7: Anlisis no lineal (Grandes desplazamientos) en el SAP 2000.
Fuerza axial en los elementos
Verificando el equilibrio para los resultados obtenidos
El desplazamiento del nudo 2 es de 0.17746 m, lo que indica que el seno del ngulo que para
la posicin deformada de la estructura es:
Sen( ) =
( z + )
L
El valor de la tensin en el elemento es de T = 57.3083 tn y para que se cumpla el equilibrio
debe de cumplirse que:
P = 2TSen( )
El valor de 2T Sen() debe ser igual a 14 tn.
El valor de 2T Sen() = 13.99 tn. y 13.99 14 , se puede decir que se cumple con el
equilibrio.
2.- Aplicacin al anlisis de cables colgante
El anlisis de cables colgantes presenta no linealidad geomtrica debida a grandes
desplazamientos.
Para ilustrar esto se har una introduccin a la teora de cables flexibles y posteriormente se
realizara un anlisis no lineal con el SAP 2000, comparando los resultados obtenidos con la
solucin analtica a fin de evaluar las ventajas que presenta esta versin del programa para
este tipo de anlisis.
En el SAP 2000 hay dos opciones de calculo para el anlisis no lineal, estas dependen del
comportamiento no lineal de las estructuras y se pueden identificar por:
- El efecto P- delta (grandes esfuerzos).- Este tipo de no linealidad se presenta cuando las
fuerzas y momentos presentes en la estructura son grandes, y las ecuaciones de equilibrio que
se plantean en la geometra inicial y la deformada difieren incluso si las deformaciones sean
muy pequeas.
- El Efecto de los grandes desplazamientos: Este tipo de no linealidad se presenta cuando la
estructura se somete a grandes deformaciones ( en particular grandes traslaciones y
rotaciones), las formulaciones ingenieriles de esfuerzo y deformacin ya no son aplicables, el
equilibrio se evala en la geometra deformada.
2.1.- Teora de los cables flexibles
La ecuacin diferencial de los cables flexibles es:
dy
p
=
dx
T
Gerencia XIV CONIC: ICG Instituto de la0 Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe / [email protected]
Ver indice
11
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Cuando es constante la intensidad de la carga p la descripcin se aproxima a la de un puente
colgante, obtenindose la forma parablica del cable, la cual se expresa e como:
y=
px
2T0
Figura 8: Cable colgante sometido una carga constante por unidad de longitud
Tomando el sistema de coordenadas indicado en la figura 9, las expresiones con las que se
trabaja son:
y=
4f
(lx x 2 )
2
l
8 f 2
S = L +
3 L
H=
S =
wl 2
8f
wL2 L 4 f
+
2 EA 4 f 3L
Figura 9: Puente colgante cuyas torres se encuentran al mismo nivel
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
12
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Donde S es la longitud del cable y S es la elongacin del mismo.
2.2.- Comportamiento frente a cargas uniformemente distribuidas
Figura 10: Posicin de montaje (plomo) y posicin final (azul)
En la figura 10 se observa que inicialmente el cable tiene una configuracin de montaje, la cual
despus de aplicar las cargas se deforma hasta llegar a la configuracin final.
A diferencia de los anlisis convencionales, en el caso de cables colgantes, la incgnita es la
posicin de montaje, puesto que la posicin final es de a cuerdo a los requerimientos en
campo.
2.2.1.- Calculo analtico del cable suspensor
L = 145.5 m.
f = 14.55 m.
Figura 11: Configuracin final de un puente colgante
La figura muestra la posicin final, lo que nos interesa conocer es la posicin de montaje, en
otras palabras la flecha de montaje (f1).
Para ello se va emplear la condicin de compatibilidad, la cual requiere que para diferentes
estados de carga, la longitud final del cable menos su elongacin, es igual a la longitud inicial
del cable (S0)
S1 S1 = S 2 S 2 = S 0
............(4)
Denominaremos a la posicin final como estado (2), el cual esta sometido a una carga
uniformemente distribuida, equivalente al peso del tablero (w2) y a la posicin de montaje como
estado (1), el cual esta sometido a una carga uniformemente distribuida equivalente al peso
propio del cable suspensor (w1)
La siguiente tabla muestra los datos de un ejemplo ilustrativo.
L=
145.5 m
f2 =
14.55 m
w2 =
w1=
E=
A =
772.37 kg/m
29 kg/m
1.61E+10 kg/m
0.0027226 m2
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
13
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Tabla 6: Datos del ejemplo ilustrativo para el anlisis del cable colgante
8 f 2
S = L +
3 L
S =
wL2 L 4 f
+
2 EA 4 f 3L
H=
wl 2
8f
.....(5)
.....(6)
.....(7)
Empleando las expresiones (4),(5) y (6) calculamos S2, S2 y S0
S2 =
S2 =
So =
149.380 m
0.491 m
148.889 m
Y reemplazando las expresiones (5), (6) y (7) en la ecuacin de compatibilidad (4), se obtiene
lo siguiente:
8 3 16w1 L 2
3w1 L3
f1
f1 ( S 0 L) f1
=0
3L
24 EA
24 EA
.....(8)
La expresin (8) es una cbica, donde la incgnita es f1 , que es la flecha de montaje y la que
define la posicin inicial, que deseamos conocer.
El valor de la flecha de montaje f1, se obtiene resolviendo la ecuacin (8); as como la
componente horizontal de la tensin del cable H ( la cual es constante y la denominamos
empuje) se pueden obtener empleando la expresin (7).
Estos valores son:
f1 = 13.64 m.
H2 = 140474.99 kg.
H1 = 5627.45 kg.
Este ejemplo se trabajara en el SAP 2000 versin 8 Haremos empleo de la versatilidad que
presenta este programa para interactuar con Excel y obtener de esta manera la posicin de
montaje, por medio de un proceso iterativo.
Figura 12: Numeracin de los nudos de la estructura
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
14
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
La siguientes tablas muestran las coordenadas finales y las coordenada correspondientes a la
posicin de montaje, despus del proceso iterativo, el calculo efectuado es el correspondiente
a una etapa de carga a la cual he denominado Etapa 1 y corresponde a la accin de cargas
uniformemente distribuidas tales como el peso propio del tablero del puente colgante, la etapa
2 corresponde a la accin de una carga puntual.
El modelo no cuenta con una viga de rigidez, esto porque queremos evaluar nicamente el
comportamiento del cable colgante ante la accin de cargas distribuidas y puntuales.
La figura 13 muestra los desplazamientos para la Etapa 1, la posicin de montaje esta
representada de plomo y la posicin final o deformada esta representada por azul, el
desplazamiento del nudo central del cable suspensor (nudo 16, ver figura 12), es de 0.7471 m.
16
Figura 13: Desplazamientos producidos por la carga uniformemente distribuida
Verificando el equilibrio en el nudo 16 - Etapa1 ( carga uniformemente distribuida)
De lo obtenido en el SAP 2000 se tiene:
H = 140421.6 kg
z = 14.551702 m.
Donde H es la componente horizontal de la tensin en el cable suspensor, z es la distancia
vertical entre el nudo 16 y el apoyo del cable suspensor.
El momento flector en el cable es cero, as que el momento flector en el nudo central (nudo 16)
debe ser cero.
Hz
wL2
=0
8
....(9)
Evaluando en la ecuacin (9) la diferencia no es cero, es 536.88, valor que representa un
porcentaje de error del 0.026% con respecto al valor del momento flector, y para fines del
calculo es aceptable.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
15
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
TABLE: Joint
Joint
Text
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
Coordinates
XorR
m
3.6E-15
4.85
9.7
14.55
19.4
24.25
29.1
33.95
38.8
43.65
48.5
53.35
58.2
63.05
67.9
72.75
77.6
82.45
87.3
92.15
97
101.85
106.7
111.55
116.4
121.25
126.1
130.95
135.8
140.65
145.5
4.85
9.7
14.55
19.4
24.25
29.1
33.95
38.8
43.65
48.5
53.35
58.2
63.05
67.9
72.75
77.6
82.45
87.3
92.15
97
101.85
106.7
111.55
116.4
121.25
126.1
130.95
135.8
140.65
0
145.5
Z
m
14.65
12.77467
11.02867
9.412
7.92467
6.56667
5.338
4.23867
3.26867
2.428
1.71667
1.13467
0.682
0.35867
0.16467
0.1
0.16467
0.35867
0.682
1.13467
1.71667
2.428
3.26867
4.23867
5.338
6.56667
7.92467
9.412
11.02867
12.77467
14.65
-1.8E-15
-3.6E-15
-3.6E-15
-7.1E-15
-7.1E-15
-8.9E-15
-8.9E-15
-8.9E-15
-1.24E-14
-1.24E-14
-1.42E-14
-1.42E-14
-1.78E-14
-1.78E-14
-1.78E-14
-2.13E-14
-2.13E-14
-2.31E-14
-2.31E-14
-2.66E-14
-2.66E-14
-2.66E-14
-2.84E-14
-2.84E-14
-3.2E-14
-3.2E-14
-3.38E-14
-3.38E-14
-3.38E-14
0
-3.55E-14
Tabla 7: Configuracin final
Desplazamientos y coord. De montaje iteracion 2 U1
U3
X
Joint
1
0.000
0.000
0.000
2
-0.019
-0.089
4.869
3
-0.035
-0.173
9.735
4
-0.047
-0.252
14.597
5
-0.056
-0.327
19.456
6
-0.061
-0.395
24.311
7
-0.064
-0.459
29.164
8
-0.064
-0.517
34.014
9
-0.062
-0.568
38.862
10
-0.058
-0.614
43.708
11
-0.051
-0.654
48.551
12
-0.043
-0.686
53.393
13
-0.033
-0.713
58.233
14
-0.023
-0.732
63.073
15
-0.012
-0.742
67.912
16
0.000
-0.746
72.750
17
0.012
-0.742
77.588
18
0.023
-0.732
82.427
19
0.033
-0.713
87.267
20
0.043
-0.686
92.107
21
0.051
-0.654
96.949
22
0.058
-0.614
101.792
23
0.062
-0.568
106.638
24
0.064
-0.517
111.486
25
0.064
-0.459
116.336
26
0.061
-0.395
121.189
27
0.056
-0.327
126.044
28
0.047
-0.252
130.903
29
0.035
-0.173
135.765
30
0.019
-0.089
140.631
31
0.000
0.000
145.500
32
-0.003
-0.098
4.853
33
-0.003
-0.181
9.703
34
-0.002
-0.259
14.552
35
-0.002
-0.333
19.402
36
-0.001
-0.401
24.251
37
-0.001
-0.463
29.101
38
-0.001
-0.520
33.951
39
-0.001
-0.571
38.801
40
0.000
-0.617
43.650
41
0.000
-0.656
48.500
42
0.000
-0.688
53.350
43
0.000
-0.714
58.200
44
0.000
-0.733
63.050
45
0.000
-0.743
67.900
46
0.000
-0.746
72.750
47
0.000
-0.743
77.600
48
0.000
-0.733
82.450
49
0.000
-0.714
87.300
50
0.000
-0.688
92.150
51
0.000
-0.656
97.000
52
0.000
-0.617
101.850
53
0.001
-0.571
106.699
54
0.001
-0.520
111.549
55
0.001
-0.463
116.399
56
0.001
-0.401
121.249
57
0.002
-0.333
126.098
58
0.002
-0.259
130.948
59
0.003
-0.181
135.797
60
0.003
-0.098
140.647
61
-0.004
0.000
0.004
62
0.004
0.000
145.496
Etapa 1
Z
14.650
12.864
11.202
9.664
8.251
6.962
5.797
4.755
3.837
3.042
2.370
1.821
1.395
1.090
0.907
0.846
0.907
1.090
1.395
1.821
2.370
3.042
3.837
4.755
5.797
6.962
8.251
9.664
11.202
12.864
14.650
0.098
0.181
0.259
0.333
0.401
0.463
0.520
0.571
0.617
0.656
0.688
0.714
0.733
0.743
0.746
0.743
0.733
0.714
0.688
0.656
0.617
0.571
0.520
0.463
0.401
0.333
0.259
0.181
0.098
0.000
0.000
Tabla 8: Configuracin de montaje, obtenida en la 2da iteracin
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
16
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Una vez obtenida la posicin de montaje, procedemos a calcular los desplazamientos
producidos por la carga correspondiente a la Etapa 1, esta posicin debera ser igual a la
posicin final, pero por efectos de calculo es una posicin muy aproximada.
La tabla 9 muestra los desplazamientos para la Etapa 1 as como la posicin deformada y la
variacin de esta ultima con respecto a la posicin final.
Una vez realizado el anlisis para cargas uniformemente distribuidas, se calcula los
desplazamientos para una Etapa 2, la cual corresponde a la accin de una carga puntual P =
10000 kg, aplicada en el nudo 10, que corresponde a una distancia de x1 = 43. 65 m.
La tabla 10 muestra los desplazamientos para la Etapa 2 as como la posicin deformada.
Al final de ambos anlisis se evaluara el equilibrio en las posiciones deformadas, con la
finalidad de verificar si se cumple este ultimo y ver cuan preciso es el anlisis no lineal de
estructuras sometidas a grandes desplazamientos con el SAP 2000 versin 8.
La figura 14 muestra los desplazamientos correspondientes a la etapa2.
10
Figura 14: Desplazamientos producidos por la carga uniformemente distribuida mas una carga
puntual.
Verificando el equilibrio en el nudo 10 - Etapa2 (carga uniformemente distribuida + puntual)
Para ello en el nudo 10 que es el nudo de aplicacin de la carga debe de cumplirse que el
momento debe se cero, tal como lo indica la siguiente expresin:
2
( H + h) ( z + w)
P( L x1 ) x1 w( x1 ) wL x1
+
=0
L
2
2
......(10 )
Donde (H+h) es la componente horizontal de la tensin ene le cable suspensor, debida a la
carga uniformemente distribuida (H) y a la carga puntual (h).
El valor de (z+w) corresponde a la distancia vertical entre el nudo 10 y el apoyo del cable
suspensor.
De lo obtenido en el SAP 2000 tenemos:
(H+h) = 155803.3 kg.
(z+w) = 12.949252 m.
Evaluando los valores en la ecuacin (10) se obtiene una diferencia de 220.62, la cual
representa un error de 0.010%, que para fines prcticos es aceptable.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
17
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
ord. De monta
Joint
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
Desplazamientos y coordenadas - Etapa 1
Desplazamientos
Coordenadas desplazadas Coomparacion con P. F
X
Z
U1
U3
X
Z
0.0000
0.0000
0.0000
14.6500
0.0000
0.0000
-0.0193
-0.0896
4.8499
12.7740
0.0001
0.0007
-0.0350
-0.1745
9.6998
11.0273
0.0002
0.0013
-0.0471
-0.2544
14.5498
9.4101
0.0002
0.0019
-0.0560
-0.3291
19.3997
7.9222
0.0003
0.0025
-0.0618
-0.3985
24.2497
6.5637
0.0003
0.0030
-0.0646
-0.4623
29.0997
5.3345
0.0003
0.0035
-0.0647
-0.5204
33.9497
4.2348
0.0003
0.0039
-0.0624
-0.5726
38.7997
3.2644
0.0003
0.0042
-0.0579
-0.6186
43.6497
2.4235
0.0003
0.0045
-0.0514
-0.6584
48.4998
1.7119
0.0002
0.0048
-0.0432
-0.6915
53.3498
1.1297
0.0002
0.0050
-0.0337
-0.7177
58.1998
0.6769
0.0002
0.0051
-0.0230
-0.7367
63.0499
0.3536
0.0001
0.0051
-0.0116
-0.7471
67.9000
0.1600
0.0000
0.0046
0.0000
-0.7506
72.7500
0.0956
0.0000
0.0044
0.0116
-0.7471
77.6000
0.1600
0.0000
0.0046
0.0230
-0.7367
82.4501
0.3536
-0.0001
0.0051
0.0337
-0.7177
87.3002
0.6769
-0.0002
0.0051
0.0432
-0.6915
92.1502
1.1297
-0.0002
0.0050
0.0514
-0.6584
97.0002
1.7119
-0.0002
0.0048
0.0579
-0.6186
101.8503
2.4235
-0.0003
0.0045
0.0624
-0.5726
106.7003
3.2644
-0.0003
0.0042
0.0647
-0.5204
111.5503
4.2348
-0.0003
0.0039
0.0646
-0.4623
116.4003
5.3345
-0.0003
0.0035
0.0618
-0.3985
121.2503
6.5637
-0.0003
0.0030
0.0560
-0.3291
126.1003
7.9222
-0.0003
0.0025
0.0471
-0.2544
130.9502
9.4101
-0.0002
0.0019
0.0350
-0.1745
135.8002
11.0273
-0.0002
0.0013
0.0193
-0.0896
140.6501
12.7740
-0.0001
0.0007
0.0000
0.0000
145.5000
14.6500
0.0000
0.0000
-0.0040
-0.0991
4.8493
-0.0007
0.0007
0.0007
-0.0033
-0.1827
9.6995
-0.0013
0.0005
0.0013
-0.0026
-0.2614
14.5496
-0.0019
0.0004
0.0019
-0.0021
-0.3351
19.3997
-0.0025
0.0003
0.0025
-0.0016
-0.4036
24.2497
-0.0030
0.0003
0.0030
-0.0012
-0.4666
29.0998
-0.0035
0.0002
0.0035
-0.0009
-0.5241
33.9499
-0.0039
0.0001
0.0039
-0.0006
-0.5756
38.7999
-0.0043
0.0001
0.0043
-0.0004
-0.6211
43.6499
-0.0046
0.0001
0.0046
-0.0002
-0.6604
48.5000
-0.0049
0.0000
0.0049
-0.0001
-0.6931
53.3500
-0.0051
0.0000
0.0051
-0.0001
-0.7190
58.2000
-0.0052
0.0000
0.0052
0.0000
-0.7377
63.0500
-0.0052
0.0000
0.0052
0.0000
-0.7476
67.9000
-0.0047
0.0000
0.0047
0.0000
-0.7506
72.7500
-0.0044
0.0000
0.0044
0.0000
-0.7476
77.6000
-0.0047
0.0000
0.0047
0.0000
-0.7377
82.4500
-0.0052
0.0000
0.0052
0.0001
-0.7190
87.3000
-0.0052
0.0000
0.0052
0.0001
-0.6931
92.1500
-0.0051
0.0000
0.0051
0.0002
-0.6604
97.0000
-0.0049
0.0000
0.0049
0.0004
-0.6211
101.8501
-0.0046
-0.0001
0.0046
0.0006
-0.5756
106.7001
-0.0043
-0.0001
0.0043
0.0009
-0.5241
111.5501
-0.0039
-0.0001
0.0039
0.0012
-0.4666
116.4002
-0.0035
-0.0002
0.0035
0.0016
-0.4036
121.2503
-0.0030
-0.0003
0.0030
0.0021
-0.3351
126.1003
-0.0025
-0.0003
0.0025
0.0026
-0.2614
130.9504
-0.0019
-0.0004
0.0019
0.0033
-0.1827
135.8005
-0.0013
-0.0005
0.0013
0.0040
-0.0991
140.6507
-0.0007
-0.0007
0.0007
-0.0050
0.0000
-0.0008
0.0000
0.0008
0.0000
0.0050
0.0000
145.5008
0.0000
-0.0008
0.0000
Tabla 9: Desplazamientos, posicin deformada y variacin de esta ultima respecto a la posicin
final Etapa 1
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
18
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Desplazamientos y coordenadas desplazadas - Etapa 2
rd. De mont Desplazamientos Etapa 2 Coord. desplazadas Etapa 2
U1
U3
X
Z
Joint
1
0.000
0.000
0.000
14.650
2
-0.030
-0.121
4.839
12.743
3
-0.060
-0.249
9.675
10.953
4
-0.089
-0.383
14.508
9.281
5
-0.117
-0.525
19.339
7.727
6
-0.144
-0.673
24.167
6.289
7
-0.170
-0.829
28.994
4.967
8
-0.194
-0.993
33.821
3.762
9
-0.215
-1.165
38.647
2.672
10
-0.234
-1.341
43.474
1.701
11
-0.206
-1.227
48.345
1.143
12
-0.181
-1.115
53.212
0.707
13
-0.160
-1.008
58.074
0.386
14
-0.142
-0.906
62.931
0.184
15
-0.126
-0.808
67.786
0.099
16
-0.113
-0.715
72.637
0.131
17
-0.102
-0.628
77.487
0.279
18
-0.092
-0.546
82.335
0.544
19
-0.085
-0.469
87.182
0.926
20
-0.078
-0.397
92.029
1.424
21
-0.073
-0.331
96.876
2.039
22
-0.068
-0.271
101.725
2.771
23
-0.063
-0.217
106.575
3.620
24
-0.059
-0.169
111.427
4.586
25
-0.054
-0.127
116.282
5.669
26
-0.048
-0.092
121.140
6.871
27
-0.042
-0.062
126.002
8.189
28
-0.034
-0.038
130.869
9.627
29
-0.025
-0.020
135.741
11.182
30
-0.013
-0.007
140.617
12.857
31
0.000
0.000
145.500
14.650
32
-0.099
-0.130
4.754
-0.032
33
-0.099
-0.256
9.604
-0.075
34
-0.099
-0.390
14.454
-0.130
35
-0.098
-0.531
19.303
-0.198
36
-0.099
-0.678
24.153
-0.278
37
-0.099
-0.834
29.002
-0.371
38
-0.100
-0.996
33.851
-0.476
39
-0.101
-1.167
38.699
-0.596
40
-0.103
-1.348
43.548
-0.731
41
-0.105
-1.228
48.395
-0.573
42
-0.107
-1.116
53.243
-0.428
43
-0.109
-1.009
58.091
-0.295
44
-0.110
-0.907
62.940
-0.175
45
-0.112
-0.809
67.788
-0.066
46
-0.113
-0.715
72.637
0.031
47
-0.113
-0.628
77.487
0.115
48
-0.114
-0.547
82.336
0.185
49
-0.114
-0.470
87.186
0.244
50
-0.114
-0.398
92.036
0.290
51
-0.114
-0.333
96.886
0.323
52
-0.114
-0.273
101.735
0.343
53
-0.114
-0.220
106.586
0.351
54
-0.114
-0.173
111.436
0.347
55
-0.113
-0.132
116.286
0.331
56
-0.113
-0.097
121.136
0.304
57
-0.113
-0.068
125.986
0.265
58
-0.112
-0.045
130.835
0.215
59
-0.112
-0.028
135.685
0.154
60
-0.112
-0.016
140.535
0.082
61
-0.100
0.000
-0.096
0.000
62
-0.112
0.000
145.384
0.000
Tabla 10: Desplazamientos, posicin Etapa 2
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
19
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
Conclusiones:
-
Las estructuras de cables, particularmente en el caso de puentes colgantes sin viga
de rigidez, estn sometidas a grandes desplazamientos, y cuyo anlisis requiere
tomar en cuenta los efectos de no linealidad.
La opcin que presenta el SAP 2000 para anlisis no lineal
por grandes
desplazamientos es una herramienta verstil, de gran utilidad y que facilita el anlisis
de estructuras de cables y de puentes colgantes en general.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
20
XIV CONGRESO NACIONAL DE INGENIERIA CIVIL - IQUITOS 2003
Captulo de Ingeniera Civil del Consejo Departamental de Loreto del Colegio de Ingenieros del Per
BIBLIOGRAFA
-
William McGoire, Richard H. Gallagher y Ronald D. Ziemian.
analysis John Wiley & Sons, Inc, New York., 2000.
Matrix structural
Max Irvine, Cable Structures, Dover Publications, Inc. New York.,1992.
Jack Lpez Jara, Anlisis de puentes atirantados, El puente atirantado Yanango,
Universidad Nacional de Ingeniera, 1996.
Larry Segerlind, Applied finite element analysis, John Wiley & Sons, Inc, New
York.1984.
Gerencia XIV CONIC: ICG Instituto de la Construccin y Gerencia
Calle Nueve 1056 Urb. Corpac San Isidro, LIMA PERU / (51 1) 225-9066 / www.construccion.org.pe /
[email protected]Ver indice
21
