Macroscopic Young ’ s Elastic Modulus Model of Particle Packing Rock Layers

Based on the Hertzian granular contact mechanics model, the paper built up a Macroscopic Young’s Elastic Modulus of particle/granular packing rock layers, and built up a ties to connecting Young’s Elastic Modulus of sand particle in Meso and the Macroscopic Young’s Modulus of granular packing rock layers. The Macroscopic Young’s Modulus of granular packing rock layers is far less than the Young’s Modulus of sand particle. The Macroscopic Young’s Modulus of granular packing rock layers is proportioned to the powers of 1/3 of the vertical contact force of sand particles. The Macroscopic Young’s Modulus is inversely proportional to particle diameter. The paper calculated the vertical contact force of five types aligning mode of the particles. When equal stress, the increased of the coordination number lead to the decrease of the contact force fn, this lead to the coordination number is an inverse proportion to Macroscopic Young’s Modulus. But the larger coordination number change only means very little Macroscopic Young’s Modulus change.


Introduction
Young's Modulus of rocks or soil is the important parameters in calculating seismic velocity, liquefaction of sand-soil, seismic exploration and hydrocarbon prediction, earthquake engineering, tunnel excavating, deformation of sand body, dynamic diagenesis of sandstone.Scholars always study the macroscopic structure of porous media used the relation of wave velocity and density of material experientially.But there was a positive correlation between wave velocity and Young's Modulus theoretically and this logical accords with the principle of wave propagation.The difference between experiential methods and theoretical methods make the engineering application method and entirely pure research model each go its own way.Young's Modulus obtained by experiment mostly.But few studies focused on Young's Modulus of granular mixtures or granular packing rock layers.Liu Xu (2002) [1] used multi-phase medium micromechanical model to educe out some macroscopic elastic parameter.But Liu's method was based on most hypothesis on factor of porosity [2].Zhong Xiao-xiong (1992) [3] set up the relationship between fabric tensor and contact density distribution functions, and the relationship between fabric tensor and stress tensor are analyzed.Those model on stress tensor and arrangements of particle is too fussy to used on geologic and macroscopic mechanics of rocks and soils.
The Macroscopic Young's Modulus of dry particle or granular packing rock layers is an important method to distinguish elastic deformation and plastic deformation of sedimentary rock layer."From a grain sand can we find a world?".The paper effort to set up a ties model between the macroscopic elastic parameter of granular packing rock layers and the elastic parameter of those sand particles.And effort to make out a series of simple theoretical logics and catch hold of those dominant factors of those logics.Maybe this can help those investigators and engineers to qualitatively and semi-quantitatively grasp the macroscopic elastic parameter of unconsolidated sandstone.

Model of Elastic Contact
Based on the mechanics of spherical particles contact, classical Hertzian contact theory (Johnson, 1985 [4]; Oda, 1999 [5]; Sun Qi-cheng, 2009) [6] consider that the contact points transformed into interface after the elastic deformation.As Figure 1(a) shows.
The Hertzian contact theory as formula 1 to 7 shows.Among Figure 1, a is the round radius of the contact interface and its distribution as


is the relative approach of the centroids of the two sphere in contact.
Among them, 1 and 2 is the material's Young's Modulus of spherical particles.and is the radius of particles 1 and particles 2, 1  and 2  is the Max normal deformation amount of particles 1 and 2. So, the total normal force is defined as formula 5.
  The inverse ratios of Normal Contact Stiffness n K can express the normal direction deformation amount between two particles center by unit normal direction force.
Sun Qi-cheng [6] thought that if the normal deformation amount  is little, according to as formula 9, the force may be calculated by Hooke's Law.And some studies looks the E  in formula 6 as the Shearing Young's Modulus also mean a great errors.

Macroscopic Young's Modulus
The Macroscopic Young's Modulus is the Young's Modulus of particle packing rock layers.It is smaller than the material's Young's Modulus of spherical particles.
Because of the different of material's and rock layer's, the wave velocity of various depth or petrofabric of rock layers is different.The percentage of sedimentary rock of the rock cover the earth surface was 75% and formed in sand particle packing layer.The compaction and diagenesis in particle packing layer is also the process of pore evolution and change of Macroscopic Young's Modulus.Liu Yu (2010) [7] and Xia Tang-dai (2011) [2] put forward the concept of "effective shear modulus" and expression the relation between the material's Young's Modulus of spherical particles and the Macroscopic Young's Modulus of particle packing layers.The "effective shear modulus" has expression the macroscopic modulus partially but not considered the affect of diameters of the spherical particles in those papers.When we study the sand and soil then looks them as granular mixtures, the "force chain" often as a important mechanism.The affect of the "force chain" is used in the lower stress and we should neglects it in the higher tension.We should attend to that the E  in formula 6 is only a process parameters and a constant, not the real Young's Modulus of particle packing rock layers and no some good application meaning.In fact the Normal Contact Stiffness K n is proportional to Macroscopic Young's Modulus g E when only single contact point.The formula 5 be transformed into formula 7.
When we know the stress of external force which the particle packing rock layers suffered, calculating the derivative of the formula 7, My paper get the formula 8.
The Macroscopic Young's Modulus g E as formula 9 shows.

 
The formula 9 is the theoretical Macroscopic Young's Modulus for arbitrary spherical particles contact for small deformation respectively.If that is equal diameter, equal Young's Modulus, equal Poisson's Ratio spherical particles, formula 9 can be transformed into formula 10.
From those component in formula 10, we can find that the theoretical Macroscopic Young's Modulus of particle packing rock layers mainly be effected by following parameters: 1) The material's Young's Modulus of spherical particles.The theoretical Macroscopic Young's Modulus of particle packing rock layers is proportional to the material's Young's Modulus of spherical particles and inversely proportional to the material's Poisson's Ratio of spherical particles.
2) g E is proportional to 1/3 n f , so proportional to the depth and stress and can not be calculated by Hooke's Law.
3) The material's Young's Modulus of spherical particles also be affected by temperature and pressure.Such as the material's Young's Modulus of sand particles (quartz) in the sandstone will decrease with higher temperature.

How the Arrangements of Particle Affected
In the unconsolidated sandstone formation, the normal contact external force f n of equal diameter spherical particles is related to the corresponding depth, principal stress, diameter of particles.The f n is the forces be decomposed by f in the Figure 2. Their has five types familiar arrangements modes of particles [6,8], their names and coordination number as Table 1 shows and their figures of modes shows as Figure 2. The particle packing rock layers tend to transform to make the direction of the direction of the maximum principal stress 1 Max ability of resistance contact force consistent with  (parallel to Z-axis in Figure 2).We may begins th the example that the v wi  is the maximum principal stress 1  in a rock layer a fault basin.The i in  is the angle force-decomposition.The force-decomp sition of five types familiar arrangements modes as Figure 2 shows.
The relations between corresponding total normal con of o tact force f n and 1  (parallel to Z-axis) of the particles arrangements modes i as formula 11 shows.
That the relations between th M e Macroscopic Young's odulus E g and 1  of the particles arrangements modes i as formula 13 shows.
The is the coefficient that can tell the relations be i k tween maximum principal stress 1  and its corresponding normal contact force f ni of a single particle under the affect of the particles arrangements modes i.
The k n is the area coefficient which express a single pa rticle undertake the stress/theoretical area.k ni is the area coefficient k n of the particles arrangements modes i .  We can find that the f in Figure 2 acc 15 ord with formula .
The k m is the coefficient whic be o the force de h express contact angle tween the two layers particles, The k mi is the coefficient k m of the particles arrangements modes i.The k z is a coefficient which relate to the coordination numbers N. The rhi v is the undertake volume of a single particle of the p cles arrangements modes i.The h i is the layer height of the particles arrangements modes i.
The k mi should be calculated according t Besides the particles arrangements mode 1-the Simple Cubic Mode, the zi k ( 2 3 4 5 , , ,  1 shows.If we sum and average the five types familiar particles arrangements modes simply.The sum and average can not representative the fact on the particle packing rock layers.But in a larger scale, the five familiar particles arrangements modes must all exist and closer to the average.The sum and average not means a bigger error.From the Table 1, the coordination number 1 N and vertical   R .Because differentiation of those material's Young's Modulus E of spherical particles is often not bi the particle p ing g, but ack diffe lay rentiation of the radius of particles R is often difference by several magnitudes.So the R is often the principal influencing factors of E g and maybe lead to the E g of sandstone bigger than that of mud rock.
In populous sedimentary basin and its most sedimentary stage, for example the loose sand ers in quaternary, in those layers, the principal stress in vertical direction v  is the Max principal stress 1  , that is say that v  parallel to Z-axis, the stress field is

Discussion
Young's Modulus of quartz bulk parent rock is often 40 In fact the quartz bulk parent rocks is d its Young's Modulus is less than that h than several hundred mete so nd h dulus of particle packing rock layers an GPa -100 GPa. also have pore an of material's of spherical sand particles.Other material's Young´s Modulus of material's has hardly been reported before.Because of the existing of some other mineral particles and clay particles, Young's Modulus of mineral particles maybe less than the quartz mineral particles, but maybe not less than the quartz bulk parent rock.LiuYu (2010) [7] thought that 20 GPa -80 GPa is a reasonable range of Young's Modulus of mineral particles E. Some data from Tian Jia-ning (1988) [9] tell that those Young's Modulus of rocks distributed from 5 GPa -60 GPa and very discrete.Those Young's Modulus of some rocks of clay particles is only 0.17 Ga, this less two magnitudes than that of sandstones, and maybe related to the spatial structure of those clay particles.Some shear wave velocity of particles packing rock layers which not deep than 15 m is often only 60 -200 m/s (Elnashai, 2008) [10] and far less than that of average value of sedimentary rock on the upper crust (2800 m/s -3500 m/s), also less than the compressional wave velocity of water (1400 m/s).Some data from An-Ou (1992 [11]) also tell that the Young's Modulus of rocks distributed from 0.1 GPa -20 GPa and very discrete.Those bigger differentiation of Young's Modulus and wave velocity verify that the Young's Modulus tend to approaching to zero when close to earth surface synchronously.
The former model on Macroscopic Young's Modulus of particle packing rock layers is only suitable for the rock layers not dept rs or me loose sandstone with rapid sedimentation.Plastic deformation maybe the dominating mode of the change of Macroscopic Young's Modulus of particle packing rock layers in a depth or an older stratum of rocks.In this paper, E g is proportional to 2/3  R , but some data express that is inverse proportion, for example An-Ou (1992: p. 34) [11], The paper think that because of those data was come from rocks in depth a ad undertake long-term plastic deformation and dynamic diagenesis.The plastic deformation velocity of little diameter particles is greater than that of the bigger diameter particles.The plastic deformation velocity model of particle packing rock layers need to be developed.We should also attend to that the instantaneous Young's Modulus will determine the plastic deformation velocity in a time future .The elastic deformation will runs through modern and the elastic deformation model will be the basement to distinguish the elastic deformation and plastic deformation.

Conclusions
1) The paper set up a model to calculate the Macro-scopic Young's Mo d the model includes the model that can tell how the arrangements of particle affected.
2) g E is proportional to 1/3 n f and proportional to 1/3 1  , can not be calculated by Hooke's Law.The Macro pic dry sco Young's Modulus of sand particle packing k layers which close to the earth surface is approaching to zero in 0 m depth, and the too little Macroscopic Young's Modulus is the main cause of earthquake site effect. 3)

Figure 1 (
b) shows.The leads to normal displacements u over the contact area.The distribution of normal deformation amounts u of different spherical particles after the normal contact force f n .So, the deformation amounts of v of different points in a contact interface as formula 3 shows.The   p r

Figure 1 .
Figure 1.Hertzian contact model [5].(a) Geometry of hertzian contact between two dissimilar, unequal sized elatstic spheres subjected to an applied normal force; (b) Normal contact force distribution of Hertzian contact.

Figure 2 .
Figure 2. Five types familiar arrangements modes of particles (coordinate system B: The Z-axis perpendicular to the paper surface,

arti composition angles i 
which based on the decomposition principle for the rces.The decomposition angles i fo  shows as Figure 2.
the Simple Cubic Mode is smallest and the bi st f the Regular Tetrahedron Mode.Not all the coordination number undertake the maximum principal stress than the ch e of the coordination numbe i The model has not finished ang l in above r N . .Because m pa ode ragraphs focused on the contact force f ni which parallel to the maximum principal stress 1  .When the normal contact force f ni not parallel to th maximum principal stress 1 e  , we need to convert the normal displacements  to parallel to the maximum principal stress 1  .So the gi E transform into the formula 16.

E 4 )
roc g is proportional to 2/3 R , the differentiation of the radius of particles R is often difference by several ma it gn es.So the R is ofte e When equal stress, the increased of the coordination number lead to the decrease of the normal contact force n f , this lead to the coordination number is an inverse proportion to Macroscopic Young's Modulus.But the larger coordination number change only means very little Macroscopic Young's Modulus change.

Table 1 . Coefficients of the five types familiar particles arra ents modes (the data of N, h i , v rh wa m [6]).
the H