Issue loading Shell Element Models

[johonea4]

Hello. I'm having an issue when loading RST files with dpf post, and I believe the issue persists in dpf-core as well. When loading a verification test case RST file generated from Ansys 19 (Have checked prior versions with similar behavior) parts of the model are missing when doing a mesh.plot(). Below are images of what I'm expecting vs. what I get. Any insight? I am able to share the RST as well if wanted.

^^^ Generated from reading RST with in-house C/C++ reader

^^^ DPF Post mesh.plot result

[johonea4]

Additional Info: I opened the same file with the legacy ansys-mapdl-reader and the plotted result matches expectation (1st image above). I haven't been able to find any cause of why dpf-post is producing the 2nd image above. I will also mention that the node counts match in all readers, but dpf-post's element count is reduced by the number of missing elements. It's strange.

pyansys-reader:

\>>>reader.read_binary("v19_104.rst")
PyMAPDL Result
Units       : User Defined
Version     : 19.0
Cyclic      : False
Result Sets : 10
Nodes       : 173
Elements    : 96

Available Results:
EMS : Miscellaneous summable items (normally includes face pressures)
ENF : Nodal forces
ENS : Nodal stresses
ENG : Element energies and volume
EEL : Nodal elastic strains
ETH : Nodal thermal strains (includes swelling strains)
EUL : Element euler angles
EMN : Miscellaneous nonsummable items
EPT : Nodal temperatures
NSL : Nodal displacements
RF  : Nodal reaction forces

dpf-post

\>>>sol = post.load_solution("v19_104.rst")
\>>> print(sol)
Modal solution object.

Data Sources
------------------------------
DPF  DataSources:
  Result files:
     result key: rst and path: \\devsrv02\public\Models\EDAS\ANSYS\v190\v19_104.rst
  Secondary files:

DPF Model
------------------------------
DPF Result Info
  Analysis: modal
  Physics Type: mecanic
  Unit system: SI: m, kg, N, s, V, A, K
  Available results:
    U Displacement :nodal displacements
    RF Force :nodal reaction forces
    ENF Element nodal Forces :element nodal forces
    S Stress :element nodal component stresses
    ENG_VOL Volume :element volume
    ENG_SE Energy-stiffness matrix :element energy associated with the stiffness matrix
    ENG_AHO Hourglass Energy :artificial hourglass energy
    ENG_TH thermal dissipation energy :thermal dissipation energy
    ENG_KE Kinetic Energy :kinetic energy
    ENG_CO co-energy :co-energy (magnetics)
    ENG_INC incremental energy :incremental energy (magnetics)
    EPEL Strain :element nodal component elastic strains
    ETH Thermal Strains :element nodal component thermal strains
    ETH_EQV Thermal Strains eqv :element nodal equivalent component thermal strains
    ETH_SWL Swelling Strains :element nodal swelling strains
    BFE Temperature :element structural nodal temperatures
\------------------------------
DPF  Meshed Region:
  173 nodes
  72 elements
  Unit:
  With shell (2D) elements
\------------------------------
DPF  Time/Freq Support:
  Number of sets: 10
Cumulative     Frequency (Hz) LoadStep       Substep
1              11.972045      1              1
2              18.994242      1              2
3              40.068307      1              3
4              45.658840      1              4
5              47.868608      1              5
6              53.224182      1              6
7              69.962842      1              7
8              86.567796      1              8
9              91.567291      1              9
10             99.643575      1              10

[rlagha]

Hi @johonea4 , Thanks for finding this, please could you share the rst file? Thanks Ramdane

[johonea4]

v190.zip

Here are 2 sets of RST/CDB files this happens on.

[johonea4]

I can certainly supply more if needed. Please let me know.

[anslpa]

Hi @johonea4, The issue was that the shell43 and shell63 elements (used in the provided file) were not supported. A fix is found and will be available for the next Ansys 2022 R1 release. @rlagha Thank you.