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

geom = channel_section(d=200, b=75, t_f=12, t_w=6, r=12, n_r=8)
geom.create_mesh(mesh_sizes=[20])

section = Section(geometry=geom)
section.calculate_geometric_properties()
section.calculate_warping_properties()
section.calculate_plastic_properties()

section.plot_centroids()