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

timber0 = Material(
    name="Timber0",
    elastic_modulus=9.5e3,
    poissons_ratio=0.35,
    density=4.4e-7,
    yield_strength=5.5,
    color="burlywood",
)
timber90 = Material(
    name="Timber90",
    elastic_modulus=317,
    poissons_ratio=0.35,
    density=4.4e-7,
    yield_strength=5.5,
    color="orange",
)

geom = clt_rectangular_section(
    d=[40, 40, 40],
    layer_mat=[timber0, timber90, timber0],
    b=1000
)

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