How can I get Rotational velocity boundary condition from the RST?

Find more posts tagged with

PyMAPDL

Comments

Using the*XPL command from PyMAPDL we can get the omegacgvalues (position 20, 21 and 22 in the record) as shown in the Results Description File.

Ansys Help Link

from ansys.mapdl.core import launch_mapdl

mapdl = launch_mapdl()

rst = r"D:\…\file.rst"

mapdl.post1()
mapdl.file(rst)

xpl = mapdl.xpl
xpl.open(rst)
xpl.step('DSI')

xpl.step("SET1")
step_data = xpl.read("DPHEAD")
omega_x = step_data[19]
omega_y = step_data[20]
omega_z = step_data[21]
print(omega_x, omega_y, omega_z)
xpl.close()
mapdl.exit()

Naveen Kumar Begari

@Ayush Kumar Hi…

Is it possible to get all the Boundary Conditions from .rst using PyDPF.?

Jim Kosloski

Thanks for posting this Ayush. This has always been a struggle. I wanted to point out a few things for clarification:

You are getting the CGOMEGA values. This is what you want if you are using Mechanical because Mechanical writes the Rotational Velocity as a CGOMEGA command, but if you were running MAPDL it is likely that you would want positions 13, 14 and 15 to get the OMEGA values
You issue a mapdl.set('last') which I am not sure is necessary, but it is a little confusing, because in the XPL command you access DSI::SET1 this is actually the first set on the RST not the last. You can use DSI::SET# to access the results set desired.

Since this uses the *XPL command, this can be done directly in MAPDL as well. I have attached a macro that does that.

getomega.txt

Ayush Kumar

@Naveen Kumar Begari no, you can't. The above code is using PyMAPDL.

Ayush Kumar

@Jim Kosloski thanks for the hint - set, last is indeed isn't necessary, corrected.