h2lib-tests 13.2.801__py3-none-any.whl → 13.2.901__py3-none-any.whl

h2lib_tests/test_distributed_sections.py

@@ -1,12 +1,48 @@
+from pathlib import Path
+
 from h2lib_tests import tfp
+from numpy import testing as npt
 import pytest
+from scipy.spatial.transform import Rotation as R
+from wetb.hawc2.htc_file import HTCFile
+from wetb.hawc2.st_file import StFile
 
-from h2lib._h2lib import H2LibThread, H2Lib
-from h2lib.distributed_sections import DistributedSections, LinkType
-import numpy as np
-from numpy import testing as npt
+from h2lib import H2Lib
+from h2lib.distributed_sections import LinkType
+from h2lib_tests.dtu10mw import DTU10MWSimple
 import matplotlib.pyplot as plt
-from wetb.hawc2.htc_file import HTCFile
+import numpy as np
+
+
+class Plot():
+    def __init__(self, h2, ds_lst, ax=None):
+        self.h2 = h2
+        self.ds_lst = ds_lst
+        self.ax = ax or plt.figure().add_subplot(projection='3d')
+
+    def __call__(self, label, coo_index=[], h2=None, ds_lst=None, mainbody_coo=0):
+        self.h2 = h2 or self.h2
+        self.ds_lst = ds_lst or self.ds_lst
+        for ds in self.ds_lst:
+            sec_pos, sec_tsg = self.h2.get_distributed_section_position_orientation(ds, mainbody_coo_nr=mainbody_coo)
+            x, y, z = sec_pos.T
+            self.ax.plot(y, x, z, label=f'{label}, {ds.name}')
+
+            for (x, y, z), tsg in zip(sec_pos[coo_index], sec_tsg[coo_index]):
+                for (ex, ey, ez), c in zip(tsg.T * 10, 'rgb'):
+                    plt.plot([y, y + ey], [x, x + ex], [z, z + ez], c)
+
+    def show(self, show):
+        if show:
+            plt.axis('equal')
+            self.ax.set_zlim([10, -220])
+            plt.legend()
+            self.ax.set_xlabel('x')
+            self.ax.set_ylabel('y')
+            self.ax.set_zlabel('z')
+            plt.show()
+        else:
+            plt.close('all')
 
 
 def test_distributed_sections():
@@ -15,22 +51,11 @@ def test_distributed_sections():
         model_path = f"{tfp}DTU_10_MW/"
         htc_path = f"{tfp}DTU_10_MW/htc/DTU_10MW_RWT.htc"
         h2.init(htc_path=htc_path, model_path=model_path)
+        ds_lst = [h2.get_distributed_sections(LinkType.BLADE, link_id=i) for i in [1, 2, 3]]
+        plot = Plot(h2, ds_lst)
+        plot('', slice(None))
+        plot.show(0)
 
-        ax = plt.figure().add_subplot(projection='3d')
-        for i in [1, 2, 3]:
-            # ds = DistributedSections(name='', link_type=1, link_id=i, nsec=50)
-            ds = h2.get_distributed_sections(LinkType.BLADE, link_id=i)
-
-            sec_pos, sec_tsg = h2.get_distributed_section_position_orientation(ds, mainbody_coo_nr=0)
-            x, y, z = sec_pos.T
-            ax.plot(y, x, -z, label=ds.name)
-            for (x, y, z), tsg in zip(sec_pos, sec_tsg):
-                for (ex, ey, ez), c in zip(tsg.T * 10, 'rgb'):
-                    plt.plot([y, y + ey], [x, x + ex], [-z, -z - ez], c)
-        plt.legend()
-        plt.axis('equal')
-        if 0:
-            plt.show()
         ds = h2.get_distributed_sections(LinkType.BLADE, link_id=2)
         sec_pos, sec_tsg = h2.get_distributed_section_position_orientation(ds, mainbody_coo_nr=0)
         npt.assert_array_almost_equal(sec_pos[40], [70.58938502, -16.98280589, -77.95555399], 6)
@@ -46,21 +71,9 @@ def test_distributed_sections():
                    for mbdy_name in mbdy_name_lst}
         h2.initialize_distributed_sections()
         ds_dict['blade1_aero'] = h2.get_distributed_sections(LinkType.BLADE, 1)
-
-        ax = plt.figure().add_subplot(projection='3d')
-        for name, ds in ds_dict.items():
-            sec_pos, sec_tsg = h2.get_distributed_section_position_orientation(ds, mbdy_name_dict['hub2'])
-            x, y, z = sec_pos.T
-            ax.plot(y, x, -z, label=name)
-            for (x, y, z), tsg in zip(sec_pos, sec_tsg):
-                for (ex, ey, ez), c in zip(tsg.T * 10, 'rgb'):
-
-                    plt.plot([y, y + ey], [x, x + ex], [-z, -z - ez], c)
-
-        plt.axis('equal')
-        if 0:
-            plt.show()
-
+        plot = Plot(h2, ds_dict.values())
+        plot('test', slice(None), mainbody_coo=mbdy_name_dict['hub2'])
+        plot.show(0)
         # print(ds_lst[6].name, ds_lst[3].name)
         # mb_pos, b1_mb_tbg, b1_sec_pos, sec_tsb = h2.get_distributed_section_position_orientation(ds_lst[6])
         # hub_id = h2.get_mainbody_name_dict()['hub1']
@@ -92,16 +105,11 @@ def test_distributed_sections():
         frc[:, 1] = 100000
         h2.set_distributed_section_force_and_moment(ds_dict['tower'], sec_frc=frc, sec_mom=np.zeros((4, 3)))
 
+        plot = Plot(h2, ds_dict.values())
+        plot('t=0')
         h2.run(2)
-        for ds in ds_dict.values():
-            sec_pos, sec_tsb = h2.get_distributed_section_position_orientation(ds)
-            x, y, z = sec_pos.T
-            ax.plot(y, x, -z, '--', label=ds.name)
-        plt.legend()
-        plt.axis('scaled')
-        if 0:
-            plt.show()
-        plt.close('all')
+        plot('t=2')
+        plot.show(0)
         np.set_printoptions(linewidth=200)
         sec_pos, sec_tsb = h2.get_distributed_section_position_orientation(ds_dict['tower'])
         # print(np.round(mb_pos + [mb_tbg@sp for sp in sec_pos], 1).tolist())
@@ -112,6 +120,61 @@ def test_distributed_sections():
                                                 [0.0, 3.5, -115.6]], 1)
 
 
+def test_distributed_section_positions():
+
+    mp = Path(tfp) / 'DTU_10_MW/'
+    htc_path = mp / 'htc/DTU_10MW_RWT_only_tower.htc'
+    st = StFile(mp / 'data/DTU_10MW_RWT_Tower_st.dat')
+    st_new = st.main_data_sets[1][1].copy()
+    st_new[:, -2] = .4
+    st_new[:, -1] = .3
+    st.main_data_sets[1][1] = st_new
+    st.save(mp / 'data/DTU_10MW_RWT_Tower_st_tmp.dat')
+    htc = HTCFile(htc_path)
+    htc.new_htc_structure.main_body.timoschenko_input.filename = str(mp / 'data/DTU_10MW_RWT_Tower_st_tmp.dat')
+    htc.set_name('DTU_10MW_RWT_only_tower_tmp.htc')
+    htc.save()
+
+    with H2Lib() as h2:
+        h2.init(htc_path=htc.filename, model_path=htc.modelpath)
+        ds = h2.add_distributed_sections(mainbody_name='tower', section_relative_position=[0, .1, .5, 1])
+        h2.initialize_distributed_sections()
+
+        ax = plt.figure().add_subplot(projection='3d')
+        pos = h2.get_mainbody_nodes_state(1, 'pos')
+        npt.assert_array_equal(pos[:, :2].T, [[-0.4, -0.4, -0.4, -0.4, -0.4, -0.4, -0.4, -0.4, -0.4, -0.4, -0.4],
+                                              [0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3]])
+        ax.plot(*pos.T, '.-', label='initial nodes position')
+
+        sec_pos, sec_tsg = h2.get_distributed_section_position_orientation(ds)
+        npt.assert_array_almost_equal(sec_pos[:, :2], 0, 10)
+        ax.plot(*sec_pos.T, '.-', label='initial distributed sections, c2def')
+
+        st_new[:, -2] = 0.8
+        st_new[:, -1] = 0.9
+
+        h2.set_st(main_body_nr=1, st=st_new)
+
+        pos = h2.get_mainbody_nodes_state(1, 'pos')
+        npt.assert_array_equal(pos[:, :2].T, [[-0.8, -0.8, -0.8, -0.8, -0.8, -0.8, -0.8, -0.8, -0.8, -0.8, -0.8],
+                                              [0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9, 0.9]])
+        ax.plot(*pos.T, '.-', label='exy=1, nodes position')
+
+        sec_pos, sec_tsg = h2.get_distributed_section_position_orientation(ds)
+        ax.plot(*sec_pos.T, '.-', label='exy=1 distributed sections, c2def')
+
+        plt.legend()
+        ax.set_xlabel('x')
+        ax.set_ylabel('y')
+        ax.set_zlabel('z')
+        ax.set_zlim([0, -120])
+        if 0:
+            plt.show()
+        plt.close('all')
+
+        npt.assert_array_equal(sec_pos[:, :2], 0, 10)
+
+
 def test_distributed_sections_static_solver():
     with H2Lib() as h2:
         model_path = f"{tfp}DTU_10_MW/"
@@ -137,16 +200,17 @@ def test_distributed_sections_static_solver():
         frc[:, 1] = 100000
         h2.set_distributed_section_force_and_moment(ds, sec_frc=frc, sec_mom=frc * 0)
         h2.solver_static_run(reset_structure=True)
-        draw('set frc + static solver', [0.0, 1.8295, -115.6123])
+        draw('set frc + static solver', [0.0, 1.8293, -115.6123])
         c2_def = np.concatenate([pos, pos[:, :1]], 1)
         c2_def[:, 0] = np.r_[np.arange(0, 60, 10), np.arange(50, 0, -10)]
-        h2.set_c2_def('tower', c2_def)
+        tower_id = h2.get_mainbody_name_dict()["tower"]
+        h2.set_c2_def(tower_id, c2_def)
         draw('set_c2_def', [10.0, 0.1184, -115.6296])
         h2.solver_static_run(reset_structure=True)
-        draw('static solver', [10.2777, 2.8085, -115.5066])
+        draw('static solver', [10.2695, 2.8081, -115.5122])
         h2.set_distributed_section_force_and_moment(ds, sec_frc=-frc, sec_mom=frc * 0)
         h2.solver_static_run(reset_structure=True)
-        draw('set -frc + static solver', [10.2859, -2.8064, -115.5011])
+        draw('set -frc + static solver', [10.2695, -2.8081, -115.5122])
         if 0:
             plt.legend()
             plt.show()
@@ -160,9 +224,11 @@ def test_set_distributed_section_force_and_moment_coo():
         htc = HTCFile(model_path + "htc/DTU_10MW_RWT_only_tower.htc")
         htc.new_htc_structure.orientation.base.body_eulerang = 0, 0, 90
         htc.save(model_path + "htc/DTU_10MW_RWT_only_tower_rot90.htc")
-        print(htc)
         h2.init(htc_path=htc.filename, model_path=model_path)
         ds = h2.add_distributed_sections(mainbody_name='tower', section_relative_position=[0, .5, 1])
+        with pytest.raises(AssertionError, match='Call initialize_distributed_sections before get_distributed_section_force_and_moment'):
+            h2.get_distributed_section_force_and_moment(ds)
+
         h2.initialize_distributed_sections()
 
         ax = plt.figure().add_subplot(projection='3d')
@@ -170,7 +236,7 @@ def test_set_distributed_section_force_and_moment_coo():
         def draw(label, ref):
             h2.solver_static_run()
             pos = h2.get_mainbody_nodes_state(mainbody_nr=1, state='pos')
-            print(label, np.round(pos[-1], 4).tolist())
+            # print(label, np.round(pos[-1], 4).tolist())
             ax.plot(*pos.T, marker='.', label=label)
             npt.assert_array_almost_equal(pos[-1], ref, 4)
 
@@ -178,10 +244,10 @@ def test_set_distributed_section_force_and_moment_coo():
         frc = np.zeros((3, 3))
         frc[:, 1] = 2e6
         h2.set_distributed_section_force_and_moment(ds, sec_frc=frc, sec_mom=frc * 0)
-        draw('frc_y_global', [-0.0, 34.7544, -109.7119])
+        draw('frc_y_global', [0.0, 33.8576, -110.0939])
 
         h2.set_distributed_section_force_and_moment(ds, sec_frc=frc, sec_mom=frc * 0, mainbody_coo_nr=1)
-        draw('frc_y_tower', [-34.7545, 0.0, -109.7119])
+        draw('frc_y_tower', [-33.8576, -0.0, -110.0939])
 
         if 0:
             ax.axis('equal')
@@ -195,3 +261,145 @@ def test_set_distributed_section_force_and_moment_coo():
             plt.show()
         else:
             plt.close('all')
+
+
+def test_set_frc_blade():
+    with H2Lib(suppress_output=0) as h2:
+        htc = DTU10MWSimple(rotor_speed=.6, pitch=0, nbodies=1)
+        htc.set_name('DTU_10MW_RWT_fixed_rotspeed')
+        htc.aero.aerocalc_method = 0
+        htc.save()
+        h2.init(htc_path=htc.filename, model_path=htc.modelpath)
+
+        ds_lst = [h2.get_distributed_sections(LinkType.BLADE, link_id=i) for i in [1, 2, 3]]
+        h2.solver_static_run()
+        plot = Plot(h2, ds_lst)
+        plot('initial', [-1])
+        ini_tip_pos = [h2.get_distributed_section_position_orientation(ds, mainbody_coo_nr=0)[0][-1] for ds in ds_lst]
+        ini_tip3_tsg = h2.get_distributed_section_position_orientation(ds_lst[2], mainbody_coo_nr=9)[1][-1]
+        zeros = np.zeros((50, 3))
+        h2.set_distributed_section_force_and_moment(ds=ds_lst[0], sec_frc=zeros + [10000, 0, 0], sec_mom=zeros)
+        h2.set_distributed_section_force_and_moment(ds=ds_lst[1], sec_frc=zeros + [0, 10000, 0], sec_mom=zeros)
+        h2.set_distributed_section_force_and_moment(ds=ds_lst[2], sec_frc=zeros, sec_mom=zeros + [0, 0, 100000])
+        h2.solver_static_run()
+
+        plot('deflected', [-1])
+        def_tip_pos = [h2.get_distributed_section_position_orientation(ds, mainbody_coo_nr=0)[0][-1] for ds in ds_lst]
+        def_tip3_tsg = h2.get_distributed_section_position_orientation(ds_lst[2], mainbody_coo_nr=9)[1][-1]
+        npt.assert_array_almost_equal(np.array(def_tip_pos) - ini_tip_pos,
+                                      [[7.545511, -0.270957, 0.189638],  # mainly x
+                                       [0.942178, 11.676871, 2.980586],  # mainly y
+                                       [0.596675, -6.677523, -2.850951]],  # complex due to bend-twist coupling
+                                      6)
+
+        npt.assert_array_almost_equal(R.from_matrix(ini_tip3_tsg @ def_tip3_tsg).as_euler('xyz', 1),
+                                      [30.762956, 0.310381, 44.774262], 6)  # rotation around z (but also x due to bend-twist coupling
+        plot.show(0)
+
+
+@pytest.mark.parametrize('coo_nr', [-1, 0, 1])
+def test_get_set_frc_blade(coo_nr):
+    '''
+    1 Run static solver with aerodynamic and extract position, frc and mom of distributed sections
+    2 New simulation, set frc and mom and run static solver without aerocalc
+    3 compare position of distributed sections
+    '''
+    with H2Lib(suppress_output=0) as h2:
+        htc = DTU10MWSimple(rotor_speed=.6, pitch=0, nbodies=1)
+        htc.set_name('DTU_10MW_RWT_fixed_rotspeed')
+        htc.simulation.convergence_limits.delete()
+        htc.wind.wsp = 10
+        htc.aero.tiploss_method = 0  # tiploss makes the static solver unstable
+        htc.save()
+        h2.init(htc_path=htc.filename, model_path=htc.modelpath)
+
+        aero_ds_lst = [h2.get_distributed_sections(LinkType.BLADE, link_id=i) for i in [1, 2, 3]]
+
+        mb_dict = h2.get_mainbody_name_dict()
+        p = h2.get_distributed_section_position_orientation(aero_ds_lst[0], mb_dict['blade1'])[0]
+        r = np.r_[0, np.cumsum(np.sqrt((np.diff(p, 1, 0)**2).sum(1)))]  # curve-length radius
+        rR = r / r[-1]
+        for ds in aero_ds_lst:
+            h2.add_distributed_sections(ds.name, rR)
+        h2.initialize_distributed_sections()
+
+        body_ds_lst = [h2.get_distributed_sections(link_type=LinkType.BODY, link_id=i) for i in [1, 2, 3]]
+
+        plot = Plot(h2, aero_ds_lst)
+        h2.solver_static_run()
+        plot('aerosections after static solver')
+        # plot('bodysections after static solver', ds_lst=body_ds_lst)
+        pos_coo = np.maximum(coo_nr, 0)
+        ref_pos = [h2.get_distributed_section_position_orientation(ds, pos_coo)[0] for ds in aero_ds_lst]
+        aero2body = [np.array(h2.get_distributed_section_position_orientation(aero_ds, pos_coo)[0]) -
+                     h2.get_distributed_section_position_orientation(body_ds, pos_coo)[0]
+                     for aero_ds, body_ds in zip(aero_ds_lst, body_ds_lst)]
+        frc, mom = zip(*[h2.get_distributed_section_force_and_moment(ds, mainbody_coo_nr=coo_nr)
+                       for ds in aero_ds_lst])
+        if coo_nr == 0:
+            aero_frc = h2.get_aerosections_forces()
+            aero_mom = h2.get_aerosections_moments()
+            aero_pos = h2.get_aerosections_position()
+            npt.assert_array_equal(aero_frc, frc)
+            npt.assert_array_almost_equal(aero_mom, mom)
+
+            npt.assert_array_equal(ref_pos, aero_pos)
+
+    with H2Lib(suppress_output=0) as h2:
+        htc = DTU10MWSimple(rotor_speed=.6, pitch=0, nbodies=1)
+        htc.simulation.convergence_limits.delete()
+        htc.set_name('DTU_10MW_RWT_fixed_rotspeed')
+        htc.aero.aerocalc_method = 0
+        htc.save()
+        h2.init(htc_path=htc.filename, model_path=htc.modelpath)
+
+        aero_ds_lst = [h2.get_distributed_sections(LinkType.BLADE, link_id=i) for i in [1, 2, 3]]
+
+        # set_distributed_section_force_and_moment maps the moment from section moments to moments around
+        # elastic axis. This mapping depends on the deflection, so we need to run a couple of times to
+        # get the right deflection and thereby the right mapping
+        for _ in range(5):
+            for sec_frc, sec_mom, ds in zip(frc, mom, aero_ds_lst):
+                h2.set_distributed_section_force_and_moment(ds, sec_frc, sec_mom, mainbody_coo_nr=coo_nr)
+            h2.solver_static_run()
+
+        pos = [h2.get_distributed_section_position_orientation(ds, np.maximum(coo_nr, 0))[0] for ds in aero_ds_lst]
+        plot = Plot(h2, aero_ds_lst, ax=plot.ax)
+        plot('set frc')
+        plot.show(0)
+        npt.assert_array_almost_equal(ref_pos, pos)
+
+    if coo_nr >= 0:
+        # for global and blade1 coordinates:
+        # apply forces and moments (mapped to c2df) to distributed sections (at c2def)
+        # and compare position of aerodynamic sections
+        with H2Lib(suppress_output=0) as h2:
+            htc = DTU10MWSimple(rotor_speed=.6, pitch=0, nbodies=1)
+            htc.simulation.convergence_limits.delete()
+            htc.set_name('DTU_10MW_RWT_fixed_rotspeed')
+            htc.aero.aerocalc_method = 0
+            htc.save()
+            h2.init(htc_path=htc.filename, model_path=htc.modelpath)
+
+            aero_ds_lst = [h2.get_distributed_sections(LinkType.BLADE, link_id=i) for i in [1, 2, 3]]
+
+            for ds in aero_ds_lst:
+                h2.add_distributed_sections(ds.name, rR)
+            h2.initialize_distributed_sections()
+
+            body_ds_lst = [h2.get_distributed_sections(link_type=LinkType.BODY, link_id=i) for i in [1, 2, 3]]
+
+            # set_distributed_section_force_and_moment maps the moment from section moments to moments around
+            # elastic axis. This mapping depends on the deflection, so we need to run a couple of times to
+            # get the right deflection and thereby the right mapping
+            for _ in range(5):
+                for sec_frc, sec_mom, ds, a2b in zip(frc, mom, body_ds_lst, aero2body):
+                    body_mom = sec_mom + np.cross(sec_frc, -a2b)
+                    h2.set_distributed_section_force_and_moment(ds, sec_frc, body_mom, mainbody_coo_nr=coo_nr)
+                h2.solver_static_run()
+            pos = [h2.get_distributed_section_position_orientation(ds, np.maximum(coo_nr, 0))[0] for ds in aero_ds_lst]
+            plot = Plot(h2, aero_ds_lst, ax=plot.ax)
+            plot('set frc')
+            plot.show(0)
+            # not sure why the accuracy is much worse. Maybe the distributed
+            npt.assert_array_almost_equal(ref_pos, pos, 2)

h2lib_tests/test_topology_h2lib.py

@@ -239,27 +239,35 @@ def test_set_orientation_base_speed(h2_dtu_10mw_only_blade):
 
 
 def test_set_c2_def_too_few_sections(h2_dtu_10mw_only_blade):
+    blade_id = h2_dtu_10mw_only_blade.get_mainbody_name_dict()["blade1"]
     with pytest.raises(ValueError, match="TOO_FEW_SECTIONS_IN_C2DEF"):
-        h2_dtu_10mw_only_blade.set_c2_def("blade1", np.zeros((1, 4)))
+        h2_dtu_10mw_only_blade.set_c2_def(blade_id, np.zeros((1, 4)))
 
 
 def test_set_c2_def_wrong_number_of_columns(h2_dtu_10mw_only_blade):
+    blade_id = h2_dtu_10mw_only_blade.get_mainbody_name_dict()["blade1"]
     with pytest.raises(ValueError, match="WRONG_NUMBER_OF_COLUMNS"):
-        h2_dtu_10mw_only_blade.set_c2_def("blade1", np.zeros((2, 6)))
+        h2_dtu_10mw_only_blade.set_c2_def(blade_id, np.zeros((2, 6)))
 
 
 def test_set_c2_def_main_body_not_found(h2_dtu_10mw_only_blade):
     with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
         h2_dtu_10mw_only_blade.set_c2_def(
-            "tower",
+            -10,
+            np.array([[0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 1.0, 0.0]]),
+        )
+    with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
+        h2_dtu_10mw_only_blade.set_c2_def(
+            123,
             np.array([[0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 1.0, 0.0]]),
         )
 
 
 def test_set_c2_def_different_nsec(h2_dtu_10mw_only_blade):
+    blade_id = h2_dtu_10mw_only_blade.get_mainbody_name_dict()["blade1"]
     with pytest.raises(ValueError, match="DIFFERENT_NSEC"):
         h2_dtu_10mw_only_blade.set_c2_def(
-            "blade1",
+            blade_id,
             np.array([[0.0, 0.0, 0.0, 0.0], [0.0, 0.0, 1.0, 0.0]]),
         )
 
@@ -267,19 +275,21 @@ def test_set_c2_def_different_nsec(h2_dtu_10mw_only_blade):
 def test_set_c2_def_uniform_node_distribution(
     h2_dtu_10mw_only_blade_uniform_node_distribution,
 ):
+    blade_id = h2_dtu_10mw_only_blade_uniform_node_distribution.get_mainbody_name_dict()["blade1"]
     with pytest.raises(
         NotImplementedError,
     ):
         h2_dtu_10mw_only_blade_uniform_node_distribution.set_c2_def(
-            "blade1", np.ones((27, 4), dtype=np.float64, order="F")
+            blade_id, np.ones((27, 4), dtype=np.float64, order="F")
         )
 
 
 def test_set_c2_def_beam_too_short(h2_dtu_10mw_only_blade):
     c2def = np.zeros((27, 4), dtype=np.float64, order="F")
     c2def[:, 2] = np.linspace(0.0, 1e-9, c2def.shape[0])
+    blade_id = h2_dtu_10mw_only_blade.get_mainbody_name_dict()["blade1"]
     with pytest.raises(ValueError, match="BEAM_TOO_SHORT"):
-        h2_dtu_10mw_only_blade.set_c2_def("blade1", c2def)
+        h2_dtu_10mw_only_blade.set_c2_def(blade_id, c2def)
 
 
 def test_set_c2_def_blade_static_back_to_original(
@@ -295,6 +305,7 @@ def test_set_c2_def_blade_static_back_to_original(
     # Get the clamped DTU 10 MW blade subjected to high gravity loading.
     _, blade_c2def, _ = write_dtu10mw_only_blade_high_gravity
     h2 = h2_dtu10mw_only_blade_high_gravity
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Sensors 1 to 6 are the displacement and rotation of the blade tip.
     #   h2.get_sensor_info(1)
@@ -307,7 +318,7 @@ def test_set_c2_def_blade_static_back_to_original(
     c2def_new = blade_c2def.copy()
     for factor in (0.6, 0.7, 0.8, 0.9, 1.1, 1.2, 1.3, 1.4):
         c2def_new[:, 2] = factor * blade_c2def[:, 2]
-        h2.set_c2_def("blade1", c2def_new)
+        h2.set_c2_def(blade_id, c2def_new)
         # Since we are smoothly changing c2_def, it makes sense to start the static solver from the last converged solution.
         h2.solver_static_run(reset_structure=False)
 
@@ -319,7 +330,7 @@ def test_set_c2_def_blade_static_back_to_original(
             npt.assert_allclose(blade_tip_actual, blade_tip_desired)
 
     # Restore blade c2_def.
-    h2.set_c2_def("blade1", blade_c2def)
+    h2.set_c2_def(blade_id, blade_c2def)
     h2.solver_static_run(reset_structure=True)
     blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
     npt.assert_allclose(blade_tip_actual, blade_tip_desired)
@@ -338,7 +349,8 @@ def test_set_c2_def_blade_static_deformed(
     # Set c2_def in the original blade, thus making it equivalent to the deformed one.
     # Then, solve again the static problem.
     _, c2def_deformed = write_dtu10mw_only_blade_high_gravity_deformed
-    h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_deformed)
+    blade_id = h2_dtu10mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_dtu10mw_only_blade_high_gravity.set_c2_def(blade_id, c2def_deformed)
     h2_dtu10mw_only_blade_high_gravity.solver_static_run(reset_structure=True)
     blade_tip_actual = h2_dtu10mw_only_blade_high_gravity_deformed.get_sensor_values((1, 2, 3, 4, 5, 6))
 
@@ -346,7 +358,7 @@ def test_set_c2_def_blade_static_deformed(
 
     # Restore c2_def.
     _, c2def_original, _ = write_dtu10mw_only_blade_high_gravity
-    h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_original)
+    h2_dtu10mw_only_blade_high_gravity.set_c2_def(blade_id, c2def_original)
 
 
 @pytest.mark.skip(reason="The system eigenanalysis cannot be called more than once.")
@@ -358,6 +370,7 @@ def test_set_c2_def_blade_eig_back_to_original(
     # We use the fixture with high gravity because it also returns c2_def.
     _, blade_c2def, _ = write_dtu10mw_only_blade_high_gravity
     h2 = h2_dtu10mw_only_blade_high_gravity
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Solve the eigenvalue problem.
     h2.structure_reset()
@@ -368,7 +381,7 @@ def test_set_c2_def_blade_eig_back_to_original(
     c2def_new = blade_c2def.copy()
     for factor in (0.6, 0.7, 0.8, 0.9, 1.1, 1.2, 1.3, 1.4):
         c2def_new[:, 2] = factor * blade_c2def[:, 2]
-        h2.set_c2_def("blade1", c2def_new)
+        h2.set_c2_def(blade_id, c2def_new)
         h2.structure_reset()
         h2.linearize()
         natural_frequencies_actual, damping_ratios_actual = h2.do_system_eigenanalysis(n_modes=6, include_damping=True)
@@ -380,7 +393,7 @@ def test_set_c2_def_blade_eig_back_to_original(
             npt.assert_allclose(damping_ratios_actual, damping_ratios_desired)
 
     # Restore blade c2_def.
-    h2.set_c2_def("blade1", blade_c2def)
+    h2.set_c2_def(blade_id, blade_c2def)
     h2.structure_reset()
     h2.linearize()
     natural_frequencies_actual, damping_ratios_actual = h2.do_system_eigenanalysis(n_modes=6, include_damping=True)
@@ -396,6 +409,7 @@ def test_set_c2_def_blade_inertia_back_to_original(
     # We use the fixture with high gravity because it also returns c2_def.
     _, blade_c2def, _ = write_dtu10mw_only_blade_high_gravity_1_body
     h2 = h2_dtu10mw_only_blade_high_gravity_1_body
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Get the inertia properties.
     h2.structure_reset()
@@ -405,7 +419,7 @@ def test_set_c2_def_blade_inertia_back_to_original(
     c2def_new = blade_c2def.copy()
     for factor in (0.6, 0.7, 0.8, 0.9, 1.1, 1.2, 1.3, 1.4):
         c2def_new[:, 2] = factor * blade_c2def[:, 2]
-        h2.set_c2_def("blade1", c2def_new)
+        h2.set_c2_def(blade_id, c2def_new)
         inertia_actual = h2.body_output_mass(0)
 
         # Must differ from the desired ones.
@@ -414,7 +428,7 @@ def test_set_c2_def_blade_inertia_back_to_original(
                 npt.assert_allclose(inertia_actual[i], inertia_desired[i])
 
     # Restore blade c2_def.
-    h2.set_c2_def("blade1", blade_c2def)
+    h2.set_c2_def(blade_id, blade_c2def)
     inertia_actual = h2.body_output_mass(0)
     for i in range(4):  # Loop over tuple of arrays.
         npt.assert_allclose(inertia_actual[i], inertia_desired[i])
@@ -432,7 +446,8 @@ def test_set_c2_def_blade_inertia_deformed(
 
     # Set c2_def in the original blade, thus making it equivalent to the deformed one.
     _, c2def_deformed = write_dtu10mw_only_blade_high_gravity_deformed
-    h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_deformed)
+    blade_id = h2_dtu10mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_dtu10mw_only_blade_high_gravity.set_c2_def(blade_id, c2def_deformed)
 
     # Compare inertia properties for all bodies.
     nbdy, _ = h2_dtu10mw_only_blade_high_gravity_deformed.get_number_of_bodies_and_constraints()
@@ -445,7 +460,7 @@ def test_set_c2_def_blade_inertia_deformed(
 
     # Restore c2_def.
     _, c2def_original, _ = write_dtu10mw_only_blade_high_gravity
-    h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_original)
+    h2_dtu10mw_only_blade_high_gravity.set_c2_def(blade_id, c2def_original)
 
 
 def test_set_c2_def_blade_element_deformed(
@@ -460,7 +475,8 @@ def test_set_c2_def_blade_element_deformed(
 
     # Set c2_def in the original blade, thus making it equivalent to the deformed one.
     _, c2def_deformed = write_dtu10mw_only_blade_high_gravity_deformed
-    h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_deformed)
+    blade_id = h2_dtu10mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_dtu10mw_only_blade_high_gravity.set_c2_def(blade_id, c2def_deformed)
 
     # Compare element matrices for all bodies.
     nelem = h2_dtu10mw_only_blade_high_gravity_deformed.get_number_of_elements()
@@ -475,31 +491,36 @@ def test_set_c2_def_blade_element_deformed(
 
     # Restore c2_def.
     _, c2def_original, _ = write_dtu10mw_only_blade_high_gravity
-    h2_dtu10mw_only_blade_high_gravity.set_c2_def("blade1", c2def_original)
+    h2_dtu10mw_only_blade_high_gravity.set_c2_def(blade_id, c2def_original)
 
 
 def test_set_st_wrong_number_of_columns(h2_dtu_10mw_only_blade):
+    blade_id = h2_dtu_10mw_only_blade.get_mainbody_name_dict()["blade1"]
     with pytest.raises(ValueError, match="WRONG_NUMBER_OF_COLUMNS"):
-        h2_dtu_10mw_only_blade.set_st("blade1", np.empty((2, 6)))
+        h2_dtu_10mw_only_blade.set_st(blade_id, np.empty((2, 6)))
 
 
 def test_set_st_main_body_not_found(h2_dtu_10mw_only_blade):
     with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
-        h2_dtu_10mw_only_blade.set_st("tower", np.empty((2, 19)))
+        h2_dtu_10mw_only_blade.set_st(-3, np.empty((2, 19)))
+    with pytest.raises(ValueError, match="MAIN_BODY_NOT_FOUND"):
+        h2_dtu_10mw_only_blade.set_st(123, np.empty((2, 19)))
 
 
 def test_set_st_z_not_continuously_increasing(h2_dtu_10mw_only_blade):
+    blade_id = h2_dtu_10mw_only_blade.get_mainbody_name_dict()["blade1"]
     with pytest.raises(ValueError, match="ST_Z_NOT_CONTINUOUSLY_INCREASING"):
-        h2_dtu_10mw_only_blade.set_st("blade1", np.zeros((3, 19)))
+        h2_dtu_10mw_only_blade.set_st(blade_id, np.zeros((3, 19)))
 
 
 def test_set_st_uniform_node_distribution(
     h2_dtu_10mw_only_blade_uniform_node_distribution,
 ):
+    blade_id = h2_dtu_10mw_only_blade_uniform_node_distribution.get_mainbody_name_dict()["blade1"]
     with pytest.raises(NotImplementedError):
         st = np.zeros((2, 19))
         st[1, 0] = 1.0
-        h2_dtu_10mw_only_blade_uniform_node_distribution.set_st("blade1", st)
+        h2_dtu_10mw_only_blade_uniform_node_distribution.set_st(blade_id, st)
 
 
 def test_set_st_classic_timoshenko_blade_static_back_to_original(
@@ -515,6 +536,7 @@ def test_set_st_classic_timoshenko_blade_static_back_to_original(
     # Get the clamped DTU 10 MW blade subjected to high gravity loading.
     _, _, blade_st = write_dtu10mw_only_blade_high_gravity
     h2 = h2_dtu10mw_only_blade_high_gravity
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Sensors 1 to 6 are the displacement and rotation of the blade tip.
     #   h2.get_sensor_info(1)
@@ -530,7 +552,7 @@ def test_set_st_classic_timoshenko_blade_static_back_to_original(
         st_new = factor * blade_st.main_data_sets[1][1]
         # We do not change z, to prevent ST_Z_NOT_CONTINUOUSLY_INCREASING.
         st_new[:, 0] = blade_st.main_data_sets[1][1][:, 0]
-        h2.set_st("blade1", st_new)
+        h2.set_st(blade_id, st_new)
         h2.solver_static_run(reset_structure=True)
 
         # Get new blade displacement.
@@ -541,7 +563,7 @@ def test_set_st_classic_timoshenko_blade_static_back_to_original(
             npt.assert_allclose(blade_tip_actual, blade_tip_desired)
 
     # Restore blade c2_def.
-    h2.set_st("blade1", blade_st.main_data_sets[1][1])
+    h2.set_st(blade_id, blade_st.main_data_sets[1][1])
     h2.solver_static_run(reset_structure=True)
     blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
     npt.assert_allclose(blade_tip_actual, blade_tip_desired)
@@ -553,6 +575,7 @@ def test_set_st_classic_timoshenko_blade_static_changed_st(
     h2_dtu10mw_only_blade_high_gravity,
     h2_dtu10mw_only_blade_high_gravity_changed_st,
 ):
+    blade_id = h2_dtu10mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
     # Solve the static problem with the changed ST loaded directly by HAWC2.
     h2_dtu10mw_only_blade_high_gravity_changed_st.solver_static_run(reset_structure=True)
     blade_tip_desired = h2_dtu10mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
@@ -560,7 +583,7 @@ def test_set_st_classic_timoshenko_blade_static_changed_st(
     # Set ST in the original blade, thus making it equivalent to the changed one.
     # Then, solve again the static problem.
     _, st_changed = write_dtu10mw_only_blade_high_gravity_changed_st
-    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    h2_dtu10mw_only_blade_high_gravity.set_st(blade_id, st_changed.main_data_sets[1][1])
     h2_dtu10mw_only_blade_high_gravity.solver_static_run(reset_structure=True)
     blade_tip_actual = h2_dtu10mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
 
@@ -568,7 +591,7 @@ def test_set_st_classic_timoshenko_blade_static_changed_st(
 
     # Restore c2_def.
     _, _, st_riginal = write_dtu10mw_only_blade_high_gravity
-    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_riginal.main_data_sets[1][1])
+    h2_dtu10mw_only_blade_high_gravity.set_st(blade_id, st_riginal.main_data_sets[1][1])
 
 
 def test_set_st_classic_timoshenko_inertia_back_to_original(
@@ -579,6 +602,7 @@ def test_set_st_classic_timoshenko_inertia_back_to_original(
     # We use the fixture with high gravity because it also returns st.
     _, _, st = write_dtu10mw_only_blade_high_gravity_1_body
     h2 = h2_dtu10mw_only_blade_high_gravity_1_body
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Get the inertia properties.
     h2.structure_reset()
@@ -593,7 +617,7 @@ def test_set_st_classic_timoshenko_inertia_back_to_original(
         # Make the blade tip heavier to change center of gravity.
         density_new *= np.linspace(0.8, 2.5, density_new.size)
         st.set_value(mset=1, set=1, m=density_new)
-        h2.set_st("blade1", st.main_data_sets[1][1])
+        h2.set_st(blade_id, st.main_data_sets[1][1])
         inertia_actual = h2.body_output_mass(0)
 
         # Must differ from the desired ones.
@@ -603,7 +627,7 @@ def test_set_st_classic_timoshenko_inertia_back_to_original(
 
     # Restore ST.
     st.set_value(mset=1, set=1, m=mass_original)
-    h2.set_st("blade1", st.main_data_sets[1][1])
+    h2.set_st(blade_id, st.main_data_sets[1][1])
     inertia_actual = h2.body_output_mass(0)
     for i in range(4):  # Loop over tuple of arrays.
         npt.assert_allclose(inertia_actual[i], inertia_desired[i])
@@ -621,7 +645,8 @@ def test_set_st_classic_timoshenko_blade_inertia_changed_st(
 
     # Set ST in the original blade, thus making it equivalent to the changed one.
     _, st_changed = write_dtu10mw_only_blade_high_gravity_changed_st
-    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    blade_id = h2_dtu10mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_dtu10mw_only_blade_high_gravity.set_st(blade_id, st_changed.main_data_sets[1][1])
 
     # Compare inertia properties for all bodies.
     nbdy, _ = h2_dtu10mw_only_blade_high_gravity_changed_st.get_number_of_bodies_and_constraints()
@@ -634,7 +659,7 @@ def test_set_st_classic_timoshenko_blade_inertia_changed_st(
 
     # Restore ST.
     _, _, st_original = write_dtu10mw_only_blade_high_gravity
-    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+    h2_dtu10mw_only_blade_high_gravity.set_st(blade_id, st_original.main_data_sets[1][1])
 
 
 def test_set_st_classic_timoshenko_blade_element_changed_st(
@@ -649,7 +674,8 @@ def test_set_st_classic_timoshenko_blade_element_changed_st(
 
     # Set ST in the original blade, thus making it equivalent to the changed one.
     _, st_changed = write_dtu10mw_only_blade_high_gravity_changed_st
-    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    blade_id = h2_dtu10mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_dtu10mw_only_blade_high_gravity.set_st(blade_id, st_changed.main_data_sets[1][1])
 
     # Compare element matrices for all bodies.
     nelem = h2_dtu10mw_only_blade_high_gravity_changed_st.get_number_of_elements()
@@ -664,7 +690,7 @@ def test_set_st_classic_timoshenko_blade_element_changed_st(
 
     # Restore ST.
     _, _, st_original = write_dtu10mw_only_blade_high_gravity
-    h2_dtu10mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+    h2_dtu10mw_only_blade_high_gravity.set_st(blade_id, st_original.main_data_sets[1][1])
 
 
 def test_set_st_fpm_blade_static_back_to_original(
@@ -680,6 +706,7 @@ def test_set_st_fpm_blade_static_back_to_original(
     # Get the clamped IEA 22 MW blade subjected to high gravity loading.
     _, _, blade_st = write_iea22mw_only_blade_high_gravity
     h2 = h2_iea22mw_only_blade_high_gravity
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Sensors 1 to 6 are the displacement and rotation of the blade tip.
     #   h2.get_sensor_info(1)
@@ -695,7 +722,7 @@ def test_set_st_fpm_blade_static_back_to_original(
         st_new = factor * blade_st.main_data_sets[1][1]
         # We do not change z, to prevent ST_Z_NOT_CONTINUOUSLY_INCREASING.
         st_new[:, 0] = blade_st.main_data_sets[1][1][:, 0]
-        h2.set_st("blade1", st_new)
+        h2.set_st(blade_id, st_new)
         h2.solver_static_run(reset_structure=True)
 
         # Get new blade displacement.
@@ -706,7 +733,7 @@ def test_set_st_fpm_blade_static_back_to_original(
             npt.assert_allclose(blade_tip_actual, blade_tip_desired)
 
     # Restore blade c2_def.
-    h2.set_st("blade1", blade_st.main_data_sets[1][1])
+    h2.set_st(blade_id, blade_st.main_data_sets[1][1])
     h2.solver_static_run(reset_structure=True)
     blade_tip_actual = h2.get_sensor_values((1, 2, 3, 4, 5, 6))
     npt.assert_allclose(blade_tip_actual, blade_tip_desired)
@@ -725,7 +752,8 @@ def test_set_st_fpm_blade_static_changed_st(
     # Set ST in the original blade, thus making it equivalent to the changed one.
     # Then, solve again the static problem.
     _, st_changed = write_iea22mw_only_blade_high_gravity_changed_st
-    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    blade_id = h2_iea22mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_iea22mw_only_blade_high_gravity.set_st(blade_id, st_changed.main_data_sets[1][1])
     h2_iea22mw_only_blade_high_gravity.solver_static_run(reset_structure=True)
     blade_tip_actual = h2_iea22mw_only_blade_high_gravity_changed_st.get_sensor_values((1, 2, 3, 4, 5, 6))
 
@@ -733,7 +761,7 @@ def test_set_st_fpm_blade_static_changed_st(
 
     # Restore c2_def.
     _, _, st_riginal = write_iea22mw_only_blade_high_gravity
-    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_riginal.main_data_sets[1][1])
+    h2_iea22mw_only_blade_high_gravity.set_st(blade_id, st_riginal.main_data_sets[1][1])
 
 
 def test_set_st_fpm_inertia_back_to_original(
@@ -744,6 +772,7 @@ def test_set_st_fpm_inertia_back_to_original(
     # We use the fixture with high gravity because it also returns st.
     _, _, st = write_iea22mw_only_blade_high_gravity_1_body
     h2 = h2_iea22mw_only_blade_high_gravity_1_body
+    blade_id = h2.get_mainbody_name_dict()["blade1"]
 
     # Get the inertia properties.
     h2.structure_reset()
@@ -758,7 +787,7 @@ def test_set_st_fpm_inertia_back_to_original(
         # Make the blade tip heavier to change center of gravity.
         density_new *= np.linspace(0.8, 2.5, density_new.size)
         st.set_value(mset=1, set=1, m=density_new)
-        h2.set_st("blade1", st.main_data_sets[1][1])
+        h2.set_st(blade_id, st.main_data_sets[1][1])
         inertia_actual = h2.body_output_mass(0)
 
         # Must differ from the desired ones.
@@ -768,7 +797,7 @@ def test_set_st_fpm_inertia_back_to_original(
 
     # Restore ST.
     st.set_value(mset=1, set=1, m=mass_original)
-    h2.set_st("blade1", st.main_data_sets[1][1])
+    h2.set_st(blade_id, st.main_data_sets[1][1])
     inertia_actual = h2.body_output_mass(0)
     for i in range(4):  # Loop over tuple of arrays.
         npt.assert_allclose(inertia_actual[i], inertia_desired[i])
@@ -786,7 +815,8 @@ def test_set_st_fpm_blade_inertia_changed_st(
 
     # Set ST in the original blade, thus making it equivalent to the changed one.
     _, st_changed = write_iea22mw_only_blade_high_gravity_changed_st
-    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    blade_id = h2_iea22mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_iea22mw_only_blade_high_gravity.set_st(blade_id, st_changed.main_data_sets[1][1])
 
     # Compare inertia properties for all bodies.
     nbdy, _ = h2_iea22mw_only_blade_high_gravity_changed_st.get_number_of_bodies_and_constraints()
@@ -799,7 +829,7 @@ def test_set_st_fpm_blade_inertia_changed_st(
 
     # Restore ST.
     _, _, st_original = write_iea22mw_only_blade_high_gravity
-    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+    h2_iea22mw_only_blade_high_gravity.set_st(blade_id, st_original.main_data_sets[1][1])
 
 
 def test_set_st_fpm_blade_element_changed_st(
@@ -814,7 +844,8 @@ def test_set_st_fpm_blade_element_changed_st(
 
     # Set ST in the original blade, thus making it equivalent to the changed one.
     _, st_changed = write_iea22mw_only_blade_high_gravity_changed_st
-    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_changed.main_data_sets[1][1])
+    blade_id = h2_iea22mw_only_blade_high_gravity.get_mainbody_name_dict()["blade1"]
+    h2_iea22mw_only_blade_high_gravity.set_st(blade_id, st_changed.main_data_sets[1][1])
 
     # Compare element matrices for all bodies.
     nelem = h2_iea22mw_only_blade_high_gravity_changed_st.get_number_of_elements()
@@ -829,7 +860,7 @@ def test_set_st_fpm_blade_element_changed_st(
 
     # Restore ST.
     _, _, st_original = write_iea22mw_only_blade_high_gravity
-    h2_iea22mw_only_blade_high_gravity.set_st("blade1", st_original.main_data_sets[1][1])
+    h2_iea22mw_only_blade_high_gravity.set_st(blade_id, st_original.main_data_sets[1][1])
 
 
 def test_set_orientation_relative_main_body_not_found(
@@ -964,5 +995,5 @@ def test_set_orientation_relative_static(
 # %% Main.
 
 if __name__ == "__main__":
-    # pytest.main([__file__])
-    pytest.main([__file__, "-k test_set_st_fpm_blade_element_changed_st"])
+    pytest.main([__file__])
+    # pytest.main([__file__, "-k test_set_st_fpm_blade_element_changed_st"])

{h2lib_tests-13.2.801.dist-info → h2lib_tests-13.2.901.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.4
 Name: h2lib-tests
-Version: 13.2.801
-Summary: Tests and test_files for test h2lib (13.2.8)
+Version: 13.2.901
+Summary: Tests and test_files for test h2lib (13.2.9+4-gfd77867)
 Download-URL: 
 Author: Mads M Pedersen
 Author-email: 

{h2lib_tests-13.2.801.dist-info → h2lib_tests-13.2.901.dist-info}/RECORD

@@ -2,7 +2,7 @@ h2lib_tests/__init__.py,sha256=VjSqfGg8BzdmSjfSFhJh4hZbYZ_cME7xp9EWFKHQphA,61
 h2lib_tests/conftest.py,sha256=r4p868zCfJ5MMd2hs6XvGaiqu3jYlSR5G0OKAk5WX7Y,29629
 h2lib_tests/dtu10mw.py,sha256=a7SXfyDwDQPastYKb5CgghOQcYfgO1eTwGrd-H3Enok,4374
 h2lib_tests/test_calc.py,sha256=VNLfr2J9R2Jy9xTbdZ9dfbQ4dCwr7H7nbZRI3yP69fQ,2152
-h2lib_tests/test_distributed_sections.py,sha256=-YqB2GgRWmZMophwubRjRh8y5X6ptAm4xqor9sgDQ30,9042
+h2lib_tests/test_distributed_sections.py,sha256=xYcQwNGL3IQyVaguHGZrm0EZGYIMw4LRq9q5xwird1E,19413
 h2lib_tests/test_ellipsys_couplings.py,sha256=mfSTk1IamaKETU6Be8p8W9vA-yostWJl17m5sFlEK3Q,3345
 h2lib_tests/test_h2lib.py,sha256=t2jwmXzibej9DQoAkldwWX1MmEmpcV0kp2hLDRV5dbQ,13836
 h2lib_tests/test_h2rotor.py,sha256=cF3KMdzTyGoZGyihhHCJqdH_n_t2vjS1BPHDaUCk9pU,13918
@@ -10,7 +10,7 @@ h2lib_tests/test_lin.py,sha256=yydbsMX44ym_MSqR1cbQkqil8qC4icSW_JccPSElGac,6292
 h2lib_tests/test_mpi.py,sha256=CTT160yc6uZFRr_QNFWxyOwJ-y0qHiIp8jPYs1MQpyg,7111
 h2lib_tests/test_multiprocessinterface.py,sha256=h2o4havtK6IuMXsplNjGUa3VxOnbpEYGxdrrAKQilj0,1470
 h2lib_tests/test_static_solver.py,sha256=IuUkKrEoz_EKLCmGQ94CLFgjKPvO20iVx4773JTKO-8,6503
-h2lib_tests/test_topology_h2lib.py,sha256=g6H4E5OR4b2NlCDTZkOkUds2NyA-9XSG9Vu9Cy7TJw8,40348
+h2lib_tests/test_topology_h2lib.py,sha256=8SCrrjd8VZh_1pX9KDm_b2aCCyCZ-YbjwRgBRCOKNwQ,42472
 h2lib_tests/test_files/__init__.py,sha256=9e6ZUPb42e0wf2E1rutdcTM8hROcWFRVPXtZriU3ySw,50
 h2lib_tests/test_files/my_test_cls.py,sha256=7ZDsFkxrLfOY6q00U5Y-daxfuhATK-K5H04RP-VmQdE,850
 h2lib_tests/test_files/DTU_10_MW/control/dtu_we_controller.dll,sha256=C5T_CuAFtIuDgCXSYAoNu24yKPwj2nWOeORacJbLN9s,1134592
@@ -91,7 +91,7 @@ h2lib_tests/test_files/minimal/res/minimal_mann_turb.hdf5,sha256=Q3cs3bZyplZjBpo
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_u,sha256=byiorJmXDL6uKFbyfXthHTjJdm6ELvLR2lS202KrhRI,1048576
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_v,sha256=cxK5Rfgfm3gyJsEYi_KlmYY8DIIl_G0aizN2jt18Glc,1048576
 h2lib_tests/test_files/minimal/turb/hawc2_mann_l33.6_ae0.1000_g3.9_h0_512xd32xd16_2.000x3.00x4.00_s0001_w,sha256=xs61jAwhP3fIR1P5Oa8ovEt2baLoF8uCNs6pKIT8L4o,1048576
-h2lib_tests-13.2.801.dist-info/METADATA,sha256=c2VITTY0PxlaMfuXZrhOsoTeB8MVqNudNxuApGiQAPI,441
-h2lib_tests-13.2.801.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-h2lib_tests-13.2.801.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
-h2lib_tests-13.2.801.dist-info/RECORD,,
+h2lib_tests-13.2.901.dist-info/METADATA,sha256=tsWn93wEcdivKQ20oO9utThtj-uQRsQWDZmGR_khWj0,452
+h2lib_tests-13.2.901.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+h2lib_tests-13.2.901.dist-info/top_level.txt,sha256=WufAL3LO35YJBhWg1AfgTjSld-6l_WuRkXAkNKczUrM,12
+h2lib_tests-13.2.901.dist-info/RECORD,,