コード例 #1
0
        //Итеративный расчет содержания PO (общая функция)
        private (double, double[]) CalculateStrength(double waterContent, double acnContent, float densityFromMassFT, ref Density inputDensity, bool isIncrement)
        {
            //isDecrement = true - итерация в сторону увеличения; false - итерация в сторону уменьшения
            var newDens   = isIncrement ? 0.0 : 10000.0;
            var POContent = 0.0;
            var percArray = new double[4] {
                waterContent, acnContent, 0.0, 0.0
            };                                                                    //ACN, Water, P, PO

            inputDensity.PercArray = percArray;

            if (inputDensity == null)
            {
                return(-1.0, new double[] { -1.0, -1.0, -1.0, -1.0 });
            }

            int i = 0;

            //Формирование условия выхода из цикла в зависимости от направления итерации
            bool GetExitCondition()
            {
                return(isIncrement ? newDens > densityFromMassFT * 10.0 : newDens < densityFromMassFT * 10.0);
            }

            while (true)
            {
                //newDens = inputDensity.CalculateDensity();
                if (GetExitCondition() || i++ > 10000)
                {
                    POContent = inputDensity.PercArray[3];
                    break;
                }
                inputDensity.PercArray[3] += 0.01;                                                                                       //PO
                inputDensity.PercArray[1]  = (100.0 - inputDensity.PercArray[3] - inputDensity.PercArray[2]) * acnContent * 0.01;        //ACN
                inputDensity.PercArray[0]  = 100.0 - inputDensity.PercArray[1] - inputDensity.PercArray[2] - inputDensity.PercArray[3];  //Water
                newDens = inputDensity.CalculateDensity();
            }
            return(Math.Min(100.0, POContent), inputDensity.PercArray);
        }
コード例 #2
0
        internal void CalculateMassOfPropylene()
        {
            double propyleneMass;

            if (!isInitPropyleneSuccess)
            {
                return;
            }

            //Проверка на наличие минимального расхода ACN и воды, а также, чтобы расход из P05 был выше, чем из D02
            if ((singleTagCreator.AcnWaterMassFlow < 408.85 * densityCreator.DensityList.FirstOrDefault(d => d.TagName == "S11_A01_FC02_DENS").ValHmi * 0.0001) ||
                singleTagCreator.averReactFlow < singleTagCreator.AcnWaterMassFlow / (densityCreator.DensityList.FirstOrDefault(d => d.TagName == "S11_A01_FC02_DENS").ValHmi * 0.0001) + 1)
            {
                propyleneMass = 0.0;
            }
            else
            {
                var waterMass = singleTagCreator.AcnWaterMassFlow * acnWaterDensity.PercArray[0] * 0.01;
                var acnMass   = singleTagCreator.AcnWaterMassFlow * acnWaterDensity.PercArray[1] * 0.01;

                //Инициализурем расчитываемую плотность стартовыми значениями (без пропилена - только вода и ACN
                tempDensity.PercArray[0] = acnWaterDensity.PercArray[0];
                tempDensity.PercArray[1] = acnWaterDensity.PercArray[1];
                tempDensity.PercArray[2] = 0.0;

                var dens = tempDensity.CalculateDensity();

                var massAverReactFlowSave = singleTagCreator.averReactFlow * dens * 0.0001;
                var massAverReactFlow     = 0.0;
                var diff = 0.0;
                //Подставляем для расчета плотности температуру после P05
                tempDensity.Temperature = acnWaterPropyleneTemperature;
                var i = 0;
                do
                {   //Рассчитываем массовый расход реакционной смеси и % компонентов в ней
                    if (massAverReactFlowSave != 0)
                    {
                        tempDensity.PercArray[0] = Math.Max(0.0, Math.Min(100.0, waterMass * 100.0 / massAverReactFlowSave));
                        tempDensity.PercArray[1] = Math.Max(0.0, Math.Min(100.0, acnMass * 100.0 / massAverReactFlowSave));
                        tempDensity.PercArray[2] = Math.Max(0.0, 100.0 - tempDensity.PercArray[1] - tempDensity.PercArray[0]);
                    }


                    //Считаем плотность с новым содержанием компонентов
                    dens = tempDensity.CalculateDensity();
                    if (dens != -1)
                    {
                        massAverReactFlow = singleTagCreator.averReactFlow * dens * 0.0001;
                        propyleneMass     = massAverReactFlow * tempDensity.PercArray[2] * 0.01;
                        diff = Math.Abs(massAverReactFlowSave - massAverReactFlow);
                        massAverReactFlowSave = massAverReactFlow;
                        i++;
                    }
                    else
                    {
                        return;
                    }
                }while (diff > 1.0 || i < 10);
            }

            singleTagCreator.PropyleneMass = (short)(propyleneMass * 10.0);
        }