PyMAPDL Reader - Legacy Binary and Archive File Reader
This is the legacy module for reading in binary and ASCII files generated from MAPDL.
This Python module allows you to extract data directly from binary ANSYS v14.5+ files and to display or animate them rapidly using a straightforward API coupled with C libraries based on header files provided by ANSYS.
The ansys-mapdl-reader module supports the following formats:
*.rst - Structural analysis result file
*.rth - Thermal analysis result file
*.emat - Element matrix data file
*.full - Full stiffness-mass matrix file
*.cdb or *.dat - MAPDL ASCII block archive and Mechanical Workbench input files
Please see the PyMAPDL-Reader Documentation for the full documentation.
Installation
Installation through pip:
pip install ansys-mapdl-reader
You can also visit pymapdl-reader to download the source or releases from GitHub.
Python 3.10 Extra Instructions
PyMAPDL-Reader requires the VTK library which, at the moment, is not available for Python 3.10 in their official channel.
If you wish to install PyMAPDL-Reader in Python 3.10, you can still do it by using the unofficial VTK wheel from PyVista using --find-links. This tells pip to look for vtk at wheels.pyvista.org. Use this with:
pip install ansys-mapdl-reader --find-links https://wheels.pyvista.org/
Please visit Unofficial VTK Wheels for Python 3.10 for further details.
Examples
Loading and Plotting a MAPDL Archive File
ANSYS archive files containing solid elements (both legacy and modern), can be loaded using Archive and then converted to a vtk object.
from ansys.mapdl import reader as pymapdl_reader
from ansys.mapdl.reader import examples
# Sample *.cdb
filename = examples.hexarchivefile
# Read ansys archive file
archive = pyansys.Archive(filename)
# Print raw data from cdb
for key in archive.raw:
print("%s : %s" % (key, archive.raw[key]))
# Create a vtk unstructured grid from the raw data and plot it
grid = archive.parse_vtk(force_linear=True)
grid.plot(color='w', show_edges=True)
# write this as a vtk xml file
grid.save('hex.vtu')
# or as a vtk binary
grid.save('hex.vtk')
You can then load this vtk file using pyvista or another program that uses VTK.
# Load this from vtk
import pyvista as pv
grid = pv.UnstructuredGrid('hex.vtu')
grid.plot()Loading the Result File
This example reads in binary results from a modal analysis of a beam from ANSYS.
# Load the reader from pyansys
from ansys.mapdl import reader as pymapdl_reader
from ansys.mapdl.reader import examples
# Sample result file
rstfile = examples.rstfile
# Create result object by loading the result file
result = pyansys.read_binary(rstfile)
# Beam natural frequencies
freqs = result.time_values>>> print(freq)
[ 7366.49503969 7366.49503969 11504.89523664 17285.70459456
17285.70459457 20137.19299035]Get the 1st bending mode shape. Results are ordered based on the sorted node numbering. Note that results are zero indexed
>>> nnum, disp = result.nodal_solution(0)
>>> print(disp)
[[ 2.89623914e+01 -2.82480489e+01 -3.09226692e-01]
[ 2.89489249e+01 -2.82342416e+01 2.47536161e+01]
[ 2.89177130e+01 -2.82745126e+01 6.05151053e+00]
[ 2.88715048e+01 -2.82764960e+01 1.22913304e+01]
[ 2.89221536e+01 -2.82479511e+01 1.84965333e+01]
[ 2.89623914e+01 -2.82480489e+01 3.09226692e-01]
...Plotting Nodal Results
As the geometry of the model is contained within the result file, you can plot the result without having to load any additional geometry. Below, displacement for the first mode of the modal analysis beam is plotted using VTK.
# Plot the displacement of Mode 0 in the x direction
result.plot_nodal_solution(0, 'x', label='Displacement')
Results can be plotted non-interactively and screenshots saved by setting up the camera and saving the result. This can help with the visualization and post-processing of a batch result.
First, get the camera position from an interactive plot:
>>> cpos = result.plot_nodal_solution(0)
>>> print(cpos)
[(5.2722879880979345, 4.308737919176047, 10.467694436036483),
(0.5, 0.5, 2.5),
(-0.2565529433509593, 0.9227952809887077, -0.28745339908049733)]Then generate the plot:
result.plot_nodal_solution(0, 'x', label='Displacement', cpos=cpos,
screenshot='hexbeam_disp.png',
window_size=[800, 600], interactive=False)Stress can be plotted as well using the below code. The nodal stress is computed in the same manner that ANSYS uses by to determine the stress at each node by averaging the stress evaluated at that node for all attached elements. For now, only component stresses can be displayed.
# Display node averaged stress in x direction for result 6
result.plot_nodal_stress(5, 'Sx')
Nodal stress can also be generated non-interactively with:
result.plot_nodal_stress(5, 'Sx', cpos=cpos, screenshot=beam_stress.png,
window_size=[800, 600], interactive=False)Animating a Modal Solution
Mode shapes from a modal analysis can be animated using animate_nodal_solution:
result.animate_nodal_solution(0)If you wish to save the animation to a file, specify the movie_filename and animate it with:
result.animate_nodal_solution(0, movie_filename='/tmp/movie.mp4', cpos=cpos)
Reading a Full File
This example reads in the mass and stiffness matrices associated with the above example.
# Load the reader from pyansys
from ansys.mapdl import reader as pymapdl_reader
from scipy import sparse
# load the full file
fobj = pyansys.FullReader('file.full')
dofref, k, m = fobj.load_km() # returns upper triangle only
# make k, m full, symmetric matrices
k += sparse.triu(k, 1).T
m += sparse.triu(m, 1).TIf you have scipy installed, you can solve the eigensystem for its natural frequencies and mode shapes.
from scipy.sparse import linalg
# condition the k matrix
# to avoid getting the "Factor is exactly singular" error
k += sparse.diags(np.random.random(k.shape[0])/1E20, shape=k.shape)
# Solve
w, v = linalg.eigsh(k, k=20, M=m, sigma=10000)
# System natural frequencies
f = np.real(w)**0.5/(2*np.pi)
print('First four natural frequencies')
for i in range(4):
print '{:.3f} Hz'.format(f[i])First four natural frequencies
1283.200 Hz
1283.200 Hz
5781.975 Hz
6919.399 HzDeveloping on Windows
This package is designed to be developed on Linux, and if you need to develop on Windows you will need to install your own C++ compiler. We recommend:
- Install Visual C++
See here for a list of which Python versions correspond to which Visual C++ version
Only Python <= 3.8 appears to be supported at the moment.
- Install the development version of pymapdl-reader to your Python environment
Navigate to the project’s top level (the same directory as this README)
run pip install -e .
To get the package up and running.
License and Acknowledgments
The ansys-mapdl-reader module is licensed under the MIT license.
1 maintainer
Metadata
1 maintainer
Release files for ansys-mapdl-reader 0.51.10
For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.
Built distributions (wheels)
Total release size: 19.4 MB
Release files / ansys_mapdl_reader-0.51.10-cp310-cp310-win_amd64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp310-cp310-win_amd64.whl |
|---|---|
| Size | 1.1 MB |
| Tags | CPython 3.10 Windows x86-64 |
|
SHA-256 checksum How to use checksums |
6165ee46693790d3917f15f3b0893d8b0208013b7ecb5f93e8cd5217166d8870
|
|
BLAKE2b-256 checksum How to use checksums |
58d574a409de70c7d0b793934eb51e5ba5ec8656c139d6e48e2868fc5f74419d
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl |
|---|---|
| Size | 3.8 MB |
| Tags | CPython 3.10 Linux glibc 2.17+ x86-64 |
|
SHA-256 checksum How to use checksums |
f02645586365477083af6f2f46e81d96c7a10f1afac19c34a211e0d674f7081f
|
|
BLAKE2b-256 checksum How to use checksums |
1566fc6dea7fc038867d8154492eed426d73240683139c76ae2b44e967914514
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp39-cp39-win_amd64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp39-cp39-win_amd64.whl |
|---|---|
| Size | 1.1 MB |
| Tags | CPython 3.9 Windows x86-64 |
|
SHA-256 checksum How to use checksums |
6fc00f36948ae3cfbe1132242643dc15cd87c737528d163a1cbcfac1c98d6987
|
|
BLAKE2b-256 checksum How to use checksums |
9b9597fcbbe459f0c3fef3ba2c2e9a5792978f8fef69b3d63c25a94d676eef6e
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp39-cp39-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp39-cp39-manylinux_2_17_x86_64.manylinux2014_x86_64.whl |
|---|---|
| Size | 3.8 MB |
| Tags | CPython 3.9 Linux glibc 2.17+ x86-64 |
|
SHA-256 checksum How to use checksums |
7b5f825977339af246eadc5849ec4d32c9fd68a7adbd049fee816517353f8c4a
|
|
BLAKE2b-256 checksum How to use checksums |
1e505f6a4983011002db030b41366d848af9ad3192dd8f65b3b4b162d7469240
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp38-cp38-win_amd64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp38-cp38-win_amd64.whl |
|---|---|
| Size | 1.1 MB |
| Tags | CPython 3.8 Windows x86-64 |
|
SHA-256 checksum How to use checksums |
f214b721c89c2177f089f7a08167d956f660e58d66d1633911ee42a67cc0e118
|
|
BLAKE2b-256 checksum How to use checksums |
2d82945a6f894313f948c7f8efc8838d2e9f93827e35b8f771ad0e121a48de89
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp38-cp38-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp38-cp38-manylinux_2_17_x86_64.manylinux2014_x86_64.whl |
|---|---|
| Size | 3.9 MB |
| Tags | CPython 3.8 Linux glibc 2.17+ x86-64 |
|
SHA-256 checksum How to use checksums |
84d7bd4308401a717924e399a17fab0b5898ac4e1ee40899a26f5ae87803add7
|
|
BLAKE2b-256 checksum How to use checksums |
718f0f01a50a899c6338ab765675995fa048be8be5729f004a0702f0283beac1
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp37-cp37m-win_amd64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp37-cp37m-win_amd64.whl |
|---|---|
| Size | 1.1 MB |
| Tags | CPython 3.7 CPython 3.7 pymalloc Windows x86-64 |
|
SHA-256 checksum How to use checksums |
b4b7518b3a0b65e198530113144cccd77077be0cacf57af447b67c13b6110e32
|
|
BLAKE2b-256 checksum How to use checksums |
6859d2ca32d1d20b13b94ef00233faf957eb5c3294eca9f95b40446849d809f5
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
Release files / ansys_mapdl_reader-0.51.10-cp37-cp37m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl
| Download URL | ansys_mapdl_reader-0.51.10-cp37-cp37m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl |
|---|---|
| Size | 3.6 MB |
| Tags | CPython 3.7 CPython 3.7 pymalloc Linux glibc 2.17+ x86-64 |
|
SHA-256 checksum How to use checksums |
6ee95055715f2d105f6ccebc04cb11d1b6f73b8a86e9aed4fa2df112f4e1e6a4
|
|
BLAKE2b-256 checksum How to use checksums |
260521bdd591874ac67f9a9b28975ce8e313665e68b3b9358a6183248793a1ce
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/3.8.0 pkginfo/1.8.2 readme-renderer/32.0 requests/2.27.1 requests-toolbelt/0.9.1 urllib3/1.26.8 tqdm/4.62.3 importlib-metadata/4.11.1 keyring/23.5.0 rfc3986/2.0.0 colorama/0.4.4 CPython/3.9.10
|
