⚡ Pluto.jl ⚡

Example: Transport systems

Transport systems are typical candidates for systems. They can be described by the one-dimensional partial differential equation (or simply transport equation)

x(t,z)t+νx(t,z)x=u(t,x)for z(0,L),t>0

with velocity ν>0, left boundary condition x(t,0)=u0(t) for t>0 and the initial data

x(0,z)=xz(0)for z[0,L].

Function u(t,x) is the distributed control input along the complete length and function u0(t) is the boundary input on the left side (x=0).

For simplicity, it is assumed that control input u(t,x) is constant along the whole length.

Solution in Laplace-domain

Further the output is assumed at the right side of the system (x=L) as y(t)=x(t,L). A Laplace transform of the PDE and further calculations lead to

y^(s)=exp(s Tt) u^0(s)+1s[1exp(s Tt)] u^(s)

with Tt=Lν.

Test cases

Two scenarios are considered next:

  • Fluid flow in a pipe with single boundary input u0(t)>0, thus u(t,x)=0 and

  • Conveyor system with distributed input u(t)=u(t,x)>0, thus u0(t)=0.

4.9 ms

Fluid flow in a pipe

The transfer function is noted as

g^0=y^(s)u^0(s)=exp(s Tt) with Tt=Lν.

9.9 μs
32.1 s
L
1.0
1.9 μs
ν
0.2
2.4 μs
Tt
5.0
5.6 μs
g0
DelayLtiSystem{Float64,Float64}

P: StateSpace{Continuous,Float64,Array{Float64,2}}
D = 
 0.0  1.0
 1.0  0.0

Continuous-time state-space model

Delays: [5.0]
161 ms

The ideal delay is approximated by

g^approx=as+a exp(s Tt)

with a>0.

10.0 μs
18.2 ms
1
823 ns
g_approx
DelayLtiSystem{Float64,Float64}

P: StateSpace{Continuous,Float64,Array{Float64,2}}
A = 
 -1.0
B = 
 0.0  1.0
C = 
 1.0
 0.0
D = 
 0.0  0.0
 1.0  0.0

Continuous-time state-space model

Delays: [5.0]
115 μs

Step response

6.0 μs
Tf
10
1.6 μs
2.8 ms

Bode plot

25.9 μs
ω
0.010000000000000002:0.010000000000000002:99.99000000000002
128 ms
168 ms
10.9 ms

Nyquist plot

3.9 μs
110 ms

Conveyor system

The transfer function is noted as

g^0=y^(s)u^0(s)=1s[1exp(s Tt)] with Tt=Lν.

Values for L and ν are the same as above.

6.7 μs
g1
DelayLtiSystem{Float64,Float64}

P: StateSpace{Continuous,Float64,Array{Float64,2}}
A = 
 0.0  0.0
 0.0  0.0
B = 
 1.0  0.0
 0.0  1.0
C = 
 1.0  -1.0
 0.0   0.0
D = 
 0.0  -0.0
 1.0   0.0

Continuous-time state-space model

Delays: [5.0]
4.5 s

Step plot

3.9 μs
35.7 s

Bode plot

3.4 μs
13.7 s

Nyquist plot

3.3 μs
1.8 s

Zuletzt geändert: Montag, 18. Januar 2021, 17:18