An Introduction to Piecewise Smooth Dynamics Paul Glendinning 2024 scribd download
An Introduction to Piecewise Smooth Dynamics Paul Glendinning 2024 scribd download
com
https://textbookfull.com/product/an-introduction-to-
piecewise-smooth-dynamics-paul-glendinning/
OR CLICK BUTTON
DOWNLOAD NOW
https://textbookfull.com/product/an-introduction-to-complex-systems-
society-ecology-and-nonlinear-dynamics-1st-edition-paul-fieguth-auth/
textboxfull.com
https://textbookfull.com/product/an-introduction-to-modeling-neuronal-
dynamics-1st-edition-christoph-borgers-auth/
textboxfull.com
https://textbookfull.com/product/an-introduction-to-symbolic-dynamics-
and-coding-2nd-edition-douglas-lind/
textboxfull.com
https://textbookfull.com/product/approaching-photography-an-
introduction-to-understanding-photographs-3rd-edition-paul-hill/
textboxfull.com
Introduction To Lagrangian Dynamics Aron Wolf Pila
https://textbookfull.com/product/introduction-to-lagrangian-dynamics-
aron-wolf-pila/
textboxfull.com
https://textbookfull.com/product/market-research-in-practice-an-
introduction-to-gaining-greater-market-insight-3rd-edition-paul-hague/
textboxfull.com
https://textbookfull.com/product/practical-programming-an-
introduction-to-computer-science-using-python-3-6-3rd-edition-paul-
gries/
textboxfull.com
https://textbookfull.com/product/introduction-to-software-testing-2nd-
edition-paul-ammann/
textboxfull.com
https://textbookfull.com/product/introduction-to-averaging-dynamics-
over-networks-1st-edition-fabio-fagnani/
textboxfull.com
Advanced Courses in Mathematics CRM Barcelona
Paul Glendinning
Mike R. Jeffrey
An Introduction
to Piecewise
Smooth Dynamics
Advanced Courses in Mathematics
CRM Barcelona
Managing Editor:
Enric Ventura
An Introduction to Piecewise
Smooth Dynamics
Editors for this volume:
Elena Bossolini, Technical University of Denmark
J. Tomás Lázaro, Universitat Politècnica de Catalunya
Josep M. Olm, Universitat Politècnica de Catalunya
Paul Glendinning Mike R. Jeffrey
School of Mathematics Department of Engineering Mathematics
University of Manchester University of Bristol
Manchester, UK Bristol, UK
Mathematics Subject Classification (2010): 34A36, 34C23, 37E05, 37G10, 37E35, 37G35
This book is published under the imprint Birkhäuser, www.birkhauser-science.com by the registered company
Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Foreword
v
Contents
Introduction 1
1 Piecewise-smooth Flows
by Mike R. Jeffrey 3
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 History and applications . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3 Inclusions and combinations . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4 Types of dynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.5 Switching layers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
1.6 Multiple switches . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.7 Codimension r switching layers . . . . . . . . . . . . . . . . . . . . . . 32
1.8 Codimension r sliding . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
1.9 Boundary equilibrium bifurcations . . . . . . . . . . . . . . . . . . . . 36
1.10 Stability, equivalence, and bifurcation . . . . . . . . . . . . . . . . . . 38
1.11 Discontinuity-induced phenomena . . . . . . . . . . . . . . . . . . . . . 42
1.12 Determinacy-breaking . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
1.13 Hidden attractors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
1.14 Hidden bifurcations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
1.15 Moving forward . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
2 Piecewise-smooth Maps
by Paul Glendinning 55
2.1 Introduction to maps . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.2 Smooth theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
2.3 Piecewise-smooth maps of the interval . . . . . . . . . . . . . . . . . . 69
2.4 Lorenz maps and rotations . . . . . . . . . . . . . . . . . . . . . . . . . 83
2.5 Gluing bifurcations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
2.6 Piecewise-smooth maps of the plane . . . . . . . . . . . . . . . . . . . 97
2.7 Periodic orbits and resonance . . . . . . . . . . . . . . . . . . . . . . . 102
2.8 Robust chaos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
2.9 Two-dimensional attractors . . . . . . . . . . . . . . . . . . . . . . . . . 114
2.10 Challenges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
Bibliography 123
vii
Introduction
The study of nonsmooth dynamics has a long history, but has been recently en-
livened by the introduction of new techniques, the discovery of new phenomena,
and the explosion of new practical disciplines applying nonsmooth dynamical mod-
eling. The more formal name for our field is piecewise-smooth dynamics, concern-
ing the theory and applications of how dynamical systems that are smooth almost
everywhere behave when affected by discontinuities at isolated thresholds.
This volume presents an informal course at graduate level, aimed at mathe-
maticians, scientists, and engineers, studying models that involve a discontinuity,
or studying the theory of nonsmooth systems for its own sake. These notes are de-
rived from lectures given at the Advanced School on Piecewise-Smooth Dynamical
Systems, organized at the Centre de Recerca Matemàtica (Bellaterra, Catalonia),
during April 11–15, 2016. The School was organized with support from the Engi-
neering and Physical Sciences Research Council (EPSRC), the Societat Catalana
de Matemàtiques (SCM), WIRIS, and the Institute for Mathematics and its Appli-
cations (IMA). The course comes in two parts: flows and maps (or continuous and
discrete time systems). Each part introduces the applications and main theoretical
techniques, of piecewise-smooth dynamics, some long-established and others newly
emerging.
An introduction to the dynamics of piecewise-smooth flows is authored by
Mike Jeffrey. The chapter starts with a few key points in the history and applica-
tions of discontinuities from mechanics and control systems, and some of the key
figures in setting up a general theory for differential equations with ‘discontinuous
right-hand sides’. The methods of inclusions versus combinations are discussed,
after which we introduce the elementary dynamics of crossing and sliding at a
discontinuity surface. The analytical methods of switching layers and layer vari-
ables are presented for one switch and then for multiple switches. The concepts of
discontinuity-induced phenomena and determinacy breaking are introduced, along
with the definitions of stability, equivalence, and bifurcation in piecewise-smooth
flows. We end by introducing some of the novel attractors and bifurcations that
offer a hint at what remains to be discovered, including the present state of zoology
of singularities in the plane.
An introduction to the dynamics of piecewise-smooth maps is authored by
Paul Glendinning. After a discussion of why piecewise-smooth maps are inter-
1
2 Introduction
esting, the course moves into their phenomenology, and reviews some techniques
from smooth theory. The main topics are then: basic stability analysis, piece-
wise monotonic maps of the interval, rotation-like maps, gluing bifurcations (aka
Big Bang bifurcations and period-adding), an introduction to renormalization,
decomposition theorems and a brief guide to kneading theory, piecewise-smooth
maps of the plane including the Lozi and border collision normal forms, piecewise
isometries, bounding regions, periodic orbits and resonance, robust chaos, and
two-dimensional attractors. The course concludes with a discussion of challenges
in higher dimensions, particularly concerning periodic orbits, N -dimensional at-
tractors, and analogies with smooth cases.
Piecewise-smooth Flows
1.1 Introduction
This course is about the geometry of piecewise-smooth dynamical systems. The
solutions of a system of ordinary differential equations, such as
ẋ = f (x ), (1.1)
where x = (x1 , x2 , . . . , xn ) is some n-dimensional vector or variable, and f is an
n-dimensional vector field, can be pictured as trajectories (or orbits) in space (for
example, Rn or some subset of it). Those trajectories are organized by various
singularities, separatrices, and invariant sets, whose geometry can be studied in
great generality. A loss of continuity in the differential equations greatly adds to
the richness of that geometry.
.
or x vector field
parameters x2
dx = f(x,a,b,...)
dt state, dependent variable x1
time, independent variable
switching
surface
Figure 1.1: The vector field f tells us the velocity with which some flow evolves
through a state x . If the flow is non-differentiable, then f is discontinuous.
The mathematician Alexei Fedorovich Filippov set out methods for solving
piecewise-smooth differential equations, and these have been adopted as standard.
Recently, however, with the discovery of new singularities and singular phenomena,
we have understood that more is needed. This series of lectures will provide readers
with the tools to delve more deeply into the world of nonsmooth dynamics.
We will focus on:
geometry of piecewise-smooth vector fields,
general methods for solving and analysing them,
their key notions of stability and bifurcations.
Some important current topics that we will not cover, but you may wish to
look up, include:
special cases: e.g., piecewise linear or continuous non-differentiable vector
fields, hybrid/impact systems;
modeling non-idealities: e.g., smoothing, noise, delay, etc., and the various
types of ‘regularization’, a word that comes up a lot as an open problem in
current piecewise-smooth systems.
You may also want to look into simulation methods: piecewise-smooth sys-
tems require special consideration when simulating. There are event detection
routines built into Matlab and Mathematica. There are ways of making continua-
tion tools like AUTO or MatCont work with discontinuities (often by smoothing
them out), including the AUTO-derived TcHat. Filippov’s solving methods (we’ll
discuss these in the course) have even been built into tools like Mathematica. But
these are not up-to-date in the many advances in theory that we have seen in the
last decade. Use them all with care and critical judgement.
f2 f1 f1
f2
f− f+
f3 f4 f4
f3
Figure 1.2: Discontinuities along one threshold in two dimensions, two thresholds
in two dimenions, and two thresholds in three dimensions.
Locally, solutions take certain simple forms. Away from the switching surface
they will be smooth unique curves, thanks to the existence and uniqueness of
solutions of differentiable dynamical systems. At the surface, however, they might
cross through the discontinuity, or they might slide along it, following a sliding
vector field.
crossing sliding
The switching surface has a lower dimension than the surrounding space, so
if sliding occurs, the space the solutions occupy changes, see fig. 1.4. This results
in non-uniqueness. When sliding is:
attractive, solutions stick to the switching surface, and then many solutions
will all evolve onto the same trajectory in forward time. The history of any
point (shown in the left part of fig. 1.4) in attractive sliding is therefore
non-unique — this is common in physics as mechanical ‘sticking’.
repulsive, solutions escape a switching surface, and then the sliding mode has
many possible future trajectories (shown in the right part of fig. 1.4). The
future of any point in repulsive sliding is therefore non-unique — determinacy
is broken.
The only other things we must add to these are the elementary singularities.
In differentiable vector fields, the commonly encountered singularity is a steady
state or ‘equilibrium’, where f = 0 in a system ẋ = f . At a discontinuity we
encounter a new kind of steady state, a sliding equilibrium, fig. 1.5, (called a
6 Chapter 1. Piecewise-smooth Flows
sticking determinacy-breaking
Figure 1.4: Sliding along the switching surface leads to sticking (if attractive) or
determinacy-breaking (if repulsive).
Keeping these few elements in mind will be of great assistance in gaining some
intuition for the (at first strange) terrain of what we have come to call informally:
Nonsmoothland.
1.2. History and applications 7
in 1785 Coulomb added clarity and depth to the former theories, especially
those of da Vinci and Amontons;
in 1950 Bowden and Tabor clarified the role of contact area and surface
roughness, showing that friction does depend on the true contact area, which
is often less than the apparent contact area (i.e., the full contact surface);
1950 and into the present, the field of Tribology remains active, studying
speed-dependence and friction memory, the different friction characteristics
at different scales of length, time, or speed, and remains important, for ex-
ample, to the health, automotive, and energy industries.
This is a very fleeting description, meant only to highlight the complexity of a
contact force that we use and interact with every moment of every day. The most
enduring part of this story — that friction depends principally on load only — is
encapsulated in the now commonly adopted law
. .
x x
right left
μNFN slip slip μNFN
v v
Figure 1.7: The switching force of friction.
.
x
incoming
.
−rx
recoil
Figure 1.8: The discontinuity in velocity resulting from impact with a hard wall.
Because this mixes up continuous and discrete time evolution, it crosses the
divide between the two chapters of this course, and opens the way to study of
the incredibly general field of sytems that are hybrids of maps and flows (see, e.g.,
[12, 21, 99]), which we will not delve in to here (another example of a hybrid
system is cellular mitosis: a cellular organism grows continuously until it triggers
a discrete change in mass corresponding to mitosis; see, e.g., [27]).
However, usually a map just represents a jump through some fast continuous
interval of motion. In the case above, the restitution map ẋ ↦ −rẋ represents a
jump through a continuous impact phase. Alternatively we could model this as
⎧
⎪
⎪1 if x > c,
mẍ = f − k step(x − c), step(x − c) = ⎨ (1.5)
⎪
⎪0 if x < c,
⎩
introducing a large wall stiffness k. This now has continuous time solutions that
are non-differentiable at x = c, and falls back under our topic of piecewise-smooth
flows.
ẋ = f i (x ) on x ∈ Ri
1 We have written “1990 onwards” but one of these references is earlier. Before this Teixeira
phrased this work as ‘divergent diagrams’ or ‘pairings of fields and functions’ to circumvent the
early skepticism towards discontinuous dynamical systems. Teixeira now leads one of the most
fervent and successful communities in geometric piecewise-smooth dynamical theory.
1.3. Inclusions and combinations 11
or, in components,
(ẋ1 , ẋ2 , . . . , ẋn ) = (f1i , f2i , . . . , fni ) . (1.6)
The index i is taken from a set of labels identified with the regions Ri . The
boundary between these regions is the switching surface Σ, such that the full
space is given by Σ ∪ R1 ∪ R2 ∪ ⋯
Example 1.3.1. For two regions we could take labels i = 1, 2, or i = +, −, so that,
e.g., ẋ = f + for x ∈ R+ and ẋ = f − for x ∈ R− (fig. 1.9, left). For four regions we
may take labels i = 1, 2, 3, 4, or alternatively i = ++, +−, −+, −−.
The figure shows a system with: two regions, two regions where the boundary
between them is a corner, or four regions.
In each open region Ri the vector fields f i are smooth, so we can apply
standard dynamical systems methodology. This course is about what then happens
at the discontinuity boundary — the switching surface Σ — between them.
ẋ ∈ F where f i (x ) ∈ F ∀i ∶ lim f (x ′ ) = f i (x )
x ′ →x
f+ f+
R+
F
Σ
f−
f−
R−
Figure 1.10: The set F at a switching surface Σ between vector fields f + on region
R+ and f − on region R− .
12 Chapter 1. Piecewise-smooth Flows
.
x
kx
.
cx
FN
μNFN v
mg
Figure 1.11: A block on a moving belt.
ẋ = y, ẏ = − mc y − k
m
x − λ FmN . (1.8)
∂(ẋ, ẏ) 1 0 m
= ( ),
∂(x, y) m −k −c
2 Upper semicontinuity is an extension of the notion of continuity for sets. It says
supb∈F (x ) ρ(b, a) → 0 as p′ → p for a ∈ F (p′ ) and b ∈ F (p) where ρ(b, a) = inf a∈A, b∈B ∣a − b∣
with ∣a − b∣ the distance between a and b.
1.3. Inclusions and combinations 13
√
has eigenvalues (−c ± c2 − 4km)/2m with negative real part (giving a focus for
small damping c2 < 4km and a node otherwise); these are basic results using the
local stability theory of smooth dynamical systems (see, e.g., [6, 69]). Because
there is only a single equilibrium, at y = 0, the region in which it exists (in right
slip ẋ > v or left slip ẋ < v) depends on the sign of the constant v.
. v<0
x
right
stick slip
v
left
slip
x
. v>0
x right
slip
stick
v
left
slip
x
Figure 1.12: Dynamics of the oscillating block shown in the phase plane (x, ẋ).
As we trace out trajectories in this phase space, we find that the block can go
from left slip to right slip or vice versa, and it can go from right or left slip to stick,
then back to slip. All trajectories are eventually attracted to the equilibrium.
Sticking is represented by motion along ẋ = v, which means that the block
and surface have matching speeds and are therefore stuck together.
We need a method to describe this sticking motion. This comes from solving
the inclusion.
⎧
⎪
⎪f + (x ) if σ(x ) > 0,
ẋ = ⎨ −
⎪
⎪f (x ) if σ(x ) < 0,
⎩
f + (x ) = f (x ; +1),
i.e., i = ± ⇔ λ = ±1. (1.9)
f − (x ) = f (x ; −1),
We then define
λ = sign (σ(x )) for σ(x ) ≠ 0,
(1.10)
λ ∈ (−1, +1) for σ(x ) = 0.
This last definition is extremely important to remember. The multiplier λ takes
values ±1 outside the switching surface, and lies inside the open interval (−1, +1) on
the switching surface. We shall see later how to fix the value inside this interval.
We will often write for shorthand just λ = sign(σ), understanding the value of
sign(σ) to lie in (−1, +1) for σ = 0.
So the switching of the multiplier λ causes the vector field f (x ; λ) to switch
between functional forms f + (x ) and f − (x ). The dependence of f on λ tells us
about how the jump between those modes occurs.
The function f (x , λ) combines f + and f − into a single expression where the
switching is described by λ, so we may refer to f (x , λ) as a combination.
1.3.5 Combinations
The combination f (x ; λ) is differentiable with respect to x and λ (by which we
mean that the partial derivaties ∂x∂
f and ∂λ
∂
f exist). The discontinuity is now
encoded in the multiplier λ.
E.g., convex combination (see Filippov [29, Def. 4a, pp. 50–52])
ẋ = 12 (1 + λ) f + (x ) + 12 (1 − λ) f − (x ); (1.11)
ẋ = 1
2
(1 + λ) f + (x ) + 12 (1 − λ) f − (x ) + (λ2 − 1)g (x , λ). (1.12)
(λ2−1)g
f+ f+
R+ f(fil)
Σ
f−
f− ∇σ
R−
Figure 1.13: The Filippov convex combination (dotted line) written f (fil) =
1
2
(1 + λ) f + + 12 (1 − λ) f − , and a non-convex combination (dashed curve) formed
by hidden term (λ2 − 1)g .
f+
(λ2−1)(0,−1)
f+
f(fil)
f−
f− x1+x2=0
Figure 1.14: The Filippov convex combination (dotted line) and a non-convex
combination (dashed curve) for theorem 1.3.4 with quadratic dependence on λ.
f+
R+
Σ
Σ
f−
R− R+
f−
R−
Figure 1.15: A crossing solution (left) corresponding to the inclusion from fig. 1.10
and fig. 1.13. We can concatenate over many discontinuities (right).
d
x (t) = Φt (x0 ) ∶ Φt (x0 ) = f (Φt (x0 ); λ) , Φ0 (x0 ) = x0 .
dt
We may write the constituent flows in the regions R+ , R− , as Φ+t , Φ−t , where
d ±
Φ (x0 ) = f (Φt (x0 ); ±1) .
dt t
d
x (t) = ΦΣ
t (x0 ) ∶ σ(x (t)) = 0 on Σ (1.15)
dt
In the convex combination (1.11), since σ̇ = ẋ ⋅ ∇σ = f ⋅ ∇σ, sliding
requires that both vector fields f ± point towards (or both away from) Σ,
therefore sliding occurs where
(f + ⋅ ∇σ)(f − ⋅ ∇σ) < 0 on Σ, (1.16)
and satisfies
1+λ + 1−λ −
0 = σ̇ = ẋ ⋅ ∇σ = { f (x ) + f (x )} ⋅ ∇σ (1.17)
2 2
(f − + f + ) ⋅ ∇σ
⇒ λ = λΣ ≡ (1.18)
(f − − f + ) ⋅ ∇σ
(f − ⋅ ∇σ)f + − (f + ⋅ ∇σ)f −
⇒ ẋ = f Σ ≡ . (1.19)
(f − − f + ) ⋅ ∇σ
In three dimensions this is compactly written as
∇σ × (f + × f − )
fΣ ≡ . (1.20)
(f − − f + ) ⋅ ∇σ
For sliding to exist we must have λΣ ∈ (−1, +1) by (1.10), and indeed it is easy
to show that (1.18) satisfies this if and only if (1.16) is satisfied. Sliding can
occur outside (1.16) for a nonlinear switching system like (1.12), occurring
wherever σ̇ = σ = 0 can be satisfied for [∈] (−1, +1).
Later it will be useful to observe that, from (1.19), we have
fΣ =f+ if f + ⋅ ∇σ = 0, (1.21)
Σ − −
f =f if f ⋅ ∇σ = 0, (1.22)
meaning that when the sliding vector field becomes equal to f or f − at +
points where their normal components vanish, i.e., where they lie tangent to
Σ. These two scenarios correspond to λΣ = +1 and λΣ = −1 by (1.18), and
therefore also apply to (1.12).
18 Chapter 1. Piecewise-smooth Flows
Defining sliding dynamics therefore generally requires solving for the value
of λ that provides motion along σ = 0:
Definition 1.4.2. The sliding dynamics is given by
ẋ = f Σ (x ) on σ(x ) = 0 (1.23)
where f Σ (x ) is the sliding vector field, given by
f Σ (x ) ≡ f (x ; λΣ ) such that f (x ; λΣ ) ⋅ ∇σ = 0 and λΣ ∈ (−1, +1) (1.24)
on σ(x ) = 0.
Example 1.4.3. (i) The orbit
⎧
⎪
⎪Φ− (x0 ) if t < τ,
x (t) = ⎨ t+ (1.25)
⎪
⎪ Φ (x ) if t > τ,
⎩ t−τ 1
crosses from R− to R+ at time τ ;
(ii) The orbit
⎧
⎪
⎪Φ− (x0 ) if t < τ,
x (t) = ⎨ tΣ (1.26)
⎪
⎪ Φ (x ) if t > τ,
⎩ t−τ 1
starts in R− , then sticks to Σ at time τ , sliding thereafter; where Φ±t are the
flows of ẋ = f ± , and x1 = φ−τ (x0 ) ∈ Σ.
+
Φt
cross x1 x1 Σ
Φt
− slide
Φt −
Φt
x0 x0
Figure 1.16: Examples of crossing and sliding.
f−=(1,1) f+=(3,1)
R− x1=0 R+
Figure 1.17: A piecewise constant example of crossing. We will revisit this example
later to see it is not as trivial as it first appears.
fΣ=(0,1)
f−=(1,1) f+=(−1,1)
R− x1=0 R+
Figure 1.18: A piecewise constant example of sliding.
These are stationary points of the sliding flow on the switching surface.
(ii) Tangencies of the vector field to Σ are points where:
f + ⋅ ∇σ = 0 ⇒ ẋ = f + (1.36)
or similarly for f − . These are stationary points with respect to the switching
surface, of the flows outside the surface. It is easy seen using (1.11) or (1.12)
that they correspond to points where f (x ; ±1) = 0, and so for example λΣ =
±1 in (1.18).
Example 1.4.6 (Sliding equilibria). In two dimensions consider λ = sign x1 with
Clearly the λ = +1 system in x1 > 0 and the λ = −1 system in x1 < 0 are non-
vanishing. In both cases, to find sliding modes we solve 0 = σ̇ = ẋ1 = 12 (1 + λ)(−1) +
1
2
(1 − λ)(+1) = λ. Hence in (1.37)
λΣ = 0 ⇒ ẋ2 = 12 (b − x2 ),
λΣ = 0 ⇒ (ẋ2 , ẋ3 ) = 12 (b − x2 + x3 , c − x2 ),
x1 x1 x3
x2 x2
Figure 1.19: Sliding equilibria in two or three dimensions. For (1.37) to (1.38) the
x1 direction is vertical.
Example 1.4.7 (Tangencies). Figure 1.20 below shows quadratic tangencies of types
we call ‘visible’ (curving away from Σ) or ‘invisible’ (curving towards Σ).
In higher dimensions, sets of visible and invisible tangencies meet at higher
order tangencies, like the cubic tangency, called a ‘cusp’, illlustrated in fig. 1.21.
The examples shown in fig. 1.20 and fig. 1.21 are from the equations
(ẋ1 , ẋ2 ) = 1
2
(1 + λ)(−1, −x1 ) + 12 (1 − λ)(0, 1), λ = sign(x2 ),
(ẋ1 , ẋ2 ) = 1
2
(1 + λ)(+1, −x1 ) + 12 (1 − λ)(0, 1), λ = sign(x2 ),
(ẋ1 , ẋ2 , ẋ3 ) = 1
2
(1 + λ)(1, 0, x21 + x2 ) + 12 (1 − λ)(0, 0, 1), λ = sign(x3 ),
but we could also have cubic, quartic, etc. order.
visible
invisible
Σ Σ
Figure 1.20: Visible and invisible tangencies in the upper vector field. On the
switching surface Σ, the full line indicates the sliding region found using a convex
combination, but nonlinear dependence on λ can let this can bleed out as indicated
by the dashed region, while the remainder is crossing (dotted).
At the simple (‘fold’) tangencies in fig. 1.20 we have that the upper (λ = +1
vector field) satisfies ẋ2 = −x1 which vanishes at x1 , creating a tangency. The sign
of the second derivative determines its curvature, in the first example ẍ2 = −ẋ1 = +1
implies the flow is curving upwards, creating a visible tangency, and in the second
example ẍ2 = −ẋ1 = −1 implies the flow is curving downwards, creating an invisible
tangency
At the cusp, for the upper vector field we have ẋ3 = x21 + x2 which vanishes
√
on x3 = 0 along a curve x1 = ± −x2 (existing for x2 < 0). The second derivative is
√ √
ẍ3 = 2ẋ1 x1 + ẋ2 = 2x1 = ±2 −x2 , meaning the half-curve x1 = + −x2 is a family
√
of visible folds while the half-curve x1 = − −x2 is a family of invisible folds. The
point x1 = x2 = x3 = 0 where ẍ3 = 0 is the cusp.
22 Chapter 1. Piecewise-smooth Flows
vis.
inv.
cusp
Σ
Figure 1.21: A cusp in the upper vector field. On the switching surface Σ, the dark
shaded region indicates the sliding region found using a convex combination, but
nonlinear dependence on λ can let this can bleed out as indicated by the light
shaded region, while the remainder is crossing (unshaded).
Above we gave a necessary condition for crossing, but it does not guarantee
crossing (it is not necessary and sufficient). The region indicating by the full line
in fig. 1.20 and by shading fig. 1.21 on Σ shows where both vector fields point
towards the surface, f + ⋅ ∇σ < 0 < f − ⋅ ∇σ, so sliding must occur. In the remaining
region trajectories may cross through (and in Filippov’s convention they do cross
through), but in nonlinear combinations it is possible for the sliding region to
‘bleed out’ into the crossing region, indicated by the dashed/light-shaded region
on Σ.
We’ve hinted above that crossing or sliding on a switching surface depends
on convexity, i.e., whether the dependence on λ is linear. Most courses would
assume linearity and hasten onward. Let’s hold back a little and consider nonlinear
dependence a little further.
so σ(x ) = x1 . Let τ > 0. The obvious crossing solution is again x1 (t) = (t − τ )(2 +
sign(t − τ )), x2 (t) = t + x20 , but turns out to be wrong in the presence of the
nonlinear term. Let’s look for sliding modes by solving σ̇ = σ = 0 on Σ (x1 = 0):
Thus, there are two sliding modes. It turns out that either modes gives sliding
dynamics (ẋ1 , ẋ2 ) = (0, 1) for this simple example.
Discovering Diverse Content Through
Random Scribd Documents
about donations to the Project Gutenberg Literary Archive
Foundation.”
• You comply with all other terms of this agreement for free
distribution of Project Gutenberg™ works.
1.F.
1.F.4. Except for the limited right of replacement or refund set forth
in paragraph 1.F.3, this work is provided to you ‘AS-IS’, WITH NO
OTHER WARRANTIES OF ANY KIND, EXPRESS OR IMPLIED,
INCLUDING BUT NOT LIMITED TO WARRANTIES OF
MERCHANTABILITY OR FITNESS FOR ANY PURPOSE.
Please check the Project Gutenberg web pages for current donation
methods and addresses. Donations are accepted in a number of
other ways including checks, online payments and credit card
donations. To donate, please visit: www.gutenberg.org/donate.
Most people start at our website which has the main PG search
facility: www.gutenberg.org.
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
textbookfull.com