A micromechanics approach to homogenizing elasto-viscoplastic heterogeneous materials

The variational asymptotic method for unit cell homogenization (VAMUCH) has emerged as a general-purpose micromechanics code capable of predicting the effective properties of heterogeneous materials and recovering the local fields. The objective of this paper is to propose a micromechanics approach enabling VAMUCH to homogenize elasto-viscoplastic heterogeneous materials. An affine formulation of the constitutive relations for an elasto-viscoplastic constituent, which exhibits viscoplastic anisotropy and combined isotropic–kinematic hardening, is derived. The weak form of the problem is derived using an asymptotic method, discretized using finite elements, and implemented into VAMUCH. The new features of VAMUCH are validated with examples such as homogenizing binary, fiber-reinforced, and particle-reinforced composites. VAMUCH is found to be capable of handling various microstructure, complex material models, complex loading conditions, and complex loading paths. More sophisticated material models can be implemented into it.     

Dynamic contact problem with adhesion and damage between thermo–electro–elasto-viscoplastic bodies

We study of a dynamic contact problem between two thermo–electro–elasto-viscoplastic bodies with damage and adhesion. The contact is frictionless and is modeled with normal compliance condition. We derive variational formulation for the model and prove an existence and uniqueness result of the weak solution. The proof is based on arguments of evolutionary variational inequalities, parabolic inequalities, differential equations, and fixed point theorem.     

Solution of the plane Hertz problem

The paper considers the problem of onesided frictionless compression of plane elastic bodies that are initially in contact with each other at a point. The first terms of an asymptotic solution of the problem are constructed by the method of joined asymptotic expansions. Determination of the approach of the bodies as a function of the pressing force reduces to calculating socalled of local compliance. The problems of contact of an elastic ring and elastic circular disks with punches and an elastic disk compressed between two elastic strips are considered. An asymptotic model for the quasistatic collision of plane elastic bodies is proposed.     

The tilted shallow wedge problem

The problem of indentation of a half-plane by a titled, shallow wedge, under frictionless conditions, is studied, using half-plane theory. The contact law itself, together with the pressure distribution and internal state of stress are found. Particular attention is then given to the nature of the pressure and state of stress at the apex, and an asymptotic form appropriate to all angular features within contacts, and therefore of widespread application, is introduced.     

Asymptotic approach to analyze singular stress state in anisotropic multi-material: Application to the rivets

The problem of stress singularities due to multi-material junctions is described in this paper. A semi-analytical approach is proposed based on an asymptotic model in the case of anisotropic linear elasticity in three dimensions. The advantage of the present method is the quasi-explicit knowledge of the stress and displacement fields around the junctions. Thus, after a brief explanation of the approach leading to the mechanical fields around the singular line at the junction between different materials, numerical results are presented concerning various configurations of an assembly with rivets included in a bi-layer composite that illustrate this method.     

A rheology model of high damping rubber bearings for seismic analysis: Identification of nonlinear viscosity

The mechanical behavior of high damping rubber bearings (HDRBs) is investigated under horizontal cyclic shear deformation with a constant vertical compressive load. On the basis of experimental observations, an elasto-viscoplastic rheology model of HDRBs for seismic analysis is developed. In this model, the Maxwell model is extended by adding a nonlinear elastic spring and an elasto-plastic model (spring-slider) in parallel. In order to identify constitutive relations of each element in the rheology model, an experimental scheme comprised of three types of tests, namely a cyclic shear (CS) test, a multi-step relaxation (MSR) test, and a simple relaxation (SR) test, are carried out at room temperature. HDRB specimens with the standard ISO geometry and three different high damping rubber materials are employed in these tests. A nonlinear viscosity law of the dashpot in the Maxwell model is deduced from the experimental scheme, and incorporated into the rheology model to reproduce the nonlinear rate dependent behavior of HDRBs. Finally, numerical simulation results for sinusoidal loading are presented to illustrate capability of the proposed rheology model in reproducing the mechanical behavior of HDRBs.     

Equivalent multipolar point-source modeling of small spheres for fast and accurate electromagnetic wave scattering computations

In this paper, we develop reduced models to approximate the solution of the electromagnetic scattering problem in an unbounded domain which contains a small perfectly conducting sphere. Our approach is based on the method of matched asymptotic expansions. This method consists in defining an approximate solution using multi-scale expansions over outer and inner fields related in a matching area. We make explicit the asymptotics up to the second order of approximation for the inner expansion and up to the fifth order for the outer expansion. We validate the results with numerical experiments which illustrate theoretical orders of convergence for the asymptotic models requiring negligible computational cost.     

Frictionless contact of two parallel congruent rigid cylindrical surfaces coated with thin elastic transversely isotropic incompressible layers: an analytic solution

The paper deals with a frictionless contact problem of two parallel rigid cylindrical surfaces, one encased in the other, coated with thin elastic transversely isotropic and incompressible layers. The coatings of the two circular cylinders may differ. A simplifying approximation for the displacement in the coating enables the problem to be formulated using stress and strain averaged through the coating thickness (for the method, see [Matthewson, M.J., 1981. Axi-symmetric contact on thin compliant coatings. J. Mech. Phys. Solids 29, 89–113]). Analytical results are obtained for the contact width and contact stress distribution. Given this contact stress distribution, an asymptotic analytical solution for the displacement in the coating is then obtained. The results are applied to the human ankle joint and generalized for articular cartilage with depth-dependent properties.     

Asymptotic approximations for Bloch waves and topological mode steering in a planar array of Neumann scatterers

We study the canonical problem of wave scattering by periodic arrays, either of infinite or finite extent, of Neumann scatterers in the plane; the characteristic lengthscale of the scatterers is considered small relative to the lattice period. We utilise the method of matched asymptotic expansions, together with Fourier series representations, to create an efficient and accurate numerical approach for finding the dispersion curves associated with Floquet–Bloch waves through an infinite array of scatterers. The approach lends itself to direct scattering problems for finite arrays and we illustrate the flexibility of these asymptotic representations on topical examples from topological wave physics.     

Response of a wire rope strand to axial and torsional loads: Asymptotic modeling of the effect of interwire contact deformations

The refined discrete mathematical model of a simple helical wire rope strand is developed. The effect of the transverse contraction of the wire strand through Poisson’s ratio and also through local contact deformations (wire flattening) has been studied in detail. In order to express the interwire contact deformation in terms of the parameters describing the strand deformation, we formulate a two-dimensional model interwire contact problem. The interwire contact interaction is treated as a frictionless unilateral plain strain problem. The nonlinear model interwire contact problem has been solved by the method of matched asymptotic expansions. The constitutive equations for a helical wire rope strand, which take into account both the Poisson’s ratio effect and the effect of contact deformation, are obtained in a closed form.     

Axisymmetric frictionless contact of functionally graded materials

The main interest of this study is a new method to solve the axisymmetric frictionless contact problem of functionally graded materials (FGMs). Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into a series of sub-layers with shear modulus varying linearly in each sub-layer and continuous at the sub-interfaces. With this model, the axisymmetric frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Hankel integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically. An erratum to this article can be found at     

INTERACTION OF THREE PARALLEL CRACKS IN RECTANGULAR PLATE UNDER CYCLIC LOADS

This paper presents a numerical approach for modeling the interaction between multiple cracks in a rectangular plate under cyclic loads. It involves the formulation of fatigue growth of multiple crack tips under ruixed-mode loading and an extension of a hybrid displacement discontinuity method （a boundary element method） to fatigue crack growth analyses. Because of an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single-region formulation. In the numerical simulation, remeshing of existing boundaries is not necessary for each increment of crack extension. Crack extension is conveniently modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. As an example, the numerical approach is used to analyze the fatigue growth of three parallel cracks in a rectangular plate. The numerical results illustrate the validation of the numerical approach and can reveal the effect of the geometry of the cracked plate on the fatigue growth.     

Dynamics of elasto-viscoplastic plates and shells

Summary In this paper, the problem of the elasto-viscoplastic dynamic behaviour of geometrically non-linear plates and shells is studied under the assumption of small strains and moderate rotations. The Chaboche and Bodner-Partom models were chosen among several types of constitutive laws. To avoid the calculation of the stiffness matrix, an effective procedure using the central difference method of solving the equations of motion was applied. The trapezoidal method was used to integrate the constitutive viscoplastic laws. A nine-node isoparametric shell element was utilised for the finite element algorithm.     

Construction of the solution of the two-dimensional problem for a complex medium with an unknown boundary

A small-parameter method is developed for solving the problem of the elasto-viscoplastic state of the material of a thick plate (plane deformation) with a hole of almost regular polygonal shape under biaxial tension. Translated from Prikladnaya Mekhanika, Vol. 34, No. 11, pp. 66–77, November, 1998.     

用于设计灵敏度分析的弹粘塑性边界元法

基于一致切线算子概念的弹粘塑性隐式边界元方法,进行了弹粘塑性设计灵敏度分析.采用了Perzyna经典粘塑性本构模型,针对包含各向同性硬化和运动硬化的混合硬化模型,导出了弹粘塑性灵敏度分析的径向返回算法和一致切线算子.利用直接微分的方法,建立了设计灵敏度分析的弹粘塑性边界元增量方程,导出了弹粘塑性径向返回的灵敏度公式.给出了在不同粘塑性流动参数下的三个典型算例的结果,与ANSYS结果相吻合,证明了方法是正确的.     

Calculation for path-domain independentJ integral with elasto-viscoplastic consistent tangent operator concept-based boundary element methods

This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals. The project supported by National Natural Science Foundation of China (9713008) and Zhejiang Natural Science Foundation Special Funds No. RC.9601     

The equilibrium state of a structure subject to frictional contact

A structure in frictional contact subject to static loads has not, in general, a unique static equilibrium state. This is because the state, displacements and contact forces, depend on the load history of the structure.In cases where the exact load history is not known it would be of interest to find a state that is in some sense likely and define this as the equilibrium state. In this paper, it is assumed that the state with the smallest potential energy is the most likely one. The implication of this definition of likely state is analysed and shows that the resulting problem basically can be seen as a generalization of the frictionless contact problem to structures where no frictionless state is possible, i.e. structures where non-zero friction forces are necessary to satisfy force equilibrium.The results of several numerical experiments are given. The structures in the experiments are trusses and structures modelled by the finite element method. Both a sequential quadratic programming method and an enumeration method are used to solve the likely-state problem.     

单侧接触约束下桁架拓扑优化的若干结果

对考虑应力和无摩擦单侧接触位移约束的桁架拓扑优化问题进行了研究,采用解析方法,通过对典型算例不同初始接触间隙下的拓扑优化解的讨论,发现这一优化问题中存在解不唯一现象,研究也表明,接触状态对结构拓扑有一定影响,不同初始接触间隙可能导致不同的拓扑结构,对利用数值方法该问题进行了研究,为克服结构拓扑优化中奇异最优解问题,采用ε松弛方法建立桁架拓扑优化模型列式,并利用双层优化方法求解本拓扑优化问题,第一层,通过求解与原接触问题等价的二次规划问题来进行结构接触分析;第二层,利用第一层得到的位移及应力,采用ε-松弛方法进行结构拓扑优化。     

Bearing capacity of strip footings with horizontal confinementLa capacité portante des fondations superficielles en présence de parois rigides

For a strip footing under axial loading, the bearing capacity is influenced by the presence of rigid walls confining the foundation soil. This problem is investigated within the framework of the theory of yield design, considering both a perfectly rough and a frictionless contact condition at the interfaces with the walls in the case of a purely cohesive soil. Upper bounds for the correction factor to be applied to the classical value of the bearing capacity are determined, as functions of the non-dimensional geometric parameter of the problem, through the kinematic approach, implementing virtual velocity fields inspired from the solution to the problem of inverted extrusion. In the perfectly rough case, it appears that the new upper bound is a significant improvement of those already available. A very simple relationship is established, which derives the upper bound for the frictionless walls from the upper bound for the rough walls. A general conclusion of the analysis is that, for the values of the geometric parameter that are likely to be encountered in practice, the increase in the bearing capacity due to the presence of the rigid walls remains very small. To cite this article: J. Salençon, C. R. Mecanique 331 (2003).     

Three-dimensional asymptotic stress field in the vicinity of the circumference of a penny shaped discontinuity

An eigenfunction expansion method is presented to obtain three-dimensional asymptotic stress fields in the vicinity of the front of a penny shaped discontinuity, e.g., crack, anticrack (infinitely rigid lamella), etc., subjected to the far-field torsion (mode III), extension/bending (mode I) and sliding shear/twisting (mode II) loadings. Five different discontinuity-surface boundary conditions are considered: (i) penny shaped crack, (ii) penny shaped anticrack or perfectly bonded thin rigid inclusion, (iii) penny shaped thin transversely rigid inclusion (frictionless planar slip permitted), (iv) penny shaped thin rigid inclusion in part perfectly bonded, the remainder with frictionless slip, and (v) penny shaped thin rigid inclusion alongside penny shaped crack. The computed stress singularity for a penny shaped anticrack is the same as that of the corresponding crack. The main difference is, however, that all the stress components at the circular tip of an anticrack depend on Poisson’s ratio under modes I and II.