Merge branch 'master' of http://47.103.154.90:83/r/HStation/XHS.V1.0

lixiaojun

2024-11-08 a1fe0840ce2022ca711994ba10ba3512ae3c43be

 WinFrmUI/Yw.WinFrmUI.Phart.Core/00-core/01-calculate/PumpParallelConnectionHelper.cs

@@ -8,8 +8,8 @@
    public class PumpParallelConnectionHelper
    { 
       
        private List<Yw.Geometry.CubicSpline2d> _cubic_spline_qh_list = new List<Yw.Geometry.CubicSpline2d>();
        private List<Yw.Geometry.CubicSpline2d> _cubic_spline_qp_list = new List<Yw.Geometry.CubicSpline2d>();
        private List<Yw.Geometry.CubicSpline2d> _qh_list = new List<Yw.Geometry.CubicSpline2d>();
        private List<Yw.Geometry.CubicSpline2d> _qp_list = new List<Yw.Geometry.CubicSpline2d>();

        // 延长率, 默认不延长
        private double _extend_min_ratio = 1;
@@ -26,29 +26,29 @@
        /// <summary>
        /// 添加
        /// </summary>
        /// <param name="cubic_spline_qh"></param>
        /// <param name="cubic_spline_qp"></param>
        public void Add(Yw.Geometry.CubicSpline2d cubic_spline_qh, Yw.Geometry.CubicSpline2d cubic_spline_qp)
        /// <param name="qh"></param>
        /// <param name="qp"></param>
        public void Add(Yw.Geometry.CubicSpline2d qh, Yw.Geometry.CubicSpline2d qp)
        {
            if (cubic_spline_qh != null)
                _cubic_spline_qh_list.Add(cubic_spline_qh);
            if (qh != null)
                _qh_list.Add(qh);

            if (cubic_spline_qp != null)
                _cubic_spline_qp_list.Add(cubic_spline_qp);
            if (qp != null)
                _qp_list.Add(qp);
        }

        /// <summary>
        /// 添加
        /// </summary>
        /// <param name="pt_qh_list"></param>
        /// <param name="pt_qp_list"></param>
        public void Add(List<Yw.Geometry.Point2d> pt_qh_list, List<Yw.Geometry.Point2d> pt_qp_list)
        /// <param name="qh_pt_list"></param>
        /// <param name="qp_pt_list"></param>
        public void Add(List<Yw.Geometry.Point2d> qh_pt_list, List<Yw.Geometry.Point2d> qp_pt_list)
        {
            if (pt_qh_list != null && pt_qh_list.Count > 3)
                _cubic_spline_qh_list.Add(new Yw.Geometry.CubicSpline2d(pt_qh_list));
            if (qh_pt_list != null && qh_pt_list.Count > 3)
                _qh_list.Add(new Yw.Geometry.CubicSpline2d(qh_pt_list));

            if (pt_qp_list != null && pt_qp_list.Count > 3)
                _cubic_spline_qp_list.Add(new Yw.Geometry.CubicSpline2d(pt_qp_list));
            if (qp_pt_list != null && qp_pt_list.Count > 3)
                _qp_list.Add(new Yw.Geometry.CubicSpline2d(qp_pt_list));
        }

        /// <summary>
@@ -56,41 +56,41 @@
        /// </summary>
        public void Clear()
        {
            _cubic_spline_qh_list.Clear();
            _cubic_spline_qp_list.Clear();
            _qh_list.Clear();
            _qp_list.Clear();
        }
        #endregion

        /// <summary>
        /// 计算并联
        /// </summary>
        /// <param name="pt_qh_list"></param>
        /// <param name="pt_qe_list"></param>
        /// <param name="pt_qp_list"></param>
        /// <param name="qh_pt_list"></param>
        /// <param name="qe_pt_list"></param>
        /// <param name="qp_pt_list"></param>
        /// <returns></returns>
        public bool CalculateParallel(out List<Yw.Geometry.Point2d> pt_qh_list, out List<Yw.Geometry.Point2d> pt_qe_list, out List<Yw.Geometry.Point2d> pt_qp_list)
        public bool CalculateParallel(out List<Yw.Geometry.Point2d> qh_pt_list, out List<Yw.Geometry.Point2d> qe_pt_list, out List<Yw.Geometry.Point2d> qp_pt_list)
        {
            pt_qh_list = pt_qe_list = pt_qp_list = null;
            qh_pt_list = qe_pt_list = qp_pt_list = null;
            //求H的最大最小值
            if (!CalculateMaxMinH())
                return false;

            //计算并联
            return GetCurveParallel(out pt_qh_list, out pt_qe_list, out pt_qp_list);
            return GetCurveParallel(out qh_pt_list, out qe_pt_list, out qp_pt_list);
        }

        /// <summary>
        /// 从并联曲线中分割出单泵曲线
        /// </summary>
        /// <param name="connect_pt_qh_list"></param>
        /// <param name="separate_pt_qh_list"></param>
        /// <param name="connect_qh_pt_list"></param>
        /// <param name="separate_qh_pt_list"></param>
        /// <returns></returns>
        public bool SeparateParallelQH(List<Yw.Geometry.Point2d> connect_pt_qh_list, out List<Yw.Geometry.Point2d> separate_pt_qh_list)
        public bool SeparateParallelQH(List<Yw.Geometry.Point2d> connect_qh_pt_list, out List<Yw.Geometry.Point2d> separate_qh_pt_list)
        {
            separate_pt_qh_list = new List<Yw.Geometry.Point2d>();
            separate_qh_pt_list = new List<Yw.Geometry.Point2d>();

            _max_head = (from x in connect_pt_qh_list select x.Y).Max();
            _min_head = (from x in connect_pt_qh_list select x.Y).Min();
            _max_head = (from x in connect_qh_pt_list select x.Y).Max();
            _min_head = (from x in connect_qh_pt_list select x.Y).Min();

            int insert_num = 30;
            double space = (_max_head - _min_head) / (insert_num - 1);
@@ -99,18 +99,18 @@
            for (j = 0; j < insert_num; j++)
            {
                //计算流量
                var q_list = connect_pt_qh_list.GetInterPointsX(head)?.ToList();
                var q_list = connect_qh_pt_list.GetInterPointsX(head)?.ToList();
                if (q_list == null || q_list.Count() == 0)
                {
                    continue;
                }
                var total_q = q_list.Last();
                var is_insect = true;
                for (i = 0; i < _cubic_spline_qh_list.Count; i++)
                for (i = 0; i < _qh_list.Count; i++)
                {
                    var cubic_spline_qh = _cubic_spline_qh_list[i];
                    var qh = _qh_list[i];
                    //计算流量
                    var q_list2 = cubic_spline_qh.GetPointsX(head)?.ToList();
                    var q_list2 = qh.GetPointsX(head)?.ToList();
                    if (q_list2 == null || q_list2.Count == 0)
                    {
                        is_insect = false;
@@ -123,44 +123,44 @@
                //添加到列表
                if (is_insect)
                {
                    separate_pt_qh_list.Add(new Yw.Geometry.Point2d(total_q, head));
                    separate_qh_pt_list.Add(new Yw.Geometry.Point2d(total_q, head));
                }

                head = head + space;
            }

            if (separate_pt_qh_list.Count < 3)
            if (separate_qh_pt_list.Count < 3)
            {
                separate_pt_qh_list = null;
                separate_qh_pt_list = null;
                return false;
            }
             
            separate_pt_qh_list = new Yw.Geometry.CubicSpline2d(separate_pt_qh_list).GetPointList(20);
            separate_qh_pt_list = new Yw.Geometry.CubicSpline2d(separate_qh_pt_list).GetPointList(20);
            return true;
        }


        //计算串联
        public bool CalculateSeries(
            out List<Yw.Geometry.Point2d> pt_qh_list,
            out List<Yw.Geometry.Point2d> pt_qe_list,
            out List<Yw.Geometry.Point2d> pt_qp_list)
            out List<Yw.Geometry.Point2d> qh_pt_list,
            out List<Yw.Geometry.Point2d> qe_pt_list,
            out List<Yw.Geometry.Point2d> qp_pt_list)
        {
            pt_qh_list = pt_qe_list = pt_qp_list = null;
            qh_pt_list = qe_pt_list = qp_pt_list = null;
            //求Q的最大最小值
            if (!CalculateMaxMinQ())
                return false;

            //计算并联
            return GetCurveSeries(out pt_qh_list, out pt_qe_list, out pt_qp_list);
            return GetCurveSeries(out qh_pt_list, out qe_pt_list, out qp_pt_list);
        }
         

        private bool CalculateMaxMinH()
        {
            if (_cubic_spline_qh_list == null || _cubic_spline_qh_list.Count() == 0)
            if (_qh_list == null || _qh_list.Count() == 0)
                return false;
            foreach (var curves in _cubic_spline_qh_list)
            foreach (var curves in _qh_list)
            {
                if (curves == null)
                    continue;
@@ -174,9 +174,9 @@

        private bool CalculateMaxMinQ()
        {
            if (_cubic_spline_qh_list == null || _cubic_spline_qh_list.Count() == 0)
            if (_qh_list == null || _qh_list.Count() == 0)
                return false;
            foreach (var cubic_spline in _cubic_spline_qh_list)
            foreach (var cubic_spline in _qh_list)
            {
                if (cubic_spline == null)
                    continue;
@@ -192,19 +192,19 @@
        }

        //计算并联
        private bool GetCurveParallel(out List<Yw.Geometry.Point2d> pt_qh_list, out List<Yw.Geometry.Point2d> pt_qe_list, out List<Yw.Geometry.Point2d> pt_qp_list)
        private bool GetCurveParallel(out List<Yw.Geometry.Point2d> qh_pt_list, out List<Yw.Geometry.Point2d> qe_pt_list, out List<Yw.Geometry.Point2d> qp_pt_list)
        {
            if (_cubic_spline_qh_list == null || !_cubic_spline_qh_list.Any())
            if (_qh_list == null || !_qh_list.Any())
            {
                pt_qh_list = null;
                pt_qe_list = null;
                pt_qp_list = null;
                qh_pt_list = null;
                qe_pt_list = null;
                qp_pt_list = null;
                return false;
            }

            pt_qh_list = new List<Yw.Geometry.Point2d>();
            pt_qe_list = new List<Yw.Geometry.Point2d>();
            pt_qp_list = new List<Yw.Geometry.Point2d>();
            qh_pt_list = new List<Yw.Geometry.Point2d>();
            qe_pt_list = new List<Yw.Geometry.Point2d>();
            qp_pt_list = new List<Yw.Geometry.Point2d>();

            int insert_num = 30;
            double space = (_max_head - _min_head) / (insert_num - 1);
@@ -216,15 +216,15 @@
                double total_q = 0;
                double total_p = 0;
                bool is_insect = true;
                for (i = 0; i < _cubic_spline_qh_list.Count(); i++)
                for (i = 0; i < _qh_list.Count(); i++)
                {
                    Yw.Geometry.CubicSpline2d cubic_spline_qh = _cubic_spline_qh_list[i];
                    Yw.Geometry.CubicSpline2d cubic_spline_qp = null;
                    if (_cubic_spline_qp_list != null && _cubic_spline_qp_list.Count > i)
                        cubic_spline_qp = _cubic_spline_qp_list[i];
                    Yw.Geometry.CubicSpline2d qh = _qh_list[i];
                    Yw.Geometry.CubicSpline2d qp = null;
                    if (_qp_list != null && _qp_list.Count > i)
                        qp = _qp_list[i];

                    //计算流量
                    var listQ = cubic_spline_qh.GetPointsX(head);
                    var listQ = qh.GetPointsX(head);
                    if (listQ == null || listQ.Length == 0)
                    {
                        is_insect = false;
@@ -234,9 +234,9 @@
                    total_q = total_q + Q;

                    //计算功率
                    if (cubic_spline_qp != null)
                    if (qp != null)
                    {
                        double p = cubic_spline_qp.GetPointY(Q);
                        double p = qp.GetPointY(Q);
                        total_p = total_p + p;
                    }
                }
@@ -244,47 +244,47 @@
                //添加到列表
                if (is_insect)
                {
                    pt_qh_list.Add(new Yw.Geometry.Point2d(total_q, head));
                    qh_pt_list.Add(new Yw.Geometry.Point2d(total_q, head));
                    if (total_p > 0.001)
                    {
                        pt_qp_list.Add(new Yw.Geometry.Point2d(total_q, total_p));
                        qp_pt_list.Add(new Yw.Geometry.Point2d(total_q, total_p));

                        double E = PumpCalcHelper.CalculateE(total_q, head, total_p);
                        pt_qe_list.Add(new Yw.Geometry.Point2d(total_q, E));
                        qe_pt_list.Add(new Yw.Geometry.Point2d(total_q, E));
                    }
                }

                head = head + space;
            }

            if (!pt_qh_list.Any())
            if (!qh_pt_list.Any())
            {
                return false;
            }
            //根据流量排序
            Point2dComparer.Sort featPointComp = new Point2dComparer.Sort(Point2dComparer.eSortType.X);
            pt_qh_list.Sort(featPointComp);
            if (pt_qe_list != null && pt_qe_list.Count > 2)
                pt_qe_list.Sort(featPointComp);
            if (pt_qp_list != null && pt_qp_list.Count > 2)
                pt_qp_list.Sort(featPointComp);
            qh_pt_list.Sort(featPointComp);
            if (qe_pt_list != null && qe_pt_list.Count > 2)
                qe_pt_list.Sort(featPointComp);
            if (qp_pt_list != null && qp_pt_list.Count > 2)
                qp_pt_list.Sort(featPointComp);

            return true;
        }

        //计算串联
        private bool GetCurveSeries(out List<Yw.Geometry.Point2d> pt_qh_list, out List<Yw.Geometry.Point2d> pt_qe_list, out List<Yw.Geometry.Point2d> pt_qp_list)
        private bool GetCurveSeries(out List<Yw.Geometry.Point2d> qh_pt_list, out List<Yw.Geometry.Point2d> qe_pt_list, out List<Yw.Geometry.Point2d> qp_pt_list)
        {
            if (_cubic_spline_qh_list.Count() == 0)
            if (_qh_list.Count() == 0)
            {
                pt_qh_list = null;
                pt_qe_list = null;
                pt_qp_list = null;
                qh_pt_list = null;
                qe_pt_list = null;
                qp_pt_list = null;
                return false;
            }
            pt_qh_list = new List<Yw.Geometry.Point2d>();
            pt_qe_list = new List<Yw.Geometry.Point2d>();
            pt_qp_list = new List<Yw.Geometry.Point2d>();
            qh_pt_list = new List<Yw.Geometry.Point2d>();
            qe_pt_list = new List<Yw.Geometry.Point2d>();
            qp_pt_list = new List<Yw.Geometry.Point2d>();

            int insert_num = 30;
            double space = (_max_flow - _min_flow) / (insert_num - 1);
@@ -296,21 +296,21 @@
                double total_h = 0;
                double total_p = 0;
                bool isInsect = true;
                for (i = 0; i < _cubic_spline_qh_list.Count(); i++)
                for (i = 0; i < _qh_list.Count(); i++)
                {
                    Yw.Geometry.CubicSpline2d cubic_spline_qh = _cubic_spline_qh_list[i];
                    Yw.Geometry.CubicSpline2d cubic_spline_qp = null;
                    if (_cubic_spline_qp_list != null && _cubic_spline_qp_list.Count > i)
                        cubic_spline_qp = _cubic_spline_qp_list[i];
                    Yw.Geometry.CubicSpline2d qh = _qh_list[i];
                    Yw.Geometry.CubicSpline2d qp = null;
                    if (_qp_list != null && _qp_list.Count > i)
                        qp = _qp_list[i];

                    //计算H
                    double h = cubic_spline_qh.GetPointY(Q);
                    double h = qh.GetPointY(Q);
                    total_h = total_h + h;

                    //计算功率
                    if (cubic_spline_qp != null)
                    if (qp != null)
                    {
                        double p = cubic_spline_qp.GetPointY(Q);
                        double p = qp.GetPointY(Q);
                        total_p = total_p + p;
                    }
                }
@@ -318,13 +318,13 @@
                //添加到列表
                if (isInsect)
                {
                    pt_qh_list.Add(new Yw.Geometry.Point2d(Q, total_h));
                    qh_pt_list.Add(new Yw.Geometry.Point2d(Q, total_h));
                    if (total_p > 0.001)
                    {
                        pt_qp_list.Add(new Yw.Geometry.Point2d(Q, total_p));
                        qp_pt_list.Add(new Yw.Geometry.Point2d(Q, total_p));

                        double E = PumpCalcHelper.CalculateE(Q, total_h, total_p);
                        pt_qe_list.Add(new Yw.Geometry.Point2d(Q, E));
                        qe_pt_list.Add(new Yw.Geometry.Point2d(Q, E));
                    }
                }

@@ -333,11 +333,11 @@

            //根据流量排序
            Point2dComparer.Sort featPointComp = new Point2dComparer.Sort(Point2dComparer.eSortType.X);
            pt_qh_list.Sort(featPointComp);
            if (pt_qe_list != null && pt_qe_list.Count > 2)
                pt_qe_list.Sort(featPointComp);
            if (pt_qp_list != null && pt_qp_list.Count > 2)
                pt_qp_list.Sort(featPointComp);
            qh_pt_list.Sort(featPointComp);
            if (qe_pt_list != null && qe_pt_list.Count > 2)
                qe_pt_list.Sort(featPointComp);
            if (qp_pt_list != null && qp_pt_list.Count > 2)
                qp_pt_list.Sort(featPointComp);

            return true;
        }