serial_example.cc

#include <iostream>
#include <stdexcept>
#include <cmath>
#include <limits>
#include "dynamic/dynamic.h"
using namespace std;
const double pi = 3.14159265;
const double g = 9.81; //ms-2
const double weight_1 = 5; //kg
const double weight_2 = 2; //kg
const double rigidity_1 = 100; //Nm-1
const double rigidity_2 = 15; //Nm-1
const double damping_1 = 10; //kgs-1
const double damping_2 = 1; //kgs-1
const double start_deflection_1 = 0.0; //m
const double start_deflection_2 = 0.0; //m
const double start_speed_1 = 0.0; //ms-1
const double start_speed_2 = 0.0; //ms-1

const double delta = 0.5; //s (this isn't computing delta time but only delta for output)
const double time_max = 50.0;//s
const double amplitude = 20; //N
const double frequency = 0.05; //Hz
const double shift_phase = 0; //rad
double Q1(dynamic::Object* p_object)
{
        //return amplitude*sin(2*pi*frequency*p_object->time()+shift_phase)+weight_1*g;
        return weight_1 * g;
}
double Q2(dynamic::Object* p_object)
{
        return weight_2 * g;
}
int main(int argC, char* argV[])
{
        //Output header of table
        cout << endl;
        cout << "Dynamic library by Jakub Krasa, 2006,2007" << endl;
        cout << "example dynamic five masses with start deflection and speed." << 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::mdo* p_mass_1 = new dynamic::mdo(weight_1,rigidity_1,damping_1,
                                start_deflection_1,start_speed_1);
                dynamic::mdo* p_mass_2 = new dynamic::mdo(weight_2,rigidity_2,damping_2,
                                start_deflection_2,start_speed_2);
                dynamic::Serial* p_system = new dynamic::Serial();
                dynamic::mdo* p_mass_3 = new dynamic::mdo(weight_1,rigidity_1,damping_1,
                                start_deflection_1,start_speed_1);
                dynamic::mdo* p_mass_4 = new dynamic::mdo(weight_2,rigidity_2,damping_2,
                                start_deflection_2,start_speed_2);
                dynamic::mdo* p_mass_5 = new dynamic::mdo(weight_2,rigidity_2,damping_2,
                                start_deflection_2,start_speed_2);
                dynamic::Serial* p_sub_system = new dynamic::Serial();
                p_mass_1->name.set("mass 1");
                p_mass_2->name.set("mass 2");
                p_mass_3->name.set("mass 3");
                p_mass_4->name.set("mass 4");
                p_mass_5->name.set("mass 5");
                p_system->name.set("system");
                p_sub_system->name.set("sub_system");
                p_sub_system->push_front(p_mass_3);
                p_sub_system->push_front(p_mass_4);
                p_system->push_front(p_mass_5);
                p_system->push_front(p_mass_2);
                p_system->push_front(p_sub_system);
                p_system->push_front(p_mass_1);
                p_mass_1->load(Q1);
                p_mass_2->load(Q2);
                p_mass_3->load(Q2);
                p_mass_4->load(Q2);
                p_mass_5->load(Q2);
                double time = 0;
                do
                {
                        time += delta;
                        p_system->time_to(time);
                        cout << p_system->time.get() << "\t\t"  
                                << p_mass_1->defle() << "\t" << p_mass_1->speed() << "\t" << p_mass_1->accel()
                                << "\t\t"       
                                << p_mass_2->defle() << "\t" << p_mass_2->speed() << "\t" << p_mass_2->accel()
                                << endl;
                }
                while((((time + delta/2)-time_max)/time_max) <= numeric_limits<double>::epsilon());
                delete p_system;
                delete p_mass_1;
                delete p_mass_2;
                delete p_mass_3;
                delete p_mass_4;
                delete p_mass_5;
                delete p_sub_system;
        }
        catch (exception& e)
        {
                cout << "sds exception: " << e.what() << endl;
                return EXIT_FAILURE;
        }
        return EXIT_SUCCESS;
}

---