Next Article in Journal
Advanced Bimodal Skew-Symmetric Distributions: Methodology and Application to Cancer Cell Protein Data
Previous Article in Journal
BDAC: Boundary-Driven Approximations of K-Cliques
Previous Article in Special Issue
A Non-Second-Gradient Model for Nonlinear Electroelastic Bodies with Fibre Stiffness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Linear Elastic Beam Deformations with Four-Parameter Timoshenko Beam Element Considering Through-the-Thickness Stretch Parameter and Reduced Integration

by
Nasser Firouzi
1,* and
Ahmed S. M. Alzaidi
2
1
Institute of Structural Mechanics, Bauhaus-University Weimar, 99423 Weimar, Germany
2
Department of Mathematics and Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(8), 984; https://doi.org/10.3390/sym16080984
Submission received: 14 July 2024 / Revised: 27 July 2024 / Accepted: 30 July 2024 / Published: 2 August 2024
(This article belongs to the Special Issue Symmetry in Finite Element Modeling and Mechanics)

Abstract

:
In this work, non-linear elastic deformations of beams are investigated. The kinematics of the beam is derived based on an element with four-parameter containing a through-the-thickness stretch parameter to avoid Poisson locking. Moreover, the Kirchhoff-Saint Venant model is used to derive kinetic quantities. Next, a non-linear FE formula in Total Lagrangian form is obtained, and three-node beam element with two-node reduced integration is employed to avoid shear locking. Finally, to evaluate the performance of the derived formulations, some examples are provided. The results prove that the current formulation is in very good agreement with those available in the literature. More importantly, the formulation is capable of predicting the experimental results with high accuracy.

1. Introduction

1.1. Importance and Necessity of the Problem

Beams, fundamental structural elements, are very important piece across various engineering disciplines. Their applications span aerospace structures soaring through the atmosphere, high-rise buildings reaching for the sky, bridges connecting landscapes, and even the miniaturized marvels of micro- and nano-electromechanical systems (MEMS/NEMS) [1,2,3,4,5,6]. Biomechanics, too, finds application for beams in analyzing the intricate skeletal structures of living organisms [7,8]. Within this broad spectrum of applications, the analysis of beams subjected to large deformations has gained increasing importance. This growing focus stems from the recognition that real-world structures encounter a diverse range of loading scenarios. Static loads, dynamic forces, and even environmental factors can all induce significant deformations in beams [9,10]. Developing a more accurate framework to analyze these large deformations becomes crucial for designing safer and more reliable structures. While traditional beam theories, such as the Euler-Bernoulli model, have provided valuable tools for engineers, their limitations become apparent when dealing with scenarios involving substantial deformations. This necessitates the exploration of more sophisticated analytical approaches to accurately predict the behavior of beams under such conditions. Classical beam theories often assume linear material behavior. However, many engineering materials exhibit nonlinear stress-strain relationships. These limitations often stem from assumptions of small strains and linear material behavior, which can lead to inaccurate predictions in scenarios involving significant geometric changes. The Timoshenko beam theory, which accounts for shear deformation and rotational effects, provides a more comprehensive approach to modeling such scenarios [11,12].
The Timoshenko beam theory offers a more refined approach compared to the Euler-Bernoulli model that with additional considerations become particularly important for stockier beams or those experiencing substantial shear forces [13,14]. Within the domain of beam mechanics, achieving a nuanced understanding of the behavior exhibited by beams undergoing significant deflections necessitates the strategic deployment of two key elements: nonlinear constitutive models and geometrically exact formulations [15,16]. Hyperelasticity, a theoretical framework adept at describing the elastic response of materials even under substantial strains, emerges as a powerful tool for this endeavor [17]. This framework excels at capturing the nonlinear stress-strain relationships observed in a diverse array of materials, encompassing rubbers, polymers, and even the intricate biological tissues found within living organisms [18,19]. When applied to beam analysis, hyperelastic constitutive laws enable the accurate prediction of the complex interplay between material and geometric nonlinearities. However, directly solving the governing equations arising from the Timoshenko beam theory, especially when coupled with Hyperplastic and nonlinear material behavior, can be mathematically challenging.
The finite element method (FEM) offers a powerful tool for numerical analysis, enabling the exploration of complex deformation patterns in beams [20,21]. By discretizing the beam into smaller elements, the FEM transforms the complex problem into a system of algebraic equations that can be efficiently solved using numerical methods. This approach offers unparalleled flexibility in handling complex geometries, boundary conditions, and material properties, making it ideally suited for analyzing the large deformation response of beams even when considering nonlinear material and Hyperplastic material properties [22]. By incorporating the FEM into the framework, we can efficiently solve the governing equations arising from the Timoshenko beam theory with hyperplastic material considerations.

1.2. A Review on Literature

Several studies that have explored the finite element analysis of Timoshenko beams. Meier et al. [23] compared two geometrically exact finite element formulations for slender beams, namely the Kirchhoff-Love theory and the Simo-Reissner theory. Their study provided valuable insights into the accuracy and computational efficiency of these formulations for analyzing large deformations in beams. They highlighted the importance of considering the geometric nonlinearity in the finite element analysis of slender beams, which is particularly relevant to the study of Timoshenko hyperplastic beams. Cazzani et al. [24] performed an analytical assessment of finite element and isogeometric analyses of the whole spectrum of Timoshenko beams. Their study provided insights into the comparative performance of finite element and isogeometric analyses for modeling the behavior of Timoshenko beams under various loading conditions. The findings of this study are particularly relevant to the exploration of large deformations in Timoshenko hyperplastic beams, as they offer valuable information on the capabilities and limitations of different analysis methods. Sakai et al. [25] present a 3D elastic beam model for form-finding and analysis of elastic gridshells subjected to bending deformation at the self-equilibrium state. The accuracy of the proposed method utilizing the dynamic relaxation method is confirmed when compared to the results obtained by finite element analysis. In a study by Numanoğlu and his research group [26], a new eigenvalue problem solver was developed for the thermo-mechanical vibration of Timoshenko nanobeams using an innovative nonlocal finite element method. Finite element analyses of micro- and nanobeams were the topic of several studies [27,28,29,30]. Firouzi et al. [31] studied non-linear free vibration of Timoshenko beams using direct integration scheme. Moreover, Firouzi et al. [32] investigated large deformation of electro-visco-hyperelastic beams under electric charge. Wang et al. [33] used homotopy analysis method to analyze the deformation of a cantilever under a point load. Moreover, Beléndez et al. [34] performed an experiment on large deformation of a cantilever beam.
Cai et al. [35] proposed modeling of an active metabeam for broadband non-reciprocal wave suppression. Sahmani et al. [36] studied nonlinear in-plane buckling of FG porous microbeams by isogeometric analysis. Geometrically nonlinear stability of restrained nanobeams utilizing Stokes’ transformation was analyzed by Uzun et al. [37]. The role of homogenisation models in theoretically analysing the vibrations of three-layered microbeam was investigated by Yee et al. [38]. Settimi andLenci [39] analyzed free wave propagation of periodic flexural waves on an infinite elastic beam. Żur et al. [40] studied large deformations of hyperelastic beams using neo-Hookean material model. Bensaid et al. [41] studied free vibration and stability response of bidirectional graded material beams. Study of multiple-input multiple-output active vibration control of a composite sandwich beam by fractional order positive position feedback was done by Hameury et al. [42]. Liu et al. [43] studied free vibration of composite beams using Carrera unified formulation model. Pagani et al. [44] studied vibration of thin-walled slender structures by higher-order beam finite element.
While significant progress has been made in exploring such deformations through finite element analysis, knowledge gaps persist and necessitate further investigation. One potential future research direction is the development of a model that can accurately capture the behavior of Timoshenko beams under large deformations. Furthermore, the validation of finite element simulations through experimental testing of Timoshenko beams would enhance the confidence in the numerical results and contribute to the advancement of this research area.
This study demonstrates the progress and potential for further advancements in the finite element exploration of large deformations in Timoshenko beams. By adopting this framework, we establish a more accurate representation of the beam’s response under diverse loading conditions. These findings contribute significantly to our comprehension of the mechanical behavior exhibited by Timoshenko beams. Furthermore, they underscore the necessity for ongoing research endeavors to bridge the remaining knowledge gaps within this domain. Future investigations could strategically target free-form Timoshenko curved beams, with a particular emphasis on elucidating the nonlinear effects arising from large deformations. The nonlinear finite element analysis is developed, and through-the-thickness parameter is used to circumvent Poisson locking. Moreover, reduced integration is utilized to avoid shear locking. The experimental validation is also investigated.
The rest of the paper is organized as follows. In Section 2, the kinematics and kinetics of the four-parameter beam element are derived. Section 3 is allocated to derive FE formulation. In Section 4, some examples are provided to explore the applicability of the formulation. Finally, conclusion is discussed in Section 5.

2. Kinematics of the Four-Parameter Beam

An initially straight beam of which length, thickness and width are L , h 0 and b is considered. A fixed Cartesian coordinate system { X , Z , Y } with basis vectors { e 1 , e 2 , e 3 } is located in the neutral surface of the beam, as displayed in Figure 1.
The reference position of each material point can be described as
X = X   e 1 + Z   e 3 ,
Moreover, the position of every arbitrary point after deformation can be expressed as [32]:
x = x 0 ( X ) + [ Z + Z 2 κ ( X ) ] r ( X ) ,
where x 0 , κ and r are, respectively, the position in the centerline, the through-the-thickness stretch parameter, and the director vector. It should be mentioned that, the director vector is a unit vector in terms of rotation of the cross-section, φ is given by
r = sin φ   e 1 + cos φ   e 2 .
The deformation gradient is obtained via the following relation:
F = Grad ( x ) = x X F 0 + Z F 1
Considering { u , v , φ , κ } as unknown parameters, this expression is a four-parameter Timoshenko beam element. Moreover, the right Cauchy-Green tensor can be derived as follows:
C = C 0 + Z C 1 ,           C 0 = F 0 T F 0 ,         C 1 = F 0 T F 1 +   F 1 T F 0 .    
Finally, the Green-Lagrange tensor can be expressed via the following relation:
E = 1 2 ( C I )   E 0 + Z E 1 ,             E 0 = 1 2 ( C 0 I ) ,         E 1 = 1 2 C 1 .
In this work, the well-known Kirchhoff-St Venant material model is used as follows:
ψ = λ   [ tr ( E ) ] 2 + μ   tr E 2 ,
where Lame’ constants, λ and μ are defined by
  λ = E υ ( 1 + υ ) ( 1 2 υ ) ,                 μ = E 2 ( 1 + υ ) .
In the latter relation, E and υ are, respectively, the Young’s modulus and the Poisson’s ratio. In this case, the 2nd Piola-Kirchhoff stress tensor is expressed as:
Π = λ ( tr E ) I + 2 μ E .
In what follows, the stress is written in the linear format as follows:
Π = Π 0 + Z Π 1 ,
in which Π 0 and Π 1 are calculated as
Π 0 = λ ( tr E 0 ) I + 2 μ E 0 ,         Π 1 = λ ( tr E 1 ) I + 2 μ E 1 .
Using Equations (10) and (11), the fourth-order elasticity tensors are derived as
= 0 + Z 1 ,
0 = λ I I + μ ( E 0 ¯ E 0 + E 0 _ E 0 ) ,     1 = λ I I + μ ( E 1 ¯ E 1 + E 1 _ E 1 ) ,
where ¯ and _ are non-standard tensor product given by
[ ( ) ¯ ( ) ] i j k l = ( ) i k ( ) j l ,         [ ( ) _ ( ) ] i j k l = ( ) i l ( ) j k .
It should be mentioned that 0 and 0 are the linearized tensors of the total fourth-order elasticity tensor.

3. Finite Element Analysis

A finite element formulation for elastic deformation of 4-parameter beam in nonlinear regime is derived in this section. Based on principle of virtual work, the variation of internal energy is equal to the variation of work done by virtual external force
D U int = D U ext ,
where D U int and D U ext are the variations of internal energy and external work, respectively, and can be derived as follows:
D U int = V 0 D U   d V 0 = L A Π : D E   d A 0 d X ,
D U ext = L q 1 D u + q 2 D v   d X + p 1 D u + p 2 D v ,
where V 0 and A 0 are the referential volume and cross-sectional area, respectively. Besides, q 1 and q 2 are the components of external traction acting on the length, and p 1 and p 2 are applied point forces in X - and Z -directions. Equation (16) can be re-written as
D U int = L A D E T Π d A 0 d X = L A ( D E 0 T + Z D E 1 T ) ( Π 0 + Z Π 1 ) d A 0 d X                             = L A ( D E 0 T Π 0 + Z D E 0 T Π 1 + Z D E 1 T Π 0 + Z 2 D E 1 T Π 1 ) d A 0 d X ,
The approximation of the generalized displacement field is defined as
Q = u , v , φ , κ T = I = 1 N N H I ξ u I , v I , φ I , κ I T ,
in which Q is displacement vector, N N is total nodes per element, and H I are the interpolation functions. In this paper, three-node element is used. Therefore, the interpolation functions are defined as follows:
H 1 = 1 2 ξ 1 ξ ,         H 2 = 1 ξ 1 + ξ ,               H 3 = 1 2 ξ 1 + ξ .
Using the latter relation, the virtual Green-Lagrange strain vector can be expressed via the following relation
D E r Δ E r = I = 1 N N B r I D Q Δ Q ,         r = 0 , 1
Substituting Equations (19) and (21) into (18) one can obtain:
D U int = 𝒜 e = 1 n e D Q T f int = 𝒜 e = 1 n e I = 1 N N D Q I T f I int ,    
in which 𝒜 e = 1 n e is the operator for assemblage and n e is total number of elements. Moreover, f I int is the vector of internal force expressed as
f I int = L ( B 0 I T F + B 1 I T m )   d X .
Similarly, the variation of the external work is given by
D U ext = 𝒜 e = 1 n e D Q T f ext = 𝒜 e = 1 n e I = 1 N N D Q I T f I ext ,
where f I ext is external force vector as follows:
f I ext = L H I q   d X + H I p 1 + H I p 2 .
By back-substitution of Equations (22) and (24) in the principle of virtual work i.e., Equation (15), one can obtain that:
𝒜 e = 1 n e D Q T ( f int f ext ) = 0 .
It is noted that, in the latter relation, D Q can take arbitrary values, and therefore, the residual vector is defined as
= f int f ext = 0 .
Equation (27) is a set of non-linear equations. Hence, an iterative Newton-Raphson method is employed to deal with it. To linearize the non-linear relation, one can write:
Δ D U int = L A 0 ( D E T Δ Π + Δ D E T Π ) d A 0 d X = L A 0 D E T Δ Π d A 0 d X + L A 0 Δ D E T Π d A 0 d X = Φ + Ω ,
where
Φ = L A 0 D E T Δ Π d A 0 d X ,           Ω = L A 0 Δ D E T Π d A 0 d X .
Therefore, the linearization of variation of internal stored energy is derived as
Δ D U int = Φ + Ω .
Finally, utilizing Equation (21), the Equation (30) can be re-written as
Δ D U int = 𝒜 e = 1 n e I = 1 N N J = 1 N N D Q I T K I J Δ Q J = D Q T K Δ Q ,
in which K I J is the stiffness matrix in the IJ’th block per element expressed via the following relation
K I J = K I J mat + K I J geo .
In the latter relation K I J mat is the material stiffness matrix given by
K I J mat = L { A B 0 I T 0 B 0 J + I ˜ ( B 0 I T 1 B 1 J + B 1 I T 0 B 1 J + B 1 I T 1 B 0 J ) } d X .
Moreover, K I J geo is the geometric stiffness matrix expressed as
K I J geo = L ( A Δ D E 0 T Π + I ˜ Δ δ E 1 T Π ) d X .

4. Demonstrative Examples

The kinematics and kinetics of the four-parameter Timoshenko beam are derived in the previous sections. Moreover, the finite element formulation for large deformation of Timoshenko beam is also derived. In this section, to evaluate the performance and ability of the formulations, some examples are provided. In all examples, three-node elements are utilized. Furthermore, two-node Gauss-Legendre reduced integration is employed to calculate all integrals to avoid shear locking phenomenon.

4.1. Analysis of a Cantilever Subjected to End Force

In this example, deflection of a clamped-free beam under a tip force is studied. As shown in the Figure 2a, a cantilever with length L is applied to the point force P at the free end. The load is so that the load parameter P L 2 / E I goes to 5, where E is the elasticity modulus and I is the moment of inertia. Our convergence analysis shows that a mesh of ten elements is enough to get converged results. The deformation-force curves are depicted in Figure 2a,b. It is deduced that, the result for vertical deflection ratio v / L and horizontal deflection ratio u / L in terms of load parameter P L 2 / E I are in
Very good agreement with the results reported in Wang et al. [33]. Moreover, the deformed shapes of the beam in ten load steps are displayed in Figure 3a. Besides, the convergence plot is shown in Figure 3b.

4.2. Analysis of a Clamped-Supported Beam under Distributed Load

In this example, finite elastic deformation of a clamped-supported beam is investigated. The beam is clamped in one end, and is supported in the other end so that it can move in x-direction. The beam is applied to uniform distributed load q (see Figure 4b). Our convergence study reveals that a mesh of 20 elements is enough to get converged results. The curves for large vertical deformations in X = L / 4 , X = L / 2 and X = 3 L / 4 in terms of load parameter p * = q L 2 / E I are shown in Figure 4. It can be deduced from Figure 4a that the deflection ratio v / L in X = 3 L / 4 is larger than X = L / 4 , and obviously, the deformation is not symmetric. Furthermore, the deformed shapes of the clamped-supported beam are portrayed in Figure 5a in ten load step. Finally, the convergence diagrams of deformation in the middle of the beam in terms of number of elements are depicted in Figure 5b.

4.3. Validation with Experiment

In the final example, verification of the present model with experimental results is performed. To do so, the experiment done by Beléndez et al. [34] is considered. According to their experiment, a cantilever of which dimensions are L = 40   cm , b = 2.5   cm and h = 0.04   cm are considered. The cantilever is under its own weight W = 0.3032   N which is equal to distributed load w = W / L = 0.758   N / m . Convergence analysis shows that a mesh of 5 elements is sufficient for solving this problem. The deformed shape in of the beam is depicted in Figure 6a. As can be observed, the final shape of the simulation in the present work is coincident with the final deformed shape of the experiment reported in [34]. This proves that the present model can predict the experimental results very good. Moreover, the curves of deformations in terms of loading are shown in Figure 6b.

5. Conclusions

In this paper, a non-linear finite element formulation for large elastic deflections of the Timoshenko beam is derived. The formulation is free of Poisson and shear locking, as it employs through-the-thickness stretch parameter and reduced integration. The demonstrative examples show very good agreement with the results reported in the literature. In Section 4.1, the present formulation can get the results for large deformations of a cantilever, and is coincident with the previous study in the literature. In Section 4.2, the deformation of a clamped-supported beam is investigated, and large deflections at three points of the beam are obtained. Besides, the formulations enable to capture the experimental results reported in the literature with completely convincing accuracy. In Section 4.3, the large defection of the beam under its own weight is investigated. A comparison to the experiment proves that the present model can predict the experiment very well.

Author Contributions

N.F. Conceptualization, Methodology, Software, Investigation, Writing—Original draft preparation. A.S.M.A. Conceptualization, Methodology, Software, Resources, Funding acquisition, Writing—Review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Taif University, Saudi Arabia, Project No. (TU-DSPP-2024-257).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors have declared no conflicts of interest.

Nomenclature

ParameterDefinition
L Length of the beam
b Width of the beam
h 0 Thickness of the beam
X , Y , Z Cartesian Coordinate system
e 1 , e 2 , e 3 Orthonormal basis vectors
X Position of each material point in undeformed beam
x Position of each material point in deformed beam
x 0 Position of centerline after deformation
κ Through-the-thickness stretch parameter
r Director vector
φ Rotation of the cross-section of the beam
F Deformation gradient tensor
F 0 , F 1 Linearized forms of deformation gradient tensor
u Displacement in X-direction
v Displacement in Z-direction
C Right Cauchy-Green deformation tensor
C 0 , C 1 Linearized parts of right Cauchy-Green deformation tensor
E Green-Lagrange strain tensor
E 0 , E 1 Linearized parts of the Green-Lagrange tensor
ψ Kirchhoff-St Venant strain energy density function
μ , λ Lame constants
E Young’s modulus
υ Poisson’s ratio
Π Second Piola-Kirchhoff stress tensor
Π 0 , Π 1 Linearized parts of the second Piola-Kirchhoff stress tensor
I Second-order identity tensor
Fourth-order elasticity tensor
0 , 1 Linearized form of fourth-order elasticity tensor
D U int Variation of external force
q 1 , q 2 Components of external traction
p 1 , p 2 Applied point forces
Q Generalized displacement vector
H I Interpolation function
𝒜 e = 1 n e Assemble operator
f I int Internal force
f I ext External force
Residual
K I J Total stiffness matrix
K I J mat Material stiffness matrix
K I J geo Geometric stiffness matrix
A Area of cross-section of the beam
I ˜ Moment of inertia of the beam

References

  1. Zozulya, V.; Saez, A. A high-order theory of a thermoelastic beams and its application to the MEMS/NEMS analysis and simulations. Arch. Appl. Mech. 2016, 86, 1255–1272. [Google Scholar] [CrossRef]
  2. Mabsout, M.E.; Tarhini, K.M.; Frederick, G.R.; Tayar, C. Finite-element analysis of steel girder highway bridges. J. Bridge Eng. 1997, 2, 83–87. [Google Scholar] [CrossRef]
  3. Alzaidi, A.S.M.; Kaplunov, J.; Prikazchikova, L. The edge bending wave on a plate reinforced by a beam. J. Acoust. Soc. Am. 2019, 146, 1061–1064. [Google Scholar] [CrossRef]
  4. Nefske, D.; Sung, S. Power flow finite element analysis of dynamic systems: Basic theory and application to beams. J. Vib. Acoust. Stress Reliab. 1989, 111, 94–100. [Google Scholar] [CrossRef]
  5. Lazarus, A.; Thomas, O.; Deü, J.-F. Finite element reduced order models for nonlinear vibrations of piezoelectric layered beams with applications to NEMS. Finite Elem. Anal. Des. 2012, 49, 35–51. [Google Scholar] [CrossRef]
  6. Wu, G.; He, X.; Pai, P.F. Geometrically exact 3D beam element for arbitrary large rigid-elastic deformation analysis of aerospace structures. Finite Elem. Anal. Des. 2011, 47, 402–412. [Google Scholar] [CrossRef]
  7. Alruthea, M.S. Biomechanical Analysis of the Zirconia and Graphene-based CAD-CAM Dental Bridges at Different Pontic Length: A Finite Element Analysis. J. Clin. Diagn. Res. 2020, 14, ZF01–ZF05. [Google Scholar] [CrossRef]
  8. Huang, L.-S.; Huang, Y.-C.; Yuan, C.; Ding, S.-J.; Yan, M. Biomechanical evaluation of bridge span with three implant abutment designs and two connectors for tooth-implant supported prosthesis: A finite element analysis. J. Dental Sci. 2023, 18, 248–263. [Google Scholar] [CrossRef] [PubMed]
  9. Capozucca, R. Static and dynamic response of damaged RC beams strengthened with NSM CFRP rods. Compos. Struct. 2009, 91, 237–248. [Google Scholar] [CrossRef]
  10. Lin, Y.-H.; Trethewey, M. Finite element analysis of elastic beams subjected to moving dynamic loads. J. Sound Vib. 1990, 136, 323–342. [Google Scholar] [CrossRef]
  11. Wang, X.; So, R. Timoshenko beam theory: A perspective based on the wave-mechanics approach. Wave Motion 2015, 57, 64–87. [Google Scholar] [CrossRef]
  12. Yu, W.; Hodges, D.H. Generalized Timoshenko theory of the variational asymptotic beam sectional analysis. J. Am. Helicopter Soc. 2005, 50, 46–55. [Google Scholar] [CrossRef]
  13. Mosavi, A.; Benkreif, R.; Varkonyi-Koczy, A.R. Comparison of Euler-Bernoulli and Timoshenko beam equations for railway system dynamics. In Recent Advances in Technology Research and Education: Proceedings of the 16th International Conference on Global Research and Education Inter-Academia 2017; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  14. Carrera, E.; Giunta, G.; Petrolo, M. Beam Structures: Classical and Advanced Theories; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  15. Bessa, M.A.; Elkhodary, K.I.; Liu, W.K.; Belytschko, T.; Moran, B. Nonlinear Finite Elements for Continua and Structures; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  16. Mata, P.; Oller, S.; Barbat, A. Dynamic analysis of beam structures considering geometric and constitutive nonlinearity. Comput. Methods Appl. Mech. Eng. 2008, 197, 857–878. [Google Scholar] [CrossRef]
  17. Khaniki, H.B.; Ghayesh, M.H.; Chin, R.; Amabili, M. Hyperelastic structures: A review on the mechanics and biomechanics. Int. J. Non-Linear Mech. 2023, 148, 104275. [Google Scholar] [CrossRef]
  18. Wagner, M.; Schaeffer, J. Rubbers and polymer melts: Universal aspects of nonlinear stress–strain relations. J. Rheol. 1993, 37, 643–661. [Google Scholar] [CrossRef]
  19. Read, B.; Dean, G. Modelling non-linear stress–strain behaviour of rubber toughened plastics. Plast. Rubber Compos. 2001, 30, 328–336. [Google Scholar] [CrossRef]
  20. Oñate, E. Structural Analysis with the Finite Element Method. Linear Statics: Volume 2: Beams, Plates and Shells; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  21. Litewka, P. Finite Element Analysis of Beam-to-Beam Contact; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010; Volume 53. [Google Scholar]
  22. Holzapfel, G.A. Nonlinear Solid Mechanics: A Continuum Approach for Engineering; John Wiley & Sons: Hoboken, NJ, USA, 2000. [Google Scholar]
  23. Meier, C.; Popp, A.; Wall, W.A. Geometrically exact finite element formulations for slender beams: Kirchhoff–love theory versus Simo–reissner theory. Arch. Comput. Methods Eng. 2019, 26, 163–243. [Google Scholar] [CrossRef]
  24. Cazzani, A.; Stochino, F.; Turco, E. An analytical assessment of finite element and isogeometric analyses of the whole spectrum of Timoshenko beams. Z. Angew. Math. Mech. 2016, 96, 1220–1244. [Google Scholar] [CrossRef]
  25. Sakai, Y.; Ohsaki, M.; Adriaenssens, S. A 3-dimensional elastic beam model for form-finding of bending-active gridshells. Int. J. Solids Struct. 2020, 193–194, 328–337. [Google Scholar] [CrossRef]
  26. Numanoğlu, H.M.; Ersoy, H.; Akgöz, B.; Civalek, Ö. A new eigenvalue problem solver for thermo-mechanical vibration of Timoshenko nanobeams by an innovative nonlocal finite element method. Math. Methods Appl. Sci. 2022, 45, 2592–2614. [Google Scholar] [CrossRef]
  27. Aria, A.I.; Rabczuk, T.; Friswell, M.I. A finite element model for the thermo-elastic analysis of functionally graded porous nanobeams. Eur. J. Mech. A/Solids 2019, 77, 103767. [Google Scholar] [CrossRef]
  28. Karamanli, A.; Vo, T.P.; Civalek, Ö. Finite element formulation of metal foam microbeams via modified strain gradient theory. Eng. Comput. 2023, 39, 751–772. [Google Scholar] [CrossRef]
  29. Civalek, Ö.; Uzun, B.; Yaylı, M.Ö. Finite element formulation for nano-scaled beam elements. ZAMM J. Appl. Math. Mech. 2022, 102, e202000377. [Google Scholar] [CrossRef]
  30. Aria, A.I.; Friswell, M.I.; Rabczuk, T. Thermal vibration analysis of cracked nanobeams embedded in an elastic matrix using finite element analysis. Compos. Struct. 2019, 212, 118–128. [Google Scholar] [CrossRef]
  31. Firouzi, N.; Lenci, S.; Amabili, M.; Rabczuk, T. Nonlinear free vibrations of Timoshenko–Ehrenfest beams using finite element analysis and direct scheme. Nonlinear Dyn. 2024, 112, 7199–7213. [Google Scholar] [CrossRef]
  32. Firouzi, N.; Rabczuk, T.; Bonet, J.; Żur, K.K. A computational framework for large strain electromechanics of electro-visco-hyperelastic beams. Comput. Methods Appl. Mech. Eng. 2024, 426, 116985. [Google Scholar] [CrossRef]
  33. Wang, J.; Chen, J.-K.; Liao, S. An explicit solution of the large deformation of a cantilever beam under point load at the free tip. J. Comput. Appl. Math. 2008, 212, 320–330. [Google Scholar] [CrossRef]
  34. Beléndez, T.; Neipp, C.; Beléndez, A. Numerical and experimental analysis of a cantilever beam: A laboratory project to introduce geometric nonlinearity in mechanics of materials. Int. J. Eng. Educ. 2003, 19, 885–892. [Google Scholar]
  35. Cai, R.; Jin, Y.; Pennec, Y.; Djafari-Rouhani, B.; Rabczuk, T.; Zhuang, X. Broadband non-reciprocal wave suppression and frequency conversion by active metabeams. Mech. Syst. Signal Process. 2024, 220, 111656. [Google Scholar] [CrossRef]
  36. Sahmani, S.; Safaei, B.; Rabczuk, T. Nonlinear in-plane buckling of small-curved and large-curved FG porous microbeams via strain gradient-based isogeometric collocation formulations. Compos. Struct. 2024, 334, 117969. [Google Scholar] [CrossRef]
  37. Uzun, B.; Civalek, Ö.; Yaylı, M.Ö. Nonlinear stability analysis of embedded restrained nanobeams using the Stokes’ transformation. Mech. Based Des. Struct. Mach. 2024, 52, 2504–2531. [Google Scholar] [CrossRef]
  38. Yee, K.; Ong, O.Z.S.; Ghayesh, M.H.; Amabili, M. Various homogenisation schemes for vibration characteristics of axially FG core multilayered microbeams with metal foam face layers based on third order shear deformation theory. Appl. Math. Model. 2024, 125, 189–217. [Google Scholar] [CrossRef]
  39. Settimi, V.; Lenci, S. Periodic wave propagation in nonlocal beams resting on a bilinear foundation. Theor. Appl. Mech. 2023, 50, 117–131. [Google Scholar] [CrossRef]
  40. Żur, K.K.; Firouzi, N.; Rabczuk, T.; Zhuang, X. Large deformation of hyperelastic modified Timoshenko–Ehrenfest beams under different types of loads. Comput. Methods Appl. Mech. Eng. 2023, 416, 116368. [Google Scholar] [CrossRef]
  41. Bensaid, I.; Saimi, A.; Civalek, Ö. Effect of two-dimensional material distribution on dynamic and buckling responses of graded ceramic-metal higher order beams with stretch effect. Mech. Adv. Mater. Struct. 2024, 31, 1760–1776. [Google Scholar] [CrossRef]
  42. Hameury, C.; Ferrari, G.; Buabdulla, A.; Silva, T.M.P.; Balasubramanian, P.; Franchini, G.; Amabili, M. Multiple-input multiple-output active vibration control of a composite sandwich beam by fractional order positive position feedback. Mech. Syst. Signal Process. 2023, 200, 110633. [Google Scholar] [CrossRef]
  43. Liu, X.; Pagani, A.; Carrera, E.; Liu, X. Free vibration analysis of composite beams and laminated reinforced panels by refined dynamic stiffness method and CUF-based component-wise theory. Compos. Struct. 2024, 337, 118058. [Google Scholar] [CrossRef]
  44. Pagani, A.; Chiaia, P.; Carrera, E. Vibration of solid and thin-walled slender structures made of soft materials by high-order beam finite elements. Int. J. Non-Linear Mech. 2024, 160, 104634. [Google Scholar] [CrossRef]
Figure 1. Schematic of the beam before and after deformation.
Figure 1. Schematic of the beam before and after deformation.
Symmetry 16 00984 g001
Figure 2. The curves of deflection ratio versus load parameter P L 2 / E I , (a) vertical deflection ratio v / L , (b) horizontal deflection ratio u / L [33].
Figure 2. The curves of deflection ratio versus load parameter P L 2 / E I , (a) vertical deflection ratio v / L , (b) horizontal deflection ratio u / L [33].
Symmetry 16 00984 g002
Figure 3. (a) Deformed shapes of the cantilever under end force, (b) convergence plot in terms of number of elements.
Figure 3. (a) Deformed shapes of the cantilever under end force, (b) convergence plot in terms of number of elements.
Symmetry 16 00984 g003
Figure 4. Curves of deformation ratio in terms of load parameter p * = q L 2 / E I at three different positions of the clamped-supported beam (a) v / L , (b) u / L .
Figure 4. Curves of deformation ratio in terms of load parameter p * = q L 2 / E I at three different positions of the clamped-supported beam (a) v / L , (b) u / L .
Symmetry 16 00984 g004
Figure 5. (a) Deformed shapes of the clamped-supported beam under uniform distributed loading in ten load steps, (b) convergence in the middle of the beam in terms of number of elements.
Figure 5. (a) Deformed shapes of the clamped-supported beam under uniform distributed loading in ten load steps, (b) convergence in the middle of the beam in terms of number of elements.
Symmetry 16 00984 g005
Figure 6. (a) The final deformed shape of the cantilever under weight and comparison with experiment, (b) the deformation-loading curves of the cantilever [34].
Figure 6. (a) The final deformed shape of the cantilever under weight and comparison with experiment, (b) the deformation-loading curves of the cantilever [34].
Symmetry 16 00984 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Firouzi, N.; Alzaidi, A.S.M. Non-Linear Elastic Beam Deformations with Four-Parameter Timoshenko Beam Element Considering Through-the-Thickness Stretch Parameter and Reduced Integration. Symmetry 2024, 16, 984. https://doi.org/10.3390/sym16080984

AMA Style

Firouzi N, Alzaidi ASM. Non-Linear Elastic Beam Deformations with Four-Parameter Timoshenko Beam Element Considering Through-the-Thickness Stretch Parameter and Reduced Integration. Symmetry. 2024; 16(8):984. https://doi.org/10.3390/sym16080984

Chicago/Turabian Style

Firouzi, Nasser, and Ahmed S. M. Alzaidi. 2024. "Non-Linear Elastic Beam Deformations with Four-Parameter Timoshenko Beam Element Considering Through-the-Thickness Stretch Parameter and Reduced Integration" Symmetry 16, no. 8: 984. https://doi.org/10.3390/sym16080984

APA Style

Firouzi, N., & Alzaidi, A. S. M. (2024). Non-Linear Elastic Beam Deformations with Four-Parameter Timoshenko Beam Element Considering Through-the-Thickness Stretch Parameter and Reduced Integration. Symmetry, 16(8), 984. https://doi.org/10.3390/sym16080984

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