Advances in Building Technologies and Construction Materials 2016
View this Special IssueResearch Article  Open Access
The Static and Dynamic Mechanical Properties of Magnetorheological Silly Putty
Abstract
A novel magnetorheological material defined as magnetorheological Silly Putty (MRSP) is prepared by dispersing soft magnetic particles into Silly Putty matrix with shear stiffening property. Static mechanical properties including creep and stress relaxation and dynamic rheological properties of MRSPs are tested by rheometer. The experimental results indicate that the external magnetic field exerts significant influence on the creep and relaxation behaviors. Moreover, the storage modulus of MRSPs increases sharply in response to the external stimuli of increasing angular frequency automatically and can be enhanced by external magnetic field. Besides, temperature plays a key role in shear stiffening and magnetorheological effect of MRSPs. Furthermore, considering the obstruction to the particle chains formation induced by Silly Putty matrix, a nonperforative particle aggregated chains model is proposed. The model curve is in consistency with experimental data, which means it can describe magnetoinduced behavior of MRSPs well.
1. Introduction
Silly Putty is an intelligent material that exhibits obvious shear stiffening property. It can be easily kneaded, much like dough, into any shape desired. However, if a shock or impulsive load is applied to the putty, it will stiffen and even shatter. In rheological terms the experimenter is distorting the Silly Putty over a range of Deborah numbers. This nondimensional parameter describes the ratio of the fluid relaxation time scale to that of the experimental time scale. A high Deborah number therefore corresponds to a fast experiment in which the load or impulse is applied over a very short time scale. Due to its unique rheological property, Silly Putty is becoming a promising material for a special damper called shock transmission unit (STU) widely used in multispan bridge [1].
Shear stiffening property of Silly Putty has been investigated by many researchers. Cross [2] proposed a threeelement parameter model to describe the compression of Silly Putty at different compression speeds; Goertz et al. [3] developed a mechanical model for Silly Putty based on three Maxwell elements in parallel to study the influence of the temperature on the rheological property and found that the mechanical properties shifted significantly with temperature. Wang et al. [4, 5] prepared a novel shear stiffening gel (STG) with polydimethylsiloxane (PDMS), boric acid, benzoyl peroxide (BPO), and different additives such as calcium carbonate and carbon nanotube. The stiffening mechanism of STG was also explained in respect of chemical bond.
Magnetorheological materials are a class of smart materials that have usually been prepared by dispersing magnetic particles into different polymer matrices. The mechanical properties of magnetorheological materials can be controlled by an externally imposed magnetic field. Due to the controllable rheological property, magnetorheological materials have thus been widely applied in dampers, isolators, and different vibration controllers. However, most solidstate magnetorheological materials cannot present considerable tunability. Interestingly, on the basis of the fabricated STG, Wang et al. [6, 7] firstly developed a magnetically responsive shear stiffening gel with excellent shear stiffening performance and magnetorheological effect. This novel magnetically responsive shear stiffening gel could provide credible tunability with external stimuli.
Besides, some theoretical models have also been developed to describe the magnetoinduced behavior of magnetorheological materials. Jolly et al. [8] proposed a dipole model based on the magnetic interactions between two adjacent particles. Shen et al. [9] presented a mathematical model to represent the stressstrain relationship of magnetorheological elastomer. The model considered all the dipole interactions in a chain and the nonlinear properties of the host composite. Chen et al. [10] developed a finitecolumn model to describe the relationships between the microstructures and the viscoelastic properties. However, for magnetorheological Silly Putty (MRSP) of soft magnetic particles dispersing in Silly Putty matrix, when external magnetic field strength increases, the internal structure changes unlike that of magnetorheological elastomers, and the chain structures are not pulling through like magnetorheological fluids.
In this work, several MRSP samples were prepared by dispersing carbonyl iron particles into Silly Putty matrix. Both static and dynamic mechanical properties were tested by rheometer. The static mechanical testing indicated that MRSPs exhibit obvious creep and stress relaxation behaviors and external magnetic field strength plays a key role in creep compliance and relaxation modulus. The dynamic mechanical testing revealed that MRSPs display both good shear stiffening property and high magnetorheological effect; furthermore, influence factors of mechanical properties for MRSPs including temperature and particle concentration were analyzed. Finally, a nonperforative particle aggregated chains model considering the obstruction to the particle chains formation induced by Silly Putty matrix was proposed to describe magnetoinduced behavior of MRSPs. Due to the tunable shear stiffening and magnetorheological effect, MRSPs are candidates for applications in dampers and STUs widely used in vibration control of buildings and bridges.
2. Experimental Section
2.1. Materials
Dow Corning 3179 Dilatant Compound purchased from Dow Corning Co., Ltd. was used as the Silly Putty matrix. The composition of the silicone material was as outlined in Table 1. The soft magnetic particles carbonyl iron (CI) of average size 3.5 μm was purchased from Jiangsu Tianyi UltraFine Metal Powder Co., Ltd., Xuyi, China.

2.2. Preparing and Testing the MRSPs
The dilatant compound as matrix and different volume fractions of CI particles as fillers were homogeneously mixed by a tworoll mill (Nantong Hailite Rubber Machinery Inc., China, model XK160) at room temperature. For mechanical mixing method, up to six different volume fractions had been considered in this study: 0, 6.98, 10.11, 15.84, 20.80, and 27.29% and the prepared samples were marked as MRSP0, MRSP1, MRSP2, MRSP3, MRSP4, and MRSP5 in sequence. For particle filled materials like MRSPs, dispersion quality is crucial in the rheological properties [11]. Therefore, in this work, Hitachi S4800 scanning electron microscope (SEM) was used to observe the internal microstructure of MRSPs. MRSP sample and the SEM image are displayed in Figure 1. It can be obtained that CI particles are uniformly dispersed in Silly Putty matrix approximately and MRSP sample presents isotropic feature due to great mechanical force provided by milling roller XK160.
Static mechanical properties including creep and stress relaxation and dynamic rheological properties were carried out using a commercial rheometer (Physica MCR 302, Anton Paar Co., Austria). During the testing procedure, a parallel plate PP20 with a diameter of about 20 mm was used and the gap of 1 mm was kept all the time. At the same time, a controllable magnetic field was generated by an external coil. Besides, all the testing samples were kept at approximately the same volume value. Most of all, as a kind of ratedependent material, MRSPs exhibit obvious transition from polymeric fluids on the condition of static state or low shear rate to solidlike polymers under high shear rate. Furthermore, the applied shear strain was very small, especially for dynamic testing, the shear strain was only 0.1%. So, no matter for polymeric fluid state or solidlike state of MRSPs, the fracture or wallslip effects were not considered in this work.
3. Results and Discussion
3.1. Creep and Relaxation Behavior of MRSPs
Figures 2(a) and 2(b) show the creep characteristic of MRSP4 on the condition of different initial stress and zero magnetic field at room temperature. Creep phase (0–300 s) and recovery phase (300–900 s) are explicitly demonstrated in Figure 2(a). MRSP presents more obvious deformation (maximum strain 20.6% of 300 Pa) but lower creep compliance with higher initial stress as shown in Figure 2(b). In the creep phase, creep compliance of MRSP increases with time gradually.
(a)
(b)
Figures 3(a) and 3(b) show the effect of different external magnetic field strength on creep characteristic of MRSP4 with the initial stress 300 Pa. It proves that external magnetic field strength can substantially confine creep deformation and creep compliance of MRSP. With an external magnetic field of 23.9 kA/m, the maximum strain sharply decreases to only 0.375% comparing with the value 20.6% of zero magnetic field. That is because CI particles arrange in chain structure with external magnetic field, thus inducing the generation of magnetoinduced modulus which contributes to the reciprocal of creep compliance .
(a)
(b)
Relaxation modulus of MRSP4 as a function of time under different initial strain and external magnetic field strength is displayed in Figures 4(a) and 4(b), respectively. Without external magnetic field, MRSP exhibits obvious relaxation behavior especially under larger initial strain. Taking initial strain of 3.5% as example, relaxation modulus falls from 507 kPa to 95.8 kPa in less than one second. Besides, as shown in Figure 4(b), external magnetic field can improve relaxation modulus but delay the starting time of relaxation behavior. Interestingly, with the increasing external magnetic field strength, relaxation modulus increases with time firstly and then decreases to generate relaxation. When the external magnetic field strength reaches 156 kA/m, the relaxation would not generate until time of 1 s.
(a)
(b)
3.2. Shear Stiffening Property of MRSPs
The MRSP exhibits apparent temperature dependent property. The response of storage modulus to a change in shear angular frequency for MRSP4 at different temperature is shown in Figure 5. The angular frequency is applied from 1 to 100 rad/s. When the angular frequency is increased, the shear rate on the MRSP also increases. As the temperature decreases, the storage modulus increases more rapidly within a certain angular frequency (0–20 rad/s) that demonstrates excellent shear stiffening property. To quantify the shear stiffening property, absolute shear stiffening effect (ASTe) and relative shear stiffening effect (RSTe) referring to characterization of magnetorheological effect [6] are introduced in (1). is the maximum storage modulus of MRSP4 induced by the shear angular frequency and is the initial storage modulus. The details about MRSP4 at different temperature are listed in Table 2:

MRSP4 displays maximum ASTe 0.805 MPa at 10°C and maximum RSTe 6600% at 60°C as shown in Table 2. Due to the larger initial storage modulus at lower temperature, RSTe of MRSP increases with increasing temperature, but the ASTe is on the contrary.
Figure 6(a) shows the shear stiffening property of different CI contents of MRSPs at temperature 25°C. The storage modulus of all the testing samples increases with increasing angular frequency from 0.1 rad/s to 100 rad/s which indicates that MRSPs become stiffer under application of external shear stress. On the condition of no external magnetic field, CI particles only play a role in reinforcing the Silly Putty matrix and enhance the storage modulus of MRSPs. The storage modulus of the MRSP4 as a function of angular frequency at different magnetic field strength is shown in Figure 6(b). The shear stiffening property can be controlled through the external magnetic field due to the formation of magnetic chains of CI particles in MRSP. All the relevant details to quantify the shear stiffening effect of MRSP from Figures 6(a) and 6(b) are listed in Table 3.

(a)
(b)
Pure Silly Putty MRSP0 displays maximum RSTe and MRSP5 with CI volume fraction of 27.29% exhibits maximum ASTe as shown in Table 3. Due to the reinforcing effect of CI particles, more soft magnetic fillers can induce larger initial storage modulus at 0.1 rad/s which contribute to the descending RSTe. Meanwhile, when the external magnetic field strength is 23.9 kA/m, the ASTe of MRSP4 with CI volume fraction of 20.8% reaches to 1.01 MPa. But as the magnetic field strength is increased, both ASTe and RSTe decrease. Therefore, it can be concluded that the state of larger magnetic field strength is not conducive to the shear stiffening property of MRSP, but lower magnetic field strength can improve ASTe of MRSP.
Overall, MRSP is a dilatant material where the viscosity increases faster than the strain rate. There are two mechanisms (and hence two characteristic time scales) at work in this material [1, 12]. The high molecular weight PDMS has a characteristic polymeric relaxation time . However due to the boric acid there are also transient boron mediated “crosslinks” arising from associating boron linkages. These act to give the MRSP a behavior more like an elastic solid than a liquid. However since these “crosslinks” are dynamic (with a characteristic time , i.e., much shorter than ) the material is not permanently locked in place and can consequently flow under the correct conditions. Therefore at time scales longer than the MRSP behaves like a high molecular weight polymeric fluid (with a characteristic relaxation of ). Over time scales much shorter than MRSP behaves like a crosslinked elastic solid.
The mechanisms of shear stiffening behind MRSPs are very complex, but our goal here is to discover a relatively simple function to accurately fit the storage modulus versus angular frequency at different temperature, CI contents, and external magnetic field strength. GalindoRosales et al. [13, 14] used an apparent viscosity function with 11 parameters to describe the three characteristic regions of shear thickening fluid (STF), including the slight shear thinning at low shear rates, followed by a sharp increase in viscosity over the critical shear rate, and a subsequent, but pronounced, shear thinning region at high shear rates. Furthermore, Tian et al. [15] proposed a modified model shown in (2) to overcome the drawback of 11 parameters function that if any , , or was a decimal, the value in the bracket must be positive; otherwise the equation failed:
A new function is proposed in (3) through imitating part of (2) to describe the shear stiffening property of MRSPs. In (3), represents the storage modulus at initial experimental angular frequency ; represents the maximum storage modulus at final experimental angular frequency ; parameter refers to the slope of the shear stiffening behavior in loglog plot; parameter is responsible for the ratio of changed storage modulus induced by frequency to the initial storage modulus (RSTe); finally, parameter is introduced to overcome the shortcoming that the final experimental point () in curves cannot be fitted:
Figure 7 shows the fitting curves of all the angular frequency scanning testing data on the condition of different temperature, CI contents, and external magnetic field strength. It is obvious that (3) fits to the experimental results perfectly. LevenbergMarquardt algorithm (LMA) is applied in all the fitting procedures. This algorithm combines the GaussNewton and steepest descent methods to obtain the values of the parameters by an iterative chisquare minimisation technique. Fitting is considered to have converged when the difference between the values obtained in two successive iterations is smaller than a given tolerance, chosen here as 10^{−10}.
Table 4 lists all the fitting values of parameters in (3). It is observed from related coefficient that a high fitting precision is obtained in each fitting result. From each section, parameter is inversely proportional to the value of RSTe; interestingly, when there is an external magnetic field applied (Section III), due to higher initial storage modulus and lower changed storage modulus, parameter increases with the value magnetic field strength obviously. Parameter reflects the slope of shear stiffening behavior. In loglog plot, the slope difference in each section is not very distinct which leads to tiny change or even the same values of parameter . Besides, parameter is relatively independent of experimental data that means the value changes a little in each section.

3.3. Magnetorheological Property of MRSPs
MRSP is also a novel magnetorheological material whose mechanical properties are highly dependent on the externally applied magnetic field. Figure 8 shows the response of storage modulus to a change in magnetic flux density for MRSP4 at different temperature. The angular frequency applied on MRSP4 is 10 rad/s and magnetic flux density increases from 0 T to 0.894 T. Similar to its shear stiffening property, magnetorheological effect of MRSP is also influenced by temperature. Meanwhile, absolute magnetorheological effect (AMRe) and relative magnetorheological effect (RMRe) are introduced in (4) to quantify the magnetorheological property. is the maximum storage modulus of MRSP4 induced by the external magnetic field and is the initial storage modulus. The details obtained from Figure 8 are listed in Table 5:

The overall trends displayed in Table 5 are similar to those presented in Table 2. When the temperature is increased, RMRe of MRSP4 enhances significantly, but the AMRe is on the contrary. Higher temperature can soften the Silly Putty matrix, thus inducing the lower initial storage modulus and less resistance to the regular arrangement of CI particles.
Figures 9(a) and 9(b) show the magnetic fielddependent storage modulus of MRSPs with different CI contents at angular frequency 20 rad/s and 100 rad/s, respectively. These experimental curves demonstrate that there is a tendency for the storage modulus of MRSPs to increase with the increase of the magnetic flux density firstly and then saturate with further increases of the magnetic flux density. CI content is observed from Table 6 to have direct influence on AMRe of MRSPs, which means that higher CI content induces greater AMRe. The magnetoinduced storage modulus is mainly attributed to the particle chains formed under external magnetic field. However, greater AMRe is not certain to keep up with higher RMRe. For example, when the applied angular frequency is 20 rad/s, AMRe of MRSP4 is 1.614 MPa less than that of MRSP5, but RMRe of MRSP4 reaches to maximum value of 216.35%. Moreover, it can also be concluded that magnetorheological effect of MRSPs on lower shear rate state is greater than that of higher shear rate state.

(a)
(b)
3.4. Modelling on Magnetoinduced Behavior
Before applying external magnetic field, CI particles are homogeneously dispersed in Silly Putty matrix, but when there is an external magnetic field applied, particles are arranged in chain structure instantaneously as shown in Figure 10. In the chain structure, CI particles are simplified to spherical particles with the radius and the additional acting force between two particles caused by magnetic field leads to the increment of the storage modulus. However, the viscosity of Silly Putty matrix is much larger than general liquids which can cause tremendous obstruction to the formation of pullingthrough particle chains. So a nonperforative particle aggregated chains model is displayed in Figure 10 to explain the magnetoinduced mechanism. In the model, MRSP is divided into several square columns with the length of side and the height . At zero magnetic field strength, particles are randomly dispersed in each square column and MRSP exhibits isotropy. When a low magnetic field is applied, particles aggregate to several apart chains along the direction of magnetic field in each square column with different length and the total length of the chains included in a square column can be marked as . The chains in parallel are sparse due to the low magnetic field strength. However, when the magnetic field strength increases, the aggregated chains in each square column turn unstable and a small disturbance (shearing) may drift some particles in a square column away to aggregate to another apart chain in a new square column. So the chains in parallel become close and even to form particle columns extremely, but the total chains length in a square column becomes short.
For the square columns, the effective permeability can be calculated by using Maxwell Garnett mixing rule as [16, 17]where and are the relative permeability of CI particles and Silly Putty matrix, respectively. is the particles volume fraction in each square column. In the granular basic unit, the volume fraction can be calculated as . So can be expressed as below:where is relevant to external magnetic field strength.
Besides, if is the volume fraction of the square columns in MRSP and can be expressed as where is the particles volume fraction in MRSP.
So, the effective permeability along the direction of square columns can be calculated as follows according to the parallel connection rule [18, 19] and (5)–(7):
When shear strain is applied on MRSP under the external magnetic field, the relative permeability of MRSP changes which can induce the generation of additional force, that is, magnetoinduced shear stress. The magnetoinduced shear stress can be expressed aswhere is the vacuum permeability, is the shear strain, and is the external magnetic field. The shear strain can be marked asThe magnetoinduced shear stress can be obtained by substituting (8) and (10) into (9) as
The magnetoinduced shear modulus can be calculated as
Therefore, the storage modulus can be calculated as below:
Furthermore, the relative permeability of CI particles is the function of external magnetic field strength, which can be obtained from FrohlichKennelly equation [20]:where is the relative permeability of CI particles at zerofield strength and is the magnitude of the saturation magnetization.
Besides, the variable can be seen as the function of external magnetic field strength, and the possible function can be expressed as below, and , , and are the parameters:
From (13) to (15), the vacuum permeability is H/m and the relative permeability of Silly Putty is 1. Besides, for CI particles, the relative permeability at zerofield strength is 70 and the magnitude of the saturation magnetization reaches 1360 kA/m [21]. Furthermore, the applied shear strain is 0.1%. Therefore, according to the nonperforative particle aggregated chains model, taking MRSP2, MRSP3, and MRSP4 for example ( = 20 rad/s), the curves are displayed in Figure 11(a) comparing with the experimental data. The model curves agree well with experimental data which indicates that singlechain aggregation model can describe magnetoinduced mechanical behavior well. From Figure 11(a), the variation trend of as a function of magnetic field strength is shown in Figure 11(b). It is obvious that larger particle volume fraction induces larger value of due to more particles in each square column. However, variable is a decreasing function of magnetic field strength which is in consistency with the description in the front.
(a)
(b)
4. Conclusions
In this work, the multifunctional composite MRSPs with shear stiffening property and magnetorheological property were prepared. Both static mechanical behaviors including creep and relaxation and dynamic rheological behavior were investigated by rheometer. Firstly, MRSPs exhibited obvious creep and relaxation behavior, but a low external magnetic field strength (23.9 kA/m) could substantially confine creep deformation and creep compliance of MRSPs (strain of MRSP4 from 20.6% to 0.375%); moreover, external magnetic field could improve relaxation modulus but delay the starting time of relaxation behavior. Secondly, in dynamic experiments, the storage modulus of MRSPs changed sharply under the stimuli of shear forces with increasing angular frequency, which presented good shear stiffening effect; furthermore, the storage modulus could be enhanced by stimuli of external magnetic field, which exhibited typical magnetorheological effect; besides, both shear stiffening effect and magnetorheological effect were influenced by temperature, CI volume fraction, and external magnetic field strength. A modified parameter model and a nonperforative particle aggregated chains model which agreed well with experimental data were established, respectively, to describe the shear stiffening and the magnetoinduced behaviors perfectly.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by Research and Innovation Program of Jiangsu Province for Graduate Students (Grant no. B14042); National Natural Science Foundation of China (Grant no. 51508237); Natural Science Foundation of Jiangsu Province (Grant no. BK20140560); and Research Foundation for Advanced Talents of Jiangsu University (Grant no. 14JDG161).
References
 J. Liang and X.H. Zhang, “Rheological properties of SP in shock transmission application,” Journal of Materials in Civil Engineering, vol. 27, no. 9, Article ID 04014250, 2015. View at: Publisher Site  Google Scholar
 R. Cross, “Elastic and viscous properties of Silly Putty,” American Journal of Physics, vol. 80, no. 10, pp. 870–875, 2012. View at: Publisher Site  Google Scholar
 M. P. Goertz, X.Y. Zhu, and J. E. Houston, “Temperature dependent relaxation of a ‘solidliquid’,” Journal of Polymer Science Part B: Polymer Physics, vol. 47, no. 13, pp. 1285–1290, 2009. View at: Publisher Site  Google Scholar
 S. Wang, S. Xuan, W. Jiang et al., “Ratedependent and selfhealing conductive shear stiffening nanocomposite: a novel safeguarding material with force sensitivity,” Journal of Materials Chemistry A, vol. 3, no. 39, pp. 19790–19799, 2015. View at: Publisher Site  Google Scholar
 S. Wang, S. Xuan, Y. Wang et al., “Stretchable polyurethane sponge scaffold strengthened shear stiffening polymer and its enhanced safeguarding performance,” ACS Applied Materials & Interfaces, vol. 8, no. 7, pp. 4946–4954, 2016. View at: Publisher Site  Google Scholar
 S. Wang, W. Jiang, W. Jiang et al., “Multifunctional polymer composite with excellent shear stiffening performance and magnetorheological effect,” Journal of Materials Chemistry C, vol. 2, no. 34, pp. 7133–7140, 2014. View at: Publisher Site  Google Scholar
 Y. Wang, S. Wang, C. Xu, S. Xuan, W. Jiang, and X. Gong, “Dynamic behavior of magnetically responsive shearstiffening gel under high strain rate,” Composites Science and Technology, vol. 127, pp. 169–176, 2016. View at: Publisher Site  Google Scholar
 M. R. Jolly, J. D. Carlson, and B. C. Muñoz, “A model of the behaviour of magnetorheological materials,” Smart Materials and Structures, vol. 5, no. 5, pp. 607–614, 1996. View at: Publisher Site  Google Scholar
 Y. Shen, M. F. Golnaraghi, and G. R. Heppler, “Experimental research and modeling of magnetorheological elastomers,” Journal of Intelligent Material Systems and Structures, vol. 15, no. 1, pp. 27–35, 2004. View at: Publisher Site  Google Scholar
 L. Chen, X. L. Gong, and W. H. Li, “Microstructures and viscoelastic properties of anisotropic magnetorheological elastomers,” Smart Materials and Structures, vol. 16, no. 6, pp. 2645–2650, 2007. View at: Publisher Site  Google Scholar
 F. J. GalindoRosales, P. Moldenaers, and J. Vermant, “Assessment of the dispersion quality in polymer nanocomposites by rheological methods,” Macromolecular Materials and Engineering, vol. 296, no. 34, pp. 331–340, 2011. View at: Publisher Site  Google Scholar
 Cambridge Polymer Group, “The cambridge polymer group silly putty ‘Egg’,” 2002, http://www.campoly.com/files/7213/9937/8983/CS_CPG_Silly_Putty_Egg_Rheology.pdf. View at: Google Scholar
 F. J. GalindoRosales, F. J. RubioHernández, and A. Sevilla, “An apparent viscosity function for shear thickening fluids,” Journal of NonNewtonian Fluid Mechanics, vol. 166, no. 56, pp. 321–325, 2011. View at: Publisher Site  Google Scholar
 F. J. GalindoRosales, F. J. RubioHernández, A. Sevilla, and R. H. Ewoldt, “How Dr. Malcom M. Cross may have tackled the development of ‘An apparent viscosity function for shear thickening fluids’,” Journal of NonNewtonian Fluid Mechanics, vol. 166, no. 2324, pp. 1421–1424, 2011. View at: Publisher Site  Google Scholar
 T. Tian, G. Peng, W. Li, J. Ding, and M. Nakano, “Experimental and modelling study of the effect of temperature on shear thickening fluids,” Korea Australia Rheology Journal, vol. 27, no. 1, pp. 17–24, 2015. View at: Publisher Site  Google Scholar
 X. Dong, N. Ma, J. Ou, and M. Qi, “Predicating magnetorheological effect of magnetorheological elastomers under normal pressure,” Journal of Physics: Conference Series, vol. 412, no. 1, Article ID 012035, pp. 759–773, 2013. View at: Publisher Site  Google Scholar
 J. C. M. Garnett, “Colours in metal glasses and in metallic films,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 203, pp. 385–420, 1904. View at: Publisher Site  Google Scholar
 K. K. Kärkkäinen, A. H. Sihvola, and K. I. Nikoskinen, “Effective permittivity of mixtures: numerical validation by the FDTD method,” IEEE Transactions on Geoscience & Remote Sensing, vol. 38, no. 3, pp. 1303–1308, 2000. View at: Publisher Site  Google Scholar
 O. Wierner, Die Theorie des Mischkörpers für das Feld der stationären Strömung. Erste Abhandlung: Die Mittelwertsätze für Kraft, Polarisation und Energie, 1912.
 Y. M. Shkel and D. J. Klingenberg, “Magnetorheology and magnetostriction of isolated chains of nonlinear magnetizable spheres,” Journal of Rheology, vol. 45, no. 2, pp. 351–368, 2001. View at: Publisher Site  Google Scholar
 A. Dorfmann and R. W. Ogden, “Magnetoelastic modelling of elastomers,” European Journal of Mechanics. A. Solids, vol. 22, no. 4, pp. 497–507, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2016 Fei Guo 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.