Example #1
0
        public string[] Calculation(string sql)
        {
            using (MySqlConnection conn = new MySqlConnection(connString))
            {
                MySqlCommand cmd = new MySqlCommand(sql, conn);
                conn.Open();
                MySqlDataAdapter MydataAdapter = new MySqlDataAdapter(cmd);

                MydataAdapter.TableMappings.Add("Table", "typicalprofile");
                MydataAdapter.TableMappings.Add("Table1", "enginelifecon");
                //MydataAdapter.TableMappings.Add("Table2", "failurerateconInformation");
                DataSet ds = new DataSet();
                MydataAdapter.Fill(ds);//将典型剖面表、装备任务表填充


                int K = 9;                                      //飞行阶段列数

                List <double> taskfp = new List <double>();     //典型剖面飞行任务阶段列表Flight phase,包括地面慢车、有地效悬停、起飞增速	、斜向爬升、上升转弯、巡航、盘旋飞行、下滑、着陆

                List <double> typicalfp  = new List <double>(); //任务剖面飞行任务阶段Flight phase,包括地面慢车、有地效悬停、起飞增速	、斜向爬升、上升转弯、巡航、盘旋飞行、下滑、着陆
                List <double> typicalttt = new List <double>(); //典型剖面任务任务总时间列
                List <double> typicaldl  = new List <double>(); ///典型剖面设计寿命列DesignLife
                List <string> ENGINE     = new List <string>(); //任务表中发动机型号
                List <string> PN         = new List <string>(); //任务表中基准剖面名称
                List <string> MTID       = new List <string>(); //任务表中任务标识
                List <double> SDL        = new List <double>(); //任务剖面设计寿命列
                List <string> ep         = new List <string>(); // 环境参数列表,包括地域、季节、海拔、其他4列
                List <string> TEMP       = new List <string>(); //温度参数
                List <string> HUMP       = new List <string>(); //湿度参数


                foreach (DataRow row in ds.Tables["enginelifecon"].Rows)//将任务表中的飞行阶段的列添加到列表
                {
                    for (int i = 2; i < 11; i++)
                    {
                        if (row[i] == System.DBNull.Value)               //判断任务飞行阶段是否空
                        {
                            ds.Tables["enginelifecon"].Rows[0][i] = "0"; //如果为空,飞行阶段时间设为0
                            MySqlCommandBuilder MysqlCommandBuilder = new MySqlCommandBuilder(MydataAdapter);
                            MydataAdapter.Update(ds);
                        }
                        taskfp.Add(double.Parse(row[i].ToString()));
                    }
                }

                //任务剖面飞行任务阶段时间矩阵Ts的创建
                double[] task = taskfp.ToArray(); //将list转化为一维数组
                double[,] Ts = new double[K, 1];  //将一维数组task转化为矩阵Ts
                for (int i = 0; i < K; i++)
                {
                    Ts[i, 0] = task[i];
                }



                foreach (DataRow row1 in ds.Tables["typicalprofile"].Rows)//将典型剖面表中的飞行阶段的列添加到列表typicalfp
                {
                    for (int i = 2; i < 11; i++)
                    {
                        typicalfp.Add(double.Parse(row1[i].ToString()));
                    }
                }

                //将典型剖面飞行阶段列表typicalf转化为一维数组typical
                double[] typical = typicalfp.ToArray();
                int      N       = typical.Length / K; //典型剖面表的行数
                double[,] P = new double[N, K];        //将一维数组typical转化为矩阵P
                for (int i = 0; i < N; i++)
                {
                    for (int j = 0; j < K; j++)
                    {
                        P[i, j] = typical[i * (K) + j];
                    }
                }


                //提取出典型剖面任务总时间的列添加到列表typicalttt
                foreach (DataRow row2 in ds.Tables["typicalprofile"].Rows)
                {
                    typicalttt.Add(double.Parse(row2["TTD"].ToString()));
                }

                double[] ttt = typicalttt.ToArray(); //将典型剖面中任务总时间列表转化为一维数组ttt
                double[,] TTT = new double[N, 1];    //将一维数组ttt转化为矩阵TTT
                for (int i = 0; i < N; i++)
                {
                    TTT[i, 0] = ttt[i];
                }

                string    tmpcmd       = "SELECT TEMP,HUMP FROM tempandhump WHERE MAREA LIKE '%" + ds.Tables["enginelifecon"].Rows[0]["MAREA"].ToString() + "%'AND MSEA='" + ds.Tables["enginelifecon"].Rows[0]["MSEA"].ToString() + "';";
                DataTable tmpdatatable = FillTable(tmpcmd);
                TEMP.Add(tmpdatatable.Rows[0][0].ToString());
                HUMP.Add(tmpdatatable.Rows[0][1].ToString());

                //典型剖面设计寿命列
                foreach (DataRow row3 in ds.Tables["typicalprofile"].Rows)//将典型剖面中设计寿命列添加到typicaldl列表
                {
                    typicaldl.Add(double.Parse(row3["DLF"].ToString()));
                }

                double[]  DL             = typicaldl.ToArray();//将typicaldl列表转化为一维数组 DL
                string    tmpcmdDL       = "SELECT DLF FROM typicalprofile WHERE TPNAME ='" + ds.Tables["enginelifecon"].Rows[0]["DAP"] + "' ";
                DataTable tmpdatatableDL = FillTable(tmpcmdDL);
                SDL.Add(double.Parse(tmpdatatableDL.Rows[0][0].ToString()));

                foreach (DataRow row5 in ds.Tables["enginelifecon"].Rows) // //提取出任务表中地域、季节、海拔、其他4列添加到列表ep
                {
                    for (int i = 13; i < 17; i++)
                    {
                        ep.Add(row5[i].ToString());
                    }
                }
                foreach (DataRow row6 in ds.Tables["enginelifecon"].Rows)//提取发动机型号列添加到列表ENGINE
                {
                    ENGINE.Add(row6[0].ToString());
                }
                foreach (DataRow row7 in ds.Tables["enginelifecon"].Rows)//提取任务表中剖面名称列添加到列表PN
                {
                    PN.Add(row7[12].ToString());
                }
                foreach (DataRow row8 in ds.Tables["enginelifecon"].Rows) //提取任务表中任务标识列添加到列表MTID
                {
                    MTID.Add(row8[1].ToString());
                }


                conn.Close();

                double   E  = 1;                               //初始环境系数
                double[] EP = new double[4];                   //温度,湿度,海拔,人因影响系数4个环境参数系数

                EP[0] = Calculater.TCC(double.Parse(TEMP[0])); //通过温度参数值计算得到温度系数
                EP[1] = Calculater.HCC(double.Parse(HUMP[0])); //通过湿度参数值计算得到湿度系数
                EP[2] = Calculater.ACC(double.Parse(ep[2]));   //在环境参数列表ep中提出海拔参数值并计算得到海拔系数
                EP[3] = double.Parse(ep[3]);                   //在环境参数列表ep中提出人因影响系数



                string[] str1 = new string[5];         //建立一个字符串数组,用于接收输出结果,其中第四个元素接收计算结果
                str1[0] = ENGINE[0];                   //将发动机型号填充到数组
                str1[1] = PN[0];                       //将基准剖面名称填充到数组
                str1[2] = MTID[0];                     //将任务标识填充到数组
                str1[4] = Convert.ToString(SDL[0]);    //将任务表中设计寿命填充到数组

                double TaskLifeConsumptionResults = 0; //开始计算前寿命消耗设为0


                double[,] TransP = Calculater.Transposition(P);                                                               //将P转置
                double[,] PTP    = Calculater.Multiplication(P.GetLength(0), P.GetLength(1), TransP.GetLength(1), P, TransP); //P和其转置的乘积

                double[,] PMI = Calculater.MatrixInversion(PTP);                                                              //将PTP求逆



                if (PMI == null)//如果矩阵PMI为空即矩阵PTP的逆不存在,否则就计算寿命消耗值
                {
                    MessageBox.Show("求逆错误!寿命消耗折算无法正常计算!");

                    str1[3] = "-99999"; //将寿命计算结果设为 - 99999,说明寿命消耗不能计算;
                    return(str1);
                }
                else
                {
                    double[,] MUL = Calculater.Multiplication(PMI.GetLength(0), PMI.GetLength(1), P.GetLength(1), PMI, P);   //PMI与P的乘积
                    double[,] L   = Calculater.Multiplication(MUL.GetLength(0), MUL.GetLength(1), Ts.GetLength(1), MUL, Ts); //剖面合成权重系数L

                    //飞行剖面折合系数
                    double[,] ccofp = Calculater.CCOFP(DL);

                    //典型剖面第一个剖面为基准剖面时寿命消耗折算结果
                    double CalculationResults = Calculater.LifeConsumptionConversion(TTT, L, ccofp, E, EP);
                    //任务表中基准剖面下寿命消耗折算结果
                    TaskLifeConsumptionResults = (typicaldl[0] * CalculationResults / SDL[0]);
                    str1[3] = TaskLifeConsumptionResults.ToString("0.########");
                }


                return(str1);
            }
        }