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)
mdl.write_input_file(software='sofistik')
# Run (Abaqus/OpenSees)
#mdl.analyse_and_extract(software='abaqus', fields=['u', 'ur', 'rf'])
mdl.analyse_and_extract(software='opensees', fields=['u', 'ur', 'rf'])
rhino.plot_data(mdl, step='step_loads', field='um', scale=50)
rhino.plot_data(mdl, step='step_loads', field='urm', scale=50)
print(mdl.get_nodal_results(step='step_loads', field='rfx', nodes='nset_pins'))
print(mdl.get_nodal_results(step='step_loads', field='rfz', nodes='nset_pins'))
section='sec_beam',
elset='elset_beams'))
# Displacements
mdl.add(FixedDisplacement(name='disp_fixed', nodes='nset_support'))
# Loads
mdl.add(PointLoad(name='load_point', nodes='nset_load', z=600))
# Steps
mdl.add([
GeneralStep(name='step_bc', displacements=['disp_fixed']),
GeneralStep(name='step_load', loads=['load_point']),
])
mdl.steps_order = ['step_bc', 'step_load']
# Summary
mdl.summary()
# Run
mdl.analyse_and_extract(software='abaqus', fields=['u', 'sf', 'sm'])
rhino.plot_data(mdl, step='step_load', field='uz', radius=1)
print(mdl.get_nodal_results(step='step_load', field='uz', nodes='nset_load'))
# Displacements
mdl.add(FixedDisplacement(name='disp_pins', nodes=[0, 2, 8, 6]))
# Loads
mdl.add(PointLoad(name='load_v', nodes=[4], z=-1000))
# Steps
mdl.add(
GeneralStep(name='step_bc_loads',
displacements=['disp_pins'],
loads=['load_v'],
nlgeom=False))
mdl.steps_order = ['step_bc_loads']
# Summary
# mdl.summary()
# Run (Abaqus)
# mdl.analyse_and_extract(software='abaqus', fields=['u', 'rf'], license='research')
# Run (Ansys)
mdl.analyse_and_extract(software='ansys', fields=['u'], license='research')
print(mdl.get_nodal_results(step='step_bc_loads', field='um', nodes=[3, 4, 5]))
# Loads
mdl.add(PointLoad(name='load_top', nodes='nset_top', x=2000, y=1000,
z=-100000))
# Steps
mdl.add([
GeneralStep(name='step_bc', displacements='disp_pinned'),
GeneralStep(name='step_load', loads='load_top'),
])
mdl.steps_order = ['step_bc', 'step_load']
# Summary
mdl.summary()
# Run
mdl.analyse_and_extract(software='abaqus', fields=['u', 'cf', 'rf', 'sf', 's'])
# blender.plot_data(mdl, step='step_load', field='um')
# blender.plot_data(mdl, step='step_load', field='smises')
# blender.plot_reaction_forces(mdl, step='step_load', scale=0.05)
# blender.plot_concentrated_forces(mdl, step='step_load', scale=0.05)
# Print results
print(mdl.get_nodal_results(step='step_load', field='um', nodes='nset_top'))
print(mdl.get_nodal_results(step='step_load', field='rfm', nodes='nset_pins'))
# print(mdl.loads['gravity'])
mdl.add([
GeneralStep(name='bc', displacements='pinned'),
GeneralStep(name='loads', loads=['pointloads', 'gravity'], factor=1.5),
])
mdl.steps_order = ['bc', 'loads']
# print(mdl.steps['bc'])
# print(mdl.steps['loads'])
mdl.summary()
mdl.analyse_and_extract(software='abaqus', fields=['u', 's', 'cf', 'rf'])
print(mdl.get_nodal_results(step='loads', field='rfm', nodes='supports'))
# load / show .inp and .odb files
# show raw data, discuss methods to extract data
rhino.plot_data(mdl,
step='loads',
field='um',
radius=0.1,
scale=10,
cbar_size=0.3)
rhino.plot_data(mdl, step='loads', field='smises', radius=0.1, cbar_size=0.3)
rhino.plot_reaction_forces(mdl, step='loads', scale=0.1)
rhino.plot_concentrated_forces(mdl, step='loads', scale=0.1)
# edit parameters and geometry