def setUp(self) -> None:
self.gef = ParseGEF("files/example.gef")
self.layer_table = pd.read_csv("files/layer_table.csv")
self.calc = PileCalculationSettlementDriven(
self.gef,
self.d_eq,
self.circum,
self.area,
self.layer_table,
pile_load=1500,
soil_load=10,
pile_system="soil-displacement",
ocr=1,
elastic_modulus_pile=30e3,
settlement_time_in_days=int(1e10),
alpha_s=0.01,
gamma_m=1,
alpha_p=1,
beta_p=1,
pile_factor_s=1,
)
# Assert that tipping points are as yet none
self.assertIsNone(self.calc.negative_friction_tipping_point_nap)
self.assertIsNone(self.calc.positive_friction_tipping_point_nap)
def setUp(self) -> None:
self.cpts = []
self.layer_tables = []
basedir = "data/cpt-grid/"
for cpt_path in os.listdir(basedir):
self.cpts.append(ParseGEF(basedir + cpt_path))
self.layer_tables.append(
pd.read_csv("data/layer_tables_grid/" +
cpt_path.replace("gef", "csv")))
def setUp(self) -> None:
self.gef = ParseGEF("files/test_gef2.gef")
self.layer_table = pd.read_csv("files/test_layer_table2.csv")
self.calc = PileCalculationSettlementDriven(
self.gef,
self.d_eq,
self.circum,
self.area,
self.layer_table,
pile_load=1000,
soil_load=1,
)
from pygef import ParseGEF, nap_to_depth
from anapile.pressure import bearing
from anapile.geo import soil
import pandas as pd
import numpy as np
# read cpt and layers
gef = ParseGEF("../tests/files/example.gef")
layer_table = pd.read_csv("../tests/files/layer_table.csv")
# we know the groundwater level is 1 m NAP
gef.groundwater_level = 1
# backfill cpt values with soil parameters
df = soil.join_cpt_with_classification(gef, layer_table)
# pile tip level -13.5 m
ptl = nap_to_depth(gef.zid, -13.5)
pile_width = 0.25
circum = pile_width * 4
# find the index of the pile tip
idx_ptl = np.argmin(np.abs(gef.df.depth.values - ptl))
ptl_slice = slice(0, idx_ptl)
# assume the tipping point is on top of the sand layers.
tipping_point = soil.find_last_negative_friction_tipping_point(
df.depth.values[ptl_slice], df.soil_code[ptl_slice])
idx_tp = np.argmin(np.abs(df.depth.values - tipping_point))
negative_friction_slice = slice(0, idx_tp)
def setUp(self) -> None:
self.gef = ParseGEF("files/example.gef")
self.soil_prop = pd.read_csv("files/soil_prop_example.csv")
self.layer_table = pd.read_csv("files/layer_table.csv")