Пример #1
0
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _timenow = _propTimenow.GetState(context);
            double timenow  = Convert.ToDouble(_timenow.StateValue);
            IState _salida  = _propSalida.GetState(context);
            double salida   = Convert.ToDouble(_salida.StateValue);
            IState _aux     = _propAux.GetState(context);
            double aux      = Convert.ToDouble(_aux.StateValue);

            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc = v.MipYc;
                vectores.MipYv = v.MipYv;
                vectores.MipTc = v.MipTc;
                vectores.MipTv = v.MipTv;

                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.PlYc          = v.PlYc;
                vectores.PlYv          = v.PlYv;
                vectores.PlTc          = v.PlTc;
                vectores.PlTv          = v.PlTv;
                vectores.PlTcmalo      = v.PlTcmalo;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }


            /*
             *              vectores.PlTc[30, 0]; #flujo pala 1
             *              vectores.PlTc[30, 1]; #flujo pala 2
             *              etc etc
             *
             *              vectores.Din1[x, 4] = pala1; #Ultimo tiempo que pala 1 tuvo una asignacion
             *              vectores.Din1[x + 30, 4] = pala2; #Ultimo tiempo que pala 2 tuvo una asignacion
             *              vectores.Din1[x + 60, 4] = pala3; etc
             *              vectores.Din1[x + 90, 4] = pala4;
             *              vectores.Din1[x + 120, 4] = pala5;
             *              vectores.Din1[x + 150, 4] = pala6;
             *              vectores.Din1[x + 180, 4] = pala7;
             *              vectores.Din1[x + 210, 4] = pala8;
             *              vectores.Din1[x + 240, 4] = pala9;
             *              vectores.Din1[x + 270, 4] = pala10;
             *              vectores.Din1[x + 300, 4] = pala11;
             *              vectores.Din1[x + 330, 4] = pala12;
             *              vectores.Din1[x + 360, 4] = pala13;
             *              vectores.Din1[x + 390, 4] = pala14;
             *              vectores.Din1[x + 420, 4] = pala15;
             *              vectores.Din1[x + 450, 4] = pala16;
             *              vectores.Din1[x + 480, 4] = pala17;
             *              vectores.Din1[x + 510, 4] = pala18;
             *              vectores.Din1[x + 540, 4] = pala19;
             *              vectores.Din1[x + 570, 4] = pala20;
             *
             *              timenow; #tiempo actual
             *
             *
             *
             */
            double[] tij_op   = new double[20];
            double[] desv_tij = new double[20];
            double[] Puntaje  = new double[20] {
                0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
            };
            double Pmax     = -10000000000;
            double Destino  = 0;
            double botadero = 0;

            for (int i = 0; i < 20; i++)
            {
                if (vectores.Din1[30 * i, 4] > 0)
                {
                    tij_op[i]   = 1 / ((vectores.PlTc[30, i] + 1) / 300);
                    desv_tij[i] = (timenow - vectores.Din1[30 * i, 4]) - tij_op[i];
                    Puntaje[i]  = desv_tij[i] / tij_op[i];
                }
            }

            for (int j = 0; j < 20; j++)
            {
                if (vectores.Din1[30 * j, 4] > 0)
                {
                    {
                        if (Puntaje[j] > Pmax)
                        {
                            Pmax    = Puntaje[j];
                            Destino = j + 1;
                        }
                    }
                }
            }

            if (Destino == 4 || Destino == 5 || Destino == 6 || Destino == 10 || Destino == 17 || Destino == 18)
            {
                botadero = 1;
            }

            else if (Destino == 7 || Destino == 8 || Destino == 12)
            {
                botadero = 4;
            }

            else if (Destino == 1 || Destino == 2 || Destino == 3 || Destino == 11)
            {
                botadero = 8;
            }

            else
            {
                Destino  = Pmax;
                botadero = 3;
            }
            _salida.StateValue = Destino;
            _aux.StateValue    = botadero;

            sr.serializa(vectores);

            sr.closeStream();
            return(ExitType.FirstExit);
        }
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _pala0  = _propPala0.GetState(context);
            double pala0   = Convert.ToDouble(_pala0.StateValue);
            IState _pala1  = _propPala1.GetState(context);
            double pala1   = Convert.ToDouble(_pala1.StateValue);
            IState _pala2  = _propPala2.GetState(context);
            double pala2   = Convert.ToDouble(_pala2.StateValue);
            IState _pala3  = _propPala3.GetState(context);
            double pala3   = Convert.ToDouble(_pala3.StateValue);
            IState _pala4  = _propPala4.GetState(context);
            double pala4   = Convert.ToDouble(_pala4.StateValue);
            IState _pala5  = _propPala5.GetState(context);
            double pala5   = Convert.ToDouble(_pala5.StateValue);
            IState _pala6  = _propPala6.GetState(context);
            double pala6   = Convert.ToDouble(_pala6.StateValue);
            IState _pala7  = _propPala7.GetState(context);
            double pala7   = Convert.ToDouble(_pala7.StateValue);
            IState _pala8  = _propPala8.GetState(context);
            double pala8   = Convert.ToDouble(_pala8.StateValue);
            IState _pala9  = _propPala9.GetState(context);
            double pala9   = Convert.ToDouble(_pala9.StateValue);
            IState _pala10 = _propPala10.GetState(context);
            double pala10  = Convert.ToDouble(_pala10.StateValue);
            IState _pala11 = _propPala11.GetState(context);
            double pala11  = Convert.ToDouble(_pala11.StateValue);
            IState _pala12 = _propPala12.GetState(context);
            double pala12  = Convert.ToDouble(_pala12.StateValue);
            IState _pala13 = _propPala13.GetState(context);
            double pala13  = Convert.ToDouble(_pala13.StateValue);
            IState _pala14 = _propPala14.GetState(context);
            double pala14  = Convert.ToDouble(_pala14.StateValue);
            IState _pala15 = _propPala15.GetState(context);
            double pala15  = Convert.ToDouble(_pala15.StateValue);
            IState _pala16 = _propPala16.GetState(context);
            double pala16  = Convert.ToDouble(_pala16.StateValue);
            IState _pala17 = _propPala17.GetState(context);
            double pala17  = Convert.ToDouble(_pala17.StateValue);
            IState _pala18 = _propPala18.GetState(context);
            double pala18  = Convert.ToDouble(_pala18.StateValue);
            IState _pala19 = _propPala19.GetState(context);
            double pala19  = Convert.ToDouble(_pala19.StateValue);
            IState _pala20 = _propPala20.GetState(context);
            double pala20  = Convert.ToDouble(_pala20.StateValue);



            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc    = v.MipYc;
                vectores.MipYv    = v.MipYv;
                vectores.MipTc    = v.MipTc;
                vectores.MipTv    = v.MipTv;
                vectores.PlYc     = v.PlYc;
                vectores.PlYv     = v.PlYv;
                vectores.PlTc     = v.PlTc;
                vectores.PlTv     = v.PlTv;
                vectores.PlTcmalo = v.PlTcmalo;


                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }

            if (pala0 == 0)
            {
                for (int x = 0; x < 30; x = x + 1)
                {
                    vectores.Din1[x, 4]       = pala1;
                    vectores.Din1[x + 30, 4]  = pala2;
                    vectores.Din1[x + 60, 4]  = pala3;
                    vectores.Din1[x + 90, 4]  = pala4;
                    vectores.Din1[x + 120, 4] = pala5;
                    vectores.Din1[x + 150, 4] = pala6;
                    vectores.Din1[x + 180, 4] = pala7;
                    vectores.Din1[x + 210, 4] = pala8;
                    vectores.Din1[x + 240, 4] = pala9;
                    vectores.Din1[x + 270, 4] = pala10;
                    vectores.Din1[x + 300, 4] = pala11;
                    vectores.Din1[x + 330, 4] = pala12;
                    vectores.Din1[x + 360, 4] = pala13;
                    vectores.Din1[x + 390, 4] = pala14;
                    vectores.Din1[x + 420, 4] = pala15;
                    vectores.Din1[x + 450, 4] = pala16;
                    vectores.Din1[x + 480, 4] = pala17;
                    vectores.Din1[x + 510, 4] = pala18;
                    vectores.Din1[x + 540, 4] = pala19;
                    vectores.Din1[x + 570, 4] = pala20;
                }
            }


            //filas
            else if (pala0 == 1)
            {
                vectores.Fila[0]  = pala1;
                vectores.Fila[1]  = pala2;
                vectores.Fila[2]  = pala3;
                vectores.Fila[3]  = pala4;
                vectores.Fila[4]  = pala5;
                vectores.Fila[5]  = pala6;
                vectores.Fila[6]  = pala7;
                vectores.Fila[7]  = pala8;
                vectores.Fila[8]  = pala9;
                vectores.Fila[9]  = pala10;
                vectores.Fila[10] = pala11;
                vectores.Fila[11] = pala12;
                vectores.Fila[12] = pala13;
                vectores.Fila[13] = pala14;
                vectores.Fila[14] = pala15;
                vectores.Fila[15] = pala16;
                vectores.Fila[16] = pala17;
                vectores.Fila[17] = pala18;
                vectores.Fila[18] = pala19;
                vectores.Fila[19] = pala20;
            }



            //numero alocados

            else if (pala0 == 3)
            {
                for (int x = 0; x < 30; x = x + 1)
                {
                    vectores.Din1[x, 5]       = pala1;
                    vectores.Din1[x + 30, 5]  = pala2;
                    vectores.Din1[x + 60, 5]  = pala3;
                    vectores.Din1[x + 90, 5]  = pala4;
                    vectores.Din1[x + 120, 5] = pala5;
                    vectores.Din1[x + 150, 5] = pala6;
                    vectores.Din1[x + 180, 5] = pala7;
                    vectores.Din1[x + 210, 5] = pala8;
                    vectores.Din1[x + 240, 5] = pala9;
                    vectores.Din1[x + 270, 5] = pala10;
                    vectores.Din1[x + 300, 5] = pala11;
                    vectores.Din1[x + 330, 5] = pala12;
                    vectores.Din1[x + 360, 5] = pala13;
                    vectores.Din1[x + 390, 5] = pala14;
                    vectores.Din1[x + 420, 5] = pala15;
                    vectores.Din1[x + 450, 5] = pala16;
                    vectores.Din1[x + 480, 5] = pala17;
                    vectores.Din1[x + 510, 5] = pala18;
                    vectores.Din1[x + 540, 5] = pala19;
                    vectores.Din1[x + 570, 5] = pala20;
                }
            }

            else if (pala0 == 4) //Si es de mineral 0, si es de esteril 1
            {
                vectores.Mine[0]  = pala1;
                vectores.Mine[1]  = pala2;
                vectores.Mine[2]  = pala3;
                vectores.Mine[3]  = pala4;
                vectores.Mine[4]  = pala5;
                vectores.Mine[5]  = pala6;
                vectores.Mine[6]  = pala7;
                vectores.Mine[7]  = pala8;
                vectores.Mine[8]  = pala9;
                vectores.Mine[9]  = pala10;
                vectores.Mine[10] = pala11;
                vectores.Mine[11] = pala12;
                vectores.Mine[12] = pala13;
                vectores.Mine[13] = pala14;
                vectores.Mine[14] = pala15;
                vectores.Mine[15] = pala16;
                vectores.Mine[16] = pala17;
                vectores.Mine[17] = pala18;
                vectores.Mine[18] = pala19;
                vectores.Mine[19] = pala20;
            }

            else if (pala0 == 5) //numero de camiones que salieron de cada pala hacia el chancador
            {
                vectores.Mineralocado[0]  = pala1;
                vectores.Mineralocado[1]  = pala2;
                vectores.Mineralocado[2]  = pala3;
                vectores.Mineralocado[3]  = pala4;
                vectores.Mineralocado[4]  = pala5;
                vectores.Mineralocado[5]  = pala6;
                vectores.Mineralocado[6]  = pala7;
                vectores.Mineralocado[7]  = pala8;
                vectores.Mineralocado[8]  = pala9;
                vectores.Mineralocado[9]  = pala10;
                vectores.Mineralocado[10] = pala11;
                vectores.Mineralocado[11] = pala12;
                vectores.Mineralocado[12] = pala13;
                vectores.Mineralocado[13] = pala14;
                vectores.Mineralocado[14] = pala15;
                vectores.Mineralocado[15] = pala16;
                vectores.Mineralocado[16] = pala17;
                vectores.Mineralocado[17] = pala18;
                vectores.Mineralocado[18] = pala19;
                vectores.Mineralocado[19] = pala20;
            }

            else if (pala0 == 6) //numero de camiones que salieron de cada pala el esteril
            {
                vectores.Lastrealocado[0]  = pala1 * vectores.Mine[0];
                vectores.Lastrealocado[1]  = pala2 * vectores.Mine[1];
                vectores.Lastrealocado[2]  = pala3 * vectores.Mine[2];
                vectores.Lastrealocado[3]  = pala4 * vectores.Mine[3];
                vectores.Lastrealocado[4]  = pala5 * vectores.Mine[4];
                vectores.Lastrealocado[5]  = pala6 * vectores.Mine[5];
                vectores.Lastrealocado[6]  = pala7 * vectores.Mine[6];
                vectores.Lastrealocado[7]  = pala8 * vectores.Mine[7];
                vectores.Lastrealocado[8]  = pala9 * vectores.Mine[8];
                vectores.Lastrealocado[9]  = pala10 * vectores.Mine[9];
                vectores.Lastrealocado[10] = pala11 * vectores.Mine[10];
                vectores.Lastrealocado[11] = pala12 * vectores.Mine[11];
                vectores.Lastrealocado[12] = pala13 * vectores.Mine[12];
                vectores.Lastrealocado[13] = pala14 * vectores.Mine[13];
                vectores.Lastrealocado[14] = pala15 * vectores.Mine[14];
                vectores.Lastrealocado[15] = pala16 * vectores.Mine[15];
                vectores.Lastrealocado[16] = pala17 * vectores.Mine[16];
                vectores.Lastrealocado[17] = pala18 * vectores.Mine[17];
                vectores.Lastrealocado[18] = pala19 * vectores.Mine[18];
                vectores.Lastrealocado[19] = pala20 * vectores.Mine[19];
            }


            sr.serializa(vectores);
            sr.closeStream();
            return(ExitType.FirstExit);
        }
Пример #3
0
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _timenow = _propTimenow.GetState(context);
            double timenow  = Convert.ToDouble(_timenow.StateValue);
            IState _salida  = _propSalida.GetState(context);
            double salida   = Convert.ToDouble(_salida.StateValue);
            IState _aux     = _propAux.GetState(context);
            double aux      = Convert.ToDouble(_aux.StateValue);

            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc    = v.MipYc;
                vectores.MipYv    = v.MipYv;
                vectores.MipTc    = v.MipTc;
                vectores.MipTv    = v.MipTv;
                vectores.PlYc     = v.PlYc;
                vectores.PlYv     = v.PlYv;
                vectores.PlTc     = v.PlTc;
                vectores.PlTv     = v.PlTv;
                vectores.PlTcmalo = v.PlTcmalo;

                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }

            /* programa aca*/



            int[,] Resultado = new int[20 + 1, 3];
            int[] VectorResultado = new int[20 + 1];
            Parallel.For(1, 3, m =>
            {
                VectorResultado = Simulacion(m, timenow);
                for (int i = 0; i < 20 + 1; i++)
                {
                    Resultado[i, m - 1] = VectorResultado[i];
                }
            }
                         )
            ;



            sr.serializa(vectores);
            sr.closeStream();
            return(ExitType.FirstExit);
        }
Пример #4
0
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _pala0  = _propPala0.GetState(context);
            double pala0   = Convert.ToDouble(_pala0.StateValue);
            IState _pala1  = _propPala1.GetState(context);
            double pala1   = Convert.ToDouble(_pala1.StateValue);
            IState _pala2  = _propPala2.GetState(context);
            double pala2   = Convert.ToDouble(_pala2.StateValue);
            IState _pala3  = _propPala3.GetState(context);
            double pala3   = Convert.ToDouble(_pala3.StateValue);
            IState _pala4  = _propPala4.GetState(context);
            double pala4   = Convert.ToDouble(_pala4.StateValue);
            IState _pala5  = _propPala5.GetState(context);
            double pala5   = Convert.ToDouble(_pala5.StateValue);
            IState _pala6  = _propPala6.GetState(context);
            double pala6   = Convert.ToDouble(_pala6.StateValue);
            IState _pala7  = _propPala7.GetState(context);
            double pala7   = Convert.ToDouble(_pala7.StateValue);
            IState _pala8  = _propPala8.GetState(context);
            double pala8   = Convert.ToDouble(_pala8.StateValue);
            IState _pala9  = _propPala9.GetState(context);
            double pala9   = Convert.ToDouble(_pala9.StateValue);
            IState _pala10 = _propPala10.GetState(context);
            double pala10  = Convert.ToDouble(_pala10.StateValue);
            IState _pala11 = _propPala11.GetState(context);
            double pala11  = Convert.ToDouble(_pala11.StateValue);
            IState _pala12 = _propPala12.GetState(context);
            double pala12  = Convert.ToDouble(_pala12.StateValue);
            IState _pala13 = _propPala13.GetState(context);
            double pala13  = Convert.ToDouble(_pala13.StateValue);
            IState _pala14 = _propPala14.GetState(context);
            double pala14  = Convert.ToDouble(_pala14.StateValue);
            IState _pala15 = _propPala15.GetState(context);
            double pala15  = Convert.ToDouble(_pala15.StateValue);
            IState _pala16 = _propPala16.GetState(context);
            double pala16  = Convert.ToDouble(_pala16.StateValue);
            IState _pala17 = _propPala17.GetState(context);
            double pala17  = Convert.ToDouble(_pala17.StateValue);
            IState _pala18 = _propPala18.GetState(context);
            double pala18  = Convert.ToDouble(_pala18.StateValue);
            IState _pala19 = _propPala19.GetState(context);
            double pala19  = Convert.ToDouble(_pala19.StateValue);
            IState _pala20 = _propPala20.GetState(context);
            double pala20  = Convert.ToDouble(_pala20.StateValue);
            IState _pala21 = _propPala21.GetState(context);
            double pala21  = Convert.ToDouble(_pala21.StateValue);
            IState _pala22 = _propPala22.GetState(context);
            double pala22  = Convert.ToDouble(_pala22.StateValue);
            IState _pala23 = _propPala23.GetState(context);
            double pala23  = Convert.ToDouble(_pala23.StateValue);
            IState _pala24 = _propPala24.GetState(context);
            double pala24  = Convert.ToDouble(_pala24.StateValue);
            IState _pala25 = _propPala25.GetState(context);
            double pala25  = Convert.ToDouble(_pala25.StateValue);
            IState _pala26 = _propPala26.GetState(context);
            double pala26  = Convert.ToDouble(_pala26.StateValue);
            IState _pala27 = _propPala27.GetState(context);
            double pala27  = Convert.ToDouble(_pala27.StateValue);
            IState _pala28 = _propPala28.GetState(context);
            double pala28  = Convert.ToDouble(_pala28.StateValue);
            IState _pala29 = _propPala29.GetState(context);
            double pala29  = Convert.ToDouble(_pala29.StateValue);
            IState _pala30 = _propPala30.GetState(context);
            double pala30  = Convert.ToDouble(_pala30.StateValue);
            IState _pala31 = _propPala31.GetState(context);
            double pala31  = Convert.ToDouble(_pala31.StateValue);
            IState _pala32 = _propPala32.GetState(context);
            double pala32  = Convert.ToDouble(_pala32.StateValue);
            IState _pala33 = _propPala33.GetState(context);
            double pala33  = Convert.ToDouble(_pala33.StateValue);
            IState _pala34 = _propPala34.GetState(context);
            double pala34  = Convert.ToDouble(_pala34.StateValue);
            IState _pala35 = _propPala35.GetState(context);
            double pala35  = Convert.ToDouble(_pala35.StateValue);
            IState _pala36 = _propPala36.GetState(context);
            double pala36  = Convert.ToDouble(_pala36.StateValue);
            IState _pala37 = _propPala37.GetState(context);
            double pala37  = Convert.ToDouble(_pala37.StateValue);
            IState _pala38 = _propPala38.GetState(context);
            double pala38  = Convert.ToDouble(_pala38.StateValue);
            IState _pala39 = _propPala39.GetState(context);
            double pala39  = Convert.ToDouble(_pala39.StateValue);
            IState _pala40 = _propPala40.GetState(context);
            double pala40  = Convert.ToDouble(_pala40.StateValue);
            IState _pala41 = _propPala41.GetState(context);
            double pala41  = Convert.ToDouble(_pala41.StateValue);
            IState _pala42 = _propPala42.GetState(context);
            double pala42  = Convert.ToDouble(_pala42.StateValue);
            IState _pala43 = _propPala43.GetState(context);
            double pala43  = Convert.ToDouble(_pala43.StateValue);
            IState _pala44 = _propPala44.GetState(context);
            double pala44  = Convert.ToDouble(_pala44.StateValue);
            IState _pala45 = _propPala45.GetState(context);
            double pala45  = Convert.ToDouble(_pala45.StateValue);
            IState _pala46 = _propPala46.GetState(context);
            double pala46  = Convert.ToDouble(_pala46.StateValue);
            IState _pala47 = _propPala47.GetState(context);
            double pala47  = Convert.ToDouble(_pala47.StateValue);
            IState _pala48 = _propPala48.GetState(context);
            double pala48  = Convert.ToDouble(_pala48.StateValue);
            IState _pala49 = _propPala49.GetState(context);
            double pala49  = Convert.ToDouble(_pala49.StateValue);
            IState _pala50 = _propPala50.GetState(context);
            double pala50  = Convert.ToDouble(_pala50.StateValue);
            IState _pala51 = _propPala51.GetState(context);
            double pala51  = Convert.ToDouble(_pala51.StateValue);
            IState _pala52 = _propPala52.GetState(context);
            double pala52  = Convert.ToDouble(_pala52.StateValue);
            IState _pala53 = _propPala53.GetState(context);
            double pala53  = Convert.ToDouble(_pala53.StateValue);
            IState _pala54 = _propPala54.GetState(context);
            double pala54  = Convert.ToDouble(_pala54.StateValue);
            IState _pala55 = _propPala55.GetState(context);
            double pala55  = Convert.ToDouble(_pala55.StateValue);
            IState _pala56 = _propPala56.GetState(context);
            double pala56  = Convert.ToDouble(_pala56.StateValue);
            IState _pala57 = _propPala57.GetState(context);
            double pala57  = Convert.ToDouble(_pala57.StateValue);
            IState _pala58 = _propPala58.GetState(context);
            double pala58  = Convert.ToDouble(_pala58.StateValue);
            IState _pala59 = _propPala59.GetState(context);
            double pala59  = Convert.ToDouble(_pala59.StateValue);
            IState _pala60 = _propPala60.GetState(context);
            double pala60  = Convert.ToDouble(_pala60.StateValue);
            // Example of how to get the value of a step property.


            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc    = v.MipYc;
                vectores.MipYv    = v.MipYv;
                vectores.MipTc    = v.MipTc;
                vectores.MipTv    = v.MipTv;
                vectores.PlYc     = v.PlYc;
                vectores.PlYv     = v.PlYv;
                vectores.PlTc     = v.PlTc;
                vectores.PlTv     = v.PlTv;
                vectores.PlTcmalo = v.PlTcmalo;


                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }

            /*
             * //tiempos
             *
             * vectores.Din1[0, 7] = pala0;
             * vectores.Din1[30, 7] = pala1;
             * vectores.Din1[60, 7] = pala2;
             * vectores.Din1[90, 7] = pala3;
             * vectores.Din1[120, 7] = pala4;
             * vectores.Din1[150, 7] = pala5;
             * vectores.Din1[180, 7] = pala6;
             * vectores.Din1[210, 7] = pala7;
             * vectores.Din1[480, 7] = pala8;
             * vectores.Din1[510, 7] = pala9;
             *
             * vectores.Din1[3, 7] = pala10;
             * vectores.Din1[33, 7] = pala11;
             * vectores.Din1[63, 7] = pala12;
             * vectores.Din1[93, 7] = pala13;
             * vectores.Din1[123, 7] = pala14;
             * vectores.Din1[153, 7] = pala15;
             * vectores.Din1[183, 7] = pala16;
             * vectores.Din1[213, 7] = pala17;
             * vectores.Din1[483, 7] = pala18;
             * vectores.Din1[513, 7] = pala19;
             *
             * vectores.Din1[7, 7] = pala20;
             * vectores.Din1[37, 7] = pala31;
             * vectores.Din1[67, 7] = pala32;
             * vectores.Din1[97, 7] = pala33;
             * vectores.Din1[127, 7] = pala34;
             * vectores.Din1[157, 7] = pala35;
             * vectores.Din1[187, 7] = pala36;
             * vectores.Din1[217, 7] = pala37;
             * vectores.Din1[487, 7] = pala38;
             * vectores.Din1[517, 7] = pala39;
             *
             * vectores.Din1[17, 7] = pala40;
             * vectores.Din1[47, 7] = pala41;
             * vectores.Din1[77, 7] = pala42;
             * vectores.Din1[107, 7] = pala43;
             * vectores.Din1[137, 7] = pala44;
             * vectores.Din1[167, 7] = pala45;
             * vectores.Din1[197, 7] = pala46;
             * vectores.Din1[227, 7] = pala47;
             * vectores.Din1[497, 7] = pala48;
             * vectores.Din1[527, 7] = pala49;
             */
            for (int x = 0; x < 30; x = x + 1)
            {
                vectores.Din1[x, 4]       = pala1;
                vectores.Din1[x + 30, 4]  = pala2;
                vectores.Din1[x + 60, 4]  = pala3;
                vectores.Din1[x + 90, 4]  = pala4;
                vectores.Din1[x + 120, 4] = pala5;
                vectores.Din1[x + 150, 4] = pala6;
                vectores.Din1[x + 180, 4] = pala7;
                vectores.Din1[x + 210, 4] = pala8;
                vectores.Din1[x + 240, 4] = pala9;
                vectores.Din1[x + 270, 4] = pala10;
                vectores.Din1[x + 300, 4] = pala11;
                vectores.Din1[x + 330, 4] = pala12;
                vectores.Din1[x + 360, 4] = pala13;
                vectores.Din1[x + 390, 4] = pala14;
                vectores.Din1[x + 420, 4] = pala15;
                vectores.Din1[x + 450, 4] = pala16;
                vectores.Din1[x + 480, 4] = pala17;
                vectores.Din1[x + 510, 4] = pala18;
                vectores.Din1[x + 540, 4] = pala19;
                vectores.Din1[x + 570, 4] = pala20;
            }
            vectores.Fila[0]  = pala21;
            vectores.Fila[1]  = pala22;
            vectores.Fila[2]  = pala23;
            vectores.Fila[3]  = pala24;
            vectores.Fila[4]  = pala25;
            vectores.Fila[5]  = pala26;
            vectores.Fila[6]  = pala27;
            vectores.Fila[7]  = pala28;
            vectores.Fila[8]  = pala29;
            vectores.Fila[9]  = pala30;
            vectores.Fila[10] = pala31;
            vectores.Fila[11] = pala32;
            vectores.Fila[12] = pala33;
            vectores.Fila[13] = pala34;
            vectores.Fila[14] = pala35;
            vectores.Fila[15] = pala36;
            vectores.Fila[16] = pala37;
            vectores.Fila[17] = pala38;
            vectores.Fila[18] = pala39;
            vectores.Fila[19] = pala40;

            for (int x = 0; x < 30; x = x + 1)
            {
                vectores.Din1[x, 5]       = pala41;
                vectores.Din1[x + 30, 5]  = pala42;
                vectores.Din1[x + 60, 5]  = pala43;
                vectores.Din1[x + 90, 5]  = pala44;
                vectores.Din1[x + 120, 5] = pala45;
                vectores.Din1[x + 150, 5] = pala46;
                vectores.Din1[x + 180, 5] = pala47;
                vectores.Din1[x + 210, 5] = pala48;
                vectores.Din1[x + 240, 5] = pala49;
                vectores.Din1[x + 270, 5] = pala50;
                vectores.Din1[x + 300, 5] = pala51;
                vectores.Din1[x + 330, 5] = pala52;
                vectores.Din1[x + 360, 5] = pala53;
                vectores.Din1[x + 390, 5] = pala54;
                vectores.Din1[x + 420, 5] = pala55;
                vectores.Din1[x + 450, 5] = pala56;
                vectores.Din1[x + 480, 5] = pala57;
                vectores.Din1[x + 510, 5] = pala58;
                vectores.Din1[x + 540, 5] = pala59;
                vectores.Din1[x + 570, 5] = pala60;
            }


            sr.serializa(vectores);
            sr.closeStream();
            return(ExitType.FirstExit);
        }
Пример #5
0
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _timenow = _propTimenow.GetState(context);
            double timenow  = Convert.ToDouble(_timenow.StateValue);

            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc = v.MipYc;
                vectores.MipYv = v.MipYv;
                vectores.MipTc = v.MipTc;
                vectores.MipTv = v.MipTv;

                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.PlYc          = v.PlYc;
                vectores.PlYv          = v.PlYv;
                vectores.PlTc          = v.PlTc;
                vectores.PlTv          = v.PlTv;
                vectores.PlTcmalo      = v.PlTcmalo;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }

            int    i = 30;
            int    j = 20;
            double TA;

            TA = 0.0158333333; //tiempo de aculatado
            double expr;

            expr = 0; //variable auxiliar

            GRBEnv env = new GRBEnv();

            GRBModel disp = new GRBModel(env);

            GRBLinExpr L = 320; //Numero auxiliar de velocidad de descarga
            GRBLinExpr M = 110; //Numero max camiones

            double[,] Tcvij = new double[30, 20];
            double[,] Tccij = new double[30, 20];

            for (int a = 0; a < i; a = a + 1) // saca tiempo de viaje total
            {
                for (int b = 0; b < j; b = b + 1)
                {
                    Tcvij[a, b] = vectores.TCV[a, b] + vectores.TDi[a];

                    if (vectores.Destij[a, b] != 0)
                    {
                        Tccij[a, b] = vectores.TCC[a, b] + vectores.TCj[b] + TA;
                    }
                    else
                    {
                        Tccij[a, b] = 9999;
                    }
                }
            }



            for (int b = 0; b < j; b = b + 1)
            {
                expr = 0;

                for (int a = 0; a < i; a = a + 1)
                {
                    expr = expr + vectores.Destij[a, b];
                }
                if (expr == 0)
                {
                    vectores.Uj[b] = 0;
                }
            }


            GRBVar[,] Ycij = new GRBVar[30, 20]; //Camiones cargados de i a j
            GRBVar[,] Yvij = new GRBVar[30, 20]; //Camiones vacios de i a j
            GRBVar[,] Tcij = new GRBVar[30, 20]; //Throughtput cargado de i a j
            GRBVar[,] Tvij = new GRBVar[30, 20]; //Throughtput vacio de i a j

            for (int a = 0; a < i; a = a + 1)
            {
                for (int b = 0; b < j; b = b + 1)
                {
                    Ycij[a, b] = disp.AddVar(0.0, 1000, 0.0, GRB.CONTINUOUS, "Yc" + a + "_" + b);
                    Yvij[a, b] = disp.AddVar(0.0, 1000, 0.0, GRB.CONTINUOUS, "Yv" + a + "_" + b);
                    Tcij[a, b] = disp.AddVar(0.0, 100000, 0.0, GRB.CONTINUOUS, "Tc" + a + "_" + b); //tonelaje cargado pasando hacia j en toneladas
                    Tvij[a, b] = disp.AddVar(0.0, 100000, 0.0, GRB.CONTINUOUS, "Tv" + a + "_" + b); //tonelaje vacio pasando hacia j en toneladas
                }
            }

            for (int a = 0; a < i; a = a + 1)
            {
                for (int b = 0; b < j; b = b + 1)
                {
                    disp.AddConstr(Tcij[a, b] == ((Ycij[a, b] / Tccij[a, b]) * L), "tonelaje_minimo" + a);
                    disp.AddConstr(Tvij[a, b] == ((Yvij[a, b] / Tcvij[a, b]) * L), "tonelaje_maximo" + b);
                }
            }

            for (int b = 0; b < j; b = b + 1)
            {
                GRBLinExpr P = 0;

                for (int a = 0; a < i; a = a + 1)
                {
                    P = P + Tcij[a, b];
                }
                disp.AddConstr((vectores.Dj[b] * vectores.Uj[b]), GRB.LESS_EQUAL, P, "tonelaje_minimo_de_palas" + b);
            }



            for (int a = 0; a < i; a = a + 1)
            {
                GRBLinExpr P = 0;

                for (int b = 0; b < j; b = b + 1)
                {
                    P = P + Tcij[a, b];
                }
                disp.AddConstr(P, GRB.LESS_EQUAL, vectores.RUi[a] * vectores.Ui[a], "tonelaje_minimo_a_basureros" + a);
            }

            //Restriccion del chancador independiente

            for (int a = 1; a < i; a = a + 1)
            {
                GRBLinExpr P = 0;

                for (int b = 0; b < j; b = b + 1)
                {
                    P = P + Tcij[a, b];
                }
                disp.AddConstr(vectores.RLi[a] * vectores.Ui[a], GRB.LESS_EQUAL, P, "tonelajecargado2" + a);
            }

            for (int a = 0; a < i; a = a + 1)
            {
                GRBLinExpr P = 0;

                for (int b = 0; b < j; b = b + 1)
                {
                    P = P + Tcij[a, b];
                }
                disp.AddConstr(vectores.Ui[a] * (vectores.RLi[a] - (vectores.PlTcmalo[0, 3] * (1 - vectores.Uj[3]) + vectores.PlTcmalo[0, 4] * (1 - vectores.Uj[4]) + vectores.PlTcmalo[0, 5] * (1 - vectores.Uj[5]) + vectores.PlTcmalo[0, 16] * (1 - vectores.Uj[16]) + vectores.PlTcmalo[0, 17] * (1 - vectores.Uj[17]))), GRB.LESS_EQUAL, P, "tonelajecargado2" + a);
            }



            GRBLinExpr Qonda = 0;

            for (int a = 0; a < i; a = a + 1)
            {
                for (int b = 0; b < j; b = b + 1)
                {
                    Qonda = (Qonda + Ycij[a, b] + Yvij[a, b]);
                }
            }

            disp.AddConstr(Qonda, GRB.LESS_EQUAL, M, "maximocamiones");



            for (int b = 0; b < j; b = b + 1)
            {
                GRBLinExpr P1 = 0;
                GRBLinExpr P2 = 0;

                for (int a = 0; a < i; a = a + 1)
                {
                    P1 = P1 + Tcij[a, b];
                    P2 = P2 + Tvij[a, b];
                }
                disp.AddConstr(P1, GRB.EQUAL, P2, "conservaciondeflujo" + i);
            }


            for (int a = 0; a < i; a = a + 1)
            {
                GRBLinExpr P1 = 0;
                GRBLinExpr P2 = 0;

                for (int b = 0; b < j; b = b + 1)
                {
                    P1 = P1 + Tcij[a, b];
                    P2 = P2 + Tvij[a, b];
                }

                disp.AddConstr(P1, GRB.EQUAL, P2, "conservaciondeflujo2" + i);
            }


            GRBLinExpr Pobj = 0;

            for (int a = 0; a < i; a = a + 1)
            {
                for (int b = 0; b < j; b = b + 1)
                {
                    Pobj = Pobj + Ycij[a, b] + Yvij[a, b];
                }
            }

            disp.Update();

            disp.SetObjective(Pobj, GRB.MINIMIZE);

            disp.Optimize();


            double sumayc = 0;
            double sumayv = 0;
            double sumatc = 0;
            double sumatv = 0;

            for (int a = 0; a < i; a = a + 1)
            {
                double sumajyc = 0;
                double sumajyv = 0;
                double sumajtc = 0;
                double sumajtv = 0;

                for (int b = 0; b < j; b = b + 1)
                {
                    vectores.PlYc[a, b] = Convert.ToDouble(Ycij[a, b].X);
                    vectores.PlYv[a, b] = Convert.ToDouble(Yvij[a, b].X);
                    vectores.PlTc[a, b] = Convert.ToDouble(Tcij[a, b].X);
                    vectores.PlTv[a, b] = Convert.ToDouble(Tvij[a, b].X);
                    sumayc  = sumayc + vectores.PlYc[a, b];             //suma todos los flujos llenos
                    sumayv  = sumayv + vectores.PlYv[a, b];             //suma todos los flujos vacios
                    sumatc  = sumatc + vectores.PlTc[a, b];             //suma todos los camiones llenos
                    sumatv  = sumatv + vectores.PlTv[a, b];             //suma todos los camiones vacios
                    sumajyc = sumajyc + Convert.ToDouble(Ycij[a, b].X); //suma todo el flujo cargado al botadero "a"
                    sumajyv = sumajyv + Convert.ToDouble(Yvij[a, b].X); //suma todo el flujo vacio al botadero "a"
                    sumajtc = sumajtc + Convert.ToDouble(Tcij[a, b].X); //suma todos los camiones cargados yendo a botadero "a"
                    sumajtv = sumajtv + Convert.ToDouble(Tvij[a, b].X); //suma todos los camiones vacios yendo a botadero "a"
                }

                vectores.PlYc[a, 20] = sumajyc; // deja total
                vectores.PlYv[a, 20] = sumajyv; //
                vectores.PlTc[a, 20] = sumajtc; //
                vectores.PlTv[a, 20] = sumajtv; //
            }

            for (int b = 0; b < j; b = b + 1)
            {
                double sumaiyc = 0;
                double sumaiyv = 0;
                double sumaitc = 0;
                double sumaitv = 0;

                for (int a = 0; a < i; a = a + 1)
                {
                    sumaiyc = sumaiyc + Convert.ToDouble(Ycij[a, b].X);
                    sumaiyv = sumaiyv + Convert.ToDouble(Yvij[a, b].X);
                    sumaitc = sumaitc + Convert.ToDouble(Tcij[a, b].X);
                    sumaitv = sumaitv + Convert.ToDouble(Tvij[a, b].X);
                }

                vectores.PlYc[30, b] = sumaiyc;
                vectores.PlYv[30, b] = sumaiyv;
                vectores.PlTc[30, b] = sumaitc;
                vectores.PlTv[30, b] = sumaitv;
            }

            vectores.PlYc[30, 20] = sumayc;
            vectores.PlYv[30, 20] = sumayv;
            vectores.PlTc[30, 20] = sumatc;
            vectores.PlTv[30, 20] = sumatv;

            _timenow.StateValue = vectores.PlTc[0, 20];
            sr.serializa(vectores);
            sr.closeStream();

            disp.Dispose();
            env.Dispose();



            return(ExitType.FirstExit);
        }
Пример #6
0
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _timenow = _propTimenow.GetState(context);
            double timenow  = Convert.ToDouble(_timenow.StateValue);
            IState _salida  = _propSalida.GetState(context);
            double salida   = Convert.ToDouble(_salida.StateValue);
            IState _aux     = _propAux.GetState(context);
            double aux      = Convert.ToDouble(_aux.StateValue);

            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc = v.MipYc;
                vectores.MipYv = v.MipYv;
                vectores.MipTc = v.MipTc;
                vectores.MipTv = v.MipTv;

                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.PlYc          = v.PlYc;
                vectores.PlYv          = v.PlYv;
                vectores.PlTc          = v.PlTc;
                vectores.PlTv          = v.PlTv;
                vectores.PlTcmalo      = v.PlTcmalo;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }

            //i Botaderos. j Palas. k Camiones.

            double[,] Din12 = (double[, ])vectores.Din1.Clone();

            double[,] Tij = (double[, ])vectores.PlTc.Clone();//I Tonelaje por turno


            double[,] xij = new double[30, 20]; //I Tonelaje cumplido hasta el momento de activación del dinámico
            double[,] yij = new double[30, 20]; //Tonelaje requerido para alcazar R
            double[,] Rij = new double[30, 20]; //Porcentaje de cumplimienteo del tonelaje
            double[,] dij = new double[30, 20]; //Desviación del porcentaje de cumplimiento ideal en el momento de activación
            double[,] Mij = new double[30, 20]; //Demanda de camiones por ruta
            double[,] rij = new double[30, 20]; //I Matriz de tiempos de viajes

            double[] tk = new double[81];       //I Tiempo que le toma al camión k llegar al botadero/chancador al que va (si ya está en un botadero este tiempo es 0)
            double[] dk = new double[81];       //I Tiempo esperado de espera del camión k en el botadero/chancador al que va
            double[] ej = new double[30];       //I Tiempo promedio de descarga en el botadero j
            double[] Di = new double[20];       //Demanda de camiones por pala (int?)
            double[] Li = new double[20];       //I Tiempo promedio de carga en la pala i
            double[] Ni = new double[20];       //I Número de camiones en la pala i (int?)
            double[] Ei = new double[20];       //I Camiones en ruta a la pala i (int?)
            double[] Sk = new double[81];       //I Arreglo de variables binarias que indican si el camión tiene una tarea asignada
            int[]    Bk = new int[81];          //Botadero al que está llendo el camión k

            int    numcam       = 81;           //Número total de camiones
            double capacidad    = 300;          //Capacidad de los camiones
            int    numPalas     = 20;           //Palas activas?
            int    numBotaderos = 30;           //Botaderos o chancador activos?
            double R            = 0;

            //Se cálcula el porcentaje de cumplimiento de cuota
            for (int i = 0; i < numBotaderos; i++)
            {
                for (int j = 0; j < numPalas; j++)
                {
                    if (Tij[i, j] != 0)
                    {
                        Rij[i, j] = xij[i, j] / Tij[i, j];
                        R         = R + Rij[i, j];
                    }
                }
            }

            //Se calcula de déficit o superávit
            R = R / (numPalas * numBotaderos);
            for (int i = 0; i < numBotaderos; i++)
            {
                for (int j = 0; j < numPalas; j++)
                {
                    dij[i, j] = Rij[i, j] - R;
                }
            }
            //Se calcula la demanda adeudada para alcanzar R
            double mij_aux;
            int    mij;

            for (int i = 0; i < numBotaderos; i++)
            {
                for (int j = 0; j < numPalas; j++)
                {
                    yij[i, j] = R * Tij[i, j] - xij[i, j];
                    mij_aux   = yij[i, j] / capacidad;
                    mij       = Convert.ToInt32(mij_aux);
                    Mij[i, j] = mij + 1;
                }
            }
            double DTotal = 0;
            double di;

            for (int j = 0; j < numPalas; j++)
            {
                di = 0;
                for (int i = 0; i < numPalas; i++)
                {
                    di = Mij[i, j] + di;
                }

                Di[j]  = di;
                DTotal = DTotal + Di[j];
            }
            //En teoría se debería hacer un chequeo para que la demanda total de camiones no supere
            //el total disponible. De ocurrir, se deben limitar los flujos del PL (*).

            //Criterio de asignación de camiones
            //Se crea el arreglo Wik que contiene los tiempos que demora al camion k llegar a ser
            //atendido por la pala i. El tamaño de este arreglo depende de los camiones disponibles.
            double[,] Wik = new double[numPalas, numcam];//Matriz de tiempos esperados del camión k para ser atendido por la pala i

            for (int k = 0; k < numcam; k++)
            {
                if (Sk[k] != 0) //Sólo para los camiones sin tareas asignadas
                {
                    for (int i = 0; i < numPalas; i++)
                    {
                        Wik[i, k] = Li[i] * (Ni[i] + Ei[i]) - (tk[k] + dk[k] + ej[Bk[k]] + rij[Bk[k], i]);
                    }
                }
            }

            //Modelo
            GRBEnv   env   = new GRBEnv();
            GRBModel model = new GRBModel(env);

            model.ModelName = "Asignación Dinámica";

            GRBVar[,] X = new GRBVar[numPalas, numcam]; //Binaria que define si el camión k se asigna a la pala i

            for (int i = 0; i < numPalas; i++)
            {
                for (int k = 0; k < numcam; k++)
                {
                    X[i, k] = model.AddVar(0.0, 1.0, 0.0, GRB.BINARY, "X" + i + "_" + k);
                }
            }

            model.Update();

            //Restricciones

            for (int k = 0; k < numcam; k++)
            {
                GRBLinExpr suma_Xik = 0;

                for (int i = 0; i < numPalas; i++)
                {
                    suma_Xik = suma_Xik + X[i, k];
                }
                model.AddConstr(suma_Xik, GRB.LESS_EQUAL, Sk[k], "S" + k);
            }

            for (int i = 0; i < numPalas; i++)
            {
                GRBLinExpr suma_Xik = 0;

                for (int k = 0; k < numcam; k++)
                {
                    suma_Xik = suma_Xik + X[i, k];
                }
                model.AddConstr(suma_Xik, GRB.GREATER_EQUAL, Di[i], "D" + i);
            }

            model.Update();

            GRBLinExpr FO = 0;

            for (int i = 0; i < numPalas; i++)
            {
                for (int k = 0; k < numcam; k++)
                {
                    FO = FO + Wik[i, k] * X[i, k];
                }
            }
            model.SetObjective(FO, GRB.MINIMIZE);

            model.Update();

            model.Optimize();

            sr.serializa(vectores);

            sr.closeStream();
            return(ExitType.FirstExit);
        }
Пример #7
0
        /// <summary>
        /// Method called when a process token executes the step.
        /// </summary>
        public ExitType Execute(IStepExecutionContext context)
        {
            IState _timenow   = _propTimenow.GetState(context);
            double timenow    = Convert.ToDouble(_timenow.StateValue);
            IState _salida    = _propSalida.GetState(context);
            double salida     = Convert.ToDouble(_salida.StateValue);
            IState _aux       = _propAux.GetState(context);
            double aux        = Convert.ToDouble(_aux.StateValue);
            IState _aux2      = _propAux2.GetState(context);
            double aux2       = Convert.ToDouble(_aux2.StateValue);
            IState _numcam    = _propNumCam.GetState(context);
            double numcam     = Convert.ToDouble(_numcam.StateValue);
            IState _needtime  = _propNeedtime.GetState(context);
            double needtime   = Convert.ToDouble(_needtime.StateValue);
            IState _losttones = _propLosttones.GetState(context);
            double losttones  = Convert.ToDouble(_losttones.StateValue);
            IState _needch    = _propNeedch.GetState(context);
            double needch     = Convert.ToDouble(_needch.StateValue);
            // Example of how to get the value of a step property.


            SerializarElement sr = new SerializarElement();

            if (File.Exists(sr.SerializationFile))
            {
                Vect v = sr.deserializa();
                vectores = sr.deserializa();

                vectores.MipYc    = v.MipYc;
                vectores.MipYv    = v.MipYv;
                vectores.MipTc    = v.MipTc;
                vectores.MipTv    = v.MipTv;
                vectores.PlYc     = v.PlYc;
                vectores.PlYv     = v.PlYv;
                vectores.PlTc     = v.PlTc;
                vectores.PlTv     = v.PlTv;
                vectores.PlTcmalo = v.PlTcmalo;


                //Dinamico
                vectores.Din1   = v.Din1;
                vectores.DesCam = v.DesCam;
                vectores.Fila   = v.Fila;

                //LP
                vectores.TCV           = v.TCV;
                vectores.TCC           = v.TCC;
                vectores.TCj           = v.TCj;
                vectores.TDi           = v.TDi;
                vectores.Destij        = v.Destij;
                vectores.Dj            = v.Dj;
                vectores.Uj            = v.Uj;
                vectores.Ui            = v.Ui;
                vectores.RLi           = v.RLi;
                vectores.RUi           = v.RUi;
                vectores.Mine          = v.Mine;
                vectores.Mineralocado  = v.Mineralocado;
                vectores.Lastrealocado = v.Lastrealocado;
            }

            int    y;                              //para entrar al loop
            int    palaneedy        = 0;           //numero de la pala needy
            double caminoipalaneedy = 0;
            double min;                            //var aux para encontrar pala needy

            double[,] Din22 = new double[200, 11]; //matriz con datos para encontrar camion con menos lost_tones
            double[] need_time = new double[600];  //matriz para encontrar al camino mas needy


            for (int x = 0; x < numcam; x++)    //construye lo que falta de la matriz de los camiones (lo que depende de la pala needy)
            {
                if (vectores.DesCam[x, 3] == 1) // indicador si necesita pala
                {
                    Din22[x, 8] = 1;
                }
            }


            for (int vuelta = 0; vuelta < 20; vuelta++)
            {
                for (int x = 0; x < 30; x++)
                {
                    vectores.Din1[x + 30 * (vuelta), 1] = vectores.PlTc[30, vuelta]; //rescatar el haulage por hora de cada camino
                }
            }

            int desesperacion;
            int desesperacion2;

            desesperacion  = 0;
            desesperacion2 = 0;
            for (int x = 0; x < 600; x++) //recupera los tonelajes del PL
            {
                vectores.Din1[x, 0] = vectores.PlTc[desesperacion, desesperacion2];
                desesperacion       = desesperacion + 1;
                if (desesperacion > 29)
                {
                    desesperacion2 = desesperacion2 + 1;
                    desesperacion  = 0;
                }
            }

            for (int x = 0; x < 200; x++)
            {
                Din22[x, 0] = x + 1;                                                 //numero camion
                Din22[x, 1] = 320;                                                   //capacidad
                Din22[x, 2] = vectores.PlTc[30, 20];                                 //total rate shovels
                Din22[x, 3] = (vectores.PlYc[30, 20] + vectores.PlYv[30, 20]) * 320; //total trucks en toneladas
            }

            int newaux      = 0;
            int camionpiola = 100;

            min = 999999;//var aux para encontrar pala needy
            double min3;

            double[,] Din12 = (double[, ])vectores.Din1.Clone(); //matriz con datos para encontrar camion con menos lost_tones

            do
            {
                min3 = 0;
                min  = 999999;
                y    = 0;

                int indicemin1 = 900000000; //var aux para encontrar pala needy

                for (int x = 0; x < 600; x++)
                {
                    if (Din12[x, 0] != 0)
                    {
                        need_time[x] = Din12[x, 4] + ((Din12[x, 0] / Din12[x, 1]) * (Din12[x, 5] * 300 - Din12[x, 1]) / (Din12[x, 0]));

                        //need_time[x] = Din12[x, 4] + ((Din12[x, 0] / Din12[x, 1]) * (Din12[x, 5] * 300 - Din12[x, 1] * timenow)) / (Din12[x, 7] * 300 / (timenow));

                        if (need_time[x] < min) //encuentra el menor y lo guarda en palaneedy
                        {
                            indicemin1 = x;

                            palaneedy = Convert.ToInt32(Din12[x, 3]);

                            min = need_time[x];

                            caminoipalaneedy = Din12[x, 2];
                        }
                    }
                }

                //if (min < timenow)
                //{
                //Din12[indicemin1, 7] = Din12[indicemin1, 7] + 1;


                for (int x = 0; x < numcam; x++)     //construye lo que falta de la matriz de los camiones (lo que depende de la pala needy)

                {
                    /*
                     * sacar la fila bien
                     *
                     */

                    vectores.Ui[0] = vectores.Ui[0];

                    Din22[x, 4] = vectores.Fila[palaneedy - 1] * vectores.TCj[palaneedy - 1] + vectores.DesCam[x, 4];

                    if (Din22[x, 8] == 1)
                    {
                        Din22[x, 5] = vectores.TCV[Convert.ToInt32(vectores.DesCam[x, 2] - 1), palaneedy - 1];     //recoje tiempo de donde esta a la pala needy
                    }

                    Din22[x, 6] = vectores.PlTc[30, palaneedy - 1]; //tonelaje total asignado a palaneedy

                    if (vectores.Fila[palaneedy - 1] == 0)          //si no hay fila
                    {
                        if (Din22[x, 8] == 1)
                        {
                            Din22[x, 7] = vectores.TCV[Convert.ToInt32(vectores.DesCam[x, 2] - 1), palaneedy - 1];     //el shovel idle rate es Tij
                        }
                    }
                }


                double min2 = 90000000;
                camionpiola = 100;     //que sea mayor que la ultima linea
                double[] lost_tones = new double[200];
                newaux = 0;

                for (int x = 0; x < numcam; x++)     //calcula lost_tones para los importantes
                {
                    if (Din22[x, 8] == 1)
                    {
                        lost_tones[x] = ((Din22[x, 1] * Din22[x, 2] / Din22[x, 3]) * (Din22[x, 4] + Din22[x, 5]) + (Din22[x, 6] * Din22[x, 7]));
                        //lost_tones[x] = (Din22[x, 2] / Din22[x, 3]) * (Din22[x, 3] + Din22[x, 4]) + (Din22[x, 5] * Din22[x, 6]);

                        if (lost_tones[x] < min2)     //recupera el menor y actualiza el valor
                        {
                            min2 = lost_tones[x];

                            //camionpiola = Convert.ToInt32(Din22[x, 0]);

                            camionpiola = x + 1;
                        }
                    }

                    else

                    {
                        newaux = newaux + 1;
                    }
                }


                Din22[camionpiola - 1, 9]  = palaneedy; // asignar camino a camion
                Din22[camionpiola - 1, 8]  = 0;         //lo saca de la lista
                Din22[camionpiola - 1, 10] = caminoipalaneedy;

                for (int x = 0; x < numcam; x++)
                {
                    if (Din22[x, 8] == 1)
                    {
                        y = 1;
                    }
                }

                for (int x = 0; x < 30; x++)
                {
                    Din12[x + (palaneedy - 1) * 30, 4] = timenow;
                    Din12[x + (palaneedy - 1) * 30, 5] = Din12[x + (palaneedy - 1) * 30, 5] + 1;
                }
                newaux = indicemin1;
//                }

                /*else
                 * {
                 *  min = 959;
                 *  int palapiola = 100; //que sea mayor que la ultima linea
                 *  double min2 = 90000000;
                 *  for (int palanedy = 1; palanedy < 20; palanedy++)
                 *  {
                 *
                 *      if (vectores.PlTc[30, palanedy - 1] > 0)
                 *      {
                 *
                 *          Din22[Convert.ToInt32(aux2 - 1), 4] = vectores.Fila[palanedy - 1] * vectores.TCj[palanedy - 1] + vectores.DesCam[Convert.ToInt32(aux2 - 1), 4];
                 *
                 *          Din22[Convert.ToInt32(aux2 - 1), 5] = vectores.TCV[Convert.ToInt32(vectores.DesCam[Convert.ToInt32(aux2 - 1), 2] - 1), palanedy - 1]; //recoje tiempo de donde esta a la pala needy
                 *
                 *          Din22[Convert.ToInt32(aux2 - 1), 6] = vectores.PlTc[30, palanedy - 1]; //tonelaje total asignado a palaneedy
                 *
                 *          Din22[Convert.ToInt32(aux2 - 1), 7] = 0;
                 *
                 *          if (vectores.Fila[palanedy - 1] == 0) //si no hay fila
                 *          {
                 *              Din22[Convert.ToInt32(aux2 - 1), 7] = vectores.TCV[Convert.ToInt32(vectores.DesCam[Convert.ToInt32(aux2 - 1), 2] - 1), palanedy - 1]; //el shovel idle rate es Tij                            }
                 *          }
                 *
                 *          double lost_tones;
                 *
                 *          lost_tones = ((Din22[Convert.ToInt32(aux2 - 1), 1] * Din22[Convert.ToInt32(aux2 - 1), 2] / Din22[Convert.ToInt32(aux2 - 1), 3]) * (Din22[Convert.ToInt32(aux2 - 1), 4] + Din22[Convert.ToInt32(aux2 - 1), 5]) + (Din22[Convert.ToInt32(aux2 - 1), 6] * Din22[Convert.ToInt32(aux2 - 1), 7]));
                 *          //lost_tones[x] = (Din22[x, 2] / Din22[x, 3]) * (Din22[x, 3] + Din22[x, 4]) + (Din22[x, 5] * Din22[x, 6]);
                 *
                 *          if (lost_tones < min2) //recupera el menor y actualiza el valor
                 *          {
                 *
                 *              min2 = lost_tones;
                 *
                 *              palapiola = palanedy;
                 *
                 *          }
                 *      }
                 *  }
                 *
                 *  double min32 = 0;
                 *
                 *  for (int l = 0; l < 30; l++)
                 *  {
                 *      double last_tones = 0;
                 *      if (vectores.PlTc[l, palapiola - 1] > 0)
                 *      {
                 *          last_tones = vectores.PlTc[l, palapiola - 1] / (Din12[(palapiola - 1) * 30 + l, 7] * 300 / timenow);
                 *          if (last_tones > min32)
                 *          {
                 *
                 *              min32 = last_tones;
                 *              Din22[Convert.ToInt32(aux2 - 1), 9] = palapiola;
                 *              Din22[Convert.ToInt32(aux2 - 1), 10] = l + 1;
                 *
                 *          }
                 *      }
                 *      else
                 *      {
                 *          if (min2 == 0)
                 *          {
                 *              Din22[Convert.ToInt32(aux2 - 1), 9] = 21;
                 *          }
                 *      }
                 *
                 *
                 *  }
                 *  y = 0;
                 * }*/
            }while (y == 1);

            //for (int x = 0; x < numcam; x++) //DesCam suelta el dato como importante
            //{
            //    vectores.DesCam[x, 3] = 0;
            //}

            _salida.StateValue = Din22[Convert.ToInt32(aux2 - 1), 9];
            if (Din22[Convert.ToInt32(aux2 - 1), 9] == 4 || Din22[Convert.ToInt32(aux2 - 1), 9] == 5 || Din22[Convert.ToInt32(aux2 - 1), 9] == 6 || Din22[Convert.ToInt32(aux2 - 1), 9] == 10 || Din22[Convert.ToInt32(aux2 - 1), 9] == 17 || Din22[Convert.ToInt32(aux2 - 1), 9] == 18)
            {
                Din22[Convert.ToInt32(aux2 - 1), 10] = 1;
            }

            else if (Din22[Convert.ToInt32(aux2 - 1), 9] == 7 || Din22[Convert.ToInt32(aux2 - 1), 9] == 8 || Din22[Convert.ToInt32(aux2 - 1), 9] == 12)
            {
                Din22[Convert.ToInt32(aux2 - 1), 10] = 4;
            }

            else if (Din22[Convert.ToInt32(aux2 - 1), 9] == 1 || Din22[Convert.ToInt32(aux2 - 1), 9] == 2 || Din22[Convert.ToInt32(aux2 - 1), 9] == 3 || Din22[Convert.ToInt32(aux2 - 1), 9] == 11)
            {
                Din22[Convert.ToInt32(aux2 - 1), 10] = 8;
            }

            _aux.StateValue = Din22[Convert.ToInt32(aux2 - 1), 10];                                                                                   //
            //_aux2.StateValue = camionpiola;
            _needtime.StateValue  = min;                                                                                                              //
            _losttones.StateValue = ((Din12[33, 0] / Din12[33, 1]) * (Din12[33, 5] * 300 - Din12[33, 1] * timenow)) / (Din12[33, 7] * 300 / timenow); //(Din12[33, 7] * 300 / timenow);//
            _needch.StateValue    = ((Din12[37, 0] / Din12[37, 1]) * (Din12[37, 5] * 300 - Din12[37, 1] * timenow)) / (Din12[37, 7] * 300 / timenow); //

            sr.serializa(vectores);
            sr.closeStream();
            return(ExitType.FirstExit);
        }