Beispiel #1
0
    for element in S.elements.values():
        faces.append(element.nodes)

    try:
        clear_layer(layer_output)
    except:
        create_layer(layer_output)

    vertices = S.nodes_xyz()

    xdraw_mesh(name='welded_mesh',
               vertices=vertices,
               faces=faces,
               layer=layer_output)


# ==============================================================================
# Debugging
# ==============================================================================

if __name__ == "__main__":

    from compas_fea.structure import Structure

    # mdl = Structure.load_from_obj(filename='/home/al/compas/compas_fea/data/_workshop/example_tets.obj')
    # plot_voxels(mdl, step='step_load', field='smises', vdx=0.100)

    mdl = Structure.load_from_obj(filename='C:/Temp/block_tets.obj')
    plot_voxels(mdl, step='step_load', field='um', vdx=0.010)
Beispiel #2
0
                result = postprocess(self.nodes, self.elements, self.ux,
                                     self.uy, self.uz, data, 'element', 1,
                                     cbar, 255, iptype, nodal)
                toc, _, cnodes, *_ = result

                self.update_vertices_colors(
                    {i: j
                     for i, j in enumerate(cnodes)})
                self.update_statusbar('Plotting: {0:.3f} s'.format(toc))

        except:

            pass


# ==============================================================================
# Main
# ==============================================================================

if __name__ == "__main__":

    from compas_fea.structure import Structure

    # fnm = '/home/al/compas/compas_fea/data/_workshop/example_truss.obj'
    # fnm = '/home/al/compas/compas_fea/data/_workshop/example_shell.obj'
    fnm = '/home/al/compas/compas_fea/data/_workshop/example_tets.obj'

    mdl = Structure.load_from_obj(fnm)
    mdl.view()
Beispiel #3
0
        dtype = 'element'

    # Postprocess

    result = postprocess(nodes, elements, ux, uy, uz, data, dtype, 1, cbar, 1,
                         iptype, nodal)

    try:
        toc, U, cnodes, fabs, fscaled, celements, eabs = result
        U = array(U)
        print('\n***** Data processed : {0:.3f} s *****'.format(toc))
    except:
        print(
            '\n***** Error encountered during data processing or plotting *****'
        )

    plotvoxels(values=fscaled, U=U, vdx=vdx, plot=plot)


# ==============================================================================
# Debugging
# ==============================================================================

if __name__ == "__main__":

    from compas_fea.structure import Structure

    mdl = Structure.load_from_obj(filename='/home/al/temp/block_deepbeam.obj')

    plot_voxels(mdl, step='step_load', cbar=[0, 2], vdx=0.05, plot='vtk')