纳米压痕过程的三维有限元数值试验研究

采用有限元方法模拟了纳米压痕仪的加、卸载过程，三维有限元模型考虑了纳米压痕仪的标准Berkovich压头．介绍了有限元模型的几何参数、边界条件、材料特性与加载方式，讨论了摩擦、滑动机制、试件模型的大小对计算结果的影响，进行了计算结果与标准试样实验结果的比较，证实了模拟的可靠性．在此基础上，重点研究了压头尖端曲率半径对纳米压痕实验数据的影响．对比分析了尖端曲率半径r＝0与r＝100nm两种压头的材料压痕载荷—位移曲线．结果表明，当压头尖端曲率半径r≠0时，基于经典的均匀连续介质力学本构理论、传统的实验手段与数据处理方法，压痕硬度值会随着压痕深度的减小而升高．     

Influence of Vickers tip imperfection on depth sensing indentation tests

In the last decade, the development of depth sensing indentation equipment has permitted two of the most conventional mechanical properties of materials to be easily determined: hardness and Young’s modulus. Some care is needed to accurate results from the experimental determination of the aforementioned mechanical properties. In this study, numerical simulations were performed on two well-known materials (Bk7 glass, AISI M2 steel) and on fictitious materials with a wide range of mechanical properties, using Vickers indenters with different sizes of tip imperfections. The purpose is to estimate the influence of the size of the defect on the hardness and Young’s modulus results obtained by ultramicro and nanoindentation tests.     

Comparison between Berkovich,Vickers and conical indentation tests: A three-dimensional numerical simulation study

Three-dimensional numerical simulations of Berkovich, Vickers and conical indenter hardness tests were carried out to investigate the influence of indenter geometry on indentation test results of bulk and composite film/substrate materials. The strain distributions obtained from the three indenters tested were studied, in order to clarify the differences in the load–indentation depth curves and hardness values of both types of materials. For bulk materials, the differentiation between the results obtained with the three indenters is material sensitive. The indenter geometry shape factor, β, for evaluating Young’s modulus for each indenter, was also estimated. Depending on the indenter geometry, distinct mechanical behaviours are observed for composite materials, which are related to the size of the indentation region in the film. The indentation depth at which the substrate starts to deform plastically is sensitive to indenter geometry.     

Two new expanding cavity models for indentation deformations of elastic strain-hardening materials

Two expanding cavity models (ECMs) are developed for describing indentation deformations of elastic power-law hardening and elastic linear-hardening materials. The derivations are based on two elastic–plastic solutions for internally pressurized thick-walled spherical shells of strain-hardening materials. Closed-form formulas are provided for both conical and spherical indentations, which explicitly show that for a given indenter geometry indentation hardness depends on Young’s modulus, yield stress and strain-hardening index of the indented material. The two new models reduce to Johnson’s ECM for elastic-perfectly plastic materials when the strain-hardening effect is not considered. The sample numerical results obtained using the two newly developed models reveal that the indentation hardness increases with the Young’s modulus and strain-hardening level of the indented material. For conical indentations the values of the indentation hardness are found to depend on the sharpness of the indenter: the sharper the indenter, the larger the hardness. For spherical indentations it is shown that the hardness is significantly affected by the strain-hardening level when the indented material is stiff (i.e., with a large ratio of Young’s modulus to yield stress) and/or the indentation depth is large. When the indentation depth is small such that little or no plastic deformation is induced by the spherical indenter, the hardness appears to be independent of the strain-hardening level. These predicted trends for spherical indentations are in fairly good agreement with the recent finite element results of Park and Pharr.     

Characterization of the Grain-Level Mechanical Behavior of Eglin Sand by Nanoindentation

Characterization was made on the structure and grain-level mechanical behavior of Eglin sand (Quikrete #1961 sand quarried in Pensacola, FL). The as-received assorted sand was sorted to six grain sizes: 0.60 mm, 0.50 mm, 0.42 mm, 0.30 mm, 0.212 mm, and 0.15 mm. The sand chemical constituents and crystalline structures were determined using energy dispersive X-ray spectroscopy, X-ray diffraction and transmission electron microscopy. The Young’s modulus and hardness were determined using nanoindentation with a Berkovich tip, and the fracture toughness was measured using a cube-corner tip. The median Young’s modulus, hardness and fracture toughness were determined as 90.4 GPa, 12.8 GPa and 2.32 MPa?m^0.5, respectively. The mechanical properties were analyzed statistically and the parameters of the Weibull distribution were determined. The grains show highly ductile behavior under nanoindentation due to confinement by high pressure induced by Berkovich tip. An inverse problem solving approach using finite element method (FEM) with the consideration of the Ramberg-Osgood model was used to determine the stress–strain relationship for individual sand grains.     

Nonlinearly Viscoelastic Nanoindentation of PMMA Under a Spherical Tip

The Oliver-Pharr method has been well established to measure Young’s modulus and hardness of materials without time-dependent behavior in nanoindentation. The method, however, is not appropriate for measuring the viscoelastic properties of materials with pronounced viscoelastic effects. One well-known phenomenon is the formation of unloading “nose” or negative stiffness during unloading that often occurs during slow loading-unloading in nanoindentation on a viscoelastic material. Most methods in literature have only considered the loading curve for analysis of viscoelastic nanoindentation data while the unloading portion is not analyzed adequately to determine the nonlinearly viscoelastic properties. In this paper, nonlinearly viscoelastic effects are considered and modeled using the nonlinear Burgers model. Nanoindentation was conducted on poly-methylmethacrylate (PMMA) using a spherical indenter tip. An inverse problem solving approach is used to allow the finite element simulation results to agree with the nanoindentation load–displacement curve during the entire loading and unloading stage. This approach has allowed the determination of the nonlinearly viscoelastic behavior of PMMA at submicron scale. In addition, the nanoindentation unloading “nose” has been captured by simulation, indicating that the negative stiffness in the viscoelastic material is the result of memory effect in time-dependent materials.     

An expanding cavity model incorporating strain-hardening and indentation size effects

An expanding cavity model (ECM) for determining indentation hardness of elastic strain-hardening plastic materials is developed. The derivation is based on a strain gradient plasticity solution for an internally pressurized thick-walled spherical shell of an elastic power-law hardening material. Closed-form formulas are provided for both conical and spherical indentations. The indentation radius enters these formulas with its own dimensional identity, unlike that in classical plasticity based ECMs where indentation geometrical parameters appear only in non-dimensional forms. As a result, the newly developed ECM can capture the indentation size effect. The formulas explicitly show that indentation hardness depends on Young’s modulus, yield stress, strain-hardening exponent and strain gradient coefficient of the indented material as well as on the geometry of the indenter. The new model reduces to existing classical plasticity based ECMs (including Johnson’s ECM for elastic–perfectly plastic materials) when the strain gradient effect is not considered. The numerical results obtained using the newly developed model reveal that the hardness is indeed indentation size dependent when the indentation radius is very small: the smaller the indentation, the larger the hardness. Also, the indentation hardness is seen to increase with the Young’s modulus and strain-hardening level of the indented material for both conical and spherical indentations. The strain-hardening effect on the hardness is observed to be significant for materials having strong strain-hardening characteristics. In addition, it is found that the indentation hardness increases with decreasing cone angle of the conical indenter or decreasing radius of the spherical indenter. These trends agree with existing experimental observations and model predictions.     

Measurement of Temperature-Dependent Young’s Modulus at a Strain Rate for a Molding Compound by Nanoindentation

The mechanical properties of a molding compound on a packaged integrated circuit (IC) were measured by spherical nanoindentation using a 50 μm radius diamond tip. The molding compound is a heterogeneous material, consisting of assorted diameters of glass beads embedded in an epoxy. Statistical analysis was conducted to determine the representative volume element (RVE) size for a nanoindentation grid. Nanoindentation was made on the RVE to determine the effective viscoelastic properties. The relaxation functions were converted to temperature-dependent Young’s modulus at a given strain rate at several elevated temperatures. The Young’s modulus values at a given strain rate from nanoindentation were found to be in a good agreement with the corresponding data obtained from tensile samples at or below 90 °C. However, the values from nanoindentation were significantly lower than the data obtained from tensile samples when the temperature was near or higher than 110 °C, which is near the glass transition. The spatial distribution of the Young’s modulus at a given strain rate was determined using nanoindentation with a Berkovich tip. The spatial variation of the Young’s modulus at a given strain rate is due to the difference in nanoindentation sites (glass beads, epoxy or the interphase region). A graphical map made from an optical micrograph agrees reasonably well with the nanoindentation results.     

A procedure for estimating Young’s modulus of textured polycrystalline materials

In this study, a procedure for estimating Young’s modulus of textured and non-textured polycrystalline materials was examined based on finite element analyses, which were performed using three-dimensional polycrystalline finite element models of a random structure, generated using the Voronoi tessellation. Firstly, the local stress/strain distribution and its influence on macroscopic elastic properties were evaluated. Then, the statistical relationship between Young’s modulus obtained from the finite element analyses and averaged Young’s modulus of all grains evaluated based on Voigt’s or Reuss’ model was investigated. It was revealed that the local stress/strain in the polycrystalline body is affected by crystal orientation and deformation constraint caused by adjacent grains, whereas only the crystal orientation affects Young’s modulus of the polycrystalline body when the number of grains is large enough. It was also shown that Young’s modulus correlates well with the averaged Young’s modulus of all grains, in which the size of grains is considered in the averaging. Finally, a procedure for estimating Young’s modulus of textured and non-textured materials was proposed. Young’s modulus of various materials can be estimated from the elastic constants of single crystal and the distribution of crystal orientation and size of grains, which can be obtained by using electron backscatter diffraction (EBSD).     

考虑压头曲率半径的J2形变理论压痕有限元分析

采用经典J2形变理论,在考虑压头曲率半径的前提下,针对几组微压痕实验进行了有限元数值模拟,给出理论计算结果并与实验进行了系统比较．结果发现考虑了压头曲率半径后的经典J2形变理论得到的整段计算曲线不能和实验曲线吻合,可以推理出即使考虑压头曲率半径的的影响,经典J2形变理论也不能很好地解释实验现象．     

Determination of Mechanical Properties of Sand Grains by Nanoindentation

Determination of the mechanical properties of individual sand grains by conventional material testing methods at the macroscale is somewhat difficult due to the sizes of the individual sand particles (a few μm to mm). In this paper, we used the nanoindentation technique with a Berkovich tip to measure the Young’s modulus, hardness, and fracture toughness. An inverse problem solving approach was adopted to determine the stress-strain relationship of sand at the granular level using the finite element method. A cube-corner indenter tip was used to generate radial cracks, the lengths of which were used to determine the fracture toughness. Scatter in the data was observed, as is common with most brittle materials. In order to consider the overall mechanical behavior of the sand grains, statistical analysis of the mechanical properties data (including the variability in the properties) was conducted using the Weibull distribution function. This data can be used in the mesoscale simulations.     

Elastic modulus measurements of vertically aligned multi walled carbon nanotubes carpets by using the nanoindentation technique

Electrical, thermal and mechanical properties of Vertically Aligned Multi Walled Carbon NanoTubes (VA-MWCNT) make them an ideal candidate to replace some of conventional materials in micro and nano-electronic components. Integrating this material in micro components requires a good knowledge of their properties. As the electrical and thermal properties, the MWCNT mechanical properties are difficult to assess. Several techniques have been developed to estimate the CNT Young's modulus and the obtained results cover a large range of scale. In this study, we propose an indirect technique for MWCNT carpet Young's modulus measurements by using the nanoindentation technique. Nanoindentation tests are performed on a metallic film deposited on MWCNT. The measured equivalent reduced modulus takes into account the elastic properties of the metallic thin film and those of the MWCNT substrate. Bec et al. model, introduced in 2006, is used to separate elastic properties, and thus determine the MWCNT reduced Young’s modulus which is estimated between 329 and 352 GPa. Knowing the indenter mechanical properties, we estimate the Young’s modulus in the 461 to507 GPa range.     

Numerical evaluation of the thermomechanical effective properties of a functionally graded material using the homogenization method

When the stresses of the functionally graded materials (FGMs) are discussed under thermal and/or mechanical loading conditions, the different thermomechanical effective properties are needed. For the steady state thermal analyses, these properties include the Young’s modulus, Poisson’s ratio, thermal expansion coefficient and thermal conductivity. For the transient analyses of the heat conduction problem, on the other hand, the density and heat capacity should be added to the aforementioned properties. The homogenization method (HM) based on the finite element method (FEM) is used as it has advantages, such as it is appropriate for estimating the effective properties of composites with a given periodic fiber distribution and complicated geometries. For a periodic composite structure, it is not necessary to study the whole structure but only a representative volume element (RVE) or a unit cell (UC). As the overall behavior of composites depends on the arrangement of the reinforcements, the corresponding UCs of two different arrangements of the fibers are analyzed; namely the square and hexagonal arrangements. It is found that the square arrangement predicts higher values of the Young’s modulus than the hexagonal one but with small difference. In order to verify the computed values of the properties, the results are compared with previous experimental measurements and results of analytical and numerical methods, and good agreement is achieved.     

血红细胞力学性能的纳米压痕实验和有限元模拟

采用纳米压痕技术和有限元方法研究了血红细胞的生物力学性能. 进行了血红细胞的纳米压痕实验, 得到了血红细胞的材料参数和变形形貌; 在实验基础上, 建立了血红细胞的三维有限元模型, 模拟了血红细胞的压痕载荷-位移曲线, 并考虑了参数效应. 数值模拟结果和实验数据符合很好. 通过改变压头与材料之间的摩擦系数和压头曲率半径等参数, 比较了载荷-位移曲线的变化情况. 研究表明摩擦系数对压痕载荷-位移曲线和应力分布影响很小, 而压头曲率对载荷-位移曲线的影响明显.     

考虑压头曲率半径和应变梯度的微压痕分析

在压头尖端曲率半径取100nm的前提下，采用Chen和Wang的应变梯度理论，对微压痕实验进行了系统的数值分析. 首先通过拟合载荷-位移实验曲线的后半段来确定材料的屈服应力和幂硬化指数值,然后用有限元方法数值模拟压痕实验，并将计算得到的整段载荷-位移曲线及硬度-位移曲线和实验结果进行了比较. 结果表明应变梯度理论所预测的计算结果和实验结果很好地符合,包括压痕深度在亚微米和微米范围内的整段曲线.     

Simulations of nanowire bend tests for extracting mechanical properties

Mechanical properties of nickel nanowires are characterized based on the numerical simulations of bend tests performed with a customized atomic force microscope (AFM) and scanning electron microscope (SEM). Nickel nanowire specimens are subjected to bending loads by the tip of the AFM cantilever. The experimental force versus bending displacement curves are compared against simulations from finite element analysis and peridynamic theory, and the mechanical properties are extracted based on their best correlations. Similarly, SEM images of fractured nanowires are compared against peridynamic failure simulations. The results of this study reveal that nickel nanowires have significantly higher strengths than their bulk counterparts, although their elastic modulus values are comparable to bulk nickel modulus values.     

Local mechanical properties of the 6061-T6 aluminium weld using micro-traction and instrumented indentation

The local mechanical properties of a weld zone, in a 6061-T6 aluminium alloy subjected to the modified indirect electric arc technique have been studied. The mechanical properties of the base metal, the weld metal and the heat affected zone were determined by means of usual and instrumented indentation testing, as well as micro-traction testing. To analyse the heat input effect resulting from the welding process, the evolution of the weld zone size was evaluated by means of classical indentation under a constant applied load. The results were presented using a Vickers hardness map representation. This allows monitoring exact hardness variation while leading to the identification of the different zones of the welded joint. Instrumented indentation testing was carried out to determine the local mechanical properties, such as the yield stress, the bulk modulus and the strain-hardening exponent. Obtained results are compared to those derived from tensile tests conducted on micro-specimen cuts taken from the weld zone. It was observed that yield stress values are directly comparable for indentation and micro-traction experiments. As for the elastic properties, no comparison was possible since the bulk modulus is measured by indentation, whereas it is the Young’s modulus by tensile test. The micro-traction testing seems to be more sensitive to represent the work hardening of a material since the corresponding exponent is found to be constant by instrumented indentation.     

Numerical analysis of effective elastic properties of geomaterials containing voids using 3D random fields and finite elements

The purpose of the study is to investigate the influence of porosity and void size on effective elastic geotechnical engineering properties with a 3D model of random fields and finite element. The random field theory is used to generate models of geomaterials containing spatially random voids with controlled porosity and void size. A “tied freedom” analysis is developed to evaluate the effective Young’s modulus and Poisson’s ratio in an ideal block material of finite elements. To deliver a mean and standard deviation of the elastic parameters, this approach uses Monte-Carlo simulations and finite elements, where each simulation leads to an effective value of the property under investigation. The results are extended to investigate an influence of representative volume element (RVE). A comparison of the effective elastic stiffness of 2D and 3D models is also discussed.     

Analytical and numerical techniques to predict carbon nanotubes properties

In this paper, two different approaches for modeling the behaviour of carbon nanotubes are presented. The first method models carbon nanotubes as an inhomogeneous cylindrical network shell using the asymptotic homogenization method. Explicit formulae are derived representing Young’s and shear moduli of single-walled nanotubes in terms of pertinent material and geometric parameters. As an example, assuming certain values for these parameters, the Young’s modulus was found to be 1.71 TPa, while the shear modulus was 0.32 TPa. The second method is based on finite element models. The inter-atomic interactions due to covalent and non-covalent bonds are replaced by beam and spring elements, respectively, in the structural model. Correlations between classical molecular mechanics and structural mechanics are used to effectively model the physics governing the nanotubes. Finite element models are developed for single-, double- and multi-walled carbon nanotubes. The deformations from the finite element simulations are subsequently used to predict the elastic and shear moduli of the nanotubes. The variation of mechanical properties with tube diameter is investigated for both zig-zag and armchair configurations. Furthermore, the dependence of mechanical properties on the number of nanotubules in multi-walled structures is also examined. Based on the finite element model, the value for the elastic modulus varied from 0.9 to 1.05 TPa for single and 1.32 to 1.58 TPa for double/multi-walled nanotubes. The shear modulus was found to vary from 0.14 to 0.47 TPa for single-walled nanotubes and 0.37 to 0.62 for double/multi-walled nanotubes.     

Extraction of Mechanical Properties with Second Harmonic Detection for Dynamic Nanoindentation Testing

In this article, a new method based on the detection of the second harmonic of the displacement signal to determine mechanical properties of materials from dynamic nanoindentation testing, is presented. With this technique, the Young’s modulus and hardness of homogeneous materials can be obtained at small penetration depths from the measurement of the second harmonic amplitude. With this innovative method, the measurement of the normal displacement is indirectly used, avoiding the need for very precise contact detection. Moreover, the influence of the tip defect and thermal drift on the measurements are reduced. This method was used for dynamic nanoindentation tests performed on fused silica and on an amorphous polymer (PMMA) because these materials are supposed not to exhibit an indentation size effect at small penetration depths. The amplitude of the second harmonic of the displacement signal was correctly measured at small depths, allowing to calculate the Young’s modulus and the hardness of the tested materials. The mechanical properties calculated with this method are in good agreement with values obtained from classical nanoindentation tests.