Las Ecuaciones Diferenciales y Sus Aplic

LAS ECUACIONES DIFERENCIALES
Y SUS APLICACIONES EN LA
INGENIERÍA
INDICE:
-Generalidades. Pg (1-4)
-Etapas de resolución del problema científico. Pg (5)
.
Formulación matemática del problema científico.
.
Solución de las ecuaciones.
.
Interpretación científica de la solución.
-Aplicaciones de ecuaciones diferenciales de primer orden y simples de

orden superior. Pg (7-30)
1. Aplicaciones a la mecánica:
1.1 Introducción.
1.2 Las leyes del movimiento de Newton.
2. Aplicaciones a los circuitos eléctricos:
2.1 Introducción.
2.2 La ley de Kirchhoff.
3. Aplicaciones a flujo de calor en estado estacionario.
4. Aplicaciones a problemas combinados de crecimiento y decrecimiento.
5. El cable colgante.
6. La deflexión de vigas.
-Aplicaciones de ecuaciones diferenciales lineales. Pg (31-50)
1. Movimiento vibratorio de sistemas mecánicos:

1.1 El resorte vibrante (movimiento armónico simple).
1.2 El resorte vibrante con amortiguamiento (movimiento
amortiguado).
1.3 El resorte con fuerzas externas.
1.4 La resonancia mecánica.
Las Ecuaciones Diferenciales y sus Aplicaciones en la Ingeniería
LAS ECUACIONES DIFERENCIALES

Y SU APLICACIÓN A LA INGENIERÍA
GENERALIDADES:
El descubrimiento de Newton y Leibniz en el siglo diecisiete sobre las ideas básicas

del cálculo integral fue crucial para el avance que sufrieron las matemáticas, y más
importante fue, si cabe, la relación que encontraron entre el cálculo integral y el
diferencial, ya que consiguieron fundirlos en uno solo. Una de las aplicaciones de este
descubrimiento fue la física aplicada, dícese, la Ingeniería.
El maestro de Newton, Isaac Barrow, conocía ya la existencia de la relación entre la
tangente en un punto a una curva (derivada) y el área de una región limitada de una curva
(Integral Definida), pero fueron Newton y Leibniz los que comprendieron la importancia
de esa relación.
La derivada se utilizó, en principio, para el cálculo de la tangente en un punto, y pronto se
vió que también servía para el cálculo de velocidades, y en consecuencia para el estudio
de la variación de una función.
Desde los primeros pasos en el cálculo diferencial, de todos es conocido que dada una
= f ´(x) , en forma de diferencial de una función de una

dy
función y = f(x), su derivada
dx
sola variable, es también una función que se puede encontrar mediante ciertas reglas como
el Teorema Fundamental del Cálculo Integral, que nos muestra la vinculación entre la
derivada de una función y la integral de dicha función ; si F(x) es la función integral que
debe ser integrable en el intervalo [a,x] para cada x de [a,b], siendo c tal que a ≤ c ≤ b
tenemos que F ( x) = ∫ f (t )dt si a ≤ x ≤ b , existe entonces F´(x) en cada punto x del

x
intervalo abierto (a,b), en el que f es continua, y para tal x tenemos F´(x) = f ( x)

quedando demostrado la relación entre Integral y Derivada.
1
-La Derivada de la Integral de una función es la propia función:

F´(x) = f ( x)
-La Integral de la Derivada de una función es la propia función:
f ( x) = ∫ f ´(x)dx
x
Con lo antes mencionado, a lo que se une La Regla de Barrow (que no es más que la
aplicación del teorema fundamental), es posible conseguir la función primitiva de la
= f ´(x) mediante la integración de dicha función, que es lo que

dy
función derivada
dx
necesitamos para poder resolver las ecuaciones diferenciales, pero antes debemos
definirlas.
Hay una gran variedad de problemas en los cuales se desea conocer un elemento variable a
partir de su coeficiente de variación, o dicho de otra forma, queremos conocer cómo varía
dicho elemento en función de una o varias variables.
En definitiva, lo que se pretende es determinar una función desconocida mediante datos
relacionados por una ecuación que contiene, por lo menos, una de las derivadas de la
función desconocida.
Estas ecuaciones se denominan ecuaciones diferenciales y su estudio por parte de Newton,
Leibniz y los Bernouilli para resolver algunas de las ecuaciones diferenciales sencillas que
se presentaron en geometría y mecánica, llevaron al conocimiento sobre la resolución de
ciertos tipos de ecuaciones diferenciales; se conoce mediante la práctica que es difícil
obtener teorías matemáticas de gran generalidad para la resolución de estas ecuaciones
diferenciales, salvo para algunos tipos, como las ecuaciones lineales, muy extendidas para
problemas de tipo científico.
Definimos:
-Ecuación diferencial (E.D.) a una ecuación que relaciona una función (o variable
dependiente), su variable o variables (variables independientes), y sus derivadas.
Si la ecuación contiene derivadas respecto a una sola variable independiente, entonces se
dice que es una ecuación diferencial ordinaria (E.D.O); y si contiene las derivadas
parciales respecto a dos o más variables independientes, se llama ecuación en derivadas
2
parcialales (E.D.P.).
Otro tipo de ecuaciones son las ecuaciones diferenciales de retraso (o retardo) que están
caracterizadas por la presencia de un desplazamiento en el argumento o variable (x-x0).
-Se llama orden de la ecuación diferencial al orden de la derivada o derivada parcial

más alta que aparece en la ecuación.
Se dice que una ecuación diferencial (de orden n) está expresada en forma implícita
cuando tiene la forma F ( x, y, y´,...., y ( n ) ) = 0 , siendo F una función F : Ω ⊂ R n+2 → R
siendo Ω un subconjunto (generalmente abierto) de Rn+2
Se dice que una ecuación diferencial (de orden n) está expresada en forma explícita
cuando tenemos y(n)= f(x,y,y´,….,y(n-1)) con f : D ⊂ R n+1 → R siendo la función definida
en el subconjunto D (generalmente abierto) de Rn+1 .
-Se dice que una ecuación diferencial es lineal si tiene la forma

d n−1 y
+ + ... + a1 ( x) + a0 ( x) y = g ( x) y se llama lineal homogénea si
dny dy
an ( x ) n
a n−1 ( x ) n −1
dx dx dx
además g(x) = 0.
-Se dice que una función y = φ(x) definida en un intervalo I es solución de una
diferencial en el intervalo si, sustituida en dicha ecuación, la reduce a una identidad.
Una E. D. se dice resoluble (o integrable) por cuadraturas si su solución es expresable
mediante integrales.
En general, la solución de la ecuación diferencial de orden n dependerá de n parámetros.
Pero incluso de esta forma pueden no obtenerse todas las soluciones de una E. D. Por
ejemplo, cuando tenemos una familia uniparamétrica de soluciones de una E. D., una
sencilla interpretación geométrica nos muestra que también la envolvente de la familia de
curvas (si existe) es solución de la E. D.
-Se define como problema de valor inicial y problemas de valor frontera a

aquellos en que la ecuación diferencial se resuelve sujeta a unas condiciones dadas que la
función desconocida debe satisfacer.
3
Problema de valor inicial:

Es un problema que busca determinar una solución a una ecuación diferencial sujeta
a condiciones sobre la función desconocida y sus derivadas especificadas en un
valor de la variable independiente. Tales condiciones se llaman condiciones
iniciales.
Problemas de valor frontera:

Es un problema que busca determinar una solución a una ecuación diferencial sujeta
a condiciones sobre la función desconocida, especificadas en dos o más valores de
la variable independiente. Tales condiciones se llaman condiciones de frontera.
-La función primitiva resultante, o función solución de una ecuación diferencial,

puede tener por las condiciones iniciales o de frontera diversos valores, diferenciándose una
solución de otra en el parámetro, definiéndose este conjunto de soluciones familia de
soluciones de un parámetro (en el caso de existir sólo un parámetro) o familia de
soluciones de dos o más parámetros (en el caso de existir más de un parámetro)
4
ETAPAS DE RESOLUCIÓN DEL PROBLEMA CIENTÍFICO:
1) FORMULACIÓN MATEMÁTICA DEL PROBLEMA CIENTÍFICO:
Las leyes científicas, que por supuesto están basadas en experimentos u

observaciones, se traducen en ecuaciones matemáticas. En cada caso las ecuaciones
diferenciales representan una simplificación idealizada del problema físico con el que nos
encontramos, llamándose esta idealización Modelo Matemático.
Cada modelo es una aproximación a la realidad del problema físico, su aproximación y uso
del modelo sólo depende de los criterios impuestos a cada problema para su resolución.
Si la intuición o la evidencia del experimento coinciden con los resultados obtenidos por
medio del modelo podremos determinar cuan útil es ese modelo.
2) SOLUCIÓN DE LAS ECUACIONES:
Las ecuaciones formuladas en la etapa anterior necesitan ser resueltas, sujetas a

condiciones obtenidas del problema para determinar la incógnita o incógnitas involucradas.
Los procedimientos usados pueden producir una solución exacta o, en casos donde
soluciones exactas no se pueden obtener, soluciones aproximadas. Frecuentemente para
elaborar los cálculos numéricos se recurre al uso de la informática. El proceso de obtener
soluciones a menudo conduce a preguntas de naturaleza puramente matemática que
propician y propiciaron el avance de las susodichas matemáticas.
3) INTERPRETACIÓN CIENTÍFICA DE LA SOLUCIÓN:
Con el uso de las soluciones conocidas, el matemático o físico puede ser capaz de
interpretar lo que está sucediendo desde el punto de vista aplicado. Puede hacer
interpretaciones gráficas y tablas para poder comparar la teoría con lo obtenido de los
experimentos. Puede, incluso, basar una investigación posterior en las interpretaciones de
experimentos previos. Por supuesto que, si encuentra que los experimentos u observaciones
no están de acuerdo con la teoría, debe revisar el modelo matemático y su formulación
matemática hasta que se consiga un resultado cuyo margen de error lo marque la persona o
personas encargadas de los experimentos. Cada una de estas etapas es importante en la
solución final de un problema aplicado.
5
6
APLICACIONES DE ECUACIONES DIFERENCIALES DE PRIMER ORDEN Y

SIMPLES DE ORDEN SUPERIOR:
Tipos de aplicaciones:
1.1 Introducción.
2.1 Introducción.
1.1 Introducción:
La física trata de la investigación de las leyes que gobiernan el comportamiento del
universo físico. Por universo físico entendemos la totalidad de objetos a nuestro alrededor,
no sólo las cosas que observamos sino tambien las que no observamos, tales como los
átomos y moléculas. El estudio del movimiento de los objetos en nuestro universo es una
rama de la mecánica llamada dinámica formulada mediante las leyes del movimiento de
Newton. Para los objetos que se mueven muy rápido, cerca de la velocidad de la luz, no
podemos usar las leyes de Newton. En su lugar debemos usar una versión revisada de estas
leyes, desarrolladas por Einstein y conocidas como mecánica relativista, o mecánica de la
relatividad. Para objetos de dimensiones atómicas las leyes de Newton tampoco son
válidas. De hecho, para obtener descripciones precisas del movimiento de objetos de
dimensiones atómicas, necesitamos establecer un conjunto de leyes denominadas mecánica
cuántica. La mecánica cuántica y la relativista son muy complicadas, no siendo objeto de
estudio en este trabajo.
7
Las tres leyes del movimiento primero desarrolladas por Newton son:
1. Un cuerpo en reposo tiende a permanecer en reposo, mientras que un cuerpo en

movimiento tiende a persistir en movimiento en una línea recta con velocidad constante a
menos que fuerzas externas actúen sobre él.
2. La tasa de variación del momentum de un cuerpo en función del tiempo es proporcional
a la fuerza neta que actúa sobre el cuerpo teniendo la misma dirección de la fuerza,
(entendiéndose por momentum de un objeto al producto de su masa m multiplicado por su
velocidad v).
3. A cada acción existe una reacción igual y opuesta.
La segunda ley nos proporciona una relación importante, conociéndose como la ley de
Newton.
La tasa de cambio o variación en el momentum en el tiempo es así d (mv) /dt.
Si por F entendemos a la fuerza neta que actúa sobre el cuerpo, por la segunda ley tenemos:
d
(mv) α F siendo α el símbolo que indica proporcionalidad. Introduciendo la constante
dt
de proporcionalidad k, obtenemos:
d
(mv ) = kF
dt
= kF o ma = kF, dónde a = dv/dt es la aceleración.
dv
Si m es una constante, m
dt
F=
ma
Así vemos que donde el valor de k depende de las unidades que deseemos usar.
k
Para el sistema CGS (o sistema Centímetro, Gramo, Segundo), k = 1 siendo la ley F = ma.
En la simbología del cálculo podemos escribir las leyes de Newton en formas diferentes, al
notar que la aceleración a puede expresarse como la primera derivada de la velocidad v
(esto es, dv/dt), o como la segunda derivada de v de un desplazamiento s (esto es, d2s/dt2 ).
F =m =m 2
dv d 2s
dt dt
8
Una vez conocido el problema físico podemos aplicar estos conocimientos para obtener las
formulaciones matemáticas de varios problemas de la mecánica clásica que involucran los
conceptos anteriores, y la solución e interpretación de tales problemas.
Ejemplo aclaratorio:
Una masa de m gramos cae verticalmente hacia abajo, bajo la influencia de la gravedad,
partiendo del reposo, siendo despreciable la resistencia del aire. Vamos a establecer la
ecuación diferencial y las condiciones asociadas que describen el movimiento y a
solventarla.
Diagrama de fuerzas:
y P
A t=0
x
mg
Pi t
Tierra
Formulación matemática:
Sea A en la figura la posición de la masa m en el tiempo t = 0, y sea Pi la posición de m en
cualquier tiempo posterior t. En cualquier problema de física que involucre cantidades
vectoriales tales como fuerza, desplazamiento, velocidad y aceleración, las cuales
necesariamente requieren un conocimiento de dirección, es conveniente establecer un
sistema de coordenadas, junto con la asignación de direcciones positivas y negativas. En
este ejemplo observamos que la variación se realiza respecto del eje x.
9
La velocidad instantánea en P es v = dx/dt. La aceleración instantánea en P es a = dv/dt o

a = d2x/dt2. La fuerza que actúa es el peso, siendo su magnitud P= mg. Por la ley de
Newton tenemos:
= mg =g
dv dv
m o
dt dt
Puesto que la masa cae desde el reposo, vemos que v = 0 cuando t = 0, o en otras palabras
v(0) =0 .
Nuestra formulación matemática es el problema de valor inicial v(t)
=g
dv( t )
v(0) =0
dt
Aquí tenemos una ecuación de primer orden y su condición requerida.
Otra manera de formular el problema es escribir:
= mg =g
d 2x d 2x
m o
dt 2 dt 2
En tal caso tenemos una ecuación de segundo orden en las variables x y t, y necesitamos
dos condiciones para determinar x. Una de ellas es v = 0 o dx/dt = 0 en t = 0. La segunda
puede obtenerse al notar que x = 0 en t = 0 (puesto que escogimos el origen de nuestro
sistema de coordenadas en A).
La formulación matemática es:
= g x=0 =0
d 2x dx
y en t = 0
dt 2 dt
Cuando establezcamos ecuaciones diferenciales para describir algún fenómeno o ley,
siempre las acompañaremos de suficientes condiciones necesarias para la determinación de
las constantes arbitrarias en la solución general.
Solución:
= g (separación de variables) obtenemos por integración v = gt + c1.

dv
Empezando con
dt
= gt .
dx
Puesto que v=0 cuando t = 0, c1 = 0, ó v = gt, esto es,
dt
10
Otra integración produce de la anterior ecuación x = gt + c2 . Puesto que x= 0 en t = 0,

1 2
2
c2 = 0. Por tanto x =
1 2
gt . Podríamos haber llegado al mismo resultado al empezar con
2
= mg =g.
d 2x d 2x
m o
dt 2 dt 2
El signo más indica que el objeto se está moviendo en la dirección positiva, esto es, hacia
abajo. Se debería tener en cuenta que si hubiéramos tomado la dirección positiva hacia
= −g o = −g
dv d 2x
arriba la ecuación diferencial hubiera sido m(dv/dt) = - mg, esto es,
dt dt 2
Esto conduciría a resultados equivalentes a los obtenidos.
Para otros problemas similares la forma de actuar es la misma.
2.1 Introducción.
Así como la mecánica tiene como base fundamental las leyes de Newton, la
electricidad también tiene una ley que describe el comportamiento de los circuitos
eléctricos, conocida como la ley de Kirchhoff. Realmente, la teoría de la electricidad está
gobernada por un cierto conjunto de ecuaciones conocidas en la teoría electromagnética
como las ecuaciones de Maxwell. La ley de Kirchhoff es adecuada para estudiar las
propiedades simples de los circuitos eléctricos. El circuito eléctrico más simple es un
circuito en serie, en el cual tenemos una fem (fuerza electromotriz), la cual actúa como una
fuente de energía tal como una batería o generador, y una resistencia, la cual consume o usa
energía, tal como una bombilla eléctrica, tostador, u otro electrodoméstico.
En física elemental encontramos que la fem está relacionada con el flujo de corriente en el
circuito. En forma simple, la ley dice que la corriente instantánea I (en un circuito que
contiene sólo una fem E y una resistencia R) es directamente proporcional a la fem.
Simbólicamente: IαE o I α E de donde, E = IR donde R es una constante de
proporcionalidad llamada el coeficiente de resistencia o simplemente, resistencia. La
ecuación anterior es conocida bajo el nombre de la ley de Ohm.
11
Circuitos más complicados, pero para muchos casos más prácticos, son circuitos que
contienen otros elementos distintos a resistencias. Dos elementos importantes son
inductores y condensadores. Un inductor se opone a cambios en corriente. Tiene un efecto
de inercia en electricidad de la misma manera que una masa tiene un efecto de inercia en
mecánica. Un condensador es un elemento que almacena energía.
En física hablamos de una caída de voltaje a través de un elemento.
En la práctica podemos determinar esta caída de voltaje, o como se llama comúnmente,

caída de potencial o diferencia de potencial, por medio de un instrumento llamado
voltímetro.
Experimentalmente las siguientes leyes se cumplen:
1. La caída de voltaje a través de una resistencia es proporcional a la corriente que

pasa a través de la resistencia.
Si ER, es la caída de voltaje a través de una resistencia e I es la corriente, entonces ER α I o

ER = IR donde R es la constante de proporcionalidad llamada el coeficiente de resistencia o
simplemente resistencia.
2. La caída de voltaje a través de un inductor es proporcional a la tasa de tiempo

instantánea de cambio de la corriente.
Si EL es la caída de voltaje a través del inductor, entonces ELα o EL = L

dI dI
donde L
dt dt
es la constante de proporcionalidad llamada el coeficiente de inductancia o simplemente
inductancia.
3. La caída de voltaje a través de un condensador es proporcional a la carga

eléctrica instantánea en el condensador.
Si Ec es la caída de voltaje a través del condensador y Q la carga instantánea, entonces
Ec =
Q
Ec α Q ó donde hemos tomado 1/C como la constante de proporcionalidad,
C
C se conoce como el coeficiente de capacitancia o simplemente capacitancia.
12
El enunciado es uno de los de la ley de Kirchhoff:
La suma algebraica de todas las caídas de voltaje alrededor de un circuito eléctrico es

cero. [Otra manera de enunciar esto es decir que el voltaje suministrado (fem) es igual a la
suma de las caídas de voltaje.]
Se acostumbra indicar los diferentes elementos de un circuito como se ilustra:
Generador o batería
Resistencia
Inductor
Condensador
Interruptor
Como un ejemplo, considere un circuito eléctrico consistente en una fuente de voltaje E

(batería o generador), una resistencia R, y un inductor L, conectados en serie como se
muestra en la figura:
K
-
E I
Adoptamos la siguiente convención: la corriente fluye del lado positivo (+ ) de la batería o

generador a través del circuito hacia el lado negativo (- ).
13
Puesto que, por la ley de Kirchhoff, la fem suministrada (E) es igual a la caída de voltaje a
dI
través del inductor ( L ) más la caída de voltaje a través de la resistencia (RI), tenemos
dt
como la ecuación diferencial requerida para el circuito:
+ RI = E
dI
L
dt
Como otro ejemplo, suponga que nos dan un circuito eléctrico consistente en una batería o
generador de E en serie con una resistencia de R y un condensador de C.
E I
Aquí la caída de voltaje a través de la resistencia es RI y la caída de voltaje a través del
condensador es Q/C, de modo que por la ley de Kirchhoff RI + = E tal como aparece
Q
C
esto no es una ecuación diferencial. Sin embargo al notar que la corriente es la variación de
la carga con el tiempo, esto es, I = dQ/dt, RI + = E se convierte en R + = E , la

Q dQ Q
C dt C
cual es una ecuación diferencial para la carga instantánea. Acompañando a las ecuaciones
diferenciales obtenidas están las condiciones que se derivan, por supuesto, del problema
específico considerado.
Un generador con una fem se conecta en serie con una resistencia y un inductor. Si el
interruptor K se cierra en tiempo t = 0, establezca una ecuación diferencial para la corriente
y determine la corriente en tiempo t.
14
Llamando a I la corriente o intensidad de corriente que fluye según el primer circuito
descrito, tenemos, por la ley de Kirchhoff, E = RI + L

dI
dt
Puesto que el interruptor se cierra en t = 0, debemos tener I= 0 en t = 0.
Solución:
La ecuación diferencial anterior E = RI + L

dI
es una ecuación de primer orden lineal
dt
exacta; buscando un factor integrante obtenemos μ (t ) = e 2

R
t
. Multiplicado por este factor
R
= RIe + Le =
R R R R t
2
t
2
t
2
t dI 2
t d ( Ie 2 )
la ecuación, da Ee , es decir Ee integrando
dt dt
R
= +c
R t
2
t Ee 2
Ie
10
Puesto que I= 0 en t =0 , podemos con estas condiciones obtener la constante c.
Otro método. La ecuación E = RI + L

dI
puede también resolverse por separación de
dt
variables. Los problemas de este tipo se resuelven todos de la misma forma.
Considere una pieza de material de longitud indefinida acotada por dos planos
paralelos A y B, según la figura, suponiendo que el material es uniforme en todas sus
propiedades, por ejemplo, calor específico, densidad, etc.
15
50ºC 75ºC 90ºC 100ºC

A B
Considerando que los planos A y B se mantienen a 50º C y 100ºC, respectivamente, todo

punto en la región entre A y B alcanza cierta temperatura que no cambia posteriormente.
Así todos los puntos en el plano C entre A y B estarán a 75ºC,y en el plano E a 90ºC.
Cuando la temperatura en cada punto de un cuerpo no varía con el tiempo decimos que
prevalecen las condiciones de estado estacionario o que tenemos un flujo de calor en
estado estacionario.
Como otro ejemplo se considera un tubo de material uniforme, cuyo corte transversal
aparece en la figura.
80ºC
60ºC
40ºC
Se supone que la parte exterior se mantiene a 80ºC y la interna a 40ºC.
Habrá una superficie (línea punteada) en la cual cada punto estará a 60ºC. Sin embargo,
ésta no está en la mitad entre las superficies interna y externa.
Líneas paralelas a A y en un plano perpendicular a A (figura de la pared) se llaman líneas

isotérmicas.
16
La curva punteada del tubo es una curva isotérmica.

Los planos correspondientes de la pared y los cilindros se llaman superficies isotérmicas.
En el caso general, las curvas isotérmicas no serán líneas o círculos, como en los ejemplos,
pero pueden ser una familia de curvas como se muestra en la siguiente figura (curvas
punteadas).
Las trayectorias ortogonales de la familia se llaman líneas de flujo.
Considere pequeñas porciones de dos superficies isotérmicas contiguas separadas por una
distancia Δn.
An
S2
S1
Considerando que la temperatura correspondiente a la superficie S1 es T1 y la

correspondiente a S2 es T2. Llamando a la diferencia de temperatura T1 - T2= ΔT.
Experimentalmente se encuentra que la cantidad de calor que fluye de S1a S2 por unidad de
area y por unidad de tiempo es aproximadamente proporcional a ΔT/Δn. La aproximación
llega a ser más precisa a medida que Δn (y desde luego ΔT) se hace más pequeño. En el
17
caso límite, a medida que ΔnÆ 0, ΔT/ΔnÆdT/dn lo cual se llama el gradiente de T

(variación de T en la dirección normal a la superficie o curva isotérmica). Si q es la
cantidad de flujo de calor por unidad de área y unidad de tiempo, tomamos como nuestra
ley física:
qα ó q=k
dT dT
dn dt
Si deseamos considerar a T como una cantidad vectorial (teniendo dirección y magnitud), el
razonamiento es el siguiente. Considere como positiva la dirección de S1 a S2. Si dT/dn es
positiva, entonces T aumenta y, por tanto, debemos tener T2 < T1. Así, el calor realmente
fluye de S1 a S2 (de mayor a menor temperatura); esto es, el flujo de calor está en la
dirección negativa.
De modo similar, si dT/dn es negativa, T disminuye, T2 > T1, y el flujo es de S2 a S1; esto es,
el flujo de calor está en la dirección positiva.
La dirección del flujo de calor puede tenerse en cuenta mediante un signo menos en
q=k
dT
esto es, cantidad de calor por unidad de tiempo que fluye a través de un área A,
dt
está dada por q (cantidad vectorial) = − kA
dT
(proviene de la teoría de campos).
dn
La anterior constante de proporcionalidad k depende del material usado y se llama
conductividad térmica.
Un tubo largo de acero, de conductividad térmica k, tiene un radio interior ri y un radio

exterior re. La superficie interna se mantiene a Ti y la superficie exterior se mantiene a Te.
(a) Definir la temperatura como una función de la distancia r del eje común de los cilindros
concéntricos.
(b) Encuentre la temperatura en r.
(c) ¿Cuánto calor se pierde por minuto en una parte del tubo de L de largo?
Es claro que las superficies isotérmicas son cilindros concéntricos con los cilindros dados.
El área de tal superficie con radio ri<r<re y longitud L es 2πrL. La distancia dn en este
18
caso es dr. Así, la ecuación puede escribirse q = − k (2πrL)

dT
.En esta ecuación, q es una
dr
constante. Las condiciones iniciales son: k , L, Te en re y Ti en ri.
Solución:
Ecuación diferencial de variables separadas, integrando se obtiene:
-k2πLT=qLnr+c
Usando las condiciones iniciales obtendremos el valor de c.
Otras ideas sobre transmisión de calor:

En el caso de conducción de calor, el calor fluye de lugares de más alta temperatura a
lugares de más baja temperatura. Físicamente, cuando un extremo de una barra aumenta de
temperatura, el movimiento aleatorio de las moléculas en este extremo se aumenta con un
incremento resultante en velocidad y número de colisiones entre moléculas. Las moléculas
con más alta velocidad tienden a moverse hacia el otro extremo de la barra, dando lugar a
colisiones adicionales y a un consecuente incremento gradual en temperatura en el resto de
la barra.
Podemos considerar una conducción de calor como una difusión de moléculas. Tal difusión,
sin embargo, no está limitada a la conducción de calor.
= ay nos dice que la variación con el tiempo de una cantidad y

dy
La ecuación diferencial
dt
es proporcional a y. Si la constante de proporcionalidad a es positiva siendo y positivo,
entonces dy/dt es positivo e y aumenta. En este caso hablamos de que y crece, y el problema
es de crecimiento. Por otro lado, si a es negativo siendo y positivo, entonces dy/dt es
negativo e y decrece. Aquí el problema es uno que involucra decrecimiento.
= ay identificada como una ecuación de variables separadas

dy
Puesto que la solución de
dt
está dada por la función exponencial y = ceat , resolviendo mediante integración,
19
= ay como la ley de crecimiento exponencial si

dy
definiéndose la ecuación diferencial
dt
a > 0 y la ley de decrecimiento exponencial si a <0.
Se calienta agua a la temperatura del punto de ebullición de Te. El agua se agita y se guarda
en un cuarto el cual está a una temperatura constante de Tc. Después de ti la temperatura
del agua es Tf.
(a) Encuentre la temperatura del agua después de tf siendo tf>ti.
(b) ¿Cuándo la temperatura del agua será de Ti? Siendo Te>Tf>Ti>Tc
La diferencia de temperatura entre el agua y el cuarto es Te-Tc=ΔT. La variación en T es

dT/dt. Tomando como base en la experiencia, uno espera que la temperatura cambie más
rápidamente cuando (ΔT) es grande y más lentamente cuando (ΔT) es pequeño.
Desarrollemos un experimento en el cual tomamos temperaturas en varios intervalos de

tiempo, siendo ΔT el cambio en temperatura y Δt el tiempo para producir este cambio.
Tomando a Δt pequeño esperamos que ΔT / Δt será muy cercano a dT/dt.
Si hacemos una gráfica representando ΔT / Δt y ΔT, podríamos producir un gráfico similar

al de esta figura.
20
Los puntos marcados están determinados por el experimento. Puesto que el gráfico es,
aproximadamente, una línea recta, asumimos que dT/dt es proporcional a ΔT, esto es:
= a (ΔT ) donde a es una constante de proporcionalidad. Ahora dT/dt es negativo

dT
dt
cuando (ΔT) es positivo, y así escribiremos a = -h donde h > 0. La ecuación es
= −h(ΔT ) . Esta ecuación se conoce en física como la ley de enfriamiento de Newton y

dT
dt
es de importancia en muchos problemas de temperatura. Realmente, es sólo una
aproximación a la situación física. Las condiciones que acompañan esta ecuación se
obtienen de las condiciones iniciales dispuestas en el enunciado del ejemplo.
Solución:
Resolviendo la ecuación por separación de variables tenemos:
∫ ΔT = − ∫ hdt Æ ln (ΔT)= -ht + c1 Æ ΔT= ce-ht de la cual teniendo las condiciones

dT
iniciales podemos calcular las constantes h y c, pudiendo dar contestación al problema

planteado.
Por métodos experimentales similares a los indicados en el problema anterior de

temperatura obtenemos la siguiente ley:
Ley de desintegración radioactiva: La velocidad de desintegración de una sustancia

radioactiva es proporcional, en cualquier instante, a la cantidad de sustancia que está
presente.
Antes de formular matemáticamente esta ley, consideremos el fenómeno de radioactividad

con algún detalle. Cuando un elemento radioactivo como el radio o el uranio se desintegran,
emiten partículas de una manera aleatoria. Cada una de estas partículas tiene una masa
definida, la cual es pequeña. Si empezamos con una masa de 1 g del material radioactivo y
consideramos lo que sucede cuando se emiten las partículas, encontramos una situación
similar a la que muestra en la figura.
21
x
Pérdida de
partículas
1u
Aquí, x es la cantidad de sustancia que queda después de tiempo t, asumiendo que

empezamos con 1u (unidad) en t = 0. Cada vez que hay una baja en el valor de x significa
que se han emitido partículas; cuanto más grande sea la baja, mayor será el número de
partículas emitidas. Así, la cantidad de la sustancia radioactiva es, en realidad, una función
discontinua en t. Entonces, ¿qué se quiere decir con dx/dt?
Para obviar esta dificultad matemática, aproximamos la curva verdadera por una curva
continua (punteada en la figura de arriba). Así, no cometemos mucho error, y al mismo
tiempo aseguramos tener un gráfico para el cual dx/dt existirá en todo el dominio. Aquí
estamos construyendo una abstracción matemática de una situación física. Las ideas
presentadas aquí ocurren frecuentemente en física debido al tamaño finito, aun de la
partícula más pequeña, en otras palabras, debido a la teoría atómica. Aun en los problemas
de circuitos eléctricos ocurre la abstracción matemática.
Sea A la cantidad de elemento radiactivo presente después de t años. Entonces dA/dt (la
cual existe según la abstracción matemática anterior) representa la tasa o velocidad de
αA ó = αA , donde α es una
dA dA
desintegración del compuesto. De acuerdo con la ley
dt dt
constante de proporcionalidad. Puesto que A > 0 y decreciente, entonces dA/dt<0 y
= − kA . Sea A0 la cantidad del

dA
vemos que α debe ser positiva. Escribiendo α = -k,
dt
compuesto presente inicialmente.
Condiciones iniciales A = A0 en t = 0, A = A1 en t = t1
22
Solución:
Separando las variables, tenemos al integrar:
= − kA Æ ∫ Æ ln A = -kt + c1 Æ antilnÆ A = e − kt +c1

dA
dt
Introduciendo las condiciones iniciales en la función solución de la ecuación diferencial

podremos calcular las constantes k y c1.
Considere un cable o una cuerda que cuelga de dos puntos A y B (según la figura), no
necesariamente al mismo nivel. Suponiendo que el cable es flexible de modo que debido a
su carga (la cual puede ser debida a su propio peso, o a fuerzas externas actuantes, o a una
combinación de éstas) toma la forma como en la figura. Siendo C la posición más baja del
cable, escogiendo los ejes x e y como en la figura, donde el eje y pasa por C.
y
A
P
(x,y)
23
Ø
P
(x,y)
H
C
W (x)
Considere aquella parte del cable entre el punto más bajo y cualquier punto P en el cable
con coordenadas (x, y). Esta parte estará en equilibrio debido a la tensión T en P (según la
figura siguiente), la fuerza horizontal H en C, y la carga vertical total en el segmento CP
del cable denotada por W(x), la cual asumimos que actúa en algún punto Q, no
necesariamente en el centro del arco CP. Para el equilibrio, la suma algebraica de las
fuerzas en la dirección x (u horizontal) debe ser igual a cero, y la suma algebraica de
fuerzas en el eje y (o vertical) debe ser igual a cero. Otra forma de indicar lo mismo es que
la suma de fuerzas hacia la derecha debe ser igual a la suma de las fuerzas hacia la
izquierda, y la suma de fuerzas hacia arriba debe ser igual a la suma de fuerzas hacia abajo.
Descomponemos la tensión T en dos componentes (líneas punteadas en la figura), la

componente horizontal con magnitud TcosØ, y la componente vertical con magnitud
TsenØ.
Las fuerzas en la dirección x son H hacia la izquierda y TcosØ hacia la derecha, mientras
que las fuerzas en la dirección y son W hacia abajo y TsenØ hacia arriba. De donde,
haciendo equilibrio de acciones o fuerzas en las direcciones de los ejes tenemos:
TsenØ = W, TcosØ = H
Dividiendo, y usando el hecho de que la tangente = dy/dx = pendiente en P, tenemos
=
dy W
dx H
En esta ecuación, H es una constante, puesto que es la tensión en el punto más bajo, pero W
puede depender de x. Derivando esta última ecuación con respecto a x, tenemos:
=
d 2 y 1 dW
dx 2 H dx
24
Ahora dW/dx representa el incremento en W por unidad de incremento en x; esto es, la

carga por unidad de distancia en la dirección horizontal. La ecuación diferencial anterior es
fundamental.
Un cable flexible de poco peso (despreciable) soporta un puente uniforme.
Determine la forma del cable. (Este es el problema de determinar la forma del cable en un
puente colgante, el cual es de gran uso en la construcción de puentes).
=
d 2 y 1 dW
La ecuación se cumple aquí y nos resta determinar dW/dx, la carga por
dx 2 H dx
unidad de incremento en la dirección horizontal. En este caso dW/dx es una constante,
llamada el peso por unidad de longitud del puente. Llamando a esta constante W, tenemos
=
d2y W
que
dx 2 H
Denotando por b la distancia del punto más bajo del cable desde el puente, tenemos y = b
= 0 donde x = 0, la segunda condición debido a que el punto donde x = 0

dy
donde x = 0 ,
dx
es un punto mínimo.
25
=
d2y W
Integrando dos veces la ecuación y haciendo uso de las condiciones dadas,
dx 2 H
encontramos que y = + b , siendo esta la ecuación de una parábola.

Wx 2
2H
Considere una viga horizontal AB según la figura. Se supone que la viga es uniforme en su
sección transversal y de material homogéneo. El eje de simetría se encuentra en el plano
medio indica por la zona sombreada.
A B
Cuando está sometida a fuerzas, las cuales suponemos que están en un plano que contiene
el eje de simetría, la viga, debido a su elasticidad, puede distorsionarse en su forma como se
muestra en la siguiente figura.
A B
Estas fuerzas pueden ser debidas al peso de la viga, a cargas aplicadas externamente, o a
una combinación de ambas. El eje de simetría distorsionado resultante, situado en el plano
medio distorsionado de la segunda figura, se llama la curva elástica. La determinación de
esta curva es de importancia en la teoría de elasticidad.
Hay muchas maneras de apoyar vigas.
26
Vigas en voladizo: una viga en la cual el extremo A está rígidamente fijo, mientras que el
extremo B está libre, para moverse.
Viga simplemente apoyada: la viga está apoyada en los dos extremos A y B.
Hay más formas y más condiciones para la deflexión que serán aplicadas a cada tipo de
problema.
Así como hay diferentes maneras de apoyar vigas, también hay diferentes maneras de
aplicar fuerzas de carga externa.
Carga uniformemente distribuida sobre toda la viga.
Carga variable sobre toda la viga o sólo en una parte de ella.
Carga puntual o concentrada.
Considere la viga horizontal OB de la figura siguiente. Colocando el eje de simetría (línea

punteada) en el eje X tomado como positivo a la derecha y con origen en 0. Escoja el eje Y
como positivo hacia abajo.
F1 F2 F3
M(x)
A
B
X
0
x=0 x=L
x
F1
F2
A F3
X
0
y
x=0 x=L
B
x
Debido a la acción de las fuerzas externas F1 y F2 (y si es apreciable el peso de la viga) el

eje de simetría se distorsiona en la curva elástica que se muestra punteada en la figura de
abajo donde hemos tomado la viga como fija en 0. El desplazamiento y de la curva elástica
desde el eje X se llama la deflexión o flecha de la viga en la posición x. Así, si
determinamos la ecuación de la curva elástica, se conocerá la deflexión de la viga.
Para poder formular la ecuación debemos saber:
27
Sea M(x) el momento flector en una sección transversal vertical de la viga en x. Este
momento flector se define como la suma algebraica de los momentos de esas fuerzas que
actúan sobre un lado de x, los momentos se toman sobre una línea horizontal en la sección
transversal en x. Al calcular los momentos adoptaremos la convención de que fuerzas hacia
arriba producen momentos negativos y fuerzas hacia abajo producen momentos positivos,
asumiendo por supuesto que el eje y se toma hacia abajo como se mencionó antes. No
importa cuál lado de x se tome puesto que los momentos flectores calculados desde
cualquier lado son iguales. El momento flector en x está simplemente relacionado con el
radio de curvatura de la curva elástica en x, siendo la relación:
[1 + ( y' ) ]
= M ( x)
y' '
EI
2 32
Donde E es el módulo de elasticidad de Young y depende del material usado en el diseño

de la viga, e I es el momento de inercia de la sección transversal de la viga en x con
respecto a una línea horizontal que pasa por el centro de gravedad de esta sección
transversal. El producto EI se llama la rigidez y se considerará como una constante.
Si asumimos que la viga se dobla sólo levemente, lo cual es válido para muchos propósitos
prácticos, la pendiente y’ de la curva elástica es tan pequeña que su cuadrado es
despreciable comparado con 1, y la ecuación se puede remplazar por la buena
aproximación:
EIy' ' = M ( x)
Una viga horizontal, simplemente apoyada, de longitud L se dobla bajo su propio peso, el
cual es w por unidad de longitud. Encuentre la ecuación de su curva elástica.
En la figura se muestra la curva elástica de la viga (línea punteada) relativa a un conjunto

de ejes coordenados con origen en 0 y direcciones positivas indicadas; puesto que la viga
está simplemente soportada en 0 y en B, cada uno de estos soportes lleva la mitad del peso
de la viga, o sea wL/2.
28
wL/2 wL/2
A
B
X
0
x=0 wx P w(L-x) x=L
x L-x
El momento flector M(x) es la suma algebraica de los momentos de estas fuerzas actuando
a un lado del punto P.
Escogiendo el lado derecho de P, actuarían dos fuerzas:
1. La fuerza hacia abajo w (L - x), a una distancia (L -x)/2 de P, produciendo un momento

positivo.
2. La fuerza hacia arriba wL/2, a una distancia L-x de P, produciendo un momento negativo.
En este caso el momento flector es:
L−x
M ( x) = w( L − x)( ) − ( L − x) = −
wL wx 2 wLx
2 2 2 2
Con el valor de M(x), la ecuación fundamental es:
EIy ' ' = −

wx 2 wLx
2 2
Dos condiciones son necesarias para determinar y. Estas son, y = 0 en x = 0, y en x = L,

puesto que la viga no tiene deformación en los extremos o apoyos.
Solución:
Integrando dos veces EIy ' ' = − se obtiene EIy = − + c1 x + c2

wx 2 wLx wx 4 wLx 3
2 2 24 12
Puesto que y = 0 cuando x = 0, tenemos c2 = 0. De donde EIy = − + c1 x

wx 4 wLx 3
24 12
Puesto que y = 0 cuando x = L, c1 = wL3 /24 y tenemos, finalmente:
29
y= ( x 4 − 2 Lx 3 + L3 x)
w
24 EI
Como la ecuación requerida de la curva elástica. Es de interés práctico usar la solución
final y = ( x 4 − 2 Lx 3 + L3 x) para hallar la máxima deflexión. De la simetría o por el

w
24 EI
cálculo, el máximo ocurre en x = L/2, de donde la flecha máxima será:
ymax =
5wL4
384 EI
30
APLICACIONES DE ECUACIONES DIFERENCIALES LINEALES:
Tipos de aplicaciones:

1.1 El resorte vibrante (movimiento armónico simple).
1.2 El resorte vibrante con amortiguamiento (movimiento amortiguado).
El sistema más simple disponible para estudiar el movimiento vibratorio consiste en

un resorte de peso despreciable [figura(a)] suspendido verticalmente de un soporte fijo.
Suponga que un peso W se cuelga del resorte [figura(b)].
Cuando el peso está en reposo describimos su posición como la posición de equilibrio. Si el

peso se desplaza ejerciendo una fuerza vertical y hacia abajo y se deja de ejercer cuando
tiene una cierta distancia respeto de la posición de equilibrio, estará bajo un movimiento
vibratorio alrededor de la posición de equilibrio [figura(c)].Queremos averiguar el
movimiento que realiza el cuerpo en su desplazamiento respecto de la posición de
equilibrio. Para conseguir este propósito, tendremos que conocer las fuerzas que actúan
sobre el peso durante su movimiento. Por la experiencia vemos que hay una fuerza que
31
hace regresar o restaurar un peso desplazado a su posición de equilibrio. Esta fuerza se

llama la fuerza restauradora.
La ley que gobierna esta fuerza es un caso especial de la ley generalizada de Hooke. Nos
referiremos a este caso especial como la ley de Hooke, la cual se enuncia como sigue:
La fuerza ejercida por un resorte, tendiente a restaurar el peso W a la posición de

equilibrio, es proporcional a la distancia de W a la posición de equilibrio. (“la fuerza es
proporcional al alargamiento”).
Denotamos la magnitud de la fuerza restauradora por | f |, y sea x la posición de W medida

desde la posición de equilibrio. Se supone la dirección positiva hacia abajo, de modo que x
es positivo cuando W está por debajo de la posición de equilibrio y negativo cuando W esté
por encima de esta posición.
De acuerdo a la ley de Hooke,
f α x o f = k x donde k > 0 es una constante de proporcionalidad, que depende de la
dureza del resorte y se llama constante del resorte o módulo de elasticidad del resorte. Para
determinar la dirección de la fuerza, observamos que cuando x > 0 la fuerza está dirigida
hacia arriba y por tanto es negativa. Cuando x < 0 la fuerza está dirigida hacia abajo y es
por tanto positiva. Esto se puede satisfacer, sólo si la fuerza está dada tanto en magnitud
como dirección por - kx, de modo que la ley de Hooke es:
f=-kx
Cuando el peso W se coloca en el resorte, se estira una distancia s como se ve en la anterior

figura en la posición (b), según la ley de Hooke, la tensión T en el resorte es proporcional al
elongamiento, y así T1 = ks; puesto que el resorte y el peso están en equilibrio se tiene que
T1 = ks= W
Cuando el peso se baja más y se suelta, su posición en cualquier tiempo se muestra en la

figura (c). La tensión T, en el resorte en este tiempo es, de acuerdo a la ley de Hooke,
T2 = k(s + x)
Sigue que la fuerza neta en la dirección positiva está dada por:
32
W-T2=W-ks-kx=-kx
Por la ley de Newton, la ecuación diferencial que describe el movimiento es:
= −kx
W d 2x
g dt 2
Así, la fuerza neta es simplemente la fuerza restauradora y no depende del peso W.
Un peso P estira un resorte x1 unidades de longitud. Si el peso se halla x2 unidades de

longitud por debajo de la posición de equilibrio y se suelta:
(a) Establecer la ecuación diferencial y condiciones asociadas que describan el movimiento.
(b) Encontrar la posición del peso como una función del tiempo.
(c) Determinar la posición, velocidad y aceleración del peso t tiempo después de haberse
soltado.
Por la ley de Hooke:
f = k x1 y f =P k=
P
La ecuación diferencial que describe el movimiento
x1
= − kx
P d 2x
es por tanto
g dt 2
Puesto que inicialmente (t = 0) el peso está x2 por debajo de la posición de equilibrio,

tenemos x = x2 en t = 0 .
También, puesto que el peso se suelta (esto es, tiene velocidad cero) en t = 0,
= 0 en
dx
t=0.
dt
33
= − kx con las condiciones

P d 2x
La respuesta a la cuestión (a) está dada por la ecuación
g dt 2
= 0 en
dx
x = x2 en t = 0 , t=0
dt
Solución:
Por el teorema de existencia y unicidad, para ecuaciones lineales, sacamos la función

auxiliar que nos permite llegar a la general.
= −kx llamando a D = , x = e mt y a E=
P d 2x d kg
La ecuación auxiliar sale de
g dt 2 dt P
+ Ex = 0 y D 2 = m 2e mt .
d 2x
tenemos 2
dt
+ Ex = 0 la anterior tendremos: ( D 2 + E )e mt = 0 de donde sacamos

d 2x
Sustituyendo en 2
dt
que (m 2 + E )e mt = 0 para la cual m2 + E = 0 y tiene raíces m =± Ei. De donde la
ecuación diferencial tiene la solución de la forma: x = c1e Eit + c2e − Eit
De otra forma podemos calcular la solución general de esta ecuación haciendo a x´= v en
+ Ex = 0 , + Ex = 0 , + Ex = 0 ,
d 2x dv dv dx
la ecuación x´´ + Ex=0 tendremos que para 2
dt dt dx dt
+ Ex = 0 separado variables e integrando tendremos ∫vdv+E∫xdx=0

dv
v entonces
dx
v + Ex = c de donde c2 = 2c despejando v = ± − Ex
1 2 1 2 2 2
c como x´= v
2 2
⎛x⎞
∫ = ± ∫ dt con lo cual acsen⎜⎜ ⎟⎟ = ± Et + c2
dx
c 2 − Ex 2 ⎝ c1 ⎠
integrando en la anterior
x = c1 (± Et + c2 ) = c1senc2 cos Et ± c1 cos c2 senEt que se puede escribir como :
x = A sen Et + B cos Et siendo A y B constantes.
34
Como la ecuación x = A sen Et + B cos Et, y la ecuación x = c1e Eit + c2e − Eit ,son las dos
ecuaciones generales, debe haber una correlación entre ambas de tal forma que:
A sen Et + B cos Et ≡ c1e Eit + c2e − Eit si consideramos que en t = 0 e ± Eit = 1 tendremos
B = c1 + c2
[ ]
diferenciando en ambos lados de la ecuación
A sen Et + B cos Et ≡ c1e Eit + c2e − Eit Æ E(A cos Et − B sen Et ) ≡ Ei c1e Eit − c2e − Eit
E(A ) ≡ Ei[c1 − c2 ] con lo que A = i[c1 − c2 ]
para
y, B = c1 + c2 , A = i[c1 − c2 ] ,tendremos:
t=0 con los valores obtenidos antes de A
c1 (cosEt + isen Et) + c2 (cosEt − isen Et) ≡ c1e Eit + c2e − Eit con lo cual e Eit = cosEt + isen Et .
De la condición x = x2 en t = 0, encontramos A, así que sustituyendo en
=0
dx
x = A cos Et + B sen Et y derivando respecto de t e introduciendo la condición
dt
en t = 0 podemos encontrar el valor de B y con lo cual la x para las condiciones dadas.
Para calcular la velocidad y aceleración solo habrá que derivar una vez para la velocidad y
dos para la aceleración, respecto del tiempo, en la ecuación de la posición:
x = A sen Et + B cos Et
= v(t ) = a (t )
dx d 2x
dt dt 2
En el gráfico de la figura siguiente se obtiene dando valores a la función posición,

x = A sen Et + B cos Et, se ve que el peso empieza en x = A donde t = 0, luego pasa por la
posición de equilibrio en x =0 a la posición x =-A y luego regresa de nuevo a la posición
de equilibrio, la pasa y vuelve a x = A. Este ciclo se repite una y otra vez. Físicamente el
gráfico describe el movimiento periódico hacia arriba y abajo del resorte, el cual se llama
movimiento armónico simple.
En general, cualquier movimiento descrito por la ecuación 2 = −ax ,

d 2x
dt
35
donde a > 0 es una constante, será movimiento armónico simple. Físicamente esta ecuación
nos dice que la aceleración es directamente proporcional al desplazamiento pero en
dirección opuesta (indicado por el signo menos).
Llamamos al desplazamiento máximo del peso de su posición de equilibrio (esto es, x = 0)

amplitud.
El tiempo para un ciclo completo se llama periodo. Del gráfico vemos que el período es
π π
= T ; otra forma de ver que el período es ,sin el gráfico, es determinando cuándo el
4 4
peso está en un extremo de su trayectoria (esto es, ya sea el punto más alto, o más bajo).
Cuando tomamos el punto más bajo dado por x = A cos Et, esto ocurrirá cuando cos Et =
1, esto es, Et = 0, 2 π, 4 π, 6π, . . o t = 0; π /4, 2 π /4, 3 π /4, . . .la primera vez que x =
max es cuando t = 0, la segunda cuando t = π /4, la tercera cuando t = 2 π /4, etc. La
diferencia entre tiempos sucesivos es π /4, la cual es el período. El número de ciclos por
segundos se llama frecuencia.
Período (T) = Número unidades de tiempo por ciclo = π /4
=
π π
1 4
Frecuencia (f) = Número de ciclos por unidad de tiempo =
f =
1
En general si T es el período, la frecuencia f está dada por
T
Si ahora de la solución general x = A cos Et + B sen Et damos las condiciones iniciales:
≠0 x≠0
dx
en t=0.
dt
36
Despejando las constantes A y B obtenemos una ecuación de la forma:
x = a cos Et + b sen Et llamando a E = ω tendremos x = a cos ω t + b sen ω t esta

ecuación la podemos poner de la forma:
a cos ω t + bsen ω t = a 2 + b 2 sen (ω t + φ )
obtenida de a 2 + b 2 sen(ωt + φ ) = a 2 + b 2 ( senωt cos φ + cos ωtsenφ ) =
⎡ ⎛ ⎞ ⎛ ⎞⎤
= a 2 + b 2 ⎢( senωt )⎜⎜ ⎟ + (cos ωt )⎜
⎟ ⎜ 2
⎟⎥ = a cos ωt + bsenωt
⎟
b a
⎢⎣ ⎝ a +b ⎠ ⎝ a +b ⎠⎥⎦
2 2 2
El ángulo ø a menudo se llama el ángulo de fase.
Si un movimiento se puede describir de la forma:
1 ω 2π
x = Asen(ωt + φ ) f = = T=
T 2π ω
entonces
obteniendo la relación ω = 2πf la cual a menudo es útil.
El movimiento armónico simple se da en muchos otros casos además de las vibraciones de

resortes, como en el movimiento del péndulo de un reloj, el balanceo de un barco o un
avión, etc.
37
1.2 El resorte vibrante con amortiguamiento (movimiento amortiguado).
Los resortes considerados en el apartado anterior son modelos idealizados fuera de los
casos reales, puesto que las oscilaciones no disminuían, como uno esperaría por
experiencia, sino por el contrario, se mantenían constantes en el tiempo. En la práctica,
fuerzas de fricción y otras (tales como la resistencia del aire) actúan para decrecer la
amplitud de las oscilaciones y finalmente traer el sistema al reposo. Una manera para
obtener una mejor aproximación a la realidad, es suponiendo una fuerza amortiguadora. La
ley exacta para esta fuerza no se conoce, puesto que depende de muchas variables, pero se
ha encontrado mediante la experimentación, que para velocidades pequeñas, la magnitud de
la fuerza amortiguadora es aproximadamente proporcional a la velocidad instantánea del
peso en el resorte.
La magnitud por tanto está dada por:
β
dx
dt
donde β es la constante de proporcionalidad llamada la constante de amortiguamiento.
La fuerza amortiguadora se opone al movimiento, de modo, que cuando el peso va bajando

la fuerza amortiguadora actúa hacia arriba, mientras que actúa hacia abajo cuando el peso
va subiendo. Asumiendo hacia abajo como la dirección positiva, vemos que la fuerza
amortiguadora debe ser negativa cuando dx/dt sea positiva, y debe ser positiva cuando dx/dt
sea negativa. Así, con β > 0, es claro que la fuerza amortiguadora debe estar dada tanto en
magnitud como en dirección por -β dx/dt .
Cuando se tienen en cuenta las fuerzas restauradoras ya consideradas, según la ley de

Newton, la ecuación diferencial del movimiento es:
= −β − kx +β + kx = 0
W d 2x dx W d 2x dx
2
o 2
g dt dt g dt dt
38
Un peso W estira un resorte x1 unidades de longitud. Si el peso se halla x2 unidades de

longitud por debajo de la posición de equilibrio y se suelta:
Teniendo en cuenta la fuerza amortiguadora
+β + kx = 0
W d 2x dx
La ecuación del movimiento será 2
g dt dt
= 0 en
dx
Las condiciones iniciales son: x = x2 en t = 0 , t=0
dt
Solución:
Por el teorema de existencia y unicidad para ecuaciones lineales sacamos la función

auxiliar que nos permite llegar a la general.
+β + kx = 0 llamando a D = , x = e mt y a
W d 2x dx d
La ecuación auxiliar sale de 2
g dt dt dt
βg
E= F= +F + Ex = 0 y D = me mt D 2 = m 2e mt .
kg d 2x dx
tenemos 2
W W dt dt
+F + Ex = 0 la anterior tendremos: ( D 2 + DF + E )e mt = 0 de
d 2x dx
Sustituyendo en 2
dt dt
donde sacamos que (m 2 + mF + E )e mt = 0 para la cual (m 2 + mF + E ) = 0 y tiene raíces
39
− F ± F 2 − 4E
m= De donde la ecuación diferencial tiene la solución de la forma:
2
− F + F 2 − 4E − F − F 2 − 4E
x = c1e m1t + c2e − m2t donde m1 = m2 =
[ ]
2 2
De otra forma m = −a ± bi si F 2 − 4 E < 0 x = c1et ( − a+bi ) + c2et ( − a−bi ) = e − ta c1etbi + c2e −tbi
[ ]
siendo la expresión x = c1etbi + c2e − tbi cuya solución general debe ser como la ecuación
x = A sen bt + B cos bt , son las dos ecuaciones generales. Debe haber una correlación
entre ambas, de tal forma que:
A sen bt + B cos bt ≡ c1etbi + c2e −tbi si consideramos que en t=0 e ±tbi = 1 tendremos
B = c1 + c2 diferenciando en ambos lados de la ecuación A sen bt + B cos bt ≡ c1etbi + c2e −tbi
Æ b(A cos bt − B sen bt ) ≡ bi[c1ebit − c2e −bit ] para b(A ) ≡ bi[c1 − c2 ] con lo que
A = i[c1 − c2 ] A = i[c1 − c2 ] ,
t=0
con los valores obtenidos antes de A y B , B = c1 + c2 ,
tendremos , c1 (cosbt + isen bt) + c2 (cosbt − isen bt) ≡ c1ebit − c2e − bit con lo cual
ebit = cosbt + isen bt .
De la condición x = x2 en t = 0, encontramos A, así que sustituyendo en
x = A sen bt + B cos bt y derivando respecto de t e introduciendo la condición = 0 en

dx
dt
t = 0 podemos encontrar el valor de B y con lo cual la x para las condiciones dadas.
Por analogía de las ecuaciones:
A sen bt + B cos bt ≡ c1etbi + c2e −tbi
[
A sen bt + B cos bt ≡ e − ta c1etbi + c2e − tbi ]
tendremos que la segunda ecuación se convierte en x = e −ta c1etbi + c2e − tbi [ ] donde su
ecuación solución será:
x = e −ta [A cos bt + Bsenbt ]
40
Determinando las constantes A y B sujetas a las condiciones citadas anteriormente,

encontramos la ecuación de la posición en función del tiempo, que determina una gráfica de
este tipo:
Si hacemos uso de la identidad a cos ω t + bsen ω t = a 2 + b 2 sen (ω t + φ )
y de la ecuación x = e −ta [A cos bt + Bsenbt ] se puede escribir:
x = e − ta c 2 + d 2 [sen(ωt + φ )]
El gráfico de la figura está entre los gráficos de líneas punteadas puesto que el seno varía
entre - 1 y + 1. El punto P representa un mínimo relativo, mientras que el punto Q está
sobre la curva punteada. La diferencia constante en tiempos entre máximos (o mínimos)
sucesivos se llama cuasi período, aunque se puede referenciar éste como el período. El
adjetivo cuasi se usa, puesto que los valores funcionales no se repiten como lo harían si
realmente hubiera periodicidad.
El movimiento descrito en este ejemplo se llama movimiento oscilatorio amortiguado o

movimiento amortiguado.
Se debería notar que x = e − ta c 2 + d 2 [sen(ωt + φ )] tiene la forma
x = α (t )[sen(ωt + φ )]
41
2π
El cuasi período está dado por: T =
ω
Por analogía con el caso no amortiguado, α(t) se llama la amplitud, o más exactamente la
amplitud tiempo variante. Se ve que la amplitud decrece con el tiempo, estando así de
acuerdo con lo experimentado. Un hecho destacado es que la frecuencia con
amortiguamiento es menor que aquella sin amortiguamiento, esto es posible ya que se
esperaría oposición al movimiento para incrementar el tiempo para un ciclo completo. La
frecuencia sin amortiguamiento, esto es, con β = 0, a menudo se llama la frecuencia
natural. Esta es de gran importancia en conexión con el fenómeno de resonancia a ser
discutido posteriormente.
La fuerza amortiguadora puede ser demasiado grande comparada con la fuerza restauradora
para permitir el movimiento oscilatorio.
En los casos anteriores vimos el problema de un resorte donde solo se consideraron las
fuerzas restauradora y amortiguadora. Consideramos ahora casos donde pueden actuar otras
fuerzas externas que varían con el tiempo. Tales fuerzas pueden ocurrir, por ejemplo,
cuando el soporte que sostiene el resorte se mueve arriba y abajo en una manera
especificada tal como en movimiento periódico, o cuando al peso se le da un pequeño
empuje cada vez que alcanza la posición más baja. Si denotamos la fuerza externa por F(t),
la ecuación diferencial para el movimiento es:
= −β − kx + F (t ) +β + kx = F (t )
W d 2x dx W d 2x dx
2
o 2
g dt dt g dt dt
+ b + cx = F (t )
d 2x dx
donde la última ecuación puede escribirse: a 2
dt dt
(donde a = W/g, b= β, c = k), a menudo llamada la ecuación de vibraciones forzadas.
42
Un peso W estira un resorte x1 unidades de longitud. Si el peso se halla x2 unidades de

longitud por debajo de la posición de equilibrio y se suelta actuando sobre él una fuerza
externa periódica, dada por la ecuación F (t ) = A cos Bt :
La ecuación diferencial es:
+ b + cx = A cos Bt
d 2x dx
a 2
dt dt
= 0 en
dx
las condiciones iniciales son: x = x2 en t = 0 , t=0
dt
Solución:
Estudiando la primera parte de la ecuación diferencial
0=a + b + cx
d 2x dx
2
dt dt
escribiéndola en función de D = , x = e mt , D 2 = m 2e mt se vería como

d
dt
(aD 2 + bD + c)e mt = 0 o (aD 2 + bD + c) x = 0 donde en está ultima ecuación a,b,c,

serían constantes no todas iguales a cero de tal forma que usando la definición de
dependencia lineal:
43
Un conjunto de funciones y1 (x), y2 (x),…, yn (x) denotadas como y1, y2,…, yn , se dice que
son linealmente dependientes en un intervalo si existe un conjunto de n constantes, no todas
cero, tales que en el intervalo α1 y1 + α 2 y2 + ... + α n yn ≡ 0 en caso contrario se dice que el
conjunto es linealmente independiente.
La dependencia o independencia se averigua mediante el Wronskiano.
Teorema: Si y1, y2 son linealmente dependientes en un intervalo, y si sus derivadas y´1, y´2
existen en el intervalo, entonces el Wronskiano de y1, y2 dado por
W ( y1 , y2 ) = det y1 , y2 , y´1 , y´2 = y1 y´2 − y2 y´1
es idénticamente cero (esto es, W= 0) en el intervalo.
Este teorema puede enunciarse en términos de independencia lineal como sigue:
Si el Wronskiano de y1, y2 no es idénticamente cero (esto es, W≠ 0) en un intervalo,

entonces y1, y2 son linealmente independientes en el intervalo.
Esto es cierto puesto que si y1, y2 fueran linealmente dependientes en el intervalo su

Wronskiano sería idénticamente cero en el intervalo según el teorema anterior. Esta
contradicción muestra que las funciones no son linealmente independientes, esto es, son
linealmente dependientes en el intervalo. Con lo cual, si W = y1 y´2 − y2 y´1 ≠ 0 tenemos
α1 = 0 α 2 = 0 .
Retomando la ecuación (aD 2 + bD + c) x = 0 y por el teorema de existencia-unicidad para

n = 2, tomando como soluciones de la ecuación a
(ax1´´+bx1´+cx1 ) = 0 y (ax2´´+bx2´+cx2 ) = 0
multiplicando a la primera ecuación por x2, a la segunda por x1, restando las dos tendremos
a( x1 x2´´− x2 x1´´) + b( x1 x2´− x2 x1´) = 0
Puesto que el Wronskiano está dado por W = x1 x´2 − x2 x´1 y ya que
44
= ( x1 x´2 − x2 x´1 ) = x1 x´´2 − x2 x´´1 se convierte en a + bW = 0

dW d dW
dt dt dt
Resolviendo, esta ecuación diferencial de variables separadas
dW − b −b − ∫ tdt +c1 − ∫ tdt − ∫ tdt

= dt Æ integrando Æ lW = t + c1 Æ W = e a = e a e c1 = ce a
b b b
W a a
relación importante conocida como la identidad de Abel, dada por W = ce−tdt
∫ a tdt
Puesto que la función exponencial en W = ce
−
b
nunca es cero, vemos que el Wronskiano
W debe ser idénticamente cero en el intervalo dado, en cuyo caso c = 0, o nunca cero en el
intervalo, en cuyo caso c ≠ 0. No puede haber nada intermedio.
Con todo lo anteriormente dicho, tendremos el siguiente teorema:
Teorema: Sea x1,x2 dos soluciones de la ecuación diferencial (aD 2 + bD + c) x = 0
donde a ≠ 0, b,c son funciones continuas de x en algún intervalo. Entonces el Wronskiano

de x1,x2 está dado por la identidad de Abel.
W = det x1 , x2 , x´1 , x´2 = ce − tdt y W es ya sea idénticamente cero en el intervalo o nunca
cero en el intervalo.
Usando este teorema, podemos probar ahora el siguiente teorema para el caso de que x1,x2
sean soluciones de la ecuación (aD 2 + bD + c) x = 0
Teorema: Sean x1,x2 soluciones de la ecuación diferencial (aD 2 + bD + c) x = 0 en algún

intervalo. Entonces
(a) x1,x2 son linealmente dependientes si y sólo si W= 0 en el intervalo.
(b) x1,x2 son linealmente independientes si y sólo si W≠ 0 en el intervalo.
De los dos Teoremas anteriores podemos obtener el siguiente teorema útil:
Teorema: Sea x1 una solución a la ecuación diferencial (aD 2 + bD + c) x = 0 . Entonces
45
una solución linealmente independiente está dada por W = det x1 , x2 , x´1 , x´2 = ce − tdt Æ
− ∫ tdt − ∫ tdt
Æ ( )+k = ∫
b b
x1 Æ ( )= dt Æ
2 d x2 ce a x2 ce a
dividendo en los dos termino por
dt x1 x12 x1 x12
∫ a tdt
integrando y tomando la constante de integración c = 1, k = 0 Æ x2 = x1 ∫
−
b
e
dt
x12
Otros teoremas importante nos dicen:
-Si x1,x2 son soluciones linealmente independientes de la ecuación (aD 2 + bD + c) x = 0 ,
entonces x = c1 x1 + c2 x2 es una solución de la ecuación para cualesquiera constantes c1 ,c2 .
-Si x1,x2 son soluciones linealmente independientes de (aD 2 + bD + c) x = 0 y xp es una
solución particular de (aD 2 + bD + c) x = F (t ) entonces x = c1 x1 + c2 x2 + x p es una solución
de la ecuación diferencial (aD 2 + bD + c) x = F (t ) para cualesquiera constantes c1,c2 .
De estos dos últimos teoremas se puede desprender:
Se ha definido la solución general de una ecuación diferencial de orden n como aquella

solución que involucra n constantes arbitrarias. Todas las soluciones que no se podían
obtener de esta solución general por ninguna selección de estas constantes se denominan
soluciones singulares. En la mayoría de los problemas de naturaleza práctica la solución
general es la que proporciona la solución significativa después de determinar las constantes
de las condiciones dadas. Si surgen soluciones singulares, estas tienen “generalmente” poco
o ningún significado práctico. De estos teoremas vemos que no hemos encontrado la
solución que involucra dos constantes arbitrarias (el mismo número como el orden de la
ecuación diferencial), la cual hemos llamado la solución general, sino en realidad todas las
otras soluciones que de existir son soluciones particulares, esto es, casos especiales de la
solución general obtenidos por la selección apropiada de las constantes. Así, estos
Teoremas garantizan que no haya soluciones singulares.
Esta propiedad de no tener soluciones singulares es peculiar de las ecuaciones diferenciales

lineales, pero no de todas las ecuaciones diferenciales no lineales.
46
Como F (t ) = A cos Bt contiene una función de los términos cosBt con lo que podemos
utilizar el método de los coeficientes indeterminados para determinar la solución particular.
Ensayando como solución particular x p = c 1 senBt + c 2 cos Bt de
(aD 2 + bD + c) x = F (t )
x p = c1senBt + c2 cos Bt x´ p = Bc1 cos Bt − Bc2 senBt x´´p = − B 2 c1senBt − B 2 c2 cos Bt
Sustituyendo en la ecuación ax´´+b´x´+cx = F (t ) calculamos el valor de las constantes

c1 ,c2 y por tanto la solución particular.
La solución complementaria se obtiene de la misma forma que en el apartado anterior 1.2

con lo que sería
x c = e − ta [ A cos Bt + BsenBt ] con lo que la solución general será de la forma:
x = xc + x p = e − ta [ A cos Bt + BsenBt ] + c1senBt + c2 cos Bt
Usando las condiciones iniciales podríamos determinar el gráfico siguiente:
47
Cuando la frecuencia de una fuerza externa periódica aplicada a un sistema mecánico está
relacionada de una manera sencilla con la frecuencia natural del sistema, puede ocurrir el
fenómeno de la resonancia mecánica, la cual eleva las oscilaciones tales que el sistema
puede desplomarse o colapsarse. Una compañía de soldados marchando en fila a través de
un puente puede de esta manera hacer que el puente colapse. En una manera análoga, puede
ser posible que una nota musical de una frecuencia característica propia estalle un cristal.
Debido a los grandes daños que pueden ocurrir, la resonancia mecánica es en general algo
que necesita ser evitado, especialmente por el ingeniero al diseñar estructuras o
mecanismos vibrantes.
La resonancia mecánica puede también servir para propósitos útiles. Por ejemplo, si un
automóvil se quedara atascado en la nieve (o barro), puede, al “mecerlo”, ser puesto a
vibrar con su frecuencia natural. Entonces, al aplicarle una fuerza con esta misma
frecuencia, la resonancia mecánica resultante puede ser algunas veces suficiente para
liberarlo.
Un peso P estira un resorte x1 unidades de longitud. Si el peso se halla x2 unidades de

longitud por debajo de la posición de equilibrio y se suelta aplicándole una fuerza externa
dada por F (t ) = A cos Bt .
Describir el movimiento que resulta si se supone que inicialmente el peso está en la

posición de equilibrio (x = 0) y que su velocidad inicial es cero.
La ecuación diferencial que describe el movimiento es por tanto
= − kx + F (t )
P d 2x
g dt 2
48
Las Ecuaciones Diferenciales ysus Aplicaciones en la Ingeniería
Puesto que inicialmente (t = 0) el peso está x2 por debajo de la posición de equilibrio,

tenemos x = x2 en t = 0 .
También, puesto que el peso se suelta (esto es, tiene velocidad cero) en t = 0,
= 0 en
dx
t=0.
dt
Solución:
Como F (t ) = A cos Bt contiene una función de los términos cosBt, con lo que podemos
utilizar el método de los coeficientes indeterminados para determinar la solución particular.
Ensayando como solución particular x p = c 1 senBt + c 2 cos Bt de
D + k ) x = F (t )
P 2
(
g
x p = c1senBt + c2 cos Bt x´ p = Bc1 cos Bt − Bc2 senBt x´´p = − B 2 c1senBt − B 2 c2 cos Bt
x´´+ Kx = F (t ) calculamos el valor de las constantes c1 ,c2 y

P
Sustituyendo en la ecuación
g
por tanto la solución particular.
= −kx + F (t ) se obtiene de la misma forma que en el

P d 2x
La solución complementaria de
g dt 2
apartado anterior 1.1 con lo que sería xc = A cos Bt + BsenBt con lo que la solución
general será de la forma:
x = xc + x p = A cos Bt + BsenBt + c1senBt + c2 cos Bt
Usando las condiciones iniciales, podemos observar una solución será cuando A = B = 0 de
donde podemos sacar que:
x = c1senBt + c2 cos Bt
podríamos determinar el gráfico siguiente:
49
El gráfico de la ecuación x = c1senBt + c2 cos Bt está entre los gráficos de x = t y x = - t

como se muestra en la figura. Se puede ver en el gráfico que las oscilaciones van creciendo
sin límite. Naturalmente, el resorte está limitado a romperse dentro de un corto tiempo.
En este ejemplo, el amortiguamiento fue ignorado, ocurriendo la resonancia porque la

frecuencia de la fuerza externa aplicada fue igual a la frecuencia natural del sistema no
amortiguado. Esto es un principio general.
En el caso donde ocurre amortiguamiento las oscilaciones no crecen sin límite, pero sin
embargo pueden llegar a ser muy grandes, la resonancia en este caso ocurre cuando la
frecuencia de la fuerza externa aplicada es ligeramente menor que la frecuencia natural del
sistema.
50

Las Ecuaciones Diferenciales y Sus Aplic

Cargado por

Copyright:

Formatos disponibles

Las Ecuaciones Diferenciales y Sus Aplic

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

Las Ecuaciones Diferenciales y Sus Aplic

Cargado por

Copyright:

Formatos disponibles

LAS ECUACIONES DIFERENCIALES

-Etapas de resolución del problema científico. Pg (5)

-Aplicaciones de ecuaciones diferenciales de primer orden y simples de

-Aplicaciones de ecuaciones diferenciales lineales. Pg (31-50)

1. Movimiento vibratorio de sistemas mecánicos:

LAS ECUACIONES DIFERENCIALES

El descubrimiento de Newton y Leibniz en el siglo diecisiete sobre las ideas básicas

= f ´(x) , en forma de diferencial de una función de una

tenemos que F ( x) = ∫ f (t )dt si a ≤ x ≤ b , existe entonces F´(x) en cada punto x del

intervalo abierto (a,b), en el que f es continua, y para tal x tenemos F´(x) = f ( x)

-La Derivada de la Integral de una función es la propia función:

= f ´(x) mediante la integración de dicha función, que es lo que

-Se llama orden de la ecuación diferencial al orden de la derivada o derivada parcial

-Se dice que una ecuación diferencial es lineal si tiene la forma

-Se define como problema de valor inicial y problemas de valor frontera a

Problema de valor inicial:

Problemas de valor frontera:

-La función primitiva resultante, o función solución de una ecuación diferencial,

ETAPAS DE RESOLUCIÓN DEL PROBLEMA CIENTÍFICO:

1) FORMULACIÓN MATEMÁTICA DEL PROBLEMA CIENTÍFICO:

Las leyes científicas, que por supuesto están basadas en experimentos u

2) SOLUCIÓN DE LAS ECUACIONES:

Las ecuaciones formuladas en la etapa anterior necesitan ser resueltas, sujetas a

3) INTERPRETACIÓN CIENTÍFICA DE LA SOLUCIÓN:

APLICACIONES DE ECUACIONES DIFERENCIALES DE PRIMER ORDEN Y

1.2 Las leyes del movimiento de Newton.

1. Un cuerpo en reposo tiende a permanecer en reposo, mientras que un cuerpo en

La velocidad instantánea en P es v = dx/dt. La aceleración instantánea en P es a = dv/dt o

= g (separación de variables) obtenemos por integración v = gt + c1.

Otra integración produce de la anterior ecuación x = gt + c2 . Puesto que x= 0 en t = 0,

2. Aplicaciones a los circuitos eléctricos:

En física hablamos de una caída de voltaje a través de un elemento.

En la práctica podemos determinar esta caída de voltaje, o como se llama comúnmente,

Experimentalmente las siguientes leyes se cumplen:

1. La caída de voltaje a través de una resistencia es proporcional a la corriente que

Si ER, es la caída de voltaje a través de una resistencia e I es la corriente, entonces ER α I o

2. La caída de voltaje a través de un inductor es proporcional a la tasa de tiempo

Si EL es la caída de voltaje a través del inductor, entonces ELα o EL = L

3. La caída de voltaje a través de un condensador es proporcional a la carga

Si Ec es la caída de voltaje a través del condensador y Q la carga instantánea, entonces

2.2 La ley de Kirchhoff.

El enunciado es uno de los de la ley de Kirchhoff:

La suma algebraica de todas las caídas de voltaje alrededor de un circuito eléctrico es

Se acostumbra indicar los diferentes elementos de un circuito como se ilustra:

Como un ejemplo, considere un circuito eléctrico consistente en una fuente de voltaje E

Adoptamos la siguiente convención: la corriente fluye del lado positivo (+ ) de la batería o

Aquí la caída de voltaje a través de la resistencia es RI y la caída de voltaje a través del

la carga con el tiempo, esto es, I = dQ/dt, RI + = E se convierte en R + = E , la

Llamando a I la corriente o intensidad de corriente que fluye según el primer circuito

descrito, tenemos, por la ley de Kirchhoff, E = RI + L

Puesto que el interruptor se cierra en t = 0, debemos tener I= 0 en t = 0.

La ecuación diferencial anterior E = RI + L

exacta; buscando un factor integrante obtenemos μ (t ) = e 2

Puesto que I= 0 en t =0 , podemos con estas condiciones obtener la constante c.

Otro método. La ecuación E = RI + L

3. Aplicaciones a flujo de calor en estado estacionario.

50ºC 75ºC 90ºC 100ºC

Considerando que los planos A y B se mantienen a 50º C y 100ºC, respectivamente, todo

Se supone que la parte exterior se mantiene a 80ºC y la interna a 40ºC.