Indentation responses of piezoelectric films ideally bonded to an elastic substrate

The influences of elastic substrate on the indentation force, contact radius, electric potential and electric charge responses of piezoelectric film/substrate systems are investigated by the integral transform method. The film is assumed to be ideally bonded to the substrate and the contact interaction between the indenter and the film is assumed to be frictionless, with three kinds of axisymmetric insulating and conducting indenters (i.e., punch, cone and sphere) considered. Obtained results show that when the ratio of the contact radius to the film thickness is close to zero, the influences of the elastic substrate disappear and the indentation behaviors converge to the piezoelectric half space solutions while the indentation responses approach the corresponding ones of elastic half space as the ratio gets to infinity. The transition between the piezoelectric and the elastic half space indentation solutions for the film/substrate system is quantified in terms of the film thickness and the elasticity of the substrate. Finite element analysis on an insulating sphere indentation is conducted to verify the numerical calculations and good agreement is observed. The obtained results are believed to be useful for developing experimental techniques to extract the material properties of piezoelectric film/substrate systems.     

Mechanics of indentation for piezoelectric thin films on elastic substrate

Frictionless normal indentation problem of rigid flat-ended cylindrical, conical and spherical indenters on piezoelectric film, which is either in frictionless contact with or perfectly bonded to an elastic half-space (substrate), is investigated. Both conducting and insulating indenters are considered. With Hankel transform, the general solutions of the homogeneous governing equations for the piezoelectric layer and the elastic half-space are presented. Using the boundary conditions for a vertical point force or a point electric charge, and the boundary conditions on the film/substrate interface, the Green’s functions can be obtained by solving sets of simultaneous linear algebraic equations. The solution of the indentation problem is obtained by integrating these Green’s functions over the contact area with unknown surface tractions or electric charge distribution, which will be determined from the boundary conditions on the contact surface between the indenter and the film. The solution is expressed in terms of dual integral equations that are converted to a Fredholm integral equation of the second kind and solved numerically. Numerical examples are also presented. The comparison between two film/substrate bonding conditions is made. It shows that the indentation rigidity of the film/substrate system is lower when the film is in frictionless contact with the substrate. The effects of the Young’s modulus and Poisson’s ratio of the elastic substrate, indenter electrical condition and indenter prescribed electric potential on the indentation responses are presented.     

Contact with stick zone between an indenter and a thin incompressible layer

The dominant asymptotic term for the indentation of a thin elastic incompressible layer by an axisymmetric rigid indenter is considered. Complete adhesion is supposed everywhere in the contact area or else in a given inner region surrounded by an annular frictionless zone. Both the problems are formulated in the form of systems of coupled dual integral equations. Using operators transforming kernels of the Hankel transform into kernels of the Weber–Orr transform, the dual integral equations are reduced to systems of Fredholm integral equations of the second kind whose structures permit deriving asymptotic solutions. Simple expressions for the contact stresses, the penetration depth, and the contact radius in the case of an unknown contact area are obtained. Explicit formulae, derived for the flat and power law indenter profiles, allow us to analyze how stick and frictionless zones affect mechanical characteristics. Results manifest that the punch penetration exhibits strong sensitivity to contact conditions inspite of the fact that the radial traction is small. A conical indenter is less sensitive than flat-ended and spherical indenters.     

Orientation effect on the Boussinesq indentation of a transversely isotropic piezoelectric material

Three-dimensional finite element analysis was used to study the effect of the angle between the loading direction and the axisymmetric direction on the indentation behavior of a transversely isotropic piezoelectric half-space by a cylindrical indenter of flat end. Two cases were considered in the analysis, which included (a) the indentation by an insulating indenter, and (b) the indentation by a conducting indenter. Both the indentation load and the indentation-induced potential were found to be proportional to the indentation depth. Using the simulation results and the analytical relationship for the indentation by a rigid, insulating indenter, semi-analytical relationships were developed between the indentation load and the indentation depth and between the indentation-induced potential on the indenter and the indentation depth for the conducting indenter, respectively. The proportionality between the indentation-induced potential and the indentation depth is only a function of the angle between the loading direction and the poling direction, independent of the type of indenters, which may be used to measure the relative direction of the loading axis to the axisymmetric axis of transversely piezoelectric materials from the indentation test.     

基于能量密度等效的超弹性压入模型与双压试验方法

针对超弹性材料压入问题, 本文基于能量密度中值等效原理, 提出了描述球、平面、锥3类压头独立压入下载荷、深度、压头几何尺寸和Mooney-Rivlin本构关系参数之间关系的半解析超弹性压入模型(semi-theoretical hyperelastic-material indentation model, SHIM), 进而提出了球、平面、锥压入组合的双压试验方法(indentation method due to dual indenters, IMDI). 正向验证表明, 基于系列超弹性材料的本构关系参数, 由SHIM分别预测的球、平面、锥3类压入下的载荷-位移曲线与有限元分析(finite element analysis, FEA)结果之间密切吻合; 反向验证表明, 基于系列超弹性材料的FEA条件本构关系下3类压入的载荷-位移曲线, 由双压试验方法预测的Mooney-Rivlin本构关系与FEA条件本构关系密切吻合. 针对3种超弹性橡胶, 完成了球、平面、锥压入试验, 应用双压试验方法获得的3组Mooney-Rivlin本构关系均与单轴拉伸试验结果吻合良好.      

Spherical indentation of freestanding circular thin films in the membrane regime

We present theoretical and experimental results to describe the mechanics of indentation of a clamped circular membrane with a frictionless spherical indenter. Analytical expressions and numerical simulations are presented for the relationships between contact radius, finite indentation strains (and stresses), pre-stretch, loads and deflection. These closed-form solutions are contrasted with point-load models that neglect the contact size (i.e. classical Schwerin-type solutions), and lead to important differences in the indentation strain and load-deflection response. The accuracy of these closed form expressions is illustrated by comparisons with detailed numerical results and experiments on thin elastomer films. We show that the closed-form solutions can be used to extract mechanical properties from indentation testing of freestanding films, with important implications for developing new tests on nanoscale films and/or compliant materials such as polymers and biological substances.     

Stiffness relations for piezoelectric indentation of flat and non-flat punches of arbitrary planform: Applications to probing nanoelectromechanical properties of materials

Stiffness relations for voltage-dependent contact mechanics of piezoelectric material are derived for an indenter of arbitrary planform under normal force, centrally or non-centrally applied, and electric charge distribution at the base. Relations between indentation depth, indentation force, electric potential and electric charge are explicitly given in terms of indenter's geometry and piezoelectric material constants. The analysis covers indenters with non-flat base approximated by a second-order surface; elliptic paraboloid is considered as an example. In the case of the elliptic non-flat planform, the derived stiffness relations are exact; otherwise, they are approximate and are shown to have good accuracy. The stiffness relations are given in elementary functions and are obtained by utilizing the recently established principle of correspondence between the piezoelectric and purely elastic problems. Besides contributing to extension of Hertzian mechanics to piezoelectric materials, these results are essential for quantitative interpretation of the scanning probe microscopy and piezoelectric nanoindentation data on piezoelectric, ferroelectric, and multiferroic materials.     

Small scale,grain size and substrate effects in nano-indentation experiment of film–substrate systems

The mechanical properties of film–substrate systems have been investigated through nano-indentation experiments in our former paper (Chen, S.H., Liu, L., Wang, T.C., 2005. Investigation of the mechanical properties of thin films by nano-indentation, considering the effects of thickness and different coating–substrate combinations. Surf. Coat. Technol., 191, 25–32), in which Al–Glass with three different film thicknesses are adopted and it is found that the relation between the hardness H and normalized indentation depth h/t, where t denotes the film thickness, exhibits three different regimes: (i) the hardness decreases obviously with increasing indentation depth; (ii) then, the hardness keeps an almost constant value in the range of 0.1–0.7 of the normalized indentation depth h/t; (iii) after that, the hardness increases with increasing indentation depth. In this paper, the indentation image is further investigated and finite element method is used to analyze the nano-indentation phenomena with both classical plasticity and strain gradient plasticity theories. Not only the case with an ideal sharp indenter tip but also that with a round one is considered in both theories. Finally, we find that the classical plasticity theory can not predict the experimental results, even considering the indenter tip curvature. However, the strain gradient plasticity theory can describe the experimental data very well not only at a shallow indentation depth but also at a deep depth. Strain gradient and substrate effects are proved to coexist in film–substrate nano-indentation experiments.     

A Circular Indenter on a Piezoelectric Layer

This paper investigates a piezoelectric layer with a rigid indenter on its surface. Exact solution is given for a piezoelectric medium whose thickness is considerably larger than the diameter of the indenter. Different electrical boundary conditions that employ conducting or insulating indenters are presented. Effect of the permittivity of air (which surrounds the piezoelectric medium) is considered and is found to be negligible. Expressions for the singular mechanical and electric fields near the indenter front are established. Those expressions are useful for investigating the possible failure behavior of piezoelectric material near the indenter front. In addition, a numerical solution technique for an indentured piezoelectric layer of finite thickness is also given.     

The interaction of two parallel Mode-I limited-permeable cracks in a functionally graded piezoelectric material

In this paper, the interaction of two parallel Mode-I limited-permeable cracks in a functionally graded piezoelectric material was investigated by using the generalized Almansi's theorem. In the analysis, the electric permittivity of the air inside the crack was considered. The problem was formulated through Fourier transform into two pairs of dual integral equations, in which unknown variables are jumps of displacements across the crack surface. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. The solution of the present paper shows that the singular stresses and the singular electric displacements at the crack tips in functionally graded piezoelectric materials carry the same forms as those in homogeneous piezoelectric materials; however, the magnitudes of intensity factors depend on the electric permittivity of the air inside the crack and the gradient parameter of functionally graded piezoelectric material properties. It was also revealed that the crack shielding effect is also present in functionally graded piezoelectric materials.     

压电涂层-功能梯度压电界面层-基底结构的二维接触问题研究

在多层压电元件中,由于界面处材料成分和性质的突变,常常导致界面处应力集中,使得界面处出现开裂或蠕变现象,从而大大缩短了压电元件的使用寿命。功能梯度压电材料作为界面层,可有效的缓解界面材料不匹配导致的破坏。本文主要研究利用功能梯度压电材料界面层连接压电涂层和基底,分析三层结构在圆柱型压头作用下的力电响应。利用傅里叶积分变换技术,本文将压电涂层-功能梯度压电层-基底结构在刚性圆柱压头作用下的二维平面应变接触问题转化为带有柯西核的奇异积分方程。运用高斯-切比雪夫积分公式,将奇异积分方程转化为线性方程组并对其进行数值求解,得到压电涂层-功能梯度压电层-基底结构在圆柱形压头作用下的应力分布和电位移分布。数值结果表明,梯度压电材料参数的变化对结构中的力电响应具有重要的影响。本文研究结果对于利用功能梯度压电界面层消除界面处的应力不连续导致的界面破坏具有重要的理论指导意义,研究结果可为功能梯度压电材料界面层的设计提供帮助。     

A numerical approach to indentation technique to evaluate material properties of film-on-substrate systems

Spherical indentation approach (Lee et al., 2005, Lee et al., 2010) for the evaluation of bulk material properties is extended to that for elastic–plastic properties of film-on-substrate systems. Our interest focuses on single isotropic, metallic, and elastic–plastic film on a substrate, and we do not consider the size effects in plasticity behavior. We first determine the optimal data acquisition location, where the strain gradient is the least and the effect of friction is negligible. Dimensional analysis affords the mapping parameters as functions of normalized indentation variables. An efficient way is further introduced to reduce both the number of analyses and the regression order of mapping functions. The new numerical approach to the film indentation technique is then proposed by examining the finite element solutions at the optimal point. With the new approach, the values of elastic modulus, yield strength, and strain-hardening exponent of film materials are successfully obtained from the spherical indentation tests. We have shown that the effective property ranges such as indenter properties, substrate modulus, and E/E_s ratio can be extended without additional simulations and even loss of accuracy. For other ranges of variables or other properties, which are not dealt with in this study, this methodology is applicable through resetting FEA variables and finding proper normalized parameters.     

金属材料的强度与应力-应变关系的球压入测试方法

压入法获取材料单轴应力-应变关系和抗拉强度对服役结构完整性评价有重要的基础意义.假定材料均匀连续、各向同性、应力应变关系符合Hollomon律,基于能量等效假定,即代表性体积单元(representativevolume element, RVE)的vonMises等效和有效变形域内能量中值等效假定,本文提出了关联材料载荷、深度、球压头直径和Hollomon律的四参数半解析球压入(semi-analyticalspherical indentation,SSI)模型.通过球压入载荷-深度试验关系获得材料的应力-应变关系和抗拉强度.考虑压入过程中的损伤效应,针对金属材料提出了用于球压入测试的材料弹性模量修正模型.对11种延性金属材料完成了球压入试验,采用本文提出的球压入试验方法测到的弹性模量、应力-应变关系和抗拉强度与单轴拉伸试验结果吻合良好.      

Basic solution of a mode-I limited-permeable crack in functionally graded piezoelectric materials

In this paper, the basic solution of a mode-I crack in functionally graded piezoelectric materials was investigated by using the generalized Almansi’s theorem. In the analysis, the electric permittivity of air inside the crack were considered. To make the analysis tractable, it was assumed that the shear modulus, piezoelectric constants and dielectric constants vary exponentially with coordinate parallel to the crack. The problem was formulated through Fourier transform into two pairs of dual integral equations, in which the unknown variables are jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. The solution of the present paper shows that the effects of the electric boundary conditions on the electric displacement fields near the crack tips can not be ignored. Simultaneously, the solution of the present paper will revert to a closed form one when the functionally graded parameter equals to zero.     

Basic solution of two parallel mode-I permeable cracks in functionally graded piezoelectric materials

The basic solution of two parallel mode-I permeable cracks in functionally graded piezoelectric materials was studied in this paper using the generalized Almansi’s theorem. To make the analysis tractable, it was assumed that the shear modulus varies exponentially along the horizontal axis parallel to the crack. The problem was formulated through a Fourier transform into two pairs of dual integral equations, in which unknown variables are jumps of displacements across the crack surface. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. The solution of the present paper shows that the singular stresses and the singular electric displacements at the crack tips in functionally graded piezoelectric materials carry the same forms as those in homogeneous piezoelectric materials; however, the magnitudes of intensity factors depend on the gradient of functionally graded piezoelectric material properties. It was also revealed that the crack shielding effect is also present in functionally graded piezoelectric materials.     

Strip electric-magnetic breakdown model in a magnetoelectroelastic medium

Extending the polarization saturation model [Gao et al., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491-510] and the dielectric breakdown (DB) model [Zhang et al., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311-327] in piezoelectric materials, the Strip Electric-Magnetic Breakdown (SEMB) model is proposed for electrically and magnetically impermeable crack in a magnetoelectroelastic medium to study the effect of the nonlinear character of electric field and magnetic field on fracture of magnetoelectroelastic materials. In the SEMB model, the electric field in the strip of the electric breakdown zone ahead of the crack tip is equal to the electric breakdown strength, while the magnetic filed in the strip of the magnetic breakdown zone is equal to the magnetic breakdown strength. By using the extended Stroh formalism and the extended dislocation modeling of a crack, the Griffith crack problem under the electrically and magnetically elastic-plastic condition in a magnetoelectroelastic medium is reduced to a set of dual integral equations. The sizes of the electric breakdown zone and the magnetic breakdown zone, the extended intensity factors and the local J-integral are obtained. The effect of the combined mechanical-electric-magnetic loadings on the local J-integral is studied.     

The equivalent axisymmetric model for Berkovich indenters in power-law hardening materials

Nix and Gio [Nix, W.D., Gao, H.J., 1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. Journal of the Mechanics and Physics of Solids 46, 411–425] established an important relation between the micro-indentation hardness and indentation depth for axisymmetric indenters. For the Berkovich indenter, however, this relation requires an equivalent cone angle. Qin et al. [Qin, J., Huang, Y., Xiao, J., Hwang, K.C., 2009. The equivalence of axisymmetric indentation model for three-dimensional indentation hardness. Journal of Materials Research 24, 776–783] showed that the widely used equivalent cone angle from the criterion of equal base area leads to significant errors in micro-indentation, and proposed a new equivalence of equal cone angle for iridium. It is shown in this paper that this new equivalence holds for a wide range of plastic work hardening materials. In addition, the prior equal-base-area criterion does not hold because the Berkovich indenter gives much higher density of geometrically necessary dislocations than axisymmetric indenter. The equivalence of equal cone angle, however, does not hold for Vickers indenter.     

An analysis of nanoindentation in linearly elastic solids

The conventional method to extract elastic properties in the nanoindentation of linearly elastic solids relies primarily on Sneddon’s solution (1948). The underlying assumptions behind Sneddon’s derivation, namely, (1) an infinitely large incompressible specimen; (2) an infinitely sharp indenter tip, are generally violated in nanoindentation. As such, correction factors are commonly introduced to achieve accurate measurements. However, little is known regarding the relationship between the correction factors and how they affect the overall accuracy. This paper first proposes a criterion for the specimen’s geometry to comply with the first assumption, and modifies Sneddon’s elastic relation to account for the finite tip radius effect. The relationship between the finite tip radius and compressibility of the specimen is then examined and a composite correction factor that involves both factors, derived. The correction factor is found to be a function of indentation depth and a critical depth is derived beyond which, the arbitrary finite tip radius effect is insignificant. Techniques to identify the radius of curvature of the indenter and to decouple the elastic constants (E and ν) for linear elastic materials are proposed. Finally, experimental results on nanoindentation of natural latex are reported and discussed in light of the proposed modified relation and techniques.     

Basic solution of two parallel non-symmetric permeable cracks in piezoelectric materials

The behavior of two parallel non-symmetric cracks in piezoelectric materials subjected to the anti-plane shear loading was studied by the Schmidt method for the permeable crack electric boundary conditions. Through the Fourier transform, the present problem can be solved with two pairs of dual integral equations ip which the unknown variables are the jumps of displacements across crack surfaces. To solve the dual integral equations, the jumps of displacements across crack surfaces were directly expanded in a series of Jacobi polynomials. Finally, the relations between electric displacement intensity factors and stress intensity factors at crack tips can be obtained. Numerical examples are provided to show the effect of the distance between two cracks upon stress and electric displacement intensity factors at crack tips. Contrary to the impermeable crack surface condition solution, it is found that electric displacement intensity factors for the permeable crack surface conditions are much smaller than those for the impermeable crack surface conditions. At the same time, it can be found that the crack shielding effect is also present in the piezoelectric materials.     

SCATTERING OF THE HARMONIC STRESS WAVE BY CRACKS IN FUNCTIONALLY GRADED PIEZOELECTRIC MATERIALS

The present paper considers the scattering of the time harmonic stress wave by a single crack and two collinear cracks in functionally graded piezoelectric material （FGPM）. It is assumed that the properties of the FGPM vary continuously as an exponential function. By using the Fourier transform and defining the jumps of displacements and electric potential components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement and electric potential components across the crack surface are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the influences of material properties on the dynamic stress and the electric displacement intensity factors.