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 r
5 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
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
and
respectively. Subtracting Equation (2) from Equation (3) gives
Loop OCDE:
For the initial and
jth positions, the loop closure equations can be written as
and
respectively. Subtracting Equation (5) from Equation (6) gives
Separating the real and imaginary parts of Equations (4) and (7) leads to
and
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
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
, where
di is the degree of the
ith polynomial [
35]. Hence, the total degree of the synthesis equations for 9 precision points is 2
56 ≅ 7.21 × 10
16. 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 cos
2(
ϕ1 + ∆
ϕj)+sin
2(
ϕ1 + ∆
ϕj) = 1 to obtain
Similarly, the unknown ∆
δj can be eliminated from Equations (10) and (11) as
Let the
x- and
y-components of the initial positions of the links
r1,
r2,
r3,
r4, and
r5 be
Substituting Equation (15) into Equation (13) and gathering the coefficients of cos∆
αj and sin∆
αj gives
where
Similarly, substituting Equation (15) into Equation (14) and combining the coefficients of cos∆
αj and sin∆
αj yields
where
In order to eliminate ∆
αj from Equations (16) and (17), solve the two equations for cos∆
αj and sin∆
αj to obtain
Squaring both sides of Equation (18), adding the result, and simplifying using the trigonometric identity cos
2(∆
αj) + sin
2(∆
αj) = 1 gives
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 8
8 = 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
and
respectively. Subtracting Equation (20) from Equation (21) gives
Loop OADE:
For the initial and
jth positions, the loop closure equations can be written as
and
respectively. Subtracting Equation (23) from Equation (24) gives
Separating the real and imaginary parts of Equations (22) and (25) yields
and
respectively. Among the unknowns in Equations (26) through (29), ∆
ϕj, ∆
αj, and ∆
δj should satisfy trigonometric identities
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
and the unknown ∆
δj can be eliminated from Equations (28) and (29) as
Let the
x- and
y-components of the initial positions of the links
r1,
r2,
r3,
r4, and
r5 be
Substituting Equation (33) into Equations (31) and (32), and gathering the coefficients of cos∆
ϕj and sin∆
ϕj gives
where
and
where
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
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:
For the Stephenson III six-bar slider-crank mechanism:
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 4
82
4 = 1,048,576, which is 1/16 of 8
8 = 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.
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
of the following equation [
35].
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
Hence, the multi-homogenous Bézout number of the system, the coefficient of the term
in Equation (41), can be obtained by using the binomial theorem as (2
88!/4!4!)2
22
2 = 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.
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].
where
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 +
r2x −
r3x and
r0sin
η =
r1y +
r2y −
r3y, 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 = (
r02cos
2η +
r02sin
2η)
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]
where
and
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
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
Now, the position of the moving pivot
E is
where
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
where
β is the angle between link
r2 and link
r4 that can be determined using the solution of the synthesis equations as
Then, the coordinates of the moving pivot
E are
where
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.
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.