// Create a steam turbine system with high-, low-pressure turbine and propeller
// Number of branches: 	3
// Number of masses:	11
// Components: Two turbines and one propeller
// Engine type: Turbines


// Create inertias
dataModel.Add(M1 = new SteamTurbine() { X=551; Y=99 });
dataModel.Add(M2 = new Inertia() { X=629; Y=99 } );
dataModel.Add(M3 = new Inertia() { X=629; Y=204 } );
dataModel.Add(M4 = new Inertia() { X=683; Y=204 } );
dataModel.Add(M5 = new Inertia() { X=683; Y=311 } );
dataModel.Add(M6 = new Inertia() { X=683; Y=408 } );
dataModel.Add(M7 = new Inertia() { X=630; Y=408 } );
dataModel.Add(M8 = new Inertia() { X=630; Y=519 } );
dataModel.Add(M9 = new SteamTurbine() { X=551; Y=519 } );
dataModel.Add(M10 = new Inertia() { X=738; Y=311 } );
dataModel.Add(M11 = new Propeller() { X=793; Y=311 } );

// Add stiffness
dataModel.Add(S1 = new Stiffness(M1, M2));
dataModel.Add(S2 = new Mesh(M2, M3));
dataModel.Add(S3 = new Stiffness(M3, M4));
dataModel.Add(S4 = new Mesh(M4, M5));
dataModel.Add(S5 = new Mesh(M6, M5));
dataModel.Add(S6 = new Stiffness(M7, M6));
dataModel.Add(S7 = new Mesh(M8, M7));
dataModel.Add(S8 = new Stiffness(M9, M8));
dataModel.Add(S9 = new Stiffness(M5, M10));
dataModel.Add(S10 = new Stiffness(M10, M11));

// Set propeller as reference mass
M11.IsReferenceMass = true;

// Group into gear - LP 1st
Gear1 = new Gear();
Gear1.Add(M2);
Gear1.Add(M3);
dataModel.Add(Gear1);

// Group into gear - Wheel
Gear2 = new Gear();
Gear2.Add(M4);
Gear2.Add(M5);
Gear2.Add(M6);
dataModel.Add(Gear2);

// Group into gear - HP 1st
Gear3 = new Gear();
Gear3.Add(M7);
Gear3.Add(M8);
dataModel.Add(Gear3);

// Mass properties
M1.MomentOfInertia=992.50 kg*m^2;
M2.MomentOfInertia=24.01 kg*m^2;
M3.MomentOfInertia=4119 kg*m^2;
M4.MomentOfInertia=316.50 kg*m^2;
M5.MomentOfInertia=60330 kg*m^2;
M6.MomentOfInertia=290.10 kg*m^2;
M7.MomentOfInertia=5775 kg*m^2;
M8.MomentOfInertia=8.61 kg*m^2;
M9.MomentOfInertia=43.51 kg*m^2;
M10.MomentOfInertia=2412 kg*m^2;
M11.InertiaInAir=198000 kg*m^2;
M11.InertiaOfEntrainedWater=62970 kg*m^2;

// Stiffness properties
S1.TorsionalStiffness=6.779E+006 N*m/rad;
S2.TorsionalStiffness=1.000E+009 N*m/rad;
S3.TorsionalStiffness=4.291E+007 N*m/rad;
S4.TorsionalStiffness=1.000E+009 N*m/rad;
S5.TorsionalStiffness=1.000E+009 N*m/rad;
S6.TorsionalStiffness=9.444E+006 N*m/rad;
S7.TorsionalStiffness=1.000E+009 N*m/rad;
S8.TorsionalStiffness=2.091E+006 N*m/rad;
S9.TorsionalStiffness=8.023E+007 N*m/rad;
S10.TorsionalStiffness=2.372E+008 N*m/rad;

// Outer diameter on shafts
S1.DiamOuter = 0.0 m;
S2.DiamOuter = 0.0 m;
S3.DiamOuter = 0.0 m;
S4.DiamOuter = 0.0 m;
S5.DiamOuter = 0.0 m;
S6.DiamOuter = 0.0 m;
S7.DiamOuter = 0.0 m;
S8.DiamOuter = 0.0 m;
S9.DiamOuter = 0.605 m;
S10.DiamOuter = 0.610 m;

// Shaft damping
S1.DynamicMagnifier = 180;
S2.DynamicMagnifier = 100;
S3.DynamicMagnifier = 180;
S4.DynamicMagnifier = 100;
S5.DynamicMagnifier = 100;
S6.DynamicMagnifier = 180;
S7.DynamicMagnifier = 100;
S8.DynamicMagnifier = 180;
S9.DynamicMagnifier = 180;
S10.DynamicMagnifier = 180;

// Propeller properties
M11.PowerMCR = 27300 kW;
M11.BladeNumber = 4;
M11.Type = "Fixed pitch"
M11.DampingModel = "Schwaneke";
M11.EntrainedWaterModel = "Manual";
M11.DiameterBoss = 1400 mm;
M11.DiameterPropeller = 8800 mm;
M11.ExpandedAreaRatio = 0.59;
M11.PitchRatio = 0.84;
M11.FirstOrderTorqueAmplitude = 0.06;
M11.SecondOrderTorqueAmplitude = 0.02;

// MCR power and speed for turbines and propeller
M1.PowerMCR = 13760 kW;
M1.SpeedMCR = 3153 rpm;
M9.PowerMCR = 13540 kW;
M9.SpeedMCR = 4777 rpm;
M11.PowerMCR = 27300 kW;
M11.SpeedMCR = 83 rpm;

// Set correct rpm ratio for the 1st stage gear wheels
M3.RpmRatio = 6.375;
M7.RpmRatio = 6.375;

// Operating conditions and load cases
dataModel.AddOperatingMode(OPM1 = new OperatingMode() { Name="Full pitch" });
dataModel.AddLoadCase(OPM1,LC1 = new LoadCase() { Name="Load case 1"; MaxSpeed=90; MinSpeed=20; StepSpeed=1; Type="Propeller law - full pitch"});
