Elastic Modulus of Rock Calculator
Calculate Young's modulus of rock from stress-strain data or seismic wave velocities. Includes Poisson ratio, shear modulus, and rock classification.
Elastic Modulus of Rock Calculator
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: E = Stress / Strain (static) | E = 2G(1 + nu) (dynamic)
Worked example โ E = 66.67 GPa | Classification: Hard Rock
Formula
E = Stress / Strain (static) | E = 2G(1 + nu) (dynamic)
For the static method, divide axial stress (MPa) by axial strain (dimensionless). For the dynamic method using seismic velocities, first compute Poisson ratio and shear modulus from P-wave and S-wave velocities with rock density.
Worked Examples
Example 1: Laboratory Compression Test on Granite
Problem:A granite core sample is tested in uniaxial compression. At 50% of peak strength, the axial stress is 80 MPa and the axial strain is 0.0012. Calculate the elastic modulus.
Solution:Young's Modulus E = Stress / Strain E = 80 MPa / 0.0012 E = 66,667 MPa = 66.667 GPa Converted: 66.667 x 145.038 = 9,669 ksi Classification: Hard Rock (granite range: 40-70 GPa)
Result:E = 66.67 GPa | Classification: Hard Rock
Example 2: Dynamic Modulus from Seismic Survey
Problem:A seismic survey measures P-wave velocity of 5,000 m/s and S-wave velocity of 2,800 m/s in limestone with density 2,650 kg/m3.
Solution:Vp/Vs ratio = 5000/2800 = 1.786 Poisson's ratio = (5000^2 - 2x2800^2) / (2x(5000^2 - 2800^2)) = (25e6 - 15.68e6) / (2x(25e6 - 7.84e6)) = 9.32e6 / 34.32e6 = 0.2716 Shear modulus G = 2650 x 2800^2 = 20.776 GPa E = 2 x 20.776 x (1 + 0.2716) = 52.83 GPa
Result:E = 52.83 GPa | nu = 0.272 | G = 20.78 GPa
Frequently Asked Questions
What is the elastic modulus of rock and why is it important?
The elastic modulus (Young's modulus) of rock measures its stiffness, defined as the ratio of stress to strain within the elastic deformation range. It is one of the most fundamental mechanical properties used in geotechnical engineering, mining, tunneling, and petroleum engineering. A higher elastic modulus indicates a stiffer rock that deforms less under applied loads. Engineers need this value for designing foundations, analyzing slope stability, predicting ground settlement, planning tunnel support systems, and modeling subsurface reservoir behavior. The elastic modulus typically ranges from less than 1 GPa for very weak rocks to over 100 GPa for extremely hard crystalline rocks like quartzite.
How is elastic modulus determined from laboratory stress-strain tests?
In laboratory testing, cylindrical rock core samples are subjected to uniaxial or triaxial compression tests. The applied axial stress (force per unit area in MPa) is plotted against the measured axial strain (deformation divided by original length, dimensionless). Young's modulus is calculated from the slope of the linear portion of the stress-strain curve. Different conventions exist: tangent modulus uses the slope at 50% of ultimate stress, secant modulus draws a line from origin to a specific stress point, and average modulus calculates the slope of the best-fit line through the linear region. ISRM and ASTM standards provide detailed procedures for these measurements to ensure consistency and reproducibility.
What is the difference between static and dynamic elastic modulus?
Static elastic modulus is measured by physically deforming rock samples in compression tests, while dynamic elastic modulus is calculated from seismic wave velocities propagating through the rock. Dynamic values are typically 10% to 80% higher than static values because seismic waves involve very small, rapid strains that do not activate microcracks and other imperfections that reduce stiffness in static tests. The ratio of static to dynamic modulus depends on rock type, porosity, and degree of fracturing. Empirical correlations exist to convert between the two, such as E_static = 0.7 x E_dynamic for many sedimentary rocks. Dynamic testing is non-destructive and can be performed in the field using seismic surveys.
How do seismic velocities relate to elastic properties of rock?
P-wave (compressional) and S-wave (shear) velocities are directly related to rock elastic properties and density. Poisson's ratio is calculated as nu = (Vp^2 - 2Vs^2) / (2(Vp^2 - Vs^2)). The shear modulus G = density x Vs^2, and Young's modulus E = 2G(1 + nu). The bulk modulus K = density x (Vp^2 - 4/3 x Vs^2). These relationships assume an isotropic, homogeneous, linearly elastic medium. Typical P-wave velocities range from 2,000 m/s in weak sedimentary rocks to 7,000 m/s in dense igneous rocks. S-wave velocities are always lower than P-wave velocities, typically by a factor of about 1.7 for most rocks.
What factors affect the elastic modulus of rock?
Numerous factors influence rock elastic modulus. Mineralogy is primary: quartz-rich rocks like quartzite have high moduli (60-100 GPa), while clay-rich rocks like shale are much lower (1-20 GPa). Porosity inversely affects stiffness; increasing porosity from 5% to 25% can reduce modulus by 50% or more. Confining pressure increases modulus by closing microcracks. Water saturation generally increases dynamic modulus but can decrease static modulus in clay-bearing rocks. Temperature increases tend to reduce modulus. Weathering and alteration progressively decrease stiffness. Anisotropy due to bedding, foliation, or fracture sets means modulus varies with measurement direction, sometimes by factors of 2 to 5 in strongly foliated metamorphic rocks.
What is the difference between intact rock modulus and rock mass modulus?
The value returned from a core plug is the intact modulus E_i, the stiffness of unfractured material between joints. A tunnel or dam foundation loads a jointed rock mass, whose deformation modulus E_rm also includes closure and shear along bedding planes, joints and shear zones, so it is almost always lower than E_i. In heavily jointed ground the rock mass modulus can be a small fraction of the intact value, while in massive unjointed rock the two converge. Engineers bridge the gap with empirical correlations keyed to rock mass classification indices such as RMR, Q or GSI, or measure E_rm directly with plate-loading, dilatometer or flat-jack tests in situ.
Why must the S-wave velocity be lower than the P-wave velocity in the dynamic method?
P-waves travel through the P-wave modulus K + 4G/3, which always exceeds the shear modulus G that carries S-waves, so Vp is greater than Vs in any stable isotropic solid and the calculator rejects an input where Vs equals or exceeds Vp. There is a second, tighter check worth applying by eye: Poisson's ratio is positive only when Vp exceeds Vs times the square root of two, about 1.414, and most rocks sit around a Vp/Vs of 1.6 to 2.0. A ratio close to 1.4 suggests a mispicked shear arrival or a gas-charged interval, and a ratio above about 2.5 points to unconsolidated or fluid-rich material.
How does the elastic modulus control tunnel convergence and foundation settlement?
Elastic deformation scales inversely with stiffness, so halving the modulus doubles the predicted movement for the same load. For a circular opening in an elastic ground under a roughly hydrostatic stress p, the inward radial displacement of the wall is on the order of p times the radius times (1 + nu) divided by E, which is why a tunnel in 10 GPa mudstone converges several times more than the same tunnel in 60 GPa granite. Shallow foundation settlement follows the same inverse relationship, being proportional to bearing pressure and footing width divided by E. Because the movement is directly proportional to 1/E, an error in the modulus propagates one-for-one into the predicted displacement.
How are the P-wave and S-wave velocities for the dynamic method actually measured?
In the laboratory an ultrasonic pulse-transmission rig places piezoelectric transducers on the ends of a jacketed core, sends a pulse at several hundred kilohertz, and divides sample length by the first-arrival travel time; separate compressional and shear transducers give Vp and Vs. In the borehole a sonic or dipole shear log records the same quantities as slowness in microseconds per foot or per metre, with the dipole source needed because a monopole tool cannot excite a shear wave in slow formations. At the site scale, crosshole, downhole and surface-wave methods such as MASW return Vs profiles. All three scales need a matching bulk density, from core weighing or a density log, before a modulus can be computed.
Why is rock stiffness reported in GPa when the test stress is in MPa?
Modulus has the same dimensions as stress because strain is dimensionless, but rock strains elastically by only a few thousandths, so dividing a stress of tens of MPa by a strain near 0.001 lands three orders of magnitude higher. A stress of 80 MPa at a strain of 0.0012 gives 66,667 MPa, which is far more readable as 66.7 GPa. For conversion, 1 MPa equals 145.038 psi, so 1 GPa equals about 145,038 psi or 145.0 ksi, and the 66.7 GPa granite above is roughly 9,670 ksi. Always confirm whether a quoted figure is MPa or GPa before using it, since the thousandfold difference is the most common unit error in deformation analysis.
References
Background & Theory
History
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer ยท Editorial policy
Related Calculators
๐งฎElastic Moduli Converter (e, G, K, ฮ)
Calculate elastic moduli converter (e, g, k, ฮฝ) with inputs, formulas, and instant results.
๐งฎRock Density From Porosity Calculator
Calculate rock density from porosity with inputs, formulas, and instant results.
๐งฎRock Thermal Conductivity Estimator Calculator
Calculate rock thermal conductivity estimator with inputs, formulas, and instant results.
๐งฎEarth Energy Balance Calculator
Calculate earth energy balance with inputs, formulas, and instant results.
๐งฎBouguer Correction Calculator
Calculate bouguer correction with inputs, formulas, and instant results.
๐งฎCrustal Density Calculator
Calculate crustal density with inputs, formulas, and instant results.
๐งฎCrustal Thickness Calculator (H-k)
Calculate crustal thickness from receiver functions with inputs, formulas, and instant results.
๐งฎEarthquake Magnitude to Energy Calculator
Calculate earthquake magnitude to energy with inputs, formulas, and instant results.