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