Abstract
Micromorphic theory (MMT) envisions a material body as a continuous collection of deformable particles; each possesses finite size and inner structure. It is considered as the most successful top-down formulation of a two-level continuum model, in which the deformation is expressed as a sum of macroscopic continuous deformation and internal microscopic deformation of the inner structure. In this work, the kinematics including the objective Eringen tensors is introduced. Balance laws are derived by requiring the energy equation to be form-invariant under the generalized Galilean transformation. The concept of material force and the balance law of pseudomomentum are generalized for MMT. An axiomatic approach is demonstrated in the formulation of constitutive equations for a generalized micromorphic thermoviscoelastic solid, generalized micromorphic fluid, micromorphic plasticity, and micromorphic electromagnetic–thermoelastic solid. Applications of MMT in micro/nanoscale are discussed.
1. Introduction
Continuum approaches have dominated material modeling research over the past few decades. This approach to predict material deformation and failure, by implicitly averaging atomic scale dynamics and defect evolution over time and space, is valid only for large systems [Citation1–3 Citation Citation3]. Therefore, lots of experimental observations of material behavior, especially at micro/nanoscale, cannot be explained within the continuum mechanics framework. On the other hand, molecular dynamics (MD) simulations [Citation4] have become a powerful tool for elucidating complex mechanics phenomena. But the length and time scales that can be probed in MD are still fairly limited. The pros and cons of both continuum and atomistic models have motivated the development of various methods of microscale modeling and simulation. Microcontinuum field theories, or generalized continuum theories, including Cosserat theory [Citation5], couple stress theory [Citation6], micromorphic theory [Citation7], microstructure theory [Citation8], micropolar theory [Citation9] are extensions of the classical field theory to microscopic space and time scales. Here let us briefly mention the relation of micromorphic theory to several others microscale theories. As a pioneer of the rational theories of polar continua, Cosserat obtained balance equations of momenta in the dynamic case. However, Cosserat theory is a special case of micropolar theory, where the balance law of microinertia is missing. Couple stress theory, by including high order stresses, provided a model that can support body and surface couples. However, it can be obtained from micropolar theory as a special case when the motion is constrained so that the macrorotations and microrotations coincide. Upon the assumptions of small deformation, slow motion, constant microinertia, spin isotropy and linear isotropic elasticity, micromorphic theory can be reduced to microstructure theory. However, those assumptions lead to a limitation in applications related to atomic motions, such as the thermal mechanical coupling, phase transition, anisotropic and large deformation problems. In addition, the local balance equations of energy, entropy, and microinertia are not given in microstructure theory. The above-mentioned discussion lends the credible support to claim that micromorphic theory is the most successful formulation of a microcontinuum model.
Micromorphic theory has been developed since 1964 by Eringen. It is considered to be the most successful top-down microscale model. The development of micromorphic theory was motivated by the query, “Is it possible to construct continuum theories that can predict physical phenomena on the atomic, molecular, or nano scales?” These would require supplying additional degrees of freedom to a material point, which has three degrees of freedom in classical field theory. After all, the molecules that constitute the internal structure of the material points undergo deformations and rotations arising from the displacements and rotations of their constituent atoms. Micromorphic theory envisions a material body as a continuous collection of deformable particles; each possesses finite size and inner structure. On the other hand, classical continuum mechanics considers a material body as a continuous collection of material points, each with infinitesimal size and no inner structure. The purpose of going beyond classical continuum mechanics is to take into account the microstructure of the material body in question while still keeping the advantages of continuum theory intact. A question arises: how can we reconcile the concept of the deformable particle with the continuum hypothesis? Eringen settled this question by replacing the deformable particle with a geometric point and some vectors attached to that point, which denote the orientations and intrinsic deformations of all the material points in the deformable particle [Citation7,Citation10–18 Citation Citation Citation Citation Citation Citation Citation Citation18]. This is compatible with the classical picture where a material point in a continuum is endowed with physical properties such as mass density, displacement vector, electric field, stress tensor, etc. Therefore micromorphic theory can be expected to unveil many new classes of physical phenomena that fall beyond classical field theories.
The organization of the remainder of this paper is as follows. In Section 2, we briefly present the kinematics of micromorphic theory. Balance laws and material force are derived in Section 3. Section 4 presents the axiomatic approach to formulate constitutive equations of generalized Micromorphic solid/fluid. In Sections 5 and 6, we derive the constitutive theories for micromorphic plasticity and micromorphic thermomechanics–electromagnetics, respectively. Finally, we conclude this paper with discussions in Section 7.
2. Kinematics
Micromorphic theory, developed by Eringen and Suhubi [Citation10,Citation11] and Eringen [Citation12,Citation13], constitutes extensions of the classical field theories concerned with the deformations, motions, and electromagnetic interactions of material media, as continua, in microscopic time and space scales. In micromorphic theory, a material body is considered as a continuous collection of deformable particles, each with finite size and inner structure. Geometrically, a deformable particle P is characterized by its centroid C located at X and vector Ξ relative to X of a generic point within the particle. As shown in , a generic point in the particle is represented by the vector sum of X and Ξ in the reference state at time t = 0. The material is called a micromorphic continuum if the deformation carrying P(X,Ξ,t) to p(x,ξ,t) in the deformed state at time t can be expressed as
Then there exist unique inverse motions
with
If the microgyration tensor is defined as
Let and ΔV(Δv) denote the mass density and the volume of the deformable particle in the Lagrangian (Eulerian) state and let primed quantities refer to those of the point in the particle. This leads to
2.1. Objectivity
It is intuitively clear that the material properties do not depend on the coordinate frame selected. The measurements made by an observer, whether they are in motion or not, should be the same. If this viewpoint is accepted, then the measurements made in one frame of reference are sufficient to determine the material properties in all other frames, which are in rigid motion with respect to one another. In the formulation of the response functions, it is desirable to employ quantities that are not dependent on the motions of the observer. Such quantities are called objective or material frame-indifferent.
2.1.1. Definition 1
Two motions and
are called objectively equivalent if and only if
2.1.2. Definition 2
Any tensorial quantity is said to be objective if in any two objectively equivalent motions it obeys the following tensor transformation law for all times:
Since EquationEquation (11)(11) in general holds for arbitrary X and Ξ, it can be replaced by
The generalized Lagrangian strain tensors of micromorphic theory are defined as [Citation12]
One may verify that the strain rates can be obtained as
2.2. Eringen tensors
There are three kinds of Eringen tensors; the first order Eringen tensors are ,
, and
. Eringen tensors of order n + 1 are defined as
After lengthy but straightforward derivation one may prove that Eringen tensors of any order are objective. One may also verify that
3. Balance laws
Eringen and Suhubi [Citation10,Citation11] and Eringen [Citation7] derived the laws of conservation of mass, conservation of microinertia, balance of linear momentum, balance of momentum moments, and conservation of energy for micromorphic theory by means of a “microscopic space-averaging” process. Later Eringen [Citation12] derived the balance laws in a more elegant way by starting with the following expression for the kinetic energy per unit mass, k, of a particle:
The body force and the body moment are defined as
The moment stress, a third order tensor, is defined as
Of course, under this limiting situation (size of the particle is vanishing or, say, a particle is regarded as a mathematical point), micromorphic theory is identical to classical continuum theory.
The balance laws expressed in EquationEquations (21)(21)*
Equation
(22)*
Equation
(23)*
Equation
(24)*
Equation
(25)*–
Equation(26)
(26)* are in the Eulerian description. In Lagrangian description, the balance laws can be written as [Citation19]
3.1. Material forces
The concept of material forces was first introduced by Eshelby [Citation20], elaborated and further developed by Maugin [Citation21,Citation22]. Material forces are generated by displacement, not in physical space, but on material manifold. For example, they can be generated by (i) an infinitesimal rigid displacement of a finite region surrounding a point of singularity in an elastic body [Citation20], (ii) an infinitesimal displacement of a dislocation line [Citation23], (iii) an infinitesimal increase in the length of a crack [Citation24–26 Citation Citation26]. Material forces drive the motion of defects of various dimensions in condensed matter physics, e.g. phase-transition fronts in elasticity, Bloch and Neel walls in ferromagnetism, and elastic solitons [Citation21,Citation26]. This characteristic property of material forces also leads to their christening as inhomogeneity forces. Material inhomogeneity is defined as the dependence of properties (not the solution), such as density, elastic coefficients, viscosity, plasticity threshold, on the material point. These inhomogeneities may be more or less continuous such as in metallurgically superficially treated specimens or in a polycrystal observed at a mesoscopic scale, or it may change abruptly such as in laminated composite or in a body with foreign inclusions or cavities.
Eshelbian mechanics was extended to micromorphic theory. The formulation can be easily done by adding (a) EquationEquation (23)(23)* multiplied by xl
,p
and (b) EquationEquation (24)
(24)* multiplied by
. It leads to the balance law of pseudo-momentum,
F
2 is the material force due to the inhomogeneity of and
,
K is the kinetic energy density (per unit volume in Lagrangian coordinates)
P is named as pseudo-momentum:
The detailed expressions of Eshelby stress tensor, pseudo-momentum, and material forces were derived for thermoelastic micromorphic solid. It was found that the material forces are due to (1) body force and body moment, (2) temperature gradient and (3) material inhomogeneities in density, microinertia, and elastic coefficients. The general expression of material forces due to the presence of dynamically propagating crack front was also derived. It was found that, at the crack front, material force is reduced to the J-integral in a very special and restrictive case [Citation27].
4. Constitutive equations
The fundamental laws of micromorphic continua consist of a system of 20 partial differential equations, and one inequality. Given the external loads and h, there are 67 unknowns
and θ. Clearly then, the system is highly indeterminate. Forty-seven independent additional equations are needed for the determination of motions and temperatures of a micromorphic body. This is also clear from the fact that the balance equations are valid for all micromorphic bodies irrespective of their physical constitutions. For bodies of different constitutions, the response of the body to external stimuli is very different. Thus, we must bring the constitutional nature of bodies into our formulation. For a constitutive theory to represent a material adequately, certain physical and mathematical requirements must be satisfied. The axiomatic approach to the formulation of constitutive equations was evolved in the course of the development of continuum mechanics. Axioms of causality, determinism, equipresence, objectivity, material invariance, neighborhood, memory, and admissibility are accepted as guiding principles not only in classical continuum mechanics [Citation28] but also in micromorphic theory [Citation12].
For the sake of discussion, in the following we are going to show the formulation of constitutive equations for micromorphic thermoviscoelastic solid and micromorphic fluid.
4.1. Micromorphic thermoviscoelastic solid
For a micromorphic thermoviscoelastic solid, to begin with, the independent and dependent constitutive variables are set to be
It is noticed that the constitutive equation, EquationEquation (44)(44), automatically satisfies the axiom of objectivity because it is expressed in Lagrangian forms. Substituting EquationEquation (44)
(44) into the Clausius–Duhem (CD) inequality, EquationEquation (26)
(26)**, this leads to
Since the CD inequality must be satisfied for all independent thermomechanical processes and EquationEquation (45) is linear in
, it leads to
EquationEquation (46)(46) implies
One may decompose into elastic parts and dissipative parts as
Then the CD inequality becomes
EquationEquation (49)(49)** and the CD inequality, EquationEquation (50)
(50)**, are the constitutive equations for p-th order micromorphic thermoviscoelastic solid. In the limiting case, i.e. p = 1, one has the following constitutive equations for simple micromorphic thermoviscoelastic solid:
Clearly it can be seen that {} and {
} are the thermodynamic forces and the corresponding thermodynamic fluxes. If p is zero and
is eliminated from the list of independent constitutive variables, then we end up with the following constitutive equations for a heat-conducting micromorphic elastic solid:
4.2. Micromorphic fluid
4.2.1. Definition 3
A body is called a micromorphic fluid if every configuration of the body leaving the density and microinertia unchanged can be taken as the reference configuration.
If every configuration is to be taken as a reference configuration, then one can write X = x and Ξ = ξ with ρ and i unchanged. With this, we have
Therefore for micromorphic fluid, to begin with, the independent and dependent constitutive variables can be expressed as
It should be noticed that, in the case of a solid, the dependency on the Lagrangian coordinate X indicating material inhomogeneity and, in the case of fluid, X → x according to the definition, but the Eulerian coordinate x is not objective – that is why x is not in the list of independent variables. It is also noticed that the p-th order strain rates become p-th order Eringen tensors when the micromorphic solid is “melted” into micromorphic fluid.
The constitutive equation, EquationEquation (44)(44), because it is expressed in Lagrangian forms, automatically satisfies the axiom of objectivity. Now the constitutive equation, EquationEquation (54)
(54), is expressed in Eulerian forms, axiom of objectivity requires that [Citation13]
For example,
Substituting EquationEquation (54)(54) into the Clausius–Duhem (CD) inequality, EquationEquation (26)
(26)*, leads to
Since the CD inequality must be satisfied for all independent thermomechanical processes and EquationEquation (58)(58) is linear in
, it results in
If p = 1, the simple micromorphic fluid has the following constitutive equations:
These constitutive equations for a simple micromorphic fluid, EquationEquations (63)(63)
Equation
(64)–
Equation(65)
(65), were obtained by Eringen [Citation13].
5. Micromorphic plasticity
The formulation of constitutive theory for plasticity is unique in the sense that one needs to add a set of internal variables to the list of dependent constitutive variables and, of course, one needs to supply a set of governing equations for the newly added internal variables.
For micromorphic thermo-visco-elastic-plastic (TVEP) continuum, a set of internal variables is introduced as [Citation19]
To separate the material behavior into two distinct parts: thermo-visco-elastic (TVE) part and thermo-visco-elastic-plastic (TVEP) part, a scalar-valued yield function is introduced as
For a set of fixed values of W, a hyper surface, named yield surface, is determined in the eighty-nine-dimensional space of U and V by
We also define the loading rate λ as the scalar product between the outward normal to the yield surface and the tangent vector to the trajectory in the {U, V} space, i.e.
Three distinct cases, unloading, neutral loading, and loading, can be defined by (a) f < 0, (b) f = λ = 0, and (c) f = 0 and λ > 0, respectively. The internal variables of plasticity, W, will remain unchanged in the cases of unloading and neutral loading.
Following the axiom of equipresence, the constitutive relations of a micromorphic TVEP material are initiated as
The Kuhn–Tucker conditions for general plasticity of a micromorphic continuum can now be expressed as
The CD inequality, EquationEquation (26)(26)**, now reads
Since the inequality, EquationEquation (78)(78), must be satisfied for any value of the loading rate λ, it implies
Also, it is emphasized that a TVEP state should lead to another TVEP state; in other words, the consistency condition of plasticity requires that f = 0 and λ > 0 lead to another state with f = 0, which implies
This gives another constitutive constraint to the plasticity
6. Micromorphic electromagnetic–thermoelastic solid
Eringen [Citation15] and Lee et al. [Citation29] formulated the constitutive equations for a micromorphic electromagnetic–thermoelastic solid. Lee and Chen [Citation30] studied the coupling problem of wave propagation in a micromorphic electromagnetic elastic solid. Here, in this work, we briefly go through the formulation of the constitutive theory of micromorphic electromagnetic–thermoelastic solid.
The balance laws in electromagnatics (EM) are the well-known Maxwell's equations, which can be expressed as
The polarization vector, Pk , and the magnetization vector, Mk , are defined as
These EM vectors mentioned above are all referred to a fixed laboratory frame RC . The Galilean transformations of inertial frames form a group that consists of time-independent spatial rotations and pure Galilean transforms, i.e.
The requirement of the form-invariance of the Maxwell's equations under the Galilean transformations leads to the following transformations [Citation31]:
The balance laws of linear momentum, moment of momentum, and energy of micromorphic continuum with EM interactions can be expressed as
The second law of thermodynamics, also referred to as the Clausius–Duhem inequality, is expressed as
The Helmholtz's free energy density ψ with EM interaction is now introduced as
Then the Clausius–Duhem inequality can be becomes
Since the generalized Lagrangian strains [cf. EquationEquations (16)(16)] and their material time derivatives of any order [cf. EquationEquations (19)
(19)] are objective, and hence they are suitable for being employed as independent constitutive variables in the development of a constitutive theory. In the same spirit, define the Lagrangian forms of the electric field vector and the magnetic flux vector as
We also define
Now, the Clausius–Duhem inequality, EquationEquation (107)(107), can be rewritten as
Let the independent and dependent constitutive variables be
Since the inequality is linear in ,
,
,
,
,
,
, it holds if, and only if,
EquationEquations (117)(117)
Equation
(118)–
Equation(119)
(119) are the generalized Gibbs equations for a micromorphic electromagnetic–thermoelastic solid. These constitutive relations are further subjected to the axioms of material invariance and time reversal. It may be stated that the constitutive response functionals must be form-invariant with respect to a group of transformations of the material frame of reference
and microscopic time reversal
representing the material symmetry conditions, and these transformations must leave the density and charge at (X, t) unchanged [Citation31]. The magnetic symmetry properties of solids cannot be discussed rationally by means of three-dimensional point groups only since magnetism is the result of the spin magnetic moment of electrons, which changes sign upon the time reversal. In other words, diamagnetic and paramagnetic crystals do not exhibit any orderly distribution of their spin magnetic moments, and are therefore ‘time symmetric’. The crystallographic point group is enough for the discussion of their material symmetries; on the other hand, for ferromagnetic, ferromagnetic and anti-ferromagnetic materials, which are characterized by an orderly distribution of spin magnetic moment, an additional symmetry operator is needed to take care of the time reversal. For a complete account of this subject, interested readers are referred to Shubnikov and Belov [Citation34] and Kiral and Eringen [Citation35].
7. Discussion
There is one concern in applying micromorphic theory to practical problem: it involves many material constants. How do we measure those material constants? Chen and Lee [Citation36] proposed an algorithm to determine the material constants for micromorphic elastic solids, such as single crystals, through the phonon dispersion relations obtained by atomistic calculations or experimental measurements. Later Zeng et al. [Citation37] extended the algorithm to determine material constants in non-local micromorphic theory. The phonon dispersion relations from atomistic calculations and/or experimental measurements provide the means to determine the material constants in micromorphic theory, which offers a great promise in the applications of semiconductor physics, MEMS, and other microsystems. Beyond what has been mentioned above, micromorphic theory, including its special cases microstretch theory and micropolar theory, has been applied to liquid crystals, theory of turbulence, blood, anisotropic fluids, suspensions, etc.
However, in micromorphic theory, there is an assumption that the micromotion is an affine motion [cf. EquationEquation (2)(2)]. This assumption makes it difficult to describe the atomic motion of complex crystalline materials. To enlarge the domain of applicability of micromorphic theory, Chen and Lee [Citation38] abandoned the assumption of affine motion, let the micromotions have the generality as
, and formulated a generalized continuum field theory named atomistic field theory (AFT). Since AFT is a concurrent atomistic/continuum field theory in its own right, it has the advantage in dealing with some critical phenomena, especially those related to defect nucleation and evolution in solids such as dynamic crack propagation [Citation39].
References
- Tadmor , E.B. , Ortiz , M. and Phillips , R. 1996 . Quasicontinuum analysis of defects in solids . Phil. Mag. , 73 : 1529 – 1563 .
- Knap , J. and Ortiz , M. 2001 . An analysis of the quasicontinuum method . J. Mech. Phys. Solid , 49 : 1899 – 1923 .
- Liu , W.K. , Karpov , E.G. , Zhang , S. and Park , H.S. 2004 . An introduction to computational nanomechanics and materials . Comput. Meth. Appl. Mech. Eng. , 193 : 1529 – 1578 .
- Haile , J.M. 1997 . Molecular Dynamics Simulation: Elementary Methods , New York : Wiley-Interscience .
- Cosserat , E. and Cosserat , F. 1909 . Theories des Corps Deformable , Paris : Hermann .
- Toupin , R.A. 1962 . Elastic materials with couple-stresses . Arch. Ration. Mech. Anal. , 11 : 385 – 414 .
- Eringen , A.C. 1964 . Simple micro-fluids . Int. J. Eng. Sci. , 2 : 205 – 217 .
- Mindlin , R.D. 1964 . Microstructure in linear elasticity . Arch. Ration. Mech. Anal. , 16 : 51 – 78 .
- Eringen , A.C. and Suhubi , E.S. 1965 1964 . Nonlinear theory of simple micro-elastic solids I . Int. J. Eng. Sci. , 2 : 189 – 203 .
- Eringen , A.C. Theory of micropolar continua . Proceedings of the Ninth Midwestern Mechanics Conference, Wisconsin, 16–18 August . New York : Wiley .
- Eringen , A.C. and Suhubi , E.S. 1964 . Nonlinear theory of simple micro-elastic solids II . Int. J. Eng. Sci. , 2 : 389 – 404 .
- Eringen , A.C. 1999 . Microcontinuum Field Theories-I: Foundations and Solids , New York : Springer .
- Eringen , A.C. 2001 . Microcontinuum Field Theories-II: Fluent Media , New York : Springer .
- Eringen , A.C. 2002 . Nonlocal Continuum Field Theories , New York : Springer .
- Eringen , A.C. 2003 . Continuum theory of micromorphic electromagnetic thermoelastic solids . Int. J. Eng. Sci. , 41 : 653 – 665 .
- Eringen , A.C. 2005 . A continuum theory of dense suspensions . Z. Angew. Math. Phys. , 56 : 529 – 547 .
- Eringen , A.C. 2006 . Nonlocal continuum mechanics based on distributions . Int. J. Eng. Sci. , 44 : 141 – 147 .
- Eringen , A.C. 2010 . Micromorphic theory of turbulence . Z. Angew. Math. Phys. , 61 : 119 – 132 .
- Lee , J.D. and Chen , Y. 2003 . Constitutive relations of micromorphic thermoplasticity . Int. J. Eng. Sci. , 41 : 387 – 399 .
- Eshelby , J.D. 1951 . The force on an elastic singularity . Phil. Trans. Roy. Soc. Lond. A , 244 : 87 – 112 .
- Maugin , G.A. 1993 . Material Inhomogeneities in Elasticity , London : Chapman & Hall .
- Maugin , G.A. 1995 . Material Forces: Concepts and Applications . Appl. Mech. Rev. , 48 : 213 – 245 .
- Peach , M.O. and Koehler , J.S. 1950 . Forces exerted on dislocations and the stress field produced by them . Phys. Rev. II , 80 : 436 – 439 .
- Rice , J.R. 1968 . Path-independent integral and the approximate analysis of strain concentrations by notches and cracks . Trans. ASME J. Appl. Mech. , 33 : 79 – 385 .
- Casal , P. 1978 . Interpretation of the Rice integral in continuum mechanics . Lett. Appl. Eng. Sci. , 16 : 335 – 347 .
- Maugin , G.A. 1992 . Application of an energy–momentum tensor in nonlinear elastodynamics: Pseudomomentum and Eshelby Stress in Solitonic Elastic Systems . J. Mech. Phys. Solid , 40 : 1543 – 1558 .
- Chen , Y. and Lee , J.D. 2005 . Material force for dynamic crack propagation in multiphase micromorphic materials . Theor. Appl. Fract. Mech. , 43 : 335 – 341 .
- Eringen , A.C. 1980 . Mechanics of Continua , Melbourne, FL : Robert E. Krieger Publishing Co .
- Lee , J.D. , Chen , Y. and Eskandarian , A. 2004 . A micromorphic electromagnetic theory . Int. J. Solid. Struct. , 41 : 2099 – 2110 .
- Lee , J.D. and Chen , Y. 2004 . Electromagnetic wave propagation in micromorphic elastic solids . Int. J. Eng. Sci. , 42 : 841 – 848 .
- Eringen , A.C. and Maugin , G.A. 1990 . Electrodynamics of Continua I: Foundations and Solid Media , New York : Springer-Verlag .
- Jackson , J.D. 1975 . Classical Electrodynamics , New York : John Wiley & Sons .
- De Groot , S.R. and Suttorp , L.G. 1972 . Foundations of Electrodynamics , Amsterdam : North-Hollard .
- Shubnikov , A.V. and Belov , N.V. 1964 . Coloured Symmetry , Oxford : Pergamon Press .
- Kiral , E. and Eringen , A.C. 1990 . Constitutive Equations of Nonlinear Electromagnetic–Elastic Crystals , New York : Springer-Verlag .
- Chen , Y. and Lee , J.D. 2003 . Determining material constants in micromorphic theory through phonon dispersion relations . Int. J. Eng. Sci. , 41 : 871 – 886 .
- Zeng , X. , Chen , Y. and Lee , J.D. 2006 . Determining material constants in nonlocal micromorphic theory through phonon dispersion relations . Int. J. Eng. Sci. , 44 : 1334 – 1345 .
- Chen , Y. and Lee , J.D. 2005 . Atomistic formulation of a multiscale theory for nano/micro physics . Phil. Mag. , 85 : 4095 – 4126 .
- Lee , J.D. , Wang , X. and Chen , Y. 2009 . Multiscale material modeling and its application to a dynamic crack propagation problem . Theor. Appl. Fract. Mech. , 51 : 33 – 40 .