/* Object initialization is called once after all object have been created */
int dc_dc_converter::init(OBJECT *parent)
{
//initialize variables that are used internally
//set_terminal_voltage = 240; //V
//max_current_step_size = 100; //A
Rated_kW = 1000; //< nominal power in kW
Max_P = 1000;//< maximum real power capacity in kW
Min_P = 0;//< minimum real power capacity in kW
//Max_Q = 1000;//< maximum reactive power capacity in kVar
//Min_Q = 1000;//< minimum reactive power capacity in kVar
Rated_kVA = 2500; //< nominal capacity in kVA
Rated_kV = 100; //< nominal line-line voltage in kV
Rinternal = 0.035;
Rload = 1;
Rtotal = 0.05;
//Rphase = 0.03;
Rground = 0.03;
Rground_storage = 0.05;
Cinternal = 0;
Cground = 0;
Ctotal = 0;
Linternal = 0;
Lground = 0;
Ltotal = 0;
filter_120HZ = false;
filter_180HZ = false;
filter_240HZ = false;
pf_in = 0;
pf_out = 0;
number_of_phases_in = 0;
number_of_phases_out = 0;
phaseAIn = false;
phaseBIn = false;
phaseCIn = false;
phaseAOut = false;
phaseBOut = false;
phaseCOut = false;
switch_type_choice = IDEAL_SWITCH;
//generator_mode_choice = CONSTANT_PQ;
gen_status_v = ONLINE;
filter_type_v = BAND_PASS;
filter_imp_v = IDEAL_FILTER;
dc_dc_converter_type_v = BUCK_BOOST;
power_in = DC;
power_out = DC;
P_Out = 500; // P_Out and Q_Out are set by the user as set values to output in CONSTANT_PQ mode
Q_Out = 0;
margin = 50;
V_Set = 480;
I_out_prev = 0;
I_step_max = 50;
internal_losses = 0;
C_Storage_In = 0;
C_Storage_Out = 0;
losses = 0;
efficiency = 0;
//duty_ratio = 0;
service_ratio = 2; //if greater than 1, output voltage > input voltage
//if less than 1, output voltage < input voltage.
//input_frequency = 2000;
//frequency_losses = 0;
/* TODO: set the context-dependent initial value of properties */
gl_verbose("dc_dc_converter init: initialized the variables");
struct {
complex **var;
char *varname;
} map[] = {
// local object name, meter object name
{&pCircuit_V, "V_In"},
{&pLine_I, "I_In"},
};
static complex default_line_voltage[1], default_line_current[1];
int i;
// find parent meter, if not defined, use a default meter (using static variable 'default_meter')
if (parent!=NULL && strcmp(parent->oclass->name,"meter")==0)
{
parent_string = "meter";
for (i=0; i<sizeof(map)/sizeof(map[0]); i++)
*(map[i].var) = get_complex(parent,map[i].varname);
gl_verbose("dc_dc_converter init: mapped METER objects to internal variables");
}
else if (parent!=NULL && ((strcmp(parent->oclass->name,"inverter")==0))){
gl_verbose("dc_dc_converter init: parent WAS found, is an inverter!");
parent_string = "inverter";
// construct circuit variable map to meter
/// @todo use triplex property mapping instead of assuming memory order for meter variables (residential, low priority) (ticket #139)
for (i=0; i<sizeof(map)/sizeof(map[0]); i++){
*(map[i].var) = get_complex(parent,map[i].varname);
}
gl_verbose("dc_dc_converter init: mapped INVERTER objects to local variables");
}
else if (parent!=NULL && strcmp(parent->oclass->name,"battery")==0){
gl_verbose("dc_dc_converter init: parent WAS found, is an battery!");
parent_string = "dc_dc_converter";
// construct circuit variable map to meter
/// @todo use triplex property mapping instead of assuming memory order for meter variables (residential, low priority) (ticket #139)
for (i=0; i<sizeof(map)/sizeof(map[0]); i++){
*(map[i].var) = get_complex(parent,map[i].varname);
}
gl_verbose("dc_dc_converter init: mapped BATTERY objects to internal variables");
}
else{
// construct circuit variable map to meter
/// @todo use triplex property mapping instead of assuming memory order for meter variables (residential, low priority) (ticket #139)
gl_verbose("dc_dc_converter init: mapped meter objects to internal variables");
OBJECT *obj = OBJECTHDR(this);
gl_verbose("dc_dc_converter init: no parent meter defined, parent is not a meter");
gl_warning("dc_dc_converter:%d %s", obj->id, parent==NULL?"has no parent meter defined":"parent is not a meter");
// attach meter variables to each circuit in the default_meter
*(map[0].var) = &default_line_voltage[0];
*(map[1].var) = &default_line_current[0];
// provide initial values for voltages
default_line_voltage[0] = complex(0,0);
default_line_current[0] = complex(0,0);
//default_line123_voltage[1] = complex(V_Max/sqrt(3.0)*cos(2*PI/3),V_Max/sqrt(3.0)*sin(2*PI/3));
//default_line123_voltage[2] = complex(V_Max/sqrt(3.0)*cos(-2*PI/3),V_Max/sqrt(3.0)*sin(-2*PI/3));
}
gl_verbose("dc_dc_converter init: finished connecting with meter");
if (gen_mode_v==UNKNOWN)
{
OBJECT *obj = OBJECTHDR(this);
throw("Generator control mode is not specified");
}
if (gen_status_v== dc_dc_converter::OFFLINE)
{
//OBJECT *obj = OBJECTHDR(this);
throw("Generator is out of service!");
}else
{
//initialize variables that are used internally
//need to check for parameters SWITCH_TYPE, FILTER_TYPE, FILTER_IMPLEMENTATION, GENERATOR_MODE
/*
if (Rated_kW!=0.0) SB = Rated_kW/sqrt(1-Rated_pf*Rated_pf);
if (Rated_kVA!=0.0) SB = Rated_kVA/3;
if (Rated_kV!=0.0) EB = Rated_kV/sqrt(3.0);
if (SB!=0.0) ZB = EB*EB/SB;
else throw("Generator power capacity not specified!");
double Real_Rinternal = Rinternal * ZB;
double Real_Rload = Rload * ZB;
double Real_Rtotal = Rtotal * ZB;
double Real_Rphase = Rphase * ZB;
double Real_Rground = Rground * ZB;
double Real_Rground_storage = Rground_storage * ZB;
double[3] Real_Rfilter = Rfilter * ZB;
double Real_Cinternal = Cinternal * ZB;
double Real_Cground = Cground * ZB;
double Real_Ctotal = Ctotal * ZB;
double[3] Real_Cfilter = Cfilter * ZB;
double Real_Linternal = Linternal * ZB;
double Real_Lground = Lground * ZB;
double Real_Ltotal = Ltotal * ZB;
double[3] Real_Lfilter = Lfilter * ZB;
tst = complex(Real_Rground,Real_Lground);
AMx[0][0] = complex(Real_Rinternal,Real_Linternal) + tst;
AMx[1][1] = complex(Real_Rinternal,Real_Linternal) + tst;
AMx[2][2] = complex(Real_Rinternal,Real_Linternal) + tst;
// AMx[0][0] = AMx[1][1] = AMx[2][2] = complex(Real_Rs+Real_Rg,Real_Xs+Real_Xg);
AMx[0][1] = AMx[0][2] = AMx[1][0] = AMx[1][2] = AMx[2][0] = AMx[2][1] = tst;
*/
//all other variables set in input file through public parameters
switch(dc_dc_converter_type_v){
case BUCK:
efficiency = 0.9;
break;
case BOOST:
efficiency = 0.9;
break;
case BUCK_BOOST:
efficiency = 0.8;
break;
default:
efficiency = 0.8;
break;
}
internal_switch_resistance(switch_type_choice);
filter_circuit_impact(filter_type_v, filter_imp_v);
}
return 1;
}
/* Object initialization is called once after all object have been created */
int inverter::init(OBJECT *parent)
{
OBJECT *obj = OBJECTHDR(this);
// construct circuit variable map to meter
static complex default_line123_voltage[3], default_line1_current[3];
int i;
// find parent meter or triplex_meter, if not defined, use default voltages, and if
// the parent is not a meter throw an exception
if (parent!=NULL && gl_object_isa(parent,"meter"))
{
// attach meter variables to each circuit
parent_string = "meter";
struct {
complex **var;
char *varname;
}
map[] = {
// local object name, meter object name
{&pCircuit_V, "voltage_A"}, // assumes 2 and 3 follow immediately in memory
{&pLine_I, "current_A"}, // assumes 2 and 3(N) follow immediately in memory
};
/// @todo use triplex property mapping instead of assuming memory order for meter variables (residential, low priority) (ticket #139)
NR_mode = get_bool(parent,"NR_mode");
for (i=0; i<sizeof(map)/sizeof(map[0]); i++)
*(map[i].var) = get_complex(parent,map[i].varname);
node *par = OBJECTDATA(parent, node);
number_of_phases_out = 0;
if (par->has_phase(PHASE_A))
number_of_phases_out += 1;
if (par->has_phase(PHASE_B))
number_of_phases_out += 1;
if (par->has_phase(PHASE_C))
number_of_phases_out += 1;
}
else if (parent!=NULL && gl_object_isa(parent,"triplex_meter"))
{
parent_string = "triplex_meter";
number_of_phases_out = 4; //Indicates it is a triplex node and should be handled differently
struct {
complex **var;
char *varname;
}
map[] = {
// local object name, meter object name
{&pCircuit_V, "voltage_12"}, // assumes 1N and 2N follow immediately in memory
{&pLine_I, "current_1"}, // assumes 2 and 3(N) follow immediately in memory
{&pLine12, "current_12"}, // maps current load 1-2 onto triplex load
/// @todo use triplex property mapping instead of assuming memory order for meter variables (residential, low priority) (ticket #139)
};
NR_mode = get_bool(parent,"NR_mode");
// attach meter variables to each circuit
for (i=0; i<sizeof(map)/sizeof(map[0]); i++)
{
if ((*(map[i].var) = get_complex(parent,map[i].varname))==NULL)
{
GL_THROW("%s (%s:%d) does not implement triplex_meter variable %s for %s (house:%d)",
/* TROUBLESHOOT
The Inverter requires that the triplex_meter contains certain published properties in order to properly connect
the inverter to the triplex-meter. If the triplex_meter does not contain those properties, GridLAB-D may
suffer fatal pointer errors. If you encounter this error, please report it to the developers, along with
the version of GridLAB-D that raised this error.
*/
parent->name?parent->name:"unnamed object", parent->oclass->name, parent->id, map[i].varname, obj->name?obj->name:"unnamed", obj->id);
}
}
}
else if ((parent != NULL && strcmp(parent->oclass->name,"meter") != 0)||(parent != NULL && strcmp(parent->oclass->name,"triplex_meter") != 0))
{
throw("Inverter must have a meter or triplex meter as it's parent");
/* TROUBLESHOOT
Check the parent object of the inverter. The inverter is only able to be childed via a meter or
triplex meter when connecting into powerflow systems. You can also choose to have no parent, in which
case the inverter will be a stand-alone application using default voltage values for solving purposes.
*/
}
else
{
parent_string = "none";
number_of_phases_out = 3;
struct {
complex **var;
char *varname;
}
map[] = {
// local object name, meter object name
{&pCircuit_V, "voltage_A"}, // assumes 2 and 3 follow immediately in memory
{&pLine_I, "current_A"}, // assumes 2 and 3(N) follow immediately in memory
};
gl_warning("Inverter:%d has no parent meter object defined; using static voltages", obj->id);
NR_mode = false;
// attach meter variables to each circuit in the default_meter
*(map[0].var) = &default_line123_voltage[0];
*(map[1].var) = &default_line1_current[0];
// provide initial values for voltages
default_line123_voltage[0] = complex(Rated_kV*1000/sqrt(3.0),0);
default_line123_voltage[1] = complex(Rated_kV*1000/sqrt(3.0)*cos(2*PI/3),Rated_kV*1000/sqrt(3.0)*sin(2*PI/3));
default_line123_voltage[2] = complex(Rated_kV*1000/sqrt(3.0)*cos(-2*PI/3),Rated_kV*1000/sqrt(3.0)*sin(-2*PI/3));
}
if (gen_mode_v == UNKNOWN)
{
gl_warning("Inverter control mode is not specified! Using default: SUPPLY_DRIVEN");
gen_mode_v = SUPPLY_DRIVEN;
}
if (gen_status_v == UNKNOWN)
{
gl_warning("Inverter status is unknown! Using default: ONLINE");
gen_status_v = ONLINE;
}
if (inverter_type_v == UNKNOWN)
{
gl_warning("Inverter type is unknown! Using default: PWM");
inverter_type_v = PWM;
}
//need to check for parameters SWITCH_TYPE, FILTER_TYPE, FILTER_IMPLEMENTATION, GENERATOR_MODE
/*
if (Rated_kW!=0.0) SB = Rated_kW/sqrt(1-Rated_pf*Rated_pf);
if (Rated_kVA!=0.0) SB = Rated_kVA/3;
if (Rated_kV!=0.0) EB = Rated_kV/sqrt(3.0);
if (SB!=0.0) ZB = EB*EB/SB;
else throw("Generator power capacity not specified!");
double Real_Rinternal = Rinternal * ZB;
double Real_Rload = Rload * ZB;
double Real_Rtotal = Rtotal * ZB;
double Real_Rphase = Rphase * ZB;
double Real_Rground = Rground * ZB;
double Real_Rground_storage = Rground_storage * ZB;
double[3] Real_Rfilter = Rfilter * ZB;
double Real_Cinternal = Cinternal * ZB;
double Real_Cground = Cground * ZB;
double Real_Ctotal = Ctotal * ZB;
double[3] Real_Cfilter = Cfilter * ZB;
double Real_Linternal = Linternal * ZB;
double Real_Lground = Lground * ZB;
double Real_Ltotal = Ltotal * ZB;
double[3] Real_Lfilter = Lfilter * ZB;
tst = complex(Real_Rground,Real_Lground);
AMx[0][0] = complex(Real_Rinternal,Real_Linternal) + tst;
AMx[1][1] = complex(Real_Rinternal,Real_Linternal) + tst;
AMx[2][2] = complex(Real_Rinternal,Real_Linternal) + tst;
// AMx[0][0] = AMx[1][1] = AMx[2][2] = complex(Real_Rs+Real_Rg,Real_Xs+Real_Xg);
AMx[0][1] = AMx[0][2] = AMx[1][0] = AMx[1][2] = AMx[2][0] = AMx[2][1] = tst;
*/
//all other variables set in input file through public parameters
switch(inverter_type_v)
{
case TWO_PULSE:
efficiency = 0.8;
break;
case SIX_PULSE:
efficiency = 0.8;
break;
case TWELVE_PULSE:
efficiency = 0.8;
break;
case PWM:
efficiency = 0.9;
break;
default:
efficiency = 0.8;
break;
}
internal_switch_resistance(switch_type_choice);
filter_circuit_impact(filter_type_v, filter_imp_v);
return 1;
}
/* Object initialization is called once after all object have been created */
int rectifier::init(OBJECT *parent)
{
OBJECT *obj = OBJECTHDR(this);
//initialize variables that are used internally
//set_terminal_voltage = 240; //V
//max_current_step_size = 100; //A
Rated_kW = 1000; //< nominal power in kW
Max_P = 1000;//< maximum real power capacity in kW
Min_P = 0;//< minimum real power capacity in kW
//Max_Q = 1000;//< maximum reactive power capacity in kVar
//Min_Q = 1000;//< minimum reactive power capacity in kVar
Rated_kVA = 1500; //< nominal capacity in kVA
Rated_kV = 10; //< nominal line-line voltage in kV
Rinternal = 0.035;
Rload = 1;
Rtotal = 0.05;
//XphaseA = complex(1 * 1,0);
//XphaseB = complex(1 * -0.5, 1 * 0.866);
//XphaseC = complex(1 * -0.5, 1 * -0.866);
XphaseA = complex(5 * 1,0);
XphaseB = complex(5 * 1,0);
XphaseC = complex(5 * 1,0);
Rground = 0.03;
Rground_storage = 0.05;
Vdc = 480;
Cinternal = 0;
Cground = 0;
Ctotal = 0;
Linternal = 0;
Lground = 0;
Ltotal = 0;
filter_120HZ = false;
filter_180HZ = false;
filter_240HZ = false;
pf_in = 1;
pf_out = 0;
number_of_phases_in = 3;
number_of_phases_out = 0;
phaseAIn = true;
phaseBIn = true;
phaseCIn = true;
phaseAOut = false;
phaseBOut = false;
phaseCOut = false;
V_In_Set_A = complex(360,0);
V_In_Set_B = complex(-180, 311.769);
V_In_Set_C = complex(-180,-311.769);
switch_type_choice = IDEAL_SWITCH;
//generator_mode_choice = CONSTANT_PQ;
gen_status_v = ONLINE;
filter_type_v = BAND_PASS;
filter_imp_v = IDEAL_FILTER;
rectifier_type_v = SIX_PULSE;
power_in = AC;
power_out = DC;
P_Out = 500; // P_Out and Q_Out are set by the user as set values to output in CONSTANT_PQ mode
Q_Out = 0;
margin = 50;
V_Rated = 360;
I_out_prev = 0;
I_step_max = 100;
internal_losses = 0;
C_Storage_Out = 0;
efficiency = 0;
losses = 0;
on_ratio = 0;
input_frequency = 2000;
frequency_losses = 0;
gl_verbose("rectifier init: initialized the variables");
int i;
if (parent!=NULL && gl_object_isa(parent,"meter"))
{
// TO DO: Figure out how to connect a DC device to meter.
}
if (parent!=NULL && gl_object_isa(parent,"inverter"))
{
parent_string = "inverter";
struct {
complex **var;
char *varname;
}
map[] = {
// local object name, meter object name
{&pCircuit_V, "Vdc"},
{&pLine_I, "I_In"}
};
// attach meter variables to each circuit
for (i=0; i<sizeof(map)/sizeof(map[0]); i++)
{
if ((*(map[i].var) = get_complex(parent,map[i].varname))==NULL)
{
GL_THROW("%s (%s:%d) does not implement triplex_meter variable %s for %s (inverter:%d)",
/* TROUBLESHOOT
The rectifier requires that the inverter contains certain published properties in order to properly connect. If you encounter this error, please report it to the developers, along with
the version of GridLAB-D that raised this error.
*/
parent->name?parent->name:"unnamed object", parent->oclass->name, parent->id, map[i].varname, obj->name?obj->name:"unnamed", obj->id);
}
}
}
if (parent_string == NULL){
GL_THROW("Rectifier's parent is not an inverter. The rectifier object's parent must be an inverter.");
}
/* TODO: set the context-dependent initial value of properties */
if (gen_mode_v==UNKNOWN)
{
GL_THROW("Generator control mode is not specified.");
}
else if(gen_mode_v == CONSTANT_V)
{
GL_THROW("Generator mode CONSTANT_V is not implemented yet.");
}
else if(gen_mode_v == CONSTANT_PQ)
{
GL_THROW("Generator mode CONSTANT_PQ is not implemented yet.");
}
else if(gen_mode_v == CONSTANT_PF)
{
GL_THROW("Generator mode CONSTANT_PF is not implemented yet.");
}
//need to check for parameters SWITCH_TYPE, FILTER_TYPE, FILTER_IMPLEMENTATION, GENERATOR_MODE
/*
if (Rated_kW!=0.0) SB = Rated_kW/sqrt(1-Rated_pf*Rated_pf);
if (Rated_kVA!=0.0) SB = Rated_kVA/3;
if (Rated_kV!=0.0) EB = Rated_kV/sqrt(3.0);
if (SB!=0.0) ZB = EB*EB/SB;
else throw("Generator power capacity not specified!");
double Real_Rinternal = Rinternal * ZB;
double Real_Rload = Rload * ZB;
double Real_Rtotal = Rtotal * ZB;
double Real_Xphase = Xphase * ZB;
double Real_Rground = Rground * ZB;
double Real_Rground_storage = Rground_storage * ZB;
double[3] Real_Rfilter = Rfilter * ZB;
double Real_Cinternal = Cinternal * ZB;
double Real_Cground = Cground * ZB;
double Real_Ctotal = Ctotal * ZB;
double[3] Real_Cfilter = Cfilter * ZB;
double Real_Linternal = Linternal * ZB;
double Real_Lground = Lground * ZB;
double Real_Ltotal = Ltotal * ZB;
double[3] Real_Lfilter = Lfilter * ZB;
tst = complex(Real_Rground,Real_Lground);
AMx[0][0] = complex(Real_Rinternal,Real_Linternal) + tst;
AMx[1][1] = complex(Real_Rinternal,Real_Linternal) + tst;
AMx[2][2] = complex(Real_Rinternal,Real_Linternal) + tst;
// AMx[0][0] = AMx[1][1] = AMx[2][2] = complex(Real_Rs+Real_Rg,Real_Xs+Real_Xg);
AMx[0][1] = AMx[0][2] = AMx[1][0] = AMx[1][2] = AMx[2][0] = AMx[2][1] = tst;
*/
//all other variables set in input file through public parameters
switch(rectifier_type_v){
case ONE_PULSE:
efficiency = 0.5;
break;
case TWO_PULSE:
efficiency = 0.7;
break;
case THREE_PULSE:
efficiency = 0.7;
break;
case SIX_PULSE:
efficiency = 0.8;
break;
case TWELVE_PULSE:
efficiency = 0.9;
break;
default:
efficiency = 0.8;
break;
}
internal_switch_resistance(switch_type_choice);
filter_circuit_impact((power_electronics::FILTER_TYPE)filter_type_v, (power_electronics::FILTER_IMPLEMENTATION)filter_imp_v);
gl_verbose("rectifier init: about to exit");
return 1;
}
