Eclipse SUMO - Simulation of Urban MObility
Loading...
Searching...
No Matches
MSCFModel_Rail.h
Go to the documentation of this file.
1/****************************************************************************/
2// Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
3// Copyright (C) 2012-2023 German Aerospace Center (DLR) and others.
4// This program and the accompanying materials are made available under the
5// terms of the Eclipse Public License 2.0 which is available at
6// https://www.eclipse.org/legal/epl-2.0/
7// This Source Code may also be made available under the following Secondary
8// Licenses when the conditions for such availability set forth in the Eclipse
9// Public License 2.0 are satisfied: GNU General Public License, version 2
10// or later which is available at
11// https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
12// SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
13/****************************************************************************/
18// <description missing>
19/****************************************************************************/
20#pragma once
21#include <config.h>
22
23
24#include "MSCFModel.h"
25
26
27
28class MSCFModel_Rail : public MSCFModel {
29
30public:
34 MSCFModel_Rail(const MSVehicleType* vtype);
35
36 double followSpeed(const MSVehicle* const veh, double speed, double gap2pred, double predSpeed,
37 double predMaxDecel, const MSVehicle* const pred = 0, const CalcReason usage = CalcReason::CURRENT) const;
38
39 virtual int getModelID() const;
40
41 virtual MSCFModel* duplicate(const MSVehicleType* vtype) const;
42
43
44
45 virtual ~MSCFModel_Rail();
46
47 virtual double maxNextSpeed(double speed, const MSVehicle* const veh) const;
48
49 virtual double minNextSpeed(double speed, const MSVehicle* const veh) const;
50
51 virtual double minNextSpeedEmergency(double speed, const MSVehicle* const veh = 0) const;
52
53 double getSpeedAfterMaxDecel(double v) const;
54
56
57 double finalizeSpeed(MSVehicle* const veh, double vPos) const;
58
59 double freeSpeed(const MSVehicle* const veh, double speed, double seen, double maxSpeed,
60 const bool onInsertion, const CalcReason usage = CalcReason::CURRENT) const;
61
62 bool startupDelayStopped() const {
63 // starup delay in trains is dominated by inertia + brake delay and thus applies to any kind of stopping
64 return true;
65 }
66
67private:
68
69
70 typedef std::map<double, double> LookUpMap;
71
72 struct TrainParams {
73 double weight;
74 double mf;
75 double length;
76 double decl;
77 double vmax;
78 double recovery;
79 double rotWeight;
82 };
83
84 double getInterpolatedValueFromLookUpMap(double speed, const LookUpMap* lookUpMap) const;
85
86
87public:
88 double stopSpeed(const MSVehicle* const veh, const double speed, double gap, double decel, const CalcReason usage = CalcReason::CURRENT) const;
89
90// class VehicleVariables : public MSCFModel::VehicleVariables {
91//
92// public:
93//
94// const std::string getTrainType() const { return myTrainType; };
95//
96// void setTrainType(std::string trainType) { myTrainType = trainType; }
97//
98// bool isNotYetInitialized() {
99// return notYetInitialized;
100// }
101//
102// void setInitialized() {
103// notYetInitialized = false;
104// }
105//
106// private:
107// std::string myTrainType;
108// bool notYetInitialized = true;
109//
110//
111// };
112
113private:
114
115
117
118 LookUpMap initNGT400Traction() const { // early version of NGT 400
119 LookUpMap map;
120 map[0] = 716.0;
121 map[10] = 700.0;
122 map[20] = 684.0;
123 map[30] = 669.0;
124 map[40] = 653.0;
125 map[50] = 637.0;
126 map[60] = 622.0;
127 map[70] = 606.0;
128 map[80] = 590.0;
129 map[90] = 574.0;
130 map[100] = 559.0;
131 map[110] = 543.0;
132 map[120] = 527.0;
133 map[130] = 507.0;
134 map[140] = 471.0;
135 map[150] = 439.0;
136 map[160] = 412.0;
137 map[170] = 388.0;
138 map[180] = 366.0;
139 map[190] = 347.0;
140 map[200] = 329.0;
141 map[210] = 314.0;
142 map[220] = 299.0;
143 map[230] = 286.0;
144 map[240] = 275.0;
145 map[250] = 264.0;
146 map[260] = 253.0;
147 map[270] = 244.0;
148 map[280] = 235.0;
149 map[290] = 227.0;
150 map[300] = 220.0;
151 map[310] = 213.0;
152 map[320] = 206.0;
153 map[330] = 200.0;
154 map[340] = 194.0;
155 map[350] = 188.0;
156 map[360] = 183.0;
157 map[370] = 178.0;
158 map[380] = 173.0;
159 map[390] = 169.0;
160 map[400] = 165.0;
161 map[410] = 160.6;
162 map[420] = 156.8;
163 map[430] = 153.2;
164 map[440] = 149.7;
165 map[450] = 146.4;
166 map[460] = 143.2;
167 map[470] = 140.1;
168 map[480] = 137.2;
169 return map;
170 }
171
172 LookUpMap initNGT400Resistance() const { // early version of NGT 400
173 LookUpMap map;
174 map[0] = 1.9;
175 map[10] = 2.1;
176 map[20] = 2.4;
177 map[30] = 2.9;
178 map[40] = 3.5;
179 map[50] = 4.2;
180 map[60] = 5.0;
181 map[70] = 6.0;
182 map[80] = 7.0;
183 map[90] = 8.2;
184 map[100] = 9.6;
185 map[110] = 11.0;
186 map[120] = 12.6;
187 map[130] = 14.3;
188 map[140] = 16.1;
189 map[150] = 18.0;
190 map[160] = 20.1;
191 map[170] = 22.3;
192 map[180] = 24.6;
193 map[190] = 27.0;
194 map[200] = 29.5;
195 map[210] = 32.2;
196 map[220] = 35.0;
197 map[230] = 37.9;
198 map[240] = 40.9;
199 map[250] = 44.1;
200 map[260] = 47.4;
201 map[270] = 50.8;
202 map[280] = 54.3;
203 map[290] = 58.0;
204 map[300] = 61.7;
205 map[310] = 65.6;
206 map[320] = 69.7;
207 map[330] = 73.8;
208 map[340] = 78.1;
209 map[350] = 82.4;
210 map[360] = 86.9;
211 map[370] = 91.6;
212 map[380] = 96.3;
213 map[390] = 101.2;
214 map[400] = 106.2;
215 map[410] = 111.3;
216 map[420] = 116.6;
217 map[430] = 121.9;
218 map[440] = 127.4;
219 map[450] = 133.0;
220 map[460] = 138.7;
221 map[470] = 144.6;
222 map[480] = 150.6;
223 return map;
224 }
225
227 TrainParams params;
228 params.weight = 384;
229 params.mf = 1.04;
230 params.length = 200;
231 params.decl = 0.9;
232 params.vmax = 500 / 3.6;
233 params.recovery = 0.6;
234 params.rotWeight = params.weight * params.mf;
235 params.traction = initNGT400Traction();
237 return params;
238 }
239
241 LookUpMap map;
242 map[0] = 274.5;
243 map[10] = 274.5;
244 map[20] = 274.5;
245 map[30] = 274.5;
246 map[40] = 274.5;
247 map[50] = 274.5;
248 map[60] = 274.5;
249 map[70] = 274.5;
250 map[80] = 274.5;
251 map[90] = 274.5;
252 map[100] = 274.5;
253 map[110] = 274.5;
254 map[120] = 274.5;
255 map[130] = 274.5;
256 map[140] = 274.5;
257 map[150] = 274.5;
258 map[160] = 274.5;
259 map[170] = 274.5;
260 map[180] = 274.5;
261 map[190] = 274.5;
262 map[200] = 274.5;
263 map[210] = 273;
264 map[220] = 262.3;
265 map[230] = 250;
266 map[240] = 239;
267 map[250] = 229;
268 map[260] = 222;
269 map[270] = 213;
270 map[280] = 205;
271 map[290] = 198;
272 map[300] = 191;
273 map[310] = 186;
274 map[320] = 180;
275 map[330] = 174;
276 map[340] = 169;
277 map[350] = 165;
278 map[360] = 160;
279 map[370] = 155.5;
280 map[380] = 151;
281 map[390] = 148;
282 map[400] = 144;
283 map[410] = 140;
284 map[420] = 136;
285 map[430] = 134;
286 map[440] = 131;
287 map[450] = 128;
288 map[460] = 125;
289 map[470] = 123;
290 return map;
291 }
292
294 LookUpMap map;
295 map[0] = 5.71;
296 map[10] = 6.01;
297 map[20] = 6.4;
298 map[30] = 7.0;
299 map[40] = 7.51;
300 map[50] = 8.34;
301 map[60] = 9.2;
302 map[70] = 10.4;
303 map[80] = 11.3;
304 map[90] = 12.58;
305 map[100] = 14.05;
306 map[110] = 15.5;
307 map[120] = 16.9;
308 map[130] = 18.75;
309 map[140] = 20.77;
310 map[150] = 22.9;
311 map[160] = 25.0;
312 map[170] = 26.9;
313 map[180] = 29.38;
314 map[190] = 32.0;
315 map[200] = 34.5;
316 map[210] = 37.0;
317 map[220] = 39.91;
318 map[230] = 43.11;
319 map[240] = 46.4;
320 map[250] = 49.3;
321 map[260] = 52.4;
322 map[270] = 56.09;
323 map[280] = 59.8;
324 map[290] = 63.7;
325 map[300] = 66.8;
326 map[310] = 71.0;
327 map[320] = 75.39;
328 map[330] = 79.9;
329 map[340] = 84.1;
330 map[350] = 87.9;
331 map[360] = 92.7;
332 map[370] = 97.6;
333 map[380] = 102.0;
334 map[390] = 107.0;
335 map[400] = 112.0;
336 map[410] = 117.5;
337 map[420] = 123.0;
338 map[430] = 128.3;
339 map[440] = 133.3;
340 map[450] = 139.2;
341 map[460] = 145.5;
342 map[470] = 150.0;
343 return map;
344 }
345
347 TrainParams params;
348 params.weight = 384;
349 params.mf = 1.04;
350 params.length = 200;
351 params.decl = 0.9;
352 params.vmax = 430 / 3.6;
353 params.recovery = 0.6;
354 params.rotWeight = params.weight * params.mf;
357 return params;
358 }
359
361 LookUpMap map;
362 map[0] = 400;
363 map[10] = 394;
364 map[20] = 388;
365 map[30] = 382;
366 map[40] = 377;
367 map[50] = 372;
368 map[60] = 369;
369 map[70] = 366;
370 map[80] = 363;
371 map[90] = 361;
372 map[100] = 349;
373 map[110] = 317;
374 map[120] = 290;
375 map[130] = 268;
376 map[140] = 249;
377 map[150] = 232;
378 map[160] = 217;
379 map[170] = 205;
380 map[180] = 193;
381 map[190] = 183;
382 map[200] = 174;
383 map[210] = 165;
384 map[220] = 158;
385 map[230] = 151;
386 map[240] = 145;
387 map[250] = 139;
388 return map;
389 }
390
392 LookUpMap map;
393 map[0] = 10.7;
394 map[10] = 12.3;
395 map[20] = 14.2;
396 map[30] = 16.4;
397 map[40] = 18.7;
398 map[50] = 21.3;
399 map[60] = 24.2;
400 map[70] = 27.3;
401 map[80] = 30.6;
402 map[90] = 34.1;
403 map[100] = 37.9;
404 map[110] = 41.9;
405 map[120] = 46.2;
406 map[130] = 50.6;
407 map[140] = 55.4;
408 map[150] = 60.4;
409 map[160] = 65.6;
410 map[170] = 71.1;
411 map[180] = 76.7;
412 map[190] = 82.6;
413 map[200] = 88.8;
414 map[210] = 95.2;
415 map[220] = 101.8;
416 map[230] = 108.7;
417 map[240] = 115.8;
418 map[250] = 123.1;
419 return map;
420 }
421
423 TrainParams params;
424 params.weight = 876;
425 params.mf = 1.1;
426 params.length = 358;
427 params.decl = 0.5;
428 params.vmax = 250 / 3.6;
429 params.recovery = 0.1;
430 params.rotWeight = params.weight * params.mf;
431 params.traction = initICE1Traction();
433 return params;
434 }
435
437 TrainParams params;
438 params.weight = 420;
439 params.mf = 1.04;
440 params.length = 200;
441 params.decl = 0.5;
442 params.vmax = 300 / 3.6;
443 params.recovery = 0.6;
444 params.rotWeight = params.weight * params.mf;
445 params.traction = initICE3Traction();
447 return params;
448 }
449
451 LookUpMap map;
452 map[0] = 300;
453 map[10] = 298;
454 map[20] = 297;
455 map[30] = 295;
456 map[40] = 293;
457 map[50] = 292;
458 map[60] = 290;
459 map[70] = 288;
460 map[80] = 286.5;
461 map[90] = 285;
462 map[100] = 283;
463 map[110] = 262;
464 map[120] = 240;
465 map[130] = 221;
466 map[140] = 206;
467 map[150] = 192;
468 map[160] = 180;
469 map[170] = 169;
470 map[180] = 160;
471 map[190] = 152;
472 map[200] = 144;
473 map[210] = 137;
474 map[220] = 131;
475 map[230] = 125;
476 map[240] = 120;
477 map[250] = 115;
478 map[260] = 111;
479 map[270] = 107;
480 map[280] = 103;
481 map[290] = 99;
482 map[300] = 96;
483 return map;
484 }
485
487 LookUpMap map;
488 map[0] = 7.4;
489 map[10] = 7.6;
490 map[20] = 8.0;
491 map[30] = 8.4;
492 map[40] = 9.1;
493 map[50] = 9.8;
494 map[60] = 10.7;
495 map[70] = 11.7;
496 map[80] = 12.8;
497 map[90] = 14.1;
498 map[100] = 15.5;
499 map[110] = 17.1;
500 map[120] = 18.8;
501 map[130] = 20.6;
502 map[140] = 22.6;
503 map[150] = 24.6;
504 map[160] = 26.9;
505 map[170] = 29.2;
506 map[180] = 31.7;
507 map[190] = 34.3;
508 map[200] = 37.1;
509 map[210] = 40.0;
510 map[220] = 43.1;
511 map[230] = 46.2;
512 map[240] = 49.6;
513 map[250] = 53.0;
514 map[260] = 56.6;
515 map[270] = 60.3;
516 map[280] = 64.1;
517 map[290] = 68.1;
518 map[300] = 71.8;
519 return map;
520 }
521
523 TrainParams params;
524 params.weight = 425.5;
525 params.mf = 1.1;
526 params.length = 207;
527 params.decl = 0.5;
528 params.vmax = 160 / 3.6;
529 params.recovery = 0.1;
530 params.rotWeight = params.weight * params.mf;
533 return params;
534 }
535
537 LookUpMap map;
538 map[0] = 300;
539 map[10] = 300;
540 map[20] = 300;
541 map[30] = 300;
542 map[40] = 300;
543 map[50] = 300;
544 map[60] = 300;
545 map[70] = 289;
546 map[80] = 253;
547 map[90] = 224;
548 map[100] = 202;
549 map[110] = 183;
550 map[120] = 168;
551 map[130] = 155;
552 map[140] = 144;
553 map[150] = 134;
554 map[160] = 125;
555 return map;
556 }
557
559 LookUpMap map;
560 map[0] = 8.5;
561 map[10] = 8.9;
562 map[20] = 9.5;
563 map[30] = 10.4;
564 map[40] = 11.4;
565 map[50] = 12.7;
566 map[60] = 14.1;
567 map[70] = 15.8;
568 map[80] = 17.7;
569 map[90] = 19.8;
570 map[100] = 22.6;
571 map[110] = 24.6;
572 map[120] = 27.3;
573 map[130] = 30.2;
574 map[140] = 33.3;
575 map[150] = 36.6;
576 map[160] = 40.2;
577 return map;
578 }
579
581 TrainParams params;
582 params.weight = 72.2;
583 params.mf = 1.04;
584 params.length = 46;
585 params.decl = 0.5;
586 params.vmax = 120 / 3.6;
587 params.recovery = 0;
588 params.rotWeight = params.weight * params.mf;
589 params.traction = initRB628Traction();
591 return params;
592 }
593
595 LookUpMap map;
596 map[0] = 60;
597 map[10] = 53.8;
598 map[20] = 47.6;
599 map[30] = 36.9;
600 map[40] = 28.7;
601 map[50] = 23.5;
602 map[60] = 20.0;
603 map[70] = 17.5;
604 map[80] = 15.2;
605 map[90] = 13.9;
606 map[100] = 12.8;
607 map[110] = 11.7;
608 map[120] = 10.8;
609 return map;
610 }
611
613 LookUpMap map;
614 map[0] = 1.29;
615 map[10] = 1.46;
616 map[20] = 1.73;
617 map[30] = 2.08;
618 map[40] = 2.52;
619 map[50] = 3.05;
620 map[60] = 3.66;
621 map[70] = 4.36;
622 map[80] = 5.16;
623 map[90] = 6.03;
624 map[100] = 7.00;
625 map[110] = 8.06;
626 map[120] = 9.2;
627 return map;
628 }
629
631 TrainParams params;
632 params.weight = 1440;
633 params.mf = 1.06;
634 params.length = 512;
635 params.decl = 0.3;
636 params.vmax = 120 / 3.6;
637 params.recovery = 0.05;
638 params.rotWeight = params.weight * params.mf;
639 params.traction = initFreightTraction();
641 return params;
642 }
643
645 LookUpMap map;
646 map[0] = 300;
647 map[10] = 296;
648 map[20] = 293;
649 map[30] = 289;
650 map[40] = 286;
651 map[50] = 282;
652 map[60] = 279;
653 map[70] = 275;
654 map[80] = 272;
655 map[90] = 255;
656 map[100] = 230;
657 map[110] = 209;
658 map[120] = 190;//guessed value
659 return map;
660 }
661
663 LookUpMap map;
664 map[0] = 1.9;
665 map[10] = 4.3;
666 map[20] = 8.5;
667 map[30] = 14.3;
668 map[40] = 21.7;
669 map[50] = 30.8;
670 map[60] = 41.5;
671 map[70] = 53.8;
672 map[80] = 67.8;
673 map[90] = 83.5;
674 map[100] = 110.7;
675 map[110] = 119.6;
676 map[120] = 140.2;
677 return map;
678 }
679
681 LookUpMap map;
682 map[0] = 150;
683 map[10] = 150;
684 map[20] = 150;
685 map[30] = 150;
686 map[40] = 150;
687 map[50] = 150;
688 map[60] = 140;
689 map[70] = 120;
690 map[80] = 105;
691 map[90] = 93;
692 map[100] = 84;
693 map[110] = 75;
694 map[120] = 70;
695 map[130] = 56;
696 map[140] = 52;
697 map[150] = 46;
698 map[160] = 40;
699 return map;
700 }
701
703 LookUpMap map;
704 map[0] = 2.6;
705 map[10] = 2.9;
706 map[20] = 3.3;
707 map[30] = 3.7;
708 map[40] = 4.3;
709 map[50] = 4.9;
710 map[60] = 5.7;
711 map[70] = 6.6;
712 map[80] = 7.5;
713 map[90] = 8.6;
714 map[100] = 9.7;
715 map[110] = 11.0;
716 map[120] = 12.3;
717 map[130] = 13.8;
718 map[140] = 15.3;
719 map[150] = 16.9;
720 map[160] = 18.7;
721 return map;
722 }
723
725 TrainParams params;
726 params.weight = 138;
727 params.mf = 1.04;
728 params.length = 67.5;
729 params.decl = 1.0;
730 params.vmax = 160 / 3.6;
731 params.recovery = 0.6;
732 params.rotWeight = params.weight * params.mf;
733 params.traction = initRB425Traction();
735 return params;
736 }
737
738// void initVehicleVariables(const MSVehicle *const pVehicle, MSCFModel_Rail::VehicleVariables *pVariables)const;
739
740};
741
742
TrainParams initICE3Params() const
TrainParams initNGT400_16Params() const
virtual ~MSCFModel_Rail()
LookUpMap initNGT400Traction() const
TrainParams initREDosto7Params() const
LookUpMap initRB628Traction() const
LookUpMap initRB628Resistance() const
MSCFModel::VehicleVariables * createVehicleVariables() const
Returns model specific values which are stored inside a vehicle and must be used with casting.
virtual double minNextSpeedEmergency(double speed, const MSVehicle *const veh=0) const
Returns the minimum speed after emergency braking, given the current speed (depends on the numerical ...
LookUpMap initICE3Traction() const
TrainParams initICE1Params() const
double freeSpeed(const MSVehicle *const veh, double speed, double seen, double maxSpeed, const bool onInsertion, const CalcReason usage=CalcReason::CURRENT) const
Computes the vehicle's safe speed without a leader.
TrainParams initRB628Params() const
double getInterpolatedValueFromLookUpMap(double speed, const LookUpMap *lookUpMap) const
virtual MSCFModel * duplicate(const MSVehicleType *vtype) const
Duplicates the car-following model.
LookUpMap initNGT400Resistance() const
LookUpMap initRB425Traction() const
LookUpMap initICE1Resistance() const
virtual int getModelID() const
Returns the model's ID; the XML-Tag number is used.
LookUpMap initICE1Traction() const
LookUpMap initNGT400_16Traction() const
double getSpeedAfterMaxDecel(double v) const
Returns the velocity after maximum deceleration.
TrainParams myTrainParams
virtual double minNextSpeed(double speed, const MSVehicle *const veh) const
Returns the minimum speed given the current speed (depends on the numerical update scheme and its ste...
LookUpMap initRB425Resistance() const
LookUpMap initREDosto7Resistance() const
double followSpeed(const MSVehicle *const veh, double speed, double gap2pred, double predSpeed, double predMaxDecel, const MSVehicle *const pred=0, const CalcReason usage=CalcReason::CURRENT) const
Computes the vehicle's follow speed (no dawdling)
LookUpMap initNGT400_16Resistance() const
TrainParams initRB425Params() const
TrainParams initFreightParams() const
LookUpMap initFreightResistance() const
LookUpMap initREDosto7Traction() const
virtual double maxNextSpeed(double speed, const MSVehicle *const veh) const
Returns the maximum speed given the current speed.
TrainParams initNGT400Params() const
LookUpMap initICE3Resistance() const
LookUpMap initFreightTraction() const
std::map< double, double > LookUpMap
bool startupDelayStopped() const
whether startupDelay should be applied after stopping
double finalizeSpeed(MSVehicle *const veh, double vPos) const
Applies interaction with stops and lane changing model influences. Called at most once per simulation...
double stopSpeed(const MSVehicle *const veh, const double speed, double gap, double decel, const CalcReason usage=CalcReason::CURRENT) const
Computes the vehicle's safe speed for approaching a non-moving obstacle (no dawdling)
The car-following model abstraction.
Definition MSCFModel.h:55
CalcReason
What the return value of stop/follow/free-Speed is used for.
Definition MSCFModel.h:77
@ CURRENT
the return value is used for calculating the next speed
Definition MSCFModel.h:79
Representation of a vehicle in the micro simulation.
Definition MSVehicle.h:77
The car-following model and parameter.