コード例 #1
0
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'))
コード例 #2
0
# 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'))
コード例 #3
0
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'))