Next: Intrinsic properties for Up: Exact solutions for Previous: Field fluctuation result


Equal shear modulus

 

A result that can be very useful in testing elastic programs is the equal shear modulus result, true for any microstructure [21,22]. If there are only two phases present, with equal shear moduli G but different bulk moduli K1 and K2, then the effective shear modulus of the entire system is just G, and the effective bulk modulus K of the system is, in 2-D,

img358.gif

and in 3-D,

img359.gif

where c1 and c2 are the volume fractions of the two phases. Note that the only difference between the two equations is the 4/3 factor in front of the shear modulus G. In 2-D, the result is somewhat simpler if expressed in terms of the Young's modulus E and the Poisson's ratio ν, E = c1 E1 + c2 E2, and ν = c1 ν11 + c2 ν2 [22], which actually looks like a parallel result but is valid for general two phase microstructures.

We choose the same cubic image as was used for the small contrast in conductivity case to test the equal shear moduli result. A 22 x 22 x 22 pixel cube (phase 2) is centered in a 30 x 30 x 30 unit cell, so that the volume fractions are c1 = 0.60563 and c2 = 0.39437. The two shear moduli are taken equal to unity, G1 = G2 = 1, and the two bulk moduli are K1 = 1, and K2 = 20. The exact answer, according to eq. (53), is K = 2.263250, while the numerical answer is K = 2.269684, a difference of only 0.3%.

In 2-D, a test of the equal shear moduli result can be combined with a test of the effect of digital resolution by considering the checkerboard problem. The shear moduli are both equal to 2, and there is a ratio of 10 between the bulk moduli (K1 = 1, K2 = 10) (see Fig. 6 for a picture of the microstructure). The effective bulk modulus is computed as a function of system size L x L, where each "check" is L/2 x L/2. The exact value of K is 2.8 using eq. (60). Table 9 displays the data from this test. For L > 16, the error in the effective bulk modulus is about 1% or less.


System
size
K% diff.
25.50096.4
43.27116.8
82.9374.9
162.8411.5
322.8120.4
642.8040.1
1282.80120.04
2562.80030.01
5122.80010.004

Table 9: Effective bulk modulus for checkerboard--equal shear moduli case.

Next: Intrinsic properties for Up: Exact solutions for Previous: Field fluctuation result