Sliding Modes in Filippov Systems.
|
|
- Benjamin Stewart
- 6 years ago
- Views:
Transcription
1 in Filippov Systems. School of Mathematics, Georgia Tech Joint work with Luciano Lopez (Univ. Bari) Geneva Conference In honor of Ernst Hairer 60th birthday
2 References Filippov [DEs with discontinuous RHS, 1988] Utkin [Sliding modes in control & optimization, 1991] Yakubovich Leonov Gelig [Stability of stationary sets in control systems with discontinuous nonlinearities, 2004] di Bernardo-Budd-Champneys-Kowalczyk [P.wise smooth dynamical systems, 2008] Acary-Brogliato [Numerical methods for nonsmooth DynSys, 2008] Hairer-Norsett-Wanner [v.1, 1992, pp. 196 & ss]... and many more in specialized literature... Control (Bartolini, Usai, Fridman, Johansson,...), in Optimization (Cortes, Glazos-Hui-Zak, Venets,...), in Dynamical Systems (Leine et al., Y. Kuznetsov, Mueller, Zou-Kupper-Beyn, Dercole, Rinaldi, Galvanetto,...), and Numerical: D.Stewart, Mannshardt, Shampine-Thompson, Gear, Enright, DAE,... [D. Lopez].
3 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes.
4 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes. 2. Entry and exit.
5 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes. 2. Entry and exit. 3. Projections and sliding vector field.
6 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes. 2. Entry and exit. 3. Projections and sliding vector field. 4. Sliding on intersection of p 2 surfaces.
7 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes. 2. Entry and exit. 3. Projections and sliding vector field. 4. Sliding on intersection of p 2 surfaces. (a) Case of p = 2: [DiLo] vector and [Alexander-Seidman] choice.
8 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes. 2. Entry and exit. 3. Projections and sliding vector field. 4. Sliding on intersection of p 2 surfaces. (a) Case of p = 2: [DiLo] vector and [Alexander-Seidman] choice. (b) p > 2, generalization of [DiLo].
9 Outline of Talk 1. Which non-smooth? Filippov theory, sliding modes. 2. Entry and exit. 3. Projections and sliding vector field. 4. Sliding on intersection of p 2 surfaces. (a) Case of p = 2: [DiLo] vector and [Alexander-Seidman] choice. (b) p > 2, generalization of [DiLo]. (c) p = 3: Extension of [AS]-method. p > 3....
10 Which non-smooth? Focus on piecewise smooth systems (switched systems). One surface case. x (t) = f(x(t)) = { f1 (x(t)), x S 1, f 2 (x(t)), x S 2, x(0) = x 0 R n, R n = S 1 Σ S 2, Σ = {x R n h(x) = 0}, h : R n R, S 1 = {x R n h(x) < 0}, S 2 = {x R n h(x) > 0}. Assume f 1 smooth on S 1 Σ, f 2 smooth on S 2 Σ, h is smooth on Σ and h(x) 0 for all x Σ.
11 ... Which non-smooth?... S 2 n(x) S 2 n(x) x T T x x x S 1 Σ Σ S1
12 Filippov theory Issue: f(x) not defined if x Σ!
13 Filippov theory Issue: f(x) not defined if x Σ! Work with differential inclusion: f 1 (x(t)) x S 1 x (t) F(x(t)) = {(1 α)f 1 + αf 2, α [0, 1]} x Σ. f 2 (x(t)) x S 2
14 Filippov theory Issue: f(x) not defined if x Σ! Work with differential inclusion: f 1 (x(t)) x S 1 x (t) F(x(t)) = {(1 α)f 1 + αf 2, α [0, 1]} x Σ. f 2 (x(t)) x S 2 Filippov Solution. An absolutely continuous function x : [0, τ) R n is a solution if for almost all t [0, τ) it holds that x (t) F(x(t)).
15 ... Filippov... Theorem Let f 1,2 be C 1, f 1 f 2 also C 1 (on Σ), and h be C 2. If, at any point x Σ we have that at least one of n T (x)f 1 (x) > 0 and n T (x)f 2 (x) < 0 holds, then there exists a unique Filippov solution from each initial condition.
16 ... Filippov... Theorem Let f 1,2 be C 1, f 1 f 2 also C 1 (on Σ), and h be C 2. If, at any point x Σ we have that at least one of n T (x)f 1 (x) > 0 and n T (x)f 2 (x) < 0 holds, then there exists a unique Filippov solution from each initial condition. What is there? Non-trivial case: x Σ.
17 Crossing and Sliding Transversal Intersection [n T (x)f 1 (x)] [n T (x)f 2 (x)] > 0. Leave Σ to S 1, if n T (x)f 1 (x) < 0, or S 2, if n T (x)f 1 (x) > 0.
18 Crossing and Sliding Transversal Intersection [n T (x)f 1 (x)] [n T (x)f 2 (x)] > 0. Leave Σ to S 1, if n T (x)f 1 (x) < 0, or S 2, if n T (x)f 1 (x) > 0. Sliding Mode [n T (x)f 1 (x)] [n T (x)f 2 (x)] < 0.
19 Crossing and Sliding Transversal Intersection [n T (x)f 1 (x)] [n T (x)f 2 (x)] > 0. Leave Σ to S 1, if n T (x)f 1 (x) < 0, or S 2, if n T (x)f 1 (x) > 0. Sliding Mode [n T (x)f 1 (x)] [n T (x)f 2 (x)] < 0. Attracting. Have existence and uniqueness: [n T (x)f 1(x)] > 0 and [n T (x)f 2(x)] < 0.
20 Crossing and Sliding Transversal Intersection [n T (x)f 1 (x)] [n T (x)f 2 (x)] > 0. Leave Σ to S 1, if n T (x)f 1 (x) < 0, or S 2, if n T (x)f 1 (x) > 0. Sliding Mode [n T (x)f 1 (x)] [n T (x)f 2 (x)] < 0. Attracting. Have existence and uniqueness: [n T (x)f 1(x)] > 0 and [n T (x)f 2(x)] < 0. Repulsive. Not covered by result. No uniqueness: [n T (x)f 1(x)] < 0 and [n T (x)f 2(x)] > 0.
21 Sliding f 2 S 2 n(x) n(x) S 2 f 2 x T x x T x Σ f 1 S 1 Σ f 1 S 1
22 Sliding f 2 S 2 n(x) n(x) S 2 f 2 x T x x T x Σ f 1 S 1 Σ f 1 S 1 How is attracting sliding motion taking place on Σ?
23 Sliding f 2 S 2 n(x) n(x) S 2 f 2 x T x x T x Σ f 1 S 1 Σ f 1 S 1 How is attracting sliding motion taking place on Σ? Filippov vector field: f F (x) = (1 α(x))f 1 (x) + α(x)f 2 (x) α(x) such that n T (x)f F (x) = 0: α(x) = n T (x)f 1 (x) n T (x)(f 1 (x) f 2 (x)).
24 Entry and Exit As we reach Σ, or are sliding on it, we can have several distinct possible behaviors: Transversal crossing, attractive Sliding, smooth exits. And these can occur at first order, or higher order.
25 Entry and Exit As we reach Σ, or are sliding on it, we can have several distinct possible behaviors: Transversal crossing, attractive Sliding, smooth exits. And these can occur at first order, or higher order. As far as sliding motion is concerned, entries are less troublesome than exits. Still, for the sake of this presentation, we will tacitly assume that (insofar as sliding motion) entries and exits happen in the generic way. This means that:
26 Entry and Exit As we reach Σ, or are sliding on it, we can have several distinct possible behaviors: Transversal crossing, attractive Sliding, smooth exits. And these can occur at first order, or higher order. As far as sliding motion is concerned, entries are less troublesome than exits. Still, for the sake of this presentation, we will tacitly assume that (insofar as sliding motion) entries and exits happen in the generic way. This means that: (a) Entries are non-tangential (we hit Σ);
27 Entry and Exit As we reach Σ, or are sliding on it, we can have several distinct possible behaviors: Transversal crossing, attractive Sliding, smooth exits. And these can occur at first order, or higher order. As far as sliding motion is concerned, entries are less troublesome than exits. Still, for the sake of this presentation, we will tacitly assume that (insofar as sliding motion) entries and exits happen in the generic way. This means that: (a) Entries are non-tangential (we hit Σ); (b) Exits are tangential. That is, we leave the surface with the sliding vector field being (aligned with) one of the vector fields f i s we have outside of the surface.
28 Entry and Exit As we reach Σ, or are sliding on it, we can have several distinct possible behaviors: Transversal crossing, attractive Sliding, smooth exits. And these can occur at first order, or higher order. As far as sliding motion is concerned, entries are less troublesome than exits. Still, for the sake of this presentation, we will tacitly assume that (insofar as sliding motion) entries and exits happen in the generic way. This means that: (a) Entries are non-tangential (we hit Σ); (b) Exits are tangential. That is, we leave the surface with the sliding vector field being (aligned with) one of the vector fields f i s we have outside of the surface. We will also implicitly assume that entry/exit points are isolated events (no accumulation of events).
29 Projections and tangents: Revisit sliding vector field Beginning from the end... sliding must keep us on Σ.
30 Projections and tangents: Revisit sliding vector field Beginning from the end... sliding must keep us on Σ. Projection onto T Σ (x), π Σ = [I nn T ]. Define admissible sliding vector field: f S (x) = (1 β(x))(π Σ f 1 )(x) + β(x)(π Σ f 2 )(x), 0 β 1.
31 Projections and tangents: Revisit sliding vector field Beginning from the end... sliding must keep us on Σ. Projection onto T Σ (x), π Σ = [I nn T ]. Define admissible sliding vector field: f S (x) = (1 β(x))(π Σ f 1 )(x) + β(x)(π Σ f 2 )(x), 0 β 1. Choose β of the form β(x) = a 1n T (x)f 1 (x) + b 1 n T (x)f 2 (x) + c 1 a 2 n T (x)f 1 (x) + b 2 n T (x)f 2 (x) + c 2 and find coefficients by imposing 1st Order Exit Conditions: { β(x) = 0 when n T (x)f 1 (x) = 0, β(x) = 1 when n T (x)f 2 (x) = 0, β(x) = n T (x)f 1(x) n T (x)(f 1(x) af 2(x)), a > 0 is free; a = 1 is Filippov.
32 Two and more sliding surfaces Consider now the case of two surfaces of discontinuity: Σ 1 := {x R n : h 1 (x) = 0}, Σ 2 := {x R n : h 2 (x) = 0}. Assume h 1 (x) 0, for all x Σ 1, and h 2 (x) 0, for all x Σ 2. Assume n 1 (x) and n 2 (x) linearly independent for all x Σ 1 Σ 2. Issue: To define sliding vector field on Σ 1 Σ 2.
33 ... Two or more... f 3 Σ 2 n 1 (x) f 4 T x x n (x) 2 f 1 T x Σ 1 01 f 2 Figure: Two surfaces
34 ... Two or more... Attractive sliding... w i (x) = [ ] w 1 i (x) wi 2(x) = [ ] n T 1 f i (x) n T 2 f = N i(x) T f i (x), i = 1, 2, 3, 4.
35 ... Two or more... Attractive sliding... w i (x) = [ ] w 1 i (x) wi 2(x) = Attractive sliding motion means ( ) ( ) > 0 > 0 w 1 =, w > 0 2 =, w < 0 3 = [ ] n T 1 f i (x) n T 2 f = N i(x) T f i (x), i = 1, 2, 3, 4. ( ) < 0, w > 0 4 = ( ) < 0. < 0 f 3 Σ 2 n 1 (x) f 4 T x x n (x) 2 f 1 T x Σ f 2
36 ... Two or more... Filippov vector field, ambiguity: 4 f F (x) = α i f i (x), i=1 4 α i = 1, 0 α i 1. i=1 Impose: One relation too few. n T 1 (x)f F (x) = 0, n T 2 (x)f F (x) = 0.
37 ... Two or more... Filippov vector field, ambiguity: 4 f F (x) = α i f i (x), i=1 4 α i = 1, 0 α i 1. i=1 Impose: One relation too few. n T 1 (x)f F (x) = 0, n T 2 (x)f F (x) = 0. Special cases treated in [Filippov], [Piiroinen Kuznetsov], [Utkin], [Stewart],....
38 ... Two or more... Geometric approach: Use tangent vectors. Let Π(x) be projection onto T x, x Σ 1 Σ 2 : Π = I N(N T N) 1 N T.
39 ... Two or more... Geometric approach: Use tangent vectors. Let Π(x) be projection onto T x, x Σ 1 Σ 2 : Π = I N(N T N) 1 N T. Define f = 4 λ i v i (x), i=1 with λ i 0, and 4 i=1 λ i = 1. v i (x) = Πf i (x)
40 ... Two or more... Geometric approach: Use tangent vectors. Let Π(x) be projection onto T x, x Σ 1 Σ 2 : Π = I N(N T N) 1 N T. Define f = 4 λ i v i (x), i=1 with λ i 0, and 4 i=1 λ i = 1. To determine λ i... v i (x) = Πf i (x)
41 ... Two or more... We want to satisfy exit conditions v i = f i λ i = 1, i = 1,...,4.
42 ... Two or more... We want to satisfy exit conditions v i = f i λ i = 1, i = 1,...,4. We end up proposing to take µ i λ i = 4 i=1 µ, i = 1, 2, 3, 4, i µ i = [ 4 ] 1/3 j=1,j i at j w j [ 4 j=1,j i at j w j [ 1 a 1 =, a 1] 2 = [ 1 1 ] 1/3 a T i w i i = 1,...,4 ] [ [ 1 1, a 3 =, a 1] 4 = 1],
43 ... Two or more... We want to satisfy exit conditions v i = f i λ i = 1, i = 1,...,4. We end up proposing to take µ i λ i = 4 i=1 µ, i = 1, 2, 3, 4, i µ i = [ 4 ] 1/3 j=1,j i at j w j [ 4 j=1,j i at j w j [ 1 a 1 =, a 1] 2 = [ 1 1 ] 1/3 a T i w i i = 1,...,4 ] [ [ 1 1, a 3 =, a 1] 4 = 1], Notice: a T 1 y 1 > 0, a T 2 y 2 < 0, a T 3 y 3 < 0, a T 4 y 4 < 0 and we satisfy 1st order exit conditions: v i = f i λ i = 1, i = 1,...,4.
44 ... Two or more... [Alexander-Seidman] interpolation approach. ([AS] call it sigmoid blending).
45 ... Two or more... [Alexander-Seidman] interpolation approach. ([AS] call it sigmoid blending). The idea is to obtain a vector field on the intersection, call it f B, by forming a bilinear interpolant amongst the four vector fields and then imposing the orthogonality conditions: f B (x) = (1 α(x))(1 β(x))f 1 (x) + (1 α(x))β(x)f 2 (x)+ α(x)(1 β(x))f 3 (x) + α(x)β(x)f 4 (x), where α, β [0, 1] [0, 1] must be found by solving the nonlinear system n T 1 (x)f B (x) = 0, nt 2 (x)f B (x) = 0.
46 ... Two or more... [Alexander-Seidman] interpolation approach. ([AS] call it sigmoid blending). The idea is to obtain a vector field on the intersection, call it f B, by forming a bilinear interpolant amongst the four vector fields and then imposing the orthogonality conditions: f B (x) = (1 α(x))(1 β(x))f 1 (x) + (1 α(x))β(x)f 2 (x)+ α(x)(1 β(x))f 3 (x) + α(x)β(x)f 4 (x), where α, β [0, 1] [0, 1] must be found by solving the nonlinear system n T 1 (x)f B (x) = 0, nt 2 (x)f B (x) = 0. As it turns out, under the attractivity conditions, this system is uniquely solvable (see [AS]; their proof makes heavy use of planar geometry).
47 ... Two or more... On the intersection of p surfaces (n p + 1)... Let us generalize the geometric approach of [DiLo].
48 ... Two or more... On the intersection of p surfaces (n p + 1)... Let us generalize the geometric approach of [DiLo]. Now, have p surfaces in R n, Σ 1, Σ 2,..., Σ p, with n p + 1, and the vectors h i (x) 0 are linearly independent for x p i=1 Σ i. Locally, we have 2 p vector fields f i, i = 1,...,2 p. Set N = [n 1, n 2,...,n p ] and as before let the projection onto N be Π = I N(N T N) 1 N T. Let w i = N T f i, for i = 1,...,2 p.
49 ... Two or more... On the intersection of p surfaces (n p + 1)... Let us generalize the geometric approach of [DiLo]. Now, have p surfaces in R n, Σ 1, Σ 2,..., Σ p, with n p + 1, and the vectors h i (x) 0 are linearly independent for x p i=1 Σ i. Locally, we have 2 p vector fields f i, i = 1,...,2 p. Set N = [n 1, n 2,...,n p ] and as before let the projection onto N be Π = I N(N T N) 1 N T. Let w i = N T f i, for i = 1,...,2 p. Attractive sliding means that the components of the vectors w i, i = 1,...,2 p, take all signs patterns in the 2 p combinations of [±1, ±1,..., ±1]. E.g., for appropriate ordering of the vector fields, hence of the w i s, we can assume that all components of w 1 are positive, the last component of w 2 is negative all the others being positive, etc..
50 ... Two or more... On the intersection of p surfaces (n p + 1)... Let us generalize the geometric approach of [DiLo]. Now, have p surfaces in R n, Σ 1, Σ 2,..., Σ p, with n p + 1, and the vectors h i (x) 0 are linearly independent for x p i=1 Σ i. Locally, we have 2 p vector fields f i, i = 1,...,2 p. Set N = [n 1, n 2,...,n p ] and as before let the projection onto N be Π = I N(N T N) 1 N T. Let w i = N T f i, for i = 1,...,2 p. Attractive sliding means that the components of the vectors w i, i = 1,...,2 p, take all signs patterns in the 2 p combinations of [±1, ±1,..., ±1]. E.g., for appropriate ordering of the vector fields, hence of the w i s, we can assume that all components of w 1 are positive, the last component of w 2 is negative all the others being positive, etc.. We still want to satisfy first order exit conditions (i.e., tangential).
51 ... Two or more...
52 ... Two or more... Propose to take the sliding vector field on the intersection as: f S (x) = 2 p i=1 For the µ i s we take µ i = λ i v i (x), λ i = [ 2 p j=1,j i at j w j] m [ 2 p j=1,j i at j w j] m a T i w i µ i 2 p j=1 µ, i = 1,..., 2 p. j, i = 1,...,2 p, m = 1 2 p 1, and a i = [±1, ±1,..., ±1] T, with the signs chosen so that sgn(a i ) j = sgn(w i ) j < 0 for i = 2,...,2 p, j = 1,...,p, and the signs of the entries of a 1 being the same as the signs of the entries of w 1.
53 ... Two or more... Extension of [AS] to the case of attractive sliding motion on the intersection of p = 3 surfaces: ẋ =(1 α)(1 β)(1 γ)f 1 (x) + (1 α)(1 β)(1 + γ)f 2 (x) +(1 α)(1 + β)(1 γ)f 3 (x) + (1 α)(1 + β)(1 + γ)f 4 (x) +(1 + α)(1 β)(1 γ)f 5 (x) + (1 + α)(1 β)(1 + γ)f 6 (x) +(1 + α)(1 + β)(1 γ)f 7 (x) + (1 + α)(1 + β)(1 + γ)f 8 (x)
54 ... Two or more... Extension of [AS] to the case of attractive sliding motion on the intersection of p = 3 surfaces: ẋ =(1 α)(1 β)(1 γ)f 1 (x) + (1 α)(1 β)(1 + γ)f 2 (x) +(1 α)(1 + β)(1 γ)f 3 (x) + (1 α)(1 + β)(1 + γ)f 4 (x) +(1 + α)(1 β)(1 γ)f 5 (x) + (1 + α)(1 β)(1 + γ)f 6 (x) +(1 + α)(1 + β)(1 γ)f 7 (x) + (1 + α)(1 + β)(1 + γ)f 8 (x) and α, β, γ [0, 1]: n T 1 (x)ẋ = n T 2 (x)ẋ = n T 3 (x)ẋ = 0.
55 ... Two or more... Not at all clear that this nonlinear system is solvable.
56 ... Two or more... Not at all clear that this nonlinear system is solvable. However, we recently proved that: System has a solution (α, β, γ) [0, 1] 3.
57 ... Two or more... Not at all clear that this nonlinear system is solvable. However, we recently proved that: System has a solution (α, β, γ) [0, 1] 3. Solution is isolated (locally unique).
58 ... Two or more... Not at all clear that this nonlinear system is solvable. However, we recently proved that: System has a solution (α, β, γ) [0, 1] 3. Solution is isolated (locally unique). Extend construction, and proof of existence of solution, to the case of the intersection of any number of surfaces.
59 Numerical Method Each particular state of the system is integrated with an appropriate numerical method, and the events (where structural changes in the system occur) are located accurately: Event driven method.
60 Numerical Method Each particular state of the system is integrated with an appropriate numerical method, and the events (where structural changes in the system occur) are located accurately: Event driven method. Must decide how to proceed at event points (sliding or crossing). Must check exit conditions (1st and higher order). Must decide how to slide. Difficulties as we get close to the surface(s) of discontinuity, since vector fields may not extend past surfaces,
61 Numerical Method Each particular state of the system is integrated with an appropriate numerical method, and the events (where structural changes in the system occur) are located accurately: Event driven method. Must decide how to proceed at event points (sliding or crossing). Must check exit conditions (1st and higher order). Must decide how to slide. Difficulties as we get close to the surface(s) of discontinuity, since vector fields may not extend past surfaces, Standard integration off Σ. Integration on Σ using projected RK methods. [We use Lagrange s multiplier technique and a simplified Newton iteration.]
62 Numerical Method Each particular state of the system is integrated with an appropriate numerical method, and the events (where structural changes in the system occur) are located accurately: Event driven method. Must decide how to proceed at event points (sliding or crossing). Must check exit conditions (1st and higher order). Must decide how to slide. Difficulties as we get close to the surface(s) of discontinuity, since vector fields may not extend past surfaces, Standard integration off Σ. Integration on Σ using projected RK methods. [We use Lagrange s multiplier technique and a simplified Newton iteration.] Alternative: Regularization Obvious simplifications in the theory. But, regularization may lead to losing the characteristics of the original problem.
63 Numerical Method Each particular state of the system is integrated with an appropriate numerical method, and the events (where structural changes in the system occur) are located accurately: Event driven method. Must decide how to proceed at event points (sliding or crossing). Must check exit conditions (1st and higher order). Must decide how to slide. Difficulties as we get close to the surface(s) of discontinuity, since vector fields may not extend past surfaces, Standard integration off Σ. Integration on Σ using projected RK methods. [We use Lagrange s multiplier technique and a simplified Newton iteration.] Alternative: Regularization Obvious simplifications in the theory. But, regularization may lead to losing the characteristics of the original problem. Small integration steps are usually required during the numerical simulation of regularized system due to the large derivatives that replace the changes in the structure of the system.
64 Example: Sliding and Chattering on Intersection. A linear problem [IEEE-AC, 2002] x 1 y < 0, (t) = Ax(t) + bu(t) y(t) = c T with u = sign(y) = [ 1, 1] y = 0, x(t) 1 y > A = , c = 0 0, and (i) b = 3β 3β 2, or (ii) b = 1 2β, β = 1/5. β 3 β 2
65 ... Sliding and Chattering... (i) The trajectory has a portion on the sliding set of order 1, the plane x 1 = 0. The solution slides on the plane x 1 = 0, twice. Past the transient, the behavior is periodic, and the lack of differentiability when the solution hits the sliding plane is apparent.
66 ... Sliding and Chattering... (i) The trajectory has a portion on the sliding set of order 1, the plane x 1 = 0. The solution slides on the plane x 1 = 0, twice. Past the transient, the behavior is periodic, and the lack of differentiability when the solution hits the sliding plane is apparent. (ii) The trajectory exhibits chattering (frequent switching) about the sliding set of order 2, x 1 = x 2 = 0. The solution appears to slide on the sliding set of order 2, x 1 = x 2 = 0, but it does not.
67 ... Sliding and Chattering x1 3 x x x x3 1 x x x 10 4
68 ... Sliding and Chattering x 1 x x x 10 5 x blowup x blowup x x x 10 4
69 Example: Double stick-slip [Galvanetto]. Mechanical system composed by two blocks on a moving belt. The velocity of the belt is constant, v. Each block is connected to a fixed support and to the other block by elastic springs.
70 Example: Double stick-slip [Galvanetto]. Mechanical system composed by two blocks on a moving belt. The velocity of the belt is constant, v. Each block is connected to a fixed support and to the other block by elastic springs. The surface between the blocks and the belt is rough so that the belt exerts a dry friction force on each block that sticks on the belt to the point where the elastic forces due to the springs exceed the maximum static force.
71 Example: Double stick-slip [Galvanetto]. Mechanical system composed by two blocks on a moving belt. The velocity of the belt is constant, v. Each block is connected to a fixed support and to the other block by elastic springs. The surface between the blocks and the belt is rough so that the belt exerts a dry friction force on each block that sticks on the belt to the point where the elastic forces due to the springs exceed the maximum static force. At this point the blocks start slipping and the slipping motion will continue to the point where the velocity of the block will equal that of the belt and the elastic forces will be equilibrated by the static friction force.
72 Example: Double stick-slip [Galvanetto]. Mechanical system composed by two blocks on a moving belt. The velocity of the belt is constant, v. Each block is connected to a fixed support and to the other block by elastic springs. The surface between the blocks and the belt is rough so that the belt exerts a dry friction force on each block that sticks on the belt to the point where the elastic forces due to the springs exceed the maximum static force. At this point the blocks start slipping and the slipping motion will continue to the point where the velocity of the block will equal that of the belt and the elastic forces will be equilibrated by the static friction force. Continuous repetition of this type of motion stick-slip oscillation.
73 ... stick-slip... { m1 x = k 1 x + k 12 (x y) + f k1 (x v), m 2 y = k 2 y + k 12 (y x) + f k2 (y v), The kinetic force has the form f k2 (x v) = βf k1 (x v) with: 1 δ f k1 (x 1 γ(x v) + δ + η(x 1 v) 2, for v > x, v) =. (1 δ) 1 γ(x v) δ η(x 1 v)2 for v < x, Set x 1 = x, x 2 = y, x = x 3, y = x 4. Discontinuity surfaces Σ 1 and Σ 2 : h 1 (x) = x 3 v, h 2 (x) = x 4 v.
74 ... stick-slip... Fix m 1 = m 2 = 1, k 1 = k 2 = k 12 = 1, δ = 0, γ = 3, η = 0, v = 0.295, β = Time interval [0, 80], initial point at the origin x x x
75 ... stick-slip... Periodic solution x x x Can define a monodromy, do stability analysis, etc. etc..
Nonsmooth systems: synchronization, sliding and other open problems
John Hogan Bristol Centre for Applied Nonlinear Mathematics, University of Bristol, England Nonsmooth systems: synchronization, sliding and other open problems 2 Nonsmooth Systems 3 What is a nonsmooth
More informationDynamics of Non-Smooth Systems
Dynamics of Non-Smooth Systems P Chandramouli Indian Institute of Technology Madras October 5, 2013 Mouli, IIT Madras Non-Smooth System Dynamics 1/39 Introduction Filippov Systems Behaviour around hyper-surface
More informationDerivation of border-collision maps from limit cycle bifurcations
Derivation of border-collision maps from limit cycle bifurcations Alan Champneys Department of Engineering Mathematics, University of Bristol Mario di Bernardo, Chris Budd, Piotr Kowalczyk Gabor Licsko,...
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationSzalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing.
Szalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing. Early version, also known as pre-print Link to publication record in Explore Bristol Research PDF-document University
More informationThe reduction to invariant cones for nonsmooth systems. Tassilo Küpper Mathematical Institute University of Cologne
The reduction to invariant cones for nonsmooth systems Tassilo Küpper Mathematical Institute University of Cologne kuepper@math.uni-koeln.de International Workshop on Resonance Oscillations and Stability
More informationImplicit Euler Numerical Simulation of Sliding Mode Systems
Implicit Euler Numerical Simulation of Sliding Mode Systems Vincent Acary and Bernard Brogliato, INRIA Grenoble Rhône-Alpes Conference in honour of E. Hairer s 6th birthday Genève, 17-2 June 29 Objectives
More informationOne-Sided Event Location Techniques in the Numerical Solution of Discontinuous Differential Systems
BIT manuscript No. (will be inserted by the editor) One-Sided Event Location Techniques in the Numerical Solution of Discontinuous Differential Systems Luca Dieci Luciano Lopez Received: date / Accepted:
More informationLecture Note 7: Switching Stabilization via Control-Lyapunov Function
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 7: Switching Stabilization via Control-Lyapunov Function Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio
More informationSMOOTH TO DISCONTINUOUS SYSTEMS: A GEOMETRIC AND NUMERICAL METHOD FOR SLOW-FAST DYNAMICS. Luca Dieci. Cinzia Elia
DISCRETE AND CONTINUOUS doi:.3934/dcdsb.282 DYNAMICAL SYSTEMS SERIES B Volume 23, Number 7, September 28 pp. 2935 295 SMOOTH TO DISCONTINUOUS SYSTEMS: A GEOMETRIC AND NUMERICAL METHOD FOR SLOW-FAST DYNAMICS
More informationHouston Journal of Mathematics. c 1999 University of Houston Volume 25, No. 1, 1999
Houston Journal of Mathematics c 999 University of Houston Volume 5, No., 999 SLIDING MODES IN INTERSECTING SWITCHING SURFACES, II: HYSTERESIS JAMES C. ALEXANDER AND THOMAS I. SEIDMAN Communicated by Giles
More informationStability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach
Stability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach Kamyar Mehran School of Engineering University of Warwick Coventry, CV4 7LW, UK Email: k.mehran@warwick.ac.uk
More informationTHE GEOMETRY OF GENERIC SLIDING BIFURCATIONS
THE GEOMETRY OF GENERIC SLIDING BIFURCATIONS M. R. JEFFREY AND S. J. HOGAN Abstract. Using the singularity theory of scalar functions, we derive a classification of sliding bifurcations in piecewise-smooth
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
More informationModelling Research Group
Modelling Research Group The Importance of Noise in Dynamical Systems E. Staunton, P.T. Piiroinen November 20, 2015 Eoghan Staunton Modelling Research Group 2015 1 / 12 Introduction Historically mathematicians
More informationLecture Outline Chapter 6. Physics, 4 th Edition James S. Walker. Copyright 2010 Pearson Education, Inc.
Lecture Outline Chapter 6 Physics, 4 th Edition James S. Walker Chapter 6 Applications of Newton s Laws Units of Chapter 6 Frictional Forces Strings and Springs Translational Equilibrium Connected Objects
More informationOn rigorous integration of piece-wise linear systems
On rigorous integration of piece-wise linear systems Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology 06.2011, Trieste, Italy On rigorous integration of PWL
More informationSecond Order Optimality Conditions for Constrained Nonlinear Programming
Second Order Optimality Conditions for Constrained Nonlinear Programming Lecture 10, Continuous Optimisation Oxford University Computing Laboratory, HT 2006 Notes by Dr Raphael Hauser (hauser@comlab.ox.ac.uk)
More informationKinetics of Particles
Kinetics of Particles A- Force, Mass, and Acceleration Newton s Second Law of Motion: Kinetics is a branch of dynamics that deals with the relationship between the change in motion of a body and the forces
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More information03. Simple Dynamical Systems
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 03. Simple Dynamical Systems Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationWritten Examination
Division of Scientific Computing Department of Information Technology Uppsala University Optimization Written Examination 202-2-20 Time: 4:00-9:00 Allowed Tools: Pocket Calculator, one A4 paper with notes
More informationControlling Attractivity of Friction-Induced Equilibrium Sets
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 TuA16.5 Controlling Attractivity of Friction-Induced Equilibrium
More informationLagrange Multipliers
Lagrange Multipliers (Com S 477/577 Notes) Yan-Bin Jia Nov 9, 2017 1 Introduction We turn now to the study of minimization with constraints. More specifically, we will tackle the following problem: minimize
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationUniversity of Bath DOI: / Publication date: Link to publication SIAM
Citation for published version: Jeffrey, MR & Hogan, SJ 2011, 'The geometry of generic sliding bifurcations' Siam Review, vol. 53, no. 3, pp. 505-525. https://doi.org/10.1137/090764608 DOI: 10.1137/090764608
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationAttractors near grazing-sliding bifurcations. Glendinning, Paul and Kowalczyk, Piotr and Nordmark, Arne. MIMS EPrint:
Attractors near grazing-sliding bifurcations Glendinning, Paul and Kowalczyk, Piotr and Nordmark, Arne 11 MIMS EPrint: 11.63 Manchester Institute for Mathematical Sciences School of Mathematics The University
More informationSimulating Two-Dimensional Stick-Slip Motion of a Rigid Body using a New Friction Model
Proceedings of the 2 nd World Congress on Mechanical, Chemical, and Material Engineering (MCM'16) Budapest, Hungary August 22 23, 2016 Paper No. ICMIE 116 DOI: 10.11159/icmie16.116 Simulating Two-Dimensional
More informationPreview of Period 5: Forces and Newton s Laws
Preview of Period 5: Forces and Newton s Laws 5.1 The Fundamental Forces of Nature What are the four fundamental forces of nature? How do we see their effects? 5.2 Forces and Newton s Laws What causes
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Input-Output and Input-State Linearization Zero Dynamics of Nonlinear Systems Hanz Richter Mechanical Engineering Department Cleveland State University
More informationFirst Year Physics: Prelims CP1 Classical Mechanics: DR. Ghassan Yassin
First Year Physics: Prelims CP1 Classical Mechanics: DR. Ghassan Yassin MT 2007 Problems I The problems are divided into two sections: (A) Standard and (B) Harder. The topics are covered in lectures 1
More informationA Geometric Interpretation of Newtonian Impacts with Global Dissipation Index
A Geometric Interpretation of Newtonian Impacts with Global Dissipation Index Christoph Glocker Institute B of Mechanics, Technical University of Munich, D-85747 Garching, Germany (glocker@lbm.mw.tu-muenchen.de)
More informationActually, no perfectly frictionless surface exists. For two surfaces in contact, tangential forces, called friction forces, will develop if one
Chapter 8 Friction Actually, no perfectly frictionless surface exists. For two surfaces in contact, tangential forces, called friction forces, will develop if one attempts to move one relative to the other.
More informationTheoretical Statistical Physics
Janosh Riebesell, Adrian van Kan Lecturer: Manfred Salmhofer December nd, 06 Theoretical Statistical Physics Solution to Exercise Sheet 5 Ideal gas work (3 oints Within the kinetic model of an ideal gas,
More informationChapter 3a Topics in differentiation. Problems in differentiation. Problems in differentiation. LC Abueg: mathematical economics
Chapter 3a Topics in differentiation Lectures in Mathematical Economics L Cagandahan Abueg De La Salle University School of Economics Problems in differentiation Problems in differentiation Problem 1.
More informationModeling nonlinear systems using multiple piecewise linear equations
Nonlinear Analysis: Modelling and Control, 2010, Vol. 15, No. 4, 451 458 Modeling nonlinear systems using multiple piecewise linear equations G.K. Lowe, M.A. Zohdy Department of Electrical and Computer
More information1 Trajectory Generation
CS 685 notes, J. Košecká 1 Trajectory Generation The material for these notes has been adopted from: John J. Craig: Robotics: Mechanics and Control. This example assumes that we have a starting position
More informationFINITE-DIFFERENCE APPROXIMATIONS AND OPTIMAL CONTROL OF THE SWEEPING PROCESS. BORIS MORDUKHOVICH Wayne State University, USA
FINITE-DIFFERENCE APPROXIMATIONS AND OPTIMAL CONTROL OF THE SWEEPING PROCESS BORIS MORDUKHOVICH Wayne State University, USA International Workshop Optimization without Borders Tribute to Yurii Nesterov
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
PH 105 Exam 2 VERSION A Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) Is it possible for a system to have negative potential energy? A)
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part IV: Applications ISS Consider with solutions ϕ(t, x, w) ẋ(t) =
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
PH 105 Exam 2 VERSION B Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) A boy throws a rock with an initial velocity of 2.15 m/s at 30.0 above
More informationStability and Robustness of Homogeneous Differential Inclusions
Stability and Robustness of Homogeneous Differential Inclusions Arie Levant 1, Denis Efimov 23, Andrey Polyakov 23, Wilfrid Perruquetti 2 Abstract The known results on asymptotic stability of homogeneous
More informationReduction to invariant cones for non-smooth systems
Available online at www.sciencedirect.com Mathematics and Computers in Simulation 8 () 98 995 Reduction to invariant cones for non-smooth systems T. Küpper, H.A. Hosham Mathematisches Institut, Universität
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
More informationContractionanalysisofswitchedsystemsviaregularization
Contractionanalysisofswitchedsystemsviaregularization DavideFiore a,s.johnhogan b,mariodibernardo a,b a Department of Electrical Engineering and Information Technology, University of Naples Federico II,
More informationLecture 6 Friction. Friction Phenomena Types of Friction
Lecture 6 Friction Tangential forces generated between contacting surfaces are called friction forces and occur to some degree in the interaction between all real surfaces. whenever a tendency exists for
More informationONE-PARAMETER BIFURCATIONS IN PLANAR FILIPPOV SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 3, No. 8 (23) 257 288 c World Scientific Publishing Company ONE-PARAMETER BIFURCATIONS IN PLANAR FILIPPOV SYSTEMS YU. A. KUZNETSOV, S. RINALDI and A.
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationChapter 6 Applications of Newton s Laws. Copyright 2010 Pearson Education, Inc.
Chapter 6 Applications of Newton s Laws Units of Chapter 6 Frictional Forces Strings and Springs Translational Equilibrium Connected Objects Circular Motion 6-1 Frictional Forces Friction has its basis
More informationBifurcations in Nonsmooth Dynamical Systems
SIAM REVIEW Vol. 5, No. 4, pp. 629 7 c 28 Society for Industrial and Applied Mathematics Bifurcations in Nonsmooth Dynamical Systems Mario di Bernardo Chris J. Budd Alan R. Champneys Piotr Kowalczyk Arne
More informationComputation of Local ISS Lyapunov Function Via Linear Programming
Computation of Local ISS Lyapunov Function Via Linear Programming Huijuan Li Joint work with Robert Baier, Lars Grüne, Sigurdur F. Hafstein and Fabian Wirth Institut für Mathematik, Universität Würzburg
More informationGEO E1050 Finite Element Method Autumn Lecture. 9. Nonlinear Finite Element Method & Summary
GEO E1050 Finite Element Method Autumn 2016 Lecture. 9. Nonlinear Finite Element Method & Summary To learn today The lecture should give you overview of how non-linear problems in Finite Element Method
More informationarxiv: v1 [nlin.cd] 7 Jun 2010
arxiv:006.236v [nlin.cd] 7 Jun 200 Simple models of bouncing ball dynamics and their comparison Andrzej Okninski, Bogus law Radziszewski 2 Physics Division, Department of Mechatronics and Mechanical Engineering
More informationarxiv: v2 [cs.sy] 7 Apr 2017
arxiv:6.58v [cs.sy] 7 Apr 7 Observer design for piecewise smooth and switched systems via contraction theory Davide Fiore Marco Coraggio Mario di Bernardo, Department of Electrical Engineering and Information
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationLecture I: Basic Physics
1 Velocity: Instantaneous change in position! = $% ' $& Suppose object position ( ) and constant velocity!. After time step +: ( ) + + + = ( ) + +! + ( ) = ( ) + + + ( ) + =! +.! is never constant in practice
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationA sub-optimal second order sliding mode controller for systems with saturating actuators
28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 FrB2.5 A sub-optimal second order sliding mode for systems with saturating actuators Antonella Ferrara and Matteo
More informationChapter 10: Friction A gem cannot be polished without friction, nor an individual perfected without
Chapter 10: Friction 10-1 Chapter 10 Friction A gem cannot be polished without friction, nor an individual perfected without trials. Lucius Annaeus Seneca (4 BC - 65 AD) 10.1 Overview When two bodies are
More informationThe problem of singularity in impacting systems
The problem of singularity in impacting systems Soumitro Banerjee, Department of Physics Indian Institute of Science Education & Research, Kolkata, India The problem of singularity in impacting systems
More informationHermite Interpolation
Jim Lambers MAT 77 Fall Semester 010-11 Lecture Notes These notes correspond to Sections 4 and 5 in the text Hermite Interpolation Suppose that the interpolation points are perturbed so that two neighboring
More informationAXIAL STICK-SLIP LIMIT CYCLING IN DRILL-STRING DYNAMICS WITH DELAY
ENOC-2005, Eindhoven, Netherlands, 7-12 August 2005 AXIAL STICK-SLIP LIMIT CYCLING IN DRILL-STRING DYNAMICS WITH DELAY Christophe Germay Department of Electrical Engineering and Computer Science University
More informationPhysics 8, Fall 2011, equation sheet work in progress
1 year 3.16 10 7 s Physics 8, Fall 2011, equation sheet work in progress circumference of earth 40 10 6 m speed of light c = 2.9979 10 8 m/s mass of proton or neutron 1 amu ( atomic mass unit ) = 1 1.66
More informationWhen a rigid body is in equilibrium, both the resultant force and the resultant couple must be zero.
When a rigid body is in equilibrium, both the resultant force and the resultant couple must be zero. 0 0 0 0 k M j M i M M k R j R i R F R z y x z y x Forces and moments acting on a rigid body could be
More informationThe... of a particle is defined as its change in position in some time interval.
Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle
More informationCHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS. p. 1/73
CHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS p. 1/73 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS Mixed Inequality Constraints: Inequality constraints involving control and possibly
More informationJoint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018
EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,
More informationDEFINITION AND CLASSIFICATION OF FORCES
DEINITION AND CLASSIICATION O ORCES As defined before, force is an action of one body on another. It is a vector quantity since its effect depends on the direction as well as on the magnitude of the action.
More informationDynamical Analysis of a Hydraulic Pressure Relief Valve
Proceedings of the World Congress on Engineering 9 Vol II WCE 9, July -, 9, London, U.K. Dynamical Analysis of a Hydraulic Pressure Relief Valve Gábor Licskó Alan Champneys Csaba Hős Abstract A mathematical
More informationON THE ROBUSTNESS OF PERIODIC SOLUTIONS IN RELAY FEEDBACK SYSTEMS. Mario di Bernardo Λ Karl Henrik Johansson ΛΛ;1 Ulf Jönsson ΛΛΛ Francesco Vasca Λ
Copyright IFAC 1th Triennial World Congress, Barcelona, Spain ON THE ROBUSTNESS OF PERIODIC SOLUTIONS IN RELAY FEEDBACK SYSTEMS Mario di Bernardo Λ Karl Henrik Johansson ΛΛ;1 Ulf Jönsson ΛΛΛ Francesco
More informationLecture 8: Friction Part 1 Friction Phenomena
Lecture 8: Friction Part 1 Friction Phenomena Friction is a force which is generated between. Its direction is always to that of the motion, or tendency for a motion. F mg P N Figure 1 Figure 1 shows a
More informationarxiv: v1 [math.ds] 20 Sep 2018
ON NON SMOOTH SLOW FAST SYSTEMS JAIME R. DE MORAES 1 AND PAULO R. DA SILVA 2 arxiv:1809.07612v1 [math.ds] 20 Sep 2018 Abstract. We deal with non smooth differential systems ż = X(z), z R n, with discontinuity
More informationEveryday Forces. MCHS Honors Physics Weight. Weight
Everyday Forces Weight The Force of Friction Kinetic vs. Static The Coefficient of Friction MCHS Honors Physics 2014-15 Weight The gravitational force exerted on an object by Earth, F g, is a vector quantity,
More informationNonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control
Nonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control 19/4/2012 Lecture content Problem formulation and sample examples (ch 13.1) Theoretical background Graphical
More informationCALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD
Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen
More informationOutline: Types of Friction Dry Friction Static vs. Kinetic Angles Applications of Friction. ENGR 1205 Appendix B
Outline: Types of Friction Dry Friction Static vs. Kinetic Angles Applications of Friction ENGR 1205 Appendix B 1 Contacting surfaces typically support normal and tangential forces Friction is a tangential
More informationSECTION C: CONTINUOUS OPTIMISATION LECTURE 9: FIRST ORDER OPTIMALITY CONDITIONS FOR CONSTRAINED NONLINEAR PROGRAMMING
Nf SECTION C: CONTINUOUS OPTIMISATION LECTURE 9: FIRST ORDER OPTIMALITY CONDITIONS FOR CONSTRAINED NONLINEAR PROGRAMMING f(x R m g HONOUR SCHOOL OF MATHEMATICS, OXFORD UNIVERSITY HILARY TERM 5, DR RAPHAEL
More informationThe Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction
Send Orders for Reprints to reprints@benthamscience.ae 308 The Open Mechanical Engineering Journal, 2014, 8, 308-313 Open Access The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationPhysics 331 Introduction to Numerical Techniques in Physics
Physics 331 Introduction to Numerical Techniques in Physics Instructor: Joaquín Drut Lecture 12 Last time: Polynomial interpolation: basics; Lagrange interpolation. Today: Quick review. Formal properties.
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationPIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES. 1. Introduction
PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES JOAN C. ARTÉS1, JAUME LLIBRE 1, JOAO C. MEDRADO 2 AND MARCO A. TEIXEIRA 3 Abstract. In this paper we study piecewise linear differential systems
More informationThe Two-Fold Singularity of Discontinuous Vector Fields
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 8, No. 2, pp. 624 640 c 2009 Society for Industrial and Applied Mathematics The Two-Fold Singularity of Discontinuous Vector Fields M. R. Jeffrey and A. Colombo Abstract.
More informationARCS IN FINITE PROJECTIVE SPACES. Basic objects and definitions
ARCS IN FINITE PROJECTIVE SPACES SIMEON BALL Abstract. These notes are an outline of a course on arcs given at the Finite Geometry Summer School, University of Sussex, June 26-30, 2017. Let K denote an
More informationKrylov Subspace Methods for Nonlinear Model Reduction
MAX PLANCK INSTITUT Conference in honour of Nancy Nichols 70th birthday Reading, 2 3 July 2012 Krylov Subspace Methods for Nonlinear Model Reduction Peter Benner and Tobias Breiten Max Planck Institute
More informationForce a push or a pull exerted on some object the cause of an acceleration, or the change in an objects velocity
Chapter 4 Physics Notes Changes in Motion Force a push or a pull exerted on some object the cause of an acceleration, or the change in an objects velocity Forces cause changes in velocity Causes a stationary
More information41514 Dynamics of Machinery
41514 Dynamics of Machinery Theory, Experiment, Phenomenology and Industrial Applications Ilmar Ferreira Santos 1. Recapitulation Mathematical Modeling & Steps 2. Example System of Particle 3. Example
More informationThe dynamics of a nonautonomous oscillator with friction memory
The dynamics of a nonautonomous oscillator with friction memory L. A. Igumnov 1, V. S. Metrikin 2, M. V. Zaytzev 3 1, 2 Research Institute for mechanics, National Research Lobachevsky State University
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationLecture 6 Physics 106 Spring 2006
Lecture 6 Physics 106 Spring 2006 Angular Momentum Rolling Angular Momentum: Definition: Angular Momentum for rotation System of particles: Torque: l = r m v sinφ l = I ω [kg m 2 /s] http://web.njit.edu/~sirenko/
More informationGeneralized mid-point and trapezoidal rules for Lipschitzian RHS
Generalized mid-point and trapezoidal rules for Lipschitzian RHS Richard Hasenfelder, Andreas Griewank, Tom Streubel Humboldt Universität zu Berlin 7th International Conference on Algorithmic Differentiation
More informationDry Friction Static vs. Kinetic Angles
Outline: Types of Friction Dry Friction Static vs. Kinetic Angles Applications of Friction 1 Contacting surfaces typically support normal and tangential forces Friction is a tangential force Friction occurs
More informationSecond Order Sliding Mode Control for Nonlinear Affine Systems with Quantized Uncertainty
Second Order Sliding Mode Control for Nonlinear Affine Systems with Quantized Uncertainty Gian Paolo Incremona a, Michele Cucuzzella b, Antonella Ferrara b a Dipartimento di Elettronica, Informazione e
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationChapter 6 Dynamics I: Motion Along a Line
Chapter 6 Dynamics I: Motion Along a Line Chapter Goal: To learn how to solve linear force-and-motion problems. Slide 6-2 Chapter 6 Preview Slide 6-3 Chapter 6 Preview Slide 6-4 Chapter 6 Preview Slide
More informationChapter 3 HIGHER ORDER SLIDING MODES. L. FRIDMAN and A. LEVANT. Chihuahua Institute of Technology, Chihuahua, Mexico.
Chapter 3 HIGHER ORDER SLIDING MODES L. FRIDMAN and A. LEVANT Chihuahua Institute of Technology, Chihuahua, Mexico. Institute for Industrial Mathematics, Beer-Sheva, Israel. 3.1 Introduction One of the
More information