Coarita Te PDF

UNIVERSIDAD NACIONAL DE INGENIERA
FACULTAD DE INGENIERA CIVIL
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA
TESIS
Para Optar el Ttulo Profesional de:
INGENIERO CIVIL
EVER DENNYS COARITA TINTAYA
Lima- Per
2013
Dedicatoria:
Dedico este trabajo a mis padres por haberme apoya-
do en todo momento, por sus consejos, sus valores y
su amor. A mis hermanos, por estar conmigo y apoyar-
me siempre, los quiero mucho. A todos mis amigos, por
compartir los buenos y malos momentos.

Agradecimientos:
Agradezco a mi asesor de tesis el M.Sc. lng. Leonar-
do Flores Gonzlez, por la orientacin y ayuda que me

brind para la realizacin de esta tesis.
FACULTAD DE INGENIERA-CIVIL NDICE
NDICE
NDICE IV
RESUMEN VI
LISTA DE FIGURAS VIl
LISTA DE TABLAS . IX
LISTA DE SIGLAS Y SMBOLOS X
CAPTULO 1 GENERALIDADES 1
1.1 INTRODUCCIN . . .
1.2 ESTADO DEL ARTE 2

1.3 METODOLOGA 4
1.4. OBJETIVOS . . . 4
1.4.1 Objetivos especficos . 5
1.5 HIPTESIS . .' . . . . 5
1.6 ORGANIZACIN DE LA TESIS 5
CAPTULO 11 FUNDAMENTO TERICO 6
2.1 ELEMENTO TIPO ARMADURA . 6
2.1.1 Formulacin variacional del elemento tipo armadura 6
2.2 ELEMENTO TIPO CABLE . . . . . 19
2.2.1 Hiptesis del elemento tipo cable . 20
2.2.2 Formulacin Lagrangiana del elemento tipo cable 20
2.3 ANLISIS NO LINEAL ESTTICO . 27
2.4 ANLISIS NO LINEAL DINMICO . 27
2.4.1 Matriz de masa . . . . . . . 29
2.4.2 Matriz de amortiguamiento : 29
2.5 MTODOS DE SOLUCIN DE LA ECUACIN NO LINEAL . 29
2.5.1 Mtodo de Newton . . 30
2.5.2 Mtodo de la secante . 30
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA iv

Bach. Coarita Tintaya, Ever Dennys
FACULTAD DE INGENIERA CIVIL NDICE
2.5.3 Mtodo de control de desplazamiento generalizado 31

CAPTULO 111 ALGORITMO DE SOLUCIN . . . . . 33
3.1 MATRIZ DEL ELEMENTO TIPO ARMADURA 33
3.2 MATRIZ DEL ELEMENTO TIPO CABLE 34
3.2.1 Algoritmo del elemento tipo cable 34
3.2.2 Descripcin del algoritmo . . . 39
3.3 ESTRATEGIA DE SOLUCIN 41
3.3.1 Formato de entrada del modelo 41
3.3.2 Algoritmo para el anlisis esttico no lineal 42
3.3.3 Algoritmo para el anlisis dinmico no lineal 49
CAPTULO IV APLICACIONES 52
4.1 ANLISIS ESTTICO . . 52
4.1 .1 Cable suspendido a diferente nivel . 52
4.1.2 Cable suspendido con carga puntual 55
4.1.3 Armadura bidimensional 57
4.1.4 Armadura tridimensional 60
4.1.5 Armadura tridimensional con cables atirantados 65
4.1.6 Armadura tridimensional con cables pre-esforzados 68
4.1.7 Emparrillado tridimensional de cables 71
4.2 ANLISIS DINMICO . 73
4.2.1 Armadura tridimensional 73
4.2.2 Emparrillado de cables . 75
CAPTULO V CONCLUSIONES V RECOMENDACIONES 77
5.1 CONCLUSIONES . . . . 77
5.2 RECOMENDACIONES . 78
BIBLIOGRAFA . . . . . . . . . . 80
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA V

FACULTAD DE INGENIERA CIVIL RESUMEN
RESUMEN
En esta tesis se desarrollan las matrices de rigidez no lineales de los elementos tipo cable
y armadura mediante una formulacin en coordenadas Lagrangianas. Luego, se desarrollan

algoritmos que emplean los elementos propuestos para resolver el sistema de ecuaciones
no lineales de una estructura en un sistema global de referencia, bajo la accin de fuerzas
externas estticas y dinmicas. Finalmente, se presentan algunas aplicaciones relativas a la
Ingeniera Civil.
SUMMARY
In this thesis, the nonlinear stiffness matrices of cable and truss elements using a Lagrangian
formulation are developed. Then, algorithms are developed using the proposed elements to
salve nonlinear equations of a structure in a global reference system, under action of externa!
static and dynamic torces. Finally, some Civil Engineering applications are presented.
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA vi

FACULTAD DE INGENIERA CIVIL LISTA DE FIGURAS
LISTA DE FIGURAS
Figura 2.1 Movimiento de un elemento armadura en el espacio. (Fuente: Yang,

1994) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Figura 2.2 Grados de libertad para un elemento armadura espacial. 8
Figura 2.3 Coordenadas para la catenaria elstica. . . . . . . . . . 20
Figura 2.4 Fuerzas en un segmento del perfil deformado de cable .. 21
Figura 2.5 Lnea de material definida entre dos partculas de un cuerpo en un
espacio. (Fuente: Linero y Garzn, 201 O). . . . . . . . . . 22
Figura 2.6 Mtodo de Newton modificado con incrementos de carga. 30
Figura 2.7 Mtodo de la secante con incrementos de carga. . . . . . 31
Figura 2.8 Caracterstica general de sistemas no lineales. (Fuente: Yang, 1994). 32
Figura 3.1 Diagrama de flujo del elemento cable para Lo. . 36

Figura 3.2 Diagrama de flujo del elemento cable para To. 38
Figura 3.3 Configuracin en un paso de la iteracin. . 39
Figura 3.4 Elemento cable luego de la convergencia. 39
Figura 3.5 Cable con nodos internos. . 40
Figura 3.6 Fuerzas de cada elemento. 40
Figura 3.7 Diagrama de flujo, mtodo de la secante. . 47
Figura 3.8 Diagrama de flujo, mtodo de Newton. . . 48
Figura 3.9 Diagrama de flujo de anlisis dinmico no lineal. . 51
Figura 4.1 Comparacin numrica de aplicacin con apoyos situados a distinto
nivel. . . . . . . . . . . . . . . . . . . . . 53
Figura 4.2 Coordenada y tensin mnima en el punto ms bajo. 54
Figura 4.3 Configuracin inicial de su geometra bajo peso propio. 55
Figura 4.4 Desplazamiento vertical del nodo 2 versus nmero de incrementos de

carga. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE ARMADURA vii

FACULTAD DE INGENIERA CIVIL LISTA DE FIGURAS
Figura 4.5 Tiempo de convergencia versus nmero de incrementos de carga. . 56 .

Figura 4.6 Geometra de armadura bidimensional. . . . . . . . . . . . . . 57
Figura 4.7 Diagrama carga P versus desplazamiento lateral para cada a. 59
Figura 4.8 Relacin de desplazamiento no lineal y lineal para cada valor de a. 60
Figura 4.9 Geometra de torre tridimensional. . . . . . . . . . . . . 60
Figura 4.1 O Geometra con los nodos de la estructura tridimensional. 61

Figura 4.11 Deformacin del anlisis no lineal (no a escala). . . . . . 64
Figura 4.12 Desplazamientos con nmero de incrementos para el nodo 26 de la
armadura, al considerar anlisis no lineal. . . . . . . . . . . . 64
Figura 4.13 Nmero de incrementos de carga versus tiempo de solucin. 65

Figura 4.14 Geometra de torre tridimensional arriostrada. . . . . . . . . . 66
Figura 4.15 Geometra deformada por peso propio de torre tridimensional arrios-
trada (no a escala). . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
Figura 4.16 Desplazamientos en las tres componentes del nodo 26 con diferentes
cantidades de elementos cable .. 68

Figura 4.17 Modelo tridimensional analizado. 68
Figura 4.18 Diagrama de carga versus desplazamiento lateral del nodo 2. 69
Figura 4.19 Envolvente de cada curva de tensin. 69
Figura 4.20 Relacin de las fuerzas ( ~) y desplazamientos (~) versus tensin

de los cables. . . . . . . . . . . . . 70
Figura 4.21 Relacin de rigideces asintticas. . 70
Figura 4.22 Relacin de rigidez inicial y asinttica negativa. 71

Figura 4.23 Geometra del emparrillado de cables. . . . . . 71
Figura 4.24 Configuracin inicial y deformada del emparrillado de cables. . 72

Figura 4.25 Registro de aceleraciones, El Centro 1940. 73
Figura 4.26 Geometra de la armadura tridimensional con masas concentradas. 74
Figura 4.27 Registro de desplazamiento del nodo 19 en direccin x. . 74
Figura 4.28 Emparrillado de cables con masas concentradas. . . . . 75
Figura 4.29 Registro de desplazamiento del nodo 4 en direccin vertical. 76
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA viii

FACULTAD DE INGENIERA CIVIL LISTA DE TABLAS
LISTA DE TABLAS
Tabla 3.1 Formato de entrada de geometra y masa concentrada. 41

Tabla 3.2 Formato de entrada de propiedades. . . . . . . . . . . 42
Tabla 3.3 Formato de entrada de restricciones y fuerzas nodales. 42
Tabla 4.1 Informacin para el problema .. 52

Tabla 4.2 Comparacin de resultados. . . 53
Tabla 4.3 Coordenada, tensin en el punto ms bajo y fuerzas en el apoyo l. 54
Tabla 4.4 Informacin para el cable suspendido. . . . . . . . 55
Tabla 4.5 Comparacin de desplazamientos del nodo 2 (m). 57
Tabla 4.6 Propiedades de la armadura bidimensional. . . 58
Tabla 4.7 Desplazamiento del nodo 1 en cm, con a =0,5. 58

Tabla 4.1 O Propiedades de la armadura tridimensional. 61

Tabla 4.11 Desplazamientos del anlisis lineal (cm). . . 62
Tabla 4.12 Desplazamientos del anlisis no lineal (cm). 63
Tabla 4.13 Propiedades del elemento cable. . . . . . . 65

Tabla 4.14 Desplazamientos del anlisis no lineal por peso propio (cm). 67
Tabla 4.15 Propiedades de los elementos. . . . . . . . . 69
Tabla 4.16 Puntos de inflexin en cada curva de tensin. 70
Tabla 4.17 Informacin para los cables. . . . . . . . . . . 72
Tabla 4.18 Comparacin de desplazamientos en nodo 8 (mm). 72
Tabla 4.19 Informacin de la armadura tridimensional. . . . . . 73
Tabla 4.20 Comparacin de resultados de desplazamientos mximos. 75
Tabla 4.21 Comparacin de resultados de desplazamientos mximos. 76
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE -ARMADURA ix

FACULTAD DE INGENIERA CIVIL LISTA DE SIGLAS Y SMBOLOS
LISTA DE SIGLAS Y SMBOLOS
A rea de seccin transversal
e Matriz de amortiguamiento
Matriz de amortiguamiento viscoso global
Configuracin inicial no deformada
ltima configuracin calculada
Configuracn actual deformada
Matriz de correccin en el caso de T 0 conocido
E Mdulo de elasticidad de Young
F Matriz de flexibilidad del elemento cable
Componentes horizontales en las direcciones x e y de la tensin en el nodo i
Componente vertical en la direccin z de la tensin en el nodo i
Componentes horizontales en las direcciones x e y de la tensin en el nodo j
Componente vertical en la direccin z de la tensin en el nodo j
Fuerzas internas del elemento cable
Fuerzas internas del elemento cable en el incremento de carga i

lp Fuerza axial del elemento armadura en e1
Fuerzas en el elemento en las direcciones x, y y z en C1
Vector de fuerzas
Vector de carga incremental
Vector de direccin del sismo x, y z
Vector de fuerzas externas
Vector de fuerzas internas
Vector de fuerzas internas del elemento cable
Elemento de la matriz de flexibilidad
Vector de masa en la direccin del sismo
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA X

fs Vector de fuerza residual de la estructura
fsj Vector de fuerza residual de la estructura en la iteracin j
h Sub matriz de s 3
Momento de inercia
I,J Nodos del elemento cable
Jm Masa concentrada en el nodo
K Matriz de rigidez
Ka Matriz de rigidez global
Kr Matriz de rigidez tangente
K Matriz de rigidez tangente dinmica
Matriz de rigidez global en el incremento de carga i y la iteracin j
Matriz de rigidez del cable
Matriz de rigidez tangente del cable
Matriz global de rigidez elstica
Matriz de rigidez elstica de un elemento
Matriz de rigidez geomtrica de un elemento
Rigidez inicial de la envolvente
Rigidez asinttica de la envolvente en la parte positiva
Rigidez asinttica de la envolvente en la parte negativa
Longitud del elemento
Longitud inicial del elemento
Longitud inicial del elemento en el incremento de carga i
L Longitud final del elemento cable
Proyecciones horizontales de la longitud inicial del cable

Proyecciones horizontales de la longitud del cable
Proyeccin vertical de la longitud inicial del cable
Proyeccin vertical de la longitud del cable
M Matriz de masas concentradas
Ma Matriz de masa global
m Peso por unidad de longitud
Peso por unidad de longitud en el incremento de carga i
Fuerza axial en el incremento de carga i
Ns Nmero de segmentos del cable
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE -ARMADURA xi

Nacel Nmero de datos del registro ssmico
Nmero de incrementos de carga

P-b. Efecto P-Delta
Fuerza en la inflexin de la envolvente positiva
Fuerza en la inflexin de la envolvente negativa
R Vector de carga actuante
Vector de carga actuante en la iteracin j
Trabajo virtual realizado por una carga externa en C1 y 0 2
S Incremento del esfuerzo axial del elemento armadura
Trminos de orden superior en las matrices de rigidez
T Tensin del cable
Tensin inicial del cable
Tensin del cable en los nodos i y j respectivamente
tol Tolerancia
telem Tipo de elemento
u,v,w Desplazamientos en las dire~ciones de los ejes x, y y z
Ux,Uy,Uz Desplazamientos en las direcciones x, y y z
Vector de desplazamientos en el incremento de carga i
u Vector de desplazamientos
Vector de desplazamientos en el tiempo
Vector de velocidades en el tiempo
Vector de aceleraciones en el tiempo
Ug Aceleracin ssmica en los puntos de apoyo
w Peso por unidad de longitud del cable
w Peso por unidad de longitud luego de la elongacin del cable
Frecuencia circular del modo i
Desplazamiento en la inflexin de la envolvente positiva
Desplazamiento en la inflexin de la envolvente negativa
x,y,z Coordenadas cartesianas globales
Coordenadas cartesianas globales en el incremento de carga i
Xj, Yj, Zj Coordenadas cartesianas globales temporales en la iteracin j
b.Fintj Variacin de las fuerzas internas del elemento cable en la iteracin j
b.fsj Vector incremento de fuerza residual de la estructura en la iteracin j
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA xii

FACULTAD DE INGENIERA CIVIL LISTA DE SIGLAS V SMBOLOS
l::iL Incremento de longitud del elemento cable

l::iNj Variacin de la fuerza axial en la iteracin j
l::iPi Vector de carga externa incremental en el tiempo
l::iPi Vector de carga externa incremental dinmica en el tiempo
l::iT Incremento de la tensin en el cable
l::it Intervalo de tiempo del registro
/:).U, /:).V, /:). W Variacin de desplazamientos de los nodos
l::iui Vector de variacin de desplazamientos en el tiempo
l::iui Vector de variacin de velocidades en el tiempo
l::ii Vector de variacin de aceleraciones en el tiempo
!::iu/ Vector incremento de desplazamientos en incremento de carga i e iteracin j
l::igi Variacin de la aceleraCin ssmica del registro en el tiempo t

Fraccin de la carga vertical
Factor de amortiguamiento proporcional a la masa
Factor de amortiguamiento proporcional a la rigidez
c, e, r Funcin de perturbacin
e Deformacin lineal
Deformacin axial
,/3 Parmetros del mtodo de Newmark
0 Dimetro del elemento
Deformacin no lineal
Fraccin del amortiguamiento crtico del modo i
Factor de aproximacin inicial del cable
T Esfuerzo axial inicial
lllloo Norma infinita
00 Infinito
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA xiii

FACULTAD DE INGENIERA CIVIL CAPTULO 1: GENERALIDADES
CAPTULO 1
GENERALIDADES
1.1. INTRODUCCIN
Muchas estructuras en su estado inicial tienen un comportamiento lineal, sin embargo
existen excepciones donde la relacin entre la fuerza y el desplazamiento no puede ser des-
crita correctamente con un modelo de comportamiento lineal de la estructura. Las fuentes de
no linealidad se deben al comportamiento no lineal del material, a la no linealidad geomtrica
o a un efecto combinado de las anteriores [22]. Dentro de las excepciones, los cables siempre
presentan un comportamiento no lineal, con una no linealidad debida a que slo tienen un
comportamiento tensional y no linealidad geomtrica [6]. En general, las estructuras de cable
manifiestan una flexibilidad que no se puede atribuir a la baja rigidez axial de sus elementos,
sino a la geometra de la estructura [17].
Los mtodos usuales de anlisis de estructuras de cable estn basados en la discretiza-
cin del cable en elementos ms pequeos, para formular las ecuaciones de equilibrio. Luego
se emplean mtodos numricos para la solucin de las ecuaciones resultantes o se emplean

elementos con valores del mdulo de Young adecuados. En el presente trabajo de investiga-
cin se estudian las estructuras de cable dividindolas en elementos cable con matrices de
rigidez deducidas a partir de la ecuacin de la catenaria [11 ].
En los ltimos aos, el elemento cable viene siendo empleado en estructuras tales como
puentes atirantados, torres ventadas, techos de estructuras que cubren amplias reas y otros.
Cabe mencionar que el comportamiento de este elemento siempre es no lineal tanto en el
rango elstico como el inelstico. As mismo, en muchos casos son estructuras que presentan
grandes desplazamientos, por ende se hace necesaria la profundizacin de su estudio en el
rango no lineal.
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA

Bach. Coarita Tintay, Ever Dennys
1.2. ESTADO DEL ARTE
El estudio del elemento cable ha sido de mucho inters durante aos. Es as que Galileo
en 1638 dice que la forma de un cable colgante es parablica, esto por una analoga al vuelo
de un proyectil. En 1691 los hermanos Bernoulli, Leibnitz y Huygens llegaron a la conclusin
que tiene la forma de la catenaria [14] y encontraron las ecuaciones de equilibrio para un
cable inextensible. Leibnitz uso la teora de clculo infinitesimal para derivar la ecuacin de
esta curva [25].
Los hermanos Bernoulli tambin formularon la ecuacin diferencial general de equilibrio
de un elemento cadena bajo varias cargas y que considera los efectos de estiramiento por
incorporacin de la ley de Hooke. Routh en 1891 resolvi dicha ecuacin para una catenaria
elstica suspendida simtricamente y compuesta de un material elstico lineal [1 ], mientras
que Feld extendi el trabajo anterior para el caso no simtrico, pero desafortunadamente las
coordenadas usadas en ambas soluciones hicieron difcil su aplicacin [14].
En 1981 lrvine, para evitar dicha dificultad, adopta un aproximacin Lagrangiana en la
solucin de la catenaria elstica [14] obteniendo una expresin para la matriz de rigidez tan-
gente de un cable elstico no simtrico.
Dentro de las diversas metodologas que han sido propuestas para modelar estructuras
de cable, stas son representadas empleando dos mtodos diferentes.
El primer mtodo est basado en funciones de interpolacin polinomial para describir la
forma y desplazamiento. Este es un mtodo comn en el desarrollo de elementos finitos.
En este contexto, el elemento recto de dos nodos es el elemento ms usado en el modela-

miento de cables [25). Dicho elemento slo posee rigidez axial y generalmente es aplicable
para cables pretensados. Cuando se tienen cables aflojados con una curvatura grande, la
representacin de su geometra se realiza por un gran nmero de elementos barra, anlisis
que se convierte en ineficiente mientras ms grados de libertad existan. Cuando se desea

modelar al cable aflojado, el mdulo de elasticidad se reemplaza por un mdulo equivalente,
el cual toma en cuenta la flexin del cable; esta rigidez equivalente fue propuesta por Ernst
[25]. Cuando la tensin aumenta la flexin disminuye y la rigidez axial aparente aumenta.
Generalmente en un cable que tiene relativamente altos esfuerzos y pequea longitud, el
mdulo equivalente aproximado es recomendable para obtener buenos resultados [23]. Otro
modelo es el elemento isoparamtrico de multiples nodos, el que se consigue al adicionar
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA 2

ms nodos al elemento finito; lo comn son los elementos de tres a cuatro nodos, los cua-
les usan funciones de interpolacin parablica o cbica respectivamente. stos dan buenos
resultados en cables con pequeas deflexiones; sin embargo, para grandes deflexiones se
deben usar ms elementos y la continuidad de las derivadas entre elementos no se cumple
[11 ]. La continuidad de la derivada se puede cumplir agregando grados de libertad rotacional
a los nodos. Dicho modelo se conoce como elementos curvos con grados de rotacin libre;
este elemento fue desarrollado por Gambhir y Batchelor [25] quienes usan polinomios cbi-
cos que describen el desplazamiento y la forma del cable. Por ello los elementos basados
en interpolacin polinomial son generalmente apropiados para modelar cables con pequeas
deflexiones verticales; por lo mencionado anteriormente, para elementos cable con grandes
deflexiones verticales, se emplean muchos elementos para conseguir la geometra curva del
cable, lo que demanda un gran nmero de operaciones aritmticas. Esto influye en el tiempo
de solucin del problema [11 ].
El segundo mtodo usa formulaciones analticas que toman en cuenta los efectos de car-
ga aplicada a lo largo del cable y describe el comportamiento real del cable. En este contexto,
el elemento elstico parablico es generalmente ms usado en el anlisis de estructuras de
cable debido a su forma simple en comparacin a la catenaria [25]. Otro modelo es la ca-
tenaria elstica: asociada, la cual es una extensin del modelo de lrvine, desarrollada por
Ahmadi-Kashani y Bell [6]. En este modelo, una carga constante por unidad de longitud se
considera aplicada a pesar de la elongacin del cable. La elongacin positiva es producida
por un incremento de la carga total sostenida, siendo este modelo utilizable para cables suje-
tos a presiones externas como viento, nieve u otro tipo de fuerzas externas [1 ]. En el presente
trabajo se considera un elemento elstico catenaria, el cual es exacto debido a que la con-
figuracin de equilibrio de un cable colgante es de naturaleza catenaria [23] y este modelo
es apropiado para un cable perfectamente flexible sujeto a cargas de peso propio [25]. La
formulacin del elemento elstico catenaria est basada en la solucin analtica exacta para
la catenaria elstica desarrollada por O'Brien [23]. En este mtodo, slo un elemento cate-
naria de dos nodos sin uniones internas es necesario para modelar un cable; este elemento
puede ser usado para modelar pequeas y grandes deflexiones verticales y horizontales. De
esta manera el cable puede ser representado por un slo elemento, lo que tiene ciertas ven-
tajas, tales como un reducido nmero de grados de libertad y el considerar el efecto no lineal
elstico del cable en todos los casos [11 ]. Gracias a las expresiones analticas se preserva la
continuidad a travs de los lmites del elemento.

Otro modelo usado es el elemento catenario inextensible, el cual es un elemento infini-
tamente rgido axialmente (AE --+ oo) lo cual evita incrementar su longitud, perfectamente
flexible (Eln = 0), libre de rigidez torsional y capaz de soportar slo fuerzas tensionales
[25]. Las aplicaciones basadas en este tipo de elementos encuentran algunas dificultades
debido a que la curva que relaciona fuerza y desplazamiento es asinttica, por lo que tien-
den a experimentar gran inestabilidad numrica, acusando gran dificultad o incluso imposible
convergencia [1].
Dentro de los programas ms usados en nuestro medio se conocen el programa SAP2000
(CSI), que modela al cable como un elemento cable catenaria [5] y el programa MStower V6
(Engineering Systems), que tambin usa un cable catenaria [7].
1.3. METODOLOGA
La metodologa para el desarrollo de la presente tesis es la siguiente:
Recopilacin de informacin existente sobre no linealidad de cables, no linealidad de
elementos armadura, mtodos para resolver ecuaciones no lineales y anlisis de res-
puesta no lineal en tesis, informes de suficiencia, libros, revistas, artculos de investiga-
cin.
Luego de obtener la informacin pertinente se resume toda la informacin, se elabo-
ra el fundamento terico y seguidamente se realizan algoritmos de solucin para los
elementos tipo armadura y cable.
El programa desarrollado se comprueba con modelos y ejemplos ya estudiados en la

bibliografa revisada o con programas de reconocido prestigio y uso en nuestro medio.
Con el programa desarrollado se analizan aplicaciones de estructuras tipo cable y ar-

madura, considerando un comportamiento no lineal elstico y efectos de orden superior.
1.4. OBJETIVOS
El objetivo general de la presente tesis es analizar estructuras tipo armadura y cable
considerando su no linealidad geomtrica, a partir del desarrollo de un elemento cable y un
algoritmo de convergencia global propuesto en el presente trabajo.
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE -ARMADURA 4

1.4.1. Objetivos especficos
Determinar la matriz de rigidez no lineal de un elemento cable.
Determinar la matriz de rigidez no lineal de un elemento armadura.
Elaborar un algoritmo para resolver el problema global no lineal geomtrico.
Presentar ejemplos aplicados a la Ingeniera Civil resueltos con el algoritmo desarrolla-
do en el item anterior, el cual ser implementado en Matlab.
1.5. HIPTESIS
Mediante el enfoque de las expresiones analticas exactas de la catenaria elstica, pode-
mos desarrollar un elemento cable espacial de dos nodos para realizar el anlisis no lineal
elstico de estructuras de cables y combinar con matrices de rigidez no lineales de otros
elementos como armadura. La ecuacin no lineal f = Ku del conjunto cable - armadura se
puede resolver mediante un mtodo numrico incremental e iterativo.
1.6. ORGANIZACIN DE LA TESIS
El captulo 1trata de las generalidades, estado del arte y presenta los objetivos e hiptesis
como parte central de la tesis.
El captulo 11 documenta el marco referencial de manera terica, con el que se deducen
las matrices que representan a la armadura y al cable.
En el captulo 111 se presentan los algoritmos globales y la estrategia de solucin en el
rango no lineal de cables y armaduras, as tambin se propone un algoritmo de inicio para
generar la matriz de rigidez del elemento cable.
En el captulo IV se desarrollan aplicaciones relativas a la Ingeniera Civil con la teora e
implementacin numrica desarrolladas en los captulos 11 y 111 respectivamente.
En el captulo V se indican las principales conclusiones y recomendaciones del presente
trabajo.

FACULTAD DE INGENIERA CIVIL CAPTULO 11: FUNDAMENTO TERICO
CAPTULO 11
FUNDAMENTO TERICO
2.1. ELEMENTO TIPO ARMADURA
Las armaduras son estructuras utilizadas en torres, naves industriales, coberturas en
grandes luces y otros. Los elementos que componen a una armadura trabajan exclusiva-
mente a fuerza axial, ya sea tensin o compresin.
Estos elementos pueden ser de aceros estructurales galvanizados o aleaciones de alu-
minio; materiales que se utilizan en una variedad de formas y tamaos y son seleccionados
por su resistencia, capacidad de soldadura, resistencia a la corrosin, facilidad de trabajo,

factores econmicos, entre otros factores.
2.1.1. Formulacin variacional del elemento tipo armadura
El movimiento de un cuerpo puede ser descrito por tres configuraciones tpicas, tal co-
mo se muestra en la figura (2.1 ), en un anlisis incremental no lineal: la configuracin inicial
no deformada C0 , la configuracin ltima previamente calculada C1 , y la configuracn actual

deformada C2 (deformacin final). En la formulacin Lagrangiana actualizada, la ltima confi-
guracin conocida C 1 es usada como configuracin de referencia para establecer la ecuacin
del equilibrio del cuerpo en la configuracin actual C 2 .
Las ecuaciones de equilibrio de acuerdo al principio de trabajos virtuales pueden ser
escritas para el cuerpo en forma incremental no lineal como [27]:
j v eT E e dV + j v 'TJT r dV + R =
1 2
R (2.1)

oy 1y 2y
2 c2 2x
y~
actual (deformada)
1\~
y~ X
ltima configuracin
calculada
o 1
YL-~ox
Inicial (no deformada)
Figura 2.1: Movimiento de un elemento armadura en el espacio. (Fuente: Yang, 1994).
La expresin anterior es la versin linalizada de la ecuacin incremental de equilibrio pa-
ra describir el movimiento del cuerpo de C 1 a C 2 con referencia a los ejes de la configuracin
inicial no deformada C0 , y con todos los parmetros fsicos referidos a la ltima configuracin
calculada el.
Para un elemento armadura, debemos considerar slo la componente axial de los tenso-
res de esfuerzo y deformacin. Denotando el esfuerzo axial inicial por T, el incremento de
deformacin axial por e y el incremento de esfuerzo axial por S, el cual se puede expresar
como [27]:
S=Ec (2.2)
Descomponiendo el incremento de deformacin e en su parte lineal (e) y parte no lineal
(r), la ecuacin puede ser expresada como:
Para el elemento armadura espacial como en la figura (2.2), los desplazamientos (u, v, w)
internos del elemento pueden ser relacionados al desplazamiento (ui, vi, wi) y (uj, Vj, wj) de
los dos extremos del elemento usando funciones de interpolacin lineal tales como [27]:
U = Ui ( 1 - ~) + Uj ( f) (2.4a)
V = Vi ( 1 - ~) + Vj ( f) (2.4b)
W = Wi ( 1- ~) + Wj (f) (2.4c)

ly
"' J
~L
X
/ .
Figura 2.2: Grados de libertad para un elemento armadura espacial.
Estudios hechos por Jagannathan y otros (1975), en el anlisis no lineal de armaduras
estructurales espaciales, muestran que se requiere una representacin Green-Lagrange de
la deformacin axial (E:) [19]:
2 2 2
_ dudx
1
E:--+-
2
[(du)
dx + (dv)
--dx + (dw)
-dx ] =e+r (2.5)
Cuyos componentes lineal y no lineal para el incremento de deformacin E: pueden ser
escritos como:
du tl.u
e=-=-
dx L (2.6a)
~ ~ ~ [( ~~) 2+ ( ~:) 2+ ( ~) 2] ~ ~ [(~u) 2+ ( ~v) 2+ ( !>Lw) 2] (2.6b)
Donde:
(2.7a)
(2.7b)
(2.7c)
El vector de desplazamiento nodal para el elemento armadura es:
(2.8)
Los vectores de fuerza del elemento en C1 y C2 son respectivamente:
1
JT = [ 1 Fxi 1
Fyi 1
Fzi 1
Fxj 1
Fyj 1
Fzj ] (2.9)
2 2 2
JT = [ 2 Fxi 2
Fyi
2
Fzi
2
Fxj Fyj Fzj ] (2.1 O)

Para un elemento armadura en equilibrio en C1 , todas las fuerzas de corte transversal
desaparecen, es decir:
(2.11)
Y las fuerzas axiales en los dos extremos de la barra son iguales en magnitud, pero
opuestas en signo:
(2.12)
Se supone que slo actan cargas concentradas en los extremos del elemento armadura,
mientras que las cargas distribuidas slo se consideran si son reemplazadas por fuerzas
nodales equivalentes.
Cada uno de los trminos de la ecuacon (2.3), puede ser expresado como:
J V eT E e dV = J VE ~u (j (~u) dV (2.13)
(2.14)
(2.15)
{ c5 TE dV =
Jv TJ TJ
{~E (~u
Jv 2
2
~v
L2 + L2 + L2
2 2
~w ) [~u c5 (~u)
L L
. ~v
+ L
c5 (~v)
L +
~w c5 (~w)]
L L
dV (2.16)
[ L(j
~::!.u (~::!.u) !:l.v (!:l.v) !:l.w (!:l.w)]
Jvr? rdV =
JLlp L + L(j L + L(j L dx (2.17)
(2.18)
(2.19)
Donde q se refiere a las componentes en los nodos i y j.

Teorema 2.1.1. La ecuacin variacional (2. 13) se expresa en forma matricial como:
(2.20)
Demostracin. Al desarrollar la primera ecuacin (2.13):
Es posible expresar lo anterior en funcin del diferencial de longitud, debido a que dV =

Adx y al sustituir (2.7a)
Al integrar:
Al expresar como producto escalar:
-D. u Ui- Uj
o o
o = EAUT
o
D.u L -Ui + Uj
o o
o o
Lo anterior se puede expresar como un producto entre una matriz y el vector de despla-
zamientos u:
1 0 0 -1 0 0 Ui
0 0 0 0 0 0 Vi
o o o o o o
-1 o o 1 o o Uj
0 0 0 0 0 0 Vj
0 0 0 0 0 0 Wj
Entonces:
donde:

ke es la matriz de rigidez elstica, expresada como:
1 o o -1 o o
o o o o o o
EA o o o o o o
.ke=----; (2.21)
-1 o o 1 o o
o o o o o o
o o o o o o
o
j v r? rdV = 8uTk9 u (2.22)
Demostracin. De la ecuacin (2.17):
= Llp2 [~u (~u)+ ~v (~v) + ~w (~w)] J L dx
Al integrar y sustituir (2.7a), (2.7b) y (2.7c):
lp
=L [~u (uj- ui) + ~v (vj- vi)+ ~w (wj- wi)]
Al expresar como producto escalar:
-~u Ui- Uj
-~V Vi- Vj
-~w Wi -Wj
= uTlp
~u L -Ui + Uj
~V -Vi+ Vj
~w -Wi +wj

Se puede expresar como un producto de una matriz y el vector de desplazamientos u:
1 o o -1 o o
o 1 o o -1 o
= uTlp o o 1 o o -1
L -1 o o 1 o o uJ
o -1 o o 1 o
o o -1 o o 1
Entonces:
j vr? rdV = u7k u 9
donde:
k9 es la matriz de rigidez geomtrica:
1 o o -1 o o
o 1 o o -1 o
lp o o 1 o o -1
kg=- (2.23)
L -1 o o 1 o o
o -1 o o 1 o
o o -1 o o 1
A continuacin se resuelve las expresiones de las matrices de rigidez no lineal de orden
superior:
(2.24)
Demostracin. De la ecuacin (2.14)
Al sustituir dV en relacin a Adx:

Al integrar en el dominio de la longitud y al reemplazar (2.7a):
=EA
- (6.u 2 + 6.v 2 + 6.w 2) 8(u- u)
22L 3 ~
Al acomodar convenientemente:
Si se expresa en forma de producto escalar:
- ( 6.u2 + 6.v2 + 6.w2)

o
o
( 6.u2 + 6.v2 + 6.w2)
o
o
Al sustituir (2.7a), (2.7b) y (2.7c) en una parte del producto y al multiplicar cada uno de
los trminos:
o
o
o
o
Lo anterior se expresa como producto de una matriz y el vector de desplazamientos u:
!:!..u 6.v 6.w -6.u -6.v -6.w Ui
o o o o o o Vi
=ouT EA o o o o o o Wi
2L 2 -6.u -6.v -6.w !:!..u 6.v 6.w Uj
o o o o o o Vj
o o o o o o Wj
Entonces:
JV OeT ErdV =OuT SI u

donde:
..u ..v ..w -..u -..v -..w

o o o o o o
EA o o o o o o
Sl = 2L2 (2.25)
-..u -..v -..w ..u ..v ..w
o o o o o o
o o o o o o
D
(2.26)
Al cambiar de dV a dx, mediante la relacin dV = Adx:
= EA..u
L 2 L [..u (..u) + ..v (..v) + +..w (..w)] J L dx
Al integrar y reemplazar (2.7a), (2.7b) y (2.7c):
Al multiplicar todos los trminos:
Se puede expresar en forma de producto escalar como:
-..u..v
-..u..w
6.u2
..u..v
..u..w

Al multiplicar y dividir en 2:
-2~u~v
rEA -2~u~w
= u 2L2
2~u 2
2~u~v
2~u~w
Al sustituir convenientemente (2.7a), (2.7b) y (2.7c):
2ui~U- 2uj~U
Ui~V + Vi~U- Uj~V- Vj~U

- rEA Ui~W + Wi~U- Uj~W- Wj~U
- u 2L2
-2Ui~U + 2Uj~U
-~~v-~~u+~~v+~~u
-Ui~W- Wi~U + Uj~W + Wj~U

Es posible expresar como producto de una matriz y un vector de desplazamiento u:
2~u o o -2~u o o Ui
~V ~u o -~V -~u o Vi
= UT 2P
EA
~w o ~u -~w o -..u Wi
-2~u o o 2~u o o Uj
-~v -~u o ~V ..u o Vj
-..w o -~u ..w o ~u wJ
Entonces:
j v r? EedV = uTs 2 u
donde:
2~u o o -2~u o o
~V ~u o -~v -~u o
EA ~w o ~u -~w o -~u
s2 = 2L2. (2.27)
-2~u o o 2~u o o
-~v -~u o ~V ~u o
-~w o -~u ~w o ~u

(2.28)
J 8 TE dV =
v 1J 1J
J ~E (f:l.u2 +
v2 L2
f:l.v2
L2 +
f:l.w2) [!:l.u 8
L2 L
(!:l.u)
L
!:l.v 8 (!:l.v)
+ L L
!:l.w 8 (!:l.w)] dV
+ L L
Al expresar el diferencial de volumen dV en relacin al diferencial de longitud dL como
dV = Adx, se tiene:
Al sustituir (2.7a), (2.7b) y (2.7c):
= 2EA
L 3 (b..u 2 + b..v 2 + b..w 2) [!1u5 (uj- ui) + !1v5 (vj- vi)+ !1w8 (wj- wi)]
Al multiplicar trmino a trmino:
Se expresa como producto escalar:
= 2EA
L 3 (b..u 2 + b..v 2 + b..w 2) [
5ui 8vi 8wi 8uj Vj 8wj ]
Al introducir !1u 2 + !1v 2 + !1w 2 al vector:
+ !1v 2 + !1w 2 )
-b..u (11u 2
-b..v (11u 2 + !1v 2 + !1w 2 )
-5 rEA -b..w (11u 2 + !1v 2 + !1w 2 )

- u 2L3
b..u (!1u 2 + !1v 2 + !1w 2 )
b..v (!1u 2 + !1v 2 + !1w 2 )
b..w (11u 2 + !1v 2 + !1w 2 )

Al acomodar convenientemente:
-b..u (3b..u 2 + b..v 2 + b..w 2 ) - 2b..ub..v 2 - 2b..ub..w 2

-b..v (b..u 2 + 3b..v 2 + b..w 2 ) - 2b..vb..u 2 - 2b..vb..w 2
_ r5 rEA -b..w (b..u 2 + b..v 2 + 3b..w 2 ) - 2b..wb..u2 - 2b..wb..v 2

- u 6V
b..u (3b..u 2 + b..v 2 + b..w 2 ) + 2b..ub..v 2 + 2b..ub..w 2
b..v (b..u 2 + 3b..v 2 + b..w 2 ) + 2b..vb..u2 + 2b..vb..w 2
b..w (b..u 2 + b..v 2 + 3b..w 2 ) + 2b..wb..u 2 + 2b..wb..v 2
Al reemplazar (2.7a), (2.7b) y (2.7c) en algunos casos:
+ ~v + ~w )- 2 (vj- vi) ~u~v- 2 (wj- wi) ~u~w

- (uj- ui) (3~u 2 2 2
- (vj- vi) (~u + 3~v + ~w )- 2 (uj- ui) ~u~v- 2 (wj- wi) ~v~w
2 2 2
- (wj- wi) (~u + ~v + 3~w ) - 2 (uj- ui) ~u~w- 2 (vj- vi) ~v~w
2 2 2
(uj- ui) (3~u + ~v + ~w ) + 2 (vj- vi) ~u~v + 2 (wj- Wi) ~u~w

2 2 2
(vj- Vi) (~u + 3~v + ~w ) + 2 (uj- ui) ~u~v + 2 (wj- wi) ~v~w
2 2 2
(wj- wi) (~u + ~v + 3~w ) + 2 (uj- ui) ~u~w + 2 (vj- vi) ~v~w
2 2 2
Al multiplicar y agrupar segn los desplazamientos:
u; (3b.u 2 + b.v 2 + b.w 2 )

+ v; (2b.ub.v) + w; (2b.ub.w)- Uj (3b.u2 + b.v 2 + b.w 2 ) - Vj (2b.ub.v)- Wj (2b.ub.w)
u; (2b.ub.v) + v; (b.u + 3b.v 2 + b.w 2 ) + w; (2b.vb.w)- Uj (2b.ub.v)- Vj (b.u 2 + 3b.v 2 + b.w 2 ) - Wj (2b.vb.w)
2
u; (2b.ub.w) + v; (2b.vb.w) + w; (b.u 2 + b.v 2 + 3b.w 2 ) - Uj (2b.ub.w)- Vj (2b.vb.w)- Wj (b.u 2 + b.v 2 + 3b.w 2 )
-u; (3b.u2 + b.v 2 + b.w 2 ) - v; (2b.ub.v)- w; (2b.ub.w) + Uj (3b.u2 + b.v 2 + b.w 2 ) + Vj (2b.ub.v) + Wj (2b.ub.w)
-u; (2b.ub.v)- v; (b.u 2 + 3b.v 2 + b.w 2 ) - w; (2b.vb.w) + Uj (2b.ub.v) + Vj (b.u 2 + 3b.v 2 + b.w 2 ) + Wj (2b.vb.w)
-u; (2b.ub.w)- v; (2b.vb.w)- w; (b.u 2 + b.v 2 + 3b.w 2 ) + Uj (2b.ub.w) + Vj (2b.vb.w) + Wj (b.u 2 + b.v 2 + 3b.w 2 )
Es posible expresar en forma de una multiplicacin de una matriz con el vector de despla-
zamiento u.
3Llu2 + Llv2 + Llw2 2Llullv 2Llullw - (3Llu2 + Llv 2 + Llw2 ) - (2Llullv) - (2Llullw) U
2~u..v ..u2 + 3..v2 + .w2 2Llv..w -(2Ll.uLl.v) - (Ll.u2 + 3Llv2 + Llw2 ) - (2Llvllw) ti
2Llullw 2Llvllw Llu2 + Llv 2 + 3Llw2 - (2Llullw) - (2Llvllw) - (Llu 2 + Llv2 + 3Llw2) W
= 6ur ft}
- (3Ll.u2 + Llv2 + Ll.w2 ) - (2Llullv) - (2Llullw) 3Llu2 + Llv2 + Llw 2 2Llullv 2Llullw Uj
- (Mullv) - (Llu 2
+ 3Llv + Llw
2 2
) - (Mvllw) 2Llullv Llu2 + 3Llv2 + Llw2 2Llvllw Vj
- (2Llullw) - (2Llvllw) - (Llu2 + Llv2 + 3Llw2 ) 2Llullw 2..v..w Llu2 + Llv 2 + 3Llw2 Wj
Entonces:
donde:
(2.29)
En la cual:

3~u 2 + ~v 2 + ~w 2 2~u~v
2~u~w 1
h= 2~u~v ~u 2
+ 3~v + ~w
2 2
2~v~w (2.30)
1 2~u~w 2~v~w ~u 2 + ~v + 3~w 2 2
Teorema 2.1.6. Las ecuaciones variacionales (2. 18) y (2. 19} se expresan en forma matricial
como:
(2.31)
(2.32)
Demostracin. Respecto a la ecuaciones (2.18) y (2.19):
Debido a que q se refiere a las componentes de los nodos i y j, entonces se puede

expresar como:
Al integrar en la superficie se obtiene el rea transversal: j s dS = A; adems como 1r qA

representa a la fuerza superficial 1F ' entonces se tiene:
Se puede expresar lo anterior como producto escalar:
1
Fxi
1 F:yt.
1 Fzi
1F.
XJ
1pyj
1p.
ZJ
Entonces:
1R = 8uT 1/

De manera anloga para 2 R, tenemos:
Luego de desarrollar cada trmino de la expresin {2.3), se reemplazan las ecuaciones
{2.20), (2.22), {2.24), (2.26), (2.28), {2.31) y (2.32) en dicha expresin:
Al agrupar convenientemente:
Debido a que la expresin es vlida para todo 8uT, entonces la ecuacin de equilibrio que
relaciona el estado ltimo y actual de la armadura, se puede expresar como:
{2.33)
2.2. ELEMENTO TIPO CABLE
Los cables son elementos estructurales muy flexibles y tienen un comportamiento no li-
neal bajo cargas estticas [1 0]. Mientras ms largos sean los cables, estos presentan una
inherente flexibilidad que involucra grandes desplazamientos y son dominados por una geo-
metra no lineal [19].
Este comportamiento no lineal es a causa de la variacin de su rigidez por efectos del

cambio en su tensin. Los cables slo pueden realizar trabajo en estado de tensin, con toda
la seccin en este estado, no ofrecen resistencia a la flexin y fuerzas de comprensin.
En ciertas estructuras como torres ventadas, los cables son la primera fuente de no linea-
lidad. stos sirven para proporcionar un soporte lateral a la torre y constituyen un elemento
que aporta rigidez al sistema estructural. Debido a la flexibilidad de los cables de soporte, las
torres ventadas son muy flexibles, con desplazamientos laterales relativamente grandes, lo
que hace a la torre ventada susceptible a excitaciones dinmicas del viento y efectos signifi-
cativos de segundo orden, conocidos como efectos P- .6. [24].

Para un adecuado proceso constructivo, los cables son seleccionados de acuerdo a sus
esfuerzos de fluencia. Los cables son normalmente pretensados con el fin de eliminar el
aflojamiento.
2.2.1. Hiptesis del elemento tipo cable
Las consideraciones bsicas a tomar en cuenta son:
La seccin transversal es constante antes y despus de la deformacin del elemento.
El elemento es supuesto como perfectamente flexible, con slo fuerzas de tensin.
La masa del elemento cable se conserva despus de la deformacin, variando el peso
por unidad de longitud despus de cada elongacin del cable.
2.2.2. Formulacin Lagrangiana del elemento tipo cable
Con base en la formulacin Lagr~ngiana, se toma en cuenta el efecto geomtrico no lineal

causado por un gran desplazamiento.
Se considera al cable de una longitud inicial L 0 , suspendido entre dos apoyos fijos I y J,
los cuales tienen coordenadas cartesianas (0, O, O) y (lx, ly, lz) respectivamente. As, las pro-
yecciones horizontales de su longitud son lx y ly, mientras que la proyeccin vertical relativa
entre ambos extremos es lz. Sea el punto P un punto interno del elemento cable. Cuando el
cable se deforma por accin de su peso propio (wL 0 ), este punto P pasa a ocupar una nue-
va posicin descrita por las coordenadas cartesianas (x, y, z) y la coordenada Lagrangiana
p, lo que se define como la longitud del segmento deformado comprendido entre el origen
de coordenadas y el punto sobre la geometra deformada que corresponde al punto P en la
geometra original [20] tal como se muestra en la figura (2.3).

z
Figura 2.3: Coordenadas para la catenaria elstica.

La restriccin geomtrica que debe cumplirse por equilibrio de fuerzas (2.35a), (2.35b) y
(2.35c) es:
(2.34)
En la configuracin inicial de equilibrio, el cable esta sujeto a seis fuerzas nodales F 1 , F 2 ,
F3, F4, F 5 y F5; donde las primeras tres se.localizan en el nodo I tal como se aprecia en la
figura (2.4)
Figura 2.4: Fuerzas en un segmento del perfil deformado de cable.
De la figura (2.4), el equilibrio de las fuerzas para Tes:
T (~;) = -F1 (2.35a)
T (~~) = -F2 (2.35b)
T (~;) = -F3 +ws (2.35c)
La tensin (T) del cable en la coordenada Lagrangiana s, est expresada como:
(2.36)
Para medir la deformacin, se usa el concepto de deformacin de Cauchy o deformacin
ingenieril:
(2.37)
Suponiendo que la longitud de una lnea material en un entorno diferencial es igual a ds
en la configuracin de referencia y a dp en la configuracin actual, tal como lo indica la figura
(2.5) [18].

t=O
(i---;---~-----1H')-------x
ds t=t
(_--- - t - - - " " " 1 ' - - - - ---+()1----- X
dp
Figura 2.5: Lnea de material definida entre dos partculas de un cuerpo en un espacio.
(Fuente: Linero y Garzn, 201 O).
La tensin del cable T, est relacionada con la deformacin e y a la ley de Hooke:
T = EAc = EA ( dp ~ ds) = EA ( ~~ - 1) (2.38)
Donde E es mdulo de elasticidad, A es el rea transversal constante en el perfil no

deformado.
Las condiciones en los extremos de los nodos son:
x = O, y = O, z = O, p = Oen s = O (2.39a)
X = lx, y = ly, z = lz, p = Len S = Lo (2.39b)
Donde L es la longitud del cable en el perfil deformado.
Las coordenadas Lagrangianas s y las coordenadas cartesianas se relacionan como se

indica en las ecuaciones (2.40a), (2.40b) y(2.40c):
x( s) = dx
- ds =
dxdp
- - ds (2.40a}
s ds s dp ds
y( s) = dy
- ds =
1dydp
- - ds (2.40b)
s ds s dp ds
z(s) = dz
-ds =
dzdp
--ds (2.40c)
1 8 ds s dp ds
Teorema 2.2.1. Las coordenadas cartesianas x e y en funcin de la longitud del segmento
deformado son:
y(s) = -:}ts- ~ { ln ( JF[ + Fi + (ws- F3) 2 + ws- F3) -ln ( JF[ + Fi + F}- F3)} (2.42}

Demostracin. Para desarrollar las expresiones de x e y, primero se expresa ~~ y !fls en
funcin de las otras variables de las ecuaciones (2.35a) y (2.38) respectivamente, para luego
reemplazar ambas ecuaciones en (2.40a):
x( s)
{ - F1 (
= Js T
T )
1 + EA ds
f (- F1 + -EA
= Js T
F1)
ds
x(s) = -Fl
- d s - F1-ds
S EA S T
Al reemplazar la ecuacin (2.36) en la expresin anterior e integrar la primera parte:
() -As f A d
x s = EA - Js JF'f + F:j + (ws- F3) 2 s
Si: ws - F3 = uJ F'f + F:j, entonces: wds = duJ F'f + F:j con lo que se obtiene
x(s) = -F1s _ f F1 JF'f + F:j du
EA j Dom(u) wJF'f + F:j + (uJF'f + F:j)2
x(s) = -F1s _ f F1 JF'f + F:j du

EA JDom(u) wJF'f + F:jy'(1 + u2)
x(s) -F1s
=- ---F1
EA W
JDom(u)
1
V1 + u 2
du
Si: u = senM, entonces: du = cosh(}d()
x(s) -F1s
=- ---
EA
F1
w
J Dom(O)
1
J1 + senh 2()
cosh(}d()
De la relacin hiperblica: senh 2 () + 1 = cosh2 (), se obtiene:

x(s) = -F1s _ F1 f d() = -F1s _ F1 ()
EA W } Dom(O) EA W
Al sustituir () en funcin de u:
-F1s F1
x(s) = - - - -senh- 1u
EA w
Como: senh- 1u = ln (u+ (1 + u 2 ) 112 )

-F1s F1 2 1/2
x(s) = EA -:- -:;;ln(u + (1 +u ) )
Al reemplazar u en funcin de s:
x()s =-F1s F1l

---- n ( ws-F3 +
EA w JF'f+F:j

x(s) = -F1s _ F1ln (ws- F3 + JF'f + Fi + (ws- F3)2)

EA w JF'f+Fi
Al integrar la expresin anterior de O a s:
x(s) = -:ls- ~ { ln (V F'f + Fi + (ws- F3) 2 + ws- F3) -ln (V F'f + Fi + Fl- F3)}
De manera anloga para y:
y(s) = -:ls- ~ { ln (V F'f + Fi + (ws- F3) 2 + ws- F3) -ln (V F'f + Fi + Fl- F3)}
D
Teorema 2.2.2. La coordenada cartesiana z en funcin de la longitud del segmento deforma-

do es:
(2.43)
Demostracin. Para desarrollar la expresin de z, se debe expresar ~; y ~~ en funcin de las
variables de las ecuaciones (2.35c) y (2.38) respectivamente, luego reemplazar en (2.40c):
Al sustituir la ecuacin (2.36) en la expresin anterior e integrar la primera parte:
ZS
( )
- F3s
=--+--+
EA
ws
2EA
2 S J F'f
ws - F3
+ Fi + (ws - F3) 2
d
S
Si: ws- F3 = uJF'f + Fi, entonces: wds = duJF'f + Fi, con lo que se obtiene:
-F3s ws2
z(s) = - - + - - +
~ u ( JF'f + Fir
du
EA 2EA Dom(u) wJF 12 + F 22 Vl + u 2
2
- F3s
( ) =--+--+ws J F'f + Fi ~ u d
zs u
EA 2EA W Dom(u) Vl + u 2
2
z(s) = -F3s + ws + JF'f + Fi ~ f d (1 + u
2
)
EA 2EA W 2 } Dom(u) Vl + u2
_ -F3s ws 2 JF'f+Fi. 1 2
z ()
S - EA + 2EA + w V 1+u

Al reemplazar u en funcin de s:
2
_ -F3s ws 2 )F'[ + F:j ws- F3
z (S ) - EA + 2EA + w 1
+( )F[+F:f )
_ -F3s
z (s ) - EA +
ws 2
2EA +w
1 VFr + F2 + (ws- F3)
2 2 2
Si se integra de O a s resulta:
Al incluir las condiciones de contorno (2.39) en (2.41 ), (2.42) y (2.43) se tiene:
lx = - F}J;~" - ~ { ln ( J F[ + F:j + (wL 0 - F3) 2 + wL 0 - F3) - ln ( J F'f + F:j + Fi - F3)} (2.44a)
ly = _F_~o- ~ {zn ( JF[ + F:j + (wLo- F3) 2 + wLo- F3) -ln ( JF[ + F:j + Fi- F3)} (2.44b)
(2.44c)
Como las variables lx, ly y lz son escritas en trminos de las fuerzas del extremo I
(Fr, F 2 , F 3 ), stas pueden expresarse en pequeas variaciones de las mismas mediante una
derivacin de primer orden para obtener una mejor aproximacin [15]:
(2.45a)
{2.45b)
(2.45c)
O expresadas en forma matricial:
(2.46)
Donde F es la matriz de flexibilidad incremental, cuyos elementos fij se indican a conti-
nuacin:
(2.47a)
(2.47b)

(2.47c)
(2.47d)
(2.47e)
(2.47f)
Donde T y Tj son las tensiones del cable en los nodos 1 y J respectivamente. Las
variables H. F2, F3, F4, Fs, F6, Ti y Tj estn relacionadas como:
(2.48a)
(2.48b)
(2.48c}
T = Jp12 + p22 + p32 (2.49a)
Tj = VF42+Fs2+p62 - (2.49b}
La matriz de rigidez del cable se obtiene al invertir la matriz de flexibilidad F:

-1
fn !12 !13
Kc = p- 1 = h1 !22 h3 (2.50)
1 h1 h2 h3 1
La matriz de rigidez tangente y el correspondiente vector de fuerzas internas del elemento
cable, se expresan en trminos de seis grados de libertad [11]:
(2.51)
(2.52)
La matriz de rigidez tangente Kct relaciona el vector de fuerzas internas incrementales al
vector de desplazamientos incrementales nodales a travs de la ley de Hooke:
!:lF1 !:lu
!:lF2 !:lvi
!:lF3 !:lwi
= Kct (2.53}
!:lF4 !:lu.J
!:lFs !:lvj
!:lF6 !:lwj

Luego de encontrar las fuerzas nodales F 1 , F 2 y F3 que satisfacen la tolerancia buscada
de las coordenadas del nodo J, es posible determinar la geometra total del cable al calcular
las coordenadas de puntos a lo largo del cable. Esto es posible debido a que las fuerzas F 1 ,
F 2 y F 3 son conocidas y al reemplazar alguna fraccin de L 0 en las ecuaciones (2.44) se
puede determinar la coordenada comprendida entre los nodos I y J.
Cuando se ejerce alguna fuerza externa, el cable se deforma y la elongacin final del
cable es L =Lo+ ~L, donde ~L se calcula considerando la ley de Hooke:
~L=~T (2.54)
EA
Debido a la hiptesis de conservacin de la masa del cable, el peso unitario por unidad
de longitud (w) en la configuracin deformada se actualiza respecto al peso inicial (w) [1 ].
wL 0
w=-- (2.55)
L
2.3. ANLISIS NO LINEAL ESTTICO
El anlisis comienza con la evaluacin del estado de esfuerzo inicial de los elementos
bajo cargas de peso propio y vivas existentes en la estructura antes de aplicar cargas ccli-
cas o cargas ssmicas. Estas cargas estticas son especificadas como cargas puntuales en
elementos armadura y distribuidas en elementos cable.
La carga esttica debe ser aplicada incrementalmente para obtener la redistribucin de
esfuerzos debido al comportamiento no lineal. Para esto la carga esttica se subdivide en un
razonable nmero de incrementos para encontrar la respuesta no lineal adecuada [26].
2.4. ANLISIS NO LINEAL DINMICO
Se resuelve con el mtodo de Newmark de forma incremental e iterativa. La parte iterativa
es resuelta con el mtodo de la rigidez secante (mtodo de la secante) rigidez tangente
(mtodo de Newton).
El algoritmo de Newmark es frecuentemente utilizado en el anlisis dinmico no lineal de
estructuras.

Las ecuaciones de velocidad y desplazamiento en el tiempo t + flt, se expresan tal como

se indica a continuacin:
(2.56a)
(2.56b)
Donde 'Y y (3, son parmetros que se designan segn la funcin de aceleracin sea cons-
tante o tenga una variacin lineal durante el incremento de tiempo. Para la aceleracin cons-
tante 'Y= ! y (3 = ~;para la aceleracin lineal 'Y= ! y (3 = f; [8].

Al emplear cantidades incrementales y al usar el mtodo de aceleracin constante el cual
es incondicionalmente estable [21], se obtiene la ecuacin incremental del movimiento:
(2.57)
La expresin anterior (2.57) puede ser expresada como:
(2.58)
Donde:
!lPt = flPt + (_i__M

flt + 2e) u+
t 2Mt {2.59)
A 2 4
K = K + flt e + flt2 M (2.60)
Donde M, eyK son las matrices de masa, amortiguamiento y rigidez tangente respec-
tivamente; flPi es el vector de carga externa incremental [3].
Al resolver la ecuacin {2.58) mediante el mtodo de la secante el mtodo de Newton,
se actualizan el desplazamiento, velocidad y aceleracin en el tiempo t + flt como:
{2.61 a)
(2.61 b)
{2.61 e)

2.4.1. Matriz de masa
Para el modelamiento dinmico se opta por una matriz de masas concentradas. Esta ma-
triz considera simplificadamente las propiedades de masa de cualquier elemento. Se supone
que la masa est concentrada en los nodos de los elementos cable o armadura [4]. La matriz
de masas concentradas es una matriz diagonal.
1 o o o o o
o 1 o o o o
M= mL o o 1 o o o
(2.62)
2 o o o 1 o o
o o o o 1 o
o o o o o 1
2.4.2. Matriz de amortiguamiento
Segn el modelo de amortiguamiento de Rayleigh, la matriz de amortiguamiento se define
como:
(2.63)
Donde aM y aK son factores de amortiguamiento proporcionales a la masa y a la rigidez
de la estructura respectivamente. Dependiendo del tipo de amortiguamiento que se suponga
predominante (proporcional a la masa, rigidez ambos) estos factores se calculan con (2.64),
(2.65) (2.66) [26]:
(2.64)
CiM =O; CiK = (2f,i)/wi (2.65)

2
2f,iWiW/ - 2f,jWjWi . 2~jWj - 2~iWi
Ci !VI = Ci K = ____::_:t.----;;,_----=.,,.-- (2.66)
Wj2 - Wi2 , Wj2 - Wi2
Donde ~i y wi son la fraccin de amortiguamiento crtico y la frecuencia circular del modo "i"
respectivamente.
2.5. MTODOS DE SOLUCIN DE LA ECUACIN NO LINEAL
Para resolver el sistema de ecuaciones no lineales 1 = Ku, cuando las cargas 1 son
conocidas, podemos usar los siguientes mtodos:

2.5.1. Mtodo de Newton
Este mtodo es utilizado frecuentemente para encontrar los desplazamientos a partir de
las ecuaciones de equilibrio. La solucin es obtenida como una suma sucesiva de los incre-
mentos de desplazamientos ~u, los que se calculan durante cada ciclo de iteracin.
El esquema iterativo en cada incremento de carga es ilustrado en la figura (2.6).
1--------U------1
Figura 2.6: Mtodo de Newton modificado con incrementos de carga.
En el mtodo de Newton modificado, la matriz tangente Kr se mantiene constante en
cada ciclo de iteracin hasta completar el incremento de carga (~f); a ste se le denomina
como el mtodo de Newton en el presente trabajo.
2.5.2. Mtodo de la secante
Este mtodo es un proceso numrico cuya efectividad se pone de manifiesto en los ca-
ptulos 111 y IV de la presente tesis. En esta metodologa, para cualquier iteracin la carga
incremental es aplicada directamente, para obtener la deformacin total u por iteracin [2]. El
esquema de solucin de la ecuacin, se muestra en la figura (2.7).

fl !1-----;r-:,{
l~~ ------+-~ .6.u ---1

!---- .6. u -----j
u
~------------u------------~
Figura 2.7: Mtodo de la secante con incrementos de carga.
Las fuerzas internas y desplazamientos por cada elemento y nodo respectivamente, son
actualizadas basadas en el ltimo desplazamiento calculado; luego se forma la matriz de
rigidez secante actualizada, en la siguiente iteracin. La matriz secante actualizada y la carga
total de la iteracin son usadas para estimar el siguiente desplazamiento; cuando la diferencia
del ltimo desplazamiento y el desplazamiento calculado previamente es menor a la tolerancia
requerida, se dice que el proceso alcanz la convergencia. Este mtodo tiene la ventaja de
ser sencillo, ya que slo emplea la matriz secante y no requiere de una formulacin compleja
[13].
2.5.3. Mtodo de control de desplazamiento generalizado
Dependiendo de la historia de carga, la rigidez de la estructura puede estar en ablanda-
miento (softening) o endurecimiento (stiffening), la trayectoria de equilibrio puede ser estable
o inestable, y la estructura puede estar en una etapa de carga o descarga. Todos estos fen-
menos son identificados por la ocurrencia de puntos crticos, como los puntos lmites y puntos
de retorno (snap-back points) en la trayectoria carga - deformacin [27], tal como se muestra
en la figura (2.8).
Para trazar esta trayectoria de equilibrio existen diversas tcnicas, las que se basan en
el control de carga (mtodo de la secante, mtodo de Newton) desplazamientos, cada uno
con ventajas y desventajas.

Los principales inconvenientes de los mtodos de secante o Newton se presentan al re-

solver problemas que involucren puntos lmites (punto A de la figura 2.8). Por esta razn, no
son recomendados para trazar la respuesta de la estructura ms all del punto lmite; mien-
tras que en tcnicas de control de desplazamientos, la divergencia puede ocurrir en particular
cuando un punto de control seleccionado muestra un comportamiento de punto de retorno

(punto B de la figura 2.8) [13].
Carga
Punto limite ---....A
Puntos de retorno
E
1
\
1
1
e ~r limite\ !
1
:Endurecimiento:
1 1
1 1
1 1
: : Desplazamiento
1
Estable Inestable. Estable
Figura 2.8: Caracterstica general de sistemas no lineales. (Fuente: Yang, 1994).
Uno de los mtodos numricos ms robustos para un anlisis no lineal, es el mtodo de
control de desplazamiento generalizado propuesto por Yang y Shieh [27], el cual tiene las
siguientes ventajas:
La estabilidad numrica puede ser siempre asegurada en regiones cercanas a puntos
lmites.
La variacin de rigidez de la estructura debido al comportamiento no lineal, puede ser

tomada en cuenta a travs de un parmetro, el que indica el cambio de direccin de
carga en la curva carga - desplazamiento.
Para los fines de esta tesis, los mtodos de la secante y de Newton cumplen con los
requisitos de convergencia con relacin al comportamiento no lineal (elstico, geomtrico),
por lo que se consideran suficientes para resolver las aplicaciones que se presentan en el
captulo IV.

FACULTAD DE INGENIERA CIVIL CAPTULO 111: ALGORITMO DE SOLUCIN
CAPTULO 111
ALGORITMO DE SOLUCIN
3.1. MATRIZ DEL ELEMENTO TIPO ARMADURA
La ecuacin de equilibrio de la armadura esta expresada como (2.33):
( ke + kg + SI + S2 + S3) U+ I j =
2
j
Donde:
ke es la matriz de rigidez elstica (2.21):
1 o o -1 o o
o o o o o o
k_ EA o o o o o o
e- L
-1 o o 1 o o
o o o o o o
o o o o o o
k 9 es la matriz de rigidez geomtrica (2.23):
1 o o -1 o o
o 1 o o -1 o
F o o 1 o o -1
kg = L
-1 o o 1 o o
o -1 o o 1 o
o o -1 o o 1
y las matrices si (2.25), s 2 (2.27) y s 3 (2.29) son de orden superior:

/j.u /j.v /j.w -/j.u -/j.v -/j.w
o o o o o o
EA o o o o o o
SI= 2L2
-/j.u -/j.v -/j.w /j.u /j.v /j.w
o o o o o o
o o o o o o
2/j.u o o -2/j.u o o
/j.v /j.u o -/j.v -/j.u o
EA /j.w o /j.u -/j.w o -/j.u
s2 = 2L2
-2/j.u o o 2/j.u o o
-/j.v -/j.u o /j.v /j.u o
-/j.w o -/j.u /j.w o /j.u
s3 = 6L3
EA
[:h ~h l
Donde:
3~u2 + ~v2 + ~w2

2~v~w
[ 2/j.u/j.v 2/j.u/j.w
h = 2/j.u/j.v /j.u 2 + 3fj.v 2 + /j.w 2

1
2/j.u/j.w 2/j.v/j.w /j.u 2 + /j.v 2 + 3fj.w2
3.2. MATRIZ DEL ELEMENTO TIPO CABLE
3.2.1. Algoritmo del elemento tipo cable
Para determinar la matriz de rigidez tangente Kct. se debe determinar primero las fuerzas
F 1 , .F2 y F 3 del nodo 1 [16]. Para ello, se comienza con una estimacin inicial de estas fuerzas
nodales mediante [11]:

Fl = - wlx (3.1 a)
2-X
F2 =- wly (3.1 b)
2-X
P.3 = w
2
(-z z
cosh.X
senh.X +
L0 ) (3.1 e)

Donde:
2
106 si Ux + ly 2 ) = O
2
A= 0,2 si Lo :s; lx 2 + ly 2 + lz 2 (3.2)
3 ( Lo2_z2 - 1) si Lo
2
> lx 2 + ly 2 + lz 2
lx2+ly2
Al reemplazar las fuerzas estimadas en (2.44) se determina las coordenadas Uxo, lyo, lzo)
del cable, y respecto a la coordenada indicada del nodo J (lx, ly, lz), se determina el vec-
tor de diferencia .L = { (lxo - lx) (lyo - ly) (lzo - lz) } T el cual debe ser menor a la
tolerancia buscada.
Si la tolerancia no se cumple, las fuerzas de aproximacin del nodo I se corrigen usando
el vector de diferencia como:
.F1 .lx
( .F2
.F3
)=K, .ly
.lz
) {3.3)
Luego se actualizan las fuerzas:
F1 F1 .F1
( )=( )+( )
F2
F3
F2
F3
.F2
.F3
(3.4)
donde Kc es la matriz de rigidez de la ecuacin (2.50).
A continuacin se presenta el algoritmo empleado para encontrar la matriz de rigidez de
un elemento cable.
Algoritmo 3.2.1. Cable con La conocido
Paso 1: Los datos iniciales son las coordenadas de los nodos I(xi,Yi,zi) y J(xj,yj,Zj), las
propiedades fsicas (E, A, m, La) y la tolerancia (tol).
Paso 2: Calcular las coordenadas relativas entre los nodos I y J: lxo = Xj - Xi, lyo = Yi - Yi y
Paso 3: Dar una primera aproximacin de las fuerzas F 1 , F 2 y F 3 del nudo I, al utilizar la
ecuacin (3. 1).
Paso 4: Actualizar la coordenada relativa del nodo J (lx, ly, lz), al reemplazar F1 , F2 y F3 en la
ecuacin (2.44).

Paso 5: Calcular el vector diferencia:
Paso 6: Si II.6.LIIoo < tol, calcular Kct y Fint con (2.51) y {2.52) respectivamente. Caso contrario
continuar los siguientes pasos.
Paso 7: Calcular el vector de correccin de las fuerzas del nodo I: .6.F = Kc.6.L (3.3).
Paso 8: Actualizar las fuerzas de la siguiente iteracin empleando {3.4) y retornar al Paso 4.
INICIO
~
Datos iniciales
w,E,A,Lo
I(x;,y,z) J(:t:j,Yj,ZJ)1
lxo
lyo
Calcular
= :t:j- :r
= 1/J- Jli
lzo = .!J- z;
Aproximar
F1, F2 y F3
--"'
Ir
Calcular
lx, l11 , lz
D.L

Calcular el vector diferencia
={(t~,o - lx) (lyo - ly) Uzo -U} r
4
NO SI
Calcular el vector de
correccin: D.F = KcD.L

Calcular Kct y fuerzas
nodales del elemento:
Ft.F2.F.~.F4,F5

y Fa
j_
Actualizar e FIN
rl rl rl
F2
Fa
= F2
F3
+ D.F2
D.F3
Figura 3.1: Diagrama de flujo del elemento cable para Lo.

Cuando la tensin inicial del cable T 0 en el nodo I es conocida, se emplea el siguiente
proceso de iteracin para determinar la longitud del cable no esforzado.
Algoritmo 3.2.2. Cable con T 0 conocido
Paso 1: Como datos iniciales se necesitan las propiedades fsicas w, E, A, la tensin en el nodo
I (I'0 ) y la coordenada de los nodos I ( Xi, Yi, zi) y J (Xj, Yi, Zj).
Paso 2: Calcular las coordenadas relativas entre los nodos I y J: lxo = Xj - xi, lyo = Yi - Yi y
lzo = Zj- Zi
Paso 3: Inicializar la longitud no esforzada L 0 y las fuerzas del nodo I (F1 , F 2 , F 3 ) como sigue:
Lo~ Vlxo 2 + lyo 2 + lzo 2 (3.5a)
Fl = _lxoTo (3.5b)
Lo
F2 = _lyoTo (3.5c)
Lo
F3 = _ZzoTo (3.5d)
Lo
Paso 4: Actualizar la coordenada relativa del nodo J (lx, ly, lz), al reemplazar F 1 , F 2 , F 3 y L 0
en la ecuacin (2.44).
Paso 5: Calcular el vector diferencia: b.L ={ (lxo -lx) (lyo -ly) (lzo -lz) } T, la diferen-
cia entre la tensin conocida T 0 y la estimada Ti (2.49a): b.T =T 0 - Ti
Paso 6: Si /lb.LIIoo < tol y llb.TIIoo < tol, calcular Kct y Fint con (2.51} y (2.52) respectivamen-
te. Caso contrario continuar con los siguientes pasos.
Paso 7: Calcular Ce por diferenciacin de las expresiones (2.44) y (2.49a} con respecto a F 1 , F 2 , F 3
y Lo como:
Blx Blx Blx Blx

BFr BF2 BF3 BLo
.0:J_ .0:J_ .0:J_ .E!:JL
BFr BF2 BF3 BLo
Ce= (3.6)
2h 2h 2h Blz
BFr BF2 BF3 BLo
ari ari ff!l ari
BFr BF2 BF3 BLo

Paso 8: Calcular el vector de correccin:
b..F1 b..lx
b..F2 b..ly
=Ce {3.7)
b..F3 b..lz
b..Lo b..Ti
Paso 9: Actualizar las fuerzas: F 1 = F 1 + b..F1, F 2 = F 2 + b..F2, F 3 = F 3 + b..F3, la longitud:

Lo = Lo + b..L y retornar al Paso 4.
0
Calcular
lxo = "'J- :;
lytl = YJ- .IJ;
f::o = Zj- Z
Calcular el vector diferencia

C.L ={U.ro- lx) (1 11o- 111 ) Uzo -I=J} T
Calcular
C.F}
C.F2 = Ce
{C.Ix}
C./ u
{C.F3 C./z
C.Lu C.T;
F Actualizar:
F + C.F
=
F2 = F2 + C.F2
F3 = F3 + C.F3
Lo= Lo+ C.Lo
Figura 3.2: Diagrama de flujo del elemento cable para To.

3.2.2. Descripcin del algoritmo
Con los algoritmos descritos en el acpite 3.2.1, es posible determinar la configuracin
del cable apoyado en los nodos I y J, y cuya distancia relativa es (lxo, lyo, lzo), para ello se
deben de realizar iteraciones de convergencia al nodo J tal como se ilustra en la figura (3.3).
Figura 3.3: Configuracin en un paso de la iteracin.
Luego de satisfecha la convergencia al nodo J, se pueden calcular las fuerzas (F4 , F 5 , F 6 )
en dicho nodo usando las expresiones de equilibrio (2.48}:
:r:
Figura 3.4: Elemento cable luego de la convergencia.
Las figuras (3.3} y (3.4) slo muestran un elemento cable de extremos I y J, el cual no se

puede utilizar cuando una carga externa se aplica en una parte del cable, tal como se apre-
ciar en la aplicacin 2 del captulo 4 de la presente tesis. Para ello, se divide internamente
al cable segn un cierto nmero de segmentos (N s) respecto a una longitud conocida,

empleando las expresiones (2.41 ), (2.42) y (2.43).
Figura 3.5: Cable con nodos internos.
En cada segmento es posible determinar las fuerzas en los nodos adicionados, mediante
condiciones de equilibrio. Estas fuerzas nodales y la longitud de cada segmento son emplea-
das para determinar la matriz de rigidez de cada segmento cable.
Figura 3.6: Fuerzas de cada elemento.
La matriz de cada segmento cable se agrupa en una matriz global del cable, la cual es
usada en el anlisis no lineal.

3.3. ESTRATEGIA DE SOLUCIN
Como estrategia de solucin primero se construye la matriz de rigidez de los cables; a tra-
vs de la formulacin numrica descrita en el acpite 2.2.2., y con la formulacin variacional
expuesta en el acpite 2.1.1. se construye la matriz de rigidez de los elementos armadura.
Se combinan ambos tipos de elementos, el cable tensado por su peso propio o con una ten-
sin inicial en estado no lineal y la armadura en estado lineal. La matriz de rigidez global se
ensambla a partir de las matrices indicadas en los acpites 3.1 y 3.2. El sistema de ecuacio-
nes ensambladas es resuelto con ambos elementos de acuerdo al mtodo de la secante o al
mtodo de Newton en un proceso incremental de cargas.
3.3.1. Formato de entrada del modelo
Las siguientes figuras muestran el formato de entrada para los programas del mtodo de
la secante y del mtodo de Newton, los cuales se explican en el acpite 3.3.2.
La tabla (3.1) muestra el formato de entrada de la geomtria, donde la columna 1 (Nodo)
enumera los nodos de la estructura, las columnas 2-5 (x, y, z) son las coordenadas de los
nodos y finalmente la columna 6 (Jm) es la masa concentrada del nodo. Cabe indicar que la
masa concentrada slo es necesaria para el anlisis dinmico.
Tabla 3.1: Formato de entrada de geometra y masa concentrada.
Nodo x y. z Jm
n1 x1 y1 z1 Jm1
n2 x2 y2 z2 Jm2
n3 x3 y3 z3 Jm3
nn xn yn zn Jmn
La tabla (3.2) muestra el formato de propiedades de los elementos, en donde la columna
1 (Elemento) enumera los elementos que componen la estructura, las columnas 2 y 3 (Nodo)
asignan el nodo de comienzo (ni) y el nodo final (nj) que identifica al elemento, la columna
4 (E) es el mdulo de Elasticidad, la columna 5 (A) es el rea transversal, la columna 6
(m) es el peso por unidad de longitud, la columna 7 (telem) identifica a qu_e tipo de elemento

corresponde (1: Armadura 2: Cable), la columna 8 (N s) en el caso del elemento cable divide
a ste y finalmente la columna 9 (To Lo) es la tensin inicial T 0 del nodo i la longitud no
deformada L 0 del elemento cable.
Tabla 3.2: Formato de entrada de propiedades.
Elemento Nodo E A m telem Ns To Lo
e1 n1 n2 E1 A1 m1 Ns1 Ton1 Lon1
en ni nj En An mn Nsn Tonn Lonn
La tabla (3.3) muestra las restricciones y fuerzas externas aplicadas en el nodo, donde la
columna 1 (Nodo) enumera los nodos de la estructura, las columnas 2-4 (Restricciones) son
las restricciones al movimiento del nodo en la direcciones x, y y/o z tan slo colocando (1:
Restriccin 0: Libertad) y por ltimo las columnas 5-7 (Fuerzas) son las fuerzas externas
aplicadas al nodo en las direcciones x, y y/o z.
Tabla 3.3: Formato de entrada de restricciones y fuerzas nodales.
Nodo Restricciones Fuerzas

n1 R1x R1y R1z F1x F1y F1z
nn Rnx Rny Rnz Fnx Fny Fnz
3.3.2. Algoritmo para el anlisis esttico no lineal
Algoritmo 3.3.1. Mtodo de la secante
Paso 1: Como datos iniciales se necesita la geometra del modelo (x, y, z), propiedades fsicas
(E,A,m), tipo de elemento (telem). Si es elemento cable (telem = 2) se indica el
nmero de segmentos (N s), tensin longitud inicial del cable (To Lo). Adems las

restricciones, vector de fuerzas externas(! e), nmero de incrementos de carga (Ninc)
y la tolerancia (tal).
Paso 2: Para el elemento cable se determina su configuracin de equilibrio de acuerdo a su

longitud Lo tensin iniciaiTo, y se calculan las fuerzas internas nodales de cada seg-
mento (Fint 1 ) (2.52) segn como se ha subdividido el cable (esto se puede resolver
como se indica en el acpite 3.2. 1.). Luego se actualizan los datos iniciales de geome-
tra, propiedades, restricciones y fuerza aplicada en los nodos debido a la segmentacin
del cable.
Paso 3: Condiciones iniciales: desplazamientos acumulados u 1 = {0}, incrementos de des-
plazamientos ..u 1 = {O}, fuerzas axiales N 1 = {0}, incrementos de fuerzas axiales
..N 1 = {0}, vector de fuerzas internas /i = {O} y vector de carga incrementada
fe= {0}.
Paso 4: El incremento de carga es i = 1, 2, 3 ... Ninc.
Paso 5: Para la primera iteracin (j = 1) en el incremento de carga i:
5. 1: La carga aplicada incremental es: 1e = 1e + li~e.

5.2: El vector de carga actuante: R = fe - f i.
5.3: Formar la matriz de rigidez global Kj i de la estructura:
5.3.1: Para el elemento armadura emplear las matrices de elasticidad ke (2.21 ), geo-
mtrica k9 (2.23) y las de orden superior s 1 (2.25), s 2 (2.27) y s 3 (2.29). Para
su ensamble usar las posiciones (xi, Yi, zi), propiedades fsicas, esfuerzos
normales Ni y las deformaciones ui.
5.3.2: Para el elemento cable emplear la matriz Ket (2.51) y usar para su ensam-
ble las propiedades fsicas, longitud inicial de cada elemento Loi y fuerzas
internas Finti.
Paso 6: Realizar el siguiente proceso iterativo:
6.1: Resolver el incremento de desplazamientos ..u/ de la ecuacin K/ ..u/ = R.
6.2: Actualizar temporalmente las posiciones de los nodos, deformaciones, fuerzas
normales y fuerzas internas en los elementos:
6.2.1: Actualizar los desplazamientos: ui+ 1 = ui + ..u/

6.2.2: Actualizar posicin de los nodos de (xi, Yi, zi) a (Xj, Yj, Zj)-
6.2.3: Calcular tlNj con respecto a la deformacin llu/.
6.2.4: Actualizar las fuerzas axiales: Ni+ 1 =Ni+ tlNj.
6.2.5: Para los elementos cable, calcular la variacin de fuerzas internas de cada
elemento al multiplicar cada rigidez con su respectiva deformacin y almace-
nar en tlFintj, actualizar la longitud del cable Loi+l debido a su elongacin
y adems por conservacin de la masa actualizar la densidad lineal del ele-
mento a mi+ 1
6.2.6: Actualizar las fuerzas internas de los elementos cable: Finti+ 1 = Finti +
tlFintj
6.3: Ensamblar la matriz de rigidez global K/:
6.3.1: Para el elemento armadura, respecto a las posiciones usar (xi, Yi, zi) en el
ensamble de la matriz geomtrica ke y (Xj, Yj, Zj) en las dems matrices;
adems las propiedades fsicas, esfuerzos normales Ni+ 1 y las deformacio-

nesui+1.
6.3.2: Para el elemento cable emplear la matriz Ket y usar para su ensamble las
propiedades fsicas incluido mi+ 1 , longitud Loi+l y fuerzas internas Finti+l
6.4: Calcular la fuerza residual de la estructura: f s = K/ llu/.
Paso 7: Verificar la convergencia: Si [[R- f s [[ 00 < tol ir al paso 8, en caso contrario j = j +1

y repetir paso 6.
Paso 8: Almacenar la posicin alcanzada: (xi+l, Yi+ 1, zi+ 1) = (Xj, Yj, Zj), la fuerza interna
alcanzada hasta este incremento de carga: fi = fi + f s
Paso 9: Si i ~ N in e, entonces ir al paso 4 con el siguiente incremento i i + 1, pero si
i = N inc detener el proceso.
Algoritmo 3.3.2. Mtodo de Newton
Los pasos 1, 2, 3 y 4 del mtodo de la secante se repiten. A continuacin se describen los

pasos que siguen:
Paso 5: Para la primera iteracin (j = 1) en el incremento de carga i:
5.1: La carga aplicada incremental es: fe = fe + rh~e.

5.2: Formar la matriz de rigidez global K/ de la estructura:
5.2.1: Para el elemento armadura emplear las matrices de rigidez elstica ke (2.21 },
geomtrica k 9 (2.23} y las de orden superior s 1 (2.25), s 2 (2.27) y s 3 (2.29).
Para su ensamble usar las posiciones (xi, Yi, zi), propiedades fsicas, esfuer-
zos normales Ni y las deformaciones ui.
5.2.2: Para el elemento cable emplear la matriz Kct (2.51) y usar para su ensam-
ble las propiedades fsicas, longitud inicial de cada elemento Loi y fuerzas
internas Finti.
Paso 6: Inicializar las condiciones de Newton: .6.fsj = {0}, !sj = fi, vector de fuerza residual
Rj =fe- fi, vector de desplazamiento total por incremento ui+ 1 = ui, el vector de
fuerzas axiales Ni+ 1 = Ni y las fuerzas internas Finti+I = Finti. Adems hacer
constante la matriz tangente: Kr = K/
Paso 7: Realizar el siguiente proceso iterativo:
7.1: Resolver el incremento de desplazamientos .6.uji de la ecuacin Kr.6.u/ = Rj.

7.2: Actualizar y acumular posicin de los nodos:
7.2.1: Acumular deformacin: ui+ 1 = ui+ 1 + .6.u/.

7.2.2: Actualizar posicin de los nodos: si j = 1 actualizar de (xi, Yi, zi) a (Xj, Yj, Zj),
cuando j 2:: 2 actualizar de (Xj- 1 , 1J- 1 , Zj- 1 ) a (Xj, Yj, Zj)
7.2.3: Calcular .6.Ni con respecto a la deformacin .6.uj i.
7.2.4: Acumular los esfuerzos normales: Ni+ 1 = Ni+ 1 + .6.Ni.

7.2.5: Para los elementos cable, calcular la variacin de fuerzas internas de cada
elemento al multiplicar cada rigidez con su respectiva deformacin y almace-
nar en .6.Fintj, actualizar la longitud del cable Loi+I debido a su elongacin
y adems por conservacin de la masa actualizar la densidad lineal del ele-

mento a mi+1
7.2.6: Actualizar las fuerzas internas de los elementos cable: Finti+ 1 = Finti+ 1 +
.6.Fintj.
7.3: Ensamblar la matriz de rigidez global K/:
7.3.1: Para el elemento armadura, respecto a las posiciones usar (xi, Yi, zi) en el
ensamble de la matriz geomtrica ke y (Xj, Yj, Zj) en las dems matrices;

adems las propiedades fsicas, esfuerzos normales Ni+ 1 y las deformacio-

nesui+1.
7.3.2: Para el elemento cable empleqr la matriz Kct y usar para su ensamble las
propiedades fsicas incluido mi+ l. longitud Loi+ 1 y fuerzas internas Finti+l
7.4: Calcular la fuerza ejercida por la respuesta de la estructura: b.f sj = K/ D.u/.
7.5: Actualizar el vector fuerza residual f s: f sj+I = b.f sj + f sj
7.6: Actualizar la fuerza actuante: Rj+l = Rj - b.f sj
Paso 8: Verificar la convergencia: Si IIRj+III 00 < tol ir al paso 9, en caso contrario j = j +1y
repetir paso 7.
Paso 9: Almacenar la posicin alcanzada (xi+I, Yi+l, Zi+I) (Xj, Yj, Zj) y la fuerza interna
alcanzada fi = !sj+I
Paso 10: Si i :S N inc, entonces ir al paso 4 con i = i + 1, pero si i = Ninc detener el proceso.

( INICIO)
1 Geometra, propiedades,
restricciones y fuerzas
j (j)
"

Determinar equilibrio del cable
y actualizar geometra, propiedades,

(j)
restricciones y fuerzas.
2/
,1

Condiciones iniciales:
u 1 = {0}, ..u 1 = {0},
3
N 1 = {0}, M 1 = {0},
/; = {O} Y fe= {0} . 1
..
1
--- i = 1, 2, 3 ... Ninc )' ;,
1 r=l
2
l?'Z
1/c:fc+~
1
1
1 R-Jc-/i
1
1
1
kc, k9 , St,

Matriz armadura global:
82 1 S3, 4
3;
1
.,
Emplear (x;, y;,z;), E,A, Ni y u.i.
3 @
1
1
Matriz cable global:
](el JI" lJ
~,t
1 Emplear E, m, Lo; y Fint;.
1
1 1 Matriz

de rigidez global
ambos elementos: K/
1
_1
-..
1 K/ ..u/ = R.!
+
Actualizar temporalmente:
Deformaciones: ui+t =u;+ ..u/
Posiciones: de (x;, y;, z;) a (X;, Yj, Zi)
Armadura: Ni+ 1 =Ni + ..Ni.
Cable: ..Finti, Lo;+t mi+t
Finti+l = Fint; + ..Finti .

kc emplear (x;, y;, z;)
k9 , s, s2, s3 emplear (X;, Yj, Zi) 4
adems de E,A, Ni+I y ui+I

I<ct .Lu; ...L...-b.1lj i ---J. u
Emplear E, A, mi+, Lo;+t Y Finti+t J._.u;+1___J.
1 Matriz

de rigidez global
ambos elementos: K/
1
+
1 fs = K/..u/ 1
1 j_
j=j+l
1 iR-fslloc<to
1
1"
1
(x;+, Yi+J.Z+J) = (Xj, lj, Zj) 1
1 /;=/;+/.
L ------~
FIN
Figura 3.7: Diagrama de flujo, mtodo de la secante.

3
m
CD
CD
u 1 = {0}, ..u 1 = {0},
JI
N 1 = {0}, M 1 = {O}, 3
/; = {O} Y Je= {0}. m

JI
3
4;
L "
Condiciones de Newton:
..Jsj = {0}, fsj = /;.
Rj =fe-/;.
ui+l =ui,
Ni+t =N;,
Fint;+l = Fint,
Ky =K/
Actualizar v acumular:
Deformacin: ui+l = ui+l +a u/.
Posicin: sij = 1, de (:r;,y;,z;) a (Xj,Yj,Zj)
si j ;?: 2, de (Xi_,l'j_,Zi_) a (Xj,l'J,Zi)
Armadura: Ni+l = Ni+ 1 +Mi.
Cable: ..Finti, Lo;+! m;+!
Fint;+! = Finti+l + D.Finti
k. emplear (:r;, y,z;),
J.,9 , s1 , s2 , s3 emplear (Xi, l'j, Zi),
adems de E,A, Ni+l y ui+l.
-@
Figura 3.8: Diagrama de flujo, mtodo de Newton.

3.3.3. Algoritmo para el anlisis dinmico no lineal
El siguiente algoritmo es la continuacin del anlisis esttico, del cual se obtiene la matriz
de rigidezglobal (a la que se denominar Ka), geometra deformada (xi+ 1 , Yi+l, Zi+l), pro-
piedades fsicas (E, A, mi+l, Loi+l), fuerzas externas (/e) y fuerzas internas (Ji) despus
de un incremento de cargas i = Ninc. Esta geometra, propiedades y esfuerzos finales son
considerados datos de entrada del proceso dinmico. Adems se debe ensamblar la matriz
de masa global Ma (2.62) que incluya las masas concentradas Jm.
El siguiente algoritmo de anlisis dinmico se basa en el fundamento terico desarrollado
.en el acpite 2.4.
Algoritmo 3.3.3. Anlisis dinmico no lineal
Paso 1: Los datos que se adicionan son: Registro ssmico (aceleracin 9 ), cantidad de datos
del registro (N acel), direccin del sismo (f d), intervalo de tiempo del registro (flt) y
desplazamiento inicial (u 0 ).
Paso 2: Calcular el vector de masa f m que depende de la direccin de anlisis del sismo.
Paso 3: Determinar los factores de amortiguamiento proporcionales (a. K y O'. M) segn el amor-
tiguamiento considerado, y ensamblar la matriz de amortiguamiento viscoso global Ca
(2.63).
Paso 4: Condiciones iniciales: desplazamientos u 1 u 0 , velocidad 1 {O} y aceleracin

1 = Ma- 1 (/mfl 91 - Ca- fi +fe)-
Paso 5: Calcular: a = 1tMa + 2Ca y b = 2Ma.
Paso 6: El proceso se realiza hasta completar la cantidad de datos de aceleracin ssmica de
cada flt para i = 1, 2, 3 .. .N acel - l.
Paso 7: Para j = 1 de la aceleracin i:
7.1: Calcular fl.Pi = fm(fl 9 i) + aiJ.i + bi.

7.2: Ensamblar la matriz de rigidez global K/ de la estructura:
7.2.1: Para elemento armadura emplear las matrices de rigidez elstica ke (2.21) y
geomtrica k 9 (2.23). Para su ensamble usar las posiciones (xi, Yi, zi), pro-
piedades fsicas y esfuerzos normales Ni.

7.2.2: Para elemento cable emplear la matriz Kct (2.51) y usar para su ensamble las
propiedades fsicas, longitud inicial de cada elemento Loi y fuerzas internas

Finti.
7.3: Formar la matriz de rigidez: k= K/+ ltCc + bMc.

7.4: Resolver b..u/ de b..Pi = k b..u/ mediante cualquier mtodo como el de la se-
cante el de Newton.
7.5: Actualizar:
7.5.1: Desplazamiento:ui+l = ui + b..u{

7.5.2: Velocidad:i+ 1 = -i + ltb..ui.
7.5.3: Aceleracin:i+ 1 = -i- f.ti + bD..ui.
Paso 8: Si i ~ N acel - 1, entonces ir al paso 6 con i = i + 1, pero si i = N inc- 1 detener el
proceso.
Estos algoritmos se implementan en Matlab.

INICIO
Registro sismico (aceleracin 9 ),

cantidad de datos del registro (N acel),
direccin del sismo (/ d), intervalo de registro (~t) y
desplazamiento inicial (u 0 ).
Desplazamientos u 1 = u0 , Velocidad - 1 = {O}
Aceleracin 1 = Mc- 1 Urn~91 - Cc- /;+fe)
,---
k. y kg
Emplear (x;, y;, z;), E, A y Ni

Kct
Emplear E, m;, Lo; y Fint;
1
1
1 Actualizar:
1 Desplazamiento: ui+ 1 =u;+ ~u/
1
Velocidad: i+t = -; + &~ui
1
Aceleracin i+l = -ii - :l.ii
At
+~
4 ~u;
L ______ _
Figura 3.9: Diagrama de flujo de anlisis dinmico no lineal.

FACULTAD DE INGENIERA CIVIL CAPTULO IV: APLICACIONES
CAPTULO IV
APLICACIONES
4.1. ANLISIS ESTTICO
4.1.1. Cable suspendido a diferente nivel
La aplicacin consiste en determinar la tensin del cable y representarla grficamente
junto con la geometra no deformada, como es planteada por M. Paliares y W. Rodrguez

[20]. Los datos del problema se encuentran en la tabla (4.1 ).
Tabla 4.1: Informacin para el problema.
ltem Dato
w: Peso propio del cable 0,85 kN/m

A: rea de la seccin transversal 2xlo- 4 cm2
E: Mdulo de elasticidad 1,5x107 kN/m 2

Lo: Longitud del cable 28 m
lx: Distancia horizontal entre apoyos 20 m
lz: Desnivel entre apoyos 8,5 m
El clculo se realiza mediante un elemento cable el cual converge al punto J con una
tolerancia de w- 10 . Luego se segmenta internamente en 14 partes en las cuales se determina
la tensin de los nodos y sus coordenadas, los resultados se muestran en la tabla (4.2} y la
figura (4.1 ).

Tabla 4.2: Comparacin de resultados.
Paliares y Rodrguez Presente trabajo
S x(s) z(s) Tensin S x(s) z(s) Tensin
(m) (m) (m) (kN/m) (m) (m) (m) (kN/m)
o 0,000 0,000 17,172 o 0,000 0,000 17,172

2 0,765 -1 ,860 15,600 2 0,765 -1 ,860 15,600
4 1,610 -3,683 14,058 4 1,610 -3,683 14,058
6 2,552 -5,457 12,557 6 2,552 -5,457 12,557
8 3,610 -7,163 11,112 8 3,610 -7,163 11,112
10 4,811 -8,770 9,751 10 4,811 -8,770 9,751
12 6,184 -10,231 8,513 12 6,184 -10,231 8,513
14 7,754 -11,475 7,459 14 7,754 -11,475 7,459
16 9,529 -12,397 6,676 16 9,529 -12,397 6,676
18 11,469 -12,878 6,268 18 11,469 -12,878 6,268
20 13,467 -12,831 6,308 20 13,467 -12,831 6,308
22 15,384 -12,266 6,788 22 15,384 -12,266 6,788
24 17,125 -11,279 7,625 24 17,125 -11,279 7,625
26 18,660 -9,991 8,716 26 18,660 -9,991 8,716
28 20,000 -8,500 9,980 28 20,000 -8,500 9,980
l-e- Geometrla del cable calculado --&- Geometrla del cable -+Tensin en el cable calculado .. ~ .. Tensin en el cable 1
o
Tr :: 2
. .
Qi . : -- ~:::..:.:.~:.:.::::;:~:::::..:::..:::.:::=4:_:~:L.:.::~ ~~l._ -...~. :.~. e.'.=.. 6

~ 8 . -~>7-~~r~..... [: .......... ; ~ ~-~ 8 ~
..,
Cll
,/j ~ '~., .a
~ 10 ............... - ....... ... .. .. ~ ..._....~;r:
:,..-
-~ = -~ .....~:> 10 ~
lii ..
!a 12 ................. ;.......
: . .>/
/'...~:\:...................[. 12
=/
14 ;.,/ ............. <:~ ................, .................. ,.................. ;... ................ . 14
/ .
fo" .. ~- .................. ~ .... ......................... .
16 / .. .. 16
1 ~~----~2----~4~----~6----~B~----~10~--~1~2----~~--~1~6----~1~8~--~2~0 18
Distancia hozontal entre apoyos (m)
Figura 4.1: Comparacin numrica de aplicacin con apoyos situados a distinto nivel.

De la figura (4.1) y la tabla (4.2), los valores de las coordenadas x y z as como la tensin
en los nodos del cable dividido en 14 elementos, se ajustan a lo calculado por M. Paliares y
W. Rodrguez [20].
Para estimar la ubicacin y la tensin en el punto ms bajo, se divide al cable en 1000
elementos cable, y se busca la menor coordenada z correspondiente a la tensin mnima, la
cual es igual en magnitud a la componente horizontal F 1 del apoyo l. Estos resultados se
aprecian en la figura (4.2) y se listan en la tabla (4.3).
1- Geometrla del cable calculado -Tensin en el cable calculado 1

~--~----;===~====~==~====~==~==~~~~--~o
............. T ............ 2
---~---
............... ~- ................ ~-.-........... . ... .................. ! . . . . . . . . . . .. -~ .................. ~-
............ ;............. -6
"18o~--~-----':-------:-------:8,___ _1:'::0-----:'12::----:':----:'::-------:'::----:2::':!0 "18
Distancia horizontal entre apoyos (m)
Figura 4.2: Coordenada y tensin mnima en el punto ms bajo.
Tabla 4.3: Coordenada, tensin en el punto ms bajo y fuerzas en el apoyo l.
Investigador X z Tensin F1 Fs
(m) (m) (kN) (kN) (kN)
M. Paliares y W. Rodrguez [20] 12,269 -12,925 6,229 -6,299 16,003
Presente trabajo 12,285 -12,925 6,229 -6,229 16,003

4.1.2. Cable suspendido con carga puntual
La siguiente aplicacin es tomada como referencia para validar diferentes mtodos para
simular cables. Propuesta inicialmente por Michalos y Birnstiel, luego analizada por O'Brien
y Francis (como lo indica Hu u [11 ]), Jayaraman y Knudson [15], Tibert [25], Andreu [1 ], Hu u y
Seung [11].
El problema consiste en determinar el desplazamiento del nodo 2, cuando al cable sus-
pendido en condiciones de propio peso, se le aplica una carga concentrada. La configuracin
inicial y su informacin estn en la figura (4.3) y la tabla (4.4).
r ~ 35,586 kN
~------------------------------- 3
30,48 m
121,92 m
1524m
304,8m
Figura 4.3: Configuracin inicial de su geometra bajo peso propio.
Tabla 4.4: Informacin para el cable suspendido.
ltem Dato
w: Peso propio del cable 46,1167 N/m

A: rea de la seccin transversal 5,4839 cm 2
E: Mdulo de elasticidad 13100039 Nfcm2
Lo1-2: Longitud inicial del cable 1-2 125,8470 m
Lo2-3: Longitud inicial del cable 2-3 186,8552 m
Flecha bajo peso propio en la posicin de carga puntual 29,2759 m
El modelo usado est compuesto por dos elementos cable, los cuales han sido derivados
de la configuracin inicial bajo peso propio. La carga concentrada se aplica incrementalmente
hasta completar dicha carga, siendo la diezmilsima parte de la carga total (35,586 N /1 0000)
un incremento de carga aceptable, tal como se aprecia en la figura (4.4).

g- -6
-9L---~~~~~~~--~~~~~~---L--~~~~L_--~_L_L~~~
3 5
10 10 10
Nmero de Incrementos
Figura 4.4: Desplazamiento vertical del nodo 2 versus nmero de incrementos de carga.
oQ.
E
F! 100
OL__L~~Lll~~~~~~~~~LL~~~llW
3
10 10
Figura 4.5: Tiempo de convergencia versus nmero de incrementos de carga.
De las figuras (4.4) y (4.5), el mtodo de la secante se aproxima rpidamente a la solucin
con una menor cantidad de incrementos y en un menor tiempo de convergencia comparado
al mtodo de Newton.
Los desplazamientos calculados, son aceptables al compararlos con los resultados de
diversos investigadores, tal como se aprecia en la tabla (4.5).

Tabla 4.5: Comparacin de desplazamientos del nodo 2 (m).
Investigador Tipo de elemento Desplazamiento vertical Desplazamiento horizontal
Michalos y Birnstiel (11] Recta elstica -5,472 -0,845
O'Brien [11] Catenaria elstica -5,627 -0,860
Jayaraman y Knudson (15] Catenaria elstica -5,626 -0,859
Tibert [25] Parbola elstica -5,601 -0,866
Tibert [25] Catenaria asociada -5,656 -0,860
Tibert [25] Catenaria elstica -5,626 -0,859
Andreu [1] Catenaria elstica -5,626 -0,860
Huu y Seung (11] Catenaria elstica -5,626 -0,859
Presente trabajo Catenaria elstica -5,627 -0,860
4.1.3. Armadura bidimensional
Los efectos lineales y no lineales de segundo orden (P-~) en un elemento armadura, se
determinan usando el algoritmo desarrollado y se comparan con el programa SAP2000. Esta
aplicacin es planteada por Fook Heng [9], cuyos datos se encuentran en la tabla (4.6). La
aplicacin consiste en colocar una carga vertical P en los nodos 1 y 2, y una fraccin a de
dicha carga como fuerza horizontal en el nodo 1, tal como es mostrado en la figura (4.6).
7m
2m. J
'
Figura 4.6: Geometra de armadura bidimensional.

Tabla 4.6: Propiedades de la armadura bidimensional.
ltem Dato
A: rea de seccin transversal 157 cm 2
E: Mdulo de elasticidad 205 N/cm2
A continuacin se muestran las tablas de resultados obtenidos:
Tabla 4.7: Desplazamiento del nodo 1 en cm, con o: =0,5.
SAP2000 Presente trabajo Diferencia ( %)

p (N)
Lineal No lineal Lineal No lineal Lineal No lineal
o 0,0000 0,0000 0,0000 0,0000 - -

50 7,0945 7,3605 7,0946 7,3901 0,00 -0,40
100 14,1891 15,3824 14,1891 15,4252 0,00 -0,28
150 21,2836 24,1321 21,2837 24,1987 0,00 -0,28

200 28,3781 33,7733 28,3783 33,8249 0,00 -0,15
250 35,4726 44,3400 35,4729 44,4457 0,00 -0,24
300 42,5671 56,0027 42,5674 56,2413 0,00 -0,43
350 49,6617 69,0266 49,662 69,4458 0,00 -0,61
400 56,7562 84,0064 56,7566 84,3743 0,00 -0,44
450 63,8507 10Q,7695 63,8512 101,4683 0,00 -0,69
Tabla 4.8: Desplazamiento del nodo 1 en cm, con o: =0,7.

p (N)
o 0,0000 0,0000 0,0000 0,0000 - -
50 9,9406 10,2972 9,9407 10,3629 0,00 0,63
100 19,8813 21,6187 19,8814 21,6505 0,00 0,15
150 29,8219 33,9628 29,8221 34,0036 0,00 0,12
200 39,7626 47,3779 39,7629 47,5984 0,00 0,46
250 49,7032 62,5121 49,7036 62,6606 0,00 0,24
300 59,6439 79,2348 59,6443 79,4882 0,00 0,32
350 69,5845 97,9714 69,585 98,4909 0,00 0,53
400 79,5252 119,2117 79,5257 120,2623 0,00 0,87
450 89,4658 143,5468 89,4664 145,7285 0,00 1,50

Tabla 4.9: Desplazamiento del nodo 1 en cm, con a =0,9.

p (N)
o 0,0000 0,0000 0,0000 0,0000 - -

50 12,7869 13,3407 12,7868 13,3274 0,00 O, 1
100 25,5737 27,8993 25,5735 27,8138 0,00 0,31
150 38,3606 43,8720 38,3603 43,7362 0,00 0,31
200 51,1474 61,5116 51,1471 61,3909 0,00 0,2
250 63,9343 81,1536 63,9338 80,9073 0,00 0,3
300 76,7211 103,2638 76,7206 102,7169 0,00 0,53
350 89,508 128,5274 89,5073 127,355 0,00 0,91
400 102,2948 158,039 102,2941 155,4844 0,00 1,62
450 115,0817 193,7816 115,0809 189,0041 0,00 2,47
De las tablas (4.7), (4.8) y (4.9), en el anlisis lineal no existe diferencia con el programa
SAP2000, pero en el anlisis no lineal la diferencia es menor a 2,47% cuando se aplica un
90% de la carga vertical, lo cual es aceptable.
La figura (4.7), muestra la diferencia de considerar un comportamiento lineal y no lineal
en la determinacin de desplazamientos.
. -- ! ..
:
........................ 1 - ..
! -Lineal (a=O,S)
f .. ...... : ---~No lineal (a=O~)
-Lineal (a=0,7)
....:...................:................... i. .................. --No lineal (a0,7)
. -Lineal (<>=0,9)
-- .. No lineal (a=0,9)
:: .. -:--- .......... ! ................ : .. . . . ........ --~- .. -- ........ --. -= ............. .
40 60 80 100 120 160 180 200

Desplazamiento (cm)
Figura 4.7: Diagrama carga P versus desplazamiento lateral para cada a.

En la figura (4.8} se muestra la razn de desplazamientos entre el comportamiento no
lineal y el comportamiento lineal, en la cual se aprecia que sta diferencia se incrementa con
el aumento de la carga P y por ende la carga lateral aP.
(V
~ 1.8.------.-------,-------,-------.,,------.-------.-------~,.------~
i ~ 1.s .....~ ................ ~ ............... , ................ : ............... -... .............. ~
~ 1.4 ;!:~~ : :
l '~: - l - - : - - - - r - r - -1 - - -- :~ :~
lil Q)
~e 08
so 100 150 200
i :
250
i
300
:
350 400 450
p (N)
Figura 4.8: Relacin de desplazamiento no lineal y lineal para cada valor de a.
4.1.4. Armadura tridimensional
Para verificar el algoritmo desarrollado en armaduras tridimensionales, se anliza una
seccin de torre de 3,00 m de altura, sujeta a una carga lateral de 1000 kg.
1000 kg
.//
\
\\
/1
\\
/_
'\'\ -81
-82
--~ -83
84
//
\,\
<
X L-.l_
30cm.
Figura 4.9: Geometra de torre tridimensional.

El presente modelo tambin es analizado en SAP2000, para verificar la diferencia entre
ambos programas. La configuracin geomtrica y sus propiedades fsicas son presentadas
en la figura (4.9) y en la tabla (4.1 O) respectivamente.
Tabla 4.1 O: Propiedades de la armadura tridimensional.
Seccin Tipo rea (cm 2 ) E (kg/cm 2 ) Densidad (kgfcm 3 )
81: L2 11 x2 11 x3/16 11 4,6119 2038901,9 0,007849

82: L1 5"x1 5"x3/16" 3,4022 2038901,9 0,007849
' '
83: PX1,5" 6,8912 2038901,9 0,007849
84: 05/8" 1,9793 2038901,9 0,007849
La denominacin de los nodos y la geometra tridimensional, se muestran en la figura
(4.10).
: .
...
300
-
250
:--. ~-;gl ... ............:
.. :... : :..
.
200
-
'
E'
-
u
N
w
150
iil .
100
r--r-
. ..
'8: .......:............ :
.. ..
' ... ...
! ~
50 .'
: ....... .
100
50 100 EJE Y(cm}

EJE X(cm}
Figura 4.1 O: Geometra con los nodos de la estructura tridimensional.

El primer anlisis es el lineal, para ello se resuelve la matriz global de rigidez elstica Ke
en un nico incremento de carga, y se obtienen los resultados listados en la tabla (4.11 ).
Tabla 4.11: Desplazamientos del anlisis lineal (cm).

nodo
Ux Uy Uz Ux Uy Uz Ux Uy Uz
4 0,05227 -0,01684 -0,01347 0,0523 -0,01685 -0,01348 -0,06 -0,09 -0,09
5 0,05227 -0,01684 0,02887 0,0523 -0,01685 0,02889 -0,06 -0,09 -0,09
6 0,05227 -0,01684 -0,0154 0,0523 -0,01685 -0,01541 -0,06 -0,09 -0,09
7 0,16287 -1,3E-14 -0,02694 O, 16298 -7,1 E-15 -0,02697 -0,07 43,25 -0,09
8 O, 16287 3,6E-15 0,05388 o, 16298 -2,3E-15 0,05393 -0,07 164,81 -0,09
9 O, 16287 -1 ,3E-14 -0,02694 O, 16298 -7,1E-15 -0,02697 -0,07 43,20 -0,09
10 0,32347 -0,01203 -0,03656 0,32371 -0,01204 -0,0366 -0,08 -0,09 -0,09
11 0,32347 -0,01203 0,07505 0,32371 -0,01204 0,07512 -0,08 -0,09 -0,09
12 0,32347 -0,01203 -0,03849 0,32371 -0,01204 -0,03852 -0,08 -0,09 -0,09
13 0,52573 -3,2E-14 -0,04619 0,52614 -2,3E-14 -0,04623 -0,08 27,86 -0,09
14 0,52573 -7,8E-16 0,09237 0,52614 -1 ,2E-14 0,09245 -0,08 -1428,63 -0,09
15 0,52573 -3,2E-14 -0,04619 0,52614 -2,3E-14 -0,04623 -0,08 27,86 -0,09
16 0,76132 -0,00722 -0,05196 0,76193 -0,00722 -0,05201 -0,08 -0,1 o -0,09
17 0,76132 -0,00722 O, 10584 0,76193 -0,00722 O, 10594 -0,08 -0,10 -0,09
18 0,76132 -0,00722 -0,05388 0,76193 -0,00722 -0,05393 -0,08 -0,10 -0,09
19 1,02191 -5,4E-14 -0,05773 1,02274 -4,3E-14 -0,05778 -0,08 21,60 -0,09
20 1,02191 -7,6E-15 O, 11546 1,02274 -2,7E-14 o, 11557 -0,08 -261 '15 -0,09
21 1,02191 -5,4E-14 -0,05773 1,02274 -4,3E-14 -0,05778 -0,08 21,39 -0,09
22 1,29916 -0,00241 -0,05966 1,30023 -0,00241 -0,05971 -0,08 -0,11 -0,09
23 1,29916 -0,00241 0,12124 1,30023 -0,00241 0,12135 -0,08 -0,11 -0,09
24 1,29916 -0,00241 -0,06158 1,30023 -0,00241 -0,06164 -0,08 -0,11 -0,09
25 1,58475 -7,2E-14 -0,06158 1,58606 -6,8E-14 -0,06164 -0,08 6,49 -0,09
26 1,58581 -0,00184 0,12316 1,58713 -0,00184 0,12327 -0,08 0,00 -0,09
27 1,58475 -7,2E-14 -0,06158 1,58606 -6,8E-14 -0,06164 -0,08 6,46 -0,09
En la tabla (4.11 ), se puede ver aparentemente grandes diferencias en las componentes
uy; sin embargo, los desplazamientos son no apreciables debido a que el desplazamiento de
esos nodos en la componente uy est en el orden de magnitud de 10- 14 .
El anlisis no lineal se realiza con incrementos de la diezmilsima parte de la carga total

externa (1 000 kg/1 0000), y se obtienen los desplazamientos mostrados en la tabla (4.12),
cuya deformada se muestra en la figura (4.11).
Tabla 4.12: Desplazamientos del anlisis no lineal (cm).

nodo
4 0,05191 -0,01702 -0,01353 0,05194 -0,01704 -0,01354 -0,06 -0,1 o -0,09
5 0,05235 -0,01639 0,02882 0,05238 -0,01641 0,02885 -0,06 -0,1 o -0,09
6 0,05259 -0,01702 -0,01546 0,05262 -0,01704 -0,01547 -0,06 -0,09 -0,09
7 0,16281 -0,00004 -0,02719 0,16292 -0,00004 -0,02721 -0,07 -1,20 -0,09
8 O, 16304 0,00004 0,05367 o, 16316 0,00003 0,05372 -0,07 1'11 -0,09
9 O, 16293 -0,00004 -0,02719 o, 16304 -0,00004 -0,02721 -0,07 0,59 -0,09
10 0,32325 -0,01216 -0,03717 0,32349 -0,01218 -0,0372 -0,08 -0,09 -0,09
11 0,32378 -0,01171 0,0745 0,32402 -0,01172 0,07457 -0,08 -0,09 -0,09
12 0,32375 -0,01216 -0,03909 0,32399 -0,01218 -0,03913 -0,08 -0,1 o -0,09
13 0,5257 -0,00007 -0,04735 0,52611 -0,00007 -0,04739 -0,08 0,38 -0,09
14 0,52621 0,00013 0,09128 0,52662 0,00013 0,09137 -0,08 -0,12 -0,09
15 0,52594 -0,00007 -0,04735 0,52635 -0,00007 -0,04739 -0,08 0,25 -0,1 o
16 0,7613 -0,00732 -0,05386 0,76192 -0,00733 -0,05391 -0,08 -0,1 o -0,10
17 0,76199 -0,00698 O, 10403 0,76261 -0,00699 0,10412 -0,08 -0,1 o -0,09
18 0,76168 -0,00732 -0,05579 0,76229 -0,00733 -0,05584 -0,08 -0,10 -0,09
19 1,02204 -0,00008 -0,06055 1,02288 -0,00008 -0,0606 -0,08 -0,22 -0,1 o
20 1,02278 0,00018 o, 11276 1,02362 0,00018 0,11286 -0,08 0,15 -0,09
21 1,02235 -0,00008 -0,06055 1,02319 -0,00008 -0,06061 -0,08 -0,24 -0,1 o
22 1,29946 -0,00249 -0,0635 1,30053 -0,00249 -0,06356 -0,08 -0,09 -0,1 o
23 1,30025 -0,00222 o, 11752 1,30133 -0,00222 0,11763 -0,08 -0,08 -0,09
24 1,29976 -0,00249 -0,06542 1,30084 -0,00249 -0,06549 -0,08 -0,08 -0,10
25 1,58522 -0,00008 -0,06651 1,58655 -0,00008 -0,06658 -0,08 0,47 -0,10
26 1,58711 -0,00164 O, 11837 1,58844 -0,00164 0,11847 -0,08 0,01 -0,09
27 1,58555 -0,00008 -0,06651 1,58688 -0,00008 -0,06658 -0,08 0,38 -0,1 o
En la tabla (4.12), la diferencia de los desplazamientos tiene un mximo de 1,2 %, lo que
demuestra una buena aproximacin del algoritmo propuesto con el programa SAP2000, a
pesar de usar matrices de orden superior.

.r .
. :
:
300

. . . . .. . . .... : .
~
\ :
250 i
... ..;
200
. .. ..: .... .. . .
e
u
150
w
ii3 :
100
...............
50
-:-:"'--.---------- . :-------------
100
.............. ::= .......
..
o ... o
-100 -50
o 50 100
-100
EJE Y(cm)
EJE X(cm)
Figura 4.11: Deformacin del anlisis no lineal (no a escala).
Si se analiza la aproximacin tomando como parmetro el nmero necesario de incre-
mentos de carga, el mtodo de la secante se aproxima ms rpido a la solucin, como se
observa en la figura (4.12).
X 10-a
-1.615r---~--~------,
0.119
-1.62
E
u
';;' 1.51185
E
u
-;: -1.625
0.1188
:S :S
j -1.63
~
:;
'ii. -1.635
:;
e
0.118
- 1 . 64 ~o'-:,-~~1~o2:---~~1o7'~~.......
1o4 0.1170 '-,--~---~---l
2 4
10 10' 10 1i 10' 10
Nmero de incrementos Nmero de incrementos Nmero de incrementos
1-- Secante - - Newton - - - Solucin 1
Figura 4.12: Desplazamientos con nmero de incrementos para el nodo 26 de la armadura,

al considerar anlisis no lineal.

Al considerar como parmetro el tiempo de solucin, ambos mtodos tienen similar tiempo
de convergencia, salvo en el rango de 800 a 4500 incrementos de carga, donde el mtodo de
la secante es ms rpido, como se muestra en la figura (4.13).
Figura 4.13: Nmero de incrementos de carga versus tiempo de solucin.
En la figura (4.13), se observan particulares saltos, los que indican mejor estabilidad del
mtodo de la secante para nmero de incrementos de carga menores a 4000.
4.1.5. Armadura tridimensional con cables atirantados
La siguiente aplicacin tiene por objetivo acoplar una armadura tridimensional con cables
atirantados en una sola estructura, a la cual luego de determinar sus desplazamientos en
condiciones de peso propio y cuando al cable es tensado en la parte inferior con 30 kg, se
le aplica una carga lateral de 1000 kg. Esta estructura tambin es analizada con el programa
SAP2000 para verificar los desplazamientos nodales.
Otro objetivo es verificar la cantidad de elementos cable necesarios para obtener resulta-
dos ms precisos en la solucin.
Tabla 4.13: Propiedades del elemento cable.
Seccin Tipo rea (cm 2 ) E (kgjcm 2 ) Densidad (kgjcm 3 )

85: 01/ 4" 0,3167 2038901,9 0,007849
ANLISIS NO LINEAL DE ESTRUCTURAS CABLE ARMADURA 65

Las propiedades fsicas del elemento armadura son similares a las de la aplicacin 4.1.4,
mientras que las propiedades del elemento cable se encuentran en la tabla (4.13). El modelo
usado se muestra en la figura (4.14).
300
250
200
e
.
~150
iii"
100
50
o
-300
300
Figura 4.14: Geometra de torre tridimensional arriostrada.
El anlisis no lineal por peso propio y el cable tensado, se realiza con incrementos de
la centsima parte de la carga (peso propio/1 00), obtenindose la deformacin mostrada
en la figura (4.15) y desplazamientos listados en la tabla (4.14). Cabe indicar que dichos
desplazamientos son similares cuando al cable que arriostra se subdivide en varios elementos
antes de aplicar la carga lateral.
3!ll
250
200
e
;:15)
..
iii"
100
50
o
-300 300
Figura 4.15: Geometra deformada por peso propio de torre tridimensional arriostrada (no a
escala).

Tabla 4.14: Desplazamientos del anlisis no lineal por peso propio (cm).
SAP2000 Presente trabajo Diferencias ( %)

nodo
4 -0,00029 -0,00017 -0,00013 -0,00029 -0,00017 -0,00013 -0,23 -0,05 o, 11

5 2,03E-9 0,00033 -0,00013 5,20E-9 0,00033 -0,00013 -156,04 -0,05 o, 11
6 0,00029 -0,00017 -0,00013 0,00029 -0,00017 -0,00013 -0,23 -0,05 o, 11
7 -0,00002 -0,00001 -0,00026 -0,00002 -0,00001 -0,00026 1,78 0,78 -0,06
8 1,745-9 0,00002 -0,00026 1 ,77E-9 0,00002 -0,00026 -1,55 0,76 -0,06
9 0,00002 -0,00001 -0,00026 0,00002 -0,00001 -0,00026 1,78 0,75 -0,06
10 -0,00027 -0,00015 -0,00037 -0,00027 -0,00015 -0,00037 0,02 -0,08 0,02
11 1,48E-9 0,00031 -0,00037 4,20E-9 0,00031 -0,00037 -183,63 -0,08 0,02
12 0,00027 -0,00015 -0,00037 0,00027 -0,00015 -0,00037 0,02 -0,08 0,02
13 -0,00004 -0,00002 -0,00047 -0,00004 -0,00002 -0,00047 -0,62 0,76 -0,03
14 1 ,24E-9 0,00005 -0,00047 1 ,30E-9 0,00005 -0,00047 -4,80 0,76 -0,03
15 0,00004 -0,00002 -0,00047 0,00004 -0,00002 -0,00047 -0,62 0,75 -0,03
16 -0,00025 -0,00014 -0,00057 -0,00025 -0,00014 -0,00057 -0,1 o -0,12 -0,15
17 1 ,03E-9 0,00028 -0,00057 3,32E-9 0,00028 -0,00057 -223,77 -0,12 -0,15
18 0,00025 -0,00014 -0,00057 0,00025 -0,00014 -o,opo57 -0,1 o -0,12 -0,15
19 -0,00006 -0,00004 -0,00065 -0,00006 -0,00004 -0,00065 0,19 0,76 -0,15
20 8,36E-10 0,00007 -0,00065 9,55E-10 0,00007 -0,00065 -14,25 -0,64 -0,15
21 0,00006 -0,00004 -0,00065 0,00006 -0,00004 -0,00065 0,19 0,76 -0,15
22 -0,00023 -0,00013 -0,00073 -0,00023 -0,00013 -0,00073 -0,24 -0,16 -0,08
23 6,27-1 o 0,00026 -0,00073 2,53E-9 0,00026 -0,00073 -303,11 -0,17 -0,08
24 0,00023 -0,00013 -0,00073 0,00023 -0,00013 -0,00073 -0,24 -0,17 -0,08
25 -0,00010 -0,00009 -0,00080 -0,00010 -0,00009 -0,00080 -0,30 -0,47 -0,08
26 -0,00002 0,00013 -0,00080 -0,00002 0,00013 -0,00080 0,89 -0,03 -0,08
27 0,00013 -0,00005 -0,00080 0,00013 -0,00005 -0,00080 -0,09 0,75 -0,08
En la tabla (4.14), se muestran aparentemente grandes diferencias en la componente ux,
pero los desplazamientos son no apreciables debido a que el desplazamiento de esos nodos
en la componente ux est en el orden de magnitud de 10-9 .
Cuando se le aplica la carga lateral de 1000 kg, se registra el desplazamiento del nodo 26
para diversas subdivisiones del elemento cable, tal como se muestra en la figura (4.16), con la
finalidad de determinar la cantidad necesaria de elementos cable. De esta figura se observa
una tendencia a emplear ms de 1O divisiones por cable. Cabe indicar que mientras ms

elementos cable se adicionen, este difiere ms con los resultados del programa SAP2000.
- 0.6701
!.. 0.67
~ 0.6699 -~ . ,.; ... i.... i....~---::-~
1:0;95
.uu
J,1r,,;J\,_,:_,;
.. ~----~ .... ~:-
'l,
. : : ~ ~
0.66
94
o 1 2 3 4 5 6 7 B 9 10
-7:1. L....C~~_.__.....__._~.....___,
O 1 2 3 4 5 6 7 B 9 ID om4 o~;-2~~~~~5~:~~~;~a~9~1u
fl Elementos cable # Elementos cable fl Elementos cable
1 ~Presenta trabajo SAP2000 1
Figura 4.16: Desplazamientos en las tres componentes del nodo 26 con diferentes cantida-
des de elementos cable.
4.1.6. Armadura tridimensional con cables pre-esforzados
Se analiza el aporte en rigidez y el comportamiento de los cables atirantados con dife-
rentes cargas de pre-esforzado (T0 ). Para ello se considera un elemento armadura vertical
que tiene una seccin adecuada de rigidez axial y una deformacin insignificante debida a la
tensin ejercida en los cables, al cual se le aplican diversas cargas laterales en el nodo 2 en
la direccin x, con sentido positivo y negativo.

300 ............. , ............................. : ..
2lJJ .........
Figura 4.17: Modelo tridimensional analizado.
Los resultados obtenidos para diversas tensiones y fuerzas laterales en la direccin del
eje x se muestran en la figura (4.18).

Tabla 4.15: Propiedades de los elementos.
Seccin Tipo
Armadura: 05cm 19,63 2038901,9 0,007849

Cable: 01/4" 0,3167 2038901,9 0,007849
2500
2000
1500
a
:::. 1000
;
]!
:a"'
Cl
200~
4
Figura 4.18: Diagrama de carga versus desplazamiento lateral del nodo 2.
Si se realiza un ajuste lineal de datos en cada curva de tensin, tal como se indica en
la figura (4.19}, se pueden determinar los puntos de inflexin (xi+,pi+) y (xi-,pi-) de cada
curva, tal como se lista en la tabla (4.16}.
- P.Oefonnacl6n
-l\lustallneal
p+ .....

xj
............. ~-
Figura 4.19: Envolvente de cada curva de tensin.

Tabla 4.16: Puntos de inflexin en cada curva de tensin.
p,+ X.,+
Tensin (kg) P;+ (kg) P;- (kg) p.- x;+ (cm) x;- (cm) Xi-
50 331,4 -75,9 -4,37 0,363 -0,083 -4,37

100 404,1 -119,1 -3,39 0,361 -0,107 -3,37
200 571,7 -232,6 -2,46 0,501 -0,204 -2,46
300 761,3 -322,9 -2,36 0,667 -0,283 -2,36
500 1124,4 -523,5 -2,15 0,985 -0,459 -2,15
1000 2066,6 -1042,1 -1,98 1,811 -0,913 -1,98
La figura (4.20), relaciona el punto de inflexin del lado positivo determinado por la coor-
denada (xi+,pi+), con la inflexin del lado negativo descrito por la coordenada (xi-,pi-).
donde para tensiones mayores a T =1 000 kg se hace asinttica, es decir:

.
x< ~ 2 y PP< ~ 2.
~ t
..... -
1
::e
'+3 ...................................... -~ .............. : ............ ............. : ............ ............ -
-~ ~
:1<
e
:22
M
7ii i
~ 1~--~=---~~--~~----~----~----~----~~--~~--~~--~-
i 1 i
0 100 200 300 400 600 600 700 800 900 1000
Tensin del cable (kg)
Figura 4.20: Relacin de las fuerzas ( ~; ~) y desplazamientos ( ~; ~) versus tensin de los
cables.
Comparando las rigideces asintticas del lado positivo k+ y negativo k-, como se muestra
en la figura (4.21 ), el aporte de rigidez de los tres cables en la direccin x+ es aproximada-
mente el doble de la rigidez en la direccin x-' es decir: z~ ~ 2.
3~----------~----.-----,!-----,-----.-----.------.-----,.,----~
............. .; ............... :. . . . . ..... ....... . . .. ......... -~- .............. -~- ...............:................:................: ............... ~ ... ' ........ .
. . .
<>-~~~--~---~+--------..----------------------------
............. _, ...............; .............. r .............. t ............... ;............... :.............. -~ .............. r .............. ,............. .
1
i i i i
o 100 200 300 400 500 600 700 800 900 1000
Tensfon del cable (kg)
Figura 4.21: Relacin de rigideces asintticas.

Al relacionar la rigidez inicial ki con la rigidez asinttica del lado negativo k-, tal como
se muestra en la figura (4.22}, la relacin de rigideces para una tensin T 0 ;::: 100 kg es
~ ~3,07, y para las condiciones descritas en sta aplicacin, se tiene k- ~370 kgjcm.
5.-----~----~----~----~----~----~----~----~----------~
! ! !
~ 4rJ:+> :t::t-
~ :"~;'r~~~I-~J~--]~~~~I--I~--r~-r.,-r -~:
i i 1
.
100 200 300 400 500 600 700 800 900 1000
Tensin del cable (kg)
Figura 4.22: Relacin de rigidez inicial y asinttica negativa.
Con las relaciones de rigidez y el ajuste lineal de cada curva, es posible estimar el aporte
de rigidez de un cable tridimensional y realizar un equivalente para diversos modelos que
involucren a este elemento.
4.1.7. Emparrillado tridimensional de cables
En esta aplicacin se presenta una estructura de slo elementos cable. sta consiste en
un emparrillado de cables tensados con una configuracin geomtrica mostrada en la figura
(4.23) y cuyas propiedades fsicas son listadas en la tabla (4.17).
T 3 ,,.
30:r-m~w~3 --~--~--~L-~--~~~ 7
3l~Wf!-7 __l2'-\!.L_ _ ~ .L'o.!JL--"t----'-"'"-L-~

__
~r
~ ~
}-30,48 m---+--30,48 m+-30,48 m---..,f
Figura 4.23: Geometra del emparrillado de cables.

Tabla 4.17: Informacin para los cables.
ltem Dato
w: Peso propio del cable 1,459 N/m

A: rea de seccin transversal 146,4513 mm2
E: Mdulo de elasticidad 82,737 MPa

Ton: Tensin de cables horizontales 24,2828 kN
To: Tensin de cables inclinados 23,6868 kN
Pv: Carga vertical para nodos internos 35,5858 kN
.6.z: Desnivel entre nodos 9,144 m
En la configuracin tensada de la figura (4.23), se aplican las cargas verticales Pv en los
nodos 4, 5, 8 y 9. El anlisis se realiza para incrementos de carga de la milsima parte de Pv
(Pv/1 000). Los resultados de los desplazamientos en el nodo 8 son listados en la tabla (4.18);
la configuracin inicial y su deformada, se presentan en la figura (4.24).
Tabla 4.18: Comparacin de desplazamientos en nodo 8 (mm).
Investigador Tipo de elemento u, Uy Uz
Jayaraman y Knudson[15] Catenaria elstica -39,62 -40,20 -446,32
Tibert[25] Parbola elstica -40,78 -40,78 -453,33
Tibert[25] Catenaria asociada -40,78 -40,78 -453,36
Tibert[25] Catenaria elstica -40,48 -40,39 -450,04
Huu y Seung[11] Catenaria elstica -40,13 -40,13 -446,50
SAP2000 Catenaria elstica -40,47 -40,47 -449,46
Presente trabajo Catenaria elstica -40,53 -40,53 -450,12
... : ....
.. :
:. .. : .. : ..
o ~- ...
.... : .... ~ . . ..
-2
.. : .: ..
"E -4
...........

Cl) -6 :.
iii' .. ;. ...-: ..
-8 : _.. ---
____ ... ::.
::;. __
-10
o
20
Figura 4.24: Configuracin inicial y deformada del emparrillado de cables.

4.2. ANLISIS DINMICO
El anlisis dinmico se realiza con el registro ssmico de El Centro (1940) el cual se
muestra en la figura (4.25).
OA,---------,---------,----------.---------,---------,---------,
0.3 .............................................................. ;........ .
0.2
-0.1
S
e

O
o
.
;;
e( ~.1
.0.3 .............. . ....................... ! ........................................................ .
~ 4o'----------.L5________---,~:10,----------15L,-----------,2.L0________---,1:25,-----------'30
Tiempo (seg)
Figura 4.25: Registro de aceleraciones, El Centro 1940.
4.2.1. Armadura tridimensional
Una armadura de 72 elementos es mostrada en la figura (4.26) con el objetivo de verificar
la efectividad de predecir los efectos no lineales cuando es sometida a una aceleracin ss-
mica. Esta aplicacin es mostrada por Huu [12] y los valores de las propiedades fsicas son
presentados en la tabla (4.19).
Tabla 4.19: Informacin de la armadura tridimensional.
ltem Dato
A: rea de seccin transversal 38,485 cm 2

E: Mdulo de elasticidad 200 GPa
M: Masa concentrada en cada nodo 100 kNs 2 /m
~i y~{ Fraccin de amortiguamiento crtico 0,05

Figura 4.26: Geometra de la armadura tridimensional con masas concentradas.
La figura (4.27) muestra el registro de desplazamientos del nodo 19 en la direccin x
comparado al obtenido por Huu. En la tabla (4.20) se muestran que las diferencias de los
desplazamientos mximos son menores al 0,46 %, lo que evidencia una baja discrepancia
estimada entre Huu y el presente trabajo.
150r---------r---------,---------~---------------r=~======~
-Presente trabajo
---Huu
e 1oo
..
>C
50 ......................... .
.S
'
u
f!
:S
eCll
o
e
Cll
-~ -50
~
i5.
m
Cll
e -100 ....
"150oL__ _ _ __J5L___ _ _ _--:':---------:15':------:--::'::--------::.:---------=3o

Tiempo (seg)
Figura 4.27: Registro de desplazamiento del nodo 19 en direccin x.

Tabla 4.20: Comparacin de resultados de desplazamientos mximos.
Hu u Presente trabajo Diferencia ( %)
Max 132,21 131,60 0,46
Min -124,46 -124,85 0,31
4.2.2. Emparrillado de cables
Para la misma estructura de emparrillado de cables analizada con fuerzas estticas (4.1.7.),
se realiza el anlisis dinmico cuando la aceleracin ssmica acta en direccin vertical. Para
ste modelo se asignan masas concentradas de 8,76 kN s2 /m en los nodos 4, 5, 8 y 9 tal
como se muestra en la figura (4.28), y se emplea un factor de amortiguamiento crtico de
0,05.

IMovimtento
SISffiiCO
Figura 4.28: Emparrillado de cables con masas concentradas.
El anlisis dinmico comienza cuando el emparrillado se encuentra en la condicin de
equilibrio indicada en la aplicacin 4.1.7.
El registro de desplazamiento vertical del nodo 4 se muestra en la figura (4.29), el cual es
comparado con el obtenido del programa SAP2000.

-100-------,-----,--------,-----,----=======-
-Presente trabajo
---SAP2000
eg -2oo
;N -300
e:
O
'
u
....
Q)
'6
e:
Q)
o ...... ...... ~
e:
Q)
1
~ 1
f: .t .. -
~
c.
t/1
~ -700 ........... !''. -
~o~L__ _ _ _ _sL__ _ _ _~1oL-----~1sL-----~2o-----~25----~3o
Tiempo (seg)
Figura 4.29: Registro de desplazamiento del nodo 4 en direccin vertical.
En la figura (4.29) y la tabla (4.21 ), se observa que los registros obtenidos por el progra-
ma SAP2000 y el presente trabajo se ajustan perfectamente en los primeros 1O seg, luego se
aprecia un desfase que se incrementa con el tiempo entre ambos registros, pero los despla-
zamientos mximos tienen una diferencia menor al 3,04 %.
Tabla 4.21: Comparacin de resultados de desplazamientos mximos.
Max -196,79 -202,77 -3,04
Min -728,84 -728,47 0,05

FACULTAD DE INGENIERA CIVIL CAPTULO V: CONCLUSIONES V RECOMENDACIONES
CAPTULO V
CONCLUSIONES Y RECOMENDACIONES
5.1. CONCLUSIONES
Es posible acoplar elementos armadura y cable en una estructura, y considerar la no

linealidad de dichos elementos.
Las ecuaciones variacionales del elemento armadura, se pueden expresar en forma
matricial con reacomodos algebraicos como los indicados en el Captulo 2.
La formulacin Lagrangiana proporciona un mtodo eficaz para analizar el comporta-
miento de un elemento armadura o cable. En este trabajo para el elemento armadura
se emplea la formulacin variacional no lineal propuesta por Yang [27], la que permite
deducir matrices que incluyan los efectos elsticos, geomtricos y de orden superior;
mientras que para el elemento cable, a partir de ecuaciones analticas debidamente
comprobadas, permite deducir matrices no lineales.
El anlisis de una estructura cable - armadura compuesta por elementos rgidos y flexi-
bies, se puede resolver considerando las siguientes etapas:
- Comenzar con el ensamble de los elementos desarrollados a partir de la formu-

lacin que emplea coordenadas Lagrangianas; si se trata de un elemento cable,
obtener la posicin de uno de los nodos a partir del nodo identificado y caracteri-
zado, con fines de tener la configuracin inicial de equilibrio por peso propio.
- Para lograr una mejor discretizacin de un cable, segmentar el elemento en ele-
mentos ms pequeos, y las coordenadas de los nodos intermedios se obtienen
a partir de la metodologa propuesta en el acpite 3.2.2.

FACULTAD DE INGENIERA CIVIL CAPTULO V: CONCLUSIONES Y RECOMENDACIONES
- Finalmente, emplear un mtodo incremental iterativo para resolver el sistema de
ecuaciones no lineales en coordenadas globales.
En el anlisis de estructuras compuestas por elementos cable, el mtodo de la secante
es 1,2 veces ms rpido que el mtodo de Newton (figura 4.5), adems se aproxima
mejor a la solucin cuando se tiene menores incrementos de carga (figura 4.4); mien-
tras que en estructuras conformadas por elementos armadura, el tiempo es similar en
ambos mtodos (figura 4.13). Sin embargo, el mtodo de la secante se aproxima ms
a la solucin cuando los incrementos de carga son cada vez ms pequeos.
La relacin entre los desplazamientos no lineales y desplazamientos lineales, se incre-
mentan con el aumento de la carga lateral. En el caso de la armadura bidimensional
(aplicacin 4.1.3.) esta razn es mayor o igual a 1 para cargas de P=50 N.
Es posible determinar una curva de rigidez simplificada de la relacin Carga - Des-
plazamiento de una estructura atirantada. En la armadura tridimensional con cables
pre esforzados (aplicacin 4.1.6.) esto ocurre para tensiones (T0 ) del cable mayores o
iguales a 50 kg.
El registro tiempo historia de estructuras compuestas por elementos armadura y estruc-
turas compuestas por elementos cable, tiene desplazamientos mximos y mnimos del
mismo orden de magnitud que los obtenidos por el investigador Huu [12] y el programa
[5] respectivamente.
Los resultados obtenidos a partir de los elementos y la metodologa propuesta, mues-
tran concordancia con los resultados obtenidos con programas de uso profesional y
con aquellos resultados reportados por otros investigadores, tal como se observa en el
Captulo 4.
5.2. RECOMENDACIONES
Se recomienda emplear la formulacin Lagrangiana para el anlisis inicial de un ele-
mento cable, para considerar el efecto de su peso propio.
Se recomienda hacer un estudio no lineal elstico del elemento cable, considerando la
variacin del rea transversal del cable producto de la elongacin de ste.

FACULTAD DE INGENIERA CIVIL CAPTULO V: CONCLUSIONES Y RECOMENDACIONES
Usar el mtodo de control de desplazamiento generalizado propuesto por Yang y Shieh
[27], para trazar la trayectoria no lineal de la estructura.
Se recomienda usar el mtodo de Newton el mtodo de la secante para el estudio no
lineal elstico que considere las hiptesis de la presente tesis.
Se debe analizar el elemento cable considerando por lo menos un comportamiento no
lineal elstico.

FACULTAD DE INGENIERA CIVIL BIBLIOGRAFA
BIBLIOGRAFA
[1] Andreu A., Gil L. y Roca P., "A new deformable catenary element for the analysis of cable
net structures", Computer and Structures, 84, p. 1882-1890, 2006.
[2] Baron F. y Venkatesan M., "Nonlinear analysis of cable and truss structures", American
Society of Civil Engineers, 624.071, p. 679-71 O, 1971.
[3] Chopra A., "Dynamics of structures", Segunda edicin, Prentice Hall, N. Y., 1995.
[4] Clough R. y Penzien J., "Dynamics of structures", Tercera edicin, McGraw-Hill, N. Y.,
1995.
[5] Computer and Structures, "CSI Analysis Reference Manual", Berkeley, 201 O.
[6] Cuomo M. y Greco L., "A finite element cable for the analysis of cable nets", Faculty of
Engineering, University of Catania, Ancona, 2009.
[7] Engineering Systems, "MStower V6 User's manual", Londres, 2008.
[8] Flores L., "Mtodos numricos aplicados a la Ingeniera Civil", Universidad Nacional de
Ingeniera, Lima, 201 O.
[9] Fook N., "Geometric nonlinear analysis of truss structure', Bachelor Thesis, University of
Technology, Malasia, 201 O.
[1 O] Haijian S., "Nonlinear finite element modeling and characterization of guyed towers un-
der severe loading", PhD. Thesis, Faculty of the Graduate School, Missouri University,
Columbia, 2007.
[11] Huu T. y Seung K., "Nonlinear static and dynamic of cable structures", Finite Elements in
Analysis and Design, 47, p. 237-246, 2011.
[12] Huu T. y Seung K., "Practica! advanced analysis software for nonlinear dynamic analysis
of steel structures", Journal of Construcctional Steel Research, 67, p. 453-461, 2011.

[13] Huu T., "Practica! nonlinear inelastic static and dynamic analysis of space steel structu-
res", PhD. Thesis, Department of Civil and Environmental Engineering, Hanyang Univer-
sity, Seoul, 2009.
[14] lrvine M., "Cable structure", MIT Press, Cambridge, 1982.
[15] Jayaraman H. y Knudson W., "A curved element for the analysis of cable structure",
Computer and Structures, 14, 3-4, p. 325-333, 1981.
[16] Karoumi R., "Sorne modeling aspects in the nonlinear finite element analysis of cable
supported bridges", Computer and Structures, 71, p. 397-412, 1999.
[17] Kwan A. "A new aproach to geometric nonlinearity of cable structures", Computer and
Structures, 67, p. 243-352, 1998.
[18] Linero S. y Garzn A., "Elementos de la mecnica del medio continuo para cuerpos s-
lidos", Departamento de Ingeniera Civil y Agricola, Facultad de Ingeniera, Universidad
Nacional de Colombia, Bogot, 201 O.
[19] Madugula M., "Dynamic response of lattice towers and guyed masts", Report of American
Society of Civil Engineers, Virginia, 2002.
[20] Paliares M. y Rodriguez C., "Validacin de la formulacin numrica de la catenaria els-
tica con Ansys ", Escuela Regional de Matemticas, Universidad del Valle, Cali, XVI, 1,
p. 63-85, 2008.
[21] Paz M. "Dinmica estructural", Reverte, Tercera edicin, Estados Unidos, 1989.
[22] Pintea A., "Comparison between the linear and nonlinear responses of cable structures 1
- Static loading", Civil Engineering and Architecture, 55, 2, p. 182-188, 2012.
[23] Ren W., Huang M. y Hu W. "A parabolic cable element for static analysis of cable struc-
tures", lnternational Journal for Computer- Aided Engineering and Software, 25, 4, p.
366-384, 2008.
[24] Rohit B., "Dynamic analysis of guyed towers subjected to wind loads incorporating nonli-
nearity of the guys", Master Thesis, Texas Tech College, Lubbock, 1999.
[25] Tibert G., "Numerical analyses of cable roofs structures", Licentiate Thesis, Department
of Structural Engineering, Royal lnstitute of Technology, Estocolmo, 1998.
'
[26] Valles R., Reinhorn A. y Kunnath S., "IDARC 20 Version 4.0: A Program for the lnelastic
Damage Analysis of Buildings- Technical Report NCEER-96-001 0", University at Buffalo,
N. Y., 1996.
[27] Yang B., "Theory and analysis of nonlinear framed structures", Prentice Hall, Singapore,
1994.


Coarita Te PDF

Cargado por

Copyright:

Formatos disponibles

Coarita Te PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Coarita Te PDF

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE INGENIERA

FACULTAD DE INGENIERA CIVIL

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA

EVER DENNYS COARITA TINTAYA

do en todo momento, por sus consejos, sus valores y

su amor. A mis hermanos, por estar conmigo y apoyar-

me siempre, los quiero mucho. A todos mis amigos, por

compartir los buenos y malos momentos.

do Flores Gonzlez, por la orientacin y ayuda que me

LISTA DE FIGURAS VIl

LISTA DE SIGLAS Y SMBOLOS X

1.2 ESTADO DEL ARTE 2

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA iv

2.5.3 Mtodo de control de desplazamiento generalizado 31

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA V

y armadura mediante una formulacin en coordenadas Lagrangianas. Luego, se desarrollan

no lineales de una estructura en un sistema global de referencia, bajo la accin de fuerzas

externas estticas y dinmicas. Finalmente, se presentan algunas aplicaciones relativas a la

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA vi

Figura 2.1 Movimiento de un elemento armadura en el espacio. (Fuente: Yang,

Figura 2.8 Caracterstica general de sistemas no lineales. (Fuente: Yang, 1994). 32

Figura 3.1 Diagrama de flujo del elemento cable para Lo. . 36

Figura 4.1 Comparacin numrica de aplicacin con apoyos situados a distinto

Figura 4.2 Coordenada y tensin mnima en el punto ms bajo. 54

Figura 4.3 Configuracin inicial de su geometra bajo peso propio. 55

Figura 4.4 Desplazamiento vertical del nodo 2 versus nmero de incrementos de

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE ARMADURA vii

Figura 4.5 Tiempo de convergencia versus nmero de incrementos de carga. . 56 .

Figura 4.7 Diagrama carga P versus desplazamiento lateral para cada a. 59

Figura 4.8 Relacin de desplazamiento no lineal y lineal para cada valor de a. 60

Figura 4.9 Geometra de torre tridimensional. . . . . . . . . . . . . 60

Figura 4.1 O Geometra con los nodos de la estructura tridimensional. 61

armadura, al considerar anlisis no lineal. . . . . . . . . . . . 64

Figura 4.13 Nmero de incrementos de carga versus tiempo de solucin. 65

trada (no a escala). . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66

cantidades de elementos cable .. 68

Figura 4.21 Relacin de rigideces asintticas. . 70

Figura 4.22 Relacin de rigidez inicial y asinttica negativa. 71

Figura 4.24 Configuracin inicial y deformada del emparrillado de cables. . 72

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE- ARMADURA viii

Tabla 3.1 Formato de entrada de geometra y masa concentrada. 41

Tabla 4.1 Informacin para el problema .. 52

Tabla 4.7 Desplazamiento del nodo 1 en cm, con a =0,5. 58

Tabla 4.8 Desplazamiento del nodo 1 en cm, con a =0,7. 58

Tabla 4.1 O Propiedades de la armadura tridimensional. 61

Tabla 4.13 Propiedades del elemento cable. . . . . . . 65

Tabla 4.19 Informacin de la armadura tridimensional. . . . . . 73

Tabla 4.20 Comparacin de resultados de desplazamientos mximos. 75

Tabla 4.21 Comparacin de resultados de desplazamientos mximos. 76

ANLISIS NO LINEAL DE ESTRUCTURAS CABLE -ARMADURA ix

LISTA DE SIGLAS Y SMBOLOS

A rea de seccin transversal

Matriz de amortiguamiento viscoso global

Configuracin inicial no deformada

ltima configuracin calculada

Configuracn actual deformada

Matriz de correccin en el caso de T 0 conocido

E Mdulo de elasticidad de Young

F Matriz de flexibilidad del elemento cable