Extract eigen values and eigen vectors from d3eigv files in pydpf

Adya

Extract eigen values and eigen vectors from d3eigv files in pydpf

Adya

from ansys.dpf import core as dpf
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
from pxr import Vt

# Load the datasource from a d3eigv file
# only the first d3eigv file need to be loaded
ds = dpf.DataSources()
ds.set_result_file_path(r'\path\to\file\d3eigv', 'd3plot')
model = dpf.Model(ds)
print(model)


# Extract displacements for all Load steps from d3eigv
U = model.results.displacement.on_all_time_freqs.eval()

# Get the displacement field for a desired loadstep
# Here 7 refers  to the 7th loadstep in d3eigv
print(U)
u = U.get_field({'time': 6})

"""
# Plot the extracted state
# Get the camera position once so as to set the view
cpos = u.plot(deform_by=u,show_edges=False,
       screenshot=r'\path\to\image\state7.png',
       off_screen=False, return_cpos=True)
print(cpos)
"""

cpos = [(-8739.018807515353, 4113.571013794808, 2869.855509287625),
 (-2424.0205557548243, -648.3712183048058, 229.55202115622455),
 (0.23537139638660276, -0.21357474711933727, 0.9481487927303034)]

cpos = u.plot(deform_by=u,show_edges=False,
       screenshot=r'\path\to\image\state7.png',
       off_screen=True, return_cpos=False,cpos=cpos)
# Extract the x,y,z displacements from this loadstep
disp = u.data

# Convert that to a pandas dataframe of numpy ararys
df = pd.DataFrame({'NodID':np.array(u.scoping.ids),
                   'x-disp':disp[:,0],'y-disp':disp[:,1],
                   'z-disp':disp[:,2]})

print(df.head())

# Convert the displacements to VtArray
def convert_np_to_vt(my_array: np.ndarray) -> Vt.Vec3fArray:
    return Vt.Vec3fArray.FromNumpy(my_array)

from_numpy: Vt.Vec3fArray = convert_np_to_vt(disp)
print(from_numpy[:5])