def shear_force(self, factor=None, figsize=None, verbosity=0, scale=1, offset=(0, 0), show=True, gridplot=False):
self.plot_structure(figsize, 1, scale=scale, offset=offset, gridplot=gridplot)
if factor is None:
max_force = max(map(lambda el: np.max(np.abs(el.shear_force)), self.system.element_map.values()))
factor = det_scaling_factor(max_force, self.max_val_structure)
for el in self.system.element_map.values():
if math.isclose(el.node_1.Ty, 0, rel_tol=1e-5, abs_tol=1e-9) and \
math.isclose(el.node_2.Ty, 0, rel_tol=1e-5, abs_tol=1e-9) and el.q_load is None:
# If True there is no bending moment and no shear, thus no shear force, so no need for plotting.
continue
axis_values = plot_values_shear_force(el, factor)
shear_1 = el.shear_force[0]
shear_2 = el.shear_force[-1]
if verbosity == 0:
node_results = True
else:
node_results = False
self.plot_result(axis_values, shear_1, shear_2, node_results=node_results)
if show:
self.plot()
else:
return self.fig
def shear_force(
self,
factor=None,
figsize=None,
verbosity=0,
scale=1,
offset=(0, 0),
show=True,
gridplot=False,
include_structure=True,
figure=None,
):
if include_structure:
self.plot_structure(
figsize, 1, scale=scale, offset=offset, gridplot=gridplot, figure=figure
)
con = len(self.system.element_map[1].shear_force)
if factor is None:
max_force = max(
map(
lambda el: np.max(np.abs(el.shear_force)),
self.system.element_map.values(),
)
)
factor = det_scaling_factor(max_force, self.max_val_structure)
color_counter = 0
for el in self.system.element_map.values():
if (
math.isclose(el.node_1.Ty, 0, rel_tol=1e-5, abs_tol=1e-9)
and math.isclose(el.node_2.Ty, 0, rel_tol=1e-5, abs_tol=1e-9)
and math.isclose(el.all_qp_load[0], 0, rel_tol=1e-5, abs_tol=1e-9)
and math.isclose(el.all_qp_load[1], 0, rel_tol=1e-5, abs_tol=1e-9)
):
# If True there is no bending moment and no shear, thus no shear force, so no need for plotting.
continue
iteration_factor = np.linspace(0, 1, con)
x = iteration_factor * el.l
eq = np.polyfit(x, el.bending_moment, 3)
x2, x = symbols("x² x")
section_function = str(round(-eq[0] * 3, 2) * x2 + round(-eq[1] * 2, 2) * x + round(-eq[2], 2))
axis_values = plot_values_shear_force(el, factor)
shear_1 = el.shear_force[0]
shear_2 = el.shear_force[-1]
self.plot_result(
axis_values, shear_1, shear_2, node_results=not bool(verbosity),
color=self.cs.color_list[color_counter], label=section_function
)
color_counter += 1
if color_counter > len(self.cs.color_list) - 1:
color_counter = 0
if show:
self.plot()
else:
return self.fig
def shear_force(self, factor=None, figsize=None, verbosity=0, scale=1, offset=(0, 0), show=True, gridplot=False):
self.plot_structure(figsize, 1, scale=scale, offset=offset, gridplot=gridplot)
if factor is None:
max_force = max(map(lambda el: np.max(np.abs(el.shear_force)), self.system.element_map.values()))
factor = det_scaling_factor(max_force, self.max_val_structure)
for el in self.system.element_map.values():
if math.isclose(el.node_1.Ty, 0, rel_tol=1e-5, abs_tol=1e-9) and \
math.isclose(el.node_2.Ty, 0, rel_tol=1e-5, abs_tol=1e-9) and el.q_load is None:
# If True there is no bending moment and no shear, thus no shear force, so no need for plotting.
continue
axis_values = plot_values_shear_force(el, factor)
shear_1 = el.shear_force[0]
shear_2 = el.shear_force[-1]
if verbosity == 0:
node_results = True
else:
node_results = False
self.plot_result(axis_values, shear_1, shear_2, node_results=node_results)
if show:
self.plot()
else:
return self.fig