def setUp(self):
self.profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND']
})
self.profile.title = "Test"
def test_applyfunction(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Total unit weight [kN/m3]': [19, 18, 17, 20],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, np.nan, 10, 40],
'qc to [MPa]': [2, np.nan, 10, 50]
})
profile.applyfunction(
function=voidratio_bulkunitweight,
outputkey="Void ratio [-]",
resultkey="e [-]",
parametermapping={'bulkunitweight': 'Total unit weight [kN/m3]'})
self.assertAlmostEqual(profile.loc[1, "Void ratio [-]"], 1.0625, 4)
profile.applyfunction(function=liquefactionprobability_moss,
outputkey='Liquefaction probability [pct]',
resultkey='Pl [pct]',
parametermapping={'qc': 'qc [MPa]'},
sigma_vo_eff=100,
Rf=0.4,
CSR=0.2,
CSR_star=0.2)
self.assertAlmostEqual(
profile.loc[2, 'Liquefaction probability from [pct]'], 27, 0)
self.assertAlmostEqual(
profile.loc[2, 'Liquefaction probability to [pct]'], 27, 0)
def test_calculate_parameter_center(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [1, 3, 10, 40]
})
profile.calculate_parameter_center('qc [MPa]')
self.assertEqual(profile.loc[2, 'qc center [MPa]'], 10)
def test_soilparameter_series(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 3, 20, 50]
})
self.assertEqual(profile.soilparameter_series("qc [MPa]")[1][2], 3)
self.assertEqual(profile.soilparameter_series("Dr [%]")[1][3], 60)
def test_cut_profile(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 4, 20, 50]
})
pf = profile.cut_profile(top_depth=1.5, bottom_depth=19)
self.assertEqual(pf.min_depth, 1.5)
self.assertEqual(pf.max_depth, 19)
self.assertEqual(pf.loc[0, "qc from [MPa]"], 3.125)
self.assertEqual(pf["qc to [MPa]"].iloc[-1], 49)
def test_parameter_mapping(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 3, 20, 50]
})
mapped_df = profile.map_soilprofile(np.linspace(0, 20, 21))
self.assertEqual(mapped_df.loc[3, 'Soil type'], 'SILT')
self.assertRaises(ValueError, profile.map_soilprofile,
(np.linspace(-1, 20, 22), ))
self.assertEqual(mapped_df.loc[15, 'qc [MPa]'], 45)
def test_wrong_linearvariation(self):
# Test incomplete linear parameter variation
self.assertRaises(IOError, sp.SoilProfile, ({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'qc from [MPa]': [1, 3, 10, 40]
}))
# Check that complete parameter variation passes
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [1, 3, 10, 40]
})
def test_convert_to_constant(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Total unit weight from [kN/m3]': [19, 18, 17, 20],
'Total unit weight to [kN/m3]': [19, 18, 19, 20],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 4, 20, 50]
})
profile.convert_to_constant("qc [MPa]")
profile.convert_to_constant("Total unit weight [kN/m3]", rule='min')
self.assertEqual(profile['qc [MPa]'].iloc[0], 1.5)
self.assertEqual(profile['Total unit weight [kN/m3]'].iloc[2], 17)
def test_calculate_overburden(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Total unit weight [kN/m3]': [19, 18, 17, 20],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 4, 20, 50]
})
profile.calculate_overburden(waterlevel=4)
self.assertEqual(profile['Total vertical stress to [kPa]'].iloc[-1],
376)
self.assertEqual(
profile['Effective vertical stress to [kPa]'].iloc[-1], 216)
self.assertEqual(profile['Hydrostatic pressure to [kPa]'].iloc[-1],
160)
def test_depth_integration(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Unit weight [kN/m3]': [9, 8, 7, 10],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 4, 20, 50]
})
profile.depth_integration(
parameter='Unit weight [kN/m3]',
outputparameter='Total vertical stress [kPa]')
self.assertEqual(profile['Total vertical stress from [kPa]'].iloc[0],
0)
self.assertEqual(profile['Total vertical stress to [kPa]'].iloc[0], 9)
self.assertEqual(profile['Total vertical stress to [kPa]'].iloc[1], 41)
def test_remove_parameter(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Relative density': ['Loose', 'Medium dense', None, 'Dense'],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [2, 3, 20, 50]
})
profile.remove_parameter('qc [MPa]')
def test_func():
return profile.loc[0, "qc from [MPa]"]
self.assertRaises(KeyError, test_func)
profile.remove_parameter('Dr [%]')
def test_func_2():
return profile.loc[0, "Dr [%]"]
self.assertRaises(KeyError, test_func_2)
def test_soilparameter_retrieval(self):
profile = sp.SoilProfile({
'Depth from [m]': [0, 1, 5, 10],
'Depth to [m]': [1, 5, 10, 20],
'Soil type': ['SAND', 'SILT', 'CLAY', 'SAND'],
'Dr [%]': [40, 60, np.nan, 80],
'qc from [MPa]': [1, 3, 10, 40],
'qc to [MPa]': [1, 3, 10, 40]
})
self.assertIn('qc [MPa]', profile.soil_parameters())
self.assertIn('Soil type', profile.soil_parameters())
self.assertIn('qc from [MPa]',
profile.soil_parameters(condense_linear=False))
self.assertIn('qc [MPa]', profile.numerical_soil_parameters())
self.assertNotIn('Soil type', profile.numerical_soil_parameters())
self.assertIn('qc from [MPa]',
profile.numerical_soil_parameters(condense_linear=False))
self.assertNotIn(
'Soil type',
profile.numerical_soil_parameters(condense_linear=False))
self.assertNotIn('qc from [MPa]', profile.string_soil_parameters())
self.assertIn('Soil type', profile.string_soil_parameters())
self.assertTrue(profile.check_linear_variation('qc [MPa]'))
self.assertFalse(profile.check_linear_variation('Dr [%]'))
