Simulation model of seat with implemented pneumatic spring

David Cirkl1 , Tien Tran Xuan2

1, 2Department of Applied Mechanics, Faculty of Mechanical Engineering, Technical University of Liberec, Liberec, Czech Republic

1Corresponding author

Vibroengineering PROCEDIA, Vol. 7, 2016, p. 154-159.
Received 2 August 2016; accepted 4 August 2016; published 31 August 2016

Copyright © 2016 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License

This article deals with description of simulation model of a seat with pneumatic element inserted into its cushioning. Suggested solution comprises built-in pneumatic spring with feedback control which influences the hardness of cushioning. The seat with variable cushioning properties should contribute to the increase of passenger’s comfort. In order to investigate properties of this seat-mass system the mathematical model was created and solved numerically. This model describes dynamic behavior of the electro-pneumatic mechanical system and calculates the displacement and acceleration of mass, pressure inside pneumatic spring, flow rate in controlled solenoid pneumatic valves and needed coil current in case of kinematic excitation of the seat, as well as the change of setting of desired pressure in pneumatic spring.

Keywords: car seat, vibration isolation, polyurethane foam, pneumatic spring.

1. Introduction

Increasing passenger’s comfort is the task inherently connected with using seats in different means of transport or working machines. The travel conditions are varied in dependence on quality of the road surface in case of moving vehicle, its speed, and other circumstances. One of the ways how to increase passenger’s comfort in case of dynamic loading by vibrations might be the use of the cushioning with variable hardness. The derivation of simulation model of such seat is presented in this article.

The seat with the implemented pneumatic spring contains an active vibroisolation element in accordance with the patented solution [1]. This system allows to change the stiffness of seat cushioning by a computer controlled electro-pneumatic feedback circuit. This electro-pneumatic part of the system comprises a compressor, reservoir of compressed air, 4 pneumatic valves, spring element, 2 pressure sensors and electronic control unit (Fig. 1). The pressure in the spring element is changed on user’s demand to chosen value or it is controlled by computer in accordance with chosen time course, for instance harmonic or triangular with predefined mean, amplitude and frequency.

Principally two different modes of operation are possible: constant stiffness and constant pressure. In the simpler constant stiffness mode, the pressure in the pneumatic spring denoted ps is set to chosen desired pressure pd at initial time. After that, the inlet and outlet valves are closed and volume of air inside the spring is then constant even if the system is mechanically loaded or not. In this mode the controlling system allows to set the fixed stiffness characteristics of cushioning. In the constant pressure mode, the electro-pneumatic feedback circuit tries to keep the pressure ps equal to pd for all time of the loading process.

There are two kinds of external loading taken into consideration. The first of them is sitting down of the passenger on the seat, and the second one is subsequent start of kinematical excitation of loaded seat by harmonic signal of displacement.

2. Model of the system

For purpose of creating the simulation model the system is divided in two parts in accordance with Fig. 1. The first one is represented by the controlled electro-pneumatic feedback circuit and the second one, the mechanical system, is represented by cushioning with implemented pneumatic spring.

Fig. 1. The scheme of seat with implemented pneumatic spring

 The scheme of seat with implemented pneumatic spring

2.1. Model of electro-pneumatic feedback circuit

Electro-pneumatic system of pneumatic circuit consists of two kinds of pneumatic valves. The proportional valves denoted V1 and V2 in Fig. 1 are SMC – PVQ13-6M-08-M5-A. They are characterized by dependence of flow rate on pressure difference p between pressure on input p1 and output p2 and by input coil current i. This characteristics (Fig. 2) is taken from the datasheet [2] and by fitting algorithm it was transformed into two-parametric Eq. (1):

q s j = k 00 + k 10 i Δ p + k 01 Δ p + k 20 i 2 + k 11 i . Δ p + k 02 Δ p 2 + k 30 i 3
            + k 21 i 2 Δ p + k 12 i Δ p 2 + k 03 Δ p 3 ,

where j= 1, 2. In the pneumatic system there are also implemented two digital valves V3 and V4 of type SMC-S070B-6A. The function of these additional valves is to increase air flow in or out of pneumatic spring in case the control error between instantaneous pressure ps and desired value pd is too high. According with [3], flow characteristics of digital valve includes sonic conductance C= 0.083 l/( and critical pressure ratio b= 0.28. According with [4], the formula for the flow rate calculation of digital valve is:

Choked flow:

p 2 p 1 b q s j = 600 . C . p 1 . 293 273 + T         l / min .

Subsonic flow:

p 2 p 1 b q s j = 600 . C . p 1 . 1 - p 2 p 1 - b 1 - b 2 293 273 + T         l / min ,

where j= 3, 4 and p1 represents upstream pressure and p2 downstream pressure respectively. The total air flow qs is then given by addition of flow rates of all individual valves qsj as it is expressed by equation:

q s = q s 1 - q s 2 + q s 3 - q s 4 .

For the control of coil current, the PID algorithm is used in form Eq. (5):

i t = K P e t + K I 0 t e τ d τ + K D d e t d t ,

where it is the current controlling the proportional valve, et is the control error given by expression:

e t = p s t p d t .

Fig. 2. Characteristics of proportional valve, from [5]

 Characteristics of proportional valve, from [5]

Fig. 3. Free body diagram of mechanical part

 Free body diagram of mechanical part

2.2. Model of mechanical part of the system

This part of the system is a simplified representation of seat cushioning with implemented pneumatic spring. According with [5], mechanical properties of flexible cushioning are described by restoring force given by nonlinear progressive function (7), and linear damping force given by Eqs. (10-12). Pneumatic spring is in a simplified way represented by closed air cylinder with internal pressure ps and the mass of passenger by value m. Value of pressure inside the cylinder is given by differential Eq. (8) and finally with respect to free body diagram in Fig. 3 it is possible to set up the equation of motion of mass m in form Eq. (9):

F k + F b + F p - m g = m x ¨ ,
p ˙ s = κ q s R T V - κ p s V ˙ V ,
V = S h + z - x V ˙ = S z ˙ - x ˙ ,
F k = k 1 z - x + k 2 z - x 2 + k 3 z - x 3 ,
F b = b 1 z ˙ - x ˙ ,
F p = S p s - p o u t .

Used physical quantities in the system of Eq. (7-12) are: z – excitation, xdisplacement of mass m, Fk – restoring force, Fb – damping force, Fp – force of pneumatic spring element, S – the contact area between mass and cylinder (S= 0.2 m2), ps – pressure inside cylinder, pin – pressure in the reservoir, pout – outside pressure (atmospheric pressure pout= 0.1 MPa), h0= 0.1 m – the height of cylinder at the initial position, V – volume of cylinder, T – temperature of gas inside cylinder, it is assumed to be constant (T= 297 °K), κ= 1.4 – adiabatic exponent, R= 287 J kg-1 K-1 – gas constant.

3. Numerical simulation

The system of equations derived above make up the mathematical model of combined electro-pneumatic mechanical system of the seat. The simulations are performed in Matlab software for two modes of seat operation, first in the mode of constant stiffness and second in the mode of constant pressure. For the purpose of this numerical experiment the excitation zt is given by harmonic function. Simulated process takes 20 seconds in total. At initial time the desired pressure in pneumatic spring element is set to pd= 10 kPa, after that at time 6 s the seat is loaded by mass m= 50 kg which simulates sitting down of passenger, at time 10 s the external harmonic excitation zt is applied as a source of vibration load, finally, but only in case of constant pressure mode at time 16 s the desired value of pressure inside air cylinder is changed to pd= 25 kPa under continuing vibration.

Values of constant used for calculation were k1= 10 000 N/m, k2= 0 N/m2, k3= 1 000 N/m3 in case of restoring force Eq. (10), and b1= 500 N.s/m in case of damping force Eq. (11). Outside pressure pout is equal to atmospheric pressure 100 kPa, pressure inside the air reservoir pin=pout+ 200 kPa. Frequency of excitation zt is 1 Hz and amplitude is 10 mm. Parameters of PID controller were set to KP= 7, KI= 0 and KD= 0, sampling time was dt= 0.01 s in case of constant pressure mode. Then the system is in mode of active control. In case of constant stiffness all pneumatic valves V1-4 are closed therefore total air flow qs is zero and the system is considered as passive one with no feedback control.

The results of simulations are depicted in Fig. 7 and Fig. 8 and comprises instantaneous pressure ps inside pneumatic spring, displacement xt and acceleration at of mass m, total air flow in/out the spring element qs, restoring and damping force of cushioning, force of pneumatic spring element and their summation.

The transmission of acceleration is calculated as ratio of amplitude of acceleration of mass and amplitude of acceleration of excitation (Ax/Az) when excitation is harmonic function. When changing frequency of excitation f, the transmission curve is obtained as in Fig. 9. It shows transmission curves in cases of desired pressure pd {1, 5, 10, 15, 20, 25} kPa in dependence on frequency of excitation in interval [1, 11] Hz.

Fig. 7. Response of the system in constant pressure mode

 Response of the system in constant pressure mode

Fig. 8. Response of the system in constant stiffness mode

 Response of the system in constant stiffness mode

Fig. 9. Transmission of acceleration

 Transmission of acceleration

a) Constant pressure mode

 Transmission of acceleration

b) Constant stiffness mode

4. Conclusions

In this article the simulation model of seat with implemented pneumatic spring is introduced which allows to influence the hardness of the cushioning. The solution which comprises electro-pneumatic circuit with feedback control is described by mathematical model. This model allows to investigate the response of the system to kinematic excitation or the change of the desired pressure inside pneumatic spring. This pressure can be set to constant value or defined by deterministic course. The analysis of the response helps to find suitable model parameters which represent properties of real parts used for seat construction. These parameters are for example size and type of used pneumatic valves or pressure of air inside reservoir and others. The influence of constants of PID controller can be analyzed as well. The model will be used for evaluation of vibration isolation effect of the system which is possible to represent e.g. in form of the curves of transmission of acceleration.


This article was written at the Technical University of Liberec, Faculty of Mechanical Engineering with the support of the Institutional Endowment for the Long Term Conceptual Development of Research Institutes, as provided by the Ministry of Education, Youth and Sports of the Czech Republic in the year 2016.


  1. Cirkl David Seat. Patent No. 303163, 2012. [Search CrossRef]
  2. Compact Proportional Solenoid Valve. Series PVQ, SMC Catalog. [Search CrossRef]
  3. 3 Port Solenoid Valve. Series S070, SMC Catalog. [Search CrossRef]
  4. Kulakowski Bohdan T., Gardner John F., Shearer J. Lowen Dynamic Modeling and Control of Engineering Systems. Cambridge University Press, 2007, p. 219-243. [Search CrossRef]
  5. Cirkl David Mechanical Properties of Polyurethane Foam. Ph.D. Thesis, Technical University of Liberec, 2005, p. 134-136. [Search CrossRef]