Inputted .anf file has no elements

[ay8tb]

Hi,
I'm facing an interesting issue where any .anf file with a void (whether full or partial) either fails to mesh (in the case of a full void) or fails to have any elements. I don't have this issue with other .anf files that have no voids or cuts. SYS-30 is an example of one that fails. Geom_trial_mm is an example of one that doesn't.
Geom_trial_mm.zip
SYS-30.zip
UserWarning: There are no elements to plot. warnings.warn("There are no elements to plot.")

As far as I can tell when opening it in Notepad, the file isn't empty and is exported the same way through Spaceclaim.

[mikerife]

Hi @ay8tb
SYS-30.anf did not fail to open and be volume plotted for me. Before I forget you may want to open the anf in a text editor and comment out this first /out line (with a !) just in case. Since I did not have that exact path on my system it threw an error, and a Fortran error!, but continued on.

Here is the script I used:

from ansys.mapdl.core import launch_mapdl
import os
path = os.getcwd()
mapdl = launch_mapdl(run_location=path, additional_switches='-smp')
print(mapdl)
input = mapdl.input('SYS-30.anf')
mapdl.prep7()
mapdl.vplot()

I'm running MAPDL 2023R2 and PyMAPDL 0.68.1.
Mike

[ay8tb]

Hi @mikerife,
I could see that on my end thanks. I'm still having an issue with being unable to see the element plot. Here's the code:
mapdl.esize(0.5) mapdl.et(1,185) mapdl.vsweep( "ALL") mapdl.eplot(color="grey", background='b') mapdl.eplot(vtk=True, cpos="xy", show_edges=True, show_axes=True, plot_bc=True, plot_bc_legend=True) mapdl.show("png") mapdl.eplot(vtk=False) mapdl.show("close") mapdl.esize() mapdl.et(1, "SOLID185") mapdl.vsweep("all") mapdl.eplot(vtk=True, show_edges=True, show_axes=False, line_width=2, background="w")
And this is the error:
C:\Users\####\anaconda3\envs\####\lib\site-packages\ansys\mapdl\core\mapdl_extended.py:1219: UserWarning: There are no elements to plot. warnings.warn("There are no elements to plot.") Traceback (most recent call last): File "C:\Program Files\JetBrains\PyCharm 2022.3\plugins\python\helpers\pydev\pydevd.py", line 1535, in _exec pydev_imports.execfile(file, globals, locals) # execute the script File "C:\Program Files\JetBrains\PyCharm 2022.3\plugins\python\helpers\pydev\_pydev_imps\_pydev_execfile.py", line 18, in execfile exec(compile(contents+"\n", file, 'exec'), glob, loc) File "C:\Users\####\PycharmProjects\####\####.py", line 50, in <module> mapdl.eplot(vtk=True, cpos="xy", show_edges=True, show_axes=True, plot_bc=True, plot_bc_legend=True) File "C:\Users\####\anaconda3\envs\####\lib\site-packages\ansys\mapdl\core\mapdl_extended.py", line 1220, in eplot return general_plotter([], [], [], mapdl=self, **kwargs) File "C:\Users####ayehi\anaconda3\envs\####\lib\site-packages\ansys\mapdl\core\plotting.py", line 921, in general_plotter pl = bc_plotter( File "C:\Users\####\anaconda3\envs\####\lib\site-packages\ansys\mapdl\core\plotting.py", line 997, in bc_plotter bc_glyph_size = get_bounding_box(mapdl.mesh.nodes) File "C:\Users\####\anaconda3\envs\####\lib\site-packages\ansys\mapdl\core\misc.py", line 551, in get_bounding_box min_ = np.min(nodes_xyz, axis=0) File "<__array_function__ internals>", line 180, in amin File "C:\Users\####\anaconda3\envs\####\lib\site-packages\numpy\core\fromnumeric.py", line 2916, in amin return _wrapreduction(a, np.minimum, 'min', axis, None, out, File "C:\Users\####\anaconda3\envs\####\lib\site-packages\numpy\core\fromnumeric.py", line 86, in _wrapreduction return ufunc.reduce(obj, axis, dtype, out, **passkwargs) ValueError: zero-size array to reduction operation minimum which has no identity

I'm running MAPDL 2022R2 and PyMAPDL 0.68.dev0