def setUp(self): self.pts = [[0.,0.,0.], [1.,0.,0.], [0.,1.,0.], [0.,0.,1], [0.5,0.,0.], [0.5, 0.5, 0.], [0., 0.5, 0.], [0., 0., 0.5], [0.5, 0., 0.5], [0., 0.5, 0.5]] self.tet = quad_tet(1, self.pts, range(10)) self.vtk = create_vtk([-1.,0,1]) self.vtk2 = create_vtk([-1,0,1],[5]*27) self.vtk3 = create_vtk([0,1,2],[0, 5, 10]*9)
def create_quad_tet_part(): p = part('Part1','C3D10','quad_tet') pts = [[-1., 1., -1.], [-1., -1., 1.], [1., -1., -1.], [1., 1. , 1.], [-1., 0., 0.], [0., -1., 0.], [0., 0., -1.], [0., 1., 0.], [0., 0., 1.], [1., 0., 0.]] nodes = [repr(n+1) for n in range(10)] p.add_element(quad_tet(1, pts, nodes)) p.moduli = [1.] return [p]