WEIS v2.1.2

.github/workflows/CI_WEIS.yml

@@ -17,7 +17,7 @@ jobs:
       fail-fast: false
       matrix:
         os: ["ubuntu-latest", "macOS-14", "windows-latest"] # mac13 is intel, mac-14 arm
-        python-version: ["3.11","3.12","3.13"]
+        python-version: ["3.12","3.13","3.14"]
 
     steps:
       - uses: actions/checkout@v5
@@ -50,7 +50,7 @@ jobs:
       # - name: Add dependencies windows specific
       #   if: contains( matrix.os, 'windows')
       #   run: |
-      #     conda install -y gfortran libpython
+      #     conda install -y gfortran gcc 
       #     pip install --no-deps control
 
 
@@ -62,12 +62,12 @@ jobs:
 
       # Don't delete after this
 
-
+      # Note- could add mpich through conda, but takes forever to resolve the environment in testing
       - name: Add dependencies ubuntu specific
         if: contains( matrix.os, 'ubuntu')
         run: |
-            sudo apt install -y libglu1-mesa
-            conda install -c conda-forge -y petsc4py mpi4py pyoptsparse mpich
+            sudo apt install -y libglu1-mesa mpich
+            conda install -y petsc4py mpi4py pyoptsparse
             pip install pygmsh==7.1.17
             pip install https://github.com/LHEEA/meshmagick/archive/master.zip  
             which -a mpiexec
@@ -156,10 +156,10 @@ jobs:
       - name: Test model creation notebooks and drivers
         if: contains( matrix.os, 'ubuntu') && contains( github.event_name, 'pull_request')
         run: |
-          cd examples/11_model_creation_process
-          treon 0_notebooks/chapter1.ipynb
-          treon 0_notebooks/chapter2.ipynb
-          treon 0_notebooks/chapter3.ipynb
+          cd examples/11_model_creation_process/0_notebooks
+          treon chapter1.ipynb
+          treon chapter2.ipynb
+          treon chapter3.ipynb
 
       # Run parallel script calling OpenFAST
       - name: Run parallel cases

README.md

@@ -70,7 +70,7 @@ The installation instructions below use the environment name, "weis-env," but an
 
 3. If you are on Kestrel, first load some modules and then install:
 
-        module load intel-oneapi-compilers intel-oneapi-mpi intel-oneapi-mkl conda
+        module load intel-oneapi-compilers intel-oneapi-mpi intel-oneapi-mkl
         pip install --no-deps -e . -v
 
 **NOTE:** To use WEIS again after installation is complete, you will always need to activate the conda environment first with `conda activate weis-env` (or `source activate weis-env`). On Kestrel, make sure to reload the necessary modules

environment.yml

@@ -5,36 +5,35 @@ channels:
 dependencies:
   - bs4
   - dash
+  - dash-bootstrap-components
+  - dash-mantine-components
   - dill
   - jsonmerge
   - mat4py
+  - moorpy==1.2.1
   - nlopt
   - numpydoc
   - numpy
   - openfast==4.2
   - openfast-io==4.2
   - pyopenfast==4.2
   - openraft>=2.0.3
-  - openmdao==3.42
-  - osqp  
+  - openmdao
+  - osqp
   - pcrunch>=2.1.5
   - pip
   - pyhams>=1.3
+  - pyvista
   #- pyoptsparse
+  - pyvista
   - rosco>=2.10.1
   - trimesh
-  - wisdem>=4.1
+  - wisdem>=4.2
   - pip:
     - control
-    - dash-bootstrap-components
-    - dash-mantine-components
     - dash-vtk
-    - pyvista
     - dearpygui
-    - orbit-nrel
     - smt
-    - windIO>=2.0.1
-    - wombat>=0.13.1
     # - git+https://github.com/NLRWindSystems/RAFT.git@dev
 # Needs to be done outside of environment file:
 #  - m2w64-toolchain # [win]

examples/03_design_with_openfast/iea22_ptfm_opt_analysis.yaml

@@ -73,7 +73,7 @@ constraints:
             lower_bound: 0.37
     floating:
         survival_heel:
-            upper_bound: 0.17453292519943295 # 10 deg
+            upper_bound: 10.0 # deg
         metacentric_height:
             flag: False
             lower_bound: 15.0

examples/04_frequency_domain_analysis_design/iea22_raft_opt_analysis.yaml

@@ -19,9 +19,6 @@ design_variables:
                   lower_bound: -30.0
                   upper_bound: -12.0
                   flag: True
-                # - names: [col1_freeboard, col2_freeboard, col3_freeboard]
-                #   lower_bound: 10.0
-                #   upper_bound: 20.0
             r_coordinate:
                 - names: [col1_keel, col1_freeboard, col2_keel, col2_freeboard, col3_keel, col3_freeboard]
                   lower_bound: 50.0
@@ -35,18 +32,7 @@ design_variables:
                       lower_bound: 12.0
                       upper_bound: 16.0
                       constant: True
-                #   ballast:
-                #       lower_bound: 0
-                #       upper_bound: 1000
-                # - names: [main_column]
-                #   ballast:
-                #       lower_bound: 0
-                #       upper_bound: 1000
-                # - names: [Y_pontoon_lower1, Y_pontoon_lower2, Y_pontoon_lower3]
-                #   diameter:
-                #       lower_bound: 10.0
-                #       upper_bound: 10.8
-                #       constant: True
+
 constraints:
     tower:
         height_constraint:
@@ -73,7 +59,7 @@ constraints:
             lower_bound: 0.37
     floating:
         survival_heel:
-            upper_bound: 0.17453292519943295 # 10 deg
+            upper_bound: 10.0 # deg
         metacentric_height:
             flag: True
             lower_bound: 15.0
@@ -91,14 +77,6 @@ constraints:
             flag: True
         buoyancy:
             flag: False
-        stress:
-            flag: False
-        global_buckling:
-            flag: False
-        shell_buckling:
-            flag: False
-        mooring_heel:
-            flag: False
         freeboard_margin:   # keep freeboard from being submerged below water during survival_heel, largest wave
             flag: True
         draft_margin:   # keep draft from raising above water line during survival_heel, largest wave
@@ -109,9 +87,6 @@ constraints:
         Max_PtfmPitch:
             flag: True
             max: 6.0
-        Std_PtfmPitch:
-            flag: False
-            max: 1.25  # Same as IEA-15MW with same DLCs
         nacelle_acceleration:
             flag: True
             max: 2.0

examples/11_model_creation_process/0_notebooks/chapter2.ipynb

@@ -502,7 +502,7 @@
     "constraints:\n",
     "    floating:\n",
     "        survival_heel:\n",
-    "            upper_bound: 0.17453292519943295 # 10 deg\n",
+    "            upper_bound: 10.0 # deg\n",
     "        metacentric_height:\n",
     "            flag: True\n",
     "            lower_bound: 1.0 # 15.0 --> Dan's experience\n",

examples/11_model_creation_process/1_drivers/stage-3-semisub/stage-3-semisub_of_analysis.yaml

@@ -53,7 +53,7 @@ constraints:
             lower_bound: 0.37
     floating:
         survival_heel:
-            upper_bound: 0.17453292519943295 # 10 deg
+            upper_bound: 10.0 # deg
         metacentric_height:
             flag: True
             lower_bound: 15.0

examples/11_model_creation_process/1_drivers/stage-3-semisub/stage-3-semisub_raft_analysis.yaml

@@ -53,7 +53,7 @@ constraints:
             lower_bound: 0.37
     floating:
         survival_heel:
-            upper_bound: 0.17453292519943295 # 10 deg
+            upper_bound: 10.0 # deg
         metacentric_height:
             flag: True
             lower_bound: 15.0

examples/11_model_creation_process/1_drivers/stage-3.5-semisubCCD/stage-3.5-semisub_of_analysis.yaml

@@ -72,7 +72,7 @@ constraints:
             lower_bound: 0.37
     floating:
         survival_heel:
-            upper_bound: 0.17453292519943295 # 10 deg
+            upper_bound: 10.0 # deg
         metacentric_height:
             flag: True
             lower_bound: 15.0

weis/aeroelasticse/openmdao_openfast.py

@@ -625,6 +625,7 @@ def setup(self):
 
         # Floating outputs
         self.add_output('Max_PtfmPitch', val=0.0, desc='Maximum platform pitch angle over a set of OpenFAST simulations')
+        self.add_output('Mean_PtfmPitch', val=0.0, units='deg', desc='Maximum (across cases) mean (of each case) platform pitch angle over a set of OpenFAST simulations')
         self.add_output('Std_PtfmPitch', val=0.0, units='deg', desc='standard deviation of platform pitch angle')
         self.add_output('Max_Offset', val=0.0, units='m', desc='Maximum distance in surge/sway direction')
 
@@ -1200,9 +1201,18 @@ def update_FAST_model(self, fst_vt, inputs, discrete_inputs):
                 if not np.any(inputs[f'{fass}_modes'][idir,:]):
                     logger.warning(f'WARNING: {fass} tower shape coefficients are zero which will cause errors in using ElastoDyn')
         fst_vt['ElastoDynTower']['TwFAM1Sh'] = inputs['fore_aft_modes'][0, :]  / np.sum(inputs['fore_aft_modes'][0, :])
-        fst_vt['ElastoDynTower']['TwFAM2Sh'] = inputs['fore_aft_modes'][1, :]  / np.sum(inputs['fore_aft_modes'][1, :])
         fst_vt['ElastoDynTower']['TwSSM1Sh'] = inputs['side_side_modes'][0, :] / np.sum(inputs['side_side_modes'][0, :])
+
+        # Since the 2nd tower modes are sometimes problematic, if the DOF is not enabled, let's give it a safe, dummy value that won't cause errors in ElastoDyn
+        fst_vt['ElastoDynTower']['TwFAM2Sh'] = inputs['fore_aft_modes'][1, :]  / np.sum(inputs['fore_aft_modes'][1, :])
+        if not fst_vt['ElastoDyn']['TwFADOF2']:
+            fst_vt['ElastoDynTower']['TwFAM2Sh'] = np.zeros_like(inputs['fore_aft_modes'][1, :])
+            fst_vt['ElastoDynTower']['TwFAM2Sh'][0] = 1.0
+
         fst_vt['ElastoDynTower']['TwSSM2Sh'] = inputs['side_side_modes'][1, :] / np.sum(inputs['side_side_modes'][1, :])
+        if not fst_vt['ElastoDyn']['TwSSDOF2']:
+            fst_vt['ElastoDynTower']['TwSSM2Sh'] = np.zeros_like(inputs['side_side_modes'][1, :])
+            fst_vt['ElastoDynTower']['TwSSM2Sh'][0] = 1.0
 
         # Calculate yaw stiffness of tower (springs in series) and use in servodyn as yaw spring constant
         k_tow_tor = inputs['tor_stff'] / np.diff(inputs['tower_z'])
@@ -3450,6 +3460,7 @@ def get_floating_measures(self, inputs, outputs):
         calculate floating measures:
             - Std_PtfmPitch (max over all dlcs if constraint, mean otheriwse)
             - Max_PtfmPitch
+            - Mean_PtfmPitch
 
         given:
             - sum_stats : pd.DataFrame
@@ -3463,6 +3474,7 @@ def get_floating_measures(self, inputs, outputs):
             outputs['Std_PtfmPitch'] = np.mean(sum_stats['PtfmPitch']['std'])
 
         outputs['Max_PtfmPitch']  = np.max(sum_stats['PtfmPitch']['max'])
+        outputs['Mean_PtfmPitch'] = np.max(sum_stats['PtfmPitch']['mean'])
 
         # Max platform offset        
         outputs['Max_Offset'] = np.max(sum_stats['PtfmOffset']['max'])

weis/glue_code/gc_PoseOptimization.py

@@ -260,7 +260,14 @@ def set_constraints(self, wt_opt):
                 raise Exception('Please turn on the call to OpenFAST or RAFT if you are trying to optimize Max_PtfmPitch constraints.')
             wt_opt.model.add_constraint(f'{self.floating_solve_component}.Max_PtfmPitch',
                 upper = control_constraints['Max_PtfmPitch']['max'])
-
+
+        # Mean platform pitch
+        if control_constraints['Mean_PtfmPitch']['flag']:
+            if not any(self.level_flags):
+                raise Exception('Please turn on the call to OpenFAST or RAFT if you are trying to optimize Mean_PtfmPitch constraints.')
+            wt_opt.model.add_constraint(f'{self.floating_solve_component}.Mean_PtfmPitch',
+                upper = control_constraints['Mean_PtfmPitch']['max'])
+
         # Platform pitch motion
         if control_constraints['Std_PtfmPitch']['flag']:
             if not any(self.level_flags):

weis/inputs/analysis_schema.yaml

@@ -382,9 +382,21 @@ properties:
                                 minimum: 0.0
                                 maximum: 30.0
                                 unit: deg
+                    Mean_PtfmPitch:
+                        type: object
+                        description: The maximum mean platform pitch displacement over all cases. Can be computed in both RAFT and OpenFAST.  The higher fidelity option will be used when active.
+                        default: {}
+                        properties:
+                            flag: *flag
+                            max:
+                                type: number
+                                default: 3.0
+                                minimum: 0.0
+                                maximum: 30.0
+                                unit: deg
                     Std_PtfmPitch:
                         type: object
-                        description: Maximum platform pitch standard deviation over all cases. Can be computed in both RAFT and OpenFAST.  The higher fidelity option will be used when active. 
+                        description: Maximum platform pitch standard deviation over all cases. Can be computed in both RAFT and OpenFAST.  The higher fidelity option will be used when active.
                         default: {}
                         properties:
                             flag: *flag

weis/inputs/validation.py

@@ -106,4 +106,40 @@ def write_analysis_yaml(instance, foutput):
     sfx_str = "-analysis.yaml"
     wisval.write_yaml(instance, foutput + sfx_str)
     return foutput + sfx_str
-
+
+
+def make_paths_absolute(data, base_dir=None):
+    """Recursively convert relative paths in a nested dict/list to absolute paths.
+
+    Any string value that contains ``/`` or ``\\`` and is not already an
+    absolute path is joined with *base_dir* and resolved to an absolute path.
+    This is useful when loading YAML configuration files that contain
+    relative references to other files or directories.
+
+    Parameters
+    ----------
+    data : dict, list, or str
+        The data structure to process (typically the result of
+        ``yaml.safe_load``).
+    base_dir : str, optional
+        The directory that relative paths are relative to.  Defaults to
+        the current working directory.
+
+    Returns
+    -------
+    data : dict, list, or str
+        A copy of the input with relative paths replaced by absolute paths.
+    """
+    if base_dir is None:
+        base_dir = os.getcwd()
+
+    if isinstance(data, dict):
+        return {k: make_paths_absolute(v, base_dir) for k, v in data.items()}
+    elif isinstance(data, list):
+        return [make_paths_absolute(item, base_dir) for item in data]
+    elif isinstance(data, str):
+        if ("/" in data or "\\" in data) and not os.path.isabs(data):
+            return os.path.abspath(os.path.join(base_dir, data))
+        return data
+    else:
+        return data