fluid = getJsonData(folderName + "fluid.json")
solid = getJsonData(folderName + "solid.json")
h_upper = 6.0
h_lower = 4.0
return tz.Solution(h_upper, h_lower, tao, solid, fluid)
folderName = "results//FIM//"
# folderName = "results//FSS//"
print(folderName + "p.txt")
res_p = ReadResults(folderName + "p.txt")
res_u = ReadResults(folderName + "u.txt")
L = res_p.coord[-1]
ic = getJsonData(folderName + "IC.json")
g = ic.get("Gravity")
bound_u = getJsonData(folderName + "BC_u.json")
top_stress = bound_u.get("TOP").get("Value")
terza = getTerza(folderName, L, top_stress, g)
np = 50
z_terza = terza.getPositionValues(np)
kPa = 1e-3
mm = 1e-3
plt.figure(figsize=(15, 6))
plt.subplots_adjust(left=0.05, right=0.95, bottom=0.1, top=0.95)
times = [1, 70, -1]
def getTerza(folderName, H, tao, g):
fluid = getJsonData(folderName + "fluid.json")
solid = getJsonData(folderName + "solid.json")
h_upper = 6.0
h_lower = 4.0
return tz.Solution(h_upper, h_lower, tao, solid, fluid)
def getTerza(folderName, H, tao, g):
fluid = getJsonData(folderName + "fluid.json")
solid = getJsonData(folderName + "solid.json")
return tz.Solution(H, tao, solid, fluid, -g)