# Chapter6-Soil Compaction¶

## Ex2-pg127¶

In [3]:
import math
#calculate maximum dry density and optimum moisture content
##initialisation of variables
G= 2.6
LL= 20.
P= 20.
##calclations
R= (4804574.*G-195.55*(LL)**2+156971*(P)**0.5-9527830)**0.5
n= (1.195e-4)*((LL)**2)-1.964*G-(6.617e-5)*(P)+7.651
w= math.e**n
##results
print'%s %.1f %s'% ('maximum dry density = ',R,' kg/m^3 ')
print'%s %.2f %s'%('optimum moisture content = ',w,' ')

maximum dry density =  1894.2  kg/m^3
optimum moisture content =  13.34


## Ex3-pg143¶

In [3]:
import math
#calculate dry unit weight of compaction in the field and dry unit weight of compaction in the field
##initialisation of variables
do= 1570. ##kg/m^3
mo= 0.545 ##kg
M1= 7.59 ##kg
M2= 4.78 ##kg
M3= 3.007 ##kg
w= 0.102 ##
dmax= 19. ##KN/m^3
##calculations
Ms= M1-M2
Mc= Ms-mo
Vh= Mc/do
Dc= M3/Vh
Du= Dc*9.81/1000.
f= Du/(1.+w)
Rc= f*100./dmax
##results
print'%s %.2f %s'% ('dry unit weight of compaction in the field = ',f,' kN/m^3 ')
print'%s %.1f %s'% ('relative compaction in the field = ',Rc,'')
#calculate the value of gamma and plot the graph
%matplotlib inline
import warnings
warnings.filterwarnings('ignore')
import math
from math import log
import numpy
from math import tan
import matplotlib
from matplotlib import pyplot
#given
p=numpy.array([6,8,9,11,12,14])
e=numpy.array([14.80,17.45,18.52,18.9,18.5,16.9])

#calculations

#results

pyplot.plot(p,e)
pyplot.xlabel('gamma')
pyplot.ylabel('weight ,w')
pyplot.title('Graph of gamma vs w')
pyplot.show()
print('look at the axis reverse in text book')

dry unit weight of compaction in the field =  18.55  kN/m^3
relative compaction in the field =  97.7

look at the axis reverse in text book


## Ex4-pg155¶

In [1]:
import math
#calculate sustainabilty number
##initialisation of variables
D1= 0.36 ##mm
D2= 0.52 ##mm
D5= 1.42 ##mm
##calculations
Sn= 1.7*(math.sqrt((3./(D5)**2)+(1./(D2)**2)+(1./(D1)**2)))
##results
print'%s %.1f %s'% ('sustainabilty number = ',Sn,' ')

sustainabilty number =  6.1