PROGRAMA DE INGENIERÍA MECÁNICA

LABORATORIO DE MOTORES TÉRMICOS

Cinemática del motor.
Excelente informe, buen José Cervantes Rivera, Nereida Llinás Giraldo,
trabajo Nota:5.0 Jassir Olivares Jiménez, Jesús Sinning Mercado
Profesor: Roberto Algarín. 14-09-2023
Universidad Del Atlántico, Barranquilla

Resumen. Datos Fabricante

Tipo de combustible Gasolina
En el presente documento se encuentra el análisis
cinemático del motor Subaru FB25, el cual cuenta Configuración H4-Boxer
con una configuración tipo Bóxer de cuatro Diámetro Pistón (mm) 94
cilindros, 4 tiempos. Para realizar dicho análisis se Carrera (mm) 90
utilizó el programa Matlab y, a su vez, la Cilindrada (cm3) 2498
información brindada por el fabricante, la cual fue
Relación compresión 10:01
complementada con datos obtenidos por medio de
cálculos numéricos. Potencia Max a 5800 rpm (hp) 170
Torque Max a 4100 rpm (lb ft) 174
Palabras claves: Motor Bóxer, Subaru, posición, Orden de encendido 1-2-3-4
velocidad, aceleración. Imagen 1: Información del motor Subaru FB25[1]

I. Objetivo. Para complementar la información obtenida

anteriormente, se calcularon unos parámetros
Realizar el análisis de posiciones, velocidades y necesarios para realizar el análisis de dicho motor y
aceleraciones de un motor reciprocante cumplir con el objetivo.
configuración tipo bóxer.
Los parámetros faltantes son:
II. Metodología.
• Volumen desplazado:
Para comenzar el análisis del motor Bóxer con cuatro 2
cilindros, primero se buscó un motor comercial que 𝜋 2 𝜋 94 90
𝑉𝐷 = 𝐷 𝑆 = ( 𝑐𝑚) ( 𝑐𝑚)
contara con esta configuración. Esto, con el fin de 4 4 10 10
obtener su ficha técnica y encontrar datos como
potencia, cilindrada, relación de compresión,
diámetro del pistón, carrera, longitud de la biela y 𝑉𝐷 = 624.580 𝑐𝑚3
orden de encendido.

El motor que se eligió fue el Subaru FB25 y en la • Volumen cámara de combustión:

imagen 1 se encuentran los datos suministrados por
el fabricante.
𝑉𝐷 624.580 𝑐𝑚3
𝑉𝑐𝑐 = =
𝑟−1 10 − 1

𝑉𝑐𝑐 = 65.3977 𝑐𝑚3

1
 • Longitud de la manivela:
𝐿1 = 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒𝑙 𝑣𝑒𝑐𝑡𝑜𝑟 𝑚𝑎𝑛𝑖𝑣𝑒𝑙𝑎.
𝑆 90 𝑚𝑚 𝐿2 = 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒𝑙 𝑣𝑒𝑐𝑡𝑜𝑟 𝑏𝑖𝑒𝑙𝑎.
𝐿𝑚 = = 𝜃1 = á𝑛𝑔𝑢𝑙𝑜 𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 𝑣𝑒𝑐𝑡𝑜𝑟 𝑚𝑎𝑛𝑖𝑣𝑒𝑙𝑎.
2 2
𝜃2 = á𝑛𝑔𝑢𝑙𝑜 𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 𝑣𝑒𝑐𝑡𝑜𝑟 𝑏𝑖𝑒𝑙𝑎.
𝜃3 = á𝑛𝑔𝑢𝑙𝑜 𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 𝑝𝑖𝑠𝑡ó𝑛
𝐿𝑚 = 45 𝑚𝑚 𝜔1 = 𝑣𝑒𝑙. 𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑚𝑎𝑛𝑖𝑣𝑒𝑙𝑎.
𝜔2 = 𝑣𝑒𝑙. 𝑎𝑛𝑔𝑢𝑙𝑎𝑟 𝑏𝑖𝑒𝑙𝑎.
n = 𝑟𝑒𝑣𝑜𝑙𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑑𝑒𝑙 𝑚𝑜𝑡𝑜𝑟
𝑡 = 𝑣𝑒𝑐𝑡𝑜𝑟 𝑡𝑖𝑒𝑚𝑝𝑜
• Longitud de la biela:
Ecuaciones de posición:
Para este cálculo se consideró 𝛌 = 𝟎, 𝟑
𝑥1 = 𝑦1 = 0
𝐿𝑚 45 𝑚𝑚
𝐿𝑏 = = 𝑥2 = 𝐿1 ∗ cos 𝜃1
λ 0.3 𝑦2 = 𝐿1 ∗ sin 𝜃1
𝑥3 = 𝑥2 + 𝐿2 ∗ cos 𝜃2
𝑦3 = 𝑦2 − 𝐿2 ∗ sin 𝜃2
𝐿𝑏 = 150 𝑚𝑚
𝜔1 = nπ/30
• Ángulo de desfase entre pistones:
𝜃1 = 𝜔1 ∗ 𝑡
720 720
𝜃= =
𝑍 4 Aplicando ley de senos:

𝜃 = 180° 𝜃2 = sin−1[(𝐿1 ∗ sin 𝜃1 )/𝐿2]

Ecuaciones de velocidad:
➢ Análisis cinemático de motores
reciprocantes tipo Bóxer. V1 = V1𝑥 = V1𝑦 = 0
V2 = V1 + V2/1 = 𝐿1 ∗ 𝜔1
Cada segmento del motor conformado por V2𝑥 = −𝐿1 ∗ 𝜔1 ∗ sin 𝜃1
cigüeñal, biela y cilindro-pistón, se V2𝑦 = 𝐿1 ∗ 𝜔1 ∗ cos 𝜃1
comportan como un mecanismo manivela-
biela- corredera (imagen 2.) V3 = V2 + V3/2 = V2 + 𝐿2 ∗ 𝜔2

Para las componentes x, y del vector velocidad

2 de 3:

V3𝑥 = V2𝑥 − 𝐿2 ∗ 𝜔2 ∗ sin 𝜃2

V3𝑦 = V2𝑦 + V3/2𝑦 = 0

3 De V3𝑦 = 0 podemos despejar 𝜔2

1
𝜔2 = (𝐿1 ∗ 𝜔1 ∗ cos 𝜃1 )/(𝐿2 ∗ cos 𝜃2 )
Imagen 2. Referencia Mec. manivela-biela-
corredera (fuente: [Link])

2
 Ecuaciones de aceleración: ➢ Análisis cinemático cilindro-pistón 2

𝑎1 = 𝑎1𝑥 = 𝑎1𝑦 = 0 𝜃1,2 = 𝜃1 + 2𝜋

𝑎2 = 𝑎1 + 𝑎2/1 = 𝜃2,2 = sin−1[(𝐿1 ∗ sin 𝜃1,2 )/𝐿2]
𝑎2/1 = 𝑎2/1𝑛 + 𝑎2/1𝑡 𝜔1 = 𝜔12 = 𝜔13 = 𝜔14
𝑎2/1𝑛 = 𝐿1 ∗ 𝜔1 2
Posición:
Debido a que trabajamos con revoluciones
constantes, la velocidad angular será constante y 𝑥2,2 = 𝐿1 ∗ cos 𝜃1,2
por lo tanto 𝛼1 = 0: 𝑦2,2 = 𝐿1 ∗ sin 𝜃1,2
𝑥3,2 = 𝑥2,2 + 𝐿2 ∗ cos 𝜃2,2
𝑎2/1𝑡 = 𝐿1 ∗ 𝛼1 = 0 𝑦3,2 = 𝑦2,2 − 𝐿2 ∗ sin 𝜃2,2
𝑎2 = 𝐿1 ∗ 𝜔1 2
𝑎2𝑥 = −𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1
𝑎2𝑦 = −𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1 Velocidad:
𝑎3 = 𝑎2 + 𝑎3/2 V2,2𝑥 = −𝐿1 ∗ 𝜔1 ∗ sin 𝜃1,2
𝑎3/2 = 𝑎3/2𝑛 + 𝑎3/2𝑡 V2,2𝑦 = 𝐿1 ∗ 𝜔1 ∗ cos 𝜃1,2
𝑎3/2𝑛 = 𝐿2 ∗ 𝜔2 2 V3,2𝑥 = V2,2𝑥 − 𝐿2 ∗ 𝜔2,2 ∗ sin 𝜃2,2
𝑎3/2𝑡 = 𝐿2 ∗ 𝛼2 𝜔2,2 = (𝐿1 ∗ 𝜔1 ∗ cos 𝜃1,2 )/(𝐿2 ∗ cos 𝜃2,2 )

Descomponiendo la aceleración 3 en sus Aceleración:

componentes x, y:
𝑎2,2𝑥 = −𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1,2
𝑎3𝑥 = 𝑎2𝑥 + 𝑎3/2𝑥 𝑎2,2𝑦 = −𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1,2
𝑎3𝑥 = −(𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1 ) 𝑎3,2𝑥 = −(𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1,2 )
− ( 𝐿2 ∗ 𝜔2 2 ∗ cos 𝜃2 )
− ( 𝐿2 ∗ 𝜔2,2 2 ∗ cos 𝜃2,2 )
− (𝐿2 ∗ 𝛼2 ∗ sin 𝜃2 )
− (𝐿2 ∗ 𝛼2,2 ∗ sin 𝜃2,2 )
Como el pistón se mueve en una sola dirección
(x) podemos considerar que la aceleración en la [−𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1,2 + 𝐿2 ∗ 𝜔2,2 2 ∗ sin 𝜃2,2 ]
𝛼2,2 =
otra dirección (y) se cancela o se hace 0 𝐿2 ∗ cos 𝜃2,2

𝑎3𝑦 = 0 ➢ Análisis cinemático cilindro-pistón 3

De la ecuación cinemática de 𝑎3𝑦 podemos 𝜃1,3 = 𝜃1 + 3𝜋

obtener la aceleración angular 𝛼2 : 𝜃2,3 = sin−1[(𝐿1 ∗ sin 𝜃1,3 )/𝐿2]

[−𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1 + 𝐿2 ∗ 𝜔2 2 ∗ sin 𝜃2 ] Posición:

𝛼2 =
𝐿2 ∗ cos 𝜃2
𝑥2,3 = 𝐿1 ∗ cos 𝜃1,3
Este análisis cinemático anterior es para un solo 𝑦2,3 = 𝐿1 ∗ sin 𝜃1,3
mecanismo manivela-biela-pistón, por lo cual se 𝑥3,3 = 𝑥2,3 + 𝐿2 ∗ cos 𝜃2,3
debe repetir las ecuaciones para los otros 𝑦3,3 = 𝑦2,3 − 𝐿2 ∗ sin 𝜃2,3
pistones sumando el ángulo de desface entre
pistones.

3
 Velocidad: [−𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1,4 + 𝐿2 ∗ 𝜔2,4 2 ∗ sin 𝜃2,4 ]
𝛼2,4 =
𝐿2 ∗ cos 𝜃2,4
V2,3𝑥 = −𝐿1 ∗ 𝜔1 ∗ sin 𝜃1,3
V2,3𝑦 = 𝐿1 ∗ 𝜔1 ∗ cos 𝜃1,3 Dicho análisis se llevó a cabo haciendo uso del
V3,3𝑥 = V2,3𝑥 − 𝐿2 ∗ 𝜔2,3 ∗ sin 𝜃2,3 programa Matlab, en el cual se realizó un código
𝜔2,3 = (𝐿1 ∗ 𝜔1 ∗ cos 𝜃1,3 )/(𝐿2 ∗ cos 𝜃2,3 ) (Anexo 1) que modela el movimiento de cada pistón
y permitió evaluar la posición, velocidad y
Aceleración: aceleración de cada uno de estos; a su vez, permitió
evaluar la velocidad de rotación del cigüeñal.
𝑎2,3𝑥 = −𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1,3 El análisis de dicho motor se realizó para distintas
revoluciones y se obtuvieron gráficas de resultados
𝑎2,3𝑦 = −𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1,3
para cada una de ellas.
𝑎3,3𝑥 = −(𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1,3 )
− ( 𝐿2 ∗ 𝜔2,3 2 ∗ cos 𝜃2,3 ) III. Resultados.
− (𝐿2 ∗ 𝛼2,3 ∗ sin 𝜃2,3 )
Como resultado del código que se realizó, y con base
2 2
[−𝐿1 ∗ 𝜔1 ∗ sin 𝜃1,3 + 𝐿2 ∗ 𝜔2,3 ∗ sin 𝜃2,3 ] a los parámetros del motor, obtenidos de la ficha
𝛼2,3 = técnica. Se realizaron las siguientes gráficas a
𝐿2 ∗ cos 𝜃2,3
distintas revoluciones para observar su
➢ Análisis cinemático cilindro-pistón 4 comportamiento:

𝜃1,4 = 𝜃1 + 2𝜋 - Posición del pistón vs ángulo del cigüeñal.

𝜃2,4 = sin−1[(𝐿1 ∗ sin 𝜃1,4 )/𝐿2]
- Velocidad del pistón vs ángulo del cigüeñal.
Posición:
- Aceleración del pistón vs ángulo del
𝑥2,4 = 𝐿1 ∗ cos 𝜃1,4 cigüeñal.
𝑦2,4 = 𝐿1 ∗ sin 𝜃1,4
𝑥3,4 = 𝑥2,4 + 𝐿2 ∗ cos 𝜃2,4 Entonces, se tiene que:
𝑦3,4 = 𝑦2,4 − 𝐿2 ∗ sin 𝜃2,4
• Para 1000 rpm:
Velocidad:

V2,4𝑥 = −𝐿1 ∗ 𝜔1 ∗ sin 𝜃1,4

V2,4𝑦 = 𝐿1 ∗ 𝜔1 ∗ cos 𝜃1,4
V3,4𝑥 = V2,4𝑥 − 𝐿2 ∗ 𝜔2,4 ∗ sin 𝜃2,4
𝜔2,4 = (𝐿1 ∗ 𝜔1 ∗ cos 𝜃1,4 )/(𝐿2 ∗ cos 𝜃2,4 )

Aceleración:

𝑎2,4𝑥 = −𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1,4

𝑎2,4𝑦 = −𝐿1 ∗ 𝜔1 2 ∗ sin 𝜃1,4
𝑎3,4𝑥 = −(𝐿1 ∗ 𝜔1 2 ∗ cos 𝜃1,4 )
− ( 𝐿2 ∗ 𝜔2,4 2 ∗ cos 𝜃2,4 )
− (𝐿2 ∗ 𝛼2,4 ∗ sin 𝜃2,4 )
Imagen 3: Gráficas para cada pistón a 1000 rpm.

4
 • Para 4000 rpm:
• Para 2000 rpm:

Imagen 6: Gráficas para cada pistón a 4000 rpm.

Imagen 4: Gráficas para cada pistón a 2000 rpm. • Para 5000 rpm:

• Para 3000 rpm:

Imagen 7: Gráficas para cada pistón a 5000 rpm.

• Para 5800 rpm-Potencia max.

Imagen 5: Gráficas para cada pistón a 3000 rpm.

Imagen 8: Gráficas para cada pistón a 5800 rpm-

potencia max.

5
 V. Referencias.

• Para 4100 rpm-Torque max. [1] Motorreviewer, “Subaru FB25 2.5L Engine
Specs, Problems & Reliability” 2016, Engine
professional. [En línea] Disponible en:
[Link]
gine_id=73 [último acceso: 13 de septiembre,
2023]
[2] D. Hagen, “Subaru FB engines” 2015,
Engine professional. [En línea] Disponible en:
[Link]
Q215_12-[Link] [último acceso: 13 de
septiembre, 2023]

Imagen 9: Gráficas para cada pistón a 4100 rpm-

torque max.

IV. Conclusiones.

Se puede observar en las gráficas de ángulo de

cigüeñal vs posición del pistón, como
efectivamente la carrera de los pistones no varía
al aumentar las RPM del cigüeñal. Y que el
pistón sigue un movimiento armónico simple. Se
puede concluir de las gráficas de ángulo del
cigüeñal vs velocidad del pistón, que a medida
que aumentan las RPM del cigüeñal aumenta la
velocidad máxima que alcanzan los pistones. De
forma análoga al observar las gráficas de ángulo
del cigüeñal vs aceleración del pistón se puede
inferir que al aumentar las rpm del cigüeñal,
aumenta la aceleración máxima que alcanzan los
pistones.

6
VI. Anexos.

1. Código de Matlab:

%Motor boxer 4cilindros

clc; clear; close all;

L1=input('Longitud de la manivela(mm): ') ;

L2=input('longitud de la biela (mm): ' );
n=input('revoluciones del cigueñal (rpm): ');
D=input('Diametro del piston (mm): ');

%calculos preliminares
w=n*pi/30;%velocidad angular (rad/s)
t0=2*pi/w;%periodo (s)
t=linspace(0,10*t0,1000)';%vectortiempo
thetha=w*t;%vector angulo del cigueñal(rad/s)
thetha1=thetha*180/pi;%angulo en grados

%como las x y las y cambian de posicion con el tiempo:

for i=1:1000
%cilindro1
thetha0=thetha(i,1);
[x21 y21 x31 y31 v3x1 a3x1]=pistone(L1,L2,thetha0,n);

X31(i,1)=x31;%vector posision piston

V3x1(i,1)=v3x1;%vector velocidad piston
A3x1(i,1)=a3x1;%vector aceleración piston
%cilindro2
thetha0=thetha(i,1)+2*pi;
[x22 y22 x32 y32 v3x2 a3x2]=pistone(L1,L2,thetha0,n);

X32(i,1)=x32;%vector posision piston

V3x2(i,1)=v3x2;%vector velocidad piston
A3x2(i,1)=a3x2;%vector aceleración piston
%cilindro3
thetha0=thetha(i,1)+3*pi;
[x23 y23 x33 y33 v3x3 a3x3]=pistone(L1,L2,thetha0,n);

X33(i,1)=x33;%vector posision piston

V3x3(i,1)=v3x3;%vector velocidad piston
A3x3(i,1)=a3x3;%vector aceleración piston
%cilindro4
thetha0=thetha(i,1)+pi;

7
 [x24 y24 x34 y34 v3x4 a3x4]=pistone(L1,L2,thetha0,n);

X34(i,1)=x34;%vector posision piston

V3x4(i,1)=v3x4;%vector velocidad piston
A3x4(i,1)=a3x4;%vector aceleración piston

%VISTA FRONTAL
%cilindro1
subplot(2,1,1)
hold off
plot([0 x21 x31],[0 y21 y31],'b')%manibela biela
hold on
plot([x31-0.25*D x31-0.25*D x31+0.25*D x31+0.25*D x31-0.25*D],[-0.5*D 0.5*D 0.5*D -
0.5*D -0.5*D],'b')%PISTON
plot([L2-L1-0.3*D L2+L1+0.3*D L2+L1+0.3*D L2-L1-0.3*D],[-0.5*D -0.5*D 0.5*D
0.5*D],'k')%CILINDRO
%cilindro3
plot([0 -x22 -x32],[0 -y22 -y32],'r')%manibela biela
plot([-x32+0.25*D -x32+0.25*D -x32-0.25*D -x32-0.25*D -x32+0.25*D],[0.5*D -0.5*D -
0.5*D 0.5*D 0.5*D],'r')%PISTON
plot([-L2+L1+0.3*D -L2-L1-0.3*D -L2-L1-0.3*D -L2+L1+0.3*D],[0.5*D 0.5*D -0.5*D -
0.5*D],'k')%CILINDRO
%cilindro2
plot([0 x23 x33],[0 y23 y33],'y')%manibela biela
plot([x33-0.25*D x33-0.25*D x33+0.25*D x33+0.25*D x33-0.25*D],[-0.5*D 0.5*D 0.5*D -
0.5*D -0.5*D],'y')%PISTON
%cilindro4
plot([0 -x24 -x34],[0 -y24 -y34],'m')%manibela biela
plot([-x34+0.25*D -x34+0.25*D -x34-0.25*D -x34-0.25*D -x34+0.25*D],[0.5*D -0.5*D -
0.5*D 0.5*D 0.5*D],'m')%PISTON
axis([-4*L1-0.5*L2-0.5*D 2*L1+L2+0.5*L2+0.5*D,-4*L1 2*D])
pause(0.01)

%VISTA SUPERIOR
%cilindro1
subplot(2,1,2)
hold off
plot([0 0 x31-L2 x31-L2 0 0 -x32+L2 -x32+L2 0 0 x33-L2 x33-L2 0 0 -x34+L2 -x34+L2 0
0],[-0.5*D -0.17*D -0.17*D 0.16*D 0.16*D 0.58*D 0.58*D 0.91*D 0.91*D 1.33*D 1.33*D
1.66*D 1.66*D 2.08*D 2.08*D 2.41*D 2.41*D 2.74*D],'k')%manibela
hold on
plot([x31-L2 x31],[0 0],'b')%biela
plot([x31-0.25*D x31-0.25*D x31+0.25*D x31+0.25*D x31-0.25*D],[-0.5*D 0.5*D 0.5*D -
0.5*D -0.5*D],'b')%PISTON

8
 plot([L2-L1-0.3*D L2+L1+0.3*D L2+L1+0.3*D L2-L1-0.3*D],[-0.5*D -0.5*D 0.5*D
0.5*D],'k')%CILINDRO
%cilindro3
plot([-x32+L2 -x32],[0.75*D 0.75*D],'r')%biela
plot([-x32+0.25*D -x32+0.25*D -x32-0.25*D -x32-0.25*D -x32+0.25*D],[0.5*D -0.5*D -
0.5*D 0.5*D 0.5*D]+0.75*D,'r')%PISTON
plot([-L2+L1+0.3*D -L2-L1-0.3*D -L2-L1-0.3*D -L2+L1+0.3*D],[0.5*D 0.5*D -0.5*D -
0.5*D]+0.75*D,'k')%CILINDRO
%cilindro2
plot([x33-L2 x33],[1.5*D 1.5*D],'y')%biela
plot([x33-0.25*D x33-0.25*D x33+0.25*D x33+0.25*D x33-0.25*D],[-0.5*D 0.5*D 0.5*D -
0.5*D -0.5*D]+1.5*D,'y')%PISTON
plot([L2-L1-0.3*D L2+L1+0.3*D L2+L1+0.3*D L2-L1-0.3*D],[-0.5*D -0.5*D 0.5*D
0.5*D]+1.5*D,'k')%CILINDRO
%cilindro4
plot([-x34+L2 -x34],[2.25*D 2.25*D],'m')%biela
plot([-x34+0.25*D -x34+0.25*D -x34-0.25*D -x34-0.25*D -x34+0.25*D],[0.5*D -0.5*D -
0.5*D 0.5*D 0.5*D]+2.25*D,'m')%PISTON
plot([-L2+L1+0.3*D -L2-L1-0.3*D -L2-L1-0.3*D -L2+L1+0.3*D],[0.5*D 0.5*D -0.5*D -
0.5*D]+2.25*D,'k')%CILINDRO
axis([-4*L1-0.5*L2-0.5*D 2*L1+L2+0.5*L2+0.5*D,-2*L1 2*L1+L2+D])
pause(0.01)
end

subplot(3,1,1)
plot(thetha1(1:100,1),X31 (1:100,1),'b')
hold on
plot(thetha1(1:100,1),X32 (1:100,1),'r')
plot(thetha1(1:100,1),X33 (1:100,1),'y')
plot(thetha1(1:100,1),X34 (1:100,1),'m')
title('angulo del cigueñal vs posicion del piston')
legend('piston1','piston3','piston2','piston4')

subplot(3,1,2)
plot(thetha1 (1:100,1) ,V3x1 (1:100,1),'b')
hold on
plot(thetha1 (1:100,1) ,V3x2 (1:100,1),'r')
plot(thetha1 (1:100,1) ,V3x3 (1:100,1),'y')
plot(thetha1 (1:100,1) ,V3x4 (1:100,1),'m')
title('angulo del cigueñal vs velocidad del piston')
legend('piston1','piston3','piston2','piston4')

subplot(3,1,3)
plot(thetha1 (1:100,1) ,A3x1 (1:100,1),'b')
hold on

9
 plot(thetha1 (1:100,1) ,A3x2 (1:100,1),'r')
plot(thetha1 (1:100,1) ,A3x3 (1:100,1),'y')
plot(thetha1 (1:100,1) ,A3x4 (1:100,1),'m')
title(['angulo del cigueñal vs aceleraciòn del piston'])
legend('piston1','piston3','piston2','piston4')

1.1 Función:

function[x2 y2 x3 y3 v3x a3x]=pistone(L1,L2,thetha,n)

w=n*pi/30;%velocidad angular (rad/s)
thetha2=asin(L1*sin(thetha)/L2);%Angulo de la biela
w2=L1*cos(thetha)*w/L2/cos(thetha2); %vel angular de la biela
alfa2=(-L1*w^2*sin(thetha)+L2*w2^2*sin(thetha2))/L2/cos(thetha2);

x2=L1*cos(thetha);
y2=L1*sin(thetha);
x3=L1*cos(thetha)+L2*cos(thetha2);
y3=L1*sin(thetha)-L2*sin(thetha2);

v3x=-L1*w*sin(thetha)-L2*w2*sin(thetha2);%velocidad del piston

a3x=-L1*w*w*cos(thetha)-L2*w2*w2*cos(thetha2)-L2*alfa2*sin(thetha2); %Acel del piston

2. Vista frontal y superior del motor bóxer obtenida del código de Matlab.

4
2
3
1

Imagen 10: Grafica frontal y superior del movimiento del motor bóxer.

10