from sectionproperties.pre.library import angle_section
from sectionproperties.analysis import Section

# create geometry, mesh and section
geom = angle_section(d=150, b=90, t=12, r_r=10, r_t=5, n_r=8)
geom.create_mesh(mesh_sizes=20)
sec = Section(geometry=geom)

# conduct analyses
sec.calculate_geometric_properties()
sec.calculate_warping_properties()
stress = sec.calculate_stress(
    n=50e3, mxx=-5e6, m22=2.5e6, mzz=0.5e6, vx=10e3, vy=5e3
)

# plot stress contour
stress.plot_stress(stress="vm", cmap="viridis", normalize=False)