# Nonlinear Classical Elasticity Model for Materials with Fluid and Matrix Phases

Authors: Rafael Muñoz; Juan Melchor
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Materials with fluid and matrix phases present different acoustic responses in each phase. While longitudinal waves propagate in both phases, shear waves do it only through the solid matrix. Longitudinal waves are mainly described by volumetric propagation and shear waves by deviatoric processes. In the case of nonlinear propagation cross effects occur between both components. This paper presents a new classical nonlinear model proposing a constitutive equation that separates volumetric and deviatoric effects. Four nonlinear constants of third order are defined. The formulation is compared to constitutive equations with Landau constants for weakly elasticity and both types of nonlinear constants related. Some reinterpretation of the Landau’s constants arises in terms of parallel or cross nonlinear effects between volumetric and deviatoric components.


---

## Body

## 1. Introduction

The study of the coefficients of acoustic nonlinearity has recently been recovered again [1–6]. The classical nonlinear behavior of a material can be characterized using different parameters: first-second nonlinear parameters, Landau and Murnaghan constants, and stiffness matrices of second and third order [7–9]. An example of this interest could be the advances in early damage detection [10–13] of materials with complex structures in Nondestricutive Evaluation field, allowing new paradigms that explore the connection between the micromechanical scales.On the mathematical description of nonlinear effects, it is important to separate the geometrical nonlinearity from the material-dependent constitutive nonlinearity. The nonlinear analysis is based on Taylor series expansions of the strain energy, where the Landau [14] and Murnaghan [15] coefficients surge and only account for the material-dependent constitutive nonlinearity. The dynamic problem and its wave equation deduction bring several calculations that introduce geometrical nonlinearities. The first one arises when the stress tensor is calculated from the strain energy expansion, where several operations involving the deformation gradient tensor generate new nonlinear terms (see (13), terms that are independent of Landau constants). The second, the kinematic relations add the well-known cross derivative term of the displacements, which is independent on the material nonlinear behavior.The understanding of the nonlinear constants meaning (Landau, Murnaghan,β, δ, or stiffness coefficients) turns tough, and their interpretation does not appear as straight forward as it could seem. Probably, the most intuitive interpretation comes from the first nonlinear parameter β and the Finite Amplitude technique [16], where the parameter is a measurement of the growth on the second harmonic amplitude with distance. However, this parameter must be further carefully interpreted, since different definitions of them can be found in literature while nomenclature is maintained the same [10, 17]. As some examples in this vain, first and second nonlinear parameters (β and δ), defined directly as wave equation coefficients [18], usually accounts for both geometrical and constitutive nonlinearity. There are simplifications on some developments that define only constitutive versions of the coefficients β and δ, directly applying a Taylor series expansion on the stress as a function of strain that are used for unidimensional plane wave propagation cases and neglecting geometrical nonlinearity [19, 20]. Anyway, all these theoretical proposals have definitely contributed to understand the main nonlinear experimental phenomena. But new approaches can better interpret different situations and open new visions.This paper proposes deepening on the study of classical nonlinear acoustics, proposing a new mathematical formulation with the aim of presenting a different classical nonlinear equation. It splits the general first-order nonlinearity (usually measured byβ or by Landau’s constants) into four specific nonlinear phenomena (measured with four new nonlinear constants), due to the interaction of the deviatoric and volumetric components of the deformation and stress tensors. This paradigm turns specially useful to explain shear wave nonlinear propagation in materials constituted by two phases, as solid fibers embedded in a quasi-fluid matrix. This is the case of soft tissues, where the propagation of shear waves in the liquid phase can be neglected against that occurring in the solid structure.Nonlinear first- and second-order constants are usually referred to the nonlinear order of the terms in the wave equation, or in other cases at the constitutive equation, while other nonlinear constants are referred to the order of the Taylor series expansion of the Energy, like Landau constants as third order. As the new four nonlinear constants are compared to Landau constants in this paper, they will be referred as third-order nonlinear constants.The concept would be scalable to the third and higher harmonics [21–23]. Non-classical nonlinear effects are not considered at the moment, in order to ease the analysis, and viscosity, just partially, but they could be considered in future works, since they could be noteworthy in some cases. The connection of these new four nonlinear constants with the Landau third-order elasticity constants is also explored.


## 2. Methods

The general method of this paper is the deduction of the nonlinear wave equation, where the new proposed constitutive equation is the starting point.(1)σij=σijv,Dijwhere v and Dij are the volumetric and deviatoric parts of the deformation tensor εij. This equation defines the new four β nonlinear constants. The constitutive equation is transformed into σij=σij(εij) and compared with the similar expression given by Landau formulation for weakly elasticity to relate the four β constants to Landau nonlinear constants.Finally, the nonlinear equation of motion as a function of the new paremeters is deducted,(2)ρu¨i=σij,j


### 2.1. Theoretical Preliminars

The dynamic elastic problem comprises the momentum balance equation, compatibility equations, and the kinematic relationships. For wave propagation (body forces neglected), the first and third equations would be(3)ρu¨i+γρu˙i=σij,j(4)εij=12ui,j+uj,i+uk,iuk,jwhere Table 1 shows the symbols of the magnitudes.Table 1
List of variables.
Quantity Symbol Units Displacement u i m Stress σ i j Pa Strain ε i j - Space x i m Time t s Lamé constants λ, μ Pa Density ρ kg/m3 Kinematic viscosity γ=ηρ s - 1The compatibility equations can be considered separating them into linear, viscous, and nonlinear stress components,σijL, σijV, and σijNL.(5)σij=σijL+σijV+σijNLThe nonlinear term can be deducted, following a similar concept of series expansion put forth by Landau [24, 25], where the volumetric and deviatoric decomposition of the stress tensor was considered, and only the volumetric part was detailed in terms of the nonlinear parameter β, in the following expression:(6)-p=-3Kv+9Kβv2-3ηvv˙. where p and v are the volumetric components of the stress and deformation tensors, respectively, so that(7)σij=-pδij︸volumetric+τij︸deviatoric,p=-13σkk(8)εij=-vδij+Dij,v=-13εkkwith τij and Dij being the deviatoric parts of the stress and strain tensors.


### 2.2. Proposal of a New Constitutive Equation

However, a more generalized constitutive relationship can be hypothesized using a combination of four of nonlinear parameters of third-orderβ that may explain a different scenario of calculations. These combinations could be expanded as exploring the whole set of combinations by quadratic terms as follows:(9)σij=-3Kvδij︸pressure+2μDij︸shear︸σijLLinear-3ηvv˙δij︸pressure+2ηD˙ij︸shear︸σijVViscous+9Kβvpv2δij+9KβdpDkpDpkδij︸pressure+4μβdsDikDkj+4μβcsvDij︸shear︸σijNLNonlinearThe constants K and μ accompanying nonlinear parameters have been chosen in accordance with (6), as the quadratic power expansion. Four nonlinear parameters of third order (third, in terms of Energy expansion) have been defined and their terms can be interpreted, as follows:(i)
Term withβvp: nonlinear volumetric stress generated from the behavior of the volumetric component of the deformation(ii)
Term withβdp: nonlinear volumetric stress generated from the behavior of the deviatoric component of the deformation(iii)
Term withβds: nonlinear stress generated from the behavior of the deviatoric component of the deformation(iv)
Term withβcs: nonlinear deviatoric stress generated from the interaction between the volumetric and deviatoric components of the deformationTo complete all possible combinations of quadratic terms, there are two additional combinations which have been removed, since the deviatoric traceDkk is always null.(i)
A term withvDkkδij: it would be a nonlinear volumetric stress generated from the volumetric part of the deformation(ii)
A term withDkk2δij: with similar interpretation than the previous oneThis expression of the stress presents a similar structure to the expression of the stress with Third-Order Elastic Constants in the Landau form(A,B,C), shown below in (13), so it will be considered later to obtain the differential equation in a parallel deduction.Notwithstanding, there is another interesting expression fully separating the nonlinear volumetric and deviatoric components of the stress. The stress term withβds is dependent on the square DikDkj of the deviatoric component of the deformation Dij and can be split into a volumetric and deviatoric part(10)DikDkj=13DpkDkpδij+dij(11)dij=deviatoricDikDkjThus,(12)σij=-3Kvδij︸pressure+2μDij︸shear︸σijLLinear-3ηvv˙δij︸pressure+2ηD˙ij︸shear︸σijVViscous+9Kβvpv2δij+9Kβdp+4μβdsDkpDpkδij︸pressure+4μβdsdij+4μβcsvDij︸shear︸σijNLNonlinear


### 2.3. Relationship to Formulations with Landau Constants

As aforementioned, a quite similar expression of the Cauchy stress tensor can be found as a function of the Landau constants and the deformation tensor, for weakly elasticity:(13)σij=λδijεkk+2μεij︸σijL+2ηε˙ij-23ηδijε˙kk+ηvδijε˙kk︸σijV+A+4μεikεkj+Bεkpεpkδij+2B+K-53μεkkεij+C-K+23μεkk2δij︸σijNLIn order to obtain the relationship between the Landau constants and the new four nonlinear parameters β, (13) is separated into volumetric and deviatoric components:(14)εij=-vδij+Dij(15)εkk=-vδkk+Dkk=-3v(16)εik=-vδik+Dik(17)εkj=-vδkj+Dkj(18)εikεkj=v2δikδkj-vDikδkj-vDkjδik+DikDkj(19)εkp=-vδkp+Dkp(20)εpk=-vδpk+Dpk(21)εpkεkp=v2δpkδkp-vDpkδkp-vDkpδpk+DpkDkp(22)εkk2=9v2 The nonlinear terms of (13) result in the following expressions:(23)A+4μεikεkj=A+4μv2δikδkj-vDikδkj-vδikDkj+DikDkj(24)Bεkpεpkδij=B3v2+DkpDpkδij(25)2B+K-53μεkkεij=2B+K-53μ3v2δij-3vDij(26)C-K+23εkk2δij=C-K+23μ9v2δij The above analysis is also valid by combining the nonlinear part of the stress with v and Dij in the constitutive equation,(27)σijNL=A+9B+9C-3Kv2δij+BDkpDpkδij+A+4μDikDkj+-2A+6B-6K+2μvDijComparing this expression with (9),(28)βvp=A+9B+9C-3K9K(29)βdp=B9K(30)βds=A+4μ4μ(31)βcs=-2A-6B-6K+2μ4μ(32)A=4μβds-1(33)B=9Kβdp(34)C=Kβvp-49μβds-9Kβdp+13KFor the viscous components, it is shown that the establish definition in (6) matches Landau’s one:(35)2ηε˙ij=2η-v˙δij+D˙ij(36)-23ηε˙kkδij=2ηv˙δij(37)ηvε˙kkδij=-3ηvv˙δij(38)σijV=-3ηvv˙δij+2ηD˙ij


### 2.4. Differential Wave Equation

Adopting the acoustic nonlinear constitutive equation presented in (9), in terms of deviatoric a volumetric parts of the deformation, it is possible to obtain the three dimensional nonlinear equation of motion in terms of the new four parameters β. Previously, the kinematic relations must be introduced in the constitutive equation.(39)εij=12∂ui∂xj+∂uj∂xi+∂uk∂ui∂uk∂uj The wave equation is,(40)
ρ
∂
2
u
i
∂
t
2
=
K
∂
2
u
k
∂
x
k
∂
x
i
+
∂
u
l
∂
x
k
∂
2
u
l
∂
x
k
∂
x
i
+
μ
∂
2
u
i
∂
x
j
2
+
∂
2
u
j
∂
x
i
∂
x
j
+
∂
2
u
k
∂
x
j
2
∂
u
k
∂
x
i
+
∂
2
u
k
∂
x
i
∂
x
j
∂
u
k
∂
x
j
-
2
3
∂
2
u
k
∂
x
k
∂
x
i
+
∂
u
l
∂
x
k
∂
2
u
l
∂
x
k
∂
x
i
+
2
K
β
v
p
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
9
K
β
d
p
1
2
∂
2
u
l
∂
x
k
∂
x
i
∂
u
l
∂
x
k
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
l
∂
x
k
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
l
∂
x
k
-
2
3
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
4
μ
β
d
s
1
4
∂
2
u
i
∂
x
k
∂
x
j
∂
u
k
∂
x
j
+
∂
2
u
k
∂
x
i
∂
x
j
∂
u
k
∂
x
j
+
∂
2
u
i
∂
x
k
∂
x
j
∂
u
j
∂
x
k
+
∂
2
u
k
∂
x
i
∂
x
j
∂
u
j
∂
x
k
+
4
μ
β
d
s
1
4
∂
2
u
k
∂
x
j
∂
x
j
∂
u
i
∂
x
k
+
∂
2
u
j
∂
x
k
∂
x
j
∂
u
i
∂
x
k
+
∂
2
u
k
∂
x
j
∂
x
j
∂
u
k
∂
x
i
+
∂
2
u
j
∂
x
k
∂
x
j
∂
u
k
∂
x
i
+
4
μ
β
d
s
-
1
3
∂
2
u
k
∂
x
k
∂
x
j
∂
u
i
∂
x
j
+
∂
2
u
k
∂
x
k
∂
x
j
∂
u
j
∂
x
i
+
∂
2
u
i
∂
x
j
2
∂
u
k
∂
x
k
+
∂
2
u
j
∂
x
i
∂
x
j
∂
u
k
∂
x
k
+
2
9
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
4
μ
β
d
s
-
1
6
∂
2
u
l
∂
x
k
∂
x
i
∂
u
l
∂
x
k
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
l
∂
x
k
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
l
∂
x
k
+
2
9
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
4
μ
β
c
s
-
1
6
∂
2
u
k
∂
x
k
∂
x
j
∂
u
i
∂
x
j
+
∂
2
u
k
∂
x
k
∂
x
j
∂
u
j
∂
x
i
+
∂
2
u
i
∂
x
j
2
∂
u
k
∂
x
k
+
∂
2
u
j
∂
x
i
∂
x
j
∂
u
k
∂
x
k
+
2
9
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
lwhere, as it was defined, K is the Bulk modulus, μ is the shear modulus, ρ is the density, and βvp,βdp,βds, and βcs are the four nonlinear parameters of third order already explained in the constitutive expression.


## 2.1. Theoretical Preliminars

The dynamic elastic problem comprises the momentum balance equation, compatibility equations, and the kinematic relationships. For wave propagation (body forces neglected), the first and third equations would be(3)ρu¨i+γρu˙i=σij,j(4)εij=12ui,j+uj,i+uk,iuk,jwhere Table 1 shows the symbols of the magnitudes.Table 1
List of variables.
Quantity Symbol Units Displacement u i m Stress σ i j Pa Strain ε i j - Space x i m Time t s Lamé constants λ, μ Pa Density ρ kg/m3 Kinematic viscosity γ=ηρ s - 1The compatibility equations can be considered separating them into linear, viscous, and nonlinear stress components,σijL, σijV, and σijNL.(5)σij=σijL+σijV+σijNLThe nonlinear term can be deducted, following a similar concept of series expansion put forth by Landau [24, 25], where the volumetric and deviatoric decomposition of the stress tensor was considered, and only the volumetric part was detailed in terms of the nonlinear parameter β, in the following expression:(6)-p=-3Kv+9Kβv2-3ηvv˙. where p and v are the volumetric components of the stress and deformation tensors, respectively, so that(7)σij=-pδij︸volumetric+τij︸deviatoric,p=-13σkk(8)εij=-vδij+Dij,v=-13εkkwith τij and Dij being the deviatoric parts of the stress and strain tensors.


## 2.2. Proposal of a New Constitutive Equation

However, a more generalized constitutive relationship can be hypothesized using a combination of four of nonlinear parameters of third-orderβ that may explain a different scenario of calculations. These combinations could be expanded as exploring the whole set of combinations by quadratic terms as follows:(9)σij=-3Kvδij︸pressure+2μDij︸shear︸σijLLinear-3ηvv˙δij︸pressure+2ηD˙ij︸shear︸σijVViscous+9Kβvpv2δij+9KβdpDkpDpkδij︸pressure+4μβdsDikDkj+4μβcsvDij︸shear︸σijNLNonlinearThe constants K and μ accompanying nonlinear parameters have been chosen in accordance with (6), as the quadratic power expansion. Four nonlinear parameters of third order (third, in terms of Energy expansion) have been defined and their terms can be interpreted, as follows:(i)
Term withβvp: nonlinear volumetric stress generated from the behavior of the volumetric component of the deformation(ii)
Term withβdp: nonlinear volumetric stress generated from the behavior of the deviatoric component of the deformation(iii)
Term withβds: nonlinear stress generated from the behavior of the deviatoric component of the deformation(iv)
Term withβcs: nonlinear deviatoric stress generated from the interaction between the volumetric and deviatoric components of the deformationTo complete all possible combinations of quadratic terms, there are two additional combinations which have been removed, since the deviatoric traceDkk is always null.(i)
A term withvDkkδij: it would be a nonlinear volumetric stress generated from the volumetric part of the deformation(ii)
A term withDkk2δij: with similar interpretation than the previous oneThis expression of the stress presents a similar structure to the expression of the stress with Third-Order Elastic Constants in the Landau form(A,B,C), shown below in (13), so it will be considered later to obtain the differential equation in a parallel deduction.Notwithstanding, there is another interesting expression fully separating the nonlinear volumetric and deviatoric components of the stress. The stress term withβds is dependent on the square DikDkj of the deviatoric component of the deformation Dij and can be split into a volumetric and deviatoric part(10)DikDkj=13DpkDkpδij+dij(11)dij=deviatoricDikDkjThus,(12)σij=-3Kvδij︸pressure+2μDij︸shear︸σijLLinear-3ηvv˙δij︸pressure+2ηD˙ij︸shear︸σijVViscous+9Kβvpv2δij+9Kβdp+4μβdsDkpDpkδij︸pressure+4μβdsdij+4μβcsvDij︸shear︸σijNLNonlinear


## 2.3. Relationship to Formulations with Landau Constants

As aforementioned, a quite similar expression of the Cauchy stress tensor can be found as a function of the Landau constants and the deformation tensor, for weakly elasticity:(13)σij=λδijεkk+2μεij︸σijL+2ηε˙ij-23ηδijε˙kk+ηvδijε˙kk︸σijV+A+4μεikεkj+Bεkpεpkδij+2B+K-53μεkkεij+C-K+23μεkk2δij︸σijNLIn order to obtain the relationship between the Landau constants and the new four nonlinear parameters β, (13) is separated into volumetric and deviatoric components:(14)εij=-vδij+Dij(15)εkk=-vδkk+Dkk=-3v(16)εik=-vδik+Dik(17)εkj=-vδkj+Dkj(18)εikεkj=v2δikδkj-vDikδkj-vDkjδik+DikDkj(19)εkp=-vδkp+Dkp(20)εpk=-vδpk+Dpk(21)εpkεkp=v2δpkδkp-vDpkδkp-vDkpδpk+DpkDkp(22)εkk2=9v2 The nonlinear terms of (13) result in the following expressions:(23)A+4μεikεkj=A+4μv2δikδkj-vDikδkj-vδikDkj+DikDkj(24)Bεkpεpkδij=B3v2+DkpDpkδij(25)2B+K-53μεkkεij=2B+K-53μ3v2δij-3vDij(26)C-K+23εkk2δij=C-K+23μ9v2δij The above analysis is also valid by combining the nonlinear part of the stress with v and Dij in the constitutive equation,(27)σijNL=A+9B+9C-3Kv2δij+BDkpDpkδij+A+4μDikDkj+-2A+6B-6K+2μvDijComparing this expression with (9),(28)βvp=A+9B+9C-3K9K(29)βdp=B9K(30)βds=A+4μ4μ(31)βcs=-2A-6B-6K+2μ4μ(32)A=4μβds-1(33)B=9Kβdp(34)C=Kβvp-49μβds-9Kβdp+13KFor the viscous components, it is shown that the establish definition in (6) matches Landau’s one:(35)2ηε˙ij=2η-v˙δij+D˙ij(36)-23ηε˙kkδij=2ηv˙δij(37)ηvε˙kkδij=-3ηvv˙δij(38)σijV=-3ηvv˙δij+2ηD˙ij


## 2.4. Differential Wave Equation

Adopting the acoustic nonlinear constitutive equation presented in (9), in terms of deviatoric a volumetric parts of the deformation, it is possible to obtain the three dimensional nonlinear equation of motion in terms of the new four parameters β. Previously, the kinematic relations must be introduced in the constitutive equation.(39)εij=12∂ui∂xj+∂uj∂xi+∂uk∂ui∂uk∂uj The wave equation is,(40)
ρ
∂
2
u
i
∂
t
2
=
K
∂
2
u
k
∂
x
k
∂
x
i
+
∂
u
l
∂
x
k
∂
2
u
l
∂
x
k
∂
x
i
+
μ
∂
2
u
i
∂
x
j
2
+
∂
2
u
j
∂
x
i
∂
x
j
+
∂
2
u
k
∂
x
j
2
∂
u
k
∂
x
i
+
∂
2
u
k
∂
x
i
∂
x
j
∂
u
k
∂
x
j
-
2
3
∂
2
u
k
∂
x
k
∂
x
i
+
∂
u
l
∂
x
k
∂
2
u
l
∂
x
k
∂
x
i
+
2
K
β
v
p
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
9
K
β
d
p
1
2
∂
2
u
l
∂
x
k
∂
x
i
∂
u
l
∂
x
k
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
l
∂
x
k
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
l
∂
x
k
-
2
3
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
4
μ
β
d
s
1
4
∂
2
u
i
∂
x
k
∂
x
j
∂
u
k
∂
x
j
+
∂
2
u
k
∂
x
i
∂
x
j
∂
u
k
∂
x
j
+
∂
2
u
i
∂
x
k
∂
x
j
∂
u
j
∂
x
k
+
∂
2
u
k
∂
x
i
∂
x
j
∂
u
j
∂
x
k
+
4
μ
β
d
s
1
4
∂
2
u
k
∂
x
j
∂
x
j
∂
u
i
∂
x
k
+
∂
2
u
j
∂
x
k
∂
x
j
∂
u
i
∂
x
k
+
∂
2
u
k
∂
x
j
∂
x
j
∂
u
k
∂
x
i
+
∂
2
u
j
∂
x
k
∂
x
j
∂
u
k
∂
x
i
+
4
μ
β
d
s
-
1
3
∂
2
u
k
∂
x
k
∂
x
j
∂
u
i
∂
x
j
+
∂
2
u
k
∂
x
k
∂
x
j
∂
u
j
∂
x
i
+
∂
2
u
i
∂
x
j
2
∂
u
k
∂
x
k
+
∂
2
u
j
∂
x
i
∂
x
j
∂
u
k
∂
x
k
+
2
9
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
4
μ
β
d
s
-
1
6
∂
2
u
l
∂
x
k
∂
x
i
∂
u
l
∂
x
k
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
l
∂
x
k
∂
x
i
∂
u
k
∂
x
l
+
∂
2
u
k
∂
x
l
∂
x
i
∂
u
l
∂
x
k
+
2
9
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
l
+
4
μ
β
c
s
-
1
6
∂
2
u
k
∂
x
k
∂
x
j
∂
u
i
∂
x
j
+
∂
2
u
k
∂
x
k
∂
x
j
∂
u
j
∂
x
i
+
∂
2
u
i
∂
x
j
2
∂
u
k
∂
x
k
+
∂
2
u
j
∂
x
i
∂
x
j
∂
u
k
∂
x
k
+
2
9
∂
2
u
k
∂
x
k
∂
x
i
∂
u
l
∂
x
lwhere, as it was defined, K is the Bulk modulus, μ is the shear modulus, ρ is the density, and βvp,βdp,βds, and βcs are the four nonlinear parameters of third order already explained in the constitutive expression.


## 3. Discussion

Considering the new nonlinear constitutive equation (12), it can be applied to the case of a pure plane horizontal shear wave. Disregarding the viscous term, the deformation tensor εij will have all components null but ε12, yielding Dij=εij, v=-(1/3)εkk=0, and DkpDpk=2ε122,(41)σij=02με1202με1200000︸linear+9Kβdp2+4μβcsε1220009Kβdp2+4μβcsε1220009Kβdp2ε122︸nonlinearIn this case, the linear part of the stress continues with the same shear structure while the nonlinear part is a diagonal matrix meaning that there exists a volumetric component in it that would contribute to posterior volumetric deformations. Only two constants survive reducing the mathematical complexity for shear waves.If it is a compressional wave in thex1 direction, with εij=0∀ij≠11,(42)σij=Kε11δij+2μ23ε11000-13ε11000-13ε11︸linear+Kβvp+6Kβdp+4μβcsε112δij+4μβds-βcs29ε112000-19ε112000-19ε112︸nonlinear, where now the four nonlinear coefficients must be considered, indicating the existence of all kind of cross and parallel effects, although pure volumetric plus deviatoric nonlinear stress components are the result, allowing the separation of both nonlinear components.Landau constants only account for the constitutive nonlinear behavior of the material, expressed in all the nonlinear terms of (13) that include AB and C. As mentioned above, there are several additional nonlinear terms in this expression without these constants, but with K and μ. They represent the geometrical nonlinearity originated in the geometrical transformation from the Piola tensor Pij as a function of the strain energy U to get the Cauchy stress tensor.(43)Pij=∂U∂Fij=Fij∂U∂εij(44)σij=1det⁡FPijFji(45)σij=1detFFij∂U∂εijFjiIt must be noted that the four new β constants would include constitutive and geometrical effects, with the exception of βdp that is only of constitutive nature (see (29) and (28), (30), and (31)).Considering the relations among the four newβ parameters and Landau constants in (32), (33), (34), (28), (29), (30), and (31), some conclusions can be extracted.(i)
As the Landau constantB is only dependent on βdp (see (33)), B is a measure of the intensity on the cross effect from the deviatoric deformation into a nonlinear volumetric stress.(ii)
As the Landau constantA is only dependent on βds (see (32)), A is a measure of the intensity on the nonlinear effect from the deviatoric deformation into the stress, which can be split into both volumetric and deviatoric parts (see (12)).(iii)
As the Landau constantC is only dependent on βvp (see (34)), C is the only one that measures the intensity on the nonlinear effect from the volumetric deformation into the volumetric component of the stress.(iv)
As expected,βcs is linearly dependent on βds and βdp:(46)βcs=-2βds-27K2μβdp-12-3K2μ


## 4. Conclusions

A new approach is proposed in the field of nonlinear acoustics introducing the volumetric and deviatoric separation in deformation and stress. The relevance of this expression is directly linked to the possibility of separating longitudinal and shear waves and selecting the responses of interest. Additional interpretations on the Landau constants are deducted in terms of the cross and parallel effects between volumetric and deviatoric components.



---

*Source: 5049104-2018-11-27.xml*