Journal cover
Journal topic
**Mechanical Sciences**
An open-access journal for theoretical and applied mechanics

Journal topic

**Research article**
11 Jan 2021

**Research article** | 11 Jan 2021

Kinematic synthesis method and eccentricity effects of a Stephenson mechanism

^{1}U. A. Joldasbekov Institute of Mechanics and Mechanical Engineering, Almaty, Kazakhstan^{2}Turkish Machine Theory Association, Turkey^{3}Department of Mechanics and Mathematics, Al-Farabi Kazakh National University, Almaty, Kazakhstan

^{1}U. A. Joldasbekov Institute of Mechanics and Mechanical Engineering, Almaty, Kazakhstan^{2}Turkish Machine Theory Association, Turkey^{3}Department of Mechanics and Mathematics, Al-Farabi Kazakh National University, Almaty, Kazakhstan

**Correspondence**: Recep Halicioglu (recephalicioglu@gmail.com)

**Correspondence**: Recep Halicioglu (recephalicioglu@gmail.com)

Abstract

When implementing the technological process on crank presses, it is necessary to provide a predetermined working cycle of the slider motion: fast lifting, dwell, and slow lowering. The cycle cannot be realized without controlling the motor. In addition, using controllable motors increases the manufacturing cost. Due to the geometric and kinematic capabilities of the mechanism, changing the kinematics of the working link is the best choice. Thanks to the use of the Stephenson II mechanism, the slider skew is eliminated due to the parallel connecting rods and the increased area of slider contact. This study presents a numerical method for kinematic synthesis of the Stephenson mechanism that has kinematic advantages. The method is based on mean square deviation which is the minimizing of an objective function. Thanks to the proposed synthesis method, approximate dwell movement can be performed when the slider is on the bottom dead center. In this study, values of the crank length and parallel connecting rods' lengths, angular coordinates of the crank and connecting rods, and the eccentricity of the guide slider relative to the crank rotation axis were obtained. It is observed that eccentricity affects the lower forward and higher backward speed of the slider. The kinematic results of the slider movement are comparatively presented in this article.

Download & links

How to cite

How to cite.

Tuleshov, A., Halicioglu, R., Shadymanova, A., and Kuatova, M.: Kinematic synthesis method and eccentricity effects of a Stephenson mechanism, Mech. Sci., 12, 1–8, https://doi.org/10.5194/ms-12-1-2021, 2021.

1 Introduction

The main disadvantage of a four-bar slider-crank press scheme is the impossibility of ensuring a long dwell of a working slide and the irregularity of its stroke (Bocharev, 2008). This is not possible because there is no intermediate link between the drive crank and the connecting rod bonded with a slider. Jomartov et al. (2020) showed that a crank press based on the Stephenson mechanism allows solving many problems characterized by crank presses based on four-bar mechanisms. For example, to correct a slider deformation caused by parallel connecting rods, the contact area of the slider needs to be increased. By using a simplified structure, they tried to achieve a dwell motion on slider (Jomartov et al., 2019; Tuleshov et al., 2019; Drakunov et al., 2018).

The development of new designs of machinery mechanisms, including crank presses (Bocharev, 2008), begins with the solution of analysis and synthesis problems based on mathematical modeling. When implementing the technological process in crank presses, it is necessary to provide a predetermined cycle of the working slider movement: fast lifting (reverse stroke), dwell, and slow lowering (working stroke). In studies conducted on crank presses, there are two ways to achieve these goals. First, these properties can be realized by modifying the kinematic chain in a mechanism with one degree of freedom (1 DOF). Tuleshov and colleagues (Jomartov and Tuleshov, 2019; Jomartov et al., 2019; Tuleshov et al., 2019; Jomartov and Tuleshov, 2018) presented an application of the four-class kinematic chain (structure group) to perform the synthesis of the crank press mechanism. The method of kinematic analysis of two connecting rod crank presses was developed by Jomartov et al. (2020) based on four-bar mechanism (Jomartov and Tuleshov, 2019). The second way is to solve this problem through the additional freedom of kinematic chain, which is called the hybrid press system (Kütük and Dülger, 2016; Tokuz, 1992; Dülger and Kireçci, 1995, 1999). Kütük and Dülger (2016) presented a scientific literature review on the analysis and synthesis of the hybrid mechanisms of the crank press. To provide a technological cycle, they studied a seven-bar mechanism that had 2 DOF and whose one degree was driven by a constant speed direct current power motor (for the implementation of the main technological process) and whose second degree was driven by a servomotor. In many other studies, hybrid press systems based on five-bar and seven-bar mechanisms with two degrees of freedom have been examined. A similar technology was carried out by Tokuz, Dülger, and their colleagues in a hybrid configuration (Tokuz, 1992; Dülger and Kireçci, 1995, 1999). The constant speed motor and servomotor were combined by a differential gearbox, which subsequently drove the crank mechanism (Tokuz, 1992).

Yuan et al. (2005) investigated two different mechanisms whose DOF was two. Ouyang et al. (2004) proposed a five-bar linkage mechanism, which consisted of a constant velocity alternating current motor and a speed controller, a brushless alternating current servomotor and a gear servo amplifier, a shift encoder, a flywheel, and a belt. Zhang also proposed a hybrid five-bar mechanism with an adjustable link (Zhang, 2006).

Connor et al. (1995) conducted studies on the synthesis of five-bar hybrid mechanisms by using genetic algorithms. Dülger et al. (2003) carried out a study on modeling and kinematic analysis of a hybrid actuator of the seven-bar mechanism with a controllable crank. Yu (2006) performed a study by using the hybrid machine (HM) system with a five-bar mechanism. Li and Zhang (2010) applied a seven-bar configuration through kinematic analysis and optimal design of the hybrid system. Li and Tso (2008) introduced a seven-bar mechanism that was later used to study the performance of stamping and energy distribution between a servo motor and a flywheel with various motion inputs (Tso and Li, 2008). In addition, Tso (2010) also developed a control system for the seven-bar mechanism by using iterative training management and feedback control techniques.

In all above-mentioned mechanisms, the issue of providing the necessary slider motion cycle was solved by controlling two or more motors, and accordingly, the tasks of dynamic synthesis of drive control functions were achieved. The implementation of the technological cyclogram by a mechanism with one degree of freedom requires a significant complication of the structure linkage kinematic chain which is called the Assur group (Yevgrafov et al., 2006; Huang and Ding, 2020).

Metal-forming equipment is commonly based on the four-bar structural group mechanisms. However, six-bar structural group mechanisms, such as Stephenson II mechanisms, have broader functionality. In this study, an analytical kinematic synthesis method for a Stephenson mechanism that provides balancing was developed for metal-forming process. This synthesis was carried out based on the mean-square minimization of the objective function which was obtained taking into account the geometric features of the Stephenson mechanism structure. The effect of eccentricity was also investigated. Kinematic results were obtained by using the mathematical relations (position, speed, and acceleration). It was assessed that the eccentricity affected the possibility of ensuring the approximate slider dwell and the increase of the average slider velocity coefficient.

2 System definition and kinematic model

Figure 1 demonstrates the kinematic scheme of a simplified six-bar
Stephenson mechanism. The feature of the mechanism is that the variable
contour BB^{′}C^{′}C is a parallelogram and a triangle ABB^{′} is equilateral. In
this figure, *r*, *a*, and *l* indicate crank length, ABB^{′} triangle height, and parallel connecting
rods length of BC = B^{′}C^{′}, respectively. *φ*, *ψ*, *S*, *e* and *b* show the angle coordinate of
a crank, angle coordinate of two connecting rods, linear coordinate of a
slider, slider eccentricity (i.e., the deviation of a trajectory of a slider
center gravity from the axis O_{y}), and distance between joint C and slider center
along an axis O_{x}, respectively.

The kinematics of the Stephenson crank press mechanism have the following Eq. (1) (Jomartov et al., 2020):

$$\begin{array}{}\text{(1)}& \left\{\begin{array}{l}r\mathrm{cos}\mathit{\phi}+l\mathrm{cos}\mathit{\psi}=e\\ r\mathrm{sin}\mathit{\phi}-a+l\mathrm{sin}\mathit{\psi}=-S.\end{array}\right.\end{array}$$

A clear solution of Eq. (1) in relation to *S*=*S*(*φ*), *ψ*=*ψ*(*φ*) is obtained in Eq. (2) (the ± signs
correspond to various assemblies of the mechanism):

$$\begin{array}{}\text{(2)}& \left\{\begin{array}{l}S=a-r\mathrm{sin}\mathit{\phi}\pm \sqrt{{l}^{\mathrm{2}}-(e-r\mathrm{cos}\mathit{\phi}{)}^{\mathrm{2}}}\\ \mathit{\psi}=\pm \mathrm{arccos}\left[\frac{\mathrm{1}}{l}(e-r\mathrm{cos}\mathit{\phi})\right].\end{array}\right.\end{array}$$

According to the scheme on Fig. 1, rods of 4 and 5 rotate around the
joints C and C^{′} (*l* – radius of joints). The trajectories of joints
B and B^{′} are a circle in the interval $\mathit{\phi}\in \left({\mathit{\phi}}^{\left(\mathrm{1},H\right)},{\mathit{\phi}}^{\left(\mathrm{1},K\right)}\phantom{\rule{0.125em}{0ex}}\right)$. It is necessary for
the slider to work well. Equation (3) is a vector form of the mechanism
scheme in Fig. 1. The scheme of the Stephenson press mechanism is based on
vector Eq. (3) that is presented in Fig. 2.

$$\begin{array}{}\text{(3)}& \mathit{S}+\mathit{e}=\mathit{r}+\mathit{a}+\mathit{l}\end{array}$$

3 Synthesis method

The origin of the *x*, *y* coordinate system is in the center of the rotation of
the crank O. The slider stroke *h*_{x} is considered as given. According to
Fig. 2, when rod and crank are on parallel lines, the slider is in extreme
positions. Considering the geometric positions, the following statements can
be written:

$$\begin{array}{c}{\displaystyle}{S}_{max}=r\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{1}\right)}+l\mathrm{sin}{\mathit{\psi}}^{\left(\mathrm{1}\right)}+a\\ {\displaystyle}{S}_{min}=r\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{2}\right)}+l\mathrm{sin}{\mathit{\psi}}^{\left(\mathrm{2}\right)}+a.\end{array}$$

By subtracting the first equation from the second, we obtain the slider stroke:

$$\begin{array}{}\text{(4)}& \begin{array}{rl}{h}_{x}& ={S}_{max}-{S}_{min}=r(\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{1}\right)}-\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{2}\right)})\\ & +l(\mathrm{sin}{\mathit{\psi}}^{\left(\mathrm{1}\right)}-\mathrm{sin}{\mathit{\psi}}^{\left(\mathrm{2}\right)}).\end{array}\end{array}$$

Given that *φ*^{(1)}=*ψ*^{(1)} and ${\mathit{\phi}}^{\left(\mathrm{2}\right)}+\mathit{\pi}={\mathit{\psi}}^{\left(\mathrm{2}\right)}$ from Eq. (1), we have the following:

$$\begin{array}{}\text{(5)}& {\displaystyle}\begin{array}{rl}& {h}_{x}=r(\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{1}\right)}-\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{2}\right)})+l(\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{1}\right)}-\mathrm{sin}({\mathit{\phi}}^{\left(\mathrm{2}\right)}+\mathit{\pi}\left)\right),\\ & {S}_{max}=(r+l)\mathrm{sin}{\mathit{\phi}}^{\left(\mathrm{1}\right)}+a,\\ & r\mathrm{cos}{\mathit{\phi}}^{\left(\mathrm{2}\right)}+l\mathrm{cos}\left({\mathit{\phi}}^{\left(\mathrm{2}\right)}+\mathit{\pi}\right)=e,\\ & r\mathrm{cos}{\mathit{\phi}}^{\left(\mathrm{1}\right)}+l\mathrm{cos}{\mathit{\phi}}^{\left(\mathrm{1}\right)}=e,\\ & \mathrm{cos}{\mathit{\phi}}^{\left(\mathrm{1}\right)}=e/\left(r+l\right),\phantom{\rule{0.125em}{0ex}}\mathrm{cos}{\mathit{\phi}}^{\left(\mathrm{2}\right)}=e/\left(-r+l\right).\end{array}\end{array}$$

Based on Fig. 2, considering the geometric positions, the following statements can be written:

$$\begin{array}{c}{\displaystyle}{x}_{\mathrm{A}}^{\mathrm{2}}+{y}_{\mathrm{A}}^{\mathrm{2}}={r}^{\mathrm{2}},\\ {\displaystyle}{x}_{\mathrm{O}}={y}_{\mathrm{O}}=\mathrm{0},\\ {\displaystyle}{\left({x}_{\mathrm{C}}-{x}_{\mathrm{B}}\right)}^{\mathrm{2}}+{\left({y}_{C}-{y}_{B}\right)}^{\mathrm{2}}={l}^{\mathrm{2}},\\ {\displaystyle}{x}_{\mathrm{B}}={x}_{\mathrm{A}},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{y}_{\mathrm{B}}-{y}_{\mathrm{A}}=a,\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{x}_{\mathrm{C}}=e,\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{y}_{\mathrm{C}}=S.\end{array}$$

Then

$${\left(e-{x}_{\mathrm{A}}\right)}^{\mathrm{2}}+{\left(S-{y}_{\mathrm{A}}-a\right)}^{\mathrm{2}}={l}^{\mathrm{2}}$$

or

$${e}^{\mathrm{2}}-\mathrm{2}e{x}_{\mathrm{A}}+{x}_{\mathrm{A}}^{\mathrm{2}}+{\left(S-a\right)}^{\mathrm{2}}-\mathrm{2}\left(S-a\right){y}_{\mathrm{A}}+{y}_{\mathrm{A}}^{\mathrm{2}}={l}^{\mathrm{2}}.$$

When statements of ${x}_{\mathrm{A}}=r\mathrm{cos}\mathit{\phi},\phantom{\rule{0.125em}{0ex}}{y}_{\mathrm{A}}=r\mathrm{sin}\mathit{\phi}$ are written on their places, Eq. (6) is obtained:

$$\begin{array}{}\text{(6)}& {l}^{\mathrm{2}}={e}^{\mathrm{2}}+{r}^{\mathrm{2}}+{\left(S-a\right)}^{\mathrm{2}}-\mathrm{2}r\left(S-a\right)\mathrm{sin}\mathit{\phi}-\mathrm{2}re\mathrm{cos}\mathit{\phi}.\end{array}$$

To find the parameters *r*, *a*, *l*, *e*, and *φ*_{0}, which satisfy Eq. (5)
and $S={S}_{max}={S}_{m}$ for all values, the following can be written:

$$\begin{array}{}\text{(7)}& {\displaystyle}\begin{array}{rl}& \mathit{\phi}={\mathit{\phi}}_{\mathrm{0}}+i\mathrm{\Delta}\mathit{\phi},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}i=-\mathrm{2},-\mathrm{1},\mathrm{0},\mathrm{1},\mathrm{2}\\ & {l}^{\mathrm{2}}={e}^{\mathrm{2}}+{r}^{\mathrm{2}}+{\left({S}_{m}-a\right)}^{\mathrm{2}}-\mathrm{2}r\left({S}_{m}-a\right)\mathrm{sin}\left({\mathit{\phi}}_{\mathrm{0}}+i\mathrm{\Delta}\mathit{\phi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.125em}{0ex}}-\mathrm{2}re\mathrm{cos}\left({\mathit{\phi}}_{\mathrm{0}}+i\mathrm{\Delta}\mathit{\phi}\right),\\ & i=-\mathrm{2},-\mathrm{1},\mathrm{1},\mathrm{2}\\ & {l}^{\mathrm{2}}={e}^{\mathrm{2}}+{r}^{\mathrm{2}}+{\left({S}_{m}-a\right)}^{\mathrm{2}}-\mathrm{2}r\left({S}_{m}-a\right)\mathrm{sin}{\mathit{\phi}}_{\mathrm{0}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.125em}{0ex}}-\mathrm{2}re\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}}.\end{array}\end{array}$$

After the transformations of Eq. (7),

$$\begin{array}{rl}& {\mathrm{\Delta}}_{i}\left({\mathit{\phi}}_{i},{S}_{m}\right)={e}^{\mathrm{2}}+{r}^{\mathrm{2}}+{\left({S}_{m}-a\right)}^{\mathrm{2}}-{l}^{\mathrm{2}},\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-\mathrm{2}r\left({S}_{m}-a\right)\mathrm{sin}\left({\mathit{\phi}}_{\mathrm{0}}+i\mathrm{\Delta}\mathit{\phi}\right)-\mathrm{2}re\mathrm{cos}\left({\mathit{\phi}}_{\mathrm{0}}+i\mathrm{\Delta}\mathit{\phi}\right),\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}i=-\mathrm{2},-\mathrm{1},\mathrm{0},\mathrm{1},\mathrm{2},\end{array}$$

or

$$\begin{array}{}\text{(8)}& {\displaystyle}\begin{array}{rl}& {\mathrm{\Delta}}_{i}\left({\mathit{\phi}}_{i},{S}_{m}\right)={e}^{\mathrm{2}}+{r}^{\mathrm{2}}+{\left({S}_{m}-a\right)}^{\mathrm{2}}-{l}^{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-\mathrm{2}r\left({S}_{m}-a\right)\phantom{\rule{0.125em}{0ex}}\left[\mathrm{sin}{\mathit{\phi}}_{\mathrm{0}}\mathrm{cos}i\mathrm{\Delta}\mathit{\phi}+\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}}\mathrm{sin}i\mathrm{\Delta}\mathit{\phi}\right]\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-\mathrm{2}re\left[\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}}\mathrm{cos}i\mathrm{\Delta}\mathit{\phi}-\mathrm{sin}{\mathit{\phi}}_{\mathrm{0}}\mathrm{sin}i\mathrm{\Delta}\mathit{\phi}\right]\phantom{\rule{0.125em}{0ex}},\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}i=-\mathrm{2},-\mathrm{1},\mathrm{0},\mathrm{1},\mathrm{2},\end{array}\end{array}$$

which can be presented as

$$\begin{array}{}\text{(9)}& {\mathrm{\Delta}}_{i}={P}_{\mathrm{0}}{f}_{\mathrm{0}i}+{P}_{\mathrm{1}}{f}_{\mathrm{1}i}+{P}_{\mathrm{2}}{f}_{\mathrm{2}i}+{P}_{\mathrm{3}}{f}_{\mathrm{3}i}+{P}_{\mathrm{4}}{f}_{\mathrm{4}i}-{F}_{i}.\end{array}$$

Here the following signs are inputted:

$$\begin{array}{c}{\displaystyle}{P}_{\mathrm{0}}={e}^{\mathrm{2}}+{r}^{\mathrm{2}}-{l}^{\mathrm{2}},\\ {\displaystyle}{P}_{\mathrm{1}}=-\mathrm{2}r\left({S}_{m}-a\right)\mathrm{sin}{\mathit{\phi}}_{\mathrm{0}},\\ {\displaystyle}{P}_{\mathrm{2}}=-\mathrm{2}r\left({S}_{m}-a\right)\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}},\\ {\displaystyle}{P}_{\mathrm{3}}=-\mathrm{2}re\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}},\\ {\displaystyle}{P}_{\mathrm{4}}=\mathrm{2}re\mathrm{sin}{\mathit{\phi}}_{\mathrm{0}},\\ {\displaystyle}{f}_{\mathrm{0}i}=\mathrm{1},\\ {\displaystyle}{f}_{\mathrm{1}i}=\mathrm{cos}i\mathrm{\Delta}\mathit{\phi},\\ {\displaystyle}{f}_{\mathrm{2}i}=\mathrm{sin}i\mathrm{\Delta}\mathit{\phi},\\ {\displaystyle}{f}_{\mathrm{3}i}=\mathrm{cos}i\mathrm{\Delta}\mathit{\phi},\\ {\displaystyle}{f}_{\mathrm{4}i}=\mathrm{sin}i\mathrm{\Delta}\mathit{\phi},\\ {\displaystyle}{F}_{\mathrm{4}i}={\left({S}_{m}-a\right)}^{\mathrm{2}}.\end{array}$$

When the statements of $\mathrm{sin}{\mathit{\phi}}_{\mathrm{0}}=\frac{({S}_{m}-a)}{(r+l)}$ , $\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}}=\frac{e}{\left(r+l\right)}$ are written on their place, the following statements are obtained:

$$\begin{array}{c}{\displaystyle}{P}_{\mathrm{0}}={e}^{\mathrm{2}}+{r}^{\mathrm{2}}-{l}^{\mathrm{2}},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{P}_{\mathrm{1}}=-\mathrm{2}r{\displaystyle \frac{({S}_{m}-a{)}^{\mathrm{2}}}{(r+l)}},\\ {\displaystyle}{P}_{\mathrm{2}}=-\mathrm{2}r\left({S}_{m}-a\right){\displaystyle \frac{e}{\left(r+l\right)}}\\ {\displaystyle}{P}_{\mathrm{3}}=-\mathrm{2}r{\displaystyle \frac{{e}^{\mathrm{2}}}{\left(r+l\right)}},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{P}_{\mathrm{4}}=\mathrm{2}re{\displaystyle \frac{({S}_{m}-a)}{(r+l)}}\\ {\displaystyle}{f}_{\mathrm{0}i}=\mathrm{1},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{f}_{\mathrm{1}i}=\mathrm{sin}i\mathrm{\Delta}\mathit{\phi},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{f}_{\mathrm{2}i}=\mathrm{cos}i\mathrm{\Delta}\mathit{\phi},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{f}_{\mathrm{3}i}=\mathrm{cos}i\mathrm{\Delta}\mathit{\phi},\\ {\displaystyle}{f}_{\mathrm{4}i}=\mathrm{sin}i\mathrm{\Delta}\mathit{\phi},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{F}_{i}={\left({S}_{m}-a\right)}^{\mathrm{2}}.\end{array}$$

As a result, we get these ratios:

$$\begin{array}{}\text{(10)}& {\displaystyle}\begin{array}{rl}& {\displaystyle \frac{{P}_{\mathrm{1}}}{{P}_{\mathrm{2}}}}=-{\displaystyle \frac{{P}_{\mathrm{4}}}{{P}_{\mathrm{3}}}}={\displaystyle \frac{{S}_{m}-a}{e}}=\mathit{\rho},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{P}_{\mathrm{4}}=-{P}_{\mathrm{2}};\\ & {P}_{\mathrm{1}}={\mathit{\rho}}^{\mathrm{2}}{P}_{\mathrm{3}},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathrm{where}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathit{\rho}=\mathrm{tan}{\mathit{\phi}}_{\mathrm{0}}.\end{array}\end{array}$$

So, four required parameters of the crank press linkage mechanism $r,\phantom{\rule{0.125em}{0ex}}{\mathit{\phi}}_{\mathrm{0}},\phantom{\rule{0.125em}{0ex}}l,\phantom{\rule{0.125em}{0ex}}e$ determined by five parameters are *P*_{0}, *P*_{1}, *P*_{2}, *P*_{3}, and *P*_{4}.

For compatibility condition of the solution, the following constraints must be applied:

$$\begin{array}{}\text{(11)}& {\displaystyle \frac{({S}_{m}-a{)}^{\mathrm{2}}}{(r+l{)}^{\mathrm{2}}}}+{\displaystyle \frac{{e}^{\mathrm{2}}}{{\left(r+l\right)}^{\mathrm{2}}}}=\mathrm{1}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathrm{or}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{\displaystyle \frac{{\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)}^{\mathrm{2}}}{{P}_{\mathrm{3}}^{\mathrm{2}}+{P}_{\mathrm{4}}^{\mathrm{2}}}}=\mathrm{1}+{\displaystyle \frac{{P}_{\mathrm{1}}^{\mathrm{2}}}{{P}_{\mathrm{2}}^{\mathrm{2}}}}\text{(12)}& {\displaystyle}\mathit{\varphi}\phantom{\rule{0.125em}{0ex}}\left({P}_{i}\right)={P}_{\mathrm{2}}^{\mathrm{2}}{\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)}^{\mathrm{2}}-({P}_{\mathrm{3}}^{\mathrm{2}}+{P}_{\mathrm{4}}^{\mathrm{2}})\left({P}_{\mathrm{2}}^{\mathrm{2}}+{P}_{\mathrm{1}}^{\mathrm{2}}\right)=\mathrm{0}.\end{array}$$

According to Eq. (11), there are constraints between unknowns and we can use the method of an objective function of conditional optimization in Eq. (9) (Bertsekas, 1982). Then, by using Eq. (12), the Lagrange function can be obtained as an objective function as follows:

$$\begin{array}{}\text{(13)}& \begin{array}{rl}L& =({P}_{\mathrm{0}}{f}_{\mathrm{0}i}+{P}_{\mathrm{1}}{f}_{\mathrm{1}i}+{P}_{\mathrm{2}}{f}_{\mathrm{2}i}+{P}_{\mathrm{3}}{f}_{\mathrm{3}i}\\ & +{P}_{\mathrm{4}}{f}_{\mathrm{4}i}-{F}_{i}{)}^{\mathrm{2}}-\mathit{\lambda}\mathit{\varphi}\phantom{\rule{0.125em}{0ex}}\left({P}_{i}\right).\end{array}\end{array}$$

For determination of *P*_{0}, *P*_{1}, *P*_{2}, *P*_{3}, *P*_{4}, and *λ*, the
quadratic approximation method, which determines the sum of squared
deviations minimum for *N* mechanism positions, is used (Jomartov and Tuleshov, 2018):

$$\begin{array}{}\text{(14)}& L\underset{min}{\Rightarrow}\mathrm{0}.\end{array}$$

By taking the derivative of *L* with respect to *P*_{j}, the minimum conditions
are given below:

$$\begin{array}{}\text{(15)}& {\displaystyle}\begin{array}{rl}& {\displaystyle \frac{\partial L}{\partial {P}_{j}}}=\mathrm{2}({P}_{\mathrm{0}}{f}_{\mathrm{0}i}+{P}_{\mathrm{1}}{f}_{\mathrm{1}i}+{P}_{\mathrm{2}}{f}_{\mathrm{2}i}+{P}_{\mathrm{3}}{f}_{\mathrm{3}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{P}_{\mathrm{4}}{f}_{\mathrm{4}i}-{F}_{i}){f}_{ji}-\mathit{\lambda}{\displaystyle \frac{\partial \mathit{\varphi}}{\partial {P}_{i}}}=\mathrm{0},\\ & i=-\mathrm{2},-\mathrm{1},\mathrm{0},\mathrm{1},\mathrm{2},\\ & {\displaystyle \frac{\partial L}{\partial \mathit{\lambda}}}=\phantom{\rule{0.125em}{0ex}}{P}_{\mathrm{2}}^{\mathrm{2}}{\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)}^{\mathrm{2}}-({P}_{\mathrm{3}}^{\mathrm{2}}+{P}_{\mathrm{4}}^{\mathrm{2}})\left({P}_{\mathrm{2}}^{\mathrm{2}}+{P}_{\mathrm{1}}^{\mathrm{2}}\right)=\mathrm{0}.\end{array}\end{array}$$

Equation (16) presents a nonlinear equation set related to the required parameters:

$$\begin{array}{}\text{(16)}& {\displaystyle}\begin{array}{rl}& {P}_{\mathrm{0}}\sum {f}_{\mathrm{0}i}^{\mathrm{2}}+{P}_{\mathrm{1}}\sum {f}_{\mathrm{1}i}{f}_{\mathrm{0}i}+{P}_{\mathrm{2}}\sum {f}_{\mathrm{2}i}{f}_{\mathrm{0}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{P}_{\mathrm{3}}\sum {f}_{\mathrm{3}i}{f}_{\mathrm{0}i}+{P}_{\mathrm{4}}\sum {f}_{\mathrm{4}i}{f}_{\mathrm{0}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\phantom{\rule{0.125em}{0ex}}\sum {f}_{\mathrm{0}i}{\left({S}_{m}-a\right)}^{\mathrm{2}},\\ & {P}_{\mathrm{0}}\sum {f}_{\mathrm{0}i}{f}_{\mathrm{1}i}+{P}_{\mathrm{1}}\sum {f}_{\mathrm{1}i}^{\mathrm{2}}+{P}_{\mathrm{2}}\sum {f}_{\mathrm{2}i}{f}_{\mathrm{1}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{P}_{\mathrm{3}}\sum {f}_{\mathrm{3}i}{f}_{\mathrm{1}i}+{P}_{\mathrm{4}}\sum {f}_{\mathrm{4}i}{f}_{\mathrm{1}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\sum {f}_{\mathrm{1}i}{\left({S}_{m}-a\right)}^{\mathrm{2}}+\phantom{\rule{0.125em}{0ex}}\mathrm{2}\mathit{\lambda}[{P}_{\mathrm{2}}^{\mathrm{2}}\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-({P}_{\mathrm{3}}^{\mathrm{2}}+{P}_{\mathrm{4}}^{\mathrm{2}}){P}_{\mathrm{1}}],\\ & {P}_{\mathrm{0}}\sum {f}_{\mathrm{0}i}{f}_{\mathrm{2}i}+{P}_{\mathrm{1}}\sum {f}_{\mathrm{1}i}{f}_{\mathrm{2}i}+{P}_{\mathrm{2}}\sum {f}_{\mathrm{2}i}^{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{P}_{\mathrm{3}}\sum {f}_{\mathrm{3}i}{f}_{\mathrm{2}i}+{P}_{\mathrm{4}}\sum {f}_{\mathrm{4}i}{f}_{\mathrm{2}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\sum {f}_{\mathrm{2}i}{\left({S}_{m}-a\right)}^{\mathrm{2}}+\mathrm{2}\mathit{\lambda}{P}_{\mathrm{2}}{\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)}^{\mathrm{2}},\\ & {P}_{\mathrm{0}}\sum {f}_{\mathrm{0}i}{f}_{\mathrm{3}i}+{P}_{\mathrm{1}}\sum {f}_{\mathrm{1}i}{f}_{\mathrm{3}i}+{P}_{\mathrm{2}}\sum {f}_{\mathrm{2}i}{f}_{\mathrm{3}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{P}_{\mathrm{3}}\sum {f}_{\mathrm{3}i}^{\mathrm{2}}+{P}_{\mathrm{4}}\sum {f}_{\mathrm{4}i}{f}_{\mathrm{3}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\phantom{\rule{0.125em}{0ex}}\sum {f}_{\mathrm{3}i}{\left({S}_{m}-a\right)}^{\mathrm{2}}+\phantom{\rule{0.125em}{0ex}}\mathrm{2}\mathit{\lambda}[{P}_{\mathrm{2}}^{\mathrm{2}}\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-({P}_{\mathrm{2}}^{\mathrm{2}}+{P}_{\mathrm{1}}^{\mathrm{2}}){P}_{\mathrm{3}}],\\ & {P}_{\mathrm{0}}\sum {f}_{\mathrm{0}i}{f}_{\mathrm{4}i}+{P}_{\mathrm{1}}\sum {f}_{\mathrm{1}i}{f}_{\mathrm{4}i}+{P}_{\mathrm{2}}\sum {f}_{\mathrm{2}i}{f}_{\mathrm{4}i}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{P}_{\mathrm{3}}\sum {f}_{\mathrm{3}i}{f}_{\mathrm{4}i}+{P}_{\mathrm{4}}\sum {f}_{\mathrm{4}i}^{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\phantom{\rule{0.125em}{0ex}}\sum {f}_{\mathrm{4}i}{\left({S}_{m}-a\right)}^{\mathrm{2}}-\phantom{\rule{0.125em}{0ex}}\mathrm{2}\mathit{\lambda}({P}_{\mathrm{2}}^{\mathrm{2}}+{P}_{\mathrm{1}}^{\mathrm{2}}){P}_{\mathrm{4}},\\ & {P}_{\mathrm{2}}^{\mathrm{2}}{\left({P}_{\mathrm{1}}+{P}_{\mathrm{3}}\right)}^{\mathrm{2}}-({P}_{\mathrm{3}}^{\mathrm{2}}+{P}_{\mathrm{4}}^{\mathrm{2}})\left({P}_{\mathrm{2}}^{\mathrm{2}}+{P}_{\mathrm{1}}^{\mathrm{2}}\right)=\mathrm{0}.\end{array}\end{array}$$

The nonlinear set (16) has to satisfy minimum conditions (15). It is solved
in relation to ${P}_{\mathrm{0}},{P}_{\mathrm{1}},{P}_{\mathrm{2}},{P}_{\mathrm{3}},{P}_{\mathrm{4}}$, and *λ*. From the solution
set, the equation of ${P}_{\mathrm{4}}=-{P}_{\mathrm{2}}$ is selected.

The required geometric parameters are found by using the following equations:

$$\begin{array}{}\text{(17)}& {\displaystyle}\begin{array}{rl}& \mathrm{tan}{\mathit{\phi}}_{\mathrm{0}}={\displaystyle \frac{{P}_{\mathrm{1}}}{{P}_{\mathrm{2}}}};\\ & e={\displaystyle \frac{{S}_{m}-a}{\mathrm{tan}{\mathit{\phi}}_{\mathrm{0}}}};\\ & r={\displaystyle \frac{-{P}_{\mathrm{3}}}{\mathrm{2}e\mathrm{cos}{\mathit{\phi}}_{\mathrm{0}}}};\\ & l=\sqrt{{e}^{\mathrm{2}}+{r}^{\mathrm{2}}-{P}_{\mathrm{0}}}.\end{array}\end{array}$$

4 Results and discussion

The mechanism's geometric parameters were obtained by using the proposed
Stephenson II synthesis method. The Eq. (16) was solved for the
mechanism's five positions, and the geometric parameters of the mechanism
were determined by using Eq. (17). As a result, the mechanism
dimensions that could meet the desired conditions were found with the
proposed method. Table 1 shows the mechanism dimensions of the Stephenson
crank press synthesized in this study. The kinematic scheme of the mechanism
is presented in Fig. 3a, with the size of links: *r*=60 mm, *φ*_{0}=260^{∘}, *l*=160 mm, $e=-\mathrm{40}$ mm, when set ${S}_{m}-a=\mathrm{216.65}$ mm and *a*=41 mm. By using these synthesized dimensions, the effects of
dimension *e* (eccentricity) were studied additionally.

Numerical modeling programs were developed by the
MAPLE^{®} software. The kinematic results (position,
speed, and acceleration results) were verified by
ASIAN-2014^{®} software.

Figure 3b presents a graph of slider movement. After crank angle passes
200^{∘}, a movement close to dwell is seen in the interval of
$\mathit{\phi}\in (\mathrm{251}{}^{\circ},\phantom{\rule{0.125em}{0ex}}\mathrm{271}{}^{\circ})$ with an accuracy of
Δ=0.063. Such accuracy was obtained by solving the
nonlinear set (Eq. 16) in the first approximation.

Figure 4 shows the three different eccentricity designs with *e* dimensions of
−40, 0, and +40, respectively. The crank moves with a constant velocity,
and Fig. 5 shows the crank angle versus time. This input was used for each
simulation. Figure 6 presents the slider position versus time that changes
with different eccentricity dimensions of *e*.

The bottom dead center of the slider was set to zero position. Figure 7 shows graphs of the set position, velocity, and acceleration of the working slider. It is observed that changing of eccentricity affects the kinematics of the Stephenson mechanism. The eccentricity of −40 provides slower forwarding and faster backwarding. While the similar maximum slider points are shown for each eccentricity, their start and end points are different as given in Fig. 7a.

An important parameter of the mechanism is the coefficient of increase of the average slider velocity:

$$\begin{array}{}\text{(18)}& K={\displaystyle \frac{{\mathit{\phi}}_{m}}{\mathrm{2}\mathit{\pi}-{\mathit{\phi}}_{m}}},\end{array}$$

where *φ*_{m} is crank rotation angle during the slider stroke; the
crank rotates counterclockwise.

The velocity graph analysis in Fig. 7b allows determining the *φ*_{m} values at different eccentricity values: ${e}_{\mathrm{1}}=-\mathrm{40}$ mm, *e*_{2}=0, *e*_{3}=40 mm. They are equal to ${\mathit{\phi}}_{{m}_{\mathrm{1}}}=\mathrm{194.8}$^{∘}, ${\mathit{\phi}}_{{m}_{\mathrm{2}}}=\mathrm{180}$^{∘}, ${\mathit{\phi}}_{{m}_{\mathrm{3}}}=\mathrm{165.2}$^{∘}. When the values are substituted in Eq. (18), the
average slider velocity coefficient values corresponding to eccentricity
values are found as *K*_{1}=1.18; *K*_{2}=1.0; *K*_{3}=0.847.

A slow lowering in the working stroke and fast lifting in the reverse stroke
of the slider is significant in press tools. As shown in figures, this
process is realized when ${e}_{\mathrm{1}}=-\mathrm{40}$ mm and *K*_{1}=1.18. The graph
analysis in Fig. 7c shows that the slider acceleration at the forward
stroke beginning is almost unchanged and noticeably accelerates at the end
of the working process. This is common in all cases. The maximum velocity of
47 mm s^{−1} is reached at a maximum acceleration of 33 mm s^{−2} in 2 s, which is the cause of the shock load between the slider and the
forging.

5 Conclusions

In this study, an analytical method for kinematic synthesis of linkage mechanisms based on the mean-square minimization of the objective function was developed. By taking into account the structural features of the mechanism, a method for deriving this objective function was proposed. Based on this method, algorithms were obtained, and Stephenson's press mechanism was synthesized. The crank length, two parallel connecting rods length, and the eccentricity of the slider relative to the axis of the crank rotation were found in an appropriate way. Eccentricity effects were also compared. The approximate slider dwell movement and the regulated coefficient of increase of the average slider velocity of the Stephenson II mechanism were achieved. It was seen that the coefficient of increase of the average slider velocity allowed regulating the technological process character and providing the desired cycle of the slider movement.

Data availability

Data are available from the authors.

Author contributions

AT, RH, AS, and MK discussed and decided on the methodology in the study and prepared the manuscript. AT conceived and developed the kinematic synthesis of the mechanism. Kinematic modeling and analysis were performed by RH, who is the corresponding author. AS participated in the kinematic synthesis and analysis. MK participated in developing the kinematic scheme and coordinated helped to draft the manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

The authors would like to thank the Ministry of Education and Science of the Republic of Kazakhstan, the U. A. Joldasbekov Institute of Mechanics and Mechanical Engineering, and LLP “Semizbay-U” for their support.

Financial support

This research has been supported by LLP “Semizbay-U” (agreement no. 478227-2/2020/1).

Review statement

This paper was edited by Daniel Condurache and reviewed by three anonymous referees.

References

Bertsekas, D.: Constrained Optimization and Lagrange Multiplier Methods, Elsevier, Academic Press, USA, https://doi.org/10.1016/C2013-0-10366-2, 1982.

Bocharev, I. A.: Press forging equipment, Akademiia, Moscow, 2008 (in Russian).

Connor, A. M., Douglas, S. S., and Gilmartin, M. J.: The synthesis of hybrid five-bar path generating mechanisms using genetic algorithms, 1st International Conference on Genetic Algorithms in Engineering Systems: Innovations and Applications, 313–318, Galesia, 1995.

Drakunov, Y., Tuleshov, A., Jomartov, A., Tuleshov, E. A., Seidakhmet, A. Zh., Kayim, T. T., Ibrayev, S. M., Jamalov, N. K., Drakunov, A. Yu., and Drakunov, S. U.: Development of methods and technology of designing of power press-machines based on new crank operating mechanisms, National Center of Science and Technology Evaluation, Report No. 0218PK00926, 1–160, 2018 (in Russian).

Dülger, L. and Kireçci, A.: A study on a hybrid actuator, Mech. Mach. Theory, 35, 1141–1149, 1999.

Dülger, L. C. and Kireçci, A.: Hybrid Manipulator Modelling and Simulation, 7th National Machine Thery Symposium, YTU, Istanbul, 1995 (in Turkish).

Dülger, L. C., Kireçci, A., and Topalbekiroglu, M.: Modeling and simulation of a hybrid actuator, Mech. Mach. Theory, 38, 395–407, 2003.

Huang, P. and Ding, H.: Structural synthesis of Assur groups with up to 12 links and creation of their classified databases, Mech. Mach. Theory, 145, 103668, https://doi.org/10.1016/j.mechmachtheory.2019.103668, 2020.

Jomartov, A. and Tuleshov, A.: Vector method for kinetostatic analysis of planar linkages, J. Braz. Soc. Mech. Sci., 40, 1–9, https://doi.org/10.1007/s40430-018-1022-y, 2018.

Jomartov, A. and Tuleshov, A.: Modeling Dynamics of Planetary Gears of Crank Press on SimulationX, in: Advances in Italian Mechanism Science, edited by: Carbone, G. and Gasparetto, A., IFToMM ITALY 2018, Mechanisms and Machine Science, Vol. 68, Springer, Chamber, 41–48, https://doi.org/10.1007/978-3-030-03320-0_5, 2019.

Jomartov, A., Tuleshov, A., and Kuatova, M.: Designing of a Crank press on the Basis of High Class Planar Linkages, in: Advances in Mechanism and Machine Science, edited by: Uhl, T., IFToMM WC 2019. Mechanisms and Machine Science, Vol. 73, Springer, Chamber, 3027–3036, https://doi.org/10.1007/978-3-030-20131-9_298, 2019.

Jomartov, A., Tuleshov, A., Jamalov, N., Kuatova, M., and Kaimov, A.: Designing of the Stephenson II six-link linkage actuator for servo mechanical press, International Journal of Mechanical and Production Engineering Research and Development, 10, 501–512, https://doi.org/10.24247/ijmperdapr202053, 2020.

Kütük, M. E. and Dülger, L. C.: A hybrid press system: Motion design and inverse kinematics issues, Engineering Science and Technology, an International Journal, 19, 846–856, 2016.

Li, C. H. and Tso, P. L.: Experimental study on a hybrid driven press using iterative learning control, International Journal of Machine Tools and Manufacture, 48, 209–219, 2008.

Li, H. and Zhang, Y.: Seven bar mechanical press with hybrid driven mechanism for deep drawing; part 1: kinematics analysis and optimum design, part 2: dynamic modeling and simulation, J. Mech. Sci. Technol., 11, 2153–2160, 2161–2167, 2010.

Ouyang, P. R., Li , Q., and Zhang, W. J.: Design, modelling and control of a hybrid machine system, Mechatronics, 14, 1197–1217, 2004.

Tokuz, L. C.: Hybrid machine modelling and control, PhD thesis, Liverpool Polytechnic, 1992.

Tso, P. L.: Optimal Design of a Hybrid-Driven Servo Press and Experimental Verification, J. Mech. Design, 132, 034503, https://doi.org/10.1115/1.4000213, 2010.

Tso, P. L. and Li, C. H.: Study of servo press with a flywheel, J. Adv. Mech. Des. Syst., 2, 1–11, 2008.

Tuleshov, A. K., Jomartov, A. A., and Kuatova M. Zh.: The Model of The Movement of The Crank Press Based on The Lever Mechanism of the 4th Class, Annotations of reports of the 12th All-Russian Congress on fundamental problems of theoretical and applied mechanics, Russian Federation, Ufa, 2019 (in Russian).

Yevgrafov, M. Z., Semyonov, Y. A., and Slousch, A. V.: Theory of mechanisms and machines: a training manual, Academia, Moscow, 2006 (in Russian).

Yu, H.: Modeling and control of hybrid machine systems – a five-bar mechanism case, International Journal of Automation and Computing, 3, 235–243, https://doi.org/10.1007/s11633-006-0235-1, 2006.

Yuan, Z., Gilmartin, M.J., and Douglas, S.S.: Design of hybrid machines for nonuniform motion production, Proceedings of the Institution of Mechanical Engineers Pt. C, 219, 491–499, https://doi.org/10.1243/095440605X16992, 2005.

Zhang, K.: Optimization dynamics design of hybrid driving mechanisms, Proceedings of the 2006 IEEE International Conference on Mechatronics and Automation, 1914–1919, 2006.

Mechanical Sciences

An open-access journal for theoretical and applied mechanics