Singular integral equation method for contact problem for rigidly connected punches on elastic half-plane

In the contact problem of a rigid flat-ended punch on an elastic half-plane, the contact stress under punch is studied. The angle distribution for the stress components in the elastic medium under punch is achieved in an explicit form. From obtained singular stress distribution, the punch singular stress factor (abbreviated as PSSF) is defined. A fundamental solution for the multiple flat punch problems on the elastic half-plane is investigated where the punches are disconnected and the forces applied on the punches are arbitrary. The singular integral equation method is suggested to obtain the fundamental solution. Further, the contact problem for rigidly connected punches on an elastic half-plane is considered. The solution for this problem can be considered as a superposition of many particular fundamental solutions. The resultant forces on punches are the undetermined unknowns in the problem, which can be evaluated by the condition of relative descent between punches. Finally, the resultant forces on punches can be determined, and the PSSFs at the corner points can be evaluated. Numerical examples are given.     

Dynamic behavior of a moving frictional punch over the surface of anisotropic materials

This paper investigates dynamic, frictional contact of a moving punch over the surface of anisotropic materials. An eigenvalue analysis of the governing equations is performed. The application of the complex function theory produces a singular integral equation exhibiting a non-square-root or unconventional singularity. Numerical tests demonstrate that both the friction coefficient and the moving velocity contribute to the contact behaviors under a moving punch with a flat or cylindrical profile. Furthermore, the present results illustrate that the surface in-plane stress possesses singularity and discontinuation at both edges of the flat punch and has a tensile spike at one edge of the cylindrical punch, which may account for the fatigue and fracture under the contact loading.     

Static frictional indentation of an elastic half-plane by a rigid unsymmetrical punch

Static rigid 2-D indentation of a linearly elastic half-plane in the presence of Coulomb friction which reverses its sign along the contact length is studied. The solution approach lies within the context of the mathematical theory of elastic contact mechanics. A rigid punch, having an unsymmetrical profile with respect to its apex and no concave regions, both slides over and indents slowly the surface of the deformable body. Both a normal and a tangential force may, therefore, be exerted on the punch. In such a situation, depending upon the punch profile and the relative magnitudes of the two external forces, a point in the contact zone may exist at which the surface friction changes direction. Moreover, this point of sign reversal may not coincide, in general, with the indentor's apex. This position and the positions of the contact zone edges can be determined only by first constructing a solution form containing the three problem's unspecified lengths, and then solving numerically a system of non-linear equations containing integrals not available in closed form.The mathematical procedure used to construct the solution deals with the Navier-Cauchy partial differential equations (plane-strain elastostatic field equations) supplied with boundary conditions of a mixed type. We succeed in formulating a second-kind Cauchy singular integral equation and solving it exactly by analytic-function theory methods.Representative numerical results are presented for two indentor profiles of practical interest—the parabola and the wedge.     

Analysis and numerical solution of a piezoelectric frictional contact problem

We consider a mathematical model which describes the frictional contact between an electro-elastic–visco-plastic body and a conductive foundation. The contact is modelled with normal compliance and a version of Coulomb’s law of dry friction, in which the stiffness and the friction coefficients depend on the electric potential. We derive a variational formulation of the problem and we prove an existence and uniqueness result. The proof is based on a recent existence and uniqueness result on history-dependent quasivariational inequalities obtained in [15]. Then we introduce a fully discrete scheme for solving the problem and, under certain solution regularity assumptions, we derive an optimal order error estimate. Finally, we present some numerical results in the study of a two-dimensional test problem which describes the process of contact in a microelectromechanical switch.     

Non-local theory solution of two collinear mode-I cracks in piezoelectric materials

In this paper, the basic solution of two collinear cracks in a piezoelectric material plane subjected to a uniform tension loading is investigated by means of the non-local theory. Through the Fourier transform, the problem is solved with the help of two pairs of integral equations, in which the unknown variables are the jumps of displacements across the crack surfaces. To solve the integral equations, the jumps of displacements across the crack surfaces are directly expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the interaction of two cracks, the materials constants and the lattice parameter on the stress field and the electric displacement field near crack tips. Unlike the classical elasticity solution, it is found that no stress and electric displacement singularities are present at crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to using the maximum stress as a fracture criterion in piezoelectric materials.     

Fracture of a piezoelectric fiber/elastic matrix composite

A piezoelectric fiber/elastic matrix system subjected to axially symmetric mechanical and electric loads is considered. The fiber contains a penny-shaped crack located at its center perpendicularly to the fiber. By using the Fourier and Hankel transforms, the problem is reduced to the solution of an integral equation. Numerical solutions for the crack tip fields are obtained for various crack sizes and different fiber volume fractions. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 42, No. 3, pp. 301–318, May–June, 2006.     

Uniform stabilization of a thermo piezoelectric/piezomagnetic model

We consider an evolution system describing an electro/magneto/thermoelastic phenomenon in a bounded domain of R³${\mathbb{R}}^3$. The resulting anisotropic model is a coupling between two hyperbolic equations, two elliptic equations and one parabolic equation. Assuming that a dissipative mechanism is present on the boundary and suitable boundary conditions are given we prove that the total energy decays exponentially as t?+∞$t?+\infty$ provided the region satisfies a geometric condition.     

Homogenization of a viscoelastic matrix in linear frictional contact

The paper is devoted to study of acoustic wave propagation in a partially consolidated composite material containing loose particles. Friction of particles against the consolidated part of the material causes mechanical energy dissipation. This situation is modelled by assuming that the medium has a periodic microstructure changing rapidly on the small scale ε. Each of the periodic microscopic cells is composed of a viscoelastic matrix containing a rigid particle in frictional contact with the matrix. We use the methods of two‐scale convergence to obtain effective acoustic equations for the homogenized material. The effective equations are history‐dependent and contain the body force term, reminiscent of the well‐known Stokes drag force. Copyright © 2004 John Wiley & Sons, Ltd.     

Local solutions to a model of piezoelectric materials

A local existence theorem is proved for a non‐linear coupled system modelling the electromechanical motion of a one‐dimensional piezoelectric body with domain switching. The system is composed by a heat equation describing the behaviour of the number of electric dipoles and by a wave equation governing the dynamic of the electric displacement. The main coupling in the system appears in the time‐dependent velocity of the waves depending on the number of electric dipoles. The proof of the result relies on a time decay estimate satisfied by the number of electric dipoles and an uniform estimate of the solution of the regularized wave equation. Copyright © 2004 John Wiley & Sons, Ltd.     

Galerkin solution of a singular integral equation with constant coefficients

Galerkin methods are used to approximate the singular integral equation

with solution φ having weak singularity at the endpoint −1, where a, b≠0 are constants. In this case φ is decomposed as φ(x)=(1−x)^α(1+x)^βu(x), where β=−α, 0<α<1. Jacobi polynomials are used in the discussions. Under the conditions fH^μ[−1,1] and k(t,x)H^μ,μ[−1,1]×[−1,1], 0<μ<1, the error estimate under a weighted L² norm is O(n^−μ). Under the strengthened conditions f^″H^μ[−1,1] and , 2α−<μ<1, the error estimate under maximum norm is proved to be O(n^2α−−μ+), where , >0 is a small enough constant.     

Elementary solution of contact problems for a transversely isotropic elastic layer bonded to a rigid foundation

The contact problem for an arbitrary punch acting on a transversely isotropic elastic layer bonded to a rigid foundation is solved by the generalized images method developed by the author earlier. The problem is reduced to that of an electrostatic problem of infinite row of coaxial charged disks in the shape of the domain of contact. The solution can be obtained by the method of iteration, collocations or any other standard procedure for solving integral equations. Exact inversion can be obtained in the case of a circular domain of contact. The mean value theorem can be used for estimation of the resultant force and tilting moment acting on a punch of arbitrary shape and circular domain of contact. A limiting case of general solution gives the solution for an isotropic layer. (Received: August 11, 2003)     

Pure parametric excitation of a micro cantilever beam actuated by piezoelectric layers

This paper deals with the stability analysis of transverse motions of a cantilever microbeam sandwiched by two piezoelectric layers located on the lower and upper surfaces of the microbeam. Application of same DC and AC voltages to the upper and lower piezoelectric layers creates an axial force with steady and time-varying components. The eigenfunction expansion of the transverse motion equation leads to the creation of a Mathieu type parametric equation which is mostly seen in the stability analysis of the structures in the literature; using Floquet theory for single degree of freedom systems the stable and unstable regions of the problem are investigated. The effect of viscous damping and DC voltage on the stability region of the problem is also studied. The results show the stabilizing effect of the viscous damping and positive DC voltage on the behavior of the microbeam. The achieved results are finally compared with those reported in the literature.     

直线上带平移的奇异积分方程的解

该文讨论了直线上带实平移或复平移的奇异积分方程的解，获得了方程的可解条件，证明了方程解的唯一性，给出了方程有解的情况下解的积分或级数表达式。     

一类非线性奇异积分方程的离散近似解及其误差估计

当函数f(x,t,u)满足一些[1]中常假定的条件时,我们可借助算子S~(N)和不等式证明非线性奇异积分方程有唯一的离散近似解,这个解可用关于距离的逐次逼近法得到.     

A solution method for the sub-surface stresses and local deflection of a semi-infinite inhomogeneous elastic medium

This paper proposes analytical Fourier series solutions (based on the Airy stress function) for the local deflection and subsurface stress field of a two-dimensional graded elastic solid loaded by a pre-determined pressure distribution. We present a selection of numerical results for a simple sinusoidal pressure which indicates how the inhomogeneity of the solid affects its behaviour. The model is then adapted and used to derive an iterative algorithm which may be used to solve for the contact half width and pressure induced from contact with a rigid punch. Finally, the contact of a rigid cylindrical stamp is studied and our results compared to those predicted by Hertzian theory. It is found that solids with a slowly varying elastic modulus produce results in good agreement with those of Hertz whilst more quickly varying elastic moduli which correspond to solids that become stiffer below the surface give rise to larger maximum pressures and stresses whilst the contact pressure is found to act over a smaller area.     

Effects of electric and mechanical loads on the morphological evolution of a void in piezoelectric films

This paper reports the result of an investigation into the effect of electric and mechanical loads on the morphological evolution of a void in piezoelectric materials based on a model for the morphological evolution of a void, the thermodynamics potential and energy principle. Thus, the path and the bifurcation condition of the morphological evolution of the void in piezoelectric materials are described, which gives some insight into the reliability of piezoelectric films under electric and mechanical loads.     

Exact solutions of functionally graded piezoelectric material sandwich cylinders by a modified Pagano method

The three-dimensional (3D) coupled analysis of simply-supported, functionally graded piezoelectric material (FGPM) circular hollow sandwich cylinders under electro-mechanical loads is presented. The material properties of each FGPM layer are regarded as heterogeneous through the thickness coordinate, and obey an exponent-law dependent on this. The Pagano method is modified to be feasible for the study of FGPM sandwich cylinders. The modifications are as follows: a displacement-based formulation is replaced by a mixed formulation; a set of the complex-valued solutions of the system equations is transferred to the corresponding set of real-valued solutions; a successive approximation method is adopted to approximately transform each FGPM layer into a multilayered piezoelectric one with an equal and small thickness for each layer in comparison with the mid-surface radius, and with the homogeneous material properties determined in an average thickness sense; and a transfer matrix method is developed, so that the general solutions of the system equations can be obtained layer-by-layer, which is significantly less time-consuming than the usual approach. A parametric study is undertaken of the influence of the aspect ratio, open- and closed-circuit surface conditions, and material-property gradient index on the assorted field variables induced in the FGPM sandwich cylinders.     

Solvability of a dynamic contact problem between a piezoelectric body and a conductive foundation

We consider a mathematical model which describes the dynamic process of contact between a piezoelectric body and an electrically conductive foundation. We model the material’s behavior with a nonlinear electro-viscoelastic constitutive law; the contact is frictionless and is described with the normal compliance condition and a regularized electrical conductivity condition. We derive a variational formulation for the problem and then, under a smallness assumption on the data, we prove the existence of a unique weak solution to the model. We also investigate the behavior of the solution with respect the electric data on the contact surface and prove a continuous dependence result. Then, we introduce a fully discrete scheme, based on the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. We treat the contact by using a penalized approach and a version of Newton’s method. We implement this scheme in a numerical code and, in order to verify its accuracy, we present numerical simulations in the study of two-dimensional test problems. These simulations provide a numerical validation of our continuous dependence result and illustrate the effects of the conductivity of the foundation, as well.     

Inverse prediction of frictional heat flux and temperature in sliding contact with a protective strip by iterative regularization method

In this study, an inverse algorithm based on the conjugate gradient method and the discrepancy principle is applied to estimate the unknown time-dependent frictional heat flux at the interface of two semi-spaces, one of them is covered by a strip of coating, during a sliding-contact process from the knowledge of temperature measurements taken within one of the semi-space. It is assumed that no prior information is available on the functional form of the unknown heat generation; hence the procedure is classified as the function estimation in inverse calculation. Results show that the relative position between the measured and the estimated quantities is of crucial importance to the accuracy of the inverse algorithm. The current methodology can be applied to the prediction of heat generation in engineering problems involving sliding-contact elements.     

Transient response of two collinear dielectric cracks in a piezoelectric solid under inplane impacts

The paper is focused on the dynamic analysis of two collinear dielectric cracks in a piezoelectric material under the action of in-plane electromechanical impacts. Considering the dielectric permeability of crack interior, the electric displacements at the crack surfaces are governed by the jumps of electric potential and crack opening displacement across the cracks. The permeable and impermeable crack models are the limiting cases of the general one. The Laplace and Fourier transform techniques are further utilized to solve the mixed initial-boundary-value problem, and then to obtain the singular integral equations with Cauchy kernel, which are solved numerically. Dynamic intensity factors of stress, electric displacement and crack opening displacement are determined in time domain by means of a numerical inversion of the Laplace transform. Numerical results for PZT-5H are calculated to show the effects of the dielectric permeability inside the cracks, applied electric loadings and the geometry of the cracks on the fracture parameters in graphics. The observations reveal that based on the COD intensity factor, a positive electric field enhances the dynamic dielectric crack growth and a negative one impedes the dynamic dielectric crack growth in a piezoelectric solid.