Open access peer-reviewed chapter

# Optimization of Kinematics of Inclined Swinging Pin

By Stanislav Barton

Submitted: November 15th 2018Reviewed: April 29th 2019Published: June 19th 2019

DOI: 10.5772/intechopen.86589

## Abstract

The inclined swinging pin simply transforms the rotary motion into a rotational oscillation. It consists of three components, a rotating inclined tenon, a crossbeam, and a stirrup. The pitch angle of the inclined tenon relative to the x-axis of its rotation is decisive for the kinematics of this joint. Too small inclination angle will result in small amplitude of oscillation; too much inclination will lead to an impulse dynamic stress that can damage the pin. The optimal angle of inclination can be determined from the mathematical model, created in the Maple environment, which uses linear algebra resources to describe the behavior of the joint. Vectors of coordinates of the critical points are created for each joint component. Furthermore, transformation matrices are created which describe the behavior of the corresponding kinematic chain in relation to the inclination angle of the inclined tenon and time.

### Keywords

• mathematical modeling in Maple
• transformation matrix
• angular velocity
• angular acceleration
• optimization

## 1. Introduction

The inclined swinging pin is used for the simple conversion of rotational motion into oscillation in many agricultural machines. For example, in combine harvesters, it is used to drive a mowing bar that cuts stems. It is composed of three basic components: an inclined tenon—red, a crossbeam—green, and a stirrup—blue, see Figure 1. We study the movement of individual components in a rectangular coordinate system, see Figure 1.

Similar mechanisms are described, for example, in the book [1] that is available in an electronic version [2]. Many animations can be found on the Web, if you use keywords rotation to oscillation mechanism, Google will offer you over 20,000 links.

## 2. The basic components

### 2.1 Tenon

The pin is a basic component. Its basic part is an oblique nose. Its axis of symmetry forms an angle ϕwith the axis of rotation of the whole pin, which is identical to the axis xof the coordinate system. The axis of symmetry of the nose passes through the beginning of the coordinate system. The spacer serves to secure the nose to the support shaft, see Figure 2.

### 2.2 Crossbeam and stirrup

The crossbeam is the middle member of the kinematic chain. The tenon nose is pushed through its center hole. As the nose rotates, the crossbar turns. The stirrup is attached to its side pins, and therefore the axis of the side pins is constantly in the plane z = 0, see Figure 3.

Depending on how the crossbeam swings, the stirrup oscillates. The driven device connects to the vertical axis of the stirrup, see Figure 4.

## 3. Deriving of transform matrices

We will derive the transformation matrices in the Maple13 program environment. The decisive element for which the transformation matrix will be derived is the crossbeam.

Allow us to specify axes of the base coordinate system in which tenon—Figures 1 and 2—are displayed as [x, y, z]. The axes of the moving coordinate system in which we plot the actual position of the crossbeam are denoted as [X, Y, Z]. Let us denote the transformation matrix from system [x, y, z] into system [X, Y, Z] as R_xyz. The individual rows of the R_xyzmatrix are the x, y, zcoordinates of the unit vectors in the X, Y, Zdirection.

The direction of the X-axis is entered as the direction of the tenon’s nose axis. The axis X, which is equal to the axis of the nose, creates the angle ϕwith the x-axis and generally the entire tenon is rotated by an angle ψaround the x-axis.

A detailed description of the properties of transformation matrices can be found, for example, in [3].

>restart; with(LinearAlgebra): with(plots): #Maple Start+Libraries

>R_z:=Matrix([[cos(phi),sin(phi),0],

[sin(phi),cos(phi),0],[0,0,1]]): #R_z - Rotation around z

>R_x:=Matrix([[1,0,0],[0,cos(psi),sin(psi)],

[0,-sin(psi),cos(psi)]]):]]): #R_x - Rotation around x

>R_xz:=R_x.R_z;

R_xzcosφsinφ0cosψsinφcosψcosφsinψsinψsinφsinψcosφcosψE1

>X:=convert(R_xz.<1,0,0>,list);

Xcosφcosψsinφsinψsinφ.E2

Equation (2) specifies the direction of the nose axis, thus X-axis vector in the coordinates [x, y, z]. Coordinates Yare calculated easily. We know that |Y| = 1 and Yzmust always be equal to 0 and that Xand Ymust be perpendicular to each other, i.e., X·Y = 0, the scalar product Xand Ymust be equal to 0. Therefore, we choose the vector coordinates Y = [Yx, Yy, 0] and Yxand Yyare calculated from the following equations:

>Y:=[Yx,Yy,0];

>simplify(allvalues(solve({E1,E2},{Yx,Yy}))[2],symbolic);

>Y:=subs(%,Y);

Ycosψsinφcosψ2cosψ2cosφ2+cosϕ2cosφcosψ2cosψ2cosφ2+cosϕ20.E3

Equations E1and E2have two solutions. In order for the XYZsystem to be oriented counterclockwise, a second solution has to be chosen. Coordinates Zare calculated from the conditions X·Z = 0, Y·Z = 0, and |Z| = 1.

>Z:=[Zx,Zy,Zz];

>simplify(allvalues(solve({E1,E2,E3},{Zx,Zy,Zz}))[1],symbolic);

>Z:=subs(%,Z);

Zsinψcosφsinφ%1sinψcosψsinφ2%1%1,where%1=cosψ2+cosψ2cosφ2+cosφ2.E4

From the vectors X, Y, Z, it is possible to create a transformation matrix M_XYZ, which will describe the orientation of the crossbeam depending on the angle ϕand the angle of rotation ψ.

>M_XYZ:=Matrix([X,Y,Z]);

M_XYZcosφcosψsinφsinψsinφcosψsinφ%1cosφ%10sinψsinφcosφ%1sinψcosψsinφ2%1%1.E5

Finally, it is possible to create a graph showing the curves described by the endpoints of the X, Y, and Zvectors. For instance, if ϕ=2π/15and ψ02πare selected, then the resulting graph is Figure 5. The corresponding animation is stored in the file An_Axes.avi, available to download in the Videos section.

## 4. Optimization

The orientation of the Yaxis can be derived from the matrix M_XYZ, see Eq. (5). The orientation of this axis indicates the angle of rotation of the stirrup about the zaxis.

>Y:=evalm([0,1,0].M_XYZ);

Ycosψsinφ%1cosφ%10.E6

The angle of rotation of the stirrup is equal to the angle between the yand Yaxes. Let us denote this angle as θ.

>theta:=arctan(rz[1]/rz[2]);

θarctancosψsinφcosφ.E7

The angle ψindicates the angle of rotation of the pin about the x-axis. We can therefore substitute:

>su:=psi=2*Pi*tau:

suψ=2πτ.E8

>theta:=subs(su,sin(phi)=cos(phi)*tan(phi),theta);

θarctancos2πτtanφ,E9

where τ = t/T = dimensionless time, t = time, and T = duration of one revolution = one period.

It is now possible to determine the angular velocity of the stirrup oscillations, ω=dθτ/, and the angular acceleration of the stirrup oscillations, ε=d2θτ/2.

>omega:=diff(theta,tau);

ω2sin2πτπtanφT1+cos2πτ2tanφ2E10

>epsilon:=normal(diff(omega,tau));

ε4cos2πτπ2tanφ1+cos2πt2tanφ2+2sin2πt2tanφ21+cos2πt2tanφ22E11

Now we plot the graph of ε(τ) for τ01.The values of ϕwill be chosen between 5 and 45° with a step of 5°, see Figure 6.

>EP:=plot([seq(epsilon,phi=expand(Pi/36*[$1..9]))], tau=0..1,color=red): EP; It is quite clear from Figure 6 that with an increase in angle ϕ—that is, the angle between the pin nose and the x-axis—the course of the angular acceleration of the oscillations begins to deviate sharply from the sinusoidal shape. Thus, large variations in angular acceleration will occur for large ϕvalues, resulting in rapid wear of the entire mechanism. This can lead to mechanical damage, because high angular acceleration values require high torque forces. Conversely, small values of the angle ϕwill lead only to a low oscillation amplitude, which may not be sufficient for the proper operation of the driven device. It is therefore necessary to find the optimum angle value ϕ. After studying the graphs in Figure 6, the value at which the angular acceleration for τ = 1/2 reaches the local maximum possible value appears optimal. If the value of ϕfurther increases, then this local maximum becomes the local minimum. This means that the optimum value of the angle ϕmust be the one at which the transition from local maximum to local minimum occurs. From a mathematical point of view, the following condition must be satisfied: d2ετφdτ2τ=12=0.E12 The value of the angle ϕcan be easily determined from Eq. (12) using Maple. >Ett:=normal(diff(epsilon,tau,tau)); Ett16cos2πτπ4tanφ(1+24sin2πτ4tanφ6cos2πτ2+8sin2πτ2tanφ4cos2πτ220sin2πτ2tanφ2+9cos2πτ4tanφ4+5cos2πτ6tanφ624sin2πτ4tanφ4+3cos2πτ2tanφ2+28sin2πτ2tanφ6cos2πτ4)/1+cos2πτ2tanφ24E13 >Ett2:=eval(subs(tau=1/2,Ett)); Ett216π4tanφ1+9tanφ4+5tanφ6+3tanφ21+tanφ24E14 >Sol:=[solve(Ett2,phi)]; Solarctan55arctan550E15 > Sold:=map(u->evalf(convert(u,degrees)),Sol); Sold24.09484255degrees24.09484255degrees0degreesE16 > Phi:=Sol[1]; Φarctan55.E17 > Epsilon:=simplify(subs(phi=Phi,epsilon),symbolic); Ε´4π25cos2πτ7+cos2πτ225+10cos2πτ2+cos2πτ4.E18 It is now possible again to plot curves describing the angular acceleration dependence on τand ϕ, and highlight the curve for ϕ = Φ, see Figure 7. > EO:=plot(subs(phi=Phi,epsilon),tau=0..1,color=blue,thickness=3): > display({EP,EO}); To make sure that the angle Φ is optimal, we will plot the angular velocity depending on τ, see Figure 8. > Omega:=simplify(subs(phi=Phi,T=1,omega),symbolic); Ω2π5sin2πτ5+cos2πτ2.E19 > OP:=plot([seq(omega,phi=expand(Pi/36*[$1..9]))],

tau=0..1,color=red):

> OO:=plot(Omega,tau=0..1,color=blue,thickness=3):

> display({OP,OO});

And yet it is possible to plot the dependence of the angle θon τ, see Figure 9.

> Theta:=simplify(subs(phi=Phi,theta));

Θarctan55cos2πτ.E20

> TP:=plot([seq(theta,phi=expand(Pi/36*[$1..9]))], tau=0..1,color=red): > TO:=plot(Theta,tau=0..1,color=blue,thickness=3): > display({TP,TO}); Advertisement ## 5. Plotting of the dependencies in phase spaces Further evidence of the choice of the optimum angle value Φ, that is, the choice of a suitable compromise between the oscillation amplitude and smooth oscillation over time, can be given by plotting the phase diagrams, which are graphs in which one kinematic variable is plotted against another kinematic variable. In this case, we will present three phase diagrams, where we will plot the dependence of the oscillation angular velocity on the current value of the oscillation angle—ω(θ(τ)), see Figure 10, then the course of the angular acceleration of the oscillations against the current value of the oscillation angle—ε(θ(τ)), see Figure 11 and finally the course of the angular acceleration of the oscillations versus the angular velocity of the oscillations—ε(ω(τ)), see Figure 12. > TW:=plot([seq([theta,omega,tau=0..1], phi=expand(Pi/36*[$1..9]))],color=red,numpoints=500):

> TWO:=plot([Theta,Omega,tau=0..1],color=blue,thickness=3):

> display({TW,TWO});

> TE:=plot([seq([theta,epsilon,tau=0..1],

phi=expand(Pi/36*[$1..9]))],color=red,numpoints=500): > TEO:=plot([Theta,Epsilon,tau=0..1],color=blue,thickness=3): > display({TE,TEO}); > OE:=plot([seq([omega,epsilon,tau=0..1], phi=expand(Pi/36*[$1..9]))],color=red,numpoints=500):

> OEO:=plot([Omega,Epsilon,tau=0..1],color=blue,thickness=3):

> display({OE,OEO});

## 6. Results and conclusions

All of the presented calculations and graphs show that the optimal angle of the nose inclination is ϕ = arctan(5−1/2) = 0.4205343352 rad = 24° 5′ 41.43″. If the angle ϕis less than the optimum value, then the amplitude of the resulting oscillations will be low and gears will be required. This will lead to a significant increase in the price of the whole pin. Otherwise, when the value of the angle ϕis higher than the optimum value, there is a sharp increase in the instantaneous angular acceleration of the oscillations. High angular acceleration values mean that the entire pin must transmit high torque values to which the entire pin must be dimensioned. This will mean either increasing the size of the pin or using materials with satisfactory mechanical properties. This solution again leads to an increase in the price of the whole pin. If the pin is not appropriately sized for the transmitted torque, it may be mechanically damaged or blocked. The condition for the optimum angle value ϕin Eq. (12), thus, appears to be the optimal compromise between the oscillation amplitude and the mechanical stress of the entire joint.

Of course, the optimum shape of the oscillations should have a sinusoidal waveform, i.e., θ(t) = a sin(2πt/T), where θ = oscillation angle, t = time, a = oscillation amplitude, and T = oscillation period. The corresponding dependencies of angular velocity and angular acceleration over time are also in the form of sine or cosine functions. If we draw phase diagrams for these oscillations, then the graph of the dependencies [θ(τ), ω(τ)] = ω(θ(τ)) and [ω(τ), ε(τ)] = ε(ω(τ)) is the shape of an ellipse and the dependency [θ(τ), ε(τ)] = ε(θ(τ)) is the shape of a line segment. The amplitude of the oscillations is then equal to the length of the main half-axis of the ellipse, which displays ε(θ(τ)).

If we examine the graphs in Figures 1012, we find that the highlighted curve that corresponds to the optimum angle ϕare an excellent compromise between the shape of the curve (ellipse or line segment) and the dimensions of the curve that correspond to the oscillation amplitude. Similarly, the graphs of functions θ(τ), ω(τ), and ε(τ), which correspond to the optimal value of the angle ϕ, are the best compromise between the desired function shape and the oscillation amplitude, see the graphs in Figures 79.

Figure 1 shows the entire inclined swinging pin, while Figure 2 shows the tenon plotted for this angle. Several bottom views of the entire joint during one quarter of the working period are shown in Figure 13.

The main advantages of an inclined swinging pin compared to similar mechanisms are:

1. Simplicity: The whole pin consists of only three main components that can be forged.

2. Compactness: All components are similar in dimensions.

3. Durability: All parts of the mechanism are subjected to approximately the same forces. No component will receive more wear than others.

For these reasons, it is used in a wide variety of agricultural machines.

An animation (An_Joint.avi) that displays the entire device's activity can be downloaded from the Videos section.

## 7. Videos

All videos referenced in this chapter are available from: https://bit.ly/2vZVGXR

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Stanislav Barton (June 19th 2019). Optimization of Kinematics of Inclined Swinging Pin, Kinematics - Analysis and Applications, Joseph Mizrahi, IntechOpen, DOI: 10.5772/intechopen.86589. Available from:

Next chapter

By Jiří Ondrášek