Added a tutorial on reversepropellerbladex

Author: DavideOberto

bladex/__init__.py

@@ -17,6 +17,10 @@
 from .deform import Deformation
 from .params import ParamFile
 from .ndinterpolator import RBF, reconstruct_f, scipy_bspline
+from .reversepropeller import ReversePropellerInterface
+from .reversepropeller import BaseReversePropeller
 from .reversepropeller import ReversePropeller
+from .reversepropeller import ReversePropellerBladeX
+from .reversepropeller.reversepropeller import ReversePropeller
 from .shaft.cylinder_shaft import CylinderShaft
 from .intepolatedface import InterpolatedFace

tutorials/tutorial-7-blades-reversepropeller.ipynb

@@ -0,0 +1,145 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "c2ec2fd3",
+   "metadata": {},
+   "source": [
+    "# BladeX"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "ad979fa2",
+   "metadata": {},
+   "source": [
+    "## Tutorial 7: Retrieve blade's parameters from iges file using reversepropeller class"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "c0ab2647",
+   "metadata": {},
+   "source": [
+    "The goal of this tutorial is to show how to start from an iges file containing a blade to obtain its parameters, both for the sections (camber, thickness) and for the whole blade (chord lengths, pitch, rake, skew angles)."
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "5f21ab8e",
+   "metadata": {},
+   "source": [
+    "We start from some useful imports, as usual."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 5,
+   "id": "4362678f",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "from bladex import Blade, NacaProfile, ReversePropellerBladeX, ReversePropeller"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "45d4e8f4",
+   "metadata": {},
+   "source": [
+    "Then, we create the same blade as for Tutorial 6 and store it as iges file."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "aa8f973f",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Create sections with NACAProfile\n",
+    "n_sections = 10\n",
+    "sections = np.asarray([NacaProfile(digits='4412', n_points=50) for i in range(n_sections)])\n",
+    "\n",
+    "# Define blade parameters\n",
+    "radii = np.arange(1.0, 11.0, 1.0)\n",
+    "chord_lengths = np.array([0.05, 2.5, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7])\n",
+    "pitch = np.arange(1., 11.)\n",
+    "rake = np.arange(0.1, 1.1, 0.1)\n",
+    "skew_angles = np.arange(1., 21, 2.)\n",
+    "\n",
+    "# Create the Blade\n",
+    "blade = Blade(sections=sections,\n",
+    "        radii=radii,\n",
+    "        chord_lengths=chord_lengths,\n",
+    "        pitch=pitch,\n",
+    "        rake=rake,\n",
+    "        skew_angles=skew_angles)\n",
+    "\n",
+    "# We build the blade to be able to store its iges version\n",
+    "# The method build() internally calls apply_transformations()\n",
+    "blade.build()\n",
+    "\n",
+    "# Storing the blade in iges format\n",
+    "iges_filename='data/blade_reversepropeller.igs'\n",
+    "blade.export_iges(iges_filename)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "61aabd42",
+   "metadata": {},
+   "source": [
+    "Now that we have our starting iges file, we can retrieve the blades parameters. In BladeX you can either use the `ReversePropellerBladeX`, if the blade has been generated using BladeX, or `ReversePropeller` otherwise. The former is robust and more accurate for BladeX blades because it exploits known structure of the iges file."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "78b71d18",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# We need to specify: \n",
+    "#   i) the name of the iges file; \n",
+    "#   ii) values of the radii in mm at which the parameters are to be extracted;\n",
+    "#   iii) how many points to use to reconstruct the sections for each radius value\n",
+    "# In this case we use the same radii values (extrema excluded) of the blade's generation for a better comparison\n",
+    "radii_reverse = np.arange(2.0, 10., 1.)\n",
+    "# This operation may take few minutes\n",
+    "reverse = ReversePropellerBladeX(iges_filename, 1000*radii_reverse, num_points_top_bottom=100)\n",
+    "\n",
+    "# Now we can retrieve the blade's parameters\n",
+    "print(reverse.chord_length_list)\n",
+    "print(reverse.pitch_list)\n",
+    "print(reverse.rake_list)\n",
+    "print(reverse.skew_angles_list)\n",
+    "\n",
+    "# Finally, we can store the parameters to a csv file\n",
+    "reverse.save_global_parameters('data/blade_parameters.csv')"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "marinai",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.10.13"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}