Análisis dinámico del pistón en motores de

combustión interna utilizando un enfoque

práctico con Python
Esdiver González Patiño
Emmanuel Rivera
Juan Fernando Vallejo

Estudiantes de pregrado Universidad de Antioquia

October 26, 2024

Abstract
En este artı́culo se busca asimilar el uso de variables en el programa
Python para ası́ resolver un problema de ingenierı́a mecánica. Con la
ayuda de esta herramienta se facilita el proceso de resolución de prob-
lemas. En este desarrollo se utilizaron gráficas para ilustrar el cambio
de volumen en el interior de una cámara de combustión para un motor
genérico de N-max. De forma esquemática este sistema se presenta en
la figura 1; donde “ r “ es el radio del cigüeñal, “ l ” es la longitud de
la biela, ”e” es la excentricidad entre el centro del pistón y el centro del
cigüeñal, ”θ” es el ángulo recorrido por el cigüeñal y ”s” es el desplaza-
miento vertical del pistón. El objetivo es poder resolver este anterior
problema relacionando el lenguaje de programación, se utilizarán datos
de análisis y cálculos.
Palabras claves: Pistón, Motor, Biela, Cigueñal, Cámara de combustión,
Excentricidad, Puntos muertos del pistón

1 Introducción
La dinámica de los motores de combustión interna es la que permite identificar
cuales son las causas del movimiento, esto es fundamental para optimizar su
desempeño en aplicaciones mecánicas tales como la transformación de la energı́a
con el fin de efectuar un trabajo. Para el desarrollo de un planteamiento que
permita entender cómo varı́a la posición del pistón, su velocidad, su aceleración
y el volumen de la cámara de combustión a lo largo del ciclo del motor, teniendo
en cuenta cada una de las variables asociadas a la ecuación 1. que describe el
desplazamiento vertical del pistón , se utiliza un software de programación que
facilita la lectura de los resultados y las conclusiones a las que se puedan llegar.
Este artı́culo está basado en el desarrollo de un código en Spyder para la
simulación del movimiento del pistón en un motor de combustión interna, el
motor que se ha tomado como referencia es el de una N-max para hallar su
velocidad, aceleración y volumen respecto a la posición en función del ángulo θ.

1
 A continuación, se hará la discusión del paso a paso para la construcción
del código que ayuda a resolver este problema, utilizando estructuras de ciclos y
funciones del lenguaje Python, además se presentan los resultados de ejecutar el
programa, los cuales son gráficas posición-ángulo, velocidad-ángulo, aceleración-
ángulo y volumen-ángulo. Estas se utilizarán en el análisis para entender cómo
cambia cada una de las variables respecto del ángulo y cuál es su significado.
La figura 1, representa el problema descrito para el movimiento del pistón.

Figure 1: Esquema del problema

2 Metodologı́a
La metodologı́a implementada consistió en el uso de variables, bloques de deci-
siones y subprogramas que proporciona Python para analizar el comportamiento
dinámico de un pistón en el interior de una cámara de combustión de un motor
genérico, en este caso el de una motocicleta N-max. Los resultados se encuentran
expresados en términos de la posición del pistón para cada ángulo, la velocidad
del pistón para cada ángulo, la cual se puede calcular como el cambio de posición
por unidad de tiempo, la aceleración del pistón para cada ángulo, que se puede
calcular como el cambio de velocidad por unidad de tiempo y finalmente el vol-
umen de la cámara de combustión para cada ángulo y esto será más fácil de
interpretarlo a través de gráficas que expresan los resultados de cada valor en

2
incrementos de 0.1° para cualquier ángulo desde 0° hasta 720°.
La estructura del código desarrollado consta de usar la definición de al-
gunas funciones, de las cuales se logra hacer determinados cálculos como por
ejemplo, (s): desplazamiento vertical del pistón, volumen de la cámara de com-
bustión, Listas de posiciones, velocidades, aceleraciones y volúmenes desde i=0
y el gráfico correspondiente a cada situación. Posteriormente se define la función
principal, donde se almacenan los valores de entrada, se digita las condiciones
para cada entrada por medio de condicionales if y while, también se utiliza la
estructura de decisión break, que permite salir de un bucle, raise para declarar
un error en caso de que no se cumpla con las condiciones adecuadas para que
las funciones se ejecuten correctamente y obtener los valores de las variables
nombradas en cada subprograma.
A continuación se presenta con más detalle el procedimiento que se llevó a
cabo con cada fracción de código y su respectiva documentación.

Figure 2: Definir funciones

3
 Figure 3: Definir función

Figure 4: Condiciones para la simulación del motor

4
Figure 5: Plots posición, velocidad, aceleración respecto de θ

Figure 6: Plot volumen respecto de θ

5
 Figure 7: Función principal y condicionales

Figure 8: Bucles a través de condicionales para solicitar entradas

Esta es la representación del código que sirve para la solución de cualquier

problema que el usuario quiera abordar, conociendo las dimensiones geométricas
de un motor de cualquier referenncia en el cual la velocidad angular siempre es
1000 rpm.

3 Resultados y discusión
Esta tabla describe las diferentes variables que se utilizaron en este trabajo.

6
 Parámetro Magnitud
Longitud de la biela 130 mm
Radio del cigueñal 48 mm
Excentricidad 0
Volumen de la cámara de combustión 20000 mm3
Diametro del pistón 58 mm
Velocidad angular 1000 rpm

Table 1: Valores de entrada

La expresión del comportamiento del desplazamiento vertical del pistón (s)

se indica en la siguiente ecuación:

 s 
 2
l e − r cos θ
s = r sin θ + 1− 
r l

El comportamiento del pistón, respecto a su posición, velocidad y aceleración

en función del ángulo del cigüeñal en un motor de combustión interna para
una motocicleta de marca comercial, la cual corresponde a una N-max, además
del volumen respecto del ángulo, se visualiza por medio de las gráficas que se
presentan en este artı́culo, ası́ como el análisis detallado de cada una de los
resultados obtenidos para relacionar los conceptos de la ingenierı́a mecánica con
la relevancia de los métodos de análisis basados en programación para la toma
de decisiones y consideraciones que se puedan hacer en el futuro con el objetivo
de buscar mayor eficiencia de los diseños en motores de combustión interna.

Figure 9: Posición-Ángulo

De la figura 8 podemos observar que el pistón sigue un movimiento a través

de una curva similar a la que se forma con la función seno. Esta primera

7
observación es crucial para entender que la posición del pistón es oscilatoria en
función del ángulo, es decir, que el movimiento del pistón es hacia arriba y abajo,
pero siempre se mueve en valores de posición positivos. De las propiedades del
seno, se sabe que esta función es periódica y este mismo hallazgo se puede
deducir de la posición del pistón conforme el cigüeñal se mueve. El pistón
adquiere puntos máximos cuando el ángulo que forma el cigüeñal con el pistón
es aproximadamente 100° ( 95 π radianes) y 450° ( 52 π radianes). El pistón adquiere
puntos mı́nimos en 280° ( 14 7
9 π radianes) y en 630° ( 2 π radianes).
El comportamiento es oscilatorio, regular y periódico en la simulación de
este motor y este hecho se debe en parte a que la velocidad angular (ω
) es 1000 rpm y es constante, si las revoluciones por minuto no fuesen constantes,
entonces la curva de la Figura 8. podrı́a variar considerablemente.

Figure 10: Velocidad-Ángulo

La velocidad del pistón con respecto al ángulo en la figura 9 también sigue

un patrón, alcanzando valores máximos positivos y negativos. La velocidad del
pistón es cero cuando se encuentra en su posición más alta, esto se conoce como
punto muerto superior y más baja en su punto muerto inferior, alrededor de los
100, 280, 450 y 630 grados.
Se encuentra que el pistón tiene su máxima velocidad en la mitad de su
recorrido, cuando pasa por el centro del cilindro, es decir cuando ha completado
una revolución. A medida que se acerca al punto muerto superior e inferior, la
velocidad disminuye hasta tal punto de ser nula, lo que concuerda con el ciclo de
movimiento de un pistón en el motor. La cuerva tiene caracterı́sticas similares
a la de posición, por lo tanto se puede concluir que la velocidad es oscilatoria,
puede ser positiva o negativa y también es periódica.

8
 Figure 11: Aceleración-Ángulo

La aceleración muestra varios cambios significativos, con picos positivos y

negativos. En particular, se observan valores de aceleración muy altos cuando la
velocidad cambia repentinamente (con grandes intervalos de aceleración), como
cuando el pistón cambia de dirección en el punto muerto superior y en el punto
muerto inferior.
Se considera que este comportamiento indica que la aceleración es mayor
cuando el pistón está en el punto muerto superior e inferior, donde actualmente
se detiene antes de cambiar de dirección. Por lo tanto, el comportamiento es
muy sensible a la velocidad, ya que la aceleración cambia drásticamente al frenar
o acelerar repentinamente el pistón. Se puede considerar la gráfica de la posición
con respecto al ángulo tetha, entonces se encuentra que la información obtenida
es concordante; puesto que si esta gráfica se analiza aplicando el criterio de
primera y segunda deriva es posible deducir en que punto la velocidad es cero,
es máxima o mı́nima y esto mismo para la aceleración.
En los puntos correspondientes a la aceleración máxima, los cuales son de
200 a 350 grados aproximadamente, la aceleración trata de estabilizarse, por lo
tanto se podrı́a decir que en la curva hay una aceleración marcada y esto sucede
para cada revoluión, es decir, cada vez que se complete un ciclo de 360 grados.

9
 Figure 12: Volumen-Ángulo

El volumen en la cámara de combustión cambia periódicamente y alcanza

el volumen máximo cuando el pistón está en el punto muerto inferior (dentro
del rango especificado de 280 grados y 630 grados) y cuando el pistón está en el
punto muerto superior (aproximadamente 100 y 450 grados) llega a su volumen
mı́nimo.
Esto corresponde a un ciclo del motor donde el volumen cambia a medida
que se mueve el pistón. Cuando el pistón está en su posición más alta está en
su volumen más pequeño y entonces se produce la combustión. El volumen es
mayor cuando el pistón está en su posición más baja, lo que permite que entre
aire junto con la mezcla de combustible, es decir, en el momento de admisión,
compresión, explosión y escape.
Para llegar a estos resultados únicamente se consideraron variables asociadas
con la geometrı́a del motor y con su funcionamiento interno, sin embargo, no se
tuvo en cuenta otros aspectos que serı́a importante revisar, como la fricción en-
tre piezas mecánicas, las pérdidas de energı́a por calor, tolerancias en el diseño,
entonces la simulación de un motor de combustión interna es un problema com-
plejo y a la hora de poder encontrar aspectos por mejorar serı́a importante tener
en cuenta cada uno de estos parámetros.

4 Conclusión
Existen factores importantes en el estudio de un motor que afectan la posición,
velocidad y aceleración del pistón, ası́ como el volumen de la cámara de com-
bustión, estos factores están relacionados con la geometrı́a del motor y la rotación
del cigüeñal, los cuales deben considerarse al momento de diseñar. El uso de
programación como Python, para el manejo de datos de ingenierı́a, ası́ como de
los cálculos, es una ventaja en el análisis del movimiento del pistón de un motor
realizado en este estudio se puede generalizar como un método eficiente para la
simulación de problemas y la solución de los mismos, en este caso, se pueden
identificar momentos en los que el motor experimenta cambios más extremos
debido a la fuerza que se produce durante los ciclos del motor, lo cual es clave
para mejorar el rendimiento del motor y reducir el desgaste mecánico. Si se

10
conoce al detalle el comportamiento de los sistemas en el motor, entonces se
puede considerar cómo reducir pérdidas de energı́a en forma de calor utilizando
materiales especı́ficos, el desarrollo de lubricantes para aumentar la vida de
funcionamiento de las piezas ası́ como la fı́sica de cada parte, subensamble y
ensamble. Por medio de la implementación del código desarrollado es que se
hacen los hallazgos, identificando patrones y el comportamiento del sistema.

References
[1] S. W. Kwak, J. K. Shim, and Y. K. Mo, “Kinematic conceptual design of in-
line four-cylinder variable compression ratio engine mechanisms considering
vertical second harmonic acceleration,” Applied Sciences, vol. 10, no. 11, p.
3765, 2020.

[2] H. Nigus, “Kinematics and load formulation of engine crank mechanism,”

Mechanics, Materials Science & Engineering Journal, 2015.
[3] M. Ranjbarkohan, M. Rasekh, A. H. Hoseini, K. Kheiralipour, and M. R.
Asadi, “Kinematics and kinetic analysis of the slider-crank mechanism in
otto linear four cylinder z24 engine,” Journal of mechanical engineering re-
search, vol. 3, no. 3, pp. 85–95, 2011.

[1] [2] [3]

11