Next Article in Journal
Parameter Identification of Fiber Orientation Models Based on Direct Fiber Simulation with Smoothed Particle Hydrodynamics
Next Article in Special Issue
Variable Stiffness Composites: Optimal Design Studies
Previous Article in Journal
Metal–Organic Framework (MOF) through the Lens of Molecular Dynamics Simulation: Current Status and Future Perspective
Previous Article in Special Issue
Surface Modification of Flax Fibers for Manufacture of Engineering Thermoplastic Biocomposites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Peridynamic Mindlin Plate Formulation for Functionally Graded Materials

PeriDynamics Research Centre, Department of Naval Architecture, Ocean and Marine Engineering, University of Strathclyde, 100 Montrose Street, Glasgow G4 0LZ, UK
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2020, 4(2), 76; https://doi.org/10.3390/jcs4020076
Submission received: 1 June 2020 / Revised: 14 June 2020 / Accepted: 16 June 2020 / Published: 19 June 2020
(This article belongs to the Special Issue Feature Papers in Journal of Composites Science in 2020)

Abstract

:
In this study, a new peridynamic Mindlin plate formulation is presented which is suitable for the analysis of functionally graded materials. The governing equations of peridynamic formulation are obtained by using Euler-Lagrange equations in conjunction with Taylor’s expansion. To validate the new formulation, three different numerical benchmark problems are considered for a Mindlin plate subjected to simply supported, fully clamped and mixed (clamped-simply supported) boundary conditions. Peridynamic results are compared against results from finite element analysis and a good agreement is observed between the two methods.

1. Introduction

As the manufacturing technology advances, it becomes possible to design materials that have better properties with respect to traditional materials including metals and fibre-reinforced composite materials. Usage of fibre-reinforced composite materials is increasing in different industries including automotive, aerospace and marine fields. Although fibre-reinforced composite materials are light, have good impact properties and corrosion resistance, they suffer from delamination damage occurring between neighbouring plies which significantly reduces the load carrying capacities of these materials. Delamination damage can be avoided by continuously varying material properties in the form of functionally graded materials.
There are currently numerous studies in the literature focusing on plate formulations of functionally graded materials. Amongst these Vel and Batra [1] provided three-dimensional exact solution for the vibration of functionally graded rectangular plates. Shen [2] obtained nonlinear bending response of functionally graded plates subjected to transverse loading by using Reddy’s higher order shear deformation plate theory. Zenkour [3] presented generalised shear deformation theory for bending analysis of functionally graded plates. Bian et. al. [4] derived analytical solutions for functionally graded plates under cylindrical bending using first- and third-order shear deformation theories. Carrera et. al. [5] investigated the effect of thickness stretching in functionally graded plates by using Carrera’s Unified Formulation. Ferreira et. al. [6] used meshless method and third-order shear deformation theory to analyse functionally graded plates. Kashtalyan [7] utilised Plevako general solution of the equilibrium equations for inhomogeneous isotropic media to obtain three-dimensional elasticity solution for functionally graded simply supported plates. Xiang and Kang [8] performed bending analysis of functionally graded plates by using nth-order shear deformation theory and meshless global collocation method based on the thin plate spline radial basis function.
In this study, an alternative formulation is presented using peridynamics [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29] to analyse functionally graded plates based on Mindlin plate theory by taking into account transverse shear deformation. Peridynamics (PD) is a new continuum mechanics formulation. PD equations are in the form of integro-differential equations which do not contain spatial derivatives which makes it convenient for analysing cracks. Moreover, it has a length scale parameter that can allow analysing problems that cannot be represented by using classical formulations. On the other hand, partial differential equations are used in Refs. [1,2,3,4,5,6,7,8] and they do not contain length scale parameters. Since its introduction, there has been a rapid progress on peridynamics especially in the recent years. PD has been used to analyse different material systems such as metals [30], composites [31], functionally graded materials [32], concrete [33], graphene [34], etc. PD is not limited to elasticity, but can represent plastic [35], viscoelastic [36] and viscoplastic [37] material behaviour. Moreover, in addition to structural analysis, PD can also be used for the analysis of other fields including heat transfer [38], moisture diffusion [39], porous flow [40], etc. Simplified structures such as beams, plates and shells can also be represented in PD framework. Taylor and Steigmann [41] developed PD formulation for thin plates. Non-ordinary state-based Euler beam [42] and Kirchhoff plate [43] formulations were presented by O’Grady and Foster. To analyse relatively thick plates, Diyaroglu et. al. [44] introduced PD Timoshenko and Mindlin plate formulations by taking into account transverse shear deformations. This study only considers isotropic materials and Mindlin plate formulation is limited to a constant Poisson’s ratio of 1/3. Chowdhury et al. [45] introduced PD formulation for linear elastic shells. Note that there are fundamental differences between peridynamics and other non-local continuum mechanics approaches such as Eringen’s non-local elasticity formulation. Please refer to Madenci and Oterkus [46] for more information about the fundamentals of peridynamic theory. An extensive review of peridynamic research is given in Javili et al. [47]. To obtain governing equations of peridynamic Mindlin plate formulation for functionally graded materials Euler-Lagrange equations in conjunction with Taylor’s expansion are utilised. The formulation does not have any limitation on material constants as in bond-based peridynamics. To validate the current formulation, several benchmark problems are considered and peridynamic results are compared against results from finite element analysis.

2. Classical Mindlin Plate Formulation

Mindlin plate formulation was developed to analyse relatively thick plates by taking into account transverse shear deformations which is neglected in Kirchhoff plate theory suitable for thin plates. According to Mindlin plate theory, the displacement of a material point can be expressed in terms of displacement and rotation fields of the material points along the mid-plane (xy-plane)
u ( x , y , z , t ) = u ( x , y , 0 , t ) + z θ x ( x , y , 0 , t )
v ( x , y , z , t ) = v ( x , y , 0 , t ) + z θ y ( x , y , 0 , t )
w ( x , y , z , t ) = w ( x , y , 0 , t )
where θ x ( x , y , 0 , t )   and θ y ( x , y , 0 , t ) denote the rotation of the material points on the mid-plane about positive y-direction and negative x-direction, respectively. Hereafter, u ( x , y , 0 , t ) , v ( x , y , 0 , t ) , w ( x , y , 0 , t ) , θ x ( x , y , 0 , t ) and θ y ( x , y , 0 , t ) are written simply as u , v , w , θ x and θ y , respectively. The positive set of the degrees-of-freedom is shown in Figure 1.
Thus, the strain-displacement relationship can be written as
ε x x = u x + z θ x x
ε y y = v y + z θ y y
ε x y = ε y x = 1 2 [ ( u y + u x ) + z ( θ x y + θ y x ) ]
γ x z = κ s ( θ x + w x )
γ y z = κ s ( θ y + w y )
ε z z = 0
which can be also expressed in terms of indicial notation as
ε I J = 1 2 [ ( u I x J + u J x I ) + z ( θ I x J + θ J x I ) ]
γ I 3 = κ s ( θ I + w x I )
where κ s is introduced as shear coefficient. Note that the subscript indices, I ,   J ,   = 1 ( = x ) , 2 ( = y ) , and this convention will be applied throughout this study.
For planar isotropic materials, the stress–strain relationships can be written as:
σ x x = E ( z ) 1 ν ( z ) 2 ( ε x x + ν ε y y )
σ y y = E ( z ) 1 ν ( z ) 2 ( ε y y + ν ε x x )
σ x y = E ( z ) 2 ( 1 + ν ( z ) ) ( ε x y + ε y x )
τ x z = G ( z ) γ x z
τ y z = G ( z ) γ y z
where G = E ( z ) 2 ( 1 + ν ( z ) ) is the shear modulus.
Note that the transverse normal stress, σ z z , is considered to be small compared to in-plane stresses. Thus, it is discarded from the stress components set and this simplifies the 3-dimensional Hooke’s law and makes it a 2-dimensional plane-stress material constitutive law. The stress components can also be expressed in indicial notation as:
σ I J = C I J K L ε K L
τ I 3 = G ( z ) γ I 3
where in-plane stiffness tensor can be defined as
C I J K L = E ( z ) 1 ν ( z ) 2 [ 1 ν ( z ) 2 ( δ I L δ J K + δ I K δ J L ) + ν ( z ) δ I J δ K L ]
The strain energy per unit area of the plate can be casted as:
W = h 2 h 2 1 2 ( σ I J ε I J + τ I 3 γ I 3 ) d z
where h represents the thickness of the plate. Inserting Equations (3), (5) and (6) into (7) and rearranging indices yields
W = 1 2 h 2 h 2 ( E ( z ) 1 ν ( z ) 2 [ ( 1 ν ( z ) ) ε I J ε I J + ν ( z ) ε J J ε I I ] + G ( z ) γ I 3 γ I 3 ) d z
or
W = h 2 h 2 G ( z ) 2 d z ( u I x J u I x J + u I x J u J x I + u I x I u J x J ) + h 2 h 2 G ( z ) 2 z 2 d z ( θ I x J θ I x J + θ I x J θ J x I + θ I x I θ J x J ) + h 2 h 2 G ( z ) z d z ( u I x J θ I x J + θ I x J u J x I + u I x I θ J x J ) + ( h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] d z u I x I u J x J + h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] z 2 d z θ I x I θ J x J + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z u I x I θ J x J ) + h 2 h 2 κ s 2 G ( z ) 2 d z ( θ I + w x I ) ( θ I + w x I )

3. Peridynamic Mindlin Plate Formulation

Peridynamics (PD) is a non-local continuum mechanics formulation. Therefore, material points inside the solution domain can interact with each other material points in a non-local manner. The range of non-local interactions is defined as ‘horizon’, H . The equations of motion of peridynamics can be written as
ρ ( x ) u ¨ ( x , t ) = H ( t t ) d V + b ( x , t )
or in discrete form for the material point k as
ρ ( k ) u ¨ ( k ) = j = 1 N ( t ( k ) ( j ) t ( j ) ( k ) ) V ( j ) + b ( k )
where t ( k ) ( j ) and t ( j ) ( k ) are the peridynamic force densities between two material points k and j as shown in Figure 2 and N is the number of material points inside the horizon. Note that the PD equations of motion given in Equations (9) and (10) have the same form as in classical continuum mechanics except the integration or summation term.
The PD equations of motion can be derived by utilising Euler-Lagrange’s equation:
d d t L u ˙ ( k ) L u ( k ) = 0
where L = T U is the Lagrangian. The kinetic energy per unit area, T ¯ , can be expressed as
T ¯ = 1 2 h 2 h 2 ρ u ˙ ( x , y , z , t ) · u ˙ ( x , y , z , t ) d z
where u is the generalised displacement vector, which in this study can be defined as
u = ( u 1 u 2 θ 1 θ 2 w ) T
Substituting Equation (1) into Equation (12) results in
T ¯ = 1 2 h 2 h 2 ρ [ ( u ˙ 1 2 + u ˙ 2 2 ) + z 2 ( θ ˙ 1 2 + θ ˙ 2 2 ) + 2 z ( u ˙ 1 θ ˙ 1 + u ˙ 2 θ ˙ 2 ) + w ˙ 2 ] d z = 1 2 ρ [ h ( u ˙ 1 2 + u ˙ 2 2 + w ˙ 2 ) + h 3 12 ( θ ˙ 1 2 + θ ˙ 2 2 ) ]
The total kinetic energy of the body can be casted by integrating Equation (14) over the whole mid-plane as
T = A T ¯ d A = 1 2 A ρ [ h 3 12 ( θ ˙ 1 2 + θ ˙ 2 2 ) + h ( u ˙ 1 2 + u ˙ 2 2 + w ˙ 2 ) ] d A
which can be written in discretised form as
T = k T ¯ ( k ) A ( k ) = 1 2 k ρ ( k ) [ h 3 12 ( θ ˙ 1 ( k ) ) 2 + h 3 12 ( θ ˙ 2 ( k ) ) 2 + h ( ( u ˙ 1 ( k ) ) 2 + ( u ˙ 2 ( k ) ) 2 + w ˙ ( k ) 2 ) ] A ( k )
where A ( k ) is the area that contains the material point k .
The PD strain energy density function has a non-local form such that the strain energy of a certain material point k depends on both its displacement and all other material points in its family, which can be expressed as
W ( k ) = W ( k ) ( u ( k ) ,   u ( 1 k ) ,   u ( 2 k ) ,   u ( 3 k ) ,   )
where u ( k ) is the displacement vector of material point k and u ( i k ) ( i = 1 ,   2 ,   3 ,   ) is the displacement vector of the i th material point within the horizon of the material point k .
Similar to Equation (15), the total potential energy stored in the body can be obtained by summing potential energies of all material points including strain energy and energy due to external loads as
U = k W ( k ) ( u ( k ) ,   u ( 1 k ) ,   u ( 2 k ) ,   u ( 3 k ) ,   ) A ( k ) k b ( k ) u ( k ) A ( k )
where b is the body force density vector, which in this study has the following components
b = ( b u 1 b u 2 b θ 1 b θ 2 b z ) T
with b u , b θ and b z correspond to in-plane loads, moments and transverse body loads, respectively.
The first term of the Euler-Lagrange’s equation can be obtained by inserting Equations (15b) and (17) into Equation (11) as
d d t L u ˙ ( k ) = d d t T u ˙ ( k ) = ρ ( k ) { h u ¨ 1 ( k ) h u ¨ 2 ( k ) h 3 12 θ ¨ 1 ( k ) h 3 12 θ ¨ 2 ( k ) h w ¨ ( k ) } A ( k )
Similarly, the second term of the Euler-Lagrange’s equation becomes
L u ( k ) = U u ( k ) = u ( k ) n W ( n ) ( u ( n ) ,   u ( 1 n ) ,   u ( 2 n ) ,   u ( 3 n ) ,   ) A ( n ) u ( k ) n b ( n ) u ( n ) A ( n ) = n W ( n ) u ( k ) ( δ n k + δ n i k ) A ( n ) n b ( n ) δ n k A ( n ) = ( W ( k ) u ( k ) A ( k ) + j W ( j ) u ( k ) A ( j ) ) b ( k ) A ( k )
Inserting Equations (19) and (20) into the Euler-Lagrange’s equation yields:
ρ ( k ) { h u ¨ 1 ( k ) h u ¨ 2 ( k ) h 3 12 θ ¨ 1 ( k ) h 3 12 θ ¨ 2 ( k ) h w ¨ ( k ) } A ( k ) = { W ( k ) u 1 ( k ) A ( k ) + j W ( j ) u 1 ( k ) A ( j ) W ( k ) u 2 ( k ) A ( k ) + j W ( j ) u 2 ( k ) A ( j ) W ( k ) θ 1 ( k ) A ( k ) + j W ( j ) θ 1 ( k ) A ( j ) W ( k ) θ 2 ( k ) A ( k ) + j W ( j ) θ 2 ( k ) A ( j ) W ( k ) W ( k ) A ( k ) + j W ( j ) W ( k ) A ( j ) } + { b u 1 ( k ) b u 2 ( k ) b θ 1 ( k ) b θ 2 ( k ) b z ( k ) } A ( k )
In order to write the non-local form of strain energy function of the material point k , Equation (8b), it is necessary to transform all the local terms into an equivalent PD form by also considering PD strain energy expression given in Equation (16). As derived in Appendix A, the strain energy density function of the material point k and its family member j can be expressed as
W ( k ) = 12 π δ 3 ( h 2 h 2 G ( z ) 2 d z i [ ( u I ( i k ) u I ( k ) ) n I ( i k ) ( k ) ] 2 ξ ( i k ) ( k ) A ( i k ) + h 2 h 2 G ( z ) 2 z 2 d z i [ ( θ I ( i k ) θ I ( k ) ) n I ( i k ) ( k ) ] 2 ξ ( i k ) ( k ) A ( i k ) + h 2 h 2 G ( z ) z d z i ( u I ( i k ) u I ( k ) ) ( θ J ( i k ) θ J ( k ) ) ξ ( i k ) ( k ) n I ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) ) + 3 π δ 3 h 2 h 2 κ s 2 G ( z ) 2 d z i ( w ( i k ) w ( k ) + θ I ( k ) + θ I ( i k ) 2 ξ ( i k ) ( k ) n I ( i k ) ( k ) ) 2 ξ ( i k ) ( k ) A ( i k ) + ( 2 π δ 2 ) 2 ( h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] d z ( i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] z 2 d z ( i θ I ( i k ) θ I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) i θ J ( i k ) θ J ( k ) ξ ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) )
and
W ( j ) = 12 π δ 3 ( h 2 h 2 G ( z ) 2 d z i [ ( u I ( i j ) u I ( j ) ) n I ( i j ) ( j ) ] 2 ξ ( i j ) ( j ) A ( i j ) + h 2 h 2 G ( z ) 2 z 2 d z i [ ( θ I ( i j ) θ I ( j ) ) n I ( i j ) ( j ) ] 2 ξ ( i j ) ( j ) A ( i j ) + h 2 h 2 G ( z ) z d z i ( u I ( i j ) u I ( j ) ) ( θ J ( i j ) θ J ( j ) ) ξ ( i j ) ( j ) n I ( i j ) ( j ) n J ( i j ) ( j ) A ( i j ) ) + 3 π δ 3 h 2 h 2 κ s 2 G ( z ) 2 d z i ( w ( i j ) w ( j ) + θ I ( j ) + θ I ( i j ) 2 ξ ( i j ) ( j ) n I ( i j ) ( j ) ) 2 ξ ( i j ) ( j ) A ( i j ) + ( 2 π δ 2 ) 2 ( h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] d z ( i u I ( i j ) u I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] z 2 d z ( i θ I ( i j ) θ I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z i u I ( i j ) u I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) i θ J ( i j ) θ J ( j ) ξ ( i j ) ( j ) n J ( i j ) ( j ) A ( i j ) )
with n 1 = cos φ , n 2 = sin φ and φ is the orientation of the peridynamic interaction with respect to horizontal x axis .
Inserting Equations (22a) and (22b) into Equation (21) results in complete PD equations of motion for functionally graded Mindlin plates as
ρ ( k ) h u ¨ L ( k ) = 24 π δ 3 ( h 2 h 2 G ( z ) z 2 d z j u I ( j ) u I ( k ) ξ ( j ) ( k ) n I ( j ) ( k ) n L ( j ) ( k ) A ( j ) + h 2 h 2 G ( z ) z d z i θ J ( i k ) θ J ( k ) ξ ( i k ) ( k ) n L ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) ) + ( 2 π δ 2 ) 2 [ h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) d z j n L ( j ) ( k ) ξ ( j ) ( k ) ( i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) + i u I ( i j ) u I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) A ( j ) + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z j n L ( j ) ( k ) ξ ( j ) ( k ) ( i θ J ( i k ) θ J ( k ) ξ ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) + i θ J ( i j ) θ J ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) A ( j ) ] + b u L ( k )
ρ ( k ) h 3 12 θ ¨ L ( k ) = 24 π δ 3 ( h 2 h 2 G ( z ) z 2 d z j θ I ( j ) θ I ( k ) ξ ( j ) ( k ) n I ( j ) ( k ) n L ( j ) ( k ) A ( j ) + h 2 h 2 G ( z ) z d z j u I ( j ) u I ( k ) ξ ( j ) ( k ) n I ( j ) ( k ) n L ( j ) ( k ) A ( j ) ) 3 π δ 3 κ s 2 h 2 h 2 G ( z ) d z j ( w ( j ) w ( k ) + θ I ( k ) + θ I ( j ) 2 ξ ( j ) ( k ) n I ( j ) ( k ) ) n L ( j ) ( k ) A ( j ) + ( 2 π δ 2 ) 2 [ h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z 2 d z j n L ( j ) ( k ) ξ ( j ) ( k ) ( i θ I ( i k ) θ I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) + i θ I ( i j ) θ I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) A ( j ) + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z j n L ( j ) ( k ) ξ ( j ) ( k ) ( i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) + i u I ( i j ) u I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) A ( j ) ] + b θ L ( k )
ρ ( k ) h w ¨ ( k ) = 6 π δ 3 h κ s 2 h 2 h 2 G ( z ) d z i ( w ( j ) w ( k ) ξ ( j ) ( k ) + θ J ( k ) + θ J ( j ) 2 n J ( j ) ( k ) ) V ( j ) + b z ( k )
In particular, when Poisson’s ratio, ν ( z ) = 1 3 , Equation (23) can be reduced to PD bond-based formulation as
ρ ( k ) h u ¨ L ( k ) = 24 π δ 3 ( h 2 h 2 G ( z ) d z j u I ( j ) u I ( k ) ξ ( j ) ( k ) n I ( j ) ( k ) n L ( j ) ( k ) A ( j ) + h 2 h 2 G ( z ) z d z i ( θ J ( i k ) θ J ( k ) ) ξ ( i k ) ( k ) n L ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) ) + b u L ( k )
ρ ( k ) h 3 12 θ ¨ L ( k ) = 24 π δ 3 ( h 2 h 2 G ( z ) z 2 d z j θ I ( j ) θ I ( k ) ξ ( j ) ( k ) n I ( j ) ( k ) n L ( j ) ( k ) A ( j ) + h 2 h 2 G ( z ) z d z j u I ( j ) u I ( k ) ξ ( j ) ( k ) n I ( j ) ( k ) n L ( j ) ( k ) A ( j ) ) 3 π δ 3 κ s 2 h 2 h 2 G ( z ) d z j ( w ( j ) w ( k ) + θ I ( k ) + θ I ( j ) 2 ξ ( j ) ( k ) n I ( j ) ( k ) ) n L ( j ) ( k ) A ( j ) + b θ L ( k )
ρ ( k ) h w ¨ ( k ) = 6 π δ 3 κ s 2 h 2 h 2 G ( z ) d z i ( w ( j ) w ( k ) ξ ( j ) ( k ) + θ J ( k ) + θ J ( j ) 2 n J ( j ) ( k ) ) A ( j ) + b z ( k )

4. Numerical Results

To verify the validity of the PD formulation for functionally graded Mindlin plates, the PD solutions are compared with the corresponding finite element (FE) analysis results. In this study, the functionally graded material properties are chosen as Young’s Modulus, E ( z ) and shear modulus G ( z ) and they are assumed to vary linearly through the thickness as
E ( z ) = ( E t E b ) z h + 1 2 ( E t + E b )     ( GPa )
G ( z ) = E ( z ) 2 ( 1 + 0.3 )
ν ( z ) = 0.3
where E t and E b denote the Young’s modulus of the top and bottom surfaces of the plate, and h denotes the total thickness of the plate. The shear correction coefficient is chosen as κ s 2 = π 2 12 .
In the following three numerical cases, a square plate with length and width of L = W = 1   m and thickness of h = 0.15   m is considered. The plate is subjected to different boundary conditions. The Young’s modulus of the top and bottom surface are chosen as E t = 200   GPa and E t = 100   GPa .
The PD models are discretized into one single row of material points through the thickness and 65 × 65 material points throughout the x y plane. Thus, the distance between two adjacent material points is Δ x = 1 65   m and the area attached on each material points is Δ A = Δ x 2 . A fictitious region is introduced outside the edges as the external boundaries with a width of 6 Δ x to apply boundary conditions as explained in Appendix B. The horizon size can be approximately chosen as δ = 3 Δ x .
The corresponding FE models are created in ANSYS by using SHELL181 elements with 50 × 50 elements throughout the plate. In order to obtain the functionally graded character, the model is divided into 50 layers with varying homogeneous material properties throughout the thickness. The Young’s modulus varies gradually over the thickness from the first layer E 1 = 101   GPa to the last layer E 50 = 199   GPa , as shown in Figure 3. The Poisson’s ratio, ν = 0.3 , is applied in ANSYS.

4.1. Simply Supported Functionally Graded Mindlin Plate

In the first example case, a simply supported functionally graded Mindlin plate subjected to a distributed load of p = 100,000   N / m through the central line is taken into consideration (see Figure 4). For the PD model, the load is transformed into a body load of b z = pW 65 Δ A = 6,500,000   N / m 2 and it is imposed on a row of material points through the central line of the plate as shown in Figure 5.
Assuming the coordinate system is attached at the geometric centre of the plate, the following boundary conditions are applied in ANSYS:
u ( L 2 ,   y ,   0 ) = u ( L 2 ,   y ,   0 ) = u ( x , W 2 ,   0 ) = u ( x ,   W 2 ,   0 ) = 0
v ( L 2 ,   y ,   0 ) = v ( L 2 ,   y ,   0 ) = v ( x , W 2 ,   0 ) = v ( x ,   W 2 ,   0 ) = 0
w ( L 2 ,   y ,   0 ) = w ( L 2 ,   y ,   0 ) = w ( x , W 2 ,   0 ) = w ( x ,   W 2 ,   0 ) = 0
θ y ( L 2 ,   y ,   0 ) = θ y ( L 2 ,   y ,   0 ) = 0
θ x ( x , W 2 ,   0 ) = θ x ( x ,   W 2 ,   0 ) = 0
The PD results for in-plane and transverse displacements, and rotations are obtained and compared with FEA results for the material points located along central x- and y- axes. As depicted in Figure 6, PD and FEA results agree very well with each other.

4.2. Fully Clamped Functionally Graded Mindlin Plate

In the second example case, the functionally graded Mindlin plate considered in the previous example is subjected to fully clamped boundary condition (see Figure 7). In ANSYS, the clamped boundary condition is achieved by constraining all degrees of freedom along the external boundaries.
Based on the comparison between the PD and FEA results as shown in Figure 8, it can be concluded that current PD formulation can also provide accurate results for fully clamped boundary condition.

4.3. Functionally Graded Mindlin Plate Subjected to Mixed Boundary Conditions

The last numerical case aims to verify the current PD formulation for mixed boundary conditions, i.e., clamped—simply supported. As shown in Figure 9, edges along the horizontal direction are subjected to clamped boundary conditions whereas the remaining two edges are subjected to simply supported boundary conditions.
Assuming the coordinate system is attached at the geometric centre of the plate, the following boundary conditions are applied in ANSYS:
u ( L 2 ,   y ,   0 ) = u ( L 2 ,   y ,   0 ) = u ( x , W 2 ,   0 ) = u ( x ,   W 2 ,   0 ) = 0
v ( L 2 ,   y ,   0 ) = v ( L 2 ,   y ,   0 ) = v ( x , W 2 ,   0 ) = v ( x ,   W 2 ,   0 ) = 0
w ( L 2 ,   y ,   0 ) = w ( L 2 ,   y ,   0 ) = w ( x , W 2 ,   0 ) = w ( x ,   W 2 ,   0 ) = 0
θ x ( L 2 ,   y ,   0 ) = θ x ( L 2 ,   y ,   0 ) = θ x ( x , W 2 ,   0 ) = θ x ( x ,   W 2 ,   0 ) = 0
θ y ( L 2 ,   y ,   0 ) = θ y ( L 2 ,   y ,   0 ) = 0
As depicted in Figure 10, also for this mixed-boundary conditions case, PD and FEM results agree well with each other.

5. Conclusions

In this study, a new peridynamic formulation was presented for functionally graded Mindlin plates. The governing equations were obtained by using Euler-Lagrange equations. To validate the formulation, three different benchmark problems were considered for simply supported, fully clamped and mixed (clamped-simply supported) boundary conditions. Peridynamic results were compared against finite element analysis results and a good agreement is observed for both in-plane and transverse displacements, and rotations. Hence, it can be concluded that the present peridynamic formulation can be a suitable alternative for the analysis of functionally graded Mindlin plates subjected to different types of boundary conditions.

Author Contributions

Conceptualisation, Z.Y., E.O. and S.O.; methodology, Z.Y. and E.O.; software, Z.Y.; validation, Z.Y.; writing—original draft preparation, Z.Y.; writing—review and editing, E.O. and S.O.; visualisation, Z.Y. and S.O.; supervision, E.O. and S.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

As explained above, according to classical continuum mechanics, the strain energy per unit area of functionally graded Mindlin plates can be written as
W = h 2 h 2 G ( z ) 2 d z ( u I x J u I x J + u I x J u J x I + u I x I u J x J ) + h 2 h 2 G ( z ) 2 z 2 d z ( θ I x J θ I x J + θ I x J θ J x I + θ I x I θ J x J ) + h 2 h 2 G ( z ) z d z ( u I x J θ I x J + θ I x J u J x I + u I x I θ J x J ) + ( h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] d z u I x I u J x J + h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] z 2 d z θ I x I θ J x J + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z u I x I θ J x J ) + h 2 h 2 κ s 2 G ( z ) 2 d z ( θ I + w x I ) ( θ I + w x I )
In order to obtain the strain energy function in PD form, it is necessary to transform each local term in Equation (A1) into their equivalent non-local form. This can be achieved by using Taylor expansion.
Figure A1. Peridynamic interaction between two material points.
Figure A1. Peridynamic interaction between two material points.
Jcs 04 00076 g0a1
As shown in Figure A1, the in-plane displacement function u ( x , y ) , can be Taylor expanded up to 1st order terms about point x :
u I ( x + ξ ) u I ( x ) = u I ( x ) x J ξ n J
u K ( x + ξ ) u K ( x ) = u K ( x ) x L ξ n L
where ξ = | ξ | , and unit direction vector n is defined as
n = { n 1 n 2 } = { cos φ sin φ }
Multiplying Equation (A2a) with Equation (A2b) gives
[ u I ( x + ξ ) u I ( x ) ] [ u K ( x + ξ ) u K ( x ) ] ξ = u I ( x ) x J u K ( x ) x L ξ n J n L
Multiplying both sides of Equation (A4) twice by directional vector yields
[ u I ( x + ξ ) u I ( x ) ] [ u K ( x + ξ ) u K ( x ) ] ξ n R n S = u I ( x ) x J u K ( x ) x L ξ n J n L n R n S
Considering x as a fixed point, integrating both sides of Equation (A5) over a circular domain with centre of x and radius of δ yields:
0 2 π 0 δ [ u I ( x + ξ ) u I ( x ) ] [ u K ( x + ξ ) u K ( x ) ] ξ n R n S ξ d ξ d φ = u I ( x ) x J u K ( x ) x L 0 2 π 0 δ ξ n J n L n R n S ξ d ξ d φ = u I ( x ) x J u K ( x ) x L ( δ J L δ R S + δ J R δ L S + δ J S δ L R ) π δ 3 12 = π δ 3 12 ( u I ( x ) x J u K ( x ) x J δ R S + u I ( x ) x R u K ( x ) x S + u I ( x ) x S u K ( x ) x R )
Multiplying both sides of Equation (A6) by δ R I δ S K results in:
0 2 π 0 δ [ u I ( x + ξ ) u I ( x ) ] [ u K ( x + ξ ) u K ( x ) ] ξ n R n S δ R I δ S K ξ d ξ d φ = π δ 3 12 ( u I ( x ) x J u K ( x ) x J δ R S + u I ( x ) x R u K ( x ) x S + u I ( x ) x S u K ( x ) x R ) δ R I δ S K
Rearranging the dummy indices gives:
u I ( x ) x J u I ( x ) x J + u I ( x ) x I u J ( x ) x J + u I ( x ) x J u J ( x ) x I = 12 π δ 3 0 2 π 0 δ [ u I ( x + ξ ) u I ( x ) ] [ u J ( x + ξ ) u J ( x ) ] ξ n I n J ξ d ξ d φ
which can be written in the discretized form as
u I ( k ) x J u I ( k ) x J + u I ( k ) x I u J ( k ) x J + u I ( k ) x J u J ( k ) x I = 12 π δ 3 i ( u I ( i k ) u I ( k ) ) ( u J ( i k ) u J ( k ) ) ξ ( i k ) ( k ) n I ( i k ) ( k ) n J ( i k ) ( k ) A ( i k )
If we make an analogy with Equation (A9a), following expressions can be obtained for the rotation field:
θ I ( k ) x J θ I ( k ) x J + θ I ( k ) x I θ J ( k ) x J + θ I ( k ) x J θ J ( k ) x I = 12 π δ 3 i ( θ I ( i k ) θ I ( k ) ) ( θ J ( i k ) θ J ( k ) ) ξ ( i k ) ( k ) n I ( i k ) ( k ) n J ( i k ) ( k ) A ( i k )
and
u I ( k ) x J θ I ( k ) x J + u I ( k ) x I θ J ( k ) x J + u I ( k ) x J θ J ( k ) x I = 12 π δ 3 i ( u I ( i k ) u I ( k ) ) ( θ J ( i k ) θ J ( k ) ) ξ ( i k ) ( k ) n I ( i k ) ( k ) n J ( i k ) ( k ) A ( i k )
Recalling Equation (A2):
u I ( x + ξ ) u I ( x ) = u I ( x ) x J ξ n J
and multiplying Equation (A10) by a directional vector gives:
u I ( x + ξ ) u I ( x ) ξ n K = u I ( x ) x J n J n K
Considering x as a fixed point, integrating both sides of Equation (A11) over a circular domain with centre of x and radius of δ yields:
0 2 π 0 δ u I ( x + ξ ) u I ( x ) ξ n K ξ d ξ d φ = u I ( x ) x J 0 2 π 0 δ n J n K ξ d ξ d φ = u I ( x ) x J π δ 2 2 δ J K = u I ( x ) x K π δ 2 2
which results in
u I ( x ) x K = 2 π δ 2 0 2 π 0 δ u I ( x + ξ ) u I ( x ) ξ n K ξ d ξ d φ
Multiplying both sides of Equation (A13) by δ I K gives
u I ( x ) x I = 2 π δ 2 0 2 π 0 δ u I ( x + ξ ) u I ( x ) ξ n I ξ d ξ d φ
Rewriting Equation (A14a) with a different index gives:
u J ( x ) x J = 2 π δ 2 0 2 π 0 δ u J ( x + ξ ) u J ( x ) ξ n J ξ d ξ d φ
Multiplying Equation (A14a) with (A14b) yields:
u I ( x ) x I u J ( x ) x J = ( 2 π δ 2 ) 2 0 2 π 0 δ u I ( x + ξ ) u I ( x ) ξ n I ξ d ξ d φ 0 2 π 0 δ u J ( x + ξ ) u J ( x ) ξ n J ξ d ξ d φ
which can be written in discretized form as
u I ( k ) x I u J ( k ) x J = ( 2 π δ 2 ) 2 i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) i u J ( i k ) u J ( k ) ξ ( i k ) ( k ) n J ( i k ) ( k ) A ( i k )
Similarly, following expressions can also be obtained
θ I ( k ) x I θ J ( k ) x J = ( 2 π δ 2 ) 2 i θ I ( i k ) θ I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) i θ J ( i k ) θ J ( k ) ξ ( i k ) ( k ) n J ( i k ) ( k ) A ( i k )
u I ( k ) x I θ J ( k ) x J = ( 2 π δ 2 ) 2 i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) i θ J ( i k ) θ J ( k ) ξ ( i k ) ( k ) n J ( i k ) ( k ) A ( i k )
Similar to Equation (A2), the following relationship can be established for the transverse displacements:
w ( x + ξ ) w ( x ) = w ( x ) x I ξ n I
The rotations of material point x can be estimated by calculating the average rotation of point x and its family member point x + ξ as
θ I ( x + ξ ) + θ I ( x ) 2 = θ I ( x )
If Equation (A18) is multiplied by ξ n I :
θ I ( x + ξ ) + θ I ( x ) 2 ξ n I = θ I ( x ) ξ n I
and added with (A17) yields:
w ( x + ξ ) w ( x ) + θ I ( x + ξ ) + θ I ( x ) 2 ξ n I = w ( x ) x I ξ n I + θ I ( x ) ξ n I
Rewriting Equation (A20a) with a different index results in:
w ( x + ξ ) w ( x ) + θ J ( x + ξ ) + θ J ( x ) 2 ξ n J = w ( x ) x J ξ n J + θ J ( x ) ξ n J
Multiplying Equation (A20a) with (A20b) and then dividing each term by ξ yields:
[ w ( x + ξ ) w ( x ) + θ J ( x + ξ ) + θ J ( x ) 2 ξ n J ] [ w ( x + ξ ) w ( x ) + θ I ( x + ξ ) + θ I ( x ) 2 ξ n I ] ξ = ( w ( x ) x I + θ I ( x ) ) ( w ( x ) x J + θ J ( x ) ) ξ n I n J
Considering x as a fixed point, integrating both sides of Equation (A21) over a circular domain with centre of x and radius of δ results in
0 2 π 0 δ [ w ( x + ξ ) w ( x ) + θ J ( x + ξ ) + θ J ( x ) 2 ξ n J ] [ w ( x + ξ ) w ( x ) + θ I ( x + ξ ) + θ I ( x ) 2 ξ n I ] ξ ξ d ξ d φ = ( w ( x ) x I + θ I ( x ) ) ( w ( x ) x J + θ J ( x ) ) 0 2 π 0 δ ξ n I n J d ξ d φ = ( w ( x ) x I + θ I ( x ) ) ( w ( x ) x J + θ J ( x ) ) π δ 3 3 δ I J = π δ 2 3 ( w ( x ) x I + θ I ( x ) ) 2
which gives
( w ( x ) x I + θ I ( x ) ) 2 = 3 π δ 2 0 2 π 0 δ [ w ( x + ξ ) w ( x ) + θ J ( x + ξ ) + θ J ( x ) 2 ξ n J ] [ w ( x + ξ ) w ( x ) + θ I ( x + ξ ) + θ I ( x ) 2 ξ n I ] ξ ξ d ξ d ξ d ξ d φ
Equation (A22b) can be written in discretized form as
( w ( k ) x I + θ I ( k ) ) 2 = 3 π δ 3 i ( w ( i k ) w ( k ) + θ I ( k ) + θ I ( i k ) 2 ξ ( i k ) ( k ) n I ( i k ) ( k ) ) ( w ( i k ) w ( k ) + θ J ( k ) + θ J ( i k ) 2 ξ ( i k ) ( k ) n J ( i k ) ( k ) ) ξ ( i k ) ( k ) A ( i k )
Finally, combining Equations (A9), (A16), (A23) and (A1) gives the strain energy density of the material point k in PD form as
W ( k ) = 12 π δ 3 ( h 2 h 2 G ( z ) 2 d z i [ ( u I ( i k ) u I ( k ) ) n I ( i k ) ( k ) ] 2 ξ ( i k ) ( k ) A ( i k ) + h 2 h 2 G ( z ) 2 z 2 d z i [ ( θ I ( i k ) θ I ( k ) ) n I ( i k ) ( k ) ] 2 ξ ( i k ) ( k ) A ( i k ) + h 2 h 2 G ( z ) z d z i ( u I ( i k ) u I ( k ) ) ( θ J ( i k ) θ J ( k ) ) ξ ( i k ) ( k ) ( n I ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) ) + 3 π δ 3 h 2 h 2 κ s 2 G ( z ) 2 d z i ( w ( i k ) w ( k ) + θ I ( k ) + θ I ( i k ) 2 ξ ( i k ) ( k ) n I ( i k ) ( k ) ) 2 ξ ( i k ) ( k ) A ( i k ) + ( 2 π δ 2 ) 2 ( h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] d z ( i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] z 2 d z ( i θ I ( i k ) θ I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z i u I ( i k ) u I ( k ) ξ ( i k ) ( k ) n I ( i k ) ( k ) A ( i k ) i θ J ( i k ) θ J ( k ) ξ ( i k ) ( k ) n J ( i k ) ( k ) A ( i k ) )
Regarding the strain energy for the material point j , a similar form will hold if we replace the index k with j as
W ( j ) = 12 π δ 3 ( h 2 h 2 G ( z ) 2 z 2 d z i [ ( θ I ( i j ) θ I ( j ) ) n I ( i j ) ( j ) ] 2 ξ ( i j ) ( j ) A ( i j ) + h 2 h 2 G ( z ) 2 z 2 d z i [ ( θ I ( i j ) θ I ( j ) ) n I ( i j ) ( j ) ] 2 ξ ( i j ) ( j ) A ( i j ) + h 2 h 2 G ( z ) z d z i ( u I ( i j ) u I ( j ) ) ( θ J ( i j ) θ J ( j ) ) ξ ( i j ) ( j ) n I ( i j ) ( j ) n J ( i j ) ( j ) A ( i j ) ) + 3 π δ 3 h 2 h 2 κ s 2 G ( z ) 2 d z i ( w ( i j ) w ( j ) + θ I ( j ) + θ I ( i j ) 2 ξ ( i j ) ( j ) n I ( i j ) ( j ) ) 2 ξ ( i j ) ( j ) A ( i j ) + ( 2 π δ 2 ) 2 ( h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] d z ( i u I ( i j ) u I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 2 [ 1 ν ( z ) ] z 2 d z ( i θ I ( i j ) θ I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ) 2 + h 2 h 2 G ( z ) 3 ν ( z ) 1 1 ν ( z ) z d z i u I ( i j ) u I ( j ) ξ ( i j ) ( j ) n I ( i j ) ( j ) A ( i j ) ( i θ J ( i j ) θ J ( j ) ξ ( i j ) ( j ) n J ( i j ) ( j ) A ( i j ) )

Appendix B

Application of boundary conditions in PD theory is different than the classical theory. Boundary conditions can be applied by creating fictitious domain, R c , outside of the actual solution domain, R as shown in Figure A2. The length of this layer can be chosen as the double size of the horizon if ν ( z ) 1 / 3 , or the size of the horizon if ν ( z ) = 1 / 3 . The application procedure for two common types of boundary conditions, i.e., clamped and simply supported, is given below.
Figure A2. The fictitious domain, R c , outside of the actual solution domain, R .
Figure A2. The fictitious domain, R c , outside of the actual solution domain, R .
Jcs 04 00076 g0a2

Clamped Boundary Condition

Figure A3. Application of clamped boundary conditions in peridynamic theory.
Figure A3. Application of clamped boundary conditions in peridynamic theory.
Jcs 04 00076 g0a3
To implement the clamped boundary condition, a fictitious boundary layer is created outside the actual material domain (see Figure A3). The horizon size can be approximately chosen as δ = 3 Δ x in which the discretisation size is Δ x . The size of the fictitious region can be specified to be equal to 6 Δ x for ν ( z ) 1 / 3 and 3 Δ x for ν ( z ) = 1 / 3 .
According to classical theory, the clamped boundary condition imposes zero in-plane displacements, zero transverse displacement and zero rotations for the material point adjacent to the clamped end as
u I = 0
w = 0
θ I = 0
In this study, these conditions can be achieved by enforcing mirror image of the transverse displacement field for the material points adjacent to the clamped end and anti-symmetric image of in-plane and rotational displacements fields as following
u I ( n ) ( i ) = u I ( n ) ( i * )
w ( n ) ( i ) = w ( n ) ( i * )
w ( n ) ( 1 ) = w ( n ) ( 1 * ) = 0
θ I ( n ) ( i ) = θ I ( n ) ( i * )   for   { i = 1 , 2 , , 6               if   ν ( z ) 1 / 3 i = 1 , 2 , 3                           if   ν ( z ) = 1 / 3

Simply Supported Boundary Condition

Figure A4. Application of simply supported boundary condition in peridynamic theory.
Figure A4. Application of simply supported boundary condition in peridynamic theory.
Jcs 04 00076 g0a4
To implement the simply supported boundary condition, the fictitious layer is introduced outside the actual material domain (see Figure A4), whose size is again chosen to be equal to 6 Δ x for ν ( z ) 1 / 3 and 3 Δ x for ν ( z ) = 1 / 3 . According to classical theory, the simply supported boundary condition can be achieved by imposing zero in-plane displacements, zero transverse displacements and zero curvature for the material point adjacent to the boundaries as
{ u ( x , y ) = v ( x , y ) = w ( x , y ) = 0 θ y ( x , y ) = 0   for   x = constant
and
{ u ( x , y ) = v ( x , y ) = w ( x , y ) = 0 θ x ( x , y ) = 0   for   y = constant
In our study, the in-plane and transverse displacements boundary conditions can be applied by enforcing anti-symmetrical displacement fields to the material points in the fictitious region with respect to the actual region as
{ u I ( n ) ( i ) = u I ( n ) ( i * ) w ( n ) ( i ) = w ( n ) ( i * )   for   { i = 1 , 2 , , 6               if   ν ( z ) 1 / 3 i = 1 , 2 , 3                           if   ν ( z ) = 1 / 3
The curvature condition can be satisfied by enforcing the following conditions adjacent to the edge x I = c o n s t a n t as
θ I ( n ) ( i ) = θ I ( n ) ( i * )   and   θ J ( n ) ( i ) = θ J ( n ) ( i * ) ( I J ) for   { i = 1 , 2 , , 6               if   ν ( z ) 1 / 3 i = 1 , 2 , 3                           if   ν ( z ) = 1 / 3

References

  1. Vel, S.; Batra, R.C. Three-dimensional exact solution for the vibration of functionally graded rectangular plates. J. Sound Vib. 2004, 272, 703–730. [Google Scholar] [CrossRef]
  2. Shen, H.-S. Nonlinear bending response of functionally graded plates subjected to transverse loads and in thermal environments. Int. J. Mech. Sci. 2002, 44, 561–584. [Google Scholar] [CrossRef]
  3. Zenkour, A. Generalized shear deformation theory for bending analysis of functionally graded plates. Appl. Math. Model. 2006, 30, 67–84. [Google Scholar] [CrossRef] [Green Version]
  4. Bian, Z.; Chen, W.; Lim, C.W.; Zhang, N. Analytical solutions for single- and multi-span functionally graded plates in cylindrical bending. Int. J. Solids Struct. 2005, 42, 6433–6456. [Google Scholar] [CrossRef] [Green Version]
  5. Carrera, E.; Brischetto, S.; Cinefra, M.; Soave, M. Effects of thickness stretching in functionally graded plates and shells. Compos. Part B Eng. 2011, 42, 123–133. [Google Scholar] [CrossRef] [Green Version]
  6. Ferreira, A.; Batra, R.C.; Roque, C.; Qian, L.; Martins, P. Static analysis of functionally graded plates using third-order shear deformation theory and a meshless method. Compos. Struct. 2005, 69, 449–457. [Google Scholar] [CrossRef]
  7. Kashtalyan, M. Three-dimensional elasticity solution for bending of functionally graded rectangular plates. Eur. J. Mech. A Solids 2004, 23, 853–864. [Google Scholar] [CrossRef]
  8. Xiang, S.; Kang, G.-W. A nth-order shear deformation theory for the bending analysis on the functionally graded plates. Eur. J. Mech. A Solids 2013, 37, 336–343. [Google Scholar] [CrossRef]
  9. Silling, S. Reformulation of elasticity theory for discontinuities and long-range forces. J. Mech. Phys. Solids 2000, 48, 175–209. [Google Scholar] [CrossRef] [Green Version]
  10. Silling, S.A.; Epton, M.; Weckner, O.; Xu, J.; Askari, E. Peridynamic States and Constitutive Modeling. J. Elast. 2007, 88, 151–184. [Google Scholar] [CrossRef] [Green Version]
  11. Jenabidehkordi, A.; Rabczuk, T.; Ma, B.; Dui, G.; Yang, S.; Xin, L. The Multi-Horizon Peridynamics. Comput. Model. Eng. Sci. 2019, 121, 493–500. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, X.; Kulkarni, S.S.; Tabarraei, A. Concurrent coupling of peridynamics and classical elasticity for elastodynamic problems. Comput. Methods Appl. Mech. Eng. 2019, 344, 251–275. [Google Scholar] [CrossRef]
  13. Ni, T.; Zaccariotto, M.; Zhu, Q.-Z.; Galvanetto, U. Static solution of crack propagation problems in Peridynamics. Comput. Methods Appl. Mech. Eng. 2019, 346, 126–151. [Google Scholar] [CrossRef]
  14. Chowdhury, S.R.; Roy, P.; Roy, D.; Reddy, J.N. A modified peridynamics correspondence principle: Removal of zero-energy deformation and other implications. Comput. Methods Appl. Mech. Eng. 2019, 346, 530–549. [Google Scholar] [CrossRef]
  15. Liu, S.; Fang, G.; Liang, J.; Fu, M.; Wang, B. A new type of peridynamics: Element-based peridynamics. Comput. Methods Appl. Mech. Eng. 2020, 366, 113098. [Google Scholar] [CrossRef]
  16. Diehl, P.; Prudhomme, S.; Lévesque, M. A Review of Benchmark Experiments for the Validation of Peridynamics Models. J. Peridynamics Nonlocal Model. 2019, 1, 14–35. [Google Scholar] [CrossRef]
  17. Katiyar, A.; Agrawal, S.; Ouchi, H.; Seleson, P.; Foster, J.T.; Sharma, M.M. A general peridynamics model for multiphase transport of non-Newtonian compressible fluids in porous media. J. Comput. Phys. 2020, 402, 109075. [Google Scholar] [CrossRef]
  18. Song, X.; Khalili, N. A peridynamics model for strain localization analysis of geomaterials. Int. J. Numer. Anal. Methods Géoméch. 2018, 43, 77–96. [Google Scholar] [CrossRef] [Green Version]
  19. Ozdemir, M.; Kefal, A.; Imachi, M.; Tanaka, S.; Oterkus, E. Dynamic fracture analysis of functionally graded materials using ordinary state-based peridynamics. Compos. Struct. 2020, 244, 112296. [Google Scholar] [CrossRef]
  20. Oterkus, E.; Madenci, E. Peridynamics for failure prediction in composites. In Proceedings of the 53rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference 20th AIAA/ASME/AHS Adaptive Structures Conference 14th AIAA, Honolulu, HI, USA, 23–26 April 2012; p. 1692. [Google Scholar]
  21. Yang, Z.; Oterkus, E.; Nguyen, C.T.; Oterkus, S. Implementation of peridynamic beam and plate formulations in finite element framework. Contin. Mech. Thermodyn. 2018, 31, 301–315. [Google Scholar] [CrossRef] [Green Version]
  22. Diyaroglu, C.; Oterkus, S.; Oterkus, E.; Madenci, E. Peridynamic Modeling of Diffusion by Using Finite-Element Analysis. IEEE Trans. Components Packag. Manuf. Technol. 2017, 7, 1–9. [Google Scholar] [CrossRef] [Green Version]
  23. Zhu, N.; De Meo, D.; Oterkus, E. Modelling of Granular Fracture in Polycrystalline Materials Using Ordinary State-Based Peridynamics. Materials 2016, 9, 977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Imachi, M.; Tanaka, S.; Bui, T.Q.; Oterkus, S.; Oterkus, E. A computational approach based on ordinary state-based peridynamics with new transition bond for dynamic fracture analysis. Eng. Fract. Mech. 2019, 206, 359–374. [Google Scholar] [CrossRef] [Green Version]
  25. Gao, Y.; Oterkus, S. Ordinary state-based peridynamic modelling for fully coupled thermoelastic problems. Contin. Mech. Thermodyn. 2018, 31, 907–937. [Google Scholar] [CrossRef] [Green Version]
  26. Wang, H.; Oterkus, E.; Oterkus, S. Predicting fracture evolution during lithiation process using peridynamics. Eng. Fract. Mech. 2018, 192, 176–191. [Google Scholar] [CrossRef] [Green Version]
  27. Alpay, S.; Madenci, E. Crack growth prediction in fully-coupled thermal and deformation fields using peridynamic theory. In Proceedings of the 54th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Boston, MA, USA, 8–11 April 2013; p. 1477. [Google Scholar]
  28. Diyaroglu, C.; Oterkus, E.; Oterkus, S. An Euler–Bernoulli beam formulation in an ordinary state-based peridynamic framework. Math. Mech. Solids 2017, 24, 361–376. [Google Scholar] [CrossRef]
  29. Zhao, J.; Chen, Z.; Mehrmashhadi, J.; Bobaru, F. A stochastic multiscale peridynamic model for corrosion-induced fracture in reinforced concrete. Eng. Fract. Mech. 2020, 229, 106969. [Google Scholar] [CrossRef] [Green Version]
  30. De Meo, D.; Diyaroglu, C.; Zhu, N.; Oterkus, E.; Siddiq, M.A. Modelling of stress-corrosion cracking by using peridynamics. Int. J. Hydrogen Energy 2016, 41, 6593–6609. [Google Scholar] [CrossRef] [Green Version]
  31. Ghajari, M.; Iannucci, L.; Curtis, P. A peridynamic material model for the analysis of dynamic crack propagation in orthotropic media. Comput. Methods Appl. Mech. Eng. 2014, 276, 431–452. [Google Scholar] [CrossRef] [Green Version]
  32. Cheng, Z.; Zhang, G.; Wang, Y.; Bobaru, F. A peridynamic model for dynamic fracture in functionally graded materials. Compos. Struct. 2015, 133, 529–546. [Google Scholar] [CrossRef]
  33. Oterkus, E.; Guven, I.; Madenci, E. Impact damage assessment by using peridynamic theory. Open Eng. 2012, 2, 523–531. [Google Scholar] [CrossRef]
  34. Liu, X.; He, X.; Wang, J.; Sun, L.G.; Oterkus, E. An ordinary state-based peridynamic model for the fracture of zigzag graphene sheets. Proc. R. Soc. A Math. Phys. Eng. Sci. 2018, 474, 20180019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Madenci, E.; Oterkus, S. Ordinary state-based peridynamics for plastic deformation according to von Mises yield criteria with isotropic hardening. J. Mech. Phys. Solids 2016, 86, 192–219. [Google Scholar] [CrossRef] [Green Version]
  36. Weckner, O.; Mohamed, N.A.N. Viscoelastic material models in peridynamics. Appl. Math. Comput. 2013, 219, 6039–6043. [Google Scholar] [CrossRef]
  37. Foster, J.T.; Silling, S.A.; Chen, W.W. Viscoplasticity using peridynamics. Int. J. Numer. Methods Eng. 2009, 81, 1242–1258. [Google Scholar] [CrossRef]
  38. Oterkus, S.; Madenci, E.; Agwai, A. Peridynamic thermal diffusion. J. Comput. Phys. 2014, 265, 71–96. [Google Scholar] [CrossRef] [Green Version]
  39. Diyaroglu, C.; Oterkus, S.; Oterkus, E.; Madenci, E.; Han, S.; Hwang, Y. Peridynamic wetness approach for moisture concentration analysis in electronic packages. Microelectron. Reliab. 2017, 70, 103–111. [Google Scholar] [CrossRef] [Green Version]
  40. Ouchi, H.; Katiyar, A.; York, J.; Foster, J.; Sharma, M.M. A fully coupled porous flow and geomechanics model for fluid driven cracks: A peridynamics approach. Comput. Mech. 2015, 55, 561–576. [Google Scholar] [CrossRef]
  41. Taylor, M.; Steigmann, D. A two-dimensional peridynamic model for thin plates. Math. Mech. Solids 2013, 20, 998–1010. [Google Scholar] [CrossRef]
  42. O’Grady, J.; Foster, J. Peridynamic beams: A non-ordinary, state-based model. Int. J. Solids Struct. 2014, 51, 3177–3183. [Google Scholar] [CrossRef] [Green Version]
  43. O’Grady, J.; Foster, J. Peridynamic plates and flat shells: A non-ordinary, state-based model. Int. J. Solids Struct. 2014, 51, 4572–4579. [Google Scholar] [CrossRef] [Green Version]
  44. Diyaroglu, C.; Oterkus, E.; Oterkus, S.; Madenci, E. Peridynamics for bending of beams and plates with transverse shear deformation. Int. J. Solids Struct. 2015, 69, 152–168. [Google Scholar] [CrossRef] [Green Version]
  45. Chowdhury, S.R.; Roy, P.; Roy, D.; Reddy, J.N. A peridynamic theory for linear elastic shells. Int. J. Solids Struct. 2016, 84, 110–132. [Google Scholar] [CrossRef]
  46. Madenci, E.; Oterkus, E. Peridynamic Theory and Its Applications; Springer: New York, NY, USA, 2014; Volume 17. [Google Scholar]
  47. Javili, A.; Morasata, R.; Oterkus, E.; Oterkus, S. Peridynamics review. Math. Mech. Solids 2018, 24, 3714–3739. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The positive set of the degrees-of-freedom for Mindlin plate formulation.
Figure 1. The positive set of the degrees-of-freedom for Mindlin plate formulation.
Jcs 04 00076 g001
Figure 2. Peridynamics (PD) force density vector and horizon [46].
Figure 2. Peridynamics (PD) force density vector and horizon [46].
Jcs 04 00076 g002
Figure 3. Variation of material properties through the thickness direction for the finite element model.
Figure 3. Variation of material properties through the thickness direction for the finite element model.
Jcs 04 00076 g003
Figure 4. Simply supported functionally graded Mindlin plate.
Figure 4. Simply supported functionally graded Mindlin plate.
Jcs 04 00076 g004
Figure 5. Numerical discretization, loading and fictitious region.
Figure 5. Numerical discretization, loading and fictitious region.
Jcs 04 00076 g005
Figure 6. Comparison of PD and FEA results along the central x- and y-axes.
Figure 6. Comparison of PD and FEA results along the central x- and y-axes.
Jcs 04 00076 g006
Figure 7. Fully clamped functionally graded Mindlin plate.
Figure 7. Fully clamped functionally graded Mindlin plate.
Jcs 04 00076 g007
Figure 8. Comparison of PD and FEA results along the central x- and y-axes.
Figure 8. Comparison of PD and FEA results along the central x- and y-axes.
Jcs 04 00076 g008aJcs 04 00076 g008b
Figure 9. Functionally graded Mindlin plate subjected to mixed boundary conditions.
Figure 9. Functionally graded Mindlin plate subjected to mixed boundary conditions.
Jcs 04 00076 g009
Figure 10. Comparison of PD and FEA results along the central x- and y-axes.
Figure 10. Comparison of PD and FEA results along the central x- and y-axes.
Jcs 04 00076 g010

Share and Cite

MDPI and ACS Style

Yang, Z.; Oterkus, E.; Oterkus, S. Peridynamic Mindlin Plate Formulation for Functionally Graded Materials. J. Compos. Sci. 2020, 4, 76. https://doi.org/10.3390/jcs4020076

AMA Style

Yang Z, Oterkus E, Oterkus S. Peridynamic Mindlin Plate Formulation for Functionally Graded Materials. Journal of Composites Science. 2020; 4(2):76. https://doi.org/10.3390/jcs4020076

Chicago/Turabian Style

Yang, Zhenghao, Erkan Oterkus, and Selda Oterkus. 2020. "Peridynamic Mindlin Plate Formulation for Functionally Graded Materials" Journal of Composites Science 4, no. 2: 76. https://doi.org/10.3390/jcs4020076

APA Style

Yang, Z., Oterkus, E., & Oterkus, S. (2020). Peridynamic Mindlin Plate Formulation for Functionally Graded Materials. Journal of Composites Science, 4(2), 76. https://doi.org/10.3390/jcs4020076

Article Metrics

Back to TopTop