Abstract

In this paper, the mathematical model of a kind of two-robot security system with an early warning function is studied. By using strongly continuous operator semigroup theory and Volterra integral equation theory, the properties of the semigroup of the system operator, the existence and uniqueness of nonnegative solution, and the well-posedness of solution are discussed, respectively. Under the assumption that the failure rate and repair rate of the system are constants, the equations of the robot system are transformed into an ordinary differential equation group, and then the instantaneous reliability and stable-state reliability of the system are obtained. The reliability and zero-state controllability of the system are proved. Finally, the numerical solution of the system model is obtained by using MATLAB mathematical software, and the corresponding numerical simulation diagram is given. The results show that the conclusions of numerical calculation and numerical simulation are in accordance with the results of reliability theory, thereby the reliability of the robot safety system is verified.

1. Introduction

With the progress and development of the times, robots are becoming more and more widely used. Robots are mechanical devices that automatically perform work, replacing or assisting human work, such as manufacturing, construction, or hazardous work. Therefore, the stability and reliability of the robot security system has become a hot issue. Early warning function refers to the ability to provide early warning and timely issuance of alarm signals before a system or component malfunctions, allowing sufficient time to predict various impending disasters and reduce economic losses caused to humans by disasters. The reliability [1] of a system means the ability or property of the system to accomplish the specified functions under the specified time and the specified using conditions. It is also one of the important characteristics of a repairable system. The robot security system with early warning function can accurately detect the abnormal condition of the parts of the robot working system according to the early warning signal and carry out advanced control, continue to work, or carry out maintenance, so as to reduce the economic losses. In order to ensure the security and reliability of the system and avoid the occurrence of unexpected accidents, the early warning function is introduced into the repairable robot security system in this paper. Based on the abstract Cauchy problem theory, the system reliability is analyzed.

The repairable system is an important system in the research of reliability mathematics, and reliability is one of the important contents in the research of the repairable system. Many scholars and researchers at home and abroad have done a lot of research on this kind of system, and achieved fruitful results. Some of the literatures adopt the method of qualitative analysis; for example, Wang et al. used linear operator semigroup theory in literatures [29] to study the semigroup properties of the main operator of repairable systems composed of faulty components and repairmen and discussed the well-posedness of the system model solution by using functional analysis methods. In literatures [1015], Kamranfar et al. studied the stability of repairable systems. On the basis of discussing the asymptotic stability of the solutions of a repairable system with two different components connected in parallel, Wang et al. studied the controllability of the system in the zero state and the criterion of using the controllability of the system in literatures [1624], and obtained the corresponding optimal control element of the system by using the minimization sequence, thus solving the optimal control problem of the system solution. In other literatures, the reliability index of the system is obtained by quantitative analysis. For example, Marco et al. obtained the reliability index of the system in the literatures [2532] by using the Laplace transformation inversion method, MATLAB mathematical software calculation, linear differential equation with constant coefficients satisfying initial conditions, and so on. In addition, some literatures combine the quantitative analysis method with computer technology to perform the numerical calculation and simulation of the repairable system. For example, Gao et al. studied a repairable system with early warning function in literature [33], utilizing that when the risk coefficient α ⟶ ∞, the warning system approximates a new model with weak solution—the model of a nonearly warning system; the corresponding conclusions were drawn and numerical simulation examples were given by computer, which greatly enriches the theory and practice of the repairable system with an early warning function; Zhou et al. studied the reliability of the system in literatures [3437] and simulated the graphs of transient and steady-state reliability of the system by using Maple software. In summary, in the existing literatures and materials both domestically and internationally, the method of combining qualitative analysis with quantitative analysis is less effective than that of using qualitative or quantitative analysis alone. In this article, a kind of two-robot security system with an early warning function is taken as the research object, and attempts are made to combine qualitative analysis with quantitative analysis to study the reliability of the system.

The main contents of this article are the modeling of the robot safety system, the well-posedness and controllability of the system solution, the reliability index of the system, and the numerical experiment of the reliability theory results. However, with the development and application of system reliability, when the system models become more complex, new methods and theories are needed to guide them, which is an important development problem of system reliability research in the future.

The reliability of repairable system has achieved certain achievements in both theory and application. Nowadays, in the 21st century, the research on reliability has been elevated both domestically and internationally to the high level of understanding of saving resources and energy. Products (or equipment) are developing in the direction of gradually improving reliability and gradually reducing maintenance time, maintenance personnel, and maintenance costs. At the same time, reliability technology is the result of the joint development of multiple fields and technologies. It belongs to a comprehensive basic industry, and the current development trend is moving towards comprehensiveness, usability, automation, informatization, virtualization, and intelligence. Thus, higher economic benefits and stronger competitiveness can be achieved.

2. Mathematical Model

A robot is a complex system including mechanical, electronic, electrical, hydraulic, pneumatic, computer, and other types of components and control software. It is relatively complicated to study its reliability and safety [15, 16]. In order to consider this issue more clearly, this article assumes that the robot security system consists of two robots、security devices and a repairman [17]. Assuming that the entire system does not consider the repair and replacement process and starts to run at the time t = 0. At the time t = 0, the two robots and the safety device are brand new, the system starts to operate normally, and the maintenance personnel goes on holiday. If the system fails, the repairman immediately terminates the vacation and repairs the faulty system immediately. If N(t) represents the state at the moment t, the system has the following situations:State 0. One robot and safety device work, and one robot is in hot standby stateState 1. One robot and safety device are working, and one robot is in a malfunction stateState 2. The state in which both robots are malfunctioningState 3. The state of the system when the safety device failsState c. The state of the normal faulty system

The state-transfer chart of the system is shown in Figure 1.

In order to facilitate the modeling and model analysis, the following general assumptions can be made according to the state transition diagram of the repairable system:(1)Various faults are independent of each other in the statistical sense(2)Only when both robots fail (or the safety device fails due to conventional reasons), the whole system is in a fault state(3)The failure rate of the system is constant, and the repair rate after system failure is nonconstant (4)The normal working time of the robot security system follows the exponential distribution function , , and (5)The repair time of faulty parts of robot system follows the general distribution function , , , and (6)The two robots are exactly the same and repaired as new

Since the time distribution of transitions between states of the system does not completely obey the negative exponential distribution, it can be known from the above hypothesis that represents the state at the moment t, so is not a Markov process. But we can make it a high-dimensional Markov process by using the method of supplementary variable, assuming that the fault system is in a state of maintenance, the supplementary variable represents the maintenance time from the beginning of maintenance to the present, and let

Then, constitutes a Markov process; that is, is a continuous-time two-dimensional Markov process; that is, t at any time, given the concrete values of and , then the probability law of the process after time t has nothing to do with the history of the process before time t. Let denote the probability that the system is in state i at time t, represents the probability density of the time t that the system is in state i and the faulty part has been repaired x, i.e.,

Here, it should be noted that although is only defined as , but for the sake of discussion, can be defined according to the actual physical background of the system and is extended on , that is, supplementary definition . At the same time, the system is out of state of the risk function; that is, the fix quotiety of the faulty component in the system in the state can be defined as follows by the conditional probability [21, 22]:

And from the actual physical meaning of , the following reasonable assumptions can be made:

The following is a discussion of the system state transition after time. Therefore, let represents the damage rate of the running system robot caused by its own reasons, indicates the normal failure rate of the system in the state , indicates the human fault quotiety of the operating system, and represents the damage quotiety of hot standby robot, indicates the constant fix quotiety of the running robot, indicates the constant fix quotiety of the running system, and indicates the probability that the system is in state i at the moment , represents the probability that the system is in the state i and the repaired time x at the moment t, in . represents the fix quotiety when the system is in state i and the repaired time x. Then, it is deduced from the formula of total probability and the properties of Markov process (for convenience, it is assumed that is the same as ):

(at t when the system is in state 0, the system remains in state 0) + P (at t when the system is in state 1, the system is fixed to state 0) + P (at t when the system is in state c, the system is fixed to state 0) + P (at t when the system is in state 2, the system is fixed to state 0) + P (at t when the system is in state 3, the system is fixed to state 0):

From formula (5) and the definition of derivative,

(at t the system is in state 1, the faulty robot in has not been repaired, and the other robot has not failed, and the system in has not left state 1) + P (at t the system is in state 0, and one robot in has failed):

From formula (7) and the definition of derivative,

(at t when the system is in state c, the system routine fault has not been repaired, and the system has not left the state c) + P (at t when the system is in state 0, the system routine fault) + P (at t when the system is in state 1, the system routine fault):

From formula (9) and the definition of derivative,

Similarly, there are

Derived from formula (11) and the definition of partial derivative,

The boundary conditions and initial conditions of the system are discussed below [20].

Because represents the probability that the system just enters the state 2 at t, that is, the probability that the system just leaves the state 1 at t, i.e.,

From formula (13) and the definition of derivative,

In addition, represents the probability that the system just enters state 3 at t, that is, the probability that the system just leaves state 0 or state 1 at t, i.e.,

From formula (15) and the definition of derivative

Assuming that time , both parts are good, that is, the initial condition is

At this time, from probability analysis formulas (5) to (17) [18, 19], the integro-differential equations of the mathematical model of the two-robot security system with early warning function can be described in the Figure 1 are

3. Well-Posedness of the System Solution

Since the repairable robot system with early warning function contains both integral and differential, it is difficult to directly deal with it. Therefore, it is necessary to perform the necessary conversion before the reliability analysis [23]. In order to facilitate the subsequent discussion of the adaptability of the system solution, the system equations (18)–(23) are translated into an abstract Cauchy problem in Banach space [24, 25]. For this definition, make .

Take state spaceamong which . At this point, it can be proved is a Banach space. The following operators A and B and their domains are defined as follows:

For , the range of the operator A is defined as follows:

And for , the range of the operator B is defined as follows:

Therefore, the system equations (18)–(23) can be rewritten as an abstract Cauchy problem in Banach space:

Several lemmas are given below, which play a significant part in describing the semigroup characteristics of system operators.

Lemma 1 (see [26]). When , have .

Proof. For any , discuss the equation , which is equivalent toTherefore,So, from equations (29) and (30) and Fubini theorem, we haveAmong them, and since , is a bounded linear operator. and .

Lemma 2 (see [27]). Dense in X.

Proof. Suppose and there is a constant , so that for any , there are . Construct a setAt this time, using the knowledge of functional analysis, it is not difficult to verify that L is dense in Banach space X. Therefore, to verify that is dense in X, just prove is dense in L. For this, take ; then, for any , there is a constant ; when , there is always, so when , , there is always . LetAmong themAt this point, it is easy to verify , andThe above mentioned formula indicates that is dense in L, so is dense in X.

Lemma 3 (see [28]). Operator creates positive compressed hemigroup .

Proof. We split the proof this theorem into three steps.
The first step is to verify that creates hemigroup .
It can be verified by the demi-closed operator theorem that A is a closed linear operator; from the definition of B operator, B is a bounded linear operator, andFrom the above mentioned proof and the Hille–Yosida generation theorem, we know that the operator A creates a hemigroup, and then from the perturbation theorem of the semigroup, the operator creates a hemigroup, denoted as .
The second step is to verify that creates a hemigroup that is a positive semigroup.
In fact, when is a nonnegative vector, P is a nonnegative vector, so is a positive operator. At the same time, it is known from the meaning of B that B is also a positive operator. Hence,So, when , there is, which exists and is bounded becauseTherefore, is also a positive operator, so when , is a positive operator, so the semigroup created by is a positive semigroup.
The third step is to prove that creates positive compressed hemigroup .
In fact, whatever , letAmong them , .
Therefore,Then, is a diffusion operator. From the above mentioned proof and Pillips theorem, creates a positive contraction hemigroup . Therefore, by the uniqueness theorem of a hemigroup, is a positive contraction hemigroup.
The well-posedness of the system solution is discussed below; because the system equations (18)–(23) are complex equations composed of differential, partial differential and integral, it is difficult to solve them directly. Therefore, in order to discuss the existence and uniqueness of the nonnegative solution of the system, system (18)–(23) is converted into the form of convolution Volterra integral equation. To this end, notation is introduced:For convenience, is still represented by in the following system equations. can be obtained from the partial differential equations (21):First, substituting (42) into (18), then make as a variable to obtainSince the initial condition  = 1, , can be obtained as follows:Among them,Similarly, and can be solved by (19) and (20), respectively:Among them,Substituting formulas (44) and (46) into formula (22), respectively:Put expressions (44), (46), (47), (49), and (50) together to form the convolution Volterra integral equations:And write the Volterra integral equation in vector formAmong them,According to the above mentioned equation, the following theorem holds.

Theorem 4. The existence and uniqueness of nonnegative solutions for the system (18)–(23) are a necessary and sufficient condition for the existence and uniqueness of nonnegative solutions for the Volterra integral equation (52).
For any , from the expressions of and , each component of and each of are nonnegative bounded functions, according to documents [36, 37] the following theorem holds.

Theorem 5. The Volterra integral equation (52) has a unique nonnegative solution on .
According to Theorem 5 and the specific expression of the solution of , we know that has a unique nonnegative solution on . So, the Volterra integral (52) has a unique nonnegative solution on , that is, exists and is unique on .
To sum up, the following theorem holds.

Theorem 6. The system (18)–(23) has a unique nonnegative solution on .
The main results of this paper and the well-posedness conclusions of the system solution are given below.

Theorem 7 (see [29]). Abstract Cauchy problem has a sole nonnegative solution which satisfies

Proof. From the above Lemmas 13, it can be seen that the abstract Cauchy problem has a sole nonnegative solution which satisfies , also ; therefore,At the same time, since satisfies the system model, if is regarded as time t function, there areTherefore, .

Theorem 8 (see [30]). 0 is the simple eigenvalue of the operator .

Proof. Discuss the equation , that is,Solve the above equations to getAssume , from the system of equations (57) and (58),It can be seen that, , so is the eigenvector of 0 eigenvalue corresponding to the operator .
Take , then we have . For arbitrary, have ; thus, 0 is the simple eigenvalue of .

4. Reliability of the System Solution

The reliability of the system is one of the significant contents of the repairable system model. In the above system (18)–(23), because the size of indicates the probability that both robots are in good condition and can work normally at , the larger the is, the closer the system is to working properly; therefore, the size of determines the probability that the system will work properly and thus gets one of the important factors affecting the reliability of the system.

So as to discuss the reliability of the system, several definitions are given at first.

Definition 9. is called the transient reliability of system (18)–(23).

Definition 10. If exists, then is called the stable-state reliability of the system (18)–(23).

Definition 11. If , the system is said to be reliable [31].

Theorem 12. The system (18)–(23) has a time-dependent asymptotically steady stable-state solution , namely,where is called the stable-state solution of the system (18)–(23).

Theorem 13. When the fault quotiety and fix quotiety are constants, the system (18)–(23) has reliability.

Proof. Let the fault quotiety and fix quotiety be constants, i.e.,Simultaneously letThen, the actual physical background of the system (18)–(23) is , which can be transformed into formula (18):A system of ordinary differential equations can be converted from the upper system (18)–(23):At this point, ifThen, the first three equations of system (18)–(23) can be converted into an abstract Cauchy problem:The solution of the system (66) is obtained by using ordinary differential equation theory and advanced algebra knowledge, including the following four steps:Step 1: Find all eigenvalues of matrix AThat is, , then the eigenvalues of the matrix A are .Step 2: Seek , letAmong them,Substitute the formulas (69) and (70) into the formula (68) and sort them outStep 3: Find . From the knowledge of linear algebra, it is easy to find the inverse matrix of A, that is,Step 4: Finding the solution of the system (66), we haveSubstitute the above formulas and into the formula (73) and sort them outSo, the transient reliability of the system (18)−(23) isThus, by definition of 10, the stable-state reliability of system (18)−(23) isFrom the above mentioned discussion, it is obvious that you can get , Therefore, according to the definition of 11, the system (18)–(23) is reliable.

5. Controllability of System Zero State

Using the method of functional analysis to find a control element , to study the controllability of , in the system model can be transferred to the specified state in a finite time , and the allowable control set is selected as follows:

Theorem 14. Let the failure rate , the repair rate , and let be the probability that the system reaches the desired state at finite time and fits , then there is such that .

Proof. By the formula of (75),Considering as a function of the variable , we derivesinceand thereforeThus, we haveThis formula shows that is a monotonically increasing function about the variable , Note that , . Therefore, for any , according to the intermediate value theorem, there is , such thatThe above mentioned proves that the zero state of the system is controllable.
Theorem 14 indicates that the zero state of the system is controllable, in other words, the initial state of the system is controllable, but it does not mean that other states of the system are controllable, that is, the conclusion of controllability of the whole system is not necessarily derived from zero-state controllability, Similarly, to prove that the system is completely controllable, it is also necessary to prove other states of the system, such as , etc. are controllable.

6. Numerical Simulation

According to Theorem 7, there is a unique nonnegative time-dependent solution for the model of a repairable robot system with early warning function. In addition, the reliability of the two-robot security system with the function of early warning is discussed by using the theory and method of ordinary differential equation, On this basis, the controllability of the zero state of the system is proved by using the method of functional analysis. The following is a numerical simulation of the above results by using the numerical calculation method, which verifies the correctness of the results of the system reliability theory [32]. Assuming that the system failure rate and repair rate are constants , respectively, i.e.,

To do this, as shown by the above proof of reliability Theorem 13, the system is converted into an ordinary differential equation system (64), solve the linear equations, we have

At this time, if , the numerical solution of system (18)–(23) can be obtained by using Matlab mathematical software, and the result is shown in Figure 2.

It can be seen from Figure 2 that the system equations (18)–(23) has a time-dependent asymptotically steady stable-state solution , which is consistent with the conclusion of the Theorem 12, at the same time, Theorem 12 indicates that the system equations (18)–(23) is steady. It should be pointed out here that the importance of the Theorem 12 lies in that it not only proves the asymptotic stability [33] of the solution of the system equations (18)−(23), but also proves the existence of the stable-state solution of the system related to the stable-state reliability of the two-robot safety system with early warning function.

In addition, by selecting , respectively, the 3 groups of instantaneous reliability and stable-state reliability of the system equations (18)–(23) can be obtained by using MATLAB software. For more intuitive comparison, put the charts of and together. The result is shown in Figure 3:

It can be seen from the Figure 3 that the transient reliability curves of the system equations (18)–(23) are all above the stable-state reliability curves [34], when , can be obtained by mathematical limit analysis, which is consistent with the reliability of the conclusion system equations (18)–(23) of the Theorem 13. To some extent, the above results show that the numerical method can reflect and depict that each state of the system tends to be stable with the change of time, and the results are in line with the actual situation of the robot security system.

7. Concluding Remarks

In this paper, the mathematical model of a two-robot safety system with an early warning function is studied, the semigroup characteristics of system operators are discussed by using linear operator semigroup theory, and the well-posedness of the solution of the system is proved. Under the assumption that the fault quotiety and fix quotiety of the system are constants, the early warning model equations are converted into ordinary differential equations, the transient reliability and stable-state reliability of the system are obtained, and the reliability and zero-state controllability of the system are proved. Finally, the numerical solution of the ordinary differential equation set of (64) is obtained by using Matlab mathematical software, and the graphs of the transient reliability and stable-state reliability of the system of equations (18)–(23) are simulated, which shows that the results obtained by numerical calculation and numerical simulation are consistent with the proof of the above theory [35].

The research on the repairable system model with an early warning function is mainly carried out by qualitative analysis, quantitative analysis, and the combination of qualitative analysis and quantitative analysis. Most of the existing literatures adopt qualitative analysis methods, but few use quantitative analysis methods. And the application of the qualitative analysis combined with the quantitative analysis method in the repairable system model is rarely reported in the repairable system model. First, the innovation of this paper is to combine the qualitative analysis with the quantitative analysis method in order to perfect and enrich the theory and method of repairable systems. Second, in the early warning system, when the warning prompt is invalid, the robot system model with an early warning function approaches to a repairable system model without an early warning function; the relationship between the early warning system and the nonearly warning system should be further studied, and the relationship between their steady-state solutions should be discussed. Finally, when studying the important indexes of system reliability, transient reliability, and steady-state reliability because of the relative complexity of the solution of the system (66) model, this paper uses the method of solving linear differential equations with constant coefficients that satisfy initial conditions, which is computationally complex. Whether the solution process can be simplified by the methods of Laplace transformation inversion, MATLAB mathematical software coding and compiling, and so on remains to be further studied.

Data Availability

All data generated or analyzed during this study are included in this article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Fund project: Innovation Team of Universities in Henan Province (21IRTSTHN014).