PROYECTO FINAL

APLICACIÓN DE METODOS NUMERICOS EN LA INGENIERIA CIVIL

PRESENTADOR POR:
JHON JAMES LASTRA PRADO

PRESENTADO AL DOCENTE:

OSWALDO ROJAS CAMACHO

UNIVERSIDAD NACIONAL DE COLOMBIA

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA CIVIL
BOGOTA D.C.
 1. PROBLEMA
Aplicar métodos numéricos para resolver una estructura mediante el análisis no lineal de la
deformación de estructuras para conocer el desplazamiento total en un punto específico sobre la
longitud del mismo mediante el uso del método Newton- Raphson para ecuaciones no lineales y
compararlo con otro método de solución con análisis estructural y así demostrar la importancia
de los métodos numérico.
2. MARCO TEÓRICO
Todo elemento estructural sometido a la acción de cargas se deforma. En la mayoría de los casos
las deformaciones son imperceptibles, y deben ser medidas con el auxilio de instrumentos.
Podemos definir la deformación de un elemento estructural como una alteración del estado físico
debido a una fuerza mecánica externa, a una variación de temperatura, a un cedimiento de
apoyos, etc.
La deformación puede ser elástica, cuando desaparece al cesar la acción que lo produce, por lo
que las partículas elementales del cuerpo vuelven a su posición inicial; y permanente, cuando
persiste después de desaparecer la causa que lo ha producido. El elemento en el primer caso tiene
un comportamiento elástico y en el segundo un comportamiento plástico o anelástico. (Morales,
2013)
2.1. ANÁLISIS NO LINEAL PARA DEFORMACIONES EN ESTRUCTURAS
Cuando se hace un estudio de deformaciones de un material es común encontrarse con el término
“análisis lineal” el cual hace referencia a el análisis de deformaciones que ocurren dentro del
rango elástico del material, las cuales generalmente constan de una magnitud pequeña, se
denomina “análisis lineal” debido a que estas deformaciones abarcadas dentro del rango elástico
presentan comportamiento lineal, este comportamiento lineal facilita en gran medida el estudio
de estas deformaciones ya que muchas características fundamentales del material como lo son
por ejemplo su módulo de elasticidad presentan valores constantes. Pero este tipo de análisis solo
aplica para la región elástica del material en el cual ocurren deformaciones de escaza magnitud, a
pesar de que al momento de efectuar el diseño o análisis de una estructura se busca que las
deformaciones tengan la menor magnitud posible, pero existen ciertos tipos de estructuras que
debido a su naturaleza están expuestos a desplazamientos mayores, lo cual genera deformaciones
que no pueden ser analizadas linealmente.
A continuación se presenta un ejemplo de una curva de esfuerzo contra deformación de un material
elástico, la cual muestra el comportamiento de cualquier material ante la aplicación de una fuerza externa
que le cause una deformación en su composición física.
 Gráfica esfuerzo vs deformación.
Un ejemplo de este tipo de análisis podría ser el de las deformaciones que se presentan en las estructuras
que contienen cables, como es el caso de los puentes colgantes que no presentan vigas de rigidez, estas
estructuras debido a su naturaleza carecen de una linealidad geométrica gracias a que los cables de la
estructura presentan deformaciones tanto en su longitud como en su rotación y curvatura, para efectos
prácticos la deformación en la curvatura se puede despreciar, ya que presenta magnitudes grandes. Por
otra parte se tiene que este tipo de materiales presentan grandes deformaciones longitudinales y
rotacionales, que impiden que las ecuaciones de equilibrio de la estructura se planteen en la geometría
lineal como pasa en las vigas con apoyos simple, ubicándose de esta manera en la geometría deformada
que ocasiona que las ecuaciones de deformación y desplazamiento contengan términos no lineales que
requieran una solución medianamente compleja para hallar así los valores requeridos. A continuación se
muestra una imagen con un ejemplo de construcción real de este tipo de estructuras en la ingeniería civil:
De acuerdo a lo explicado anteriormente, una de las situaciones que permiten comprender el
grado de dificultad que puede llegar a tener un análisis no lineal de las deformaciones consiste en
la estructura mostrada en la 2, que consiste en una estructura conformada por los cables que se
encuentran sosteniendo un peso muerto en la longitud L/2 tal y como se muestra a continuación:

L/2 L/2

α1
α2 Li sin α1
Li
x
Lf P
𝜐
T T
P

A partir de esta estructura presentada anteriormente se realizará el desarrollo del planteamiento de una
aplicación de los métodos numéricos en la ingeniería civil, el cual consiste en encontrar el desplazamiento
máximo en el presentado en 𝐿/2. Para evaluar la estructura se deben usar las ecuaciones de equilibrio de la
siguiente manera:
2 T sin α 2=P
T: tensión de los cables.
P: Peso muerto
De donde se obtiene que T es igual a:
P Lf P
T= = (1)
2∗sin α 2 2∗¿ ¿
Lf : Longitud
Li : Longitud
Donde Lf se puede expresar como:
T Li
Lf =Li + Según las relaciones constitutivas (2)
EA
Lf =√ (Li sin α 1 + v ¿ ¿ )2 +( Li cos α 2)2 Según las relaciones de compactibilidad (3)

E: Modulo de elasticidad
A: Área
De las ecuaciones 1, 2,3 se obtiene la siguiente expresión:
2 EAx 2 EAx
− 2 =P
Li √
x + ( L cos α )
2
i 2

x:
Donde x se expresa como:
x=Li sin α 1 + v
Como se puede evidenciar, las fórmulas obtenidas mediante el análisis no lineal no son de fácil solución
mediante los métodos convencionales, debido a esto la implementaciones de métodos numéricos es de
vital importancia para la solución de este tipo de problemas, para el caso expuesto anteriormente se
optó por el método Newton – Raphson para encontrar la raíz de ecuaciones no lineales.

2.2. MÉTODOS NÚMERICOS PARA ECUACIONES NO LINEALES

Los métodos numéricos para determinar configuraciones de equilibrio de una estructura con
comportamiento no lineal requiere la capacidad de resolver sistemas de ecuaciones no lineales.
Por ende, los procedimientos de análisis no lineal de estructuras están directamente basados en
métodos numéricos de resolución de sistemas de ecuaciones no lineales. En varios casos, las
aplicaciones prácticas de resolución de estructuras han sido las precursoras de los que
posteriormente serían métodos numéricos aplicados a diversas ramas de la Ingeniería, como el
método de elementos finitos, pero de acuerdo al curso de métodos numéricos para ingeniería se
usará el método de Newton- Raphson a ver cómo se comportan en estos casos de análisis
estructural no lineal.
El desarrollo de ecuaciones no lineales tiene muchos métodos con procedimientos distintos para el
momento de abarcarse en cuanto a la búsqueda de una solución, sin embargo, el caso general más usado
para encontrar raíces es por medio de Newton – Raphson que es un método iterativo que se basa en la
repetición de acciones análogas.
Este método permite hallar raíces de funciones, conocido un valor numérico cercano a la raíz ( x i=¿ x ¿).0

Es un método abierto e iterativo, en general de rápida convergencia, muy útil para el cálculo de raíces
cuadradas y de mayor grado.
'
Sabemos por el teorema de Taylor que para un x 0 ∈ (𝑎, 𝑏) tal que F (x 0 ) es diferente de cero y además si
se cumple que tanto la función como la derivada son continuas en el intervalo (a, b), entonces;
'
f ( x )=f ( x 0 ) + ( x 1−x 0 ) f ( x0 ) + …
Como se busca el punto donde f ( x 1 ) =0, entonces:
0=f ( x )=f ( x0 ) + ( x1 −x 0 ) f ' ( x 0 ) + …
Por lo tanto se puede concluir que;
f ( x0 ) f ( xn )
x 1 ≈ x 0− ' =¿ x n+1 ≈ xn − ' En forma general.
f ( x0) f ( xn)
Al fijar el valor 𝜀𝑦 (error permitido en “y”) y el valor 𝜀𝑥 (error permitido en “x”) la función (𝑥𝑖+1) y los
puntos (𝑥𝑖+1) tienen que cumplir las siguientes condiciones:
| 𝒇(𝒙𝒊+𝟏)| ≤ 𝜺𝒚 𝒚 | (𝒙𝒊+𝟏 − 𝒙𝒊)| ≤ 𝜺𝒙
 Descripción grafica del método Newton-Raphson
2.2.1. Deducción de la formula
La fórmula se deduce con base a la función tangente, dado lo siguiente :
a '
tanθ= Asumiendo que : tan θ=f ( x i ) a=f ( x i ) ; b=x i−x i +1
b
f ( x i)
f ' ( x i )=
x i−x i+1

f ( xi )
x i+1=x i− '
f ( xi )
El método Newton – Raphson es compatible para hallar raíces complejas, teniendo en cuenta que el valor
inicial 𝑥0 en el camp complejo siempre y cuando la parte compleja e imaginaria difiera de “0”.

2
 3. Resultados

Para el desarrollo del ejercicio planteado se tienen los siguientes datos de entrada:

Y a partir de esto se calcula lo siguiente:

Y se plantean las siguientes ecuaciones:

Con lo que al reemplazar los datos de entrada se tendrían las siguientes ecuaciones para resolver por el
método de Newton – Raphson y así encontrar el x necesario para la determinación del desplazamiento:

Que aplicando el método de Newton – Raphson mediante el algoritmo creado en el software Matlab que
se encuentra anexo en este documento.
Teniendo como valores de entrada un valor inicial de 0,5, un número máximo de iteraciones de 20 y una
tolerancia de 0,00001 que arrojan que en 5 iteraciones se encuentra el valor de convergencia igual a
0,61497, es decir, que el valor encontrado para la x en la ecuación no lineal planteada para la resolución
de la estructura es de x = 0,61497 m que es el valore con el que se procedería a encontrar el
desplazamiento del nudo inferior del cable de la siguiente forma:

𝝊 = 𝐱 − 𝑳𝒊 𝐬𝐢𝐧 𝜶𝟏 = 𝟎, 𝟔𝟏𝟒𝟗𝟕 − 𝟎, 𝟒𝟑𝟕𝟒𝟒 = 𝟎, 𝟏𝟕𝟕𝟓𝟑 𝐦 ↓

Se tiene un valor de desplazamiento total del nudo es de 0,17753 metros, que se proceden a comprobar
mediante un modelado con ayuda del software de análisis estructural SAP2000 para obtener una
estructura de la siguiente forma:

Que al cargarla y hacer trabajar el modelo quedaría una estructura que presenta una deformación como la
mostrada a continuación.
Donde se tiene que al revisar los resultado de la deflexión total en el punto medio de la estructura se
encuentra un valor de 0,172 metros, tal y como se muestra a continuación:

El valor exacto de la deflexión que se encontró en la imagen anterior es muy cercano al obtenido
mediante el método de Newton – Raphson que fue de 0,177 metros, lo que comprueba el valor de los
métodos numéricos en cualquier ámbito de la ingeniería.
 4. Análisis de los resultados

El problema del desplazamiento vertical en elementos cargados axialmente y trabajados a tensión es una
situación que se presenta a menudo en el desarrollo de proyectos de ingeniería civil, como en puentes
colgantes o es cables eléctricos, por citar algunos ejemplos.

En consecuencia a esto el ejercicio estudiado anteriormente muestra las ventajas de tener a la mano
métodos iterativos que permitan resolver ecuaciones no lineales, ya que, problemas de desplazamiento
verticales desarrollados mediante el método de compatibilidad de deformaciones representan una gran
complejidad al momento de solucionar la ecuación ahí descrita.

Los datos de entrada fueron los siguientes:

Se debe tener en cuenta que el elemento es de sección circular y de Acero, por lo que el diámetro a usar es
de 0,035 m y el Módulo de elasticidad es de 16100000 Tn/m2. Para la solución del problema
analíticamente se hizo uso de ecuaciones mostradas anteriormente de las cuales en su mayoría son
sacadas haciendo uso de la geometría de la estructura estudiada y las restantes de bibliografía empleada.
El desplazamiento en el punto central, analíticamente calculado fue de 0,1775 m, es decir, 17,75 cm, un
valor bastante razonable a los datos iniciales o, de entrada. Si bien, el desplazamiento teórico o usado para
el diseño en Ing. Civil se calcula haciendo uso de software especializados, para nuestro caso SAP 2000,
no hay que dejar de lado que el método aproximado de newton Raphson desarrollado en Matlab para
calcular este valor X se realizó con un error de convergencia de 0,00001, por lo que, los valores ahí
encontrados se pueden considerar aceptables. Para realizar un análisis más profundo, fue necesario el
desarrollo de la estructura en SAP 2000, en donde se halló que el desplazamiento en medio de la luz es de
0,172606 m, un valor muy parecido al encontrado analíticamente.

Si calculamos un porcentaje de error entre ambas respuestas tenemos lo siguiente.

Por lo anterior se puede concluir que la implementación de métodos numéricos ayuda a aproximar las
respuestas de los problemas a tal grado que los errores entre lo esperado y lo que se tiene es menor al 3%,
además a esto, herramientas como Matlab ayudan que estos métodos se pueden optimizar de tal manera
que la solución a ecuaciones no lineales, por detallar un ejemplo, se puede realizar de manera instantánea.