mdl.add([
PinnedDisplacement(name='disp_pins', nodes='pins'),
RollerDisplacementXZ(name='disp_roller', nodes='middle'),
])
# Loads
mdl.add(PointLoad(name='load_middle', nodes='middle', z=-500))
# Steps
mdl.add([
GeneralStep(name='step_bc', displacements=['disp_pins', 'disp_roller']),
GeneralStep(name='step_load', loads='load_middle'),
])
mdl.steps_order = ['step_bc', 'step_load']
# Summary`
mdl.summary()
# Run
mdl.analyse_and_extract(software='opensees', fields=['u', 'spf'], ndof=3)
rhino.plot_data(mdl, step='step_load', field='um', radius=0.02)
rhino.plot_data(mdl, step='step_load', field='spfx', radius=0.02)
print(mdl.get_element_results(step='step_load', field='spfx'))
# Loads
mdl.add_load(PointLoad(name='load_top', nodes='nset_top', z=-100000))
# Steps
mdl.add_steps([
GeneralStep(name='step_bc', displacements='disp_pinned'),
GeneralStep(name='step_load', loads='load_top', nlmat=False)
])
mdl.steps_order = ['step_bc', 'step_load']
# Summary
mdl.summary()
# Run (Abaqus)
exe = '/home/al/abaqus/Commands/abaqus cae '
mdl.analyse_and_extract(software='abaqus', exe=exe, fields=['u', 'rf', 's'])
blender.plot_data(mdl, step='step_load', field='um', layer=3)
print(mdl.get_nodal_results(step='step_load', field='um', nodes='nset_top'))
blender.plot_data(mdl, step='step_load', field='sxx', layer=4)
print(mdl.get_element_results(step='step_load', field='sxx', elements=[10]))
blender.plot_data(mdl, step='step_load', field='rfm', layer=5)
print(mdl.get_nodal_results(step='step_load', field='rfm', nodes='nset_pins'))
mdl.add([
PointLoad(name='load_v', nodes='nset_load_v', z=-15500),
PointLoad(name='load_h', nodes='nset_load_h', x=5000),
GravityLoad(name='load_gravity', elements=['elset_diag', 'elset_main']),
])
# Steps
mdl.add([
GeneralStep(name='step_bc', displacements=['disp_pinned']),
GeneralStep(name='step_loads', loads=['load_v', 'load_h', 'load_gravity'], factor=1.5, increments=300),
])
mdl.steps_order = ['step_bc', 'step_loads']
# Summary
mdl.summary()
# Run
mdl.analyse_and_extract(software='abaqus', fields=['u', 's', 'sf', 'cf', 'rf'], ndof=3)
rhino.plot_data(mdl, step='step_loads', field='um', radius=0.1, scale=10, cbar_size=0.3)
rhino.plot_data(mdl, step='step_loads', field='sxx', radius=0.1, cbar_size=0.3) # abaqus:sxx opensees:sf1
rhino.plot_reaction_forces(mdl, step='step_loads', scale=0.05)
rhino.plot_concentrated_forces(mdl, step='step_loads', scale=0.2)
print(mdl.get_nodal_results(step='step_loads', field='um', nodes='nset_load_v'))
print(mdl.get_nodal_results(step='step_loads', field='rfm', nodes='nset_pins'))
print(mdl.get_element_results(step='step_loads', field='sxx', elements='elset_main'))