# Properties
ep = Properties(name='ep', material='mat_steel', section='sec_pipe', elsets='elset_lines')
mdl.add_element_properties(ep)
# Displacements
mdl.add_displacements([
PinnedDisplacement(name='disp_pins', nodes='nset_pins'),
RollerDisplacementXZ(name='disp_rollers', nodes='nset_rollers')])
# Loads
mdl.add_loads([
PointLoad(name='load_h', nodes='nset_load_h', x=4000),
PointLoad(name='load_v', nodes='nset_load_v', z=-6000),
])
# Steps
mdl.add_steps([
GeneralStep(name='step_bc', displacements=['disp_pins', 'disp_rollers']),
GeneralStep(name='step_loads', loads=['load_h', 'load_v'], iterations=50)])
mdl.steps_order = ['step_bc', 'step_loads']
# Summary
mdl.summary()
# Run (Sofistik)
elsets='elset_ties'),
Properties(name='ep_concrete',
material='mat_concrete',
section='sec_solid',
elsets='elset_tets')
])
# Displacements
mdl.add_displacement(
RollerDisplacementY(name='disp_rollers', nodes='nset_supports'))
# Loads
mdl.add_loads([
GravityLoad(name='load_gravity', elements='elset_tets'),
PointLoad(name='load_point', nodes='nset_load', z=-5)
])
# Steps
mdl.add_steps([
GeneralStep(name='step_bc', displacements=['disp_rollers']),
GeneralStep(name='step_loads',
loads=['load_gravity', 'load_point'],
factor=1.35)
])
mdl.steps_order = ['step_bc', 'step_loads']
# Summary
mdl.summary()
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
mesh = mesh_from_guid(Mesh(), rs.ObjectsByLayer('load_mesh')[0])
mdl.add_loads([
GravityLoad(name='load_gravity', elements='elset_concrete'),
PrestressLoad(name='load_prestress', elements='elset_ties', sxx=50*10**6),
TributaryLoad(mdl, name='load_tributary', mesh=mesh, z=-2000)])
# Steps
mdl.add_steps([
GeneralStep(name='step_bc', displacements=['disp_roller', 'disp_pinned', 'disp_xdof']),
GeneralStep(name='step_prestress', loads=['load_prestress']),
GeneralStep(name='step_loads', loads=['load_gravity', 'load_tributary'], factor=1.1)])
mdl.steps_order = ['step_bc', 'step_prestress', 'step_loads']
# Summary
mdl.summary()
# Run (Abaqus)
Properties(name='ep_stays',
material='mat_steel',
section='sec_stays',
elsets='elset_stays')
])
# Displacements
mdl.add_displacement(PinnedDisplacement(name='disp_pinned', nodes='nset_pins'))
# Loads
# Note: GravityLoad doesnt activate for OpenSees
mdl.add_loads([
PointLoad(name='load_pl_v', nodes='nset_load_v', z=-15500),
PointLoad(name='load_pl_h', nodes='nset_load_h', x=5000),
GravityLoad(name='load_gravity', elements=['elset_diag', 'elset_main'])
])
# Steps
mdl.add_steps([
GeneralStep(name='step_bc', displacements=['disp_pinned']),
GeneralStep(name='step_loads',
loads=['load_pl_v', 'load_pl_h', 'load_gravity'],
factor=1.5,
increments=200)
])
mdl.steps_order = ['step_bc', 'step_loads']
# Summary