• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Optimizing reliability,maintainability and testability parameters of equipment based on GSPN

    2015-11-14 06:52:28TingpengLiYueLiYanlingQianandYongchengXu

    Tingpeng Li,Yue Li,Yanling Qian,and Yongcheng Xu

    Science and Technology on Integrated Logistics Support Laboratory,National University of Defense Technology,Changsha 410073,China

    1.Introduction

    Reliability,maintainability and testability(RMT)are important design parameters of equipment.RMT parameters affect not only the property and capability of the equipment,but also the life cycle costs(LCC).

    In recent years,some progress has been made on determination of RMT parameters.Liu[1]proposed a method to determinate reliability indexes,which integrates experience of similar equipment,requirements of users and capability of suppliers.In[2],a method of determination of reliability parameters based on the constraint of LCC was presented.Quality function deployment(QFD)was introduced in[3]to determine maintainability parameters.Li[4]established a maintainability evaluation model to evaluate schemes of maintainability based on grey relational analysis.With the results of evaluation and repeated attempts,designers can get proper testability parameters.The genetic algorithm(GA)was proposed to determine system testability indexes in[5,6].Zhang[7]proposed a method to select testability parameters on the basis of modified analytic hierarchy process,whereas this method did not discuss how to determine the values of parameters.Qian[8]introduced generalized stochastic Petri nets(GSPN)to define system-level testability parameters.However,he only established the relationship between operational availability and testability parameters without other constraints(e.g.cost constraint)considered.On the other hand,only testability parameters can be obtained through repeated analysis on the relation curves of availability and testability parameters.Ding[9]proposed a system of RMT parameters which can be divided into three levels:tactics,techniques and engineering design.However,the method for how to determine parameters was not discussed.

    Jonge presented a method to determine the maintenance strategy based on lifetime distribution,which is assumed to obey Weibull distribution[10].Saranga studied the application of probabilistic maintenance models to determine the optimal inspection rates considering the tradeoff between reliability and cost[11].Tian developed a physical programming based approach to tradeoff between cost and reliability[12].Kim estimated the condition-based maintenance parameter based on the continuous time Markov process and the EM algorithm[13].Jeang optimized the design and maintenance parameter(interval time)considering LCC of the product[14].

    However,most of these existing works mentioned above only focus on one of the RMT properties(e.g.[1–6]),and very few consider all of them synthetically.Besides,though the method of selection of RMT parameters is introduced in some papers,there is no discussion on how to determine the values of the parameters(e.g.[7,9]).Moreover,most of the existing methods are based on the history data of similar equipment or experts’experiences,but seldom consider the effects of RMT on the operation process(e.g.[3,5]).

    In this paper,a new method is proposed to determine RMT parameters.This method considers not only the history data of similar equipment and experts’experiences,but also the effects of RMT on the operation process.

    2.Problem description

    RMT parameters are proposed in the design phase,and have significant effects on the operation process of equipment.For example,mean time between faults(MTBF),an important reliability parameter,affects the continuous working time without faults.Therefore,this paper takes the operation process of equipment into account to determine RMT parameters.

    Generally,equipment can be categorized into two types:continuous working equipment and intermittent working equipment.This paper focuses mainly on continuous working equipment whose work time is much longer than MTBF,such as communication equipment(e.g.radar).In addition,the equipment,whose down time without faults is very short,could also be viewed as continuous working equipment.

    2.1 The operation process of equipment

    Continuous working equipment has only two states:working and under maintenance.Maintenance falls into preventative maintenance and corrective maintenance.

    Preventive maintenance consists of on-condition maintenance and scheduled maintenance,and corrective maintenance consists of accurate maintenance and fuzzy maintenance.Accurate maintenance means motivated maintenance with successful results of faults detection and diagnosis,while fuzzy maintenance means exploratory maintenance when faults detection or diagnosis fails.

    The operation process of continuous working equipment is shown in Fig.1.The equipment works normally for most of the time.However,whenever faults occur,the equipment has to be stopped to be detected and diagnosed.If faults detection and diagnosis are successful,the equipment will recover to work immediately after accurate maintenance.Otherwise,fuzzy maintenance must be executed.In addition,on-condition maintenance will be conducted when the states of the equipment satisfy the conditions of preventive maintenance.

    Fig.1 Operation process of continuous working equipment

    2.2 RMT parameters of equipment

    Reliability is the ability to complete required functions under specified conditions and within stipulated time[15].Equipment reliability consists of basic reliability and mission reliability[16].Basic reliability stands for the failurefree probability under specific conditions within a period of time.Mission reliability is the ability to complete prescribed tasks.In the design phase,designers mainly consider the basic reliability of equipment.Therefore,the most important parameter of basic reliability,i.e.MTBF,is the mainly considered reliability parameter in this paper.Generally,MTBF can be divided into preventive mean time between faults(PMTBF)and never preventive mean time between faults(NPMTBF).PMTBF is the mean time between faults which can be maintained preventively,and NPMTBF means the mean time between faults which can only be maintained after occurring.If the failure rate of equipment obeys the exponential distribution,the relationship among the three parameters is shown as follows:

    whereTBFis MTBF,TBFPis PMTBF,andTBFNPis NPMTBF.

    Maintainability means the ability to be maintained or repaired in a specified state according to prescribed procedures and methods[15].Mean time of accurate maintenance(TAM),mean time of fuzzy maintenance(TFM)and mean time of preventive maintenance(TPM)are main parameters of maintainability.

    Testability is a design characteristic which represents the ability to determine its states timely and accurately and to diagnose its faults effectively[17].This paper chooses four most important factors as testability parameters:faults detection rate(Pd),fault isolation rate(PI),faults detection time(Td)and fault isolation time(TI).The faults detection rate is the ratio between the number of detected faults and the total number of faults in a period.The fault isolation rate is the ratio between the number of isolation faults and the number of detected faults.Faults detection time means the mean time of faults detection per time.And fault isolation time is the mean time of faults isolation per time.

    RMT parameters of equipment considered in this paper are shown in Table 1.

    Table 1 RMT parameters of equipment

    2.3 Assumptions

    The operation conditions of equipment are usually complex,and different equipment has different operation characteristics.To focus on the key point,some assumptions are made as follows:

    (i)MTBF,faults detection time,fault isolation time and maintenance time are assumed to obey the exponential distribution.

    (ii)The mean time between preventive maintenance is assumed to be as the same as PMTBF,which means the lead time of preventive maintenance is ignored.

    (iii)All the human factors are ignored.

    3.Solution

    As mentioned before,most design characteristics will affect the operation process of equipment in the operation phase.Therefore,taking the operation process into consideration is beneficial for determination of RMT parameters.Thus,we firstly establish the GSPN model of the operation process based on the operation flow of equipment.Afterwards,by solving the GSPN model,the quantitative relationship between operational availability and RMT parameters is obtained.Besides,the cost model can be built based on the life time and the history data of similar equipment.With costs,operational availability and available ranges of parameters considered,the determination and optimization model of RMT parameters can be established.Finally,the appropriate RMT parameters will be calculated though solving the model with some calculation tools,such as Matlab.The detail flow of the solution is shown as Fig.2.

    Fig.2 Workflow of optimizing RMT parameters of equipment

    3.1 The GSPN model of operation process

    The concept of GSPN is derived from stochastic Petri nets(SPN),and it extends timed and immediate transition conceptions[18].GSPN can be defined algebraically as follows[19]:

    A GSPN is an 8-tuple,GSPN=(P,T,A,V,S,B,W,M0),whereP={P1,P2,...,Pm}is a finite set of place;T={T1,T2,...,Tn}is a finite set of transitions,including timed transitions and immediate transitions;A?(P×T)∪(T×P)is a set of arcs;Vis a forbidden arc;Sis the fire rate of timed transitions;Bis activation probability of immediate tensions;Wis a weight function that takes values 1,2,3,...,andM0is the initial marking.

    Fig.3 GSPN model of the operation process of continuous working equipment

    GSPN is used in many application fields for the validation and evaluation of distributed systems which characterized by concurrency,synchronization and congestion[21,22].Compared with SPN,GSPN not only needs fewer calculation costs but also has the ability to model indeterminacy processes[23–25].In the operation process of equipment,some tasks,such as faults detection and diagnosis,are indeterminacy processes.Therefore,GSPN can be chosen to model the operation process of equipment.

    According to the operation flow of continuous working equipment which is shown as in Fig.1,the GSPN model of the operation process is established and shown in Fig.3.

    In the GSPN model shown as Fig.3,the places denote the states of the equipment and the physical meaning of each place is shown in Table 2.

    Table 2 Physical meanings of places

    In the GSPN model shown in Fig.3,immediate transitions mainly model the results of work which contains two different possibilities.Timed transitions mainly model physical work,such as faults detection,diagnosis and maintenance.

    Most of the working time in this operation process is stochastic variable.For example,the accurate maintenance time is determined by the occurring faults modes,while the faults modes occur randomly.In addition,for most continuous working equipment,it is reasonable to assume that the working time obeys exponential distribution and the parameter of the exponential distribution is the inverse of the mean time parameter,described as

    whereF(t)is the distribution function of the working time,andTMis the mean time of working per time.

    The physical meanings and parameters of all transitions of the GSPN model are shown in Table 3.

    In Table 2,the parameters of immediate transitions mean the probability of being fired,and the parameters of timed transitions are the inverse of the mean time.

    Table 3 Physical meanings and parameters of all transitions

    3.2 The solution of GSPN model

    Simulation and mathematical analysis are the two main methods to analyze the GSPN model.The simulation method employs computer technology to analyze the statistic results of the GSPN model,e.g.stochastic Petri net package(GPNP)is one of the widely used simulation anal-ysis software which was developed by America Duke University.The mathematical analysis method employs mathematical theory to get properties of the GSPN model.For example,the qualitative relationships between parameters and probability of stability can be obtained based on continuous time Markov process(CTMP)theory[23].Since we need to establish an optimization model,the qualitative relationship is crucial to achieve this.The mathematical analysis method is chosen to solve the GSPN model in this paper.

    3.2.1 Reachability graph and state-transition matrix

    All mathematical analysis methods are based on the theorem that the reachability graph of GSPN is isomorphic to a CTMP[26],and the reachability graph is the key point for mathematical analysis.Therefore,the reachability graph of the GSPN model should be obtained at first.

    The reachability graph of the GSPN model consists of states and the relationships among states.Based on the GSPN model of continuous working equipment(shown in Fig.3),the states of the model are defined and shown in Table 4,whereMi(i=0,1,...,8)are the reachable states of the GSPN model andPi(i=0,1,...,8)are places of the GSPN model.The relationships among states can be directly got from the GSPN model.Thus we can obtain the reachability graph of the continuous working equipment,as is shown in Fig.4.

    Table 4 States of GSPN model of continuous working equipment

    In the reachability graph of the GSPN model,states can be classified into virtual states and actual states[26].Virtual states refer to the states which connect with other states via immediate transitions.The remaining time of virtual states is zero.In order to simplify the calculation of CTMP,the virtual states should be removed from the reachability graph.According to the definition of virtual states and the GSPN model,M2andM4are virtual states.The simplified reachability graph is shown in Fig.5.

    Fig.4 Reachability graph of the continuous working equipment

    Fig.5 Simplified reachability graph of continuous working equipment

    According to the transition parameters(shown in Table 1)and the simplified reachability graph(shown in Fig.5),the state-transition matrix(P)can be calculated directly and is shown as(4).Pijdenotes the mean speed of transition from stateito statej.For example,the mean speed of transition fromM0toM1is 1/TBFNP.In addition,since the transition probability of each state is 1 and the transition speed obeys exponent distribution,the sum of elements of each row is zero:

    3.2.2 Probability of stability

    According to the basic characteristic of CTMP[27],probabilities of stability of all states must satisfy the following equation:whereΠ=[η0,η1,η3,η5,η6,η7,η8]contains the probabilities of stability ofMi(i=0,1,3,5,6,7,8)andPdenotes the state-transition matrix of CTMP.

    By solving(5),all the probabilities of stability can be obtained and shown as(6).

    3.3 Optimization model of RMT parameters

    To build an optimization model,it needs to determine constraint conditions and objective functions.Operational availability is important to evaluate the integrated property of equipment.Therefore,operational availability is a constraint condition.Besides,every RMT parameter should be limited in a reasonable range to make the results of optimization acceptable.As always,high quality and low price are two aspects of desired properties in the design phase.However,higher quality usually comes with higher cost.So,we choose the costs of equipment as the objective function.

    3.3.1 Operational availability

    The operational availabilityA0is the ratio between the available working time and the total life time.A0is an integrated parameter to evaluate the properties of the equipment,especially RMT[27].The operational availability of continuous working equipment is the probability of working states,which means,the operational availability is the probability stability of stateM0.

    3.3.2 Available range of RMT parameters

    RMT are the important design properties of equipment,and each RMT parameter has a special physical meaning.In order to ensure the optimization results of parameters to be acceptable,every parameter should be limited in a certain range,as shown in Table 5.Generally,the ranges of parameters are determined according to that of similar equipment and experts’suggestions.

    3.3.3 Cost model

    LCC contains all costs from the design phase to the decommission phase.Statistics data show that design costs,manufacture costs and maintenance costs are the main parts of LCC.Recently,some methods have been proposed to establish the cost model of equipment.In[28–30],the gray system theory was introduced to predict the costs ofequipment.The costs prediction model was developed based on support vector machine(LVM)in[31].In[32],the adaptive network-based fuzzy inference system(ANFIS)was used to estimate the costs of avionics equipment.In[33,34],the data processing for the cost model is studied.The parametric method and the analogy method are the two main methods to estimate the costs.The analogy method estimates the costs only based on similar equipment and experts suggestions.The parametric method calculates the costs based on the qualitative relationship between costs and design parameters[35].As the objective function and optimization variable parameter of the optimization model,the qualitative relation of costs and RMT parameters needs to be estimated.Therefore,we choose the parametric method to calculate the costs of the equipment.

    Table 5 The ranges of RMT parameters

    (i)Design and manufacture costs

    No general model for design and manufacture costs exists(CDM),since different equipment has different characteristics.Generally,the qualitative relationship between design parameters and the costs of design and manufacture is established through analyzing the history data of similar equipment.It requires enough history data on design parameters and the costs of design and manufacture.Based on the history data of similar equipment,we choose the least square algorithm(LSA)to estimate the design and manufacture costs.

    (ii)Maintenance costs

    Statistics show that maintenance costs consume over 60%of LCC[17].Therefore,reducing maintenance costs is significant to reduce LCC.The maintenance costs of continuous working equipment comprise accurate maintenance costs,fuzzy maintenance costs and predictive maintenance costs.

    whereCSis maintenance costs,CSAMis accurate maintenance costs,CSFMis fuzzy maintenance costs andCSPMis predictive maintenance costs.

    The mean costs per time of maintenance can be calculated based on the history data of similar equipment or experts’suggestions,and the maintenance costs can be calculated with the following equation.

    whereKAMis the times of accurate maintenance,KFMis the times of fuzzy maintenance andKPMis the times of preventive maintenance in all the life time.CAMis the costs per time of accurate maintenance,CFMis the costs per time of fuzzy maintenance andCPMis the costs per time of preventive maintenance.

    The times of maintenance can be calculated through the following equations:

    whereTALAMis the total time of accurate maintenance,TALFMis the total time of fuzzy maintenance andTALFMis the total time of preventive maintenance.Once the design life of equipment is determined,the time of maintenance can be calculated by the following equations:

    wherePAMis the probability of stability of the accurate maintenance state,PFMis the probability of stability of the fuzzy maintenance state andPPMis the probability of stability of the preventive maintenance state.According to the results of the analysis of GSPN in Section 3.2,the three important probabilities of stability can be calculated with the following equations:

    Therefore,

    3.4 Solving determination and optimization model

    Taking costs of equipment as the objective function and operation availability and available ranges of RMT parameters as constraint conditions,the determination and optimization model of RMT parameters can be established and is shown as(14).In this model,the constraint condition(A0min≤A0≤A0max)and the objective function(CDM+CS)are nonlinear functions.The general method to solve this nonlinear programming problem is the itera-tion algorithm.The first step of the iteration algorithm is to determine an appropriate start point.A better point can be found based on a certain searching algorithm.By iterative computation,the ideal parameters will be found in the end.Generally,the calculated quantity of the solution of the determination and optimization model is rather huge.However,many optimization tools,such as Matlab,Lingo and Mathematic[36–38],have the capability to alleviate it.

    4.Application

    To validate the effectiveness and practicability of the proposed approach,we apply it to an example of determination and optimization of RMT parameters of certain equipment.

    The operation flow of the equipment is shown in Fig.1 and the GSPN model is shown in Fig.3.

    4.1 Costs of the equipment

    (i)Design and manufacture costs

    The equipment is an improved version of certain equipment.The costs data of the original equipment can be used to modify the design and manufacture costs of the improved equipment.The history data of the original equipment are shown in Table 6.

    Table 6 Design and manufacture costs of similar equipment with different RMT parameters

    Based on the characteristics of similar equipment and experts’suggestions,the cost model is presented as follows:

    whereCDMis the design and manufacture costs,andλi(i=1,2,...,9)denotes undermined parameters.

    All the undermined parameters can be determined with the Matlab functionnlinfit.The result is shown as follows.

    (ii)Maintenance costs

    With the statistic results of the costs of accurate maintenance,fuzzy maintenance and preventive maintenance of the similar equipment,the costs per maintenance can be obtained and shown as in Table 7.The design life time of the equipment is 5 000 h.

    Table 7 The costs per maintenance

    The maintenance cost model can be calculated with(10)and the data are shown in Table 7.The result is shown as follows:

    4.2 Optimization model of RMT parameters

    (i)Available ranges of RMT parameters

    According to the task requirement and suggestions of experts,the available ranges of RMT parameters can be obtained as shown in Table 8.

    (ii)Operational availability

    Operational availability is an integrated index of RMT properties.According to the task requirement,the operational availability should be over 98%.

    Table 8 The available ranges of RMT parameters

    (iii)Determination and optimization model

    Taking costs of equipment as the objective function and considering operation availability and available ranges,the determination and optimization model can be established as follows:

    4.3 Result of optimization

    Many optimization tools can be used to optimize the model.We choose the Matlab optimization tool box to optimize the RMT parameters.With the nonlinear optimization tool and the interior point optimization algorithm,the optimization results are calculated as shown in Table 9.

    Table 9 The results of RMT parameters

    The operational availability,costs and other maintenance data are shown in Table 10.

    Table 1 0 Operational availability,costs and other maintenance data

    As shown in Table 9 and Table 10,all the RMT parameters are determined and the values satisfy operational availability constraint and ranges constraints.Therefore,the method proposed in this paper can be used to obtain the appropriate RMT parameters.In addition,the results of costs can be chosen as an important reference for expense budget of equipment.

    Compared with the classical method,such as[1,3,8],the method proposed in this paper has some advantages and are shown as follows:

    (i)The influence of RMT parameters on the operation process is adequately considered when determining parameters.

    Most of the methods(e.g.[3,6]),proposed to calculate these parameters,are only based on experts’suggestions or experience and do not consider the impact of the parameters on the operation process.While the impact is very important for determining appropriate RMT parameters.

    (ii)Reliability parameters,maintainability parameters and testability parameters can be determined and optimized synthetically.

    The three properties of equipment have strong dependency on each other and separately determining any of them is difficult to get the accurate results.

    (iii)In contrast to qualitative analysis methods,this paper builds a mathematic optimization model and the values of parameters can be obtained,which has significant meaning for engineering projects.

    5.Conclusions

    RMT are important design properties for equipment.The determination of RMT parameters has impact on not only the capability but also supportability and LCC of the equipment.In addition,the three properties interact with each other deeply and only synthetic determination can get appropriate parameters.

    In this paper,a determination and optimization model of RMT parameters is proposed based on the GSPN model of continuous working equipment.Moreover,the effectiveness and practical utility of the model is validated by an example of application.Although the main aim of the model is to determine and optimize the RMT parameters,the intermediate results,such as maintenance costs,can also be an important reference for expense budget of equipment.A limitation of this approach is the assumption that all the random variables obey exponential distribution.To make it more general,in the future this approach will be tested out with different distributions.

    [1]J.H.Liu.Reliability indexes determination and proof scheme for wideband tactical radio.Radio Communications Technology,2012,38(1):7–9.

    [2]X.J.Ma,H.Q.Zhang.Research on the contribution of reliability index adjustment to LCC of weapon equipment.Equipment Manufacturing Technology,2013,2013(7):52–86.

    [3]P.F.Guo,H.Z.Zhang,C.F.Yu,et al.Design requirements for construction machinery maintainability based on QFD.Science & Technology Review,2011,29(34):49–53.

    [4]F.Li,Y.Yang,J.J.Yang.Integrative evaluation on parameters of maintainability design of mechanical equipment.Journal of Engineering Design,2008,15(2):11–15.

    [5]X.M.Lu,K.L.Huang,G.Y.Lian,et al.Research on deciding testability figures of complicated equipment.Computer Measurement & Control,2008,16(3):357–359.

    [6]J.X.Li,G.Liu,P.Di.Research oil testiability index allocation method of warship equipment.Ship Electronic Engineering,2012,32(6):108–109.

    [7]Y.S.Zhang,K.L.Huang,G.Y.Lian.Research on selection of testability parameters of missile based on modified AHP.Computer Measurement & Control,2011,2012(2):412–414.

    [8]Y.L.Qian,J.Qiu,X.S.Wen.Generalized stochastic petri network for defining the systme-level testability parameters and figure of figure of merit.Systems Engineering and Electronics,2002,24(5):4–7.(in Chinese)

    [9]D.H.Ding,J.Lu.Discussion on the system of reliability,maintainability and supportability parameters.Journal of China Academy of Electronics and Information Technology,2011,06(2):171–174.(in Chinese)

    [10]B.de Jonge,W.Klingenberg,R.Teunter,et al.Optimum maintenance strategy under uncertainty in the lifetime distribution.Reliability Engineering and System Safety,2015,(133):59–67.

    [11]K.Saranga.Application of probabilistic maintenance models for selecting optimal inspection rates considering reliability and cost tradeoff.Transactions on Power Delivery,2014,29(1):178–186.

    [12]Z.G.Tian,D.M.Lin,B.R.Wu.Condition based maintenance optimization considering multiple objectives.Journal of International Manufacture,2012,(23):333–340.

    [13]M.J.Kim,V.Makis,R.Jiang.Parameter estimation in a condition-based maintenance model.Statistics and Probability Letters,2010,(80):1633–1639.

    [14]A.Jeang.Optimal parameters design and maintenance interval for a product with quality and cost considerations.Proc.of the Institution of Mechanical Engineers,Part B:Journal of Engineering Manufacture,2014,(223):737–747.

    [15]General Reserve Department of PLA.Terms of reliability,maintainability and testability,2005.

    [16]Y.X.Qing.Overview of reliability,maintainability and testability.Beijing:National Defense Industry Press,2002.

    [17]J.Qiu.Model and design technology of equipment testability.Beijing:Science Press,2012.

    [18]H.T.Wang,F.He,H.G.Xiong.Modeling of avionics blueprint architecture based on GSPN and LP.Aerospace Science and Technology,2013,(26):111–119.

    [19]G.Thangamani.Availability analysis of a Lube oil system using generalized stochastic Petri net.Proc.of the IEEE International Conference on Quality and Reliability,2011:186–189.

    [20]E.G.Amparore.A new great SPN GUI for GSPN editing and CSLTAmodel checking.Proc.of the 11th International Conference on Quantitative Evaluation of Systems,2014.

    [21]C.Eisentraut,H.Hermanns,J.P.Katoen,et al.A semantics for every GSPN.Proc.of the 34th International Conference of Petri Nets,2013:90–109.

    [22]G.Balbo,M.Beccuti,M.De Pierro,et al.First passage time computation in tagged GSPNs with queue places.The Computer Journal,2011,54(5):653–673.

    [23]R.Y.Al-Jaar.Performance evaluation of automated manufacturing systems using generalized stochastic Petri nets.New York,1989.

    [24]H.Zhan.Study of the normal generalized stochastic Petri nets and its application in testing system.Proc.of the Instrumentation and Measurment Technology Conference,2006.

    [25]S.A.Reveliotis,Y.C.Jin.A generalized stochastic Petri net model for performance analysis and control of capacitated reentrant lines.IEEE Trans.on Robotics and Automation,2003,19(3):474–480.

    [26]J.Y.Zhang.Modularized modeling and reliability analysis of manufacture systems based on GSPN & ESPN.Shenyang:Northeastern University,2007.

    [27]Z.C.Xu.Engineer and management of equipment supportability.Beijing:National Defense Industry Press,2006.

    [28]C.J.Wang.Research of costs of missile based on gray predictive model.Ship Electronic Engineering,2007,27(2):134–136.

    [29]M.Y.Yang,M.Z.Shen.Study on the estimation of the engine development cost based on the gray combined model.Mathematics in Practice and Theory,2006,36(10):161–166.

    [30]Y.Wang,F.Liu,P.Y.Zhang,et al.Prediction model of airforce equipment maintenance and support cost based on gray system theory.Mathematics in Practice and Theory,2012,42(9):137–141.

    [31]J.Y.Zhu,X.B.Zhang,H.X.Zhang,et al.Multi-parameter equipment cost estimation based on a support vector machine.Systems Engineering and Electronics,2003,25(6):701–703.(in Chinese)

    [32]J.Y.Zhu,J.F.Yu,H.X.Zhang.The application of ANFIS network in the study of the cost estimation of avionics equipment.Systems Engineering and Electronics,2002,24(3):57–60.(in Chinese)

    [33]Y.J.Zhao,M.Liu.Data diagnosis and recognition of costs parameter model of weapon system.Aerospace Shanghai,2000,2000(2):43–45.(in Chinese)

    [34]Y.J.Zhao.Data processing for the weapon types’cost parameter model.Aerospace Shanghai,2001,(5):18–20.(in Chinese)

    [35]S.Zhang,Y.G.Chen.Three development costs estamitation models based on paramter method.Journal of Sichuan Armaments Factories,2010,31(3):53–56.

    [36]M.L.Zhang.Computer implementation of solving of nonli-near programming problem.Journal of Armed Police College of Engineering,2006,22(2):5–9.

    [37]C.T.Yu,F.C.Li.Operations research.Beijing:Science Press,2006.

    [38]W.H.Chao,Z.Guo.Optimization method and MATLAB implementation.Beijing:Chemical Industry Press,2005.

    国产中年淑女户外野战色| 国产免费一级a男人的天堂| 国产精品av视频在线免费观看| 校园人妻丝袜中文字幕| 国产视频内射| 国产在线男女| 久久精品人妻少妇| 国产黄a三级三级三级人| 禁无遮挡网站| 午夜久久久久精精品| 少妇高潮的动态图| 精品一区二区免费观看| 免费av毛片视频| 国精品久久久久久国模美| 国产69精品久久久久777片| 波多野结衣巨乳人妻| 精品一区二区三区人妻视频| 七月丁香在线播放| 久久久精品欧美日韩精品| 久久久久久伊人网av| 真实男女啪啪啪动态图| 国产精品国产三级国产专区5o| 能在线免费观看的黄片| 蜜桃亚洲精品一区二区三区| 最近2019中文字幕mv第一页| 精品99又大又爽又粗少妇毛片| 国产在线男女| 精品熟女少妇av免费看| 国产一区二区在线观看日韩| 久久久久久久久久久免费av| 日韩欧美一区视频在线观看 | 成人亚洲欧美一区二区av| 永久免费av网站大全| 日韩欧美精品v在线| 床上黄色一级片| 亚洲国产最新在线播放| 久热久热在线精品观看| 国产一区有黄有色的免费视频 | 日韩欧美精品v在线| 精品久久国产蜜桃| 最近手机中文字幕大全| 丝瓜视频免费看黄片| 亚洲在久久综合| 久久热精品热| 日日啪夜夜撸| 欧美日韩在线观看h| 亚洲激情五月婷婷啪啪| 成年女人看的毛片在线观看| 国产精品不卡视频一区二区| 日韩一区二区视频免费看| 国产成人精品一,二区| 日日啪夜夜爽| 亚洲人成网站在线播| 国产一级毛片在线| 啦啦啦啦在线视频资源| a级毛色黄片| 亚洲精品第二区| 国产麻豆成人av免费视频| 精品一区二区三区视频在线| 99久国产av精品国产电影| 性色avwww在线观看| 性插视频无遮挡在线免费观看| 欧美zozozo另类| 一二三四中文在线观看免费高清| 亚洲最大成人手机在线| 18禁裸乳无遮挡免费网站照片| 能在线免费观看的黄片| 伦精品一区二区三区| 深夜a级毛片| 一夜夜www| 欧美成人精品欧美一级黄| 亚洲,欧美,日韩| 久久国产乱子免费精品| 六月丁香七月| 国产成人精品一,二区| 我的老师免费观看完整版| 又粗又硬又长又爽又黄的视频| 日韩亚洲欧美综合| 日韩欧美一区视频在线观看 | 亚洲国产成人一精品久久久| 熟女电影av网| 97在线视频观看| 国产成人aa在线观看| 秋霞伦理黄片| 91aial.com中文字幕在线观看| 嫩草影院新地址| 亚洲av不卡在线观看| 国产精品一区二区三区四区免费观看| 最近最新中文字幕免费大全7| 亚洲最大成人手机在线| 午夜视频国产福利| 久久久午夜欧美精品| or卡值多少钱| 国产91av在线免费观看| 亚洲欧美精品自产自拍| 狂野欧美激情性xxxx在线观看| 国产视频内射| 亚洲成色77777| 色综合亚洲欧美另类图片| 国语对白做爰xxxⅹ性视频网站| 直男gayav资源| 国产有黄有色有爽视频| 一级毛片电影观看| 99久国产av精品| 网址你懂的国产日韩在线| 亚洲国产av新网站| 久久国内精品自在自线图片| 欧美高清成人免费视频www| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 欧美日韩亚洲高清精品| 国产免费福利视频在线观看| 日韩av免费高清视频| 最近的中文字幕免费完整| 国产伦理片在线播放av一区| 亚洲av成人精品一二三区| 成人无遮挡网站| 亚洲在线自拍视频| 只有这里有精品99| 欧美xxxx黑人xx丫x性爽| 日本一本二区三区精品| 亚洲乱码一区二区免费版| 国产av国产精品国产| 直男gayav资源| 国内精品美女久久久久久| 色哟哟·www| 建设人人有责人人尽责人人享有的 | 日韩成人av中文字幕在线观看| 欧美极品一区二区三区四区| 日韩三级伦理在线观看| 老司机影院成人| 亚洲精品一二三| av免费观看日本| 欧美成人a在线观看| 日本一本二区三区精品| 亚洲自拍偷在线| 一二三四中文在线观看免费高清| 成人欧美大片| 亚洲欧美成人综合另类久久久| 日韩欧美精品免费久久| 久久亚洲国产成人精品v| 国产不卡一卡二| 欧美日韩国产mv在线观看视频 | 精品人妻一区二区三区麻豆| 91午夜精品亚洲一区二区三区| 超碰av人人做人人爽久久| 成年av动漫网址| 亚洲国产高清在线一区二区三| 国国产精品蜜臀av免费| 蜜桃久久精品国产亚洲av| 国产午夜福利久久久久久| 一级毛片电影观看| 欧美高清性xxxxhd video| 国产熟女欧美一区二区| 久久综合国产亚洲精品| 日日啪夜夜爽| 久久99蜜桃精品久久| 亚洲婷婷狠狠爱综合网| 国产精品人妻久久久久久| 日韩中字成人| 人妻少妇偷人精品九色| 天天一区二区日本电影三级| 国产男人的电影天堂91| a级毛色黄片| 成年版毛片免费区| or卡值多少钱| 国产高清三级在线| 亚洲成人中文字幕在线播放| 国产精品一区www在线观看| 亚洲性久久影院| 草草在线视频免费看| 国产片特级美女逼逼视频| 国产精品不卡视频一区二区| 亚洲av免费在线观看| 免费黄色在线免费观看| 全区人妻精品视频| 一区二区三区四区激情视频| 午夜免费观看性视频| 国产 一区精品| 男插女下体视频免费在线播放| 久久久久久久久中文| 在线观看一区二区三区| 最近2019中文字幕mv第一页| 熟女电影av网| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| 在线免费观看的www视频| 亚洲精品第二区| 一个人看的www免费观看视频| 中文精品一卡2卡3卡4更新| 啦啦啦啦在线视频资源| 日韩欧美 国产精品| 日本爱情动作片www.在线观看| 在线观看人妻少妇| 亚洲国产精品sss在线观看| 男插女下体视频免费在线播放| 亚洲精品一二三| 亚洲伊人久久精品综合| 超碰97精品在线观看| 99视频精品全部免费 在线| av一本久久久久| 狠狠精品人妻久久久久久综合| av在线蜜桃| 最近中文字幕高清免费大全6| 免费黄频网站在线观看国产| 欧美日韩精品成人综合77777| 国内精品宾馆在线| 国产白丝娇喘喷水9色精品| 又粗又硬又长又爽又黄的视频| 91精品国产九色| av在线蜜桃| 免费观看无遮挡的男女| av在线亚洲专区| 高清日韩中文字幕在线| 如何舔出高潮| 婷婷色麻豆天堂久久| 97热精品久久久久久| 亚洲精华国产精华液的使用体验| 亚洲精品色激情综合| 日本黄色片子视频| 国产爱豆传媒在线观看| 久久久久九九精品影院| 黄色一级大片看看| 久久综合国产亚洲精品| 波多野结衣巨乳人妻| www.色视频.com| 日韩一区二区三区影片| 久久久久久九九精品二区国产| 欧美激情在线99| 久久这里有精品视频免费| 最新中文字幕久久久久| 久久99热6这里只有精品| 国产熟女欧美一区二区| 男女视频在线观看网站免费| 特级一级黄色大片| 成人特级av手机在线观看| 丝瓜视频免费看黄片| 国产极品天堂在线| 日本一二三区视频观看| 国产精品1区2区在线观看.| 国产精品爽爽va在线观看网站| 午夜久久久久精精品| 国产精品无大码| 大香蕉97超碰在线| 久久久久久久久久人人人人人人| 欧美zozozo另类| 看黄色毛片网站| 色吧在线观看| 国产在线一区二区三区精| 亚洲精品,欧美精品| 亚洲四区av| 99久久九九国产精品国产免费| 97在线视频观看| 天堂影院成人在线观看| 又大又黄又爽视频免费| 女人被狂操c到高潮| 精品一区二区三区人妻视频| 最近手机中文字幕大全| 三级男女做爰猛烈吃奶摸视频| 青春草亚洲视频在线观看| 国产精品麻豆人妻色哟哟久久 | 亚洲在线自拍视频| 久久久久久久久大av| 如何舔出高潮| 深夜a级毛片| 国产精品美女特级片免费视频播放器| 成人亚洲欧美一区二区av| 深夜a级毛片| 国产精品不卡视频一区二区| 人体艺术视频欧美日本| 久久久久久久久大av| 国产淫片久久久久久久久| 国产精品久久久久久久久免| 一个人看的www免费观看视频| 精华霜和精华液先用哪个| 在线观看免费高清a一片| av线在线观看网站| 久久久久久久久久久免费av| 麻豆国产97在线/欧美| 免费看日本二区| 99re6热这里在线精品视频| 亚洲av.av天堂| av在线播放精品| 国产单亲对白刺激| 国产伦一二天堂av在线观看| 91av网一区二区| 亚洲高清免费不卡视频| 国产v大片淫在线免费观看| 麻豆久久精品国产亚洲av| 最近视频中文字幕2019在线8| 内地一区二区视频在线| 精品熟女少妇av免费看| 听说在线观看完整版免费高清| 国产乱人偷精品视频| 91午夜精品亚洲一区二区三区| 日韩欧美精品免费久久| 赤兔流量卡办理| 国产av在哪里看| 中国国产av一级| 日本三级黄在线观看| 国内精品美女久久久久久| 国产综合精华液| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 久久人人爽人人爽人人片va| 午夜精品在线福利| 免费看光身美女| 国产精品人妻久久久影院| 久久久久性生活片| 黄色一级大片看看| 黑人高潮一二区| 国产精品一二三区在线看| 中国美白少妇内射xxxbb| 日本色播在线视频| 国产免费视频播放在线视频 | 国产精品综合久久久久久久免费| 舔av片在线| 国产久久久一区二区三区| 插逼视频在线观看| 国产国拍精品亚洲av在线观看| 久久久久久久大尺度免费视频| 精品国产一区二区三区久久久樱花 | a级毛片免费高清观看在线播放| 午夜爱爱视频在线播放| 亚洲av二区三区四区| 色综合站精品国产| 亚洲av电影不卡..在线观看| 成人毛片60女人毛片免费| 夫妻午夜视频| 久久精品国产鲁丝片午夜精品| 亚洲精品视频女| 午夜福利视频精品| 国产午夜精品一二区理论片| 丰满人妻一区二区三区视频av| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 91aial.com中文字幕在线观看| 一区二区三区乱码不卡18| 波多野结衣巨乳人妻| 亚洲熟女精品中文字幕| 亚洲国产最新在线播放| 国产成人福利小说| 久久韩国三级中文字幕| 啦啦啦韩国在线观看视频| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 久久精品久久久久久久性| 最近手机中文字幕大全| 欧美变态另类bdsm刘玥| 美女大奶头视频| 久久6这里有精品| 中文字幕av在线有码专区| 久久国产乱子免费精品| 精品一区二区三区视频在线| 人人妻人人看人人澡| 亚洲av成人av| 直男gayav资源| 美女主播在线视频| a级毛色黄片| 午夜福利网站1000一区二区三区| 丰满少妇做爰视频| 久久久久久久久久成人| 中文精品一卡2卡3卡4更新| 色尼玛亚洲综合影院| 免费观看a级毛片全部| 日韩 亚洲 欧美在线| 久久久久国产网址| 国产成人福利小说| 亚洲人成网站在线播| 亚洲真实伦在线观看| 久久久色成人| 久久鲁丝午夜福利片| 久久久久久国产a免费观看| 国产在视频线在精品| 日本与韩国留学比较| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 夫妻午夜视频| 最近最新中文字幕大全电影3| 亚洲国产精品sss在线观看| 国产成人aa在线观看| www.av在线官网国产| 乱人视频在线观看| 国产在线一区二区三区精| 哪个播放器可以免费观看大片| 久99久视频精品免费| 成人高潮视频无遮挡免费网站| 少妇的逼好多水| 直男gayav资源| 91av网一区二区| 亚洲国产精品成人久久小说| 久久亚洲国产成人精品v| 久久久久久国产a免费观看| 中文字幕亚洲精品专区| 久久久亚洲精品成人影院| 男人和女人高潮做爰伦理| 欧美激情久久久久久爽电影| 久久久久久国产a免费观看| 成人一区二区视频在线观看| 97超视频在线观看视频| 晚上一个人看的免费电影| 蜜臀久久99精品久久宅男| 亚洲四区av| 欧美人与善性xxx| 最近视频中文字幕2019在线8| 别揉我奶头 嗯啊视频| 97人妻精品一区二区三区麻豆| 99热6这里只有精品| 午夜福利在线观看免费完整高清在| 久久99热这里只频精品6学生| 高清视频免费观看一区二区 | 亚洲精品中文字幕在线视频 | eeuss影院久久| 建设人人有责人人尽责人人享有的 | 亚洲av福利一区| 蜜桃久久精品国产亚洲av| 美女大奶头视频| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载 | 搡老乐熟女国产| 青春草视频在线免费观看| 日本wwww免费看| 国产精品蜜桃在线观看| 欧美日韩在线观看h| 久久99热6这里只有精品| 黄色日韩在线| or卡值多少钱| 成人综合一区亚洲| 成人国产麻豆网| 国产黄色视频一区二区在线观看| 成人二区视频| 夜夜看夜夜爽夜夜摸| 一二三四中文在线观看免费高清| 国产伦在线观看视频一区| 久久久国产一区二区| 伊人久久国产一区二区| 精品人妻视频免费看| 狂野欧美白嫩少妇大欣赏| 国内精品美女久久久久久| 亚洲av一区综合| kizo精华| 九色成人免费人妻av| 久久精品熟女亚洲av麻豆精品 | 超碰av人人做人人爽久久| 人人妻人人看人人澡| 少妇的逼好多水| 黄色欧美视频在线观看| 久久久久久久久大av| 免费黄网站久久成人精品| 男女边摸边吃奶| 亚洲综合精品二区| 国产一区二区三区综合在线观看 | 最近2019中文字幕mv第一页| 女人久久www免费人成看片| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 日韩中字成人| 麻豆av噜噜一区二区三区| 亚洲乱码一区二区免费版| 国产视频内射| 国产老妇伦熟女老妇高清| 精品一区二区三卡| 日日撸夜夜添| 久久精品国产鲁丝片午夜精品| 日韩中字成人| 亚洲国产成人一精品久久久| 亚洲欧美日韩无卡精品| 久久久久性生活片| 男人舔奶头视频| 精品久久久久久久末码| 免费av不卡在线播放| 视频中文字幕在线观看| 欧美高清成人免费视频www| av一本久久久久| 精品一区二区三卡| 亚洲国产欧美在线一区| 女人被狂操c到高潮| 人妻夜夜爽99麻豆av| 在现免费观看毛片| 成年女人看的毛片在线观看| 久久国内精品自在自线图片| 九九在线视频观看精品| 欧美xxxx黑人xx丫x性爽| 成人高潮视频无遮挡免费网站| 国产精品女同一区二区软件| 尾随美女入室| 在线观看美女被高潮喷水网站| 欧美激情国产日韩精品一区| 午夜亚洲福利在线播放| 久久精品国产亚洲av涩爱| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 国产精品一及| 国产在线一区二区三区精| 男女下面进入的视频免费午夜| 国产69精品久久久久777片| av在线观看视频网站免费| 美女高潮的动态| 深夜a级毛片| 麻豆精品久久久久久蜜桃| 美女大奶头视频| 久久久久久久大尺度免费视频| 精品久久久久久久久亚洲| 久久久久久久国产电影| 国内揄拍国产精品人妻在线| 精品亚洲乱码少妇综合久久| 日韩中字成人| 99热这里只有是精品在线观看| 日韩av不卡免费在线播放| 最近的中文字幕免费完整| 亚洲精品乱码久久久久久按摩| 一区二区三区免费毛片| 大话2 男鬼变身卡| 久久久久精品性色| 日日啪夜夜撸| 久久草成人影院| 日韩欧美精品v在线| 亚洲av不卡在线观看| 亚洲av成人精品一二三区| 伊人久久国产一区二区| 久久久久久久久久黄片| 日本一本二区三区精品| 有码 亚洲区| 亚洲欧美成人精品一区二区| 国产亚洲最大av| 国产又色又爽无遮挡免| 日本猛色少妇xxxxx猛交久久| 亚洲丝袜综合中文字幕| 高清av免费在线| 国产成人午夜福利电影在线观看| 青青草视频在线视频观看| 伦理电影大哥的女人| 91精品伊人久久大香线蕉| 国产老妇女一区| 高清欧美精品videossex| 午夜福利高清视频| 国产视频内射| 欧美日韩精品成人综合77777| 别揉我奶头 嗯啊视频| 日日撸夜夜添| 观看免费一级毛片| 亚洲精品色激情综合| 亚洲在线观看片| 色综合站精品国产| 国产成人午夜福利电影在线观看| 日本免费a在线| 久久精品国产亚洲av涩爱| 超碰av人人做人人爽久久| 波野结衣二区三区在线| 99久久精品国产国产毛片| 男女边摸边吃奶| 韩国高清视频一区二区三区| 亚洲精品成人久久久久久| 男人狂女人下面高潮的视频| 哪个播放器可以免费观看大片| 久久久久久久久久久丰满| av又黄又爽大尺度在线免费看| 小蜜桃在线观看免费完整版高清| 日韩成人av中文字幕在线观看| 精品一区二区三卡| 97在线视频观看| 中文字幕av在线有码专区| 青青草视频在线视频观看| 日韩,欧美,国产一区二区三区| 嫩草影院入口| 欧美日韩一区二区视频在线观看视频在线 | 一级毛片黄色毛片免费观看视频| 中文乱码字字幕精品一区二区三区 | 国产亚洲最大av| 成人高潮视频无遮挡免费网站| 亚洲欧美清纯卡通| 日韩一本色道免费dvd| 国产精品伦人一区二区| 最新中文字幕久久久久| 又黄又爽又刺激的免费视频.| 激情五月婷婷亚洲| 日韩成人伦理影院| 69人妻影院| 插逼视频在线观看| 日本一二三区视频观看| 日韩av在线免费看完整版不卡| 麻豆av噜噜一区二区三区| 久久精品国产自在天天线| 中文在线观看免费www的网站| 国产国拍精品亚洲av在线观看| 人妻一区二区av| 国产精品久久久久久av不卡| 三级国产精品欧美在线观看| 在线免费十八禁| av免费观看日本| 国产高清三级在线| 亚洲av电影在线观看一区二区三区 | 伊人久久国产一区二区| 一级a做视频免费观看| 高清视频免费观看一区二区 | 免费大片黄手机在线观看| 一区二区三区四区激情视频| 在线天堂最新版资源| 少妇人妻精品综合一区二区| 精华霜和精华液先用哪个| 国产乱人偷精品视频| 天堂√8在线中文| 免费人成在线观看视频色| 久久久国产一区二区| 国产亚洲精品av在线| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久久久人人人人人人| 精品一区二区三区视频在线| 黄片无遮挡物在线观看| 亚洲av在线观看美女高潮| 五月玫瑰六月丁香| 久久久欧美国产精品| 啦啦啦中文免费视频观看日本| h日本视频在线播放| 亚洲四区av|