You will be provided with a reference and some statements. Please determine whether each statement is 'supported', 'unsupported', or 'unknown' with respect to the reference. Please note:
First, assess whether the reference contains any valid content. If the reference contains no valid information, such as a 'page not found' message, then all statements should be considered 'unknown'.
If the reference is valid, for a given statement: if the facts or data it contains can be found entirely or partially within the reference, it is considered 'supported' (data accepts rounding); if all facts and data in the statement cannot be found in the reference, it is considered 'unsupported'.

You should return the result in a JSON list format, where each item in the list contains the statement's index and the judgment result, for example:
[
    {
        "idx": 1,
        "result": "supported"
    },
    {
        "idx": 2,
        "result": "unsupported"
    }
]

Below are the reference and statements:
<reference>
Specified-duration shapers for suppressing residual vibrations | PLOS One

Input-shaping control has received considerable research attention for suppressing residual vibrations. Although numerous studies have been conducted on designing input shapers with arbitrary robustness to modeling errors, no studies have focused on the design of input shapers with arbitrarily specified shaping times. In this study, a specified-duration (SD) shaper, which is an input shaper with an arbitrarily specified shaping time, and a systematic method to design an SD shaper using impulse vectors are proposed. As the specified shaping time increases, the SD shaper increases the number of impulses one by one according to the number of added derivative constraints, thereby improving robustness to modeling errors. The performance of the SD shaper was evaluated for a second-order system through computer simulations. The simulation results revealed that the SD shaper suppresses residual vibrations of the vibratory system at the specified shaping time. The validity of the SD shaper was experimentally verified using a horizontal beam vibration apparatus. The results of this study provide insight into the development of vibration suppression strategies with input shaping control.

Skip to main content

Advertisement

PLOS One

Publish

Submissions

Getting Started

Submission Guidelines

Figures

Tables

Supporting Information

LaTeX

What We Publish

Preprints

Revising Your Manuscript

Submit Now

Calls for Papers

Policies

Best Practices in Research Reporting

Human Subjects Research

Animal Research

Competing Interests

Disclosure of Funding Sources

Licenses and Copyright

Data Availability

Complementary Research

Materials, Software and Code Sharing

Ethical Publishing Practice

Authorship

Corrections, Expressions of Concern, and Retractions

Manuscript Review and Publication

Criteria for Publication

Editorial and Peer Review Process

Editor Center

Guidelines for Reviewers

Accepted Manuscripts

Comments

Submit Your Manuscript

Discover a faster, simpler path to publishing in a high-quality journal.
PLOS ONE
promises fair, rigorous peer review,
broad scope, and wide readership – a perfect fit for your research every time.

Learn More

Submit Now

About

Why Publish with PLOS ONE

Journal Information

Staff Editors

Editorial Board

Section Editors

Advisory Groups

Find and Read Articles

Publishing Information

Publication Fees

Press and Media

Contact

Browse

Search

Search

advanced search

Browse Subject Areas

?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click

here
.

Loading metrics

Open Access

Peer-reviewed

Research Article

Specified-duration shapers for suppressing residual vibrations

Brian Byunghyun Kang

Roles

Writing – original draft,

Writing – review & editing

* E-mail:

brianbkang@sejong.ac.kr

Affiliation

School of Intelligent Mechatronics Engineering, Sejong University, Seoul, Korea

https://orcid.org/0000-0002-1370-2867

⨯

Specified-duration shapers for suppressing residual vibrations

Brian Byunghyun Kang

x

Published: November 28, 2022

https://doi.org/10.1371/journal.pone.0276669

Article

Authors

Metrics

Comments

Media Coverage

Peer Review

Reader Comments

Figures

Figures

Abstract
Input-shaping control has received considerable research attention for suppressing residual vibrations. Although numerous studies have been conducted on designing input shapers with arbitrary robustness to modeling errors, no studies have focused on the design of input shapers with arbitrarily specified shaping times. In this study, a specified-duration (SD) shaper, which is an input shaper with an arbitrarily specified shaping time, and a systematic method to design an SD shaper using impulse vectors are proposed. As the specified shaping time increases, the SD shaper increases the number of impulses one by one according to the number of added derivative constraints, thereby improving robustness to modeling errors. The performance of the SD shaper was evaluated for a second-order system through computer simulations. The simulation results revealed that the SD shaper suppresses residual vibrations of the vibratory system at the specified shaping time. The validity of the SD shaper was experimentally verified using a horizontal beam vibration apparatus. The results of this study provide insight into the development of vibration suppression strategies with input shaping control.

Citation:
Kang BB (2022) Specified-duration shapers for suppressing residual vibrations. PLoS ONE 17(11):
e0276669.

https://doi.org/10.1371/journal.pone.0276669
Editor:
Ning Cai, Beijing University of Posts and Telecommunications, CHINA
Received:
June 9, 2022;
Accepted:
October 12, 2022;
Published:
November 28, 2022
Copyright:
© 2022 Brian Byunghyun Kang. This is an open access article distributed under the terms of the
Creative Commons Attribution License
, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability:
https://github.com/BBKangIRL/SDshaperExp.git
.
Funding:
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2021R1G1A109521912) and the faculty research fund of Sejong University in 2021.
Competing interests:
The authors have declared that no competing interests exist.

Introduction
In engineering systems, it is often beneficial to be able to rapidly reduce vibrations to avoid system performance degradation. For example, the forced lateral vibration of a building during earthquakes should be minimized to avoid structural damage [
1
]. Another example is that when cranes [
2
–
4
], disk drives [
5
,
6
], sloshing liquid containers [
7
–
10
], industrial robots [
11
], and soft robots [
12
,
13
] are operated, residual vibrations should be mitigated swiftly to reduce working hours and increase working efficiency.

Input-shaping control has attracted considerable research attention for removing residual vibration. Input-shaping control was first introduced as posicast control in the 1950s by O. J. M. Smith [
14
]. Since the 1980s, N. Singer, J. Hyde, W. Singhose, and T. Singh have considerably improved input-shaping control by adding derivative constraints [
15
–
19
]. This control exhibits considerable residual vibration removal efficiency without an additional feedback loop in a vibratory system.

The block diagram of a general input-shaping control system is displayed in
Fig 1
. The vibratory system in
Fig 1
includes a feedback loop and a feedback controller. Input-shaping control does not affect the stability of the vibratory system because the input shaper is connected in series outside the feedback loop as shown in
Fig 1
. The input shaper consists of a sequence of impulse functions. The convolution of this impulse sequence (i.e., input shaper) with the command input of the system can remove residual vibration caused by the command input [
20
].

Download:

PNG
larger image
TIFF
original image
Fig 1.
Block diagram of an input-shaping control system.
The input shaper is connected in series with a vibratory system and located outside feedback loops. Thus, the shaper does not affect the stability of the vibratory system.

https://doi.org/10.1371/journal.pone.0276669.g001
Kang [
21
] has recently introduced an impulse vector as a mathematical tool for designing an input shaper. Residual vibrations can be suppressed regardless of the damping condition by designing input shapers using impulse vectors such that the sum of the impulse vectors in the polar coordinate plane becomes zero.

In the 1990s, specified-insensitivity (SI) shapers were developed to be robust to modeling errors with arbitrarily specified insensitivity [
22
–
25
]. The SI shaper was designed to ensure that the sensitivity function
V
of the input shaper was within a tolerance value
V
tol
for a specified range of modeling errors. The sensitivity function
V
indicates the relative magnitude of residual vibration for modeling errors in natural frequency and damping ratio. Thus, to design an SI shaper with, for example,
V
tol
= 0.05, we first specify an arbitrary 5% insensitivity
I
0.05
, which forms the range of modeling errors such that the relative residual vibration is within 5%, and then design the input shaper to satisfy this
I
0.05
condition. The 5% insensitivity
I
0.05
represents the robustness to modeling errors of the input shaper when
V
tol
= 0.05.

The shaping time (or shaper duration) as well as the robustness to modeling errors are important factors in the design of input shapers. Because the shaping time represents the time to remove residual vibrations, smaller shaping times are preferred. As the robustness to modeling errors (e.g., 5% insensitivity) increases, the shaping time also increases nonlinearly in general [
26
]. Moreover, the shaping time can be reduced to a certain limit because the smaller the shaping time is, the greater the actuator power is required. Although the shaping time is an important factor in the design of input shapers [
27
], no studies have been conducted on designing an input shaper that can satisfy an arbitrarily specified shaping time.

In this paper, a novel specified-duration (SD) shaper that satisfies an arbitrarily specified shaping time with the largest 5% insensitivity to modeling errors when the magnitude of the last impulse is taken as a free parameter is proposed. The SD shaper is an input shaper that is designed to remove residual vibration during a given specified shaping time. As the specified shaping time increases, the SD shaper increases the number of impulses one by one according to the number of added derivative constraints, thereby improving robustness to modeling errors. A systematic method is presented to design the SD shaper with an arbitrarily specified shaping time by using impulse vectors. The performance of the proposed SD shaper is evaluated through computer simulations. Furthermore, the performance of the SD shaper is verified experimentally using a horizontal beam vibration apparatus.

This paper is organized as follows. Section 2 presents a systematic method to design an SD shaper with an arbitrarily shaping time by using impulse vectors. In Section 3, the performance of the SD shaper is discussed via MATLAB simulations using sensitivity curves and step responses. Section 4 demonstrates the validity of the SD shaper experimentally. The conclusion of the paper is presented in Section 5.

Design of an SD shaper
In this section, the design of an SD shaper with an arbitrarily specified shaping time
t
s
is described. Vibratory system dynamics is assumed to be an undamped or underdamped second-order system
, where
ω
n
is the undamped natural frequency and
ζ
is the damping ratio of the system. An impulse vector
I
i
(
i
= 1, 2, ⋯,
N
) with magnitude
I
i
and angle
θ
i
in the polar coordinate plane was proposed by Kang [
21
], which was as follows:

(1)

(2)

where
A
i
is the impulse magnitude of the
i
th
impulse function,
t
i
is the impulse time of the
i
th
impulse function, and
ω
d
is the damped natural frequency of the system, which is defined as

(3)

We can design an SD shaper with an arbitrarily specified shaping time
t
s
by using the impulse vector. The shaping time
t
s
is equal to the last impulse time
t
N
of the impulse sequence in which
N
represents the number of impulses of the input shaper. First, for an arbitrarily specified shaping time
t
s
=
t
N
, the angle
θ
N
of the last impulse vector
I
N
is determined using the impulse vector definition,
Eq (2)
, and this angle is fixed for the SD shaper design.

(4)
Thus, the given shaping time
t
s
determines the angle
θ
N
of
I
N
, and the magnitude
I
N
of
I
N
determines
A
N
. Here,
I
N
is a value between 0 and
because 0 <
A
N
< 1.

The damped period
T
d
of the second-order system with a damping ratio
ζ
is expressed as follows:

(5)

The dimensionless shaping time
T
s
is defined as
t
s
divided by the damped period
T
d
.

(6)
Then, the actual shaping time
t
s
can be expressed as
t
s
=
T
d
T
s
and the angle of the final impulse vector
I
N
of the SD shaper can be expressed as
θ
N
=
ω
d
T
d
T
s
.

It is well known that the ZV, ZVD and ZVD
n
shapers have the shaping times of a half period, one period, and (
n
+ 1)/2 period, which are
T
s
= 0.5,
T
s
= 1, and
T
s
= (
n
+ 1)/2, respectively [
21
]. Also, it is well known that the residual vibration is removed if the resultant of the impulse vectors in an impulse vector diagram is zero [
21
]. Moreover, we need some flexibility in order to achieve robustness to modeling errors of the SD shaper, which can be accomplished by adding an impulse and letting its magnitude serve as a free parameter. Observing these facts, the rule of deciding the number of impulses of SD shapers are determined as follows.

If the dimensionless shaping time
T
s
is between 0.5 and 1 (that is, between half and a period), the SD shaper consists of three impulses. If
T
s
is between 1 and 1.5, the SD shaper consists of four impulses, and if it is between 1.5 and 2, it consists of five impulses, and so on. Thus, we have the following expression:

(7)

Eq (7)
can be explained as follows. If
T
s
is between 0.5 and 1, the angle of the final impulse vector is in
π
<
θ
N
≤ 2
π
, and thus, an impulse vector diagram with zero resultant can be drawn using three impulse vectors, as shown in
Fig 2
. Therefore, three impulses are required. Three impulses have four unknowns
A
1
,
A
2
,
A
3
,
t
2
since
t
3
is given from the specified shaping time. These four unknowns can be obtained by simultaneously solving two residual vibration constraints and one normalization constraint with
A
3
chosen in advance by the designer. Without losing generality, we set
t
1
= 0. Further, if
T
s
is between 1 and 1.5, the angle of the final impulse vector is in 2
π
<
θ
N
≤ 3
π
, and an impulse vector diagram with zero resultant can be drawn using four impulse vectors, as shown in
Fig 3
, which corresponds to an input shaper with four impulses. The four impulses can be obtained by simultaneously solving two residual vibration constraints, two derivative constraints, and a normalization constraint (six unknowns
A
1
,
A
2
,
A
3
,
A
4
,
t
2
,
t
3
with
A
4
chosen in advance by the designer). Similarly, if
T
s
is between 1.5 and 2,
Eq (7)
states that we should use five impulses; these five impulses can be obtained by simultaneously solving two residual vibration constraints, two derivative constraints, two second-order derivative constraints, and a normalization constraint (eight unknowns
A
1
,
A
2
,
A
3
,
A
4
,
A
5
,
t
2
,
t
3
,
t
4
, with
A
5
chosen in advance by the designer) [
21
–
24
]. For an example that provides a partial justification for using five impulses in this case, refer to the last paragraph of the next section, Analysis of SD shapers.

Download:

PNG
larger image
TIFF
original image
Fig 2.
Impulse vector diagram of an SD shaper with a shaping time 0.5 <
T
s
≤ 1.
Because the shaping time
T
s
is in 0.5 <
T
s
≤ 1, the SD shaper is composed of three impulses.

https://doi.org/10.1371/journal.pone.0276669.g002

Download:

PNG
larger image
TIFF
original image
Fig 3.
Impulse vector diagram of an SD shaper with a shaping time 1 <
T
s
≤ 1.5.
Because the shaping time
T
s
is between 1 <
T
s
≤ 1.5, the SD shaper is composed of four impulses.

https://doi.org/10.1371/journal.pone.0276669.g003
Using the final impulse magnitude
A
N
as a free parameter to minimize the sensitivity to modeling errors, we can maximize the 5% insensitivity by increasing
A
N
discretely by 0.01 from 0.01 to 0.99 in the design of SD shapers.

The design procedure of SD shapers can be summarized as follows:

Determine the number of impulses
N
from the specified shaping time using
Eq (7)
.

Obtain the constraint equations for the SD shaper with
N
impulses.

Find impulse magnitudes and times
A
i
,
t
i
by solving the constraint equations analytically or numerically, in which the last impulse magnitude
A
N
is a free parameter (discrete values with 0.01 intervals) to maximize the 5% insensitivity to modeling errors.

Here, the design of an SD shaper with 0.5 <
T
s
≤ 1 is described. If
T
s
is between 0.5 and 1, we obtain
N
= 3, as determined by
Eq (7)
, and the impulse vector diagram of the SD shaper can be drawn as shown in
Fig 2
. Here,
θ
3
is determined from the given shaping time
t
3
(=
t
s
), and
θ
2
,
I
1
,
I
2
,
I
3
are unknowns. To ensure residual vibration
V
is zero when there is no modeling error, the sum of impulse vectors must be zero. Therefore, the SD shaper for 0.5 <
T
s
≤ 1 should satisfy the following constraints:

(8)

(9)

(10)

Eqs (
8
)–(
10
) have three unknowns, namely
A
1
,
A
2
and
t
2
because
θ
3
,
I
3
are given.
θ
3
is obtained from the specified shaping time, and
I
3
is given from the free parameter
A
3
that is used for maximizing the 5% insensitivity to modeling errors. Thus, Eqs (
8
)–(
10
) have a solution. Here,
I
1
=
A
1
,
,
and
θ
2
=
ω
d
t
2
. For a given
θ
3
, we assume a specific value of the final impulse magnitude
A
3
(0 <
A
3
< 1), and find
A
1
,
A
2
,
t
2
that satisfy Eqs (
8
)–(
10
). Then, we vary the value of
A
3
, and find
A
1
,
A
2
,
t
2
that satisfy Eqs (
8
)–(
10
), repeatedly.

SD shapers are obtained separately for two cases with
ζ
= 0 and 0 <
ζ
< 1.

In case of
ζ
= 0, Eqs (
8
) and (
9
) can be expressed as follows:

(11)

(12)

By rearranging Eqs (
10
) and (
11
), we obtain the following expression:

(13)

and, from
Eq (12)
, we have the following expression:

(14)

Dividing
Eq (14)
by
Eq (13)
results in the following expression:

(15)

By selecting a specific value of
A
3
in
Eq (15)
,
b
is determined, and
Eq (15)
can be expressed as follows:

(16)

Eq (16)
is a quadratic equation for cos
θ
2
. Because cos
θ
2
< 1, the solution of
Eq (16)
is expressed as follows:

(17)

Thus, the angle
θ
2
of impulse vector
I
2
is obtained as follows:

(18)

When
A
3
is specified,
θ
2
is determined by
Eq (18)
. Because
θ
3
is given, the impulse magnitudes
A
1
,
A
2
can be obtained as follows:

(19)

From Eqs (
18
) and (
19
),
A
1
,
A
2
,
θ
2
can be calculated repeatedly by increasing
A
3
by 0.01 from 0.01 to 0.99, and, for each case, an input shaper
A
1
δ
(
t
) +
A
2
δ
(
t
−
t
2
) +
A
3
δ
(
t
−
t
3
) can be obtained, in which
δ
(
t
) implies Dirac delta function. For each input shaper
A
1
δ
(
t
) +
A
2
δ
(
t
−
t
2
) +
A
3
δ
(
t
−
t
3
), a sensitivity curve is drawn and 5% insensitivity
I
0.05
is obtained. Among these sensitivity curves, the input shaper
A
1
δ
(
t
) +
A
2
δ
(
t
−
t
2
) +
A
3
δ
(
t
−
t
3
) with the largest
I
0.05
is selected as the SD shaper. The iterative procedure of finding numerical solutions can be performed using numeric computing platforms (e.g., MATLAB).

The sensitivity curve is a plot of the sensitivity function
V
with respect to modeling errors
, in which
represents the modeled value of natural frequency, and
ω
n
is the actual value of natural frequency. In the sensitivity curves,
is assumed to be fixed, and
ω
n
is a varying parameter. The sensitivity function
V
represents the relative magnitude of residual vibration and is defined as follows [
21
]:

(20)

where

(21)

In Eqs (
20
) and (
21
),
ω
n
,
ω
d
,
ζ
are the varying actual values of the system, whereas
t
i
,
A
i
are the constant values obtained from the fixed modeled values
. The 5% insensitivity
I
0.05
is the distance between two points where
V
= 0.05 line and the sensitivity curve intersect.

As an example, let us determine an SD shaper with a shaping time
t
s
= 0.3s for an undamped system
with a natural frequency
ω
n
of 2 Hz (= 4
π
rad/s). In this case, the dimensionless shaping time is
T
s
= 0.6 according to
Eq (6)
. Thus, the SD shaper is composed of three impulses according to
Eq (7)
. The SD shaper can be obtained using Eqs (
1
), (
2
), (
15
), (
18
) and (
19
) as follows:

(22)

In this SD shaper, the angle
θ
3
of the final impulse vector is 1.2π rad by
Eq (4)
.

Next, an SD shaper with 0.5 <
T
s
≤ 1 is obtained for an underdamped system with a damping ratio 0 <
ζ
< 1. If 0 <
ζ
< 1, we cannot determine a closed-form solution such as
Eq (19)
for given
A
3
because Eqs (
8
) and (
9
) include exponential functions. However, a numerical solution for Eqs (
8
)–(
10
) can be obtained because nonlinear simultaneous equations can be solved numerically along with the impulse vector definitions using the fsolve() function in MATLAB. Next, by increasing
A
3
by 0.01 from 0.01 to 0.99, a numerical solution can be determined for
A
1
,
A
2
and
θ
2
. Then, a sensitivity curve can be drawn for each input shaper. Among these sensitivity curves, the input shaper
A
1
δ
(
t
) +
A
2
δ
(
t
−
t
2
) +
A
3
δ
(
t
−
t
3
) that maximizes 5% insensitivity
I
0.05
is selected as the SD shaper for the system with 0 <
ζ
< 1.

As an example, let us determine an SD shaper with a shaping time
t
s
= 0.3s for an underdamped system
with a natural frequency
ω
n
of 2 Hz and a damping ratio
ζ
of 0.1. In this case, the dimensionless shaping time is
T
s
= 0.597 according to
Eq (6)
. Thus, the SD shaper is composed of three impulses according to
Eq (7)
. The SD shaper obtained numerically using MATLAB is expressed as follows:

(23)

From
Eq (4)
, the angle
θ
3
of the final impulse vector in this SD shaper is 1.194π rad.

Next, we obtain an SD shaper for the case when
T
s
is between 1 and 1.5. Even when
T
s
is between 1 and 1.5, an input shaper with three impulses can be obtained from constraints Eqs (
8
)–(
10
). However, to maximize 5% insensitivity, we obtain an SD shaper with four impulses as displayed in
Fig 3
by adding derivative constraints. Therefore, when
T
s
is greater than 1 and less than or equal to 1.5, the SD shaper with
N
= 4 should satisfy the following constraints:

(24)

Two derivative constraints, the third and fourth equations of
Eq (24)
, are obtained by differentiating two residual vibration constraints, the first and second equations of
Eq (24)
, with respect to
ω
n
. The first and second constraints of
Eq (24)
can be expressed as functions of
ω
n
as follows:

Then, we can obtain two equations,
dC
/
dω
n
= 0 and
dS
/
dω
n
= 0, by differentiating
C
(
ω
n
) and
S
(
ω
n
) with respect to
ω
n
. These two equations are equivalent to the third and fourth constraints of
Eq (24)
, respectively.

Regardless of whether
ζ
is 0, or greater than 0 in this case, by increasing
A
4
by 0.01 from 0.01 to 0.99, we can determine
A
1
,
A
2
,
A
3
,
t
2
, and
t
3
numerically from
Eq (24)
and impulse vector definitions (1), (2). The sensitivity curve is then drawn for each input shaper with
A
1
,
A
2
,
A
3
,
t
2
,
t
3
. The input shaper
A
1
δ
(
t
) +
A
2
δ
(
t
−
t
2
) +
A
3
δ
(
t
−
t
3
) +
A
4
δ
(
t
−
t
4
) that maximizes 5% insensitivity, is selected as the SD shaper for 1 <
T
s
≤ 1.5.

As an example, let us determine an SD shaper with a shaping time
t
s
= 0.6 s for the aforementioned underdamped system
with a natural frequency
ω
n
of 2 Hz, and a damping ratio
ζ
of 0.1. The dimensionless shaping time is
T
s
= 1.194 according to
Eq (6)
. Thus, the SD shaper is composed of four impulses according to
Eq (7)
. The SD shaper can be obtained numerically using MATLAB as follows:

(25)

Next, we obtain an SD shaper for the case when
T
s
is between 1.5 and 2. When
T
s
is between 1.5 and 2, to maximize the 5% insensitivity, we obtain the SD shaper with five impulses by adding second-order derivative constraints. Therefore, when
T
s
is greater than 1.5 and less than or equal to 2, the SD shaper with
N
= 5 should satisfy the following constraints:

(26)

Two second-order derivative constraints, the fifth and sixth equations of
Eq (26)
, are obtained by differentiating two derivative constraints, the third and fourth equations of
Eq (26)
, with respect to
ω
n
. Regardless of whether
ζ
is 0 or greater than 0, by increasing
A
4
by 0.01 from 0.01 to 0.99, we can determine
A
1
,
A
2
,
A
3
,
A
4
,
t
2
,
t
3
and
t
4
numerically from
Eq (26)
and impulse vector definitions (1), (2). The sensitivity curve is then drawn for each input shaper with
A
1
,
A
2
,
A
3
,
A
4
,
t
2
,
t
3
and
t
4
. The input shaper
A
1
δ
(
t
) +
A
2
δ
(
t
−
t
2
) +
A
3
δ
(
t
−
t
3
) +
A
4
δ
(
t
−
t
4
) +
A
5
δ
(
t
−
t
5
) that maximizes 5% insensitivity is selected as the SD shaper for 1.5 <
T
s
≤ 2.

As an example, let us determine an SD shaper with a shaping time
t
s
= 0.85 s for the aforementioned underdamped system
with a natural frequency
ω
n
of 2 Hz, and a damping ratio
ζ
of 0.1. The dimensionless shaping time is
T
s
= 1.6915 according to
Eq (6)
. Thus, the SD shaper is composed of five impulses according to
Eq (7)
. The SD shaper with
t
s
= 0.85 s can be obtained numerically using MATLAB as follows:

(27)

When the dimensionless shaping time
T
s
is greater than 2, the SD shaper can be designed by repeating the aforementioned numerical analysis procedure by increasing the number of impulses by one. However, if the shaping time of the SD shaper is longer than two natural periods, the usefulness of the SD shaper decreases because the time for removing residual vibration increases correspondingly. Therefore, SD shapers with more than two natural periods were not considered. Moreover, if the dimensionless shaping time is less than 0.5, we require negative impulses to remove residual vibrations. Thus, SD shapers with
T
s
< 0.5 were not considered in this study.

Note that an SD shaper satisfying all nonlinear constraints always exists because an impulse vector diagram satisfying all nonlinear constraints of the SD shaper can always be drawn geometrically; then, the SD shaper can be obtained from the impulse vector diagram.

Analysis of SD shapers
The performance of the SD shaper obtained in the previous section is evaluated through sensitivity curves and unit-step responses. For the undamped system
modeled with a natural frequency of 2 Hz (i.e.,
), the SD shaper with a shaping time
t
s
= 0.3s is expressed by
Eq (22)
. The sensitivity curve of the SD shaper (22) is shown in
Fig 4
. For comparison, the sensitivity curves of the well-known ZV (zero vibration) shaper and ZVD (zero vibration and derivative) shaper [
18
] are also shown in
Fig 4
. The robustness of the SD shaper to modeling errors shown in
Fig 4
is between the ZV and ZVD shapers if the actual natural frequency is less than approximately 150% of the modeled value; however, the robustness of the SD shaper to modeling errors is better than that of the ZV and ZVD shapers if the actual natural frequency is 150% larger than the modeled value. The 5% insensitivity of the SD shaper is 0.073, which is between those of the ZV shaper (
I
0.05
= 0.063) and ZVD shaper (
I
0.05
= 0.287). The sensitivity curves of ZV and ZVD shapers for undamped systems are symmetric about the axis
, but the sensitivity curves of SD shapers are asymmetric even for undamped systems, as shown in
Fig 4
. In
Fig 4
,
is a fixed modeled value, and
ω
n
is a varying actual value.

Download:

PNG
larger image
TIFF
original image
Fig 4.
Sensitivity curve of the SD shaper with the shaping time
t
s
= 0.3s for an undamped system
with modeled natural frequency
.
Because the dimensionless shaping time is 0.6, the SD shaper is composed of three impulses. The 5% insensitivity
I
0.05
of this SD shaper is 0.073.

https://doi.org/10.1371/journal.pone.0276669.g004
Fig 5
shows the unit-step responses of the input-shaping control system with an SD shaper when
ζ
= 0.
Fig 5A
shows the unit-step responses when no modeling error (
) exists, and
Fig 5B
shows the unit-step responses when 20% modeling error exists in natural frequency (
ω
n
= 2.4 Hz,
). If no modeling error exists, the SD shaper has a shaping time of exactly 0.3 s as designed in
Eq (22)
. If 20% modeling error exists in natural frequency, the magnitude of the residual vibration is between those of the ZV and ZVD shaper as shown in
Fig 5B
, which is consistent with the sensitivity curve trend in
Fig 4
.

Download:

PNG
larger image
TIFF
original image
Fig 5.
Unit-step responses when the SD shaper with
t
s
= 0.3s is used for an undamped system
with
.
(A) No modeling error (
), (B) 20% modeling error in natural frequency (
ω
n
= 2.4 Hz,
).

https://doi.org/10.1371/journal.pone.0276669.g005
For the underdamped system
modeled with a natural frequency of 2 Hz and a damping ratio of 0.1, SD shapers with a shaping time of 0.3, 0.6, and 0.85 s are obtained as Eqs (
23
), (
25
) and (
27
), with three, four, and five impulses, respectively. The sensitivity curves for these SD shapers are shown in
Fig 6
, in which the robustness of SD shapers to modeling errors improves as the shaping time increases. The 5% insensitivities
I
0.05
of these SD shapers calculated using MATLAB were 0.088, 0.452, and 1.133, respectively.
Fig 6
shows the sensitivity curves of the ZV and ZVD shapers for comparison. In case of a shaping time of 0.3 s,
I
0.05
of the undamped system is 0.073, whereas
I
0.05
of the underdamped system with
ζ
= 0.1 is increased to 0.088. The
I
0.05
increased because the residual vibration
V
decreased owing to increased damping.

Download:

PNG
larger image
TIFF
original image
Fig 6.
Sensitivity curves of the SD shapers with shaping time
t
s
= 0.3, 0.6, 0.85s for an underdamped system
with
and
ζ
= 0.1.
Because the dimensionless shaping times are 0.597, 1.194, and 1.692, SD shapers are composed of three, four, and five impulses, respectively, and 5% insensitivities
I
0.05
are 0.088, 0.452, and 1.133, respectively.

https://doi.org/10.1371/journal.pone.0276669.g006
Fig 7
shows unit-step responses when these SD shapers are applied to the underdamped system with
ζ
= 0.1.
Fig 7A
shows unit-step responses when no modeling error exists (
), and
Fig 7B
shows unit-step responses when 20% modeling error exists in natural frequency (
ω
n
= 2.4 Hz,
). If no modeling error exists, the SD shapers exhibit the exact shaping time of 0.3, 0.6 and 0.85 s as designed in Eqs (
23
), (
25
) and (
27
), respectively. If 20% modeling error occurs in natural frequency, the magnitude of residual vibration decreases as the shaping time increases as shown in
Fig 7B
, which reveals the same trend as the sensitivity curve in
Fig 6
.

Download:

PNG
larger image
TIFF
original image
Fig 7.
Unit-step responses when SD shapers with shaping time
t
s
= 0.3, 0.6, 0.85s are used for an underdamped system
with
and
ζ
= 0.1.
(A) No modeling error (
), (B) 20% modeling error in natural frequency (
ω
n
= 2.4 Hz,
).

https://doi.org/10.1371/journal.pone.0276669.g007
If
t
s
= 0.85 s for the underdamped system with
and
ζ
= 0.1, the dimensionless shaping time is
T
s
≈ 1.692 > 1.5, and thus, the SD shaper is composed of five impulses according to (7). However, considering the impulse vector diagram for the case with
T
s
> 1.5, an input shaper with three, four or six impulses can also be designed because the sum of the corresponding impulse vectors can be made to be zero.
Table 1
shows the SD shaper with five impulses and three input shapers with three, four and six impulses obtained using design procedures (ii) and (iii) and excluding procedure (i) in the previous section. However, this input shaper with three, four or six impulses is not a satisfactory solution because of its high sensitivity to modeling errors compared to the SD shaper.
Fig 8
shows that the SD shaper with five impulses for
t
s
= 0.85 s is less sensitive to modeling errors than the input shaper with three, four or six impulses for
t
s
= 0.85 s since the
V
values of the SD shaper are smaller than the
V
values of the three input shapers around
in
Fig 8
. The main reason for the improved robustness of the SD shapers than input shapers with less impulses for the same shaping time is that SD shapers satisfy more derivative constraints than input shapers with less impulses. However, it may not be possible to obtain an input shaper with six impulses because all nine constraints [
Eq (26)
plus two more third-order derivative constraints] cannot be satisfied for the case with
t
s
= 0.85 s. Even if we find an input shaper with six impulses for
t
s
= 0.85 s numerically, it may be an approximate solution with less robustness such that the sum of
A
i
is close to one (not exact one) as can be checked in the last input shaper of
Table 1
. Note that the sum of
A
i
of the first three input shapers in
Table 1
are exactly one, but the sum of
A
i
of the last input shaper is 1.0024.

Download:

PNG
larger image
TIFF
original image
Fig 8.
Sensitivity curves of the SD shaper with five impulses and input shapers with three, four or six impulses for the case with the shaping time
t
s
= 0.85 s for an underdamped system with
and
ζ
= 0.1.
The SD shaper with five impulses is more robust to modeling errors than the input shaper with three, four or six impulses for
t
s
= 0.85 s.

https://doi.org/10.1371/journal.pone.0276669.g008

Download:

PNG
larger image
TIFF
original image
Table 1.
The SD shaper with five impulses and input shapers with three, four and six impulses for the shaping time
t
s
= 0.85s.

https://doi.org/10.1371/journal.pone.0276669.t001

Experimental demonstrations
A horizontal beam vibration apparatus that can control a flexible horizontal cantilever beam vertically was used to verify the validity of the SD shapers (
Fig 9
) [
28
]. In the apparatus, vertical motions of the flexible beam were generated by controlling a servo motor with a timing belt. A laser sensor was placed under the tip of the beam to collect motion data. Note that the sensor data were not used to control the flexible beam because the input shaper is the open-loop control method. To enhance the effect of the residual vibration, masses were attached at the end of the flexible beam.

Download:

PNG
larger image
TIFF
original image
Fig 9.
Horizontal beam vibration apparatus.
The apparatus can generate vertical movements with a flexible horizontal cantilever beam with added masses at the end. The system has natural frequency
ω
n
of 16.7 Hz and damping ratio
ζ
of 0.002.

https://doi.org/10.1371/journal.pone.0276669.g009
The block diagram of the horizontal beam vibration apparatus is shown in
Fig 10
. As shown in
Fig 10
, the major three components of the block diagram are a proportional-and-derivative (PD) controller, motor dynamics, and beam dynamics. The PD feedback controller
with
K
p
= 0.18 and
T
d
= 0.075 was used to control the vertical motion of the moving block. The motor dynamics was modeled as a first-order system with a time constant of 16 ms. Here,
K
LG
represents the conversion factor of the linear guide with the timing belt, which generates linear motion from the motor rotation. The flexible horizontal beam dynamics was modeled as a second-order underdamped system
.

Download:

PNG
larger image
TIFF
original image
Fig 10.
Block diagram of the horizontal beam vibration apparatus.
The system includes a PD controller, motor dynamics, and beam dynamics.

https://doi.org/10.1371/journal.pone.0276669.g010
To validate the design methodology and performance of the SD shaper, three SD shapers with three shaping times, namely 0.2, 0.5, and 0.7 s, were used. To obtain SD shapers, determining the natural frequency
ω
n
and damping ratio
ζ
of the system, which is the horizontal beam vibration apparatus in the experiments, is crucial. Based on 10 s of the step response of the flexible beam with no input shaper, the natural frequency and damping ratio were modeled to be 16.7 Hz and 0.002, respectively. With these modeled values, the damped period
T
d
was calculated to be 0.3762. The black dotted line in
Fig 11
represents the experimental step response of the flexible beam with no input shaper.

Download:

PNG
larger image
TIFF
original image
Fig 11.
Step responses through simulations and experiments.
The command input of step responses is 100 mm step. Black: No input shaper. Blue: SD shaper with 0.2 s shaping time. Red: SD shaper with 0.5 s shaping time. Brown: SD shaper with 0.7 s shaping time.

https://doi.org/10.1371/journal.pone.0276669.g011
Based on
Eq (6)
, the SD shaper with a shaping time of 0.2 s has the dimensionless shaping time
T
s
of 0.53, and then, according to
Eq (7)
, three impulses are required (
N
= 3). In the same manner, the SD shaper with a shaping time of 0.5 s requires four impulses (
N
= 4) because
T
s
is 1.33, and the SD shaper with a shaping time of 0.7 s requires five impulses (
N
= 5) because
T
s
is 1.86. Based on these parameters, the SD shaper for each case was obtained by solving numerically using MATLAB.
Table 2
summarizes the obtained three SD shapers.

Download:

PNG
larger image
TIFF
original image
Table 2.
Obtained SD shapers with three shaping times.

https://doi.org/10.1371/journal.pone.0276669.t002
To validate the performance of the obtained SD shapers, both simulations and experiments were performed with a 100 mm step input displacement, and step responses were compared. The three obtained SD shapers were simulated using Simulink, and the results are shown in
Fig 11
. The simulations were performed according to the block diagram displayed in
Fig 10
. The solid blue, red, and green lines in
Fig 11
depict the simulation results for SD shapers with shaping times of 0.2, 0.5, and 0.7 s, respectively. The SD shapers were experimented using the apparatus, and the results are plotted in
Fig 11
for comparison. The dotted blue, red, and green lines in
Fig 11
denote the experimental results of the SD shapers with 0.2, 0.5, and 0.7 s shaping times, respectively. As the specified shaping time is increased, the actual shaping time is also increased in both simulation and experimental results, as expected. In
Fig 11
, the experimental results of the SD shapers have the same trend as the simulation results with very little discrepancy.

Fig 12
shows the simulation results to compare the effect of the motor feedback dynamics for the obtained SD shapers. For both solid and dashed lines in
Fig 12
, red, blue, and green colors represent the results of SD shapers with shaping times 0.2, 0.5, and 0.7 s, respectively. The dashed lines represent the results of SD shapers when the system model does not include motor feedback dynamics, and thus, the system model only consists of the SD shaper and beam dynamics. The solid lines represent the results of SD shapers when the system model includes motor feedback dynamics, which is the system model based on the experimental apparatus. The solid lines in
Fig 12
are the same as the simulation results of SD shapers in
Fig 11
. The dashed lines in
Fig 12
reveal that the SD shapers can suppress residual vibration in a planned specified shaping time. However, the solid lines in
Fig 12
show that the shaping times of the SD shapers are slightly longer than the specified shaping times because of the motor feedback dynamics. The experimental and simulation results reveal that SD shapers with arbitrarily specified shaping times can eliminate the residual vibrations of a flexible beam as expected.

Download:

PNG
larger image
TIFF
original image
Fig 12.
Simulation results for three SD shapers with and without motor feedback dynamics.

https://doi.org/10.1371/journal.pone.0276669.g012

Conclusion
This paper introduced an SD shaper that can suppress residual vibrations at an arbitrarily specified shaping time with the largest 5% insensitivity to modeling errors when the magnitude of the last impulse is taken as a free parameter. From the specified shaping time, the angle of the last impulse vector was calculated, and the number of impulses was determined from the dimensionless shaping time. Then, by solving the constraint equations of the SD shaper analytically or numerically, the desired SD shaper was obtained. As the specified shaping time increases, the SD shaper increases the number of impulses one by one according to the number of added derivative constraints to improve the robustness to modeling errors.

The performance of the SD shapers for specified shaping times were evaluated through simulations by plotting sensitivity curves for modeling errors and step responses of the second-order system. The simulations revealed that the SD shapers suppress residual vibrations of the vibratory system at the specified shaping times, and the robustness of the SD shapers to modeling errors improves as the shaping time increases. Furthermore, the validity of the SD shapers was experimentally verified using a horizontal beam vibration apparatus.

For further study, the SD shaper proposed in this paper could be extended to a negative SD shaper or a multi-mode SD shaper.

Acknowledgments

The author acknowledges Prof. Manh-Tuan Ha and Prof. Chul-Goo Kang for lending the horizontal beam vibration apparatus for experimental works.

References
1.

Tanveer M, Usman M, Khan IU, Ahmad S, Hanif A, Farooq SH. Application of tuned column ball damper (TLCBD) for improved vibration control performance of multi-storey structure. PLoS ONE. 2019;14(10): e0224436. pmid:31648266

View Article

PubMed/NCBI

Google Scholar

2.

Hong KT, Huh CD, Hong KS. Command shaping control for limiting the transient sway angle of crane systems. International Journal of Control, Automation and Systems. 2003;1(1):43–53.

View Article

Google Scholar

3.

Sorensen KL, Singhose W, Dickerson S. A controller enabling precise positioning and sway reduction in bridge and gantry cranes. Control Engineering Practice. 2007;15:825–837.

View Article

Google Scholar

4.

Alghanim KA, Alhazza KA, Masoud ZN. Discrete-time command profile for simultaneous travel and hoist maneuvers of overhead cranes. Journal of Sound and Vibration. 2015;345:47–57.

View Article

Google Scholar

5.

Runghimmawan W, Chatlatanagulchai W. Track-seeking control using input-shaping method to reduce vibration in HDD. Proc. 13
th
International Symposium on Communications and Information Technologies. 2013;1:576–579.
6.

Singer N, Tanquary M, Pasch K. System for removing selected unwanted frequencies in accordance with altered settings in a user interface of a data storage device. United States patent US 6314473B1. 2001 Nov 6.
7.

Guagliumi L, Berti A, Monti E, Carricato M. Antisloshing trajectories for high-acceleration motions in automatic machines. Journal of Dynamic Systems, Measurement, and Control. 2022;144(7):071006.

View Article

Google Scholar

8.

Pridgen B, Bai K, Singhose W. Shaping container motion for multimode and robust slosh suppression. Journal of Spacecraft and Rockets. 2013;50(2):440–448.

View Article

Google Scholar

9.

Alshaya A, Alghanim K. Command shaping for sloshing suppression of a suspended liquid container. Journal of Dynamic Systems, Measurement, and Control. 2020;142(12):121003.

View Article

Google Scholar

10.

Alshaya A, Alshayji A. Robust multi-steps input command for liquid sloshing control. Journal of Vibration and Control. 2021;

View Article

Google Scholar

11.

Park J, Chang PH, Park HS, Lee E. Design of learning input shaping technique for residual vibration suppression in an industrial robot. IEEE/ASME Transactions on Mechatronics. 2006;11(1):55–65.

View Article

Google Scholar

12.

Sahoo PP, Singhose W. Effects of Input shaping on impact loads during collisions involving flexible robots. IEEE 14
th
International Conference on Control and Automation. 2018:1004–1009.
13.

Kim D, Kim SH, Kim T, Kang BB, Lee M, Park W, et al. Review of machine learning methods in soft robotics. PLoS One. 2021;16(2): e0246102. pmid:33600496

View Article

PubMed/NCBI

Google Scholar

14.

Smith OJM. Posicast control of damped oscillatory systems. Proceedings of the IRE. 1957;45:1249–1255.

View Article

Google Scholar

15.

Singer N. Residual vibration reduction in computer-controlled machines [dissertation]. Massachusetts Institute of Technology; 1989.
16.

Singer N, Seering W. Preshaping command inputs to reduce system vibration. Journal of Dynamic Systems, Measurement and Control. 1990;112:76–82.

View Article

Google Scholar

17.

Hyde JM, Seering WP. Using input command pre-shaping to suppress multiple mode vibration. Proc. 1991 IEEE International Conference on Robotics and Automation. 1991;2604–2609.
18.

Singhose W. Command generation for flexible systems [dissertation]. Massachusetts Institute of Technology; 1997.
19.

Singh T, Vadali SR. Robust time-delay control of multimode systems. International Journal of Control. 1995;62(6):1319–1339.

View Article

Google Scholar

20.

Xie X, Huang J, Liang Z. Vibration reduction for flexible systems by command smoothing. Mechanical Systems and Signal Processing. 2013;39(1–2):461–470.

View Article

Google Scholar

21.

Kang CG. Impulse vectors for input-shaping control: A mathematical tool to design and analyze input shapers. IEEE Control Systems Magazine. 2019;39(4):40–55.

View Article

Google Scholar

22.

Singer NC, Seering WP. An extension of command shaping methods for controlling residual vibration using frequency sampling, Proc. 1992 IEEE International Conference on Robotics and Automation. 1992;800–805.
23.

Singhose W, Seering W, Singer N. Input shaping for vibration reduction with specified insensitivity to modeling errors. Proc. of Japan-USA Symposium on Flexible Automation. 1996;1:307–313.
24.

Vaughan J, Yano A, Singhose W. Comparison of robust input shapers. Journal of Sound and Vibration. 2008;315:797–815.

View Article

Google Scholar

25.

Singhose W. Command shaping for flexible systems: A review of the first 50 years. International Journal of Precision Engineering and Manufacturing. 2009;10(4):153–168.

View Article

Google Scholar

26.

Singhose W, Kim D, Kenison M. Input shaping control of double-pendulum bridge crane oscillations. Journal of Dynamic Systems, Measurement, and Control. 2008;130:034504.

View Article

Google Scholar

27.

Boscariol P, Richiedei D. Robust point-to-point trajectory planning for nonlinear underactuated systems: Theory and experimental assessment. Robotics and Computer-Integrated Manufacturing. 2018;50:256–265.

View Article

Google Scholar

28.

Ha MT, Hassan R, Kang CG. Development of a benchmark testbed for input shaping control, Proc. of International Conference on Control, Automation and Systems. 2018;1:1370–1375.

Download PDF

Citation

XML

Print

Share

Reddit

Facebook

LinkedIn

Mendeley

Bluesky

Email

Advertisement

Subject Areas
?

For more information about PLOS Subject Areas, click

here
.

We want your feedback.
Do these Subject Areas make sense for this article? Click the target next to the incorrect Subject Area and let us know. Thanks for your help!

Vibration engineering

Is the Subject Area
"Vibration engineering"
applicable to this article?

Yes

No

Thanks for your feedback.

Vibration

Is the Subject Area
"Vibration"
applicable to this article?

Yes

No

Thanks for your feedback.

Control systems

Is the Subject Area
"Control systems"
applicable to this article?

Yes

No

Thanks for your feedback.

Computer modeling

Is the Subject Area
"Computer modeling"
applicable to this article?

Yes

No

Thanks for your feedback.

Earthquake engineering

Is the Subject Area
"Earthquake engineering"
applicable to this article?

Yes

No

Thanks for your feedback.

Simulation and modeling

Is the Subject Area
"Simulation and modeling"
applicable to this article?

Yes

No

Thanks for your feedback.

System stability

Is the Subject Area
"System stability"
applicable to this article?

Yes

No

Thanks for your feedback.

Soft robotics

Is the Subject Area
"Soft robotics"
applicable to this article?

Yes

No

Thanks for your feedback.

Publications

PLOS Aging and Health

PLOS Biology

PLOS Climate

PLOS Complex Systems

PLOS Computational Biology

PLOS Digital Health

PLOS Ecosystems

PLOS Genetics

PLOS Global Public Health

PLOS Medicine

PLOS Mental Health

PLOS Neglected Tropical Diseases

PLOS One

PLOS Pathogens

PLOS Sustainability and Transformation

PLOS Water

Home

Blogs

Collections

Give feedback

LOCKSS

Privacy Policy

Terms of Use

Advertise

Media Inquiries

Contact

PLOS is a nonprofit 501(c)(3) corporation, #C2354500, based in California, US
</reference>

<statements>
1. To eliminate residual structural vibration without introducing closed-loop phase lags during high-throughput stop-and-strobe acquisition, Zero Vibration Derivative (ZVD) input shaping filters the trajectory setpoints
2. Convolving f_ZVD(t) with the reference command profile r(t) produces a shaped input that cancels residual structural vibrations at ω_n, eliminating motion-induced image blur during high-speed step-and-scan indexing
3. In uncompensated indexing stages, settling delays (>140 ms) dominate the duty cycle to prevent motion-induced blur. Implementing ZVD input shaping suppresses dominant torsional resonances at ω_n, shortening settling times to 36 ms and permitting strobe triggering during smooth deceleration phases
</statements>

Begin the assessment now. Output only the JSON list, without any conversational text or explanations.