Unified elastoplastic finite difference and its application

Two elastoplastic constitutive models based on the unified strength theory (UST) are established and implemented in an explicit finite difference code, fast Lagrangian analysis of continua (FLAC/FLAC3D), which includes an associated/non-associated flow rule, strain-hardening/softening, and solutions of singularities. Those two constitutive models are appropriate for metallic and strength-different (SD) materials, respectively. Two verification examples are used to compare the computation results and test data using the two-dimensional finite difference code FLAC and the finite element code ANSYS, and the two constitutive models proposed in this paper are verified. Two application examples, the large deformation of a prismatic bar and the strain-softening behavior of soft rock under a complex stress state, are analyzed using the three-dimensional code FLAC3D. The two new elastoplastic constitutive models proposed in this paper can be used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials. The effect of the intermediate principal stress on metallic or SD material structures under complex stress states, including large deformation, three-dimensional and non-association problems, can be analyzed easily using the two constitutive models proposed in this paper.     

Electromechanical behaviour of a finite piezoelectric layer under a flat punch

This paper solves the problem of a smooth and frictionless punch on a piezoelectric ceramic layer. Different electrical boundary conditions that employ conducting or insulating punches are presented. The stress and electric displacement intensity factors are used to characterize the electromechanical fields at the punch tip. The field intensity factors are obtained numerically for finite layer thickness. Effects of the thickness of the piezoelectric layer on the stress and electric displacement, and the stress and electric displacement intensity factors at the punch tip are discussed. Solution technique for two identical and collinear surface punches on the piezoelectric layer is also provided and the effect of relative distance between the punches is investigated. Numerical results for some interesting special cases, such as conducting punch and insulating punch, and infinite piezoelectric layer thickness, are presented.     

Elastic analysis of some punch problems for a layered medium

The problems of flat-ended cylindrical, quadrilateral, and triangular punches indenting a layered isotropic elastic half-space are considered. The former two are analyzed using a basis function technique, while the latter problem is analyzed via a singular integral equation. Solutions are obtained numerically. Load-deflection relations are obtained for a series of values of the ratio of Young's modulus in the layer and substrate, and for a variety of punch sizes. These solutions provide an accurate basis for the estimation of Young's modulus of thin films from the initial unloading compliance observed in indentation tests, and are specifically relevant to axisymmetric, Vicker's, and triangular indenters. The results should also be of interest in foundation engineering.     

Axi-symmetric bodies of equal material in contact under torsion or shift

Summary Two axi-symmetric bodies are pressed together, so that their axes of symmetry coincide with the contact normal and the normal force is held constant. A small torque about the contact normal or a small tangential force is applied. For bodies of equal material, the normal and tangential stress states are uncoupled, and can solved separately. The surfaces of the bodies are thought as a superposition of infinitesimal rigid flat-ended punches. Consequently, the normal stress distribution can be calculated as a summation of differential flat punch solutions. A formula results, which is identical with the solution of Green and Collins. After application of a torque an annular sliding area forms at the border of the contact area. For reasons of symmetry, the common displacement of the inner stick area must be a rigid body rotation. Similarly to the normal problem, the solution can be thought as a superposition of rigid punch rotations. The tangential solution can be derived analogically, in form of a superposition of rigid punch displacements. The present method also solves the problem of simultanous normal and torsional or tangential loading with complete adhesion. As an example, Steuermann's problem for polynomial surfaces of the formA _2nr²ⁿis solved. The solutions for constant normal forces can be used as basic functions for loading histories with varying normal and tangential forces.     

Self-similar problems of elastic contact for non-convex punches

Self-similar problems of contact for non-convex punches are considered. The non-convexity of the punch shapes introduces differences from the traditional self-similar contact problems when punch profiles are convex and their shapes are described by homogeneous functions. First, three-dimensional Hertz type contact problems are considered for non-convex punches whose shapes are described by parametric-homogeneous functions. Examples of such functions are numerous including both fractal Weierstrass type functions and smooth log-periodic sine functions. It is shown that the region of contact in the problems is discrete and the solutions obey a non-classical self-similar law. Then the solution to a particular case of the contact problem for an isotropic linear elastic half-space when the surface roughness is described by a log-periodic function, is studied numerically, i.e. the contact problem for rough punches is studied as a Hertz type contact problem without employing additional assumptions of the multi-asperity approach. To obtain the solution, the method of non-linear boundary integral equations is developed. The problem is solved only on the fundamental domain for the parameter of self-similarity because solutions for other values of the parameter can be obtained by renormalization of this solution. It is shown that the problem has some features of chaotic systems, namely the global character of the solution is independent of fine distinctions between parametric-homogeneous functions describing roughness, while the stress field of the problem is sensitive to small perturbations of the punch shape.     

Two Kinds of Contact Problems in Dodecagonal Quasicrystals of Point Group 12 mm

In this paper,two kinds of contact problems in 2-D dodecagonal quasicrystals were discussed using the complex variable function method:one is the finite frictional contact problem,the other is the adhesive contact problem.The analytic expressions of contact stresses in the phonon and phason fields were obtained for a flat rigid punch,which showed that:(1) for the finite frictional contact problem,the contact stress exhibited power-type singularities at the edge of the contact zone;(2) for the adhesive contact problem,the contact stress exhibited oscillatory singularities at the edge of the contact zone.The distribution regulation of contact stress under punch was illustrated;and the low friction property of quasicrystals was verified graphically.     

Development of a solution method for frictional contact problems with complex loading

In this work, solution methods for frictional contact problems are extended to the case of moving punches and to the external loading history-dependent system states. To solve the frictional contact problems in the contact area, an iterative method is developed and implemented. Solutions of two-dimensional problems are constructed using the boundary element method. Numerical analysis is aimed at the quantitative study of effects such as the interaction of contact pressure and friction forces, estimates of the friction force differences due to the differences in the choice of local basis for the calculation of normal pressure and friction forces, and evaluation of the effects of complex loading (rotation of the rigid punch after its preliminary penetration into the solid). We find that, for the same definition of the friction force, different initial approximations lead to the same solution. At the same time, the friction forces defined either as projections onto the common tangent plane or as projections onto the plane tangent to the punch can differ quite substantially. Similar conclusions are derived for the solutions corresponding to single or multiple loading steps. The work relies on the variational principle for the solution of contact problems and numerical algorithms developed for the problems with one-sided constraints. The variational principle was first applied by Signorini [1] to the determination of the stress-strain state in a linearly deformed body in a rigid smooth shell. The modern view of the problem and its generalizations to the frictional problems and some other problems involving unilateral constraints in given in the monograph [2]. Finite difference and finite element methods in application to the problems with unilateral constraints are described in [3]. Analytical solution methods are developed in the monographs [4–6].     

闭孔泡沫铝的高温局部压入力学响应

通过不同形状(平头和半球头)的压头在不同温度下对闭孔泡沫铝材料进行塑性压入实验,研究不同温度下闭孔泡沫铝的压入变形模式及载荷响应特性。并基于闭孔泡沫铝在高温下的准静态塑性压入载荷响应的实验结果,结合多种分析方法,(如量纲分析和有限元计算等),探索既考虑温度影响也包含压入深度影响的预测闭孔泡沫铝平头和半球头压入力学响应的经验公式。结果表明,本文得到的两种压头情况下的经验公式都能够较好地预测闭孔泡沫铝在不同温度下的压入力学响应。     

On indentation of a pair of rigid punches connected by an elastic beam into an elastic half-plane with regard to the friction and adhesion forces in the contact region

The contact problem of indentation of a pair of rigid punches with plane bases connected by an elastic beam into the boundary of an elastic half-plane is considered under the conditions of plane strain state. The external load is generated by lumped forces applied to the punches and a uniformly distributed normal load acting on the beam.It is assumed that the contact between the punch and the elastic half-plane can be described by L. A. Galin’s statement, i.e., it is assumed that the adhesion acts in the interior part of each of the contact regions and the tangential stresses obeying the Coulomb law act on their boundaries.With the symmetry taken into account, the problem is stated only for a single punch, and solving this problem is reduced to a system of four singular integral equations for the tangential and normal stresses in the adhesion region and the contact pressure in the sliding zones. The solution of the constitutive system together with three conditions of equilibrium of the system of punches connected by a beam is constructed by direct numerical integration by the method of mechanical quadratures.As a result of the numerical analysis, the contact stress distribution functions were constructed and the values of the sliding zones and the punch rotation angle were determined for various values of the geometric, elastic, and force characteristics.     

A unified treatment of axisymmetric adhesive contact problems using the harmonic potential function method

A unified treatment of axisymmetric adhesive contact problems is provided using the harmonic potential function method for axisymmetric elasticity problems advanced by Green, Keer, Barber and others. The harmonic function adopted in the current analysis is the one that was introduced by Jin et al. (2008) to solve an external crack problem. It is demonstrated that the harmonic potential function method offers a simpler and more consistent way to treat non-adhesive and adhesive contact problems. By using this method and the principle of superposition, a general solution is derived for the adhesive contact problem involving an axisymmetric rigid punch of arbitrary shape and an adhesive interaction force distribution of any profile. This solution provides analytical expressions for all non-zero displacement and stress components on the contact surface, unlike existing ones. In addition, the newly derived solution is able to link existing solutions/models for axisymmetric non-adhesive and adhesive contact problems and to reveal the connections and differences among these solutions/models individually obtained using different methods at various times. Specifically, it is shown that Sneddon’s solution for the axisymmetric punch problem, Boussinesq’s solution for the flat-ended cylindrical punch problem, the Hertz solution for the spherical punch problem, the JKR model, the DMT model, the M-D model, and the M-D-n model can all be explicitly recovered by the current general solution.     

A Mixed Elastic Problem for an Isotropic Half-Plane with Circular Openings

A mixed problem is solved for a multiply connected half-plane with circular openings. Punches rigidly mated to the half-plane act on the rectilinear boundary. By using an analytic continuation through the unloaded parts of the rectilinear boundary and solving the obtained linear-conjunction problem for the slits of the multiply connected domain, the general representation of the complex potential containing unknown functions is found. These functions are holomorphic outside the openings and determined from the boundary conditions on the opening periphery and some additional equilibrium conditions for the punches. The indicated boundary conditions are satisfied with the help of the least-squares method. In the case where a punch acts on the boundary of a half-plane with one opening, the effects of the punch width and the relative position of the punch and the opening on the stress concentration and distribution are numerically evaluated     

Conformal contact between layered foundations and punches

We study the contact interaction between rigid punches and viscoelastic foundations with thin coatings for the cases in which the punch and coating surfaces are conformal (mutually repeating). Such problems can arise, for example, when the punch immerses into a solidificating coating before its complete solidification; as a result, the surface takes the shape of the punch base. Examples of such coatings can be a layer of glue, concrete at its young age, many polymeric materials. We consider plane contact problems for inhomogeneous aging viscoelastic basements in the case of their conformal contact with rigid punches. We present the statements of the problems and derive their basic mixed integral equation. The solution of this equation is constructed by using the generalized projection method. We present numerical computations of model problems, including the problem in which the shape of the punch base is described by a rapidly oscillating function.     

A four-part boundary value problem in elasticity — indentation by a discontinuously concave punch

A solution is given for the frictionless indentation of an elastic half-space by a flat-ended cylindrical punch with a central circular recess, when the load is large enough to establish a circular region of contact in the recess. The problem is reduced to two simultaneous Fredholm equations using the method of complex potentials due to Green and Collins. Results are presented for the relationship between load, contact radius and penetration for various punch geometries.     

Thermoelastic State of a Transversely Isotropic Layer between Two Annular Punches

A technique is developed to solve contact problems for annular punches interacting with a transversely isotropic layer. The contact problem for two heated annular punches interacting with a layer is solved. The formulas for the contact stresses under the punches are derived, and the effect of the shape of the punches on the magnitude and distribution of these stresses is analyzed     

Galin’s problem for a spatial elastic wedge

We study a three-dimensional contact problem on the indentation of an elliptic punch into a face of a linearly elastic wedge. The wedge is characterized by two parameters of elasticity and its edge is subjected to the action of an additional concentrated force. The other face wedge is free from stresses. The problem is reduced to an integral equation for the contact pressure. An asymptotic solution of this equation is obtained which is effective for a given contact region fairly remote from the edge. Calculations are performed that allow one to evaluate the effect of a force applied outside the contact region on the contact pressure distribution. The problem under study is a generalization of L. A. Galin’s problem on a force applied outside a circular punch on an elastic half-space [1, 2]. In a special case of a wedge with an opening angle of 180° and zero contact ellipse eccentricity, the obtained asymptotic relation coincides with the expansion of Galin’s exact solution in a series. Problems of indentation of an elliptic punch into a spatial wedge with the face not loaded outside the contact region have been studied previously. For example, the paper [3] dealt with the case of a known contact region (asymptotic method) and the paper [4] considered the case of an unknown contact region (numerical method). The solution of Galin’s problem allowed the authors of [2] to reduce the contact problem on the interaction of several punches applied to a half-space to a system of Fredholm integral equations of the second kind (Andreikin-Panasyuk method). A topical direction in contact mechanics is the model of discrete contact as well as related problems on the interaction of several punches [2, 5–8]. The interaction of several punches applied to a face of a wedge can be treated in a similar manner and an asymptotic solution can be obtained for the case where a concentrated force is applied at an arbitrary point of this face beyond the contact region rather than on the edge.     

Uniform indentation of the elastic half-space by a rigid rectangular punch

The problem considered is that of a rigid flat-ended punch with rectangular contact area pressed into a linear elastic half-space to a uniform depth. Both the lubricated and adhesive cases are treated. The problem reduces to solving an integral equation (or equations) for the contact stresses. These stresses have a singular nature which is dealt with explicitly by a singularity-incorporating finite-element method. Values for the stiffness of the lubricated punch and the adhesive punch are determined: the effect of adhesion on the stiffness is found to be small, producing an increase of the order of 3%.     

Indentation of punches into a piezoceramic body: two-dimensional contact problem of electroelasticity

The paper establishes the relationship between the solutions of the static contact problems of elasticity (no friction) for an isotropic half-plane and problems of electroelasticity for a transversely isotropic piezoelectric half-plane with the boundary perpendicular to the polarization axis. This allows finding the contact characteristics in the electroelastic case from the known elastic solution, without the need to solve the electroelastic problem. The contact problems of electroelasticity for different types of wedge-shaped punches (flat punch with rounded one or two edges, half-parabolic punch, and a periodic system of punches) are solved as examples Translated from Prikladnaya Mekhanika, Vol. 44, No. 11, pp. 55–70, November 2008.     

A generalized velocity field for plane strain backward extrusion through punches of any shape

In this paper, the process of plane strain backward extrusion process through arbitrarily curved punches, by means of the upper bound method and the finite element method is investigated. A generalized velocity field is developed and by calculating of the internal, shear and frictional powers, the extrusion force is estimated. Then, by using the developed analytical model, optimum punch lengths which minimize the required extrusion forces, are determined for a wedge shaped punch and a streamlined punch shape. The corresponding results for those two punch shapes are also determined by using a finite element code and compared with the upper bound results. This comparison shows that the upper bound predictions are in good agreement with the FE results.     

Elliptical crack and punch problems solved by integral equation method for piezoelectric media

The problem of an elliptical crack embedded in an unbounded transversely isotropic piezoelectric media with the crack-plane parallel to the plane of isotropy of the media and subjected to remote normal mechanical as well as electric loading is considered first. The problem has been successfully reduced to a pair of coupled integral equations that are suitable for the application of an integral equation method developed earlier for three-dimensional problems of LEFM. Solution to the mechanical displacement and electric potentials are obtained for prescribed uniform loadings and expressions for corresponding intensity factors and crack opening displacement are deduced. The above method has further been applied to solve the problem of a rigid flat-ended elliptical punch indenting a transversely isotropic piezoelectric half-space surface with the plane of isotropy parallel to the surface. Solutions to mechanical stress and electric displacement are obtained for prescribed constant normal displacement and constant electric potential interior to the elliptical region and expression for the total force required to maintain a prescribed indentation is deduced.     

Interaction of an Ellipsoidal Inclusion with an Elliptic Crack in an Elastic Material under Triaxial Tension

The interaction of an elastic ellipsoidal inclusion with an elliptic crack in an infinite elastic medium under triaxial loading is analyzed. The stress state in the elastic space is represented as a superposition of the principal state and perturbed states, which are due to the presence and interaction of the inclusion and the crack. The analytical solution of the problem is found using the method of equivalent inclusion, the potential of an inhomogeneous ellipsoid, and a system of harmonic functions for an elliptic crack. The effect of triaxial loading on the stress intensity factors is analyzed