Next Article in Journal
A Systems Thinking Archetype to Understand, Analyze, and Evaluate the Evolution of International Political Crises
Previous Article in Journal
Using Systems Thinking and Modelling: Ecological Land Utilisation Efficiency in the Yangtze Delta in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementation and Validation of a Simple Direct Identification Method for an Experimental Multi-Span Web Transport System

by
Nicola Ivan Giannoccaro
1,* and
Tetsuzo Sakamoto
2
1
Department of Innovation Engineering, University of Salento, 73100 Lecce, Italy 2 Department of Control Engineering, Kyushu Institute of Technology, Kitakyushu 804-8550, Japan
2
Department of Control Engineering, Kyushu Institute of Technology, Kitakyushu 804-8550, Japan
*
Author to whom correspondence should be addressed.
Systems 2022, 10(1), 17; https://doi.org/10.3390/systems10010017
Submission received: 26 December 2021 / Revised: 5 February 2022 / Accepted: 11 February 2022 / Published: 13 February 2022
(This article belongs to the Special Issue Industrial Applications of System Identification)

Abstract

:
The industrial processes that require the use of the web require a control system which allows for preserving the properties of the web unaltered, avoiding the risk of wrinkling, tearing, breakage and other defects. This control generally takes place by detecting the tension and the speed in certain points of the system since these variables determine the stress state on the web, which, if altered beyond certain ranges, can lead to the defects mentioned above. The problem of tension and web speed control is very demanding because the system’s dynamic is a function of many process variables that often vary over a wide range. In this study, an experimental system consisting of 12 rollers, four motorised, was analysed. This system was divided into four subsystems according to the logic of decentralised control. The tension of the initial and final subsystems and the speeds of the two central subsystems were monitored. This study proposes estimating continuous-time transfer functions using experimental time-domain data. A nonlinear least-squares search-based method minimises a weighted prediction error norm for directly identifying the mathematical model used to describe the web transport system. To test the performance of the proposed strategy, experimental data were collected in an open-loop configuration with constant voltage given to the four servo motors. The collected data were subsequently processed to define an extremely simple system model composed of a very limited number of parameters representing the system through transfer functions. The model was further validated by comparing the results obtained through simulations with the experimental data obtained with different inputs, and was also validated with some closed-loop tests.

1. Introduction

The term web is used to describe thin materials manufactured and processed in a continuous and flexible strip. Web materials cover a wide spectrum, from extremely thin plastics to paper, fabrics, metals and composite materials. Today many industrial processes employ webs, stored in the form of reels and moved by rollers, to mass produce a wide variety of products made from materials that are characterized by a continuous strip. The possibility of modelling and controlling web handling systems has been studied for a long time (i.e., [1,2,3,4,5,6,7,8]) to limit possible web damage during the transport system’s continuous operation. Accurate modelling may be useful for designing the control system or estimating some parameters by comparing them with the experimental tests. For building a dynamical model, lumped parameter expressions may designate a web section between two adjacent drive rolls. Nevertheless, there is also the necessity of introducing the property of viscoelasticity to the web [1] by using an opportune expression. A multi-span web tension system’s physical lumped model must mainly consider three equations: the conservation mass, torque balance and viscoelasticity (Voigt–Kelvin approach) [9].
The problem of tension and web speed control for systems using the web is very demanding and by no means trivial because the system’s dynamic is a function of many process variables that often vary over a wide range. Traditional linear control systems cannot adequately address this problem, even using robust methods. To solve the problem, decentralised control strategies were formulated which break the system down into subsystems and consider only the interactions between the nearest subsystems ([2,5,7]). Through these particular strategies, a MIMO (multiple-input, multiple-output) system can be traced back to a set of SISO subsystems (single-input single-output) controlled separately with PI controllers turned to use other algorithms to evaluate mutual interactions.
Optimal control of a system by means of PI(D) controllers implies optimising each controller’s proportional and integral parameters (and derivative if required). A very effective solution for this optimisation can be obtained by simulating the process through the estimated system model [10]. In this sense, the dynamic modelling of the system is ineffective since the calibration of this model involves numerous parameters linked to the viscoelastic nature of the material under examination, which are not constant but are a function of the state of stress, temperature, deformation speed and of other system and process variables. Recent studies propose specific robust control strategies to deal with control problems of unmodeled components of the web transport system. In [11], a robust controller based on dynamic surface control (DSC) is proposed to ensure the stability of the controlled system in the face of a mathematical model of the web transport system that is affected by bounded uncertainties. In [12], a robust decentralised control scheme is proposed to divide the control input for each subsystem into two parts to compensate for the dynamic model parameters variation and error. In [13], a decentralised guaranteed constant H∞ control strategy is proposed for large-scale web-winding systems with uncertain and time-varying parameters. In [14], an observer-based decentralised output feedback control method is presented for possible application to interconnected discrete systems, such as the multi-spans web transport platform.
Given the extreme difficulty of getting to a dynamic model of origin, in industrial practice, alternative approaches to identify the model, called Identification for Control (I4C), are used, reaching extremely simple models [15]. These models do not necessarily faithfully represent the system in question, but they do it in a limited way according to the purpose of the model, i.e., to control the system [16,17]. These approaches consider models consisting of an extremely limited number of parameters (generally two or three) and evaluate the validity of a model by comparing, for example, the response of the model and the real system to a step input. These modelling methodologies use experimental data to directly identify a model with few parameters and are named data-driven techniques for identification. Several studies can be found on the development and the applications of these techniques. In [18], the data-driven technique controls smart power generation systems. In [19], the data-driven technique is used to analyse the fault identifiability performance of a dynamic system.
From the previous considerations, considering the development of robust control techniques applicable to dynamic models with uncertain parameters and the possibility of modelling in a simplified way by identifying the model from the experimental behaviour of the system, the present study dealt with experimentally testing the validity of models simplified for a multi-roller web transport platform.
In this work, an approach of this type was applied to identify the model of an experimental system consisting of twelve rollers (four motorised). This system was divided into four subsystems in the logic of decentralised control. The tension of the initial and final subsystems and the speeds of the two central subsystems were monitored. Therefore, this system is a MIMO system, which is traced back to four subsystems. The model of the system was estimated in this work using the data of a single test in an open-loop that allows for simulating the system’s behaviour with step inputs in a fairly wide range of operations. For each subsystem, the algorithm for estimating a continuous-time transfer function starting from time-domain data uses nonlinear least-squares search-based updates to minimise a weighted prediction error norm [20,21].
The estimated model was also used to simulate the closed-loop system’s operation. The developed application allows for changing the reference set-points to impose and modify the controller parameters, i.e., the gains of the PI regulators and the smoothing constant of the reference signal. Therefore, once the model’s validity in a certain field of operation is verified, the application becomes a powerful tool to carry out an unlimited number of tests, acting on different degrees of freedom and observing the system response variation. The proposed work’s main novelty consists in having tested the applicability of data-driven identification techniques of the dynamic model to an experimental platform with multiple rolls and different connected sections. The analysis was carried out with an extensive experimental campaign by moving the web on the platform with different speeds and forces in open-loop and closed-loop control modes. The positive results obtained are particularly interesting considering the complexity of a classic physical model of the considered system and the related uncertainties and parameters that are difficult to measure, as shown. The data-driven model obtained can allow for managing the experimental system with good accuracy. This paper describes all the theoretical and experimental details and the results obtained.

2. The Experimental Web Transport System and Its Model

The realised system, developed in the last years and introduced in [22,23], comprises four main sections, each driven by one servomotor. The sections are strongly interlaced, and a total of 12 rolls, placed on a mechanical frame at different heights, make it up. The realised platform is intended to represent a large transport system similar to many industrial ones. The system was completely renewed at the end of 2015, with its performance improved by substituting all the rolls and their bearings with new ones (low weight and low friction). The platform is depicted in Figure 1, and a scheme is shown in Figure 2. The four servomotors that drive the transport system have 750 (W) of power 2.39 [Nm] of maximum torque. They are set in torque control mode by using an internal control of the servomotor to guarantee a precise force delivered to the system. The first servomotor is dedicated to the unwinder section, the second to the lead section, the third to the draw-roll section, the last to the winder section, as shown in Figure 2.
Two pairs of tension sensors (one for each side of the web) were placed after the unwinder roll (first pair of tension sensors) and before the winder roll (second pair of tension sensors). The first pair of tension sensors was followed by a pair of guide rolls (lead section) wrapped to maximize the contact area between the web and the drive rolls, reducing the possibility of slippage between the web and the rolls during the web transport.
The servomotors were driven by voltage input signals of Ui, i = 1, …, 4, which, in torque control mode, are proportional to the torques that the motors exert on the corresponding roller. The tension sensor signals related to the unwinder roll and the winder roll were acquired by the A/D board, considering the average value of the two corresponding sides for measuring the tension after the unwinder (T1) and before the winder (T4). The four motor encoder signals fed a digital counter: in the proposed identification strategy, the speed of the lead section (v2) and the draw-roll section (v3) are the control variables. The controller’s CPU receives signals through A/D boards and counters, performs the control algorithm and generates the command signals in real-time, driving the servomotors via D/A boards with a sampling time of 0.01 [s]. Moreover, a set of filtering circuits was inserted; they were hardware filter banks (model FN270D) connected to the servo packs driving the four motors. Moreover, a digital moving average filter was applied to the acquired data.
A functional diagram of the transport line indicating the four servomotors (type SGMAS 08ACA21), their driving servo packs, CPU controller, analogue (A/D Board PCI 3180) and digital (Digital Counter PCI 632206 TK) acquisition boards, sensors, filters banks (type FN207D), their connections (electrical, digital and analogue signals) and power supply is shown in Figure 3.
As depicted in Figure 1, the whole system was mounted on a mechanical frame designed to support the system’s components. Each of the four motors was connected to the respective roller using specific low-friction connection joints. The power cables were placed far away concerning the signal cables, and all the cables were covered with insulating tape. The system comprises four drive rolls divided into three sections with lengths named L1, L2, and L3. The [8] geometric and mass properties of the four drive rolls are shown in Table 1.
Table 2 sums up the main material properties of the Bi-Oriented Polypropylene (OPP) film used in the experimental web handling system, and Figure 4 shows a stock image of the film roll (a) and the rolls mounting on the experimental system (b).

2.1. Lumped Model for the Experimental Web Tension Control System

The model of the web transport systems is based on three laws applied at each section between two consecutive rolls [1]:
-
Conservation of mass: the law of conservation of mass for the web section for evaluating the relationship between the speeds of two adjacent rolls and the strain in the web;
-
Torque balance: the angular velocity of the kth roll can be obtained through a torque balance as a function of the tension forces Tk+1 and Tk applied to the roll from the web;
-
Voigt model: for taking into account the linear viscoelasticity of web-material expressed as a combination of linear springs and dampers, the Voigt model was considered.
The torque balance equation is expressed in Equation (1) for the generic kth section. Jk and rk are respectively the total inertia and the roll’s radius of the kth roll, uk is the motor torque applied at the kth roll, Ck and kfk are respectively the dry friction torque and the viscous friction coefficient. The conservation of mass equation is expressed in Equation (2) for the kth and kth+1 section with longitudinal velocities vk and vk+1 where ε is elasticity coefficient, and Lk is the length between the kth and the kth+1 roll. The Voigt model equation is shown in Equation (3), where σk is the stress applied at the kth section, E is the material Young’s modulus, η is the material viscosity coefficient.
d J k · ω k d t = r k · T k + 1 T k + 1 + u k C k k f k · ω k
ε = 1 L k · ( v k + 1 v k ) d t
σ k = E · ε + η · d ε d t
Equations (1)–(3) may be applied for each of the four sections of the system and then transformed in the Laplace domain considering starting conditions equal to 0 for tension force Tk and velocity ωk of each kth section. Substituting Equation (2) in Equation (3) in the Laplace domain, Equation (3) may be transformed in a relation between Tk(s) and ε(s) as expressed in Equation (4), where A is the web section area. The last equality of Equation (4) introduces a polynomial, named P(s), depending only on the physical characteristics of the web.
T k s = A · η · 1 + η E · s η E · s · s · ε s = P s · s · ε s
Equation (1), transformed in the Laplace domain, permits the calculation of the longitudinal velocities vk(s) for the kth section; Equations (5)–(8) show the expression of the longitudinal velocities for the four sections of the considered platform, considering the orientation of the motor for determining the sign of the motor torque uk applied to the kth section.
v 1 s = r 1 s · J 1 · T 1 s · r 1 u 1 s C 1 k f 1 · ω 1 s   ,
v 2 s = r 2 s · J 2 · ( T 2 s T 1 s ) · r 2 + u 2 s C 2 k f 2 · ω 2 s   ,
v 3 s = r 3 s · J 3 · ( T 4 s T 2 s ) · r 3 + u 3 s C 3 k f 3 · ω 3 s   ,
v 4 s = r 4 s · J 4 · u 4 s T 4 s · r 4 C 4 k f 4 · ω 4 s   ,
The inputs of the system are the motor torque values u1, u2, u3, u4, the outputs of the system are the forces T1 and T4 and the longitudinal roll speed values are v2 and v3 of Section 2 and Section 3. Equations (5)–(8) may be integrated with the Voigt–Kelvin viscoelasticity law (Equation (5) applied at each section) allowing for the determination of the values of T1, T2 and T4. It is not difficult to show how complex the system is, as interconnections exist between the input and outputs of the four subsystems [6,22].
To have an idea of the structure of the mathematical physical model in a manner that simplifies Equations (5)–(8) by neglecting the frictional effects expressed by the coefficients Ck and kfk, which are very difficult to estimate experimentally, the resulting model without frictional terms is shown in Figure 5a in the Laplace domain. The model structure shown in Figure 5a permits highlighting the inputs (u1, u2, u3, u4), the outputs (the forces T1 and T4 and the longitudinal roll speed v2 and v3) and the four subsystems corresponding to the system sections.
This type of mathematical description does not seem to be effective for the following reasons:
  • There is no way to appreciate the interactions that affect the individual subsystems. This is a MIMO system, and each input influences individual outputs, as is clearly depicted in Figure 5a.
  • The mathematical model used to describe the system’s behaviour cannot provide any indication of the value of the dissipative parameters; if not appropriately identified, viscous friction coefficients make the mathematical model wrong and unable to describe the reality of the physical phenomenon.
  • It is impossible to study the system’s stability in any way and to predict its behaviour by varying the operating conditions. Most control techniques require formulating the problem as a “transfer function” to control and operate on the poles and zeros that characterise the system considered.
For all these reasons, a new approach was proposed to analyse the problem considering exclusively a direct identification of the dynamical MIMO model shown in Figure 5b based on the experimental data obtained evaluating the inputs and outputs of the system.

2.2. Development of a Procedure for Model Identification

The input–output relation of the system model shown in Figure 5b is expressed with input u1, u2, u3, u4 and output T1, v2, v3, T4 and may be expressed in a generic form with 16 transfer functions Gij as expressed in Equation (9):
( T 1 v 2 v 3 T 4 ) = ( G 11 G 12 G 13 G 14 G 21 G 22 G 23 G 24 G 31 G 32 G 33 G 34 G 41 G 42 G 43 G 44 ) · ( u 1 u 2 u 3 u 4 )
The procedure here proposed for model identification considers that for a system such as the one analysed, the MIMO model Equation (9) includes transfer functions that, for the physical constitution of the system, have a different weight and importance for the system behaviour. In particular, the extra-diagonal transfer functions referring to the effects of the external subsystems on the ith subsystem may be considered an uncertainty of the estimation procedure, a disturbance leading to a model constituted by 4 SISO systems, which are characterised by the estimated SISO transfer functions Giie. The estimated identified simplified model is expressed by Equation (10).
( T 1 v 2 v 3 T 4 ) = ( G 11 e 0 0 0 0 G 22 e 0 0 0 0 G 33 e 0 0 0 0 G 44 e ) · ( u 1 u 2 u 3 u 4 )

2.3. Collecting Experimental Data for System Identification

Model identification requires a preliminary collection of input and output data from experimental tests. In this phase, it is important to demonstrate, through statistical inference methods, that the system outputs show a certain repeatable behaviour. When these conditions are satisfied, it is possible to accept such inputs to the bench machine.
When the web tension platform works in an open loop, choosing appropriate servomotor voltage input values becomes important for the experimental campaign. Each motor of the system is set up for torque control modality, so that it is possible to set the values of voltage input signals (in Volts) for each servo motor. This information is converted in torque through the motor constant Kk = 0.8 [Nm/V], so that torque [Nm] = tension signal [V] x motor constant [Nm/V]. The block diagram of the open-loop tests is schematised in Figure 6.
The following sets of voltage inputs (Table 3) were assigned respectively to unwinder, lead-section, draw-roll and winder servomotors (u1, u2, u3, u4). For each combination of inputs, five tests were carried out to have detailed information about the statistical repeatability of the tests.
The five repetitions of every test show a similar trend: in Figure 7, the outputs (T1, v2, v3, T4) of five repetitions of Test E (Table 3). The temporal length of all the tests was 15 s.
The proposed identification procedure is based on a complete routine developed in the Matlab environment (using the System Identification Toolbox app) [21] with an integrated graphical user interface that permits the user first to import the data of the experimental tests to be analysed. All the different repetitions of the same test are used to identify and validate the model, making the validation more reliable. The following step consists of estimating the transfer function models; it is necessary to define the number of poles, the number of zeros and, eventually, the delay. Using the System Identification Toolbox app [21] makes it possible to compare with the same experimental data different types of transfer functions, having different combinations of zeros, poles and delay. In Figure 8, the comparison of three different estimated transfer functions (tf1, tf2 and tf3) with the experimental data referred to T1 in one of the five repetitions for Test E. The transfer functions shown in Figure 8 have respectively 2 poles, 1 zero, 0.1 delay (tf1), 2 poles, 0 zeros, 0 delay (tf2) and 2 poles, 0 zeros and 0.1 delay for tf3.

3. Results

In the industrial application of system identification, an important objective is often to obtain the least complex model within the limits of the required model accuracy. In our case, after a few attempts, an optimal number of parameters was determined for web platform system identification: 2 poles, 1 zero and 0 delays for all the transfer functions to estimate G 11 e , G 22 e , G 33 e , G 44 e The realised routine, after loading the experimental data of the considered test, permits the insertion of the characteristics of the transfer functions to estimate the number of poles, number of zeros and delay, and then estimates the transfer functions, explaining the estimated poles, zeros and delays and, finally, produces a comparative plot between the experimental data related to the variables T1, v2, v3 and T4 and the same variables obtained by the estimated transfer functions with the same input of the experimental data, also giving an index of performance in the percentage of the similarity of behaviours. Figure 9 shows the identification results for the first repetition of Test E (in all the Figures of the manuscript, T1 and T4 are expressed in Newton and v2 and v3 in m/s), and the poles and zeros of the estimated transfer functions are reported in Table 4.
The poles and zeros are properties of the transfer function and, therefore, the differential equation describing the input–output system dynamics. They characterise the differential equation and provide a complete description of the system, including interesting information about the system stability.
Moreover, the estimated model was validated using a different repetition of a test in the same condition used for identification. Figure 10 shows the validation results obtained comparing the simulated response from Repetition 3 of Test E (named E-OL-R3) and the experimental response from Test E Repetition 1 (OL-E-R1). It is evident, also considering the performance parameter, that the transfer functions estimated from E-OL-R3 overlap very well with the experimental data of another test (OL-E-R1), validating the proposed procedure with different repetitions of the same combination of inputs.

System Simulation Using the Measured Input

The purpose of the proposed procedure is to verify that the estimated model fulfils the modelling requirements according to subjective and objective criteria of good model approximation. The model of the system estimated can be used to simulate the system’s behaviour with a certain combination of torque to the motors, different from those used for the estimation. The outputs simulated can be compared with the measured ones to validate the model in different conditions. If this comparison gives good results, then it is possible to predict the system’s behaviour, during both the transient and the steady state, at a work point different from the one of identification. The validation procedure is explained in the diagram shown in Figure 11, where a generic open-loop experimental test (i.e., Test B Repetition 1 as indicated in Figure 11, referring to the terminology in Table 3) is loaded. The measured outputs (the variables T1, v2, v3, T4) are compared with the simulated outputs obtained by the model identified from a different test (i.e., Test E) and then analysed with the same input used for the experimental data. If this comparison gives good results, it is possible to predict the system’s behaviour, during both the transient and the steady state, at a work point different from the one of identification.
The comparison was carried out for all the experimental tests indicated in Table 3; for each combination of input, the average of the five repetitions was considered, and the identified model was obtained from Test E (Repetition 3). The results are depicted in Figure 12, Figure 13 and Figure 14, where the validations with the tests A, B, C, D, F and G are shown for each control variable of each test.
The results illustrated in Figure 12, Figure 13 and Figure 14 clearly show that the model identified by a single experimental test with a specific combination of inputs (third repetition of typology E) has good behaviour even for different and various combinations of inputs that vary in an extended range. The performance indicator of the comparison always assumes very high values (close to 100%) for the speed v2 and v3, while it assumes somewhat lower values for the voltages T1 and T4. However, by analysing all the corresponding diagrams, it can be noted that the difference between the experimental trends of T1 and T4 and those generated by the simple model proposed is mainly due to the fluctuations of the experimental values, which may depend on various factors and measurement errors. It can therefore be concluded that the model obtained with the proposed methodology is inserted in the variability field of the same repetitions (see Figure 7) and allows, in all the input combinations analysed, for predicting with good accuracy the behaviour of the complex experimental platform. Considering that the proposed identification and modelling procedure does not need the implementation of any equation, as illustrated above, but only some preliminary experimental tests, it is believed that it may be useful for industrial applications, even complex ones, in which the identification of a sufficiently accurate model must be made quickly and easily.
Finally, to test the identified model’s performance in the Identification for Control (I4C) target, the identified model was used by completing a closed-loop model in Simulink environment [24] with a classical PI controller applied for each of the four subsystems. The scheme of the experimental test in the closed loop is shown in Figure 15, and the feedback is ensured by the tension sensors (T1 and T4) and by the servomotor encoders (v2 and v3).
In Figure 16, the comparison is between the experimental data and the data from the identified model with a target value for the PI controller equal to 10 Newton for T1 and T4 (T1setpoint and T4setpoint in Figure 15) and 0.8 m/s for v2 and v3 (v2setpoint and v3setpoint in Figure 15). The comparison in Figure 16 shows a good correlation between the experimental data and the prediction of the identified model. Similar results were carried out with other experimental tests in the closed loop.

4. Discussion and Conclusions

This work shows how a MIMO experimental web handling system can be identified through simple transfer functions, composed of a few parameters, in the Identification for control (I4C) perspective. The model of the system is estimated in Section 3, using the data of a single test (Test E Repetition 3) in an open loop, allows for simulating the system’s behaviour with step inputs in a wide range of operations. The validation depicted in Figure 12, Figure 13 and Figure 14 demonstrates that the identified model can describe the experimental behaviour of the web transport platform from a steady-state velocity of about 0.5 m/s (Tests G) to the velocity of about 1.4 m/s (Tests A). Therefore, it is possible to predict the system’s behaviour, during both the transient and the steady state, with a certain combination of torque to the motors. Therefore, once the model’s validity in a certain field of operation is verified, the application becomes a powerful tool to carry out an unlimited number of tests, acting on different degrees of freedom and observing the variations in the system’s response. There is a need to underline the extreme simplicity of the procedure and the absence of any complex system of equations to solve, which is different from the classical approach for modelling this kind of system.
Using the MATLAB Simulink environment [17], the estimated model simulated the closed-loop system’s operation. The developed application allows for changing the reference set-points to impose and modify the controller parameters, i.e., the gains of the PI regulators and the smoothing constant of the reference signal. Therefore, once the model’s validity in a certain field of operation is verified, the application becomes a powerful tool to conduct simulated tests, acting on different degrees of freedom and observing variation in the system’s response.
Simulating numerous tests in the same conditions as those performed in the laboratory and comparing the simulated and measured data, it was found that, as for the open loop, even the estimated model in this case provides a good approximation of the behaviour of the system, as long as the operative conditions do not change significantly from the point of operation of the test used for identification. As a possible future development, it is suggested that a more complex system model be implemented with parameters that depend on the work point. However, it is not required to act in a wide field of operation in most cases, so that this safe methodology has great chances for industrial applications, which often demand fixed work points.
In the future, an application that performs tuning of the controller could be implemented, which, once the model is identified and set-points are chosen, establishes the gains that allow for reaching the set-points quickly but without any overshoots that could damage the web.

Author Contributions

Conceptualization, N.I.G. and T.S.; methodology, N.I.G. and T.S.; software, N.I.G. and T.S.; validation, N.I.G. and T.S.; formal analysis, N.I.G. and T.S.; investigation, N.I.G. and T.S.; resources, N.I.G. and T.S.; data curation, N.I.G. and T.S.; writing—original draft preparation, N.I.G.; writing—review and editing, N.I.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Lucrezia Petruzzo for her important contribution to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Akatsuka, M. Models of Longitudinal Web behaviour incorporating Visco-Elasticity. Mitshubishi Juko Giho 1991, 28, 364–367. [Google Scholar]
  2. Benlatreche, A.; Knittel, D.; Osterbag, E. Robust Decentralised control strategies for large-scale web handling systems. Control. Eng. Pract. 2008, 16, 736–750. [Google Scholar] [CrossRef]
  3. Koc, H.; Knittel, D.; de Mathelin, M.; Abba, G. Modeling and robust control of winding systems for elastic webs. IEEE Trans. Control. Syst. Technol. 2002, 10, 197–208. [Google Scholar] [CrossRef]
  4. Laroche, E.; Knittel, D. An improved linear fractional model for robustness analysis of a winding system. Control. Eng. Pract. 2005, 13, 659–666. [Google Scholar] [CrossRef]
  5. Pagilla, P.R.; Siraskar, N.B.; Dwivedula, R.V. Decentralized control of web processing lines. IEEE Trans. Control Syst. Technol. 2007, 15, 106–117. [Google Scholar] [CrossRef]
  6. Giannoccaro, N.I.; Messina, A.; Sakamoto, T. Updating of a lumped model for an experimental web tension control system using a multivariable optimisation method. Appl. Math. Model. 2010, 34, 671–683. [Google Scholar] [CrossRef] [Green Version]
  7. Giannoccaro, N.I.; Nishida, T.; Sakamoto, T. Decentralized Control Performances of an Experimental Web Handling System. Int. J. Adv. Robot. Syst. 2012, 12, 141. [Google Scholar] [CrossRef]
  8. Giannoccaro, N.I.; Sakamoto, T.; Uchitomi, I. Performance of a web transport system with tension control subsystem using speed as control input. In Proceedings of the IECON 2015-41st Annual Conference of the IEEE Industrial Electronics Society, Yokohama, Japan, 9–12 November 2015; pp. 1–6.
  9. Roisum, D.R. The Mechanism of Web Handling; Tappi Press: Peachtree Corners, GA, USA, 1998; pp. 1–4. [Google Scholar]
  10. Ljung, L. Identification for Control: Simple Process Model. In Proceedings of the 41st IEEE Conference on Decision and Control, Las Vegas, NV, USA, 10–13 December 2002; Volume 4, pp. 4652–4657. [Google Scholar]
  11. Thi, L.T.; Mahn, C.N.; Nguyen, T.L. A control approach to web speed and tension regulation of web transport systems based on dynamic surface control. J. Control. Autom. Electr. Syst. 2021, 32, 573–581. [Google Scholar] [CrossRef]
  12. Hou, H.; Nian, X.; Xu, S.; Sun, M.; Xiong, H. Robust decentralised control for large-scale web-winding systems: A linear matrix inequality approach. Trans. Inst. Meas. Control. 2017, 39, 953–964. [Google Scholar] [CrossRef]
  13. Hou, H.; Nian, X.; Xu, S. Decentralized guaranteed cost control with H∞ performance for large scale web-winding system. Asian J. Control. 2022, 24, 459–473. [Google Scholar] [CrossRef]
  14. Yu, X.; Liao, F.; Li, L. H∞ decentralised output feedback preview tracking control via state observer for a class of nonlinear disturbed interconnected discrete systems. Int. J. Control. 2020. [Google Scholar] [CrossRef]
  15. Kubohira, K.; Sakamoto, T.; Isamu, Y. Closed-Loop Identification and Instability of Web Transport System. In Proceedings of the IEEE International Conference on Industrial Technology (ICIT), Seville, Spain, 17–19 March 2015. [Google Scholar]
  16. Sano, A.; Sun, L. Algorithm and Implementation of Direct Closed loop Identification via Output Over-Sampling Scheme. In Proceedings of the 3rd International Conference on Innovative Computing Information and Control, Dalian, China, 18–20 June 2008. [Google Scholar]
  17. Sano, A.; Sun, L. Output over-sampling approach to direct closed-loop identification and its performance. IFAC Proc. Vol. 2009, 15, 687–692. [Google Scholar]
  18. Sun, L.; You, F. Machine learning and data-driven techniques for the control of smart power generation systems: An uncertainty handling perspective. Engineering 2020, 7, 1239–1247. [Google Scholar] [CrossRef]
  19. Fu, F.; Wang, D.; Li, W.; Li, F. Data-driven fault identifiability for discrete-time dynamic systems. Int. J. Syst. Sci. 2020, 51, 404–412. [Google Scholar] [CrossRef]
  20. Garnier, H.; Mensler, M.; Richard, A. Continuous-time domain model identification from sample data: Implemented issues and performance evaluation. Int. J. Control 2003, 76, 1336–1357. [Google Scholar] [CrossRef]
  21. The Mathworks, Inc. System Identification Toolbox for Use with Matlab; The Mathworks, Inc.: Natick, MA, USA, 2020. [Google Scholar]
  22. Giannoccaro, N.I.; Manieri, G.; Martina, G.; Sakamoto, T. Implementation and validation of a lumped model for an experimental multi-span web transport system. In Proceedings of the Australian and New Zealand Control Conference-ANZCC, Gold Cost, Australia, 17–20 December 2017; pp. 1389–1394. [Google Scholar]
  23. Giannoccaro, N.I.; Manieri, G.; Martina, G.; Sakamoto, T. Genetic algorithm for decentralised PI controller tuning of a multi-span web transport system based on overlapping decomposition. In Proceedings of the Australian and New Zealand Control Conference-ANZCC, Gold Cost, Australia, 17–20 December 2017; pp. 994–998. [Google Scholar]
  24. The Mathworks, Inc. Simulink for Use with Matlab; The Mathworks, Inc.: Natick, MA, USA, 2020. [Google Scholar]
Figure 1. The experimental web transport system.
Figure 1. The experimental web transport system.
Systems 10 00017 g001
Figure 2. Scheme of the experimental web transport system with an indication of the different sections and measured parameters.
Figure 2. Scheme of the experimental web transport system with an indication of the different sections and measured parameters.
Systems 10 00017 g002
Figure 3. Functional diagram of the transport line.
Figure 3. Functional diagram of the transport line.
Systems 10 00017 g003
Figure 4. OPP film roll: (a) image of the film roll; (b) details of the experimental system.
Figure 4. OPP film roll: (a) image of the film roll; (b) details of the experimental system.
Systems 10 00017 g004
Figure 5. Resulting model in the Laplace domain without frictional terms: (a) model and subsystems; (b) dynamical model as a black box between inputs and outputs.
Figure 5. Resulting model in the Laplace domain without frictional terms: (a) model and subsystems; (b) dynamical model as a black box between inputs and outputs.
Systems 10 00017 g005
Figure 6. Scheme of open-loop tests.
Figure 6. Scheme of open-loop tests.
Systems 10 00017 g006
Figure 7. Outputs of Test E: 5 repetitions.
Figure 7. Outputs of Test E: 5 repetitions.
Systems 10 00017 g007
Figure 8. Output of the estimation of transfer function models compared to real outputs.
Figure 8. Output of the estimation of transfer function models compared to real outputs.
Systems 10 00017 g008
Figure 9. Simulated response and experimental response (third repetition Test E).
Figure 9. Simulated response and experimental response (third repetition Test E).
Systems 10 00017 g009
Figure 10. Simulated response (identified from third repetition Test E) and experimental response (third repetition Test E).
Figure 10. Simulated response (identified from third repetition Test E) and experimental response (third repetition Test E).
Systems 10 00017 g010
Figure 11. Diagram of validation of the proposed identification procedure.
Figure 11. Diagram of validation of the proposed identification procedure.
Systems 10 00017 g011
Figure 12. Validation of the identified model with the experimental test combinations of Inputs A (a) and B (b).
Figure 12. Validation of the identified model with the experimental test combinations of Inputs A (a) and B (b).
Systems 10 00017 g012
Figure 13. Validation of the identified model with the experimental test combinations of Inputs C (a) and D (b).
Figure 13. Validation of the identified model with the experimental test combinations of Inputs C (a) and D (b).
Systems 10 00017 g013
Figure 14. Validation of the identified model with the experimental test combinations of Inputs F (a) and G (b).
Figure 14. Validation of the identified model with the experimental test combinations of Inputs F (a) and G (b).
Systems 10 00017 g014
Figure 15. Scheme of the closed-loop experimental test.
Figure 15. Scheme of the closed-loop experimental test.
Systems 10 00017 g015
Figure 16. Validation of the identified model in a closed-loop (CL) test.
Figure 16. Validation of the identified model in a closed-loop (CL) test.
Systems 10 00017 g016
Table 1. Description of the physical properties of the platform.
Table 1. Description of the physical properties of the platform.
SectionSymbolValue
UNWINDER
Radiusr13.26 × 10−2 [m]
Moment of InertiaJ14.42 × 10−3 [kg m2]
Span of the sectionL10.75 [m]
LEADING
Radiusr22.5 × 10−2 [m]
Moment of InertiaJ26.10 × 10−3 [kg m2]
Span of the sectionL21.2 [m]
DRAW ROLL
Radiusr31.5 × 10−2 [m]
Moment of InertiaJ36.1 × 10−3 [kg m2]
Span of the sectionL31.25 [m]
WINDER
Radiusr43.22 × 10−2 [m]
Moment of InertiaJ44.39 × 10−3 [kg m2]
Table 2. Physical properties of the OPP film.
Table 2. Physical properties of the OPP film.
PropertiesSymbolValue
Widthw0.3 [m]
ThicknessTh40 [μ]
Cross-sectional areaA1.2 × 10−5 [m2]
Densityρ910 [kg/m3]
Young modulusE9.8 × 109 [N/m2]
Tensile strengthσR32 × 106 [N/m2]
Yield strengthσY22 × 106 [N/m2]
Viscosity coefficientη1.5 × 109 [Ns/m2]
Table 3. Experimental test executed in open loop.
Table 3. Experimental test executed in open loop.
Name of Testu1 [V]u2 [V]u3 [V]u4 [V]
A0.350.080.080.5
B0.20.060.060.4
C0.20.080.080.3
D0.30.080.080.4
E0.40.080.080.55
F0.20.060.060.3
G0.10.060.060.18
Table 4. Poles and zeros of the estimated transfer functions (first repetition Test E), 2 poles, 1 zero.
Table 4. Poles and zeros of the estimated transfer functions (first repetition Test E), 2 poles, 1 zero.
Transfer FunctionZeroPole 1Pole 2
G 11 e −0.0146−4.28330
G 22 e −0.0158−0.2194−0.1041
G 33 e −0.0162−0.2185−0.1053
G 44 e −0.2666−5.8543−0.1943
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Giannoccaro, N.I.; Sakamoto, T. Implementation and Validation of a Simple Direct Identification Method for an Experimental Multi-Span Web Transport System. Systems 2022, 10, 17. https://doi.org/10.3390/systems10010017

AMA Style

Giannoccaro NI, Sakamoto T. Implementation and Validation of a Simple Direct Identification Method for an Experimental Multi-Span Web Transport System. Systems. 2022; 10(1):17. https://doi.org/10.3390/systems10010017

Chicago/Turabian Style

Giannoccaro, Nicola Ivan, and Tetsuzo Sakamoto. 2022. "Implementation and Validation of a Simple Direct Identification Method for an Experimental Multi-Span Web Transport System" Systems 10, no. 1: 17. https://doi.org/10.3390/systems10010017

APA Style

Giannoccaro, N. I., & Sakamoto, T. (2022). Implementation and Validation of a Simple Direct Identification Method for an Experimental Multi-Span Web Transport System. Systems, 10(1), 17. https://doi.org/10.3390/systems10010017

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