# Estimation of Elastic Moduli of Particulate Composites by New Models and Comparison with Moduli Measured by Tension, Dynamic, and Ultrasonic Tests

Authors: G. Bourkas; I. Prassianakis; V. Kytopoulos; E. Sideridis; C. Younis
Journal: Advances in Materials Science and Engineering
 (2010)

---

## Abstract

The elastic constants of particulate composites are evaluated employing a theoretical cube-within-cube formation. Two new models of four and five components, respectively, formed by geometrical combination of three-component models existing in the literature, are used as Representative Volume Elements. Using the governing stress and strain equations of the proposed models, two new equations providing the static elastic and shear moduli of particulate composites are formulated. In order to obtain the dynamic elastic and shear moduli, the correspondence principle was applied successively to components connected in series and/or in parallel. The results estimated by the proposed models were compared with values evaluated from existing formulae in the literature, as well as with values obtained by tensile, dynamic, and ultrasonic experiments in epoxy/iron particulate composites. They were found to be close to values obtained by static and dynamic measurements and enough lower compared with values obtained from ultrasonic experiments. The latter is attributed to the high frequency of ultrasonics. Since measurements from ultrasonic's and from dynamic experiments depend on the frequency, the modulus of elasticity estimated by ultrasonic's is compared with that (storage modulus) estimated by dynamic experiments.


---

## Body

## 1. Introduction

Metal particles added in polymer matrices produce composite of greater density, improved electrical conductivity, better thermal conductivity, and consequently, improved behaviour at high operating temperature, and above all, highly improved mechanical properties. In general rigid fillers increase the elastic and shear moduli, and many theories have been developed to explain this effect.Epoxy resins are the most suitable polymers for composite matrices and extensive research has been carried out on their rheological behavior [1–3] and their mechanical properties [4, 5].A rigorous description of a composite system consisting of a matrix in which filler particles have been dispersed is not an easy task. In fact, a great number of geometrical, topological, and mechanical parameters are necessary, the majority of which varies statistically or is simply unknown.Theoretical treatments usually attempt to exploit as much as possible readily available information, which in most cases consists of the mechanical properties of the matrix and the filler and the volume fraction of the latter, while suitable assumptions cover missing data. The best approximation appears to be the determination of upper and lower bounds for the effective moduli of the composite, based on variational principles of mechanics, developed by Hashin [6].Analytical solutions are valid up to some fairly low filler volume fraction as they have to ignore, for reasons of efficiency, any mechanical interaction between neighbouring inclusions [7]. Referring in particular to the moduli, the existing expressions are resulted from the elasticity theory or express a kind of law of mixtures or are simply an attempt to match theoretical curves to experimental data [8–22]. In most of them a perfect adhesion between matrix and filler was assumed as existing between phases of the composite. In [23] the shear moduli of particulate composites are obtained by means of two cube-within-cube models. In [24] a concept of the interaction between the fillers using different distribution of the inclusions into the volume of the matrix is presented. In [25–28] the effect of the interphase on the values of the elastic modulus of the particulate composites is examined.There is also a sizeable volume of literature dealing with dynamic behaviour of heterogeneous systems, where the dispersed phase is a relatively rigid inclusion [29–32]. In [6] a correspondence principle is developed, by means of which effective complex moduli of viscoelastic composites can be determined on the basis of analytical expressions for effective elastic moduli of composites. It is generally found that the storage modulus increases by increasing frequency, for constant temperature, where for high frequency values tend to the value of static elastic modulus of the particulate composite.In this paper by considering the cube-within-cube formation the elastic constants of particulate-filled composites, by using two models of four and five parts, were evaluated. These models are considered geometrically hybridique in relation with the three-part models existing in the literature [18–20]. The assumptions of [18–23] were used, for the evaluation of the static elastic and shear moduli. According to [6] the correspondence principle has been used for the prediction of the dynamic elastic constants of the particulate composite. Tensile experimental results, dynamic results, and results from ultrasonic tests in epoxy/iron particulate composites were compared to the derived theoretical results. Also theoretical results predicted from existing models in the literature, whose equations are in Appendices A and B, were used for comparison. In order to investigate the effect of frequency, the results obtained from ultrasonic experiments were compared with those obtained from dynamic experiments carried out in epoxy/iron particulate composites, since frequency is a common characteristic.


## 2. Theoretical Considerations

The theoretical analysis will be based on the following assumptions.(i)
Particles are perfectly cubic in shape.(ii)
The distribution of the matrix volume to each individual inclusion is also perfectly cubic in shape. The respective sides of two cubes are parallel.(iii)
There are many filler particles and their distribution is uniform, so that the composite may be regarded as a quasihomogeneous and isotropic material.(iv)
The matrix and the filler are elastic, isotropic, and homogeneous.(v)
The volume fraction of filler is sufficiently small for the interaction between fillers to be neglected.(vi)
The deformation applied to the composite is small enough to maintain linearity of the stress-strain relations.(vii)
There is no transverse variation of the strains in the components which are connected in parallel and have the same length in the load direction.(viii)
The stresses do not vary in the direction of the applied load in the components which are connected in series and have the same cross-sections.


### 2.1. Static Elastic Constants

#### 2.1.1. Elastic Moduli

According to [18–20], the models shown in Figures 1(a) and 1(b), named model 1 and 2 respectively, are three-part composites. From Figure 1 the filler volume fraction is given by(1)υf=α3c3.
For uniaxial load in the figures direction, their elastic moduli are given, respectively, by [18–20](2)Ec(1)=Em(1+(m-1)υf2/31+(m-1)(υf2/3-υf)),(3)Ec(2)=Em(1+υfm/(m-1)-υf1/3),
where m=Ef/Em and Ef, Em, and Ec are the elastic moduli of the filler, matrix and composite, respectively. The indices (1) and (2) refer to model 1 and 2, respectively.The three-part models. (a) Paul model, model 1,   (b) Issay-Cohen model, model 2.
(a)(b)The model presented in Figure2(a) according to [18–20] named model 3, for uniaxial loading in the figure direction, has four components. The components (1) and (2) are in parallel, and both are in series with component (3). The above components are in parallel with component (4). For uniaxial loading in Figure 2(a) direction, from forces equilibrium and elongations equality, we can write the following equations:(4)σc=σ3υf1/3+σ4(1-υf1/3),(5)σ3υf1/3=σ1υf2/3+σ2(υf1/3-υf2/3),(6)ε1=ε2,(7)ε4=εc,(8)εc=ε1υf1/3+ε3(1-υf1/3),
where the indices 1, 2, 3, 4, and c correspond to the parts (1), (2), (3), and (4) and the composite, respectively.The proposed models. (a) Four-part model, model 3,   (b) Five-part model, model 4.
(a)(b)The constitutive equations are given from the Hooke law as follows:(9)σ1=ε1Ef,(10)σj=εjEm,j=2,3,4,(11)σc=εcEc.By combining (4)–(11) the following expression for the elastic modulus of particulate composite is obtained:(12)Ec(3)=Em(1+υf1/(m-1)+υf1/3-υf2/3),
where m=Ef/Em.Considering now the model presented in Figure2(b) called model 4 one can observe that it consists of five components and that the coupling of components (1), (2), (3) and (4) is the same as in the case of model 1. The element consisting of these four components is connected in series with component (5). For uniaxial load along the direction shown in Figure 2(b) the force equilibrium and the strain compatibility give the following equations(13)σc=σ5,σc=σ3υf1/3+σ4(1-υf1/3),σ3υf1/3=σ1υf2/3+σ2(υf1/3-υf2/3),ε1=ε2,2εc=ε4(1+υf1/3)+ε5(1-υf1/3),ε4(1+υf1/3)=2ε1υf1/3+ε3(1-υf1/3).
The constitutive equations are given from Hook’s law (14)σ1=ε1Ef,σj=εjEmj=2,3,4,5,σc=εcEc.For the case of model 4 which is evaluated using the superposition principle, the following constitutive equation is also used:(15)σa=εaEa
where index a corresponds to the composite consisting of components 1, 2, and 3 of model 4 with(16)σa=σ3,εa=ε4.
Solving for Ea and Ec the respective equation of the model one obtains (17)Ec(4)=Em·(1+2υf2/3(mυf1/3+1)/(m-1)-υf2/3),Ec(4)=Em(1+υf1/3+υf2/32/(n-1)+υf1/3-υf2/3),
where m=Ef/Em, and n=Eα/Em.


#### 2.1.2. Shear Moduli

In [23] models 1 and 2 under shear loading are presented in state of deformation. According to [23] models 3 and 4 in state of  deformation under shear loading in the figure directions are presented in Figures 3(a) and 3(b). The couplings remain the same as in the case of uniaxial loading in Figures 2(a) and 2(b). The same remark has been made for models 1 and 2 in [23]. Thus, analogous equations for the shear, as in the case of uniaxial loading are obtained. From this fact for the shear moduli identical equations are founded as in the case of elastic modulus where in the place of Em, Ef, and Ec, Gm, Gf, and Gc, respectively, are placed, where G is the shear modulus. Therefore, it follows  that(18)Gc(1)=Gm(1+(ρ-1)υf2/31+(ρ-1)(υf2/3-υf)),(19)Gc(2)=Gm(1+υfρ/(ρ-1)-υf1/3),(20)Gc(3)=Gm(1+υf1/(ρ-1)+υf1/3-υf2/3),(21)Ga(4)=Gm(1+2υf2/3(ρ·υf1/3+1)/(ρ-1)-υf2/3),(22)Gc(4)=Gm(1+υf1/3+υf2/32/(q-1)+υf1/3-υf2/3),
where ρ=Gf/Gm and q=Ga(4)/Gm.The proposed models in state of shear deformation. (a) model 3, (b) model 4.
(a)(b)For comparison, from the existing equations in the literature (A.1), (A.2), (A.3), (A.4), (A.8), and (A.11) are examined where in the place of Ec, Em, and Ef, G, Gm, and Gf, respectively, are placed.


#### 2.1.3. Poisson Ratios

The Poisson ratios can be evaluated from the following two different modes and then by comparing the theoretical results among each other and with experimental results.(a)
Considering the particulate composite as homogeneous and isotropic material and having evaluated as previously the respective elastic and shear moduli, the Poisson ratio of each model may be found from the relation(23)Gc=Ec2(νc+1).(b)
From the inverse law of mixture, which seems to fit fairly well the experiment, for the case of particulate composites [25], we write
(24)1νc=υfνf+υmνm.


### 2.2. Dynamic Elastic Constants

#### 2.2.1. Dynamic Elastic Moduli

Models 1, 2, 3, and 4 consist of components connected in parallel and/or in series. Consequently the static elastic and shear moduli could also be obtained by successive application of the law of mixtures and the inverse law of mixtures. This procedure is used for the evaluation of the dynamic elastic and shear moduli.Considering a two-component composite, the complex modulusEj* of each component is given by(25)Ej*=Ej′+iEj′′,
where j=1,2, Ej′ is the storage modulus, and Ej′′ is the loss modulus.According to the correspondence principle when these two components are connected in parallel it comes out that(26)E1,2*=E1*υ1+E2*υ2,
where υ1,υ2 are the volume fractions of components (1) and (2), and E1*,E2* and E1,2* are the complex moduli of components (1) and (2) and the two-component composite, respectively, and(27)E1,2*=E1,2′+iE1,2′′.From (25), (26), and (27) it comes out that(28)E1,2′=E1′υ1+E2′υ2,E1,2′′=E1′′υ1+E2′′υ2,
where E1,2′ and E1,2′′ are the storage modulus and the loss modulus of the two-component composites.Using the correspondence principle when the components of a two-component composite are connected in series, one obtains(29)1E1,2*=υ1E1*+υ2E2*.Introducing (25) and (26) into (29) one obtains(30)tanδ1,2=E1′′(E2′2+E2′′2)υ1+E2′′(E1′2+E2′2)υ2E1′(E2′2+E2′′2)υ1+E2′(E1′2+E2′2)υ2,1E1,2′(1+tanδ1,2)=E1′υ1E1′2+E1′′2+E2′υ2E1′2+E1′′2,tanδ1,2=E1,2′′E1,2′
from which one obtains E1,2′,E1,2′′, where δ1,2 is the phase lag between strain and stress.The dynamic modulus of the composite consisting of a   number components, as in the presented models, can now be obtained by a successive application of (28) and (30).


#### 2.2.2. Dynamic Shear Moduli

As in the case of the static elastic moduli, the dynamic shear moduli can be obtained from (28) and (30) by replacing E1′,E1′′,E2′,E2′′,E1,2′ and E1,2′′ by G1′,G1′′,G2′,G2′′,G1,2′, and G1,2′′, respectively.


#### 2.2.3. Dynamic Poisson Ratios

Since in the most polymers the same phase lag is observed in tension and shear, one has(31)Em′′Em′=Gm′′Gm′.Applying the correspondence principle in the equation(32)Gm=Em2(1+νm)
with(33)νm*=νm′-iνm′′,Gm*=Gm′+iGm′′,Em*=Em′+iEm′′,
it comes out that(34)νm′=Em′Gm′+Em′′Gm′′2(Gm′2+Gm′′2)-1,νm′′=0.The dynamic Poisson ratios of the composite are now obtained applying the correspondence principle into equation(35)Gc=Ec2(1+νc)
which gives(36)νc*=Ec*2Ec*-1
from which one obtains(37)νc′=Ec′Gc′+Ec′′Gc′′2(Gc′2+Gc′′2)-1,(38)νc′′=Ec′Gc′′-Ec′′Gc′2(Gc′2+Gc′′2).
When the inverse law of mixture is considered for the Poisson ratio of the composite it comes out that(39)1νc′=υfνf+υmνm′,νc′′=0
which gives(40)Ec′′Ec′=Gc′′Gc′.


### 2.3. Ultrasonic Equipment and Measurement Procedures

Energy pulse propagation through the structure at frequencies above the audible range can be related to the material properties. The velocity propagation can be measured since modulus = density× (velocity)2. However, the main aim of the ultrasonic testing of materials which contain discontinuities is to determine the effects of interaction between sound waves and material properties. The basic parameters required for all ultrasonic measuring methods are sound velocities and sound attenuation through the material in which the sound wave travels. Sound velocities cℓ and ct of the longitudinal and transverse waves, respectively, and the density ρc of the material are used for the evaluation of the elastic modulus Ec, the Poisson ratio νc, and the shear modulus Gc via he following relationships:(41)Ec=(1+νc)(1-2νc)(1-νc)ρc·cℓ2,νc=1/2(cℓ/ct)2-1(cℓ/ct)2-1,Gc=ρcct2.
Figure 4 shows a schematic diagram of the ultrasonic pulse-echo measuring system used. The system consists of a broad band (0.5–15 MHz) ultrasonic pulser-receiver flaw detector (Krautkramer) which can generate and receive electric pulses up to 15MHz. K2G and K2N probes were used as transmitting-receiving transducers of sound waves, producing ultrasounds of 2 and 4MHz, respectively. A simple machine oil was used as the transducer/specimen interface couplant. A contact load for both probes of 9.88N was applied to the transducer/specimen interface.Figure 4
Schematic diagram of the used ultrasonic pulse-echo measuring system.The pulser section produces and injects ultrasonic pulses into the specimen through the transducer, and the reflected signals produced are amplified by the receiver section of the equipment and displayed on the oscilloscope.The sound velocitycℓ* of the longitudinal waves of each specimen was evaluated using the relationship(42)cℓ*=cℓdxdy,
where cℓ is the sound velocity of the reference block, dx is the real specimen thickness, and dg is the equivalent thickness of the specimen, which is measured on the screen of the oscilloscope.


## 2.1. Static Elastic Constants

### 2.1.1. Elastic Moduli

According to [18–20], the models shown in Figures 1(a) and 1(b), named model 1 and 2 respectively, are three-part composites. From Figure 1 the filler volume fraction is given by(1)υf=α3c3.
For uniaxial load in the figures direction, their elastic moduli are given, respectively, by [18–20](2)Ec(1)=Em(1+(m-1)υf2/31+(m-1)(υf2/3-υf)),(3)Ec(2)=Em(1+υfm/(m-1)-υf1/3),
where m=Ef/Em and Ef, Em, and Ec are the elastic moduli of the filler, matrix and composite, respectively. The indices (1) and (2) refer to model 1 and 2, respectively.The three-part models. (a) Paul model, model 1,   (b) Issay-Cohen model, model 2.
(a)(b)The model presented in Figure2(a) according to [18–20] named model 3, for uniaxial loading in the figure direction, has four components. The components (1) and (2) are in parallel, and both are in series with component (3). The above components are in parallel with component (4). For uniaxial loading in Figure 2(a) direction, from forces equilibrium and elongations equality, we can write the following equations:(4)σc=σ3υf1/3+σ4(1-υf1/3),(5)σ3υf1/3=σ1υf2/3+σ2(υf1/3-υf2/3),(6)ε1=ε2,(7)ε4=εc,(8)εc=ε1υf1/3+ε3(1-υf1/3),
where the indices 1, 2, 3, 4, and c correspond to the parts (1), (2), (3), and (4) and the composite, respectively.The proposed models. (a) Four-part model, model 3,   (b) Five-part model, model 4.
(a)(b)The constitutive equations are given from the Hooke law as follows:(9)σ1=ε1Ef,(10)σj=εjEm,j=2,3,4,(11)σc=εcEc.By combining (4)–(11) the following expression for the elastic modulus of particulate composite is obtained:(12)Ec(3)=Em(1+υf1/(m-1)+υf1/3-υf2/3),
where m=Ef/Em.Considering now the model presented in Figure2(b) called model 4 one can observe that it consists of five components and that the coupling of components (1), (2), (3) and (4) is the same as in the case of model 1. The element consisting of these four components is connected in series with component (5). For uniaxial load along the direction shown in Figure 2(b) the force equilibrium and the strain compatibility give the following equations(13)σc=σ5,σc=σ3υf1/3+σ4(1-υf1/3),σ3υf1/3=σ1υf2/3+σ2(υf1/3-υf2/3),ε1=ε2,2εc=ε4(1+υf1/3)+ε5(1-υf1/3),ε4(1+υf1/3)=2ε1υf1/3+ε3(1-υf1/3).
The constitutive equations are given from Hook’s law (14)σ1=ε1Ef,σj=εjEmj=2,3,4,5,σc=εcEc.For the case of model 4 which is evaluated using the superposition principle, the following constitutive equation is also used:(15)σa=εaEa
where index a corresponds to the composite consisting of components 1, 2, and 3 of model 4 with(16)σa=σ3,εa=ε4.
Solving for Ea and Ec the respective equation of the model one obtains (17)Ec(4)=Em·(1+2υf2/3(mυf1/3+1)/(m-1)-υf2/3),Ec(4)=Em(1+υf1/3+υf2/32/(n-1)+υf1/3-υf2/3),
where m=Ef/Em, and n=Eα/Em.


### 2.1.2. Shear Moduli

In [23] models 1 and 2 under shear loading are presented in state of deformation. According to [23] models 3 and 4 in state of  deformation under shear loading in the figure directions are presented in Figures 3(a) and 3(b). The couplings remain the same as in the case of uniaxial loading in Figures 2(a) and 2(b). The same remark has been made for models 1 and 2 in [23]. Thus, analogous equations for the shear, as in the case of uniaxial loading are obtained. From this fact for the shear moduli identical equations are founded as in the case of elastic modulus where in the place of Em, Ef, and Ec, Gm, Gf, and Gc, respectively, are placed, where G is the shear modulus. Therefore, it follows  that(18)Gc(1)=Gm(1+(ρ-1)υf2/31+(ρ-1)(υf2/3-υf)),(19)Gc(2)=Gm(1+υfρ/(ρ-1)-υf1/3),(20)Gc(3)=Gm(1+υf1/(ρ-1)+υf1/3-υf2/3),(21)Ga(4)=Gm(1+2υf2/3(ρ·υf1/3+1)/(ρ-1)-υf2/3),(22)Gc(4)=Gm(1+υf1/3+υf2/32/(q-1)+υf1/3-υf2/3),
where ρ=Gf/Gm and q=Ga(4)/Gm.The proposed models in state of shear deformation. (a) model 3, (b) model 4.
(a)(b)For comparison, from the existing equations in the literature (A.1), (A.2), (A.3), (A.4), (A.8), and (A.11) are examined where in the place of Ec, Em, and Ef, G, Gm, and Gf, respectively, are placed.


### 2.1.3. Poisson Ratios

The Poisson ratios can be evaluated from the following two different modes and then by comparing the theoretical results among each other and with experimental results.(a)
Considering the particulate composite as homogeneous and isotropic material and having evaluated as previously the respective elastic and shear moduli, the Poisson ratio of each model may be found from the relation(23)Gc=Ec2(νc+1).(b)
From the inverse law of mixture, which seems to fit fairly well the experiment, for the case of particulate composites [25], we write
(24)1νc=υfνf+υmνm.


## 2.1.1. Elastic Moduli

According to [18–20], the models shown in Figures 1(a) and 1(b), named model 1 and 2 respectively, are three-part composites. From Figure 1 the filler volume fraction is given by(1)υf=α3c3.
For uniaxial load in the figures direction, their elastic moduli are given, respectively, by [18–20](2)Ec(1)=Em(1+(m-1)υf2/31+(m-1)(υf2/3-υf)),(3)Ec(2)=Em(1+υfm/(m-1)-υf1/3),
where m=Ef/Em and Ef, Em, and Ec are the elastic moduli of the filler, matrix and composite, respectively. The indices (1) and (2) refer to model 1 and 2, respectively.The three-part models. (a) Paul model, model 1,   (b) Issay-Cohen model, model 2.
(a)(b)The model presented in Figure2(a) according to [18–20] named model 3, for uniaxial loading in the figure direction, has four components. The components (1) and (2) are in parallel, and both are in series with component (3). The above components are in parallel with component (4). For uniaxial loading in Figure 2(a) direction, from forces equilibrium and elongations equality, we can write the following equations:(4)σc=σ3υf1/3+σ4(1-υf1/3),(5)σ3υf1/3=σ1υf2/3+σ2(υf1/3-υf2/3),(6)ε1=ε2,(7)ε4=εc,(8)εc=ε1υf1/3+ε3(1-υf1/3),
where the indices 1, 2, 3, 4, and c correspond to the parts (1), (2), (3), and (4) and the composite, respectively.The proposed models. (a) Four-part model, model 3,   (b) Five-part model, model 4.
(a)(b)The constitutive equations are given from the Hooke law as follows:(9)σ1=ε1Ef,(10)σj=εjEm,j=2,3,4,(11)σc=εcEc.By combining (4)–(11) the following expression for the elastic modulus of particulate composite is obtained:(12)Ec(3)=Em(1+υf1/(m-1)+υf1/3-υf2/3),
where m=Ef/Em.Considering now the model presented in Figure2(b) called model 4 one can observe that it consists of five components and that the coupling of components (1), (2), (3) and (4) is the same as in the case of model 1. The element consisting of these four components is connected in series with component (5). For uniaxial load along the direction shown in Figure 2(b) the force equilibrium and the strain compatibility give the following equations(13)σc=σ5,σc=σ3υf1/3+σ4(1-υf1/3),σ3υf1/3=σ1υf2/3+σ2(υf1/3-υf2/3),ε1=ε2,2εc=ε4(1+υf1/3)+ε5(1-υf1/3),ε4(1+υf1/3)=2ε1υf1/3+ε3(1-υf1/3).
The constitutive equations are given from Hook’s law (14)σ1=ε1Ef,σj=εjEmj=2,3,4,5,σc=εcEc.For the case of model 4 which is evaluated using the superposition principle, the following constitutive equation is also used:(15)σa=εaEa
where index a corresponds to the composite consisting of components 1, 2, and 3 of model 4 with(16)σa=σ3,εa=ε4.
Solving for Ea and Ec the respective equation of the model one obtains (17)Ec(4)=Em·(1+2υf2/3(mυf1/3+1)/(m-1)-υf2/3),Ec(4)=Em(1+υf1/3+υf2/32/(n-1)+υf1/3-υf2/3),
where m=Ef/Em, and n=Eα/Em.


## 2.1.2. Shear Moduli

In [23] models 1 and 2 under shear loading are presented in state of deformation. According to [23] models 3 and 4 in state of  deformation under shear loading in the figure directions are presented in Figures 3(a) and 3(b). The couplings remain the same as in the case of uniaxial loading in Figures 2(a) and 2(b). The same remark has been made for models 1 and 2 in [23]. Thus, analogous equations for the shear, as in the case of uniaxial loading are obtained. From this fact for the shear moduli identical equations are founded as in the case of elastic modulus where in the place of Em, Ef, and Ec, Gm, Gf, and Gc, respectively, are placed, where G is the shear modulus. Therefore, it follows  that(18)Gc(1)=Gm(1+(ρ-1)υf2/31+(ρ-1)(υf2/3-υf)),(19)Gc(2)=Gm(1+υfρ/(ρ-1)-υf1/3),(20)Gc(3)=Gm(1+υf1/(ρ-1)+υf1/3-υf2/3),(21)Ga(4)=Gm(1+2υf2/3(ρ·υf1/3+1)/(ρ-1)-υf2/3),(22)Gc(4)=Gm(1+υf1/3+υf2/32/(q-1)+υf1/3-υf2/3),
where ρ=Gf/Gm and q=Ga(4)/Gm.The proposed models in state of shear deformation. (a) model 3, (b) model 4.
(a)(b)For comparison, from the existing equations in the literature (A.1), (A.2), (A.3), (A.4), (A.8), and (A.11) are examined where in the place of Ec, Em, and Ef, G, Gm, and Gf, respectively, are placed.


## 2.1.3. Poisson Ratios

The Poisson ratios can be evaluated from the following two different modes and then by comparing the theoretical results among each other and with experimental results.(a)
Considering the particulate composite as homogeneous and isotropic material and having evaluated as previously the respective elastic and shear moduli, the Poisson ratio of each model may be found from the relation(23)Gc=Ec2(νc+1).(b)
From the inverse law of mixture, which seems to fit fairly well the experiment, for the case of particulate composites [25], we write
(24)1νc=υfνf+υmνm.


## 2.2. Dynamic Elastic Constants

### 2.2.1. Dynamic Elastic Moduli

Models 1, 2, 3, and 4 consist of components connected in parallel and/or in series. Consequently the static elastic and shear moduli could also be obtained by successive application of the law of mixtures and the inverse law of mixtures. This procedure is used for the evaluation of the dynamic elastic and shear moduli.Considering a two-component composite, the complex modulusEj* of each component is given by(25)Ej*=Ej′+iEj′′,
where j=1,2, Ej′ is the storage modulus, and Ej′′ is the loss modulus.According to the correspondence principle when these two components are connected in parallel it comes out that(26)E1,2*=E1*υ1+E2*υ2,
where υ1,υ2 are the volume fractions of components (1) and (2), and E1*,E2* and E1,2* are the complex moduli of components (1) and (2) and the two-component composite, respectively, and(27)E1,2*=E1,2′+iE1,2′′.From (25), (26), and (27) it comes out that(28)E1,2′=E1′υ1+E2′υ2,E1,2′′=E1′′υ1+E2′′υ2,
where E1,2′ and E1,2′′ are the storage modulus and the loss modulus of the two-component composites.Using the correspondence principle when the components of a two-component composite are connected in series, one obtains(29)1E1,2*=υ1E1*+υ2E2*.Introducing (25) and (26) into (29) one obtains(30)tanδ1,2=E1′′(E2′2+E2′′2)υ1+E2′′(E1′2+E2′2)υ2E1′(E2′2+E2′′2)υ1+E2′(E1′2+E2′2)υ2,1E1,2′(1+tanδ1,2)=E1′υ1E1′2+E1′′2+E2′υ2E1′2+E1′′2,tanδ1,2=E1,2′′E1,2′
from which one obtains E1,2′,E1,2′′, where δ1,2 is the phase lag between strain and stress.The dynamic modulus of the composite consisting of a   number components, as in the presented models, can now be obtained by a successive application of (28) and (30).


### 2.2.2. Dynamic Shear Moduli

As in the case of the static elastic moduli, the dynamic shear moduli can be obtained from (28) and (30) by replacing E1′,E1′′,E2′,E2′′,E1,2′ and E1,2′′ by G1′,G1′′,G2′,G2′′,G1,2′, and G1,2′′, respectively.


### 2.2.3. Dynamic Poisson Ratios

Since in the most polymers the same phase lag is observed in tension and shear, one has(31)Em′′Em′=Gm′′Gm′.Applying the correspondence principle in the equation(32)Gm=Em2(1+νm)
with(33)νm*=νm′-iνm′′,Gm*=Gm′+iGm′′,Em*=Em′+iEm′′,
it comes out that(34)νm′=Em′Gm′+Em′′Gm′′2(Gm′2+Gm′′2)-1,νm′′=0.The dynamic Poisson ratios of the composite are now obtained applying the correspondence principle into equation(35)Gc=Ec2(1+νc)
which gives(36)νc*=Ec*2Ec*-1
from which one obtains(37)νc′=Ec′Gc′+Ec′′Gc′′2(Gc′2+Gc′′2)-1,(38)νc′′=Ec′Gc′′-Ec′′Gc′2(Gc′2+Gc′′2).
When the inverse law of mixture is considered for the Poisson ratio of the composite it comes out that(39)1νc′=υfνf+υmνm′,νc′′=0
which gives(40)Ec′′Ec′=Gc′′Gc′.


## 2.2.1. Dynamic Elastic Moduli

Models 1, 2, 3, and 4 consist of components connected in parallel and/or in series. Consequently the static elastic and shear moduli could also be obtained by successive application of the law of mixtures and the inverse law of mixtures. This procedure is used for the evaluation of the dynamic elastic and shear moduli.Considering a two-component composite, the complex modulusEj* of each component is given by(25)Ej*=Ej′+iEj′′,
where j=1,2, Ej′ is the storage modulus, and Ej′′ is the loss modulus.According to the correspondence principle when these two components are connected in parallel it comes out that(26)E1,2*=E1*υ1+E2*υ2,
where υ1,υ2 are the volume fractions of components (1) and (2), and E1*,E2* and E1,2* are the complex moduli of components (1) and (2) and the two-component composite, respectively, and(27)E1,2*=E1,2′+iE1,2′′.From (25), (26), and (27) it comes out that(28)E1,2′=E1′υ1+E2′υ2,E1,2′′=E1′′υ1+E2′′υ2,
where E1,2′ and E1,2′′ are the storage modulus and the loss modulus of the two-component composites.Using the correspondence principle when the components of a two-component composite are connected in series, one obtains(29)1E1,2*=υ1E1*+υ2E2*.Introducing (25) and (26) into (29) one obtains(30)tanδ1,2=E1′′(E2′2+E2′′2)υ1+E2′′(E1′2+E2′2)υ2E1′(E2′2+E2′′2)υ1+E2′(E1′2+E2′2)υ2,1E1,2′(1+tanδ1,2)=E1′υ1E1′2+E1′′2+E2′υ2E1′2+E1′′2,tanδ1,2=E1,2′′E1,2′
from which one obtains E1,2′,E1,2′′, where δ1,2 is the phase lag between strain and stress.The dynamic modulus of the composite consisting of a   number components, as in the presented models, can now be obtained by a successive application of (28) and (30).


## 2.2.2. Dynamic Shear Moduli

As in the case of the static elastic moduli, the dynamic shear moduli can be obtained from (28) and (30) by replacing E1′,E1′′,E2′,E2′′,E1,2′ and E1,2′′ by G1′,G1′′,G2′,G2′′,G1,2′, and G1,2′′, respectively.


## 2.2.3. Dynamic Poisson Ratios

Since in the most polymers the same phase lag is observed in tension and shear, one has(31)Em′′Em′=Gm′′Gm′.Applying the correspondence principle in the equation(32)Gm=Em2(1+νm)
with(33)νm*=νm′-iνm′′,Gm*=Gm′+iGm′′,Em*=Em′+iEm′′,
it comes out that(34)νm′=Em′Gm′+Em′′Gm′′2(Gm′2+Gm′′2)-1,νm′′=0.The dynamic Poisson ratios of the composite are now obtained applying the correspondence principle into equation(35)Gc=Ec2(1+νc)
which gives(36)νc*=Ec*2Ec*-1
from which one obtains(37)νc′=Ec′Gc′+Ec′′Gc′′2(Gc′2+Gc′′2)-1,(38)νc′′=Ec′Gc′′-Ec′′Gc′2(Gc′2+Gc′′2).
When the inverse law of mixture is considered for the Poisson ratio of the composite it comes out that(39)1νc′=υfνf+υmνm′,νc′′=0
which gives(40)Ec′′Ec′=Gc′′Gc′.


## 2.3. Ultrasonic Equipment and Measurement Procedures

Energy pulse propagation through the structure at frequencies above the audible range can be related to the material properties. The velocity propagation can be measured since modulus = density× (velocity)2. However, the main aim of the ultrasonic testing of materials which contain discontinuities is to determine the effects of interaction between sound waves and material properties. The basic parameters required for all ultrasonic measuring methods are sound velocities and sound attenuation through the material in which the sound wave travels. Sound velocities cℓ and ct of the longitudinal and transverse waves, respectively, and the density ρc of the material are used for the evaluation of the elastic modulus Ec, the Poisson ratio νc, and the shear modulus Gc via he following relationships:(41)Ec=(1+νc)(1-2νc)(1-νc)ρc·cℓ2,νc=1/2(cℓ/ct)2-1(cℓ/ct)2-1,Gc=ρcct2.
Figure 4 shows a schematic diagram of the ultrasonic pulse-echo measuring system used. The system consists of a broad band (0.5–15 MHz) ultrasonic pulser-receiver flaw detector (Krautkramer) which can generate and receive electric pulses up to 15MHz. K2G and K2N probes were used as transmitting-receiving transducers of sound waves, producing ultrasounds of 2 and 4MHz, respectively. A simple machine oil was used as the transducer/specimen interface couplant. A contact load for both probes of 9.88N was applied to the transducer/specimen interface.Figure 4
Schematic diagram of the used ultrasonic pulse-echo measuring system.The pulser section produces and injects ultrasonic pulses into the specimen through the transducer, and the reflected signals produced are amplified by the receiver section of the equipment and displayed on the oscilloscope.The sound velocitycℓ* of the longitudinal waves of each specimen was evaluated using the relationship(42)cℓ*=cℓdxdy,
where cℓ is the sound velocity of the reference block, dx is the real specimen thickness, and dg is the equivalent thickness of the specimen, which is measured on the screen of the oscilloscope.


## 3. Material and Experimental Work

### 3.1. Testing Material in Tension and Ultrasonic Measurements

The specimen used consisted of a matrix material, which was a cold setting system based on a diglycidyl ether of bisphenol. A resin having an epoxy equivalent of 185–192, a viscosity of15Nsm-2 at 25°C, and molecular mass between 370 and 384, was cured with 8wt-% triethylenetetramine filled with iron particles of average radius 75μm. The elastic moduli of the matrix and filler were 3, 5, and 210GNm-2, respectively, whereas the Poisson ratios were 0, 36, and 0, 29, respectively.


#### 3.1.1. Tensile Experiments

Dogbone specimens with constant dimensions of measuring area6×3mm and length 45mm were used during the tensile tests which were carried out with an Instron type testing machine at room temperature. The specimens were tested at a rate of extension of 1mmmin-1. Five filler volume fractions υf and five specimens for each volume fraction were used and the values given correspond to their arithmetic mean value. For the obtention of the stress-strain diagrams strain gauges (KYOWA type, gauge factor k=1.99) were located on the specimen to measure the strains.


#### 3.1.2. Ultrasonic Experiments

The NDE technique used in the present work was the ultrasonic pulse-echo technique. When ultrasonic pulses are introduced into a specimen, they reflect on a discontinuity or on the back wall of the specimen. The magnitude of the echo reflections depends on the changes in the impedance across the specimen.To determine the velocities of longitudinal and transverse waves, five specimens from each volume fraction of the composite material were tested ultrasonically at ambient temperature. During each experiment the quantities obtained from the oscilloscope screen were the equivalent thicknessdg of the particle-filled composite and the echo heights. Measurements at three different points in each of the five specimens were carried out. From these quantities and using (42), the velocity cℓ was evaluated. A suitable probe for the longitudinal waves with frequency 4MHz was used. For the evaluation of the velocity ct, a suitable probe for the transverse waves with frequency 2MHz was used. From the analogous (42) this velocity was calculated(43)Ct*=Ctdxdy.


### 3.2. Dynamic Experiments

The material of the matrix was the same as the material used in the tension experiments and ultrasonic measurements. Four filler volume fractions (5, 10, 15, and 20%) were used for the study of the effect of filler content on the dynamic properties.A Dynastat and Dynalyzer apparatus was used for the measurement of the moduliEc′ and Ec′′. This apparatus could apply a sinusoidal load of maximum amplitude 100 N on a specimen 50 mm in length and 3.5 mm wide. The specimen was mounted between a long upper rod connected to a load cell and a short lower rod coupled to a displacement transducer and connected to a motor, which was a coil suspended in the gap of a permanent pole magnet. By passing a servocontrolled current through the coil, the specimen could be subjected to various sinusoidal loads of prescribed amplitude and frequency.By taking into account the rigidity of the load cell and the type and dimensions of the specimen, the storage and loss moduli were calculated. The measurements were performed at frequencies from 0.1 to 100 Hz at ambient temperature (20–22°C).


## 3.1. Testing Material in Tension and Ultrasonic Measurements

The specimen used consisted of a matrix material, which was a cold setting system based on a diglycidyl ether of bisphenol. A resin having an epoxy equivalent of 185–192, a viscosity of15Nsm-2 at 25°C, and molecular mass between 370 and 384, was cured with 8wt-% triethylenetetramine filled with iron particles of average radius 75μm. The elastic moduli of the matrix and filler were 3, 5, and 210GNm-2, respectively, whereas the Poisson ratios were 0, 36, and 0, 29, respectively.


### 3.1.1. Tensile Experiments

Dogbone specimens with constant dimensions of measuring area6×3mm and length 45mm were used during the tensile tests which were carried out with an Instron type testing machine at room temperature. The specimens were tested at a rate of extension of 1mmmin-1. Five filler volume fractions υf and five specimens for each volume fraction were used and the values given correspond to their arithmetic mean value. For the obtention of the stress-strain diagrams strain gauges (KYOWA type, gauge factor k=1.99) were located on the specimen to measure the strains.


### 3.1.2. Ultrasonic Experiments

The NDE technique used in the present work was the ultrasonic pulse-echo technique. When ultrasonic pulses are introduced into a specimen, they reflect on a discontinuity or on the back wall of the specimen. The magnitude of the echo reflections depends on the changes in the impedance across the specimen.To determine the velocities of longitudinal and transverse waves, five specimens from each volume fraction of the composite material were tested ultrasonically at ambient temperature. During each experiment the quantities obtained from the oscilloscope screen were the equivalent thicknessdg of the particle-filled composite and the echo heights. Measurements at three different points in each of the five specimens were carried out. From these quantities and using (42), the velocity cℓ was evaluated. A suitable probe for the longitudinal waves with frequency 4MHz was used. For the evaluation of the velocity ct, a suitable probe for the transverse waves with frequency 2MHz was used. From the analogous (42) this velocity was calculated(43)Ct*=Ctdxdy.


## 3.1.1. Tensile Experiments

Dogbone specimens with constant dimensions of measuring area6×3mm and length 45mm were used during the tensile tests which were carried out with an Instron type testing machine at room temperature. The specimens were tested at a rate of extension of 1mmmin-1. Five filler volume fractions υf and five specimens for each volume fraction were used and the values given correspond to their arithmetic mean value. For the obtention of the stress-strain diagrams strain gauges (KYOWA type, gauge factor k=1.99) were located on the specimen to measure the strains.


## 3.1.2. Ultrasonic Experiments

The NDE technique used in the present work was the ultrasonic pulse-echo technique. When ultrasonic pulses are introduced into a specimen, they reflect on a discontinuity or on the back wall of the specimen. The magnitude of the echo reflections depends on the changes in the impedance across the specimen.To determine the velocities of longitudinal and transverse waves, five specimens from each volume fraction of the composite material were tested ultrasonically at ambient temperature. During each experiment the quantities obtained from the oscilloscope screen were the equivalent thicknessdg of the particle-filled composite and the echo heights. Measurements at three different points in each of the five specimens were carried out. From these quantities and using (42), the velocity cℓ was evaluated. A suitable probe for the longitudinal waves with frequency 4MHz was used. For the evaluation of the velocity ct, a suitable probe for the transverse waves with frequency 2MHz was used. From the analogous (42) this velocity was calculated(43)Ct*=Ctdxdy.


## 3.2. Dynamic Experiments

The material of the matrix was the same as the material used in the tension experiments and ultrasonic measurements. Four filler volume fractions (5, 10, 15, and 20%) were used for the study of the effect of filler content on the dynamic properties.A Dynastat and Dynalyzer apparatus was used for the measurement of the moduliEc′ and Ec′′. This apparatus could apply a sinusoidal load of maximum amplitude 100 N on a specimen 50 mm in length and 3.5 mm wide. The specimen was mounted between a long upper rod connected to a load cell and a short lower rod coupled to a displacement transducer and connected to a motor, which was a coil suspended in the gap of a permanent pole magnet. By passing a servocontrolled current through the coil, the specimen could be subjected to various sinusoidal loads of prescribed amplitude and frequency.By taking into account the rigidity of the load cell and the type and dimensions of the specimen, the storage and loss moduli were calculated. The measurements were performed at frequencies from 0.1 to 100 Hz at ambient temperature (20–22°C).


## 4. Results and Discussion

The 3-component models 1 and 2 were used in [18–20] for the evaluation of the elastic and shear moduli of particulate composites. The proposed models 3 and 4, as it has been mentioned, are formed from geometrical combinations of the above 3-component models. In Figures 5 and 6 the elastic modulus of the composite is plotted versus the filler volume fraction. In these figures one can also observe the curves corresponding to theoretical values obtained from existing equations in the literature, as well as experimental values obtained from tensile experiments and ultrasonic measurements in epoxy/iron particulate composites. In these figures at first one can observe that the predicted values of Ec obtained from models 3 and 4 are bounded from the values predicted from models 1 and 2. In these figures one can also see that the values of Ec predicted from model 4 are above those predicted from model 3. This can be explained by the fact that model 4 consists of more elements from model 1 than model 3 does. In these figures it can also be observed that the predicted values from (12) fit fairly well to the experimental values mainly for high filler volume fraction while the values from (17) are above the experimental values. In Figure 5 it is noticeable that the values predicted by model 3 are close to those obtained by Counto’s (A.11) and that they fit more satisfactorily to the experimental values than Counto’s. In the same figure it can also be observed that the values from model 2 are close to the values predicted from Kerner (A.3) and (A.4) which are used as a lower bound of the elastic modulus of a particulate composite, and that the values of model 2 are lower than those of Kerner. It is also observed in this figure that (A.5) and (A.6) give lower values than those of (3). Another remark from this figure is that the curves (b), (c), (d), and (e) fit well to the experimental values for low filler volume fraction whereas discrepancies appear for high filler volume fraction. In Figure 6 the theoretical values of Ec according to Hashin and Shtrikman bounds are also presented. Comparing the results from Figures 5 and 6 it seems that all the presented curves are bounded by these two curves with little discrepancy in the lower bound. In Figure 6 one can also observe that the values from [16], (A.8), fit to the experimental values for low filler fraction, while discrepancies appear for high filler fraction. This behavior can be explained by the fact that in (A.8) it is assumed that the elastic modulus takes infinite value as υf→1. In the same figure one can also see that Narkis (A.7) gives values close to model 3 for the upper bound of this equation. Another point that must be mentioned by comparing these two figures is that the low bound of Hashin and Shtrikman equations gives values close to those calculated from Kerner equations. Finally from these two figures the high discrepancies between theoretical and experimental static results and ultrasonic measurements become obvious. This fact can be attributed to the high frequency that appears in ultrasonic experiments.Figure 5
Elastic modulus versus filler in iron/epoxy particulate composites volume fraction. (a) Equations (A.5) and (A.6),  (b) equation (A.1),  (c) equation (3),  (d) equation (A.4),  (e) equation (A.3), (f) equation (A.11),  (g) equation (12),  (h) equations (17),  (i) equation (2),  and (j) equation (A.2).Figure 6
Elastic modulus versus filler content in iron/epoxy particulate composites. (a) Equation (3),  (b) equation (A.9),  (c) upper bound, equation (A.7),  (d) equation (12),  (e) equations (17),  (f) equation (2),  (g) equation (A.8), and (h) equation (A.10).The comparison between the theoretical values ofEc, as they result from models 2, 3, and 4, and the experimental values of Ec, can be done taking into account the following facts. (a) In [33] it is mentioned that model 1 corresponds to high adhesion quality between matrix and filler and that model 2 corresponds to low adhesion quality while in [34] it is mentioned that models 3 and 4 correspond to an intermediate adhesion quality. (b) In [29] the concept of the interphase is introduced as a third phase whose physicochemical properties assume intermediate values between the values of the corresponding properties of the filler and those of the matrix. The interphase is a zone located between the filler and the matrix. In [29] the thickness of the interphase is also evaluated. It is observed that for low filler volume fractions the dependence of the thickness of the interphase upon the filler volume fraction is weak while for higher filler volume fractions this dependence is stronger.Based upon the above dependence of the interphase thickness upon the filler volume fraction, the observed disagreement between the theoretical and the experimental values ofEc can be explained as follows. In the presented models 1, 2, 3, and 4 the existence of interphase is not taken into account, but merely a constant adhesion quality is considered for each model, which does not vary with the filler volume fraction. Hence, it is expected that the values of Ec furnished by model 3 are higher than the experimental values for low filler volume fractions and are close to the experimental values for higher filler volume fractions. The same remark holds also for models 2 and 4.The contribution of models 3 and 4 is that they furnish values ofEc lying between the values given by models 1 and 2. These intermediate values can approximate better the experimental values corresponding to higher volume fractions. At the same time the disagreement between the experimental values of Ec and the theoretical values, as given by model 3 for low filler volume fractions, can be considered as acceptably small.It is noticed that the modulus of elasticity given by models 1, 2, 3, and 4 does not depend upon the number of components of each model but upon the geometry of each one. It seems that model 1 and model 2, which are three-component models, give an upper bound and a lower bound, respectively, for the values ofEc that are evaluated by cube-within-cube models. Model 4 furnishes higher modulus of elasticity than model 3 because it consists of more geometric elements taken by model 1, than model 3 does.In Figure7 the shear modulus Gc in epoxy/iron particulate composite is plotted versus the filler volume fraction. The values of Gc predicted from models 1, 2, 3, and 4 as well as those evaluated from existing formulae in the literature are presented. The experimental values have been obtained indirectly from tensile experiments by determining simultaneously the elastic modulus and Poisson ratio, and then the shear moduli were calculated from the well-known relation G=E/2(1+ν), assuming that the material is macroscopically isotropic. In this figure it can also be observed that the same remarks hold as in the case of uniaxial tensile loading.Figure 7
Shear modulus via filler content in iron/epoxy particulate composites. (a) The inverse law of mixtures,  (b) equation (A.1),  (c) equation (19),  (d) equation (A.3), (e) equation (A.4),  (f) equation (A.11),  (g) equation (20),  (h) equations (21) and (22),  (i) equation (18), (j) equation (A.2),  (k) equation (A.8), and (l) law of mixtures.In Figure8 the Poisson ratio values predicted from (23) and from the inverse law of mixtures are plotted versus the filler volume fraction. From this figure it can be seen that the theoretical values predicted by (23) are above those predicted by the inverse law of mixtures, which fit fairly well to the experimental results.Figure 8
Poisson ratio versus filler content in iron/epoxy particulate composites. (a) Equation (24),  (b) equation (23).In Figures9, 10, 11, and 12 the storage modulus Ec′ and the loss modulus Ec′′ are plotted versus the filler volume fraction for the frequencies f=0,1Hz and f=50Hz, respectively. Although the values of Ec predicted by models 3 and 4 approximate better the static experimental results than models 1 and 2, this is not the case for the dynamic experimental results, where the values of Ec′ and Ec′′ predicted by model 2 simulate satisfactorily the experimental results. This can be explained by a possible different behavior of the interphase in the static and in the dynamic loading. Assuming that the elastic modulus Ei of the interphase varies linearly inside the thickness of the interphase, satisfying the inequality Em≤Ei≤Ef, it comes out that the mean value of Ei is Ei¯=(Em+Ef)/2. Then it seems that models 3 and 4 correspond to a constant value of Ei=Ei¯ inside the thickness of the inerphase, by means of which it comes out that ε1=ε2 in models 3 and 4. Inversely, it seems that model 2 corresponds to continuously varying value of Ei inside the thickness of the interphase, by means of which it results that ε1≠ε3 in model 2. From the above remarks one can conclude that probably the first behavior takes place in the static experiments while the second takes place in the dynamic experiments.Figure 9
The storage modulusEc′ versus filler volume fraction for frequency f=0,1Hz in epoxy/iron particulate composite. (a) equation (B.1),  (b) model 2, equations (28) and (30),  (c) eq(B.5),  (d) model 3, equations (28) and (30),  (e) model 4, equations (28) and (30),  (f) model 1, equations (28) and (30),  (g) equation (B.3), and (h) equation (B.7).Figure 10
The loss modulusEc′′ versus filler volume fraction for frequency f=0,1Hz in epoxy/iron particulate composite. (a) equation (B.1),  (b) model 2, equations (28) and (30),  (c) eq(B.5),  (d) model 3, equations (28) and (30),  (e) model 4, equations (28) and (30),  (f) model 1, equations (28) and (30),  (g) equation (B.3), and (h) equation (B.7).Figure 11
The storage modulusE'c versus filler volume fraction for frequency f=50Hz in epoxy/iron particulate composite. (a) equation (B.1),  (b) model 2, equations (28) and (30),  (c) eq(B.5),  (d) model 3, equations (28) and (30),  (e) model 4, equations (28) and (30),  (f) model 1, equations (28) and (30),  (g) equation (B.3),  (h) equation (B.7).Figure 12
The loss modulusEc′′ versus filler volume fraction for frequency f=50Hz in epoxy/iron particulate composite. (a) Equation (B.1),  (b) model 2, equations (28) and (30),  (c) eq(B.5),  (d) model 3, equations (28) and (30),  (e) model 4, equations (28) and (30),  (f) model 1, equations (28) and (30),  (g) equation (B.3), and (h) equation (B.7).Whenm=Ef/Em≫1, the dynamic Poisson ratio predicted by (35) leads to (40). One can easily see this, because in (12) and (17), m/(m-1)≈1, m/(m-1)≈0, and the reinforcing coefficients Ec/Em and Gc/Gm are identical. In this case, in epoxy/particle systems, it has been found that the inverse law of mixture simulates satisfactorily the experimental results [25]. For lower values of m, deviations of νc from the inverse law of mixtures could lead the second member of (38) to be different from zero.In Figure13 the storage modulus Ec′ is plotted versus the filler volume fraction. By extending the experimental values of dynamic measurements of the storage modulus Ec′ in epoxy/iron particulate composites, one can reach the values predicted by ultrasonic measurements. This can be explained by the fact that the classification of the local oscillations to oscillations owed to covalent bonds between the atoms of the main chains or to intermolecular bonds probably is not the same for the low and for the high frequencies.Figure 13
Variation of dynamic storage modulus of compositeEc′ versus frequency for different values of filler volume fraction.


## 5. Conclusions

The values of the elastic and the shear moduli predicted by the new models 3 and 4 are bounded by the values predicted by models 1 and 2. The values given by model 3 are close to the experimental values mainly for high filler volume fractions, as well as to the values predicted by Counto’s equation. The values ofEc predicted by this model almost coincide with the average values of Ec given by models 1 and 2. Similarly the values of Ec given by model 4 fit satisfactorily to the experimental results. Since model 4 is composed by more geometrical elements resulting from model 1, which is the stiffer model in cube-within-cube formation, it is stiffer than model 3. The values of Poisson ratio resulting from (23) of models (3) and (4) are close to and lower than the values of the Poisson ratio of the matrix. These values have been found to be higher than those predicted by the inverse law of mixtures, which fit fairly well the experimental values.The dynamic experimental results are approximated satisfactorily by the values of dynamic moduli predicted by model 2. Probably this is due to different behavior of the interphase material in static tension than in dynamic loading.The values of the elastic modulus measured by ultrasonic measurements are higher than those measured by tensile and dynamic experiments. This can be explained by the fact that the number of the local oscillations owed to covalent bonds between the atoms of the main chains is rather greater in higher than in lower frequencies, resulting to an increase of the elastic modulus in high frequencies. Equivalently the number of the local oscillations owed to intermolecular or other bonds weaker than the covalent bonds seems to be lower in higher than in lower frequencies.



---

*Source: 891824-2010-07-12.xml*