Application of the Methodology of Multi-Valued Logic Trees with Weighting Factors in the Optimization of a Proportional Valve

: Hydraulic valves are used to determine the set values of hydraulic quantities (ﬂow rate, pressure, or pressure difference) in a hydraulic system or its part. This is achieved through the appropriate throttling of the stream ﬂowing through the valve, which is automatically set by the operator (e.g., opening the throttle valve). The procedures for determining its static and dynamic properties were described using the example of modeling a two-stage proportional relief valve. Subsequently, the importance of the design and operational parameters was determined using multi-valued logic trees. Modeling began with the determination of equations describing the ﬂow and movement of moving parts in a valve. Based on the equations, a numerical model was then created, e.g., in the Matlab/Simulink environment (R2020b). The static characteristics were obtained as the result of a model analysis of slow changes in the ﬂow rate through the valve. Various coefﬁcients of logical products have not been taken into account in the separable and common minimization processes of multi-valued logic equation systems in any available literature. The results of the model tests can be used to optimize several types of hydraulic valve constructions.


Introduction
In recent years, intensive development in the field of hydraulic valves has been observed. This development is mainly related to the integration of electronics designed to control the valves. Modern hydraulic valves-especially those controlled via the proportional technique-are often equipped with various types of sensors, e.g., an inductive spool-position sensor inside the body of a proportional valve. The integration of classic hydraulics with electronics and sensors creates new, previously unattainable possibilities for using hydraulic proportional valves [1]. The course of the control signal of proportional valves is shown in the form of a block diagram in Figure 1.
An analog electrical signal with a voltage value typically not greater than 10 V is fed to an electronic amplifier. From the electronic amplifier, the electric control signal is fed through wires with a current that usually does not exceed 1.5 A per coil of a proportional electromagnet. Depending on the type of proportional electromagnet, a force or displacement of the proportional electromagnet armature is generated. If the valve uses a proportional solenoid with an adjustable stroke, a displacement of the electromagnet armature is generated proportionally to the value of the control current. This affects the proportional valve's control element (e.g., a spool-bushing pair or seat plug), causing its displacement (x). If the valve uses a proportional solenoid with a regulated force, the F force is generated proportionally to the value of the control current on the armature of the proportional solenoid. This force is transmitted to the valve control, which is usually a poppet in the proportional pressure valve. With a change in the displacement of the valve actuator or the force acting on it, the Q flow rate or the p pressure functionally vary depending on if a proportional valve is controlling the flow rate or pressure. These parameters control the operation of the hydraulic receiver where the current determines the n or v speed of the hydraulic receiver and the pressure determines the external M or F load. An analog electrical signal with a voltage value typically not greater than 10 V is fed to an electronic amplifier. From the electronic amplifier, the electric control signal is fed through wires with a current that usually does not exceed 1.5 A per coil of a proportional electromagnet. Depending on the type of proportional electromagnet, a force or displacement of the proportional electromagnet armature is generated. If the valve uses a proportional solenoid with an adjustable stroke, a displacement of the electromagnet armature is generated proportionally to the value of the control current. This affects the proportional valve's control element (e.g., a spool-bushing pair or seat plug), causing its displacement (x). If the valve uses a proportional solenoid with a regulated force, the F force is generated proportionally to the value of the control current on the armature of the proportional solenoid. This force is transmitted to the valve control, which is usually a poppet in the proportional pressure valve. With a change in the displacement of the valve actuator or the force acting on it, the Q flow rate or the p pressure functionally vary depending on if a proportional valve is controlling the flow rate or pressure. These parameters control the operation of the hydraulic receiver where the current determines the n or v speed of the hydraulic receiver and the pressure determines the external M or F load.
There are no significant differences in terms of the mechanical design between conventionally controlled valves and proportionally controlled valves. The main difference relies on the fact that, in proportional valves, instead of a spring, hand wheel, lever, or conventional solenoid, there are one or more proportional solenoids. There are two types of proportional electromagnets: force-adjustable and stroke-adjustable. Their use depends on the type and function of the proportional valve and is determined by its characteristics. The benefits of using proportional valves include combining several functions in one valve, a smooth control of flow and pressure parameters, the ability to program the valve and the receiver controlled by it, and the reduction of dynamic surpluses [2]. Proportional valves also have some disadvantages in comparison to conventional ones, including a higher price, stricter requirements regarding purity of the working liquid, and sensitivity to operating conditions (moisture, salinity of the environment, and external mechanical vibrations). Technically speaking, proportional valves were initially a bridge between conventionally (mechanically or electrically) controlled valves and servo valves. Presently, the latest proportional control valves have dynamic parameters equal to those of servo valves and sometimes even surpass them [3]. For example, it can be stated that the limit frequency of a two-stage servo valve is 240-270 Hz depending on the manufacturer, and the limit frequency of the latest-generation proportional valve with VCD (Voice Coil Drive) technology is 350 Hz. Several years ago, this frequency was 6-10 Hz for a singlestage proportional valve [4]. In many industries, particularly mechanical engineering, proportional relief valves with one-or more often, two-stages are widespread. An example of the use of proportional valves is their use in the hydraulic system for lifting and lowering loads with significant masses [5] such as agricultural machinery [6], CNC machine tools, hydraulic presses, wheel loaders, and ships. There are no significant differences in terms of the mechanical design between conventionally controlled valves and proportionally controlled valves. The main difference relies on the fact that, in proportional valves, instead of a spring, hand wheel, lever, or conventional solenoid, there are one or more proportional solenoids. There are two types of proportional electromagnets: force-adjustable and stroke-adjustable. Their use depends on the type and function of the proportional valve and is determined by its characteristics. The benefits of using proportional valves include combining several functions in one valve, a smooth control of flow and pressure parameters, the ability to program the valve and the receiver controlled by it, and the reduction of dynamic surpluses [2]. Proportional valves also have some disadvantages in comparison to conventional ones, including a higher price, stricter requirements regarding purity of the working liquid, and sensitivity to operating conditions (moisture, salinity of the environment, and external mechanical vibrations). Technically speaking, proportional valves were initially a bridge between conventionally (mechanically or electrically) controlled valves and servo valves. Presently, the latest proportional control valves have dynamic parameters equal to those of servo valves and sometimes even surpass them [3]. For example, it can be stated that the limit frequency of a two-stage servo valve is 240-270 Hz depending on the manufacturer, and the limit frequency of the latest-generation proportional valve with VCD (Voice Coil Drive) technology is 350 Hz. Several years ago, this frequency was 6-10 Hz for a single-stage proportional valve [4]. In many industries, particularly mechanical engineering, proportional relief valves with one-or more often, two-stages are widespread. An example of the use of proportional valves is their use in the hydraulic system for lifting and lowering loads with significant masses [5] such as agricultural machinery [6], CNC machine tools, hydraulic presses, wheel loaders, and ships.

Related Work
The research to determine the importance of hydraulic valves' design and/or operational parameters is still ongoing. For example, study [7] optimized the relief valve by minimizing partial multi-valued logic functions. Multi-valued logical equations which constituted design guidelines for the entire series of types of such valves were used. The analysis of the stability of hydraulic elements based on the systems of multi-valued logical equations and the method of multi-valued logical trees, taking into account weighting factors, allows for the consideration of the conditions of global stability. The most favorable result is the relationship specification, which binds the design and operational parameters limitations. In addition, the conditions that limit the parameters of the valve and the system are brought to a simple analytical and graphical relationship. Overall, it is limited to a The modeling of hydraulic systems usually uses ordinary differential equations. From these equations a system of equations is created, initial conditions are assumed, and, after parameterization of the equations, the equations are solved to obtain time courses of the relevant parameters of the hydraulic system, e.g., pressure as a function of time or the velocity of the receiver as a function of time. These are models with focused parameters. One hydraulic system sometimes exhibits wave phenomena that lead to hydraulic resonance. This is referred to as a hydraulic long line. In such cases, partial differential equations are used for modeling and the method of characteristics (MOC) is used to solve them. This paper considers a system in which no wave phenomena occur and uses ordinary differential equations to describe the valve state. Decision-support systems are also applied to hydraulic and pneumatic systems [8][9][10][11][12][13][14][15][16]. The paper [8] mainly presents related methods, from classical clustering and classification topics to database methods (e.g., association), and from database methods (e.g., association rules, data cubes) to newer and more advanced topics (e.g., SVD/PCA, wavelets, support vector machines). The work of [9][10][11][12] focused on concepts for integrating decision-support systems of poorly structured data with a data warehouse based on relational or multidimensional structures. In [13], a framework was developed to evaluate different rainwater-discharge options for urban areas in arid regions. The modeling of rainfall runoff was carried out using the Hydrological-Engineering-Centre and Hydrological-Modelling-System (HEC-HMS). Hydraulic modeling was carried out using SewerGEM to evaluate the effectiveness of the various alternatives for a given design flood [14,15]. The authors of [14,15] presented further applications of multi-criteria decision support methods. In particular, in the work [15] of the Geospatial Information System (GIS), a multi-criteria decision-making system (MCDM) was applied to logic. The decision-making Trial and Evaluation Laboratory (DEMATEL) approach was used to create a network of relationships between criteria. The author of [17] described a model-driven decision-support system (software tool) implementing a model-based online leak-detection and localization methodology that is useful for a large class of water distribution networks.
The present work presents the use of multi-valued logical trees with multivalent weighting factors in the analysis of a two-stage proportional relief valve and a nozzleaperture preliminary stage [18][19][20]. A significant amount of literature exists on the applications of decision trees in decision-making systems. However, there are only a small number of publications on their application in design methodology. Cognitive decision theories seek sufficient and effective solutions for so-called real-world problems and welldefined problems. There are a number of decision-support methods that are familiar to the authors and, in particular, have already been used by the authors to solve a number of problems in decision-support areas, e.g., in the use of special types of parametric dependency graphs [21,22]; inductive decision trees [23,24]; and in particular multi-valued logic trees [25]. Specifically, the recent paper has shown how methods based on multivalued logic trees can be very beneficial when other methods are ineffective. However, multi-valued logic tree methods have plenty of advantages in design methodology and are still being developed. The advantage of the method of multi-valued logical trees is that the measurement data can be recorded by means of appropriate formal notations and it is even possible to combine complex quantitative and qualitative features with different degrees of detail according to the rules of the multi-valued morphological array. The canonical alternative normal form (KAPN) of a bivariate or multi-valued logical function describes all variants, i.e., true (realizable) solutions of a given problem obtained according to the rules of the morphological table, as the full array of combinations of values of logical variables describes all theoretical variants. As a result of minimization (after applying the Quine-McCluskey algorithm), one obtains from the realizable solutions the true sub-solutions as a shortened alternative normal form of the SAPN of the logical function. In this way, the real sub-solutions of the problem are appropriately grouped and therefore the computational time required to obtain the most important real sub-solutions is reduced.

The Tested Object
The tested object is a two-stage proportional relief valve with a preliminary nozzleaperture stage ( Figure 2) [23]. Figure 2 demonstrates a two-stage proportional relief valve. The main stage is pressure controlled, while the pilot is controlled by a proportional solenoid. Changing the pressure in the chamber above the main-stage spool is possible by throttling the fluid flowing out of the pilot. This throttling is altered by changing the position of the diaphragm driven by the proportional solenoid.
ical alternative normal form (KAPN) of a bivariate or multi-valued logical function describes all variants, i.e., true (realizable) solutions of a given problem obtained according to the rules of the morphological table, as the full array of combinations of values of logical variables describes all theoretical variants. As a result of minimization (after applying the Quine-McCluskey algorithm), one obtains from the realizable solutions the true sub-solutions as a shortened alternative normal form of the SAPN of the logical function. In this way, the real sub-solutions of the problem are appropriately grouped and therefore the computational time required to obtain the most important real sub-solutions is reduced.

The Tested Object
The tested object is a two-stage proportional relief valve with a preliminary nozzleaperture stage ( Figure 2) [23]. Figure 2 demonstrates a two-stage proportional relief valve. The main stage is pressure controlled, while the pilot is controlled by a proportional solenoid. Changing the pressure in the chamber above the main-stage spool is possible by throttling the fluid flowing out of the pilot. This throttling is altered by changing the position of the diaphragm driven by the proportional solenoid.      The receiver in the analyzed system is a throttle valve whose performance characteristics are described as follows:    p ≤ 1 MPa : Q odb = 1.2446666 · 10 −10 p, 1 MPa ≤ p ≤ 6 MPa : Q odb = 0.3533333 · 10 −10 p + 0.8913333 · 10 −4 , p ≤ 6 MPa : Q odb = 0.2425893 · 10 −10 p + 14.55 · 10 −5 . (1) where Q odb is the hydraulic actuator flow rate.
In order to describe the flow through a proportional valve it is necessary to consider the value of the loss factor as a function of the displacement of the moving element. The actual course is similar to the solution of a second-order differential equation with a variable throttling factor. This relationship is described in the following form [26]: The following course was used in the control stage (∆y < 0): The force generated by the electromagnetic transducer used in the valve is described as follows: where T m = 15 ms when i increases U 18 − i > 0 and T m = 7.5 ms when i decreases

Mathematical Model of the Tested Valve
The mathematical model of the valve under consideration was built on the basis of ordinary differential equations of the second order. The first equation of the system of equations is the flow rate balance equation, which takes into account the compressibility of the working fluid (its capacitance).
The flow balance of the drive system can be written as [23]: The flow balance through the main valve stage is described as: The flow through the nozzle is described as: The flow balance through the control stage is described as: In addition, the flow rate is distinguished in the main stage as: and at the control stage as: The flow rate through the valve is represented as: • through the main stage: • through the control stage: Ultimately, the flow rates are represented as: • through the nozzle D 1 : • through the nozzle D 3 : An additional equation described is the equilibrium equation of the forces acting on the valve control element (according to d'Alembert's principle) on the main stage and the secondary stage. This equation takes into account the forces of inertia, spring stiffness, frictional force, and the hydrodynamic reaction force associated with the change in momentum of the fluid stream.
Forces in the valve: Dynamic loads: In the main stage: where the following values indicate: • 0.0675-the mass of the main stage slider; • 0.008-spring mass; • 0.00439-the mass of the associated liquid.
In the control stage: Sticky friction: • forces in the main stage: • forces in the control stage: Forces of the hydrodynamic reaction are described as follows: • of the main stage: • of the nozzle-aperture pair: The dynamic equations of the proportional valve forces are described at any point in the transient state after the introduction of the step function: of the main stage: of the control stage: The feedback loop equation is written as follows: when U z − U p − e 0 < 0, if none of these conditions are met: The output equations for the computer simulation of the operation of the hydraulic part are shown in Appendix A.

Methodology of Multi-Valued Logic Trees with Weight Coefficients as Discrete Optimization
The methodology presented is based on two algorithms:
The Quine-McCluskey method of minimising partial multi-valued logical functions with multi-valued weighting factors.

Quine-Mc Cluskey Algorithm for the Minimization of Partial Multi-Valued Logical Functions
In the case of logic trees, the logical values of the variables are encoded on the branches of the tree. There can only be one Boolean variable per level of the tree, with the number of floors being equal to the number of independent variables of a given Boolean function. Representing a given Boolean function written in canonical alternative normal form (KAPN) on a logic tree involves encoding the individual canonical products on a tree path from the root to the end vertex. An individual path on the tree (from root to vertex) is a component of the singularity of the logical function, describing the realization of one possible solution. On the contrary, the set of paths is the set of all possible solutions. Figure 4 shows a logic tree in which a fixed Boolean function of three variables is encoded.    A minimized form of the output function (with a minimum number of literals) is subsequently obtained. However, given that so-called isolated branches exit, this is not the minimum decision form, meaning that there is no continuity between the root and the vertices. In the case of multi-valued logical functions-as in Boolean functions-the notions of incomplete gluing and elementary absorption, which are applied to the APN of a given logical function, play a fundamental role in the search for prime implicants.  A minimized form of the output function (with a minimum number of literals) is subsequently obtained. However, given that so-called isolated branches exit, this is not the minimum decision form, meaning that there is no continuity between the root and the vertices. In the case of multi-valued logical functions-as in Boolean functions-the notions of incomplete gluing and elementary absorption, which are applied to the APN of a given logical function, play a fundamental role in the search for prime implicants.
A gluing operation is called a transformation: and A denotes an elementary partial product, of which, variables of the individual literals belong to the set 1 , … , − , + , … , . An incomplete gluing operation is called a transformation:  A gluing operation is called a transformation: where r = 1, . . . , n and A denotes an elementary partial product, of which, variables of the individual literals belong to the set {x 1 , . . . , x r−i , x r+i , . . . , x n }. An incomplete gluing operation is called a transformation: where r = 1, . . . , n and A denotes a partial product of which the variables of the individual literals belong to a set of {x 1 , . . . , x r−i , x r+i , . . . , x n }. An elementary absorption operation is called a transformation: where 0 ≤ u ≤ m r − 1, 1 ≤ r ≤ n, and A denotes a partial product of which the variables of the individual literals belong to a set of {x 1 , . . . , x r−1 , x r+1 , . . . , x n }. If the above equation holds, then A absorbs Aj u (x r ). Signs (v) denote that a given partial product of the elementary, written using the digits of the system (m 1 , . . . , m n )-positional, takes part in the gluing with those products that have a sign (v) in the same column. The notation marks of the gluing operation are entered separately in the columns and not in a single column as was the case in previous literature studies of bivalent cases. In the case of equal, multivalued variables x 1 , . . . , x n of a given logical function, the set of first implicants is obtained as a special case from different multi-valued variables.
Example 1. Using the relationship: The successive steps of minimizing a multi-valued logical function can be represented as follows: The rank of importance of successive decision variables is determined using complex alternative normal forms through the swapping of floors in logical decision trees. The swapping of logical tree floors in complex, multi-valued logical functions establishes the rank of importance of logical variables from the most important (at the root) to the least important (at the top). There is a generalization of a bivariate quality indicator to a multivariate one; (C k − k i m i ) + (k i + K i ), where C k represents the number of branches of the k-th floor, k i is the simplification factor on the k-th floor of the m i -value variable, and K i represents the number of branches (k − 1)-th floors from which the non-simplifying branches of the k-th floor are formed. In this way, it is possible to obtain the minimum complexity alternative normal form (MZAPN) of a given logical function without isolated branches on the decision tree and with a concomitant minimum number of real (realizable) branches, which in particular can be considered to be elementary design guidelines. All transformations refer to the so-called Quine-McCluskey algorithm for minimising individual partial multi-valued logical functions.  Figure 6 shows all possible ZKAPNs of a given multi-valued logical function.  Figure 6 shows all possible ZKAPNs of a given multi-valued logical function. 1. The first stage of minimization due to x3: 1.

Generalization of the Quine-Mc Cluskey Algorithm for Minimization of Partial Multi-Valued Logical Functions for Multi-Valued Weighting Factors
In multi-valued logical functions with weighted products it is possible to apply the Quine-McCluskey algorithm for the minimization of multi-valued functions. As with the minimization of multi-valued logical functions without weighting coefficients, in the algorithm the elementary products are written as numbers in the corresponding positional systems. Additional elements and operations are introduced to account for the weighting coefficients.
In partial data of multi-valued logical functions f i (x 1 , . . . , x n ) n of variables (m 1 , . . . , m n ), value-added gluing and pseudo-gluing operations should include weighting factors (w n , w n−1 , w n−2 , . . . , w 1 ) assigned to the corresponding multi-valued logical products.
The Quine-McCluskey algorithm for minimizing multi-valued logical functions is built from n columns with (w 1 , . . . , w n ) weighting factors.
Symbols indicating pseudo-gluing (V) and gluing (v) sequentially relative to groups of indices differing by one are placed in the columns corresponding to the values of the weighting factors for the corresponding logical products.
Given multi-valued weighting coefficients, individual (parallel) pseudo-bonding operations sequentially against groups of indices, differing by at least one, and containing at most (m i − 1) elements can proceed in canonical products with different weighting coefficients.
The characters appear in different columns. In addition, they may be in columns with a corresponding coefficient (w n , w n−1 , w n−2 , . . . , w 1 ). Therefore, the columns with (w 1 , . . . , w n ) weighting coefficients introduce position numbers p i , with i = 1,...,n, which is useful for calculating the quality of the minimization in further stages.
Definitions of 'pure' and 'impure' gluing are introduced for gluing operations of individual partial multi-valued logical functions with weighted coefficients. A pure gluing operation is a transformation of: where r = 1, . . . , n and A represents a partial product of which the variables of the individual literals belong to a set of {x 1 , . . . , x r−i , x r+i , . . . , x n }. In n m-value variables, the weighting factor before the partial canonical product takes values in the interval w 1 , . . . , w n , with w j = w j−1 + w j−2 + . . . + w 1 and j = 2, . . . , n. The impure gluing operation for multi-valued canonical elementary products is performed with respect to the weighting factor with the smallest value, i.e., min{w 1 , . . . , w n }.
An impure bonding operation is a transformation: where r = 1, . . . , n, w s > min{w 1 , . . . , w m r −1 }, and A denotes a partial product of which the variables of the individual literals belong to a set of {x 1 , . . . , x r−i , x r+1 , . . . , x n }. In n(m 1 , . . . , m n )-value variables, the weighting factor w i before the partial canonical product takes values in the interval w 1 , . . . , w n , with w 1 where j = 2, . . . , n.

Definition 3.
An incomplete gluing operation is a transformation that retains the original records to be glued after the algorithm has been executed in the result.
Given that there is an isomorphic interpretation of logical transformations, the Quine-McCluskey algorithm for minimizing individual partial multi-valued logical functions can be considered with the weighting factors mentioned, which is important for describing the rank validity of design guidelines.
The following are the successive steps in the minimisation of logical functions due to given decision variables: Tree interpretation. Figure 7 shows the MZAPN of the multi-valued logical function from Example 3.
Tree interpretation. Figure 7 shows the MZAPN of the multi-valued logical function from Example 3. The proposed methodology can be described by the flow chart shown in Figure 8. The proposed methodology can be described by the flow chart shown in Figure 8. The structuring of the described problem takes into account the methodology of multi-valued logical trees, allowing for the introduction of appropriate formal notations and even making it possible to combine complex quantitative and qualitative features with different degrees of detail according to the principles of a multidimensional morphological array. Therefore, there is no need to extend the generation process to sub-arrays when using a multidimensional morphological array as all information about the varieties of the main and detailed features and their numerous modifications can be immediately The structuring of the described problem takes into account the methodology of multivalued logical trees, allowing for the introduction of appropriate formal notations and even making it possible to combine complex quantitative and qualitative features with different degrees of detail according to the principles of a multidimensional morphological array. Therefore, there is no need to extend the generation process to sub-arrays when using a multidimensional morphological array as all information about the varieties of the main and detailed features and their numerous modifications can be immediately stored in this array and marked on the variant tree.
In addition, the morphological and decision tables can be encoded analytically and numerically according to the definitions and theorems of the logic of multi-valued decision processes. This enables a variant way of identifying and classifying information in computer science terms when seeking and modifying solutions in the design process.
In such a situation it is possible to introduce CAD, e.g., for the generation of all theoretical variants of the designed system, selection for realizability, search for realizable solutions and-most importantly-realizable sub-solutions, etc. In order to ensure the stable operation of the actual system, model tests are carried out on the basis of which of the relevant parameters are selected. The phenomena occurring during the flow of a medium are quite often not precisely defined, so it is necessary to identify an analytical model when carrying out such studies.

Application of the Methodology of Multi-Valued Logic Trees with Weighting Factors in the Optimization of a Proportional Valve
Tests have already been carried out for valves of the direct-acting UPZ type [4], which are designed to regulate the upstream pressure of steam and non-flammable, chemically inert gases and liquids regardless of the pressure at their outlet. Multivalent weighting factors were not considered in the tested valve class. For this reason, it was decided that an improved Quine-McCluskey algorithm with weighting factors for the hydraulic proportional valve would be used. Therefore, three 'novelties' are presented in this paper.
One of the optimal methods presented are multi-valued logic algorithms. For example, in [24] the authors presented a description of the dynamics of molecular states caused by a sequence of laser pulses using multi-valued logic. In turn, the authors of [25] used multivalued logical schemes to calculate significance measures based on incompletely-defined data. This method is based on the definition of a mathematical model of an analyzed system in the form of a structure function that determines the correlation of the system reliability and the states of its components.
In [27], the authors described the historical and technical background of MVL, as well as the areas of present and future applications of quadrivalent logic. It was also intended to serve as a guide for non-specialists. The wide application of multi-valued logic in particular in these microelectronic circuits is presented in [28]. Additionally, there are many original works describing the practical application of multi-valued logic trees.
In addition, there are other works in which multi-valued decision trees and logic algorithms have been applied. For example, the authors of [29] presented the applications of machine learning and classification and regression trees (CART) in medicine. Specifically, they presented the concept of a gradient-boosting algorithm. The authors of [30] presented the application of a rotation forest with decision trees as a base classifier and a new ensemble model in the spatial modeling of groundwater potential. The use of fault-tree analysis to calculate system-failure probability bounds from qualitative data in an intuitive, fuzzy environment is presented in paper [31]. Meanwhile, in paper [32] the authors adopted component fault trees (CFTs) to support fault tree analysis, failure mode, and effect analysis as extensions of SysML models. Boolean decision support methods were presented in paper [33]. A very modern optimization method was proposed by the authors of [34]: the use of root trees. The root-tree algorithm was used for high-order sliding mode control using a super-twist algorithm based on the DTC scheme for DFIG.
The initial conditions of a differential equation can be determined by entering dx i dt = 0. The simulations were performed using the Matlab/Simulink package: Assuming that: It can be obtained that: The Importance of the Design and/or Operational Parameters of a Hydraulic Proportional Valve In the optimization process, the changed parameters of the proportional valve while observing the Q flow rate and p pressure are represented by the regulator K p1 ·K p2 gain (as a complex variable), the Q obd receiver flow rate (depending on the impulse input of the U z control voltage), and the F m magnetic force.
The arithmetic values of the tested parameters were selected for the analysis. They were coded by the authors of this work with logical decision variables: Q rz = 24 → 12 dm 3 /min ∼ 1; Q rz = 36 → 12 dm 3 /min ∼ 2; (54) In the operation of the relief valve, the authors introduced restrictions on the Q and p design parameters in terms of the stabilization time t w :t w < 0.48 t 0 . Subsequently, dynamic calculations of the valve were carried out, resulting in the t w :t w < 0.48 t 0 limitation. Following the dynamic calculations, 23 charts were selected. The code changes of the K p1 , K p2 , Q rz and F m design parameters are presented in Table 1.
Furthermore, in the code changes of the K p1 ·K p2 , Q rz and F m design parameters multivalued w i weighting factors are introduced, similar to the relief valve. The greater the weighting number, the faster the Q and p functions reach a stable state (t t > t j ).
The following weighting factors were adopted in the t w < 0.48 t 0 limitation: Table 1 presents the code changes of the K p1 ·K p2 , Q rz and F m design parameters, taking into account the multi-valued weighting factors and the t w < 0.48 t 0 limitation.
Notably, the value of the weighting factor for changes in the code K p1 ·K p2 , Q rz and F m design parameters in Table 1 is minimal among the coefficients defined separately in the Q and p function. If one of the functions stabilizes faster than the other, then the canonical product for the same code changes of K p1 ·K p2 , Q rz and F m parameters should be assigned a smaller weighting factor.  2  2  1  2  2  2  3  2  2  2  2  1  2  2  2  2  2  1  2  1  3  3  0  2  1  1  0  2  2  0  2  1  1  0  2  2  2  0  2  0  3  0  3  1  2  0  3  2  3  0  1  1  1  0  1  2  3  0  1  0  2  0  0  1  1  0  0  2  1  1  2  2  2  1  1  2  1  1  3  2  1  3  2  2  1  3  1  2  3 3 3 2 In the system of multi-valued logic functions with weighting factors, weighting factors are assigned separately for each of the functions.  2  1  1  2  2  2  1  1  2  1  1  3  2  1  3  2  2  1  3  1  2  3 3 3 2 In the system of multi-valued logic functions with weighting factors, weighting factors are assigned separately for each of the functions.     The multi-valued logical trees with the weighting factors from Table 1 are shown in Figure 14. The multi-valued logical trees with the weighting factors from Table 1 are shown in Figure 14.  For the t w < 0.48 t 0 criterion of limitation, one optimal multi-valued logical tree is presented in Figure 14. For a hydraulic proportional valve, the most crucial parameter is the Q odb flow rate of the receiver (depending on the U z step function of the control voltage).
One of the issues presented in this paper is the application of Boolean equations in the optimization of machine systems. This paper generalizes the Quine-McCluskey algorithm for minimizing multi-valued logical functions with multi-valued weight coefficients. In addition, a procedure for the combinatorial solution of weighted multi-valued systems of logical equations describing design guidelines in terms of morphological analysis with the Rosser-Turguette axioms is discussed.
The application of the methodology of multi-valued logic trees with weighting coefficients for relief valves allows for the determination of alternative sets of design guidelines to find the most crucial design guidelines in any fixed design and/or operational parameters while ensuring that the constraints and extremes of the criterion are met. In particular, the novelty presented in this paper is: • the development of a Quine-McCluskey algorithm for minimizing logical functions with weighted coefficients; • the development of algorithms for multi-valued Boolean equations, allowing for the obtainment of the most important design guidelines in the process of optimization of machine systems; and • carrying out calculations on a real object (relief valve) to confirm the developed algorithms.

Conclusions
This paper presents the use of multi-valued logical trees with weighting factors to determine the importance of the constructional and operational parameters of a two-stage proportional relief valve. As has been demonstrated by this research, relief valves do not keep up with a pressure increase in the system, react with a certain delay, and can vibrate under fixed operating conditions.
The above incorrect response of the valves usually occurs during the transition period. Hence, it is necessary to carry out model tests of valves in the transition state and to determine the importance of the operational parameters directly affecting their dynamics. Model tests aim to select essential parameters to ensure the stability of the real system. It is crucial to determine the importance of design and/or operational parameters during the model verification and subsequently select the appropriate optimization procedure.
This work discusses the procedure of a combinatorial solution for weight-multi-valued systems of logic equations describing the design guidelines in terms of morphological analysis with the preservation of Rosser-Turguette axioms. It has been shown that, in general, the minimization of logic functions with weight coefficients may be the same as without weight coefficients. However, a better reflection of the physical models of hydraulic relief systems was obtained through mathematical models. The literature shows that various coefficients of logical products have not been taken into account in the separable and common minimization of systems of multi-valued logic equations.
Three following 'novelties' are presented in this paper: 1.
The Quine-McCluskey algorithm was applied to completely new proportional valves previously not tested using such methods; 2.
Multivariate weighting factors were taken into account, which allowed for the introduction of multi-valued weighting factors to determine the most important design guidelines. This made it possible to apply the method of multi-valued decision trees to solve the technological problems of the studied proportional valve series. It also made it possible to use systems of logical equations as a formal decision-making description of the study of the importance rank of changes in design and operating parameters as a method for optimizing relief valves and other hydraulic systems; 3.
For the optimization process of hydraulic valve design, an approach using CFD computer simulation methods is encountered. However, any results obtained from simulation solutions should be, at least to some extent, verified by reliable experimental studies.

Limits of the Methodology Used
Each of the KAPN products should be assigned corresponding discrete changes in parameter values. Therefore, it is not possible to fully apply the developed methods in continuous linear and non-linear optimisation. However, for construction/engineering purposes, the use of discrete analysis is preferable in the opinion of the authors (who have computational experience). If one were to change the numerical values of the input variables in a mathematical model, one would obtain changes in the values of the output variables. In order to obtain a different planned behaviour of a system (component), one can often make many changes to the numerical values of the input variables. Concerns related to change include: which values might be changed, how the change might be made (by increasing values, keeping them unchanged, or decreasing them), or in what order the variables might be changed, etc. Such conjecturing is akin to subjectively (according to a given designer) changing the numerical values in a mathematical model. This means that another designer, according to their own experience, may subjectively redesign the layout (element) quite differently for new work conditions that are identical to those of the previous designer.
The multi-valued weighting system of logical equations describing the design guidelines can be minimized separately or together with logical equivalence. Still, even in the bivalent (Boolean) case, the common minimization is not inferior to the separate minimization in terms of literal multiplicity. Increasing, reducing, or keeping the numerical values unchanged in the process of redesigning a system for other operating conditions can be coded using multi-valued logic while sets of design guidelines can be presented as sums of multi-valued logical products.
Model tests are particularly important in the design of new valves. These design parameters, which significantly affect the dynamics of valves, cannot be selected randomly (depending on the assumptions and experience of the designer). Their values should be closely related to the permissible peak overload of the controlled signal, operation speed, time constant, and eliminating vibration. Model tests will be more useful if the described valve is outlined more accurately in the transition state. Thus, building a correct analytical equation that presents a given valve in a dynamic course determines the sense of any theoretical considerations.
In further research, it will also be necessary to take into account the modified methodology of multi-valued logic trees as parametrically-playing out graphs, i.e., a heuristic simulation method for solving linear-dynamic decision models for relief valves. In the instance analyzed, a number of simplifications were additionally taken into account; for example, the impact of the closing element hitting the valve seat was not considered and the effect of the valve-wall compressibility was not taken into account. These factors will be considered in further papers. Additionally, a new control system using optimised proportional-directional control valves, throttling valves, and flow controllers will be proposed in further studies.

MZAPN
minimal complex alternative normal form; ZKAPN complex canonical alternative normal form; f (x 1 , x 2 , x 3 ) multi-valued logic function of three variables; w i weighting factor for multi-valued logical products; A partial elementary product; k vx , k vy loss factors in the control stage of a hydraulic proportional valve; F m electromagnetic force; p operating pressure; Q p actual pump capacity; Q D1 , Q D2 , Q D3 flow rate through the D 1 , D 2 and D 3 nozzles in a proportional valve; Q odb receiver flow rate in a proportional valve; K p1 K p2 regulator boost in a proportional valve; Q t theoretical pump capacity (for a fixed hydraulic system); Q zp flow rate through a relief valve; R p0 the resultant internal leakage resistance in a system; R pp the leakage resistance of a pump determined by means of volumetric efficiency; R pz the leakage resistance of a relief valve determined by means of the slope of the valve static characteristic for pressures below the opening pressure; R r the leakage resistance of a distributor; R s the leakage resistance of a motor determined by means of volumetric efficiency.