CC BY
59
mat. nauk [Issues of numerical implementation of the method consistent perturbation parameters when calculating shell designs : Dr. phys. and math. sci. diss.]. Saratov, 2005. (in Russian).
5. Kuznecov V. N. Metod posledovatel'nogo vozmush-cheniia parametrov v prilozhenii k raschetu dina-micheskoi ustoichivosti tonkostennykh obolochechnykh konstruktsii : Diss. kand. tehn. nauk [The method of successive perturbations of the parameters in the application to the calculation of dynamic stability of thin-walled shell structures : Dr. techn. sci. diss.]. Saratov, 2000 (in Russian).
6. Bessonov L. V. Chislennaia realizatsiia algoritma
spektral'nogo kriteriia lokal'noi poteri ustoichivosti obolochechnoi konstruktsi [Numerical realization of algorithm of spectral criterion of local buckling of shell structures]. Issledovaniia po algebre, teorii chisel, funktsional'nomu analizu i smezhnym voprosam : mezhvuz. sb. nauch. tr. [Research in algebra, number theory, functional analysis and related issues : inter-university collection of scientific papers], Saratov, Saratov Univ Press, 2012, iss. 7, pp. 3-9. 7. Pisarenko G. S., Jakovlev A. P., Matveev V. V. Spravochnik po soprotivleniiu materialov [The reference resistance of materials]. Kiev, Naukova dumka, 1988, 734 p.
ON WEAK DISCONTINUITIES AND JUMP EQUATIONS ON WAVE SURFACES IN MICROPOLAR THERMOELASTIC CONTINUA
V. A. Kovalev1, E. V. Murashkin2, Yu. N. Radayev3
1 Moscow City Government University of Management Moscow, 28, Sretenka str., 107045, Moscow, Russia, vlad_koval@mail.ru 2National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), 31, Kashirskoe shosse, 115409, Moscow, Russia, murashkin@ipmnet.ru, evmurashkin@gmail.com
3Institute for Problems in Mechanics of RAS, 101-1, Vernadskogo ave., 119526, Moscow, Russia, murashkin@ipmnet.ru, evmurashkin@gmail.com, radayev@ipmnet.ru, y.radayev@gmail.com
The present study is devoted to problem of propagating surfaces of weak and strong discontinuities of translational displacements, microrotations and temperature in micropolar (MP) thermoelastic (TE) continua. Problems of propagation of weak discontinuities in type-I MPTE continua are discussed. Geometrical and kinematical compatibility conditions due to Hadamard and Thomas are used to study possible wave surfaces of weak discontinuities. Weak discontinuities are discriminated according to spatial orientations of the discontinuities polarization vectors (DPVs). It is shown that the surfaces of weak discontinuities can propagate exist without weak discontinuities of the temperature field. Second part of the paper is concerned the discussions of the propagating surfaces of strong discontinuities of field variables in type-II MPTE continua. Constitutive relations for hyperbolic thermoelastic type-II micropolar continuum is derived by the field theory. The special form of the first variation of the action integral is used in order to obtained 4-covariant jump conditions on wave surfaces. Three-dimensional form of the jump conditions on the surface of a strong discontinuity of thermoelastic field are derived from 4-covariant form.
Key words: micropolar thermoelasticity, type-I continuum, type-II continuum, weak discontinuity, strong discontinuity, shock wave, longitudinal wave, transverse wave, compatibility condition, jump.
1. PRELIMINARY REMARKS
Problems of micropolar continua take its origin from the classical E. and F. Cosserat's paper [1]. Micropolar (MP) continuum theories include not only translational displacements but also additional degrees of freedom. These degrees of freedom are coupled with changes in reper (three directors) associated with microvolume. Such changes may be described by a rotation vector when reper associated with microvolume are rigid rotated. In contrary to conventional elasticity a continuum with microstucture is described by the asymmetric strain and stress tensors known from many previous discussions. Thus the asymmetric elastic theory is characterized by a comparatively large number of constitutive elastic constants need to be determined from the experimental observations. There are several phenomena (for example, the anomalous piezoelectric effect in quartz, the dispersion of elastic waves, as well as a number of other experimentally observed elastic properties of the pure crystals) being beyond the scope of the conventional thermoelasticity (CTE) and piezoelectroelasticity. That is why a development of complex theories seems to be actual.
The type-I micropolar thermoelastic (MPTE-I) continuum may be described in terms of field formalism, for example, from positions of the Green-Naghdi thermoelasticity (GN-theory). Now such mathematical frameworks of the thermoelastic behavior of solids are rapidly refined [2,3]. They are based on different modifications of the classical Fourier law of heat conduction. The refinements aim at derivations of hyperbolic partial differential equations of coupled thermoelasticity. Those are to
simultaneously fulfill the following conditions: 1.) Finiteness of the heat signal propagation velocity, and 2.) The ability of the spatial propagation of the thermoelastic waves without attenuation, and 3.) Existence of distortionless wave forms akin to the classical d'Alembert type waves. GN-thermoelasticity theory itself can be subdivided into the three different types. Type-I theory (GNI/CTE thermoelasticity) is based on the classical Fourier law of heat conduction with an infinite velocity of propagation of an exponentially decaying heat signal. Type-II theory (GNII, hyperbolic thermoelasticity) is characterised by the energy conservation and the finite propagation velocity of thermal waves known as second sound waves.
In-depth study of plane harmonic type-I thermoelastic waves is given in [4]. It is shown that dispersion equation has exactly two complex wavenumbers for a given frequency. Moreover their real and imaginary parts are strictly positive. Plane waves in type-II continua are studied in [5]. Plane thermoelastic wave or second sound wave is characterized by four real wavenumbers. Two of them are strictly positive. In [6] problems concerning plane harmonic wavenumbers of coupled type-III thermoelastic waves are discussed. In all studies [4-6] the linear symmetrical thermoelasticity is employed. They are summarized in the monograph [7].
The present study is devoted to problem of propagating surfaces of weak and strong discontinuities of translational displacements, microrotations and temperature in micropolar (MP) thermoelastic (TE) continua. In first part of the paper is discussed the problems of the propagation of weak discontinuities in type-I MPTE continua. After preliminary remarks (Sec. 1) in Sec. 2 the requisite equations for type-I linear micropolar thermoelastic continuum are considered. Those includes the equation of motion, conservation energy principle and entropy balance equations.
In Sec. 3 compatibility conditions for weak discontinuities propagation surfaces of displacements, microrotations and temperature in the type-I micropolar thermoelastic continuum are given. These conditions are derived from the requisite equations (Sec. 2) and the geometrical and kinematical conditions due to Hadamard and Thomas [8]. The admissible normal velocities of the propagating surfaces of weak discontinuities are found.
In Sec. 4 the distinguished cases of propagating wave surfaces in MPTE-I continuum are separately considered. Propagating surfaces of weak discontinuities are discriminated depending on the spatial orientations of DPVs.
In Sec. 5 the constitutive equations for hyperbolic thermoelastic type-II micropolar continuum is derived by the field theory. The special form of the first variation of the action integral is used in order to obtained 4-covariant jump conditions on wave surfaces.
In Sec. 6 the problems of the propagating surfaces of strong discontinuities of field variables in type-II MPTE continua are considered. Three-dimensional form of the jump conditions on the surface of a strong discontinuity of thermoelastic field are derived from 4-covariant form.
2. REMINDER OF BASIC EQUATIONS OF MPTE-I CONTINUUM
A framework of MPTE-I continuum (see [9] for details) is used throughout the paper. Such continua permit the existence of rotational degrees of freedom (a rotation vector 0) along with translational degrees of freedom (the translational displacement vector u). For MPTE-I continuum, equations of motion of micropolar medium are written in direct tensor representation for the case of the absence of mass forces and mass moments
iV'a=piU' (i)
IV ■ m + e -a = J0,
where p is the mass density tensor, J is a tensor of inertia, e is the three-dimensional Levi-Civita symbol (permutation symbol, antisymmetric symbol, or alternating symbol), V is the three-dimensional Hamiltonian differential operator (the nabla symbol), dot over a symbol denotes partial differentiation with respect to time at fixed spatial coordinates. For p and J tensors the symmetry conditions pij = pji, Jij = Jji are obviously valid. where. Another form of the tensor representation of equations (1) is (see [10])
V ■ a = pU,
(2)
V ■ m + ax I = J0, v '
wherein I is the second-order unit tensor.
Equations (1) and (2) in rectangular co-ordinate system can be deduced in form
(dj Gji = pij Uj
[dj mji + eijk ajk = Jij
where dj is partial derivative on spatial coordinate xj. If the equation of entropy balance are taken in form
' = <-*•(!). (3)
where s denotes the entropy referred to the unit of volume, h is the heat flux vector, 9 is the actual temperature, £ denotes the internal entropy product, then the equations of energy balance may be written as ( )
p = tr(a ■ e) + tr(m ■ r) - s9 - 9 ^h ■ V - ^ , (4)
where p is the free energy per unit of volume.
The asymmetric strain tensor e and bending-torsion tensor r are associated with translational displacements u and microrotations 0 by formulae
e = V ® u - e ■ 0, r = V ® 0. (5)
The equation (5) are represented in a rectangular co-ordinate system as
eji dj Ui Cjik 0k, j dj 0i.
Since the free energy p = p(e, r,9) is the function of the independent variables e, r, 9 then
/" dp dp • dp ■
p = dj j + drjir ji+ ds9 (6)
Therefore the following equations can be obtained from the comparison of equations (4) and (6)
dp dp dp t i u V9
' "ji = dr" jiji
j = deji - "j- = drj-• s = -§• £+h ■ % = o- <7>
The second law of thermodynamics is correct when the £ ^ 0. Whence the Clausius - Duhem inequality can be written in form
V9
-h ■ -§r ^ 0. (8)
The inequality (8) is satisfied by using the Fourier's law of heat conduction, which states the proportionality of the heat flux vector h and negative spatial temperature gradient 9
h = -A ■ V9, (9)
where A is the material thermal conductivity tensor (thermal diffusion tensor). A is a positive definite matrix.
The equation (3) and takin into account the equations (7) and (9) may be transformed as
S9 = -V ■ h = V ■ A ■ V9. (10)
Let us expand the free energy p into the Taylor series in the vicinity of the natural state e = r = 0, 9 = 0, disregarding the terms of higher order than the second one. The following form of the expansion is obtained for isotropic, homogeneous and centrosymmetric bodies
" + n + " - n + A + y + g _ _ + Y - g
2 ejieji + 2 ejieij + 2 ekkenn + 2 1 i + 2
P A* 2 + 2rkkTnn - aekk9 - ?rkk9--— 9 •
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 1 The constitutive equations of MPTE-I continuum can be derived according to equations (7)
a = (д + n)e + (д — n)eT + (Atre — a0)I, m = (7 + e)r + (7 — е)Гт + (£trr — ?0)I, s = atre + ?trr + Л*0,
(11)
wherein A, n,Y, are isothermal constitutive constants of type-I micropolar thermoelastic continuum; a, £ are constitutive constants providing coupling of equations of motion and heat conduction. Constants a, £ depend not only on the mechanical properties of the continuum, but also depend on the thermal properties.
Rewrite (11) as follows
aij = + + - n)ej + (Atrekk - aV)5j, mj = (y + + (y - £)r*j + (^trrfcfc - ,
s = aekk + £ Tkk + A*V.
Note that
ds\ k
r,e =
where k is the heat capacity (per unit volume) at constant (zero) strains. Thus the equation is derived
ds = adtre + £dtrr + KdV.
V
Integrating this equation for entropy in assumption that s = 0 for natural state is obtained
s = atre + £trr + kln ^. (12)
Expand the logarithm in (12) under condition of |V/T0| ^ 1 allows rewrite entropy equation in form
s = atre + £trr + ^V. (13)
To
Consequently the heat conduction equation can be derived by transforming equations (10) and (13) as
V ■ A ■ VV - kv - aV ■ U - £V ■ ( = 0.
After substituting the stress tensors a and m from the formulas (11) in the equations of motion (1) and taking account of (5) the system of coupled partial differential equations of motion and heat conduction for a linear isotropic type-I micropolar thermoelastic continuum in the absence of mass forces, moments, and heat sources can be written as [9, 11]:
(A + ^ - n)VV ■ u + (^ + n)V ■ Vu + 2nV x ( - aVV - pU = 0, (^ + y - £)VV ■ ( + (y + e) V ■ V( - 4n( + 2nV x u - £VV - Jp = 0, (14)
v2 v - fV - a v ■ * - AT V ■ ( = o.
Hereafter V is the temperature increment over the referential temperature, k is the heat capacity (per unit volume) at constant (zero) strains.
As pointed out the nonzero constitutive constants a, £ provide coupling of micropolar thermoelasticity equations. It is usually assumed that £ = 0 (see [9]). This constitutive constant is kept in all further considerations for the completeness analysis .
The scalar equation in the system (14) is called a generalized heat conduction equation conjugate to the equations of motion (the first and the second equations in (14)).
з. PROPAGATING WEAK DISCONTINUITIES OF THERMOMECHANICAL FIELDS IN MPTE-I CONTINUUM
This section is devoted to processes of propagating weak discontinuities of translational displacements
и, microrotations 0 and temperature 9 in the MPTE-I continuum. System of partial differential equations (14) includes partial derivative order not higher than the second. Let a wavefront (wave surface £) of weak discontinuities displacements u, microrotations 0 and temperature 9 be propagated in three-dimensional space with normal velocity G. Unit normal vector to that wave surface is denoted by n. If surface £ is parameterized by
x = P' (y\y2,t),
then components of first gaß and second baß fundamental forms of surface £ are computed by formulas
dp'1 dp'1 0 ^ d 2p' ö-WS = da' X dß' x , baß = 7;-^ «
dya dyß a ß aß dyadyß
Kinematical and geometrical compatibility conditions of the second and the first order due to Hadamard and Thomas valid for an arbitrary physical field p read
9aß = тга ТГв = да' X3ß' X, baß = в П% = да' dß' Xпг, ЭаПг = -QßlbßaOy X
[didjj 0]] = Dn-nj + gaVda' C(n-x3 + Uj x-) - gaVgaTCbaa3^ x-3t' x3,
[di 0]] = (-Dk G + ¿4 C )n- + gaV da' (CG)dp' x-, (15)
[0] = (DG2 - ¿4C)G + CS4G, ¡840] = -CG, [d-vl = Cn-,
where i,j = 1, 2,3; a,P,r,a = 1, 2; the prime in superscripts are denoted by the Greek indices related to surface Gaussian coordinates; D = [didj0]n-nj; S4 is the delta-derivative operator, [ ■ ] denote 3-jumps of the field variables..
If, in particular, the first derivative of the function 0 is continuous across the surface £ the relations (15) can be reduced to
[d-dj 0]] = Dn-Uj, [d- 0] = -DGn-, [0] = DG2.
The obtained expressions are easy to rewrite in direct tensor notations
[V ® V ® u]] = n ® n ® A, [V ® V ® 0] = n ® n ® S,
[V ® u]] = -Gn ® A, [V ® 00]] = -Gn ® S, (16)
[u]] = G2A, [0]= G2S, ( )
[V ® V9] = Bn ® n,
where square brackets denote jump across surface of weak discontinuities. B, A, S are fields defined on this surface, and the equalities B = 0, A = 0, S = 0 cannot be satisfied simultaneously at any point of the surface, if the surface is indeed the surface of weak discontinuities. Equations (14) and (16) give the following relations between the jumps of partial derivatives of the second order across the wave surface:
' (pG2 - (" + n))A - (A + " - n)n(n ■ A) = 0,
(JG2 - (y + g))S - (P + y - g)n(n ■ S) = 0, (17)
„ aG „ qG „ B + — n ■ A + V-n ■ S = 0. A* A*
The polarization vectors of the wave A, S can be splited into a sum of projections onto the tangent plane and on the normal direction to the wave surface:
A = A± t + Ayn, S = S±T + Syn,
A± = A ■ t, Ay = A ■ n, (18)
S± = S ■ t, Sy = S ■ n,
where t is the tangential unit vector and n is the normal unit one of the wave surface respectively. Taking account of equations (18) the system (17) is transformed into
'(pG2 - (" + n))(A±t + Ayn) - (A + " - n)Ayn = 0,
' (JG2 - (y + g))(S±t + Syn) - (P + y - g)Syn = 0, (19)
b+f Ay + AG s, = 0.
After rearrangements in equations (19) in view of (18) the following system is derived:
' (pG2 - (p + n))A±t + (pG2 - (A + 2p)) A||n = 0,
(JG2 - (7 + s))sxt + (JG2 - (y + 2£))S||n = 0, „ aG , ^
+ a: a + A: s«=°-
(20)
The third equation in (20) allows us to couple a weak discontinuity of temperature with a normal projection of weak discontinuities of translational displacements and microrotations by
B = -- iGS,
A,
A,
This new equation is obtained by restrictions ç = 0.
By the aid of orthogonality of the vectors t and n the following equations are obtained
(pG2 - (p + p))Ax = 0, (pG2 - (A + 2p))A|| =0,
(JG2 - (y + e))Sx =0, (JG2 - (y + 20))S|| =0.
(21)
4. A CLASSIFICATION OF WEAK DISCONTINUITIES IN THE MPTE-I CONTINUUM
The 16 cases can be discriminated according to (5). These cases are gathered into the following Table. We proceed by considering the discriminated cases separately.
Cases Intensity vectors projections of weak discontinuities of translational displacements Intensity vectors projections of weak discontinuity of microrotations Intensity of weak discontinuity of temperature
I a = 0 a± =0 s = 0 s± = 0 b = 0
II a = 0 a± =0 s = 0 s± = 0 B = 0
III a =0 a± =0 s =0 s± = 0 b = - ^an Л* 11
IV a =0 a± =0 s = 0 s± = 0 b = 0
V a =0 a± =0 s = 0 s± = 0 b = - ^s„ л* 11
VI a =0 a± =0 s =0 s± = 0 b = -o—an - i^sn л* 11 л* 11
VII a =0 a± =0 s = 0 s± = 0 b = 0
VIII a =0 a± =0 s =0 s± = 0 b = - a—an л* 11
IX a =0 a± =0 s = 0 s± = 0 b = - i—sii л* 11
X a =0 a± =0 s =0 s± = 0 b = - a—an л* 11
XI a =0 a± =0 s = 0 s± = 0 b = 0
XII a =0 a± =0 s = 0 s± = 0 в = -a—an - i—sii л* 11 л* 11
XIII a =0 a± =0 s = 0 s± = 0 b = - i—sii л*
XIV a =0 a± =0 s = 0 s± = 0 в=-ai - i—sn л* л*
XV a =0 a± =0 s =0 s± = 0 b = - a—ay л*
XVI a =0 a± =0 s = 0 s± = 0 в=-ai - i—sn л* i л* i
Evidently in the first case A = 0, S = 0 and returning to (17) the scalar equation in (17) is satisfied identically, so the surface £ is actually not a surface of weak discontinuities.
In the second case, the first equation of system (21) is valid only on a wave surface propagating with normal velocity
G = c||, (22)
where
X + 2^
c»=v
In this case the normal velocity of a weak discontinuity surface propagation G, (22) is the well known velocity of longitudinal elastic wave cy. In addition, if wave surface is bypassed, then the derivative of temperature field is broke according to the third equation in (20):
^ aG 4
B = -17 ■
In the third case a weak discontinuity of translational displacements exists only on the surface propagating with the velocity
G = '"_L
where
G = cf , (23)
4 = ; ^
It can be elucidated that the velocity G (23) exactly coincides with the one of transverse elastic wave if the constitutive micropolar constant n equals to zero c^ = \f^Tp-
The fourth case implies existence of a weak discontinuity of microrotation. Then the third equation of system (20) is satisfied only on the surface of weak discontinuities propagating with normal velocity
G = cW = c|| ,
where
P + 2Y J
is the wave velocity of torsion. As in the previous case, a weak discontinuity of microrotation is associated with a weak discontinuity of a temperature field according to third equation of system (20):
B — - x:Sy-
Thus, the constitutive constants a, ç characterize the intensity of thermal weak discontinuity. In the fifth case, the fourth equation of the system (20) allows to compute the propagation velocity of a weak discontinuities of microrotation
G = c^,
where
c=
Y + £
J
Focus attention that a weak discontinuity of temperature is not associated with the tangential projections of polarization vectors of weak discontinuities A± and S±.
In other cases (from the VI to the XVI) the propagation of wave surfaces of weak discontinuities displacements, microrotations and temperature is impossible if the constitutive characteristic of the type-I micropolar thermoelastic continuum does not satisfied certain limitations as determined by (21).
In a 6-th case the factors at Ay and Sy should be vanished simultaneously, it is impossible if cy — c^. Clearly, then the considered surface will not be a wave one. If the continuum can be described constitutive equation of cy — cj^, then the wave surface is a wave propagating with a velocity
G — cy — cjj,
and weak discontinuity of temperature field is out of condition
B = -fA» - £S»- <24>
The above Eqs (24) implies a possibility of nondissipative propagation of longitudinal waves when B = 0, i.e.
a = _ Si
In the seventh case, the factors at A± and S± can not be equal to zero simultaneously, if c^ = c^. In contrary, if the continuum characterizes matching of the transverse waves velocities c^ = then the velocity of propagation of a wave surface is equal
G = c^ = cf.
Weak discontinuity of temperature field is again identified from the condition (24) and it disappears in the absence of normal polarization vectors projections .
Underline that in the 3-th and 5-th cases the surface of weak discontinuities in MPTE-I continuum can propagated along with the absence of weak discontinuities of temperature field, which means non-dissipative mechanism for the such surfaces propagating.
5. ACTION. COMPATIBILITY CONDITIONS FOR JUMPS
A field theory formalism involves mathematical description of physical fields by integral action functional. A general form of action within a variable domain of 4-spacetime with the elementary volume of d4X = dXMX2dX3dX4 is
I = J L (X 0,da / )d4 X,
where is the physical fields array, l — the Lagrangian density.
The least action principle states that the actual field is realized in the spacetime in a way that the action of (1) is minimum, i. e. for any admissible variations of physical fields and non-variable coordinates X0 are Si = 0. Then the classical Euler-Lagrange equations are held:
dL dL
E (l) da d(da^k) 0
In general, a conservation law has the following form
Op J0 = 0,
where the vector J0 is the vector 4-current. By finite variation Sv = S/e the 4-current can be obtained as
d L
d(d07ky r > v>d(dp^)
^ = дтд^к)^ + - )d(Ilk^) X
if the variational symmetries of the action are known. The finite variation of the action is (see [7])
Svi = J (dp J0)d4X. (25)
Assuming the temperature field is continuous and temperature gradient of the first order can be discontinuous by passing through some bilateral surface £± propagating with the normal velocity G and normal unit 4-vector n0 in 4-spacetime.
We replace the integral over 4-volume in equation (25) on sum of surface integrals. Then only two surface integrals over surface £ are remained when the variations Sv^k and SvX0 are fixed on the outer boundary of the field:
Svi = y (J0n0)+d£ _ J (J0n0)-d£ = J [J0]n0d£. (26)
E+ E" E
Hereafter square brackets denote 4-jumps.
The equation S Im =0 is valid for the actual field and a variations of Spk and SX3 are continuous in passing through the surface £. Therefore following 4-covariant compatibility equations for strong discontinuities is obtained from (26):
д L
д (др )_
= 0,
l^ - (дарк) dl
д (др pk )_
= 0. (27)
One can see, that the compatibility conditions on strong discontinuities surfaces contain the jumps of the energy-momentum 4-tensor and the Piola-Kirchhoff 4-tensor
43 = - a£k) ■ T3= ls3 - ) a&) ■ (28)
The compatibility conditions of jumps of the energy-momentum 4-tensor in view of (28) is rewritten in form
N3 [T a ] = + [3apk S3k ]N = 0. (29)
4
Relations (29) can be transformed by the compatibility conditions of the strong discontinuities of the Piola - Kirchhoff 4-tensor into
[L]N + [da pk ]S3kN3 = 0.
4
The only significant equation from the compatibility conditions for jumps of the energy-momentum 4-tensor by using the Hadamard- Thomas geometric first order compatibility conditions [8] [Oapk] = NaFk is obtained
[L] + S3k N3 Fk =0. (30)
4
The three-dimensional form of the compatibility conditions on the surface of strong discontinuity of the field is derived from the resulting 4-covariant form (27), (30):
[l] + (—G44k + n.S^k x-G^ ]] + п,[д^к ]]) = 0
Ф • Ц _
4 4 (31)
-G[S4k ]+ n,[S»k ]=0 (А, д =1, 2,3). ( )
44
Hereafter n^ stand for normal unit 3-vector.
The compatibility conditions for strong discontinuities are complemented by well known three-dimensional geometrical and kinematic Hadamard - Thomas compatibility conditions of the first and the second order [8] valid for an arbitrary field pk. Those are also due to Rankine and Hugoniot [12,13].
6. ACTION FOR MPTE-II CONTINUA
The theory GNII is the only thermodynamically correct theory which satisfies the principles listed in Sec. 1. The action taking into account of the polar microstructure can be adopted in the following form [7]
l = L(Xa, xj ,dj, ti, дАх3 ,д4dj ,д4$, дрxj ,дрdj ,дрti) dX1 dX2 dX3 dX4, (32)
a a a
where Xa (a = 1, 2,3) are the Lagrangian coordinates; xj (j = 1,2,3) are the Eulerian coordinates; dj (a = 1, 2,3) are micropolar directors associated with microvolume; ti is the thermal displacement field.
a
Action density (32) we define in the form
1 , 1 . ab
l = -(д4 xk )pkj (д4 xj ) + -(д4 d% ) I ij (д4 dj ) - P(X a,xj ,dj ,ti,д4 $,дах! , дад? , да ti),
2 2 a b a a
ab
wherein I ij — the microinertion tensor; pij — the mass density tensor; ф — volume density of the
ab ab ab
Helmholtz's free energy. For tensors I ij and pij the symmetry conditions pij = pji, I j = I ji are stated.
The field equations in this case read:
c)L
• - д4Pj = -(a = 1, 2,3, j = 1, 2,3), Механика 87
в
в
daMa + Aj - d4(Qj) = 0 (а =1, 2, 3, а = 1, 2, 3, j = 1, 2,3),
dajR + d4 S =
д L
дд
and are supplemented by the constitutive equations:
д L
S a- = — d
•j = д(daxj ) '
= d l j d(d4xj)
a
Qj =
д L
д (d4 dj)'
(а = 1, 2, 3),
д L
a dL
aj = ddj, a
S =
d l
д(d4д),
jR =
d l
д(daд) ■
Compatibility conditions on strong discontinuity surface propagating in MPTE-II continuum, according eqs. (30), (31), take the forms:
ab ab , a ,, , , 1 ,
-[d4diUijd4dj]] + (Glikd4d1 _ n.m.)(G[d4dk]] _ n%[d%dk]]) + -[d4xkpkj^4xj]] +
2 a b a a a 2
+(Gpikd4x1 _ n%S^k)(G[d4xk]] _ n%[d%xk]]) _ M + (Gs _ n.jR)(GM _ n.M) = 0, (33) Gpki[d4xk ]] = n%[S%l, G:Jki[d4 dk ]] = n^mt]], G[s] = nj ],
a
1 , 1 .ab
[l]] = -[d4xkPkjd4Xj]] + -[d4diUijd4dj] _ M (j,m, = 1, 2, 3).
2 2 a b
In the case of the propagating surface of strong discontinuity through unperturbed continuum the first of the compatibility equations (33) can be transformed to
1 ab ab a
-[d4di]] Iij[d4dj] + (Glik^d1 _ n%M%)(G[d4dk]] _ n%[d%dk]])+
2 a b a a a
+ -[d4Xk]]pkj[d4Xj]] + (Gpikd4X1 _ n%S%)(G[d4xk] _ n%[d%xk]])_ _M + (Gs _ n.jR)(GM _ n%[d%tf]]) = 0 (j, k, l, A, p = 1, 2,3).
The present work was partially supported by the Russian Foundation for Basic Research (project no. 13-01-00139 "Hyperbolic thermal waves in solids with microstructure") and the Ministry of Education and Science of Russia Federation grant given to Samara State Technical University (project no. 16.2518.2014/K).
References
1. Cosserat E., Cosserat F. Theories of Deformable Bodies. Paris, Scientific Library A. Hermann et Fils, 1909, 242 p.
2. Green A. E., Naghdi P. M. On Undamped Heat Waves in an Elastic Solid. J. Therm. Stress, 1992, vol. 15, pp. 253-264.
3. Green A. E., Naghdi P. M. Thermoelasticity without Energy Dissipation. J. Elasticity, 1993, vol. 31, pp. 189208.
4. Radayev Yu. N., Semenov D. A. Harmonic Coupled CTE-Thermoelastic Waves in a Free Cylindrical Waveguide. Vestn. Samar. Gos. Univ. Natural sciences ser., 2008, no. 8/1(67), pp. 411-459 (in Russian).
5. Kovalev V. A., Radayev Yu. N., Semenov D. A. Coupled Dynamic Problems in Hyperbolic Thermoelasticity. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2009, vol. 10, iss. 4, pt. 2, pp. 94-128 (in Russian).
6. Kovalev V. A., Radayev Yu. N. Wavenumbers of Plane
GNIII-Thermoelastic Waves and Inequality, Providing Their Normality. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2010, vol. 10, iss. 3, pp. 46—53 (in Russian).
7. Kovalev V. A., Radayev Yu. N. Wave Problems of Field Theory and Thermomechanics. Saratov, Saratov. Univ. Press, 2010, 328 p. (in Russian).
8. Thomas T. Y. Plastic Flow and Fracture in Solids. New York, Academic Press, 1961, 271 p.
9. Nowacki W. Theory of Asymmetric Elasticity. Oxford, Pergamon Press, 1986, 384 p.
10. Mindlin R. D., Tiersten H. F. Effects of couple-stresses in linear elasticity. Arch. for Rat. Mech. and Anal., 1962, vol. 11, iss. 1, pp. 415-448.
11. Eringen A. C. Microcontinuum field theories. Vol. 1. Foundations and Solids. Berlin, Heidelberg, New York, Springer, 1999, 325 p.
12. Rankine W. J. M. On the thermodynamic theory
a
a
a
of waves of finite longitudinal disturbance. Proc. of the Royal Society of London. London, The Royal Society, 1870, vol. 18, pp. 80-83.
13. Hugoniot P. H. Sur la propagation du mouvement dans les corps et spécialement dans les gaz parfaits. J. Ecole Polytechnique, 1887, vol. CLVII, pp. 3-98.
УДК 539.3
О слабых разрывах и уравнениях скачков на волновых поверхностях в микрополярных термоупругих континуумах
В. А. Ковалев1, Е. В. Мурашкин2, Ю. Н. Радаев3
1 Доктор физико-математических наук, профессор кафедры управления проектами и инвестициями, Московский городской университет управления Правительства Москвы, vlad_koval@mail.ru
2Кандидат физико-математических наук, доцент, Национальный исследовательский ядерный университет МИФИ, Москва, murashkin@ipmnet.ru, evmurashkin@gmail.com
3Доктор физико-математических наук, профессор, ведущий научный сотрудник, Институт проблем механики им. А. Ю. Ишлинского РАН, Москва, radayev@ipmnet.ru, y.radayev@gmail.com
Настоящее исследование посвящено проблеме распространения поверхностей сильных и слабых разрывов поступательных перемещений, микровращений и температуры в микрополярных (МП) термоупругих континуумах. В первой части статьи обсуждаются проблемы распространения слабых разрывов в MPTE континуумах первого типа. Геометрические и кинематические условия совместимости Адамара и Томаса используются для изучения возможных волновых поверхностей слабых разрывов. Слабые разрывы классифицируются в соответствии с пространственной ориентировкой векторов поляризации разрывов (DPVs). Показано, что поверхности слабых разрывов могут распространяться без слабых разрывов температурного поля. Вторая часть работы посвящена распространению поверхностей сильных разрывов полевых переменных в MPTE континуумах второго типа. Определяющие соотношения для гиперболических термоупругих микрополярных континуумов второго типа получены с помощью формализма теории поля. Специальная форма первой вариации интеграла действия позволяет получить 4-ковариантные условия скачков на волновых поверхностях. Трехмерная форма условий скачков на поверхности сильного разрыва поля выводится из ее четырехмерной ковариантной формы.
Ключевые слова: микрополярная термоупругость, континуум первого типа, континуум второго типа, слабый разрыв, сильный разрыв, ударная волна, продольная волна, поперечная волна, условие совместимости, скачок.
Работа выполнена при частичной финансовой поддержке РФФИ (проект № 13-01-00139 «Гиперболические тепловые волны в твердых телах с микроструктурой») и Минобрнауки РФ в рамках проектной части гос. задания ФГБОУ ВПО «СамГТУ» (проект № 16.2518.2014/(К)).
Библиографический список
1. Cosserat E. et F. Theorie des corps deformables. Paris : Librairie Scientifique A. Hermann et Fils, 1909. 226 p.
2. Green A. E., Naghdi P. M. On undamped heat waves in an elastic solid // J. Therm. Stress. 1992. Vol. 15. P. 253-264.
3. Green A. E., Naghdi P. M. Thermoelasticity without energy dissipation // J. Elasticity. 1993. Vol. 31. P. 189208.
4. Радаев Ю. Н., Семенов Д. А. Гармонические связанные СТЕ-термоупругие волны в свободном цилиндрическом волноводе // Вестн. Самар. гос. ун-та. Есте-ственнонауч. сер. 2008. № 8/1(67). С. 411-459.
5. Ковалев В. А., Радаев Ю. Н., Семенов Д. А. Связанные динамические задачи гиперболической термоупругости // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2009. Т. 9, вып. 4, ч. 2. С. 94-128.
6. Ковалев В. А., Радаев Ю. Н. Волновые числа плоских GNIII-термоупругих волн и неравенства, обеспечивающие их нормальность // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2010. Т. 10, вып. 3. C. 46-53.
7. Ковалев В. А., Радаев Ю. Н. Волновые задачи теории поля и термомеханика. Саратов : Изд-во Сарат. ун-та, 2010. 328 с.
8. Thomas T. Y. Plastic Flow and Fracture in Solids. N. Y. : Academic Press, 1961. 271 p.
9. Nowacki W. Theory of Asymmetric Elasticity. Oxford : Pergamon Press, 1986. 384 p.
10. Mindlin R. D., Tiersten H. F. Effects of couple-stresses in linear elasticity // Arch. for Rat. Mech. and Anal. 1962. Vol. 11, iss. 1. P 415-448.
11. Eringen А. С. Microcontinuum field theories. Vol. 1. Foundations and Solids. Berlin ; Heidelberg ; N.Y. : Springer, 1999. 325 p.
12. Rankine W. J. M. On the thermodynamic theory of waves of finite longitudinal disturbance // Proc. of the Royal Society of London. London : The Royal Society, 1870. Vol. 18. P. 80-83.
13. Hugoniot P. H. Sur la propagation du mouvement dans les corps et specialement dans les gaz parfaits // J. Ecole Polytechnique. 1887. Vol. CLVII. P. 3-98.
