IMPORTANCIA DE LAS ECUACIONES

DIFERENCIALES EN INGENIERIA MECANICA:

INTRODUCCIÓN

En este documento se hablara sobre la importancia de las ecuaciones

diferenciales en la ingeniería mecánica, partiendo de lo general a lo particular, es
decir, desde el concepto matemático de una ecuación diferencial, hasta algunas
de sus aplicaciones en el área. Sin embargo, se debe destacar que en particular
para la mecánica, las ecuaciones diferenciales son fundamentales, sencillamente
porque los principios de la mecánica se apoyan en ecuaciones diferenciales.

La ingeniería mecánica es una rama de la ingeniería, que aplica las ciencias

exactas, específicamente los principios físicos termodinámica, mecánica, ciencia
de materiales, mecánica de fluidos y análisis estructural, entre otros, para el
diseño y análisis de diversos elementos usados en la actualidad, además de
solucionar diferentes problemas de las mismas áreas. Apoyándose en la definición
de ingeniería mecánica o simplemente al atender al nombre, uno puede ver que la
mecánica es la ciencia elemental de esta ingeniería, por eso vale la pena dedicar
un poco de espacio al análisis de la mecánica.

La mecánica, es la ciencia física que estudia el movimiento de los cuerpos, en la

actualidad se cuenta con diferentes enfoques, mecánica clásica, mecánica
relativista, mecánica frecuentista. Afortunadamente, para estudiar el movimiento
de los objetos que encontramos en nuestra vida diaria, objetos que ni alcanzan
velocidades cercanas a la de la luz ni objetos con dimensiones atómicas, no
necesitamos mecánica cuántica o relativista. La mecánica tiene clásica tiene como
base 4 principios que son: El principio de Stevin, y las 3 leyes de Newton, cosa
que debemos tener en mente para esta lectura.

1
La mecánica clásica, se divide en mecánica de fluidos, de cuerpos deformables y
de cuerpos rígidos, es evidente que en cada una de estas divisiones podemos
encontrar aplicación de las ecuaciones diferenciales, pero me gustaría que nos
concentráramos en la mecánica de cuerpos rígidos y en sus sub divisiones, que
son la estática y la dinámica.

La estática es la encargada del estudio de los sistemas de fuerzas en equilibrio

mientras la dinámica, describe el movimiento de los cuerpos y los relaciona con
sus causas, y se apoyan principalmente en las leyes de newton, para la dinámica
la segunda ley de newton toma el papel protagónico. Newton escribió:

LEX II. Mutationem motus proportionalem esse vi motrici impressæ, et fieri

secundum
lineam rectam qua vis illa imprimitur.

Cuya traducción es como sigue: “El cambio del movimiento es proporcional a la

fuerza motriz impresa, y ocurre según la línea recta en la cual se imprime esa
fuerza. “
En esta ley podemos leer ya la presencia de la ecuación diferencial, al leer en
cambio en la cantidad de movimiento, sabiendo que una ecuación diferencial es
aquella que incluye derivadas, y la aplicación física de la derivada la podemos
decir cómo, la tasa de cambio de un elemento A con respecto a un elemento B. Se
tiene entonces el cambio de la cantidad de movimiento, nuestro elemento A y al
tiempo como nuestro elemento B. A continuación se presenta una explicación
detallada.

2
La cantidad de movimiento ( momentum ) de un objeto se define como su masa
“m” multiplicada por su velocidad “v”. La tasa de cambio en momentum en el
tiempo es así:

Si denotamos por “F” la fuerza neta que actúa sobre el cuerpo la segunda ley dice
que:

Donde el símbolo denota proporcionalidad. Introduciendo la constante de

proporcionalidad obtenemos

Si “m” es una constante y a = dv/dt’ es la aceleración.

Así vemos que

El valor de k depende de las unidades que deseemos usar.

En el sistema CGS la longitud se mide en centímetros (cm), la masa en gramos
(g), y el tiempo en segundos (s). El valor más simple para k es k = 1, de modo
que :

3
Con esta ley se pueden resolver un sin número de problemas donde se
relacionan, velocidades, distancias, aceleraciones donde cada una de ellas se
puede expresar como tasas de cambio.

De lo anterior podemos afirmar la importancia de las ecuaciones diferenciales, al

encontrarlas en uno de los fundamentos de la mecánica, eso solo por mencionar
un caso, aunque bien podríamos hablar del enfoque clásico de la termodinámica
donde encontramos la ecuación fundamental de la hidrostática, los balances de
masa y energía, entre otros y así seguir con cada una de las ciencias para
descubrir que están íntimamente ligadas.

¿QUÉ ES UNA ECUACIÓN DIFERENCIAL?

Una ecuación diferencial es aquella ecuación que contiene derivadas o

diferenciales. Las derivadas pueden ser ordinarias, en cuyo caso a ecuación
diferencial será ordinaria; o bien derivadas parciales y la ecuación diferencial será
en derivadas parciales.

Las ecuaciones diferenciales aparecen en modelos matemáticos que tratan de

describir situaciones de la vida real. Así, muchas leyes naturales pueden ser
traducidas al lenguaje matemático mediante ecuaciones que envuelven derivadas,
como en física, donde la velocidad y la aceleración aparecen como derivadas

En diversos modelos matemáticos, para obtener una ecuación diferencial que

describa un problema real, se asume que la situación se rige por leyes simples.
Una vez que se construye el modelo en forma de ecuación diferencial, lo que
viene es solucionarla y con estas soluciones, se hacen predicciones relativas al
comportamiento del problema en cuestión.

4
APLICACIONES EN DIFERENTES ÁREAS DE LA INGENIERÍA
MECÁNICA

EL MOVIMIENTO ARMÓNICO SIMPLE

Supongamos que tenemos un resorte dispuesto en forma vertical, con el

extremo superior fijado a una superficie horizontal y el otro extremo libre al
cual se lefia una carga de masa m, la fuerza de recuperación del resorte esta
dada por la Ley de Hooke:

F = ks

donde
S: elongación del resorte.
K: constante elástica del resorte.

Deduzcamos la E.D. que rige el movimiento de la masa sujeta al resorte

Por la segunda Ley de Newton tenemos

5
Llamemos

ED segundo orden coeficientes ctes.

Ecuación Característica

Condiciones Iniciales si :

Si α > 0: el cuerpo esta por debajo de la P.E. (Posición de Equilibrio).

Si α < 0: el cuerpo esta por encima de la P.E.

Si α = 0: el cuerpo esta en la P.E.

Si β> 0 el cuerpo inicia su movimiento con velocidad hacia abajo.

Si β< 0 el cuerpo inicia su movimiento con velocidad hacia arriba.

Si β= 0 el cuerpo inicia su movimiento desde el reposo.

6
Periodo de Vibraciones Libres T se define : T = 2¼

Frecuencia f de vibraciones libres se define como el inverso del periodo,

es decir:

F=1/T

MOVIMIENTO AMORTIGUADO

Cuando el cuerpo sujeto al resorte se mueve en un medio que produce fricción sobre el
cuerpo, entonces decimos que el movimiento se afecta con amortiguamiento,
supongamos que el amortiguamiento es directamente proporcional a la velocidad. Por la
segunda ley de Newton tenemos:
La constante A se le llama la amplitud del Movimiento Armónico Simple

(M.A.S.). y el ángulo Φ se le llama ángulo de Fase.

Como

Donde β constante de fricción o

7
constante de amortiguamiento y dividiendo entre m.

Con carga en movimiento

Donde

Ecuación característica:

Cuando > 0, entonces a este movimiento se le llama sobre amortiguado en

este caso p1 y p2 son negativos, y como

por lo tanto

Cuando = 0, a este movimiento se le llama críticamente amortiguado

8
por lo tanto

Cuando < 0 o lo que es lo mismo > 0, a este movimiento se le llama

subamortiguado

9
Donde Φ se le llama el ángulo de fase.

Nota: respecto al movimiento sobre amortiguado y críticamente amortiguado debemos

anotar que no son movimientos oscilatorios y cruzan a lo sumo una vez la posición de
equilibrio como se ve en las primeras dos graficas.

Péndulo Simple

En la naturaleza hay muchos movimientos que se repiten a intervalos iguales de tiempo,

estos son llamados movimientos periódicos. En Física se ha idealizado un tipo de
movimiento oscilatorio, en el que se considera que sobre el sistema no existe la acción
de las fuerzas de rozamiento, es decir, no existe disipación de energía y el movimiento

10
se mantiene invariable, sin necesidad de comunicarle energía exterior a este. Este
movimiento se llama MOVIMIENTO ARMÖNICO SIMPLE (MAS), este es un
movimiento vibratorio bajo la acción de una fuerza recuperadora elástica, proporcional al
desplazamiento y en ausencia de todo rozamiento.

Un ejemplo de este movimiento se puede encontrar a partir del desplazamiento de un

punto cualquiera alrededor de toda la longitud de una circunferencia.

Cuando un punto (P) recorre una circunferencia con velocidad uniforme, su proyección
(Q) sobre cualquiera de los diámetros de esta, realiza un tipo de movimiento armónico
simple. Cada vez que el punto se encuentre en uno de los cuatro cuadrantes de la
circunferencia, se trazará una perpendicular desde el punto a un diámetro fijo de la
circunferencia. A medida que el punto escogido se mueve a velocidad uniforme, el punto
proyectado en el diámetro, realizará un movimiento oscilatorio rectilíneo.

Para representar gráficamente (en una función) el movimiento armónico simple de un

punto, se toman como abscisas los tiempos medidos como fracciones del período (T/12,
T/6, T/4...) que es el tiempo que este punto tarda en dar una vuelta completa a la
circunferencia; y como a ordenadas las sucesivas prolongaciones del mismo. La
resultante es una sinusoide, ya que la variación del tiempo t, se traduce como una
variación del sin x, donde x es el ángulo que forma el radio con el semi-eje positivo de
abscisas (x es proporcional al tiempo).

11
Elementos:

1. Oscilación o vibración: es el movimiento realizado desde cualquier posición hasta

regresar de nuevo a ella pasando por las posiciones intermedias.

2. Elongación: es el desplazamiento de la partícula que oscila desde la posición de

equilibrio hasta cualquier posición en un instante dado.

3. Amplitud: es la máxima elongación, es decir, el desplazamiento máximo a partir de la

posición de equilibrio.

4. Periodo: es el tiempo requerido para realizar una oscilación o vibración completa. Se

designa con la letra "t".

5. Frecuencia: es el número de oscilación o vibración realizadas en la unidad de tiempo.

6. Posición de equilibrio: es la posición en la cual no actúa ninguna fuerza neta sobre la

partícula oscilante.

Relación entre el M.A.S. y el Movimiento Circular Uniforme

El M.A.S. de un cuerpo real se puede considerar como el movimiento de la "proyección"

(sombra que proyecta) de un cuerpo auxiliar que describiese un movimiento circular
uniforme (M.C.U.) de radio igual a la amplitud A y velocidad angular ω, sobre el díametro
vertical de la circunferencia que recorre.

En lo siguiente podrás visualizar dicha relación.

Vamos a establecer una relación entre un movimiento vibratorio armónico simple y el

12
movimiento circular uniforme. Esto nos va a permitir dos cosas:

- Hallar la ecuación del (MAS )sin tener que recurrir a cálculos matemáticos complejos.

- Conocer de donde vienen algunos de los conceptos que usamos en el (MAS), como
frecuencia angular o el desfase.

Observando tememos inicialmente el resorte azul, que oscila verticalmente. En la

circunferencia tienes un punto negro que gira con movimiento circular uniforme,
ocupando en cada instante una posición en la circunferencia. Traza mentalmente la
proyección de esa posición sobre el diámetro vertical de la circunferencia. En cada
momento, la masa que cuelga del resorte ocupa una posición determinada. Observa
que la posición de la masa del resorte coincide exactamente con la proyección de la
posición del objeto sobre el diámetro, que verás en forma de línea azul en el diámetro
vertical.

Es decir, como resumen, cuando un objeto gira con movimiento circular uniforme en una
trayectoria circular, el movimiento de la proyección del objeto sobre el diámetro es un
movimiento armónico simple.

Lo mismo podríamos decir del resorte amarillo y la proyección sobre el diámetro

horizontal, que verás como un trazo amarillo sobre dicho diámetro.

Los vectores azul y amarillo, que varían en el applet, corresponden al valor de la

velocidad del resorte, azul para diámetro vertical y amarillo para el horizontal. Observa
su variación y comprobarás que la velocidad es máxima en el centro de equilibrio del
resorte y mínima en los extremos, en los puntos de mínima y máxima elongación.

13
Observa también como el vector rojo de la gráfica de la derecha, la velocidad del (MAS),
coincide con el vector azul, la velocidad de la proyección sobre el diámetro vertical, lo
que supone una prueba más de lo que hemos afirmado anteriormente.

Ecuaciones del Movimiento Armónico Simple

x = A . cos . w . t

x = elongación

r = A = radio

t = tiempo

w = velocidad angular

Vx = - V . sen Ø

V=w.r

h=w.t

w . t = V = Vector representativo de la velocidad lineal.

Vx = proyección de "Y" sobre el eje "X"

h = ángulo

Vx = -2 . F . A . sen (2 .)

Vx = + w " A2 - x2

Ax = - w2 . A . cos. w . t

Ax = - Ac . cos Ø

Ac = proyección de aceleración sobre el eje horizontal

Ac = w2 . x

Ac = aceleración centrípeta

t=2 " mk

T = periodo

14
Péndulo simple

Es llamado así porque consta de un cuerpo de masa m, suspendido de un hilo largo de

longitud l, que cumple las condiciones siguientes:

 el hilo es inextensible
 su masa es despreciable comparada con la masa del cuerpo
 el ángulo de desplazamiento que llamaremos 0 debe ser pequeño

Vemos pues que, considerando únicamente el desplazamiento tangente a la trayectoria,

es decir, el arco que se está recorriendo, podemos poner

Donde no hemos hecho si no aplicar la ley de Newton. Esto se puede ver considerando
que el arco es y , como l es la longitud del hilo y es constante , la aceleración será

Por otra parte aplicando , en este caso la fuerza es solo la de la

gravedad , mg que se descompone en la trayectoria marcada por el arco ;

15
Esta ecuación no es nada fácil de resolver , es por eso que recurrimos a decir

quedándonos entonces

Que a veces también se expresa como .

Esta ecuación es absolutamente análoga a la de un movimiento armónico simple, y por

tanto su solución también será la misma teniendo, únicamente, la precaución de sustituir
el valor de antiguo por el que tiene ahora para un péndulo

A partir de aquí se pueden extraer todas las demás relaciones para un péndulo simple,
el periodo, frecuencia, etc.

Período de un Péndulo

Período: Se define como el tiempo que se demora en realizar una oscilación completa.
Para determinar el período se utiliza la siguiente expresión T/ N° de Osc. ( Tiempo
empleado dividido por el número de oscilaciones).

1) El periodo de un péndulo es independiente de su amplitud. Esto significa que si se

tienen 2 péndulos iguales (longitud y masa), pero uno de ellos tiene una amplitud de
recorrido mayor que el otro, en ambas condiciones la medida del periodo de estos
péndulos es el mismo.

2) El periodo de un péndulo es directamente proporcional a la raíz cuadrada de su

longitud. Esto significa que el periodo de un péndulo puede aumentar o disminuir de
acuerdo a la raíz cuadrada de la longitud de ese péndulo.

Un ejemplo sencillo y aplicable en ingeniería mecánica es una grúa de demolición:

16
CIRCUITOS ELECTRICOS

Las ecuaciones diferenciales también se pueden aplicar en electricidad , como a los

transitorios de circuitos :

Régimen transitorio en circuitos RCL

Régimen transitorio en corriente continua

En el siguiente circuito, al cerrar la llave L se producirá un fenómeno transitorio que

hemos de estudiar.

En el instante t = 0 en que se cierra la llave L, la intensidad i , variable y función del

tiempo, será cero ya que la inductancia en ese instante ha de actuar como una llave
abierta o una resistencia de valor infinito, cayendo toda la tensión de la fuente en ella.

17
Cuando el tiempo tienda a infinito la intensidad i, variable y función del tiempo, también
será cero ya que el capacitor actuará como una llave abierta, cayendo en él toda la
tensión de la fuente.

En medio tenemos el régimen transitorio, y en él, las caidas de tensión en cada

elemento serán:

Si aplicamos la ley de las mallas de Kirchoff al circuito RLC de la figura, obtendremos la

siguiente expresión:

si derivamos esta última expresión:

reordenando:

y esta última es una ecuación lineal de 2º orden con coeficientes constantes y

homogénea, la misma la resolvemos ya hemos explicado en el punto 3, es decir:

reemplazando:

18
y si hacemos:

entonces: y
Se nos presentarán los tres casos que estudiamos y eso dependerá del valor del
discriminante:

Caso I] B> 0 entonces r1 y r2 raíces reales y distintas.

Solución:

Para determinar las constantes aplicamos las condiciones iniciales:

a. en t = 0 el valor de i = 0
b. en t = 0 la tensión de la fuente cae en la inductancia

de a: es decir (a)

de b:

(b)
(a) y (b) forman un sistema de ecuaciones que nos permitirá determinar el valor de las
constantes, por ejemplo en (b) reemplazamos C1 por –C2 y hallamos C2 y luego con (a)
hallamos C1.

19
Caso II] B = 0 entonces r1 y r2 raíces reales e iguales.

Solución:
Para determinar las constantes aplicamos las condiciones iniciales de la misma forma
que en el caso I.

Caso III] B < 0 entonces r1 y r2 raíces complejas conjugadas.

Solución:
Para determinar las constantes aplicamos las condiciones iniciales de la misma forma
que en el caso I.

Resolución de regímenes transitorios para el circuito RLC del punto anterior

Basaremos, por supuesto, la resolución en todo lo desarrollado en el punto anterior.

Caso I] Sea el siguiente circuito:

20
 B> 0 entonces r1 y r2 raíces reales y distintas

(1)
Condiciones iniciales:

para t = 0 la i = 0 y V = vL

reemplazando en (1): es decir

reemplazando

y
con lo que:

La representación gráfica que nos muestra el desarrollo en el tiempo de este efecto o

régimen transitorio es el siguiente:

Caso II] Sea el siguiente circuito:

21
 B = 0 entonces r1 y r2 raíces reales e iguales

(2)
Condiciones iniciales:

para t = 0 la i = 0 y V = vL

reemplazando en (2): .0 es decir

y
con lo que:

22
 La representación gráfica muestra el desarrollo en el tiempo de este efecto

Caso III] Sea el siguiente circuito:

B < 0 entonces r1 y r2 raíces complejas conjugadas

(3)
Condiciones iniciales:

para t = 0 la i = 0 y V = vL

23
reemplazando en (3): es decir

con lo que:

La representación gráfica que nos muestra el desarrollo en el tiempo de este efecto o

régimen transitorio es el siguiente:

Régimen transitorio en corriente alterna.

24
Al cerrar la llave L la fuente aplica una tensión v variable en el tiempo de forma
sinusoidal.

Aplicando la ley de las mallas de Kirchoff tenemos:

derivando:

reordenado la expresión anterior y dividiendo por L miembro a miembro:

(1)
ecuación diferencial de 2º orden, con coeficiente constantes y no homogénea.

Para resolver esta ecuación diferencia y estudiar el régimen transitorio respectivo,

emplearemos uno de los métodos de resolución explicados en los apartados anteriores,
por ejemplo, el método de los coeficientes indeterminados.

Determinación de la función complementaria yh :

Resolvemos
Esto ya lo sabemos hacer, de manera que obtendremos una de las siguientes
soluciones:

25
Determinación de la solución particular yp :

Partes variables:

por lo tanto:

reemplazando en (1):

en esta última expresión agrupamos los términos demejantes:

(a)

(b)
reordenando las dos últimas expresiones:

resolviendo este sistema hallaremos los coeficientes A y B.

Si planteamos el determinante del coeficiente de las incógnitas y lo resolvemos,

obtendremos:

26
operando:
por otro lado desarrollando el determinante de la incógnita A:

dividiendo numerador y denominador por :

pero: y

reemplazando en la expresión de A:

por otro lado desarrollando el determinante de la incógnita B:

con lo cual:

27
dividiendo ambos miembros por y operando nos encontraremos con un
denominador idéntico al de la resolución de la incógnita A, por tanto:

y como :

pero sabemos que donde Z es la impedancia del circuito y que

donde x es la reactancia resultante del circuito, por lo tanto:

donde y donde  es el ángulo de fase de la impedancia de carga.

28
y por trigonometría:

Las soluciones completas podrán ser, según la solución de la homogénea:

Para que se entienda bien lo anterior mostraremos un :

Ejemplo.

29
Resolvemos la homogénea:

y obtenemos y

Hallamos la solución particular:

donde

donde

La intensidad completa será:

luego, para t = 0 la i = 0 y , por lo tanto

aplicamos a la solución general la primera condición para t = 0 la i = 0

por otro lado:

30
resolviendo este sistema obtenemos: y

Regímenes transitorios en circuitos de dos mallas.

Estudiaremos el siguiente circuito de dos mallas y en el instante en que para t = = se

cierra la llave L y veremos que ocurre con las intensidades en cada una de las ramas del
mismo.

31
Para t = 0 las corientes i1 e i2 tienen valor cero ya que las inductancia actúan como una
llave abierta; luego de pasado mucho tiempo, la inductancia de la rama central actuará
como un cortocircuito y la intensidad i1 = 3 A mientras que la otra corriente será cero.

Pero veamos que ocurre en la transición, planteando las ecuaciones de malla:

Si a la expresión la llamamos o representamos por nuestro sistema quedará:

si a ambas ecuaciones las dividimos m.a.m por el valor de las inductancias, con el fin de
dejar a sin un factor que lo multiplique:

En este sistema de ecuaciones lo podemos representar:

. =
donde para calcular cada una de las corrientes aplicamos la Regla de Cramer:

. i1 = desarrollando:

pero por lo tanto y volviendo a la notación normal:

32
donde la solución de la homogénea nos da valores y por tanto:

y la solución particular, vista al principio de este punto es ; luego la solución

general será:

donde las constantes las determinamos a partir de las condicione iniciales:

para t = 0 la corriente i1 = 0 y toda la tensión de la fuente cae en la inductancia:

, aplicando esta condiciones a la expresión anterior de la intensidad 1,

determinamos las constantes y obtenemos:

. i2 = desarrollando:

pero por lo tanto y volviendo a la notación normal:

donde la solución de la homogénea será la misma que para la corriente anterior por
tener la misma ecuación característica:

y la solución particular, vista al principio de este punto es ; luego la solución

general será:

donde las constantes las determinamos a partir de las condicione iniciales:

para t = 0 la corriente i2 = 0 y toda la tensión de la fuente cae en la inductancia:

33
 , aplicando esta condiciones a la expresión anterior de la intensidad 2,
determinamos las constantes y obtenemos:

y la intensidad en la inductancia central será:

LEY DE ENFRIAMIENTO DE NEWTON

Debido a que en ingeniería Mecánica trabajamos con diversas sustancias y
objetos (placas metálicas, motores, turbinas etc.) a temperaturas variables. Las
ecuaciones diferenciales nos ayudan a modelar el comportamiento y la variación
de temperatura de estos cuerpos y sustancias, estas por lo general son

34
ecuaciones diferenciales en derivadas parciales como es el caso de la ecuación
del calor de Newton;

Cuando la diferencia de temperaturas entre un cuerpo y su medio ambiente no es

demasiado grande, el calor transferido en la unidad de tiempo hacia el cuerpo o
desde el cuerpo por conducción, convección y radiación es aproximadamente
proporcional a la diferencia de temperatura entre el cuerpo y el medio externo

Donde a es el coeficiente de intercambio de calor y S es el área del cuerpo.

Si la temperatura T del cuerpo es mayor que la temperatura del medio ambiente

Ta, el cuerpo pierde una cantidad de calor dQ en el intervalo de tiempo
comprendido entre t y t+dt, disminuyendo su temperatura T en dT.

dQ=-m·c·dT

donde m=r V es la masa del cuerpo (r es la densidad y V es el volumen), y c el

calor específico.

La ecuación que nos da la variación de la temperatura T del cuerpo en función del

tiempo es

o bien,

Integrando esta ecuación con la condición inicial de que en el instante t=0, la

temperatura del cuerpo es T0.

35
Obtenemos la relación lineal siguiente.

ln(T-Ta)=-k·t +ln(T0-Ta)

Despejamos T

CONDUCCION DEL CALOR

Otra ecuación diferencial que podemos usar en Ingeniería mecánica es la llamada

ecuación diferencial de conducción de calor, la cual se obtiene a partir de la Ley
de Fourier.

Ley de Fourier

36
Sea J la densidad de corriente de energía (energía por unidad de área y por
unidad de tiempo), que se establece en la barra debido a la diferencia de
temperaturas entre dos puntos de la misma. La ley de Fourier afirma que hay una
proporcionalidad entre el flujo de energía J y el gradiente de temperatura.

Siendo K una constante característica del material denominada conductividad

térmica.

Consideremos un elemento de la barra de longitud dx y sección S. La energía que

entra en el elemento de volumen en la unidad de tiempo es JS, y la que sale es
J’S. La energía del elemento cambia, en la unidad de tiempo, en una cantidad
igual a la diferencia entre el flujo entrante y el flujo saliente.

Esta energía, se emplea en cambiar la temperatura del elemento. La cantidad de

energía absorbida o cedida (en la unidad de tiempo) por el elemento es igual al
37
producto de la masa de dicho elemento por el calor específico y por la variación de
temperatura.

Igualando ambas expresiones, y teniendo en cuenta la ley de Fourier, se obtiene

la ecuación diferencial que describe la conducción térmica

Solución analítica

Supongamos una barra metálica de longitud L, conectada por sus extremos a dos
focos de calor a temperaturas Ta y Tb respectivamente. Sea T0 la temperatura
inicial de la barra cuando se conectan los focos a los extremos de la barra.

Al cabo de cierto tiempo, teóricamente infinito, que en la práctica depende del tipo
de material que empleamos, se establece un estado estacionario en el que la
temperatura de cada punto de la barra no varía con el tiempo. Dicho estado está
caracterizado por un flujo J constante de energía. La ley de Fourier establece que
la temperatura variará linealmente con la distancia x al origen de la barra.

Para describir el estado transitorio buscamos una solución de la forma T(x,

t)=F(x)·G(t), variables separadas

38
El signo negativo asegura el carácter transitorio.

Integramos la primera ecuación diferencial

Integramos la segunda ecuación diferencial

Es una ecuación diferencial similar a la de un Movimiento Armónico Simple, cuya

solución es a·sen(ωx+δ)

La temperatura en cualquier punto x a lo largo de la barra, en un instante

determinado, T(x, t) es la solución de la ecuación diferencial, que es una
combinación de dos términos, la que corresponde al régimen permanente más la
del régimen transitorio.

Condiciones de contorno

 En x=0, T(0, t)=Ta, temperatura fija del extremo izquierdo de la barra

39
  En x=L, T(L, t)=Tb, temperatura fija del extremo derecho de la barra

El régimen variable general de temperaturas de la barra es

Distribución inicial de temperaturas

Solamente, queda por determinar los coeficientes an, identificando esta solución
con la distribución inicial de temperaturas en la barra T(x, 0)=T0 en el instante t=0.

Más abajo, se proporcionan los detalles del cálculo de los coeficientes an del
desarrollo en serie al lector interesado.

La temperatura en cualquier punto de la barra x, en un instante t, se compone de

la suma de un término proporcional a x, y de una serie rápidamente convergente
que describe el estado transitorio.

40
El valor de α=K/(ρc) nos da una medida de la rapidez con la que el sistema
alcanza el estado estacionario. Cuanto mayor sea α antes se alcanza el estado
estacionario

VIBRACIONES MECANICAS:

El aumento permanente de las potencias en máquinas, junto con una disminución

simultánea de gasto de materiales, y la alta exigencia de calidad y productividad
industrial, hacen que el análisis dinámico de las vibraciones mecánicas en
máquinas e instalaciones industriales sea cada vez más exacto.

Por ello en ingeniería Mecánica se debe ser capaz de trabajar sobre vibraciones,
calcularlas, medirlas, analizar el origen de ellas y aplicar correctivos, para ello se
usan modelos matemáticos que nos ayuden a calcular dichas vibraciones, tales
modelos se logran a partir de la aplicación de algunas ecuaciones diferenciales.

Básicamente las vibraciones se encuentran estrechamente relacionadas con

tolerancias de mecanización, desajustes, movimientos relativos entre superficies
en contacto, desbalances de piezas en rotación u oscilación, etc. con lo que
notamos que en Ingeniería Mecánica es común ver su estudio.

Los fenómenos anteriormente mencionados producen casi siempre un

desplazamiento del sistema desde su posición de equilibrio estable originando una
vibración mecánica.

La mayor parte de vibraciones en máquinas y estructuras son indeseables porque

aumentan los esfuerzos y las tensiones y por las pérdidas de energía que las
acompañan. Además, son fuente de desgaste de materiales, de daños por fatiga y
de movimientos y ruidos molestos.

“Todo sistema mecánico tiene características elásticas, de amortiguamiento y

de oposición al movimiento; unas de mayor o menor grado a otras; pero es

41
 debido a que los sistemas tienen esas características lo que hace que el sistema
vibre cuando es sometido a una perturbación ".

"Toda perturbación se puede controlar, siempre y cuando anexemos bloques

de control cuya función de transferencia sea igual o invertida a la función de
transferencia del sistema".

"Si la perturbación tiene una frecuencia igual a la frecuencia natural del sistema, la
amplitud de la respuesta puede exceder la capacidad física del mismo,
ocasionando su destrucción".

La ecuación general de las vibraciones es:

…………………(1)

Donde Y es la magnitud que sufre variaciones periódicas temporales, P(t) la

variable de reforzamiento o fenómeno incidente de la vibración; a, b, y k son las
constantes características del sistema.

Utilizando transformada de Laplace, tenemos que:

……………………………………….(2)

……………………………………………………(3)

Observamos que la Ecuación 3 nos define la función de transferencia general de

nuestro sistema de vibración, lo cual facilita el modelamiento y elaboración de
simulaciones.

42
MODELAMIENTO GENERAL DEL MOVIMIENTO VIBRACIONAL:

Un fenómeno incidente cualquiera puede provocar vibraciones en un sistema,

por lo tanto, en términos generales y universales, el movimiento vibracional de
un sistema afectado por un fenómeno incidente desencadenador de vibraciones
puede ser:

…………………………………………………………..(4)

En donde a constituye las propiedades del sistema, P(t) la función de forzamiento, Y la

magnitud que varía periódicamente, i identifica la magnitud del
sistema.
Sacando transformada de Laplace tenemos:

……………………………………………………..(5)

……………………………….(6)
La Ecuación 6 identifica la función de transferencia universal de las vibraciones

……………………………………………………………..(7)

43
Así mismo, la ecuación 7 identifica la Ecuación característica universal de las
vibraciones
Cuando la estructura de la vibración es lineal, tenemos que la solución a la
función general de transferencia es

……………………………………………………(8)

…………………………..………………………..(9)
Como:

Tenemos que;

Entonces:

44
EL CABLE COLGANTE
Considere un cable o una cuerda que cuelga de dos ,puntos, A y B (Figura
3.29), no necesariamente al mismo nivel. Asuma 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. Sea C la posición más baja del cable, y escoja los ejes x y
y como en la figura, donde el eje y pasa por C.

45
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 (Figura 3.30), la fuerza horizontal H en C, y la carga
vertical total en el segmento CP del cable denotada por W(X) [Link] se
asume 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 manera de decir 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.

46
Descompongamos la tensión T en dos componentes (líneas punteadas en
la Figura 3.30), la componente horizontal con magnitud T cos , y la componente
vertical con magnitud T sen .
Las fuerzas en la dirección x son H hacia la izquierda y T cose hacia
la derecha, mientras que las fuerzas en la dirección y son W hacia abajo y
T sen hacia arriba. De donde,

…1
Dividiendo, y usando el hecho que tan e = dy/dx = a la pendiente de la tangente en
P tenemos:

…2
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. Por diferenciación de (2) con respecto
a x,

…3
Ahora dW/dx representa el incremento en W por unidad de incremento en
x; ésto es, la carga por unidad de distancia en la dirección horizontal.
La ecuación (3) es fundamental.

47
ALGUNOS PUNTOS DE VISTA DE PROFESIONALES

José Luis Ruiz Ameca

Ingeniero en electrónica y comunicaciones.
Instituto Politécnico Nacional
Ultimo trabajo: TELMEX

Yo no recuerdo haber utilizado las ecuaciones diferenciales en mis trabajos,

porque siempre estuve en departamentos operativos. Pero los instrumentos que
usábamos para mediciones, o los dispositivos electrónicos que se usaban en los
proyectos, desde los más elementales como las compuertas, por ejemplo, hasta
otros de mayor complejidad, estaban regidos o representados por ecuaciones
diferenciales que explicaban su funcionamiento, pero como esa parte era de
diseño, no me toco a mi trabajar con ellas, sin embargo , en física, matemáticas,
aeronáutica , química etc. las ecuaciones diferenciales representan los procesos
que siguen los fenómenos naturales.

José Luis Oviedo

Maestro en ciencias en ingeniería eléctrica con especialidad en control automático
Diferentes Instituciones.
Ultimo trabajo: estudio el doctorado en ingeniería eléctrica, misma especialidad

Pues una de las aplicaciones es en control automático, en el diseño de

controladores, comprobación de estabilidad de los mismos
Aplicaciones en robótica
En el caminado de bípedos

48
Carlos Hernández Ruiz
Estudiante de Ingeniera en Eléctrica y Electrónica
UNAM
Ultimo trabajo: Servicio Social en la Facultad de ingeniería

Yo no he tenido la oportunidad de aplicarlas en el campo laboral, pero las he

usado casi siempre que se trate de flujo de corriente.

Alfonso Ruiz Ameca

Ingeniero Civil
ITESM
Ultimo trabajo: KENWORTH

Se aplican casi todo el tiempo, mientras se trate de la resolución de algún

problema físico, en el área de ingeniería civil se usan por ejemplo en el análisis de
estructuras.

Juan Antonio Oseguera Trejo

Ingeniero en electrónica
ITO
Ultimo trabajo: en varias empresas dando mantenimiento a equipos electrónicos

Pues para resolver problemas diarios que la industria se van presentando, para
calcular valores de algunos parámetros o valores de los instrumentos con los que
se están trabajando.

49
BIBLIOGRAFIA

Ecuaciones Diferenciales con Aplicaciones de Modelado

Dennis G. Zill.

Ecuaciones Diferenciales Aplicadas.

Murray Spiegel 3ed.

Serie de Fundamentos de la Mecánica clásica. Estática.

Juan Ocariz Castelazo.

Ecuaciones Diferenciales, Octava Edición.

Earl D. Rainville

FUENTES ELECTRONICAS
[Link]

[Link]

[Link]

CONCLUSIONES:
Con la realización de este trabajo nos dimos cuenta de la verdadera importancia
de las Ecuaciones diferenciales en Ingeniería, en particular en Mecánica.
Ampliamos nuestros conocimientos de la variedad de aplicaciones que pueden
tener las ecuaciones diferenciales, ya que nos ayudan a plantear un modelo
matemático que represente algún fenómeno físico o de otra índole.

50
Además, tuvimos la oportunidad de saber que piensan algunas personas acerca
de la utilización de las ecuaciones diferenciales en la vida profesional como
ingenieros, y como es que ellos han notado su utilización y su gran valor dentro de
la Ingeniería.
Esta investigación acerca de la importancia de las ecuaciones diferenciales en
Ingeniería mecánica, nos dio un panorama más amplio acerca de este tema, pues
como estudiantes, en ocasiones damos por sentado la aplicación e importancia
de cada materia pensando que tal vez algunas cosas no se ocupen en nuestra
carrera. Sin embargo al indagar un poco más en la relación que tienen las
ecuaciones diferenciales en Ingeniería Mecánica, quedamos ampliamente
convencidos que es fundamental el estudio de las mismas, para lograr en un
futuro como ingenieros resolver los diversos problemas que se nos presenten en
la vida profesional.

51