maglev.cc

Go to the documentation of this file.

#include "maglev.h"
#include <math.h>


maglev::maglev(dsim_model_setup *setup) : dsim_model(setup)
{
}

maglev::~maglev()
{
}
        
void maglev::initialize_data()
{               
    // Parameters
    
    const double pi = 3.14159265358979;
    
    double lMagnet = 0.01;                      // Magnet length
    double aMagnet = pi*pow(0.005,2);           // Magnet area
    double mu      = 5000e-6;                   // Transformer steel
    double mu0     = 4e-7*pi;                   // Permeability of free space
    double n       = 100;                       // Wire turns
    
    double a        = sqrt(2.0)*mu0*lMagnet/mu; // Scale factor
    double l0       = mu0*pow(n,2)*2*aMagnet/a; // Inductance
    double m        = 0.01;                     // Mass
    double g        = 9.806;                    // Acceleration of gravity
    double i        = 0;                        // Current
    
    double position = 0;
    double velocity = 0;
    
 
    i_dsim          = create_input("i",                     sd_type_double, &i,         1,1,    "A",        "Current");                 i_dsim.mark_telemetry();
    a_dsim          = create_parameter("a",                 sd_type_double, &a,         1,1,    "Nm/A^2",   "Scale");                   a_dsim.mark_telemetry();
    g_dsim          = create_parameter("g",                 sd_type_double, &g,         1,1,    "m/s^2",    "Acceleration of gravity"); g_dsim.mark_telemetry();
    l0_dsim         = create_parameter("l0",                sd_type_double, &l0,        1,1,    "H",        "Inductance");              l0_dsim.mark_telemetry();
    m_dsim          = create_parameter("m",                 sd_type_double, &m,         1,1,    "kg",       "Mass");                    m_dsim.mark_telemetry();
    position_dsim   = create_integrated_state("position",       sd_type_double, &position,  1,1,    "m",        "Position");                position_dsim.mark_telemetry();
    velocity_dsim   = create_integrated_state("velocity",       sd_type_double, &velocity,  1,1,    "m/s",      "Velocity");                velocity_dsim.mark_telemetry();
 
    configure_timestep(false,true,false);
}

void maglev::initialization_complete()
{
}

void maglev::initialize_timestep()
{
}

void maglev::rhs(double t,double jd)
{
    double a        = a_dsim.value_as_double();
    double l0       = l0_dsim.value_as_double();
    double m        = m_dsim.value_as_double();
    double g        = g_dsim.value_as_double();
    double i        = i_dsim.value_as_double();
    double position = position_dsim.value_as_double();
    double velocity = velocity_dsim.value_as_double();
    
    double f        = -(0.5/a)*l0*pow(i,2)/pow(1 + position/a,2);
    
    double accel    = f/m + g;
   
    position_dsim.set_derivative(velocity);
    velocity_dsim.set_derivative(accel);
}

void maglev::complete_timestep()
{
        dsim_model::complete_timestep();
}

extern "C"
{
    dsim_model *Maglev(void *setup);

        dsim_model *Maglev(void *setup)
        {
                return new maglev((dsim_model_setup *)setup);
        }
        
}