CC BY
71
UDC 539.3
ON NEW PHASE INTERMEDIATE LAYERS
IN NANOCOMPOSITES AS A SOURCE OF INCREASING THE ELASTIC MODULI
R. A. Filippov1, A.B. Freidin1, E.N. Vilchevskaya1
institute for Problems in Mechanical Engineering, RAS, Saint-Petersburg [email protected], [email protected], [email protected]
PACS 46.25.Cc, 68.35.Rh
We study an opportunity to increase elastic moduli of a nanocomposite due to stress-induced phase transformations which lead to the formation of intermediate new phase layers around nanoparticles. These layers enlarge the effective size of the particles which from now become inclusions made up of kernels (the initial nanoparticles) enclosed by new phase layers shells. Increasing the volume fraction of the inclusions can change the effective elastic moduli of the composite much more than one could expect in a case of the composite with a small volume fraction of initial nanoparticles. As an example we consider an isotropic composite with spherical particles under hydrostatic loading. We begin with considering the new phase formation around an isolated inclusion including the interface stability analysis. We show that stable two-phase states are impossible if both elastic moduli of the matrix increase due to phase transition and possible if the bulk modulus increases and the shear module decreases. Then, basing on a self-consistent approach, we describe the new phase formation around spatially distributed particles and study how external strains effects the new phase areas growth. Finally we demonstrate that the new phase layers formation can lead to increasing the effective bulk modulus of the composite.
Keywords: nanocomposites, intermediate layers, effective elastic moduli, phase transformations, stability. 1. Introduction
Carbon nanotube and fulleren reinforced nanocomposites have received increased attention in recent years because of their mechanical, electrical, optical and magnetoelectrical properties [4,15,16]. One of the surprising facts is that a small volume fraction of nano-inclusions may lead to significant effects which cannot be predicted by common mathematical model of composite materials made of a matrix and uniform inclusions. In a number of papers (see, e.g., [1,4]) it was suggested that these effects may be related to the quality of the boundaries which separate the inclusions and matrix. But the theories of composite material usually imply ideal contacts at the inclusions boundaries and even in this case cannot give the effects declared.
To explain how a small volume fraction of particles may lead to disproportionately significant changes of mechanical properties of nanocomposites we suppose that this may be a result of the formation of intermediate layers surrounding the particles. The major hypothesis is that these layers appear as a result of phase transformation in matrix around the particles and the phase interfaces must satisfy thermodynamic equilibrium conditions. The formations of new phase layers increases the "effective" sizes of the particles. As a result the effective volume concentration of inclusions increases and, thus, effective elastic moduli change.
To simplify the derivations we consider a composite with isotropic spherical inclusions in isotropic matrix. This gives us a foreseeable script of general consideration and analytical basement for further verifications of numerical simulations.
Note that stress-induced phase transformations can be induced by external all-round compression/tension as well as by thermal stresses appeared because of the difference between thermal expansion coefficients of the inclusions and surrounding material. In fact, new phase areas can be formed under combined thermo-mechanical actions at the stage of the composite manufacturing. Since from the mathematical point of view thermal stresses can be reduced to stresses produced by mechanical loading, we consider further only external all-round tension/compression.
The following questions are being discussed within the above formulated elementary problem:
1. Is it possible to form the stable intermediate layers of a new phase if at least some of elastic moduli increase?
Given average strain and temperature, equilibrium and stable two-phase states are the microstructures which correspond to minimal value of the Helmholtz free energy of a body that differs from the strain energy only by a term related with free energies of the phases in stress free states and proportional to the new phase volume. The strain energy cannot decrease if all elastic moduli of the matrix increase due to phase transformation around the particles. Thus, we have to examine the case when some elastic moduli increase but the others decrease. From the mechanical point of view this means that we have to deal with a case for which the jump of the elasticity tensor due to phase transition is a non-signdefinite tensor, and in any case we need a stability analysis that is carried out in the present paper.
2. Is it possible to provide the controlled growth of the effective elastic moduli?
This question arises because, generally speaking, the stability of the two-phase state of the matrix can be lost and the instability may lead, particularly, to a spontaneous phase transformation of the whole matrix. That is why we find the interval of external strains at which two-phase states adjust to external strains and remain stable.
3. Is it possible to reduce stress concentration in the matrix due to the formation of new phase transition layers?
To answer this question we examine how the phase transformation effects the strain energy density spatial distribution. We demonstrate that in a case of the stable interface the phase transformation leads to the energy density decrease in the particle itself and in the matrix on the side of the parent phase. We also show the opportunity to reduce the energy density in all areas of the body. The intensity of the shear stresses also decreases.
2. Constitutive equations
We consider stress-induced phase transformations of martensite type accompanied by the transformation strain and changes of elastic moduli. Two-phase deformation of an elastic body is studied by a small strain approach.
Let new phase areas V+ be bounded by interfaces r. A problem on the equilibrium two-phase configurations of an elastic body is reduced to finding position of unknown interfaces r and displacement u(x) sufficiently smooth at material points x E r, continuous across r, and satisfying equilibrium and boundary conditions including conditions on the interface:
x E r : V o- = 0, 9 = const, (1)
x E r : [u] = 0, [a] ■ n = 0, (2)
[/] - (<r> : [e] = 0, (<T> = i(<T+ + <T_) (3)
where e and a are strain and stress tensors, B is the temperature, f = f (s,9) is the Helmholtz free energy volume density, n is the unit normal to the interface r, super- or subscripts "—" and "+" identify the values being taken for a material in the "—" and "+" phase states, respectively,
the brackets [ • ] = (•)+ — (•)- denote the jump across r. Body forces, thermoelastic effects and surface energy are not taken into account. Note that from the displacement and traction continuity it follows that {a} : [e] = : [e].
The conditions (2) are common contact conditions and can be satisfied on any contact surfaces in a body, as it is in composite materials. The additional thermodynamic condition (3) (known as the Maxwell relation) is related to an additional degree of freedom caused by an unknown phase boundary. This condition is the necessary (but not sufficient) condition for the energy minimization with respect to the position of the interface. Namely the thermodynamic equilibrium condition (3) imposes the restrictions on the existence and geometry of new phase areas (see, e.g., [9,10]).
We consider a two-well strain-energy function:
f (e, e) = min { f-(e,d),f+(e,d)}
(4)
f ±(e,
f?{e) + -{e-e?±):C±:{e-e?±).
Upper and lower scripts "+" and "—" in relationships correspond to each other. ID-sketch of the free energy density dependence is shown in Fig. 1
e1'
£
a b
Fig. 1. Free energy volume density: a - case 1, b - case 2
By (4), the constitutive equations take the form:
where £
a±(e) = C± :(e — e±), e
8- = [e : v(e) > 0}, = (s : p(e) < 0}, p(e) = f+(e) — f-(e)
(5)
are the domains of definition of the phases "—" and "+" in strain space.
The parameters C±, f± and ep± are the elasticity tensors, free energy densities and strain tensors in unstressed phases "±", respectively. If ep+ = 0 or e— = 0, then [ep] = ep is the transformation eigenstrain.
£
Fig. 2. A particle enclosed by a new phase intermediate layer
We emphasize that one of the unstressed states can be only hypothetical. The strain e = £p+ may belong to the domain of the phase "—", not to the phase "+", if f+(e+) > f ~(£+). Further we say that we are dealing with the case 1, if f +(e+) < f ) (Fig. 1 a) and with the case 2, if f +(e+) >f"(e+) (Fig. 1 b).
The thermodynamic condition can be rewritten in the form that relates the strain on one side of the interface and the normal to the interface (see [7] or [2] and reference therein):
x(e±,7,n)^7 + ^[ep:C:ep] + ie±:C1:e±-e±:[C:eP]±iq±:KT(n):q± = 0, (6) where
q± = Ci : e± — [C : £p], KT(n) = {n 0 GT(n) 0 n}s, GT(n) = (n • CT • n)-1, Ci = C+ — C-, Bi = C-1 — C-1, 7 =[/c].
The parameter 7 acts as the temperature.
If the inverse tensor C-1 exists then (3) can be rewritten in terms of the tensor q taken on one side of the interface [14]:
X(q±, 7*, n) 4 7, + I q± : (C^ ± K») : q± = 0, 7* = 7 + \ [e"] : Br1 : (7)
Thermodynamic condition (6) cannot be satisfied for all strains. This leads to the concept of phase transitions zones (PTZs) [9,11,12], which are the areas formed in strain space by all strains that may exist in the body at the equilibrium phase boundaries. In other words, PTZ is formed by the strains, for which the equation (6) can be solved with respect to the unit normal n. Given material parameters and the parameter 7, the PTZ presents all possible local orientations of the equilibrium phase boundaries and strain jumps across the interfaces. Thus, the PTZ is the passport of a material which undergoes phase transformation (see, e.g., [3,9]). In this paper, the construction of PTZs is used to reject the equilibrium but unstable solutions.
Further we consider a matrix that can undergo phase transformations and inclusions, which act as stress concentrators and induce phase transformation. The material of the matrix
is piece-wise linear elastic described by the constitutive equations (4), (5). Material of the inclusion is also linear-elastic, its elasticity tensor is C and the free energy density is given by the dependence
Si (e) = \e-.Ct-.e. (8)
3. Spherically-symmetric two-phase deformations induced by a spherical inclusion
We consider an unbounded medium with a spherical inclusion Vi of the radius Ri under all-round tension/compression. Let e— = 0, e+ = ep = (§p/3)E, and e0 = ($0/3)E be the external strain. Suppose that there exists a two-phase deformation with one spherical interface (Fig. 2) and find the equilibrium interface radius Rr in dependence of external strain and material parameters.
Radial displacement u in areas Vi} V+, V— is given by the formula
D
u = Ar + — (9)
where dependencies of the integration constants A = [Ai7 A+,A—} and D = {Di} D+, D—} on e0 and the inverse of the relative interface radius
P = f^, pe[ 0,1]
can be found from the contact conditions at the particle boundary and the interface, the condition at the vicinity, and the boundedness of the displacement in the center:
r = Ri, r = Rr : [u] = 0, [ar] = 0
^r V3 = u |r=0<
and ar, er, e^ are stress and strains in spherical coordinates.
It can be found that the integration constants are given by the following relationships
^ = + A = 0, = ~z)t
Z ' 1 Z
A. = e0, ^ = D_=D£_e0 + Dt
K-n
D£_ 3(dtk+z - djkx) Z Dt = d+( 1 - z)t Z (10)
= A£+£Q + A+, = d+ dZ Z /11? (dt ~ - z)t Z
D+ Rf = D+ = D+£o + D+, D£+ = ^d-kii ' Z ' + d~t Z
where
Z = d+ d+ — k+
li
z ^ p , k\ ^ k+ k—, k+i ^ k+ kî, ^ P—, ^ ^
d+ = 3k± + 4p,+ , df = 3k± + df = 3ki + 4^f
The mechanical equilibrium equation (1)1 and contact conditions (2) are satisfied at any interface radius p if the solution is taken in the form (9) with integration constants given by formulae (10). The dependence of the thermodynamically equilibrium inverse interface radius
p* on external strain e0 is determined by the thermodynamic condition (7) that, after substituting strains obtained by the differentiation of the displacement (9), takes the form of the algebraic equation
o(t - 3kiA+ + A^izD(+_ ou a )2 ^
27, + -A--tL + t* + = 0
d-
1* =1
ki
k k. i i)1' 2k\
(11)
r« 0
So
Fig. 3. The inverse relative equilibrium interface radius p in dependence on e0,
¡i- > ¡i+.
The dependence of the inverse relative equilibrium interface radius p* on e0 is shown in Fig. 3 for the case > 0, ¡i1 < 0. New phase nucleation starts at the particle boundary at e0 = es. The value es can be found from (11) if substitute p =1 into the formulae for D+ and A+.
The whole matrix transforms into the phase state "+" at e0 = e?. Substituting p = 0 into (11) leads to to the formula
1
(12)
If < 0,kl > 0 then the sign "-" is to be taken in (12).
4. Energy changes due to new phase areas growth and the interface stability
The interaction energy is determined as
En
np - n-
(13)
where np is the Gibbs free energy of a body with an inclusion that has the interface inverse radius p embedded into two-phase matrix, n- is the Gibbs free energy of a body with an inclusion embedded into one-phase matrix in the phase state "—" at the same boundary conditions as np is calculated.
f
f
In the case under consideration
EP = f (h(e) - h(e0))dV + [ (f+(e) - f-(e-)) dV + f (f-(e) - f-(e-)) dV
IVi JV+ JV—
where e is a current strain field in a case of two-phase matrix with the interface inverse radius p, £0 is a strain inside the inclusion when the matrix is in the phase state "—", e- is a strain in one-phase matrix in the "—" phase state.
Following [2] we derive from (4) and (8) that in a case of an arbitrary interface radius
¿3 = (t + - 3e0) + ¿(3he0 - ^ + 2»1D+D°_(1 - z),
50 = ^ Do = 3eo(k_ - h)RÏ ' " 3ki +
The equilibrium interface radius p*, i.e. the radius satisfying the thermodynamic condition, corresponds to the extremum point on the dependence Ep(p \ e0).
Examples of the dependencies of the interaction energy on the interface radius are given in Fig. 4 b1 (case 1) and Fig. 4 62(case 2) (/_ > /+). The results for two values of the external strains e0 = e01 (solid lines) and £0 = £02 (dotted lines) are shown, £02 < £01 < 0. Minimum points in Fig. 4 1, 2 demonstrate that the equilibrium solutions are stable with respect to interface radius perturbations if > /+. If < ¡i+ then similar calculations show that the dependence of the interaction energy has a maximum at the new phase spherical layer equilibrium radius, and this means instability of the interface. That is why further we discuss in detail phase transformation, resulting in decreasing the shear modulus but the increasing the bulk modulus.
So, if the bulk modulus increases, while the shear modulus decreases due to phase transformation, i.e. k1 > 1, ¡i1 < 0, then the interface is stable with respect to perturbations of the interface radius. However, when certain radius is reached, the interface becomes unstable with respect to axisymmetric perturbations. This instability can be detected with making use of an additional stability analysis basing on the procedure developed in [10,13]. The instability of equilibrium two-phase deformations have not been detected if the strain at the interface corresponded to the outer boundaries of the PTZs. The fact that strains at the stable interface must belong to the PTZ boundaries also follows from the recent paper [17].
Fig. 4a1, a2 shows how strains at the interface are related with the PTZ at different external strains. Black dots denote to the local strains which take place at the equilibrium interface at the external strain e01. These strains are located at the external boundaries of PTZ. Crosses mark the local strains at the external strain £02 (e02 < £01 < 0). One of these strains is located at the internal boundaries of PTZ.
In the case under consideration the strains belong to the outer PTZ boundaries if
(£- - £-) \r > 0 (14)
This inequality allows us to estimate the interval of external strains at which two-phase states remain stable. Straightforward calculations show that the transformation interval \ £ ? - £ s\ is larger in the case 2 that is if the value ^(9) is chosen so that f+(sp) > f~(sp).
Note that the stable homogeneous spherical new phase areas can exist in a uniform media under external all-round tension/compresion only if ¡i+ > /j,_ [2,5,10,13]. Therefore, depending on the relationship between shear moduli of the phases, two types of new phase nucleation are possible. If < then inhomogeneities (stress concentrators) are the points of the new phase growth. If > then new phase areas can appear in the domains of homogeneity of the material.
5. Free energy density spatial re-distribution and energy release due to phase transformations
Two-phase deformations are always accompanied with some stress concentration at the interface. To clarify why two-phase deformations can be preferential from the energy minimization point of view we study how the free energy density is redistributed in a body if the phase transformation takes place around the particles.
Fig. 4 c1 and c2 shows the spatial distributions of the difference
' U(e<>) - fi(e®), 0 <r<Rz ^(r) = f-(ep) - f-(e-), <r<Rr U(e') - f-(e-), r>Rr
between the free energy density in the case of the equilibrium two-phase state of a matrix with embedded inclusion and the free energy density in the case of one-phase state of the matrix with embedded inclusion at the same boundary conditions for the cases 1 and 2. As in the previous section, solid lines and dotted lines in Fig. 4 c1 and c2 correspond to external strains e01 and e02,
£02 < £01 < 0.
£a = £.
E
ip(r)
\
- — — ,
\
\
\
s
0
0
£ Z D r
a2 &2 ^ C2
Fig. 4. a1 and a2 - the PTZ and local strains at the interfaces, b1 and b2 - the dependencies of the interaction energy E on z = p3 (the cube of the inverse interface radius), c1 and c2 - the spatial redistribution of the free energy density in the cases 1 and 2, correspondingly.
Negative values of ^(r) mean that the local energy release takes place due to equilibrium new phase zone formation around the inclusion. In cases 1 and 2 the free energy density decreases both in initial inclusion and outside the coated one. In case 2 the free energy density is also decreased in intermediate layer and the relative contribution of this decreasing in body energy is very significant.
e* (p, so)
Fig. 5. The effective field method.
Note that if the strain energy density decreases in the vicinity of the inclusion then the stress concentration decreases. In this case the phase transformation acts as a reinforcing factor.
6. Composites with a matrix undergoing phase transformations
We consider a composite material with spherical inclusions. As in the previous sections the material of matrix can undergo phase transformations. The constitutive relationships of the material of the matrix and the particles are the same as in the case of the isolated inclusion, (4) and (8), correspondingly.
Following a self-consistent method [6,8] we can consider every inclusion as an isolated one embedded into the homogeneous medium with elasticity tensor C_. The influence of other inclusions is taken into account by introducing the effective field e* that is the sum of the external field e0 and the field induced by surrounding inclusions (Fig. 5). The main hypothesis of the method is that the field e* is constant and the same for all the inclusions. We also assume statistical independence of the inclusion properties from their spatial locations and use ergodic properties of the functions considered.
The following equation for the effective field can be obtained [18]:
e* = eo + pA :(Po : e* - Pp0) (15)
where p is a volume fraction of coated inclusions, p = p0 + p+, p0 is a volume fraction of the particles (initial inclusions), p+ is a new phase volume fraction, the tensor A take the form [6]
1 — c6_ 5 — _ 3k_ + L—
A =-EE +-(I — -EE), ae_ h
9l_ 15l_ 3 3k_ +
where I is the isotropic unity fourth-range tensor, Iijki = 1/2(SikSji + SjkSn).
The effective field must be spherical for hydrostatic loads e* = e*E. We have
e* = P° = {klA£+P+ ~ ^^ = + " klA+p+) ^ (16)
Af, A£+,Af, A+ are determined by (10) at z = p0/p+.
If p ^ 0 then the solution of a composite problem becomes the solution for the isolated inclusion given in previous sections.
After averaging the expression for the stresses and strains we derive the formula for the effective elastic moduli of the composite. The tensor of effective elastic moduli takes the form:
C* = A:*EE + 2/j,* [ I - ^EE
(i-Iee)
where the dependencies of effective bulk and shear moduli on volume fractions of the particles and the thickness of the new phase layers is given by formulae
d-Po 3d~_Po
k* = k---7Z-rrr, P* = P- —
dZ - 3P0 ^" 2(dZ - 3P0)
The dependencies of relative effective bulk modulus k*/k°, shear modulus and
Young's modulus E^/E* on the volume fraction of layered inclusions are shown in Fig. 6. Here k°, p®, and El are the bulk, shear and Young's moduli of the initial composite without intermediate new phase layers. Fig. 7 shows the dependence of the equilibrium inclusions volume fraction p* on the external strain. Both figures are pictured for the case 2, kZ = 20, pz = 9, k+ = 140, = 6, ki = 500, ¡ii = 170, № = -0.08, 7 = 0.07, the initial volume fraction of the particles 0.02.
Fig. 6. The dependencies of the relative effective bulk modulus k*/k°, shear modulus and Young's modulus E*/E° on the volume fraction of the layered inclusions.
The solid parts of the lines correspond to the equilibrium volume fraction of layered inclusions and corresponding external strains at which the local strains at the interface of the effective inclusion belong to the external PTZ-boundaries. We assume that corresponding interfaces are stable. The dotted parts of the lines correspond to external strain and volume fractions for which at least one of the strains at the effective interface does not belong to the external PTZ-boundary. The stability of the effective interface is lost in this case. Chosen parameters k±, ki, p±, 7 and volume fraction of the particle p0 = 0.02, two-phase structure remains stable until the equilibrium volume fraction rp* of the layered inclusions reaches 0.1. The effective bulk modulus of the composite demonstrates 10 percent growth. Both Young's modulus and shear moduli decrease.
*
1
0.1
0.06
0.02
-0.023 -0.021 -0.019 £o
Fig. 7. The dependencies of the equilibrium volume fraction * on the external field o
7. Conclusion
We developed a model to describe intermediate layers in composite materials as new phase areas. Phase transformation around spherical isolated inclusion was analyzed in detail. Then the growth of new phase areas around the particles in composite materials was described using the self-consistent approach.
Due to the fact that phase transformations lead to the total energy decrease and stability reasons, simultaneous increasing of all elastic moduli is impossible but it is possible to increase some of the moduli. We demonstrated that if materials of the components are isotropic and inclusions are spherical then the new phase formation under all-round compression can lead to increasing the bulk modulus. In this case the shear modulus is to decrease. We also demonstrated that due to the energy release the phase transformation can act as a reinforcing factor.
We developed a procedure for the interface stability analysis and showed that the desirable process of the new phase intermediate layers formation may be very sensitive with respect to the temperature and loading conditions. This in turn may result in experimental data scattering and misunderstending.
8. Acknowledgement
The financial support of RFBR (Grant No 10-01-00670-a), Programs for Fundamental Research No 22,13, Grant from DFG (AL 341/35-1), the Russian Ministry of Education and Science (Contract 14.740.11.0353) is gratefully acknowledged.
References
[1] Cadek M., Coleman J. N., Barron V., Hedicke K., Blau W. J. Morphological and mechanical properties of carbon-nanotube-reinforced semicrystalline and amorphous polymer composites // Applied Physics Letters, 2002. V. 81. 27. P. 5123-5125.
[2] Freidin A. B. On new phase inclusions in elastic solids // ZAMM, 2007. V. 87. 2. P. 102-116.
[3] Freidin A.B., Y.B. Fu, Sharipova L.L., Vilchevskaya E.N. Spherically symmetric two-phase deformations ans phase transition zones // Int. J. Solids and Struct, 2006. V. 43. P. 4484-4508.
[4] George R., Kashyap K.T., Rahul R., Yamdagni S. Strengthening in carbon nanotube/aluminum (CNT/Al) composites // Scripta Materialia, 2005. V. 53. P. 1159-1163.
[5] Freidin A.B., Vilchevskaya E.N. Multiple development of new phase inclusions in elastic solids // International Journal of Engineering Science, 2009. V. 47. P. 240-260.
[6] Kanaun S.K., Levin V.M. Self-consistent methods for composites, Vol. 1: Static Problems. — Springer, 2007.
1 1 1
-0.023 -0.021 -0.019
[7] Kublanov L. B., Freidin A. B. Solid phase seeds in a deformable material // J. Appl. Math. Mech., 1988. V. 52. P. 382-389.
[8] Kunin L. A. Elastic media with Microstructure. II. Three Dimensional Models. — Springer Series in Solid State Sciences, Berlin, New York, etc. Springer-Verlag, 1983, 272 p.
[9] Morozov N. F., Freidin A. B. Phase transition zones and phase transformations of elastic solids under different stress states // Proc. Steklov Math. Inst., 1998. V. 223. P. 220-232.
[10] Yeremeyev V. A., Freidin A. B., Sharipova L. L. The stability of the equilibrium of two-phase elastic solids // J. Appl. Math. Mech., 2007. V. 71. P. 61-84.
[11] Freidin A. B., Chiskis A. M. Regions of phase transitions in nonlinear-elastic isotropic materials. Part 1: Basic relations // Mechanics of Solids, 1994. V. 29(4). P. 91-109; Mechanics of Solids, 1994. V. 29(5). P. 46-58. Translation from Izv. RAN, Mekhanika Tverdogo Tela, 1994. V. 4. P. 91-109; V. 5. P. 49-61.
[12] Freidin A. B. Small strains approach to the theory on solid-solid phase transformations under the process of deformation // Studies on Elasticity and Plasticity (St. Petersburg State University), 1999. V. 18. P. 266-290.
[13] Eremeev V.A., Freidin A.B., Sharipova L.L. Nonuniqueness and Stability in Problems of Equilibrium of Elastic Two-Phase Bodies // Doklady Physics, 2003. V. 48(7). P. 359363. Translated from Doklady Akademii Nauk, 2003. V. 391. 2. P. 189193.
[14] Morozov N.F., Nazyrov I.R., Freidin A.B. One-dimensional problem on phase transformation of an elastic ball // Dokl. Akad. Nauk, 1996. V. 346. P. 188-191.
[15] Thostenson E.T., Ren Z., Chou T.-W. Advance in the science and technology of carbon nanotubes and their composites: A review // Composites Sciece and Technology, 2001. V. 61. P. 1899-1912.
[16] Xie X.-L., Mai Y.-W., Zhou X.-P. Dispersion and alignment of carbon nanotubes in polymer matrix: A review // Materials Science and Engineering, 2005. V. 49. P. 89-112.
[17] Grabovsky Y., Truskinovsky L. Roughening instability of broken extremals // Arch. Rat. Mech. Anal., 2011. V. 200(1). P. 183202.
[18] Filippov R.A., Freidin A.B., Vilchevskaya E.N. On transition layers in nanocomposites as new phase areas // In: Proc. XXXIX Summer School-Conference Advanced problems in mechanics (APM 2011), 2011, St. Petersburg (Repino), P. 158-167.
