Add modal example

examples/00-mapdl-examples/3d_notch.py

@@ -5,9 +5,10 @@
 ----------------------------------------------------
 
 This tutorial is the 3D corollary to the 2D plane example
-:ref:`ref_plane_stress_concentration`, but This example verifies the
+:ref:`ref_plane_stress_concentration`.
+However, this example verifies the
 stress concentration factor :math:`K-t` when modeling opposite single
-notches in a finite width thin plate
+notches in a finite width thin plate.
 
 First, start MAPDL as a service and disable all but error messages.
 """

examples/00-mapdl-examples/modal_beam.py

@@ -0,0 +1,143 @@
+"""
+.. _ref_modal_beam:
+
+=================================
+MAPDL modal beam analysis example
+=================================
+
+This example demonstrate how to performa a simple modal analysis
+and animate its results.
+
+
+Objective
+~~~~~~~~~
+This example models a simple 3D elastic beam made of BEAM188 elements.
+These beams elements are made of a linear elastic material similar to steel,
+and have a rectangular section.
+
+
+Procedure
+~~~~~~~~~
+
+* Launch MAPDL instance
+* Material properties
+* Geometry
+* Finite element model
+* Boundary conditions
+* Solving the model
+* Post-processing
+* Stop MAPDL
+
+"""
+
+###############################################################################
+# Launch MAPDL instance
+# =====================
+#
+# Launch MAPDL with interactive plotting
+from ansys.mapdl.core import launch_mapdl
+
+nmodes = 10
+# start MAPDL
+mapdl = launch_mapdl()
+print(mapdl)
+
+
+###############################################################################
+# Define material
+# ===============
+#
+# Define material
+mapdl.prep7()
+mapdl.mp("EX", 1, 2.1e11)
+mapdl.mp("PRXY", 1, 0.3)
+mapdl.mp("DENS", 1, 7800)
+
+###############################################################################
+# Create geometry
+# ===============
+#
+# Create keypoints and line
+mapdl.k(1)
+mapdl.k(2, 10)
+mapdl.l(1, 2)
+mapdl.lplot()
+
+###############################################################################
+# Define finite element model
+# ===========================
+#
+# Define element type/section type - Rectangular beam section
+mapdl.et(1, "BEAM188")
+mapdl.sectype(1, "BEAM", "RECT")
+mapdl.secoffset("CENT")
+mapdl.secdata(2, 1)
+
+# Mesh the line
+mapdl.type(1)
+mapdl.esize(1)
+mapdl.lesize("ALL")
+mapdl.lmesh("ALL")
+mapdl.eplot()
+mapdl.finish()
+
+###############################################################################
+# Specify boundary conditions
+# ===========================
+#
+# Fully fixed (clamped) end.
+mapdl.solution()  # Entering the solution processor.
+mapdl.nsel("S", "LOC", "X", "0")
+mapdl.d("ALL", "ALL")
+mapdl.allsel()
+mapdl.nplot(plot_bc=True, nnum=True)
+
+###############################################################################
+# Solve the model
+# ===============
+#
+# Setting modal analysis
+mapdl.antype("MODAL")
+mapdl.modopt("LANB", nmodes, 0, 200)
+mapdl.solve()
+mapdl.finish()
+
+###############################################################################
+# Postprocess
+# ===========
+#
+# Enter the post processor (post1)
+mapdl.post1()
+output = mapdl.set("LIST")
+print(output)
+
+result = mapdl.result
+
+###############################################################################
+# Animate results
+
+mode2plot = 2
+normalizeDisplacement = 1 / result.nodal_displacement(mode2plot - 1)[1].max()
+result.plot_nodal_displacement(
+    mode2plot,
+    show_displacement=True,
+    displacement_factor=normalizeDisplacement,
+    n_colors=10,
+)
+
+result.animate_nodal_displacement(
+    mode2plot,
+    loop=True,
+    add_text=False,
+    n_frames=100,
+    displacement_factor=normalizeDisplacement,
+    show_axes=False,
+    background="w",
+)
+
+###############################################################################
+# Stop MAPDL
+# ==========
+#
+mapdl.finish()
+mapdl.exit()