Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Teoria de Matrices

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 19

Republica Bolivariana de Venezuela

Ministerio del Poder Popular para la Educacin


I.U.P SANTIAGO MARIO
Valencia-EDO-Carabobo

Integrantes:
Mayrene Anceno
C.I: 20163815

Valencia Octubre de 2016

INTRODUCCION

La Ingeniera Civil es una especialidad que nos adhiere a la


sociedad mediante el diseo y ejecucin de obras, y en el
proceso de sus acciones tambin hace uso de las matrices ya
que se utilizan para el diseo de sistemas estructurales en las
diversas reas que ocupa la Ingeniera Civil. Las matrices sirven
para resolver sistemas de ecuaciones lineales, estos a su vez
tienen mltiples aplicaciones en el rea de ingeniera dando
lugar a al ptimo manejo de recursos humanos y de materiales
monitoreados y controlados desde un sistema de diseo dando
as paso a la modernidad y a la ingeniera del futuro.
Es importante las matrices ya que nos llevan a conocer un sin fin
de mtodos que nos facilitan algunos problemas matemticos,
ya que las matices representan de forma implcita una particular
relacin evolutiva. La eleccin de una matriz determinada puede
afectar enormemente al resultado del anlisis, y por lo cual es
necesario saber utilizarlas Las matrices se usan en cualquier
comparacin de secuencias. Tambin se utilizan matrices de
sustitucin para incrementar la sensibilidad en los alineamientos
dbiles. Estas matrices se basan en observaciones.
Se utiliza para varias ramas de la ingeniera, para resolver
problemas que se encuentran en muchas dimensiones cuando
se tienen problemas que solo se pueden resolver con sistemas
de ecuaciones diferenciales se arman matrices con dichas
ecuaciones de tal manera que se pueda solucionar ese
problema, ejemplos prcticos enfrentan los ingenieros civiles
que teniendo sistemas de ecuaciones con muchas ecuaciones
las resuelven por mtodos matriciales.

Teora de Matrices
La teora de matrices es una rama de las matemticas que se centra en el estudio de
matrices. Inicialmente una rama secundaria del lgebra lineal, ha venido cubriendo tambin
los temas relacionados con la teora de grafos, el lgebra, la combinatoria y la estadstica.
Se denomina matriz a todo conjunto de nmeros o expresiones dispuestos en forma
rectangular, formando filas y columnas.

Las matrices se utilizan para mltiples aplicaciones y sirven, en particular, para representar
los coeficientes de los sistemas de ecuaciones lineales o para representar transformaciones
lineales dada una base. En este ltimo caso, las matrices desempean el mismo papel que
los datos de un vector para las aplicaciones lineales.
Pueden sumarse, multiplicarse y descomponerse de varias formas, lo que tambin las hace
un concepto clave en el campo del lgebra lineal.
Elemento de una matriz
Cada uno de los nmeros de que consta la matriz se denomina elemento.
Un elemento se distingue de otro por la posicin que ocupa, es decir, la fila y la columna a
la que pertenece.
Dimensin de una matriz
El nmero de filas y columnas de una matriz se denomina dimensin de una matriz. As,
una matriz de dimensin min es una matriz que tiene m filas y n columnas.
De este modo, una matriz puede ser de dimensin: 2x4 (2 filas y 4 columnas), 3x2 (3 filas y
2 columnas), 2x5 (2 filas y 5 columnas),...

S la matriz tiene el mismo nmero de filas que de columnas, se dice que es de orden: 2, 3,
4,...
El conjunto de matrices de m filas y n columnas se denota por Amxn o (aij).
Un elemento cualquiera de la misma, que se encuentra en la fila i y en la columna j, se
denota por aij.
Matrices iguales
Dos matrices son iguales cuando tienen la misma dimensin y los elementos que ocupan el
mismo lugar en ambas, son iguales.
Calculo con matrices
En trminos generales, existen dos procedimientos genricos en mecnica de medios
continuos de slidos 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 mtodo denominado de la
compatibilidad o de la flexibilidad. Hablando en trminos de las variables implicadas, en
este caso llegamos a formular los desplazamientos en funcin 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 mtodo
se denomina de la rigidez o del equilibrio. En la Figura 1.4.1 se esquematiza brevemente
este proceso
Suma de matrices
Dadas dos matrices de la misma dimensin, A= (aij) y B= (bij), se define la matriz suma
como: A+B= (aij+bij).
La matriz suma se obtiene sumando los elementos de las dos matrices que ocupan la misma
misma posicin.

Producto de un escalar por un matriz


Dada una matriz A= (aij) y un nmero real k R, se define el producto de un nmero real
por una matriz: a la matriz del mismo orden que A, en la que cada elemento est
multiplicado por k.
KA= (k aij)

Producto de matrices
Dos matrices A y B son multiplicables si el nmero de columnas de A coincide con
el nmero de filas de B.
Mm x n x Mn x p = M m x p
El elemento cij de la matriz producto se obtiene multiplicando cada elemento de la fila ide
la matriz A por cada elemento de la columna j de la matriz B y sumndolos.

Matriz inversa
El producto de una matriz por su inversa es igual a la matriz identidad.
A A-1 = A-1 A = I
Clculo de la matriz inversa por el mtodo de Gauss
Sea A una matriz cuadrada de orden n. Para calcular la matriz inversa de A, que
denotaremos como A-1, seguiremos los siguientes pasos:
1 Construir una matriz del tipo M = (A | I), es decir, A est en la mitad izquierda de M y la
matriz identidad I en la derecha.
Consideremos una matriz 3x3 arbitraria

La ampliamos con la matriz identidad de orden 3.

2 Utilizando el mtodo Gauss vamos a transformar la mitad izquierda, A, en la matriz


identidad, que ahora est a la derecha, y la matriz que resulte en el lado derecho ser la
matriz inversa: A-1.
F2 - F1

F3 + F2

F2 - F3

F1 + F2

(-1) F2

La matriz inversa es:

Clculo de la matriz inversa por determinantes


Rango de una matriz: es el nmero de lneas de esa matriz (filas o columnas) que son
linealmente independientes.
Una lnea es linealmente dependiente de otra u otras cuando se puede establecer una
combinacin lineal entre ellas.

Una lnea es linealmente independiente de otra u otras cuando no se puede establecer una
combinacin lineal entre ellas.
El rango de una matriz A ser simboliza: rang(A) o r(A).
Tambin podemos decir que el rango es: el orden de la mayor submatriz cuadrada no
nula. Utilizando esta definicin se puede calcular el rango usando determinantes.
Se puede calcular el rango de una matriz por dos mtodos:
Clculo del rango de una matriz por el mtodo de Gauss
Podemos descartar una lnea si:
Todos sus coeficientes son ceros.
Hay dos lneas iguales.
Una lnea es proporcional a otra.
Una lnea es combinacin lineal de otras.

F3 = 2F1
F4 es nula
F5 = 2F2 + F1
R(A) = 2.
En general consiste en hacer nulas el mximo nmero de lneas posible, y el rango ser el
nmero de filas no nulas.

F2 = F2 - 3F1
F3 = F3 - 2F1

Por tanto r(A) = 3.


Matriz de rigidez
Sea una barra, a la que denominaremos a, situada entre dos nudos i y j mediante uniones
articuladas, en una posicin arbitraria con respecto al sistema de coordenadas global del
problema < x, y >, tal y como se ilustra en la Figura 2.2.1, y a la que tambin se le asocia
un sistema de coordenadas local < x, y > tal que x est alineado con la direccin de la
barra i j. El vector de todas las fuerzas sobre dicha barra articulada contiene cuatro
escalares: las dos componentes (x e y) para cada extremo i y j de una barra. Este vector,
donde las fuerzas estn expresadas en coordenadas locales a la barra (ver Figura 2.2.2), se
denota por:

Sobre la notacin empleada, remarcar que un vector o componente cualquiera relacionado


con una barra a se escribir como a para referirse a sus coordenadas globales del
problema, mientras que en coordenadas locales se escribir a, como en la ecuacin de
arriba. En cuanto a los 4 grados de libertad (gdl) de desplazamientos, que asumimos seran
pequenos para que el mtodo de clculo nos d una aproximacin adecuada, los
denotaremos por el siguiente vector en coordenadas locales

Una vez definidas las fuerzas y desplazamientos que sufre una barra a (fa y ua,
respectivamente), se puede demostrar fcilmente que ambos vectores estn relacionados
linealmente entre s mediante una matriz de rigidez K a, de tamao 4 4, que a su vez se
compone de 4 submatrices de 2 2, de forma que:

El significado de cada una de estas submatrices queda claro si expandimos los vectores de
fuerzas y desplazamientos en sus dos elementos (para cada uno de los nudos de la barra i y
j):

Est claro ahora que cada submatriz K a establece la relacin entre los desplazamientos
en el nudo y las solicitaciones provocadas por estos en el nudo . Observando estas
ecuaciones, y por el principio de simetra en las acciones y reacciones, se puede demostrar
(ver teorema de Maxwell-Betti de resistencia de materiales) que para cualquier par i y j, las
dos submatrices involucradas son la transpuesta una de otra, es decir: K a ij = K a ji .
En el caso particular de barras articuladas es aun ms sencillo, ya que K a ij = K a ji
debido a que todos los elementos menos uno son ceros. Por lo tanto, la matriz de rigidez de
una barra particular, K a (de 4 4), se compone nicamente de dos submatrices de 2 2
distintas, dos matrices K a en las dos diagonales y una misma matriz K a en las dos
esquinas, de forma que la matriz final tambin es simtrica. La matriz de rigidez de una

barra, o ms en general de cualquier elemento o estructura, tiene una interpretacin fsica


muy ilustrativa. Para verlo, reescribiremos la Ec. 2.2.3 reemplazando la matriz de rigidez
K a por sus elementos genricos ka ij , es decir:

Entonces, se puede ver que la matriz de rigidez se puede en realidad definir como aquella
matriz cuyas columnas representan las solicitaciones asociadas a desplazamientos (y giros)
unitarios. Como ejemplo, tomemos el vector de desplazamiento unitario correspondiente al
eje x del nodo i, es decir ua = (1 0 0 0) . Usando la anterior expresin matricial
obtenemos:

Se verifica que la primera columna de la matriz define las solicitaciones asociadas a un


desplazamiento unitario en el primer grado de libertad (uix en este caso), y en general, la
nesima columna se corresponde con el nesimo grado de libertad. Es importante tener en
mente que aunque hablemos de desplazamientos y giros unitarios, en realidad se est
modelando en todo momento segn una aproximacin de pequenos desplazamientos, a
pesar de lo grandes que realmente sean los desplazamientos y giros unitarios (que si se usa
el SI, seran metros y radianes, respectivamente). Como se ilustra en la Figura 2.2.3,
plantear todos los problemas unitarios permite deducir los valores de cada uno de los
coeficientes de la matriz de rigidez. Este es, de hecho, una de las formas ms comunes de
derivar dicha matriz, y en nuestro caso concreto llegamos a:

Donde L es la longitud de la barra, E el modulo de elasticidad del material y A el rea de su


seccin transversal. Es trivial ver aqu que las solicitaciones en las direcciones
perpendiculares a la barra en ambos extremos (f a iy y f a jy, respectivamente) seran
nulas, es decir: en una barra articulada en la que solo se aplican cargas en sus extremos
(como se est asumiendo hasta ahora) solamente existirn fuerzas axiles, que darn lugar a
traccin o compresin pero nunca a esfuerzos cortantes. En el caso de que existan cargas
distribuidas sobre la barra s que aparecern esfuerzos tanto axiles como cortantes, como:

Caractersticas de la matriz de rigidez


Antes de finalizar este captulo introductorio conviene resumir aqu algunas de las
principales caractersticas de las matrices de rigidez K: La matriz de rigidez es una
propiedad del sistema estructural, no cambia en funcin del estado de cargas o de
condiciones de contorno a que se someta a la estructura. Solo se ver afectada si se
introduce algn 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 funcin de peso que representa la proporcin de contribucin 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 relacionado
Otro punto del Mtodo matricial de la rigidez.
El mtodo matricial de la rigidez es un mtodo de clculo aplicable a estructuras
hiperestticas de barras que se comportan de forma elstica y lineal. En ingls se le
denomina direct stiffness method (DSM, mtodo directo de la rigidez), aunque tambin se
le denomina el mtodo de los desplazamientos. Este mtodo est diseado para realizar
anlisis computarizado de cualquier estructura incluyendo a estructuras estticamente
indeterminadas.
El mtodo matricial se basa en estimar los componentes de las relaciones de rigidez para
resolver las fuerzas o los desplazamientos mediante un ordenador. El mtodo de rigidez
directa es la implementacin ms comn del mtodo de los elementos finitos. Las
propiedades de rigidez del material son compilados en una nica ecuacin matricial que
gobierna el comportamiento interno de la estructura idealizada.
Los datos que se desconocen de la estructura son las fuerzas y los desplazamientos que
pueden ser determinados resolviendo esta ecuacin. El mtodo directo de la rigidez es el
ms comn en los programas de clculo de estructuras (tanto comerciales como de fuente
libre).
El mtodo directo de la rigidez se origin en el campo de la aeronutica. Los investigadores
consiguieron aproximar el comportamiento estructura de las partes de un avin mediante
ecuaciones simples pero que requeran grandes tiempos de clculo. Con la llegada de los
ordenadores estas ecuaciones se empezaron a resolver de forma rpida y sencilla.
El mtodo consiste en asignar a la estructura de barras un objeto matemtico, llamado
matriz de rigidez, que relaciona los desplazamientos de un conjunto de puntos de la
estructura, llamados nodos, con las fuerzas exteriores que es necesario aplicar para lograr

esos desplazamientos (las componentes de esta matriz son fuerzas generalizadas asociadas
a desplazamientos generalizados). La matriz de rigidez relaciona las fuerzas nodales
equivalentes y desplazamientos sobre los nodos de la estructura, mediante la siguiente
ecuacin:

Donde: son las fuerzas nodales equivalentes asociadas a las fuerzas exteriores aplicadas
sobre la estructura;
son las reacciones hiperestticas inicialmente desconocidas sobre la
estructura; los desplazamientos nodales incgnita de la estructura y el nmero de
grados de libertad de la estructura.
La energa de deformacin elstica tambin puede expresarse en trminos de la matriz de
rigidez mediante la relacin:

Del teorema de Maxwell-Betti se deduce que la matriz de rigidez debe ser simtrica y por
tanto:

El mtodo matricial requiere asignar a cada barra elstica de la estructura una matriz de
rigidez, llamada matriz de rigidez elemental que depender de sus condiciones de enlace
extremo (articulacin, nudo rgido,...), la forma de la barra (recta, curvada, etc.) y las
constantes elsticas del material de la barra (mdulo de elasticidad longitudinal y mdulo
de elasticidad transversal). A partir del conjunto de matrices elementales mediante un
algoritmo conocido como acoplamiento que tiene en cuenta la conectividad de unas barras
con otras se obtiene una matriz de rigidez global, que relaciona los desplazamientos de los
nudos con las fuerzas equivalentes sobre los mismos.
Igualmente a partir de las fuerzas aplicadas sobre cada barra se construye el llamado vector
de fuerzas nodales equivalentes que dependen de las acciones exteriores sobre la estructura.
Junto con estas fuerzas anteriores deben considerarse las posibles reacciones sobre la
estructura en sus apoyos o enlaces exteriores (cuyos valores son incgnitas).
Finalmente se construye un sistema lineal de ecuaciones, para los desplazamientos y las
incgnitas. El nmero de reacciones incgnita y desplazamientos incgnita depende del
nmero de nodos: es igual a 3N para problemas bidimensionales, e igual a 6N para un

problema tridimensional. Este sistema siempre puede ser dividido en dos subsistemas de
ecuaciones desacoplados que cumplen:

Subsistema 1. Que agrupa todas las ecuaciones lineales del sistema original que slo
contienen desplazamientos incgnita.

Subsistema 2. Que agrupa al resto de ecuaciones, y que una vez resuelto el


subsistema 1 y substituido sus valores en el subsistema 2 permite encontrar los
valores de las reacciones incgnita.

Una vez resuelto el subsistema 1 que da los desplazamientos, se substituye el valor de estos
en el subsistema 2 que es trivial de resolver. Finalmente a partir de las reacciones, fuerzas
nodales equivalentes y desplazamientos se encuentran los esfuerzos en los nudos o uniones
de las barras a partir de los cuales pueden conocerse los esfuerzos en cualquier punto de la
estructura y por tanto sus tensiones mximas, que permiten dimensionar adecuadamente
todas las secciones de la estructura.
El uso de las matrices en la Ingeniera Civil es muy importante para resolver un diverso tipo
de problemas, principalmente en el rea de anlisis y diseo estructural.
El mtodo matricial de la rigidez, por ejemplo, es de gran utilidad para estudiar una
estructura, determinando su estabilidad por medio de tres tipos de ecuaciones que deben
cumplirse:

Ecuaciones de compatibilidad.

Ecuaciones constitutivas.

Ecuaciones de equilibrio.

De esta manera, vemos la importancia que tiene una materia bsica y fundamental en la
Ingeniera como es el lgebra, adems de sus aplicaciones prcticas en la vida profesional.
Ejemplo de aplicacin de matrices
El precio para los productos A, B, C y D por unidad son los siguientes: $3.80, $4.90,
$6.50, $10.80; y las cantidades que se adquieren de cada producto son: A = 500, B = 600,
C = 850, D = 720 Determina el costo total de las adquisiciones:
Solucin aplicando matrices:

P 3.80 4.90 6.50 10.80(1x 4)

500
600
(4 x1)
C
850

720

Se cumple la condicin del nmero de columnas es igual al nmero de renglones

En donde.
(3.80)(500) (4.90)(600) (6.50)(850) (10.80)( 720) 18141

PC 18141 (1x1)
Por lo tanto el Costo Total es de $18,141

Ejemplo para resolver un sistema de ecuaciones a travs de la matriz:


Sistema de ecuaciones lineales
A11x1 A12 x2 ...... A1n xn b1
A21x1 A22 x2 ...... A2 n xn b2
.
.
.
An1 x1 An 2 x2 ...... Ann xn bn
En forma matricial:
O sea

AX B

A = matriz de coeficientes numricos de las variables


X = matriz de las variables
B = matriz de resultados
Multiplicando ambos miembros de la igualdad por la matriz inversa

A1 * AX A1B

En donde
tenemos

A1 * A I matriz inversa
IX A1B X A1B

Para determinar en valor de las variables se determina primero la matriz inversa como se
indica a continuacin:

Matriz inversa
La inversa de un matriz se emplea en la resolucin de ecuaciones lineales simultneas y en
otros anlisis.
El producto de una matriz por su matriz inversa es igual a la matriz unidad
1
A1
A

Matriz inversa

nicamente las matrices cuadradas tienen inversa


Manera de obtener la inversa de una matriz aplicando el mtodo de Gauss-Jordan
3 7
(2 x 2)
2 5

A
Ejemplo: Dada la matriz

Determina la matriz inversa

Primera fase
1. El primer rengln se divide entre el trmino A11
3 7 1 0
1 7
3
3

1 0
3

2. El rengln base se multiplica por el trmino A21 con signo contrario

1 0 2 2 14
2
0
3
3
3

3. Despus el rengln que se va a modificar se suma algebraicamente al resultado


anterior

14

0 2 5 0 1 0

1 2 1
3
3

Segunda fase
1. El segundo rengln se divide entre el termino

1 2 1
3
3
0 1 2 3
1
3

A22

2. El rengln base se multiplica por el termino

1 2 3 7

14

A12

3. Despus el rengln que se va a modificar se suma algebraicamente al resultado


anterior

14

7 1 7

1 0 1 0 5 7
3

Comprobacin:

3 7 5
A. * A 1

2 5 2

7 1 0

3 0 1

En donde:

(3)(5) (7)( 2) 111


(3)( 7) (7)(3) 012
(2)(5) (5)( 2) 021
( 2)( 7) (5)(3) 122
Despus se multiplica la matriz inversa por la matriz de resultados y se obtiene el valor de
las variables.
El ejemplo anterior es de una matriz de (2 x 2) pero el procedimiento es el mismo para la
matriz de (3 x 3), el rengln base es la herramienta para modificar uno o ms renglones.

También podría gustarte