Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 590271, 6 pages

http://dx.doi.org/10.1155/2014/590271

Research Article

The Application of Mesh-Free Method in the Numerical Simulation of Beams with the Size Effect

Yuan Tian,¹ Caili Zhang,² and Yuzhou Sun²

¹School of Earth and Space Science, University of Science and Technology of China, Hefei 230026, China

²School of Civil Engineering and Architecture, Zhongyuan University of Technology, Zhengzhou 450007, China

Received 5 January 2014; Accepted 6 March 2014; Published 30 March 2014

Academic Editor: Zan Zhang

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.

Abstract.

The mesh-free method is employed to implement the numerical simulation of the bending behavior of beams with the size effect. On the basis of the classical Bernoulli-Euler beam theory, two higher-order 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 moving-least square approximation is used to construct the shape function and its second- and third-order derivatives, and the choice of the scaling factor is discussed in detail. A mesh-free scheme is built to implement numerical simulation, in which the higher-order strains are directly approximated with the nodal components due to the higher-order 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 higher-order 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 higher-order 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 micro-sized beams under various boundary conditions by some researchers [9–11]. A microstructure-dependent 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 size-dependent 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 higher-order strain. Specifically, the inclusion of the higher-order 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 mesh-free method [16–21] to simulate the materials with strain-gradients effects. The mesh-free method is a newly developed technique that has some distinct advantages. In particular, the moving-least square (MLS) approximations possess nonlocal properties and automatically satisfy the higher-order continuity requirement. This intrinsic nonlocal property leads to real rotation-free 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 mesh-free method to build the numerical scheme for beams with size effect. To produce the continuum third-order 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. Higher-Order Bernoulli-Euler 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 Bernoulli-Euler principle, the classical strain is defined as

in which denotes the displacement of the elastic line, that is, deflection.

Figure 1: A beam model and related markers.

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 cross-section.

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 cross-section 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 higher-order 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 third-order 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 as

in which

with

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 Mesh-Free Computational Scheme

In our mesh-free 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.

Table 1: The midpoint deflection for different numbers of nodes ( is fixed as 3.0).

Node number
5
9
13
17
21
25
29
33
61

0.2716
0.2545
0.2516
0.2508
0.2505
0.2503
0.2502
0.2501
0.2500

Table 2: The midpoint deflection for varying scaling factor (61 nodes are used).

2.0
2.5
3.0
3.5
4.0
4.5
5.0

Cubic basis
0.0564
0.0634
0.2500
0.2502
0.2817
0.2605
0.3544

Quadratic basis
0.2499
0.2501
0.2500
0.2502
0.3382
0.2601
0.2165

Figure 2: A simply supported beam and its deformation under a vertical force .

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.

Figure 3: The relatively increasing amount (%) of the midpoint deflection for the varying value of in comparison with the classical beam theory.

Figure 4: The relatively increasing amount (%) of the midpoint deflection for the varying value of in comparison with the classical beam theory.

6. Conclusions

The higher-order 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 mesh-free scheme has been built to implement the numerical simulation. Specifically, the second- and third-order derivatives of deflection are directly approximated with the nodal components, which largely make the discretization easy. It is proven that the MLS approximation and mesh-free method are very efficient for the problems in which the higher-order strains need to be involved.

The current research employs the cubic basis to produce the continuum third-order 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. NCET-12-0700), 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 micro-scaled beams,” International Journal of Engineering Science, vol. 49, no. 11, pp. 1268–1280, 2011.

Ö. Civalek and C. Demir, “Buckling and bending analyses of cantilever carbon nanotubes using the euler-bernoulli beam theory based on non-local continuum model,” Asian Journal of Civil Engineering, vol. 12, no. 5, pp. 651–661, 2011.

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.

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.

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.

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.

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.

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.

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.

Q. Wang and K. M. Liew, “Application of nonlocal continuum mechanics to static analysis of micro- and nano-structures,” Physics Letters A: General, Atomic and Solid State Physics, vol. 363, no. 3, pp. 236–242, 2007.

Y. Z. Sun and K. M. Liew, “Effect of higher-order deformation gradients on buckling of single-walled carbon nanotubes,” Composite Structures, vol. 109, pp. 279–285, 2014.

H. M. Ma, X. L. Gao, and J. N. Reddy, “A microstructure-dependent 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.

M. H. Kahrobaiyan, M. Asghari, M. Rahaeifard, and M. T. Ahmadian, “Investigation of the size-dependent 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.

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.

R. D. Mindlin and N. N. Eshel, “On first strain-gradient theories in linear elasticity,” International Journal of Solids and Structures, vol. 4, no. 1, pp. 109–124, 1968.

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.

Z. Tang, S. Shen, and S. N. Atluri, “Analysis of materials with strain-gradient effects: a Meshless Local Petrov-Galerkin (MLPG) approach, with nodal displacements only,” CMES: Computer Modeling in Engineering and Sciences, vol. 4, no. 1, pp. 177–196, 2003.

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.

J. Ren, K. M. Liew, and S. A. Meguid, “Modelling and simulation of the superelastic behaviour of shape memory alloys using the element-free Galerkin method,” International Journal of Mechanical Sciences, vol. 44, no. 12, pp. 2393–2413, 2002.

Y. Z. Sun and K. M. Liew, “Bending buckling of single-walled carbon nanotubes: higher order gradient continuum and mesh-free method,” Computer Method in Applied Mechanics and Engineering, vol. 197, pp. 3001–3013, 2008.

Y. Sun and K. M. Liew, “Application of the higher-order Cauchy-Born rule in mesh-free continuum and multiscale simulation of carbon nanotubes,” International Journal for Numerical Methods in Engineering, vol. 75, no. 10, pp. 1238–1258, 2008.