Constructor.
The constructor reads in the defining parameters from a configuration file.
52 : TankNumber(tank_number)
53{
54 string token, strFuelName;
55 Element* element;
56 Element* element_Grain;
57 auto PropertyManager = exec->GetPropertyManager();
58 Area = 1.0;
59 Density = 6.6;
60 InitialTemperature = Temperature = -9999.0;
61 Ixx = Iyy = Izz = 0.0;
62 InertiaFactor = 1.0;
63 Radius = Contents = Standpipe = Length = InnerRadius = 0.0;
64 ExternalFlow = 0.0;
65 InitialStandpipe = 0.0;
66 Capacity = 0.00001; UnusableVol = 0.0;
67 Priority = InitialPriority = 1;
68 vXYZ.InitMatrix();
69 vXYZ_drain.InitMatrix();
70 ixx_unit = iyy_unit = izz_unit = 1.0;
71 grainType = gtUNKNOWN;
72
73 string type = el->GetAttributeValue("type");
74 if (type == "FUEL") Type = ttFUEL;
75 else if (type == "OXIDIZER") Type = ttOXIDIZER;
76 else Type = ttUNKNOWN;
77
78 Name = el->GetAttributeValue("name");
79
80 element = el->FindElement("location");
81 if (element) vXYZ = element->FindElementTripletConvertTo("IN");
82 else cerr << el->ReadFrom() << "No location found for this tank."
83 << endl;
84
85 vXYZ_drain = vXYZ;
86
87 element = el->FindElement("drain_location");
88 if (element) {
89 vXYZ_drain = element->FindElementTripletConvertTo("IN");
90 }
91
92 if (el->FindElement("radius"))
93 Radius = el->FindElementValueAsNumberConvertTo("radius", "IN");
94 if (el->FindElement("inertia_factor"))
95 InertiaFactor = el->FindElementValueAsNumber("inertia_factor");
96 if (el->FindElement("capacity"))
97 Capacity = el->FindElementValueAsNumberConvertTo("capacity", "LBS");
98 if (el->FindElement("contents"))
99 InitialContents = Contents = el->FindElementValueAsNumberConvertTo("contents", "LBS");
100 if (el->FindElement("unusable-volume"))
101 UnusableVol = el->FindElementValueAsNumberConvertTo("unusable-volume", "GAL");
102 if (el->FindElement("temperature"))
103 InitialTemperature = Temperature = el->FindElementValueAsNumber("temperature");
104 if (el->FindElement("standpipe"))
105 InitialStandpipe = Standpipe = el->FindElementValueAsNumberConvertTo("standpipe", "LBS");
106 if (el->FindElement("priority"))
107 InitialPriority = Priority = (int)el->FindElementValueAsNumber("priority");
108 if (el->FindElement("density"))
109 Density = el->FindElementValueAsNumberConvertTo("density", "LBS/GAL");
110 if (el->FindElement("type"))
111 strFuelName = el->FindElementValue("type");
112
113
114 SetPriority( InitialPriority );
115
116 if (Capacity == 0) {
117 cerr << el->ReadFrom()
118 << "Tank capacity must not be zero. Reset to 0.00001 lbs!" << endl;
119 Capacity = 0.00001;
120 Contents = 0.0;
121 }
123 cerr << el->ReadFrom() << "Tank capacity (" << Capacity
124 <<
" lbs) is lower than the amount of unusable fuel (" <<
GetUnusable()
125 << " lbs) for tank " << tank_number
126 << "! Did you accidentally swap unusable and capacity?" << endl;
127 throw("tank definition error");
128 }
129 if (Contents > Capacity) {
130 cerr << el->ReadFrom() << "Tank content (" << Contents
131 << " lbs) is greater than tank capacity (" << Capacity
132 << " lbs) for tank " << tank_number
133 << "! Did you accidentally swap contents and capacity?" << endl;
134 throw("tank definition error");
135 }
137 cerr << el->ReadFrom() << "Tank content (" << Contents
138 <<
" lbs) is lower than the amount of unusable fuel (" <<
GetUnusable()
139 << " lbs) for tank " << tank_number << endl;
140 }
141
142 PctFull = 100.0*Contents/Capacity;
143
144
145
146 element_Grain = el->FindElement("grain_config");
147 if (element_Grain) {
148
149 string strGType = element_Grain->GetAttributeValue("type");
150 if (strGType == "CYLINDRICAL") grainType = gtCYLINDRICAL;
151 else if (strGType == "ENDBURNING") grainType = gtENDBURNING;
152 else if (strGType == "FUNCTION") {
153 grainType = gtFUNCTION;
154 if (element_Grain->FindElement("ixx") != 0) {
155 Element* element_ixx = element_Grain->FindElement("ixx");
156 if (element_ixx->GetAttributeValue("unit") == "KG*M2") ixx_unit = 1.0/1.35594;
157 if (element_ixx->FindElement("function") != 0) {
158 function_ixx = new FGFunction(exec, element_ixx->FindElement("function"));
159 }
160 } else {
161 throw("For tank "+to_string(TankNumber)+" and when grain_config is specified an ixx must be specified when the FUNCTION grain type is specified.");
162 }
163
164 if (element_Grain->FindElement("iyy")) {
165 Element* element_iyy = element_Grain->FindElement("iyy");
166 if (element_iyy->GetAttributeValue("unit") == "KG*M2") iyy_unit = 1.0/1.35594;
167 if (element_iyy->FindElement("function") != 0) {
168 function_iyy = new FGFunction(exec, element_iyy->FindElement("function"));
169 }
170 } else {
171 throw("For tank "+to_string(TankNumber)+" and when grain_config is specified an iyy must be specified when the FUNCTION grain type is specified.");
172 }
173
174 if (element_Grain->FindElement("izz")) {
175 Element* element_izz = element_Grain->FindElement("izz");
176 if (element_izz->GetAttributeValue("unit") == "KG*M2") izz_unit = 1.0/1.35594;
177 if (element_izz->FindElement("function") != 0) {
178 function_izz = new FGFunction(exec, element_izz->FindElement("function"));
179 }
180 } else {
181 throw("For tank "+to_string(TankNumber)+" and when grain_config is specified an izz must be specified when the FUNCTION grain type is specified.");
182 }
183 }
184 else
185 cerr << el->ReadFrom() << "Unknown propellant grain type specified"
186 << endl;
187
188 if (element_Grain->FindElement("length"))
189 Length = element_Grain->FindElementValueAsNumberConvertTo("length", "IN");
190 if (element_Grain->FindElement("bore_diameter"))
191 InnerRadius = element_Grain->FindElementValueAsNumberConvertTo("bore_diameter", "IN")/2.0;
192
193
194
195 switch (grainType) {
196 case gtCYLINDRICAL:
197 if (Radius <= InnerRadius) {
198 const string s("The bore diameter should be smaller than the total grain diameter!");
199 cerr << element_Grain->ReadFrom() << endl << s << endl;
200 throw BaseException(s);
201 }
202 Volume = M_PI * Length * (Radius*Radius - InnerRadius*InnerRadius);
203 break;
204 case gtENDBURNING:
205 Volume = M_PI * Length * Radius * Radius;
206 break;
207 case gtFUNCTION:
208 Volume = 1;
209 break;
210 case gtUNKNOWN:
211 {
212 const string s("Unknown grain type found in this rocket engine definition.");
213 cerr << el->ReadFrom() << endl << s << endl;
214 throw BaseException(s);
215 }
216 }
217 Density = (Capacity*lbtoslug)/Volume;
218 }
219
220 CalculateInertias();
221
222 if (Temperature != -9999.0) InitialTemperature = Temperature =
FahrenheitToCelsius(Temperature);
223 Area = 40.0 * pow(Capacity/1975, 0.666666667);
224
225
227
228 bind(PropertyManager.get());
229
230 Debug(0);
231}
static constexpr double FahrenheitToCelsius(double fahrenheit)
Converts from degrees Fahrenheit to degrees Celsius.
double ProcessFuelName(const std::string &name)
Returns the density of a named fuel type.
double GetUnusable(void) const
Returns the amount of unusable fuel in the tank.
Inheritance diagram for FGTank:
1.9.8 
