UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA

FACULTAD DE INGENIERÍA DE PROCESOS

ESCUELA PROFESIONAL DE INGENIERÍA METALÚRGICA

ECUACIONES DIFERENCIALES – GRUPO A

TRABAJO DE INVESTIGACIÓN FORMATIVA

BÚSQUEDA Y ANÁLISIS DE PROBLEMAS DE APLICACIÓN DEL CURSO DE

ECUACIONES DIFERENCIALES

TRABAJO PRESENTADA POR:

Jimenez Cardenas Percy Santiago

Ccama Bravo, Víctor Alexis
Cerpa Quicaña Sharon Nicool
Choquemamani Mansilla Alexander
Humani Flores Alexander Rodrigo
Meza Villanueva Luis Karlo
Quispe Hancco Yuliño

ASESOR: PARISACA ZAIRA, OSCAR LEONIDAS

AREQUIPA – PERÚ

JULIO – 2023
 1. INTRODUCCIÓN:
Las ecuaciones diferenciales son herramientas matemáticas fundamentales en el campo de la
ingeniería. Estas ecuaciones describen el comportamiento de sistemas dinámicos y son
utilizadas para modelar una variedad de fenómenos físicos y procesos en los diferentes
campos de la ingenierí[Link] de las aplicaciones más comunes de las ecuaciones diferenciales
en ingeniería es en la descripción y análisis del comportamiento de circuitos eléctricos. Las
leyes de Kirchhoff, que son ecuaciones diferenciales lineales, permiten modelar y resolver
circuitos eléctricos complejos, tanto en corriente continua como en corriente alterna.

Además, las ecuaciones diferenciales también se utilizan en el análisis de sistemas mecánicos.

Por ejemplo, en ingeniería civil, se aplican para modelar y predecir el comportamiento de
estructuras bajo cargas variables, como puentes, edificios y aerogeneradores. En ingeniería
aeroespacial, estas ecuaciones son fundamentales en el análisis del movimiento de naves
espaciales y cohetes.

En el área de Ingeniería Metalúrgica, las ecuaciones diferenciales también juegan un papel

importante en el análisis y diseño de diferentes procesos y fenómenos. Algunas de las
aplicaciones comunes de las ecuaciones diferenciales en la Ingeniería Metalúrgica son las
siguientes , Las ecuaciones diferenciales se utilizan para modelar y analizar la transferencia de
calor en diferentes procesos metalúrgicos, como la fundición, el temple, el enfriamiento y el
calentamiento de metales. Estas ecuaciones permiten predecir y controlar el gradiente de
temperatura, calcular los perfiles térmicos y optimizar los parámetros de diseño para
garantizar la calidad del producto final.
En la fabricación de materiales metálicos, las ecuaciones diferenciales se emplean para
modelar el crecimiento de cristales durante la solidificación de aleaciones. Estas ecuaciones
ayudan a determinar la velocidad de crecimiento de los granos y la formación de estructuras
cristalinas, lo que influye en las propiedades mecánicas y térmicas de los [Link]én
se utilizan para describir la cinética de las reacciones químicas en procesos como la refinación
de metales, la producción de acero, la recristalización y el tratamiento térmico. Estas
ecuaciones permiten predecir cómo evolucionan las concentraciones de los diferentes
elementos o compuestos a medida que ocurren las reacciones.
Las ecuaciones diferenciales también se emplean en la modelación y análisis del transporte de
masa en diferentes procesos metalúrgicos, como la lixiviación, extracción, electroobtención y
precipitación. Estas ecuaciones permiten predecir la difusión de especies químicas en
soluciones y determinar los perfiles de concentración, lo que es esencial para la optimización
de los procesos de separación y purificación.

Las ecuaciones diferenciales desempeñan un papel fundamental en diversas áreas de la

ingeniería, permitiendo modelar, analizar y predecir el comportamiento de sistemas dinámicos
complejos. Su aplicación se extiende desde el diseño de circuitos eléctricos hasta el análisis
de sistemas mecánicos, químicos y de control, lo que demuestra la importancia de esta
herramienta matemática en el campo de la ingeniería.
 ÍNDICE

1. INTRODUCCIÓN
2. OBJETIVOS
2.2. OBJETIVO GENERAL
2.3. OBJETIVOS ESPECÍFICOS
3. PROBLEMAS DE APLICACIÓN DEL CURSO DE ECUACIONES
DIFERENCIALES A INGENIERÍA METALÚRGICA
3.1. LEYES DEL MOVIMIENTO DE NEWTON
3.1.1 Ejercicio 1
3.2. CIRCUITOS ELÉCTRICOS
3.2.1 Ejercicio 2
3.3. FLUJOS VOLUMÉTRICO DE SALIDA Y ENTRADA IGUALES
3.3.1. Ejercicio 3
3.4. DEFLEXIÓN DE VIGAS
3.4.1. Ejercicio 4
3.5. Conducción de Calor en una Aleta de Enfriamiento (EDO de 2do Orden)
3.5.1. Ejercicio 5
3.5.2. Ejercicio 6
3.6. Movimiento Armónico de Resorte Vibrante
3.6.1. Ejercicio 7
3.7. Flujo de Calor en Estado Estacionario
3.7.1. Ejercicio 8
4. Conclusiones
5. Bibliografía
2. OBJETIVOS

2.2. OBJETIVO GENERAL

- Aplicar los conocimientos y técnicas de las ecuaciones diferenciales en el ámbito

de la ingeniería metalúrgica para resolver problemas complejos, permitiendo el
análisis y la optimización de procesos metalúrgicos, la predicción de fenómenos
físicos relevantes y la toma de decisiones fundamentadas en base a modelos
matemáticos, con el fin de mejorar la eficiencia, la calidad y la seguridad de los
procesos de producción metalúrgica.

2.3. OBJETIVOS ESPECÍFICOS

- Aplicar los conceptos y métodos de las ecuaciones diferenciales para modelar y

resolver problemas reales en la ingeniería metalúrgica, como la transferencia de
calor en procesos metalúrgicos, la cinética de reacciones químicas en la
metalurgia, la difusión de especies en sólidos, etc.

- Analizar y resolver problemas específicos de la ingeniería metalúrgica utilizando

técnicas y métodos de ecuaciones diferenciales, como ecuaciones diferenciales
lineales de primer y segundo orden, ecuaciones diferenciales parciales, sistemas
de ecuaciones diferenciales, entre otros.

- Desarrollar habilidades para la formulación y resolución de problemas propios de

la ingeniería metalúrgica mediante la aplicación de técnicas de ecuaciones
diferenciales, fomentando así el pensamiento crítico y analítico.

- Comunicar de manera efectiva los resultados obtenidos al resolver problemas de

aplicación del curso de ecuaciones diferenciales en la ingeniería metalúrgica,
utilizando lenguaje técnico adecuado y representaciones gráficas cuando sea
necesario.
PROBLEMAS DE APLICACIÓN DEL CURSO DE ECUACIONES DIFERENCIALES
A INGENIERÍA METALÚRGICA

3.1 LEYES DEL MOVIMIENTO DE NEWTON

Las tres leyes del movimiento primeras desarrolladas por Newton son:
1. Primera ley inercia: 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 esta.
2. Segunda ley Fundamental de la Dinámica: La tasa de cambio en movimiento de un
cuerpo en el tiempo es proporcional a la fuerza neta que actúa sobre el cuerpo y tiene
la misma dirección a la fuerza.
3. Tercera ley acción y reacción: A cada acción existe una reacción igual y opuesta.
La segunda ley nos proporciona una relación importante conocida a los estudiantes de
física elemental y nos referiremos a ella brevemente como la ley de Newton. El
movimiento de un objeto se define como su masa “ m ” multiplicada por su velocidad
La tasa de cambio en movimiento en el tiempo es así 𝑑/𝑑𝑡 (𝑚. 𝑎) . Si denotamos por
la fuerza neta que actúa sobre el cuerpo la segunda ley dice que:
𝑑/𝑑𝑡 (𝑚. 𝑎) 𝑘 𝐹
Si m es constante:
𝑚 𝑑𝑣/𝑑𝑡 = 𝑘 𝐹 ó 𝑚𝑎 = 𝑘 𝐹
Y es la aceleración donde vemos que a = dv / dt
𝐹 = 𝑚𝑎/𝑘
El valor de k , depende de las unidades que deseemos usar. Hasta el momento se usan
dos sistemas principales.
3.1.1 Ejercicio :
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:
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 involucra
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.
La velocidad instantánea en P es𝑣 = 𝑑𝑥/𝑑𝑡 . La aceleración instantánea en P es
𝑎 = 𝑑𝑣/𝑑𝑡 o
2 2
𝑎 = 𝑑 𝑥/𝑑𝑡 . La fuerza que actúa es el peso, siendo su magnitud 𝑃 = 𝑚𝑔. Por la
ley de Newton tenemos:
𝑚 𝑑𝑣/𝑑𝑡 = 𝑚𝑔 𝑜 𝑑𝑣/𝑑𝑡 = 𝑔
Puesto que la masa cae desde el reposo, vemos que v = 0 cuando t = 0, o en otras
palabras
𝑉(0) = 0
Nuestra formulación matemática es el problema de valor inicial 𝑉(𝑡)
𝑑𝑣(𝑡)/𝑑𝑡 = 𝑔 𝑣(0) = 0
Aquí tenemos una ecuación de primer orden y su condición requerida. Otra manera
de formular el problema es escribir:
2 2 2 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:

2 2
𝑑 𝑥/𝑑𝑡 = 𝑔 𝑥 = 0 𝑦 𝑑𝑥/𝑑𝑡 = 0 𝑒𝑛 𝑡 =0

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:
Empezando con 𝑑𝑣/𝑑𝑡 = 𝑔 (separación de variables) obtenemos por integración
𝑣 = 𝑔𝑡 + 𝑐1
Puesto que 𝑣 = 0 cuando 𝑡 = 0, 𝐶1 = 0, ó 𝑣 = 𝑔𝑡 , 𝑒𝑠𝑡𝑜 𝑒𝑠, 𝑑𝑥/𝑑𝑡 = 𝑔𝑡
2
Otra integración se produce de la anterior ecuación 𝑥 = 1/2𝑔𝑡 + 𝐶2 . Puesto que
2
𝑥 = 𝑡 𝑒𝑛 𝑡 = 0 𝐶2 = 0 . 𝑃𝑜𝑟 𝑡𝑎𝑛𝑡𝑜 𝑥 = 1/2𝑔𝑡 . Podríamos haber llegado al
mismo resultado al empezar con
2 2 2 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 arriba la ecuación diferencial hubiera sido
2 2
𝑚(𝑑𝑣/𝑑𝑡) =− 𝑚𝑔, 𝑒𝑠𝑡𝑜 𝑒𝑠, 𝑑𝑣/𝑑𝑡 = − 𝑔 𝑜 𝑑 𝑥/𝑑𝑥 = −𝑔
Esto conduciría a resultados equivalentes a los obtenidos. Para otros problemas
similares la forma de actuar es la misma.

3.2 CIRCUITOS ELÉCTRICOS:

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.
a. 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:
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:

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 (- ).
Puesto que, por la ley de Kirchhoff, la fem suministrada (E) es igual a la
caída de voltaje a través del inductor ( 𝐿 𝑑𝑙/𝑑𝑡 ) más la caída de voltaje a
través de la resistencia (RI), tenemos como la ecuación diferencial requerida
para el circuito:
𝐿 𝑑𝑙/𝑑𝑡 + 𝑅𝐼 = 𝐸
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.
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
𝑅𝐼 + 𝑄/𝐶 = 𝐸 tal como aparece 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
𝐼 = 𝑑𝑄/ 𝑑𝑡 𝑅𝐼 + 𝑄/𝐶 = 𝐸 se convierte en
𝑅 = 𝑑𝑄/𝑑𝑡 + 𝑄/𝐶 = 𝐸 la 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.

3.2.1 Ejercicio :
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..

Formulación matemática:

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

circuito descrito, tenemos, por la ley de Kirchhoff, 𝐸 = 𝑅𝐼 + 𝐿 𝑑𝐼/𝑑𝑡
Puesto que el interruptor se cierra en t = 0, debemos tener I= 0 en t = 0.

Solución:
La ecuación diferencial anterior 𝐸 = 𝑅𝐼 + 𝐿 𝑑𝐼/𝑑𝑡 es una ecuación de
primer orden lineal exacta; buscando un factor integrante obtenemos
 𝑅/2
𝑢(𝑡) = 𝑒 . Multiplicado por este factor la ecuación, da
𝑡 𝑅/2 𝑡 𝑅/2 𝑡 𝑅/2
𝐸𝑒 = 𝑅𝐼𝑒 + 𝐿𝑒 𝑑𝐼/𝑑𝑡, es decir, e integrando en c:
𝑡𝑅/2 𝑡𝑅/2
𝐼𝑒 = 𝐸𝑒 /10 + 𝑐
Puesto que I= 0 en t =0 , podemos con estas condiciones obtener la constante
c. Otro método.
La ecuación 𝐸 = 𝑅𝐼 + 𝐿 𝑑𝐼/𝑑𝑡 puede también resolverse por separación de
variables. Los problemas de este tipo se resuelven todos de la misma forma.

3.3 FLUJOS VOLUMÉTRICO DE SALIDA Y ENTRADA IGUALES - 1er Orden

Tasa de flujo volumétrico:

La tasa de flujo volumétrico de un fluido es el volumen de fluido que pasa a través de una
sección transversal dada por unidad de tiempo. El término sección transversal es solo una
forma elegante de describir el área a través de la que algo fluye, por ejemplo, el área circular
dentro de la línea punteada en el diagrama que se muestra a continuación.

Como la tasa de flujo volumétrico mide la cantidad de volumen que pasa a través de un área
en un tiempo dado, su ecuación se ve así:
𝑉 𝑉𝑜𝑙𝑢𝑚𝑒𝑛
𝑄= 𝑡
= 𝑡𝑖𝑒𝑚𝑝𝑜
En unidades del SI (Sistema Internacional de Unidades), la tasa de flujo volumétrico tiene
unidades de metros cúbicos por segundo. Para ello se desarrollará el siguiente ejercicio de
aplicación. En la primera parte se analizan las razones de entrada y
salida de masa en el tanque para flujos volumétricos distintos, desarrollando como solución
una ecuación algebraica solucionable en un instante de tiempo dentro de un intervalo
acotado,con la restricción de que el flujo de entrada es mayor al de salida por lo que el tanque
comenzará a derramar solución cuando se llene.
​
3.3.1 Ejercicio de aplicación:
Si se considera un tanque de mezclado por agitación uniforme,que en un tiempo inicial
t=0 contiene una solución de una sustancia conocida (salmuera) como el que se muestra en la
figura 1, y que se está llenando de la misma solución,pero con una concentración 𝐶1

constante y a un flujo volumétrico 𝑉1 igualmente constante. Se puede determinar la

cantidad de masa m(t) de soluto en el tanque en cualquier instante de tiempo

conociendo la concentración 𝐶2y el flujo volumétrico de salida 𝑉2.

El análisis comienza con una suposición simple acerca de que la masa de soluto en el tanque
𝑚(𝑡) es igual a la diferencia de la masa de soluto que entra 𝑚1 y la masa 𝑚2 que sale del
tanque, es decir:

𝑚(𝑡) = 𝑚1 − 𝑚2 (1)

El incremento ∆𝑚 durante un intervalo de tiempo ∆𝑡 se puede expresar como:

Si se considera el límite cuando ∆𝑡 ⟶ 0 lo que indica la razón de cambio instantánea de

𝑚(𝑡), la Ec.(2) se convierte en:

𝑑𝑚(𝑡) 𝑑𝑚1 𝑑𝑚2

𝑑𝑡
= 𝑑𝑡
- 𝑑𝑡
(3)

La Ec. (3) expresa que la razón de entrada de soluto menos la razón de salida de soluto es
igual a la razón de cambio de soluto en el tanque, lo cual expone la hipótesis inicial. En el
desarrollo de la Ec. (3) primero se analizará la razón de entrada y posteriormente la razón de
salida de la solución.
Si la concentración de una solución es igual a la masa del soluto por unidad de volumen de
disolvente, la masa 𝑚1se obtiene mediante la siguiente relación:

𝑚1 = 𝑐1𝑉1 (4)

Y la razón de cambio en el tiempo:

𝑑𝑚1 𝑑𝑐1𝑉1
𝑑𝑡
= 𝑑𝑡
(5)

Debido a que 𝑐1 se considera un parámetro constante en el análisis y que el flujo volumétrico

𝑉1 es:
𝑑𝑐1(𝑡)
𝑑𝑡
= 𝑉1 (6)

La Ec. (5) se reduce a:

𝑎𝑚2
𝑑𝑡
=𝑐1𝑉1 (7)

Recordando que 𝑐1 y 𝑉1 se consideran como datos constantes en el análisis. De

manera similar para la razón de salida se tiene:

𝑑𝑚2
𝑑𝑡
= 𝑐2𝑉2 (8)

Se hace la suposición de que la concentración de salida 𝑐2 de la solución es igual a la

concentración en el tanque, por lo que 𝑐2 se puede expresar como:

𝑐2(𝑡) = 𝑚(𝑡)/𝑉𝑇 (9)

Donde 𝑉𝑇 , 𝑐2(𝑡) y 𝑚(𝑡) son el volumen de solución en el tanque, la concentración y la masa

del soluto en cualquier instante de tiempo 𝑡 en el tanque respectivamente. El volumen de
solución en el tanque 𝑉𝑇 es igual al volumen inicial 𝑉0 más el volumen que entra y menos el
volumen que sale, es decir:

𝑉𝑇 = 𝑉0 + 𝑉1 − 𝑉2 (10)

Recordando que el flujo volumétrico constante es:

𝑉 = 𝑉̇ 𝑡 (11)
La ecuación (11) se convierte en:

𝑉𝑇 = 𝑉0 + (𝑉2 − 𝑉1)𝑡 (12)

Sustituyendo a la Ec. (12) en la Ec. (9):

𝑚(𝑡)
𝑐2(𝑡) = ( 𝑉 +𝑉 −𝑉 )𝑡 (13)
0 1 2

Y la Ec. (13) usando la relación de la Ec. (8):

𝑑𝑚2 𝑉2
𝑑𝑡
= 𝑉0+(𝑉1−𝑉2)𝑡
𝑚(𝑡) (14)

Finalmente reemplazando los valores de la Ec. (3) por los obtenidos en las Ecs. (7) y (14):

𝑑𝑚(𝑡) 𝑉2
𝑑𝑡
=𝑐1𝑉2- 𝑉0+(𝑉1−𝑉2)𝑡 m(t) (15)

La Ec. (15) se puede reescribir como:

𝑑𝑚(𝑡) 𝑉
𝑑𝑡
+ 𝑉 +(𝑐 2−𝑉 )𝑡 𝑚(𝑡) = 𝑐1𝑉2 (16)
0 1 2

La Ec. (16) es una ecuación diferencial ordinaria lineal de primer orden, primer grado, con
coeficientes constantes que modela el sistema de mezclado del tanque de la figura 1
expresando la cantidad de soluto en el tanque en unidad de masa.

La solución de la Ec. (16) se obtiene determinando un factor integrante apropiado, el cual es:

Multiplicado la Ec. (16) por el factor integrante obtenido en la Ec. (17):

𝑉2
𝑉1− 𝑉2
[𝑉0 + (𝑉1 − 𝑉2 )𝑡] 𝑚(𝑡) =
𝑉2
𝑉1− 𝑉2
∫[𝑉0 + (𝑉1 − 𝑉2)𝑡] . 𝑐1𝑉2𝑑𝑡 (18)

La solución de la integral en el segundo miembro de la Ec. (18), aplicando las fórmulas

básicas de integración, es:
 𝑉2
𝑉1− 𝑉2
∫[𝑉0 + (𝑉1 − 𝑉2)𝑡] . 𝑐1𝑉2𝑑𝑡 =

𝑉
2
𝑉 −𝑉
+1
1 2
𝑐1𝑉1 .[𝑉0 +(𝑉1− 𝑉2)𝑡]
𝑉 (19)
(𝑉1− 𝑉2)( 𝑉 −2𝑉 +1)
1 2

Simplificando la Ec. (19):

𝑉
2
𝑉 −𝑉
1 2

∫[𝑉0 + (𝑉1 − 𝑉2)𝑡] . 𝑐1𝑉1𝑑𝑡 =

𝑉
2
𝑉 −𝑉
1 2

𝑐1[ 𝑉0 + (𝑉1 − 𝑉2)𝑡] + 𝑐1 (20)

Sustituyendo a la Ec. (20) en la Ec. (18) y resolviendo para 𝑚(𝑡):

𝑉
2
𝑉 −𝑉
1 2

𝑐1[ 𝑉0 +(𝑉1− 𝑉2)𝑡] + 𝐶1

𝑚(𝑡) = 𝑉
2
(21)
𝑉 −𝑉
1 2

[ 𝑉0 +(𝑉1− 𝑉2)𝑡]

𝑐1[ 𝑉0 +(𝑉1− 𝑉2)𝑡]

𝑚(𝑡) = 𝑉
2
(22)
𝑉 −𝑉
1 2

[ 𝑉0 +(𝑉1− 𝑉2)𝑡]

Donde 𝐶1 en la Ec. (22) es una constante de integración y se determina considerando que en

un tiempo inicial 𝑡 =0 la masa 𝑚(𝑡) es igual a la masa inicial 𝑚0 en el tanque, la condición
inicial del sistema se expresa como:

𝑡 = 0 ⟶ 𝑚(𝑡) = 𝑚0 (23)

Sustituyendo la relación de la Ec. (23) en la Ec. (22):

𝑐1
𝑚0 = 𝑐1𝑉1 + 𝑉
2
(24)
𝑉 −𝑉
1 2
𝑉0

El valor de la constante 𝐶 es:

 𝑉2

𝑐1 = (𝑚0 − 𝑐1𝑉0 )𝑉0 𝑉1− 𝑉2

(25)

Sustituyendo a la Ec. (25) en la Ec. (22):

𝑚(𝑡) = 𝑐1[𝑉0 + (𝑉1 − 𝑉2 )𝑡] (26)

𝑉
2
𝑉 −𝑉
(𝑚0−𝑐1𝑉0)𝑉0 1 2

+ 𝑉
2
𝑉 −𝑉
1 2
[𝑉0+ (𝑉1− 𝑉2 )𝑡]

Simplificando en la Ec. (26):

𝑚(𝑡) = 𝑐1[𝑉0 + (𝑉1 − 𝑉2 )𝑡]

𝑉1
𝑉0 𝑉1− 𝑉2
+ (𝑚0 − 𝑐1𝑉0)[ 𝑉 +(𝑉 −𝑉 )𝑡 ] (27)
0 1 2

Con la Ec. (27) se puede calcular la masa del soluto 𝑚(𝑡) en el tanque de mezclado (figura 1)
en un instante de tiempo 𝑡 considerando los flujos volumétricos de entrada y salida
[Link] se considera 𝑉𝑚á𝑥 como el volumen máximo de solución que puede tener el
tanque antes de que se desborde, entonces el tiempo en el que la solución se derramará se
puede obtener mediante la Ec. (12) considerando 𝑉1> 𝑉2̇ :

𝑉𝑚á𝑥 = 𝑉0 + (𝑉1 − 𝑉2 )𝑡

𝑉𝑚𝑎𝑥−𝑉0
𝑡= 𝑉1 − 𝑉2
(28)

Sustituyendo a la Ec. (28) en la Ec. (27):

𝑚(𝑡) = 𝑐1𝑉𝑚á𝑥
𝑉1
𝑉0 𝑉1− 𝑉2
+ (𝑚0 − 𝑐1𝑉0 ) ( 𝑉𝑚á𝑥 ) (29)
La Ec. (29) se puede utilizar para conocer la cantidad de soluto cuando el tanque está
completamente lleno. La concentración en el tanque 𝑐𝑇 utilizando la Ec. (13) es:
𝑚(𝑡)
𝐶𝑡 = 𝑉𝑚𝑎𝑥
(30)

Y sustituyendo la Ec. (29) en la Ec. (30):

𝑉
1
𝑉0 𝑉 −𝑉
1 2
(𝑚0 − 𝑐1𝑉0)( 𝑉𝑚𝑎𝑥 )
𝐶𝑡=𝑐1+ 𝑉𝑚𝑎𝑥
(31)

Flujos volumétricos de salida y entrada iguales

Para el caso en el que los flujos de salida y entrada son iguales y constantes la Ec. Diferencial
(16) se expresa como:

𝑑𝑚(𝑡) 𝑉
𝑑𝑡
+ 𝑉2 𝑚(𝑡) = 𝐶1𝑉1 (32)
0

Y siguiendo el mismo método de solución, el factor de integración es:

𝑉2 𝑉2
𝑒∫ 𝑉0
𝑑𝑡 = 𝑒∫ 𝑉0
𝑡 (33)

Por lo que la solución de la Ec. (33) se expresa como:

𝑉2 𝑉2
𝑉0 𝑉0
𝑒 𝑚(𝑡) = ∫𝑒 . 𝑐1𝑉1𝑑𝑡 (34)

Resolviendo la Ec. (34):

𝑉2 𝑉
𝑉0
𝑉0𝑒 𝑉2 .𝐶1𝑉1
𝑒 𝑚(𝑡) = 0

𝑉2
+𝐶2 (35)

La constante 𝐶2 se calcula utilizando la condición inicial de la Ec. (23) como:

𝑉0.𝐶1𝑉1
𝑚0 = 𝑉2
+ 𝐶2

𝑉0.𝐶1𝑉1
𝐶2 = 𝑚0 − 𝑉2
Sustituyendo a la Ec. (36) en la Ec. (35):

𝑉2 𝑉2
𝑉0 𝑉0.𝐶1𝑉1 𝑉0
𝑒 𝑚(𝑡) = 𝑉2
(𝑒 -1)+𝑚0 (37)

Resolviendo la Ec. (37) para 𝑚(𝑡) y considerando que 𝑉1= 𝑉2̇

𝑉2 𝑉
𝑉0
− 𝑉2 𝑡
𝑚(𝑡) =𝐶1𝑉0(1 −𝑒 )+𝑚0𝑒 0
(38)

RESPUESTA: Con la Ec. (38) se puede calcular la masa del soluto 𝑚(𝑡) en el tanque
de mezclado (figura 1) en cualquier instante de tiempo 𝑡 considerando los flujos
volumétricos de entrada y salida iguales.

3.4 LA DEFLEXIÓN DE VIGAS

Considere una viga horizontal AB de la Figura 3.49(a). Se asume que la viga es

uniforme en su sección transversal y de material homogéneo. El eje de simetría se
indica por la línea punteada. Cuando está sometida a fuerzas, las cuales se asumen 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 Figura 3.49(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, punteado en la
Figura 3.49(b), se llama la curva elástica. La determinación de esta curva es de
importancia en la teoría de elasticidad y será parte del propósito de esta sección
mostrar cómo se hace esto.
Sea M(x) el momento flexionante en una sección transversal vertical de la viga en x.
Este momento flexionante 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. Como se mostrará en el Ejemplo ilustrativo 1, no
importa cuál lado de x se tome puesto que los momentos flexionantes calculados desde
cualquier lado son [Link] el momento flexionante en x está simplemente
relacionado con el radio de curvatura de la curva elástica en X, siendo la relación:
"
𝑦
𝐸𝐿 3 =M(𝑥) (1)
[1+(𝑦´)

donde E es el módulo de elasticidad de Young y depende del material usado en el

diseño de la viga, e Z 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 EZ se llama la rigidez flexural, 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 (1) se puede reemplazar por
la buena aproximación.

Ely” = M(x)

3.4.1 Ejercicio de aplicación :

Encontrar la curva elástica de una viga en voladizo uniforme de longitud L con un peso
constante 𝑤 por unidad de longitud y determine la deflexión en el extremo libre.
Formulación matemática:

La curva punteada en la Figura 3.53 es la curva elástica de la viga en voladizo. El

origen 0 del sistema de coordenadas se toma en el extremo fijo, y los ejes positivos z y
y como se muestran. Para calcular M(x) es más sencillo considerar la parte de la viga a
la derecha de P, puesto que sólo una fuerza actúa, a saber, la fuerza hacia abajo w(L -
x1, produciendo un momento positivo dado por:
2 2
𝐾 𝑤(𝐿−𝑥) 𝑤(𝐿−𝑥)
𝑀(𝑥) = 𝑤(𝐿 − 𝑥(𝐿 − 2
)= 2
de modo que 𝐸𝐼y”= 2

el cual debemos resolver sujeto a las condiciones y = y’= 6 donde x = 6,

puesto que no hay deflexión en x = 0, y puesto que la pendiente de la tangente a la

curva elástica en x = 0 es cero.

Integrando dos veces y usando las condiciones es fáciL mostrar que:

𝑤 4 2 2 2
y= 24𝐸𝐼 (𝑥 − 4𝐿𝑥 + 6𝐿 𝑥 )

Es la ecuación de la curva elástica. Colocando x = L, encontramos deflexión en el

extremo libre:
2
𝑤𝐿
= 8𝐸𝐼

3.5 CONDUCCIÓN DE CALOR EN UNA ALETA DE ENFRIAMIENTO(EDO - 2DO

ORDEN)

Observamos a diario que existen muchas aplicaciones sobre todo en el ámbito de enfriamiento
de una pared en un proceso metalúrgico. Por eso mismo existen parámetros como el aire que
rodea mediante un flujo una pared plana. Entre ellos se buscó opciones rentables para
solucionar y la más adecuada fue el uso de aletas que se extienden desde la pared hacia el
fluido de los alrededores.

En la figura se indica un esquema de una aleta rectangular sencilla.

Indicamos los siguientes parámetros más importantes como:

-La densidad de flujo de calor q, en la superficie de la aleta (en la dirección de z), viene dada
por q = h(T − Ta), donde h, coeficiente de transmisión de calor convectivo entre la superficie
de la aleta y el fluido, es constante y uniforme sobre la superficie.

-Para determinar la transferencia de calor asociada a una aleta debemos primero obtener la
distribución de temperatura T(z) a lo largo de la misma.

-Consideremos un segmento diferencial ∆z de la aleta como el que se muestra en la figura. En

este, se presentan dos tipos diferentes de transferencia de calor: una por conducción (a través
de las superficies transversales del segmento) y otra por convección (a través de la superficie
lateral del segmento).

Resumiendo obtendremos lo siguiente:

2
𝑑𝑇 ℎ𝑃
2 = 𝑘𝐴𝑡
(𝑇 − 𝑇𝑎)
𝑑𝑧

Donde k es la conductividad térmica del material con que está hecha la aleta y que se
considera constante; AT es el área de la sección transversal de la aleta (2BD en nuestro caso)
y P es el per´ımetro de la sección transversal (2D + 4B para nuestra aleta). Para el tipo de
 ℎ𝑃 ℎ
aleta que estamos analizando, donde B << D, el factor 𝑘𝐴𝑇
puede aproximarse por 𝑘𝐵
.

Teniendo en cuenta esto, la ecuación resulta

2
𝑑𝑇 ℎ
2 = 𝑘𝐵
(𝑇 − 𝑇𝑎)
𝑑𝑧

Realizamos un cambio de variable tomando que Θ = 𝑇 − 𝑇𝑎, obtendremos

2
𝑑Θ 2
2 =𝑟 Θ
𝑑𝑧

2 ℎ
donde 𝑟 = 𝑘𝐵 . Resolviendo esta ecuación y teniendo en cuenta que Θ = T − Ta la solución

general de la ecuación se expresa.

𝑟𝑧 −𝑟𝑧
𝑇(𝑧) = 𝐶1𝑒 + 𝐶2𝑒 + 𝑇𝑎

Las constantes 𝐶1 y 𝐶2 se calculan con la ayuda de condiciones aceptables. La condición T(0)

= 𝑇𝑤 no tiene discusión ya que hemos supuesto conocer la temperatura Tw de la pared. En z =

𝑑𝑇
L, puede considerarse 𝑑𝑧
(L) = 0 lo cual significa que estamos despreciando las pérdidas en el

extremo de la aleta debido a su delgadez. Ası, obtenemos

𝑇𝑤−𝑇𝑎
𝐶1 = 2𝑟𝐿
𝑒 +1

𝑇𝑤−𝑇𝑎 2𝑟𝐿
𝐶2 = 2𝑟𝐿 𝑒
𝑒 +1

El calor total transferido por la aleta, QT , se obtiene calculando la integral

𝐿
𝑄𝑇 = ∫ ℎ𝑃(𝑇 − 𝑇𝑎)𝑑𝑧
0

Tomando en cuenta todo lo anteriormente dicho, podemos resolver problemas aplicativos, por
ejemplo tendremos el siguiente:
3.5.1 Ejercicio de aplicación

1.-Imaginemos una aleta rectangular con las siguientes dimensiones: longitud igual a 35 mm,
ancho igual a 30 mm y espesor igual a 2 mm. La conductividad térmica de la aleta es k = 205
𝑊 𝑊
𝑚°𝐶
, el coeficiente de transferencia de calor es igual a 600 2 . ◦C , la temperatura en la
𝑚 °𝐶

base es igual a 135°C y la temperatura del aire ambiente es igual a 40°C. Calculemos el calor
disipado por la aleta
𝐿
𝑟𝑧 −𝑟𝑧
𝑄𝑇 = ℎ𝑃 ∫(𝐶1𝑒 + 𝐶2𝑒 )𝑑𝑧𝑑𝑦
0

Sabemos:

ℎ 600 1 1
𝑟= 𝑘𝐵
= 205*0.001 𝑚
= 54. 10 𝑚

P = 2D + 2B = 0.064m

𝑇𝑤−𝑇𝑎 135−40
𝐶1 = 2𝑛𝐿 = 45.12
°𝐶 = 2. 11°𝐶
𝑒 +1

2𝑛𝐿
𝐶2 = 𝑒 𝐶1 = 92. 66°𝐶

Integrando y sustituyendo valores se obtiene

ℎ𝑃 𝑟𝐿 −𝑟𝐿 ℎ𝑃
𝑄𝑇 = 𝑇
(𝐶1𝑒 − 𝐶2𝑒 − 𝐶1 + 𝐶2) = 𝑇
(𝐶2 − 𝐶1

𝑄𝑇 = 62. 93 𝑊

Por último, para comparar situaciones, calculemos cuál hubiese sido el valor de QT si no se
hubiera colocado la aleta. El calor disipado por la porción rectangular de la pared, de área A,
correspondiente a la base que ocupa la aleta se calcula con la expresión hA(Tw − Ta). En
nuestro ejemplo A = 2mm·30mm. Entonces, si la aleta no estuviera, la porción de superficie
de la pared ocupada por ella disipará solo 600 · 0, 00006 · 95W = 3, 42W.

Aplicaciones de E. D. de primer orden

A veces las ecuaciones diferenciales de primer orden se escriben en la forma

M(x, y) dx + N(x, y) dy = 0.

Por ejemplo, si suponemos que y representa la variable dependiente en (y − x) dx + 4x dy = 0

entonces y′ =
dy/dx, y al «dividir» ambos miembros de la ecuación por dx obtenemos la expresión
equivalente 4xy′ + y = x.
Definición
En general una ecuación diferencial de primer orden adopta la forma
𝑑𝑦
𝑓(𝑥, 𝑦, 𝑦) = 0, 𝑜 𝑏𝑖𝑒𝑛 𝑑𝑥
= 𝑓(𝑥, 𝑦)

Dónde ,en ambos casos𝑦 = 𝑦(𝑥) , la primera de las ecuaciones anteriores denota una
ecuación diferencial en forma implícita ,mientras que la segunda se dice que está en forma
explícita.
observación
En la aplicación es habitual considerar como variable independiente t en lugar de X así las
ecuaciones anteriores toman la forma
𝑑𝑦
𝑓(𝑡, 𝑦, 𝑦) = 0, 𝑜 𝑏𝑖𝑒𝑛 𝑑𝑥
= 𝑓(𝑡, 𝑦)

donde 𝑥 = 𝑥(𝑡)
definición
Se llama solución de una ecuación diferencial al primer orden a una función derivable con
derivada continua que alser sustitución en la ecuación
𝑑𝑦
𝑓(𝑥, 𝑦, 𝑦) = 0, 𝑜 𝑏𝑖𝑒𝑛 𝑑𝑥

3.5.2 Una solución de salmuera de sal fluye a razón constante de 6L/min. hacia el interior de
un depósito que inicialmente contiene 50L de solución de salmuera en la cual se disolvieron 5
kg de sal. La solución contenida en el depósito se mantiene bien agitada y fluye hacia el
exterior con la misma rapidez. Si la concentración de sal en la salmuera que entra en el
depósito es de 0.5kg/L, determinar la cantidad de sal presente en el depósito al cabo de t
minutos. ¿Cuándo alcanzará la concentración de sal en el depósito el valor de 0.3kg/L?

SOLUCIÓN

𝑑𝑥
𝑑𝑡
= 4(0, 2) − 3( 𝑥
100 + 𝑡 )
𝑥(0) = 0

3
∫ 100 + 𝑑𝑡
𝑑𝑥
𝑑𝑡
+ 3( 𝑥
100 + 𝑡 ) = 0, 8 µ(𝑡) = 𝑒 = (100 + 𝑡)
3

4
3 3 (100+𝑡)
(100 + 𝑡) 𝑥(𝑡) = 0, 8∫ (100 + 𝑡) 𝑑𝑡 + 𝑐 = 0, 8 4
+ 𝑐

𝑐
𝑥(𝑡) = 0, 2(100 + 𝑡) + 3
(100+𝑡)

𝑐 3 7
𝑥(𝑡)= 0 → 0 = 20 + ⇒ 𝑐 =− 20𝑥100 =− 2𝑥10
1003

7 7
2𝑥10 𝑥(𝑡) 2𝑥10
𝑥(𝑡) = 0, 2(100 + 𝑡) − 3
𝑐(𝑡) = 100+𝑡
= 0, 2 − 4
(100+𝑡) (100+𝑡)

7 7
2𝑥10 2𝑥10
0, 1 = 0, 2 − 4
⇒ 0, 1 = 4
⇒
(100+𝑡) (100+𝑡)

4 8
⇒ 𝑡 = 100( 2 − 1)
4
(100 + 𝑡) = 2𝑥10

3.5.3 Ejercicio de aplicación

Se sabe que un cierto material radiactivo se desintegra proporcionalmente a la cantidad

presente. Si inicialmente hay 50 miligramos de material presente y después de dos horas se
observa que el material ha perdido el 10% de su masa original, hallar:

a) una expresión para la masa de material presente en un momento t

b) la masa después de cuatro horas

c) el tiempo para el cual el material se ha desintegrado en la mitad de su masa inicial.

SOLUCIÓN

a) 𝑥(𝑡)=miligramos de material en el instante t

𝑑𝑥
𝑑𝑡
=− 𝑘𝑥(𝑡)
𝑥(0) = 50

𝑑𝑥 −𝑘𝑡
𝑥(𝑡)
=-kdt⇒ 𝑙𝑛(𝑥(𝑡)) =− 𝑘𝑡 + 𝑐 ⇒ 𝑥(𝑡) = 𝑐𝑒

𝑥(0) = 50 50 = 𝑐

−2𝑘 45
𝑥(2) = 45 45 = 50𝑒 ⇒ − 2𝑘 = 𝑙𝑛|| 50 || ⇒ 𝑘 = 0, 0526803

−0,0526803𝑡
𝑥(𝑡) = 50𝑒

−0,0526803𝑡
b) 𝑥(4) = 50𝑒 = 40, 5miligramos

−0,0526803𝑡 −0,0526803𝑡
c) 𝑥(𝑡) = 25 ⇒ 25 = 50𝑒 ⇒𝑒 = 0, 5 ⇒ 𝑡 = 13, 1576 horas

3.6 Movimiento Armónico de un Resorte Vibrante

Un resorte ordinario de peso despreciable está siendo suspendido de forma vertical de un

soporte fijo. Se coloca un peso W en el resorte y que este cuelgue. Su punto de reposo se
denomina como Posición de Equilibrio, si se jala el peso y se suelta se producirá un
movimiento vibratorio alrededor. Se quiere analizar el movimiento del peso, siendo necesario
saber las fuerza que influencian en el peso durante ese movimiento. Sabiendo que hay una
fuerza que restaura el peso desplazado a la posición de equilibrio, llamado la fuerza
restauradora. Esta se rige con la Ley de Hooke: “La fuerza ejercida que restaura el peso W a
la Posición es proporcional a la distancia de W al Punto de Equilibrio”.
Denominando la fuerza restauradora como f, x siendo la posición del peso W desde la
posición de equilibrio. x siendo positivo cuando W está por debajo de la posición y negativo
cuando está encima de esta, formando:

|𝑓|∞|𝑥| → |𝑓| = 𝑘|𝑥| k > 0

siendo una constante de proporcionalidad según la dureza del resorte, la constante del resorte.
La dirección de la fuerza depende si x > 0 entonces va hacia arriba y el valor es negativo y
viceversa si es x < 0, satisfaciendose si la fuerza está en magnitud y dirección por -kx, siendo
la fórmula:

𝑓 =− 𝑘𝑥 (1)

Al colocarse el peso W este se estirara una s distancia. Y en la Ley de Hooke, la tensión T es

proporcional al estiramiento

𝑇1 = 𝑘𝑠 (2)

Dado a que el peso y el resorte están en equilibrio da paso a que

𝑇1 = 𝑘𝑠 = 𝑊 (3)

Mientras que el peso se hala y suelta, la posición cambia en cada punto del tiempo, siendo la
Segunda Tensión este tiempo según la ley de Hooke:

𝑇2 = 𝐾(𝑠 + 𝑥) (4)

Sigue que la fuerza neta en dirección positiva está dada por:

𝑊 · 𝑇 = 𝑊 · 𝑘𝑠 · 𝑘𝑥 =− 𝑘𝑥 (5)

Y en base a la Ley de Newton y con la ecuación 3, la ecuación de movimiento es:

2
𝑤 𝑑𝑥
𝑔
· 2 =− 𝑘𝑥 La Fuerza Neta es la Restauradora y no depende del Peso W
𝑑𝑡

3.6.1. Ejercicio de aplicación :

Se encontró experimentalmente que un peso de 6 Ib estira un resorte 6 pulg. Si el peso se hala
4 pulg por debajo de la posición de equilibrio y se suelta: (a) Establezca una ecuación
diferencial y condiciones asociadas que describen el movimiento; (b) Encuentre la posición
del peso como una función del tiempo y (c) determine la posición, velocidad y aceleración del
peso 1/2 seg después de haberse soltado.

6 𝑝𝑢𝑙𝑔 = 1/2 𝑝𝑖𝑒 𝑓 = 𝑘𝑥 → 6 = 𝑘 · 1/2 𝑘 = 12

2 2
6 𝑑𝑥 𝑑𝑥
32
· 2 =− 12𝑥 ⇒ 2 + 64𝑥 = 0 (5)
𝑑𝑡 𝑑𝑡

En inicio es t=0 siendo el peso 4 pulg debajo de la posición: x = ⅓ pie en t=0 (6)

𝑑𝑥
y al soltarse el peso se forma: 𝑑𝑡
= 0 (7)

2
La Ecuación auxiliar para el 5 es 𝑚 + 64 = 0 con raíces 𝑚1/2 =± 8𝑖 dando como

solución:

1
𝑥 = 𝐴𝑐𝑜𝑠8𝑡 + 𝐵𝑠𝑒𝑛8𝑡 siendo A = ⅓ entonces: 𝑥 = 3
𝐴𝑐𝑜𝑠8𝑡 + 𝐵𝑠𝑒𝑛8𝑡

𝑑𝑥 8 1
𝑑𝑡
=− 3
𝑠𝑒𝑛8𝑡 + 8𝐵𝑐𝑜𝑠8𝑡 luego se halla B = 0, entonces: 𝑥 = 3
𝑐𝑜𝑠𝑟 (8)

x está en pies, teniendo que transformarse en pulg dándose la ecuación 𝑥 = 4𝑐𝑜𝑠8𝑡

2
𝑑𝑥 8 𝑑𝑥 64
Derivando 8: 𝑟 = 𝑑𝑡
=− 3
𝑠𝑒𝑛8𝑡 𝑎= 2 =− 3
𝑐𝑜𝑠8𝑡
𝑑𝑡

Tras pasar en t = ½ seg y que 4 𝑟𝑎𝑑 = 4(180/π) 𝑔𝑟𝑎𝑑𝑜𝑠 = 229 𝑔𝑟𝑎𝑑𝑜𝑠 Encontrando:

1 8
𝑋= 3
(− 0. 656) =− 0. 219 υ =− 3
(0. 755) =+ 2. 01
64
𝑎 =− 3
(− 0. 656) =+ 14

Tras ½ seg el peso está a 0.219 pies encima de la Posición de Equilibrio con una velocidad de
2.01 pies/seg y aceleración de 14 pies/seg2

Con el gráfico se ve que empieza con x = 1⁄3 donde t = 0, llendo de la posición de equilibrio
con x = 0 a x = -⅓, este regreso de nuevo al punto y vuelve al inicio de x = ⅓, siendo ciclico
viendose el proceso periodico en la imagen. siendo este el Movimiento Armónico Simple,
siendo descrito en la fórmula:

2
𝑑𝑥
2 =− 𝑎𝑥 a>0 (9)
𝑑𝑡

es la constante del movimiento armónico, siendo esta directamente proporcional al

desplazamiento pero en distinta dirección. El desplazamiento máximo del peso se denomina
amplitud, siendo la amplitud de ⅓ pie, siendo un periodo todo el ciclo que recorre el peso a la
amplitud, el periodo es de π/4 seg. Otro método para identificar el periodo es con la
trayectoria y si usamos el punto más bajo ⅓, haciendo que esto ocurrirá cuando:

𝑐𝑜𝑠8𝑡 = 1 8𝑡 = 0, 2π, 4π, 6π ó 𝑡 = 0, π/4, 2π/4, 3π/4

𝑃𝑒𝑟𝑖𝑜𝑑𝑜 (𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑠𝑒𝑔𝑢𝑛𝑑𝑜𝑠 𝑝𝑜𝑟 𝑐𝑖𝑐𝑙𝑜) = π/4

𝐹𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 (𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑐𝑖𝑐𝑙𝑜𝑠 𝑝𝑜𝑟 𝑠𝑒𝑔𝑢𝑛𝑑𝑜) = 4/π

1
Por lo tanto T es el periodo, entonces la frecuencia f está dada por: 𝑓 = 𝑇

3.7. Flujo de calor en estado estacionario

Tenemos un flujo de calor en estado estacionario cuando la temperatura en cada punto de un

cuerpo no varía con respecto al [Link] ejemplo de esto puede representarse en la siguiente
figura(ver Figura 1), donde asumimos que el material es uniforme en todas sus propiedades y
tiene una longitud desconocida delimitada por planos paralelos, los cuales denominan A y B,
por un lado el plano A presenta una temperatura de 50°C en este ejemplo, mientras que el
plano B presenta una temperatura de 100°C, en ambos casos la temperatura se mantiene
constante, por lo que no varía. De igual manera sucede en el plano medio C y en su sucesivo
plano E, que presenta una temperatura constante de 75°C y de 90°C respectivamente.

(Imagen extraída de : Murray Spiegel - Ecuaciones Diferenciales Aplicadas)

Otro claro ejemplo es el de un tubo de material uniforme, donde se supone que la parte
interior y la exterior se mantienen a temperatura constante, de 40°C y 80°C respectivamente.
Además habrá una superficie(línea punteada), que tendrá 60°C , pero no será el punto medio
entre las superficies internas y externas.

(Imagen extraída de : Murray Spiegel - Ecuaciones Diferenciales Aplicadas)

Debemos de tener en cuenta que las líneas paralelas al plano A y en un plano perpendicular a
A se llaman líneas isotérmicas. La curva isotérmica se representa como la curva punteada .
Los planos de las figuras anteriores se llaman superficies isotérmicas.

De manera general, las curvas isotérmicas no se representan como líneas o círculos, como en
las figuras anteriores, pero pueden ser una familia de curvas punteadas y las trayectorias
ortogonales se denominan líneas de flujo. A continuación una figura con estos elementos.
 (Imagen extraída de : Murray Spiegel - Ecuaciones Diferenciales Aplicadas)

Para entender este fenómeno debemos de saber que si tomamos porciones pequeñas de 2
superficies isometricas, las cuales están separadas por una distancia ∆𝑛 y asumimos que la
temperatura de la superficie 𝑆1 y de 𝑆2 es correspondiente a 𝑈1 y 𝑈2 respectivamente. Por lo
que la diferencia de temperaturas ∆𝑈 sería igual a 𝑈2 − 𝑈1 y de manera experimental se
determinó que la cantidad de calor que fluye de 𝑆1 a 𝑆2 es aproximadamente proporcional a
∆𝑈
∆𝑛
, esta aproximación es más precisa, cuando los valores de ∆𝑈 y∆𝑛 se hacen más
pequeños.

(Imagen extraída de : Murray Spiegel - Ecuaciones Diferenciales Aplicadas)

𝑑𝑈
Gracias a esto podemos determinar la gradiente de U, la cual se representa como 𝑑𝑛
o
también conocida como la tasa de cambio de U en la dirección normal a la superficie o curva
isoté[Link] ser H la cantidad de flujo de calor por unidad de área y unidad de tiempo,
utilizamos la siguiente ley física:
 𝑑𝑈
𝐻= 𝐾 𝑑𝑛

𝑑𝑈
Tener en cuenta que si “U” depende de otros factores además de “n” se reemplazara 𝑑𝑛
por
δ𝑈
δ𝑛
. También si consideramos a H como un vector, entonces la dirección del flujo de calor
puede tenerse en cuenta por un signo menos en la ecuación anterior.
Finalmente la cantidad de calor por unidad de tiempo que fluye a través de un área A
está dada por la siguiente ecuación:

𝑑𝑈
𝑞 =− 𝐾𝐴 𝑑𝑛

Donde la constante K, depende del material usado y se conoce como conductividad térmica.
Además cabe mencionar que la cantidad de calor deberá expresarse en calorías.

3.7.1. Ejercicio aplicativo

“Un tubo largo de acero, de conductividad térmica k = 0,15 unidades cgs, tiene un radio
interior de 10 cm y un radio exterior de 20 cm. La superficie interna se mantiene a 20”
C y la superficie exterior se mantiene a 50” C. (a) Encuentre la temperatura como una
función de la distancia r del eje común de los cilindros concéntricos. (b) Encuentre la
temperatura cuando r = 15 cm. (c)¿Cuánto calor se pierde por minuto en una parte del
tubo de 20 m de largo?” (Murray S, 1992)

Solución

Si r = 15, encontramos por sustitución que U = 114°C. Del valor anterior de 9, la

respuesta para la parte c) es:
q = 408.000 × 60 𝑐𝑎𝑙/𝑚𝑖𝑛 = 24. 480. 000 𝑐𝑎𝑙/𝑚𝑖𝑛
3.8. La resonancia mecánica

Cuando la frecuencia de una fuerza periódica aplicada a un sistema mecánico se relaciona

de forma simple con la frecuencia natural del sistema, puede ocurrir la resonancia mecánica.
Este fenómeno puede generar oscilaciones amplificadas que pueden llevar al colapso del
sistema, como en el caso de un puente que se derrumba debido al paso de una compañía de
soldados marchando en fila. Del mismo modo, una nota musical con una frecuencia
específica puede hacer estallar un cristal. Dado el potencial dañino de la resonancia
mecánica, es importante que los ingenieros eviten este fenómeno al diseñar estructuras o
mecanismos vibrantes.

Sin embargo, la resonancia mecánica también puede ser utilizada de manera útil. Por
ejemplo, cuando un automóvil queda atascado en la nieve o el barro, es posible liberarlo al
hacerlo vibrar con su frecuencia natural. Luego, al aplicar una fuerza con esa misma
frecuencia, la resonancia mecánica resultante puede ser suficiente para liberar el automóvil.

3.8.1. Ejemplo aplicativo:

Un peso 𝑃 estira un resorte 𝑥1 unidades de longitud. Si el peso se halla 𝑥2 unidades de

longitud por debajo de la posición de equilibrio y se suelta aplicándole una fuerza externa
dada por 𝐹(𝑡) = 𝐴 𝑐𝑜𝑠 𝐵𝑡 .

Describir el movimiento que resulta si se supone que inicialmente el peso está en la posición
de equilibrio (𝑥 = 0) y que su velocidad inicial es cero.

Formulación matemática

La ecuación diferencial que describe el movimiento es por tanto:

2
𝑃 𝑑𝑥
𝑔 𝑑𝑡2
=− 𝑘𝑥 + 𝐹(𝑡)

Puesto que inicialmente (𝑡 = 0) el peso está 𝑥2 por debajo de la posición de equilibrio,

tenemos 𝑥 = 𝑥2 en 𝑡 = 0.
𝑑𝑥
También, puesto que el peso se suelta (esto es, tiene velocidad cero) en 𝑡 = 0, 𝑑𝑡
= 0 en
𝑡 = 0.

Solución: Como 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. = cos)(
BtAtF Ensayando como solución particular

4. Conclusiones:
 -Se logró utilizar los diversos métodos de ecuaciones diferenciales para modelar y
resolver distintos problemas que surgen al momento de poner en práctica los
procesos metalúrgicos.

-Con la aplicación de las fórmulas dadas en el presente trabajo podemos llegar a

calcular diversos ejercicios , como es en el caso de la masa del soluto en un tanque
de mostrado .

-Desarrollar y aplicar los métodos son de gran ayuda para ver como es el
funcionamiento de cada problemática y también nos ayuda a fomentar un
pensamiento crítico.

-Los resultados obtenidos en los distintos ejercicios aplicativos fueron positivos ya

que como podemos visualizar se ha podido solucionar el problema principal.

[Link]ía:

● ¿Qué es la tasa de flujo volumétrico? (s/f). Khan Academy. Recuperado el 9 de julio 2023,
[Link]
e
● Universidad Autónoma del Estado de Hidalgo. (s/f). Vista de Aplicación de ecuaciones
diferenciales lineales de primer orden a mezclas. [Link]. Recuperado el 9 de julio de 2023,
de [Link]
● (S/[Link][Link]
rabajo%20Ec%20Diferenciales%20en%[Link]