electromotor_with_flywheel.cc

#include <iostream>
#include <stdexcept>
#include <cmath>
#include <cstdlib>
#include <limits>
#include "dynamic/dynamic.h"
using namespace std;
const double pi = 3.14159265;
const double g = 9.81; //ms-2
const double motor_inertia = 5; //
const double motor_torque_nominal = 1000; //Nm
const double motor_torque_start = 1400; //Nm
const double motor_torque_breakdown = 2100; //Nm
const double flywheel_inertia = 200; //
const double flywheel_rigidity = 10; //Nrad-1
const double flywheel_damping = 1; //Nrad-1s
const double flywheel_start_deflection = 0.0; //rad
const double flywheel_start_speed = 0.0; //raSystem-1
const double delta = 0.01;//s (this isn't computing delta time but only delta for output)
const double time_max = 5.0;//s
int main(int argC, char* argV[])
{
        //Output header of table
        cout << endl;
        cout << "Dynamic library by Jakub Krasa, 2006" << endl;
        cout << "example dynamic two rotary masses, first mass is electromotor." << endl << endl;
        cout << "\t\tMass 1 \t\t\t\tMass 2" << endl;
        cout << "Time\t\t" << "defl.\t" << "speed\t" << "accel.\t\t" << "defl.\t" << "speed\t" << "accel." << endl;
        cout.precision(3);
        cout << right << fixed;
        //computing
        try
        {
                dynamic::meacas* p_motor = new dynamic::meacas(motor_inertia,motor_torque_nominal,motor_torque_start,motor_torque_breakdown,6);
                dynamic::rdo* p_flywheel = new dynamic::rdo(flywheel_inertia,flywheel_rigidity,flywheel_damping,
                                flywheel_start_deflection,flywheel_start_speed);
                dynamic::Serial* p_system = new dynamic::Serial();
                p_system->push_back(p_motor);
                p_system->push_back(p_flywheel);
                double time = 0;
                do
                {
                        time += delta;
                        p_system->time_to(time);
                        cout << p_system->time.get() << "\t\t"  
                                << p_motor->defle() << "\t" << p_motor->speed() << "\t" << p_motor->accel()
                                << "\t\t"       
                                << p_flywheel->defle() << "\t" << p_flywheel->speed() << "\t" << p_flywheel->accel()
                                << endl;
                }
                while((((time + delta/2)-time_max)/time_max) <= numeric_limits<double>::epsilon());
                delete p_system;
                delete p_motor;
                delete p_flywheel;
        }
        catch (invalid_argument& e)
        {
                cout << "Invalid argument: " << e.what() << endl;
                //throw;
                return EXIT_FAILURE;
        }
        return EXIT_SUCCESS;
}

---