Magneto elastic stress subjected to uniform magnetic field over a thin infinite plate containing an elliptical hole with an edge crack

Two dimensional solutions of the magnetic field and magneto elastic stress are presented for a magnetic material of a thin infinite plate containing an elliptical hole with an edge crack subjected to uniform magnetic field. Using a rational mapping function, each solution is obtained as a closed form. The linear constitutive equation is used for these analyses. According to the electro-magneto theory, only Maxwell stress is caused as a body force in a plate. In the present paper, it raises a plane stress state for a thin plate, the deformation of the plate thickness and the shear deflection. Therefore the magneto elastic stress is analyzed using Maxwell stress. No further assumption of the plane stress state that the plate is thin is made for the stress analysis, though Maxwell stress components are expressed by nonlinear terms. The rigorous boundary condition expressed by Maxwell stress components is completely satisfied without any linear assumptions on the boundary. First, magnetic field and stress analyses for soft ferromagnetic material are carried out and then those analyses for paramagnetic and diamagnetic materials are carried out. It is stated that those plane stress components are expressed by the same expressions for those materials and the difference is only the magnitude of the permeability, though the magnetic fields H_x, H_y are different each other in the plates. If the analysis of magnetic field of paramagnetic material is easier than that of soft ferromagnetic material, the stress analysis may be carried out using the magnetic field for paramagnetic material to analyze the stress field, and the results may be applied for a soft ferromagnetic material. It is stated that the stress state for the magnetic field H_x, H_y is the same as the pure shear stress state. Solutions of the magneto elastic stress are nonlinear for the direction of uniform magnetic field. Stresses in the direction of the plate thickness and shear deflection are caused and the solutions are also obtained. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived and investigated for the crack length.     

Magnetic field and magneto elastic stress in an infinite plate containing an elliptical hole with an edge crack under uniform electric current

Two-dimensional solutions of the electric current, magnetic field and magneto elastic stress are presented for a magnetic material of a thin infinite plate containing an elliptical hole with an edge crack under uniform electric current. Using a rational mapping function, the each solution is obtained as a closed form. The linear constitutive equation is used for the magnetic field and the stress analyses. According to the electro-magneto theory, only Maxwell stress is caused as a body force in a plate which raises a plane stress state for a thin plate and the deformation of the plate thickness. Therefore the magneto elastic stress is analyzed using Maxwell stress. No further assumption of the plane stress state that the plate is thin is made for the stress analysis, though Maxwell stress components are expressed by nonlinear terms. The rigorous boundary condition expressed by Maxwell stress components is completely satisfied without any linear assumptions on the boundary. First, electric current, magnetic field and stress analyses for soft ferromagnetic material are carried out and then those analyses for paramagnetic and diamagnetic materials are carried out. It is stated that the stress components are expressed by the same expressions for those materials and the difference is only the magnitude of the permeability, though the magnetic fields H_x, H_y are different each other in the plates. If the analysis of magnetic field of paramagnetic material is easier than that of soft ferromagnetic material, the stress analysis may be carried out using the magnetic field for paramagnetic material to analyze the stress field, and the results may be applied for a soft ferromagnetic material. It is stated that the stress state for the magnetic field H_x, H_y is the same as the pure shear stress state. Solving the present magneto elastic stress problem, dislocation and rotation terms appear, which makes the present problem complicate. Solutions of the magneto elastic stress are nonlinear for the direction of electric current. Stresses in the direction of the plate thickness are caused and the solution is also obtained. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived and investigated for the crack length and the electric current direction.     

Characterization of a class of polycrystals whose effective elastic bulk moduli can be exactly determinedCaractérisation d'une classe de polycristaux dont les modules de compression isotrope effectifs peuvent exactement être déterminés

Necessary and sufficient conditions are established for the stress response of a linearly elastic material to an isotropic stain to be hydrostatic. In the 3D case, these conditions are satisfied not only by the isotropic and cubic materials but also by all other anisotropic materials provided appropriate restrictions are imposed. In the 2D case, only the isotropic and square materials have an isotropic stress response to an isotropic strain. Using a uniform field argument, the elastic bulk modulus of a polycrystal consisting of monocrystals compatible with the established conditions is shown to equal that of any constituent monocrystal. Thus, Hill's relevant result about a polycrystal composed of cubic monocrystals is generalized. To cite this article: Q.-C. He, C. R. Mecanique 331 (2003).     

Bounds for effective elastic moduli of disordered materials

Recently P.H. Dederichs and R. Zeller (1973) have developed a formal theory of the bounds of odd order n for the effective elastic moduli of linearly elastic, disordered materials. The bounds are established by use of statistical information given in terms of correlation functions up to order n (= 1, 3, 5,…). This theory is extended to include the bounds of even order n. It is indicated how these bounds can be made optimum under the given statistical information. The results for bounds of even and odd order are obtained in forms which resemble Neumann series, containing multiple integrals up to order (n?1). These integrals can be calculated for certain materials which are classified in terms of a gradual statistical homogeneity, isotropy and disorder. Materials which possess these properties up to the correlation functions of nth order are called overall grade n materials. The optimum bounds for overall grade 2 and grade 3 materials are given explicitly. Optimum bounds for materials which are of grade ∞ in homogeneity and isotropy (i.e. (statistically) perfectly homogeneous and isotropic) and, at the same time, disordered of grade 2 or 3 are also derived. Those for grade 2 in disorder are the Z. Hashin and S. Shtrikman's (1963) bounds. Those for grade 3 are the narrowest, explicit bounds so far derived for random elastic materials. They contain within themselves the so-called self-consistent elastic moduli.     

Bounds on the distribution of extreme values for the stress in composite materials

Suitable macroscopic quantities beyond effective elastic properties are used to assess the distribution of stress within a composite. The composite is composed of N anisotropic linearly elastic materials and the length scale of the microstructure relative to the loading is denoted by ε. The stress distribution function inside the composite λ^ε(t) gives the volume of the set where the norm of the stress exceeds the value t. The analysis focuses on the case when 0<ε?1. A rigorous upper bound on lim_ε→0λ^ε(t) is found. The bound is given in terms of a macroscopic quantity called the macro stress modulation function. It is used to provide a rigorous assessment of the volume of over stressed regions near stress concentrators generated by reentrant corners or by an abrupt change of boundary loading.     

Stress intensity factor calculations based on a conservation integral

It is observed that one of the integral conservation laws of elastostatics, the so-called M-integral conservation law, has certain special features which make it possible to apply this conservation law for a class of plane elastic crack problems in order to calculate the elastic stress intensity factor in each case without solving the corresponding boundary value problem. The main characteristics which a problem must have in order for the approach to be useful are (1) for points very near to the origin of coordinates, the known elastic stresses are 0(r^?r) where r is the radial coordinate and γ ? 1, (2) for points very far from the origin, the known elastic stresses are 0(r^?r) where γ ? 1, and (3) the boundary of the body is made up of radial lines on which certain traction and/or displacement conditions are satisfied. The approach is demonstrated by determining the stress intensity factors for four familiar elastic crack problems directly from the conservation law, and then four similar additional applications of the M-integral conservation law are discussed.     

Image singularities of Green’s functions for isotropic elastic bimaterials subjected to concentrated forces and dislocations

Green’s functions for isotropic materials in the two-dimensional problem for elastic bimaterials with perfectly bonded interface are reexamined in the present study. Although the Green’s function for an isotropic elastic bimaterial subjected to a line force or a line dislocation has been discussed by many authors, the physical meaning and the structure of the solution are not clear. In this investigation, the Green’s function for an elastic bimaterial is shown to consist of eight Green’s functions for a homogeneous infinite plane. One of the novel features is that Green’s functions for bimaterials can be expressed directly by knowing Green’s functions for the infinite plane. If the applied load is located in material 1, the solution for the half-plane of material 1 is constructed with the help of five Green’s functions corresponding to the infinite plane. However, the solution for the half-plane of material 2 only consists of three Green’s functions for the infinite plane. One of the five Green’s functions of material 1 and all the three Green’s functions of material 2 have their singularities located in the half-plane where the load is applied, and the other four image singularities of material 1 are located outside the half-plane at the same distance from the interface as that of the applied load. The nature and magnitude of the image singularities for both materials are presented explicitly from the principle of superposition, and classified according to different loads. It is known that for the problem of anisotropic bimaterials subjected to concentrated forces and dislocations, the image singularities are simply concentrated forces and dislocations with the stress singularity of order O(1/r). However, higher orders (O(1/r²) and O(1/r³)) of stress singularities are found to exist in this study for isotropic bimaterials. The highest order of the stress singularity is O(1/r³) for the image singularities of material 1, and is O(1/r²) for material 2. Using the present solution, Green’s functions associated with the problems of elastic half-plane with free and rigidly fixed boundaries, for homogeneous isotropic elastic solid, are obtained as special cases.     

Quasi-static cylindrical cavity expansion in an elastoplastic compressible Mises solid

The self-similar elastoplastic field induced by quasi-static expansion of a pressurized cylindrical cavity is investigated for Mises solids under the assumption of plane-strain. Material behavior is modeled by the elastoplastic J₂ flow theory with the standard hypoelastic version. The theory accounts for elastic-compressibility and allows for arbitrary strain-hardening (or softening) in the plastic range. A formulation of the exact governing equations is presented and analyzed in detail for the remote elastic field and for asymptotic plastic behavior near the cavity wall, along with numerical investigations for the entire deformation zone. An analytical solution was obtained under the axially-hydrostatic assumption (axial stress coincides with hydrostatic stress) within an error of about 2% or less as compared to the exact, numerically evaluated, value of cavitation pressure. Two ad-hoc compressibility approximations for cavitation pressure are suggested. These relations, which give very accurate results, appear to provide tight lower and upper bounds on the exact value of cavitation pressure within an error of less than 0.5%.     

On thermoelastostatics of composites with nonlocal properties of constituents II. Estimation of effective material and field parameters

One considers a linear thermoelastic composite medium, which consists of a homogeneous matrix containing a statistically homogeneous random set of ellipsoidal uncoated or coated heterogeneities. It is assumed that the stress–strain constitutive relations of constituents are described by the nonlocal integral operators, whereas the equilibrium and compatibility equations remain unaltered as in classical local elasticity. The general integral equations connecting the stress and strain fields in the point being considered and the surrounding points are obtained. The method is based on a centering procedure of subtraction from both sides of a known initial integral equation their statistical averages obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly exploited in the known centering methods. In a simplified case of using of the effective field hypothesis for analyzing composites with one sort of heterogeneities, one proves that the effective moduli explicitly depend on both the strain and stress concentrator factor for one heterogeneity inside the infinite matrix and does not directly depend on the elastic properties (local or nonlocal) of heterogeneities. In such a case, the Levin’s (1967) formula in micromechanics of composites with locally elastic constituents is generalized to their nonlocal counterpart. A solution of a volume integral equation for one heterogeneity subjected to inhomogeneous remote loading inside an infinite matrix is proposed by the iteration method. The operator representation of this solution is incorporated into the new general integral equation of micromechanics without exploiting of basic hypotheses of classical micromechanics such as both the effective field hypothesis and “ellipsoidal symmetry” assumption. Quantitative estimations of results obtained by the abandonment of the effective field hypothesis are presented.     

Local and overall extremal properties of time-independent materials and non-linear elastic comparison materials

For time-independent materials which undergo non-linear deformations from some given reference configuration two (dual) hypotheses are considered. Firstly it is supposed that the work done to a given state of deformation is bounded below and that the bound is attainable on a physically possible path; secondly that the complementary work to a given state of stress is bounded above and that this bound too is attainable on a physically possible path. The consequences of these assumptions are analysed, and the results of Ponter and Martin [1] in the linear theory are generalized to account for non-linear deformations, due attention being paid to questions of stability.A non-linear elastic comparison material is defined whose strain energy is equal to the work done on a minimum path for the time-independent material. Extremum principles for non-linear elastic materials given in [2] are then applied to the comparison elastic material, and bounds are thereby placed on the work and complementary-work functional of the time-independent material. Corresponding overall properties of the time-independent and elastic materials are then compared by defining respective overall constitutive laws and overall stress and deformation variables.Following the definition of strengthening (weakening) of a non-linear elastic solid given by Ogden[2] a time-independent material is said to be strengthened (weakened) when its comparison elastic material is strengthened (weakened). Local and overall aspects of this definition are examined.     

Three-phase plane composites of minimal elastic stress energy: High-porosity structures

The paper establishes exact lower bound on the effective elastic energy of two-dimensional, three-material composite subjected to the homogeneous, anisotropic stress. It is assumed that the materials are mixed with given volume fractions and that one of the phases is degenerated to void, i.e., the effective composite is porous. Explicit formula for the energy bound is obtained using the translation method enhanced with additional inequality expressing certain property of stresses. Sufficient optimality conditions of the energy bound are used to set the requirements which have to be met by the stress fields in each phase of optimal effective material regardless of the complexity of its microstructural geometry. We show that these requirements are fulfilled in a special class of microgeometries, so-called laminates of a rank. Their optimality is elaborated in detail for structures with significant amount of void, also referred to as high-porosity structures. It is shown that geometrical parameters of optimal multi-rank, high-porosity laminates are different in various ranges of volume fractions and anisotropy level of external stress. Non-laminate, three-phase microstructures introduced by other authors and their optimality in high-porosity regions is also discussed by means of the sufficient conditions technique. Conjectures regarding low-porosity regions are presented, but full treatment of this issue is postponed to a separate publication. The corresponding “G-closure problem” of a three-phase isotropic composite is also addressed and exact bounds on effective isotropic properties are explicitly determined in these regions where the stress energy bound is optimal.     

Porous materials with two populations of voids under internal pressure: I. Instantaneous constitutive relations

This study is devoted to the mechanical behavior of uranium dioxide (UO₂) which is a porous material with two populations of voids of very different size subjected to internal pressure. The smallest voids are intragranular and spherical in shape whereas the largest pores located at the grain boundary are ellipsoidal and randomly oriented. In this first part of the study, attention is focused on the effective properties of these materials with fixed microstructure. In a first step, the poro-elastic properties of these doubly voided materials are studied. Then two rigorous upper bounds are derived for the effective poro-plastic constitutive relations of these materials. The first bound, obtained by generalizing the approach of Gologanu et al. (Gologanu, M., Leblond, J., Devaux, J., 1994. Approximate models for ductile metals containing non-spherical voids-case of axisymmetric oblate ellipsoidal cavities. ASME J. Eng. Mater. Technol. 116, 290–297) to compressible materials, is accurate at high stress-triaxiality. The second one, which derives from the variational method of Ponte Castañeda (Ponte Castañeda, P., 1991. The effective mechanical properties of non-linear isotropic composites. J. Mech. Phys. Solids 39, 45–71), is accurate when the stress triaxiality is low. A N-phase model, inspired by Bilger et al. (Bilger, N., Auslender, F., Bornert, M., Masson, R., 2002. New bounds and estimates for porous media with rigid perfectly plastic matrix. C.R. Mécanique 330, 127–132), is proposed which matches the best of the two bounds at low and high triaxiality. The effect of internal pressures is discussed. In particular it is shown that when the two internal pressures coincide, the effective flow surface of the saturated biporous material is obtained from that of the drained material by a shift along the hydrostatic axis. However, when the two pressures are different, the modifications brought to the effective flow surface in the drained case involve not only a shift along the hydrostatic axis but also a change in shape and size of the surface.     

A micromechanical iterative approach for the behavior of polydispersed composites

In this paper, an iterative homogenization method is proposed in order to predict the behavior of polydispersed materials. Various families of heterogeneities according to their geometrical or mechanical properties are progressively introduced into a volume of matrix. At each step, the behavior of intermediate medium is obtained by any analytical homogenization method and is used as matrix of the following step. All homogenization methods, like dilute strain or stress approximations, Hashin’s bounds, three phases method, Mori–Tanaka’s approach or for example the N-layered inclusions method lead to the same effective behavior for the polydispersed material after convergence of the iterative process. Moreover, this convergence is obtained even for significant fractions of heterogeneities and for highly contrasted or polydispersed materials. This method is applied to various composites and validated by comparison with other modellings and experimental results.     

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

Computing lower and upper bounds on stress intensity factors in bimaterials

This paper presents a finite element approach for finding complementary bounds of stress intensity factors (SIFs) in bimaterials. The SIF is formulated as an explicit computable linear function of displacements by means of the two-point extrapolation method. An appropriate and computable form of the SIF plays a crucial role in the dual problem involved in the computing procedure of both lower and upper bounds. In our discussions, computable forms of stress intensity factors, K₀ and K_r, are derived, which are related to the energy release rate, and the ratio of the open mode and shear mode SIFs, respectively. Based on a posteriori finite element error estimation, a bounding procedure is used to compute the bounds on the two stress intensity factors. Finally, bounds on the SIFs in a bimaterial interface crack problem are provided to verify the procedure.     

Bounds on the elastic moduli of statistically isotropic multicomponent materials and random cell polycrystals

Minimum energy and complementary energy principles are used to derive the upper and lower bounds on the effective elastic moduli of statistically isotropic multicomponent materials in d ($d=2$ or 3) dimensions. The trial fields, involving harmonic and biharmonic potentials, and free parameters to be optimized, lead to the bounds containing, in addition to the properties and volume proportions of the material components, the three-point correlation information about the microgeometries of the composites. The relations and restrictions among the three-point correlation parameters are explored. The upper and lower bounds are specialized to symmetric cell materials and asymmetric multi-coated spheres, which are optimal or even converge in certain cases. New bounds for random cell polycrystals are constructed with particular results for random aggregates of cubic crystals.     

Stress state of an elastic ball in a spherically symmetric gravitational field

We consider a spherically symmetric static problem of general relativity whose solution was obtained in 1916 by Schwarzschild for a metric form of a special type. This solution determines the metric coefficients of the exterior and interior Riemannian spaces generated by a gravitating solid ball of constant density and includes the so-called gravitational radius r _g. For a ball of outer radius R=r _g, the metric coefficients are singular, and hence the radius r _g is traditionally assumed to be the radius of the event horizon of an object called a black hole. The solution of the interior problem obtained for an incompressible ideal fluid shows that the pressure at the ball center increases without bound for R=9/8r _g, which is traditionally used for the physical justification of the existence of black holes. The discussion of Schwarzschild’s traditional solution carried out in this paper shows that it should be generalized with respect to both the geometry of the Riemannian space and the elastic medium model. In this connection, we consider the general metric form of a spherically symmetric Riemannian space and prove that the solution of the corresponding static problem exists for a broad class of metric forms. A special metric form based on the assumption that the gravitation generating the Riemannian space inside a fluid ball or an elastic ball does not change the ball mass is singled out from this class. The solution obtained for the special metric form is singular with respect to neither the metric coefficients nor the pressure in the fluid ball and the stresses in the elastic ball. The obtained solution is compared with Schwarzschild’s traditional solution.     

Matched asymptotic expansions in nonlinear fracture mechanics—II. Longitudinal shear of an elastic work-hardening plastic specimen

Earlier analysis given by T.M. Edmunds and J.R. Willis (1976) is extended to deal with cracks in elastic work-hardening plastic specimens subjected to longitudinal shear loads. Solutions are expressed in terms of a set of parameters that are determined from linear elastic solutions alone. It is proved, for any specimen geometry and any loading symmetric about the plane of the crack, that a ‘plastic-zone correction’, obtained by solving a linear elastic problem for a crack which is a length r_y longer than the actual crack, provides a two-term asymptotic expansion for the J-integral, if r_y is defined suitably in terms of the linear elastic stress concentration factor and the initial slope of the work-hardening curve. The general method is applied in detail to a strip of finite width containing an edge crack, for which the effect of the work-hardening on the maximum extent of the plastic zone and on the J-integral is summarized graphically.     

A finite element investigation of the effect of crack tip constraint on hole growth under mode i and mixed mode loading

In this work, the effect of constraint on hole growth near a notch tip in a ductile material under mode I and mixed mode loading (involving modes I and II) is investigated. To this end, a 2-D plane strain, modified boundary layer formulation is employed in which the mixed mode elastic K–T field is prescribed as remote boundary conditions. A finite element procedure that accounts for finite deformations and rotations is used along with an appropriate version of J₂ flow theory of plasticity with small elastic strains. Several analyses are carried out corresponding to different values of T-stress and remote elastic mode-mixity. The interaction between the notch and hole is studied by examining the distribution of hydrostatic stress and equivalent plastic strain in the ligament between the notch tip and the hole, as well as the growth of the hole. The implications of the above results on ductile fracture initiation due to micro-void coalescence are discussed.