CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   ParaView (http://www.cfd-online.com/Forums/paraview/)
-   -   Extract Paraview data (vtkUnstructured Grid) using pvpython (http://www.cfd-online.com/Forums/paraview/125648-extract-paraview-data-vtkunstructured-grid-using-pvpython.html)

lu_siyah October 28, 2013 14:02

Extract Paraview data (vtkUnstructured Grid) using pvpython
 
I am working with large FE simulation data and I wish to extract information on only coordinates (at different timesteps, precisely a 3xN List of coordinates for every timestep considered or alternatively the displacements, which I know to be stored in a Point Array) and the basic connectivity of the grid into python arrays.
For now, I am working with Paraview in the pvpython-shell and vtk and struggling with finding and arranging the information. So far I used this code to get the coordinates in the beginning (data considered originates from LS-Dyna, hence 3dplot):
Code:

import numpy as np
from paraview import simple 
reader=simple.OpenDataFile('/home/data/d3plot')
reader.PartArrays = ['Part1 (Matl2000001)']
reader.UpdatePipeline()
# reader: <paraview.servermanager.LSDynaReader object at 0x376b5d0> 
data = simple.servermanager.Fetch(reader)
# data: (vtkMultiBlockDataSet)
grid = data.GetBlock(0)
# grid: (vtkUnstructuredGrid)
npoints = grid.GetNumberOfPoints()
coord = np.zeros((npoints,3))
for i in range(npoints):
  coord[i,:] = np.array(grid.GetPoint(i))

Any ideas on how to proceed from there? That is: Get a matrix of connectivities, and get acess to the data saved in the Point Arrays OR the coordinates at different time steps??
Any help will be appreciated as I am fairly new to this sort of programming.

wyldckat November 3, 2013 06:10

Greetings lu_siyah and welcome to the forum!

Sorry, I'm a bit confused with your description. Is what you're trying to do similar to what's described on this thread: http://www.cfd-online.com/Forums/par...-paraview.html :confused:

Because if I understand you correctly, you want to resample/probe the simulation volume you have from the FE simulation, with an array of point/probe locations.
The part that I don't understand is: why do you need a connectivity matrix?

Best regards,
Bruno


All times are GMT -4. The time now is 10:16.