from sectionproperties.pre import Material
from sectionproperties.pre.library import concrete_tee_section
from sectionproperties.analysis import Section

concrete = Material(
    name="Concrete",
    elastic_modulus=30.1e3,
    poissons_ratio=0.2,
    yield_strength=32,
    density=2.4e-6,
    color="lightgrey",
)
steel = Material(
    name="Steel",
    elastic_modulus=200e3,
    poissons_ratio=0.3,
    yield_strength=500,
    density=7.85e-6,
    color="grey",
)

geom = concrete_tee_section(
    d=900,
    b=450,
    d_f=150,
    b_f=1800,
    dia_top=24,
    area_top=450,
    n_top=12,
    c_top=42,
    dia_bot=28,
    area_bot=620,
    n_bot=5,
    c_bot=42,
    dia_side=16,
    area_side=200,
    n_side=4,
    c_side=42,
    n_circle=16,
    conc_mat=concrete,
    steel_mat=steel,
)

geom.create_mesh(mesh_sizes=[0])  # a size of zero creates a coarse mesh
Section(geometry=geom).plot_mesh()