Next Article in Journal
DEA-BSC and Diamond Performance to Support Museum Management
Previous Article in Journal
Frequency and Severity Dependence in the Collective Risk Model: An Approach Based on Sarmanov Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Analytical Model Used in Finite Element Analysis of Solids Mechanics

by
Sorin Vlase
1,2,*,
Adrian Eracle Nicolescu
3,* and
Marin Marin
4,*
1
Department of Mechanical Engineering, Transilvania University of Brașov, 500036 Brașov, Romania
2
Romanian Academy of Technical Sciences, B-dul Dacia 26, 030167 Bucharest, Romania
3
Faculty of Mechanical, Industrial and Maritime Engineering, Ovidius University of Constanta, 900527 Constanța, Romania
4
Department of Mathematics and Computer Science, Transilvania University of Brasov, 500036 Brașov, Romania
*
Authors to whom correspondence should be addressed.
Mathematics 2020, 8(9), 1401; https://doi.org/10.3390/math8091401
Submission received: 13 July 2020 / Revised: 15 August 2020 / Accepted: 18 August 2020 / Published: 21 August 2020
(This article belongs to the Section Mathematical Physics)

Abstract

:
In classical mechanics, determining the governing equations of motion using finite element analysis (FEA) of an elastic multibody system (MBS) leads to a system of second order differential equations. To integrate this, it must be transformed into a system of first-order equations. However, this can also be achieved directly and naturally if Hamilton’s equations are used. The paper presents this useful alternative formalism used in conjunction with the finite element method for MBSs. The motion equations in the very general case of a three-dimensional motion of an elastic solid are obtained. To illustrate the method, two examples are presented. A comparison between the integration times in the two cases presents another possible advantage of applying this method.

1. Introduction

The analytical methods of mechanics, clearly stated in Lagrange’s book (1788) [1], were alternative formulations to Newton-Euler’s equations of motion (EoM). Although initially the interest was more theoretical, these methods proved their usefulness when the complexity of the studied systems began to increase. Lagrangian mechanics experienced a continuous development in the 19th century and various formulations appeared in this context, such as Gibbs-Appell’s (1879, 1899) methods [2,3], Hamilton’s equations (1835) [4], Maggi’s equations (1896) [5] and, in the middle of the last century, Kane’s equations (1924) [6] (essentially a new form of writing Maggi’s equations [7,8]). The interest in these possibilities of expressing the evolution equations of a mechanical system, starting from the fundamental principles was, for a long time, a theoretical interest, without significant practical consequences.
A flexible multibody system (FMS) represents a group of rigid and deformable elements, each of which can undergo large translational and rotational movements. The components may come into contact with each other or with the environment. The classic kinematic linkages implied in engineering applications are: cylindrical and spherical joints, prismatic couplings, flat connections, screws, gears or cams. These elements can be connected in a closed loop (closed kinematic chain) as in classical mechanism theory or open (open kinematic chain) as in the case of manipulators. The flexible multibody dynamics (FMD) of an MBS impose computational strategies that are used to determine the dynamic response that involves knowing the time history of the mechanical behavior of the system, deformations, stresses due to external forces, liaisons and contour conditions.
In recent decades, however, the continuous industrialization that led to the emergence of machines and equipment operating at high speeds and accelerations, acting with intense forces, required the continuous improvement of the models and computational methods used. As a result, a reconsideration of the analytical methods previously used for the study of multibody systems followed and results were obtained which demonstrated that, in special cases, the use of one of the above methods may be more advantageous, both in terms of modeling and from the point of view of the calculation volume and, consequently, of the necessary computational time [7,8,9].
It is obvious that the need to solve practical problems has led to the development of the field and harmonized strategies. The excellent and well-documented work [9] makes a classification and ordering of the methods used. The modeling of flexible elements, within an MBS system, is one of the important steps for such a study.
In the analysis of an MBS, an inertial reference system is used as a global reference system to describe the motion of the entire system. Within the system, intermediate reference systems are linkage attached to each flexible element (component). This flexible elastic element has an average general motion in space. In this case, the movement of the flexible element in relation to the local reference system is mainly due to the elastic deformation. Researchers mainly use two types of local reference systems: floating and corotational frames. The first type, the floating frame, considers an average rigid body motion and follows this hypothetical rigid body. The second type, the corotational frame, considers an individual finite element and follows its average rigid body motion. However, it is possible to not use intermediate reference systems and the deformations reported directly to the global coordinate system. Our presentation uses the second type of reference frame.
We make a brief review of the main analytical methods used in FMD.
The use of Lagrange’s equations allows for an algorithmization of dynamics problems, the steps necessary to be completed being [10,11]:
  • choosing generalized coordinates (which do not need to be Cartesian coordinates);
  • calculating kinetic energy;
  • establishing generalized forces.
After that, some pre-established differentiation calculations are performed, finally reaching the equations of motion of the studied mechanical system, which are second-order differential equations that allow for the determination of the unknown generalized coordinates. As a conclusion, the basic notions studied using analytical mechanics are the generalized coordinates (or velocities), the kinetic energy, the potential energy and the work.
Known since the early twentieth century, the equations called Gibbs-Appell (GA) equations [10] have rarely been used by researchers to solve practical engineering problems related to the evolution of mechanical systems. They have advantages in terms of the number of the equations needed for modeling and computational effort. The GA method was established as a convenient procedure for the study of non-holonomic systems. To apply this method, a new notion is introduced, with which researchers are less familiar, namely “acceleration energy”, which will replace the Lagrangian one. The advantages of applying the method have been revealed by the latest studies on the dynamic analysis of mechanical systems that operate at high speeds and accelerations and operate with high forces. MBSs are very suitable for using the method [12]. The equations finally obtained coincide with those obtained by applying Lagrange’s equations, but there is the advantage that the number of operations to be performed is smaller. The use of the method shows the possibility of eliminating Lagrange multipliers and, consequently, reducing the number of computation operations [13].
Due to these advantages, it may become a procedure and researchers are beginning to use it with predilection in the analysis of modern mechanical systems [14].
Maggi’s equations are relatively recent, and represent an alternative form to other formalisms in analytical mechanics [10,11]. The feedback linearization technique in the case of control systems with which robots and manipulators are equipped naturally leads to the application of these types of equations. They are a simple way to obtain the equations of motion of the system in the case of non-holonomic liaisons. For large systems, defined by a large number of degrees of freedom (DOF), the application of the method allows for the shortening of the computational time required for the analysis. However, the method is not familiar to researchers and, as a result, is less used. Classical methods of obtaining EoM, such as Newton-Euler (NE) or Lagrange (LE), involve the liaisons forces. If the systems are large, this means a high number of calculus and, as a consequence, a high modeling and computation effort [15,16,17].
The method of Kane’s equations, which is relatively new, has begun to be successfully applied in recent decades in the modeling and calculation of MBSs [18]. However, if the system has a large number of bodies, with many constraints, the use of this method becomes problematic in determining the dynamic response of such a system [19,20]. In essence, the Kane equations method eliminates some of the disadvantages of classical LE and NE methods and can be used for both holonomic and non-holonomic constraint systems. Kane’s equations are an alternative form of Maggi’s equations [21]. Some interesting applications of Kane’s equations can be found in [22].
The analytical classic formalism in the study of MBSs was obtained using Lagrange’s equations as starting point. The Lagrangian methods are an efficient way to carry out the dynamics analysis of a flexible mechanical system. The algorithms based on these equations introduce, through a set of multipliers, the liaison forces in the dynamic equations—usually in a linear form. However, there are also formalisms that have been developed, starting from Hamilton’s equations. These are written in terms of generalized displacements and momenta. Hamiltonian mechanics achieve this precisely, providing a result of the modeling a system of differential equations of the first order with 2n unknowns. The first n unknowns are the generalized coordinates, the others being the generalized momenta [23,24]. The Hamilton method could have the advantage of providing us with a system of first-order equations, a system that can be used directly for numerical solving, without the need for the prior processing of obtained systems [25].

2. Preliminary

The first papers applying FEA to the study of an elastic MBS first considered a one-dimensional finite element with a two-dimensional motion [26,27,28]. The method used with predilection in all these works was the method of Lagrange’s equations. This method involves the calculation of well-known notions of mechanics (kinetic energy, potential and work). The complexity of the studied systems has increased over time, as well as the types of finite elements used. It went from the study of two-dimensional motion to the study of three-dimensional motion and increasingly complex two-dimensional and three-dimensional finite elements were used [29,30,31]. Different aspects of MBSs with elastic elements have been analyzed lately, but the method of obtaining the equations of motion was also the method of Lagrange’s equations [32,33,34,35].
Consider one single three-dimensional finite element referred to a local coordinate system, moving together with the element. The velocity of an arbitrary point of the element, M, after deformation, when it becomes M’, is [36]:
{ v M } G = { r ˙ M } G = { r ˙ O } G + [ R ˙ ] { r } L + [ R ˙ ] [ N ] { δ } L + [ R ] [ N ] { δ ˙ } L
where:
  • { v M } G = { r ˙ M } G is the velocity of the arbitrary chosen point;
  • { v O } G = { r ˙ O } G is the velocity of the origin of the mobile reference frame;
  • [ R ] is the rotation matrix (transform a vector from the mobile reference frame to the fixed reference frame);
  • [ N ] is the matrix of the shape functions depending on the type of finite element used;
  • { r } L is the position vector of a current point in the local reference frame;
  • { δ } L is the vector of the generalized coordinates;
  • { δ ˙ } L is the vector of the generalized velocity;
  • with the index L, a size is noted (vector, matrix), with the components expressed in the local reference frame, and with G, a size with the components expressed in the global reference frame.
The kinetic energy of the moving element is:
E c = 1 2 V ρ { v M } G T { v M } G d V
where ρ is denoted the density of the material. The potential energy has the well-known form:
E p = 1 2 V   { δ } L T [ k ] { δ } d V
The matrix [ k ] represents the rigidity of the element [18]. Concentrated forces { q } L give the work:
W c = { q } L T { δ } L
and the vector of volume forces { p } = { p ( x , y , z ) } , the work:
W = V { p } L T { f } L d V = ( V { p } L T [ N ] d V ) { d } L = { q } L T { d } L
Now it is possible to define the Lagrangian equation for a moving finite element:
L = E c E p + W + W c
The expressions of all these notions, in the case of the FEA of an elastic MBS, are presented in Appendix A.
The definition of the generalized momentum is:
{ p } L = { L { δ ˙ } L }
It results in:
{ p } L = L δ ˙ = { r ˙ O } L T [ m O i ] + V ρ ( { r } L T [ R ˙ ] T [ R ] [ N ] ) d V + ( V ρ [ N ] T [ R ˙ ] T [ R ] [ N ] d V ) { δ } L + [ m ] { δ ˙ } L
from where it is possible to obtain the vector of velocities { δ ˙ } L expressed as function of { p } L :
{ δ ˙ } L = [ m ] 1 ( { p } { r ˙ O } L T [ m O i ] V ρ ( { r } L T [ R ˙ ] T [ R ] [ N ] ) d V ( V ρ [ N ] T [ R ˙ ] T [ R ] [ N ] d V ) { δ } L )
which is useful in the further considerations.
The Lagrange equations, for a single finite element with three motions, are [31]:
[ m ] { δ ¨ } L + [ c ] { δ ˙ } L + ( [ k ] + [ k ( ε ) ] + [ k ( ω ) ] ) { δ } L = { q } L + { q } L { q i ( ε ) } L { q i ( ω ) } L [ m O i ] { r ¨ O } L

3. Hamiltonian Equations for a Finite Element

Hamilton’s formalism has been used by many researchers in the analysis of elastic systems and Hamilton’s equations have advantages in the case of mechanical systems with symmetries, which allow for obtaining prime integrals and, as a result, the size of the system decreases [37,38,39,40].
In our case, the Hamiltonian equation becomes:
H = { p } L T { δ ˙ } L L
where, for the Lagrangian, the rel. (6) is used.
Hamilton’s equations are [10]:
{ δ ˙ } L = { H { p } L } ; { p ˙ } L = { H { δ } L }
In Table 1, with the results presented in the Appendix A, Hamiltonian derivatives are shown. Finally, the motion equations in the Hamiltonian formulation, as a set of first-order differential systems, are obtained:
{ δ ˙ } L = [ m ] 1 ( { p } { r ˙ O } L T [ m O i ] V ρ ( { r } L T [ R ˙ ] T [ R ] [ N ] ) d V ( V ρ [ N ] T [ R ˙ ] T [ R ] [ N ] d V ) { δ } L )
{ p ˙ } L = { r ˙ O } L T [ R ] T [ R ˙ ] [ m O i ] + V ρ ( { r } L T [ R ˙ ] T [ R ˙ ] [ N ] ) d V + ( V ρ [ N ] T [ R ˙ ] T [ R ˙ ] [ N ] d V ) { δ } L + + V ρ [ N ] T [ R ˙ ] T [ R ] [ N ] d V { δ ˙ } L 1 2 [ k ] { δ } + { q } L + { q } L
Hamilton’s equations are attractive due to their simplicity and symmetry. Theoretically, these equations were analyzed from different points of view in different fields, from analytical mechanics to the geometry of vector spaces [9,23,37]. So, the way to solve it is the numerical way. In this case, the equations have an advantage because they are first-order differential equations and can be easily programmed and solved. In the case of second-order systems of differential equations (obtained if the Lagrange equation or other alternate analytical formulation is used), for the numerical solution, the transition to a first-order system of differential equations must be made [41]. Hamilton’s equations give us these equations directly.
In conclusion, Lagrange’s equations are second-order differential equations which, in order to be solved, are to be transformed into a system of first-order differential equations. Hamilton’s equations are first-order differential equations. In the form in which they are obtained, they can be used directly, without transformations, in a specific subroutine.
The literature in the studied field is vast. For Lagrange’s method, the works [42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61] can be used and, for Hamilton’s equations, [62,63,64,65,66,67,68,69].

4. Examples

An analysis of the advantages that Hamilton’s method could offer, compared to other classical methods and especially to Lagrange’s method, is difficult to do and would probably require a systematic and distinct effort. However, in the following, we will try to make a division between the two methods, applied in a few simple cases.
a. 
Rotating beam
Consider a bar of length equal to L, mass m and cross section A, which rotates about an axis (Figure 1), perpendicular to the bar at its end. The beam is divided into n finite elements, the bar nodes are numbered from 1 to n + 1. In our calculus, we study the axial deformation of the beam. Equation (10), obtained using the method of Lagrange’s equations, is used to write a program in Matlab for calculating axial deformations. In order to obtain the set of second-order differential equations applicable in this case, the parameters considered in the studied application are introduced. In order to apply Hamilton’s equations (first-order differential equations) for the studied bar, Equation (13) is used. In this case, too, a program is written in Matlab. A comparison between the results using these two programs is presented.
A sinusoidal excitation function will act at the right end of the beam. The results of the integration procedures are presented in Figure 2, Figure 3 and Figure 4. The two formalisms applied offer the same results and the integration seems to be stable. What is different is the integration time. In Figure 5, the diagram of the computing time considering the two methods is presented. In our application, the length is L = 1 m, the mass is 1 kg, the section A = 0.2 cm2 and the Young’s modulus is 210 GPa. The beam rotates with an angular speed of 15 rad/s. The excitation axial force has a period of T = π/100.
b. 
Beam in a general planar motion
An elastic beam of a mechanism with two degrees of freedom was considered. The beam had a planar motion. In order to facilitate the solution of the problem, the successive positions of some markers fixed on the bar were registered. Based on these records, the speeds and accelerations of the points where the markers were fixed were calculated and, based on them, the fields of the speed and acceleration for the beam were determined [22,35]. These results, obtained experimentally, were used as input data in Equations of motion (10) and (13). With the help of these values, integrations of these equations of motion were made.
Based on the mechanism presented in Figure 6, the positions of the markers were registered. Figure 7 shows some of the records. A coordinate system fixed to the bar was used to study the motion of different points on the beam. Figure 8 presents the acceleration field for two points on the bar (Table 2). In order to determine the computational times for performing an integration over a relatively short time interval (0.1 s), the values indicated in Table 3 were obtained.
Following the analysis of the two simple cases presented in this section, it can be concluded that Hamilton’s method seems to be more economical in terms of the computation time involved. However, the differences are not spectacular. With the increase in the studied system size, the difference between the two methods is accentuated. For this reason, it can be assumed that, for large systems, the method of Hamilton’s equations could become interesting.
A study dedicated to this problem will clarify more precisely the differences in the application of the two methods. The simple comparison, made in this paper, seems to indicate that the differences will increase if the dimensions of the systems increase, which would be consistent with the observation that computer times depend, for this type of application, exponentially on the size of the system.
c. 
Analysis of the influence of different parameters on the obtained solutions
Equations (14) and (17) were considered in order to determine the displacement of the beam presented in Figure 1.
A series of simulations of the bar rotating around a vertical axis, perpendicular to the bar axis, were made to determine the influence of different parameters on the behavior of the bar and the possibility of obtaining results. Thus, when comparing the computer time using Lagrange’s equations and Hamilton’s equations, when the accuracy of obtaining the solution is similar, the calculation time for the second method is 6–11% lower. The graphical representations of the two methods are similar. If we vary the section of the bar and consider different rotations of the bar, we obtain the graphs in Figure 9, Figure 10, Figure 11 and Figure 12.
Varying the excitation frequency, we obtain the results presented in Figure 13.
Studying the variation of the bar mass, we obtain Figure 14.

5. Conclusions and Discussions

The main step in the FEA of an elastic MBS is to obtain the equations of motion for the finite elements chosen for discretization. If the formalisms of the consecrated analytical mechanics are used, a system of second-order differential equations is obtained. Once these equations are written, the next step involves assembling the equations and solving them. These are properly done according the usual and well-known procedures. The system of equations that is finally obtained is a system of second-order differential equations, with the nodal coordinates being unknown. To solve this system, it must first be transformed into a first-order system by introducing additional unknowns. This transformation is a costly process in terms of computational time.
Hamilton’s equations offer, from the beginning, a system of first-order differential equations, with the unknowns consisting of generalized coordinates and generalized momenta. This proves to be an advantage as it is no longer necessary to transform the system of second-order equations into a system of first-order equations. As a result, the computational effort involved is smaller. The introduction of additional unknowns is done in the first step of modeling when the generalized moment is introduced. This method therefore proves to be better in terms of the simplicity of the solution.
Using Hamilton’s equations to obtain the equations of motion for a finite element (thus an approximation of the equations of motion of the continuous medium) has an advantage in terms of the number of formal operations involved and the computational time required to solve. The paper applies this formalism for writing equations of motion in the case of the three-dimensional motion of an elastic linear solid and makes some illustrative applications for applying the method to two relatively simple examples.
In addition to these obvious advantages, which make it easier to obtain systems of differential equations in the models used, the comparison we made in the paper seems to indicate that computer times, applying this formalism, are smaller, which is an additional advantage.

Author Contributions

Conceptualization, S.V., A.E.N. and M.M.; methodology, V.S., M.M. and A.E.N.; validation, A.E.N. and M.M.; formal analysis, S.V., A.E.N. and M.M.; investigation, S.V. and A.E.N.; resources, S.V.; data curation, A.E.N. and M.M.; writing—original draft preparation, S.V.; writing—review and editing, S.V., M.M. and A.E.N.; visualization, A.E.N. and M.M.; supervision, A.E.N., S.V. and M.M.; project administration, S.V. All authors have read and agreed to the published version of the manuscript and they have the same contribution.

Funding

This research received no external funding.

Acknowledgments

The second author is supported by the project ANTREPRENORDOC, in the framework of Human Resources Development Operational Programme 2014–2020, financed from the European Social Fund under the contract number 36355/23.05.2019 HRD OP /380/6/13—SMIS Code: 123847. We want to thank the reviewers who have read the manuscript carefully and have proposed pertinent corrections that have led to the improvement of our manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Kinetic Energy

The kinetic energy is:
E c = 1 2 V ρ { v M } G T { v M } G d V
Notations:
[ m ] = V ρ [ S ] T [ N ] S d V
[ m O i ] = V ρ [ S ] T d V
{ q i ( ε ) } L = V ρ [ S ] T [ ε ] L { r } L d V
{ q i ( ω ) } L = V ρ [ S ] T [ ω ] L [ ω ] L { r } L d V
[ k ( ε ) ] = V ρ [ S ] T [ ε ] [ S ] d V
[ k ( ε ) ] = V ρ [ S ] T [ ω ] L [ ω ] L [ S ] d V
[ c ] = V ρ [ S ] T [ ω ] L [ S ] d V
{ m i x } = V ρ [ S ( i ) ] T x d V
{ m i y } = V ρ [ S ( i ) ] T y d V
{ m i z } = V ρ [ S ( i ) ] T z d V

Appendix A.2. Potential Energy

From the relation:
E p = 1 2 V { σ } T { ε } d V
where
{ ε } contains the distinct component of the strain tensor;
{ σ } contains the distinct component of the stress tensor.
Using the Hooke law and some elementary transformation, we finally obtain:
E p = 1 2 V   { d } L T [ k ] { d } d V

Appendix A.3. Work

Concentrated forces acting in knots, { q } L and volume forces { p } = { p ( x , y , z ) } give a mechanical work:
W c = { q } L T { d } L
and:
W = V { p } L T { f } L d V = ( V { p } L T [ S ] d V ) { d } L = { q } L T { d } L

Appendix A.4. Lagrangian Equation

L = E c E p + W + W c
It obtains:
L = 1 2 V ρ ( { r ˙ O } L T { r ˙ O } L + 2 { r ˙ O } L T [ R ] T [ R ˙ ] { r } L + 2 { r ˙ O } L T [ R ] T [ R ˙ ] [ S ] { d } L + 2 { r ˙ O } L T [ N ] { δ ˙ } L ) d V + 1 2 V ρ ( { r } L T [ R ˙ ] T [ R ˙ ] { r } L + 2 { r } L T [ R ˙ ] T [ R ˙ ] [ N ] { δ } L + 2 { r } L T [ R ˙ ] T [ R ] [ N ] { δ ˙ } L ) d V + 1 2 V ρ ( { δ } L T [ N ] T [ R ˙ ] T [ R ˙ ] [ N ] { δ } L + 2 { δ } L T [ N ] T [ R ˙ ] T [ R ] [ N ] { δ ˙ } L + { δ ˙ } L T [ N ] T [ N ] { δ ˙ } L ) d V 1 2 V   { δ } L T [ k ] { δ } L d V + { q } L T { δ } L + { q } L T { δ } L

Appendix A.5. Generalized Momentum

The formula for the generalized momentum is:
{ p } L = { L { d ˙ } L }
It results in:
{ p } L = L δ ˙ = { r ˙ O } L T [ m O i ] + V ρ ( { r } L T [ R ˙ ] T [ R ] [ N ] ) d V + ( V ρ [ N ] T [ R ˙ ] T [ R ] [ N ] d V ) { δ } L + [ m ] { δ ˙ } L

Appendix A.6. Lagrange’s Equation

Lagrange’s equation applied for a finite element is:
d d t { L δ ˙ } L { L δ } L = 0
It obtains:
[ m ] { d ¨ } L + [ c ] { d ˙ } L + ( [ k ] + [ k ( ε ) ] + [ k ( ω ) ] ) { d } L = { q } L + { q } L { q i ( ε ) } L { q i ( ω ) } L [ m O i ] { r ¨ O } L

References

  1. Lagrange, J.L. Mecanique Analitique; Vve Desaint: Paris, France, 1788. [Google Scholar]
  2. Gibbs, J.W. On the fundamental formulae of dynamics. Am. J. Math. 1879, 2, 49–64. [Google Scholar] [CrossRef]
  3. Appell, P. Sur une forme générale des equations de la dynamique. J. Reine Angew. Math. 1900, 1900, 310–319. [Google Scholar]
  4. Hamilton, R.W. Second Essay on a General Method in Dynamics. In Phylosophical Transaction of the Royal Society, Part I; Wilkins, D.R., Ed.; Royal Society: London, UK, 1835; pp. 95–144. [Google Scholar]
  5. Maggi, G.A. Princippi della Teoria Matematica del Movimento dei Corpi: Corso di Meccanica Razionale; Ulrico Hoepli: Milano, Italy, 1896. [Google Scholar]
  6. Kane, T.T. Dynamic of Nonholonomic Systems, Transaction of ASME. J. Appl. Mech. 1961, 83, 574–578. [Google Scholar] [CrossRef]
  7. Haug, E.J. Extension of Maggi and Kane Equations to Holonomic Dynamic Systems. J. Comput. Nonlinear Dyn. 2018, 13, 121003. [Google Scholar] [CrossRef]
  8. Malvezzi, F.; Matarazzo Orsino, R.M.; Hess Coelho, T.A. Lagrange’s, Maggi’s and Kane’s equations to the dynamic modelling of serial manipulator. In DINAME 2017, Proceedings of the XVII International Symposium on Dynamic Problems of Mechanics, 5–10 March 2017; ABCM: São Sebastião, SP, Brazil, 2017. [Google Scholar]
  9. Wasfy, T.M.; Noor, A.K. Computational strategies for flexible multibody systems. ASME Appl. Mech. Rev. 2003, 56, 553–613. [Google Scholar] [CrossRef]
  10. Ursu-Fisher, N. Elements of Analytical Mechanics; House of Science Book Press: Cluj-Napoca, Romania, 2015. [Google Scholar]
  11. Negrean, I.; Kacso, K.; Schonstein, C.; Duca, A. Mechanics. Theory and Applications; UTPRESS: Cluj-Napoca, Romania, 2012. [Google Scholar]
  12. Mirtaheri, S.M.; Zohoor, H. The Explicit Gibbs-Appell Equations of Motion for Rigid-Body Constrained Mechanical System. In Proceedings of the 6th RSI International Conference on Robotics and Mechatronics (ICRoM), Tehran, Iran, 23–25 October 2018; pp. 304–309. [Google Scholar]
  13. Korayem, M.H.; Dehkordi, S.F. Motion equations of cooperative multi flexible mobile manipulator via recursive Gibbs-Appell formulation. Appl. Math. Model. 2019, 65, 443–463. [Google Scholar] [CrossRef]
  14. Vlase, S.; Negrean, I.; Marin, M.; Scutaru, M.L. Energy of Accelerations Used to Obtain the Motion Equations of a Three-Dimensional Finite Element. Symmetry 2020, 12, 321. [Google Scholar] [CrossRef] [Green Version]
  15. Amengonu, Y.H.; Kakad, Y.P. Dynamics and control for Constrained Multibody Systems modeled with Maggi’s equation: Application to Differential Mobile Robots Part II. In Proceedings of the 27th International Conference on CADCAM, Robotics and Factories of the Future, London, UK, 22–24 July 2014; Volume 65, p. 012018. [Google Scholar]
  16. Ladyzynska-Kozdras, E. Application of the Maggi equations to mathematical modeling of a robotic underwater vehicle as an object with superimposed non-holonomic constraints treated as control laws. Mechatron. Syst. Mech. Mater. 2012, 180, 152–159. [Google Scholar] [CrossRef]
  17. De Jalon, J.G.; Callejo, A.; Hidalgo, A.F. Efficient Solution of Maggi’s Equations. J. Comput. Nonlinear Dyn. 2012, 7, 021003. [Google Scholar] [CrossRef]
  18. Vlase, S.; Marin, M.; Scutaru, M.L. Maggi’s Equations Used in the Finite Element Analysis of the Multibody Systems with Elastic Elements. Mathematics 2020, 8, 399. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, J.T.; Huston, R.L. Kane’s equations with undetermined multipiers-application to constrained multibody systems. ASME J. Appl. Mech. 1987, 54, 424–429. [Google Scholar] [CrossRef]
  20. Kane, T.R.; Levinson, D.A. Formulation of Equations of Motion for Complex Spacecraft. J. Guid. Control 1980, 3, 99–112. [Google Scholar] [CrossRef]
  21. Kane, T.R.; Levinson, D.A. Multibody Dynamics. J. Appl. Mech. 1983, 50, 1071–1078. [Google Scholar] [CrossRef]
  22. Vlase, S.; Negrean, I.; Marin, M.; Nastac, S. Kane’s Method-Based Simulation and Modeling Robots with Elastic Elements, Using Finite Element Method. Mathmatics 2020, 8, 805. [Google Scholar] [CrossRef]
  23. Dowell, E. Hamilton’s principle and Hamilton’s equations with holonomic and non-holonomic constraints. Nonlinear Dyn. 2017, 88, 1093–1097. [Google Scholar] [CrossRef]
  24. Tong, M.M. Flexible Multibody Dynamics Formulation by Hamilton’s Equations. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Vancouver, BC, Canada, 12–18 November 2010; Volume 8, pp. 725–734. [Google Scholar]
  25. Qiao, Y.F.; Zhang, Y.L.; Han, G.C. Form invariance of Hamilton’s canonical equations of a nonholonomic mechanical system. Acta Phys. Sin. 2003, 52, 1051–1056. [Google Scholar]
  26. Erdman, A.G.; Sandor, G.N.; Oakberg, A. A General Method for Kineto-Elastodynamic Analysis and Synthesis of Mechanisms. J. Eng. Ind. 1972, 94, 1193–1205. [Google Scholar] [CrossRef]
  27. Bahgat, B.M.; Willmert, K.D. Finite Element Vibrational Analysis of Planar Mechanisms. Mech. Mach. Theory 1976, 11, 47–71. [Google Scholar] [CrossRef]
  28. Bagci, C. Elastodynamic Response of Mechanical Systems using Matrix Exponential Mode Uncoupling and Incremental Forcing Techniques with Finite Element Method. In Proceedings of the Sixth Word Congress on Theory of Machines and Mechanisms, New Delhi, India, 15–20 December1983; p. 472. [Google Scholar]
  29. Cleghorn, W.L.; Fenton, E.G.; Tabarrok, K.B. Finite Element Analysis of High-Speed Flexible Mechanisms. Mech. Mach. Theory 1981, 16, 407–424. [Google Scholar] [CrossRef]
  30. Vlase, S.; Marin, M.; Öchsner, A.; Scutaru, M.L. Motion equation for a flexible one-dimensional element used in the dynamical analysis of a multibody system. Contin. Mech. Thermodyn. 2019. [Google Scholar] [CrossRef]
  31. Vlase, S. Dynamical Response of a Multibody System with Flexible Elements with a General Three-Dimensional Motion. Rom. J. Phys. 2012, 57, 676–693. [Google Scholar]
  32. Deu, J.-F.; Galucio, A.C.; Ohayon, R. Dynamic responses of flexible-link mechanisms with passive/active damping treatment. Comput. Struct. 2008, 86, 258–265. [Google Scholar] [CrossRef] [Green Version]
  33. Hou, W.; Zhang, X. Dynamic analysis of flexible linkage mechanisms under uniform temperature change. J. Sound Vib. 2009, 319, 570–592. [Google Scholar] [CrossRef]
  34. Neto, M.A.; Ambrósio, J.A.C.; Leal, R.P. Composite materials in flexible multibody systems. Comput. Methods Appl. Mech. Eng. 2006, 195, 6860–6873. [Google Scholar] [CrossRef] [Green Version]
  35. Itu, C.; Oechsner, A.; Vlase, S.; Marin, M. Improved rigidity of composite circular plates through radial ribs. Proc. Inst. Mech. Eng. Part L—J. Mater. Des. Appl. 2019, 233, 1585–1593. [Google Scholar] [CrossRef]
  36. Marin, M.; Ellahi, R.; Chirilă, A. On solutions of Saint-Venant’s problem for elastic dipolar bodies with voids. Carpathian J. Math. 2017, 33, 219–232. [Google Scholar]
  37. Figueroa, H.; Carinena, J.F. Hamiltonian versus Lagrangian formulations of supermechanics. J. Phys. A Math. Gen. 1997, 30, 2705–2724. [Google Scholar]
  38. Naudet, J.; Lefeber, D.; Daerden, F.; Terze, Z. Forward Dynamics of Open-Loop Multibody Mechanisms Using an Efficient Recursive Algorithm Based on Canonical Momenta. Multibody Syst. Dyn. 2003, 10, 45–59. [Google Scholar] [CrossRef]
  39. Kim, J.; Kim, D. A quadratic temporal finite element method for linear elastic structural dynamics. Math. Comput. Simul. 2015, 117, 68–88. [Google Scholar] [CrossRef]
  40. Xue, Y.; Weng, D.W.; Chen, L.Q. Methods of analytical mechanics for exact Cosserat elastic rod dynamics. Acta Phys. Sin. 2013, 62, 044601. [Google Scholar] [CrossRef]
  41. MATLAB. Available online: https://de.mathworks.com/help/matlab/ (accessed on 29 July 2020).
  42. Dubowsky, S.; Gardner, T.N. Dynamic interactions of link elasticity and clearance connections in planar mechanical systems. ASME J. Eng. Ind. 1975, 97, 652–661. [Google Scholar] [CrossRef]
  43. Midha, A.; Erdman, A.G.; Forhib, D.A. Finite element approach to mathematical modeling of high-speed elastic linkages. Mech. Mach. Theory 1978, 13, 603–618. [Google Scholar] [CrossRef]
  44. Midha, A.; Erdman, A.G.; Forhib, D.A. A computationally efficient numerical algorithm for the transient response of high-speed elastic linkages. ASME J. Mech. Des. 1979, 101, 138–148. [Google Scholar] [CrossRef]
  45. Blejwas, T.E. The simulation of elastic mechanisms using kinematic constraints and Lagrange multipliers. Mech. Mach. Theory 1981, 16, 441–445. [Google Scholar] [CrossRef]
  46. Bricout, J.N.; Debus, J.C.; Micheau, P. A finite element model for the dynamics of flexible manipulators. Mech. Mach. Theory 1990, 25, 119–128. [Google Scholar] [CrossRef]
  47. Meirovitch, L.; Kwak, M.K. Dynamics and control of spacecraft with retargeting flexible antennas. J. Guid. Control Dyn. 1990, 13, 241–248. [Google Scholar] [CrossRef]
  48. Fattah, A.; Angeles, J.; Misra, A.K. Dynamics of a 3-DOF spatial parallel manipulator with flexible links. In Proceedings of the 1995 IEEE International Conference on Robotics and Automation, Nagoya, Japan, 21–27 May 1995; pp. 627–632. [Google Scholar]
  49. Pereira, M.S.; Ambrosio, J.A.C.; Dias, J.P. Crashworthiness analysis and design using rigid-flexible multibody dynamics with application to train vehicles. Int. J. Numer. Methods Eng. 1997, 40, 655–687. [Google Scholar] [CrossRef]
  50. Turcic, D.A.; Midha, A. Generalized equations of motion for the dynamic analysis of elastic mechanism systems. ASME J. Dyn. Syst. Meas. Control 1984, 106, 243–248. [Google Scholar] [CrossRef]
  51. Turcic, D.A.; Midha, A. Dynamic analysis of elastic mechanism systems, Part I: Applications. ASME J. Dyn. Syst. Meas. Control 1984, 106, 249–254. [Google Scholar] [CrossRef]
  52. Bakr, E.M.; Shabana, A.A. Geometrically nonlinear analysis of multibody systems. Comput. Struct. 1986, 23, 739–751. [Google Scholar] [CrossRef]
  53. Changizi, K.; Shabana, A.A. A recursive formulation for the dynamics analysis of open loop deformable multibody systems. J. Appl. Acoust. 1988, 55, 687–693. [Google Scholar] [CrossRef]
  54. Modi, V.J.; Suleman, A.; Ng, A.C.; Morita, Y. An approach to dynamics and control of orbiting flexible structures. Int. J. Numer. Methods Eng. 1991, 32, 1727–1748. [Google Scholar] [CrossRef]
  55. Pereira, M.S.; Proenca, P.L. Dynamic analysis of spatial flex ible multibody systems using joint coordinates. Int. J. Numer. Methods Eng. 1991, 32, 1799–1821. [Google Scholar] [CrossRef]
  56. Pradhan, S.; Modi, V.J.; Misra, A.K. Order N formulation for flexible multibody systems in tree topology: Lagrangian approach. J. Guid. Control Dyn. 1997, 20, 665–672. [Google Scholar] [CrossRef]
  57. Metaxas, D.; Koh, E. Flexible multibody dynamics and adaptive finite element techniques for model synthesis and estimation. Comput. Methods Appl. Mech. Eng. 1996, 136, 1–25. [Google Scholar] [CrossRef]
  58. Huang, S.J.; Wang, T.Y. Structural dynamics analysis of spatial robots with finite element method. Comput. Struct. 1993, 46, 703–716. [Google Scholar] [CrossRef]
  59. Smaili, A.A. A three-node finite beam element for dynamic analysis of planar manipulators with flexible joints. Mech. Mach. Theory 1993, 28, 193–206. [Google Scholar] [CrossRef]
  60. Meek, J.L.; Liu, H. Nonlinear dynamics analysis of flexible beams under large overall motions and the flexible manipulator simulation. Comput. Struct. 1995, 56, 1–14. [Google Scholar] [CrossRef]
  61. Pan, Y.C.; Scott, R.A.; Ulsoy, G.A. Dynamic modeling and simulation of flexible robots with prismatic joints. ASME J. Mech. Des. 1990, 112, 307–314. [Google Scholar] [CrossRef]
  62. Cavin, R.K.; Dusto, A.R. Hamilton’s principle: Finite element methods and flexible body dynamics. AIAA J. 1977, 15, 1684–1690. [Google Scholar] [CrossRef]
  63. Chang, L.W.; Hamilton, J.F. Kinematics of robotic manipulators with flexible links using and equivalent rigid link system ERLS model. ASME J. Dyn. Syst. Meas. Control 1991, 113, 48–54. [Google Scholar] [CrossRef]
  64. Serna, M.A. A simple and efficient computational approach for the forward dynamics of elastic robots. J. Rob. Syst. 1989, 6, 363–382. [Google Scholar] [CrossRef]
  65. Fung, R.F. Dynamic analysis of the flexible connecting rod of slider-crank mechanism with time-dependent boundary effect. Comput. Struct. 1997, 63, 79–90. [Google Scholar] [CrossRef]
  66. Huang, Y.; Lee, C.S.G. Generalization of Newton-Euler formulations of dynamic equations to nonrigid manipulators. ASME J. Dyn. Syst. Meas. Control 1988, 110, 308–315. [Google Scholar] [CrossRef]
  67. Iura, M.; Atluri, S.N. Dynamic analysis of planar flexible beams with finite rotations by using inertial and rotating frames. Comput. Struct. 1955, 55, 453–462. [Google Scholar] [CrossRef]
  68. Bauchau, O.A.; Damilano, G.; Theron, N.J. Numerical integration of non-linear elastic multibody systems. Int. J. Numer. Methods Eng. 1955, 38, 2727–2751. [Google Scholar] [CrossRef]
  69. Bauchau, O.A.; Choi, J.-Y.; Bottasso, C.L. On the modeling of shells in multibody dynamics. Multibody Syst. Dyn. 2002, 8, 459–489. [Google Scholar] [CrossRef]
Figure 1. Rotating beam.
Figure 1. Rotating beam.
Mathematics 08 01401 g001
Figure 2. Axial displacement of the node 3, 7, 10 (10 finite elements, 11 nodes), without damping.
Figure 2. Axial displacement of the node 3, 7, 10 (10 finite elements, 11 nodes), without damping.
Mathematics 08 01401 g002
Figure 3. Axial displacement of the node 3, 7, 10 (10 finite elements, 11 nodes) with damping.
Figure 3. Axial displacement of the node 3, 7, 10 (10 finite elements, 11 nodes) with damping.
Mathematics 08 01401 g003
Figure 4. Displacements of all nodes, with damping.
Figure 4. Displacements of all nodes, with damping.
Mathematics 08 01401 g004
Figure 5. The computational time to integrate the motion equations in the two cases.
Figure 5. The computational time to integrate the motion equations in the two cases.
Mathematics 08 01401 g005
Figure 6. The studied mechanism and the positioning of the markers on the bar.
Figure 6. The studied mechanism and the positioning of the markers on the bar.
Mathematics 08 01401 g006
Figure 7. Recording the successive positions of the markers.
Figure 7. Recording the successive positions of the markers.
Mathematics 08 01401 g007
Figure 8. Acceleration of markers 1 and 2. The measurements for one rotation of the shaft have the same color.
Figure 8. Acceleration of markers 1 and 2. The measurements for one rotation of the shaft have the same color.
Mathematics 08 01401 g008
Figure 9. Cross-section area A = 0.00001 m2.
Figure 9. Cross-section area A = 0.00001 m2.
Mathematics 08 01401 g009
Figure 10. Cross-section area A = 0.00002 m2.
Figure 10. Cross-section area A = 0.00002 m2.
Mathematics 08 01401 g010
Figure 11. Cross-section area A = 0.00004 m2.
Figure 11. Cross-section area A = 0.00004 m2.
Mathematics 08 01401 g011
Figure 12. Cross-section area A = 0.00006 m2.
Figure 12. Cross-section area A = 0.00006 m2.
Mathematics 08 01401 g012
Figure 13. Excitation frequency 50/π, 100/π, 150/π.
Figure 13. Excitation frequency 50/π, 100/π, 150/π.
Mathematics 08 01401 g013
Figure 14. Mass of the bar 1 kg, 0.5 kg, 0.25 kg.
Figure 14. Mass of the bar 1 kg, 0.5 kg, 0.25 kg.
Mathematics 08 01401 g014
Table 1. The derivatives of Hamiltonian equations.
Table 1. The derivatives of Hamiltonian equations.
Term of Hamiltonian Equation { p ˙ } = H δ = L δ
1 2 V ρ ( { r ˙ O } L T { r ˙ O } L ) d V = 1 2 m v O 2 0
V ρ ( { r ˙ O } L T [ R ] T [ R ˙ ] { r } L ) d V 0
{ r ˙ O } L T [ R ] T [ R ˙ ] ( V ρ [ N ] d V ) { δ } L = { r ˙ O } L T [ R ] T [ R ˙ ] [ m O i ] { δ } L { r ˙ O } L T [ R ] T [ R ˙ ] [ m O i ]
{ r ˙ O } L T V ρ ( [ N ] ) d V { δ ˙ } L = { r ˙ O } L T [ m O i ] { δ ˙ } L 0
1 2 V ρ ( { r } L T [ R ˙ ] T [ R ˙ ] { r } L ) d V 0
V ρ ( { r } L T [ R ˙ ] T [ R ˙ ] [ N ] ) d V { δ } L V ρ ( { r } L T [ R ˙ ] T [ R ˙ ] [ N ] ) d V
V ρ ( { r } L T [ R ˙ ] T [ R ] [ N ] ) d V { δ ˙ } L 0
1 2 { δ } L T ( V ρ [ N ] T [ R ˙ ] T [ R ˙ ] [ N ] d V ) { δ } L ( V ρ [ N ] T [ R ˙ ] T [ R ˙ ] [ N ] d V ) { δ } L
V ρ ( { δ } L T [ N ] T [ R ˙ ] T [ R ] [ N ] { δ ˙ } L ) d V V ρ [ N ] T [ R ˙ ] T [ R ] [ N ] d V { δ ˙ } L
1 2 { δ ˙ } L T ( V ρ [ N ] T [ N ] d V ) { δ ˙ } L = 1 2 { δ ˙ } L T [ m ] { δ ˙ } L 0
1 2 { δ } L T [ k ] { δ } + { q } L T { δ } L + { q } L T { δ } L 1 2 [ k ] { δ } + { q } L + { q } L
Table 2. Values of the velocity of markers 1 and 2.
Table 2. Values of the velocity of markers 1 and 2.
Node 1Node 2
v ( m m / s ) a ( m m / s 2 ) v ( m m / s ) a ( m m / s 2 )
v m i n v m a x a m i n a m a x v m i n v m a x a m i n a m a x
306.4894.219238991252.4670.214807325
263.3100617639962242.3845.8899.59851
253.81194244814,580325.91015223314,490
346.41243244513,980259.2996.5242112,610
315.51538295716,330355.21296300014,260
Table 3. Calculation times for the two methods.
Table 3. Calculation times for the two methods.
Number of elements510152025
Lagrange’s method. Time (s)3184176291460
Hamilton’s method. Time (s)2868139224322

Share and Cite

MDPI and ACS Style

Vlase, S.; Nicolescu, A.E.; Marin, M. New Analytical Model Used in Finite Element Analysis of Solids Mechanics. Mathematics 2020, 8, 1401. https://doi.org/10.3390/math8091401

AMA Style

Vlase S, Nicolescu AE, Marin M. New Analytical Model Used in Finite Element Analysis of Solids Mechanics. Mathematics. 2020; 8(9):1401. https://doi.org/10.3390/math8091401

Chicago/Turabian Style

Vlase, Sorin, Adrian Eracle Nicolescu, and Marin Marin. 2020. "New Analytical Model Used in Finite Element Analysis of Solids Mechanics" Mathematics 8, no. 9: 1401. https://doi.org/10.3390/math8091401

APA Style

Vlase, S., Nicolescu, A. E., & Marin, M. (2020). New Analytical Model Used in Finite Element Analysis of Solids Mechanics. Mathematics, 8(9), 1401. https://doi.org/10.3390/math8091401

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop