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

# create geometry 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=[0])
sec = Section(geometry=geom)

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

# plot mohr's circle
post.plot_mohrs_circles(x=10, y=88.9)