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)