Next Article in Journal
On Correspondences between Feedforward Artificial Neural Networks on Finite Memory Automata and Classes of Primitive Recursive Functions
Next Article in Special Issue
Characterization Results of Solution Sets Associated with Multiple-Objective Fractional Optimal Control Problems
Previous Article in Journal
A Variational Model for Wrapped Phase Denoising
Previous Article in Special Issue
Appearance-Based Gaze Estimation Method Using Static Transformer Temporal Differential Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Robust Tracking Control of Injection Velocity in an Injection Molding Machine

1
Guangdong Key Laboratory of IoT Information Technology, School of Automation, Guangdong University of Technology, Guangzhou 510006, China
2
Guangdong-HongKong-Macao Joint Laboratory for Smart Discrete Manufacturing, Guangdong University of Technology, Guangzhou 510006, China
3
Key Laboratory of Intelligent Detection and the Internet of Things in Manufacturing (GDUT), Ministry of Education, Guangzhou 510006, China
4
Guangdong Laboratory of Artificial Intelligence and Digital Economy (SZ), College of Mechatronics and Control Engineering, Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2619; https://doi.org/10.3390/math11122619
Submission received: 2 March 2023 / Revised: 19 April 2023 / Accepted: 5 June 2023 / Published: 8 June 2023

Abstract

:
Injection molding is a critical component of modern industrial operations, and achieving fast and stable control of injection molding machines (IMMs) is essential for producing high-quality plastic products. This paper focuses on solving an optimal tracking control problem of the injection velocity that arises in a typical nonlinear IMM. To this end, an efficient optimal robust controller is proposed and designed. The nonlinear injection velocity servo system is first approximately linearized at iteration points using the first-order Taylor expansion approach. Then, at each time node in the optimization process, the relevant algebraic Riccati equation is introduced, and the solution is used to construct an optimal robust feedback controller. Furthermore, a rigorous Lyapunov theorem analysis is employed to demonstrate the global stability properties of the proposed feedback controller. The results from numerical simulations show that the proposed optimal robust control strategy can successfully and rapidly achieve the best tracking of the intended injection velocity trajectory within a given time.

1. Introduction

Plastic products are extensively used in our daily life and industries, such as in the automobile industry, aerospace, military technology, and other fields, due to their good plasticity, chemical stability, and lightweight [1,2]. According to market statistics, the total plastic production in the world in 2021 was 399 million tons (384 million tons in 2020), an increase of 3.7% over 2020 (the same increase of 4.6% in 2020). Data shows that plastic products have become an indispensable part of all walks of life and people’s daily life. From now until 2030, the injection plastic market is anticipated to develop at a compound annual growth rate of 4.90%.
Injection molding (IM) technology has the advantages of high production speed, high efficiency, process automation, and the ability to manufacture complex shaped parts [3,4,5]. Currently, due to its dependable, high-quality performance, IM has become one of the most widely commonly applied processes for producing plastic components. IM is a sophisticated and continuous cycle operation, wherein various requirements for the IMM (as shown in Figure 1) settings for different products and stages are required. The quality of plastic products is directly connected to processing facilities, product cooling time, and subsequent processing technology. Poor management of variables such as temperature, pressure, and molding cycle will result in poor quality of the final molded plastic and affect their use [6,7,8]. However, as we are now in the field of intelligent manufacturing, the overall automation level of IM equipment is still very low, and the vast majority of process parameters still need to be adjusted manually, which often leads to unstable product quality and poor product consistency, which requires an urgent need for technological innovation to improve competitiveness. The development of related technologies, such as plastic molding, will have a profound impact on the competitiveness of the discrete manufacturing industry in the future and will also have an important radiative impact on the industrial sector. Therefore, the intelligent automation of various IMMs has become a hotspot in the study of discrete industrial manufacturing [9,10,11,12,13,14,15,16,17,18,19].
Recent years have seen a large number of pertinent studies on control and optimization methods for the temperature, pressure, and injection velocity in various IMMs [20,21,22,23,24,25,26,27,28,29,30,31,32,33]. For instance, Cho et al. [20] presented an open-loop control approach to realizing the pre-position tracking control of the melt injection flow inside an IMM. Dubay [22] developed an all-encompassing learning and self-optimizing model predictive control approach to regulate the melt temperature in an IMM. Yao et al. [23] studied a time-optimal control problem for the barrel temperature regulation in an IMM. Ruan et al. [29] proposed an intelligent temperature compensation control strategy for the injection molding process under dynamic conditions using feedforward iterative learning control and the Q-learning method. Stemmler et al. [30] proposed a norm-optimal iterative learning controller to track the desired cavity pressure value of an IMM throughout the cycle. Guo et al. [31] proposed a novel reinforcement learning framework and a self-predictive artificial neural network model to set appropriate process parameters. Tian et al. [32] proposed a double-controller scheme for the injection velocity control in thermoplastic injection molding. Yang et al. [33] proposed an improved PID method using the unsaturated integral and passivation differential to achieve injection velocity control. Although these control methods have proven good dynamic performance, on the one hand, most methods focus on open-loop control of IMMs and often have problems such as poor robustness in practical applications. On the other hand, some of the above-mentioned methods are more complex and often involve tedious calculations that are difficult to implement in practical physical environments. Therefore, it is very meaningful to design more efficient optimization and control methods to achieve quality control in the injection molding process. Motivated and inspired by the above, in this paper, we carry out further research on the velocity control approach for a typical type of injection machines, and we will consider designing a robust optimal feedback controller to realize the optimal tracking control of the intended injection velocity trajectory.
In an IMM, in addition to controlling the injection temperature and injection pressure, controlling the injection velocity has also been recognized as a key variable in the molten polymer. It indicates the rate of the polymer melt flowing to the mold cavity during the filling phase. The velocity has a significant impact on cavity pressure, which in turn may affect the quality of the product in terms of residual stress, shrinkage, impact strength, product morphology, and surface characteristics. For example, during the high-speed filling process, when the melt is injected into the cavity via the nozzle and gating system, a significant amount of frictional heat is generated, thereby increasing the melt temperature. However, if the cavity is filled at a low injection speed, the filling time will be prolonged. Conversely, if the mold cavity is filled at a high injection speed, the welding of the insert rear may be unsatisfactory, leading to a reduction in product strength. The melt injected into the cavity initially cools down rapidly and becomes more viscous, requiring higher pressure for subsequent filling [34]. Thus, selecting an appropriate injection velocity is also crucial for controlling the properties of injection-molded products.
In this paper, we study an optimal tracking control problem of the injection velocity in a typical class of IMMs with nonlinear characteristics, and an effective optimal robust state feedback controller is proposed and designed to achieve the fast-tracking of the given injection velocity at the given time. In practice, maintaining the actual injection velocity to follow a desired profile can be challenging due to the complex effects of numerous factors on velocity dynamics. Achieving high precision and robust control of injection velocity is therefore critical to improving product quality and reducing cycle-to-cycle variation. Furthermore, due to the nonlinearity of electrohydraulic actuators in the IMMs, the solutions of the related control and the proof of global stability and robustness for such methods remain a challenging task. In conventional methods, control of such systems is a non-trivial problem because of the nonlinearities and underactuation in their dynamics. In this current work, our aim mainly focuses on developing an efficient computational optimal control method using advanced control technology for the fast and robust tracking of injection velocity in a typical IMM. To address this issue, we first formulate a nonlinear tracking control problem for the injection velocity in the injection molding process. Given the nonlinear characteristics and the need for simplified controller design, we use Taylor’s first-order expansion to linearize the original nonlinear dynamical system at the iterative point. This transforms the optimal tracking control problem governed by the original continuous nonlinear system into one governed by a set of simplified linear systems. From a computational optimal control perspective, we propose an optimal feedback robust controller based on solving Riccati equations. This converts the design of a closed-loop controller for linear systems into solving a series of Riccati equations. We then use the Lyapunov theory to prove that our proposed controller has global stability properties. The optimal robust control strategy we propose can efficiently and quickly achieve optimal tracking of the desired injection velocity trajectory within a given time frame. Numerical simulations confirm the feasibility and efficiency of our feedback stabilization controller design. The proposed control method retains the advantages of classical linear optimal control and it has fast and accurate tracking of reference setpoints under moderate variations. This is a first step and quite a preliminary attempt in the intelligent injection molding process manufacturing area. Compared with traditional control methods for IMMs, such as PID and model predictive control (MPC) strategies, our proposed control strategy has some advantages. For example, unlike traditional PID control and MPC in IMMs, the proposed optimal control method is of proven global stability. Moreover, in traditional PID control, the selection of the controller’s parameters usually relies on a heuristic tuning procedure. The approach recommended in this work is simple and efficient, and can also be easily generalized to other types of industrial process optimal control systems.
In summary, the contributions in current work are listed as follows:
(1)
The dynamic mathematical modeling of the speed servo system of the variable speed pump-controlled hydraulic cylinder in a typical IMM is carried out, and the differential plane flatness characteristics of the nonlinear dynamic system are proven.
(2)
The underlying nonlinear dynamic system is high-fidelity linearized using Taylor series expansion, which can help further streamline the design of the ensuing robust controller.
(3)
A stabilizing optimal controller is proposed for the actuator, and it is demonstrated that the suggested feedback controller exhibits properties of global stability based on a rigorous Lyapunov theoretical analysis.
(4)
Sufficient numerical experimental results verify the viability and efficiency of the proposed approach.
The rest of this paper is organized as follows: In Section 2, the dynamics of an IMM is presented. The differential flatness properties of the injection hydraulic system are verified. The optimal tracking control problem of injection velocity is established. In Section 3, the dynamic model of the IMM is approximately linearized by Taylor series expansion and the optimal control problem is solved. In Section 4, the global asymptotic stability of the controller is demonstrated. In Section 5, numerical simulations are conducted to further verify the accuracy of the proposed controller. Finally, a summary and prospect of the current work are given.

2. Problem Formulation

2.1. Mathematical Model of the IMM

Figure 2 shows an organizational layout of the control hydraulic system used in a typical IMM. The entire system consists of a signal measurement control system and a hydraulic system. The hydraulic system is composed of a motor system and a loading system. The output control voltage of the signal measurement control system regulates the speed of the permanent magnet synchronous servo motor (PMSM) through the servo driver. The tiny time constant of the PMSM can typically be described as a first-order inertial element. As a result, the PMSM is reduced to a first-order inertial chain, which can be expressed as
G ( s ) = K v T d s + 1 ,
where K v is the speed gain determined by the characteristics of the servo drive, and T d is a given small constant. The PMSM has a low moment of inertia and quick acceleration and deceleration response times. In the majority of industrial operations, servo motors have been shown to be aperiodic and overdamped. Consequently, the servo motor model can be denoted as:
d n ( t ) d t = K v u ( t ) n ( t ) T d ,
where n ( t ) is the motor speed value, and u ( t ) is the input signal of the servo drive.
The gear pump is rotated by the PMSM, and the output of the gear pump is controlled by the servo motor speed. The existence of a safety valve can prevent high motor speed wear gear pump. Moreover, the input voltage signal u ( t ) complies with the following bound restrictions
u m i n u ( t ) u m a x ,
where u m i n and u m a x are given constants due to the real physical conditions.
Furthermore, the output flow of the gear pump can be formulated as:
q ( t ) = D p n ( t ) C i p F p ( t ) η n ( t ) F p ( t ) D p β e ,
where q ( t ) denotes the output flow, and D p and C i p denote the displacement and leakage coefficient of the quantitative pump, respectively. F p denotes the output pressure, η denotes the viscosity of the oil, and β e denotes the elastic modulus of the oil volume.
The pressure sensor may collect the drive system pressure F p and transmit that information to the signal measurement and the control system. The pressure drop of the high-pressure pipeline can also be disregarded because the hydraulic pipeline is composed of stainless steel. As a result, the flow change of the wheelless chamber is formulated as:
q ( t ) = A 1 y ˙ ( t ) + C e p F p ( t ) η + ( V g + A 1 y ( t ) ) F ˙ p ( t ) β e ,
where A 1 denotes the effective area of the rodless cavity of the hydraulic cylinder, y ( t ) denotes the displacement, C e p denotes the total leakage coefficient, V g denotes the total volume of the high-pressure chamber. The displacement sensor measures the displacement y ( t ) and the corresponding speed y ˙ ( t ) and transmits them to the signal measurement and control system.
The loading system that consists of a three-phase motor and an Atos proportional relief valve is an exact replica of the drive system. The electro-hydraulic loading device, which mimics the linear motion of the hydraulic cylinder, controls the output force of the loading cylinder. As a result, the hydraulic cylinder’s motion force balance can be formulated as:
F p ( t ) A 1 = m t y ¨ ( t ) + B p y ˙ ( t ) + K y + F L ,
where m t denotes the equivalent mass of the piston rod and load, B p denotes the viscous damping coefficient of the drive cylinder, K denotes the elastic coefficient of the load, and F L denotes the external load.
The hydraulic control system under consideration utilizes sensors to gather signals related to pressure, displacement, and velocity. The collected signals are then transmitted to a control system for signal measurement and analysis. The system compares the speed signal in conjunction with the required driving cylinder’s speed signal and generates a control signal accordingly. This control signal is employed to regulate the speed function of the servo driver and the PMSM, based on the input voltage. Eventually, indirect flow control for the hydraulic cylinder and the complete control of the driving hydraulic cylinder are provided by adjusting the electric motor’s rotational speed.
To summarize, we have constructed the evolution process of the speed servo system of the variable speed pump controlling the hydraulic cylinder in an IMM. Based on the above analysis, it can be deduced that the dynamic system is a complex coupled dynamics model with nonlinear characteristics. In order to simplify the subsequent optimal control problem establishment and controller design, we define [ x 1 ( t ) , x 2 ( t ) , x 3 ( t ) , x 4 ( t ) ] = [ n ( t ) , F p ( t ) , y ( t ) , y ˙ ( t ) ] as new state variables of the underlying hydraulic system. Then, the mathematical model of the variable speed pump-controlled hydraulic cylinder speed servo system can be formulated as:
{ x ˙ 1 ( t ) = K v u x 1 ( t ) T d , (7a) x ˙ 2 ( t ) = β e V g + A 1 x 3 ( t ) [ D p x 1 ( t ) ( 1 x 2 ( t ) β e ) ) ( C i p + C e p ) η x 2 ( t ) A 1 x 4 ( t ) , (7b) x ˙ 3 ( t ) = x 4 ( t ) , (7c) x ˙ 4 ( t ) = 1 m t A 1 x 2 ( t ) B p x 4 ( t ) K x 3 ( t ) F L , (7d)
Alternatively, the state-space model (7) can be also expressed as
x ˙ 1 ( t ) x ˙ 2 ( t ) x ˙ 3 ( t ) x ˙ 4 ( t ) = x 1 ( t ) T d β e V g + A 1 x 3 ( t ) D p x 1 ( t ) ( 1 x 2 ( t ) β e ) ) ( C i p + C e p ) η x 2 ( t ) A 1 x 4 ( t ) x 4 ( t ) 1 m t A 1 x 2 ( t ) B p x 4 ( t ) K x 3 ( t ) F L + K v T d 0 0 0 u ( t )
As a result, the dynamic model (8) can be expressed in the following compact form:
x ˙ ( t ) = f 1 ( x ( t ) ) + f 2 ( x ( t ) ) u ( t ) ,
where x ( t ) R 4 × 1 , f 1 ( x ( t ) ) R 4 × 1 , f 2 ( x ( t ) ) R 4 × 1 and u ( t ) R .

2.2. Objective Function

The quality of the final injection-molded products depends on the IMMs’ ability to accurately track the target injection velocity within a specific time T. Therefore, the main aim of this paper is to construct an optimal control signal u ( t ) that enables fast tracking of the target injection velocity while consuming the least amount of energy possible. As a result, the following performance cost function is defined in this study to be minimized:
m i n u ( t ) : J ( u ( t ) ) = 1 2 0 T v ( t ) v d ( t ) 2 + u 2 ( t ) d t ,
where v d ( t ) denotes desired injection velocity during injection molding process; it is a time-dependent state trajectory variable, which can be set in advance according to the actual system operation requirements. As a result, we propose the following Problem P 0 in the IMM.
Problem 1. 
Given a dynamic system (8) of an injection molding hydraulic system with strong nonlinear characteristics, the main aim is to design an optimal robust controller u ( t ) to minimize the objective function (10) to make the speed of the injecting molding v ( t ) track the desired injection velocity output trajectory v d ( t ) within a given time T.

3. Optimal Robust Controller Design

This section aims to devise an optimal robust controller to tackle Problem P 0 . The first step in this process involves a comprehensive analysis of the differential flatness characteristics of system (7). Following this, the linearization procedure is executed using the Taylor series expansion technique to streamline the original dynamic model. The last step entails designing an optimal feedback robust controller that takes into account model uncertainty and external disturbances.

3.1. Differential Flatness Properties of System (7)

In order to design the optimal controller rigorously in the subsequent section, we firstly need to prove that the hydraulic system defined in (7) is differentially flat. In dynamic model (7), note that the position of the driving hydraulic cylinder is selected as the flat output, that is, y ( t ) = x 3 ( t ) . We have
x ˙ 3 ( t ) = x 4 ( t ) x 4 ( t ) = y ˙ ( t ) ,
which means that x 4 ( t ) denotes a differential function. Then, we note
x ˙ 4 ( t ) = 1 m t A 1 x 2 ( t ) B p x 4 ( t ) K x 3 ( t ) F L y ¨ ( t ) = 1 m t A 1 x 2 ( t ) B p y ˙ K y ( t ) F L .
Next, from the dynamic model (7), one obtains that
x ˙ 2 ( t ) = β e V g + A 1 x 3 ( t ) D p x 1 ( t ) ( 1 x 2 ( t ) β e ) ) ( C i p + C e p ) η x 2 ( t ) A 1 x 4 ( t ) .
Note that (13) can be solved with respect to x 1 ( t ) , indicating that x 1 ( t ) is ultimately a differential function. Finally, from model (7), one has
u ( t ) = x ˙ 1 ( t ) T d + x 1 ( t ) K v .
This implies that u ( t ) is also a differential function. As a result, we complete the demonstration of the hydraulic system’s differential flatness qualities. The hydraulic system’s control loop setpoints can be also calculated using the differential flatness property.

3.2. Taylor Series Expansion

In our work, in order to efficiently design the feedback controller for the original nonlinear system (8), we first use the first-order Taylor series expansion method [35] to linearize and approximate the original nonlinear system (8) in this section. The linearization process is mainly based on the first-order Taylor expansion and the computation of the associated Jacobian matrices. The nonlinear system (8) undergoes an approximate linearization around the temporary operating point ( x * , u * ) , where x * is the current value and u * is the last sampled value. The modeling error due to the omission of higher order terms from the Taylor series expansion can also be seen as a perturbation. Therefore, the linearization process results in the following form:
x ˙ ( t ) = A x ( t ) + B u ( t ) + d ˜ ,
where d ˜ represents the modeling error due to the linearization process. A and B are Jacobian expansions that can be computed by the Taylor series expansion:
A = [ f 1 ( x ( t ) ) + f 2 ( x ( t ) ) u ( t ) ] x ( x * , u * ) = f 1 ( x ( t ) ) x ( x * , u * ) + f 2 ( x ( t ) ) u ( t ) x ( x * , u * ) = f 1 ( x ( t ) ) x ( x * , u * ) ,
B = [ f 1 ( x ( t ) ) + f 2 ( x ( t ) ) u ( t ) ] u ( x * , u * ) = f 2 ( x ( t ) ) u ( x * , u * ) ,
where f 1 ( x ( t ) ) x ( x * , u * ) is calculated by
f 1 ( x ( t ) ) x ( x * , u * ) = 1 T d 0 0 0 D p ( β e x 2 ) ( V g + A 1 x 3 ) D p x 1 ( C i p + C e p ) β e η V g + A 1 x 3 A 1 β e ( V g + A 1 x 3 ) 2 ( D p x 1 ( 1 x 2 β e ) ( C i p + C e p ) η x 2 A 1 x 4 ) A 1 β e V g + A 1 x 3 0 0 0 1 0 A m t k m t B p m t ,
and f 2 ( x ( t ) ) u ( x * , u * ) is computed as
f 2 ( x ( t ) ) u ( x * , u * ) = K v T d 0 0 0 .
As a result, using approximate linearization through first-order Taylor series expansion, one can avoid global linearization transformations and apply the optimal controller directly to the original nonlinear state-space system model.
Remark 1. 
The approximation system (15) is a linearized version of the original nonlinear system (8) around an equilibrium point. The relationship between the approximation system (15) and the original nonlinear system (8) is that the first approximation system approximates the behavior of the original nonlinear system near an equilibrium point. The first approximation system can also be used to analyze the stability of the original nonlinear system.

3.3. Optimal Robust Controller Design

After linearization around ( x * , u * ) by first-order Taylor series expansion, the original nonlinear dynamic system (8) can be written as follows:
x ˙ ( t ) = A x ( t ) + B u ( t ) + d ˜ 1 ,
where d ˜ 1 denotes the modeling error due to the linearization process of the system (8). Now, the main aim is to find a set of optimal inputs u * ( t ) that can derive the injection molding’s output x ( t ) to track the desired output x d ( t ) = [ x 1 d , , x 4 d ] within a specified period T. Suppose the difference between the optimal control input u * ( t ) and the current control input u ( t ) is Δ u ( t ) , that is,
u * ( t ) = u ( t ) + Δ u ( t ) ,
and we have
x ˙ d ( t ) = A x d ( t ) + B u * ( t ) + d ˜ 2 .
Then, the regulated system (17) can be formulated as
x ˙ ( t ) = A x ( t ) + B u ( t ) + B u * ( t ) B u * ( t ) + d ˜ 1 .
Assuming that the aggregate disturbance term in (20) is d ˜ 3 = d ˜ 1 B u * ( t ) , then one obtains
x ˙ ( t ) = A x ( t ) + B u ( t ) + B u * ( t ) + d ˜ 3 .
Denote the tracking error dynamics as E ( t ) = x ( t ) x d ( t ) ; subtracting (19) from (21) obtains
E ˙ ( t ) = A E ( t ) + B u ( t ) + d ˜ ,
where d ˜ = d ˜ 3 d ˜ 2 . According to (22), the objective performance function (10) can be formulated as
m i n u ( t ) m a x d ˜ : J ( t ) = 1 2 0 T E T ( t ) Q E ( t ) + u T ( t ) r u ( t ) ρ 2 d ˜ T d ˜ d t , r , ρ > 0 ,
where Q denotes a positive semi-definite symmetric matrix, and r and ρ denote the control signal and perturbation signal weighting factors, respectively. Therefore, the original Problem P 0 has been successfully transformed into the following Problem 1.
Problem 2. 
Given dynamic of the linearized injection hydraulic system (17) and the tracking error dynamics (22), in the case of minimum energy consumption, the aim is to find an optimal signal u ( t ) that minimizes objective function (23) assuming the largest disturbance d ˜ .
Assumption 1. 
It is assumed that (i) The energy transmitted by the disturbance signal d ˜ ( t ) is bounded, i.e., 0 d ˜ T ( t ) d ˜ ( t ) d t < , and (ii) Matrices [ A , B ] are stabilizable.
Lemma 1. 
For the linearized system (22) and the cost function (23), the optimal feedback control law is given by
u ( t ) = M E ( t ) ,
with M = 1 r B T H , where H denotes a positive definite symmetric matrix that can be solved by the following Riccati equation
A T H + HA + Q H ( 2 r B B T 1 ρ 2 N N T ) H = 0 .
Proof. 
The optimal control problem in closed form can be solved in the case of linear continuous-time systems with a quadratic cost function. For a linear system
E ˙ ( t ) = A ( t ) E ( t ) + B ( t ) u ( t ) ,
with a cost function that comprises quadratic terms in the form:
J ( t ) = 1 2 0 T { E T ( t ) Q ( t ) E ( t ) + u T ( t ) r u ( t ) } d t ,
then the optimal value of the cost function is
J 0 ( t ) = 1 2 E T ( t ) H ( t ) E ( t )
with H ( t ) being a symmetric positive definite matrix.
For the system mentioned above, the Hamiltonian function is defined as
G = 1 2 E T Q E + 1 2 u T r u + 1 2 E T H ( A E + B u ) + 1 2 ( A E + B u ) T H E
By taking the derivative of (29) with respect to the control input u and applying its extremum condition, one can obtain the following result:
G u = r u + B T H E = 0 ,
Solving the above equation for u ( t ) yields the optimal control input
u 0 ( t ) = 1 r B T ( t ) H ( t ) E ( t ) .
Substituting the optimal control input into the Hamiltonian function (29) yields the following result:
G 0 = min u [ 0 , T ] G = 1 2 E T Q E + 1 2 E T HB r 1 B T H E + 1 2 E T A T H E + 1 2 E T H A E E T HB r 1 B T H E = 1 2 E T ( Q + A T H + H A HB r 1 B T H ) E .
Next, differentiating the cost function of Equation (28) yields
J 0 t = 1 2 E T ( t ) d H ( t ) d t E ( t ) .
Equating (32) to (28) results in the following formulation of the Hamilton–Jacobi–Bellman:
1 2 E T ( t ) d H ( t ) d t E ( t ) = min u [ 0 , T ] G = 1 2 E T Q E + 1 2 E T HB r 1 B T H E + 1 2 E T A T H E + 1 2 E T H A E E T HB r 1 B T H E = 1 2 E T ( Q + A T H + H A HB r 1 B T H ) E .
After performing intermediate computations, one arrives at the following differential Riccati equation:
d H ( t ) d t = A T H + H A + Q HB r 1 B T H
The boundary condition for solving differential Equation (35) is derived from the following relation:
J 0 ( T ) = 1 2 E T ( T ) H ( T ) E ( T )
with H ( T ) = 0 . By integrating (35) from T to 0, one can find the matrix H ( t ) : t [ 0 , T ] . Additionally, using (31), one can determine the optimal control input to be applied to the system. Furthermore, using (27), one can find the minimum value of the system’s cost function.
The above procedure for solving the optimal control problem of a linear dynamical system results in the design of a Linear Quadratic Regulator. If matrices A and B in the state-space model of the system in (26) are time-invariant and the weight matrices Q and r in the quadratic cost function of Equation (27) are also time-invariant, and if the solution of the differential Riccati equation H ( t ) reaches a steady-state value (i.e., d H ( t ) d t = 0 ), then (35) becomes
A T H + H A + Q HB r 1 B T H = 0 .
In this case, the optimal feedback control is given by
u ( t ) = 1 r B T H E ( t ) ,
where H is now the solution of (37). As a result, the optimized cost function becomes J 0 = 1 2 E T H E . Note that the superposition principle of Bellman optimality states that u ( t ) = 1 r B T H E ( t ) can be used to solve for the optimal control u, which represents the minimal change in control input required to eliminate the tracking error of the state vector. Similarly, finding the worst-case disturbance by solving the optimal control problem for d ˜ yields d ˜ = 1 ρ 2 N T H E ( t ) . This completes the proof. □

4. Lyapunov Stability Analysis

This part provides a comprehensive theoretical analysis of Lyapunov stability, which proves the stability of the closed-loop system under our proposed optimal feedback controller. It is proved that the system based on the proposed controller (24) can achieve asymptotic convergence (global stability) to the desired output curve under reasonable conditions and guarantee the tracking performance of the injection velocity during IM.
Theorem 1 
([36]). The dynamic system described by x ˙ ( t ) = f ( x ( t ) ) is asymptotically stable in the vicinity of the equilibrium x 0 = 0 if there exists a function V ( x ) satisfying:
(i)
V ( x ( t ) ) to be continuous and to have a continuous first order derivative at x 0 = 0 .
(ii)
V ( x ( t ) ) > 0 , if x 0 and V ( 0 ) = 0 .
(iii)
V ˙ ( x ( t ) ) < 0 , x 0 .
Theorem 2 
([37]). Suppose that f ( t ) : [ 0 , ) R is uniformly continuous and that lim t + 0 t f ( τ ) d τ exists. Then lim t + f ( t ) = 0 holds.
As previously demonstrated, the tracking error of the injection velocity are stated as
E ˙ ( t ) = A E ( t ) + B u ( t ) + N d ˜ ,
where N = I R 4 × 4 , I is an identity matrix. The injection errors of the dynamic model’s uncertainty are represented by variable d ˜ . Now, we propose the following Lyapunov function
V ( t ) = 1 2 E T ( t ) H E ( t ) .
Differentiating (40) with respect to t, one obtains
V ˙ ( t ) = 1 2 E ˙ T H E + 1 2 E T H E ˙ = 1 2 A E + B u + N d ˜ T H E + 1 2 E T H A E + B u + N d ˜ = 1 2 E T A T + u T B T + d ˜ T N T T H E + 1 2 E T H A E + B u + N d ˜ = 1 2 E T A T H E + 1 2 u T B T H E + 1 2 d ˜ T N T H E + 1 2 E T H A E + 1 2 E T H B u + 1 2 E T H N d ˜ = 1 2 E T ( A T H + H A ) E + ( 1 2 u T B T H E + 1 2 E T H B u ) + ( 1 2 d ˜ T N T H E + 1 2 E T H N d ˜ )
Assumption 2. 
For a given positive definite matrix Q and coefficients r and ρ, there exists a positive definite matrix H that is the solution of the matrix problem shown below:
A T H + H A = Q + H ( 2 r B B T 1 ρ 2 N N T ) H ,
Note that the system is governed by the following state feedback law
u ( t ) = 1 r B T H E ( t ) .
Based on (42) and (43), (41) can be rewritten as
V ˙ ( t ) = 1 2 E T ( Q + H ( 2 r B B T 1 ρ 2 N N T ) H ) E + E T H B ( 1 r B T H E ) + ( 1 2 d ˜ T N T H E + 1 2 E T H N d ˜ ) = 1 2 E T Q E + 1 r E T H B B T H E 1 2 ρ 2 E T H N N T H E 1 r E T H B B T H E + ( 1 2 d ˜ T N T H E + 1 2 E T H N d ˜ ) = 1 2 E T Q E 1 2 ρ 2 E T H N N T H E + 1 2 d ˜ T N T H E + 1 2 E T H N d ˜ .
Lemma 2. 
The following inequality exists
1 2 ρ 2 E T H N N T H E + 1 2 d ˜ T N T H E + 1 2 E T H N d ˜ 1 2 ρ 2 d ˜ T d ˜
Proof. 
Considering the term ( ρ a 1 ρ b ) 2 , one can obtain
ρ 2 a 2 + 1 ρ 2 b 2 2 a b 0 ,
1 2 ρ 2 a 2 + 1 2 ρ 2 b 2 a b 0 ,
a b 1 2 ρ 2 b 2 1 2 ρ 2 a 2 ,
1 2 a b + 1 2 a b 1 2 ρ 2 b 2 1 2 ρ 2 a 2 .
Make the following substitutions: a = d ˜ and b = E T H N . Equation (46) becomes
1 2 d ˜ T N T H E + 1 2 E T H N d ˜ 1 2 ρ 2 E T H N N T H E 1 2 ρ 2 d ˜ T d ˜
Then, we complete the proof. □
Now, substituting (47) into (44) can obtain
V ˙ ( t ) 1 2 E T Q E + 1 2 ρ 2 d ˜ T d ˜ .
Integrating V ( t ) from 0 to T gives
0 T V ˙ ( t ) d t 1 2 0 T E Q 2 d t + 1 2 ρ 2 0 T d ˜ 2 d t 2 V ( T ) + 0 T E Q 2 d t 2 V ( 0 ) + ρ 2 0 T d ˜ 2 d t ,
Additionally, if there is a constant M p > 0 that
0 d ˜ 2 d t M p ,
then
0 T E Q 2 d t 2 V ( 0 ) + ρ 2 M p ,
that is
lim T + 0 T E Q 2 d t = 2 V ( 0 ) + ρ 2 M p .
Based on Theorem 2, we can obtain
lim t + E T ( t ) Q E ( t ) = 0 .
where Q is the positive definite matrix given above; due to H also being a positive definite matrix, the following formula can be derived
lim t + E T ( t ) H E ( t ) = 0 .
According to (54), we can obtain
lim t + V ( t ) = 0 ,
and
lim t + V ˙ ( t ) = 0 .
As a result, based on Theorem 1, it can be proven that the control scheme is globally asymptotically stable. As the number of iterations of the control algorithm increases, the tracking error tends to zero.

5. Numerical Simulations

In this section, the control strategy proposed in the previous sections is tested for its feasibility and effectiveness. In numerical experiments, all simulation experiments are performed on the platform of Matlab R2019a. The computer hardware platform is configured as Intel(R) Core i7-6820H CPU 2.7GHz, 16GB RAM, and 64-bit Windows 10 operating system. Table 1 shows some main characteristic parameters of the IMM model.
The initial experiment involves evaluating the tracking control of the injection velocity and other states within the IMM utilizing the optimal controller we designed. The corresponding numerical experimental outcomes are presented in Figure 3 and Figure 4. These illustrate that the optimal control scheme we developed demonstrates excellent performance in the speed-tracking control of the IMM. Note that the feedback control gain is determined at each control strategy iteration by repeatedly solving the algebraic Riccati Equation (25). Figure 3 depicts that the control input converges quickly and remains stable. Similarly, Figure 4 shows that all state variables in the IMM swiftly and precisely converge to the reference set points, enabling the desired curves to be tracked. In addition, another injection velocity tracking signal within the IMM is also evaluated, and the corresponding numerical results are presented in Figure 5 and Figure 6. Similarly, the results show that the control input converges rapidly and remains stable, and all states in the IMM converge swiftly and accurately to the preset references, further validating the effectiveness of the controller we designed.
Furthermore, during our numerical process, we also examine the developed nonlinear control method’s tracking accuracy at various reference set points. Table 2 and Table 3 show the achieved results, respectively. The control method applies a change equal to Δ α % to m t in (8) for the IMM model and the tracking performance under disturbances, as shown in Table 4, respectively. It can be noticed in Table 2 that the tracking accuracy is satisfactory despite the model perturbations that exist. As a result, the control approach suggested in this study performs satisfactorily and can precisely track the reference set points even when there are disturbances.
It is important to note that the selection of parameters r, Q, and ρ in (25) determines the transient performance of the control algorithm. The parameter r represents the control output weight factor in the objective function. A small value of r can eliminate tracking errors in the states. The parameter Q represents the state variable weight matrix in the objective function. A large value of Q allows state variables to converge quickly to their reference set points. The choice of decay coefficient ρ affects the resilience of the control strategy. By selecting a value of ρ small enough to satisfy the tracking performance criterion in (25), the global asymptotic stability of the control loop is ensured and tracking errors in state variables are eliminated.

6. Conclusions

In this paper, we propose an efficient optimal robust control strategy for the dynamic model of a variable speed pump-controlled hydraulic cylinder’s speed servo system within an IMM. The dynamic system model is augmented and approximately linearized using the first-order Taylor series. A reliable feedback controller is then designed. At each time step of the control scheme, the corresponding Riccati equation is solved iteratively to determine the reliable feedback gain of the feedback controller. The global asymptotic stability of the control strategy is demonstrated through Lyapunov analysis. Numerical experiments are conducted to validate the feasibility of the control strategy. The proposed approach can be easily adapted to other similar types of nonlinear optimal control problems that arise in discrete manufacturing processes.

Author Contributions

Methodology, G.W.; Validation, G.W. and J.L.; Formal analysis, G.W.; Investigation, J.L.; Resources, Z.W.; Writing—original draft, G.W.; Writing—review & editing, Z.R.; Supervision, Z.R.; Project administration, Z.R. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Grants of Key-Area Research and Development Program of Guangdong Province (2021B0101200005), the National Natural Science Foundation of China (62073088, U1911401), and the Open Research Fund from Guangdong Laboratory of Artificial Intelligence and Digital Economy (SZ).

Data Availability Statement

All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Zheng, R.; Tanner, R.I.; Fan, X.-J. Injection Molding: Integration of Theory and Modeling Methods; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  2. Chaoyan, L. Progress of World Plastics Industry (I), 2020; General Plastics: Nanjing, China, 2021; Volume 50, pp. 1–16. [Google Scholar]
  3. Fernandes, C.; Pontes, A.J.; Viana, J.C.; Gaspar-Cunha, A. Modeling and optimization of the injection-molding process: A review. Adv. Polym. Technol. 2018, 37, 2018. [Google Scholar] [CrossRef]
  4. Gao, H.; Zhang, Y.; Zhou, X.; Li, D. Intelligent methods for the process parameter determination of plastic injection molding. Front. Mech. Eng. 2018, 13, 85–95. [Google Scholar] [CrossRef]
  5. Fu, H.; Xu, H.; Liu, Y.; Yang, Z.; Kormakov, S.; Wu, D.; Sun, J. Overview of injection molding technology for processing polymers and their composites. ES Mater. Manuf. 2020, 8, 3–23. [Google Scholar] [CrossRef]
  6. Huang, S.; Tan, K.K.; Lee, T.H. Predictive control of ram velocity in injection molding. Polym.-Plast. Technol. Eng. 1999, 38, 285–303. [Google Scholar] [CrossRef]
  7. Dubay, R.; Pramujati, B.; Han, J.; Strohmaier, F. An investigation on the application of predictive control for controlling screw position and velocity on an injection molding machine. Polym. Eng. Sci. 2007, 47, 390–399. [Google Scholar] [CrossRef]
  8. Zhao, P.; Zhang, J.; Dong, Z.; Huang, J.; Zhou, H.; Fu, J.; Turng, L.-S. Intelligent injection molding on sensing, optimization, and control. Adv. Polym. Technol. 2020, 2020, 1–22. [Google Scholar] [CrossRef] [Green Version]
  9. Tan, K.K.; Huang, S.; Jiang, X. Adaptive control of ram velocity for the injection moulding machine. IEEE Trans. Control. Syst. Technol. 2001, 9, 663–671. [Google Scholar] [CrossRef]
  10. Reiter, M.; Stemmler, S.; Hopmann, C.; Reßmann, A.; Abel, D. Model predictive control of cavity pressure in an injection moulding process. IFAC Proc. Vol. 2014, 47, 4358–4363. [Google Scholar] [CrossRef] [Green Version]
  11. Wu, S.; Jin, Q.; Zhang, R.; Zhang, J.; Gao, F. Improved design of constrained model predictive tracking control for batch processes against unknown uncertainties. ISA Trans. 2017, 69, 273–280. [Google Scholar] [CrossRef]
  12. Zou, T.; Wu, S.; Zhang, R. Improved state space model predictive fault-tolerant control for injection molding batch processes with partial actuator faults using ga optimization. ISA Trans. 2018, 73, 147–153. [Google Scholar] [CrossRef]
  13. Froehlich, C.; Kemmetmüller, W.; Kugi, A. Control-oriented modeling of servo-pump driven injection molding machines in the filling and packing phase. Math. Comput. Model. Dyn. Syst. 2018, 24, 451–474. [Google Scholar] [CrossRef] [Green Version]
  14. Han, C.; Jia, L.; Peng, D. Model predictive control of batch processes based on two-dimensional integration frame. Nonlinear Anal. Hybrid Syst. 2018, 28, 75–86. [Google Scholar] [CrossRef]
  15. Froehlich, C.; Kemmetmüller, W.; Kugi, A. Model-predictive control of servo-pump driven injection molding machines. IEEE Trans. Control Syst. Technol. 2019, 28, 1665–1680. [Google Scholar] [CrossRef]
  16. Wu, S.; Zhang, R. A two-dimensional design of model predictive control for batch processes with two-dimensional (2D) dynamics using extended non-minimal state space structure. J. Process. Control 2019, 81, 172–189. [Google Scholar] [CrossRef]
  17. Khosravani, M.R.; Nasiri, S. Injection molding manufacturing process: Review of case-based reasoning applications. J. Intell. Manuf. 2020, 31, 847–864. [Google Scholar] [CrossRef]
  18. Liu, J.; Guo, F.; Gao, H.; Li, M.; Zhang, Y.; Zhou, H. Defect detection of injection molding products on small datasets using transfer learning. J. Manuf. Process. 2021, 70, 400–413. [Google Scholar] [CrossRef]
  19. Ren, Z.; Li, Y.; Wu, Z.; Xie, S. Deep learning-based predictive control of injection velocity in injection molding machines. Adv. Polym. Technol. 2022, 2022, 7662264. [Google Scholar] [CrossRef]
  20. Cho, Y.; Cho, H.; Lee, C.-O. Optimal open-loop control of the mould filling process for injection moulding machines. Optim. Control. Appl. Methods 1983, 4, 1–12. [Google Scholar] [CrossRef]
  21. Havlicsek, H.; Alleyne, A. Nonlinear control of an electrohydraulic injection molding machine via iterative adaptive learning. IEEE/ASME Trans. Mechatronics 1999, 4, 312–323. [Google Scholar] [CrossRef]
  22. Dubay, R. Self-optimizing MPC of melt temperature in injection moulding. ISA Trans. 2002, 41, 81–94. [Google Scholar] [CrossRef]
  23. Yao, K.; Gao, F. Optimal start-up control of injection molding barrel temperature. Polym. Eng. Sci. 2007, 47, 254–261. [Google Scholar] [CrossRef]
  24. Chen, W.-C.; Tai, P.-H.; Wang, M.-W.; Deng, W.-J.; Chen, C.-T. A neural network-based approach for dynamic quality prediction in a plastic injection molding process. Expert Syst. Appl. 2008, 35, 843–849. [Google Scholar] [CrossRef]
  25. Xia, W.; Luo, B.; Liao, X. An enhanced optimization approach based on gaussian process surrogate model for process control in injection molding. Int. J. Adv. Manuf. Technol. 2011, 56, 929C942. [Google Scholar] [CrossRef]
  26. Farahani, S.; Khade, V.; Basu, S.; Pilla, S. A data-driven predictive maintenance framework for injection molding process. J. Manuf. Process. 2022, 80, 887–897. [Google Scholar] [CrossRef]
  27. Xiao, H.; Meng, Q.-X.; Lai, X.-Z.; Yan, Z.; Zhao, S.-Y.; Wu, M. Design and trajectory tracking control of a novel pneumatic bellows actuator. Nonlinear Dyn. 2022, 111, 1–18. [Google Scholar] [CrossRef]
  28. Ren, Z.; Wu, G.; Wu, Z.; Xie, S. Hybrid dynamic optimal tracking control of hydraulic cylinder speed in injection molding industry process. J. Ind. Manag. Optim. 2023, 19, 5209–5229. [Google Scholar] [CrossRef]
  29. Ruan, Y.; Gao, H.; Li, D. Improving the consistency of injection molding products by intelligent temperature compensation control. Adv. Polym. Technol. 2019, 2019, 1591204. [Google Scholar] [CrossRef]
  30. Stemmler, S.; Vukovic, M.; Ay, M.; Heinisch, J.; Lockner, Y.; Abel, D.; Hopmann, C. Quality control in injection molding based on norm-optimal iterative learning cavity pressure control. IFAC-PapersOnLine 2020, 53, 10380–10387. [Google Scholar] [CrossRef]
  31. Guo, F.; Zhou, X.; Liu, J.; Zhang, Y.; Li, D.; Zhou, H. A reinforcement learning decision model for online process parameters optimization from offline data in injection molding. Appl. Soft Comput. 2019, 85, 105828. [Google Scholar] [CrossRef]
  32. Tian, Y.-C.; Gao, F. Injection velocity control of thermoplastic injection molding via a double controller scheme. Ind. Eng. Chem. Res. 1999, 38, 3396–3406. [Google Scholar] [CrossRef] [Green Version]
  33. Yang, A.; Guo, W.; Han, T.; Zhao, C.; Zhou, H.; Cai, J. Feedback control of injection rate of the injection molding machine based on pid improved by unsaturated integral. Shock Vib. 2021, 2021, 1–9. [Google Scholar] [CrossRef]
  34. Wang, J.; Mao, Q.; Jiang, N.; Chen, J. Effects of injection molding parameters on properties of insert-injection molded polypropylene single-polymer composites. Polymers 2021, 14, 23. [Google Scholar] [CrossRef] [PubMed]
  35. Bartle, R.G.; Sherbert, D.R. Introduction to Real Analysis; Wiley: New York, NY, USA, 2000. [Google Scholar]
  36. Nguyen, N.T. Lyapunov stability theory. In Model-Reference Adaptive Control: A Primer; Springer: Berlin/Heidelberg, Germany, 2018; pp. 47–81. [Google Scholar]
  37. Farkas, B.; Wegner, S.-A. Variations on barbălat’s lemma. Am. Math. Mon. 2016, 123, 825–830. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Structure diagram of a typical IMM.
Figure 1. Structure diagram of a typical IMM.
Mathematics 11 02619 g001
Figure 2. Structure diagram of hydraulic system in an injection molding machine.
Figure 2. Structure diagram of hydraulic system in an injection molding machine.
Mathematics 11 02619 g002
Figure 3. Optimal control input u ( t ) at tracking desired speed v d 1 .
Figure 3. Optimal control input u ( t ) at tracking desired speed v d 1 .
Mathematics 11 02619 g003
Figure 4. State variables x ( t ) at tracking desired speed v d 1 .
Figure 4. State variables x ( t ) at tracking desired speed v d 1 .
Mathematics 11 02619 g004
Figure 5. Optimal control input u ( t ) at tracking desired speed v d 2 .
Figure 5. Optimal control input u ( t ) at tracking desired speed v d 2 .
Mathematics 11 02619 g005
Figure 6. State variables x ( t ) at tracking desired speed v d 2 .
Figure 6. State variables x ( t ) at tracking desired speed v d 2 .
Mathematics 11 02619 g006
Table 1. Main parameters in the dynamic system.
Table 1. Main parameters in the dynamic system.
VariableVariable SymbolValue
Speed gain K v 100
Time constant T d 0.005
Oil viscosity η 4.6 × 10 3
Elastic modulus of oil volume β e 1.6 × 10 9
Total volume of high pressure chamber V g 3.9 × 10 3
Effective area of rodless chamber A 1 3.1 × 10 3
Displacement of quantitative pump D p 6.0 × 10 6
Leakage coefficient of quantitative pump C i p 4.5 × 10 11
Table 2. Tracking the RMSE of the desired speed v d 1 in the disturbance-free case.
Table 2. Tracking the RMSE of the desired speed v d 1 in the disturbance-free case.
RMSE x 1 RMSE x 2 RMSE x 3 RMSE x 4
s e t p o i n t 1 0.0236 102.1911 8.6393 × 10 5 1.0533 × 10 4
s e t p o i n t 2 0.0847 147.2880 2.7525 × 10 4 1.5215 × 10 4
s e t p o i n t 3 0.0991 180.0929 5.2280 × 10 4 1.8598 × 10 4
s e t p o i n t 4 0.3138 417.7483 9.5517 × 10 4 4.3165 × 10 4
s e t p o i n t 5 0.2297 464.5643 1.5000 × 10 3 4.7959 × 10 4
s e t p o i n t 6 0.2979 509.3806 2.2000 × 10 3 5.2590 × 10 4
Table 3. Tracking the RMSE of the desired speed v d 2 in the disturbance-free case.
Table 3. Tracking the RMSE of the desired speed v d 2 in the disturbance-free case.
RMSE x 1 RMSE x 2 RMSE x 3 RMSE x 4
s e t p o i n t 1 0.0236 102.1911 8.6393 × 10 5 1.0533 × 10 4
s e t p o i n t 2 0.0833 147.2576 2.7526 × 10 4 1.5211 × 10 4
s e t p o i n t 3 0.3484 541.8267 7.4757 × 10 4 5.5983 × 10 4
s e t p o i n t 4 0.3458 625.1940 1.5000 × 10 4 6.4563 × 10 4
s e t p o i n t 5 0.5635 471.6619 2.3000 × 10 4 4.8771 × 10 4
s e t p o i n t 6 0.2884 508.2908 3.1000 × 10 4 5.2452 × 10 4
Table 4. Tracking the RMSE of the desired speed v d 1 in the case of disturbances.
Table 4. Tracking the RMSE of the desired speed v d 1 in the case of disturbances.
Δ α % RMSE x 1 RMSE x 2 RMSE x 3 RMSE x 4
0 % 0.2160 464.2767 0.0015 4.5740 × 10 4
10 % 0.2188 464.8364 0.0022 4.6235 × 10 4
20 % 0.2215 465.3996 0.0029 4.6731 × 10 4
30 % 0.2269 466.5362 0.0035 4.7725 × 10 4
40 % 0.2296 467.1095 0.0042 4.8222 × 10 4
50 % 0.2350 468.2659 0.0049 4.9217 × 10 4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, G.; Ren, Z.; Li, J.; Wu, Z. Optimal Robust Tracking Control of Injection Velocity in an Injection Molding Machine. Mathematics 2023, 11, 2619. https://doi.org/10.3390/math11122619

AMA Style

Wu G, Ren Z, Li J, Wu Z. Optimal Robust Tracking Control of Injection Velocity in an Injection Molding Machine. Mathematics. 2023; 11(12):2619. https://doi.org/10.3390/math11122619

Chicago/Turabian Style

Wu, Guoshen, Zhigang Ren, Jiajun Li, and Zongze Wu. 2023. "Optimal Robust Tracking Control of Injection Velocity in an Injection Molding Machine" Mathematics 11, no. 12: 2619. https://doi.org/10.3390/math11122619

APA Style

Wu, G., Ren, Z., Li, J., & Wu, Z. (2023). Optimal Robust Tracking Control of Injection Velocity in an Injection Molding Machine. Mathematics, 11(12), 2619. https://doi.org/10.3390/math11122619

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop