Concrete Filled Steel Tube Section

[3]:

import numpy as np
from opstool.preprocessing import SecMesh, add_material, add_polygon, add_circle, offset, Rebars

[4]:

Ec = 3.45E7
Es = 2.0E8
Nus = 0.3
Nuc = 0.2
tube_d = 0.02
steel_mat = add_material(name='steel', elastic_modulus=Es, poissons_ratio=Nus)
conc_mat = add_material(
    name='conc', elastic_modulus=Ec, poissons_ratio=Nuc)
all = add_circle(xo=[0, 0], radius=0.75)
conc = add_circle(xo=[0, 0], radius=0.75 - tube_d)
tube = all - conc
conc.material = conc_mat
tube.material = steel_mat
sec = SecMesh()
sec.assign_group(dict(tube=tube, conc=conc))
sec.assign_mesh_size(dict(tube=0.1, conc=0.2))
sec.assign_group_color(dict(tube="gray", conc="green"))
sec.assign_ops_matTag(dict(tube=1, conc=2))
sec.mesh()
sec_props = sec.get_sec_props(
    Eref=Ec, display_results=False, plot_centroids=False)
# sec.rotate(45)
sec.view(fill=True, engine='p', save_html=None, on_notebook=True)

d:\AcademicSoft\Anaconda3\lib\site-packages\sectionproperties\analysis\section.py:198: ShapelyDeprecationWarning:

STRtree will be changed in 2.0.0 and will not be compatible with versions < 2.