# Sets
blender.add_nsets_from_layers(mdl, layers=['nset_supports', 'nset_load'])
# Materials
mdl.add([
Steel(name='mat_steel', fy=355),
Concrete(name='mat_concrete', fck=90),
])
# Sections
mdl.add([
ShellSection(name='sec_plate', t=0.005),
TrussSection(name='sec_tie', A=0.25 * 3.142 * 0.008**2),
SolidSection(name='sec_solid'),
])
# Properties
mdl.add([
Properties(name='ep_plate1',
material='mat_steel',
section='sec_plate',
elset='elset_top_plate'),
Properties(name='ep_plate2',
material='mat_steel',
section='sec_plate',
elset='elset_bot_plate'),
Properties(name='ep_ties',
# Sets
rhino.add_sets_from_layers(mdl, layers=['nset_pin', 'nset_roller'])
# Materials
mdl.add([
Concrete(name='mat_concrete', fck=50),
Steel(name='mat_steel', fy=460),
])
# Sections
mdl.add([
ShellSection(name='sec_planar', t=0.050),
TrussSection(name='sec_tie', A=0.0001),
])
# Properties
mdl.add([
Properties(name='ep_planar',
material='mat_concrete',
section='sec_planar',
elset='elset_mesh'),
Properties(name='ep_plate',
material='mat_steel',
section='sec_planar',
elset='elset_plates'),
Properties(name='ep_tie',
material='mat_steel',
blender.add_nodes_elements_from_layers(mdl,
line_type='TrussElement',
layers='elset_truss')
# Sets
blender.add_nsets_from_layers(mdl, layers=['nset_pins', 'nset_top'])
# Materials
mdl.add(ElasticIsotropic(name='mat_elastic', E=200 * 10**9, v=0.3, p=7850))
# Sections
mdl.add(TrussSection(name='sec_truss', A=0.0001))
# Properties
mdl.add(
Properties(name='ep_truss',
material='mat_elastic',
section='sec_truss',
elset='elset_truss'))
# Displacements
mdl.add(PinnedDisplacement(name='disp_pinned', nodes='nset_pins'))
# Loads
nodes_bot = [i for i, node in mdl.nodes.items() if node.y < ymin + 0.01]
mdl.add_set(name='nset_top', type='node', selection=nodes_top)
mdl.add_set(name='nset_bot', type='node', selection=nodes_bot)
# Materials
mdl.add([
Concrete(name='mat_concrete', fck=90),
Steel(name='mat_steel', fy=355),
])
# Sections
mdl.add([
ShellSection(name='sec_mesh', t=0.004),
TrussSection(name='sec_ties', A=0.25 * pi * 0.010**2),
RectangularSection(name='sec_ends', b=0.030, h=0.030),
])
# Properties
mdl.add([
Properties(name='ep_mesh',
material='mat_concrete',
section='sec_mesh',
elset='elset_mesh'),
Properties(name='ep_ties',
material='mat_steel',
section='sec_ties',
elset='elset_ties'),
Properties(name='ep_ends',
# Sets
layers = ['nset_corners', 'nset_corner1', 'nset_corner2']
rhino.add_sets_from_layers(mdl, layers=layers)
# Materials
mdl.add_materials([
Concrete(name='mat_concrete', fck=90),
Steel(name='mat_steel', fy=355)])
# Sections
mdl.add_sections([
ShellSection(name='sec_concrete', t=0.020),
TrussSection(name='sec_ties', A=0.0004)])
# Properties
mdl.add_element_properties([
Properties(name='ep_concrete', material='mat_concrete', section='sec_concrete', elsets='elset_concrete'),
Properties(name='ep_steel', material='mat_steel', section='sec_ties', elsets='elset_ties')])
# Displacements
mdl.add_displacements([
RollerDisplacementXY(name='disp_roller', nodes='nset_corners'),
PinnedDisplacement(name='disp_pinned', nodes='nset_corner1'),
GeneralDisplacement(name='disp_xdof', nodes='nset_corner2', x=0)])
# Loads
layers=[0])
# Sets
blender.add_elset_from_bmeshes(mdl, layer=0, name='elset_struts')
blender.add_nset_from_objects(mdl, layer=1, name='nset_pins')
blender.add_nset_from_objects(mdl, layer=2, name='nset_top')
# Materials
mdl.add_material(
ElasticIsotropic(name='mat_elastic', E=200 * 10**9, v=0.3, p=7850))
# Sections
mdl.add_section(TrussSection(name='sec_truss', A=0.00010))
# Properties
ep = Properties(name='ep_strut',
material='mat_elastic',
section='sec_truss',
elsets='elset_struts')
mdl.add_element_properties(ep)
# Displacements
mdl.add_displacement(PinnedDisplacement(name='disp_pinned', nodes='nset_pins'))
# Loads
blender.add_nsets_from_layers(mdl, layers=['nset_corner1', 'nset_corner2'])
edges = [i for i in mdl.nodes if mdl.nodes[i].z < 0.001]
mdl.add_set(name='nset_edges', type='node', selection=edges)
# Materials
mdl.add([
Concrete(name='mat_concrete', fck=90, fr=[1.16, 0.15]),
Steel(name='mat_steel', fy=355),
])
# Sections
mdl.add([
ShellSection(name='sec_floor', t=0.050),
TrussSection(name='sec_ties', A=pi * 0.25 * 0.030**2),
])
# Properties
mdl.add([
Properties(name='ep_floor',
material='mat_concrete',
section='sec_floor',
elset='elset_floor'),
Properties(name='ep_ties',
material='mat_steel',
section='sec_ties',
elset='elset_ties'),
])
line_type='TrussElement',
layers=layers)
# Sets
rhino.add_sets_from_layers(mdl,
layers=['nset_pins', 'nset_load_v', 'nset_load_h'])
# Materials
mdl.add_material(Steel(name='mat_steel', fy=355))
# Sections
mdl.add_sections([
TrussSection(name='sec_main', A=0.0008),
TrussSection(name='sec_diag', A=0.0005),
TrussSection(name='sec_stays', A=0.0001)
])
# Properties
mdl.add_element_properties([
Properties(name='ep_main',
material='mat_steel',
section='sec_main',
elsets='elset_main'),
Properties(name='ep_diag',
material='mat_steel',
section='sec_diag',
elsets='elset_diag'),
rhino.add_nodes_elements_from_layers(mdl, line_type='TrussElement', layers='elset_main', pL=A1*p)
rhino.add_nodes_elements_from_layers(mdl, line_type='TrussElement', layers='elset_diag', pL=A2*p)
rhino.add_nodes_elements_from_layers(mdl, line_type='TrussElement', layers='elset_stays', pL=A3*p)
# Sets
rhino.add_sets_from_layers(mdl, layers=['nset_pins', 'nset_load_v', 'nset_load_h'])
# Materials
mdl.add(Steel(name='mat_steel', fy=355, p=p))
# Sections
mdl.add([
TrussSection(name='sec_main', A=A1),
TrussSection(name='sec_diag', A=A2),
TrussSection(name='sec_stays', A=A3),
])
# Properties
mdl.add([
Properties(name='ep_main', material='mat_steel', section='sec_main', elset='elset_main'),
Properties(name='ep_diag', material='mat_steel', section='sec_diag', elset='elset_diag'),
Properties(name='ep_stays', material='mat_steel', section='sec_stays', elset='elset_stays'),
])
# Displacements
mdl.add(PinnedDisplacement(name='disp_pinned', nodes='nset_pins'))