Displacement Approach to Determine the Effective Tensile and Torsional Modulus of Nanowires

Surface elasticity and residual stress have a strong inﬂuence on the eﬀective properties of nanowire (NW) due to its excessively large surface area-to-volume ratio. Here, the classical displacement method is used to solve the ﬁeld equations of the core-surface layer model subjected to tension and torsion. The eﬀective Young’s modulus is deﬁned as the ratio of normal stress to axial strain, which decreases with the increase in NW radius and gradually reaches the bulk value. The positive or negative surface residual stresses will increase or decrease Young’s modulus and shear modulus due to the surface residual strains. Nonzero radial and circumferential strains enhance the inﬂuence of surface moduli on the eﬀective modulus.


Introduction
Typical nanowires (NWs) are often referred to as 1D materials with nanometer-scale diameters or perimeters and excessively large surface area-to-volume ratio. NWs have considerable potential in various applications, such as molecular electronics, nanoelectromechanical systems, and novel building materials, for disaster prevention and mitigation [1][2][3][4][5][6][7][8][9][10][11][12]. e applications of NWs into future generation nanodevices require a complete understanding of the NW mechanical properties [2]. Many direct measurements have been performed to investigate the mechanical properties of NWs [3]. Unlike the mechanical testing of bulk materials, NW testing heavily depends on the experimental setup; in particular, manipulation procedure leads to substantial challenges due to the small NW dimensions [4]. In many experiments, such as in [1][2][3][4][5][6], the measured deflections are less than the diameter which can be classified small-deflection problem. In addition to experimental endeavors, theoretical prediction can also be used in NW mechanical analysis. eoretical prediction is classified into two main categories as follows: first is the atomic modeling, which includes techniques such as tedious ab initio molecular dynamics calculations and density functional model [5], and second one is continuum mechanics modeling [6][7][8][9][10]. NWs are strongly influenced by their surface characteristics, thereby leading to distinct mechanical properties compared with their bulk counterpart. Consequently, in the continuum mechanics modeling of the mechanical properties of a solid NW, the role of surface stress must be considered. Zhang et al. analyzed the effect of surface residual stress and elasticity on the asymmetric yield strength of NWs on the basis of the potential energy method [6]. Chuang presented a simplistic theory to study the enhanced strength of a solid NW [7]. Gupta also presented a continuum formulation to investigate the finite deformation of nanorod/NWs [8]. e large-deflection deformation of NW implicates large rotational angle and infinitesimal strain. So, this situation will not be considered in our mode.
Although the classical continuum mechanics models can efficiently predict NW deformation, their applicability in identifying the surface effect on the effective modulus of NWs is tedious. erefore, a relatively simple approach that can directly characterize the mechanical properties of a solid NW should be developed. In this work, the classical displacement method is used to study the tension and torsion of NW and the influence of surface elasticity on the effective Young's moduli and shear moduli.

Model Analysis
As schematically shown in Figure 1, we consider a stretched and twisted core-shell NW model viewed as a composite, comprising a cylindrical core (bulk) and surface layer. e equilibrium equations, strain-displacement relationships, and constitutive equations for the isotropic bulk materials are expressed as follows: where σ ij , u i , and ε ij denote the stresses, displacements, and strains in core, respectively, and λ and G are the Lamé constants of the bulk. For the elastic isotropic surface layer, the linear relationship of the surface stress and elastic strain can be expressed by Gurtin-Murdoch elasticity as follows [9]: where λ s and G s are the surface modulus, τ 0 is the surface residual stress, and δ αβ denotes the surface Christoffel symbols.
According to the analysis of the mechanical equilibrium of surface layer and core, the generalized Young-Laplace equations for NW are expressed as follows [6,10]: where n j is the unit normal vector and κ s αβ is the curvature tensor of the surface ( Figure 2). e stress distribution of NW can be obtained using the classical displacement method on the basis of equations (1)∼(6), that is, the bulk stress and the surface stress can be obtained by directly substituting the bulk displacement/strain formula into the constitutive equations. e global displacements dominant governing equations are established using generalized Young-Laplace equations. e solution of the field equations in the cylinder coordinate can be simplified if the NW is subjected to tension and torsion. We assume that u z is function of z and u ρ is function of ρ. According to equation (1), the nonzero strains are After substituting equation (7) into equation (2), the stress components in the bulk are After substituting equations (8)-(11) into equation (3), it is noted that σ ρ � σ φ , and the stress components are found to be automatically satisfied the radial equilibrium equation. And, we get ε ρ � ε φ and u ρ is the linear distribution along the radial direction. e axial and circumferential equilibrium equation can be simplified into the following forms: which implies that u z and u φ are linear distribution along the axial direction of NW. Now, the displacement field can be expressed by where a 1 ∼ a 3 and b 1 ∼ b 2 are constants related to boundary conditions. e surface of NW is assumed to be characterized by the deformation of the bulk solid so that surface strains can be expressed by where R 0 is the radius of NW. Substituting equation (17) into equation (4), the surface stresses are After substituting equations (18)-(20) into equations (5) and (6), the above surface stress components automatically satisfy the axial and circumferential generalized 2 Advances in Civil Engineering Young-Laplace equations.
e generalized Young-Laplace equation along the radial direction gives σ ρ � τ s φφ /R 0 , and it is Under a constant external load P, the average normal stress on the cross section of NW is where λ p � λ s /R 0 and G p � G s /R 0 . Under a constant external torsion T, the shear stress on the cross section of NW satisfy the following equation: where I p � πR 4 0 /32 and I s � 2πR 3 0 . Combining equations (17), (22), and (23), the constants can be confirmed as follows: We next determine the distribution of the surface residual strain and then consider the effective modulus of NW. If nonzero surface residual stress is present on the NW surface and let σ � 0 in equations (25) and (26), then the relaxed surface axial strain and circumferential strain are expressed as follows: Advances in Civil Engineering e surface residual stresses of the NW are inherent and in the self-equilibrium state, that is, independent of the external load. If the NW is subjected to tension and torsion, then the strain and stress analysis of NW can be also performed using equations (7)∼(26). Let u z0 and u zl be the axial displacements of the two NW ends and l 0 be the initial length of NW; the average axial strain is e effective strains of NW are the difference between average axial strain and surface radial strain, as expressed below: (29) Consequently, the effective Young's modulus of NW is obtained as follows: e effective Poisson's ratio of NW is e effective shear modulus of NW is (32)

Results and Discussion
All the displacements, strains, bulk stresses, and surface stresses of the NW have been determined using the classical displacement approach. Figure 3 shows the variation in the effective Young's moduli of the NW compared with its radius. For example, in Al NW, the surface moduli are λ s � 6.8415 N/m and G s � −0.3755 N/m [6]. e bulk parameters are λ � 59.2 GPa and G � 25.4 GPa [6]. e solid curve on the basis of our continuum formula (equation (16)) matches Zhang et al. [6] approach (the dotted curve). e effect of surface stresses on the effective Young's modulus of NW is illustrated, where Young's modulus decreases with the increase in the NW radius (R 0 ) and gradually reaches a constant value of 69 GPa. As shown in Figure 3, the two approaches provide slightly different results with small radius (<5 nm). Our formula offers small effective Young' modulus and rapid reduction because the potential energy method adds to the influence of radial and circumferential strains in the definition of Young's modulus. Homogenization theory of nanocomposites can provide a rigorous definition to define effective properties [10]. As depicted in equation (28), the surface residual stresses (τ 0 ) also affects Young's modulus. e surface residual stresses are inherent and in the self-equilibrium state, that is, independent of the external load, due to the surface residual strain (equation [10]). Hence, the positive or negative surface residual stresses will increase (or decrease) Young's modulus of NW.
Surface moduli also have strong influence on Young's modulus. e absolute value of the NW surface elasticity is generally <1000 N/m, but it is difficult to accurately quantify. A minimal difference is observed between the experimental approach and numerical atomistic analysis even with small radius. For reference, we also present the variation of the effective modulus of NW compared with its surface moduli in Figure 4. e NW with a large radius is considered (R � 10 nm, 20 nm, 30 nm, 40 nm) solely for computational purpose. Figure 4 illustrates the increase in Young's modulus with the increase in the surface moduli (λ s ). Equations (8)- (11) show that the surface shear strains are zero, but the surface area expansion is nonzero. Equation (4) implies high surface stress values with high surface moduli, thereby resulting in large NW Young's modulus. However, the amplification of the surface moduli on Young's modulus is controlled by the NW radius, that is, the surface area-tovolume ratio increases with the decrease in NW radius. is result suggested that the effective modulus of NW is enhanced by the surface stress. Size dependence is the general characteristics of nanomaterials. Classical displacement approach can provide the surface strain and stress and simplifies the analysis of the effect of surface elasticity on Young's modulus.

Conclusions
NW can be viewed as a composite structure, with the inner core having the normal properties and the surface layer having surface elasticity on the basis of the Gurtin-Murdoch elasticity. e generalized Young-Laplace equations for NW are required in addition to the field equations for core and surface layer. Classical displacement approach in the cylinder coordinate has been used to determine the stress distribution of NW. e nominal normal stress and axial strain are defined in the initial NW configuration with initial length (l 0 ) and radius (R 0 ). e l 0 /R 0 ratio shows that the effective Young's modulus decreases with the increase in R 0 of the nanowire and gradually reaches the bulk value. e positive or negative surface residual stresses will increase (or decrease) Young's modulus of NW due to the surface residual strain. Surface moduli also have strong influence on the effective Young's modulus. Nonzero radial and circumferential strains lead to nonzero area expansion, which enhances the influence of surface moduli effective Young's modulus. e radial displacement in the NW is finite [8], and the present model can be easily extended to analyze the effective properties if true strain is used. NW torsion and bending may also be modeled if we consider a proper assumption of the displacement distribution in equation (7). In conclusion, classical displacement approach can obtain NW displacements, strains, and stress distributions and as well as its effective properties. e mechanism underlying the influence of the general characteristics of NW on Young's modulus can also be easily considered.

Data Availability
e cited data, about surface elastic constants of silver nanowire, used to support the findings of this study are included within the article. e data have been used to verify our theoretical prediction.  Advances in Civil Engineering 5