Research Article  Open Access
The Application of MeshFree Method in the Numerical Simulation of Beams with the Size Effect
Abstract
The meshfree method is employed to implement the numerical simulation of the bending behavior of beams with the size effect. On the basis of the classical BernoulliEuler beam theory, two higherorder strain components are involved in the beam model. The intrinsic bulk length and the directional surface length components are introduced into the constitutive relationship to describe the size effect, and the variation of the total potential is provided. The movingleast square approximation is used to construct the shape function and its second and thirdorder derivatives, and the choice of the scaling factor is discussed in detail. A meshfree scheme is built to implement numerical simulation, in which the higherorder strains are directly approximated with the nodal components due to the higherorder continuity of the shape function. The convergence of method is illustrated in virtue of an example of the simply supported beam, and the effect of the intrinsic bulk length and the directional surface length components are studied.
1. Introduction
The size effect has been found in the areas of micromechanics and nanomechanics, for example thin films, microelectromechanical systems, and nanoelectromechanical systems [1–3]. The scale effect can be successfully modeled by employing the higherorder continuum theory, in which the constitutive equations introduce some material length scale parameters, in addition to the classical material parameters [4–13]. Mostly the generally known higherorder theory includes the nonlocal theory of Eringen [14], strain gradient elasticity, and the couple stress theory [15]. The beam model is very efficient for studying static and dynamic behaviors of slender nanostructures such as carbon nanotubes and microtubules. Zhang et al. studied the small length effect on a nonlocal cantilever beam and carried out a detailed analysis of bending, buckling, and vibration of nanosized beams [4]. Lim gives detailed statements for static response and bending moments of nanosized beams via the nonlocal theory [5, 6]. Gao and Lei study the small scale effects on mechanical behaviors of protein microtubules based on the nonlocal elasticity theory [7]. The classical strain gradient elasticity theory and its modified version [8] have been applied for analysis of bending, vibration, and buckling of microsized beams under various boundary conditions by some researchers [9–11]. A microstructuredependent model for the Timoshenko beam has been developed by using a modified couple stress theory and Hamilton’s principle to account for the microstructural effect [12]. Kahrobaiyan et al. has also investigated the sizedependent dynamic characteristics of atomic force microscopic cantilevers based on the modified couple stress theory [13].
However, in all above studies, the microscale structural information is not incorporated, and the scale parameters need to be assumed to study their effect on the nanostructures’ behavior. On the other hand, most of studies are mainly on the theoretical analysis for some simple support and boundary conditions, and the numerical simulation is very rare due to the complexity from the introduction of the higherorder strain. Specifically, the inclusion of the higherorder gradient largely complicates the solving of the boundary value problems. For example, in the finite element method, the interpolation requires continuity [16, 17], and this leads to difficulty in the establishment of elements and the construction of the interpolation functions. Recently, researchers [17] have applied the meshfree method [16–21] to simulate the materials with straingradients effects. The meshfree method is a newly developed technique that has some distinct advantages. In particular, the movingleast square (MLS) approximations possess nonlocal properties and automatically satisfy the higherorder continuity requirement. This intrinsic nonlocal property leads to real rotationfree approximation [17], and displacements can thus be used as the only nodal freedoms. Motivated by the above fact, the present work attempts to employ the meshfree method to build the numerical scheme for beams with size effect. To produce the continuum thirdorder derivatives of shape function, the cubic basis is needed, and the choice of scaling factor is thus needed to be studied. After the convergence of the established method is tested, the effect of the buck length and surface length components is studied and discussed.
2. HigherOrder BernoulliEuler Beam Model
As shown in Figure 1, a beam is considered, in which the axis is the axis of the beam, and the axis is the deflection axis. The elastic line lies on the plane. Following the works of Akgöz and Civalek [1] and Civalek and Demir [2], we treat the strain energy density of the beam as a function of the classical strain and the double strains. In BernoulliEuler principle, the classical strain is defined as in which denotes the displacement of the elastic line, that is, deflection.
The double strains are defined as [3] The classical stress becomes and the double stresses are given by where is the elastic Young’s modulus, and and are the size effect parameters which are the intrinsic bulk length and the directional surface length, respectively.
The principle of virtual work requires in which is the work of the external forces, is the strain energy of the beam, and is the total potential. Moreover, the variation of the strain energy of the beam is defined by [3] in which denotes the length of beam, and denotes the area of the crosssection.
In addition, the variation of the work of the external forces is given by where denotes the continuously distributed loading, the shear forces at the ends of the beam, and and the bending moments and double moments of the beam.
Substituting (2)–(4) into (6), and calculating the inner integral, we can get where is the moment of inertia of the crosssection with respect to the axis.
3. The MLS Shape Functions and Their Derivatives
In the present research, the MLS approximation is used to build the numerical discretization scheme. Since the MLS approximation satisfies the higherorder continuity automatically, the components and in (8) can directly be approximated with the nodal components. The calculation of the MLS shape functions and their derivatives are illustrated as follows [18].
The deflect defined on the segment is set as where is the monomial basis function, is the number of terms in the basis, and are the coefficients of the basis function. To construct the thirdorder derivative, the following cubic basis function is used: The unknown coefficients in (9) can be determined by the minimization of the weighted discrete norm where is the weight function with compact support, and NP is the number of nodes with . The minimum of in (11) with respect to leads to a set of linear equations: where , and Thus, the unknown coefficients can be obtained from (12) as Substituting (14) into (9), we can get where the MLS shape function is defined as If we set the partial derivatives of with respect with can be obtained asin whichwith
In the present work, the weight function is chosen as where is the support size for the node, and is the cubic spline function [18]:
4. The MeshFree Computational Scheme
In our meshfree simulation, some nodes are allocated uniformly along the beam axis. To implement the numerical integrals, the interval is divided into a series of subdomains that are called integral cells in some papers [18–21]. For the calculation convenience, the integral cells are chosen to coincide with the nodal arrangements; that is, the integral cells are used if the nodes are allocated.
The second and third derivatives of the deflection can be approximated as Their institution into (8) gives the stiffness matrix in which with is the total number of nodes.
The Gaussian quadrature rule is used for each cell. In the practical calculation, the Gaussian points can be globally marked. To avoid the repeating calculation for each Gaussian point, the information that includes the nodes whose supports cover the Gaussian points, the corresponding shape functions, and their derivatives is calculated and deposited in advance.
The solution can be found by solving the equation where is the total nodal component vector, and F is the nodal force vector. It is noted that the unknown vector does not involve the rotation angle; the essential boundary condition can be imposed with the penalty function method [18, 20, 21].
5. Numerical Results
The efficiency of the proposed method is first tested. A thin beam of length is supported at two points and with an applied vertical at its midpoint (Figure 2). The dimensions are fixed asthe height , and thickness . If both and are chosen as zero, the deflection at the midpoint is according to the classical beam theory. The support size of nodes can be expressed as with denoting the scaling factor. Three point Gauss integral formula is used for each integral cell. Table 1 shows the midpoint deflection while the different numbers of nodes are used in the case of being fixed as 3.0. It can be seen that a stable convergence is achieved, and very good result can be obtained when the node number is larger than 15 (the exact solution is 0.25). Table 2 gives the results while the varying scaling factor is used if the node number is fixed as 61, which shows that has a large effect on the numerical results when the cubic basis is used, and a good choice for should be 3.0. For comparison, the results in the case of quadratic basis are also provided in Table 2. It should be pointed out that and are not necessary since and g are fixed as zero in these tests.


Next, the size effect is studied. The total 61 nodes are used, and is chosen as 3.0. The cross section is fixed as , and the boundary conditions are same with the above tests. is first fixed as 0.00125 to test the effect of on the midpoint deflection. Figure 3 plots the relatively increasing amount (%) in comparison with the results of those from the classical beam theory (). An increasing effect on the results can be found with the increasing g. To test the effect of on the deflection, g is fixed as 0.1875b. Figure 4 plots the result comparison with the classical beam theory, which shows that an obvious effect is only associated with the value of being in the certain range.
6. Conclusions
The higherorder continuum theory is used to study the bending behavior of beams, which can reflect the size effect associated with some special structures. The derived formulas provide a foundation for the theoretical analysis and numerical simulation, and the current work is concerned mainly with the numerical simulation. A meshfree scheme has been built to implement the numerical simulation. Specifically, the second and thirdorder derivatives of deflection are directly approximated with the nodal components, which largely make the discretization easy. It is proven that the MLS approximation and meshfree method are very efficient for the problems in which the higherorder strains need to be involved.
The current research employs the cubic basis to produce the continuum thirdorder derivatives of shape function, and it is found that the better convergence can be obtained, the stricter the scaling factor needs to be used in comparison with the quadratic basis. For the beams with size effect, the size effect becomes larger with the increasing buck length component if the surface length component is fixed. But, only the value of surface length component in certain range has a large effect if the buck length component is fixed.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work in this research was supported by the Natural Science Foundation of China (Grant no. 10902129), the Program for New Century Excellent Talents in University, China (Project no. NCET120700), and Program for Science and Technology Innovation Talents in Universities of Henan Province (Project no. 2012HASTIT021).
References
 B. Akgöz and Ö. Civalek, “Strain gradient elasticity and modified couple stress models for buckling analysis of axially loaded microscaled beams,” International Journal of Engineering Science, vol. 49, no. 11, pp. 1268–1280, 2011. View at: Publisher Site  Google Scholar
 Ö. Civalek and C. Demir, “Buckling and bending analyses of cantilever carbon nanotubes using the eulerbernoulli beam theory based on nonlocal continuum model,” Asian Journal of Civil Engineering, vol. 12, no. 5, pp. 651–661, 2011. View at: Google Scholar
 K. A. Lazopoulos and A. K. Lazopoulos, “Bending and buckling of thin strain gradient elastic beams,” European Journal of Mechanics, A/Solids, vol. 29, no. 5, pp. 837–843, 2010. View at: Publisher Site  Google Scholar
 Y. Y. Zhang, C. M. Wang, and N. Challamel, “Bending, buckling, and vibration of micro/nanobeams by hybrid nonlocal beam model,” Journal of Engineering Mechanics, vol. 136, no. 5, pp. 562–574, 2010. View at: Publisher Site  Google Scholar
 C. W. Lim, “Equilibrium and static deflection for bending of a nonlocal nanobeam,” Advances in Vibration Engineering, vol. 8, no. 4, pp. 277–300, 2009. View at: Google Scholar
 C. W. Lim, “On the truth of nanoscale for nanobeams based on nonlocal elastic stress field theory: equilibrium, governing equation and static deflection,” Applied Mathematics and Mechanics, vol. 31, no. 1, pp. 37–54, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Gao and F. M. Lei, “Small scale effects on the mechanical behaviors of protein microtubules based on the nonlocal elasticity theory,” Biochemical and Biophysical Research Communications, vol. 387, no. 3, pp. 467–471, 2009. View at: Publisher Site  Google Scholar
 D. C. C. Lam, F. Yang, A. C. M. Chong, J. Wang, and P. Tong, “Experiments and theory in strain gradient elasticity,” Journal of the Mechanics and Physics of Solids, vol. 51, no. 8, pp. 1477–1508, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Artan and R. C. Batra, “Free vibrations of a strain gradient beam by the method of initial values,” Acta Mechanica, vol. 223, pp. 2393–2409, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Q. Wang and K. M. Liew, “Application of nonlocal continuum mechanics to static analysis of micro and nanostructures,” Physics Letters A: General, Atomic and Solid State Physics, vol. 363, no. 3, pp. 236–242, 2007. View at: Publisher Site  Google Scholar
 Y. Z. Sun and K. M. Liew, “Effect of higherorder deformation gradients on buckling of singlewalled carbon nanotubes,” Composite Structures, vol. 109, pp. 279–285, 2014. View at: Publisher Site  Google Scholar
 H. M. Ma, X. L. Gao, and J. N. Reddy, “A microstructuredependent Timoshenko beam model based on a modified couple stress theory,” Journal of the Mechanics and Physics of Solids, vol. 56, no. 12, pp. 3379–3391, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. H. Kahrobaiyan, M. Asghari, M. Rahaeifard, and M. T. Ahmadian, “Investigation of the sizedependent dynamic characteristics of atomic force microscope microcantilevers based on the modified couple stress theory,” International Journal of Engineering Science, vol. 48, no. 12, pp. 1985–1994, 2010. View at: Publisher Site  Google Scholar
 A. C. Eringen, “On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves,” Journal of Applied Physics, vol. 54, no. 9, pp. 4703–4710, 1983. View at: Publisher Site  Google Scholar
 R. D. Mindlin and N. N. Eshel, “On first straingradient theories in linear elasticity,” International Journal of Solids and Structures, vol. 4, no. 1, pp. 109–124, 1968. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. C. Xia and J. W. Hutchinson, “Crack tip fields in strain gradient plasticity,” Journal of the Mechanics and Physics of Solids, vol. 44, no. 10, pp. 1621–1648, 1996. View at: Publisher Site  Google Scholar
 Z. Tang, S. Shen, and S. N. Atluri, “Analysis of materials with straingradient effects: a Meshless Local PetrovGalerkin (MLPG) approach, with nodal displacements only,” CMES: Computer Modeling in Engineering and Sciences, vol. 4, no. 1, pp. 177–196, 2003. View at: Google Scholar  Zentralblatt MATH
 T. Belytschko, Y. Krongauz, D. Organ, M. Fleming, and P. Krysl, “Meshless methods: an overview and recent developments,” Computer Methods in Applied Mechanics and Engineering, vol. 139, no. 1–4, pp. 3–37, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Ren, K. M. Liew, and S. A. Meguid, “Modelling and simulation of the superelastic behaviour of shape memory alloys using the elementfree Galerkin method,” International Journal of Mechanical Sciences, vol. 44, no. 12, pp. 2393–2413, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Z. Sun and K. M. Liew, “Bending buckling of singlewalled carbon nanotubes: higher order gradient continuum and meshfree method,” Computer Method in Applied Mechanics and Engineering, vol. 197, pp. 3001–3013, 2008. View at: Google Scholar
 Y. Sun and K. M. Liew, “Application of the higherorder CauchyBorn rule in meshfree continuum and multiscale simulation of carbon nanotubes,” International Journal for Numerical Methods in Engineering, vol. 75, no. 10, pp. 1238–1258, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2014 Yuan Tian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.