Computer-aided synthesis of spherical and planar 4R linkages for four specified orientations

According to the burmester theory, an infinite number of spherical or planar 4R linkages for a specific four-orientation task can be synthesized, but most of the linkage solutions calculated by this method are invalid because of motion defect, poor performance and others. In order to improve the synthesis efficiency, a program package based on Matlab is developed to find a satisfied linkage solution automatically and quickly. Firstly, the calculation on circle points and the center points based on the burmester theory in spherical problems is introduced. Secondly, the calculation methods of linkage defect discrimination, linkage type classification, linkage performance evaluation and solutions visualization based on the theory of spherical trigonometry are presented respectively. Thirdly, the synthesis calculation of program package is extended to the planar 4R linkage based on the theory of planar analytic geometry. Finally, the examples of the spherical synthesis problem and the planar synthesis problem based on solutions map are introduced to test the program package, the result proves this program package is effective and flexible.

1 Introduction

In this paper, the study work is to design a 4R linkage which can guide its coupler to pass four specified orientations smoothly and orderly. In order to accomplish this design task, graphing method and computer-aided synthesis method are used mainly at present. The theoretical basis of the graphing method is the theory of burmester, that is to say, connecting any two groups of center points and their corresponding circle points can obtain the desired linkage solution. Although the selection of center points or circle points are optional, majority of the linkage solutions have some problems such as circuit defect, branch defect, order defect or poor motion performance (Wang et al., 2010; Baskar and Bandyopadhyay, 2019a). For the linkage motion defect, the classical defect discrimination method is to study the segmentation technique of the burmester curve (Filemon, 1972; Waldron and Strong, 1978; Gupta and Beloiu, 1998). Considering the inefficiency of the graphing method and the inaccuracy of the discrimination of defective linkages, such as discarding some feasible solutions as uncertain solutions, its application have some limitations. With the development of computer-aided design, the burmester curve mentioned above can be obtained by programming with analytic geometry, but the existence of defective linkages makes it difficult to obtain feasible solutions. In order to get the most candidate solutions, one of the approaches is to allow some task variables to change in a certain range, so as to expand the number of solutions (Zhao et al., 2016; Baskar and Bandyopadhyay, 2019b; Chanekar et al., 2013; Sun et al., 2016). Although this method can increase the number of solutions, most of the increased linkage solutions still have motion defects, so it does not fundamentally eliminate the blindness of linkage synthesis process. Considering that modern computers can obtain massive solutions quickly and discriminate the defective linkage solutions in turn later (Myszka et al., 2009; Tipparthi and Larochelle, 2011; Bai, 2017), the method of type map (Han and Cao, 2019) is gradually applied in some linkage synthesis software, such as Synthetic (Su and McCarthy, 2005), Lincages (Erdman and Loftness, 2005), Sphinx (the latest version is SphinxPC) (Ruth and McCarthy, 1999; McCarthy, 2013) and the software developed by Han et al. (2013). Among them, specially, Sphinx is the first software to solve the spherical problems. The method of type map maps all the linkage solutions obtained by computer into a two-dimensional coordinate system and marks different linkage types with different colors. At the same time, the defective linkage solutions can be hidden. The blindness of linkage synthesis can be reduced by using the type map based on these synthesis software, but there are still some shortcomings: firstly, the biggest problem is that it is not flexible enough and must be designed the linkage according to the process set by the software; secondly, few software besides Sphinx have the synthesis ability of spherical linkage; finally, the type map method can only express the solutions distribution with different linkage type in the solution domain, however, in practical synthesis problems, the transmission performance and installation position of linkage solutions are more important than linkage types (Mendoza-Trejo et al., 2015). In order to make up for the shortcomings of the graphing method and computer-aided synthesis method mentioned above, a series of flexible functions for the synthesis, analysis and plotting of linkages are constructed under Matlab. These functions can be used for the synthesis calculation of spherical and planar 4R linkages. Compared with the previous studies, the main work of this study is as follows: the existing color-coded type map method is improved to express not only the discrete attribute of linkage type, but also the continuous attribute such as transmission performance; a defect discrimination method based on geometric method is proposed, and a pyramid structure of the circuit, branch and order defect discrimination are defined. This method is not only logical, but also accurate without misjudgment.

2 Principle of the computer-aided linkage synthesis

The working principle of the program package based on Matlab is shown in Fig. 1. The first step of linkage synthesis is to input design tasks, i.e. four guiding orientations on sphere or plane. The second step of linkage synthesis is to generate the center points and circle points. Based on the four given orientations, a spherical or planar 4R linkage consisting of poles and image poles is obtained. According to Burmester theory, the center points and the circle points can be obtained by kinematics analysis of the 4R linkage. The third step of linkage synthesis is to calculate the attributes of all candidate solutions. Number the center points from 1 to N in the order in which it is generated. For any center point i(i≤N), a candidate linkage solution can be formed by combining any point of the other (N−1) center points. On the basis of it, $N\times (N-\mathrm{1})$ groups of linkage solutions can be formed. For each linkage solution formed by the center point i and the circle point $j(j\le N,i\ne j)$, the linkage type, the motion defect type and the transmission performance are calculated in turn, and the calculation results are written into the data structure of the linkage solution. After above calculation, the set of candidate solutions is built. The last step of linkage synthesis is to select and analyze linkage solutions. All the candidate linkage solutions are mapped to the solutions map. Any point on the solutions map represents a linkage solution. Its color represents the linkage type and its height represents the level of transmission performance. In the solutions map, all defect solutions are hidden, which improves the efficiency of linkage synthesis. When the user uses the mouse to select a point on the solutions map, the program automatically calculates and displays the corresponding 4R linkage and its coupler curve. Users can watch the animation demonstration of the 4R linkage, or make further analysis of it.

Figure 1 The flow chart of the linkage synthesis calculation.

Specially, the synthesis on spherical 4R linkage and planar 4R linkage are both supported by this program package. Considered the calculation on spherical 4R linkage is more complex than that on planar 4R linkage, the principle of the synthesis calculation on spherical problem is mainly presented below, and the planar problem is introduced as an extension of the program package.

3 Generation of the center and circle points

Based on burmester theory for spherical problems, the circle points are calculated by kinematic analysis of 4R linkage constructed by image poles. The theory is explained by Shirazi (2005), as shown in Fig. 2. When the driving link of 4R linkage constructed by poles (which is denoted by P_ij) and image poles (which is denoted by P${}_{ij}^{\mathrm{1}}$) rotates the angle ω, the 4R linkage moves from the orientation P₁₂ P₁₃ P${}_{\mathrm{34}}^{\mathrm{1}}$ P${}_{\mathrm{24}}^{\mathrm{1}}$ to P₁₂ P₁₃ P${}_{\mathrm{34}}^{{\mathrm{1}}^{\prime }}$ P${}_{\mathrm{24}}^{{\mathrm{1}}^{\prime }}$. The circle point C₁ is determined by the intersection of sphere and two planes P_M and P_N, the plane P_M is defined by the center point of sphere (which is denoted by point O) and normal vector P${}_{\mathrm{24}}^{\mathrm{1}}$ P${}_{\mathrm{24}}^{{\mathrm{1}}^{\prime }}$, the point P_N is defined by the center point of sphere and normal vector P${}_{\mathrm{34}}^{\mathrm{1}}$ P${}_{\mathrm{34}}^{{\mathrm{1}}^{\prime }}$.

For each circle point C₁, the center point C₀ is obtained as following steps. Firstly, define a point P₁₂₃ as the reflection point of C₁ relative to the plane P${}_{{\mathrm{OP}}_{\mathrm{12}}{\mathrm{P}}_{\mathrm{13}}}$. The plane P${}_{{\mathrm{OP}}_{\mathrm{12}}{\mathrm{P}}_{\mathrm{13}}}$ is defined by the center point of sphere, the pole P₁₂ and P₁₃. And define the point P₁₂₄ as the reflection point of C₁ relative to the plane P${}_{{\mathrm{OP}}_{\mathrm{12}}{\mathrm{P}}_{\mathrm{14}}}$, the plane P${}_{{\mathrm{OP}}_{\mathrm{12}}{\mathrm{P}}_{\mathrm{14}}}$ is defined by the center point of sphere, the pole P₁₂ and P₁₄. Relative to the plane P${}_{{\mathrm{OP}}_{\mathrm{12}}{\mathrm{P}}_{\mathrm{23}}}$, the circle point at the second specified orientation C₂ is the reflection point of P₁₂₃. P${}_{{\mathrm{OP}}_{\mathrm{12}}{\mathrm{P}}_{\mathrm{23}}}$ is defined by the center point of sphere, the pole P₁₂ and P₂₃. Relative to the plane P${}_{{\mathrm{OP}}_{\mathrm{13}}{\mathrm{P}}_{\mathrm{23}}}$, the circle point at the third specified orientation C₃ is the reflection point of P₁₂₃. P${}_{{\mathrm{OP}}_{\mathrm{13}}{\mathrm{P}}_{\mathrm{23}}}$ is defined by the center point of sphere, the pole P₁₃ and P₂₃. Relative to the plane P${}_{{\mathrm{OP}}_{\mathrm{14}}{\mathrm{P}}_{\mathrm{24}}}$, the circle point at the fourth specified orientation C₄ is the reflection point of P₁₂₄. P${}_{{\mathrm{OP}}_{\mathrm{14}}{\mathrm{P}}_{\mathrm{24}}}$ is defined by the center point of sphere, the pole P₁₄ and P₂₄. Then, the center point C₀ is determined by the intersection of the sphere and any two perpendicular bisector planes of C₁C₂, C₁C₃ and C₁C₄. With the rotation of the driving link of 4R linkage P₁₂ P₁₃ P${}_{\mathrm{34}}^{\mathrm{1}}$ P${}_{\mathrm{24}}^{\mathrm{1}}$, a series of circle points or center points are generated. Note that two groups of center points and circle points are obtained for each new orientation of linkage P₁₂ P₁₃ P${}_{\mathrm{34}}^{\mathrm{1}}$ P${}_{\mathrm{24}}^{\mathrm{1}}$, and they are both valid and symmetrical for the center point of sphere.

In order to realize the calculation and visualization of linkage synthesis mentioned above, a series of functions on geometric calculation are defined in this program package, such as the calculation functions on rotation, mirror, intersection.

Figure 2 Principle of the spherical 4R linkage synthesis (reproduced from Shirazi, 2007).

4 Calculation and analysis of the linkage solution

4.1 Circuit defect discrimination

If the driving link of spherical 4R linkage rotates partially, there will be a circuit defect when driving link of given four coupler orientations located at different moving intervals. In order to discriminate this motion defect, the limited angles of driving link should be calculated firstly. Figure 3 shows a configuration of spherical triangle which is changed from spherical 4R linkage ABCD(P), where P is the guiding point on the coupler link BC; P₁–P₄ are the given four guiding orientations; α–γ and θ₁∼θ₄ represent the length and roll angle of each link of spherical linkage respectively; δ represents the angle between the link BC and BP. The point coordinates of sphere are defined by longitude θ and latitude φ, for example, the coordinates of the point A can be expressed as (θ_A, φ_A).

Figure 3 Spherical 4R linkage located at the limited orientations.

For the configuration of spherical triangle A(BC)D, define $a=\mathit{\eta }+\mathit{\beta }$, b=γ, c=α and $p=(a+b+c)/\mathrm{2}$, the limited angles of driving link with θ₄=0 is calculated by the half-angle formula

$$\begin{array}{}\text{(1)}& \left\{\begin{array}{l}\sin {\displaystyle \frac{\mathrm{\angle }A}{\mathrm{2}}}=\sqrt{{\displaystyle \frac{\sin \left(p-b\right)\sin \left(p-c\right)}{\sin b\sin c}}}\\ \cos {\displaystyle \frac{\mathrm{\angle }A}{\mathrm{2}}}=\sqrt{{\displaystyle \frac{\sin p\sin \left(p-a\right)}{\sin b\sin c}}}\end{array}\right.\end{array}$$

where ∠A represents the angle from the link AD to the link AB. In addition, the angle θ₄ can be expressed as

$$\begin{array}{}\text{(2)}& \left\{\begin{array}{l}\sin {\mathit{\theta }}_{\mathrm{4}}={\displaystyle \frac{\cos {\mathit{\varphi }}_{\mathrm{D}}\sin \left({\mathit{\theta }}_{\mathrm{A}}-{\mathit{\theta }}_{\mathrm{D}}\right)}{\sin \mathit{\gamma }}}\\ \cos {\mathit{\theta }}_{\mathrm{4}}={\displaystyle \frac{\sin {\mathit{\varphi }}_{\mathrm{D}}-\sin {\mathit{\varphi }}_{\mathrm{A}}\cos \mathit{\gamma }}{\cos {\mathit{\varphi }}_{\mathrm{A}}\sin \mathit{\gamma }}}\end{array}\right.\end{array}$$

Then, the limited angles of driving link which is denoted by ${\mathit{\theta }}_{\mathrm{1}}^{+}$ and ${\mathit{\theta }}_{\mathrm{1}}^{-}$ can be calculated as

$$\begin{array}{}\text{(3)}& \left\{\begin{array}{ll}{\mathit{\theta }}_{\mathrm{1}}^{+}& =\mathrm{2}\arctan \mathrm{2}\left(\sin {\displaystyle \frac{\mathrm{\angle }A}{\mathrm{2}}},\cos {\displaystyle \frac{\mathrm{\angle }A}{\mathrm{2}}}\right)\\ & +\arctan \mathrm{2}\left(\sin {\mathit{\theta }}_{\mathrm{4}},\cos {\mathit{\theta }}_{\mathrm{4}}\right)\\ {\mathit{\theta }}_{\mathrm{1}}^{-}& =-\mathrm{2}\arctan \mathrm{2}\left(\sin {\displaystyle \frac{\mathrm{\angle }A}{\mathrm{2}}},\cos {\displaystyle \frac{\mathrm{\angle }A}{\mathrm{2}}}\right)\\ & +\arctan \mathrm{2}\left(\sin {\mathit{\theta }}_{\mathrm{4}},\cos {\mathit{\theta }}_{\mathrm{4}}\right)\end{array}\right.\end{array}$$

If there is no real root in Eq. (1), the calculation will be invalid, and the set of the limited angles of the driving link (which is denoted by ${\mathit{\theta }}_{\mathrm{11}}^{\ast }$) will be empty, that is ${\mathit{\theta }}_{\mathrm{11}}^{\ast }=\left[\right]$. If not, ${\mathit{\theta }}_{\mathrm{11}}^{\ast }=[{\mathit{\theta }}_{\mathrm{1}}^{+},{\mathit{\theta }}_{\mathrm{1}}^{-}]$.

Similarly, for the configuration of spherical triangle A(B^′C^′)D, define $a=|\mathit{\eta }-\mathit{\beta }|$, b=γ, c=α and $p=(a+b+c)/\mathrm{2}$. If there is no real root in Eq. (1), the calculation will be invalid, and the set of the limited angles of the driving link (which is denoted by ${\mathit{\theta }}_{\mathrm{12}}^{\ast }$) will be empty, that is ${\mathit{\theta }}_{\mathrm{12}}^{\ast }=\left[\right]$. If not, ${\mathit{\theta }}_{\mathrm{12}}^{\ast }=[{\mathit{\theta }}_{\mathrm{1}}^{+},{\mathit{\theta }}_{\mathrm{1}}^{-}]$.

The final set of the limited angles of driving link denoted by ${\mathit{\theta }}_{\mathrm{1}}^{\ast }$, consider ${\mathit{\theta }}_{\mathrm{1}}^{\ast }={\mathit{\theta }}_{\mathrm{11}}^{\ast }\cup {\mathit{\theta }}_{\mathrm{12}}^{\ast }$. Obviously, the size of ${\mathit{\theta }}_{\mathrm{1}}^{\ast }$ which is denoted by $N_{\mathrm{1}}^{\ast }$ is only 0, 2 or 4. Due to the driving link can't move between two adjacent moving intervals, the circuit defect can only be checked when $N_{\mathrm{1}}^{\ast }=\mathrm{4}$. In this case, if all input angles of the driving link are corresponding to the specified coupler orientations located at the same moving interval, the linkage solution will be valid. If not, there will be a circuit defect in the motion of spherical 4R linkage. Specially, the input angle of driving link is expressed as

$$\begin{array}{}\text{(4)}& \left\{\begin{array}{l}\sin {\mathit{\theta }}_{\mathrm{1}}={\displaystyle \frac{\cos {\mathit{\varphi }}_{\mathrm{B}}\sin \left({\mathit{\theta }}_{\mathrm{A}}-{\mathit{\theta }}_{\mathrm{B}}\right)}{\sin \mathit{\alpha }}}\\ \cos {\mathit{\theta }}_{\mathrm{1}}={\displaystyle \frac{\sin {\mathit{\varphi }}_{\mathrm{B}}-\sin {\mathit{\varphi }}_{\mathrm{A}}\cos \mathit{\alpha }}{\cos {\mathit{\varphi }}_{\mathrm{A}}\sin \mathit{\alpha }}}\end{array}\right.\end{array}$$

4.2 Branch defect discrimination

Figure 4 shows two assembly driven link configurations of spherical 4R linkage ABCD and ABC^′D. Due to the linkage is difficult to move from one branch to other without jamming, a branch defect will appear once the driven link in given four coupler orientations locates at different assembly configurations. Actually, the branch of linkage is determined by the solution of motion equation of driven link.

Figure 4 Constructions of the driven link of spherical 4R linkage.

According to the coordinate of point A(θ_A,φ_A), the length of AB (that is α) and the roll angle of AB (that is θ₁), the coordinate of point B(θ_B,φ_B) is expressed as

$$\begin{array}{}\text{(5)}& \left\{\begin{array}{l}\sin {\mathit{\varphi }}_{\mathrm{B}}=\cos \mathit{\alpha }\sin {\mathit{\varphi }}_{\mathrm{A}}+\sin \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{1}}\cos {\mathit{\varphi }}_{\mathrm{A}}\\ \sin {\mathit{\theta }}_{\mathrm{B}}={\displaystyle \frac{{\displaystyle \begin{array}{c}-\sin \mathit{\alpha }\sin {\mathit{\theta }}_{\mathrm{1}}\cos {\mathit{\theta }}_{\mathrm{A}}+\cos \mathit{\alpha }\sin {\mathit{\theta }}_{\mathrm{A}}\cos {\mathit{\varphi }}_{\mathrm{A}}\\ -\sin \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{1}}\sin {\mathit{\theta }}_{\mathrm{A}}\sin {\mathit{\varphi }}_{\mathrm{A}}\end{array}}}{\cos {\mathit{\varphi }}_{\mathrm{B}}}}\\ \cos {\mathit{\theta }}_{\mathrm{B}}={\displaystyle \frac{{\displaystyle \begin{array}{c}\sin \mathit{\alpha }\sin {\mathit{\theta }}_{\mathrm{1}}\sin {\mathit{\theta }}_{\mathrm{A}}+\cos \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{A}}\cos {\mathit{\varphi }}_{\mathrm{A}}\\ -\sin \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{1}}\cos {\mathit{\theta }}_{\mathrm{A}}\sin {\mathit{\varphi }}_{\mathrm{A}}\end{array}}}{\cos {\mathit{\varphi }}_{\mathrm{B}}}}\end{array}\right.\end{array}$$

Similarly, the coordinate of point C(θ_C,φ_C) can be expressed as

$$\begin{array}{}\text{(6)}& \left\{\begin{array}{l}\sin {\mathit{\varphi }}_{\mathrm{C}}=\cos \mathit{\beta }\sin {\mathit{\varphi }}_{\mathrm{D}}+\sin \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{3}}\cos {\mathit{\varphi }}_{\mathrm{D}}\\ \sin {\mathit{\theta }}_{\mathrm{C}}={\displaystyle \frac{{\displaystyle \begin{array}{c}-\sin \mathit{\beta }\sin {\mathit{\theta }}_{\mathrm{3}}\cos {\mathit{\theta }}_{\mathrm{D}}+\cos \mathit{\beta }\sin {\mathit{\theta }}_{\mathrm{D}}\cos {\mathit{\varphi }}_{\mathrm{D}}\\ -\sin \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{3}}\sin {\mathit{\theta }}_{\mathrm{D}}\sin {\mathit{\varphi }}_{\mathrm{D}}\end{array}}}{\cos {\mathit{\varphi }}_{\mathrm{C}}}}\\ \cos {\mathit{\theta }}_{\mathrm{C}}={\displaystyle \frac{{\displaystyle \begin{array}{c}\sin \mathit{\beta }\sin {\mathit{\theta }}_{\mathrm{3}}\sin {\mathit{\theta }}_{\mathrm{D}}+\cos \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{D}}\cos {\mathit{\varphi }}_{\mathrm{D}}\\ -\sin \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{3}}\cos {\mathit{\theta }}_{\mathrm{D}}\sin {\mathit{\varphi }}_{\mathrm{D}}\end{array}}}{\cos {\mathit{\varphi }}_{\mathrm{C}}}}\end{array}\right.\end{array}$$

The relationship between point B and C is determined by the length of link BC (that is η), that is

$$\begin{array}{}\text{(7)}& \cos \mathit{\eta }=\sin {\mathit{\varphi }}_{\mathrm{B}}\sin {\mathit{\varphi }}_{\mathrm{C}}+\cos {\mathit{\varphi }}_{\mathrm{B}}\cos {\mathit{\varphi }}_{\mathrm{C}}\cos \left({\mathit{\theta }}_{\mathrm{B}}-{\mathit{\theta }}_{\mathrm{C}}\right)\end{array}$$

Substitute the Eqs. (5) and (6) into Eq. (7) and obtain following equation, that is

$$\begin{array}{}\text{(8)}& {\mathit{\tau }}_{\mathrm{1}}\sin {\mathit{\theta }}_{\mathrm{3}}\left({\mathit{\theta }}_{\mathrm{1}}\right)+{\mathit{\tau }}_{\mathrm{2}}\cos {\mathit{\theta }}_{\mathrm{3}}\left({\mathit{\theta }}_{\mathrm{1}}\right)+{\mathit{\tau }}_{\mathrm{3}}=\mathrm{0}\end{array}$$

Where

$$\begin{array}{}\text{(9)}& \left\{\begin{array}{ll}{\mathit{\tau }}_{\mathrm{1}}& ={\mathit{\lambda }}_{\mathrm{1}}\sin \mathit{\beta }\sin {\mathit{\theta }}_{\mathrm{D}}-{\mathit{\lambda }}_{\mathrm{2}}\sin \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{D}}\\ {\mathit{\tau }}_{\mathrm{2}}& =-{\mathit{\lambda }}_{\mathrm{1}}\sin \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{D}}\sin {\mathit{\varphi }}_{\mathrm{D}}-{\mathit{\lambda }}_{\mathrm{2}}\sin \mathit{\beta }\sin {\mathit{\theta }}_{\mathrm{D}}\sin {\mathit{\varphi }}_{\mathrm{D}}\\ & +{\mathit{\lambda }}_{\mathrm{3}}\sin \mathit{\beta }\cos {\mathit{\varphi }}_{\mathrm{D}}\\ {\mathit{\tau }}_{\mathrm{3}}& ={\mathit{\lambda }}_{\mathrm{1}}\cos \mathit{\beta }\cos {\mathit{\theta }}_{\mathrm{D}}\cos {\mathit{\varphi }}_{\mathrm{D}}+{\mathit{\lambda }}_{\mathrm{2}}\cos \mathit{\beta }\sin {\mathit{\theta }}_{\mathrm{D}}\cos {\mathit{\varphi }}_{\mathrm{D}}\\ & +{\mathit{\lambda }}_{\mathrm{3}}\cos \mathit{\beta }\sin {\mathit{\varphi }}_{\mathrm{D}}-\cos \mathit{\eta }\\ {\mathit{\lambda }}_{\mathrm{1}}& =\sin \mathit{\alpha }\sin {\mathit{\theta }}_{\mathrm{1}}\sin {\mathit{\theta }}_{\mathrm{A}}+\cos \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{A}}\cos {\mathit{\varphi }}_{\mathrm{A}}\\ & -\sin \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{1}}\cos {\mathit{\theta }}_{\mathrm{A}}\sin {\mathit{\varphi }}_{\mathrm{A}}\\ {\mathit{\lambda }}_{\mathrm{2}}& =-\sin \mathit{\alpha }\sin {\mathit{\theta }}_{\mathrm{1}}\cos {\mathit{\theta }}_{\mathrm{A}}+\cos \mathit{\alpha }\sin {\mathit{\theta }}_{\mathrm{A}}\cos {\mathit{\varphi }}_{\mathrm{A}}\\ & -\sin \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{1}}\sin {\mathit{\theta }}_{\mathrm{A}}\sin {\mathit{\varphi }}_{\mathrm{A}}\\ {\mathit{\lambda }}_{\mathrm{3}}& =\cos \mathit{\alpha }\sin {\mathit{\varphi }}_{\mathrm{A}}+\sin \mathit{\alpha }\cos {\mathit{\theta }}_{\mathrm{1}}\cos {\mathit{\varphi }}_{\mathrm{A}}\end{array}\right.\end{array}$$

According to above equations, the output angle ${\mathit{\theta }}_{\mathrm{3}}^{+}$ obtains two roots that is ${\mathit{\theta }}_{\mathrm{3}}^{+}$ and ${\mathit{\theta }}_{\mathrm{3}}^{-}$, which can be calculated as

$$\begin{array}{}\text{(10)}& \left\{\begin{array}{l}{\mathit{\theta }}_{\mathrm{3}}^{+}\left({\mathit{\theta }}_{\mathrm{1}}\right)=\mathrm{2}\arctan \left({\displaystyle \frac{{\mathit{\tau }}_{\mathrm{1}}+\sqrt{{\mathit{\tau }}_{\mathrm{1}}^{\mathrm{2}}+{\mathit{\tau }}_{\mathrm{2}}^{\mathrm{2}}-{\mathit{\tau }}_{\mathrm{3}}^{\mathrm{2}}}}{{\mathit{\tau }}_{\mathrm{2}}-{\mathit{\tau }}_{\mathrm{3}}}}\right)\\ {\mathit{\theta }}_{\mathrm{3}}^{-}\left({\mathit{\theta }}_{\mathrm{1}}\right)=\mathrm{2}\arctan \left({\displaystyle \frac{{\mathit{\tau }}_{\mathrm{1}}-\sqrt{{\mathit{\tau }}_{\mathrm{1}}^{\mathrm{2}}+{\mathit{\tau }}_{\mathrm{2}}^{\mathrm{2}}-{\mathit{\tau }}_{\mathrm{3}}^{\mathrm{2}}}}{{\mathit{\tau }}_{\mathrm{2}}-{\mathit{\tau }}_{\mathrm{3}}}}\right)\end{array}\right.\end{array}$$

For each location of driving link corresponded to four specified coupler orientations, if all actual output angles θ₃ are generated by the same subformula of Eq. (10), the linkage solution will be valid. If not, a branch defect will appear in the motion of spherical 4R linkage. Specially, according to the coordinates of point C and D, the actual output angles θ₃ can be calculated by the method of Eq. (4).

4.3 Order defect discrimination

The order defect occurs when the coupler link can't pass the given four-orientations with the order of 1–2–3–4 or 4–3–2–1. Obviously, if the driving link is rotated fully, the geometrical criterion ${\mathit{\theta }}_{\mathrm{21}}<{\mathit{\theta }}_{\mathrm{31}}<{\mathit{\theta }}_{\mathrm{41}}$ or ${\mathit{\theta }}_{\mathrm{21}}>{\mathit{\theta }}_{\mathrm{31}}>{\mathit{\theta }}_{\mathrm{41}}$ will be used to solve this problem according to the definition of order defect, where θ_ji represents the rotation angle of driving link revolves from orientation ABⁱ (which is corresponded to coupler orientation i) to AB^j (which is corresponded to coupler orientation j) round the axis OA, as shown in Fig. 5.

Figure 5 Driving link revolved round the pole axis of sphere.

In order to obtain rotation angle θ_ji, the normal vector of plane OABⁱ (which is denoted by n₁) and OAB^j (which is denoted by n₂) are expressed as $n_{\mathrm{1}}=\mathrm{OA}\times {\mathrm{OB}}^i$, $n_{\mathrm{2}}=\mathrm{OA}\times {\mathrm{OB}}^j$. Based on the conclusions, the normal vector of axis OA in plane OABⁱ (which is denoted by f₁) and plane OAB^j (which is denoted by f₂) are expressed as $f_{\mathrm{1}}=n_{\mathrm{1}}\times \mathrm{OA}$, $f_{\mathrm{2}}=n_{\mathrm{2}}\times \mathrm{OA}$. Define $\mathrm{OA}=[x_{\mathrm{1}},y_{\mathrm{1}},z_{\mathrm{1}}]$, and $f_{\mathrm{1}}\times f_{\mathrm{2}}=[x_{\mathrm{2}},y_{\mathrm{2}},z_{\mathrm{2}}]$. The rotation angle of driving link θ_ji is calculated as

$$\begin{array}{}\text{(11)}& \left\{\begin{array}{cc}{\mathit{\theta }}_{ji}=\arccos {\displaystyle \frac{f_{\mathrm{1}}\cdot f_{\mathrm{2}}}{\left|f_{\mathrm{1}}\right|\left|f_{\mathrm{2}}\right|}}& z_{\mathrm{1}}\cdot z_{\mathrm{2}}>\mathrm{0}\\ {\mathit{\theta }}_{ji}=\mathrm{2}\mathit{\pi }-{\mathit{\theta }}_{ji}& \mathrm{others}\end{array}\right.\end{array}$$

On the contrary, if the driving link is rotated partially, the geometrical criterion which mentioned above (that is ${\mathit{\theta }}_{\mathrm{21}}<{\mathit{\theta }}_{\mathrm{31}}<{\mathit{\theta }}_{\mathrm{41}}$ or ${\mathit{\theta }}_{\mathrm{21}}>{\mathit{\theta }}_{\mathrm{31}}>{\mathit{\theta }}_{\mathrm{41}}$) will be invalid. For example, the order 4–1–2–3 is incorrect because of the driving link can't move from one moving interval to other, as shown in Fig. 6. In this case, the geometrical criterion is modified as ${\mathit{\theta }}_{\mathrm{1}e}<{\mathit{\theta }}_{\mathrm{2}e}<{\mathit{\theta }}_{\mathrm{3}e}<{\mathit{\theta }}_{\mathrm{4}e}$ or ${\mathit{\theta }}_{\mathrm{1}e}>{\mathit{\theta }}_{\mathrm{2}e}>{\mathit{\theta }}_{\mathrm{3}e}>{\mathit{\theta }}_{\mathrm{4}e}$, where the subscript letter e represents the arbitrary limited orientations of the driving link.

Figure 6 Order defect discrimination.

4.4 Pyramid structure of the motion defect calculation

The method of defect discrimination cannot be considered in isolation. For example, the calculation on the order defect discrimination is meaningless for a 4R linkage with circle defect because the coupler of the 4R linkage cannot pass through all specified orientations. In order to solve this problem, a pyramid structure is defined in this program package. The circuit defect locates at the bottom of pyramid, the branch defect locates at the middle of pyramid, and the order defect locates at the top of pyramid. Under this definition, both the circuit defect and branch defect should be checked first before the calculation of order defect discrimination. Similarly, the circuit defect should be checked first before the calculation of branch defect discrimination.

4.5 Linkage type classification

For a number of synthesis tasks, the 4R linkage with specified linkage type should be satisfied. According to the method proposed by Murray, there are 8 kinds of linkage type for a spherical 4R linkage, as shown in Table 1, where $K_{\mathrm{1}}=\mathit{\gamma }-\mathit{\alpha }+\mathit{\eta }-\mathit{\beta }$, $K_{\mathrm{2}}=\mathit{\gamma }-\mathit{\alpha }-\mathit{\eta }+\mathit{\beta }$, $K_{\mathrm{3}}=\mathit{\eta }+\mathit{\beta }-\mathit{\gamma }-\mathit{\alpha }$, $K_{\mathrm{4}}=\mathrm{2}\mathit{\pi }-\mathit{\eta }-\mathit{\beta }-\mathit{\gamma }-\mathit{\alpha }$. Specially, at least one link of a 4R linkage will be rotated fully if K₁K₂K₃K₄>0, which is the Grashof's condition. Both calculations on the linkage type and Grashof's condition are supported by program package.

Table 1 Classification of the linkage type (Murray and Larochelle, 1998).

4.6 Performance evaluation and solutions domain analysis

A standard function interface on linkage performance calculation is provided for the designer in the program package. The input vector of this function is a data structure to describe the coordinate information of linkage solution at each specified orientation, and the output variable of this function is a real number to indicate the linkage performance. If there are more than one evaluating indicators for the linkage performance, a weight method will be used. Of course, all the functions in the program package can be called in this function.

Actually, a default function template is provided for the designer. The evaluating indicator which defined in this function template is termed as “Braking angle”. For a valid 4R linkage solution which driving link is rotated partially, if the rotation angle of driving link from the specified orientation 1 (or 4) to the nearest limited orientation is too small, the driving link will be difficult to stop before the linkage moves to the limited orientation for the reason of the high rotation speed or the inaccurate displacement control of the motor.

If N center points and circle points are generated after the synthesis calculation, N×N linkage solutions in total will be obtained. Based on the calculation on linkage performance evaluation, each linkage solution will be analyzed and the distribution of the linkage performance evaluation in the solution domain will be obtained. The calculation on the solutions domain analysis is visualized by a solutions map, which is introduced later.

5 Synthesis visualization

5.1 Task visualization

Once the four coupler orientations is specified, a figure with four orientations, the center points and the circle points plotted in a sphere are generated. The designer can adjust task parameters according to the figure. The specified coupler orientations are defined by four equilong arcs in the sphere, and the longitude and latitude of the coupler orientation are determined by the starting point of arc, and the roll angle of coupler orientation is determined by the roll angle of arc. Specially, in order to plot the arc in the surface of sphere smoothly, the arc is divided into several line segments through an interpolation method.

5.2 Solutions map

For N center points and circle points, all linkage solutions can be mapped into a three-dimensional coordinate system O−XYZ, that is the visualization of solutions domain analysis mentioned above. In this coordinate system, N×N cubes are constructed based on the O−XY plane, and each linkage solution defined by the center point (or circle point) i and the center point (or circle point) j is described by a cube $L(i,j,h)$, as shown in Fig. 7, where i is the x offset of cube L, j is the y offset of the cube L, h is the z offset of the cube L, i and j are determined by selected center points (or circle points), and h is determined by linkage performance. In order to describe the linkage type of solution, each cube in the solutions map is color-coded according to linkage type.

Figure 7 Constructions of solutions map.

5.3 Solutions visualization

The key problem of solutions visualization is the motion analysis and calculation of coupler curve. Figure 3 shows that the output angle θ₃ is calculated by the Eq. (10) aimmed at each input angle θ₁. Due to the coordinates of point A and D of spherical 4R linkage are given, the coordinates of point B and C can be obtained according to the Eqs. (5) and (6). Then, the rotation angle of link BC (that is θ₂) can be calculated as

$$\begin{array}{}\text{(12)}& \left\{\begin{array}{l}\sin {\mathit{\theta }}_{\mathrm{2}}={\displaystyle \frac{\cos {\mathit{\varphi }}_{\mathrm{C}}\sin \left({\mathit{\theta }}_{\mathrm{B}}-{\mathit{\theta }}_{\mathrm{C}}\right)}{\sin \mathit{\eta }}}\\ \cos {\mathit{\theta }}_{\mathrm{2}}={\displaystyle \frac{\sin {\mathit{\varphi }}_{\mathrm{C}}-\sin {\mathit{\varphi }}_{\mathrm{B}}\cos \mathit{\eta }}{\cos {\mathit{\varphi }}_{\mathrm{B}}\sin \mathit{\eta }}}\end{array}\right.\end{array}$$

Correspondingly, the rotation angle revolved from the link BC to BP (that is δ) is obtained by the rotation method introduced in Sect. 3.3. Obviously, the angle δ is a constant value and calculated from the coordinate information of spherical 4R linkage at any specified orientations. On this basis, the rotation angle of link BP (denoted by θ₅) is “θ₂+δ”, and the coordinate of point P is expressed as

$$\begin{array}{}\text{(13)}& \left\{\begin{array}{l}\sin {\mathit{\varphi }}_{\mathrm{P}}=\cos \mathit{\kappa }\sin {\mathit{\varphi }}_{\mathrm{B}}+\sin \mathit{\kappa }\cos {\mathit{\theta }}_{\mathrm{5}}\cos {\mathit{\varphi }}_{\mathrm{B}}\\ \sin {\mathit{\theta }}_{\mathrm{P}}={\displaystyle \frac{{\displaystyle \begin{array}{c}-\sin \mathit{\kappa }\sin {\mathit{\theta }}_{\mathrm{5}}\cos {\mathit{\theta }}_{\mathrm{B}}+\cos \mathit{\kappa }\sin {\mathit{\theta }}_{\mathrm{B}}\cos {\mathit{\varphi }}_{\mathrm{B}}\\ -\sin \mathit{\kappa }\cos {\mathit{\theta }}_{\mathrm{5}}\sin {\mathit{\theta }}_{\mathrm{B}}\sin {\mathit{\varphi }}_{\mathrm{B}}\end{array}}}{\cos {\mathit{\varphi }}_{\mathrm{P}}}}\\ \cos {\mathit{\theta }}_{\mathrm{P}}={\displaystyle \frac{{\displaystyle \begin{array}{c}\sin \mathit{\kappa }\sin {\mathit{\theta }}_{\mathrm{5}}\sin {\mathit{\theta }}_{\mathrm{B}}+\cos \mathit{\kappa }\cos {\mathit{\theta }}_{\mathrm{B}}\cos {\mathit{\varphi }}_{\mathrm{B}}\\ -\sin \mathit{\kappa }\cos {\mathit{\theta }}_{\mathrm{5}}\cos {\mathit{\theta }}_{\mathrm{B}}\sin {\mathit{\varphi }}_{\mathrm{B}}\end{array}}}{\cos {\mathit{\varphi }}_{\mathrm{P}}}}\end{array}\right.\end{array}$$

According to the calculation above, the point coordinates of the spherical 4R linkage under each input angle are calculated, and the coordinate data is applied to conduct the animation of linkage motion and plot the coupler curve. Due to there are two solutions in the Eq. (10) corresponding to the two assembly configurations of spherical 4R linkage driven link, the coupler curve is generated at twice. In order to distinguish the branch of coupler curve, different colors are used to represent different branches. Specially, for a spherical 4R linkage which driving link rotated partially, each solution of Eq. (10) should be checked because there is no real root for some input angles of driving link.

6 Expansion of the program package on planar 4R linkage

The synthesis on both spherical and planar 4R linkages is supported by the program package. The calculation on spherical problem mentioned above is mainly based on the theory of spherical trigonometry (Ratcliffe, 2006). Actually, the geometrical principle of planar 4R linkage is similar to a spherical one. As the Fig. 2 shows, the spherical 4R linkage only move along the surface of a sphere and all the axes of hinges passes through the center of sphere. If the diameter of sphere is infinite, the spherical 4R linkage ABCD will be transformed into a planar one.

For the synthesis of planar 4R linkage, the method mentioned in Sect. 2 which is presented by Shirazi (2005) is applicable, it is based on the classic burmester theory and cannot be described here.

For the discrimination of linkage circuit defect, the difference between spherical and planar problems is the calculation of limited orientations of driving link, which can be obtained by the cosine theorem of planar triangle. Define $a^{\prime }=\left|\mathrm{BC}\right|+\left|\mathrm{CD}\right|$, $b^{\prime }=\left|\mathrm{AD}\right|$, $c^{\prime }=\left|\mathrm{AB}\right|$, if the configuration of planar triangle is A(BC)D, the limited angles of driving link which is denoted by ${\mathit{\theta }}_{\mathrm{1}}^{\mathrm{1}\ast }$ and ${\mathit{\theta }}_{\mathrm{1}}^{\mathrm{2}\ast }$ could be calculated as

$$\begin{array}{}\text{(14)}& \left\{\begin{array}{l}{\mathit{\theta }}_{\mathrm{1}}^{\mathrm{1}\ast }=\arccos \left({\displaystyle \frac{{c^{\prime }}^{\mathrm{2}}+{b^{\prime }}^{\mathrm{2}}-{a^{\prime }}^{\mathrm{2}}}{\mathrm{2}{c}^{\prime }{b}^{\prime }}}\right)+\arctan \mathrm{2}\left(\sin {\mathit{\theta }}_{\mathrm{4}},\cos {\mathit{\theta }}_{\mathrm{4}}\right)\\ {\mathit{\theta }}_{\mathrm{1}}^{\mathrm{2}\ast }=-\arccos \left({\displaystyle \frac{{c^{\prime }}^{\mathrm{2}}+{b^{\prime }}^{\mathrm{2}}-{a^{\prime }}^{\mathrm{2}}}{\mathrm{2}{c}^{\prime }{b}^{\prime }}}\right)+\arctan \mathrm{2}\left(\sin {\mathit{\theta }}_{\mathrm{4}},\cos {\mathit{\theta }}_{\mathrm{4}}\right)\end{array}\right.\end{array}$$

Similarly, if the configuration of planar triangle is A(B^′C^′)D, the Eq. (14) will be applicable as $a=\left|\mathrm{BC}\right|-\left|\mathrm{CD}\right|$. If it is not all input angles of driving link corresponding to the specified coupler orientations locates at the same moving interval, the linkage solution will be invalid. Specially, the validity of Eq. (14) is determined by the validity of expression $v=(c^{\mathrm{2}}+b^{\mathrm{2}}-a^{\mathrm{2}})/\mathrm{2}cb$. If the solution of Eq. (14) is valid, the condition $v\in [-\mathrm{1},\mathrm{1}]$ should be satisfied.

For the discrimination of linkage branch defect, the Eqs. (8) and (10) are applicable, the difference between spherical and planar problems is the expression of τ₁, τ₂, τ₃, where

$$\begin{array}{}\text{(15)}& \left\{\begin{array}{ll}{\mathit{\tau }}_{\mathrm{1}}& =\mathrm{2}\left|\mathrm{AB}\right|\left|\mathrm{CD}\right|\sin {\mathit{\theta }}_{\mathrm{1}}-\mathrm{2}\left|\mathrm{AD}\right|\left|\mathrm{CD}\right|\sin {\mathit{\theta }}_{\mathrm{4}}\\ {\mathit{\tau }}_{\mathrm{2}}& =\mathrm{2}\left|\mathrm{AB}\right|\left|\mathrm{CD}\right|\cos {\mathit{\theta }}_{\mathrm{1}}-\mathrm{2}\left|\mathrm{AD}\right|\left|\mathrm{CD}\right|\cos {\mathit{\theta }}_{\mathrm{4}}\\ {\mathit{\tau }}_{\mathrm{3}}& ={\left|\mathrm{BC}\right|}^{\mathrm{2}}-{\left|\mathrm{AB}\right|}^{\mathrm{2}}-{\left|\mathrm{CD}\right|}^{\mathrm{2}}-{\left|\mathrm{AD}\right|}^{\mathrm{2}}\\ & +\mathrm{2}\left|\mathrm{AB}\right|\left|\mathrm{AD}\right|\cos \left({\mathit{\theta }}_{\mathrm{1}}-{\mathit{\theta }}_{\mathrm{4}}\right)\end{array}\right.\end{array}$$

If it was not all actual output angles θ₃ corresponding to four specified coupler orientations that could be generated by the same subformula of Eq. (14), the linkage solution would be invalid.

For the discrimination of linkage order defect, the geometrical criterion ${\mathit{\theta }}_{\mathrm{21}}<{\mathit{\theta }}_{\mathrm{31}}<{\mathit{\theta }}_{\mathrm{41}}$ or ${\mathit{\theta }}_{\mathrm{21}}>{\mathit{\theta }}_{\mathrm{31}}>{\mathit{\theta }}_{\mathrm{41}}$ for linkage with rotating fully driving link and geometrical criterion ${\mathit{\theta }}_{\mathrm{1}e}<{\mathit{\theta }}_{\mathrm{2}e}<{\mathit{\theta }}_{\mathrm{3}e}<{\mathit{\theta }}_{\mathrm{4}e}$ or ${\mathit{\theta }}_{\mathrm{1}e}>{\mathit{\theta }}_{\mathrm{2}e}>{\mathit{\theta }}_{\mathrm{3}e}>{\mathit{\theta }}_{\mathrm{4}e}$ for linkage with rotating partially driving link are applicable. The difference between spherical and planar problems is the vectors f₁ and f₂ in the Eq. (11), where f₁=ABⁱ and f₂=AB^j.

For the classification of linkage type and the condition of grashof, the method proposed by Martin and Murray (2002) is applicable to planar 4R linkage. For the performance evaluation, the “Braking angle” of planar 4R linkage is obtained from the calculation on planar vectors. Similarly, for the synthesis visualization of planar 4R linkage, the method mentioned in Sect. 4 is applicable, and the calculation is based on the principle of planar trigonometry.

Figure 8 Center points and circle points generated by the program package.

Figure 9 Solutions map generated by the program package.

Table 2 Specified coupler orientations for a spherical 4R linkage.

Table 3 Solutions information.

7 Examples

7.1 Spherical 4R linkage synthesis

In this example, it is necessary to synthesize a spherical 4R linkage to guide a camera to pass four specified orientations without circuit defect, branch defect and order defect, as shown in Table 2.

Figure 10 Solutions visualization.

Firstly, the synthesis program is called to generate center points and circle points, the result is shown in Fig. 8. Then, the solutions map is generated, the result is shown in Fig. 9. As shown in Fig. 9, 91.7 % of solutions are erased from the map because of the linkage motion defect, and there are no valid linkages with type of crank-rocker, rocker crank, grashof double rocker and 0π double rocker satisfied the synthesis task. Additionally, all linkage solutions which driving link rotate partially have a bad “Braking angle” less than 5^∘, therefore, only the solutions with a type of double crank need to be considered.

The program allows designer to select valid solutions on the solutions map directly, and this way is quick and efficient. As the Fig. 9 shows, a candidate solution is selected from region I of solutions map, and the solution information is shown in Table 3. Similarly, another three linkage solutions selected from region II–IV of solutions map are erased because of motion defect, and the solution information is shown in Table 3. In order to verify the results, the plotting program is called to generate coupler curve of four linkages, as shown in Fig. 10. Obviously, only the linkage solution I is feasible.

In this calculation, 7396 groups of linkage solutions are synthesized and analyzed. When the program is running on the computer, it takes 4.48 s to generate and display the solutions map at one time. In the process of linkage solution selection, it only takes 0.37 s to generate and display the coupler curve of the selected linkage solution each time.

7.2 Planar 4R linkage synthesis

In this example, it is necessary to synthesize a planar 4R linkage to guide the bucket of the skid steer loader to pass four specified orientations without circuit defect, branch defect and order defect, as shown in Table 4.

Table 4 Specified coupler orientations for a planar 4R linkage.

Firstly, the synthesis program is called to generate center points and circle points, the result is shown in Fig. 11. Then, the solutions map is generated, the result is shown in Fig. 12. It is noted that the installation position of the 4R linkage on the skid steer loader is limited, which becomes a strong constraint condition for the rationality of the solution. In order to provide enough candidate solutions, this example no longer calculates the transmission performance of the linkage solution, but loads the picture of the skid steer loader with reduced proportion into the drawing area to evaluate the reasonableness of the installation position of the linkage solution.

Figure 11 Center points and circle points generated by the program package.

Figure 12 Solutions map generated by the program package.

Table 5 Solutions information.

As shown in Fig. 12, 70.95 % of solutions are erased from the map because of the linkage motion defect, and there are no valid linkages with type of crank-rocker, rocker crank, double crank, 00 double rocker and ππ double rocker satisfied the synthesis task. Six candidate solutions are selected from region I–VI of solutions map, and the solution information is shown in Table 5. In order to verify the results, the plotting program is called to generate coupler curve of six linkages. As shown in Fig. 13, the installation position of one rocker in solution II interferes with the tire of the skid steer loader, and the installation position of one rocker in solution III–VI is outside the body of the skid steer loader. Obviously, only the installation position of the rocker with solution I is reasonable and meets the requirements of linkage synthesis.

Figure 13 Solutions visualization.

In this calculation, 19 600 groups of linkage solutions are synthesized and analyzed. When the program is run on the computer, it takes 12.73 s to generate and display the solutions map at one time. In the process of linkage solution selection, it only takes 0.24 s to generate and display the coupler curve of the selected linkage solution each time.

8 Concluding remarks

In this paper, a synthesis program package is developed to find a satisfied planar or spherical linkage solution automatically according to burmester theory. The main conclusions of this study are as follows:

1. The program package developed in this paper has fast calculation speed and high synthesis efficiency. The time required for synthesis calculation depends on the segmentation accuracy of buermester center curve, that is, the rotation angle increment of driving link of 4R linkage constructed by image poles when generating each center point. The more precise the burmester center curve is segmented, the more solutions are generated and the longer the calculation time is. In the two examples given in this paper, the number of linkage solutions analyzed is 7396 and 19 600 respectively, and the total time required for solutions map generation and display is 4.48 and 12.73 s respectively. The generation of the solutions map is one-time, after which the selection and analysis of linkage solution only takes less than 0.5 s each time. In addition, the solutions map proposed in this study extends the dimension of traditional type map, so that it can also express the transmission performance of linkage solutions while expressing linkage types. This method further avoids the blindness of linkage synthesis.

2. The proposed defect discrimination algorithm is not only correct but also logical. In the synthetic example of spherical 4R linkage, by comparing the result of program discrimination with that of coupler curve analysis, it can be seen that the algorithm in this paper can correctly discriminate the defect of circuit, branch and order. Moreover, this paper proposes a pyramid structure for defect discrimination. This is obviously more logical. For example, when four given orientations are located in different circuits and branches, the coupler curve of 4R linkage can not pass through all four orientations at all. It is not logical to discuss whether the coupler link can pass through four orientations orderly.

3. The program developed in this paper is more flexible than the previous software. Like the toolbox functions of MATLAB, the program contains a series of functions from the bottom graphics calculation to the generation of human-computer interaction interface. On this basis, users can not only directly use the given human-computer interaction interface for general synthesis calculation, but also reorganize these functions to solve more complex problems.

Considering that the precise orientation synthesis problem and the function synthesis problem can be transformed each other, how to extend the research method in this paper, especially the defect discrimination method, to the function synthesis field of spherical and planar 4R linkage is the next research direction of this paper.