Submitted:
07 December 2024
Posted:
11 December 2024
You are already at the latest version
Abstract
Keywords:
1. Summary of Eshelby’s Inclusion and Inhomogeneity Problem
1.1. Basic Definitions
- Inclusion: An inclusion is a region within a material that has different elastic properties (stiffness, for example) from the surrounding material (the matrix). The inclusion is assumed to be embedded in the matrix and can have its own distinct material properties.
- Inhomogeneity: An inhomogeneity refers to a region within a material where the material properties differ from the surrounding matrix. It is a broader term that includes inclusions but also refers to regions where properties such as density, thermal expansion, or other physical characteristics differ.
1.2. Eshelby’s Inclusion Problem
- Eshelby’s Tensor: At the core of Eshelby’s inclusion problem is the Eshelby tensor, which describes the elastic field inside and around an inclusion when it is subjected to an external stress or strain. This tensor is a fourth-order tensor that relates the applied strain to the strain inside the inclusion.
- Ellipsoidal Inclusions: Eshelby’s work showed that for ellipsoidal inclusions, the strain inside the inclusion is uniform and can be related to the external strain through the Eshelby tensor. This remarkable result simplifies the analysis of inclusions significantly, as it reduces the complexity of the problem.
- Inclusion vs. Matrix: The key idea is that when an inclusion is subjected to a uniform external stress or strain, the strain field inside the inclusion remains uniform, although different from the strain field in the surrounding matrix. The specific relationship between these strains is governed by the shape of the inclusion and the Eshelby tensor.
1.3. Mathematical Formulation
- Eigenstrain: The concept of eigenstrain (or transformation strain) is central to Eshelby’s analysis. Eigenstrain refers to a strain that would exist in the inclusion if it were isolated from the matrix and allowed to undergo a strain freely. When the inclusion is embedded in the matrix, the surrounding material restricts this strain, leading to an interaction between the inclusion and the matrix.
- Elastic Field Equations: The elastic field due to an inclusion is governed by the equations of elasticity. For a linear elastic material, these equations are linear partial differential equations (PDEs) involving the stress and strain fields, which are solved subject to boundary conditions at the inclusion-matrix interface.
- Eshelby’s Solution: Eshelby provided an analytical solution for the elastic field both inside and outside an ellipsoidal inclusion in an infinite medium. His solution showed that the strain inside the inclusion is constant and can be calculated using the Eshelby tensor.
1.4. Inhomogeneity Problem
- Difference from Inclusion: In the case of an inhomogeneity, the material properties of the region differ from those of the matrix, leading to a more complex interaction between the region and the surrounding material. Unlike an inclusion, where the material inside the inclusion can be imagined as having the same properties as the matrix, an inhomogeneity represents a real difference in material properties.
- Complexity: The solution to the inhomogeneity problem is more complex than the inclusion problem because the contrast in material properties must be accounted for. This typically requires solving the elasticity equations with variable material coefficients.
- Perturbation Techniques: In practice, solutions to inhomogeneity problems often involve perturbation techniques, where the problem is treated as a small deviation from the homogeneous case, or numerical methods, where the equations are solved using computational techniques.
1.5. Applications
- Material Science: Eshelby’s inclusion theory is widely used in materials science to predict how inclusions and inhomogeneities affect the mechanical properties of composites, polycrystals, and other heterogeneous materials.
- Micromechanics: The theory forms the basis for many micromechanical models that predict the behavior of materials with microstructural features such as voids, fibers, or precipitates.
- Fracture Mechanics: In fracture mechanics, Eshelby’s theory is used to understand how cracks and other defects influence the stress distribution in materials, which is crucial for predicting failure.
1.6. Extensions and Generalizations
- Non-Ellipsoidal Inclusions: While Eshelby’s original work focused on ellipsoidal inclusions, subsequent research has extended the theory to non-ellipsoidal shapes, though these cases generally require numerical solutions or approximations.
- Anisotropic Materials: The theory has also been extended to anisotropic materials, where the material properties differ in different directions, adding another layer of complexity to the problem.
1.7. Limitations and Challenges
- Finite Boundaries: Eshelby’s solution assumes an infinite medium, which is an idealization. In real-world applications, the finite size of the material can influence the stress and strain fields, requiring corrections or alternative methods.
- Nonlinear Materials: The theory is based on linear elasticity, and its application to nonlinear materials is limited. In such cases, more advanced models are needed.
2. Derivation of Double Derivative and Their Integration over the Surface §
3. Derivation of Quadruple Derivative and Their Integration over the Surface S
3.1. Method 1: With Using the Exchanging of Integration and Differentiation Operator
3.2. Method 2: With Using Exchanging of Integration and Differentiation Operator
- Assuming that the incremental solid angle measure is independent of the Differentiation operator variables .
- Assuming that the incremental volume measure (subtended by a cone emanating from a point inside the volume V to the surface which is at a distance in the direction from the point) is independent of the Differentiation operator variables .
- Assuming that the incremental surface area measure (subtended by a cone emanating from a point inside the volume V to the surface which is at a distance in the direction from the point) is independent of the Differentiation operator variables .
- Assuming that the incremental volume measure is independent of the Differentiation operator variables .
- Assuming that the incremental angular direction measure and incremental volume measure is independent of the Differentiation operator variables in the first integration and second integration respectively.
- Assuming that the incremental angular direction measure is independent of the Differentiation operator variables .
3.2.1. Assuming that the Incremental Solid Angle Measure is Constant with Respect to the Differentiation Operator Variables
3.2.2. Assuming that the Incremental Volume Measure is Independent of the Differentiation Operator Variables
3.2.3. Assuming that the Incremental Volume Measure is Independent of the Differentiation Operator Variables
3.2.4. Assuming that the Incremental Angular Direction Measure is Independent of the Differentiation Operator Variables
3.2.5. Assuming that the Incremental Angular Direction Measure and Incremental Volume Measure is Independent of the Differentiation Operator Variables in the First Integration and Second Integration Respectively
3.2.6. Assuming that the Incremental Surface Area Measure is Independent of the Differentiation Operator Variables
3.3. Method 3: Without Using Exchanging of Integration and Differentiation Operator
4. Betti’s Theorem and Reciprocity
- be the displacement field due to traction force and body force .
- be the displacement field due to traction force and body force .
- , i.e., a point force at .
- , i.e., a point force at .
5. Derivation of the Green’s Function for Isotropic Medium
6. Displacement due to Spontaneous Change of Form of Inclusion
6.1. Method 1: Simple Method using only Gauss Divergence Theorem
- is the stress tensor, which is constant throughout the volume.
- is the k-th component of the unit normal vector to the surface S.
- is the Kronecker delta.
- is the radial distance from the origin.
6.2. Method 2: Eshelby’s Method of using the Gauss Divergence Theorem and a special variation of Stokes’s theorem
7. Strain due to Spontaneous Change of Form of Inclusion
8. Discontinuities Across Inclusion Interface
9. Elastic Field in a Spherical and Ellipsoidal Inclusion
9.1. Spherical Inclusion
9.2. Ellipsoidal Inclusion
9.3. Elliptic Cylinder Inclusion
9.4. Flat Ellipsoid Inclusion
9.5. Penny Shaped Inclusion
10. Eshelby’s Tensor for 2-Dimensional Inclusions
11. Derivation of the Green’s Function for Anisotropic Medium
12. Derivation of the Auxiliary Tensor
13. Derivation of Inclusion Energy in an Infinite Solid
14. Derivation of Inclusion Energy in a Finite Solid
15. Derivation of Inclusion Energy of Finite solid with Applied Tractions
- State 1: Purely internal, generated by an eigenstrain (or some inhomogeneity) inside the solid.
- State 2: Purely applied, generated by external tractions on the surface of the solid without any internal eigenstrain.
- 1.
- is the stress field due to the applied tractions.
- 2.
- is the strain field corresponding to .
- 1.
- is the stress field due to the inclusion in an infinite medium.
- 2.
- is the averaged image stress field due to the boundary effects of the finite solid.
- 3.
- is the eigenstrain in the inclusion.
- 4.
- is the volume of the inclusion.
16. Ellipsoidal Inhomogenity
16.1. Application of Eshelby’s Inclusion Solution to Inhomogeneities
16.2. Transformed Inhomogeneity
16.3. Inhomogeneity under Uniform Applied Loads
- 1.
-
First Set: Suppose the solid with the inhomogeneity is subjected to a uniform strain . The stress fields inside the matrix and inhomogeneity are given by:However, this stress field does not satisfy equilibrium conditions unless a body force is applied on the surface of the inhomogeneity.
- 2.
- Second Set: To restore equilibrium, apply a body force on . The corresponding stress and strain fields due to this body force are and .
17. Cracks
17.1. Ellipsoidal Void
17.1.1. Stress and Strain Relations for an Ellipsoidal Void
17.1.2. Eigenstrain in the Void
17.1.3. Enthalpy Calculation
17.2. Penny-Shaped Crack
17.2.1. Eshelby’s Tensor for Penny-Shaped Crack
17.2.2. Eigenstrain and Stress in Penny-Shaped Crack
17.2.3. Limiting Behavior of Eigenstrain
17.3. Energy Considerations and Griffith Criterion
17.4. Slit-Like Crack
17.4.1. Derivation of Eshelby’s Tensor in the Slit-Like Crack Limit
17.4.2. Equivalent Eigenstrain
17.4.3. Griffith Criterion for Slit-Like Crack Growth
17.4.4. Stress Intensity Factors and Crack Tip Fields
17.5. Flat Ellipsoidal Crack
17.5.1. Eigenstrain Calculation
17.5.2. Enthalpy Change Calculation
17.5.3. Gibbs Free Energy and Griffith Criterion
17.6. Crack Opening Displacement: Rigorous Derivation
17.7. Stress Intensity Factors
- Mode I: Tensile mode, ,
- Mode II: In-plane shear mode, ,
- Mode III: Out-of-plane shear mode, .
17.8. Another Derivation of Crack Extension Force
17.9. J-integral as Driving Force
17.10. Invariance of J-Integral
- w is the strain energy density,
- is the traction vector,
- is the displacement field,
- is the spatial coordinate, and
- is a surface surrounding the crack tip.
17.11. Applications of J-Integral
- 1.
- On and , and , hence the contributions to J are zero.
- 2.
- On and , and , leading to zero contributions as well.
- 3.
- On , and .
18. Dislocations
18.1. Introduction to Dislocations
- is the shear modulus,
- b is the magnitude of the Burgers vector (which represents the magnitude of lattice distortion),
- a is a constant related to the atomic spacing.
18.2. Dislocation’s Effects on Mechanical Properties
18.2.1. Orowan’s Law
- is the mobile dislocation density (in units of ),
- b is the Burgers vector,
- v is the average dislocation velocity.
- Stage I: Immediately after yielding, plastic deformation occurs with little increase in applied stress. Dislocations primarily glide on parallel planes with minimal interaction.
- Stage II: At higher deformation, the slope of the stress-strain curve increases, indicating work hardening. Dislocations on several non-parallel slip planes interact, blocking each other’s motion and forming dense, entangled structures. The dislocation density increases significantly.
- Stage III: The hardening rate decreases as recovery mechanisms begin to annihilate dislocations, leading to a saturation in dislocation density.
18.3. Elastic Fields of a Dislocation Loop
18.4. Self Energy of a Dislocation Loop
18.5. Force on a Dislocation
18.6. Non-Singular Dislocation Model
References
- W. M. L. Challis, D. Rubin and E. Krempl, Introduction to Continuum Mechanics, (Butterworth-Heinemann, 1999).
- L. D. Landau and E. M. Lifshits, Theory of Elasticity, 2nd English ed., (Pergamon Press, Oxford, 1970).
- T. Mura, Micromechanics of Defects in Solids, 2nd rev. ed. (Kluwer Academic Publishers, 1991). [CrossRef]
- J. R. Barber, Elasticity, 2nd ed. (Kluwer Academic Publishers, 2002).
- L. Challis and L. Sheard, The Green of Green Functions, Physics Today, Dec. 2003.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
