Exemplos de Structure.get_nodal_results em Python

Exemplo n.º 1

Arquivo: beam_frame_rhino.py Projeto: ming91915/compas_fea

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'))

Exemplo n.º 2

               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'))

Exemplo n.º 3

Arquivo: shells.py Projeto: OVGULIU/compas_fea

# 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]))

Exemplo n.º 4

Arquivo: truss_tower_blender.py Projeto: OVGULIU/compas_fea

# 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'))

Exemplo n.º 5

Arquivo: example_truss.py Projeto: compas-Workshops/WS_compas_fea

# 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