Why should Yamaha to make a YZF R675?

First at all, this is a simplified exercise. Yamaha engineers have a lot of information and they are smarter than me. My conclusions are mines and perhaps my point of view have several lacks.

This began when as a issue with a friend, he said that it is better 4 cylinders and I say that it is better 3 cylinders.

I based my studies in YZF-R1 2015, and both engines share the same bore.

• YZF-R1 Bore x stroke 79.0 mm x 50.9 mm
• YZF-R6 Bore x stroke 67.0 mm x 42.5 mm
• YZF-R675 Bore x stroke 79.0 mm x 45.9 mm

Mechanical

If you have 3 cylinders you have several advantages:

• Less parts, and therefore less probability of fault
• Engine smaller, 3 cyl thiner than 4 cylinder
• More displacement engine in competition rules, allow 3 cyl 675cc against 4 cyl 600cc

The power if we could calculate with easies equations are:

Mecanical Power

Thermodynamic Power

FAR = Fuel Air ratio

The Thermodynamic Power is the power that it is providing for the fuel. If we relation with the the power measure by dynamometer and we make the relation with Thermodynamic Power we have the engine efficiency

Other part to see it is the movement piston and rod, due to the piston have speed limit and it is near of 40m/s.

But if we make a simplitication to calculate we colud approx:

The next table shows this calculus with the previous equations:

If we suppose two approximations to calculate the YZF-R675:

• The first line is with the same efficiency and the same rpm of maximum power.
• The last line is with the same piston speed for maximum power (admision and exhaust air speed) but the effiency of the YZF-R6.

The two calculates show more power tahn de current YZF-R6, but this is a exercise and this more complicated.

The next step I compared the engines with cranksahft the 3 to 120º and 4 cilinders with crossplane (90º).

The next chart shows the position, speed and acceleration of one piston, left current YZF-R6 and right the YZF-R675

YZF-R6                                                                                          YZF-R675

If we see the two movent, the 675 has more stroke therefore have more speed, this is the reason that in the Supersport championship for 4 cylinder maximun 600cc and 3 cylinder until 675cc.

In this case to 10.000rpm 

• YZF-R6   22.85m/s
• YZF-R675   24.72m/s 

Thereofore in this case, the YZF-R6 take a 7.5% more rpm than YZF-R675, therefore the this is better for 4 cylinder in the same.

Also ins the speed is more, the acceleration is more 

• YZF-R6   min -2.877e04 m/s2 max 1.783e04 m/s2
• YZF-R675   min -3.124e04 m/s2 max 1.908e04 m/s2 

Therefore in this case, the YZF-R6 take a 7.5% less vibrations than YZF-R675, therefore the this is better for 4 cylinder in the same too.

Flame speed and pressure

Into the each cylinder, the pressure is the same, the differences is the time to burn all mixture, beacause, the flame speed is the same, in this calculates, de flame speed is 100m/s, if we want all burn in the top dead centre (TDC), we need begin to burn several degrees before:

• YZF-R6   20.22º (10.000rpm)
• YZF-R675   23.81º (10.000rpm) 

Thereofore in this case, the YZF-R6 take a 15% more quick than YZF-R675, therefore the this is better for 4 cylinder.

In this study the pressure is the same in both cases (16MPa peak 1.05MPa Mean)

YZF-R6                                                                                          YZF-R675

But the force is depending of surface of piston, and therefore the cylinder.

• YZF-R6   259.67 Nm maximum and -74.87 Nm minimum peak, and 20.73 Nm mean
• YZF-R675   394.02 Nm maximum and -131.46 Nm minimum peak, and 30.54 Nm mean

Thereofore in this case, the YZF-R675 take a 47% more powerfull than YZF-R6, however the only have 3 cylinder meanwhile r6 is 4 cylinder, therefore if we add the means.

• YZF-R6   20.73 Nm x 4 = 82.92 Nm
• YZF-R675   30.54 Nm x 3 = 91.63 Nm

Therefore in this case, the advantage decrease and YZF-R675 take a 10% more powerfull than YZF-R6.

Friction in cylinders.

Friction equations, I take simplified calculus of this with a little changes web

• Nring is 3
• e is ring thickness 0.8mm
• Pressure ring is 75.000N/m2

This equations is ideal case, for example is not contemplate the force per accelerations in the piston and this affect to lateral forces, if I have any time I update the calculus.

YZF-R6                                                                                          YZF-R675

• YZF-R6   0 Nm maximum and -123.27 Nm minimum peak, and -17.32 Nm mean
• YZF-R675   0 Nm maximum and -160.41 Nm minimum peak, and -22.62 Nm mean

Thereofore in this case, the YZF-R675 take a 30% more friction than YZF-R6, however the only have 3 cylinder meanwhile r6 is 4 cylinder, therefore if we add the means.

• YZF-R6   -17.32 Nm x 4 = -69.28 Nm
• YZF-R675   -22.62 Nm x 3 = -67.87 Nm

Therefore in this case, the advantage decrease and YZF-R675 take a 2% less friction than YZF-R6, this is due to the piston is bigger and the friction of lateral force increase this.

All cylinders.

Now I show the all cylinders compound, Friction and Torque, and together (calculate to 10.000rpm)

YZF-R6                                                                                          YZF-R675

Torque

• YZF-R6   255.5 Nm maximum and -69.05 Nm minimum peak, and 82.93 Nm mean
• YZF-R675   396.61 Nm maximum and -133.77 Nm minimum peak, and 91.63 Nm mean

Friction

• YZF-R6  0 Nm maximum and -127.49 Nm minimum peak, and -69.29 Nm mean
• YZF-R675  -3.50 Nm maximum and -165.34 Nm minimum peak, and -67.87 Nm mean

Together

• YZF-R6  140.49 Nm maximum and -96.65 Nm minimum peak, and 13.63 Nm mean
• YZF-R675  248.58 Nm maximum and -175.38 Nm minimum peak, and 23.76 Nm mean

And other part important is standard deviation σ, this provide the smooth in the engine.

• YZF-R6   13.63 Nm mean and 66.42 Nm σ , and 20.83% mean/σ
• YZF-R675   23.76 Nm mean and 106.10 Nm σ , and 22.40% mean/σ

With TPS =0 and brake with engine.

• YZF-R6   -18.15 Nm mean and 29.62 Nm σ , and 61.28% mean/σ
• YZF-R675   -17.02 Nm mean and 42.19 Nm σ , and 40.35% mean/σ

The values have errors, but the important is not the values, the important is the relationship between values

The relation between Mean and σ provide the info that the 3 cylinder in acceleration is smoother than 4 cylinders (9% of difference), but the brake engine with 3 cylinder is rougher than 4 cylinders (52% of difference).

This problem is possible to solve with the slipper clutch. This clutch will be issue of other post.

Vibrations of the engine.

The axis used is show in the figure, and code for calculate is show in the end of post.

Inthe next figures show the forces in the supports A and B and the differents parts 

The parts calculates is:

• Cylinders, in the z axis
• Crankshaft, crankpins where connect with the rod
• Crankshaft, counterweights per rod
• Balancer, counter rotation axis 

And finally togethers parts. 

YZF-R6.

Cylinders (Each cylinder and all cylinders)

Crankshaft crankpins (Each crankpins and all cylinders)

Crankshaft, counterweights (Each counterweights and all cylinders)

Balancer, counter rotation axis (Two mass in the extreme of axis)

All parts together

In the next number are the last figures the maximum and minimum values are showed

• Max Force Ax =  0 N
• Min Force Ax = 0 N
• Max Force Az =  1,178.9 N
• Min Force Az = -1,212.8 N

• Max Force Bx =  0 N
• Min Force Bx = 0 N
• Max Force Bz =  1,178.9 N
• Min Force Bz = -1,212.8 N

YZF-R675.

Cylinders (Each cylinder and all cylinders)

Crankshaft crankpins (Each crankpins and all cylinders)

Crankshaft, counterweights (Each counterweights and all cylinders)

Balancer, counter rotation axis (Two mass in the extreme of axis)

All parts together

In the next number are the last figures the maximum and minimum values are showed

• Max Force Ax =  351.5 N
• Min Force Ax = -351.5 N
• Max Force Az =  636.2 N
• Min Force Az = -635.9 N

• Max Force Bx =  351.5 N
• Min Force Bx = -351.5 N
• Max Force Bz =  635.3 N
• Min Force Bz = -634.6 N

If we compares the vibrations we could  see, with balancer to double speed of crankshaft in the R6 the force in A and B would be over 600 N. Therefore the vibrarions is the same

Crankshaft inertia

This part with stroke 60mm have a inertia 0.01 kg·m2, the reason to extrapolate is k.

• YZF-R6  = I' (1 cyl) = 0,0050 kg·m2 ; I' (4 cyl) = 0,0200 kg·m2 ; 
• YZF-R675  = I' (1 cyl) = 0,0058 kg·m2 ; I' (3 cyl) = 0,0175 kg·m2 ; 

The increase mass due to the  mass calcualtes for vibrations:

• YZF-R6 = Crankshaft inertia 1.9457e-04 kg·m2
• YZF-R675 = Crankshaft inertia 7.6361e-04 kg·m2

Therefore

• YZF-R6   I' = 0,0203 kg·m2 ; 
• YZF-R675   I' = 0,0183 kg·m2 ; 

Therefore in this case, the advantage of YZF-R675 take a 10% by inertia forces than YZF-R6.

Conclusion

If  YZF-R675 used the same piston than YZF-R1

1. YZF-R6 7.5% more rpm than YZF-R675, but the YZF-R675 is 12.5% bigger than YZF-R6.
2. YZF-R675 take a 10% more powerfull than YZF-R6, to same rpm.
3. YZF-R675 take a 2% less friction than YZF-R6.
4. YZF-R675 in acceleration is smoother than YZF-R6 (9% of difference)
5. YZF-R6 in brake engine is smoother than YZF-R675 (52% of difference) this is possible to correct with Slipper Clutch.
6. YZF-R675 vibrations is smoother than YZF-R6 (mean 22% of differences without balancer)
7. YZF-R675 lenght os crankshaft is shorter than YZF-R6 (45% of difference)
8. YZF-R675 cranckshaft inertia is smaller than YZF-R6 (10% more mass, more stroke, and more axis, but less cylinders)

And other advanteges are:

• Simplifications in parts because uses de same piston and valves than YZF-R1
• Less parts (one cylinder less)
• Slipper clutch and new materials mitigate the differences between 600cc and 675

My conclusion, is 675cc, it is true that I like this engine, and I would like see the YZF-R675 in the future, (I hope in the close future).


The code is in Matlab (it has several simplifications): Engine = 1; <br /> <br /> %Selected engine for studies<br /> %Engine 1 Cross-plane YZF-R6<br /> %Engine 2 Cross-plane YZF-R1<br /> %Engine 3 Tricylinder YZF-R675<br /> <br /> TPS = 1 ;<br /> <br /> All = 0 ; Crankshaft = 0 ; Cylinder = 0 ; Balanced = 0 ; Balancer = 0 ; <br /> <br /> Crankshaft.Angle = 0 : pi()./5000 : 4*pi() ; %rd<br /> Crankshaft.Speed = 11000 ; %rpm<br /> Crankshaft.Speed1 = 4000 ; %rpm<br /> Crankshaft.Speed2 = 10000 ; %rpm<br /> Crankshaft.Acceleration = Crankshaft.Angle.*0 ; %rd/s^2<br /> Crankshaft.Time = (Crankshaft.Angle)./Crankshaft.Speed(1) ;<br /> <br /> Air.Density = 1.2754 ; % kg/m^3<br /> Air.DynViscosity = 18.27e-6 ; % Pa*s<br /> Air.Temperature = 298 ; % ºK<br /> Air.Gamma = 1.41 ; %Cp/Cv<br /> Air.BurnGamma = 1.3 ;<br /> Air.Pressure = 101325 ; % Pa<br /> Air.R = 8.314 ; % J/(ºK*mol)<br /> Air.HeatCapacity = 20 ; % J/(ºK*mol)<br /> Air.Mol = 0.0289645 ; % kg/mol<br /> Ring.Pressure = 75000 ; %N/m^2<br /> Ring.Thickness = 0.8/1000 ; %m (Ring piston)<br /> Ring.Friction = 0.04 ; %/mu Steel to Steel Lubricate<br /> Fuel.Energy = TPS.*42.4e6 ; %J/kg<br /> Fuel.Ratio = 14.6; %AFR <br /> <br /> Temp.Losses = 30 % (%) of temp losses lineal with the temperature <br /> <br /> switch Engine<br /> <br /> case 1 %Cross-plane YZF-R6<br /> Cylinder.Number = 4 ;<br /> Piston.Stroke = 42.5/1000 ; %stroke of engine one cylindrer<br /> Piston.Bore = 67/1000 ; %Bore of engine one cylindrer<br /> Piston.Height = Piston.Bore.*0.1689 ;<br /> Piston.Mass = 2700.*Piston.Height.*pi().*(Piston.Bore./2).^2 ;<br /> Crankshaft.Cylinder.Angle = [0,90,270,180].*pi()./180 ; %Degrees in crankshaft<br /> Cylinder.Angle = [0,0,0,0].*pi()./180 ;%Degrees of cylinder<br /> Cylinder.Mass = Piston.Mass.*[1,1,1,1]; %Mass associated to cylinder<br /> Crankshaft.Mass = Piston.Mass.*[1,1,1,1];<br /> Cylinder.Compression = 12.7 ;<br /> Balanced.Angle = [180,270,450,360].*pi()./180 ; %Degrees in Balanced in crankshaft<br /> Balanced.Mass = Piston.Mass.*[1.5,1.5,1.5,1.5] ; %Mass in Balanced in crankshaft<br /> Balancer.Mass = Piston.Mass.*[0.396,0,0,0.396] ; %Mass in external Balancer<br /> Balancer.Angle = [150,0,0,-30].*pi()./180 ; %Degrees in external Balancer<br /> Balancer.Position = [0,0,0,1] ; %Position in external Balancer % L crankshaft<br /> Balancer.Speed = -1 ; % facto os speed crankshaft<br /> Rod.Longitude = 3.563.*25.4/1000 ; %Rod longitude between rotation axis.<br /> Valve.Diameter.Exhaust = 23./1000 ;<br /> Valve.Diameter.Intake = 27./1000 ;<br /> Valve.Lenght.Exhaust = Piston.Bore;<br /> Valve.Lenght.Intake = Piston.Bore ;<br /> Cylinder.Volume = pi().*((Piston.Bore./2).^2).*Piston.Stroke; <br /> Cylinder.ExpCamera = Cylinder.Volume./(Cylinder.Compression-1) ;<br /> Cylinder.ExpCameraHeihgt = (Cylinder.ExpCamera./(pi().*((Piston.Bore./2).^2)));<br /> Cylinder.Spark = [0,180,270,450].*pi()./180 ;%Degrees spark of cylinder<br /> <br /> case 2 %Cross-plane YZF-R1<br /> Cylinder.Number = 4 ;<br /> Piston.Stroke = 50.9/1000 ; %stroke of engine one cylindrer<br /> Piston.Bore = 79/1000 ; %Bore of engine one cylindrer<br /> Piston.Height = Piston.Bore.*0.1689 ;<br /> Piston.Mass = 2700.*Piston.Height.*pi().*(Piston.Bore./2).^2 ;<br /> Crankshaft.Cylinder.Angle = [0,90,270,180].*pi()./180 ; %Degrees in crankshaft<br /> Cylinder.Angle = [0,0,0,0].*pi()./180 ;%Degrees of cylinder<br /> Cylinder.Mass = Piston.Mass.*[1,1,1,1]; %Mass associated to cylinder<br /> Crankshaft.Mass = -Piston.Mass.*[1,1,1,1];<br /> Cylinder.Compression = 13.0 ;<br /> Balanced.Angle = [180,270,90,0].*pi()./180 ; %Degrees in Balanced in crankshaft<br /> Balanced.Mass = -Piston.Mass.*[2,2,2,2] ; %Mass in Balanced in crankshaft <br /> Balancer.Mass = Piston.Mass.*[0,0,0,0] ; %Mass in external Balancer<br /> Balancer.Angle = [0,0,0,0].*pi()./180 ; %Degrees in external Balancer<br /> Balancer.Position = [0,0,0,1] ; %Position in external Balancer % L crankshaft<br /> Balancer.Speed = 1 ; % facto os speed crankshaft<br /> Rod.Longitude = Piston.Stroke.*2.07 ; %Rod longitude between rotation axis.<br /> Valve.Diameter.Exhaust = 26.5./1000 ;<br /> Valve.Diameter.Intake = 33./1000 ;<br /> Valve.Lenght.Exhaust = Piston.Bore;<br /> Valve.Lenght.Intake = Piston.Bore ;<br /> Cylinder.Volume = pi().*((Piston.Bore./2).^2).*Piston.Stroke; <br /> Cylinder.ExpCamera = Cylinder.Volume./(Cylinder.Compression-1) ;<br /> Cylinder.ExpCameraHeihgt = (Cylinder.ExpCamera./(pi().*((Piston.Bore./2).^2)));<br /> Cylinder.Spark = [0,180,270,450].*pi()./180 ;%Degrees spark of cylinder<br /> <br /> case 3 %Tricylinder YZF-R675<br /> Cylinder.Number = 3 ;<br /> Piston.Stroke = 45.9./1000 ; %stroke of engine one cylindrer<br /> Piston.Bore = 79./1000 ; %Bore of engine one cylindrer<br /> Piston.Height = Piston.Bore.*0.1689 ; <br /> Piston.Mass = 2700.*Piston.Height.*pi().*(Piston.Bore./2).^2 ;<br /> Crankshaft.Cylinder.Angle = [0,120,240].*pi()./180 ; %Degrees in crankshaft<br /> Cylinder.Angle = [0,0,0].*pi()./180 ;%Degrees of cylinder<br /> Cylinder.Mass = Piston.Mass.*[1,1,1]; %Mass associated to cylinder<br /> Crankshaft.Mass = Piston.Mass.*[1,1,1,1];<br /> Cylinder.Compression = 13.0 ;<br /> Balanced.Angle = [180,300,420].*pi()./180 ; %Degrees in Balanced in crankshaft<br /> Balanced.Mass = Piston.Mass.*[1.57,1.57,1.57] ; %Mass in Balanced in crankshaft<br /> Balancer.Mass = Piston.Mass.*[.25,0,.25] ; %Mass in external Balancer<br /> Balancer.Angle = [150,0,330].*pi()./180 ; %Degrees in external Balancer<br /> Balancer.Position = [0,0,1] ; %Position in external Balancer % L crankshaft<br /> Balancer.Speed = -1 ; % factor of speed crankshaft<br /> Rod.Longitude = Piston.Stroke.*2.07 ; %Rod longitude between rotation axis.<br /> Valve.Diameter.Exhaust = 26.5./1000 ;<br /> Valve.Diameter.Intake = 33./1000 ;<br /> Valve.Lenght.Exhaust = Piston.Bore;<br /> Valve.Lenght.Intake = Piston.Bore ;<br /> Cylinder.Volume = pi().*((Piston.Bore./2).^2).*Piston.Stroke; <br /> Cylinder.ExpCamera = Cylinder.Volume./(Cylinder.Compression-1) ;<br /> Cylinder.ExpCameraHeihgt = (Cylinder.ExpCamera./(pi().*((Piston.Bore./2).^2)));<br /> Cylinder.Spark = [0,240,480].*pi()./180 ;%Degrees spark of cylinder<br /> <br /> otherwise<br /> <br /> end<br /> <br /> Crankshaft.L=Piston.Bore.*Cylinder.Number + (0.03).*(Cylinder.Number-1); %Lenght crankshaft<br /> <br /> Crankshaft.Speed = (2.*pi().*Crankshaft.Speed./60).*Crankshaft.Angle.^0 ; %rd/s<br /> Crankshaft.Speed1 = (2.*pi().*Crankshaft.Speed1./60) ; %rd/s<br /> Crankshaft.Speed2 = (2.*pi().*Crankshaft.Speed2./60) ; %rd/s<br /> <br /> %Hl.Laminar = (32.*mu.*L./(density*Diameter))3.*speed ; % dP diff of pressure Lamninar<br /> %Hl.Turbulence = (32.*mu.*L./(density*Diameter))3.*speed ; % dP diff of pressure Blasius formule<br /> <br /> Cylinder.Piston.Area = ((pi().*(Piston.Bore./2).^2)) ;<br /> <br /> z.rpmCounter = 0 ;<br /> <br /> for Speed = Crankshaft.Speed1 : (Crankshaft.Speed2-Crankshaft.Speed1)./100 : Crankshaft.Speed2 ;<br /> <br /> z.rpmCounter = z.rpmCounter + 1 <br /> <br /> Crankshaft.Speed = Speed.*Crankshaft.Angle.^0 ;<br /> Crankshaft.Time = (Crankshaft.Angle)./Crankshaft.Speed(1) ;<br /> <br /> <br /> z.X = ((Rod.Longitude).^2 - (sin(Crankshaft.Angle).*Piston.Stroke./2).^2).^0.5 ;<br /> z.SIN = (Piston.Stroke./2).*sin(Crankshaft.Angle);<br /> z.COS = (Piston.Stroke./2).*cos(Crankshaft.Angle);<br /> <br /> Cylinder.Piston.Position2 = z.COS + z.X ;<br /> Cylinder.Piston.Position2 = Cylinder.Piston.Position2 - (Cylinder.Piston.Position2.^0).*(max(Cylinder.Piston.Position2)+min(Cylinder.Piston.Position2))./2 ;<br /> Cylinder.Piston.Speed = -z.SIN.*Crankshaft.Speed.*( 1 + z.COS./z.X ); %<br /> Cylinder.Piston.Acceleration = -( (z.COS.*(Crankshaft.Speed.^2) + z.SIN.*Crankshaft.Acceleration).*( 1 + z.COS./z.X ) + z.SIN.*Crankshaft.Speed.*( (-z.SIN.*Crankshaft.Speed.*z.X-z.COS.*(-z.SIN.*Crankshaft.Speed.*z.COS./z.X ))./(z.X.^2) ) );<br /> <br /> Cylinder.Piston.Volume = Cylinder.ExpCamera + Cylinder.Volume-(Cylinder.Volume./2+ pi().*((Piston.Bore./2).^2).*Cylinder.Piston.Position2);<br /> <br /> Cylinder.Piston.Pressure = 0 ;<br /> <br /> Cylinder.Piston.Force.Friction = -abs(3.*Ring.Friction.*Ring.Pressure.*Ring.Thickness.*pi().*Piston.Bore.*sin(Crankshaft.Angle)) ;<br /> <br /> Cylinder.Piston.Intake.Speed = (((Piston.Bore).^2)./((Valve.Diameter.Intake).^2)).*Cylinder.Piston.Speed/2 ;<br /> Cylinder.Piston.Intake.Reynolds = Air.Density.*Cylinder.Piston.Intake.Speed.*Valve.Diameter.Intake./Air.DynViscosity ;<br /> Cylinder.Piston.Intake.Pressure = (64./Cylinder.Piston.Intake.Reynolds).*(Valve.Lenght.Intake./Valve.Diameter.Intake).*((Air.Density.*Cylinder.Piston.Intake.Speed.^2)./2) + Air.Pressure;<br /> <br /> Cylinder.Piston.Compression.Temperature = Air.Temperature.*((Cylinder.Volume+Cylinder.ExpCamera)./Cylinder.Piston.Volume).^(Air.Gamma-1) ;<br /> %Cylinder.Piston.Compression.Pressure = Air.Pressure.*((Cylinder.Volume+Cylinder.ExpCamera)./(Cylinder.Piston.Volume)).^Air.Gamma;<br /> <br /> Air.nMol_1 = (Cylinder.Volume+Cylinder.ExpCamera).*Air.Density./Air.Mol ; %Air.Pressure.*(Cylinder.Volume+Cylinder.ExpCamera)./(Air.R.*Air.Temperature);<br /> <br /> Cylinder.Piston.Compression.Pressure = (Air.nMol_1.*Air.R./(Cylinder.Piston.Volume)).*(Air.Temperature+(Cylinder.Piston.Compression.Temperature-Air.Temperature).*(1-Temp.Losses./100)) ;<br /> <br /> Air.nMol_2 = (128./121).*Air.nMol_1 ;<br /> Energy = Fuel.Energy.*Air.Density.*(Cylinder.Volume+Cylinder.ExpCamera)./Fuel.Ratio ;<br /> <br /> Cylinder.FlameFront.LenghtTotal = ((Piston.Bore./2).^2 + (Cylinder.ExpCamera./(pi().*((Piston.Bore./2).^2))).^2).^0.5 ;<br /> Cylinder.FlameFront.Speed = 100; %(Air.Gamma.*Air.R.*Cylinder.Piston.Compression.Temperature./Air.Mol).^0.5 ;<br /> <br /> Cylinder.FlameFront.Time = Cylinder.FlameFront.LenghtTotal./Cylinder.FlameFront.Speed ;<br /> Cylinder.FlameFront.Degrees = (180.*Crankshaft.Speed(1)./pi()).*Cylinder.FlameFront.Time; %Crankshaft.Angle(2)<br /> <br /> for i = 1 : round(min((Cylinder.FlameFront.Degrees.*pi()./180)./Crankshaft.Angle(2)))<br /> Cylinder.FlameFront.Lenght(i) = i.*Cylinder.FlameFront.LenghtTotal./round(min((Cylinder.FlameFront.Degrees.*pi()./180)./Crankshaft.Angle(2))) ;<br /> <br /> if (Cylinder.ExpCameraHeihgt>Cylinder.FlameFront.Lenght(i))<br /> Cylinder.FlameFront.Volume(i) = pi().*((2./3).*Cylinder.FlameFront.Lenght(i).^3) ;<br /> end<br /> <br /> if ((Piston.Bore./2)>Cylinder.FlameFront.Lenght(i))<br /> Cylinder.FlameFront.Volume(i) = pi().*((Cylinder.FlameFront.Lenght(i).^2).*Cylinder.ExpCameraHeihgt-(Cylinder.ExpCameraHeihgt.^3)/3) ;<br /> end<br /> <br /> if ((Piston.Bore./2)<Cylinder.FlameFront.Lenght(i)) Cylinder.FlameFront.Volume(i) = (Cylinder.ExpCamera/Cylinder.ExpCameraHeihgt).*(Cylinder.FlameFront.Lenght(i).^2 -(Piston.Bore./2).^2 ).^0.5 + pi().*((Cylinder.FlameFront.Lenght(i).^2).*(Cylinder.ExpCameraHeihgt-(Cylinder.FlameFront.Lenght(i).^2 -(Piston.Bore./2).^2 ).^0.5)-(Cylinder.ExpCameraHeihgt.^3)/3) ; end end Cylinder.FlameFront.Volume(i) = Cylinder.ExpCamera ; z.Explosion1 = round(numel(Cylinder.Piston.Compression.Temperature)./2) ; z.Explosion2 = z.Explosion1 - numel(Cylinder.FlameFront.Volume) + 1; Cylinder.Piston.Explosion.Temperature = Cylinder.Piston.Compression.Temperature(z.Explosion2:z.Explosion1) + (Cylinder.FlameFront.Volume./Cylinder.ExpCamera).*(Energy./(Air.nMol_2.*Air.HeatCapacity)) ; %Cylinder.Piston.Explosion.Pressure = (Air.nMol_1.*Air.R.*Cylinder.Piston.Explosion.Temperature)./(Cylinder.Piston.Volume(numel(Cylinder.FlameFront.Volume))) - Cylinder.Piston.Compression.Pressure(numel(Cylinder.FlameFront.Volume)) + Cylinder.Piston.Compression.Pressure(numel(Cylinder.FlameFront.Volume)) ; Cylinder.Piston.Explosion.Pressure = (Air.nMol_1.*Air.R./(Cylinder.Piston.Volume(z.Explosion2:z.Explosion1))).*(Air.Temperature+(Cylinder.Piston.Explosion.Temperature-Air.Temperature).*(1-Temp.Losses./100)) ; % figure(1) % plot (Cylinder.Piston.Explosion.Pressure) % hold on; % grid on; Cylinder.Piston.Power.Pressure = max(Cylinder.Piston.Explosion.Pressure).*((Cylinder.ExpCamera)./(Cylinder.Piston.Volume)).^Air.BurnGamma ; Cylinder.Piston.Power.Temperature = max(Cylinder.Piston.Explosion.Temperature).*((Cylinder.ExpCamera)./(Cylinder.Piston.Volume)).^(Air.BurnGamma - 1 ) ; %Cylinder.Piston.Exhaust.Pressure = Cylinder.Piston.Intake.Pressure(numel(Cylinder.Piston.Pressure)+1 : numel(Crankshaft.Angle)) Cylinder.Piston.Exhaust.Temperature = min(Cylinder.Piston.Power.Temperature) ; Cylinder.Piston.Exhaust.Pressure = min(Cylinder.Piston.Power.Pressure) ; i=0; Air.nMol_3 = Air.nMol_2 ; % for j = round(3.*numel(Crankshaft.Angle)./4) : numel(Crankshaft.Angle) % i = i + 1 ; % Cylinder.Piston.Exhaust.Speed(i) = (Air.BurnGamma.*Air.R.*Cylinder.Piston.Exhaust.Temperature./Air.Mol).^0.5; % Cylinder.Piston.Exhaust.Volume(i) = Cylinder.Piston.Exhaust.Speed(i).*pi().*(Valve.Diameter.Exhaust.^2).*Crankshaft.Time(i) + Cylinder.Piston.Volume(j) ; % Cylinder.Piston.Exhaust.Pressure(i+1) = Air.nMol_3.*Air.R.*Cylinder.Piston.Exhaust.Temperature./(Cylinder.Piston.Exhaust.Volume(i)); % end Cylinder.Piston.Exhaust.Speed = (Air.BurnGamma.*Air.R.*Cylinder.Piston.Exhaust.Temperature./Air.Mol).^0.5; Cylinder.Piston.Exhaust.Volume = Cylinder.Piston.Volume(round(3.*numel(Cylinder.Piston.Volume)./4):numel(Cylinder.Piston.Volume)) + Cylinder.Piston.Exhaust.Speed.*(2.*pi().*((Valve.Diameter.Exhaust./2).^2)).*(Crankshaft.Time(round(3.*numel(Cylinder.Piston.Volume)./4):numel(Cylinder.Piston.Volume)) - Crankshaft.Time(round(3.*numel(Cylinder.Piston.Volume)./4))) ; Cylinder.Piston.Exhaust.Temperature = min(Cylinder.Piston.Power.Temperature).*(Cylinder.Piston.Exhaust.Volume(1)./Cylinder.Piston.Exhaust.Volume).^(Air.BurnGamma-1); Cylinder.Piston.Exhaust.Pressure = Air.nMol_3.*Air.R.*Cylinder.Piston.Exhaust.Temperature./(Cylinder.Piston.Exhaust.Volume); %plot(Cylinder.Piston.Exhaust.Pressure); Cylinder.Piston.Pressure = Cylinder.Piston.Intake.Pressure(1:round(numel(Crankshaft.Angle)./4)) ; Cylinder.Piston.Pressure = horzcat(Cylinder.Piston.Pressure, Cylinder.Piston.Compression.Pressure(numel(Cylinder.Piston.Pressure):round(numel(Crankshaft.Angle)./2)-numel(Cylinder.Piston.Explosion.Pressure)) ) ; Cylinder.Piston.Pressure = horzcat(Cylinder.Piston.Pressure, Cylinder.Piston.Explosion.Pressure) ; Cylinder.Piston.Pressure = horzcat(Cylinder.Piston.Pressure, Cylinder.Piston.Power.Pressure(round(2.*numel(Crankshaft.Angle)./4)+1:round(3.*numel(Crankshaft.Angle)./4))) ; Cylinder.Piston.Pressure = horzcat(Cylinder.Piston.Pressure, Cylinder.Piston.Exhaust.Pressure) ; Cylinder.Piston.Pressure = Cylinder.Piston.Pressure(1:numel(Crankshaft.Angle)) ; Cylinder.Piston.Pressure = Cylinder.Piston.Pressure-Air.Pressure ; Cylinder.Piston.Pressure(isnan(Cylinder.Piston.Pressure)) = 0; %plot(Cylinder.Piston.Pressure) Cylinder.Piston.Force.Pressure = Cylinder.Piston.Pressure.*Cylinder.Piston.Area ; Cylinder.Piston.Force.All = Cylinder.Piston.Force.Pressure;% + Cylinder.Piston.Force.Friction ; Crankshaft.Torque.Piston = sin(Crankshaft.Angle).*Cylinder.Piston.Force.All.*Piston.Stroke./2 ; Crankshaft.Torque.Piston(1)= 0; Crankshaft.Torque.Piston2 = Crankshaft.Torque.Piston ; Crankshaft.Torque.Piston3 = Crankshaft.Torque.Piston ; Crankshaft.Torque.Piston1=0; for i = 1 : numel(Crankshaft.Torque) Crankshaft.Torque.Piston1 = Crankshaft.Torque.Piston1 + Crankshaft.Torque.Piston(i).*Crankshaft.Angle(2)./max(Crankshaft.Angle) ; end Cylinder.Piston.Force.Friction_Adjust = -abs(3.*Ring.Friction.*Ring.Pressure.*Ring.Thickness.*pi().*Piston.Bore.*sin(Crankshaft.Angle)) ; Cylinder.Piston.Force.Friction_Pressure = -abs(3.*Ring.Friction.*Cylinder.Piston.Pressure.*Ring.Thickness.*pi().*Piston.Bore.*sin(Crankshaft.Angle)) ; Cylinder.Piston.Force.Friction_LateralForce = -abs((Cylinder.Piston.Pressure.*Cylinder.Piston.Area)./(1+(z.X./(z.SIN.*Ring.Friction))).*sin(Crankshaft.Angle)) ; Cylinder.Piston.Force.Friction = Cylinder.Piston.Force.Friction_Adjust + Cylinder.Piston.Force.Friction_Pressure + Cylinder.Piston.Force.Friction_LateralForce ; Crankshaft.Torque.Piston = horzcat(Crankshaft.Torque.Piston(round(numel(Crankshaft.Torque.Piston)./2):numel(Crankshaft.Torque.Piston)),Crankshaft.Torque.Piston(1:round(numel(Crankshaft.Torque.Piston)./2)-1)) ; Cylinder.Piston.Force.Friction = horzcat(Cylinder.Piston.Force.Friction(round(numel(Crankshaft.Torque.Piston)./2):numel(Crankshaft.Torque.Piston)),Cylinder.Piston.Force.Friction(1:round(numel(Crankshaft.Torque.Piston)./2)-1)) ; % Cylinder position into the crankshaft longitude for i = 1 : Cylinder.Number Cylinder.Position(i) = (Crankshaft.L)./(2.*(Cylinder.Number)) + (i-1).*(Crankshaft.L)./(Cylinder.Number) ; end for i = 1 : Cylinder.Number Crankshaft.Force.Az( i , : ) = Crankshaft.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*cos(Crankshaft.Angle + Crankshaft.Cylinder.Angle(i)).*(Crankshaft.L - Cylinder.Position(i))./(Crankshaft.L) ; Crankshaft.Force.Bz( i , : ) = Crankshaft.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*cos(Crankshaft.Angle + Crankshaft.Cylinder.Angle(i)).*(Cylinder.Position(i))./(Crankshaft.L) ; Crankshaft.Force.Ax( i , : ) = Crankshaft.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*sin(Crankshaft.Angle + Crankshaft.Cylinder.Angle(i)).*(Crankshaft.L - Cylinder.Position(i))./(Crankshaft.L) ; Crankshaft.Force.Bx( i , : ) = Crankshaft.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*sin(Crankshaft.Angle + Crankshaft.Cylinder.Angle(i)).*(Cylinder.Position(i))./(Crankshaft.L) ; z.Position = round(numel(Crankshaft.Angle).*(Crankshaft.Cylinder.Angle(i)./max(Crankshaft.Angle)))+1 ; Cylinder.Piston.Acceleration2 = horzcat(Cylinder.Piston.Acceleration(z.Position:numel(Crankshaft.Angle)),Cylinder.Piston.Acceleration(1:z.Position-1)); Cylinder.Force.Az( i , : ) = -Cylinder.Mass(i).*cos(Cylinder.Angle(i)).*Cylinder.Piston.Acceleration2.*(Crankshaft.L - Cylinder.Position(i))./(Crankshaft.L) ; Cylinder.Force.Bz( i , : ) = -Cylinder.Mass(i).*cos(Cylinder.Angle(i)).*Cylinder.Piston.Acceleration2.*(Cylinder.Position(i))./(Crankshaft.L) ; Cylinder.Force.Ax( i , : ) = -Cylinder.Mass(i).*sin(Cylinder.Angle(i)).*Cylinder.Piston.Acceleration2.*(Crankshaft.L - Cylinder.Position(i))./(Crankshaft.L) ; Cylinder.Force.Bx( i , : ) = -Cylinder.Mass(i).*sin(Cylinder.Angle(i)).*Cylinder.Piston.Acceleration2.*(Cylinder.Position(i))./(Crankshaft.L) ; Balanced.Force.Az( i , : ) = Balanced.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*cos(Crankshaft.Angle + Balanced.Angle(i)).*(Crankshaft.L - Cylinder.Position(i))./(Crankshaft.L) ; Balanced.Force.Bz( i , : ) = Balanced.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*cos(Crankshaft.Angle + Balanced.Angle(i)).*(Cylinder.Position(i))./(Crankshaft.L) ; Balanced.Force.Ax( i , : ) = Balanced.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*sin(Crankshaft.Angle + Balanced.Angle(i)).*(Crankshaft.L - Cylinder.Position(i))./(Crankshaft.L) ; Balanced.Force.Bx( i , : ) = Balanced.Mass(i).*((Piston.Stroke./2).*(Speed.^2)).*sin(Crankshaft.Angle + Balanced.Angle(i)).*(Cylinder.Position(i))./(Crankshaft.L) ; Balancer.Force.Az( i , : ) = Balancer.Mass(i).*((Piston.Stroke./2).*((Balancer.Speed.*Speed).^2)).*cos(Crankshaft.Angle.*Balancer.Speed + Balancer.Angle(i)).*(Crankshaft.L.*(1 - Balancer.Position(i)))./(Crankshaft.L) ; Balancer.Force.Bz( i , : ) = Balancer.Mass(i).*((Piston.Stroke./2).*((Balancer.Speed.*Speed).^2)).*cos(Crankshaft.Angle.*Balancer.Speed + Balancer.Angle(i)).*(Crankshaft.L.*Balancer.Position(i))./(Crankshaft.L) ; Balancer.Force.Ax( i , : ) = Balancer.Mass(i).*((Piston.Stroke./2).*((Balancer.Speed.*Speed).^2)).*sin(Crankshaft.Angle.*Balancer.Speed + Balancer.Angle(i)).*(Crankshaft.L.*(1 - Balancer.Position(i)))./(Crankshaft.L) ; Balancer.Force.Bx( i , : ) = Balancer.Mass(i).*((Piston.Stroke./2).*((Balancer.Speed.*Speed).^2)).*sin(Crankshaft.Angle.*Balancer.Speed + Balancer.Angle(i)).*(Crankshaft.L.*Balancer.Position(i))./(Crankshaft.L) ; z.PositionSpark = round(numel(Crankshaft.Angle).*(Cylinder.Spark(i)./max(Crankshaft.Angle)))+1 ; Crankshaft.Torque.All( i , : ) = horzcat(Crankshaft.Torque.Piston(numel(Crankshaft.Angle)-z.PositionSpark:numel(Crankshaft.Angle)),Crankshaft.Torque.Piston(1:numel(Crankshaft.Angle)-z.PositionSpark-1)) ; Cylinder.Piston.Position3 ( i , : ) = horzcat(Cylinder.Piston.Position2(numel(Crankshaft.Angle)-z.Position:numel(Crankshaft.Angle)),Cylinder.Piston.Position2(1:numel(Crankshaft.Angle)-z.Position-1)) ; Cylinder.Piston.Torque.Friction ( i , : ) = horzcat(Cylinder.Piston.Force.Friction(numel(Crankshaft.Angle)-z.PositionSpark:numel(Crankshaft.Angle)),Cylinder.Piston.Force.Friction(1:numel(Crankshaft.Angle)-z.PositionSpark-1)) ; Crankshaft.Torque.All( i , : ) = Crankshaft.Torque.All( i , : ) ; end Crankshaft.Force.Az( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Crankshaft.Force.Bz( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Crankshaft.Force.Ax( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Crankshaft.Force.Bx( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Cylinder.Force.Az( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Cylinder.Force.Bz( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Cylinder.Force.Ax( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Cylinder.Force.Bx( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balanced.Force.Az( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balanced.Force.Bz( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balanced.Force.Ax( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balanced.Force.Bx( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balancer.Force.Az( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balancer.Force.Bz( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balancer.Force.Ax( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Balancer.Force.Bx( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Crankshaft.Torque.All( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; Cylinder.Piston.Torque.Friction ( Cylinder.Number + 1 , : ) = 0*Crankshaft.Force.Az( 1 , : ) ; for i = 0 : Cylinder.Number - 1 Crankshaft.Force.Az( Cylinder.Number + 1 , : ) = Crankshaft.Force.Az( Cylinder.Number + 1 , : ) + Crankshaft.Force.Az( i + 1 , : ) ; Crankshaft.Force.Bz( Cylinder.Number + 1 , : ) = Crankshaft.Force.Bz( Cylinder.Number + 1 , : ) + Crankshaft.Force.Bz( i + 1 , : ) ; Crankshaft.Force.Ax( Cylinder.Number + 1 , : ) = Crankshaft.Force.Ax( Cylinder.Number + 1 , : ) + Crankshaft.Force.Ax( i + 1 , : ) ; Crankshaft.Force.Bx( Cylinder.Number + 1 , : ) = Crankshaft.Force.Bx( Cylinder.Number + 1 , : ) + Crankshaft.Force.Bx( i + 1 , : ) ; Cylinder.Force.Az( Cylinder.Number + 1 , : ) = Cylinder.Force.Az( Cylinder.Number + 1 , : ) + Cylinder.Force.Az( i + 1 , : ) ; Cylinder.Force.Bz( Cylinder.Number + 1 , : ) = Cylinder.Force.Bz( Cylinder.Number + 1 , : ) + Cylinder.Force.Bz( i + 1 , : ) ; Cylinder.Force.Ax( Cylinder.Number + 1 , : ) = Cylinder.Force.Ax( Cylinder.Number + 1 , : ) + Cylinder.Force.Ax( i + 1 , : ) ; Cylinder.Force.Bx( Cylinder.Number + 1 , : ) = Cylinder.Force.Bx( Cylinder.Number + 1 , : ) + Cylinder.Force.Bx( i + 1 , : ) ; Balanced.Force.Az( Cylinder.Number + 1 , : ) = Balanced.Force.Az( Cylinder.Number + 1 , : ) + Balanced.Force.Az( i + 1 , : ) ; Balanced.Force.Bz( Cylinder.Number + 1 , : ) = Balanced.Force.Bz( Cylinder.Number + 1 , : ) + Balanced.Force.Bz( i + 1 , : ) ; Balanced.Force.Ax( Cylinder.Number + 1 , : ) = Balanced.Force.Ax( Cylinder.Number + 1 , : ) + Balanced.Force.Ax( i + 1 , : ) ; Balanced.Force.Bx( Cylinder.Number + 1 , : ) = Balanced.Force.Bx( Cylinder.Number + 1 , : ) + Balanced.Force.Bx( i + 1 , : ) ; Balancer.Force.Az( Cylinder.Number + 1 , : ) = Balancer.Force.Az( Cylinder.Number + 1 , : ) + Balancer.Force.Az( i + 1 , : ) ; Balancer.Force.Bz( Cylinder.Number + 1 , : ) = Balancer.Force.Bz( Cylinder.Number + 1 , : ) + Balancer.Force.Bz( i + 1 , : ) ; Balancer.Force.Ax( Cylinder.Number + 1 , : ) = Balancer.Force.Ax( Cylinder.Number + 1 , : ) + Balancer.Force.Ax( i + 1 , : ) ; Balancer.Force.Bx( Cylinder.Number + 1 , : ) = Balancer.Force.Bx( Cylinder.Number + 1 , : ) + Balancer.Force.Bx( i + 1 , : ) ; Crankshaft.Torque.All( Cylinder.Number + 1 , : ) = Crankshaft.Torque.All( Cylinder.Number + 1 , : ) + Crankshaft.Torque.All( i + 1 , : ) ; Cylinder.Piston.Torque.Friction( Cylinder.Number + 1 , : ) = Cylinder.Piston.Torque.Friction( Cylinder.Number + 1 , : ) + Cylinder.Piston.Torque.Friction( i + 1 , : ) ; %Crankshaft.Torque.All1( Cylinder.Number + 1 , : ) = Crankshaft.Torque.All( Cylinder.Number + 1 , : ) + Cylinder.Piston.Torque.Friction( Cylinder.Number + 1 , : ) ; end All.Force.Az = Balanced.Force.Az + Cylinder.Force.Az + Crankshaft.Force.Az + Balancer.Force.Az; All.Force.Bz = Balanced.Force.Bz + Cylinder.Force.Bz + Crankshaft.Force.Bz + Balancer.Force.Bz; All.Force.Ax = Balanced.Force.Ax + Cylinder.Force.Ax + Crankshaft.Force.Ax + Balancer.Force.Ax; All.Force.Bx = Balanced.Force.Bx + Cylinder.Force.Bx + Crankshaft.Force.Bx + Balancer.Force.Bx; Crankshaft.Torque.All1 = 0 ; Data.Torque.Mean (z.rpmCounter) = mean(Crankshaft.Torque.All(Cylinder.Number + 1 , :)); Data.Torque.Sigma (z.rpmCounter) = std(Crankshaft.Torque.All(Cylinder.Number + 1 , :)); Data.Power.Mean (z.rpmCounter) = mean(Crankshaft.Torque.All(Cylinder.Number + 1 , :)).*Crankshaft.Speed(1)./735.5 ; Data.Power.Sigma (z.rpmCounter) = std(Crankshaft.Torque.All(Cylinder.Number + 1 , :)).*Crankshaft.Speed(1)./735.5 ; Data.Friction.Mean (z.rpmCounter) = mean(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)); Data.Friction.Sigma (z.rpmCounter) = std(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)); Data.CylSpeed (z.rpmCounter) = max(Cylinder.Piston.Speed) ; Data.IntakeSpeed.Mean (z.rpmCounter) = mean(Cylinder.Piston.Intake.Speed(1:round(numel(Cylinder.Piston.Intake.Speed)./4))); Data.IntakeSpeed.Max (z.rpmCounter) = max(Cylinder.Piston.Intake.Speed); Data.Pressure (z.rpmCounter,:) = Cylinder.Piston.Pressure ; %Data.Pressure (z.rpmCounter,:) = Cylinder.Piston.Explosion.Pressure ; Data.rpm (z.rpmCounter) = Speed./(2.*pi()./60) ; Data.SparkDegrees (z.rpmCounter) = Cylinder.FlameFront.Degrees ; Data.Torque_Piston (z.rpmCounter,:) = Crankshaft.Torque.Piston2 ; Data.Torque_Crankshaft (z.rpmCounter,:) = Crankshaft.Torque.All( i + 1 , : ) ; Data.Torque_Friction (z.rpmCounter,:) = Cylinder.Piston.Torque.Friction( Cylinder.Number + 1 , : ) ; %Data.Torque_All.All (z.rpmCounter,:) = Crankshaft.Torque.All1( Cylinder.Number + 1 , : ) ; % Data.Torque_All.Mean (z.rpmCounter) = mean(Crankshaft.Torque.All1(Cylinder.Number + 1 , :)); % Data.Torque_All.Sigma (z.rpmCounter) = std(Crankshaft.Torque.All1(Cylinder.Number + 1 , :)); % Data.Power_All.Mean (z.rpmCounter) = mean(Crankshaft.Torque.All1(Cylinder.Number + 1 , :)).*Crankshaft.Speed(1)./735.5 ; % Data.Power_All.Sigma (z.rpmCounter) = std(Crankshaft.Torque.All1(Cylinder.Number + 1 , :)).*Crankshaft.Speed(1)./735.5 ; end %Piston datas fig=1; figure (fig) hold on; subplot(2,1,1); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position2*1000,Crankshaft.Angle*180/pi(),Cylinder.Piston.Speed); title('Cylinder Position (blue) / Speed (green)') %xlabel('Crankshaft (º)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'m/s') % right y-axis grid on; subplot(2,1,2); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position2*1000,Crankshaft.Angle*180/pi(),Cylinder.Piston.Acceleration); %plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Posit); %plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Speed) %plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position2) hold off; grid on; title('Cylinder Position (blue) /Acceleration (green)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'m/s^2') % right y-axis xlabel('Crankshaft (º)') %Piston pressure and Torque fig=fig+1; figure(fig) hold on; subplot(2,1,1); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position2*1000,Crankshaft.Angle*180/pi(),Cylinder.Piston.Pressure); title('Cylinder Position (blue) / Pressure (green)') %xlabel('Crankshaft (º)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'Pa') % right y-axis grid on; subplot(2,1,2); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position2*1000,Crankshaft.Angle*180/pi(),Crankshaft.Torque.Piston3); hold off; grid on; title('Cylinder Position (blue) / Torque (green)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis xlabel('Crankshaft (º)') %Map of Pressure in one cilinder fig=fig+1; figure(fig) mesh(Crankshaft.Angle*180/pi(), Data.rpm,Data.Pressure) grid on; xlabel('Crankshaft (º)') ylabel('rpm') zlabel('Pa') %Map of Torque in one cilinder fig=fig+1; figure(fig) mesh(Crankshaft.Angle*180/pi(), Data.rpm, Data.Torque_Piston ) grid on; xlabel('Crankshaft (º)') ylabel('rpm') zlabel('Nm') %Map of Torque in one cilinder fig=fig+1; figure(fig) mesh(Crankshaft.Angle*180/pi(), Data.rpm, Data.Torque_Friction ) grid on; xlabel('Crankshaft (º)') ylabel('rpm') zlabel('Nm') %Torque and power f(rpm) fig=fig+1; figure(fig) hold on; subplot(2,1,1); [hAx,hLine1,hLine2] = plotyy(Data.rpm,Data.Power.Mean,Data.rpm,Data.Torque.Mean); title('Engine Power (blue) / Torque (green)') %xlabel('Crankshaft (º)') ylabel(hAx(1),'CV') % left y-axis ylabel(hAx(2),'Nm') % right y-axis grid on; subplot(2,1,2); [hAx,hLine1,hLine2] = plotyy(Data.rpm,Data.Power.Sigma,Data.rpm,Data.Torque.Sigma); hold off; grid on; title('Engine Power \sigma (blue) / Torque \sigma (green)') ylabel(hAx(1),'CV') % left y-axis ylabel(hAx(2),'Nm') % right y-axis xlabel('rpm') %Friction f(rpm) fig=fig+1; figure(fig) hold on; subplot(2,1,1); [hAx,hLine1,hLine2] = plotyy(Data.rpm,Data.Friction.Mean,Data.rpm,Data.Friction.Sigma); title('Engine Friction Torque (blue) / Torque \sigma (green)') %xlabel('Crankshaft (º)') ylabel(hAx(1),'Nm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis grid on; subplot(2,1,2); [hAx,hLine1,hLine2] = plotyy(Data.rpm,((Data.Torque.Mean+Data.Friction.Mean).*(2.*pi().*Data.rpm/60)./735.5),Data.rpm,(Data.Torque.Mean+Data.Friction.Mean)); hold off; grid on; title('Engine With friction Power (blue) / Torque (green)') ylabel(hAx(1),'CV') % left y-axis ylabel(hAx(2),'Nm') % right y-axis xlabel('rpm') fig=fig+1; figure(fig) hold on; subplot(2,1,1); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position3'.*1000,Crankshaft.Angle*180/pi(),Crankshaft.Torque.All(Cylinder.Number + 1 , :)+Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)); title('Cylinders Position / Torque+Friction') %xlabel('Crankshaft (º)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis grid on; subplot(2,1,2); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Crankshaft.Torque.All(Cylinder.Number + 1 , :),Crankshaft.Angle*180/pi(),Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)); hold off; grid on; title('Cylinder Torque (blue) / Friction (green)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis xlabel('Crankshaft (º)') MaxTorque = max(Crankshaft.Torque.All(Cylinder.Number + 1 , :)) MinTorque = min(Crankshaft.Torque.All(Cylinder.Number + 1 , :)) MeanTorque = mean(Crankshaft.Torque.All(Cylinder.Number + 1 , :)) MaxFriction = max(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)) MinFriction = min(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)) MeanFriction = mean(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :)) MaxTorqueFriction = max(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :) + Crankshaft.Torque.All(Cylinder.Number + 1 , :)) MinTorqueFriction = min(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :) + Crankshaft.Torque.All(Cylinder.Number + 1 , :)) MeanTorqueFriction = mean(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :) + Crankshaft.Torque.All(Cylinder.Number + 1 , :)) StdTorqueFriction = std(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :) + Crankshaft.Torque.All(Cylinder.Number + 1 , :)) MeanStdTorqueFriction = mean(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :) + Crankshaft.Torque.All(Cylinder.Number + 1 , :))./std(Cylinder.Piston.Torque.Friction(Cylinder.Number + 1 , :) + Crankshaft.Torque.All(Cylinder.Number + 1 , :)) fig=fig+1; figure(fig) hold on; subplot(2,1,1); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Cylinder.Piston.Position3'.*1000,Crankshaft.Angle*180/pi(),Crankshaft.Torque.All(Cylinder.Number - 1 , :)+Cylinder.Piston.Torque.Friction(Cylinder.Number -1 , :)); title('Cylinder Position / Torque+Friction') %xlabel('Crankshaft (º)') ylabel(hAx(1),'mm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis grid on; subplot(2,1,2); [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Crankshaft.Torque.All(Cylinder.Number - 1 , :),Crankshaft.Angle*180/pi(),Cylinder.Piston.Torque.Friction(Cylinder.Number - 1 , :)); hold off; grid on; title('Cylinder Torque (blue) / Friction (green)') ylabel(hAx(1),'Nm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis xlabel('Crankshaft (º)') fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Crankshaft.Torque.All(i , :),Crankshaft.Angle*180/pi(),Cylinder.Piston.Torque.Friction(i , :)); grid on; ylabel(hAx(1),'Nm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis title(sprintf('Cyl %d Torque(blue) Friction(green)',i)) xlabel('Crankshaft (º)') hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) [hAx,hLine1,hLine2] = plotyy(Crankshaft.Angle*180/pi(),Crankshaft.Torque.All(i+1 , :),Crankshaft.Angle*180/pi(),Cylinder.Piston.Torque.Friction(i+1 , :)); grid on; ylabel(hAx(1),'Nm') % left y-axis ylabel(hAx(2),'Nm') % right y-axis xlabel('Crankshaft (º)') title('All Cyl Torque(blue) Friction(green)') hold off %Crankshaft fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Ax(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Bx(i,:),'r') grid on; title(sprintf('Crankshaft Cyl %d Ax(blue) Bx(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Crankshaft.Force.Bx()))-0.1 max(max(Crankshaft.Force.Bx()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Ax(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Bx(i+1,:),'r') grid on; title('Crankshaft Cyl total Ax(blue) Bx(red)') axis([0 max(Crankshaft.Angle())*180./pi() min(min(Crankshaft.Force.Bx()))-0.1 max(max(Crankshaft.Force.Bx()))+0.1 ]) hold off fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Az(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Bz(i,:),'r') grid on; title(sprintf('Crankshaft Cyl %d Az(blue) Bz(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Crankshaft.Force.Bz()))-0.1 max(max(Crankshaft.Force.Bz()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Az(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Crankshaft.Force.Bz(i+1,:),'r') grid on; title('Crankshaft Cyl total Az(blue) Bz(red)') axis([0 max(Crankshaft.Angle())*180./pi() min(min(Crankshaft.Force.Bz()))-0.1 max(max(Crankshaft.Force.Bz()))+0.1 ]) hold off %Cylinder fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Ax(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Bx(i,:),'r') grid on; title(sprintf('Cylinder %d Ax(blue) Bx(red)',i)); axis([0 max(Crankshaft.Angle())*180./pi() min(min(Cylinder.Force.Bx()))-0.1 max(max(Cylinder.Force.Bx()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Ax(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Bx(i+1,:),'r') grid on; axis([0 max(Crankshaft.Angle())*180./pi() min(min(Cylinder.Force.Bx()))-0.1 max(max(Cylinder.Force.Bx()))+0.1 ]) title('Cylinder total Ax(blue) Bx(red)') hold off fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Az(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Bz(i,:),'r') grid on; title(sprintf('Cylinder %d Az(blue) Bz(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Cylinder.Force.Bz()))-0.1 max(max(Cylinder.Force.Bz()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Az(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Cylinder.Force.Bz(i+1,:),'r') grid on; title('Cylinder total Az(blue) Bz(red)') axis([0 max(Crankshaft.Angle())*180./pi() min(min(Cylinder.Force.Bz()))-0.1 max(max(Cylinder.Force.Bz()))+0.1 ]) hold off %Balanced Crankshaft fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Balanced.Force.Ax(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balanced.Force.Bx(i,:),'r') grid on; title(sprintf('Balanced %d Ax(blue) Bx(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balanced.Force.Bx()))-0.1 max(max(Balanced.Force.Bx()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Balanced.Force.Ax(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balanced.Force.Bx(i+1,:),'r') grid on; axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balanced.Force.Bx()))-0.1 max(max(Balanced.Force.Bx()))+0.1 ]) title('Balanced total Ax(blue) Bx(red)') hold off fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Balanced.Force.Az(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balanced.Force.Bz(i,:),'r') grid on; title(sprintf('Balanced %d Az(blue) Bz(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balanced.Force.Bz()))-0.1 max(max(Balanced.Force.Bz()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Balanced.Force.Az(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balanced.Force.Bz(i+1,:),'r') grid on; title('Balanced total Az(blue) Bz(red)') axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balanced.Force.Bz()))-0.1 max(max(Balanced.Force.Bz()))+0.1 ]) hold off %Balancer axis fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Balancer.Force.Ax(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balancer.Force.Bx(i,:),'r') grid on; title(sprintf('Balancer %d Ax(blue) Bx(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balancer.Force.Bx()))-0.1 max(max(Balancer.Force.Bx()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Balancer.Force.Ax(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balancer.Force.Bx(i+1,:),'r') grid on; axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balancer.Force.Bx()))-0.1 max(max(Balancer.Force.Bx()))+0.1 ]) title('Balancer total Ax(blue) Bx(red)') hold off fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),Balancer.Force.Az(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balancer.Force.Bz(i,:),'r') grid on; title(sprintf('Balancer %d Az(blue) Bz(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balancer.Force.Bz()))-0.1 max(max(Balancer.Force.Bz()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),Balancer.Force.Az(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),Balancer.Force.Bz(i+1,:),'r') grid on; title('Balancer total Az(blue) Bz(red)') axis([0 max(Crankshaft.Angle())*180./pi() min(min(Balancer.Force.Bz()))-0.1 max(max(Balancer.Force.Bz()))+0.1 ]) hold off %All forces fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),All.Force.Ax(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),All.Force.Bx(i,:),'r') grid on; title(sprintf('All %d Ax(blue) Bx(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(All.Force.Bx()))-0.1 max(max(All.Force.Bx()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),All.Force.Ax(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),All.Force.Bx(i+1,:),'r') grid on; axis([0 max(Crankshaft.Angle())*180./pi() min(min(All.Force.Bx()))-0.1 max(max(All.Force.Bx()))+0.1 ]) title('All total Ax(blue) Bx(red)') hold off fig=fig+1; figure(fig) for i = 1 : Cylinder.Number subplot(2,round(0.5 + Cylinder.Number./2),i) plot(Crankshaft.Angle*180/pi(),All.Force.Az(i,:)) hold on; plot(Crankshaft.Angle*180/pi(),All.Force.Bz(i,:),'r') grid on; title(sprintf('All %d Az(blue) Bz(red)',i)) axis([0 max(Crankshaft.Angle())*180./pi() min(min(All.Force.Bz()))-0.1 max(max(All.Force.Bz()))+0.1 ]) hold off; end subplot(2,round(0.5 + Cylinder.Number./2),i+1) plot(Crankshaft.Angle*180/pi(),All.Force.Az(i+1,:)) hold on; plot(Crankshaft.Angle*180/pi(),All.Force.Bz(i+1,:),'r') grid on; title('All total Az(blue) Bz(red)') axis([0 max(Crankshaft.Angle())*180./pi() min(min(All.Force.Bz()))-0.1 max(max(All.Force.Bz()))+0.1 ]) hold off MaxForceAx = max(All.Force.Ax(i+1,:)) MinForceAx = min(All.Force.Ax(i+1,:)) MaxForceAz = max(All.Force.Az(i+1,:)) MinForceAz = min(All.Force.Az(i+1,:)) MaxForceBx = max(All.Force.Bx(i+1,:)) MinForceBx = min(All.Force.Bx(i+1,:)) MaxForceBz = max(All.Force.Bz(i+1,:)) MinForceBz = min(All.Force.Bz(i+1,:))