ESCUELA PROFESIONAL DE INGENIERIA CIVIL

DOCENTE: RODRIGUEZ PLACENCIA

CURSO: ANALISI ESTRUCTURAL II

TEMA: ANALISIS MATRICIAL

ALUMNO: OLASCOAGA LLICAN RICARDO ALFREDO

CAJAMARCA- 2019
 RESUMEN

El proyecto que se presenta en este documento, tiene como finalidad el estudio del
comportamiento en estado estático de una estructura reticular.
El proyecto está dividido en tres grandes bloques. El primero se plantean las
ecuaciones que gobierna el método matricial de la rigidez, un método de cálculo
aplicable a estructuras hiperestáticas y reticuladas de barras que se comportan de
forma elástica y lineal.
Seguidamente en el segundo bloque se explicará la interacción necesaria para
resolver y comprender cómo se comporta la estructura estudiada en verso a una
simulación de rigidez torsional y a una máxima solicitación de rigidez en circuito y
así poder dar paso al tercer bloque.
En este último bloque y posiblemente el más importante se implementa un modelo de
optimización capaz de recalcular todas las variables y posibles configuraciones en
busca de un mínimo peso y una resistencia torsional adecuada para nuestra estructura
reticulada, en este bloque se estudiarán los resultados obtenidos para que el lector
tenga un mejor entendimiento.
La misión de este trabajo es la de unificar todas las modelizaciones matemáticas de la
realidad física
 INTRODUCCION AL CALCULO MATRICIAL
Generalidades Los métodos de cálculo matricial (CM) de estructuras son un conjunto
de métodos que tienen en común organizar toda la información en forma de matrices.
En estos métodos, todas las relaciones entre las distintas partes de una estructura dan
lugar a sistemas de ecuaciones con un alto número de variables, pero donde no se
han realizado suposiciones o simplificaciones en las que se pierda información
relevante. Esta generalidad, junto a la estructura de la información en matrices,
permite que su planteamiento y resolución pueda ser ejecutada de manera automática
por medio de programas de ordenador, lo que ha hecho que en la actualidad sean la
práctica habitual en la ingeniería.
la rigidez de cálculo matricial, aplicado a estructuras bidimensionales formadas por
barras y vigas. Este mismo esquema puede ser extendido a otras formas de discretizar
una estructura o un medio continuo. De hecho, el método de los Elementos Finitos es
la extensión del método de CM donde se trata con elementos que no son solo barras,
sino volúmenes de distintas formas geométricas que modelan un mayor número de
problemas mecánicos o físicos. En todo el desarrollo del método aceptaremos las
hipótesis generales en las que normalmente se desarrolla la Teoría de Estructuras,
esto es, comportamiento elástico y lineal del material y estado de pequeños
desplazamientos.
Métodos de cálculo matricial
En términos generales, existen dos procedimientos genéricos en mecánica de medios
continuos de solidos deformables para poder establecer el sistema completo de
ecuaciones dependiendo del orden en que las vayamos aplicando. Las ecuaciones
que podemos poner en juego son las ecuaciones de equilibrio, las de comportamiento
y las de compatibilidad del problema. Cuando partiendo de las ecuaciones de equilibrio
las utilizamos para incorporarlas a las de comportamiento y finalmente el resultado lo
introducimos en las ecuaciones de compatibilidad, estamos aplicando el método
denominado de la compatibilidad o de la flexibilidad. Hablando en términos de las
variables implicadas, en este caso llegamos a formular los desplazamientos en función
de las cargas aplicadas. Si seguimos el procedimiento inverso, inicialmente
relacionamos deformaciones y desplazamientos aplicando las ecuaciones de
compatibilidad para posteriormente aplicar las leyes de comportamiento y finalmente
las ecuaciones de equilibrio, en ese caso el método se denomina de la rigidez o del
equilibrio. En la Figura 1.4.1 se esquematiza brevemente este proceso. En cálculo
matricial, tal y como se verá a continuación en un ejemplo, es posible aplicar ambos
procedimientos. Sin embargo, tal y como se desarrollará en los capítulos siguientes,
únicamente es posible llegar a un procedimiento automático y sistematizado con el
método de la rigidez, siendo este por tanto el que se ha implantado y generalizado.
FIGURA: Los dos caminos alternativos para plantear las ecuaciones que modelan una
estructura: el método de la flexibilidad y el de la rigidez. Las funciones G, Φ, Ψ y G′ , Φ
′ , Ψ′ son nombres simbólicos de funciones diferentes cuya forma exacta no es
relevante aquí.
Ejemplo A continuación se va a desarrollar un ejemplo donde se puede ver la aplicación de
ambos métodos a un mismo problema y como ´este se puede estructurar en forma
matricial. El problema a estudiar es un problema clásico de resistencia de materiales,
se trata de una viga en voladizo como la mostrada en la Figura 1.5.1. El primer paso
que habrá que dar es la discretización del problema. Tenemos diversas opciones de
discretización, todas ellas válidas. Algunas opciones se muestran en la parte derecha
de la Figura 1.5.1. Se puede ver que incluso el número de gdl es variable. A mayor
número de gdl, mayor detalle de resultados a cambio de un problema mas voluminoso
y con mayor tiempo de cálculo.

Figura: (Izquierda) Ejemplo de problema para estudio por los métodos matriciales de flexibilidad y
de la rigidez. (Derecha) Algunas de las opciones posibles para la discretización del
problema, donde se indican los gdl elegidos para que entren en el modelo del sistema. al
desplazamiento horizontal y vertical del extremo y a su giro. Esta elección nos limita el tipo
de problemas que podremos resolver con esta discretización, por ejemplo, las cargas solo
podrán estar aplicadas en ese extremo. Una vez discretizado el problema, empezamos
obteniendo su correspondiente matriz de rigidez K, que relaciona las solicitaciones F y
desplazamientos U de manera que:
Figura: Ejemplo de cálculo de matriz de rigidez mediante problemas unitarios. Los
desplazamientos y cargas positivos se realizan en las direcciones +X y +Y, siendo los giros
y momentos positivos en sentido contrario a las agujas del reloj. donde podemos obtener
los valores de la matriz (los números Kij ) por el método de inspeccionar cada uno de los
problemas unitarios, uno por gdl. Dicho método consiste en plantear los problemas que
corresponden con un vector de desplazamiento nulo en todos los gdl menos en uno (donde
habrá un desplazamiento unitario) y calcular las solicitaciones asociadas a dicho
desplazamiento. Los valores de esfuerzos así obtenidos se colocan en la columna
correspondiente al gdl en el que se aplicó el desplazamiento, y así, columna a columna, se
puede obtener la matriz completa.
 Juntando estos tres últimos resultados obtenemos la matriz de rigidez K completa:

Si en cambio hubiésemos elegido plantear el problema mediante la matriz de flexibilidad A,

que relaciona esfuerzos y desplazamientos mediante U� = AF� el método sería similar,
planteando en este caso los distintos problemas de cargas unitarios, calculando los
desplazamientos asociados a cada uno de ellos y colocando dichos desplazamientos en cada
una de las columnas de la matriz de flexibilidad. El proceso se ilustra.

Figura: Ejemplo de cálculo de matriz de flexibilidad mediante problemas unitarios. Los

desplazamientos y fuerzas positivos se realizan en las direcciones +X y +Y, siendo los giros y
momentos positivos en sentido contrario a las agujas del reloj. en la Figura 1.5.3, y permite
crear la matriz columna a columna de esta manera:
 Una vez hayamos obtenido la matriz de rigidez o de flexibilidad de una estructura, tan simple
como la de este ejemplo o arbitrariamente compleja, podemos resolver cualquier problema
que combine gdl con desplazamientos conocidos (apoyos, asientos, etc.…) con gdl libres
cuyos desplazamientos sean incógnitas. Dependiendo de las variables podrá ser mejor aplicar
un método u otro, aunque en general los gdl conocidos y las incógnitas estarán intercalados y
por tanto siempre habrá que reorganizar la matriz para resolverla, no habiendo ventajas en
este sentido entre los métodos de la rigidez y de la flexibilidad.
Para ilustrar como se obtendría la matriz de rigidez para una estructura más compleja,
consideremos ahora el caso del pórtico de la Figura.

Figura: Ejemplo de un pórtico. Se estudiarán los gdl correspondientes a los cuatro nudos numerados
en la figura. De manera similar a los problemas de desplazamientos unitarios de la Figura
1.5.2, para este pórtico analizaríamos los 11 problemas mostrados en la Figura 1.5.5. Para
cada uno de dichas situaciones, calcularíamos las fuerzas necesarias para obtener dichos
desplazamientos y como vimos.
Características de la matriz de rigidez.

Antes de finalizar este capítulo introductorio conviene resumir aquí algunas de las principales
características de las matrices de rigidez K:

 La matriz de rigidez es una propiedad del sistema estructural, no cambia en función del estado de
cargas o de condiciones de contorno a que se someta a la estructura. Solo se verá afectada si se
introduce algún elemento adicional.
 Cada columna representa las acciones necesarias para conseguir un desplazamiento unitario en el
grado de libertad definido por el ´índice de la columna a la vez que se quedan fijados a cero el resto de
los gdl.
 Una fila es un conjunto de multiplicadores que operados sobre el vector desplazamiento completo
proporcionan el valor de la fuerza correspondiente al gdl definido por el ´índice de la fila.
 Cada termino kij se puede considerar una “función de peso” que representa la proporción de
contribución a la fuerza del gdl i debido al desplazamiento del gdl j. En caso de que su valor sea cero
significa que ambos gdl no está relacionados.

Matriz completa de una barra (6 gdl)

Matriz de rigidez
Consideremos ahora el caso más general de una barra a en la que estudiemos sus 6 grados de
libertad (el máximo posible en 2D). Por ejemplo, se deberá emplear dicho modelo con barras en
la configuración llamada “bien potrada”, es decir, con uniones rígidas en sus dos nudos extremos
i y j. Al igual que en el caso anterior, la barra podrá estar en una posición y orientación arbitraria
dentro de un sistema de coordenadas global del problema < x, y >, en el que definimos un
sistema local de coordenadas de la barra < x, ˆ y >ˆ tal que ˆx esta alineado con la dirección de la
barra i → j, como se ve en Figura.

Figura: Sistema de coordenadas local a una barra cuyos dos extremos están empotrados. Las
solicitaciones en los extremos de dicha barra constan ahora de seis escalares: las dos
componentes (x e y) de las fuerzas en cada extremo más el momento flector, también para cada
extremo.

Los seis componentes, en coordenadas locales de la barra, se denotarán por el vector � ˆ f a:

Hay que resaltar que en el caso de los momentos ma i y ma j no se ha usado la notación de
coordenadas locales (el símbolo ˆ·) al ser magnitudes invariantes con el sistema de coordenadas
empleado para estructuras planas. Los desplazamientos de los extremos de la barra en este caso
también tienen seis componentes: los desplazamientos propiamente dichos (direcciones x e y) y
los giros θ de las secciones extremas (ver Figura 2.3.2). Todos ellos se denotan mediante el
vector �una, también en coordenadas locales:

Se puede demostrar que, para pequeños desplazamientos, dichos desplazamientos �una y las
solicitaciones � ˆ fa que causan están relacionados linealmente mediante la matriz de rigidez Kˆ
a, de tamaño 6×6, formada a su vez por 4 su matriz de 3 × 3, de forma que:

Al igual que en la sección anterior, podemos expandir el producto de matrices, llegándose a una expresión
idéntico a la Ec. 2.2.5. Al igual que para el caso de las matrices articuladas en §2.2.1, se demuestra,
invocando de nuevo el teorema de Maxwell-Betti, que para cualquier par de nudos i y j las su matrices Kˆ a
ij y Kˆ a ji son traspuestas: Kˆ a ij = Kˆ a ji ⊤. Es interesante resaltar que en el caso de barras articuladas las
sus matrices correspondientes eran idénticas en lugar de traspuestas, pero sólo porque en ese caso las
submatrices son simétricas en sí mismas, y, por lo tanto, iguales a su traspuesta. De igual manera que se
hizo para una barra biarticulada, mostramos en la Figura 2.3.3 los seis problemas de desplazamientos y
giros unitarios que permiten derivar el valor de la matriz de rigidez a partir de expresiones conocidas de
resistencia de materiales. De esta forma, obtenemos las 6 ecuaciones (para coordenadas locales), que
gobiernan una barra bien potrada.

donde I y A son el momento de inercia y el ´área de su sección transversal, respectivamente, L es

la longitud de la barra y E el módulo de elasticidad del material.
Cambio de coordenadas

La relación establecida arriba asume que tanto las solicitaciones como los desplazamientos están dados en
coordenadas locales a la barra, por lo que habrá que transformarlos a coordenadas globales del problema
para poder considerar el conjunto de la estructura. Las fuerzas y los desplazamientos en coordenadas
globales serán ahora denotados como los vectores f a y una, respectivamente, como se muestra en la
Figura 2.3.4. Para el caso de una barra con sus dos nudos rígidos, definimos la siguiente matriz de
transformación de coordenadas:
 Figura: Los distintos desplazamientos unitarios (izquierda) en una barra bien potrada y las
correspondientes solicitaciones asociadas (derecha).

ANEXOS.

http://ingmec.ual.es/~jlblanco/papers/blanco2012calculo_matricial_estructuras.pdf

https://upcommons.upc.edu/bitstream/handle/2117/105804/VOLUMEN%20I%20MEMORIA.pdf?
sequence=3&isAllowed=y