`Journal of NanomaterialsVolume 2012 (2012), Article ID 860371, 7 pageshttp://dx.doi.org/10.1155/2012/860371`
Research Article

## PP/MMT Nanocomposite: Mathematic Modelling of Layered Nanofiller

1Department of Polymer Engineering, Faculty of Technology, Tomas Bata University in Zlin, Nam. T. G. Masaryka 275, 762 72 Zlin, Czech Republic
2Centre of Polymer Systems, University Institute, Tomas Bata University in Zlin, Nad Ovcirnou 3685, 760 01 Zlin, Czech Republic

Received 13 December 2011; Revised 19 January 2012; Accepted 23 January 2012

Academic Editor: Sevan P. Davtyan

Copyright © 2012 Dagmar Merinska et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The comparison of calculated data from proposed mathematic model and experimentally obtained data of PP/clay nanocomposites was done with the focus on the layered shape of MMT platelets. Based on the well-known Kerner's model and the Halpin-Tsai' equation with the use of some described presumption, the mathematic model for PP/clay nanocomposite was proposed. Data from the measurement of prepared PP/clay samples were taken and compared with the calculated ones from the proposed model. The good agreement was found.

#### 1. Introduction

Firstly, the mathematical description of mechanical and physical properties of composites based on a polymer matrix has been studied in 1950s [1, 2]. These and also the next studies are deduced only theoretically [37].

These models mainly describe the dependence of shear modulus of elasticity on the content of the filler and spherical or fibre type of filler is considered. They are named self-consistent, it means that homogenous spherical particles are placed in the homogenous continuous polymer matrix. With the increasing use of composites as construction materials in the industry models with cylindrical dispersed phase has been proposed. This was used in order to simulate short or long fibres.

The classical Kerner’s model is usually presented as:where is a bulk modulus, is shear modulus, is volume fraction, is Poson ratio index belongs to matrix, and index belongs to filler (fibre).

The equation of Halpin-Tsai is very often cited and properly it is the conversion of Hill’s equation [5, 6] for long and short fibres: where (ratio of the volume of filler (fibre) to the whole volume), is a parameter that depends on the matrix Poisson ratio and on the particular elastic property, represents the composite moduli, and, again, indexes and belong to matrix and fibre fraction.

In the literature one can find a lot of attempts of simplification and specification in earlier models [810]. The problem of these models is the fact that they are based on the simplified assumptions that do not correspond with the real states. For example, the rigidity of prepared composites does not only depend on the volume of filler and the geometry of a stress, but also on the particle size, its shape, and the tendency to agglomerate [1113]. Moreover, in case of aspherical particles it depends also on the particle orientation, especially in the case of layered clay nanoparticles where they are supposed to have the shape of platelets with significant aspect ratio.

Nanocomposite materials based on the polymer matrix have been studied more than twenty years. The main reason is the unusual combination of properties. They are unique in comparison with conventional composites because of their large interphase surface (connecting with aspect ratio) and very small distances between reinforcing particles [13]. In this way a very rigid net is created.

The behaviour of polymer chains in the neighbourhood of nanoparticles and the influence of the nanoparticles shape on the nanocomposite properties have been studied in the several next papers [1416].

#### 2. Theoretical

The proper model derivation (modelling).

At the majority of above-mentioned model is based on the mixing principle when the slope of the curve increases with the increasing filler content. Slope of curve represents the elasticity modulus—and it is the function of the filler content. This is valid till the composite modulus reaches for of the value of the filler elasticity ( is a volume percent of filler particles, see later). This behaviour is valid only in case of the composite filled with the long fibres oriented in the longitudinal direction. In all other cases at the condition , (no matrix) the composite is fully non-cohesive small pile of a filler with zero tensile modulus.

More complicated models have a disadvantage—the dependence on stress or deformational fields. In the case of our studied system PP/MMT the tensile modulus initially grows very quickly, but with the increasing filler content the slope of curve falls down. Figure 1 shows the comparison of experimental data with the Kerner’s equation. Values used for this calculation form (2) were chosen only approximately based on the authors ( 0,004, 0,008, 0,016, 0,04) in order to show how this model is far from the real data (taken from Table 3) in case of layered mineral filler in polymer matrix.

Figure 1: The dependence of reduced elasticity modulus on the filler content.

Based on these observations mentioned above we tried to deduce such an equation which at the use of the easiest mathematical apparatus would correspond with obtained experimental data and would be suitable for prognosis of future research in the use of layered nanoclays. The interval of filling was chosen between 0–10 wt. %.

The outgoing assumptions are as follows:(1)the filler has an unlimited modulus → it has zero deformability;(2)the primary deformation is observed in case of torque and strain modulus, after the created stress is measured;(3)particles are fixed in matrix with the perfect adhesion (index ), either free in the vacancy with the zero adhesion (index ), nothing else, thus their volume percent .

When the particles are freely inserted into the vacancies in the matrix, without any adhesion, the system behaves like an expanded material. With the growing number of ratio of particles the stress at the constant deformation decreases according to the equation: On the other hand, at the deformation of the system with the rigidly fixed components the deformation is higher by a part of undeformed particles and by this also the stress is higher: However, in this case of rigid fixation of both components the polymer matrix is not deformed equally. The deformation of the polymer layer at the rigid particle is almost zero. With the increasing distance from this particle the deformation increases and with the approximation to the next particle it goes to the zero again. By this the deformation gradient is created. It is similar to the speed gradient at the liquid flow along fixed plates (see Figure 2).

Figure 2: The simple scheme of the deformation gradient at the liquid flow along fixed plates.

A complete spectrum of deformation matrix from zero at the filler particle to the maximum somewhere along particles is simplified into two values. After the layer at the thickness in the nearest neighbourhood has the zero deformation and the more remote particles have the mentioned maximal value of deformation. However, the total average value of deformation must stay the same even after applying these assumptions.

The consequence of this thought is the increasing of the origin undeformable volume by the undeformable layer of the matrix. This means, for example, in case of cubic filler particles with the edge a the increasing by 6 ha2 at one particle with the volume a3. For the total volume of particles it is the value of 6 , that is, 6 . By this train of thought (5) changes into: The used constant “6” is valid only for the particles of cubic shape. It will be different for different shape of the filler particles. Moreover, this constant depends on the shape of the deformation field (strength, shear), when the rigid shell (standing only temporary during the composite deformation) is not spread evenly on the whole particle surface, thus, after simplifying, .

The next consequence based on the thought discussed above, is that the created stress and by this the rigidity of the composite depends on the size of filler particles.

At the observations it was found that the size of filler particles a depends on . The content of aggregates increases with the growing of . The reason is either imperfect mixing or after aggregation. Thus, where is the size of particle for convergent to zero

The constant will be discussed later.

After the summarization of the above-mentioned thoughts (6) changes into where .

With the growing content of filler at the same initial deformation the stress inside composite is created and increases according to (6). This causes that proportionally to this stress the part of filler particles is breaking loose (is delaminated) from the surface of the polymer matrix. Consequently, its influence on the stress inside the composite changes from the state (5) to the state (4). Thus, the is proportional to the filler content and to the existing stress in the matrix. This stress we will mark as . After where is a critical stress, at which all the filler is delaminated, thus .

Now we have all mathematical apparatus needed for the calculation of the module of elasticity:firstly, we calculate the initial stress caused applied deformation at zero delamination (according to (8)): from the previous equation we can get the ratio of no delaminated particles is (according to (9)):next we calculate the stress equal to free particles (according to (4)):after we calculate the rest of fixed particles:next the stress equal to bonded particles (according to (8)):and the resulting stress:if the deformation at measurement is unit and it is small enough (in order to be responding to Hook’s principle), after

This system of equations contains only three unknown parameters of physical merit which it is necessary to estimate for the next calculation:

—critical stress (it has connection to an adhesion degree of filler particles to the matrix);

—proportional thickness of undeformable layer of matrix in the neighbourhood of filler particle, it is also connected with the geometry of deformational field;

—the velocity of the growing of the particle aggregation degree with regard to , it decreases with the rising quality of dispergation.

In order to obtain the real result of this calculation, the particular parameters must be chosen in the values fulfilling following conditions:

From the influence of the mentioned three parameters it is possible to see from the following three-dimensional graphs (Figures 3, 4, and 5) the following.

Figure 3: The simple scheme of the deformation gradient at the liquid flow along fixed plates—with the filler in polymer melt.
Figure 4: Reduced modulus of elasticity as a function of parameters and .
Figure 5: Reduced modulus of elasticity as a function of parameters and .

The low value of the adhesion between the filler and the matrix (resp. high value of ) causes more significant validity of (4) at the cost of (5) and the rigidity of the sample after the initial increase starts to decrease. This is visible also in the next pictures.

In case of the microscopy observations of nanocomposite samples it was found that with the growing content of nanofiller the number of agglomerates increases; thus, the average size of particles increases too. The velocity of this increase is included in the constant (see above). The ideal requirement is to reach such a dispersion of filler nanoparticle platelets in the polymer matrix that the size of particles should not grow with the , the rising of the rigidity should be maximal, thus, , respectively, (Figure 6).

Figure 6: Reduced modulus of elasticity as a function of parameters and .

Parameter (Figure 5) represents the static part of the matrix in the neighbourhood of the undeformable particle relating to the size of the primary filler particle. After, the rest of the polymer matrix transfers higher stress and consequently the rigidity of the composite goes up. The value of this stress depends on the shape, orientation, and size of the surface of the particles dispersed during the mixing in the kneader.

Time and conditions of the kneading, the technology of the nanocomposite preparation are projected in the parameters and (dispersion of agglomerates, the breaking of lamellar particles bringing the change of the size and shape). depends also on the method used for the measurement of the modulus of elasticity. Parameter is connected with material.

It is necessary to underline that the presented model (like all models) is valid for the limited interval of layered nanocomposite filling only.

#### 3. Experimental

In order to verify the validity of the proposed model for the layered type of nanofiller nanocomposite samples of PP/modified MMT were prepared.

##### 3.1. Materials

Polypropylene (PP) Mosten GB 003 produced by Chemopetrol Litvínov, Czech Republic, was used as a polymer matrix. The density of chosen PP was 907 kg/m3 and melt flow index (MFI) was 3.2 g/10 min at 230°C.

The density of maleated polypropylene (PP-Ma) EXXELOR PO 1015 was 900 kg/m3, MFI was 22, and contents of maleic groups were 1 wt. % (used as a compatibilizator). Maleated polypropylene was supplied by ExxonMobil Chemical Europe, Belgium. Two types of nanofiller Dellite were used (Dellite 72T and Dellite 67G); their concentration was 2, 4, 6, and 10 wt. %. The nanofiller were supplied by Laviosa Chemical Mineraria S.p.A., Italy.

##### 3.2. Preparation of the Polymer-Clay Nanocomposites

All nanocomposites in this work were prepared by melt blending on a Brabender Plasticorder compounder at 40, 60, and 80 rpm and 220°C for 10, 20, 30, and 40 min. The samples were prepared by pressuring at 220°C at 9 min and cooling was 7 min.

##### 3.3. DMA Analysis

The measurement of E modulus was done by DMA analysis on the equipment DMA DX04T (company RMI) on FT TBU in Zlin. Values presented in this work are values at 30°C.

#### 4. Results and Discussion

Resulting data of measurement of dynamic and strength modulus of elasticity for both used nanofiller types are listed in the Tables 1 and 2.

Table 1: Dynamic modulus as a function of speed of rotation of kneader and volume ratio of nanofiller D72T.
Table 2: Summary of results of dynamic modulus and strength modulus for both nanofiller types (data calculated by the same way as in the case of Table 1).
Table 3: The summary of results from Tables 1 and 2.

Table 3 summarizes data from Tables 1–4 averaged through the speed of rotation, estimated values of parameter, and with them corresponding calculated values of elasticity modulus. In order to give the better approach the results are listed in Figure 6.

Based on Figure 6 it is possible to say that the calculated data form the proposed model correspond very well with the experimental values. It is valid especially in case of the dynamic modulus data. This is supported by the 2D graph shown in Figure 7 where the comparison of measured (D in the name) and calculated (M in the name) values for dynamic modulus is presented (data are summarised in Table 3). Modelled data were calculated with the calculated values of parameters , , and . It is possible to see that real data correspond very well with the calculated ones how it was already said. In case of dynamic modulus the data meet the calculated curve almost in 100%.

Figure 7: Reduced elasticity modulus as a function of nanofiller content of prepared samples.

The growing of the particle agglomeration with their rising content makes it possible to observe qualitatively in TEM pictures. The lowering of the parameter is crucial for the next improvement of the reinforcing effect of nanofiller.

#### 5. Conclusion

Models and mathematical descriptions of the various composite behaviours have been created for more than 50 years. Nanofillers especially due to their relatively short time of application and the specific shape and size, unique interface surface, small distance of particles, and the tendency to the agglomeration do not fit to these models.

In this work the authors tried to create a relatively simple model which with the variation of three parameters allows very nice well description of the dependency of the elasticity modulus on the filler content for the layered nanofiller type with the platelet shape of particles. This model is valid in the interval of the clay nanofiller content at least to the 10 wt. %.

Model brings the better look at the material containing particles and technological parameters influenced by the rigidity of prepared composites and the more optimal planning of the next experimental work in this field.

#### Acknowledgments

This project was supported by the Academy of Sciences of the Czech Republic (Project KAN 100400701), Board of Trade in the frame of Project FI-IM3/085 and Internal Grant Agency IGA/FT/2012/040, and by Operational Programme Research and Development for Innovations cofunded by the European Regional Development Fund (ERDF) and national budget of the Czech Republic within the framework of the Centre of Polymer Systéme Project (reg. number: CZ.1.05/2.1.00/03.0111). Some of these data were presented in ANTEC 2009.

#### References

1. E. H. Kerner, “The elastic and thermo-elastic properties of composite media radioradio,” Proceedings of the Physical Society B, vol. 69, pp. 808–813, 1956.
2. J. D. Eshelby, “The determination of the elastic field of an ellipsoidal inclusion, and related problems,” Proceedings of the Royal Society of London A, vol. 241, pp. 376–396, 1957.
3. Z. Hashin, “The elastic moduli of heterogeneous materials,” Transactions of the ASME, Journal of Applied Mechanics, vol. 29, pp. 143–150, 1962.
4. Z. Hashin and S. Shtrikman, “A variational approach to the theory of the elastic behaviour of multiphase materials,” Journal of the Mechanics and Physics of Solids, vol. 11, no. 2, pp. 127–140, 1963.
5. R. Hill, “Elasticproperties of reinforced solids: some theoreticalprinciples,” Journal of the Mechanics and Physics of Solids, vol. 11, no. 5, pp. 357–372, 1963.
6. R. Hill, “Aself-consistentmechanics of composite materials,” Journal of the Mechanics and Physics of Solids, vol. 13, no. 4, pp. 213–222, 1965.
7. B. Budiansky, “On the elastic moduli of some heterogeneous materials,” Journal of the Mechanics and Physics of Solids, vol. 13, no. 4, pp. 223–227, 1965.
8. L. E. Nielsen, “Generalized equation for the elastic moduli of composite materials,” Journal of Applied Physics, vol. 41, no. 11, pp. 4626–4627, 1970.
9. T. B. Lewis and L. E. Nielsen, “Dynamic mechanical properties of particulate-filled composites,” Journal of Applied Polymer Science, vol. 14, no. 6, pp. 1449–1471, 1970.
10. R. M. Christensen, “Solutions for effective shear properties in three phase sphere and cylinder models,” Journal of the Mechanics and Physics of Solids, vol. 27, no. 4, pp. 315–330, 1979.
11. J. Jancar, J. Kucera, and P. Vesely, “Peculiarities of mechanical response of heavily filled polypropylene composites—part I elastic modulus,” Journal of Materials Science, vol. 26, no. 18, pp. 4878–4882, 1991.
12. A. Agbossou, A. Bergeret, K. Benzarti, and N. Alberola, “Modelling of the viscoelastic behaviour of amorphous thermoplastic/glass beads composites based on the evaluation of the complex poisson's ratio of the polymer matrix,” Journal of Materials Science, vol. 28, no. 7, pp. 1963–1972, 1993.
13. Y. Bréchet, J. Y. Cavaillé, E. Chabert, et al., “Polymer based nanocomposites: effect of filler-filler and filler-matrix interactions,” Advanced Engineering Materials, vol. 3, no. 8, pp. 571–577, 2001.
14. J. Kalfus and J. Jancar, “Relaxation processes in pvac-ha nanocomposites,” Journal of Polymer Science B, vol. 45, no. 11, pp. 1380–1388, 2007.
15. J. Kalfus and J. Jancar, “Elastic response of nanocomposite poly(vinylacetate)-hydroxyapatite with varying particle shape,” Polymer Composites, vol. 28, no. 3, pp. 365–371, 2007.
16. J. Kalfus and J. Jancar, “Viscoelastic response of nanocomposite poly(vinyl acetate)-hydroxyapatite with varying particle shape—dynamic strain softening and modulus recovery,” Polymer Composites, vol. 28, no. 6, pp. 743–747, 2007.