The mathematical model of the improved system of the seat with adjustable pressure profile
Tien Tran Xuan^{1} , Dong Nguyen Phu^{2}
^{1}Department of Mechatronics, Faculty of Mechanical Engineering, Hung Yen University of Technology and Education, Hung Yen, Vietnam
^{2}Department of Vehicles and Engines, Faculty of Mechanical Engineering, Technical University of Liberec, Liberec, Czech Republic
^{1}Corresponding author
Mathematical Models in Engineering, Vol. 6, Issue 2, 2020, p. 7992.
https://doi.org/10.21595/mme.2019.21211
Received 5 December 2019; received in revised form 11 December 2019; accepted 19 December 2019; published 30 June 2020
Following a patented solution, a seat which is possible to change its stiffness was created. The seat contains an actively controlled pneumatic spring element (the PSE). For the requirement of working faster and more precisely, an improvement was applied. This article deals with derivation of mathematical model of the improved PSE system used for subsequent analysis. The model is considered as a mixed model which is a combination of singlediscipline subsystems as mechanical, electrical, fluid and control ones. The simulations are carried out for varied input parameters and both the system parameters and system characteristics are calculated. The results describe the behavior of the improved system in two modes of controller setup: constant pressure and constant stiffness under static and dynamic conditions
Keywords: pneumatic spring, mathematical model, improved system, transmission of acceleration.
1. Introduction
The seat with adjustable pressure profile is shown in Fig. 1(a). A PSE inserted inside the cushion allows changing the stiffness of cushion by a system of pneumatic actuator and a controller [1]. The PSE is inserted inside the car seat cushion (as shown in Fig. 1(b)). The mathematical model of the original system was presented in the article [2]. For the requirement of working faster and more precisely, the improved system is built. The detailed scheme of the system is in Fig. 2.
Fig. 1. The seat with a PSE inserted inside
a) The seat in reality
b) A simplified scheme
From the point of view that the PSE affects the human body locally (in contact with buttock), the model is built with the assumption that the mechanical subsystem represents a combination of foam block 100×100×50 mm^{3} and a PSE (see Fig. 3). An additional latex tube is connected to the PSE (see Fig. 4). It can be considered that the volume change is not caused by the PSE but is caused by the additional latex tube when compressed air is supplied. Therefore, the volume change in case of improved PSE is much larger than in case of the original PSE and the regulation of pressure can be enhanced.
Fig. 2. Detailed scheme of the improved system
For convenience of modeling, the improved system is divided into two parts: electropneumatic control subsystem and mechanical subsystem.
Fig. 3. The brick of foam with a PSE inserted inside
Fig. 4. The PSE connected to an additional latex tube
For the distribution of compressed air in and out of the PSE, the system of valves is divided into two pairs. Each pair (input and output) consisting of one proportional valve and one discrete valve. The first pair of input valves (${V}_{1}$, ${V}_{3}$) is used to distribute the compressed air supplied from the reservoir ${R}_{1}$ into the PSE, the second pair of output valves (${V}_{2}$, ${V}_{4}$) is to distribute the air from the PSE to the outlet. A vacuum pump generating output pressure (inside a reservoir ${R}_{2}$) less than 0 kPa is connected to the outlet valves. The proportional valves ${V}_{1}$, ${V}_{2}$ are intended for small airflow in case of small control error (the difference between desired and instant value) which is less than the threshold value set in the control software. The discrete valves ${V}_{3}$, ${V}_{4}$ for faster airflow are opened when the difference between the desired and the instant pressure value is bigger than the threshold value. Sensor PS_{1} (SMC PSE540AR06) measures pressure in the reservoir ${R}_{1}$ (denoted ${p}_{cr}$) and sensor PS_{2} (PSE543AR06) measures pressure inside the PSE (denoted ${p}_{s}$) and one more sensor PS_{3} (PSE543AR06) is added for measuring vacuum pressure inside the reservoir ${R}_{2}$_{}(denoted ${p}_{vr}$). The control software was created in Labview environment. It enables to display the courses of system values in time as desired pressure in PSE, measured unfiltered and filtered pressure in PSE, value of control error, etc. According to [3], the control software is also used for setting of system parameters as upper and lower limit value of pressure in the reservoir, constants of PID controller, selection of working mode (there are 2 modes: constant pressure and constant stiffness). The improved system can work with inlet air pressure up to 250 kPa and vacuum pressure at the output up to –60 kPa. The pressure inside the PSE can be changed in the range [0, 25] kPa.
2. Mathematical model
2.1. Model of the electropneumatic control subsystem
The components of this part are proportional and discrete valves, compressor and PID controller.
2.1.1. Model of the valves
The proportional valves denoted ${V}_{1}$ and ${V}_{2}$ are of type SMCPVQ136M08M5A. They are characterized by dependence of flow rate ${q}_{sj}$ ($j=$1, 2) on pressure difference $\u2206p$ ($\u2206p={p}_{inlet}{p}_{outlet}$) and coil current ${i}_{j}$ (which is supplied to control the proportional valve $j$). The course of these characteristics is taken from the datasheet [3] and it is transformed by means of datafit algorithm into a twoparametric Eq. (1):
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{k}_{11}{i}_{j}\mathrm{\Delta}p+{k}_{02}\mathrm{\Delta}{p}^{2}\mathrm{}+{k}_{30}{{i}_{j}}^{3}+{k}_{21}{{i}_{j}}^{2}\mathrm{\Delta}p+{k}_{12}{i}_{j}\mathrm{\Delta}{p}^{2}+{k}_{03}\mathrm{\Delta}{p}^{3},$
where: ${p}_{inlet}={p}_{cr}$, ${p}_{outlet}={p}_{s}$ for valve ${V}_{1}$ and ${p}_{inlet}={p}_{s}$, ${p}_{outlet}={p}_{vr}$ for valve ${V}_{2}$.
For the purpose of controlling the system, we define a pressure error e as a difference between desired pressure ${p}_{d}$ inside the PSE and instant internal pressure ${p}_{s}$ Eq. (2):
Beside that we define a threshold ${e}_{t}$ which addresses an insensitivity of control system. We have:
The PID controller is used for control of the coil current which is supplied to proportional valves ${V}_{1}$, ${V}_{2}$. Using the equation of PID controller the coil current is a function of the pressure error $e$ and it is presented in the form:
where: ${i}_{j}\left(t\right)$ is the current supplied for controlling the proportional valve $j$ ($j=$ 1, 2).
For cases when we need to realize considerably higher flow rates to make system work faster we need to use the pair of discrete valves ${V}_{3}$, ${V}_{4}$ (in the type of SMCS070B6A) which increase the flow rate in stepchange. According to [4], flow characteristics of the discrete valve include sonic conductance ${C}_{dv}=$0.083 l/(s.bar) and critical pressure ratio ${b}_{dv}=$0.28. The formula for the flow rate calculation of discrete valve is:
in case of choked flow, and:
in case of subsonic flow where $j=$ 3, 4.
${p}_{upstream}$, ${p}_{downstream}$ represents upstream pressure and downstream pressure of the valve respectively. In case of valve ${V}_{3}$, ${p}_{upstream}={p}_{cr}$ and ${p}_{downstream}={p}_{s}$, in case of valve ${V}_{4}$, ${p}_{upstream}={p}_{s}$ and ${p}_{downstream}={p}_{atm}$, and $T$ is the temperature inside the PSE which is assumed to be constant (297 °K).
The discrete valves ${V}_{3}$, ${V}_{4}$ are used for increasing of airflow when control error is recognized as high which is beyond the sum of the threshold (${e}_{t}$) and another parameter (denoted $\mathrm{\Delta}{p}_{s}$) which is set in the control software and called the value of pressure sensitivity. We have:
2.1.2. Mathematical model of the compressed air supply
The process of transporting air from the compressor to the reservoir and distributing compressed air from the reservoir through inlet valves (${V}_{1}$ and ${V}_{3}$) to the PSE is shown in Fig. 5. The compressor characteristics in time and pressure domain are calculated on the experimental basis. Based on the scheme, the total air flowrate ${q}_{cr}$ through the reservoir is then given by:
where: ${q}_{s1}$, ${q}_{s3}$ are flow rates of both valves ${V}_{1}$, ${V}_{3}$, ${q}_{c}$ is supplied flow rate from the compressor.
Fig. 5. The scheme of the process of transmitting compressed air to PSE
When the compressor is working, the flow rate ${q}_{c}$ can be expressed in form of a linear function of the pressure inside the reservoir ${p}_{cr}$:
The constants ${a}_{1}=$ 0.062, ${b}_{1}=$ 275.92 are experimentally determined to correspond to the results in Fig. 6(b).
Pressure ${p}_{cr}$ is measured by sensor $P{S}_{1}$ and it is kept in the interval $\left[{p}_{cr1,}{p}_{cr2}\right]$ by the control system. Values ${p}_{cr1}$, ${p}_{cr2}$ are set in control software as userdefined system parameters. The compressor is switched on in case ${p}_{cr}<{p}_{cr1}$_{}(State = ON) and switched off in case of ${p}_{cr}>{p}_{cr2}$ (State = OFF). The “State” is an internal variable that is automatically set by the controller. Normally, the values of ${p}_{cr1}$ and ${p}_{cr2}$ are set to 190 kPa and 210 kPa respectively. The function describing the flow from the compressor to the reservoir is presented as follows:
Fig. 6. The characteristics of the compressor
a) Dependence of pressure of compressed air on time
b) Dependence of flowrate air on pressure
Fig. 7. The scheme of the process of releasing compressed air in the improved system
2.1.3. Mathematical model of the vacuum pump
The mathematical model of the vacuum pump can be built similarly to the model of the air compressor because of the same operating principle. The total flow rate through the reservoir is then given by:
where ${q}_{s2}$, ${q}_{s4}$ are flow rates of valves ${V}_{2}$, ${V}_{4}$ respectively. In this case ${q}_{s2}$,_{}${q}_{s4}$ is calculated by Eqs. (6) or (7) with ${p}_{upstream}={p}_{s}$ and ${p}_{downstream}={p}_{vr}$.
When the vacuum pump works the flow rate ${q}_{vp}$ can be expressed in the form of a linear function of the pressure ${p}_{r2}$ inside the reservoir ${R}_{2}$:
Fig. 8. The characteristics of the vacuum pump
a) Pressure of vacuum pump – time diagram
b) Flow rate – pressure diagram
The constants ${a}_{2}=$ 0.23778, ${b}_{2}=$ –62.754 are experimentally determined in correspondence with the results in Fig. 8(b).
Pressure ${p}_{vr}$ is measured by sensor $P{S}_{3}$ and is kept in the interval $\left[{p}_{vr1,}{p}_{vr2}\right]$ by the control system. Values ${p}_{vr1}$, ${p}_{vr2}$ are set in control software as userdefined system parameters. The vacuum pump is switched on in case of ${p}_{vr}>{p}_{vr2}$ (State = ON) and switched off in case of ${p}_{vr}<{p}_{vr1}$ (State = OFF). The “State” is an internal variable that is automatically set by the controller. Normally the values of ${p}_{vr1}$ and ${p}_{vr2}$ are set to –60 kPa and –40 kPa respectively. The function describing the flow from the reservoir to the vacuum pump is expressed by:
2.2. Model of the mechanical subsystem
Forces acting on the mass are illustrated in the freebody diagram (see Fig. 9) where $x$, $\dot{x}$, $\ddot{x}$ represent the kinematic displacement excitation, velocity excitation and acceleration excitation, ${\overrightarrow{F}}_{foam}$, ${\overrightarrow{F}}_{p}$ represent the contact forces between the foam block, the latex air spring and the mass sequentially. It is possible to set up the equation of motion of mass $m$ in the form:
For further model derivation, we used this simplified approach.
Fig. 9. Freebody diagram
Fig. 10. Detailed scheme of the mechanical subsystem.
2.2.1. Model of polyurethane foam
Computational scheme of the rheological model the PU foam block is depicted in Fig. 10. It is the result of research taken from [6]. The concept of this lumped parameter model follows from the phenomenological approach. It comprises nonlinear restoring force ${F}_{R}$, damping force ${F}_{d}$ and frictional force ${F}_{f}$ in general force response ${F}_{foam}$. The model definition of PU foam is described in detail as below:
The parameter of PU foam model is shown in Table 1.
Table 1. The parameter of PU foam model
Force component

Parameter

Physical unit

Value


${F}_{b}$

${F}_{b0}$
${c}_{b}$

[N]
[N/m]

80
600


${F}_{dp}$

${S}_{p}$
${p}_{p}$
${n}_{p}$
${h}_{p}$

[m^{2}]
[Pa]
[1]
[m]

0.0095
100
6.2
0.06


${F}_{di}$

${S}_{0i}$
${p}_{0i}$
${n}_{0i}$
${h}_{i}$
${c}_{i}$
${n}_{i}$

[m^{2}]
[Pa]
[1]
[m]
[]
[1]

$i=$1
1.2
100
3
0.05
50
0.2

$i=$2
0.03
3.5
0.05
300
0.2

$i=$3
0.8
2
0.18
300
0.2

${F}_{f}$

${f}_{f0}$
${k}_{1}$
${k}_{2}$
${k}_{3}$

[1]
[s/m]
[–]
[1]

0.05
5000
0.2
1

2.2.2. Model of the combination of the latex air spring and the additional latex tube
The contact force ${F}_{p}$ can be considered as a function depending on the deformation ($x$$z$) and pressure ${p}_{s}$. The change of pressure, deformation and change of contact area is a complex phenomenon, therefore, the dependence of contact force on the deformation and the pressure is determined experimentally. In this way, an experiment was performed on the Instron E3000 machine as shown in Fig. 11. The contact force ${F}_{p}$ is possible to consider in the form of polynomial function:
$+{k}_{11}\left(xz\right)\left({p}_{s}{p}_{atm}\right)+{k}_{02}{\left({p}_{s}{p}_{atm}\right)}^{2}+{k}_{30}{\left(xz\right)}^{3}$
$+{k}_{21}{\left(xz\right)}^{2}\left({p}_{s}{p}_{atm}\right)+{k}_{12}\left(xz\right){\left({p}_{s}{p}_{atm}\right)}^{2}+{k}_{03}.{\left({p}_{s}{p}_{atm}\right)}^{3},$
and the 3D graph of interaction force ${F}_{p}$ is shown in Fig. 12.
According to [7], the internal pressure (${p}_{s}$) depends on the volume (${V}_{ils}$) and the flow rate of supplied air (${q}_{s}$) and is expressed in the form of the differential equation:
where: ${V}_{ils}$ is the total volume of the improved latex air spring, $T$ is the temperature of the air inside the PSE which is assumed to be constant ($T=$297 °K), $\kappa =$1.4 is adiabatic exponent, ${R}_{gas}=$287 J kg^{1 }K^{1}is the gas constant.
Fig. 11. Setup of the experiment
Fig. 12. Relationship between contact force, displacement and pressure
The additional latex tube is connected to the PSE as shown in Fig. 4. The PSE is totally covered with tape so it can be assumed that change of the PSE volume caused by internal pressure is negligible and the total change of volume is completely caused by the change of volume of the additional latex tube.
Fig. 13. Scheme of the latex tube with foam inserted inside
Fig. 14. The scheme of the additional latex tube
Firstly, we can use a simple calculation of volume of the PSE. The volume of the PSE is calculated in the way that is similar to the way and is given by:
where: $R=$0.025 m is the radius of the PSE, $l=$0.3 m is the length of the PSE, ${V}_{foam}=$ 0.15$l$ is the volume of foam (see Fig. 13).
Secondly, we consider the calculation of the volume of the additional latex tube. The additional latex tube is designed with the dimensions shown in Fig. 14. We can recognize that its structure is similar to the structure of the original PSE. The change of volume of additional latex tube depends on pressure ${p}_{s}$ and can be found experimentally. The setup of the experiment is shown in Fig. 15. At first the additional latex tube is filled with compressed air when the value of internal pressure meets the value of desired pressure. At second the additional latex tube is kept in a stable position and its surface is scanned by a 3D scanner machine.
The calculation is programmed in Matlab software. The result is a set of volume values of the additional latex tube in the dependence of pressure. So the volume (${V}_{add}$) depending on internal pressure (${p}_{s}$) can be given in form of a polynomial function Eq. (29):
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{k}_{3}{\left({p}_{s}{p}_{atm}\right)}^{3}+{k}_{4}{\left({p}_{s}{p}_{atm}\right)}^{4}{V}_{foam}.$
By using Matlab a set of coefficients is obtained as shown below:
${k}_{2}=4.111139460129753\mathrm{e}12,\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{k}_{3}=4.030675121010709\mathrm{e}16,$
${k}_{4}=1.163546770201630\mathrm{e}20.$
Fig. 15. The setup of the experiment of the scanning of the surface of the additional latex tube
Fig. 16. Example of a point cloud of the surface of the additional latex tube
Fig. 17. Volumepressure diagram and the fit function
Finally, the total volume of latex air spring in the case of the improved system can be given by:
Then the derivative of the total volume is:
In summary, the mathematical model describing the improved system is given by firstorder differential equations. For the purpose of using Matlab software, it is useful to transform this system of equations into matrix form:
where:
The predefined functions of excitation $z$, $\dot{z}$, $\ddot{z}$ are continuous in an interval of time. The variable ${q}_{s}$ is calculated when the system of equations from Eq. (1) to Eq. (15) describing the function of the electropneumatic control subsystem is being solved. The solution of this system of equations gives as the results describing system response and system characteristics as kinematics of mass ($x,\dot{x},\ddot{x}$), pressure inside the PSE (${p}_{s}$) in case of mechanical subsystem, the flow rate $\left({q}_{s},{q}_{r},{q}_{c,}{q}_{s1},{q}_{s2},{q}_{s3},{q}_{s4}\right)$, pressure $\left({p}_{cr},e,\dots \right)$, supplied coil current of proportional valves ${i}_{j}$ ($j=$1, 2) in case of electropneumatic control subsystem and forces $\left({F}_{foam},{F}_{p},\dots \right)$, volume change of PSE ($\mathrm{\Delta}{V}_{b}$), etc.
To solve the system of equations the initial and boundary conditions were provided. Initial conditions for the system of equations of the improved system are the values of the variables $x$, ${x}_{di}$, $\dot{x}$, ${p}_{s}$, ${q}_{si}$, ${q}_{c}$, ${p}_{c}$, ${V}_{ils}$ at time $t=$ 0 and boundary conditions for the system equations are the values of the parameters ${p}_{cr1}$, ${p}_{cr2}$, ${p}_{vr1}$, ${p}_{vr2}$, ${e}_{s}$, $\mathrm{\Delta}{p}_{s}$.
3. Calculation of the response of the improved system
The examples below show the behavior of the original system when it works under static conditions (without excitation $z\left(t\right)$) and dynamic conditions (with excitation $z\left(t\right)$) for both modes of operation (constant stiffness and constant pressure) and the load $m=$10 kg.
Fig. 18. The calculated response of the improved system
a) Displacement of mass $x$ and diplacement excitation $z$
b) Velocity of mass $vx$ and velocity of excitation $vz$
c) Acceleration of mass $ax$ and acceleration of excitation $az$
d) Forces
e) Volume change of latex air spring
f) Pressure responce $ps$ and desired pressure $pd$
g) The flow rate qs through the PSE
h) Supplied coil current of the valves ${V}_{1}$, ${V}_{2}$, ${V}_{3}$, ${V}_{4}$
4. Calculation of transmission of acceleration
The transmission of acceleration is an important property of efficiency of vibration isolation systems. This is a way how to evaluate the influence of the PSE under dynamic conditions. In this article, we use the signal with continual change of frequency what is advantageous for investigation of nonlinear systems. For such cases socalled sweep sine functions are usually used. According to [8] the harmonic excitation function can be replaced by a sweep sine with variable frequency. Numerous human vibration studies conducted over the past several decades have shown that the human body is sensitive to lowfrequency vibration occurring below 10 Hz [9]. Based on this the range was set to 111 Hz. The amplitude of acceleration was set to constant value of 0.1 g. The excitation function used for experiments under dynamic conditions is expressed in the form:
where: $n$ is the exponent of excitation function ($n=$2), ${f}_{start}$ is start frequency, ${f}_{stop}$ is the stop frequency, $T$ is the total time of the excitation signal, $g$ is the gravitational acceleration and the other parameters are calculated as follows:
Fig. 19 and Fig. 20 show calculated results of the transmission of acceleration in case of using constant pressure mode and constant stiffness mode of control operation. The transmission curves are shown in cases of desired pressure ${p}_{d}\in ${0, 5, 10, 15, 20, 25} kPa in dependence on the frequency of excitation in the range [1, 11] Hz. With a load of 10 kg, peaks of the transmission curves are positioned in the range of [6, 8] Hz and peak values are in range [2.5, 4].
Fig. 19. Transmission of acceleration – constant pressure mode
Fig. 20. Transmission of acceleration – constant stiffness mode
5. Conclusions
The model of the original system consists of the model of foam block which is based on the article [6] and the model of a latex air spring and the model of electropneumatic control which are found experimentally. This model allows to investigate the response of the system to kinematic excitation or the change of the desired pressure inside pneumatic spring. Calculated results describe the behavior of characteristics of the improved system in detail when the improved system works in control modes (constant pressure and constant stiffness) and under static and dynamic conditions. The model is also used for evaluation of vibration isolation effect of the system which is possible to represent in form of the curves of transmission of acceleration. The transmission of acceleration shows the same tendency to shift the peak position toward the higher frequency when the desired pressure increases. This frequency shift is followed by simultaneous increase of peak value.
References
 Cirkl D. Seat. Patent No. 303163, 2012. [Search CrossRef]
 Cirkl D., Tran Xuan T. Modelling of dynamical behavior of pneumatic springmass system. Engineering Mechanics Proceedings, Vol. 24, 2018, p. 865868. [Search CrossRef]
 Cirkl D., Tran Xuan T. Simulation model of seat with implemented pneumatic spring. Vibroengineering Procedia, Vol. 7, 2016, p. 154159. [Search CrossRef]
 Compact Proportional Solenoid Valve. Series PVQ, SMC catalog. [Search CrossRef]
 Port Solenoid Valve. Series S070, SMC catalog. [Search CrossRef]
 Cirkl D., Hrus T. Simulation model of polyurethane foam for uniaxial dynamical compression. Vibroengineering Procedia, Vol. 1, 2013, p. 5358. [Search CrossRef]
 Cirkl D. Mechanical Properties of Polyurethane Foam. Ph.D. Thesis, Technical University of Liberec, Liberec, 2005. [Search CrossRef]
 Sivčák M. Dynamics of the Vibration Isolation System with More Degrees of Freedom. Ph.D. Thesis, Technical University of Liberec, Liberec, 2009. [Search CrossRef]
 Griffin M. J. Handbook of Human Vibration. Academic Press, London, 1990. [Search CrossRef]