Printer Friendly

Nonlinear waves in dissipative microstructured two-dimensional solids/mittelineaarsed lained dissipatiivsetes mikrostruktuursetes 2-dimensionaalsetes tahkistes.


The theory of continua with microstructures (see, e.g., [1,2]) can be used to model the behaviour of real materials such as granular materials, polycrystalline solids, ceramic composites, and materials with microdefects. One of the main features of the theory is the possibility of taking into account intrinsic space scales, namely, the size of grains, distance between microcracks, etc. Nonlinearity and dispersion cannot be avoided as phenomenological effects and mathematical features intrinsic to such phenomena. To describe granular materials, we must also introduce dissipation, which is mainly due to interaction between neighbouring grains.

The model of vector microstructures ([3-5]) seems to be useful for obtaining general field equations applicable to such materials. For simplicity, we restrict our attention to plane granular media. The natural model is the Cosserat theory, where the vector microstructure is described by a triad of orthonormal vectors. In field equations conservative stresses appear, i.e. stresses related to a generalized strain energy, and dissipative stresses, which are assumed to be linear in strain velocities.

Wave propagation can be studied by means of the well-known perturbation technique and slaving principle (see [6]). This approach allows us to reduce the set of general field equations to one equation in the ruling variable; hence we may estimate the possibility of propagation of solitary waves, depending on the balance of nonlinearity, dispersion, and dissipation. As a rule, dispersion is required for the existence of the bell-shaped solitary waves in an elastic microstructured medium, while dissipation is expected to be responsible for a saturation, which prevents unbounded growth of the bell-shaped solitary wave.

Our aim is to solve analytically the two-dimensional in space nonlinear wave propagation problem, using the simplest possible but informative description of internal dissipative features. Anaytical results are of considerable interest also as a tentative fault detection test for any numerical simulation in the 1+2D problem. The numerical solution in a 1+1D problem considered in [7] was based on the pseudospectral method supported by the analytical solution for the linear case to show how dissipative effects on various scales affect the harmonic wave.

We show how the influence of dissipation in the 1+2D problem in comparison with the simple case, where only nonlinearity and dispersion are taken into account, appears through a coefficient in the final solution, having a general form in terms of the Weierstrass elliptic function. It provides periodical and solitary wave solutions as appropriate limits. Moreover, it enters essentially the soliton solution to the 1+2D problem, and it seems feasible to claim that it can affect the amplitude evolution. Physical importance of the results obtained is in getting a detectable signal--the soliton, which is elusive in experiments without preliminary analytical solution (see [8,9]).


It is well known (see [3,4,10-12]) that the Cosserat solids can be described as nonlinear elastic solids C with a vectorial microstructure constrained to be a rigid triad [d.sub.[alpha]] = [d.sub.[alpha]]([X.sup.[beta]], t), [alpha], [beta] = 1, 2, 3, where [X.sup.[beta]] are material coordinates in a reference configuration [C.sub.*], such that [d.sub.[alpha]] x [d.sub.[beta]] = [[delta].sub.[alpha][beta]], [for all]t.

In the plane case we have [alpha], [beta] = 1, 2 and we can express the directors in terms of a rotation angle [theta] only:

[d.sub.1] = d = cos [theta][e.sub.1] + sin [theta][e.sub.2], [d.sub.2] = v = -sin [theta][e.sub.1] + cos [theta][e.sub.2], (1)

where {[e.sub.i]}, i = 1, 2 is any orthonormal spatial basis, 0 [less than or equal to] [theta] [less than or equal to] 2[pi]. We have

[??] = [??]v, [d.sub.,h] = [[theta].sub.,h][degrees], (2)

and writing r = [x.sup.h]([X.sup.h], t)[e.sub.h], we obtain for kinetic and strain energy densities


If we avoid body forces, the dissipative field equations yield ([13]):


where the terms [[??].sup.h], [[??].sup.i], [??] represent the dissipation.


We study the simplest model, where the dissipation is due to the rotation [??] only, through the term [??] = -F([[theta].sub.x] + [[theta].sub.y])t.

For the sake of simplicity, we shall use the notation: [X.sup.1] = x, [X.sup.2] = y, [x.sup.1] = u, [x.sup.2] = v. Hence we consider the vector r = r(x, y, t) = u(x, y, t)[e.sub.1] + v(x, y, t)[e.sub.2] for the macrostructure and, for the microstructure, the function [theta] = [theta](x, y, t) that represents the angle of rotation of the particle with respect to the fixed basis. In the following the subscripts x, y, t will denote differentiations.

The kinetic energy density reads

T = 1/2 [[rho] ([u.sup.2.sub.t] + [v.sup.2.sub.t]) + I[[theta].sup.2.sub.t]].

The strain energy density is chosen in the form


Since we are mainly interested in the estimation of [theta], we introduce a new variable U = u + v. The field equations can be written as follows:


Consider the dimensionless form of system (5), and apply the slaving principle (for all details see [13]). In addition to dimensionless quantity [theta], for further analysis the dimensionless variables are introduced:

w = U/[w.sub.0], X = x/L, Y = y/L, T = [c.sup.2.sub.0]/L t,

where [c.sup.2.sub.0], [w.sub.0], L are physically meaningful constants (velocity, intensity, and wavelength of the initial excitation). We also need a scale for the microstructure l. Then two dimensionless parameters can be introduced: [delta] ~ [(l/L).sup.2] characterizing the relation between the microstructure and the wavelength and [epsilon] ~ ([w.sub.0]/L) accounting for elastic strain, where [delta] is the relevant characteristic length.

Following [14], we suppose I = [rho][l.sup.2][I.sup.*], C = [l.sup.2][C.sup.*], D = [l.sup.2][D.sup.*], F = [l.sup.2][F.sup.*], where [I.sup.*] is dimensionless and [C.sup.*] and [D.sup.*] have the dimension of the stress. We consider the expansion in terms of the characteristic length [delta]: [theta] = [[theta].sub.0] + [delta][[theta].sub.1] + ... and we impose the coefficients of powers of [delta] to be equal. Following the approximation for [theta] in terms of w and derivatives, we get


Finally, we obtain the governing equation for w, i.e., for a nondimensional U = u + v:


Introducing the new variable z = X + kY - cT, we consider the 4th-order ODE for function w(Z):


where we introduce the following coefficients [??], [??], [??]:


By setting w' = y, where y = y(z), integrating twice Eq. (6) and rescaling the coefficients, we get the final equation:

[(y').sup.2] + a[y.sup.2] + b[y.sup.3] + dy + e = 0, (7)

which is the Weierstrass equation, having a general solution in terms of y = [A.sub.p] + B, where p = p(z + [z.sub.0, [g.sub.2], [g.sub.3]). We define y = [A.sub.p] + B, [right arrow] y' = [A.sub.p'], since [p'.sup.2] = (4[p.sub.3] - [g.sub.2]} - [g.sub.3]), and have [(y').sup.2] = [A.sup.2](4[p.sub.3] - [g.sub.2]} - [g.sub.3]).

Replacing them into Eq. (7), we obtain

[A.sup.2](4[p.sub.3] - [g.sub.2]} - [g.sub.3]) + a([A.sup.2][p.sup.2] + [B.sup.2] + 2ABp) +b([A.sup.3][p.sub.3] + [B.sup.3] + 3[A.sup.2]B[p.sup.2] + 3A[B.sup.2]}) + d(Ap + B) + e = 0:

The coefficients A, B, [g.sub.2], [g.sub.3] are determined by making independently the coefficients of each order of p and p' equal to zero:



y = -4/b p(z + [z.sub.0], [g.sub.2], [g.sub.3]) - a/3b, (8)

where [g.sub.2] and [g.sub.3] are defined above. Eventually we have


The last step of this approach consists in the evaluation of the integral:


We consider, again, the expression (9), which is general, discontinuous, and semibounded from below or above. Solutions of this kind can be useful for tentative validation tests during numerical simulations. Of main physical interest are bounded continuous solutions, we obtain them as appropriate limits of (9) after some nontrivial algebra. Following [15], we introduce

(i) 2w and 2w' as primitive periods of p,

(ii) the discriminant [DELTA] = [g.sup.3.sub.2] - 27[g.sup.2.sub.3],

(iii) the roots [e.sub.i] = p([w.sub.i]), i = 1, 2, 3 of the equation 4[e.sup.3] - [g.sub.2]e - [g.sub.3] = 0, where [w.sub.1] = w, [w.sub.2] = w + w', [w.sub.3] = w'.

Reduction of the doubly periodic Weierstrass function p to a set of single periodic Jacobian elliptic functions is based on the following relationship between p and cn, sn with modulus M:


Moreover, (cf. [16]), if [e.sub.1], [e.sub.2], [e.sub.3] are any three numbers whose sum is zero, and if we write

y = [e.sub.3] + [e.sub.1] - [e.sub.3]/[sn.sup.2](z[square root of [e.sub.1] - [e.sub.3]]; M),

we obtain the following relation between the Weierstrass function and the Jacobian elliptic function: p(z, [g.sub.2], [g.sub.3]) = [e.sub.3] + ([e.sub.1] - [e.sub.3])[ns.sup.2] (z[square root of [e.sub.1] - [e.sub.3])]| M). It is well known that the behaviour of p depends on the sign and value of [DELTA], which allows us to extract two cases of main interest.

3.1. Two-dimensional solitary wave solution

In the case of [DELTA] = 0 one of the periods is infinite: w = [infinity] or w' = i[infinity] (the trivial case w = - iw' = 1 will be excluded). The first case w = 1 corresponds to [e.sub.1] = [e.sub.2] [not equal to] [e.sub.3]. Since [e.sub.1] + [e.sub.2] + [e.sub.3] = 0, introducing the condition [e.sub.1] = [e.sub.2] [equivalent to] E, we have:

[e.sub.3] = -2E, [g.sub.2] = 12E2, [g.sub.3] - 8E3, w0 = i1/4=p12E:

Therefore the equation y = -4/b p(z + [z.sub.0]; [g.sub.2]; [g.sub.3]) - a/3b yields


In our case M = 1 and the limiting value sn(z, 1) = tanh z. Then we obtain the exact 1+2D solitary wave solution y(z) to the original problem (9):

y(z) = -a/3b + 8E/b -12E/b [tanh.sup.-2] ([square root of 3E] (z + [c.sub.0])). (11)

Integrating (10), where y(z) is given by (11), we obtain finally

w = (-a/3b + 8E/b) z + [square root of 48E/b] coth ([square root of 3E] (z + [c.sub.0])):

The value w' = [infinitive] corresponds to [e.sub.1] [not equal to] [e.sub.2] = [e.sub.3] and does not lead to bounded solutions (however, it very often occurs in numerical simulation).

In order to get the explicit expression for [theta], we recall that


For simplicity we set


and obtain [theta] = [c.sub.1]W' + [c.sub.2]W''', namely:


3.2. Periodic bounded solutions in the case [DELTA] > 0

In this case there is a pair of primitive periods 2w and 2w' such that w is a real and w' a pure imaginary semiperiod of the p function. Assuming [DELTA] > 0, all roots [e.sub.i], i = 1, 2, 3 are real and different, [e.sub.1] > [e.sub.2] > [e.sub.3], [e.sub.1] > 0, [e.sub.3] < 0, from y = [square root of -4/b] p(z + [z.sub.0], [g.sub.2], [g.sub.3]) - a/3b we obtain the bounded periodical 1+2D cnoidal wave solution:


No bounded solutions p(z) were found from (8) when [DELTA] < 0.


We aimed to show how the influence of dissipation in the nonlinear two-dimensional problem may be encapsulated in a compact form of the coefficient in the final ODE. From that equation we found a new general solution in terms of the Weierstrass elliptic function. It provides many periodical and solitary wave solutions as appropriate limits. The approach used does not require any reduction to evolutionary equation like the KP equation, which is widely used but restrictive in both derivation and solution of a two-dimensional nonlinear wave problem. Moreover, our approach yields an essentially new soliton solution to the 1+2D problem. Thus it seems feasible to claim that the importance of our results in physics lies in the prediction of a detectable signal--the soliton, which is elusive in experiments without preliminary analytical solution.


This research was supported by the Italian MIUR-PRIN Project 2005: "Mathematical Models in Material Sciences".

Received 15 January 2007


[1.] Maugin, G. A. Nonlinear Waves in Elastic Crystals. Oxford University Press, UK, 1999.

[2.] Phillips, R. Crystals, Defects and Microstructures. Modeling Across Scales. Cambridge University Press, Cambridge, 2001.

[3.] Pastrone, F. Mathematical models of microstructured solids. Lect. Notes Mech. 4/04, Tallinn Tech. Univ. Tallinn, 2004.

[4.] Pastrone, F. Waves in solids with vectorial microstructure. Proc. Estonian Acad. Sci. Phys. Math., 2003, 52, 21-29.

[5.] Pastrone, F. Wave propagation in microstructured solids. Math. Mech. Solids, 2005, 10, 349-357.

[6.] Engelbrecht, J., Berezovski, A., Pastrone, F. and Braun, M. Waves in microstructured materials and dispersion. Phil. Mag., 2005, 85, 4127-4141.

[7.] Sillat, T. and Engelbrecht, J. Wave propagation in dissipative microstructured materials. Proc. Estonian Acad. Sci. Phys. Math., 2003, 52, 103-114.

[8.] Janno, J. and Engelbrecht, J. Solitary waves in nonlinear microstructured materials. J. Phys. A: Math Gen., 2005, 38, 5159-5172.

[9.] Giovine, P. and Oliveri, F. Dynamics and wave propagation in dilatant granular materials. Meccanica, 1995, 30.

[10.] Casasso, A. and Pastrone, F. Nonlinear waves in microstructured solids and complex structures. In Proceedings of the International Seminar "Days on Diffraction 2004", Saint Petersburg, Russia. Faculty of Physics, SpbU, 2004, 43-51.

[11.] Casasso, A. and Pastrone, F. Nonlinear waves motion in complex elastic structures. Rend. Circ. Mat. Palermo, Serie II, Suppl. 78, 2006, 45-58.

[12.] Pastrone, F. Microstructures and granular media. Rend. Sem. Mat. Univ. Polit. Torino, 2007, 65, 87-95 .

[13.] Casasso, A. and Pastrone, F. Nonlinear waves in plane granular media. In Proceedings of the International Seminar "Days on Diffraction 2005", Saint Petersburg, Russia. 2005, 30-39.

[14.] Porubov, A. V. and Pastrone, F. Nonlinear bell-shaped and kink-shaped strain waves in microstructured solids. Int. J. Non-Linear Mech., 2004, 39, 1289-1299.

[15.] Samsonov, A. M. Strain Soliton in Solids and How to Construct Them. Chapman & Hall/CRC, New York, 2000.

[16.] Whittaker, E. T. and Watson, G. N. A Course of Modern Analysis. The MacMillan Company, New York, 1946.

Alessia Casasso (a), Franco Pastrone (a), and Alexander M. Samsonov (b)

(a) Dipartimento di Matematica, Universita di Torino, Via C Alberto 10, Torino 10123, Italy,,

(b) A. F. Ioffe Physico-Technical Institute, Russian Academy of Sciences, 194021 St Petersburg, Russia,
COPYRIGHT 2007 Estonian Academy Publishers
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2007 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Casasso, Alessia; Pastrone, Franco; Samsonov, Alexander M.
Publication:Proceedings of the Estonian Academy of Science Physics/Mathematics
Date:Jun 1, 2007
Previous Article:Propagation of an austenite-martensite interface in a thermal gradient/austeniit-martensiitpiirpinna levimine soojusgradiendi toimel.
Next Article:On the propagation of localized perturbations in media with microstructure/lokaliseeritud hairituste levist mikrostruktuurses keskkonnas.

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters