Next Article in Journal
The Real-Time Optimal Attitude Control of Tunnel Boring Machine Based on Reinforcement Learning
Previous Article in Journal
Sensitivity Analysis of Structural Parameters of Unequal-Span Continuous Rigid Frame Bridge with Corrugated Steel Webs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Closed-Form Analytical Solutions for the Deflection of Elastic Beams in a Peridynamic Framework

1
Department of Mechanical Engineering, National Taiwan University, Taipei 10617, Taiwan
2
Institute of Mechanics, Otto-von-Guericke-University Magdeburg, 39106 Magdeburg, Germany
3
Department of Civil and Architectural Engineering, Harbin Far East Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(18), 10025; https://doi.org/10.3390/app131810025
Submission received: 18 August 2023 / Revised: 2 September 2023 / Accepted: 4 September 2023 / Published: 5 September 2023

Abstract

:
Peridynamics is a continuum theory that operates with non-local deformation measures as well as long-range internal force/moment interactions. The resulting equations are of the integral type, in contrast to the classical theory, which deals with differential equations. The aim of this paper is to analyze peridynamic governing equations for elastic beams. To this end, the strain energy density is formulated as a function of the non-local curvature. By applying the Lagrange principle, the peridynamic equations of motion are derived. Examples of non-local boundary conditions, including simple support, clamped edge, roller clamped edge, and free edge, are presented by introducing the interaction domain. Novel closed-form analytical solutions to integral equations are presented for beams with various boundary conditions, including clamped—simply supported, clamped–clamped, simply supported–roller-clamped, and clamped–roller-clamped beams. Furthermore, different types of loadings, including uniformly distributed load, concentrated force, and concentrated moment, are considered. The results are validated by comparing the derived solutions against solutions to the classical Bernoulli–Euler beam theory. A very good agreement between the non-local and the classical theories is observed for the case of the small horizon sizes, which shows the capability of the derived equations of motion and proposed boundary conditions.

1. Introduction

In recent years, non-local theories of solid mechanics have become of increasing interest, due to many applications such as analysis of micro- and nanometer-scale components, laminates with thin layers, composite materials with extreme differences in properties of constituents, among others.
One way to develop a non-local theory is to extend the classical continuum mechanics (CCM) by introducing higher-order strain gradients as arguments of the strain energy density [1]. By applying a variational principle, the governing equations of motion and boundary conditions can be derived. For linear-elastic material behavior, and in the case of statics, a strain gradient theory of the nth order provides a 2nth-order system of linear partial differential equations, e.g., [2]. In addition, higher-order stress tensors are introduced, and many boundary conditions are required to pose the boundary-value problem. Strain gradient theories were recently applied to the analysis of micro- and nanometer structures in, e.g., [3,4] (see also the works cited therein).
Alternatively, non-local continuum theories can be formulated within the peridynamics (PD) framework [5,6]. Here, the balance equations for long-range internal forces are introduced. The deformation gradient, higher strain gradients, or gradients of internal state variables are not required. Peridynamic theories for three-dimensional solids are discussed in [5].
Thin-walled structural components can be modeled by using classical three-dimensional continuum mechanics or alternatively by applying classical theories of beams, plates, and shells, as discussed in many textbooks, e.g., [7]. Within the framework of peridynamics, a three-dimensional theory can be used in the analysis of thin-walled structures, e.g., [8]. However, a very high computational effort is required for the numerical solution, as the thickness direction must be discretized by many node layers to to provide satisfactory accuracy. Furthermore, for components with a very low thickness-to-span-length ratio, for example, soft polymeric layers in laminated glass [9], the PD horizon size may be comparable to or even larger than the layer thickness. This would lead to non-physical long-range force interactions between the layers in a laminate. Therefore, robust peridynamic theories for beams, plates, and shells are required to analyze slender components.
Similarly to classical theories, PD theories for thin-walled structures can be derived from balance laws for a deformable line in the case of a beam, a deformable plane for a plate, and a deformable surface for a shell. These types of theories can be referred to as direct theories by analogy to the classical continuum mechanics. A direct PD beam theory proposed in [10] assumes the strain energy density to be a function of a non-local curvature. Discretized PD equations of motion with respect to the displacements of a straight beam are derived in [11]. For initially curved beams, PD equations of motion are presented in [12]. Discretized forms of PD equations of motion for plates and shells are derived in [13] based on the Lagrange variational principle. Peridynamic theories for membranes and fibers are presented in [14]. A direct PD shell theory in terms of displacements and linearized rotations is proposed in [15]. In [16], a non-linear direct PD plate theory is developed. Points of the plate midplane are assumed to have the rigid body-type degrees of freedom: displacements and finite rotations. Long-range force and moment interactions between the points of a plate are introduced. Governing balance laws are formulated directly for the plate.
Alternatively, applying the through-the-thickness approximations for the kinematic variables, a three-dimensional PD theory can be reduced to one- or two-dimensional cases. In [17], a PD plate theory is developed from the three-dimensional bond-based PD theory by applying asymptotic analysis. A peridynamic beam theory proposed in [18] is derived from three-dimensional PD equations and cross-section approximations for the deformation vector.
The governing peridynamic equations of motion are of integro-differential type, unlike CCM and strain gradient theories, where partial differential equations are introduced. To solve PD equations numerically, various approaches such as the mesh-free method [19] and discontinuous finite element formulation [20] have been proposed. Various numerical solutions are presented to demonstrate the ability of PD to represent complex phenomena such as crack initiation [8], delamination [21], and propagation of cracks [22,23].
However, only several closed-form analytical solutions to PD integral equations are derived. Examples include solutions for straight axially loaded rods [24,25,26]. For initially straight beams with specific boundary constraints, series solutions to PD equations are presented in [27]. Series solutions for rectangular plates under special boundary conditions are given in [28]. General solutions to peristatics and peridynamics for three-dimensional media are presented in [29,30].
Analytical solutions explicitly illustrate the effect of non-locality on the deformation state. In particular, the PD horizon size enters the results for the beam deformations and can be varied to analyze the convergence to the classical solution. This gives an advantage over numerical methods, where the horizon size is related to the size of the spatial discretization and can be only varied in a small range. Furthermore, different types of constraints can be introduced and analyzed, as they affect the final results, by applying the analytical solutions. Furthermore, closed-form analytical solutions are useful as reference solutions in order to validate the numerical methods developed to solve PD equations, such as the mesh-free method, and to analyze the convergence of numerical results.
The aim of this study is to derive a family of novel closed-form analytical solutions to both peristatic and peridynamic equations of shear-rigid beams. In addition to results presented in [27] for simply supported cantilever beams, various boundary conditions, including clamped–simply supported, clamped–clamped, simply supported–roller-clamped, and clamped–roller-clamped beams, will be considered. Furthermore, novel closed-form solutions for different types of loadings, including uniformly distributed load, concentrated force, and concentrated moment, will be introduced. The results will be validated by comparing the derived solutions against available solutions to the classical Bernoulli–Euler beam theory (CBT).

2. Peridynamic Beam Theory

2.1. Classical Equations and Assumptions

Let x be the coordinate for a point of the centerline of the initially straight beam, and z is the transverse coordinate of cross-section points with the origin at the centroid of the cross-section. The corresponding Cartesian base vectors will be specified by e x and e z , respectively. In what follows, let us analyze the plane-bending of the beam. In this case, the deformed centerline belongs to the x z plane. In the case of small cross-sectional rotations, the displacement vector u can be represented as follows:
u ( x , z ) = w ( x ) e z + θ ( x ) z e x ,
where w is the deflection and θ is the small angle of cross-sectional rotation.
Assuming small strains, the linear strain tensor ε is related to the displacement field according to the following equation:
ε = 1 2 u + ( u ) T , ( ) = e x ( ) x + e z ( ) z .
With Equation (1) inserted into Equation (2), the strain tensor is computed as follows:
ε = ε e x e x + 1 2 γ ( e x e z + e z e x ) ,
with the following expressions for normal strain ε and the transverse shear strain γ :
ε = z d θ d x , γ = θ + d w d x .
For the shear rigid beams, γ = 0 is assumed, resulting in
θ = d w d x , ε = z d 2 w d x 2 .
For the linear-elastic material behavior, the strain energy per unit length is defined as follows:
U = 1 2 A E ε 2 d A ,
where E is the Young’s modulus, and A is the cross-sectional area. With Equation (5), the strain energy density (6) takes the following form:
U = E I 2 κ 2 , κ = d 2 w d x 2
where I is the second area moment of the beam cross-section.

2.2. Peridynamic Equations of Motion

Unlike the CCM and the classical beam theory (CBT), where strains are defined in terms of derivatives with respect to the coordinates, non-local deformation measures are introduced in the PD framework. The PD deformations can be introduced as work conjugates to the bond force and moment densities considering the energy balance equation, e.g., [5,16]. Assuming the linear-elastic material behavior, a class of PD theories for beams and plates can be derived with the correspondence of the strain energy density from the CCM as a function of the PD deformation measures [27,28].
Let Γ be the collection of points of the beam centerline. Introduce the PD horizon as an open ball centered at point x with the PD horizon size δ defined as x Γ , i.e.,
H x : = x Γ : x x < δ .
In PD, the strain energy density of any arbitrary material point x depends upon both the displacement of itself and the displacement fields within its horizon H x , i.e.,
U ¯ ( x ) = U ¯ w ( H x ) .
To formulate the PD strain energy density, consider the following non-local curvature:
κ ¯ ( x ) = 1 δ x δ x + δ w ( x ) w ( x ) ( x x ) 2 d x .
For smooth deflection function w ( x ) , it can be shown that κ ¯ ( x ) approaches the local curvature (7) 2 as δ 0 . With Equations (7) 1 and (10), the PD bending strain energy density is formulated as follows:
U ¯ ( x ) = 1 2 E I δ 2 x δ x + δ w ( x ) w ( x ) ( x x ) 2 d x 2 .
For elastic beams, the PD equation of motion can be derived using Lagrange’s equation,
d d t L w ˙ L w = 0 ,
where L = T U ¯ + W is the Lagrangian. The kinetic energy T is computed as follows:
T = 1 2 0 L ρ A w ˙ 2 d x
in which ρ is the density and L is the length of the beam. The strain energy of the beam U ¯ and the work of external forces W are
U ¯ = 0 L U ¯ ( x ) d x , W = 0 L p ( x ) w ( x ) d x ,
where p ( x ) is the distributed lateral load.
With Equations (11)–(14), the following equation of motion is derived:
ρ A w ¨ ( x ) = E I δ 2 x δ x + δ 1 x x 2 x δ x + δ w ( ζ ) w ( x ) ζ x 2 d ζ x δ x + δ w ( ζ ) w ( x ) ζ x 2 d ζ d x + p ( x ) .
By changing the integration variables, Equation (15) can be put in the following form:
ρ A w ¨ ( x ) = E I δ 2 δ δ 1 ξ 2 δ δ w x + η w x η 2 d η δ δ w x + ξ + η w x + ξ η 2 d η d ξ + p x .
For the static case, the PD governing equation for the Bernoulli–Euler beam (16) can be expressed as follows:
E I δ 2 δ δ 1 ξ 2 δ δ w x + η w x η 2 d η δ δ w x + ξ + η w x + ξ η 2 d η d ξ = p x .

2.3. Peridynamic Boundary Conditions

The governing equation, Equation (16), is valid only if the material points are completely embedded in its PD influence domain. For material points adjacent to the boundary, it is necessary to introduce a fictitious interaction material region. For the introduced beam theory, the width of the fictitious region is double the PD horizon size, 2 δ , as shown in Figure 1.

2.3.1. Simple Support

Suppose the beam is subjected to the simply supported boundary condition on the left. In this case, the deflection and curvature must vanish, i.e.,
w 0 = 0 , d 2 w d x 2 x = 0 = 0 .
Applying the central difference for Equation (18) 2 yields
d 2 w d x 2 x = 0 w Δ x 2 w 0 + w Δ x Δ x 2 = 0 .
Replacing Δ x with ξ and associating it with Equation (18) 1 gives the PD simply supported condition as
w ξ = w ξ , for 0 ξ 2 δ .
It indicates that the simply supported end enforces in the fictitious region an anti-symmetric relation with respect to the real material domain, as shown in Figure 2.

2.3.2. Roller Clamped Edge

In the case of a roller-clamped edge, Figure 3, a rotational rigidity is required such that
d w d x x = 0 = 0 .
With the central difference, the PD boundary condition for the roller-clamped boundary condition takes the form
d w d x x = 0 w ξ w ξ 2 ξ = 0 w ξ = w ξ .
One can observe that roller-clamped boundary condition in PD enforces a mirror relation between the fictitious region and the real region.

3. Analytical Solutions to Peristatic Problems

3.1. Simply Supported Beam

Consider a simply supported beam subjected to a transverse load p ( x ) , as shown in Figure 4. In this case, the PD boundary conditions can be given as follows:
w ξ = w ξ , w L + ξ = w L ξ , 0 < ξ 2 δ .
With the BCs (23), the deflection field can be extended periodically with the period of 2 L , as shown in Figure 5.
Therefore, the deflection function can be expanded by applying the following Fourier series:
w x = n = 1 b n sin n π x L .
Substituting Equation (24) into (17) provides, after some simplifications,
E I δ 2 n = 1 b n sin n π x L δ δ 1 ξ 2 1 cos n π ξ L d ξ , 2 = p x
where the unknown coefficients can be determined by invoking the orthogonality of trigonometric functions:
b n = 2 δ 2 L E I 0 L p x sin n π x L d x δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 .
Plugging Equation (26) back into (24) gives the following explicit PD solution for the simply supported beam:
w x = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L .

3.2. Clamped–Simply Supported Beam

For the clamped–simply supported beam, Figure 6a shows the solution that can be applied. The rotational constraint of the clamped end can be released by introducing the moment M 0 , as shown in Figure 6b. By applying Castigliano’s theorem, moment M 0 can be computed as follows:
M 0 = 1 L 0 L 0 x x a p a d a d x 3 2 L 2 0 L 0 x x a p L x d a x d x + L 2 0 L p x d x
In the case of the shear rigid beam, the bending moment can be transformed into body force associated with a derivative of the Dirac delta. Indeed, a moment can be replaced with the force couple, as shown in Figure 7.
For such a case, the external load term can be formulated as follows:
p x = lim Δ x 0 F δ x x 0 Δ x 2 δ x x 0 + Δ x 2 = lim Δ x 0 M δ x x 0 + Δ x 2 δ x x 0 Δ x 2 Δ x = M δ x x 0 .
With Equations (27)–(29), the solution to the clamped–simply supported beam can be formulated as follows:
w x = 2 δ 2 L E I n = 1 0 L p x + M 0 δ x 0 sin n π x L d x δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x + 0 L M 0 δ x 0 sin n π x L d x δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x + M 0 δ x 0 sin n π x L d x 0 L 0 L n π L δ x 0 cos n π x L d x × δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x M 0 n π L δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L .

3.3. Clamped–Clamped Beam

The clamped–clamped beam, Figure 8a, can be considered as a simply supported beam subjected to bending moments M 0 and M L at the boundaries, as shown in Figure 8b, if the two rotational constrains are released. By applying Castigliano’s theorem, we obtain
M 0 = 6 L 2 0 L 0 x p a x a d a x d x 4 L 0 L 0 x p a x a d a d x M L = 0 L p x L x d x + 6 L 2 0 L 0 x p a x a d a x d x 2 L 0 L 0 x p a x a d a d x .
With the PD solution for simply supported beam (27) and (31), the PD solution for the clamped–clamped beam can be cast as
w ( x ) = 2 δ 2 L E I n = 1 0 L p x + M 0 δ x 0 + M L δ x L sin n π x L d x × × δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x + 0 L M 0 δ x 0 sin n π x L d x + 0 L M L δ x L sin n π x L d x × δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x + M 0 0 L n π L δ x 0 cos n π x L d x + M L 0 L n π L δ x L cos n π x L d x δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L = 2 δ 2 L E I n = 1 0 L p x sin n π x L d x n π L M 0 + M L 1 n δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L .

3.4. Simply Supported–Roller-Clamped Beam

A simply supported–roller-supported beam subjected to transverse loading p ( x ) is illustrated in Figure 9. The PD boundary conditions in this case admit an anti-symmetric and symmetric relation between the fictitious and real domains
w ( ξ ) = w ( ξ ) , w ( L ξ ) = w ( L + ξ ) .
Based on relation (33), the displacement field can be extended over the entire real line such that it has a periodicity of 4 L and is anti-symmetric about the origin and symmetric about x = L , as illustrated in Figure 10.
Therefore, the displacement can be expanded using the Fourier series as follows:
w x = w x w L x = w L + x w x = n = 1 b n sin 2 n 1 π x 2 L .
Following the similar approach given in Section 3.1, the PD solution to the simply supported–roller-supported beam can be expressed explicitly as follows:
w x = 2 δ 2 L E I n = 1 0 L p x sin 2 n 1 π x 2 L d x δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 sin 2 n 1 π x 2 L .

3.5. Clamped–Roller-Clamped Beam

Figure 11a illustrates a clamped–roller-clamped beam. To derive the analytical solution, the rotational constraint in the left clamped edge is removed, and the moment M 0 is introduced, as shown in Figure 11b. By applying Castigliano’s theorem, we obtain
M 0 = 1 L 0 L 0 x p a x a d a d x + L 2 0 L p x d x .
Again, the boundary moment M 0 can be absorbed as a part of the body force, and hence, based on the solution form of simply supported–roller-clamped beam, the solution for the clamped–roller-clamped beam can be given as
w x = 2 δ 2 L E I n = 1 0 L p x sin 2 n 1 π x 2 L d x M 0 2 n 1 π 2 L × δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 sin 2 n 1 π x 2 L .

3.6. Cantilever Beam

Consider a cantilever beam subjected to the transverse load p ( x ) . The solution can be identically decomposed using the following steps (see Figure 12).
  • Place a fictitious simple support at the free end.
  • Transit the origin to be level with the free end and then subtract the corresponding rigid body motion.
  • Release the vertical displacement constraint and replace it with the corresponding shear force Q.
The shear force Q must balance with the resultant of the external loading such that
Q = 0 L p x d x .
By modifying the PD solution for the simply supported–roller-clamped beam (35) with consideration of shear force Q at the left end, the explicit PD solution for the cantilever beam can be given as follows:
w x = 2 δ 2 L E I n = 1 0 L p x cos 2 n 1 π x 2 L d x + Q δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 cos 2 n 1 π x 2 L 1 = 2 δ 2 L E I n = 1 0 L p x cos 2 n 1 π x 2 L 1 d x δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 cos 2 n 1 π x 2 L 1 .

4. Analytical Solutions to Peridynamic Problems

Let us recall the PD equation of motion (16):
ρ A w ¨ ( x ) = E I δ 2 δ δ 1 ξ 2 δ δ w x + ζ w x ζ 2 d ζ δ δ w x + ξ + η w x + ξ η 2 d η d ξ + p x .
To formulate the analytical solutions to Equation (40), we apply the Laplace transform as follows:
ρ A s 2 W x , s s w 0 x v 0 x = c δ δ 1 ξ 2 δ δ W x + ζ , s W x , s ζ 2 d ζ δ δ W x + ξ + η , s W x + ξ , s η 2 d η . d ξ + P x , s ,
where
W x , s = L w x , t , P x , s = L p x , t , w 0 x = w x , 0 , v 0 x = w ˙ x , 0 .
Consider a simply supported beam subjected to some arbitrary dynamic load. By analogy to Section 2.3.1, we can assume that the solution to Equation (41) has the following form:
W x , s = n = 1 B n s sin n π x L .
After plugging Equation (43) into Equation (41), we obtain
n = 1 B n s c δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 + ρ A s 2 sin n π x L = P x , s + s ρ A w 0 x + ρ A v 0 x .
According to orthogonality, we obtain
B n s = 2 L 0 L P x , s + s ρ A w 0 x + ρ A v 0 x sin n π x L d x ρ A s 2 + c δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 = 2 L 1 ρ A 0 L P x , s sin n π x L d x + s w ^ n + v ^ n s i ω n s + i ω n ,
where
w ^ n = 0 L w 0 x sin n π x L d x , v ^ n = 0 L v 0 x sin n π x L d x , ω n = c ρ A δ δ 1 ξ 2 1 cos n π ξ L d ξ .
Applying partial fractions to Equation (45) it follows that
B n s = 2 L 1 ρ A 0 L P x , s sin n π x L d x s i ω n s + i ω n + i ω n w ^ n + v ^ n 2 i ω n 1 s i ω n + i ω n w ^ n v ^ n 2 i ω n 1 s + i ω n .
After performing the inverse Laplace transform on (47), we obtain
b n t = 2 L 1 ρ A 1 2 i ω n 0 t e i ω n u e i ω n u 0 L b x , t u sin n ^ x d x d u + i ω n w ^ n + v ^ n 2 i ω n e i ω n t + i ω n w ^ n v ^ n 2 i ω n e i ω n t = 2 L w ^ n cos ω n t + v ^ n ω n sin ω n t + 1 ρ A 1 ω n 0 t sin ω n u 0 L p x , t u sin n π x L d x d u .
In this case, the solution is
w x , t = L 1 W x , s = n = 1 b n t sin n π x L .
Coupling (48) with (49) yields the complete solution for a simply supported beam as follows:
w x , t = 2 L n = 1 w ^ n cos ω n t + v ^ n ω n sin ω n t + 1 ρ A 1 ω n 0 t sin ω n u 0 L p x , t u sin n π x L d x d u sin n π x L
where
w ^ n = 0 L w 0 x sin n π x L d x , v ^ n = 0 L v o x sin n π x L d x , ω n = c ρ A δ δ 1 ξ 2 1 cos n π ξ L d ξ .
Regarding the dynamic response of the clamped–simply supported and clamped–clamped case, we will have to add moments on the boundaries accordingly and absorb them into the body force as follows:
p ¯ x , t = M 0 t δ x + M L t δ x L + p x , t .
Here, M 0 t and M L t must balance the net effect of external loading and inertial force, and it can be computed using Castigliano’s theorem. Regarding simply supported–roller-clamped case and its ramifications, we can replace the term n π L by 2 n 1 π 2 L and follow a similar procedure.

5. Examples

5.1. Peristatic Solutions

The derived closed-form solutions can be applied to a large class of beam problems for various boundary conditions and loading cases, including distributed and concentrated loads. Below, we present plots of derived solutions for the deflection function for several examples. In addition, well-known analytical solutions to the classical beam theory will be plotted in order to demonstrate the validation of both the PD beam theory and the derived results.
Consider a beam with a rectangular cross-section with high h = 0.01 m and width b = 0.01 m. The length of the beam is assumed to be L = 1 m, and Young’s modulus of the beam material is E = 200 GPa. conditions. The PD horizon size is chosen as δ = 0.001 m. Figure 13 illustrates the plots of closed-form series solutions developed in this study. Different cases of loadings, including concentrated moment applied at different points of the beam, uniformly distributed load over the beam length, concentrated force at different points of the beam, and the load distributed over a part of the beam are considered. Furthermore, different examples of boundary conditions, including simply supported–simply supported, clamped–simply supported, clamped–clamped, simply supported–roller-clamped, clamped–roller-clamped, and clamped–free, are analyzed. For comparison, solutions according to the classical beam theory for the same loadings and boundary conditions are presented. A very good agreement between the results for the selected value of the horizon size can be observed.
Let us note that the horizon size is an important parameter in PD, providing a length-scale influence and the non-locality of the system. If the horizon size tends to zero, the PD solution converges to that of the classical beam theory in the case of elastic material behavior. The horizon size enters the derived analytical solutions in the explicit form. Plots of deflection vs. horizon size for a simply supported beam are presented in [27].
Let us summarize the closed-form analytical solutions for the six cases presented in Figure 13 as follows:
  • Simply supported–simply supported
    w x = 2 c L n = 1 n π L M cos n π 2 δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L
    based on PD and
    w x = M 24 L E I 3 L 3 + 12 L 2 x 12 L x 2 H 2 x L L 2 x + 4 x 3
    based on CBT.
  • Clamped–simply supported
    w x = 2 c L n = 1 0 L p sin n π x L d x n π L p 8 L 2 δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L
    based on PD and
    w x = p x 2 48 E I 3 L 2 5 L x + 2 x 2
    based on CBT.
  • Clamped–clamped
    w x = 2 p c L n = 1 sin n π a L 4 27 n π + n π 1 n 2 27 δ δ 1 ξ 2 1 cos n π ξ L d ξ 2 sin n π x L
    based on PD and
    w x = F 162 E I 4 x 2 3 L 5 x L 3 x 3 H 3 x L
    based on CBT.
  • Simply supported–roller-clamped
    w x = 2 p c L n = 1 L 4 3 L 4 sin 2 n 1 π x 2 L d ξ δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 sin 2 n 1 π x 2 L
    based on PD and
    w x = p 6144 E I L 4 16 L 3 x + 96 L 2 x 2 256 L x 3 + 256 x 4 H L 4 x + 81 L 4 + 432 L 3 x 864 L 2 x 2 256 x 4 + 768 L x 3 H 3 L 4 x L 3 L 1 6 x + 27 L 3 3 L 16 x + 768 L 2 x 2
    based on CBT.
  • Clamped–roller-clamped
    w x = 2 p c L n = 1 sin 2 n 1 π 4 3 16 2 n 1 π δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 sin 2 n 1 π x 2 L
    based on PD and
    w x = F 48 E I x 2 9 L 8 x L 2 x 3 H 2 x L
    based on CBT
  • Cantilever beam
    w x = 2 c L n = 1 M 2 n 1 π 2 L 1 n + 1 δ δ 1 ξ 2 1 cos 2 n 1 π ξ 2 L d ξ 2 1 cos 2 n 1 π x 2 L
    based on PD and
    w x = M 2 E I x 2
    based on CBT

5.2. Peridynamic Solutions

Consider the simply supported beam subjected to a harmonic load, as shown in Figure 14.
Let us set L = 2 m, A = 0.01 × 0.01 m 2 , E = 200 × 10 9 Pa, ρ = 7850 kg/m 3 , δ = 0.01 m. The initial conditions are assumed as follows: u x , 0 = 0.01 x x L and u ˙ x , 0 = sin π x L . The harmonic concentrated load is
p x , t = 50 sin ω 1 3 t δ x L 2 ,
where δ is the Dirac delta function. Figure 15 illustrates the transverse displacement vs. time plots for two different points along the beam axis.
We observe that for the considered value of the PD horizon, the results according to the PD beam theory and the classical beam theory almost coincide. As the next example, consider a mass drop on the midpoint of a simply supported beam; see Figure 16.
As in the previous example, we set L = 2 m, A = 0.01 × 0.01 m 2 , E = 200 × 10 9 Pa, ρ = 7850 kg/m 3 , and δ = 0.01 m. The balance of momentum for the mass yields
m v t 1 + m v t 1 = m v t 1 + + v t 1 = m v ¯ , v ¯ = v t 1 + + v t 1 .
Specifying the concentrated force acting on the beam by F and considering the interaction over Δ t , we obtain
m v ¯ H t t 1 H t t 1 + Δ t = F Δ t ,
where H is the Heaviside function. Consequently, the force is obtained as follows
F = m v ¯ H t t 1 H t t 1 + Δ t Δ t
or
F ( t ) = m v ¯ δ t t 1
as Δ t 0 . The initial conditions are u x , 0 = 0.01 x x L and u ˙ x , 0 = sin π x L The load is assumed as follows:
p x , t = F ( t ) δ x L 2 , F ( t ) = 10 δ t 0.5 .
Plots of the transverse displacement vs. time for two different points along the beam axis are shown in Figure 17.
Again, the results according to the PD beam theory and the classical beam theory agree if the horizon size is considered small.

6. Conclusions

In this paper, straight elastic beams are analyzed within the framework of peridynamics. Closed-form analytical solutions to both peristatic and peridynamic equations are presented for various examples with different types of loadings and boundary conditions. In addition to the cases analyzed in the literature, novel solutions for the clamped edge, roller-clamped edge, free edge, concentrated force, and concentrated moment are presented. The validation of PD formulations and solutions is accomplished by comparing the deflections against the corresponding results according to the Bernoulli–Euler beam theory. A very good agreement between the non-local and the classical theories is observed for the case of the small horizon sizes, which shows the capability of the derived PD equations and proposed boundary conditions.
The derived closed-form solutions include the horizon size in the explicit form such that the non-locality of the PD theory can be analyzed in a natural way. Furthermore, the results are useful to validate the numerical methods and numerical codes developed to solve PD equations.
It is worth noting that, unlike the classical beam theory, many different types of constraints can be applied within the peridynamic framework. Indeed, the deflection function can be specified over the interaction domain in many different ways, in addition to the cases discussed in this paper. Future studies should be related to aalyzing different types of PD constraints as they influence the deformed shape of the beam and deriving the corresponding analytical solutions.

Author Contributions

Conceptualization, Z.Y. and K.N.; methodology, Z.Y. and K.N.; writing—original draft preparation, Z.Y. and K.N.; writing—review and editing, Z.Y. visualization, Z.Y. and Y.C.; supervision, C.-C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bertram, A. Compendium on Gradient Materials; Springer: Cham, Switzerland, 2023. [Google Scholar] [CrossRef]
  2. Eremeyev, V.A. Strong Ellipticity and Infinitesimal Stability within Nth-Order Gradient Elasticity. Mathematics 2023, 11, 1024. [Google Scholar] [CrossRef]
  3. Cordero, N.M.; Forest, S.; Busso, E.P. Second strain gradient elasticity of nano-objects. J. Mech. Phys. Solids 2016, 97, 92–124. [Google Scholar] [CrossRef]
  4. Wang, B.; Zhao, J.; Zhou, S. A micro scale Timoshenko beam model based on strain gradient elasticity theory. Eur. J. Mech. A/Solids 2010, 29, 591–599. [Google Scholar] [CrossRef]
  5. Silling, S.A.; Lehoucq, R.B. Peridynamic theory of solid mechanics. Adv. Appl. Mech. 2010, 44, 73–168. [Google Scholar] [CrossRef]
  6. Silling, S.A. Introduction to peridynamics. In Handbook of Peridynamic Modeling; CRC Press: Boca Raton, FL, USA, 2016; pp. 63–98. [Google Scholar]
  7. Timoshenko, S.P.; Woinowsky-Krieger, S. Theory of Plates and Shells; McGraw-Hill: New York, NY, USA, 1959. [Google Scholar]
  8. Naumenko, K.; Pander, M.; Würkner, M. Damage patterns in float glass plates: Experiments and peridynamics analysis. Theor. Appl. Fract. Mech. 2022, 118, 103264. [Google Scholar] [CrossRef]
  9. Aßmus, M.; Naumenko, K.; Altenbach, H. A multiscale projection approach for the coupled global–local structural analysis of photovoltaic modules. Compos. Struct. 2016, 158, 340–358. [Google Scholar] [CrossRef]
  10. Chen, J. Nonlocal Euler-Bernoulli Beam Theories: A Comparative Study; Springer: Cham, Switzerland, 2021. [Google Scholar]
  11. Shen, G.; Xia, Y.; Hu, P.; Zheng, G. Construction of peridynamic beam and shell models on the basis of the micro-beam bond obtained via interpolation method. Eur. J. Mech. A/Solids 2021, 86, 104174. [Google Scholar] [CrossRef]
  12. Yang, Z.; Naumenko, K.; Ma, C.C.; Oterkus, E.; Oterkus, S. Peridynamic analysis of curved elastic beams. Eur. J. Mech. A/Solids 2023, 101, 105075. [Google Scholar] [CrossRef]
  13. Vazic, B.; Oterkus, E.; Oterkus, S. Peridynamic model for a Mindlin plate resting on a Winkler elastic foundation. J. Peridyn. Nonlocal Model. 2020, 2, 229–242. [Google Scholar] [CrossRef]
  14. Silling, S.A.; Bobaru, F. Peridynamic modeling of membranes and fibers. Int. J. Non-Linear Mech. 2005, 40, 395–409. [Google Scholar] [CrossRef]
  15. Chowdhury, S.R.; Roy, P.; Roy, D.; Reddy, J. A peridynamic theory for linear elastic shells. Int. J. Solids Struct. 2016, 84, 110–132. [Google Scholar] [CrossRef]
  16. Naumenko, K.; Eremeyev, V.A. A non-linear direct peridynamics plate theory. Compos. Struct. 2022, 279, 114728. [Google Scholar] [CrossRef]
  17. Taylor, M.; Steigmann, D.J. A two-dimensional peridynamic model for thin plates. Math. Mech. Solids 2015, 20, 998–1010. [Google Scholar] [CrossRef]
  18. Roy, P.; Behera, D.; Madenci, E. Peridynamic modeling of elastic instability and failure in lattice beam structures. Comput. Methods Appl. Mech. Eng. 2023, 415, 116210. [Google Scholar] [CrossRef]
  19. Silling, S.A.; Askari, E. A meshfree method based on the peridynamic model of solid mechanics. Comput. Struct. 2005, 83, 1526–1535. [Google Scholar] [CrossRef]
  20. Ren, B.; Wu, C.; Askari, E. A 3D discontinuous Galerkin finite element method with the bond-based peridynamics model for dynamic brittle failure analysis. Int. J. Impact Eng. 2017, 99, 14–25. [Google Scholar] [CrossRef]
  21. Dorduncu, M. Peridynamic modeling of delaminations in laminated composite beams using refined zigzag theory. Theor. Appl. Fract. Mech. 2021, 112, 102832. [Google Scholar] [CrossRef]
  22. Lu, W.; Oterkus, S.; Oterkus, E.; Zhang, D. Modelling of cracks with frictional contact based on peridynamics. Theor. Appl. Fract. Mech. 2021, 116, 103082. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Deng, H.; Deng, J.; Liu, C.; Yu, S. Peridynamic simulation of crack propagation of non-homogeneous brittle rock-like materials. Theor. Appl. Fract. Mech. 2020, 106, 102438. [Google Scholar] [CrossRef]
  24. Mikata, Y. Analytical solutions of peristatic and peridynamic problems for a 1D infinite rod. Int. J. Solids Struct. 2012, 49, 2887–2897. [Google Scholar] [CrossRef]
  25. Silling, S.A.; Zimmermann, M.; Abeyaratne, R. Deformation of a peridynamic bar. J. Elast. 2003, 73, 173–190. [Google Scholar] [CrossRef]
  26. Nishawala, V.V.; Ostoja-Starzewski, M. Peristatic solutions for finite one-and two-dimensional systems. Math. Mech. Solids 2017, 22, 1639–1653. [Google Scholar] [CrossRef]
  27. Yang, Z.; Naumenko, K.; Altenbach, H.; Ma, C.C.; Oterkus, E.; Oterkus, S. Some analytical solutions to peridynamic beam equations. ZAMM J. Appl. Math. Mech. 2022, 102, e202200132. [Google Scholar] [CrossRef]
  28. Yang, Z.; Naumenko, K.; Ma, C.C.; Altenbach, H.; Oterkus, E.; Oterkus, S. Some closed form series solutions to peridynamic plate equations. Mech. Res. Commun. 2022, 126, 104000. [Google Scholar] [CrossRef]
  29. Mikata, Y. Linear peridynamics for isotropic and anisotropic materials. Int. J. Solids Struct. 2019, 158, 116–127. [Google Scholar] [CrossRef]
  30. Mikata, Y. Analytical solutions of peristatics and peridynamics for 3D isotropic materials. Eur. J. Mech. A/Solids 2023, 101, 104978. [Google Scholar] [CrossRef]
Figure 1. Beam domain and interaction domain.
Figure 1. Beam domain and interaction domain.
Applsci 13 10025 g001
Figure 2. Simply supported edge with a fictitious region.
Figure 2. Simply supported edge with a fictitious region.
Applsci 13 10025 g002
Figure 3. Roller clamped edge with a fictitious region.
Figure 3. Roller clamped edge with a fictitious region.
Applsci 13 10025 g003
Figure 4. Simply supported beam.
Figure 4. Simply supported beam.
Applsci 13 10025 g004
Figure 5. Periodic extension of the deflection field for a simply supported beam.
Figure 5. Periodic extension of the deflection field for a simply supported beam.
Applsci 13 10025 g005
Figure 6. Clamped–simply supported beam: (a) original problem; (b) equivalent simply supported beam with applied moment.
Figure 6. Clamped–simply supported beam: (a) original problem; (b) equivalent simply supported beam with applied moment.
Applsci 13 10025 g006
Figure 7. Beam subjected to external moment: (a) moment applied at x 0 ; (b) force couple applied at x 0 .
Figure 7. Beam subjected to external moment: (a) moment applied at x 0 ; (b) force couple applied at x 0 .
Applsci 13 10025 g007
Figure 8. Clamped–clamped beam: (a) original problem; (b) equivalent simply supported beam with applied moments.
Figure 8. Clamped–clamped beam: (a) original problem; (b) equivalent simply supported beam with applied moments.
Applsci 13 10025 g008
Figure 9. Simply supported–roller-clamped beam.
Figure 9. Simply supported–roller-clamped beam.
Applsci 13 10025 g009
Figure 10. Periodic extension of the deflection line.
Figure 10. Periodic extension of the deflection line.
Applsci 13 10025 g010
Figure 11. Clamped–roller-clamped beam: (a) original problem; (b) equivalent simply supported beam with the applied moment.
Figure 11. Clamped–roller-clamped beam: (a) original problem; (b) equivalent simply supported beam with the applied moment.
Applsci 13 10025 g011
Figure 12. Solution representation for a cantilever beam.
Figure 12. Solution representation for a cantilever beam.
Applsci 13 10025 g012
Figure 13. Deflection vs. location of beams for different loadings and boundary conditions.
Figure 13. Deflection vs. location of beams for different loadings and boundary conditions.
Applsci 13 10025 g013
Figure 14. Simply supported beam under harmonic load.
Figure 14. Simply supported beam under harmonic load.
Applsci 13 10025 g014
Figure 15. Displacement vs. time for a simply supported beam under harmonic load. (a) Orbit of x = L / 2 , (b) orbit of x = L / 4 .
Figure 15. Displacement vs. time for a simply supported beam under harmonic load. (a) Orbit of x = L / 2 , (b) orbit of x = L / 4 .
Applsci 13 10025 g015
Figure 16. Simply supported beam under impact load.
Figure 16. Simply supported beam under impact load.
Applsci 13 10025 g016
Figure 17. Displacement vs. time for simply supported beams under impact load. (a) Orbit of x = L / 2 , (b) orbit of x = L / 4 .
Figure 17. Displacement vs. time for simply supported beams under impact load. (a) Orbit of x = L / 2 , (b) orbit of x = L / 4 .
Applsci 13 10025 g017
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

Yang, Z.; Naumenko, K.; Ma, C.-C.; Chen, Y. Closed-Form Analytical Solutions for the Deflection of Elastic Beams in a Peridynamic Framework. Appl. Sci. 2023, 13, 10025. https://doi.org/10.3390/app131810025

AMA Style

Yang Z, Naumenko K, Ma C-C, Chen Y. Closed-Form Analytical Solutions for the Deflection of Elastic Beams in a Peridynamic Framework. Applied Sciences. 2023; 13(18):10025. https://doi.org/10.3390/app131810025

Chicago/Turabian Style

Yang, Zhenghao, Konstantin Naumenko, Chien-Ching Ma, and Yang Chen. 2023. "Closed-Form Analytical Solutions for the Deflection of Elastic Beams in a Peridynamic Framework" Applied Sciences 13, no. 18: 10025. https://doi.org/10.3390/app131810025

APA Style

Yang, Z., Naumenko, K., Ma, C. -C., & Chen, Y. (2023). Closed-Form Analytical Solutions for the Deflection of Elastic Beams in a Peridynamic Framework. Applied Sciences, 13(18), 10025. https://doi.org/10.3390/app131810025

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