Hai bisogno di aiuto per modellare un sistema ferroviario con Simpy

Question

Ho bisogno di aiuto per modellare un sistema ferroviario (come un sistema di metropolitana) in Simpy, il problema è che il mio insegnante "simulazione" vuole che io uso Python + Simpy e non ho idea di come usarlo, dopo una settimana leggendo Simpy sono riuscito a capire l'esempio bancario, è facile, ma ora non ho ancora idea su come modellare il problema che l'insegnante mi ha dato ...

Il problema è grande ed esteso, ma ho solo bisogno di un piccolo aiuto, non all'intero progetto fatto, quindi, se qualcuno può aiutare, apprezzerò davvero.

Il problema ridotto che mi aiuterà a iniziare è così:

Immagina di avere 3 stazioni (a, b, c) separate da una certa distanza (diciamo 100 m) e ho 2 treni, uno su A e uno su C (treno1 che va da a a c, tran2 che va da c a) dove Ogni treno ha una velocità massima (diciamo 50m/s per entrambi) e accelerazione (diciamo 5m/s^2 e -5m/s^2 per la rottura) e tutto ciò che devono fare è fermarti su ogni stazione per alcuni Secondi (diciamo 24s) e continuano alla stazione successiva, quando raggiungono la fine, aspettano e vendali extra (cambiando i binari) e poi ricominciano da capo.

La stazione ha un limite di passeggeri (non devono essere simulati) l'unica cosa è che ogni stazione ha un numero casuale di passeggeri e quando il treno arriva ad alzarsi e altri scendono ...

Quindi, in pratica ho bisogno di un piccolo modello del treno e delle stazioni con quelle 3 stazioni ... se qualcuno può aiutarmi per favore ...

Il vero problema ha 22 stazioni, 2 stazioni extra per cambiare corsia, tempo casuale in attesa, 42 treni con diversa velocità e capacità massime, numero diverso di passeggeri generati a seconda della stazione e dell'ora del giorno, ecc ... Tutta questa roba i Riesci a gestire più tardi, ma la logica di modellazione del sistema di stazione del treno non riesco a capirlo ... grazie per il tuo aiuto!

Answer 1

Non sono affatto un esperto di Simpy, ma questo dovrebbe iniziare:

"""
Simulation of a train network
"""
from SimPy.Simulation import *
from math import sqrt
from random import randint
from itertools import cycle

def timeTo(A, maxV, d):
    """
    Given a trapezoidal velocity envelope defined by
    A       constant acceleration, m/s/s
    maxV    maximumum velocity, m/s
    return time in seconds required to travel
    d       distance, m
    """
    tA = float(maxV)/A          # time to accelerate to maxV
    dA = A*tA*tA                # distance traveled during acceleration from 0 to maxV and back to 0
    if (d < dA):                # train never reaches full speed?
        return sqrt(4.0*d/A)        # time needed to accelerate to half-way point then decelerate to destination
    else:
        return 2*tA + (d-dA)/maxV   # time to accelerate to maxV plus travel at maxV plus decelerate to destination

class Train(Process):
    def __init__(self, name, sim, accel=1.0, maxV=50.0, passengers=0, maxPassengers=400):
        super(Train, self).__init__(name, sim)
        self.accel = accel
        self.maxV  = maxV
        self.p     = passengers
        self.maxP  = maxPassengers

    def roll(self, route):
        here = route.next()     # starting location
        for dest in route:
            # travel to next station
            print "{:.1f}s: {} leaving {} for {}".format(self.sim.now(), self.name, here, dest)
            yield hold, self, timeTo(self.accel, self.maxV, here.distanceTo[dest])
            # arrive at next station
            here = dest
            print "{:.1f}s: {} at {}".format(self.sim.now(), self.name, here)
            yield hold, self, here.arrive(self)

    def getOff(self, num):
        if self.p >= num:
            print "  {} passengers got off".format(num)
            self.p -= num
        else:
            num = self.p
            print "  train is empty - only {} passengers got off".format(num)
            self.p = 0

    def getOn(self, num):
        if (self.maxP is None) or (self.p + num <= self.maxP):
            print "  {} passengers got on".format(num)
            self.p += num
        else:
            num = self.maxp - self.p
            print "  train is full - only {} passengers got on".format(num)
            self.p = self.maxp

class TrackNode(object):
    def __init__(self, name, delay=5.0):
        self.name = name
        self.delay = delay
        self.distanceTo = {}
    def arrive(self, train):
        pass
    def __str__(self):
        return self.name

class Station(TrackNode):
    def arrive(self, train):
        train.getOff(randint(1,15))
        train.getOn(randint(1,15))
        return self.delay

class Switch(TrackNode):
    def arrive(self, train):
        print("  switching tracks")
        return self.delay

class SampleRailroad(Simulation):
    def run(self, maxTime=100.0):
        self.initialize()
        # create places
        x = Switch("switch x", 20.0)
        A = Station("Station A", 24.0)
        B = Station("Station B", 27.0)
        C = Station("Station C", 25.0)
        y = Switch("switch y", 18.0)
        # distances between places
        x.distanceTo[A] = 50.0
        A.distanceTo[B] = 5000.0
        B.distanceTo[C] = 2000.0
        C.distanceTo[y] = 80.0
        y.distanceTo[C] = 80.0
        C.distanceTo[B] = 2000.0
        B.distanceTo[A] = 5000.0
        A.distanceTo[x] = 50.0
        # set up first train
        sf = Train("Santa Fe 219", self)
        self.activate(sf, sf.roll(cycle([A,B,C,y,C,B,A,x])), at=0.0)
        # set up second train
        cn = Train("Canadian National 41", self, maxPassengers=200)
        self.activate(cn, cn.roll(cycle([C,B,A,x,A,B,C,y])), at=5.0)
        # start simulating!
        self.simulate(until=maxTime)

def main():
    rr = SampleRailroad()
    rr.run(800.0)

if __name__=="__main__":
    main()