h2lib-tests 13.1.506__py3-none-any.whl → 13.1.1701__py3-none-any.whl

h2lib_tests/conftest.py

@@ -0,0 +1,373 @@
+import pytest
+import numpy as np
+from wetb.hawc2.htc_file import HTCFile
+from h2lib_tests.test_files import tfp
+from h2lib._h2lib import H2Lib
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_tower():
+    # Start from DTU_10MW_RWT and delete everything except the tower.
+    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc.set_name("DTU_10MW_RWT_only_tower")
+    for key1 in htc["new_htc_structure"].keys():
+        if key1.startswith("main_body"):
+            if "tower" not in htc["new_htc_structure"][key1]["name"].values:
+                htc["new_htc_structure"][key1].delete()
+        if key1 == "orientation":
+            for key2 in htc["new_htc_structure"]["orientation"].keys():
+                if key2.startswith("relative"):
+                    htc["new_htc_structure"]["orientation"][key2].delete()
+        if key1 == "constraint":
+            for key2 in htc["new_htc_structure"]["constraint"].keys():
+                if key2 != "fix0":
+                    htc["new_htc_structure"]["constraint"][key2].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+    # Reduce simulation time.
+    htc.simulation.time_stop = 10.0
+    # Change number of bodies in the tower.
+    htc.new_htc_structure.main_body.nbodies = 3
+    # Save the new file.
+    htc.save()
+    return htc
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_tower_rotated(write_dtu10mw_only_tower):
+    # Start from the DTU_10MW_RWT_only_tower and rotate the tower.
+    htc = write_dtu10mw_only_tower.copy()
+    htc.set_name("DTU_10MW_RWT_only_tower_rotated")
+    alpha = 30.0
+    htc.new_htc_structure.orientation.base.body_eulerang = [
+        alpha,
+        0.0,
+        0.0,
+    ]
+    htc["wind"].delete()
+    htc["output"].delete()
+    htc.save()
+    return (htc, alpha)
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_tower_encrypted(write_dtu10mw_only_tower):
+    # Start from the DTU_10MW_RWT_only_tower and then encrypt the tower.
+    htc = write_dtu10mw_only_tower.copy()
+    htc.set_name("DTU_10MW_RWT_only_tower_encrypted")
+    # Only the tower is left.
+    htc.new_htc_structure.main_body.timoschenko_input.filename = "./data/DTU_10MW_RWT_Tower_st.dat.v3.enc"
+    htc.save()
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_blade():
+    # Start from DTU_10MW_RWT and delete everything except the blade.
+    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc.set_name("DTU_10MW_RWT_only_blade")
+    for key1 in htc["new_htc_structure"].keys():
+        if key1.startswith("main_body"):
+            if "blade1" not in htc["new_htc_structure"][key1]["name"].values:
+                htc["new_htc_structure"][key1].delete()
+        if key1 == "orientation":
+            htc["new_htc_structure"][key1].delete()
+        if key1 == "constraint":
+            htc["new_htc_structure"][key1].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+
+    # Set the blade horizontal, to maximize gravity loading.
+    htc.new_htc_structure.add_section("orientation")
+    htc.new_htc_structure.orientation.add_section("base")
+    htc.new_htc_structure.orientation.base.mbdy = "blade1"
+    htc.new_htc_structure.orientation.base.inipos = [0.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base["mbdy_eulerang"] = [90.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base.mbdy_eulerang.comments = "Blade span is horizontal."
+
+    # Clamp the blade.
+    htc.new_htc_structure.add_section("constraint")
+    htc.new_htc_structure.constraint.add_section("fix0")
+    htc.new_htc_structure.constraint.fix0.mbdy = "blade1"
+
+    # Set as many bodies as elements.
+    htc.new_htc_structure.main_body__7.nbodies = 26
+
+    htc.simulation.log_deltat.delete()
+
+    # Do not use static solver, since it will be done during the test.
+    htc.simulation.solvertype = 2
+    htc.simulation.solvertype.comments = ""
+    htc.simulation.initial_condition = 1
+
+    # No output, as we will use add_sensor().
+
+    # Save the new file.
+    htc.save()
+
+    return htc
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_blade_low_max_iter(write_dtu10mw_only_blade):
+    # Start from the write_dtu10mw_only_blade and then reduce the number of max iterations,
+    # so that the static solver will not have time to converge.
+    htc = write_dtu10mw_only_blade.copy()
+    htc.set_name("DTU_10MW_RWT_only_blade_low_max_iter")
+    htc.simulation.max_iterations = 1
+    htc.save()
+
+    return htc
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_blade_rotate_base(write_dtu10mw_only_blade):
+    # Start from the write_dtu10mw_only_blade and then make it rotate by using the base command.
+    # HAWC2 will use the initial condition, but then the blade will not rotate because of the fix0 constraint.
+    # So, running the simulation will show a clamped blade that vibrates.
+    # Rotate at about 9.5 rpm.
+    htc = write_dtu10mw_only_blade.copy()
+    htc.set_name("DTU_10MW_RWT_only_blade_rotate_base")
+    speed = 1.0  # [rad/s]
+    htc.new_htc_structure.orientation.base["mbdy_ini_rotvec_d1"] = [0.0, 1.0, 0.0, speed]
+    htc.new_htc_structure.orientation.base.mbdy_ini_rotvec_d1.comments = f"= {speed * 30.0 / np.pi:.2f} rpm"
+    htc.new_htc_structure.main_body__7.gravity = 0.0
+    htc["wind"].delete()
+    htc["output"].delete()
+    htc.save()
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_blade_rotate_relative():
+    # Start from DTU_10MW_RWT and delete everything except the blade and hub.
+    # The blade now rotates because of the relative rotation.
+    # Because of the fix1 constraint, the blade will not rotate after time 0.
+    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc.set_name("DTU_10MW_RWT_only_blade_rotate_relative")
+    for key1 in ("main_body",     # tower
+                 "main_body__2",  # towertop
+                 "main_body__3",  # shaft
+                 "main_body__5",  # hub2
+                 "main_body__6",  # hub3
+                 "main_body__8",  # blade2
+                 "main_body__9",  # blade3
+                 "orientation",
+                 "constraint",
+                 ):
+        htc["new_htc_structure"][key1].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+
+    # Set the hub as vertical.
+    htc.new_htc_structure.add_section("orientation")
+    htc.new_htc_structure.orientation.add_section("base")
+    htc.new_htc_structure.orientation.base.mbdy = "hub1"
+    htc.new_htc_structure.orientation.base.inipos = [0.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base["mbdy_eulerang"] = [180.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base.mbdy_eulerang.comments = "Hub axis is up."
+
+    # Set the blade horizontal.
+    htc.new_htc_structure.orientation.add_section("relative")
+    htc.new_htc_structure.orientation.relative.mbdy1 = "hub1  last"
+    htc.new_htc_structure.orientation.relative.mbdy2 = "blade1  1"
+    htc.new_htc_structure.orientation.relative["mbdy2_eulerang"] = [-90.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.relative.mbdy2_eulerang.comments = "Blade span is horizontal."
+    speed = 1.0  # [rad/s]
+    htc.new_htc_structure.orientation.relative["mbdy2_ini_rotvec_d1"] = [0.0, 1.0, 0.0, speed]
+    htc.new_htc_structure.orientation.relative.mbdy2_ini_rotvec_d1.comments = f"= {speed * 30.0 / np.pi:.2f} rpm"
+
+    # Disable gravity.
+    htc.new_htc_structure.main_body__7.gravity = 0.0
+    htc.new_htc_structure.main_body__4.gravity = 0.0
+
+    # Clamp the hub and blade.
+    htc.new_htc_structure.add_section("constraint")
+    htc.new_htc_structure.constraint.add_section("fix0")
+    htc.new_htc_structure.constraint.fix0.mbdy = "hub1"
+    htc.new_htc_structure.constraint.add_section("fix1")
+    htc.new_htc_structure.constraint.fix1.mbdy1 = "hub1  last"
+    htc.new_htc_structure.constraint.fix1.mbdy2 = "blade1  1"
+
+    # Set as many bodies as elements.
+    htc.new_htc_structure.main_body__7.nbodies = 26
+
+    htc.simulation.log_deltat.delete()
+
+    # Do not use static solver, since it will be done during the test.
+    htc.simulation.solvertype = 2
+    htc.simulation.solvertype.comments = ""
+    htc.simulation.initial_condition = 1
+
+    # Save the new file.
+    htc.save()
+
+    return htc
+
+
+@pytest.fixture(scope="session")
+def write_dtu10mw_only_blade_rotate_bearing3():
+    # Start from DTU_10MW_RWT and delete everything except the blade and hub.
+    # The blade now rotates because of the bearing3 between the blade and hub.
+    htc = HTCFile(tfp + "DTU_10_MW/htc/DTU_10MW_RWT.htc")
+    htc.set_name("DTU_10MW_RWT_only_blade_rotate_bearing3")
+    for key1 in ("main_body",     # tower
+                 "main_body__2",  # towertop
+                 "main_body__3",  # shaft
+                 "main_body__5",  # hub2
+                 "main_body__6",  # hub3
+                 "main_body__8",  # blade2
+                 "main_body__9",  # blade3
+                 "orientation",
+                 "constraint",
+                 ):
+        htc["new_htc_structure"][key1].delete()
+    htc["wind"].delete()
+    htc["aerodrag"].delete()
+    htc["aero"].delete()
+    htc["dll"].delete()
+    htc["output"].delete()
+
+    # Set the hub as vertical.
+    htc.new_htc_structure.add_section("orientation")
+    htc.new_htc_structure.orientation.add_section("base")
+    htc.new_htc_structure.orientation.base.mbdy = "hub1"
+    htc.new_htc_structure.orientation.base.inipos = [0.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base["mbdy_eulerang"] = [180.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.base.mbdy_eulerang.comments = "Hub axis is up."
+
+    # Set the blade horizontal.
+    htc.new_htc_structure.orientation.add_section("relative")
+    htc.new_htc_structure.orientation.relative.mbdy1 = "hub1  last"
+    htc.new_htc_structure.orientation.relative.mbdy2 = "blade1  1"
+    htc.new_htc_structure.orientation.relative["mbdy2_eulerang"] = [-90.0, 0.0, 0.0]
+    htc.new_htc_structure.orientation.relative.mbdy2_eulerang.comments = "Blade span is horizontal."
+
+    # Disable gravity.
+    htc.new_htc_structure.main_body__7.gravity = 0.0
+    htc.new_htc_structure.main_body__4.gravity = 0.0
+
+    # Clamp the hub.
+    htc.new_htc_structure.add_section("constraint")
+    htc.new_htc_structure.constraint.add_section("fix0")
+    htc.new_htc_structure.constraint.fix0.mbdy = "hub1"
+
+    # Insert bearing3.
+    htc.new_htc_structure.constraint.add_section("bearing3")
+    htc.new_htc_structure.constraint.bearing3.name = "bearing"
+    htc.new_htc_structure.constraint.bearing3.mbdy1 = "hub1 last"
+    htc.new_htc_structure.constraint.bearing3.mbdy2 = "blade1 1"
+    htc.new_htc_structure.constraint.bearing3.bearing_vector = [1, 0.0, 0.0, 1.0]
+    speed = 1.0  # [rad/s]
+    htc.new_htc_structure.constraint.bearing3.omegas = speed
+    htc.new_htc_structure.constraint.bearing3.omegas.comments = f"= {speed * 30.0 / np.pi:.2f} rpm"
+
+    # Set as many bodies as elements.
+    htc.new_htc_structure.main_body__7.nbodies = 26
+
+    htc.simulation.log_deltat.delete()
+
+    # Do not use static solver, since it will be done during the test.
+    htc.simulation.solvertype = 2
+    htc.simulation.solvertype.comments = ""
+    htc.simulation.initial_condition = 1
+
+    # Save the new file.
+    htc.save()
+
+    return htc
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_tower(write_dtu10mw_only_tower):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_tower_rotated(write_dtu10mw_only_tower_rotated):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower_rotated.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_tower_encrypted(write_dtu10mw_only_tower_encrypted):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_tower_encrypted.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_blade(write_dtu10mw_only_blade):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu10mw_only_blade_low_max_iter(write_dtu10mw_only_blade_low_max_iter):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade_low_max_iter.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_blade_rotate_base(write_dtu10mw_only_blade_rotate_base):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade_rotate_base.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_blade_rotate_relative(write_dtu10mw_only_blade_rotate_relative):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade_rotate_relative.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()
+
+
+@pytest.fixture(scope="session")
+def h2_dtu_10mw_only_blade_rotate_bearing3(write_dtu10mw_only_blade_rotate_bearing3):
+    h2 = H2Lib(suppress_output=True)
+    model_path = f"{tfp}DTU_10_MW/"
+    htc_path = "htc/DTU_10MW_RWT_only_blade_rotate_bearing3.htc"
+    h2.init(htc_path=htc_path, model_path=model_path)
+    h2.stop_on_error(False)
+    yield h2
+    h2.close()

h2lib_tests/test_ellipsys_couplings.py

@@ -19,7 +19,7 @@ class Ellipsys():
         return self.time
 
     def set_fxyz(self, pos_xyz, fxyz):
-        print(self.time, np.shape(pos_xyz), np.shape(fxyz))
+        # print(self.time, np.shape(pos_xyz), np.shape(fxyz))
         pass
 
 
@@ -28,7 +28,7 @@ def test_ellipsys_dummy_workflow_1wt():
     N = 1
     with MultiH2Lib(N, suppress_output=True) as h2:
         el = Ellipsys()
-        htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT.htc')
+        htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
         for i in range(N):
             htc.set_name(f'wt{i}')
             htc.save()
@@ -41,7 +41,7 @@ def test_ellipsys_dummy_workflow_1wt():
             pos_gl_xyz = np.array(h2.get_aerosections_position(), order='F') + wt_pos[na, na, :]
             uvw = np.asfortranarray(el.get_uvw(pos_gl_xyz))
             h2.set_aerosections_windspeed(uvw)
-            h2.run(t)  # run after set_aero_windspeed requires initilalize of bem before disabling
+            h2.run(t)  # run after set_aero_windspeed requires initialize of bem before disabling
             frc_gl_xyz = h2.get_aerosections_forces()
             el.set_fxyz(pos_gl_xyz, frc_gl_xyz)
             if t == 1:
@@ -53,7 +53,7 @@ def test_ellipsys_dummy_workflow():
     N = 4
     with MultiH2Lib(N, suppress_output=True) as h2:
         el = Ellipsys()
-        htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT.htc')
+        htc = HTCFile(tfp + 'DTU_10_MW/htc/DTU_10MW_RWT_no_aerodrag.htc')
         for i in range(N):
             htc.set_name(f'wt{i}')
             htc.save()

h2lib_tests/test_files/DTU_10_MW/htc/DTU_10MW_RWT.htc

@@ -6,7 +6,7 @@ begin simulation;
   on_no_convergence continue ;
   convergence_limits 1E3 1.0 1E-7 ;
   ;logfile ./log/DTU_10MW_RWT_ver09.log ;
-  log_deltat .1;
+  log_deltat 1;
   begin newmark;
     deltat    0.01;  
   end newmark;
@@ -599,7 +599,7 @@ end dll;
 ;
 begin output;
   filename ./res/DTU_10MW_RWT_ver09 ;
-  data_format  gtsdf;
+  data_format  gtsdf64;
   buffer 500 ;
 ;
   general time;