Problem while creation of multiple volumes by VGEN command

[herttim]

Dear all,

I am struggling a bit with the creation of multiple volumes by the VGEN command.
While running the code, a failure message is shown in the console. My investigations led me to the section in which the ‘rollers’ are generated with the VGEN command around the axis of rotation of the roller bearing. This creates two rows of rollers: Inboard and Outboard.
As soon as I try to look at the volume, I get an error message (see attached).

Maybe you can tell me why I can't generate the volumes. Thank you very much! :)

General setup:

import matplotlib.pyplot as plt
import numpy as np
import math as m
import itertools
import pandas

from ansys.mapdl.core import launch_mapdl

mapdl = launch_mapdl()

mapdl.prep7()

mapdl.units("SI")

mapdl.ignore_errors = True

mapdl.et(1, "SOLID185")
mapdl.mp("EX", 1, 210e9)
mapdl.mp("DENS", 1, 7800)
mapdl.mp("NUXY", 1, 0.3)
mapdl.keyopt(1, 2, 3)

mapdl.nlgeom("ON")

#Slicing
anzahl_scheiben = 5
anzahl_scheiben_board = 20`

[...]

**Definition of Rollers:**
`#ROLLEN
#Outboard
#Keypoints
mapdl.k(60, 0, 15.11e-03+of, 8.74e-03)
mapdl.k(61, 0, 21.61e-03+of, 10.21e-03)
#mapdl.k(62, 0, 21.81e-03+of, 10.25e-03)
#mapdl.k(63, 0, 15.04e-03+of, 37.07e-03)
mapdl.k(64, 0, 14.9e-03+of, 37.04e-03)
mapdl.k(65, 0, 9.02e-03+of, 35.71e-03)

#Linien
r_ob_l60 = mapdl.l(60, 61)
r_ob_l61 = mapdl.l(61, 64) #Laufbahn Rolle - Außenring Outboard
r_ob_l62 = mapdl.l(64, 65)
r_ob_l63 = mapdl.l(65, 60) 

#Fläche Rolle Outboard
a_r_ob = mapdl.a(60, 61, 64, 65)

#Erzeugung Scheiben Rollen Outboard
mapdl.lesize(nl1=r_ob_l61, ndiv=anzahl_scheiben-1) #Outboard

#Inboard
#Keypoints
mapdl.k(70, 0, 9.02e-03+of, 79.29e-03)
mapdl.k(71, 0, 14.9e-03+of, 77.96e-03)
#mapdl.k(72, 0, 15.04e-03+of, 77.93e-03)
#mapdl.k(73, 0, 21.8e-03+of, 104.75e-03)
mapdl.k(74, 0, 21.61e-03+of, 104.79e-03)
mapdl.k(75, 0, 15.11e-03+of, 106.26e-03)

#Linien
r_ib_l70 = mapdl.l(70, 71)
r_ib_l71 = mapdl.l(71, 74) #Laufbahn Rolle - Außenring Inboard
r_ib_l72 = mapdl.l(74, 75) 
r_ib_l73 = mapdl.l(75, 70)

#Fläche Rolle Inboard
a_r_ib = mapdl.a(70, 71, 74, 75)

print("Flächen vorhanden für Outboard?", mapdl.geometry.areas)
print("Flächen vorhanden für Inboard?", mapdl.geometry.areas)
print("a_r_ob:", a_r_ob)
print("a_r_ib:", a_r_ib)

#Erzeugung Scheiben Rollen Inboard
mapdl.lesize(nl1=r_ib_l71, ndiv=anzahl_scheiben-1) #Inboard

#EXTRUDIEREN - Rollen 

#Outboard
mapdl.vsel("NONE")
mapdl.vrotat(a_r_ob, pax1=60, pax2=65, arc=360, nseg=2) #Outboard
print("Outboard:",mapdl.geometry.volumes)

mapdl.csys(1)
mapdl.vgen(itime=anzahl_rollen, nv1="ALL", dy=360/anzahl_rollen, imove=0)
mapdl.vsel("ALL")

#Inboard
mapdl.vsel("NONE")
mapdl.vrotat(a_r_ib, pax1=70, pax2=75, arc=360, nseg=2) #Inboard
print("Inboard:",mapdl.geometry.volumes)

mapdl.csys(1)
mapdl.vgen(itime=anzahl_rollen, nv1="ALL", dy=360/anzahl_rollen, imove=0)

mapdl.vsel("ALL")

mapdl.vplot()

mapdl.exit()`

Output from Console:

`C:\Users\HertingerT\Anaconda\lib\site-packages\pyvista\core\pointset.py:843: PyVistaDeprecationWarning: `PolyData` constructor parameter `n_faces` is deprecated and no longer used.
  warnings.warn(
Flächen vorhanden für Outboard? Multiblock (0x19c1dae5600)
  N Blocks    849
  X Bounds    -0.069, 0.069
  Y Bounds    -0.069, 0.069
  Z Bounds    0.000, 0.120
Flächen vorhanden für Inboard? Multiblock (0x19c1f1a80a0)
  N Blocks    849
  X Bounds    -0.069, 0.069
  Y Bounds    -0.069, 0.069
  Z Bounds    0.000, 0.120
a_r_ob: 849
a_r_ib: 850


Windows fatal exception: access violation


Thread 0x00003f50 (most recent call first):
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\launcher.py", line 721 in enqueue_output
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 953 in run
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 1016 in _bootstrap_inner
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 973 in _bootstrap



Thread 0x00001760 (most recent call first):
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\launcher.py", line 721 in enqueue_output
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 953 in run
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 1016 in _bootstrap_inner
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 973 in _bootstrap



Thread 0x00005cf4 (most recent call first):
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\launcher.py", line 721 in enqueue_output
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 953 in run
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 1016 in _bootstrap_inner
  File "C:\Users\HertingerT\Anaconda\lib\threading.py", line 973 in _bootstrap



Main thread:
Current thread 0x00001a70 (most recent call first):
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\pyvista\core\filters\__init__.py", line 40 in _update_alg
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\pyvista\core\filters\data_set.py", line 5044 in extract_cells
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\mapdl_geometry.py", line 1310 in get_volumes
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\misc.py", line 1121 in wrapper
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\mapdl_geometry.py", line 1227 in volumes
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\ansys\mapdl\core\misc.py", line 1121 in wrapper
  File "c:\users\hertingert\python programm\trb\20250227_trb_model_v5-2.py", line 247 in <module>
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\spyder_kernels\py3compat.py", line 356 in compat_exec
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 469 in exec_code
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 611 in _exec_file
  File "C:\Users\HertingerT\Anaconda\lib\site-packages\spyder_kernels\customize\spydercustomize.py", line 524 in runfile
  File "C:\Users\HertingerT\AppData\Local\Temp\ipykernel_20684\152702173.py", line 1 in <module>


Restarting kernel...

[germa89]

Hi @herttim

I am missing some parameters like of and anzahl_rollen.

But setting those to 0 and 2 respectively, I get this plot:

However, I believe your problem is related to plotting (many Windows fatal exception: access violation are). Are you on latest PyMAPDL version?

Code used to get VPLOT

import matplotlib.pyplot as plt
import numpy as np
import math as m
import itertools
import pandas

from ansys.mapdl.core import launch_mapdl

mapdl = launch_mapdl()


mapdl.clear()
mapdl.prep7()

mapdl.units("SI")

mapdl.ignore_errors = True

mapdl.et(1, "SOLID185")
mapdl.mp("EX", 1, 210e9)
mapdl.mp("DENS", 1, 7800)
mapdl.mp("NUXY", 1, 0.3)
mapdl.keyopt(1, 2, 3)

mapdl.nlgeom("ON")

#Slicing
anzahl_scheiben = 5
anzahl_scheiben_board = 20

# [...]

# **Definition of Rollers:**
#ROLLEN
#Outboard
#Keypoints
of = 0
mapdl.k(60, 0, 15.11e-03+of, 8.74e-03)
mapdl.k(61, 0, 21.61e-03+of, 10.21e-03)
#mapdl.k(62, 0, 21.81e-03+of, 10.25e-03)
#mapdl.k(63, 0, 15.04e-03+of, 37.07e-03)
mapdl.k(64, 0, 14.9e-03+of, 37.04e-03)
mapdl.k(65, 0, 9.02e-03+of, 35.71e-03)

#Linien
r_ob_l60 = mapdl.l(60, 61)
r_ob_l61 = mapdl.l(61, 64) #Laufbahn Rolle - Außenring Outboard
r_ob_l62 = mapdl.l(64, 65)
r_ob_l63 = mapdl.l(65, 60) 

#Fläche Rolle Outboard
a_r_ob = mapdl.a(60, 61, 64, 65)

#Erzeugung Scheiben Rollen Outboard
mapdl.lesize(nl1=r_ob_l61, ndiv=anzahl_scheiben-1) #Outboard

#Inboard
#Keypoints
mapdl.k(70, 0, 9.02e-03+of, 79.29e-03)
mapdl.k(71, 0, 14.9e-03+of, 77.96e-03)
#mapdl.k(72, 0, 15.04e-03+of, 77.93e-03)
#mapdl.k(73, 0, 21.8e-03+of, 104.75e-03)
mapdl.k(74, 0, 21.61e-03+of, 104.79e-03)
mapdl.k(75, 0, 15.11e-03+of, 106.26e-03)

#Linien
r_ib_l70 = mapdl.l(70, 71)
r_ib_l71 = mapdl.l(71, 74) #Laufbahn Rolle - Außenring Inboard
r_ib_l72 = mapdl.l(74, 75) 
r_ib_l73 = mapdl.l(75, 70)

#Fläche Rolle Inboard
a_r_ib = mapdl.a(70, 71, 74, 75)

print("Flächen vorhanden für Outboard?", mapdl.geometry.areas)
print("Flächen vorhanden für Inboard?", mapdl.geometry.areas)
print("a_r_ob:", a_r_ob)
print("a_r_ib:", a_r_ib)

#Erzeugung Scheiben Rollen Inboard
mapdl.lesize(nl1=r_ib_l71, ndiv=anzahl_scheiben-1) #Inboard

#EXTRUDIEREN - Rollen 

#Outboard
mapdl.vsel("NONE")
mapdl.vrotat(a_r_ob, pax1=60, pax2=65, arc=360, nseg=2) #Outboard
print("Outboard:",mapdl.geometry.volumes)

anzahl_rollen= 2

mapdl.csys(1)
mapdl.vgen(itime=anzahl_rollen, nv1="ALL", dy=360/anzahl_rollen, imove=0)
mapdl.vsel("ALL")

#Inboard
mapdl.vsel("NONE")
mapdl.vrotat(a_r_ib, pax1=70, pax2=75, arc=360, nseg=2) #Inboard
print("Inboard:",mapdl.geometry.volumes)

mapdl.csys(1)
mapdl.vgen(itime=anzahl_rollen, nv1="ALL", dy=360/anzahl_rollen, imove=0)

mapdl.vsel("ALL")

mapdl.vplot()

# mapdl.exit()