Next Article in Journal
Dual Band Antenna Design and Prediction of Resonance Frequency Using Machine Learning Approaches
Previous Article in Journal
Systematic Machine Translation of Social Network Data Privacy Policies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dimensional Synthesis of Watt II and Stephenson III Six-Bar Slider-Crank Function Generators for Nine Prescribed Positions

1
Department of Mechanical Engineering, Graduate School, Korea University, Seoul 02841, Korea
2
School of Mechanical Engineering, Korea University, Seoul 02841, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(20), 10503; https://doi.org/10.3390/app122010503
Submission received: 28 September 2022 / Revised: 12 October 2022 / Accepted: 14 October 2022 / Published: 18 October 2022
(This article belongs to the Section Mechanical Engineering)

Abstract

:
This paper proposes an efficient exact dimensional synthesis method for finding all the link lengths of the Watt II and Stephenson III six-bar slider-crank function generators, satisfying nine prescribed precision points using the homotopy continuation method. The synthesis equations of each mechanism are initially constructed as a system of 56 quadratic polynomials whose Bézout number, which represents the maximum number of solutions, is 256 ≅ 7.21 × 1016. In order to reduce the size of the system, multi-homogeneous formulation is applied to transform the system into 12 equations in 12 unknowns, and the multi-homogeneous Bézout number of the system is 286,720. The Bertini solver, based on the homotopy continuation method, is used to solve the synthesis equations to obtain the dimensions of the two mechanisms. For the arbitrarily given nine precision points, the proposed method yields 37 and 31 defect-free solutions of Watt II and Stephenson III six-bar slider-crank mechanisms, respectively, and it is confirmed that they pass through the prescribed positions.

1. Introduction

The dimensional synthesis of a mechanism determines the dimensions of the individual links that generate a desired motion. Since the development of an analytical approach to the design of four-bar linkages introduced by Freudenstein [1], numerous studies have been conducted on the dimensional synthesis, which can be classified into two categories. One is the exact synthesis method, known as the precision point approach, in which a finite number of exact positions are to be satisfied. The other is the approximate synthesis, which seeks the dimensions that approximately satisfy prescribed positions.
In general, dimensional synthesis problems are formulated as a system of nonlinear equations and require heavy computation. For this reason, various approximate synthesis approaches using numerical methods [2,3], optimization methods such as least square techniques [4,5,6], genetic algorithms [7,8,9,10], and evolutionary algorithms [11,12,13,14,15] have been used widely in recent times. However, the approximate synthesis methods require appropriate initial values that are close to the actual solutions of the synthesis equations. When unsuitable initial values are assumed, the solutions may not be found. In addition, the methods may not find all the solutions of the synthesis equations.
In contrast, since the exact synthesis method can find all the solutions without estimating an initial value, it yields diverse mechanism dimensions that satisfy a desired motion. Among the studies on the precision point approach, Huang et al. [16] and Almandeel et al. [17] used elimination methods to perform the dimensional synthesis of four-bar path generators for five precision points and slider-crank function generators for four precision points, respectively. Subbian and Flugrad [18] and Wampler et al. [19] carried out the dimensional synthesis of four-bar path generators for five precision points and nine precision points, respectively. For the dimensional synthesis of six-bar linkages with revolute joints only, Dhingra et al. [20] and Plecnik and McCarthy [21,22,23] applied the homotopy continuation method for the synthesis of Watt II, Stephenson II, and Stephenson III six-link function generators, proposing design methods for eight, nine, and eleven precision points, respectively. It is worth noting that the continuation method is developed from the parameter perturbation procedure introduced by Freudenstein and Roth [24] for solving a system of nonlinear equations starting from a known solution of a similar system.
Most studies on the precision point approach deal with four-bar linkages, four-link slider-crank mechanisms, and six-bar linkages with revolute joints only, but the exact synthesis of six-link slider-crank mechanisms with six revolute joints and one prismatic joint are extremely rare. Since the six-bar slider-crank linkage has more design parameters, it can satisfy more precision points than its four-link counterpart. Different types of six-bar slider-crank linkages are widely used in various applications such as ore crushers, mechanical presses [25], punching mechanisms [26], and quick return mechanisms [27] in machinery. These precision mechanisms are required to satisfy precise motion between input and output links. Furthermore, by adding one link and one joint, the six-link slider-crank can be modified to a two-degree-of-freedom seven-link adjustable mechanism that can be used in controllable presses for punching and metal forming [28,29], variable pumps [30], and variable compression ratio engines [31,32,33,34]. For the design of this machinery, the dimensional synthesis of the six-link slider-crank is essential.
This paper proposes a systematic exact dimensional synthesis method to find all Watt II and Stephenson III six-bar slider-crank function generators that satisfy nine precision points of the output slider with respect to prescribed input crank angles. The synthesis equations of each mechanism for nine precision points are derived as a system of 56 quadratic polynomial equations with 58 unknowns, respectively. By assuming two unknowns as free choices, the system can be solved. However, since the system size is very large, it is necessary to decrease its size so as to reduce the calculation time for pursuing an efficient dimensional synthesis. To this end, some of the unknowns are eliminated from the synthesis equations to simplify the system to eight equations in eight unknowns. Then, in order to apply the multi-homogeneous Bézout theorem [35,36], four auxiliary equations that render the system a set of homogeneous polynomials are defined, and the unknowns are divided into appropriate two-homogeneous groups. Through this process, the synthesis equations for Watt II and Stephenson III six-bar slider-crank linkages are formulated into eight quartic polynomials and four auxiliary quadratic polynomials in 12 unknowns, and the multi-homogeneous Bézout number of the system that represents the maximum number of solutions is 286,720. The synthesis equations for the two mechanisms are solved by using the Bertini software program [37] that is based on the homotopy continuation method. The computation results provided 25,630 and 36,061 nonsingular solutions, respectively. Among the nonsingular solutions, mechanisms that allow full rotation of the input crank and satisfy nine precision points on a single stroke of the slider are finally selected as feasible solutions. As a result, 37 Watt II and 31 Stephenson III six-bar slider-crank linkages without kinematic defects are determined.
The paper is arranged as follows. In Section 2, the dimensional synthesis equations for Watt II and Stephenson III six-bar slider-crank function generators are derived and a process for simplifying the synthesis equations is presented. Section 3 describes the procedure for calculating the multi-homogeneous Bézout number by applying the multi-homogeneous Bézout theorem. Section 4 then explains how to select feasible solutions that enable prescribed precision points to lie on a single stroke as well as the input link to rotate continuously. Section 5 presents examples of applying the dimensional synthesis process proposed in this paper. Finally, Section 6 concludes with a summary and discussion of the research.

2. Dimensional Synthesis Equations for Six-Bar Slider-Crank Mechanisms

Each of the dimensional synthesis equations for six-bar slider-crank function generators can be formulated by a system of polynomial equations. As the number of prescribed positions increases, the system size becomes very large and the calculation time to obtain the solution could increase tremendously. Hence, in order to carry out the dimensional synthesis efficiently, it is necessary to reduce the size of the equation system. In this section, the synthesis equations for the two mechanisms are derived and simplified by eliminating some unknowns.

2.1. Synthesis Equations of Watt II Six-Bar Slider-Crank Function Generator

A Watt II six-bar slider-crank linkage is shown in Figure 1a, in which the input link r1 rotates about the fixed pivot O, and the output slider translates vertically. A fixed frame is positioned so that its origin is at O and oriented so that its x-axis is perpendicular to the direction of slide. The horizontal distance from O to the line of path of the moving pivot E is the offset denoted by h. In the fixed frame, the input angle θ is measured from the positive x-axis and the output displacement of the slider p is the vertical distance from the x-axis to E, the center of the revolute joint attached to the slider. The angle of link r0 is η, and the angle between links r3 and r4 is β. The angular displacements of links r2, r3, r4, and r5 are denoted by ϕ, α, α + β, and δ, respectively.
Figure 1b shows the mechanism in its first and jth positions. In the Figure, Δθj, Δϕj, Δαj, and Δδj denote the rotation angles of links r1 to r5, respectively, and Δpj represents the displacement of the slider, all measured from the initial position to jth position of the mechanism. Let θ1 and p1 denote the angle of input link r1 and the position of output slider in the initial position, respectively. Then, the prescribed positions of the function generator to be synthesized are given pairwise as
Δ θ j = θ j θ 1 , Δ p j = p j p 1 ,   j = 2 , , n ,
where θj and pj are the input angle and the position of the slider in jth position, respectively, and n is the number of precision points to be satisfied.
The dimensional synthesis equations of the Watt II six-bar slider-crank function generator can be constructed by the closure equations of the two loops OABC and OCDE shown in Figure 1 as follows.
Loop OABC:
For the initial and jth positions, the loop closure equations can be written as
r 1 + r 2 = r 0 + r 3 ,
and
r 1 e i Δ θ j + r 2 e i Δ ϕ j = r 0 + r 3 e i Δ α j ,   j = 2 , , n ,
respectively. Subtracting Equation (2) from Equation (3) gives
r 1 ( e i Δ θ j 1 ) + r 2 ( e i Δ ϕ j 1 ) = r 3 ( e i Δ α j 1 ) ,   j = 2 , , n .
Loop OCDE:
For the initial and jth positions, the loop closure equations can be written as
r 0 + r 4 = h + i p 1 + r 5 ,
and
r 0 + r 4 e i Δ α j = h + i p j + r 5 e i Δ δ j ,   j = 2 , , n ,
respectively. Subtracting Equation (5) from Equation (6) gives
r 4 ( e i Δ α j 1 ) = i Δ p j + r 5 ( e i Δ δ j 1 ) ,   j = 2 , , n .
Separating the real and imaginary parts of Equations (4) and (7) leads to
r 1 { cos ( θ 1 + Δ θ j ) cos θ 1 } + r 2 { cos ( ϕ 1 + Δ ϕ j ) cos ϕ 1 } = r 3 { cos ( α 1 + Δ α j ) cos α 1 } ,   j = 2 , , n ,
r 1 { sin ( θ 1 + Δ θ j ) sin θ 1 } + r 2 { sin ( ϕ 1 + Δ ϕ j ) sin ϕ 1 } = r 3 { sin ( α 1 + Δ α j ) sin α 1 } ,   j = 2 , , n ,
and
r 4 { cos ( α 1 + β + Δ α j ) cos ( α 1 + β ) } = r 5 { cos ( δ 1 + Δ δ j ) cos δ 1 } ,   j = 2 , , n ,
r 4 { sin ( α 1 + β + Δ α j ) sin ( α 1 + β ) } = r 5 { sin ( δ 1 + Δ δ j ) sin δ 1 } + Δ p j ,   j = 2 , , n ,
respectively. Note that θ1, ϕ1, α1, and δ1 denote the angles of the corresponding links shown in Figure 1a in the initial position. Among the unknowns in Equations (8) through (11), ∆ϕj, ∆αj, and ∆δj should satisfy trigonometric identities
cos 2 Δ ϕ j + sin 2 Δ ϕ j = 1 , cos 2 Δ α j + sin 2 Δ α j = 1 , cos 2 Δ δ j + sin 2 Δ δ j = 1 ,
where j = 2, …, n.
The unknowns in the synthesis equations, Equations (8)–(12), are r1, r2, r3, r4, r5, θ1, ϕ1, α1, β, δ1, cos∆ϕj, sin∆ϕj, cos∆αj, sin∆αj, cos∆δj, and sin∆δj, where j = 2, …, n. Hence, the problem of the dimensional synthesis for the function generator that satisfies n prescribed positions is formulated as 7(n − 1) equations in 10 + 6(n − 1) unknowns. For 9 prescribed positions, n is equal to 9 and the synthesis equations consist of 56 quadratic polynomials with respect to 58 unknowns. Therefore, to match the number of unknowns and the number of equations, two free choices can be made. The total degree of a system of m polynomial equations, which is referred to as the Bézout number, is defined as i = 1 m d i , where di is the degree of the ith polynomial [35]. Hence, the total degree of the synthesis equations for 9 precision points is 256 ≅ 7.21 × 1016. Since the system size is too large to solve, it is necessary to reduce the size of the synthesis equations as explained in the following sections.

2.2. Procedure for Simplifying the Synthesis Equations of Watt II Six-Bar Slider-Crank

The system of 7(n − 1) synthesis equations, Equation (8) through (12), can be reduced to a system of n−1 polynomial equations by eliminating the unknowns ∆ϕj, ∆δj, and ∆αj, which corresponds to the elimination of 6(n − 1) unknowns cos∆ϕj, sin∆ϕj, cos∆δj, sin∆δj, cos∆αj, and sin∆αj, where j = 2, …, n. For the elimination of ∆ϕj from Equations (8) and (9), isolate the term r2cos(ϕ1 + ∆ϕj) in Equation (8) and r2sin(ϕ1 + ∆ϕj) in Equation (9) to one side. Then, square both sides of the equations, add, and simplify the resulting equation using the trigonometric identity cos2(ϕ1 + ∆ϕj)+sin2(ϕ1 + ∆ϕj) = 1 to obtain
r 2 2 = [ r 3 { cos ( α 1 + Δ α j ) cos α 1 } r 1 { cos ( θ 1 + Δ θ j ) cos θ 1 } + r 2 cos ϕ 1 ] 2 + [ r 3 { sin ( α 1 + Δ α j ) sin θ 3 } r 1 { sin ( θ 1 + Δ θ j ) sin θ 1 } + r 2 sin ϕ 1 ] 2 ,   j = 2 , , n .
Similarly, the unknown ∆δj can be eliminated from Equations (10) and (11) as
r 5 2 = [ r 4 { cos ( α 1 + β + Δ α j ) cos ( α 1 + β ) } + r 5 cos δ 1 ] 2 + [ r 4 { sin ( α 1 + β + Δ α j ) sin ( α 1 + β ) } + r 5 sin δ 1 Δ p j ] 2 ,   j = 2 , , n .
Let the x- and y-components of the initial positions of the links r1, r2, r3, r4, and r5 be
r 1 cos θ 1 = r 1 x , r 2 cos ϕ 1 = r 2 x , r 3 cos α 1 = r 3 x , r 4 cos ( α 1 + β ) = r 4 x , r 5 cos δ 1 = r 5 x , r 1 sin θ 1 = r 1 y , r 2 sin ϕ 1 = r 2 y , r 3 sin α 1 = r 3 y , r 4 sin ( α 1 + β ) = r 4 y , r 5 sin δ 1 = r 5 y .
Substituting Equation (15) into Equation (13) and gathering the coefficients of cos∆αj and sin∆αj gives
L 1 j + L 2 j cos Δ α j + L 3 j sin Δ α j = 0 ,   j = 2 , , n ,
where
L 1 j = 2 cos Δ θ j ( r 1 x 2 + r 1 y 2 + r 1 x r 2 x r 1 x r 3 x + r 1 y r 2 y r 1 y r 3 y ) + 2 sin Δ θ j ( r 1 x r 2 y + r 1 y r 2 x + r 1 x r 3 y r 1 y r 3 x ) + 2 ( r 3 x 2 + r 3 y 2 r 2 x r 3 x r 2 y r 3 y ) + 2 ( r 1 x 2 + r 1 y 2 + r 1 x r 2 x r 1 x r 3 x + r 1 y r 2 y r 1 y r 3 y ) ,
L 2 j = 2 cos Δ θ j ( r 1 x r 3 x + r 1 y r 3 y ) 2 sin Δ θ j ( r 1 x r 3 y r 1 y r 3 x ) + 2 ( r 1 x r 3 x + r 1 y r 3 y ) 2 ( r 3 x 2 + r 3 y 2 r 2 x r 3 x r 2 y r 3 y ) ,
L 3 j = 2 cos Δ θ j ( r 1 x r 3 y r 1 y r 3 x ) 2 sin Δ θ j ( r 1 x r 3 x + r 1 y r 3 y ) 2 ( r 1 x r 3 y r 1 y r 3 x ) 2 ( r 2 x r 3 y r 2 y r 3 x ) .
Similarly, substituting Equation (15) into Equation (14) and combining the coefficients of cos∆αj and sin∆αj yields
Q 1 j + Q 2 j cos Δ α j + Q 3 j sin Δ α j = 0 ,   j = 2 , , n ,
where
Q 1 j = 2 ( r 4 x 2 + r 4 y 2 r 4 x r 5 x r 4 y r 5 y ) + Δ p j 2 + 2 Δ p j ( r 4 y r 5 y ) ,
Q 2 j = 2 ( r 4 x 2 + r 4 y 2 r 4 x r 5 x r 4 y r 5 y ) 2 r 4 y Δ p j ,
Q 3 j = 2 ( r 4 x r 5 y r 4 y r 5 x ) 2 r 4 x Δ p j .
In order to eliminate ∆αj from Equations (16) and (17), solve the two equations for cos∆αj and sin∆αj to obtain
cos Δ α j = L 1 j Q 3 j + L 3 j Q 1 j L 2 j Q 3 j L 3 j Q 2 j ,   sin Δ α j = L 1 j Q 2 j L 2 j Q 1 j L 2 j Q 3 j L 3 j Q 2 j ,   j = 2 , , n .
Squaring both sides of Equation (18), adding the result, and simplifying using the trigonometric identity cos2(∆αj) + sin2(∆αj) = 1 gives
( L 1 j Q 3 j + L 3 j Q 1 j ) 2 + ( L 1 j Q 2 j L 2 j Q 1 j ) 2 ( L 2 j Q 3 j L 3 j Q 2 j ) 2 = 0 ,   j = 2 , , n .
Equation (19), the system of synthesis equations for n prescribed positions, has n − 1 eighth-degree polynomials in 10 unknowns, rix and riy, i = 1, …, 5, given in Equation (15). Hence, for 9 position synthesis, n is equal to 9 and there are 8 equations in 10 unknowns, and two unknowns can be selected as free choices. Among the 10 unknowns, r1x and r1y are chosen as the free choices in this study to reduce the calculation time as explained in Section 3.1. Now, the total degree of the 9 precision point synthesis equations is 88 = 16,777,216, which is much smaller than that of the synthesis equations derived in Section 2.1. However, the system size of the simplified synthesis equations can be further reduced by using auxiliary equations and applying the multi-homogeneous Bézout theorem. This process is described in Section 3.

2.3. Synthesis Equations of Stephenson III Six-Bar Slider-Crank Function Generator

A Stephenson III slider-crank linkage is shown in Figure 2a. The input link r1 rotates about the fixed pivot O, and the output slider translates along a line parallel to y-axis. The angle between links r2 and r4 is β. The angular displacements of links r2, r3, r4, and r5 are denoted by ϕ, α, ϕ + β, and δ, respectively. As shown in Figure 2b, Δθj, Δϕj, Δαj, and Δδj represent the rotation angles of links r1 to r5, measured from their initial positions to jth positions, respectively. The prescribed precision points are given as Equation (1).
The dimensional synthesis equations of the Stephenson III six-bar slider-crank function generator can be derived using the closure equations of the two loops, OABC and OADE, shown in Figure 2 as follows.
Loop OABC:
For the initial and jth positions, the loop closure equations can be written as
r 1 + r 2 = r 0 + r 3 ,
and
r 1 e i Δ θ j + r 2 e i Δ ϕ j = r 0 + r 3 e i Δ α j ,   j = 2 , , n ,
respectively. Subtracting Equation (20) from Equation (21) gives
r 1 ( e i Δ θ j 1 ) + r 2 ( e i Δ ϕ j 1 ) = r 3 ( e i Δ α j 1 ) ,   j = 2 , , n .
Loop OADE:
For the initial and jth positions, the loop closure equations can be written as
r 1 + r 4 = h + i p 1 + r 5 ,
and
r 1 e i Δ θ j + r 4 e i Δ α j = h + i p j + r 5 e i Δ δ j ,   j = 2 , , n ,
respectively. Subtracting Equation (23) from Equation (24) gives
r 1 ( e i Δ θ j 1 ) + r 4 ( e i Δ α j 1 ) = i Δ p j + r 5 ( e i Δ δ j 1 ) ,   j = 2 , , n .
Separating the real and imaginary parts of Equations (22) and (25) yields
r 1 { cos ( θ 1 + Δ θ j ) cos θ 1 } + r 2 { cos ( ϕ 1 + Δ ϕ j ) cos ϕ 1 } = r 3 { cos ( α 1 + Δ α j ) cos α 1 } ,   j = 2 , , n ,
r 1 { sin ( θ 1 + Δ θ j ) sin θ 1 } + r 2 { sin ( ϕ 1 + Δ ϕ j ) sin ϕ 1 } = r 3 { sin ( α 1 + Δ α j ) sin α 1 } ,   j = 2 , , n ,
and
r 1 { cos ( θ 1 + Δ θ j ) cos θ 1 } + r 4 { cos ( ϕ 1 + β + Δ α j ) cos ( ϕ 1 + β ) } = r 5 { cos ( δ 1 + Δ δ j ) cos δ 1 } ,   j = 2 , , n ,
r 1 { sin ( θ 1 + Δ θ j ) sin θ 1 } + r 4 { sin ( ϕ 1 + β + Δ α j ) sin ( ϕ 1 + β ) } = r 5 { sin ( δ 1 + Δ δ j ) sin δ 1 } + Δ p j ,   j = 2 , , n ,
respectively. Among the unknowns in Equations (26) through (29), ∆ϕj, ∆αj, and ∆δj should satisfy trigonometric identities
cos 2 Δ ϕ j + sin 2 Δ ϕ j = 1 , cos 2 Δ α j + sin 2 Δ α j = 1 , cos 2 Δ δ j + sin 2 Δ δ j = 1 .
where j = 2, …, n.
As in case of the Watt II slider-crank mechanism, the dimensional synthesis equations, Equations (26)–(30), for the Stephenson III function generator that satisfies n prescribed positions are formulated as a system of 7(n − 1) equations in 10 + 6(n − 1) unknowns which are r1, r2, r3, r4, r5, θ1, ϕ1, α1, β, δ1, cos∆ϕj, sin∆ϕj, cos∆αj, sin∆αj, cos∆δj, and sin∆δj, with j = 2, …, n. For 9 prescribed positions, n is equal to 9 and there are 56 equations in 58 unknowns. Hence, two unknowns need to be assumed as free choices in order to obtain the solutions. The total degree of the synthesis equations is 256 ≅ 7.21 × 1016. In the following section, the synthesis equations for n precision points are simplified by eliminating some of the unknowns.

2.4. Procedure for Simplifying the Synthesis Equations of Stephenson III Six-Bar Slider-Crank

The synthesis equations, Equations (26)–(30), can be reduced from 7(n − 1) to n − 1 by eliminating ∆αj, ∆δj, and ∆ϕj, j = 2, …, n. Applying the similar procedure used in Section 2.2, the unknown ∆αj can be eliminated from Equations (26) and (27) as
r 3 2 = [ r 1 { cos ( θ 1 + Δ θ j ) cos θ 1 } + r 2 { cos ( ϕ 1 + Δ ϕ   j ) cos ϕ 1 } + r 3 cos α 1 ] 2 + [ r 1 { sin ( θ 1 + Δ θ j ) sin θ 1 } + r 2 { sin ( ϕ 1 + Δ ϕ j ) sin ϕ 1 } + r 3 sin α 1 ] 2 ,   j = 2 , , n ,
and the unknown ∆δj can be eliminated from Equations (28) and (29) as
r 5 2 = [ r 1 { cos ( θ 1 + Δ θ j ) cos θ 1 } + r 4 { cos ( ϕ 1 + β + Δ ϕ j ) cos ( ϕ 1 + β ) } + r 5 cos δ 1 ] 2 + [ r 1 { sin ( θ 1 + Δ θ j ) sin θ 1 } + r 4 { sin ( ϕ 1 + β + Δ ϕ j ) sin ( ϕ 1 + β ) } + r 5 sin δ 1 Δ p j ] 2 ,   j = 2 , , n .
Let the x- and y-components of the initial positions of the links r1, r2, r3, r4, and r5 be
r 1 cos θ 1 = r 1 x , r 2 cos ϕ 1 = r 2 x , r 3 cos α 1 = r 3 x , r 4 cos ( ϕ 1 + β ) = r 4 x , r 5 cos δ 1 = r 5 x , r 1 sin θ 1 = r 1 y , r 2 sin ϕ 1 = r 2 y , r 3 sin α 1 = r 3 y , r 4 sin ( ϕ 1 + β ) = r 4 y , r 5 sin δ 1 = r 5 y .
Substituting Equation (33) into Equations (31) and (32), and gathering the coefficients of cos∆ϕj and sin∆ϕj gives
L 1 j + L 2 j cos Δ ϕ j + L 3 j sin Δ ϕ j = 0 ,   j = 2 , , n ,
where
L 1 j = 2 cos Δ θ j ( r 1 x 2 + r 1 y 2 + r 1 x r 2 x r 1 x r 3 x + r 1 y r 2 y r 1 y r 3 y ) + 2 sin Δ θ j ( r 1 x r 2 y + r 1 y r 2 x + r 1 x r 3 y r 1 y r 3 x ) + 2 ( r 2 x 2 + r 2 y 2 r 2 x r 3 x r 2 y r 3 y ) + 2 ( r 1 x 2 + r 1 y 2 + r 1 x r 2 x r 1 x r 3 x + r 1 y r 2 y r 1 y r 3 y ) ,
L 2 j = 2 cos Δ θ j ( r 1 x r 2 x + r 1 y r 2 y ) + 2 sin Δ θ j ( r 1 x r 2 y r 1 y r 2 x ) 2 ( r 1 x r 2 x + r 1 y r 2 y ) 2 ( r 2 x 2 + r 2 y 2 r 2 x r 3 x r 2 y r 3 y ) ,
L 3 j = 2 cos Δ θ j ( r 1 x r 2 y + r 1 y r 2 x ) + 2 sin Δ θ j ( r 1 x r 2 x + r 1 y r 2 y ) + 2 ( r 1 x r 2 y r 1 y r 2 x ) + 2 ( r 2 x r 3 y r 2 y r 3 x ) ,
and
Q 1 j + Q 2 j cos Δ ϕ j + Q 3 j sin Δ ϕ j = 0 ,   j = 2 , , n ,
where
Q 1 j = 2 cos Δ θ j ( r 1 x 2 + r 1 y 2 + r 1 x r 4 x r 1 x r 5 x + r 1 y r 4 y r 1 y r 5 y + r 1 y Δ p j ) + 2 sin Δ θ j ( r 1 x r 4 y + r 1 y r 4 x + r 1 x r 5 y r 1 y r 5 x r 1 x Δ p j ) + 2 ( r 1 x 2 + r 1 y 2 + r 1 x r 4 x r 1 x r 5 x + r 1 y r 4 y r 1 y r 5 y ) + 2 ( r 4 x 2 + r 4 y 2 r 4 x r 5 x r 4 y r 5 y ) + 2 Δ p j ( r 1 y + r 4 y r 5 y ) + Δ p j 2 ,
Q 2 j = 2 cos Δ θ j ( r 1 x r 4 x + r 1 y r 4 y ) + 2 sin Δ θ j ( r 1 x r 4 y r 1 y r 4 x ) 2 ( r 1 x r 4 x + r 1 y r 4 y ) 2 ( r 4 x 2 + r 4 y 2 r 4 x r 5 x r 4 y r 5 y ) 2 r 4 y Δ p j ,
Q 3 j = 2 cos Δ θ j ( r 1 x r 4 y r 1 y r 4 x ) + 2 sin Δ θ j ( r 1 x r 4 x + r 1 y r 4 y ) + 2 ( r 1 x r 4 y r 1 y r 4 x ) + 2 ( r 4 x r 5 y r 4 y r 5 x ) 2 r 4 x Δ p j ,
respectively.
For the elimination of ∆ϕj from Equations (34) and (35), solve the two equations for cos∆ϕj and sin∆ϕj, and simplify using the trigonometric identity to obtain
( L 1 j Q 3 j + L 3 j Q 1 j ) 2 + ( L 1 j Q 2 j L 2 j Q 1 j ) 2 ( L 2 j Q 3 j L 3 j Q 2 j ) 2 = 0 ,   j = 2 , , n .
As in the case of the Watt II six-link slider-crank mechanism, the synthesis equations of the Stephenson III six-link slider-crank function generator for n prescribed positions derived in Equation (36) are formulated by n − 1 eighth-degree polynomials in 10 unknowns which are given in Equation (33). For 9 position synthesis, n is equal to 9 and there are 8 equations in 10 unknowns. Hence, 2 unknowns should be assumed as free choices. In this research, the x- and y-components of link r1 in the initial position are selected as free choices. The total degree of the synthesis equations for 9 precision points is 88 = 16,777,216.

3. Additional Degree Reduction with the Multi-Homogeneous Theorem

The Bézout number of a polynomial system represents the largest number of solutions that the system can have, and is also the number of paths to be tracked to find all isolated solutions in the continuation method. According to the multi-homogeneous Bézout theorem, the number of solution paths and thus the computing time can be significantly reduced by dividing the unknowns into appropriate homogeneous groups by the multi-homogeneous formulation [35,36].
In order to apply the theorem, the 9 position synthesis equations, Equations (19) and (36), need to be converted into homogeneous structures first. For this purpose, quadratic auxiliary equations that replace particular terms in Equations (16a)–(16c), (17a)–(17c), (34a)–(34c), and (35a)–(35c) are defined as new unknowns as follows.
For the Watt II six-bar slider-crank mechanism:
M 1 = r 3 x 2 + r 3 y 2 r 2 x r 3 x r 2 y r 3 y , M 2 = r 2 x r 3 y r 2 y r 3 x , M 3 = r 4 x 2 + r 4 y 2 r 4 x r 5 x r 4 y r 5 y , M 4 = r 4 x r 5 y r 4 y r 5 x .
For the Stephenson III six-bar slider-crank mechanism:
N 1 = r 2 x 2 + r 2 y 2 r 2 x r 3 x r 2 y r 3 y , N 2 = r 2 x r 3 y r 2 y r 3 x , N 3 = r 4 x 2 + r 4 y 2 r 4 x r 5 x r 4 y r 5 y , N 4 = r 4 x r 5 y r 4 y r 5 x .
By substituting Mk in Equation (37) into Lij and Qij of Equation (19), and substituting Nk in Equation (38) into Lij and Qij of Equation (36), where i = 1, …, 3, j = 2, …, 9, and k = 1, …, 4, Equations (19) and (36) become 8 homogenous quartic polynomials, respectively. Hence, now the 9 position synthesis equations are formulated into a system of 12 polynomial equations −8 quartic polynomials, Equations (19) or (36), and 4 quadratic polynomials, Equations (37) or (38)—in 12 unknowns, rix and riy, where i = 2, …, 5, and Mk, or Nk, with k = 1, …, 4. The total degree of the synthesis equations for each mechanism is 4824 = 1,048,576, which is 1/16 of 88 = 16,777,216, the total degree of the simplified synthesis equations described in Section 2.
As the last step for reducing the system size, the multi-homogeneous Bézout theorem is applied to the synthesis equations. Then, the synthesis equations can be solved by the tracking of only 286,720 paths as explained in Section 3.1 and Section 3.2.

3.1. 2-Homogeneous Formulation for Watt II Six-Bar Slider-Crank Mechanism

For the multi-homogeneous formulation of the synthesis equations for the Watt II slider-crank mechanism, the 12 unknowns are arranged into two groups λ1 and λ2 that constitute the loop OABC and loop OCDE shown in Figure 1 as follows.
λ 1 : ( r 2 x , r 2 y , r 3 x , r 3 y , M 1 , M 2 ) ,   λ 2 : ( r 4 x , r 4 y , r 5 x , r 5 y , M 3 , M 4 ) .
In order to determine the multi-homogenous Bézout number for a polynomial system whose unknowns are arranged into m homogeneous groups λ1 to λm, let l be the number of polynomial equations to solve, di,j be the degree of the ith equation with respect to the unknowns in group λj, and kj be the number of unknowns in group λj. Then, the multi-homogenous Bézout number of the system is defined as the coefficient of the term j = 1 m λ j k j of the following equation [35].
i = 1 l ( j = 1 m d i , j λ j ) = ( d 1 , 1 λ 1 + + d 1 , m λ m ) ( d 2 , 1 λ 1 + + d 2 , m λ m ) ( d l , 1 λ 1 + + d l , m λ m ) ,
Note that λj in Equation (40) is used as a null value for calculating the multi-homogenous Bézout number.
Since the system of 12 synthesis equations for the Watt II slider-crank function generator, Equations (19) and (37) with Equation (39), is a 2-homogeneous system, substituting l = 12, m = 2, k1 = k2 = 6, and di,j shown in Table 1 into Equation (40) gives
i = 1 12 ( j = 1 2 d i , j λ j ) = ( 2 λ 1 + 2 λ 2 ) 8 ( 2 λ 1 ) 2 ( 2 λ 2 ) 2 .
Hence, the multi-homogenous Bézout number of the system, the coefficient of the term λ 1 6 λ 2 6 in Equation (41), can be obtained by using the binomial theorem as (288!/4!4!)2222 = 286,720, which agrees with the results reported by Dhingra et al. [20] for the synthesis of the Watt II six-link function generator with revolute joints only for 9 precision positions.
If other unknowns than r1x and r1y are selected as the free choices, the number of auxiliary equations to convert the system into a homogenous structure increases and 2-homogenous Bézout number becomes higher. For example, if r2x and r2y are assumed as the free choices, it is necessary to define 6 auxiliary equations and the 2-homogeneous Bézout number is (288!/4!4!)2422 = 1,146,880. Therefore, it is important to select the unknown that minimize the number of auxiliary equations as the free choices. In this case, it is efficient to assume r1x and r1y as the free choices to reduce the 2-homogeneous Bézout number in consideration of the calculation time.

3.2. 2-Homogeneous Formulation for Stephenson III Six-Bar Slider-Crank

In order to make the Stephenson III six-bar slider-crank synthesis equations a 2-homogeneous system, the 12 unknowns are divided into two groups γ1 and γ2 for the loop OABC and loop OADE shown in Figure 2 as follows.
γ 1 : ( r 2 x , r 2 y , r 3 x , r 3 y , N 1 , N 2 ) ,   γ 2 : ( r 4 x , r 4 y , r 5 x , r 5 y , N 3 , N 4 ) .
As shown in Table 2, the degrees of the individual equations with respect to the unknowns in each group are the same as those of the Watt II six-bar slider-crank function generator. Therefore, the 2-homogeneous Bézout number of the system is also 286,720, which agrees with the results of Dhingra et al. [20] for the synthesis of the Stephenson III six-bar function generator with revolute joints only for 9 precision points. For the same reason as described in Section 3.1, r1x and r1y are assumed to be the free choices in this study to reduce the 2-homogeneous Bézout number.
Consequently, the 2-homogeneous Bézout number of the synthesis equations for each mechanisms is about 1/58 of the total degree of the polynomial systems that are derived by eliminating some unknowns in Section 2.2 and Section 2.4. Hence, the process proposed in this section can significantly reduce the time to solve the synthesis equations.
In this research, the calculation of the synthesis equations is carried out by using the Bertini software program [37], which is based on the homotopy continuation method, for obtaining all the solutions to the Watt II and Stephenson III six-bar slider-crank synthesis equations that satisfy 9 prescribed precision points. Among the real solutions, however, some may have kinematic defects. The following section describes the process of sorting out feasible solutions.

4. Sorting Out Feasible Solutions

In order for the synthesized mechanisms to function as six-link slider-crank linkages, they should satisfy the following two conditions. The first is that the input link needs to be capable of complete rotation, and the second is the prescribed positions should lie on a single stroke of the slider. Hence, the solutions that do not satisfy these two conditions should be excluded. This section describes the process for screening the feasible mechanisms by means of the full rotation condition of the input link and the displacement analysis of synthesized mechanisms to check whether prescribed positions are present on a single stroke of the slider.

4.1. Full Rotatability of Input Link

The input link r1 of each mechanism is in the four-bar loop OABC shown in Figure 1 and Figure 2. Hence, it is required for the four-bar to have a full rotatable crank as the input link. Among the four-bar linkages, those that allow complete rotation of the input link are the crank-rocker and double-crank, and each mechanism satisfies the following condition, respectively [38].
Crank-rocker   m e c h a n i s m : H 1 > 0 , H 2 > 0 , H 3 > 0 , Double-crank   m e c h a n i s m : H 1 < 0 , H 2 < 0 , H 3 > 0 ,
where
H 1 = r 0 r 1 + r 2 r 3 , H 2 = r 0 r 1 r 2 + r 3 , H 3 = r 2 + r 3 r 0 r 1 .
The value of r0 in Equation (44) can be determined by the loop closure equation for the four-bar in each mechanism. From Figure 1 and Equation (2) for the Watt II slider-crank and Figure 2 and Equation (20) for the Stephenson III slider-crank, the x- and y-components of r0 can be written as r0cosη = r1x + r2xr3x and r0sinη = r1y + r2yr3y, respectively, where r1x and r1y are the prescribed values as free choices, and r2x, r2y, r3x, and r3y are obtained from the solution of the synthesis equations. Then, r0 = (r02cos2η + r02sin2η)1/2. Similarly, since rix and riy (i = 1, 2, 3) represent the x- and y-components of ri as defined in Equations (15) and (33), ri in Equation (44) can be determined by ri = (rix2 + riy2)1/2.
Based on the condition provided in Equation (43), the solutions that allow the full rotation of the input link can be found from among all the solutions of the synthesis equations.
The next section describes the displacement analysis for further discerning the mechanisms in which all the prescribed precision points are present on a single stroke of the output slider.

4.2. Displacement Analysis

If all link lengths ri, i = 0, …, 5, and the offset h are given, for a given value of the input crank angle, the Watt II and Stephenson III six-bar slider-crank mechanisms can each be assembled into four different configurations called assembly modes or branches, and some of the prescribed positions of a synthesized mechanism may lie on different assembly modes. If all the prescribed positions do not lie on a single stroke of the slider on one branch, the solution is unacceptable. Therefore, it is necessary to select the solution mechanisms that satisfy this condition by analyzing the displacement of each mechanism.
For the analysis of the six-bar slider-crank mechanism, the four-bar linkages commonly used in the two mechanisms need to be analyzed first. The angular displacements α of link r3 and ϕ of link r2 shown in Figure 1 and Figure 2 can be determined as [39]
α = 2 tan 1 ( B ± B 2 C 2 + A 2 C A ) ,   π α π ,
where
A = 2 r 0 r 3 cos η 2 r 1 r 3 cos θ , B = 2 r 0 r 3 sin η 2 r 1 r 3 sin θ , C = r 0 2 + r 1 3 + r 3 2 r 2 2 2 r 0 r 1 ( cos η cos θ + sin η sin θ ) ,
and
ϕ = tan 1 ( r 0 sin η + r 3 sin α r 1 sin θ r 0 cos η + r 3 cos α r 1 cos θ ) ,
respectively.
Notice that when only the link lengths of the four-bar mechanism and the angle of the ground link η shown in Figure 1 and Figure 2 are given, there are two sets of α and ϕ for a given crank angle θ due to the positive (α+) or negative (α) sign in Equation (45), which correspond to the two assembly modes. However, since the solution of the synthesis equations in this study yields the x- and y-components of each link in its initial position, the sign to be used in Equation (45) and the assembly mode of the four-bar can be readily identified.

4.2.1. Watt II Six-Bar Slider-Crank Mechanism

In order to analyze the slider displacement of the Watt II slider-crank mechanism, the position of the joint E in Figure 1 needs to be determined. To find the coordinates of E, the position of joint D can be expressed as
D = [ D x D y ] = [ r 0 cos η r 0 sin η ] + [ r 4 cos ( α + β ) r 4 sin ( α + β ) ] ,
where α is calculated from Equation (45), and β is the angle between link r3 and link r4 that can be determined using the solution of the synthesis equations as
β = tan 1 ( r 4 y r 4 x ) tan 1 ( r 3 y r 3 x ) .
Now, the position of the moving pivot E is
E = [ E x E y ] = [ h D y ± r 5 2 ( h D x ) 2 ] ,
where
h = r 0 cos η + r 4 x r 5 x .
Equation (50) indicates that there are two solutions for the position of the joint E. That is, the slider and link r5 can be assembled above (with the positive sign: Ey+) or below (with the negative sign: Ey) the moving pivot D. As mentioned in Section 4.1, however, since the result of this research gives the orientation of each link in its initial position, the sign to be used in Equation (50) and the assembly mode of a synthesized mechanism can be identified.

4.2.2. Stephenson III Six-Bar Slider-Crank Mechanism

The position of the moving pivot D in Figure 2 can be determined by
D = [ D x D y ] = [ r 1 cos θ r 1 sin θ ] + [ r 4 cos ( ϕ + β ) r 4 sin ( ϕ + β ) ] ,
where β is the angle between link r2 and link r4 that can be determined using the solution of the synthesis equations as
β = tan 1 ( r 4 y r 4 x ) tan 1 ( r 2 y r 2 x ) .
Then, the coordinates of the moving pivot E are
E = [ E x E y ] = [ h D y ± r 5 2 ( h D x ) 2 ] ,
where
h = r 1 x + r 4 x r 5 x .
As in the case with the Watt II slider-crank mechanism, the ± sign to be used in Equation (54) can be determined by considering the x- and y-components of link r5 of the synthesized function generator.
By applying the solution screening process explained in this section, the Watt II and Stephenson III six-bar slider-cranks without kinematic defects can be finally selected from among all the solutions of the synthesis equations. The next section provides examples of the dimensional synthesis method proposed in this research.

5. Numerical Examples

In this section, examples are given for the dimensional synthesis of the Watt II and Stephenson III six-bar slider-crank function generators that satisfy 9 precision points. The synthesis equations of the two mechanisms derived in this study were calculated by using the homotopy continuation method, through which all the solutions were obtained. Then, the feasible solutions were selected based on the two criteria described in Section 4.

5.1. Dimensional Synthesis of Watt II Six-Bar Slider-Crank Function Generator

The prescribed positions that the Watt II six-bar slider-crank function generator should satisfy are given in Table 3, which shows the relative displacement of the output slider, Δpj = pjp1, with respect to the relative input crank angle, Δθj = θjθ1, for i = 2, …, 9, each measured from the initial positions of the slider and input crank, respectively. The free choices were assumed as r1x = 0.12268 and r1y = 0.87294.
For the given precision points, Equations (19) and (37) are solved by using the Bertini software package. The computation time was 3.42 h on a single node of Intel® Xeon® W-2245 [email protected] GHz. Among 25,630 nonsingular solutions obtained, those that have less than 0.01% structural error are selected as feasible solutions. As a result, there were 37 mechanisms that have completely rotatable cranks and pass all the prescribed positions. Table 4 shows the number of the feasible mechanisms for each assembly configuration. The kinematic diagrams in their initial positions and slider displacements of two synthesized mechanisms are shown in Figure 3 and the corresponding solutions are listed in Table 5.

5.2. Dimensional Synthesis of Stephenson III Six-Bar Slider-Crank Function Generator

The prescribed positions for the Stephenson III six-bar slider-crank function generator are given in Table 6 and the free choices were assumed as r1x = 0.12859 and r1y = 1.0473.
Using the Bertini software package, the computation time to solve the synthesis equations given in Equations (36) and (38) was 4.38 h on the same work station referred to in Section 5.1, and a total of 36,061 nonsingular solutions were obtained. Among them, 31 Stephenson III six-bar slider-cranks that have less than 0.01% structural error were found, and the number of the feasible mechanisms for each assembly configuration is tabulated in Table 7. The kinematic diagrams in the initial positions and slider displacements of two synthesized mechanisms are shown in Figure 4 and the corresponding solutions are listed in Table 8.

6. Conclusions

This paper presents the dimensional synthesis of the Watt II and Stephenson III six-bar slider-crank function generators that satisfy nine precision points of the output slider for the prescribed rotational angles of the input link. In this study, the system of synthesis equations for each mechanism initially derived by 56 quadratic polynomials is simplified to eight eighth-degree polynomial equations by eliminating some unknowns. Then, by defining four auxiliary quadratic polynomials as new unknowns, the synthesis equations are converted into eight quartic and four quadratic homogeneous polynomials in 12 unknowns. By this process, the Bézout number of the system that represents the maximum number of isolated solutions is reduced from 256 ≅ 7.21 × 1016 to 1,048,576. In order to decrease the computation time of the synthesis equations further, 2-homogeneous formulation is applied by arranging the unknowns of the system into two groups and the multi-homogeneous Bézout number of the system is 286,720. In this research, the synthesis equations are solved by using the Bertini solver, which uses the homotopy continuation method.
As is known, all the solutions obtained by solving the synthesis equations do not yield feasible six-link slider-crank mechanisms. For the solution mechanism to function properly as desired without any defects, it must have a full rotatable crank and the prescribed precision points should lie on a single stroke of the slider. The screening process to select suitable linkages among the solutions is explained in detail.
The proposed method is verified by carrying out synthesis examples. For arbitrarily given nine precision points, the methods determined 37 feasible Watt II slider-crank mechanisms among 25,630 nonsingular solutions and 31 Stephenson III six-bdar slider-crank linkages out of 36,061 solutions. The structural errors of the synthesized mechanisms determined in this research are ranged from 5.42 × 10−3% and 7.95 × 10−12%.
In order to synthesize six-link slider-crank function generators for other numbers of precision points than nine, the n-1 polynomials in 10 unknowns given in Equation (19) for the Watt II and Equation (36) for the Stephenson III slider-crank mechanism would be the starting point for the synthesis. Once the number of precision points to be synthesized for is chosen, the number of free choices is determined. When selecting free choices among the 10 unknowns in Equations (15) or (33), it is efficient to choose the unknowns that can minimize the number of auxiliary equations to obtain the lower value of the multi-homogeneous Bézout number of the system. Then, define auxiliary equations that can convert the system into a homogeneous structure and apply multi-homogeneous formulation by arranging the variables into homogeneous groups. By solving the resulting system of equations and sorting out feasible solutions, the dimensions of the six-link slider-crank function generators can be determined.
In order to reduce the computation time to solve the system of synthesis equations, the use of newly developed homotopy continuation methods [40,41,42,43] can be considered. Among them, the Diagonal homotopy method [40] seeks only the nonsingular solutions for a system of nonlinear equations. Since there exist many degenerate solutions in a system of non-linear equations, using this method to solve the system may reduce the calculation time significantly.

Author Contributions

Conceptualization, J.K.S. and M.S.J.; Methodology, M.S.J., H.S.P. and J.K.S.; Software, M.S.J., H.S.P. and W.R.K.; Validation, M.S.J., J.K.S., H.S.P. and W.R.K.; Investigation, M.S.J. and W.R.K.; Data curation, M.S.J. and W.R.K.; Writing—Original Draft Preparation, M.S.J.; Writing—Review and Editing, J.K.S.; Visualization, M.S.J.; Supervision, J.K.S.; project administration, M.S.J. and W.R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Freudenstein, F. An analytical approach to the design of four-link mechanisms. Trans. ASME 1954, 76, 483–492. [Google Scholar] [CrossRef]
  2. Plecnik, M.M.; McCarthy, J.M. Five position synthesis of a Slider-Crank Function generator. In Proceedings of the 35th Mechanisms and Robotics Conference, Parts A and B, Washington, DC, USA, 28–31 August 2011; pp. 317–324. [Google Scholar] [CrossRef]
  3. McLarnan, C.W. Synthesis of Six-Link Plane Mechanisms by Numerical Analysis. J. Eng. Ind. 1963, 85, 5–10. [Google Scholar] [CrossRef]
  4. Suh, C.; Mecklenburg, A. Optimal design of mechanisms with the use of matrices and least squares. Mech. Mach. Theory 1973, 8, 479–495. [Google Scholar] [CrossRef]
  5. Akhras, R.; Angeles, J. Unconstrained nonlinear least-square optimization of planar linkages for rigid-body guidance. Mech. Mach. Theory 1990, 25, 97–118. [Google Scholar] [CrossRef]
  6. Jaiswal, A.; Jawale, H.P. Comparative study of four-bar hyperbolic function generation mechanism with four and five accuracy points. Arch. Appl. Mech. 2017, 87, 2037–2054. [Google Scholar] [CrossRef]
  7. Cabrera, J.; Simon, A.; Prado, M. Optimal synthesis of mechanisms with genetic algorithms. Mech. Mach. Theory 2002, 37, 1165–1177. [Google Scholar] [CrossRef] [Green Version]
  8. Laribi, M.A.; Mlika, A.; Romdhane, L.; Zeghloul, S. A combined genetic algorithm–fuzzy logic method (GA–FL) in mechanisms synthesis. Mech. Mach. Theory 2004, 39, 717–735. [Google Scholar] [CrossRef]
  9. Rai, A.K.; Saxena, A.; Mankame, N.D. Synthesis of Path Generating Compliant Mechanisms Using Initially Curved Frame Elements. J. Mech. Des. 2007, 129, 1056–1063. [Google Scholar] [CrossRef]
  10. Shen, J.; Wang, G.; Bi, Q.; Qu, J. A comprehensive genetic algorithm for design optimization of Z-bar loader working mechanism. J. Mech. Sci. Technol. 2013, 27, 3381–3394. [Google Scholar] [CrossRef]
  11. Acharyya, S.; Mandal, M. Performance of EAs for four-bar linkage synthesis. Mech. Mach. Theory 2009, 44, 1784–1794. [Google Scholar] [CrossRef]
  12. Shiakolas, P.; Koladiya, D.; Kebrle, J. On the optimum synthesis of six-bar linkages using differential evolution and the geometric centroid of precision positions technique. Mech. Mach. Theory 2005, 40, 319–335. [Google Scholar] [CrossRef]
  13. Gatti, G.; Mundo, D. Optimal synthesis of six-bar cammed-linkages for exact rigid-body guidance. Mech. Mach. Theory 2007, 42, 1069–1081. [Google Scholar] [CrossRef]
  14. Gogate, G.R.; Matekar, S.B. Optimum synthesis of motion generating four-bar mechanisms using alternate error functions. Mech. Mach. Theory 2012, 54, 41–61. [Google Scholar] [CrossRef]
  15. Gogate, G.R.; Matekar, S. Unified synthesis of Watt-I six-link mechanisms using evolutionary optimization. J. Mech. Sci. Technol. 2014, 28, 3075–3086. [Google Scholar] [CrossRef]
  16. Huang, X.; He, G.; Liao, Q.; Wei, S.; Tan, X. Solving a Planar Four-Bar Linkages Design Problem. In Proceedings of the International Conference on Information and Automation, Zhuhai/Macau, China, 22–24 June 2009; pp. 1586–1590. [Google Scholar] [CrossRef]
  17. Almandeel, A.; Murray, A.P.; Myszka, D.H.; Stumph, H.E. A Function Generation Synthesis Methodology for All Defect-Free Slider-Crank Solutions for Four Precision Points. J. Mech. Robot. 2015, 7, 031020. [Google Scholar] [CrossRef]
  18. Subbian, T.; Flugrad, D.R. Four-Bar Path Generation Synthesis by a Continuation Method. J. Mech. Des. 1991, 113, 63–69. [Google Scholar] [CrossRef]
  19. Wampler, C.; Morgan, A.P.; Sommese, A.J. Complete Solution of the Nine-Point Path Synthesis Problem for Four-Bar Linkages. J. Mech. Des. 1992, 114, 153–159. [Google Scholar] [CrossRef]
  20. Dhingra, A.K.; Cheng, J.C.; Kohli, D. Synthesis of Six-link, Slider-crank and Four-link Mechanisms for Function, Path and Motion Generation Using Homotopy with m-homogenization. J. Mech. Des. 1994, 116, 1122–1131. [Google Scholar] [CrossRef]
  21. Plecnik, M.M.; McCarthy, J.M. Numerical Synthesis of Six-Bar Linkages for Mechanical Computation. J. Mech. Robot. 2014, 6, 031012. [Google Scholar] [CrossRef] [Green Version]
  22. Plecnik, M.M.; McCarthy, J.M. Computational Design of Stephenson II Six-Bar Function Generators for 11 Accuracy Points. J. Mech. Robot. 2015, 8, 011017. [Google Scholar] [CrossRef]
  23. Plecnik, M.M.; McCarthy, J.M. Kinematic synthesis of Stephenson III six-bar function generators. Mech. Mach. Theory 2016, 97, 112–126. [Google Scholar] [CrossRef] [Green Version]
  24. Freudenstein, F.; Roth, B. Numerical Solution of Systems of Nonlinear Equations. J. ACM 1963, 10, 550–556. [Google Scholar] [CrossRef]
  25. Hu, J.; Sun, Y.; Cheng, Y. High mechanical advantage design of six-bar Stephenson mechanism for servo mechanical presses. Adv. Mech. Eng. 2016, 8, 1–12. [Google Scholar] [CrossRef] [Green Version]
  26. Hu, J.G.; Sun, Y.S.; Cheng, Y.Q.; Ruan, W.P. Numerical Analysis on Transmission Characteristics of Stephenson’s Six-Bar Punching Mechanism with Servo Input. Appl. Mech. Mater. 2011, 86, 623–628. [Google Scholar] [CrossRef]
  27. Podhorodeski, R.P.; Wittchen, J.D.; Nokleby, S.B. Quick-Return Mechanism Design and Analysis Projects. Int. J. Mech. Eng. Educ. 2004, 32, 100–114. [Google Scholar] [CrossRef]
  28. Li, C.-H.; Tso, P.-L. Experimental study on a hybrid-driven servo press using iterative learning control. Int. J. Mach. Tools Manuf. 2008, 48, 209–219. [Google Scholar] [CrossRef]
  29. Du, R.; Guo, W.Z. The Design of a New Metal Forming Press with Controllable Mechanism. J. Mech. Des. 2003, 125, 582–592. [Google Scholar] [CrossRef]
  30. Wilhelm, S.R.; Van De Ven, J.D. Design and Testing of an Adjustable Linkage for a Variable Displacement Pump. J. Mech. Robot. 2013, 5, 041008. [Google Scholar] [CrossRef]
  31. Hiyoshi, R.; Aoyama, S.; Takemura, S.; Ushijima, K.; Sugiyama, T. A Study of a Multiple-Link Variable Compression Ratio System for Improving Engine Performance. SAE Trans. 2006, 115, 309–317. Available online: https://www.jstor.org/stable/44687307 (accessed on 1 September 2022).
  32. Kwak, S.W.; Shim, J.K.; Mo, Y.K. Kinematic Conceptual Design of In-Line Four-Cylinder Variable Compression Ratio Engine Mechanisms Considering Vertical Second Harmonic Acceleration. Appl. Sci. 2020, 10, 3765. [Google Scholar] [CrossRef]
  33. Mo, Y.K.; Shim, J.K.; Kwak, S.W.; Jo, M.S.; Park, H.S. Type Synthesis of Variable Compression Ratio Engine Mechanisms. Appl. Sci. 2020, 10, 6574. [Google Scholar] [CrossRef]
  34. Hoeltgebaum, T.; Simoni, R.; Martins, D. Reconfigurability of engines: A kinematic approach to variable compression ratio engines. Mech. Mach. Theory 2016, 96, 308–322. [Google Scholar] [CrossRef]
  35. Tsai, L.W. Robot Analysis: The Mechanics of Serial and Parallel Manipulators, 1st ed.; Wiley-Interscience: Hoboken, NJ, USA, 1999. [Google Scholar]
  36. Wampler, C. Bezout number calculations for multi-homogeneous polynomial systems. Appl. Math. Comput. 1992, 51, 143–157. [Google Scholar] [CrossRef]
  37. Bates, D.J.; Hauenstein, J.D.; Sommese, A.J.; Wampler, C.W. Numerically Solving Polynomial Systems with Bertini; SIAM Press: Philadelphia, PA, USA, 2013; Available online: https://bertini.nd.edu/ (accessed on 1 May 2022).
  38. McCarthy, J.M. Geometric Design of Linkages; Springer: New York, NY, USA, 2000. [Google Scholar]
  39. Waldron, K.J.; Kinzel, G.L.; Agrawal, S.K. Kinematics, Dynamics, and Design of Machinery; John Wiley & Sons: New York, NY, USA, 1999. [Google Scholar]
  40. Sommese, A.J.; Verschelde, J.; Wampler, C.W. Solving Polynomial Systems Equation by Equation. In Algorithms in Algebraic Geometry; The IMA Volumes in Mathematics and its Applications; Dickenstein, A., Schreyer, F.O., Sommese, A.J., Eds.; Springer: New York, NY, USA, 2008; Volume 146. [Google Scholar] [CrossRef] [Green Version]
  41. Li, T.Y. Numerical solution of polynomial systems by homotopy continuation methods. In Handbook of Numerical Analysis; Special Volume: Foundations of Computational Mathematics; Cucker, F., Ed.; North-Holland: Harlem, The Netherlands, 2003; Volume 11, pp. 209–304. [Google Scholar] [CrossRef]
  42. Wu, T.-M. A study of convergence on the Newton-homotopy continuation method. Appl. Math. Comput. 2005, 168, 1169–1174. [Google Scholar] [CrossRef]
  43. Malinen, I.; Tanskanen, J. Homotopy parameter bounding in increasing the robustness of homotopy continuation methods in multiplicity studies. Comput. Chem. Eng. 2010, 34, 1761–1774. [Google Scholar] [CrossRef]
Figure 1. Watt II six-bar slider-crank mechanism: (a) Design parameters; (b) The link vectors drawn in its first and jth positions.
Figure 1. Watt II six-bar slider-crank mechanism: (a) Design parameters; (b) The link vectors drawn in its first and jth positions.
Applsci 12 10503 g001
Figure 2. Stephenson III six-bar slider-crank mechanism: (a) Design parameters; (b) The link vectors drawn in its first and jth positions.
Figure 2. Stephenson III six-bar slider-crank mechanism: (a) Design parameters; (b) The link vectors drawn in its first and jth positions.
Applsci 12 10503 g002
Figure 3. Watt II six-bar slider-crank function generator: (a) (α+, Ey+) with a crank and the displacement of the slider; (b) (α+, Ey+) with double cranks and the displacement of the slider.
Figure 3. Watt II six-bar slider-crank function generator: (a) (α+, Ey+) with a crank and the displacement of the slider; (b) (α+, Ey+) with double cranks and the displacement of the slider.
Applsci 12 10503 g003
Figure 4. Stephenson III six-bar slider-crank function generator: (a) (α, Ey+) with a crank and the displacement of the slider; (b) (α+, Ey+) with double cranks and the displacement of the slider.
Figure 4. Stephenson III six-bar slider-crank function generator: (a) (α, Ey+) with a crank and the displacement of the slider; (b) (α+, Ey+) with double cranks and the displacement of the slider.
Applsci 12 10503 g004aApplsci 12 10503 g004b
Table 1. Degrees of the 2-homogenous system for Watt II six-bar slider-crank synthesis equations.
Table 1. Degrees of the 2-homogenous system for Watt II six-bar slider-crank synthesis equations.
Equations (i = 1, …, 12)Equation (19) (i = 1, …, 8)M1
(i = 9)
M2
(i = 10)
M3
(i = 11)
M4
(i = 12)
di,1 (Group λ1)22200
di,2 (Group λ2)20022
Table 2. Degrees of the 2-homogenous system for the Stephenson III six-bar slider-crank synthesis equations.
Table 2. Degrees of the 2-homogenous system for the Stephenson III six-bar slider-crank synthesis equations.
Equations (i = 1, …, 12)Equation (36) (i = 1, …, 8)N1
(i = 9)
N2
(i = 10)
N3
(i = 11)
N4
(i = 12)
di,1 (Group γ1)22200
di,2 (Group γ2)20022
Table 3. Prescribed positions for the Watt II slider-crank function generator.
Table 3. Prescribed positions for the Watt II slider-crank function generator.
Precision Points (j)123456789
θj (deg)02170100124164193224298
pj0−0.49087−1.45837−1.69238−1.77397−1.77643−1.67172−1.42028−0.13685
Table 4. The number of feasible mechanisms for each configuration.
Table 4. The number of feasible mechanisms for each configuration.
ConfigurationNumber of Feasible Mechanisms
1 (Assembled by α, Ey+)0
2 (Assembled by α, Ey)1
3 (Assembled by α+, Ey+)25
4 (Assembled by α+, Ey)11
Table 5. The solutions of the Watt II six-bar slider-crank mechanism shown in Figure 3.
Table 5. The solutions of the Watt II six-bar slider-crank mechanism shown in Figure 3.
ConfigurationType of Four-Bar MechanismSolutionsLink LengthsOrientations (Deg)
3Crank-rockerr2x1.83719848448098r22.66481499912287ϕ46.4151406996215
r2y1.93026959468645
r3x2.29173698831327r32.29948451051778α−4.70463367830402
r3y−0.188601671540832
r4x2.95628516701656r43.02185450092804α + β11.9574803117514
r4y0.626085166776018
r5x2.50589035665077r53.22236485855068δ−38.9535693977896
r5y−2.02586988774379
3Double-crankr2x−0.492252951916725r20.660020417350853ϕ−138.229126201616
r2y−0.439674860151526
r3x−0.167431742713909r30.459407273662191α111.373840719186
r3y0.427810302149809
r4x−0.66096098576403r40.97212798670023α + β132.83672677491
r4y0.712855804369779
r5x−0.333010353208065r51.40620185002117δ−103.698662763378
r5y−1.36620194248844
Table 6. Prescribed positions for the Stephenson III slider-crank function generator.
Table 6. Prescribed positions for the Stephenson III slider-crank function generator.
Precision Points (j)123456789
θj (deg)03988140182225253287333
pj0−0.16691−1.08488−2.29326−2.83569−2.59666−1.93088−0.95797−0.18975
Table 7. The number of feasible mechanisms for each configuration.
Table 7. The number of feasible mechanisms for each configuration.
ConfigurationNumber of Feasible Mechanisms
1 (Assembled by α, Ey+)11
2 (Assembled by α, Ey)0
3 (Assembled by α+, Ey+)12
4 (Assembled by α+, Ey)8
Table 8. The solutions of the Stephenson III six-bar slider-crank mechanism shown in Figure 4.
Table 8. The solutions of the Stephenson III six-bar slider-crank mechanism shown in Figure 4.
ConfigurationType of Four-Bar MechanismSolutionsLink LengthsOrientations (Deg)
1Crank-rockerr2x2.54452041055486r22.99784510000088ϕ31.9204361635379
r2y1.58508394852419
r3x−0.439101961903118r33.00115483999674α98.4132035930371
r3y2.96885833961284
r4x0.265943554102272r42.55948745154112ϕ + β84.0359171635551
r4y2.54563352441547
r5x0.272684581643968r53.33598841545769δ−85.3113966753684
r5y−3.32482508216621
3Double-crankr2x−0.448150890683856r20.943859211257232ϕ−118.346779745901
r2y−0.830681039782654
r3x0.003969289812873r30.854831121637962α89.7339540343908
r3y0.854821906164786
r4x−0.725026834630641r40.729835893584505ϕ + β−173.418970576143
r4y−0.083645206854701
r5x0.650197621387553r52.41279005227063δ−74.3666995186773
r5y-2.32353155551589
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jo, M.S.; Shim, J.K.; Park, H.S.; Kim, W.R. Dimensional Synthesis of Watt II and Stephenson III Six-Bar Slider-Crank Function Generators for Nine Prescribed Positions. Appl. Sci. 2022, 12, 10503. https://doi.org/10.3390/app122010503

AMA Style

Jo MS, Shim JK, Park HS, Kim WR. Dimensional Synthesis of Watt II and Stephenson III Six-Bar Slider-Crank Function Generators for Nine Prescribed Positions. Applied Sciences. 2022; 12(20):10503. https://doi.org/10.3390/app122010503

Chicago/Turabian Style

Jo, Min Seok, Jae Kyung Shim, Ho Sung Park, and Woon Ryong Kim. 2022. "Dimensional Synthesis of Watt II and Stephenson III Six-Bar Slider-Crank Function Generators for Nine Prescribed Positions" Applied Sciences 12, no. 20: 10503. https://doi.org/10.3390/app122010503

APA Style

Jo, M. S., Shim, J. K., Park, H. S., & Kim, W. R. (2022). Dimensional Synthesis of Watt II and Stephenson III Six-Bar Slider-Crank Function Generators for Nine Prescribed Positions. Applied Sciences, 12(20), 10503. https://doi.org/10.3390/app122010503

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop