add basic static cyclic example (#1631)

Author: akaszynski

examples/00-mapdl-examples/cyclic_static.py

@@ -0,0 +1,81 @@
+"""
+.. _ref_cyclic_static_analysis:
+
+Static Cyclic Analysis
+----------------------
+
+Run a static cyclic analysis of an example rotor sector using the imperial unit
+system at 1000 RPM.
+
+
+"""
+from ansys.mapdl.reader import examples
+
+from ansys.mapdl.core import launch_mapdl
+
+# launch mapdl
+mapdl = launch_mapdl()
+
+
+###############################################################################
+# Load in the mesh
+# ~~~~~~~~~~~~~~~~
+# Load in the example sector and plot it.
+mapdl.cdread("db", examples.sector_archive_file)
+mapdl.eplot()
+
+
+###############################################################################
+# Make the rotor cyclic
+# ~~~~~~~~~~~~~~~~~~~~~
+# Enter the preprocessing routine and make the mesh cyclic.
+mapdl.prep7()
+mapdl.cyclic()
+
+
+###############################################################################
+# Set material properties
+# ~~~~~~~~~~~~~~~~~~~~~~~
+# Units are in imperial units and the material is (approximately) structural
+# steel.
+mapdl.mp("NUXY", 1, 0.31)
+mapdl.mp("DENS", 1, 4.1408e-04)
+mapdl.mp("EX", 1, 16900000)
+
+
+###############################################################################
+# Apply boundary conditions
+# ~~~~~~~~~~~~~~~~~~~~~~~~~
+# Apply a cyclic rotation at 1000 RPM and constrain the rotor at the center.
+mapdl.omega(0, 0, 1000)  # 1000 RPM
+
+mapdl.csys(1)  # enter the cyclic coordinate system
+mapdl.nsel("S", "loc", "x", 0.69, 0.71)  # radial between 0.69 - 0.71
+mapdl.d("ALL", "ALL")  # all DOF for those 8 nodes
+mapdl.allsel()
+mapdl.csys(0)  # return to cartesian coordinate system
+
+
+###############################################################################
+# Run a static analysis
+# ~~~~~~~~~~~~~~~~~~~~~
+# Run the MAPDL solver and print the output of the solution.
+mapdl.run("/SOLU")
+mapdl.antype("STATIC")
+output = mapdl.solve()
+mapdl.finish()
+print(output)
+
+
+###############################################################################
+# Plot the cyclic result
+# ~~~~~~~~~~~~~~~~~~~~~~
+# Plot the result using the legacy
+mapdl.result.plot_nodal_displacement(0)
+
+
+###############################################################################
+# Exit MAPDL
+# ~~~~~~~~~~
+# Finally, exit MAPDL.
+mapdl.exit()