from abapy.mesh import Mesh, Nodes
mesh = Mesh()
nodes = mesh.nodes
# Adding some nodes
nodes.add_node(label=1, x=0., y=0., z=0.)
nodes.add_node(label=2, x=1., y=0., z=0.)
nodes.add_node(label=3, x=1., y=1., z=0.)
nodes.add_node(label=4, x=0., y=1., z=0.)
nodes.add_node(label=5, x=2., y=0., z=0.)
nodes.add_node(label=6, x=2., y=1., z=0.)
# Adding some elements
mesh.add_element(label=1,
connectivity=(1, 2, 3, 4),
space=2,
name='QUAD4',
toset='mySet')
mesh.add_element(label=2,
connectivity=(2, 5, 6, 3),
space=2,
name='QUAD4',
toset=['mySet', 'myOtherSet'])
print(mesh)
from abapy.postproc import FieldOutput
from abapy.mesh import Mesh, Nodes
x = [0., 1., 0.]
y = [0., 0., 1.]
z = [0., 0., 0.]
labels = [1, 2, 3]
nodes = Nodes(x=x, y=y, z=z, labels=labels)
mesh = Mesh(nodes=nodes)
mesh.add_element(label=1, connectivity=[1, 2, 3], space=2,
name='tri3') # triangle element
nodeField = FieldOutput()
nodeField.add_data(data=0., label=1)
nodeField.add_data(data=10., label=2)
nodeField.add_data(data=20., label=3)
elementField = FieldOutput(position='element')
elementField.add_data(label=1, data=10.)
out = ''
out += mesh.dump2vtk()
out += nodeField.dump2vtk('nodeField')
out += elementField.dump2vtk('elementField')
f = open("FieldOutput-dump2vtk.vtk", "w")
f.write(out)
f.close()
