Modeling elasto-plastic indentation on layered materials using the equivalent inclusion method

This paper develops a fast semi-analytical model for solving the three-dimensional elasto-plastic contact problems involving layered materials using the Equivalent Inclusion Method (EIM). The analytical elastic solutions of a half-space subjected to a unit surface pressure and a unit subsurface eigenstrain are employed in this model; the topmost layer is simulated by an equivalent inclusion with fictitious eigenstrain. Accumulative plastic deformation is determined by a procedure involving an iterative plasticity loop and an incremental loading process. Algorithms of the fast Fourier transform (FFT) and the Conjugate Gradient Method (CGM) are utilized to improve the computation efficiency. An analytical elastic solution of layered body contact (O’Sullivan and King, 1988) and an indentation experiment result involving a layered substrate (Michler et al., 1999) are used to examine the accuracy of this model. Comparisons between numerical results from the present model and a commercial FEM software (Abaqus) are also presented. Case studies of a rigid ball loaded against a layered elasto-plastic half-space are conducted to explore the effects of the modulus, yield strength, and thickness of the coating on the hardness, stiffness, and plastic deformation of the composite body.     

Receptance behaviour of railway track and subgrade

Summary Vertical dynamic behaviour of a railway track on an elastic halfspace or on a layered halfspace is investigated by a frequency domain analysis. The results are compared with those for a simpler model, where ballast and subgrade are considered as a viscoelastic foundation. In the low- and medium-frequency range up to 250 Hz, great differences are observed between the results of the halfspace model and the results of the viscoelastic foundation model. This is because the damping due to wave propagation and coupling between sleepers cannot be modelled correctly by a viscoelastic foundation. Contradictions observed in the past between measured and calculated results can be explained with the new halfspace model. For frequencies higher than 250 Hz, the influence of the subgrade is negligible, so that here the simpler viscoelastic foundation model can be used. Received 5 January 1998; accepted for publication 7 April 1998     

Geometrically necessary dislocation structure organization in FCC bicrystal subjected to cyclic plasticity

Crystal plasticity finite element analysis of cyclic deformation of compatible type FCC bicrystals are performed. The model specimen used in the analysis is a virtual FCC bicrystal with an isotropic elastic property; therefore, the effect of constraint due to elastic incompatibility does not appear. The results of the analysis show the strain-amplitude-dependence of both the organization of the GND structure and the stress–strain behavior. The calculated stress–strain curve with the largest strain amplitude shows additional cyclic hardening. The microscopic mechanisms of the strain-amplitude-dependent organization of the GND structure and additional cyclic hardening behavior are discussed in terms of the activation of secondary slip system(s). Finally, the effects of the elastic anisotropy, the lattice friction stress and the interaction between dislocations are also argued.     

Structure of weak shock waves in 2-D layered material systems

Shock waves in homogeneous materials in the absence of phase transitions are understood to have a one-wave structure. However, upon loading of a layered heterogeneous material system a two-wave structure is obtained––a leading shock front followed by a complex pattern that varies with time. This dual shock-wave pattern can be attributed to material architecture through which the shock wave propagates, i.e. the impedance (and geometric) mismatch present at various length scales, and nonlinearities arising from material inelasticity and failure.The objective of the present paper is to provide a better understanding of the role of material architecture in determining the structure of weak shock waves in 2-D layered material systems. Normal plate-impact experiments are conducted on 2-D layered material targets to obtain both the precursor decay and the late-time dispersion. The particle velocity at the free surface of the target plate is measured by using a multi-beam VALYN VISAR. In order to understand the effects of layer thickness and the distance of wave propagation on elastic precursor decay and late-time dispersion several different targets with various layer and target thicknesses are employed. Moreover, in order to understand the effects of material inelasticity both elastic–elastic and elastic–viscoelastic bilaminates are utilized.The results of these experiments are interpreted by using asymptotic techniques to analyze propagation of acceleration waves in 2-D layered material systems. The analysis makes use of the Laplace transform and Floquet theory for ODE’s with periodic coefficients [Asymptotic solutions for wave propagation in elastic and viscoelastic bilaminates. In: Developments in Mechanics, Proceedings of the 14th Mid-Eastern Mechanics Conference, vol. 26, no. 8, pp. 399–417]. Both wave-front and late-time solutions for step-pulse loading on layered half-space are compared with the experimental observations. The results of the study indicate that the structure of acceleration waves is strongly influenced by impedance mismatch of the layers constituting the laminates, density of interfaces, distance of wave propagation, and the material inelasticity.     

用薄层法分析层状地基中各种基础的阻抗函数

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法.采用薄层单元和傍轴边界分别模拟层状地基和弹性半空间.利用薄层法位移基本解和容积法推导了层状地基中基础-地基动力相互作用方程及块式基础、桩基础和承台群桩基础阻抗函数的统一计算公式.通过计算半无限弹性地基中桩基础、块式基础和二层地基上基础的阻抗函数验证了方法的适用性.进而计算了某实际层状地基中承台群桩基础的阻抗函数,并与试验结果进行对比,两者吻合较好.本文方法可用于分析弹性层状地基中各种基础的阻抗函数.     

THICKNESS-SHEAR VIBRATION OF A RECTANGULAR QUARTZ PLATE WITH PARTIAL ELECTRODES

We study free vibration of a thickness-shear mode crystal resonator of AT-cut quartz. The resonator is a rectangular plate partially and symmetrically electroded at the center with rectangular electrodes. A single-mode, three-dimensional equation governing the thickness-shear displacement is used. A Fourier series solution is obtained. Numerical results calculated from the series show that there exist trapped thickness-shear modes whose vibration is mainly under the electrodes and decays rapidly outside the electrodes. The effects of the electrode size and thickness on the trapped modes are examined.     

Harmonic waves in one-, two- and three-dimensional composites: Bounds for eigenfrequencies

Harmonic waves in one-, two- and three-dimensional elastic composites with periodic structure are considered. Based on a new quotient recently proposed by Nemat-Nasser, lower and upper bounds for the eigenfrequencies are developed. For illustration waves propagating normal to the layers in layered composites, and normal to the fibers in fiber-reinforced composites, are considered. These examples show that the new quotient is very effective and yields very accurate results for the considered class of problems. While these results are of interest in their own right, they can be used to check the effectiveness of various approximate theories which recently have been proposed by various authors.     

A three-layer structure model for the effect of a soft middle layer on Love waves propagating in layered piezoelectric systems

A three-layer structure model is proposed for investigating the effect of a soft elastic middle layer on the propagation behavior of Love waves in piezoelectric layered systems, with "soft" implying that the bulk-shear-wave velocity of the middle layer is smaller than that of the upper sensitive layer. Dispersion equations are obtained for unelectroded and traction-free upper surfaces which, in the limit, can be reduced to those for classical Love waves. Systematic parametric studies are subsequently carried out to quantify the effects of the soft middle layer upon Love wave propagation, including its thickness, mass density, dielectric constant and elastic coefficient. It is demonstrated that whilst the thickness and elastic coefficient of the middle layer affect significantly Love wave propagation, its mass density and dielectric constant have negligible influence. On condition that both the thickness and elastic coefficient of the middle layer are vanishingly small so that it degenerates into an imperfectly bonded interface, the three-layer model is also employed to investigate the influence of imperfect interfaces on Love waves propagating in piezoelectric layer/elastic substrate systems. Upon comparing with the predictions obtained by employing the traditional shear-lag model, the present three-layer structure model is found to be more accurate as it avoids the unrealistic displacement discontinuity across imperfectly bonded interfaces assumed by the shearlag model, especially for long waves when the piezoelectric layer is relatively thin.     

Active control scheme for improving mass resolution of film bulk acoustic resonators

High mass resolution of sensors based on film bulk acoustic resonators (FBARs) is required for the detection of small molecules with the low concentration.An active control scheme is presented to improve the mass resolution of the FBAR sensors by adding a feedback voltage onto the driving voltage between two electrodes of the FBAR sensors. The feedback voltage is obtained by giving a constant gain and a constant phase shift to the current on the electrodes of the FBAR sensors. The acoustic energy produced by the feedback voltage partly compensates the acoustic energy loss due to the material damping and the acoustic scattering, and thus improves the quality factor and the mass resolution of the FBAR sensors. An explicit expression relating to the impedance and the frequency for an FBAR sensor with the active control is derived based on the continuum theory by neglecting the influence of the electrodes. Numerical simulations show that the impedance of the FBAR sensor strongly depends on the gain and the phase shift of the feedback voltage, and the mass resolution of the FBAR sensor can greatly be improved when the appropriate gain and the phase shift of the feedback voltage are used. The active control scheme also provides an effective solution to improve the resolution of the quartz crystal microbalance (QCM).     

Second-harmonic generation in an infinite layered structure with nonlinear spring-type interfaces

The second-harmonic generation characteristics in the elastic wave propagation across an infinite layered structure consisting of identical linear elastic layers and nonlinear spring-type interlayer interfaces are analyzed theoretically. The interlayer interfaces are assumed to have identical linear interfacial stiffness but can have different quadratic nonlinearity parameters. Using a perturbation approach and the transfer-matrix method, an explicit analytical expression is derived for the second-harmonic amplitude when the layered structure is impinged by a monochromatic fundamental wave. The analysis shows that the second-harmonic generation behavior depends significantly on the fundamental frequency reflecting the band structure of the layered structure. Two special cases are discussed in order to demonstrate such dependence, i.e., the second-harmonic generation by a single nonlinear interface as well as by multiple consecutive nonlinear interfaces. In particular, when the second-harmonic generation occurs at multiple consecutive nonlinear interfaces, the cumulative growth of the second-harmonic amplitude with distance is only expected in certain frequency ranges where the fundamental as well as the double frequencies belong to the pass bands of the layered structure. Furthermore, a remarkable increase of the second-harmonic amplitude is found near the terminating edge of pass bands. Approximate expressions for the low-frequency range are also obtained, which show the cumulative growth of the second-harmonic amplitude with quadratic frequency dependence.     

Analysis of scattering waves in an elastic layered medium by means of the complete eigenfunction expansion form of the Green’s function

A scattering problem due to an object and a plane incident wave in an elastic layered half space is presented in this paper. The complete eigenfunction expansion form of the Green’s function developed by the author and the boundary integral equation method are introduced into the analysis. First, the complete eigenfunction expansion form of the Green’s function is investigated for its application to the scattering problem. A comprehensive explanation is also given for the fact that the complex Rayleigh wave modes exhibit standing waves. Next, a method for the analysis of scattering waves by means of the Green’s function is presented. The advantage of the present method is that the formulation itself is independent of the number of layers and the scattering waves can be decomposed into the modes for the spectra defined for the layered medium. Several numerical calculations are performed to examine the efficiency of the present method as well as the properties of the scattering waves. According to the numerical results, the complete eigenfunction expansion form of the Green’s function provides accurate values for application to a boundary element analysis. The spectral structure and radiation patterns of the scattering wave are presented and investigated. The differences in directionality can be found from the radiation patterns of the scattering waves decomposed into the modes for the spectra.     

Simulation of elastic wave diffraction by multiple strip-like cracks in a layered periodic composite

The problem of numerical simulation of the steady-state harmonic vibrations of a layered phononic crystal (elastic periodic composite) with a set of strip-like cracks parallel to the layer boundaries is solved, and the accompanying wave phenomena are considered. The transfer matrix method (propagator matrix method) is used to describe the incident wave field. It allows one not only to construct the wave fields but also to calculate the pass bands and band gaps and to find the localization factor. The wave field scattered by multiple defects is represented by means of an integral approach as a superposition of the fields scattered by all cracks. An integral representation in the form of a convolution of the Fourier symbols of Green’s matrices for the corresponding layered structures and a Fourier transform of the crack opening displacement vector is constructed for each of the scattered fields. The crack opening displacements are determined by the boundary integral equation method using the Bubnov-Galerkin scheme, where Chebyshev polynomials of the second kind, which take into account the behavior of the solution near the crack edges, are chosen as the projection and basis systems. The system of linear algebraic equations with a diagonal predominance of components arising when the system of integral equations is discretized has a block structure. The characteristics describing qualitatively and quantitatively the wave processes that take place under the diffraction of plane elastic waves by multiple cracks in a phononic crystal are analyzed. The resonant properties of a system of defects and the influence of the relative positions and sizes of defects in a layered phononic crystal on the resonant properties are studied. To obtain clearer results and to explain them, the energy flux vector is calculated and the energy surfaces and streamlines corresponding to them are constructed.     

Thermal ratcheting of solder-bonded elastic and elastoplastic layers

The cyclically growing deflection of solder-bonded elastic and elastoplastic layers subjected to cyclic thermal loading is studied. Finite element analysis of a Si/Sn–95Pb/OFHC-Cu layered structure is performed by taking into account the temperature-dependent viscoplastic behavior of Sn–95Pb as well as the uniaxial ratcheting behavior of OFHC-Cu. A temperature-dependent power law is employed for the viscoplasticity of Sn–95Pb, while a combined nonlinear kinematic and isotropic hardening model is assumed for the cyclic plasticity of OFHC-Cu. It is shown that the temperature-dependent viscoplasticity of Sn–95Pb and the uniaxial ratcheting of OFHC-Cu are the controlling factors for the cyclic growth of deflection of the layered structure under temperature cycling. It is also shown that cyclic hardening of OFHC-Cu plays an important role for the cyclic growth of deflection, and that elastic stress in the Si layer cyclically develops noticeably if the cyclic growth of deflection is significant.     

Theoretical modeling of guided wave propagation in a sandwich plate subjected to transient surface excitations

A fast and efficient two-dimensional (2D) semi-analytical model based on global matrix method is developed to study the general characteristics of guided wave propagation in a honeycomb composite sandwich structure (HCSS) subjected to time-dependent transient surface excitations. The HCSS used in this study has an extremely lightweight and thick nomex honeycomb core, which is sandwiched between two thin graphite woven composite skins. The homogenized material properties of the skin and the core are considered to be elastic and quasi-isotropic in nature. Far-field time history of surface displacements are calculated for vertical and horizontal tone-burst surface excitations that are representative of thickness and radial mode of vibrations of piezoelectric transducers, respectively. Results are compared with those obtained from Finite element modeling (FEM) in LS-DYNA showing good agreement. Wavelet transform is then performed on the time-domain signal to obtain the group velocities of the propagating modes for their accurate identification on the basis of the theoretical dispersion curve. It is found that the response signal is dominated by the first anti-symmetric mode for a vertical excitation, whereas, the signal characteristics are multimodal in nature with dominating higher order symmetric and anti-symmetric modes for a horizontal excitation. The model is expected to be helpful for appropriate guided wave mode tuning and rapid analysis of data for experimental detection of disbonds in these novel structures.     

硅基薄膜电极充放电过程应力演化解析解研究

为了研究硅电极在充放电过程中的应力演化,在忽略弹性变形的情况下提出了一个考虑粘塑性的简化力学模型,分别导出了出现相分离和不出现相分离时薄膜电极内应力场的解析解,该模型的计算结果与现有的实验结果相吻合.计算结果表明,当存在相分离时,电极中的应力取决于薄膜电极的厚度和充电速率的乘积;在未出现相分离现象时,应力仅取决于充电速...     

Surface electro-elastic shear horizontal waves in a layered structure with a piezoelectric substrate and a hard dielectric layer

The existence and behaviour of surface electro-elastic shear horizontal waves in a layered structure consisting of a piezoelectric substrate of crystal class 6, 4, 6mm, or 4mm mechanically bonded at its upper surface to an elastic dielectric layer and bounded by an adjacent dielectric medium is considered when the shear bulk wave velocity in the elastic layer is greater than or equal to that in the substrate. The dispersion equation for the existence of the surface electro-elastic SH waves with respect to the phase velocity is obtained which includes all the above crystal classes i.e. the surface wave problems related to all these classes are presented in a single mathematical model. The investigation of the solutions of the dispersion equation is carried out and all the possible cases of the behaviour of the surface electro-elastic SH wave depending on the electro-mechanical coefficients of the layered structure are revealed.     

Properties of Love waves in layered piezoelectric structures

Love waves propagating in a layered structure with an elastic layer deposited on a piezoelectric substrate are analytically investigated. We present a general dispersion equation that describes the properties of Love waves in the structure. A detailed discussion regarding the dispersion equation is presented, and the parameters for Love-mode sensors are also introduced. The properties of Love waves are illustrated by means of sample results for a layered structure with an SiO₂ layer sputtered on an ST-cut 90°X-propagating quartz substrate. Interestingly, we found that a threshold-normalized layer thickness existed for the fundamental Love mode in such a structure.     

Overall dynamic constitutive relations of layered elastic composites

A method for the homogenization of a layered elastic composite is presented. It allows direct, consistent, and accurate evaluation of the averaged overall frequency-dependent dynamic material constitutive relations without the need for a point-wise solution of the field equations. When the spatial variation of the field variables is restricted by Bloch-form (Floquet-form) periodicity, then these relations together with the overall conservation and kinematical equations accurately yield the displacement or stress mode-shapes and, necessarily, the dispersion relations. The method can also give the point-wise solution of the elastodynamic field equations (to any desired degree of accuracy), which, however, is not required for the calculation of the average overall properties. The resulting overall dynamic constitutive relations are general and need not be restricted by the Bloch-form periodicity.The formulation is based on micromechanical modeling of a representative unit cell of the composite. For waves in periodic layered composites, the overall effective mass-density and compliance (stiffness) are always real-valued whether or not the corresponding unit cell (representative volume element used as a unit cell) is geometrically and/or materially symmetric. The average strain and linear momentum are coupled and the coupling constitutive parameters are always each others' complex conjugates. We separate the overall constitutive relations, which depend only on the composition and structure of the unit cell, from the overall field equations which hold for any elastic composite; i.e., we use only the local field equations and material properties to deduce the overall constitutive relations. Finally, we present solved numerical examples to further clarify the structure of the averaged constitutive relations and to bring out the correspondence of the current method with recently published results.     

Cohesive failure analysis of an array of IC chips bonded to a stretched substrate

The paper presents a mechanical model for predicting the cohesive failure of a periodic array of integrated circuit (IC) chips adhesively bonded to a stretched substrate. A unit cell of the layered structure consisting of the IC chips, adhesive layer, and substrate is modeled as an assembly of two elastic Timoshenko beams, representing the chip and substrate, connected by an elastic interface, representing the adhesive. Accordingly, the stresses and energy release rate (ERR) in the adhesive layer – responsible for the premature cracking of the adhesive and debonding of the IC chips – are identified with the corresponding quantities computed for the elastic interface. Expressions for the adhesive stresses and ERR are given in terms of geometrical dimensions and material properties, combined with integration constants obtained numerically via the multi-segment analysis method. For comparison, the stresses in the adhesive are also computed based on a finite element model, and the ERR is evaluated using the virtual crack-closure technique (VCCT). The analytical predictions and numerical results match fairly well, considering the effects of key factors, such as the distance between adjacent chips, the chip size, the material properties of adhesive and substrate. The interaction between the chips is shown to have relevant effects on the adhesive stresses. In particular, only the mode II contributes to the ERR which increases with the ratio of the chip size to the distance between the chips and with the compliance of the adhesive and substrate layers.     

An extended finite element method for modeling near-interfacial crack propagation in a layered structure

The extended finite element method (XFEM) is applied for the simulation of near-interfacial crack propagation in a metal–ceramic layered structure. An experimental evidence indicates that, in a ceramic–metal–ceramic sandwich structure, a near-interfacial crack in the ceramic layer can be drawn to or deflect away from the metal layer depending on the difference in elastic properties across the interface. To model near-interfacial fracture, only the Heaviside functions are used for the XFEM, and the vector level set method is employed for efficient evaluation of the enrichment functions. The crack propagation paths predicted by the XFEM simulation are found to be consistent with the experimental observation.