Analytical Method for Determination of Internal Forces of Mechanisms and Manipulators

This paper presents a theory for the analytical determination of internal forces in the links of planar linkage mechanisms and manipulators with statically determinate structures, considering the distributed dynamic loads. Linkage mechanisms and manipulators were divided into elements and joints. Discrete models were created for both the elements and the entire mechanism. The dynamic equations of equilibrium for the discrete model of the elements and the hinged and rigid joints, under the action of longitudinal and transverse distributed dynamic trapezoidal loads, were derived. In the dynamic equations of the equilibrium of the discrete model of the elements and joints, the connections between the components of the force vector in the calculated cross-sections and the geometric, physical, and kinematic characteristics of the element were established for its plane-parallel motion. According to the developed technique, programs were created in the Maple system, and animations of the motion of the mechanisms were produced. The links were constructed with the intensity of transverseand longitudinal-distributed dynamic loads, bending moments, and shearing and normal forces, depending on the kinematic characteristics of the links.


Introduction
One of the important problems in designing mechanisms and manipulators is ensuring the strength and stiffness of their links during full-time process.For strength and stiffness analysis, the laws of distribution of internal forces, which enable selecting the form of cross-sections and defining their linear sizes, should be considered.
To analyze the stress-strain state of linkage systems, there are different approaches: graphically analytic methods, forces methods, and displacement methods [1][2][3].When calculating using the forces method, the main sought-after values are the forces in redundant constraints.Knowledge of the forces in redundant constraints enables the use of sectioning to perform complete calculations to determine the forces that arise in cross-sections of elements in a given system.When calculating using the displacement method, the main sought-after values are the displacements of the nodal points caused by deformation of the system.Knowledge of these displacements is necessary and sufficient to determine all internal forces that arise in cross-sections of elements in the system.
As all linkage systems have distributed mass, they are always systems with a degree of freedom of equal infinity; however, in many cases, it is possible to minimize the calculation of such systems this method, computer programs are developed in the Maple system, and animations of movement of mechanism with the construction of links in terms of intensity of transverse and longitudinal inertial loads, bending moments, and shearing and normal forces, which depend on kinematic characteristics of links, were produced.

Distributed Dynamic Loads and Approximation Matrix
Let us consider the plane-parallel movement of the link k of the mechanism with constant cross-sections with respect to the fixed coordinate system OXY.The following laws of distribution of transverse and longitudinal inertial loads along the link that arise from the mass of the link are defined [21]: where: where θ k is the angle, defining position of link k with respect to the fixed coordinate system OXY; ω k , ε k are the angular velocity and angular acceleration of the link k, respectively; w x k kp and w y k kp are the components of acceleration of the point P k (pole) of the link k, are directed along and perpendicular to the axis of the link k, respectively; γ k is the specific weight of material of the link k; A k is the square of transverse cross-section of the link k; and g is the acceleration of gravity.
The found expressions show that the distribution of transverse and longitudinal inertial forces along the axis of the link with constant cross-sections is characterized by a trapezoidal law.
For the element k, which is under the influence of longitudinal distributed trapezoidal loads, as shown in Figure 1, the bending moments along the length of the element are distributed by a law of a third-order polynomial: Robotics 2018, 7, x FOR PEER REVIEW 3 of 14 inertial loads, bending moments, and shearing and normal forces, which depend on kinematic characteristics of links, were produced.

Distributed Dynamic Loads and Approximation Matrix
Let us consider the plane-parallel movement of the link k of the mechanism with constant cross-sections with respect to the fixed coordinate system OXY.The following laws of distribution of transverse and longitudinal inertial loads along the link that arise from the mass of the link are defined [21]: where: where k  is the angle, defining position of link k with respect to the fixed coordinate system OXY; The found expressions show that the distribution of transverse and longitudinal inertial forces along the axis of the link with constant cross-sections is characterized by a trapezoidal law.
For the element k, which is under the influence of longitudinal distributed trapezoidal loads, as shown in Figure 1, the bending moments along the length of the element are distributed by a law of a third-order polynomial: Now, express the bending moment in the cross-section ' k x through the required bending moments in the cross-sections demonstrated in Figure 1.For this purpose, it is enough to express the coefficients , respectively.As a result we have [22]: Now, express the bending moment in the cross-section x k through the required bending moments M k1 , M k2 , M k3 , M k4 in the cross-sections demonstrated in Figure 1.
For this purpose, it is enough to express the coefficients a 0 , a 1 , a 2 , a 3 through M k1 , M k2 , M k3 , M k4 , respectively.As a result we have [22]: Differentiating M k x k by x k results in the equation of the shearing force: Suppose, in addition to the transverse distributed load, the longitudinal distributed trapezoidal load acts on the element.In this case, the normal force in an arbitrary cross-section of the element can be expressed analogously to the previous expression by means of the normal forces in the calculated cross-sections as follows: Thus, for the element experiencing transverse and longitudinal distributed trapezoidal loads, the approximation matrix connecting the internal forces in any cross-section of the element, with the values of internal forces in the calculated cross-sections, has the form: The elements of the first row of this matrix can be seen in Equation (1); the elements of the second row can be seen in Equation (4); and the elements of the third row can be seen in Equation (5).
This expression of the approximation matrix of forces defines the relationship between the vector of forces S k x k in any cross-section of the element x k and the vector of forces in appointed cross-sections {S k }.For the element of linkage system, the approximation matrix is found exactly within the framework of the known laws of distribution of unknown forces.
We see that the equations of the bending moment, shearing force and normal forces, in Equations ( 3)- (5), respectively, which are expressed by the same values in the calculated cross-sections, demonstrate that, to define the internal loads of each element of the mechanism, it is enough to know the values of these loads in a finite number of cross-sections in each of these elements.The number of cross-sections in which it is necessary to know the values of internal forces are defined by the polynomial degrees of external actions.Thus, the internal forces of each continual link are determined unambiguously by a set of internal forces in its separate cross-sections and by the matrices of approximation.Therefore, the task is reduced to calculating internal forces in a finite number of cross-sections of the elements.Hence, we obtain a discrete model of elastic calculation of the links of linkage mechanisms and manipulators.

Discrete Models for Elastic Calculation of Elements and Entire Construction of Mechanisms and Manipulators
For elastic calculation of the linkage mechanisms based on D'Alambert's principle, all distributed dynamic and concentrated loads were attached to the links, as well as unknown driving moments (forces) were attached to the driving links, which ensured the assigned laws of their motion.Therefore, we obtained the construction with degrees of freedom equal to zero, if the rotational kinematic pair that connect the driving link and fixed base are replaced by a rigid restraint.
To define the internal forces in the links (in the elements), the construction of the mechanisms and manipulators was divided into elements and joints.The link or the part of the link could be the element, whereas the joints were the kinematic pairs, connecting the adjacent links and cross-sections where the concentrated external forces were attached.
The process of division of the construction included assignment of the calculated cross-sections of elements and their designations.When dividing the elements of the computational scheme of the construction into calculated cross-sections and joints, it was necessary to establish which internal connections between the elements should remain and which ones should be deleted.After discarding some internal connections or their combinations in the element, the element was then broken up into two elements that can rotate, move, or be removed from each other.To prevent this, the internal forces were attached in place of the discarded connections.Henceforth, these forces were considered the main unknowns [23].
Let us decompose the element of the planar linkage mechanism into three types of beams for the convenience of producing the decisive equations and determining the internal forces in the assigned cross-sections of the elements of mechanism.The first type of the element was a beam with two rigidly fixed ends, as demonstrated in Figure 2.Such beams can be rods with basic links when these rods are interconnected rigidly, as shown in Figure 3.

Discrete Models for Elastic Calculation of Elements and Entire Construction of Mechanisms and Manipulators
For elastic calculation of the linkage mechanisms based on D'Alambert's principle, all distributed dynamic and concentrated loads were attached to the links, as well as unknown driving moments (forces) were attached to the driving links, which ensured the assigned laws of their motion.Therefore, we obtained the construction with degrees of freedom equal to zero, if the rotational kinematic pair that connect the driving link and fixed base are replaced by a rigid restraint.
To define the internal forces in the links (in the elements), the construction of the mechanisms and manipulators was divided into elements and joints.The link or the part of the link could be the element, whereas the joints were the kinematic pairs, connecting the adjacent links and cross-sections where the concentrated external forces were attached.
The process of division of the construction included assignment of the calculated cross-sections of elements and their designations.When dividing the elements of the computational scheme of the construction into calculated cross-sections and joints, it was necessary to establish which internal connections between the elements should remain and which ones should be deleted.After discarding some internal connections or their combinations in the element, the element was then broken up into two elements that can rotate, move, or be removed from each other.To prevent this, the internal forces were attached in place of the discarded connections.Henceforth, these forces were considered the main unknowns [23].
Let us decompose the element of the planar linkage mechanism into three types of beams for the convenience of producing the decisive equations and determining the internal forces in the assigned cross-sections of the elements of mechanism.The first type of the element was a beam with two rigidly fixed ends, as demonstrated in Figure 2.Such beams can be rods with basic links when these rods are interconnected rigidly, as shown in Figure 3.To determine the coefficients of expressions of the bending moment in Equation ( 3), it was necessary to know the values of the bending moments in the four cross-sections.To determine the coefficients of expressions of the normal forces in Equation ( 5), it was necessary to know the values of the normal forces in the three cross-sections of the element.Therefore, in this beam, we chose four cross-sections with unknown bending moments and three cross-sections with unknown normal forces.Then, considering the conditional schemes with corresponding unknowns, we constructed the discrete model of the considered beam, shown in Figure 4.

Discrete Models for Elastic Calculation of Elements and Entire Construction of Mechanisms and Manipulators
For elastic calculation of the linkage mechanisms based on D'Alambert's principle, all distributed dynamic and concentrated loads were attached to the links, as well as unknown driving moments (forces) were attached to the driving links, which ensured the assigned laws of their motion.Therefore, we obtained the construction with degrees of freedom equal to zero, if the rotational kinematic pair that connect the driving link and fixed base are replaced by a rigid restraint.
To define the internal forces in the links (in the elements), the construction of the mechanisms and manipulators was divided into elements and joints.The link or the part of the link could be the element, whereas the joints were the kinematic pairs, connecting the adjacent links and cross-sections where the concentrated external forces were attached.
The process of division of the construction included assignment of the calculated cross-sections of elements and their designations.When dividing the elements of the computational scheme of the construction into calculated cross-sections and joints, it was necessary to establish which internal connections between the elements should remain and which ones should be deleted.After discarding some internal connections or their combinations in the element, the element was then broken up into two elements that can rotate, move, or be removed from each other.To prevent this, the internal forces were attached in place of the discarded connections.Henceforth, these forces were considered the main unknowns [23].
Let us decompose the element of the planar linkage mechanism into three types of beams for the convenience of producing the decisive equations and determining the internal forces in the assigned cross-sections of the elements of mechanism.The first type of the element was a beam with two rigidly fixed ends, as demonstrated in Figure 2.Such beams can be rods with basic links when these rods are interconnected rigidly, as shown in Figure 3.To determine the coefficients of expressions of the bending moment in Equation (3), it was necessary to know the values of the bending moments in the four cross-sections.To determine the coefficients of expressions of the normal forces in Equation ( 5), it was necessary to know the values of the normal forces in the three cross-sections of the element.Therefore, in this beam, we chose four cross-sections with unknown bending moments and three cross-sections with unknown normal forces.Then, considering the conditional schemes with corresponding unknowns, we constructed the discrete model of the considered beam, shown in Figure 4.  To determine the coefficients of expressions of the bending moment in Equation (3), it was necessary to know the values of the bending moments in the four cross-sections.To determine the coefficients of expressions of the normal forces in Equation ( 5), it was necessary to know the values of the normal forces in the three cross-sections of the element.Therefore, in this beam, we chose four cross-sections with unknown bending moments and three cross-sections with unknown normal forces.Then, considering the conditional schemes with corresponding unknowns, we constructed the discrete model of the considered beam, shown in Figure 4.

Discrete Models for Elastic Calculation of Elements and Entire Construction of Mechanisms and Manipulators
For elastic calculation of the linkage mechanisms based on D'Alambert's principle, all distributed dynamic and concentrated loads were attached to the links, as well as unknown driving moments (forces) were attached to the driving links, which ensured the assigned laws of their motion.Therefore, we obtained the construction with degrees of freedom equal to zero, if the rotational kinematic pair that connect the driving link and fixed base are replaced by a rigid restraint.
To define the internal forces in the links (in the elements), the construction of the mechanisms and manipulators was divided into elements and joints.The link or the part of the link could be the element, whereas the joints were the kinematic pairs, connecting the adjacent links and cross-sections where the concentrated external forces were attached.
The process of division of the construction included assignment of the calculated cross-sections of elements and their designations.When dividing the elements of the computational scheme of the construction into calculated cross-sections and joints, it was necessary to establish which internal connections between the elements should remain and which ones should be deleted.After discarding some internal connections or their combinations in the element, the element was then broken up into two elements that can rotate, move, or be removed from each other.To prevent this, the internal forces were attached in place of the discarded connections.Henceforth, these forces were considered the main unknowns [23].
Let us decompose the element of the planar linkage mechanism into three types of beams for the convenience of producing the decisive equations and determining the internal forces in the assigned cross-sections of the elements of mechanism.The first type of the element was a beam with two rigidly fixed ends, as demonstrated in Figure 2.Such beams can be rods with basic links when these rods are interconnected rigidly, as shown in Figure 3.To determine the coefficients of expressions of the bending moment in Equation (3), it was necessary to know the values of the bending moments in the four cross-sections.To determine the coefficients of expressions of the normal forces in Equation ( 5), it was necessary to know the values of the normal forces in the three cross-sections of the element.Therefore, in this beam, we chose four cross-sections with unknown bending moments and three cross-sections with unknown normal forces.Then, considering the conditional schemes with corresponding unknowns, we constructed the discrete model of the considered beam, shown in Figure 4.Then, the force vector in the calculated cross-sections of the discrete model for this beam is expressed by the following vector: There is a dependence between the degree of freedom of discrete model m, the number of attached external forces n, and the degree of static indeterminacy k of the computational scheme [1]: This equation conveniently simplifies determining the degrees of freedom of the discrete model.The total number of forces n in the calculated cross-sections are easily counted, and the degree of static indeterminacy of the computational scheme is found using the formula k = 3K − I I I, where K is the number of closed loops, is the number of single hinges, and k is the degree of static indeterminacy of computational scheme of mechanism.
The degrees of freedom m of the discrete model determines the number of necessary independent equations of statics.Let us define the degrees of freedom of discrete model of this beam.The number of unknowns is n = 7, the static indeterminacy of the beam is k = 3, so the degrees of freedom of the discrete model m = 4.In other words, it is possible to derive four independent equilibrium equations for this discrete beam model.
The second type of element is the beam, where one end is fixed rigidly and the other is fixed by the motionless hinge.Such beams can be the driving links of planar linkage mechanisms.
The third type of element is the beam of the intermediate links.They can be considered as beams fixed with motionless hinged supports at both ends.The discrete models for the second and third types of beams can be constructed similarly to the discrete model for the first type of beam.
The discrete model of the four-bar linkage is shown in Figure 5, where all the unknown values that define all internal forces in any cross-section of links of the mechanism are illustrated.For the first link (the second type of beam) of this mechanism, the vector of forces in the calculated cross-sections {S 1 } has the following components: Robotics 2018, 7, x FOR PEER REVIEW 6 of 14 Then, the force vector in the calculated cross-sections of the discrete model for this beam is expressed by the following vector: There is a dependence between the degree of freedom of discrete model m , the number of attached external forces n , and the degree of static indeterminacy k of the computational scheme [1]: This equation conveniently simplifies determining the degrees of freedom of the discrete model.The total number of forces n in the calculated cross-sections are easily counted, and the degree of static indeterminacy of the computational scheme is found using the formula , where K is the number of closed loops, Ш is the number of single hinges, and In other words, it is possible to derive four independent equilibrium equations for this discrete beam model.
The second type of element is the beam, where one end is fixed rigidly and the other is fixed by the motionless hinge.Such beams can be the driving links of planar linkage mechanisms.
The third type of element is the beam of the intermediate links.They can be considered as beams fixed with motionless hinged supports at both ends.The discrete models for the second and third types of beams can be constructed similarly to the discrete model for the first type of beam.
The discrete model of the four-bar linkage is shown in Figure 5, where all the unknown values that define all internal forces in any cross-section of links of the mechanism are illustrated.For the first link (the second type of beam) of this mechanism, the vector of forces in the calculated cross-sections    For the second and third links (the third type of beam) of the considered mechanism, the vector of the forces in the calculated cross-sections have the following components, respectively: For the second and third links (the third type of beam) of the considered mechanism, the vector of the forces in the calculated cross-sections have the following components, respectively: For the entire discrete model of the mechanism, the vector of the forces in the calculated cross-sections are:

Dynamic Equilibrium Equations of the Discrete Models of the Elements and Joints
Let us derive the equations for the dynamic equilibrium of the element.From the applied concentrated external loads (Q k1 , M k1 ) and the transverse distributed trapezoidal loads along the axis of the element, the bending moment arises in the arbitrary cross-section x k and is defined by Equation ( 2).The bending moment in cross-section x k of the element that is expressed through the sought-after moments in the calculated cross-sections is determined using Equation (3).
Differentiating Equations ( 2) and (3) three times by x k , equating them, and substituting the value of b kq , the first equation of dynamic equilibrium becomes: The dependence between the values of the sought-after magnitudes of the bending moments in the calculated cross-sections, and the geometric, physical, and kinematic characteristics of the element k of the mechanism are found.The second equation of dynamic equilibrium is produced by taking the sum of the moments of all acting forces on the element k relative to the center of gravity of the cross-section k4 (Figure 1).Then, the following expression is derived: where: This expression is not difficult to derive if we substitute the value x k = 0 into Equation (4).Substituting the values Q k1 , a kq and b kq into Equation ( 13), and summing the coefficients of the same known and unknown magnitudes of the equation into the right-hand side, the second equation of the dynamic equilibrium is produced: From the longitudinal distributed trapezoidal loads acting on the element, and from the force N k1 of the cross-section k1, normal force occurs in the cross-section x k of the element, which is defined by: The normal force in the cross-section x k of the element, expressed through the normal forces in the calculated cross-sections, is calculated using Equation (5).
Differentiating the Equations ( 5) and ( 15) two times by x k , equating them, and substituting the value of b kn , the third dynamic equilibrium equation is produced: Projecting all the forces acting on element k onto axle x k and substituting the values of a kn , b kn , the fourth equation of dynamic equilibrium is found: The resulting system of equations consisting of Equations ( 12), ( 14), ( 16) and (17), can be written down in matrix form: where: Let the two elements j and k of the mechanism form a rotational kinematic pair, i.e., tolerate rotational motion relative to each other.Let the length of these elements have constant cross-sections.We cut out the kinematic pair with adjacent cross-sections of the elements forming this pair from the mechanism.In this case, in the cross-sections of the elements adjacent to the joint (to the kinematic pair), the internal forces occur as shown in Figure 6.For such joints, we have two equilibrium conditions.These equilibrium equations for the joint under consideration are written as: Robotics 2018, 7, x FOR PEER REVIEW 9 of 14  The magnitudes of Q k1 , Q j4 are expressed by means of the sought-after moments in the calculated cross-sections of the discrete model of the elements.To this end, we used Equation (4) for the shearing force, and by substituting the values x k = 0 and x j = l j , we obtain, respectively: Now, substituting the values of Q k1 and Q j4 into Equation ( 19), we have the following equilibrium equations for the Rigid joints can also be cross-sections of the link where external concentrated forces are applied.The cross-sections of the links can be rigid joints if external concentrated loads are attached in this cross-section.For instance, let the concentrated loads P kx k and P ky k , and the concentrated moment M k , be attached to the cross-section G of the link k (Figure 7).Then, the link k is divided into two elements: kth and ith. Figure 7 shows the joint with adjacent cross-sections, where the arising internal forces are shown.For this joint, the following three conditions of dynamic equilibrium are expressed through the sought-after parameters of the elements: Robotics 2018, 7, x FOR PEER REVIEW 10 of 14

Decisive Equations for Determining Internal Forces
By combining the equations of the equilibrium of the elements and joints into one system, we obtain the equilibrium equations for the entire discrete model of the mechanism.They can be written in the general form: These systems of equations are sufficient for determining internal forces in links of mechanisms with statically determinate structures.The matrix of equilibrium equations for the discrete model of the mechanisms consists of the matrices of the equilibrium equations of their individual elements, as well as the equilibrium equation of their joints.The matrix of the equilibrium equations for the discrete models of mechanisms is as follows:

Decisive Equations for Determining Internal Forces
By combining the equations of the equilibrium of the elements and joints into one system, we obtain the equilibrium equations for the entire discrete model of the mechanism.They can be written in the general form: [A]{S} = {F}.
These systems of equations are sufficient for determining internal forces in links of mechanisms with statically determinate structures.The matrix of equilibrium equations for the discrete model of the mechanisms consists of the matrices of the equilibrium equations of their individual elements, as well as the equilibrium equation of their joints.The matrix of the equilibrium equations for the discrete models of mechanisms is as follows: The vector of load and the vector of the force in the calculated cross-sections of the discrete models of the mechanisms are formed from the vectors of the load and the forces in the calculated cross-sections of their individual elements.These vectors in the vector form have the following form: Determination of internal forces will be outlined using an example of a second-class six-bar linkage with one driving link, shown in Figure 8. Computer programs were developed in the Maple system to determine and construct the diagrams of inertial and internal forces on the links.The resulting inertial and internal forces are shown in Figures 8-13 for some positions of the mechanism.
in the general form: These systems of equations are sufficient for determining internal forces in links of mechanisms with statically determinate structures.The matrix of equilibrium equations for the discrete model of the mechanisms consists of the matrices of the equilibrium equations of their individual elements, as well as the equilibrium equation of their joints.The matrix of the equilibrium equations for the discrete models of mechanisms is as follows: The vector of load and the vector of the force in the calculated cross-sections of the discrete models of the mechanisms are formed from the vectors of the load and the forces in the calculated cross-sections of their individual elements.These vectors in the vector form have the following form: ., .,. , ; ., . ., ,

Results and Discussion
To verify our developed theory, we applied the theory to solving a task involving a specific six-bar linkage.We solved the kinematics problem and determined the internal forces in the links, as well as animated the motion and constructed the diagrams of the distributed dynamic loads and internal forces on the links.Thus, we wanted to show that our developed theory works and determine the dynamic loads and the internal forces in the links, depending on how the mechanism position changes.To determine the maximum values of the internal forces, it was necessary to determine internal forces in all positions of the mechanism and manipulator.The maximum values of internal forces allow, according to the appropriate strength theories, to select the shape and find linear dimensions of the cross-sections of the links.The validity of the results can be seen from the plotted diagrams, for example, in Figure 11 in the cross-section A of the link 1.The driving moment is shown and the moment in the cross-section B is zero, since there is a rotational kinematic pair.Since there are differential dependencies between the distributed transverse dynamic loads, the shearing force, and the bending moment, using these relationships all the diagrams can be checked.For example, in Figure 12, in the cross-section of link 2 where the shearing force is zero, the bending moment in the same cross-section in Figure 11 has a maximum value and so on.

Results and Discussion
To verify our developed theory, we applied the theory to solving a task involving a specific six-bar linkage.We solved the kinematics problem and determined the internal forces in the links, as well as animated the motion and constructed the diagrams of the distributed dynamic loads and internal forces on the links.Thus, we wanted to show that our developed theory works and determine the dynamic loads and the internal forces in the links, depending on how the mechanism position changes.To determine the maximum values of the internal forces, it was necessary to determine internal forces in all positions of the mechanism and manipulator.The maximum values of internal forces allow, according to the appropriate strength theories, to select the shape and find linear dimensions of the cross-sections of the links.The validity of the results can be seen from the plotted diagrams, for example, in Figure 11 in the cross-section A of the link 1.The driving moment is shown and the moment in the cross-section B is zero, since there is a rotational kinematic pair.Since there are differential dependencies between the distributed transverse dynamic loads, the shearing force, and the bending moment, using these relationships all the diagrams can be checked.For example, in Figure 12, in the cross-section of link 2 where the shearing force is zero, the bending moment in the same cross-section in Figure 11 has a maximum value and so on.

Conclusions
In this study, we established the laws of distribution of the distributed loads from inertial forces and forces of gravity, arising from the distributed weight of the links with constant cross-sections.Dependencies between the distributed loads and with geometrical, physical, and kinematic characteristics of the links were determined.The approximation matrix [H k (x)] of the internal forces of the element under the action of distributed loads with trapezoidal shape intensity was found.The approximation matrices of the internal forces define the relationship between the vector of the force S k x k in any cross-section of the element x k and the vector of forces in the calculated cross-sections {S k (x)} = [H k (x)]{S k }.The computational and discrete schemes of the linkage mechanisms for elastic calculation were developed.The dynamic equations of equilibrium for the discrete model of each element, as well as the dynamic equations of equilibrium for hinged and rigid joints, under the action of transverse and longitudinal distributed trapezoidal loads, were derived.In the dynamic equilibrium equations of the discrete model of the elements, the connections were established between the components of the vector of the forces in the calculated cross-sections and with geometric, physical, and kinematic characteristics of links with constant cross-sections, in their plane-parallel movement.Decisive equations were derived for determining internal forces in the links of the mechanisms with a statically determinate structure.In using the developed technique, programs were created in the Maple system and animations of the motion of the mechanisms were produced.The links were constructed with the intensity of the distributed transverse and longitudinal dynamic loads, the bending moments, and the shearing and normal forces, depending on the kinematic characteristics of the links.
, are the angular velocity and angular acceleration of the link k, respectively; of acceleration of the point k P (pole) of the link k, are directed along and perpendicular to the axis of the link k, respectively; k  is the specific weight of material of the link k; k A is the square of transverse cross-section of the link k; and g is the acceleration of gravity.

Figure 1 .
Figure 1.The element under the action of transverse distributed trapezoidal loads.

Figure 1 .
Figure 1.The element under the action of transverse distributed trapezoidal loads.

Figure 2 .
Figure 2. The beam in which both ends are connected rigidly (the first type of beam).

Figure 3 .
Figure 3.The basic link in which rods are interconnected rigidly.

Figure 4 .
Figure 4.The discrete model of the first type of beam under the action of distributed trapezoidal loads.

Figure 2 .
Figure 2. The beam in which both ends are connected rigidly (the first type of beam).

Figure 2 .
Figure 2. The beam in which both ends are connected rigidly (the first type of beam).

Figure 3 .
Figure 3.The basic link in which rods are interconnected rigidly.

Figure 4 .
Figure 4.The discrete model of the first type of beam under the action of distributed trapezoidal loads.

Figure 3 .
Figure 3.The basic link in which rods are interconnected rigidly.

Figure 2 .
Figure 2. The beam in which both ends are connected rigidly (the first type of beam).

Figure 3 .
Figure 3.The basic link in which rods are interconnected rigidly.

Figure 4 .
Figure 4.The discrete model of the first type of beam under the action of distributed trapezoidal loads.Figure 4. The discrete model of the first type of beam under the action of distributed trapezoidal loads.

Figure 4 .
Figure 4.The discrete model of the first type of beam under the action of distributed trapezoidal loads.Figure 4. The discrete model of the first type of beam under the action of distributed trapezoidal loads.

k
is the degree of static indeterminacy of computational scheme of mechanism.The degrees of freedom m of the discrete model determines the number of necessary independent equations of statics.Let us define the degrees of freedom of discrete model of this beam.The number of unknowns is

1 S
has the following components:

Figure 5 .
Figure 5.The discrete model of the four-bar mechanism with constant cross-sections.

Figure 5 .
Figure 5.The discrete model of the four-bar mechanism with constant cross-sections.

Figure 6 .
Figure 6.Hinge joint of the mechanism with constant cross-sections of the element.

Figure 6 .
Figure 6.Hinge joint of the mechanism with constant cross-sections of the element.

Figure 7 .
Figure 7. Rigid joint with constant cross-sections of the element, where the external concentrated loads are attached.

Figure 7 .
Figure 7. Rigid joint with constant cross-sections of the element, where the external concentrated loads are attached.


Determination of internal forces will be outlined using an example of a second-class six-bar linkage with one driving link, shown in Figure8.Computer programs were developed in the Maple system to determine and construct the diagrams of inertial and internal forces on the links.The resulting inertial and internal forces are shown in Figures8-13for some positions of the mechanism.

Figure 8 .
Figure 8.The second-class six-bar mechanism with one driving link.Figure 8.The second-class six-bar mechanism with one driving link.

Figure 8 .
Figure 8.The second-class six-bar mechanism with one driving link.Figure 8.The second-class six-bar mechanism with one driving link.Robotics 2018, 7, x FOR PEER REVIEW 11 of 14

Figure 9 .
Figure 9.The mechanism under discussion of the links showing the constructed transverse dynamic distributed loads.

Figure 9 .
Figure 9.The mechanism under discussion of the links showing the constructed transverse dynamic loads.

Figure 9 .
Figure 9.The mechanism under discussion of the links showing the constructed transverse dynamic distributed loads.

Figure 10 .
Figure 10.The mechanism of the links illustrating the constructed longitudinal dynamic distributed loads.

Figure 11 .
Figure 11.The mechanism of the links showing the constructed bending moments.

Figure 10 .
Figure 10.The mechanism of the links illustrating the constructed longitudinal dynamic distributed loads.

Figure 9 .
Figure 9.The mechanism under discussion of the links showing the constructed transverse dynamic distributed loads.

Figure 10 .
Figure 10.The mechanism of the links illustrating the constructed longitudinal dynamic distributed loads.

Figure 11 .
Figure 11.The mechanism of the links showing the constructed bending moments.Figure 11.The mechanism of the links showing the constructed bending moments.

Figure 11 . 14 Figure 12 .
Figure 11.The mechanism of the links showing the constructed bending moments.Figure 11.The mechanism of the links showing the constructed bending moments.Robotics 2018, 7, x FOR PEER REVIEW 12 of 14

Figure 12 .
Figure 12.The mechanism of the links showing the constructed shearing force.Figure 12.The mechanism of the links showing the constructed shearing force.

Figure 12 .
Figure 12.The mechanism of the links showing the constructed shearing force.

Figure 13 .
Figure 13.The mechanism of the links showing the constructed normal force.

Figure 13 .
Figure 13.The mechanism of the links showing the constructed normal force.