Errors with d3plot sph_* and part_* outputs. #67
Description
🐛 Describe the bug
Multiple issues with d3plot
- SPH data being read incorrectly.
- NRBs are read as parts
🔢 To Reproduce
1.
1.1. Read d3plot (d3plot = D3plot(filepath))
1.2. Read sph data: arr = d3plot.arrays['sph_mass']
arr.shape will be (N, ) with N being the number of states.
The same holds for many other SPH outputs (sph_n_neighbors, sph_density, ...)
sph_n_neighbors actually returns the density of individual particles. sph_density and sph_mass return some incorrect arrays (sometimes with negative values)
2.1. Read d3plot (d3plot = D3plot(filepath))
2.2. parts = d3plot.arrays['part_ids']
parts will include ids of *CONSTRAINED_NODAL_RIGID_BODY
The same happens with part_kinetic_energy, where kinetic energy of NRBs is included.
💘 Expected behavior
1.
For sph data, return the overall correct sph_mass, density etc.. or
return an array (N, X) where N is the number of states and X is the number of particles
2.
For parts, exclude nrbs from part_* outputs.
🖥️ Setup
- lasso-python version: [2.0.3]
- OS: [Windows 11]
- LS-Dyna R14.0.0