CAPÍTULO 12

Estudio de casos: ecuaciones

algebraicas lineales
El propósito de este capítulo es usar los procedimientos numéricos analizados en los
capítulos 9, 10 y 11 para resolver sistemas de ecuaciones algebraicas lineales, en algunas
aplicaciones a la ingeniería. Dichas técnicas numéricas sistemáticas tienen significado
práctico, ya que los ingenieros con mucha frecuencia se enfrentan a problemas que in-
volucran sistemas de ecuaciones que son demasiado grandes para resolverse a mano. Los
algoritmos numéricos en estas aplicaciones son particularmente adecuados para imple-
mentarse en computadoras personales.
En la sección 12.1 se muestra cómo se emplea un balance de masa para modelar un
sistema de reactores. En la sección 12.2 se le da especial énfasis al uso de la matriz in-
versa para determinar las complicadas interacciones causa-efecto entre las fuerzas en los
elementos de una armadura. La sección 12.3 constituye un ejemplo del uso de las leyes
de Kirchhoff para calcular las corrientes y voltajes en un circuito con resistores. Por úl-
timo, la sección 12.4 es una ilustración de cómo se emplean las ecuaciones lineales para
determinar la configuración en estado estacionario de un sistema masa-resorte.

12.1 ANÁLISIS EN ESTADO ESTACIONARIO DE UN SISTEMA

DE REACTORES (INGENIERÍA QUÍMICA/BIOINGENIERÍA)

Antecedentes. Uno de los principios de organización más importantes en la ingenie-

ría química es la conservación de la masa (recuerde la tabla 1.1). En términos cuantita-
tivos, el principio se expresa como un balance de masa que toma en cuenta todas las
fuentes y sumideros de un fluido que entra y sale de un volumen (figura 12.1). En un
periodo finito, esto se expresa como
Acumulación = entradas – salidas (12.1)

El balance de masa representa un ejercicio de contabilidad para la sustancia en

particular que se modela. Para el periodo en que se calcula, si las entradas son mayores
que las salidas, la masa de la sustancia dentro del volumen aumenta. Si las salidas son
mayores que las entradas, la masa disminuye. Si las entradas son iguales a las salidas,
la acumulación es cero y la masa permanece constante. Para esta condición estable, o
en estado estacionario, la ecuación (12.1) se expresa como
Entradas = salidas (12.2)

Emplee la conservación de la masa para determinar las concentraciones en estado esta-

cionario de un sistema de reactores conectados.
328 ESTUDIO DE CASOS: ECUACIONES ALGEBRAICAS LINEALES

Volumen

Salida
Entrada

Acumulación

FIGURA 12.1
Una representación esquemática del balance de masa.

Solución. Se puede usar el balance de masa para resolver problemas de ingeniería al

expresar las entradas y salidas en términos de variables y parámetros medibles. Por ejem-
plo, si se realiza un balance de masa para una sustancia conservativa (es decir, aquella
que no aumente ni disminuya debido a transformaciones químicas) en un reactor (figura
12.2), podríamos cuantificar la velocidad con la cual el flujo de la masa entra al reactor
a través de dos tuberías de entrada y sale de éste a través de una tubería de salida. Esto se
hace mediante el producto de la velocidad del fluido o caudal Q (en metros cúbicos por
minuto) por la concentración c (en miligramos por metro cúbico) en cada tubería.
Por ejemplo, en la tubería 1 de la figura 12.2, Q1 = 2 m3/min y c1 = 25 mg/m3; por lo
tanto, la velocidad con la cual la masa fluye hacia el reactor a través de la tubería 1 es
Q1c1 = (2 m3/min)(25 mg/m3) = 50 mg/min. Así, 50 mg de sustancias químicas fluyen
cada minuto hacia el interior del reactor a través de esta tubería. De forma similar, para
la tubería 2 la velocidad de masa que entra se calcula como Q 2 c 2 = (1.5 m3 /min)
(10 mg/m3) = 15 mg/min.
Observe que la concentración a la salida del reactor a través de la tubería 3 no se
especifica en la figura 12.2. Esto es así porque ya se tiene información suficiente para
calcularla con base en la conservación de la masa. Como el reactor se halla en estado
estacionario se aplica la ecuación (12.2) y las entradas deberán estar en balance con las
salidas,
Q1c1 + Q2c2 = Q3c3

Sustituyendo los valores dados en esta ecuación se obtiene

50 + 15 = 3.5c3
de la cual se despeja c3 = 18.6 mg/m3. De esta forma, hemos determinado la concen-
tración en la tercera tubería. Sin embargo, del cálculo se obtiene algo más. Como el
reactor está bien mezclado (representado por el agitador en la figura 12.2), la concentra-
ción será uniforme, u homogénea, en todo el tanque. Por lo que, la concentración en la
tubería 3 deberá ser idéntica a la concentración en todo el reactor. En consecuencia, el
balance de masa nos ha permitido calcular tanto la concentración en el reactor como en
el tubo de salida. Esta información es de gran utilidad para los ingenieros químicos y
 12.1 ANÁLISIS EN ESTADO ESTACIONARIO DE UN SISTEMA DE REACTORES 329

Q1 = 2 m3/min
FIGURA 12.2 c1 = 25 mg/m3
Un reactor en estado
Q3 = 3.5 m3/min
estacionario, completamente c3 = ?
mezclado, con dos tuberías
de entrada y una de salida. Q2 = 1.5 m3/min
Los caudales Q están en c2 = 10 mg/m3
metros cúbicos por minuto, y
las concentraciones c están
en miligramos por metro
cúbico.

Q15 = 3 Q55 = 2
c5 Q54 = 2

Q25 = 1

Q01 = 5 Q12 = 3 Q24 = 1 Q44 = 11

c1 c2 c4
c01 = 10

Q23 = 1 Q34 = 8

Q31 = 1
Q03 = 8 c3
FIGURA 12.3
c03 = 20
Cinco reactores conectados
por tuberías.

petroleros, quienes tienen que diseñar reactores que tengan mezclas de una concentración
específica.
Debido a que se utilizó álgebra simple para determinar la concentración de un solo
reactor en la figura 12.2, podría no ser obvio lo que tiene que hacer una computadora en
el cálculo de un balance de masa. En la figura 12.3 se muestra un problema donde las
computadoras no solamente son útiles, sino que son de una enorme necesidad práctica.
Debido a que hay cinco reactores interconectados o acoplados, se necesitan cinco ecua-
ciones de balance de masa para caracterizar el sistema. En el reactor 1, velocidad de la
masa que entra es
5(10) + Q31c3
y la velocidad de la masa que sale es
Q12c1 + Q15c1
330 ESTUDIO DE CASOS: ECUACIONES ALGEBRAICAS LINEALES

Como el sistema se encuentra en estado estacionario, los flujos de entrada y de salida

deben ser iguales:
5(10) + Q31c3 = Q12c1 + Q15c1
o, sustituyendo los valores de la figura 12.3,
6c1 – c3 = 50
Ecuaciones similares se obtienen para los otros reactores:

–3c1 + 3c2 = 0
–c2 + 9c3 = 160
–c2 – 8c3 + 11c4 – 2c5 = 0
–3c1 – c2 + 4c5 = 0

Se puede utilizar un método numérico para resolver estas cinco ecuaciones con las
cinco incógnitas que son las concentraciones:
{C}T = 11.51 11.51 19.06 17.00 11.51
Además, la matriz inversa se calcula como

0.16981 0.00629 0.01887 0 0

0.16981 0.33962 0.01887 0 0
[A] –1 = 0.01887 0.03774 0.11321 0 0
0.06003 0.07461 0.08748 0.09091 0.04545
0.16981 0.08962 0.01887 0 0.25000

Cada uno de los elementos aij significa el cambio en la concentración del reactor i debi-
do a un cambio unitario en la carga del reactor j. De esta forma, los ceros en la columna
4 indican que una carga en el reactor 4 no influirá sobre los reactores 1, 2, 3 y 5. Esto es
consistente con la configuración del sistema (figura 12.3), la cual indica que el flujo de
salida del reactor 4 no alimenta ningún otro reactor. En cambio, las cargas en cualquie-
ra de los tres primeros reactores afectarán al sistema completo, como se indica por la
ausencia de ceros en las primeras tres columnas. Tal información es de gran utilidad
para los ingenieros que diseñan y manejan sistemas como éste.

12.2 ANÁLISIS DE UNA ARMADURA ESTÁTICAMENTE

DETERMINADA (INGENIERÍA CIVIL/AMBIENTAL)

Antecedentes. Un problema importante en la ingeniería estructural es encontrar las

fuerzas y reacciones asociadas con una armadura estáticamente determinada. En la fi-
gura 12.4 se muestra el ejemplo de una armadura.
Las fuerzas (F) representan ya sea la tensión o la compresión sobre los componen-
tes de la armadura. Las reacciones externas (H2 , V2 y V3) son fuerzas que caracterizan
cómo interactúa dicha estructura con la superficie de soporte. El apoyo fijo en el nodo
2 puede transmitir fuerzas horizontales y verticales a la superficie, mientras que el apo-
yo móvil en el nodo 3 transmite sólo fuerzas verticales. Se observa que el efecto de la
carga externa de 1 000 lb se distribuye entre los componentes de la armadura.
 12.2 ANÁLISIS DE UNA ARMADURA ESTÁTICAMENTE DETERMINADA 331

1 000 lb

F1 908 F3

608
H2 2 308
3
F2

FIGURA 12.4
Fuerzas en una armadura V2 V3
estáticamente determinada.

F1,v

1 F1,h
308
608

F1 F3
F1 F3
F2,v F3,v
608
H2 2 308 3 F3,h
FIGURA 12.5 F2
Diagramas de fuerza F2
de cuerpo libre para los F2,h
nodos de una armadura V2 V3
estáticamente determinada.

Solución. Este tipo de estructura se puede describir como un conjunto de ecuaciones

algebraicas lineales acopladas. Los diagramas de fuerza de cuerpo libre para cada nodo
se muestran en la figura 12.5. La suma de las fuerzas en ambas direcciones, vertical y
horizontal, deben ser cero en cada nodo, ya que el sistema está en reposo. Por lo tanto,
para el nodo 1,
Σ FH = 0 = –F1 cos 30° + F3 cos 60° + F1,h (12.3)

Σ FV = 0 = –F1 sen 30° – F3 sen 60° + F1,v (12.4)

para el nodo 2,
Σ FH = 0 = F2 + F1 cos 30° + F2,h + H2 (12.5)

Σ FV = 0 = F1 sen 30° + F2,v + V2 (12.6)

332 ESTUDIO DE CASOS: ECUACIONES ALGEBRAICAS LINEALES

para el nodo 3,
Σ FH = 0 = –F2 – F3 cos 60° + F3,h (12.7)
Σ FV = 0 = F3 sen 60° + F3,v + V3 (12.8)

donde Fi,h es la fuerza horizontal externa aplicada sobre el nodo i (se considera que una
fuerza positiva va de izquierda a derecha) y Fi,v es la fuerza vertical externa que se apli-
ca sobre el nodo i (donde una fuerza positiva va hacia arriba). Así, en este problema, la
fuerza de 1 000 lb hacia abajo en el nodo 1 corresponde a F1,v = –1 000 libras. En este
caso, todas las otras Fi,v y Fi,h son cero. Observe que las direcciones de las fuerzas in-
ternas y de las reacciones son desconocidas. La aplicación correcta de las leyes de Newton
requiere sólo de suposiciones consistentes respecto a la dirección. Las soluciones son
negativas si las direcciones se asumen de manera incorrecta. También observe que en
este problema, las fuerzas en todos los componentes se suponen en tensión y actúan ti-
rando de los nodos adyacentes. Una solución negativa, por lo tanto, corresponde a com-
presión. Este problema se plantea como el siguiente sistema de seis ecuaciones con seis
incógnitas:
0.866 0 –0.5 0 0 0 F1 0
0.5 0 0.866 0 0 0 F2 –1 000
–0.866 –1 0 –1 0 0 F3 = 0
(12.9)
–0.5 0 0 0 –1 0 H2 0
0 1 0.5 0 0 0 V2 0
0 0 –0.866 0 0 –1 V3 0
Observe que, como se formuló en la ecuación (12.9), se requiere de pivoteo parcial
para evitar la división entre cero de los elementos de la diagonal. Con el uso de una es-
trategia de pivote, el sistema se resuelve mediante cualquiera de las técnicas de elimi-
nación que se analizaron en los capítulos 9 y 10. Sin embargo, como este problema es
un caso de estudio ideal, para demostrar la utilidad de la matriz inversa se utiliza la
descomposición LU para calcular
F1 = –500 F2 = 433 F3 = –866
H2 = 0 V2 = 250 V3 = 750

la matriz inversa es
0.866 0.5 0 0 0
0.25 –0.433 0 0 1 0
–1
[A] = –0.5 0.866 0 0 0 0
–1 0 –1 0 –1 0
–0.433 –0.25 0 –1 0 0
0.433 –0.75 0 0 0 –1
Ahora, observe que el vector del lado derecho representa las fuerzas horizontales y
verticales aplicadas externamente sobre cada nodo,
{F}T = F1,h F1,v F2,h F2,v F3,h F3,v (12.10)

Debido a que las fuerzas externas no tienen efecto sobre la descomposición LU, no
se necesita aplicar el método una y otra vez para estudiar el efecto de diferentes fuerzas
 12.2 ANÁLISIS DE UNA ARMADURA ESTÁTICAMENTE DETERMINADA 333

1 000 1 000

6 6

50

50
86 86

0
2 000 1 000 250 2 000 1 250 1 000

433 433 433 433

a) b)

FIGURA 12.6
Dos casos de prueba que muestran a) vientos desde la izquierda y b) vientos desde la derecha.

externas sobre la armadura. Todo lo que hay que hacer es ejecutar los pasos de sustitución
hacia adelante y hacia atrás, para cada vector del lado derecho, y así obtener de manera
eficiente soluciones alternativas. Por ejemplo, podríamos querer estudiar el efecto de
fuerzas horizontales inducidas por un viento que sopla de izquierda a derecha. Si la
fuerza del viento se puede idealizar como dos fuerzas puntuales de 1 000 libras sobre
los nodos 1 y 2 (figura 12.6a), el vector del lado derecho es
{F}T = –1 000 0 1 000 0 0 0

que se utiliza para calcular

F1 = –866 F2 = 250 F3 = –500
H2 = –2 000 V2 = –433 V3 = 433
Para un viento de la derecha (figura 12.6b), F1,h = –1 000, F3,h = –1 000, y todas las demás
fuerzas externas son cero, con lo cual resulta
F1 = –866 F2 = –1 250 F3 = 500
H2 = 2 000 V2 = 433 V3 = –433
Los resultados indican que los vientos tienen efectos marcadamente diferentes sobre la
estructura. Ambos casos se presentan en la figura 12.6.
Cada uno de los elementos de la matriz inversa tienen también utilidad directa para
aclarar las interacciones estímulo-respuesta en la estructura. Cada elemento representa
el cambio de una de las variables desconocidas a un cambio unitario de uno de los es-
–1
tímulos externos. Por ejemplo, el elemento a32 indica que la tercera incógnita (F3) cam-
biará a 0.866 debido a un cambio unitario del segundo estímulo externo (F1,v). De esta
forma, si la carga vertical en el primer nodo fuera aumentada en 1, F3 se podría aumen-
tar en 0.866. El hecho de que los elementos sean 0 indica que ciertas incógnitas no se
ven afectadas por algunos de los estímulos externos. Por ejemplo, a –1 13 = 0 significa que
F1 no se ve afectado por cambios en F2,h. Esta habilidad de aislar interacciones tiene
334 ESTUDIO DE CASOS: ECUACIONES ALGEBRAICAS LINEALES

diversas aplicaciones en la ingeniería; éstas comprenden la identificación de aquellos

componentes que son más sensibles a estímulos externos y, como una consecuencia, más
propensos a fallar. Además, esto sirve para determinar los componentes que son inne-
cesarios (véase el problema 12.18).
El procedimiento anterior resulta particularmente útil cuando se aplica a grandes
estructuras complejas. En la práctica de la ingeniería, en ocasiones es necesario resolver
estructuras con cientos y aun miles de elementos estructurales. Las ecuaciones lineales
proporcionan un medio poderoso para ganar cierta comprensión del comportamiento de
dichas estructuras.

12.3 CORRIENTES Y VOLTAJES EN CIRCUITOS

CON RESISTORES (INGENIERÍA ELÉCTRICA)

Antecedentes. Un problema común en ingeniería eléctrica es la determinación de

corrientes y voltajes en algunos puntos de los circuitos con resistores. Tales problemas
se resuelven utilizando las reglas para corrientes y voltajes de Kirchhoff. La regla para
las corrientes (o nodos) establece que la suma algebraica de todas las corrientes que
entran a un nodo debe ser cero (véase figura 12.7a), o
Σi = 0 (12.11)

donde todas las corrientes que entran al nodo se consideran de signo positivo. La regla
de las corrientes es una aplicación del principio de la conservación de la carga (recuerde
la tabla 1.1).
La regla para los voltajes (o mallas) especifica que la suma algebraica de las dife-
rencias de potencial (es decir, cambios de voltaje) en cualquier malla debe ser igual a
cero. Para un circuito con resistores, esto se expresa como
Σx – ΣiR = 0 (12.12)

donde x es la fem (fuerza electromotriz) de las fuentes de voltaje, y R es la resistencia

FIGURA 12.7 de cualquier resistor en la malla. Observe que el segundo término se obtiene de la ley de
Representaciones esquemá- Ohm (figura 12.7b), la cual establece que la caída de voltaje a través de un resistor ideal
ticas de a) la regla de las co- es igual al producto de la corriente por la resistencia. La regla de Kirchhoff para el
rrientes de Kirchhoff y b) la ley voltaje es una expresión de la conservación de la energía.
de Ohm.

i1 i3
FIGURA 12.8
Un circuito con resistores para resolverse usando ecuaciones algebraicas lineales
simultáneas.
i2 R = 10 V R=5V
3 2 1
a) V1 = 200 V

Rij
Vi Vj R=5V R = 10 V

iij V6 = 0 V
b) 4 5 6
R = 15 V R = 20 V
12.3 CORRIENTES Y VOLTAJES EN CIRCUITOS CON RESISTORES 335

3 2 1

i32 i12

i43 i52
i54 i65

4 5 6

FIGURA 12.9
Corrientes supuestas.

Solución. La aplicación de estas reglas da como resultado un sistema de ecuaciones

algebraicas lineales simultáneas, ya que las mallas que forman un circuito están conec-
tadas. Por ejemplo, considere el circuito de la figura 12.8. Las corrientes asociadas con
este circuito son desconocidas, tanto en magnitud como en dirección. Esto no presenta
gran dificultad, ya que tan sólo se supone una dirección para cada corriente. Si la solu-
ción resultante a partir de las leyes de Kirchhoff es negativa, entonces la dirección su-
puesta fue incorrecta. Por ejemplo, la figura 12.9 muestra direcciones supuestas para las
corrientes.
Dadas estas suposiciones, la regla de la corriente de Kirchhoff se aplica a cada nodo
para obtener

i12 + i52 + i32 = 0

i65 – i52 – i54 = 0
i43 – i32 = 0
i54 – i43 = 0

La aplicación de la regla de voltajes en cada una de las mallas da

–i54R54 – i43R43 – i32 R32 + i52 R52 = 0
–i65R65 – i52 R52 + i12 R12 – 200 = 0
o, sustituyendo el valor de las resistencias de la figura 12.8 y pasando las constantes al
lado derecho,

–l5i54 – 5i43 – l0i32 + l0i52 = 0

–20i65 – 10i52 + 5i12 = 200

Por lo tanto, el problema consiste en la solución del siguiente conjunto de seis ecuaciones
con seis corrientes como incógnitas:
1 1 1 0 0 0 i12 0
0 –1 0 1 –1 0 i52 0
0 0 –1 0 0 1 i32 = 0
0 0 0 0 1 –1 i65 0
0 10 –10 0 –15 –5 i54 0
5 –10 0 –20 0 0 i43 200
336 ESTUDIO DE CASOS: ECUACIONES ALGEBRAICAS LINEALES

V = 153.85 V = 169.23
V = 200

FIGURA 12.10
La solución obtenida para i = 1.5385 i = 6.1538
las corrientes y voltajes
usando un método de V=0
eliminación. V = 146.15 V = 123.08

Aunque no es práctico resolverlo a mano, este sistema se resuelve de manera sencilla

con un método de eliminación. Si se procede de esta forma, la solución es
i12 = 6.1538 i52 = –4.6154 i32 = –1.5385
i65 = –6.1538 i54 = –1.5385 i43 = –1.5385

Así, con una interpretación adecuada de los signos del resultado, las corrientes y volta-
jes en el circuito se muestran en la figura 12.10. Deben ser evidentes las ventajas de usar
algoritmos numéricos y computadoras para problemas de este tipo.

12.4 SISTEMAS MASA-RESORTE

(INGENIERÍA MECÁNICA/AERONÁUTICA)

Antecedentes. Los sistemas idealizados masa-resorte desempeñan un papel impor-

tante en la mecánica y en otros problemas de ingeniería. En la figura 12.11 se presenta
un sistema de este tipo. Después de liberar las masas, éstas son jaladas hacia abajo por
la fuerza de gravedad. Observe que el desplazamiento resultante en cada resorte de la
figura 12.11b se mide a lo largo de las coordenadas locales referidas a su posición inicial
en la figura 12.11a.
Como se mencionó en el capítulo 1, la segunda ley de Newton se emplea en conjun-
to con el equilibrio de fuerzas para desarrollar un modelo matemático del sistema. Para
cada masa, la segunda ley se expresa como
d2x
m = FD – FU (12.13)
dt 2
Para simplificar el análisis se supondrá que todos los resortes son idénticos y que se
comportan de acuerdo con la ley de Hooke. En la figura 12.12a se muestra un diagrama
de cuerpo libre para la primera masa. La fuerza hacia arriba es únicamente una expresión
directa de la ley de Hooke:
FU = kx1 (12.14)

Las componentes hacia abajo consisten en las dos fuerzas del resorte junto con la acción
de la gravedad sobre la masa,
FD = k(x2 – x1) + k(x2 – x1) = m1g (12.15)

Observe cómo la componente de fuerza de los dos resortes es proporcional al desplaza-

miento de la segunda masa, x2, corregida por el desplazamiento de la primera masa, x1.
 12.4 SISTEMAS MASA-RESORTE 337

FIGURA 12.11
k
Un sistema compuesto de
tres masas suspendidas m1 0
verticalmente por una serie m1 x1
de resortes. a) El sistema k k
antes de ser liberado, m2 0
es decir, antes de la
extensión o compresión k
m2 x2
de los resortes. b) El
m3 0
sistema después de ser
liberado. Observe que las
posiciones de las masas m3 x3
están en referencia a las
coordenadas locales con a) b)
orígenes en su posición
antes de ser liberadas.

kx1 k(x2 – x1) k(x2 – x1) k(x3 – x2)

m1 m2 m3

k(x2 – x1) m1g k(x2 – x1) m2g k(x3 – x2) m3g

a) b) c)

FIGURA 12.12
Diagramas de cuerpo libre para las tres masas de la figura 12.11.

Las ecuaciones (12.14) y (12.15) se sustituyen en la ecuación (12.13) para dar

d 2 x1
m1 = 2 k ( x 2 – x1 ) + m1g – kx1 (12.16)
dt 2

De esta forma, se ha obtenido una ecuación diferencial ordinaria de segundo orden para
describir el desplazamiento de la primera masa con respecto al tiempo. Sin embargo,
advierta que la solución no se puede obtener, ya que el modelo tiene una segunda varia-
ble dependiente, x2. En consecuencia, se deben desarrollar diagramas de cuerpo libre
para la segunda y tercera masa (figuras 12.12b y c) que se emplean para obtener
d 2 x2
m2 = k ( x3 – x 2 ) + m2 g – 2 k ( x 2 – x1 ) (12.17)
dt 2
338 ESTUDIO DE CASOS: ECUACIONES ALGEBRAICAS LINEALES

y
d 2 x3
m3 = m3 g – k ( x3 – x 2 ) (12.18)
dt 2

Las ecuaciones (12.16), (12.17) y (12.18) forman un sistema de tres ecuaciones dife-
renciales con tres incógnitas. Con las condiciones iniciales apropiadas, estas ecuaciones
sirven para calcular los desplazamientos de las masas como una función del tiempo (es
decir, sus oscilaciones). En la parte siete estudiaremos los métodos numéricos para ob-
tener tales soluciones. Por ahora, podemos obtener los desplazamientos que ocurren cuan-
do el sistema eventualmente llega al reposo, es decir, al estado estacionario. Para esto se
igualan a cero las derivadas en las ecuaciones (12.16), (12.17) y (12.18), obteniéndose
3kx1 – 2kx2 = m1g
–2kx1 + 3kx2 – kx3 = m 2g
– kx2 + kx3 = m3g
o, en forma matricial,
[K]{X} = {W}
donde [K], conocida como matriz de rigidez, es
3k –2k
[K] = –2k 3k –k
–k k

y {X} y {W} son los vectores columna de las incógnitas X y de los pesos mg, respecti-
vamente.
Solución. Aquí se emplean métodos numéricos para obtener una solución. Si m1 = 2
kg, m2 = 3 kg, m3 = 2.5 kg, y todas las k = 10 kg/s2, use la descomposición LU con el
propósito de obtener los desplazamientos y generar la inversa de [K].
Sustituyendo los parámetros del modelo se obtiene

30 –20 19.6
[K] = –20 30 –10 {W} = 29.4
–10 10 24.5
La descomposición LU se utiliza con el objetivo de obtener x1 = 7.35, x2 = 10.045 y x3 =
12.495. Estos desplazamientos se utilizaron para construir la figura 12.11b. La inversa
de la matriz de rigidez calculada es
0.1 0.1 0.1
[K] –1 = 0.1 0.15 0.15
0.1 0.15 0.25

Cada elemento de la matriz k–1 ji nos indica el desplazamiento de la masa i debido a

una fuerza unitaria impuesta sobre la masa j. Así, los valores 0.1 en la columna 1 nos
indican que una carga unitaria hacia abajo en la primera masa desplazará todas las ma-
sas 0.1 m hacia abajo. Los otros elementos se interpretan en forma similar. Por lo tanto,
la inversa de la matriz de rigidez proporciona una síntesis de cómo los componentes del
sistema responden a fuerzas que se aplican en forma externa.