Update suction pile examples + anchor and project logic

Author: Moreno

examples/05_Anchors/anchor_soil_layered.py

@@ -10,21 +10,17 @@
     soil_mode='layered',
     profile_source='inputs/GulfOfMaine_soil_profiles.yaml')
 
-for label, props in proj.soilProps.items():
-    print(f"{label}: {props}")
-
 # Step 2: Create and register an anchor at a known position in the grid
 anchor = Anchor(
     dd = {'type': 'suction', 'design': {'D': 3.5, 'L': 12.0, 'zlug': 9.67}},
     r  = [54.0, -4450.0, 0.0])
 
 # Step 3: Assign local soil profile from project (nearest neighbor lookup)
-soil_id, soil_profile = proj.getSoilAtLocation(anchor.r[0], anchor.r[1])
-anchor.soilProps = {soil_id: soil_profile}
-anchor.setSoilProfile([{'name': soil_id, 'layers': soil_profile}])  # ensures `anchor.soil_profile` is set
+proj.getSoilAtLocation(anchor.r[0], anchor.r[1])
+anchor.setSoilProfile(proj.profile_map)  
 
 # Step 4: Assign loads and line
-anchor.loads = {'Hm': 2e6, 'Vm': 1.5e6}
+anchor.loads = {'Hm': 1e6, 'Vm': 1.5e6}
 anchor.line_type = 'chain'
 anchor.d = 0.16
 anchor.w = 5000.0
@@ -52,7 +48,7 @@
     loads = None,
     lambdap_con = [3, 6],
     zlug_fix = False,
-    safety_factor = {'SF_combined': 1},
+    safety_factor = {'SF_combined': 2},
     plot = True)
 
 # Step 6: Report

examples/05_Anchors/anchor_soil_uniform.py

@@ -6,27 +6,23 @@
 proj = Project()
 proj.loadSoil(
     filename='inputs/GulfOfMaine_soil_layered_100x100.txt',
-    soil_mode='uniform')  # 'uniform' soil does not need from the profile_source yaml file
-
-for label, props in proj.soilProps.items():
-    print(f"{label}: {props}")
+    soil_mode='uniform')  # 'uniform' soil mode does not need from the profile_source yaml file
 
-# Convert to profile_map format so anchor capacity models can use it
-proj.convertUniformToLayered(default_layer=50.0)
-
-for label, props in proj.profile_map.items():
-    print(f"{label}: {props}")
-
 # Step 2: Create and register an anchor at a known position in the grid
 anchor = Anchor(
     dd = {'type': 'suction', 'design': {'D': 3.5, 'L': 12.0, 'zlug': 9.67}},
     r  = [54.0, -4450.0, 0.0])
 
-# Step 3: Assign local soil profile from project (nearest neighbor lookup)
-soil_id, _ = proj.getSoilAtLocation(anchor.r[0], anchor.r[1])
-soil_profile = proj.profile_map[soil_id]  # get compatible layered format
-anchor.soilProps = {soil_id: soil_profile}
-anchor.setSoilProfile([{'name': soil_id, 'layers': soil_profile}])  # ensures `anchor.soil_profile` is set
+# SELECT option 3a or 3b
+# ======================
+# Step 3a: Convert to profile_map format from uniform soil data
+# proj.convertUniformToLayered(default_layer=50.0)
+# profile_map = proj.profile_map[3]     # manually take relevant soil profile from the list
+# anchor.setSoilProfile(profile_map)    # automatically turns dict to list 
+
+# Step 3b: Assign local soil profile from project (nearest neighbor lookup)
+proj.getSoilAtLocation(anchor.r[0], anchor.r[1])
+anchor.assignSoilProfile(proj.profile_map)   
 
 # Step 4: Assign loads and line
 anchor.loads = {'Hm': 1e6, 'Vm': 5e4}

examples/05_Anchors/anchor_suction.py

@@ -52,7 +52,7 @@
     r = [250.0, 250.0, 000.0])
 
 # --- Step 0: Create anchor based grid CPTs ---
-anchor.interpolateSoilProfile(profile_map)
+anchor.assignSoilProfile(profile_map)
 
 # --- Step 1: Plot suction pile and soil profile ---
 # Access anchor geometrical properties
@@ -67,8 +67,8 @@
 
 # Assign loads manually
 anchor.loads = {
-    'Hm': 3.0e6,    # Horizontal mudline load (N)
-    'Vm': 1.0e6     # Vertical mudline load (N)
+    'Hm': 2.0e6,    # Horizontal mudline load (N)
+    'Vm': 8.0e5     # Vertical mudline load (N)
 }
 
 # Assign line properties manually
@@ -109,7 +109,6 @@
 anchor.getCost()  
 
 print(f"Mass: {anchor.anchorCapacity['Weight pile']/9.81:.2f} kg")
-print(f"Material unit cost: {anchor.cost['unit_cost']:.2f} USD/kg")
 print(f'Material cost: {anchor.cost["Material cost"]:.2f} USD [2024]')
 
 #%%
@@ -121,7 +120,7 @@
     loads = None,
     lambdap_con = [3, 6],
     zlug_fix = False,
-    safety_factor = {'SF_horizontal': 2, 'SF_vertical': 3},
+    safety_factor = {'SF_combined': 2},
     plot = True)
 
 print('\nFinal Optimized Anchor:')
@@ -132,8 +131,5 @@
 anchor.getCost()  
 
 print(f"Mass: {anchor.anchorCapacity['Weight pile']/9.81:.2f} kg")
-print(f"Material unit cost: {anchor.cost['unit_cost']:.2f} USD/kg")
 print(f'Material cost: {anchor.cost["Material cost"]:.2f} USD [2024]')
 
-# --- Step 7: Visualize Anchor Geometry ---
-# anchor.getCombinedPlot()

examples/05_Anchors/example_suction.ipynb

@@ -130,6 +130,14 @@ def setSoilProfile(self, profile_map):
         Assign a soil profile directly from a single CPT.
         Assumes profile_map is a list with only one entry.
         '''
+        
+        if isinstance(profile_map, dict):
+            profile_map = [profile_map]    
+        elif isinstance(profile_map, list):
+            pass  # leave it as it is
+        else:
+            raise TypeError('setSoilProfile expects a dict or a list/tuple with one dict.')
+            
         if len(profile_map) != 1:
             raise ValueError("setSoilProfile expects a profile_map with exactly one CPT.")
 
@@ -215,6 +223,42 @@ def interpolateSoilProfile(self, profile_map):
 
         if self.display > 0: print(f"[Anchor] Interpolated soil profile: {self.profile_name} with soil types {self.soil_type_list}")
 
+    def assignSoilProfile(self, profile_map):
+        '''
+        High-level soil assignment:
+
+        - If a single CPT is provided → call setSoilProfile
+        - If 4 or more CPTs are provided → call interpolateSoilProfile
+        - Otherwise → complain (not enough info for interpolation)
+
+        Accepts:
+          - dict (single CPT)
+          - list of CPT dicts
+        '''
+        # Normalize input
+        if isinstance(profile_map, dict):
+            profile_list = [profile_map]
+        elif isinstance(profile_map, list):
+            profile_list = profile_map
+        else:
+            raise TypeError('applySoilProfile expects a dict or a list of dicts.')
+
+        n = len(profile_list)
+
+        if n == 1:
+            # delegate to existing set logic
+            return self.setSoilProfile(profile_list)
+
+        elif n >= 4:
+            # delegate to existing interpolation logic
+            return self.interpolateSoilProfile(profile_list)
+
+        else:
+            # 2 or 3 CPTs: ambiguous use-case
+            raise ValueError(
+                f'applySoilProfile received {n} CPTs: '
+                'need exactly 1 for direct assignment or ≥4 for interpolation.')
+        
     def makeMoorPyAnchor(self, ms):
         '''
         Create a MoorPy anchor object in a MoorPy system.

famodel/anchors/anchor.py

@@ -1257,7 +1257,7 @@ def loadSoil(self, filename=None, yaml=None, soil_mode='uniform', profile_source
         else:
             raise ValueError("Invalid combination of filename/yaml inputs")
 
-        # --- Set defaults only for uniform mode ---
+        # --- Set defaults only for uniform mode (when values are missing) ---
         if soil_mode == 'uniform':
             for key, props in soilProps.items():
                 props['Su0']   = getFromDict(props, 'Su0',   shape=-1, dtype=list, default=[2.39])
@@ -1298,24 +1298,34 @@ def getSoilAtLocation(self, x, y):
         -------
         (str, dict or list): soil name or profile ID, and associated soil properties or layered profile
         '''
+        self.profile_map = []
+        
         if self.soil_x is not None:
             ix = np.argmin([abs(x - sx) for sx in self.soil_x])
             iy = np.argmin([abs(y - sy) for sy in self.soil_y])
             soil_id = self.soil_names[iy, ix]  # could be label or profile_id
 
             if self.soil_mode == 'uniform':
-                soil_info = self.soilProps[soil_id]
-                return soil_id, soil_info
-
+                soil_info = self.soilProps[soil_id]               
+                if not hasattr(self, 'profile_map') or not self.profile_map:
+                    self.convertUniformToLayered(default_layer=50.0)
+                    
+                # Replace with a single entry corresponding to this soil_id
+                self.profile_map = [
+                    next(e for e in self.profile_map if e['name'] == str(soil_id))]
+                self.profile_map[0]['layers']
+                
             elif self.soil_mode == 'layered':
-                profile_layers = self.soilProps[soil_id]  # list of layer dicts
-                return soil_id, profile_layers
-
+                layers = self.soilProps[soil_id]  # list of layer dicts
+                profile_entry = {'name': str(soil_id), 'layers': layers}
+                self.profile_map.append(profile_entry)
+                
             else:
                 raise ValueError(f"Unknown soil_mode: {self.soil_mode}")
 
             print(f"[DEBUG] soil_id at location ({x}, {y}) is: {soil_id}")
             print(f"[DEBUG] Available soilProps keys: {list(self.soilProps.keys())}")
+        
         else:
             raise ValueError("No soil grid defined")
 
@@ -1331,8 +1341,8 @@ def convertUniformToLayered(self, default_layer=50.0):
             name = str(name)
 
             gamma = float(props['gamma'][0])
-            Su0 = float(props['Su0'][0])
-            k = float(props['k'][0])
+            Su0   = float(props['Su0'][0])
+            k     = float(props['k'][0])
 
             layer = {
                 'soil_type': 'clay',
@@ -1341,7 +1351,7 @@ def convertUniformToLayered(self, default_layer=50.0):
                 'gamma_top': gamma,
                 'gamma_bot': gamma,
                 'Su_top': Su0,
-                'Su_bot': Su0 + k * default_layer}
+                'Su_bot': Su0 + k*default_layer}
 
             profile_entry = {'name': name, 'layers': [layer]}
             self.profile_map.append(profile_entry)