Valid Initial Pose of Any Fourbar Linkage

Stefan Gössner

August 23, 2017

Department of Mechanical Engineering, University of Applied Sciences, Dortmund, Germany.

Keywords: Fourbar Linkage; Brahmagupta's formula; cyclic quadrilateral; initial pose;


Abstract

For simulating fourbar linkages a valid initial pose is needed. Finding such a pose is sometimes a non-trivial task. In this paper an elementary method for specifying a valid initial pose of any fourbar linkage is suggested.

Introduction

Analysis and simulation of fourbar mechanisms is nothing new and done for decades. The initial pose of the fourbar linkage is usually nothing, we think about very often. Solving the task to animate one cycle of a crank-rocker, we pragmatically let the crank angle start at 0° and stop at 360° (with respect to whatever). But this implies, that we know,

  • that the fourbar satisfies the Grashof condition.
  • which link is the fully rotatable link.

McCarthy classifies eight basic types with the fourbar [1]. For each type we have to

  • check the Grashof condition.
  • identify the fully rotatable link, if any.
  • choose the driven input link.
  • specify an initial input angle for the input link.
  • set an end angle accordingly.

Now imagine a scenario, where you can graphical interactively modify the four link lengths. Then you must go through these steps with each length change. The critical step here would be to find an always valid start angle for your input link, if one is even uniquely identified.

The universal method suggested here for finding a valid initial pose is independent of above steps.

Prerequisites

The geometric shape of a fourbar mechanisms is a quadrilateral in every pose. The area enclosed by that quadrilateral becomes extreme, when it is a cyclic quadrilateral. Brahmagupta's formula helps to calculate the aera AA enclosed by that cyclic quadrilateral as well as the radius RR of the circle, all vertices are lying on [2].

Let's take as an assumption, that the fourbar's pose corresponding to a cyclic quadrilateral is always a valid pose.

The maximal area of a cyclic quadrilateral is determined by

A=(sa)(sb)(sc)(sd)A = \\sqrt{(s-a)(s-b)(s-c)(s-d)} (1)

with the semiperimeter

s=a+b+c+d2.s = \\dfrac{a+b+c+d}{2}\\,. (2)

The radius RR of the corresponding circle is then [2]

R=14A(ab+cd)(ac+bd)(ad+bc). R = \\frac{1}{4A}\\sqrt{(ab+cd)(ac+bd)(ad+bc)}\\,. (3)

Please note, that Brahmagupta's Formula is of no help to draw the cyclic quadrilatral's shape, i.e. providing the side vectors or the angles between the sides.

The qadrilateral shape corresponding to the fourbar pose is to be determined with given side lengths a,b,c,da, b, c, d only.

The vectorial method applied is making use of an orthogonal operator. The orthogonal vector to a\\bold a is denoted here as a~\\tilde\\bold a, while in other publications it is often written as a\\bold a^{\\perp} [3].

Calculate the Fourbar Pose

Fig.1: Geometry of the cyclic qadrilateral

We start analog to the proof of Brahmagupta's Formula with the sum of the area of the two triangles in Fig.1 using half area product each [2] .

A=12(a~b+c~d)andA=12(b~c+d~a) A = \\frac{1}{2}(\\tilde\\bold a \\bold b + \\tilde\\bold c \\bold d)\\quad and \\quad A = \\frac{1}{2}(\\tilde\\bold b \\bold c + \\tilde\\bold d \\bold a)(4)

Both angles in the triangles opposite to each diagonal sum up to 180180^{\\circ}, i.e.

α+γ=180,sinα=sinγ,cosα=cosγ\\alpha + \\gamma = 180^{\\circ},\\,\\sin\\alpha = \\sin\\gamma,\\,\\cos\\alpha = -\\cos\\gamma
β+δ=180,sinβ=sinδ,cosβ=cosδ\\beta + \\delta = 180^{\\circ},\\,\\sin\\beta = \\sin\\delta,\\,\\cos\\beta = -\\cos\\delta

The perp dot products of the side vectors correspond to the sines of the enclosed angles, so we can write

a~bab=c~dcdandb~cbc=d~aad\\frac{\\tilde\\bold a \\bold b}{ab} = \\frac{\\tilde\\bold c \\bold d}{cd}\\quad and \\quad \\frac{\\tilde\\bold b \\bold c}{bc} = \\frac{\\tilde\\bold d \\bold a}{ad}

Introducing these expressions to equations (4) we get for the sine's

c~dcd=2Aab+cdandd~aad=2Aad+bc\\frac{\\tilde\\bold c \\bold d}{cd} = \\frac{2A}{ab + cd}\\quad and \\quad\\frac{\\tilde\\bold d \\bold a}{ad} = \\frac{2A}{ad + bc}(5)

The dot products of the side vectors correspond to the cosines of the enclosed angles, so we can write

abab=cdcdandbcbc=daad.\\frac{\\bold a \\bold b}{ab} = -\\frac{\\bold c \\bold d}{cd}\\quad and \\quad \\frac{\\bold b \\bold c}{bc} = -\\frac{\\bold d \\bold a}{ad}\\,.(6)

Now expressing the diagonal vectors by its two triangle side vectors each

p=a+b=(c+d)andq=b+c=(d+a)\\bold p = \\bold a + \\bold b = -(\\bold c + \\bold d)\\quad and \\quad\\bold q = \\bold b + \\bold c = -(\\bold d + \\bold a)

and squaring leads us to

a2+2ab+b2=c2+2cd+d2andb2+2bc+c2=d2+2da+a2.a^2 + 2\\bold a\\bold b + b^2 = c^2 + 2\\bold c\\bold d + d^2\\quad and \\quad b^2 + 2\\bold b\\bold c + c^2 = d^2 + 2\\bold d\\bold a + a^2\\,.

Reintroducing these expression to equations (6) yields the cosines

cdcd=a2+b2c2d2ab+cdanddaad=a2+b2+c2d2ad+bc.\\frac{\\bold c \\bold d}{cd} = \\frac{a^2 + b^2 - c^2 - d^2}{ab + cd}\\quad and \\quad \\frac{\\bold d \\bold a}{ad} = \\frac{-a^2 + b^2 + c^2 - d^2}{ad + bc}\\,.(7)

With this we have all ingredients to determine the side vectors a,b,c,d\\bold a,\\bold b,\\bold c, \\bold d. Let the base vector d\\bold d be aligned with the positive x-axis. The side vectors then finally result using expressions (5) and (7).

a=a(daadexd~aadey)c=c(cdcdex+c~dcdey)d=dexb=(a+c+d) \\begin{matrix} \\bold a & = & a\\,(\\dfrac{\\bold d \\bold a}{ad} \\bold e\_x - \\dfrac{\\tilde\\bold d \\bold a}{ad} \\bold e\_y) \\\\ \\\\ \\bold c & = & c\\,(\\dfrac{\\bold c \\bold d}{cd} \\bold e\_x + \\dfrac{\\tilde\\bold c \\bold d}{cd} \\bold e\_y) \\\\ \\\\ \\bold d & = & d\\, \\bold e\_x \\\\ \\\\ {\\bold b} & = & -(\\bold a + \\bold c + \\bold d) \\end{matrix} (8)

Please note, that trigonometric functions were used here only for clarification. The resulting vector equations (8) describing the initial pose are trigonometry-free.

References

[1] J.M. McCarthy, G.S. Soh, Geometric Design of Linkages, Springer, 2010.
[2] Brahmagupta's Formula, S. Gössner.
[3] Cross Product Considered Harmful, S.Gössner.