Updated project; replaced anchor_soil example with layered/uniform variants

Author: Moreno

examples/05_Anchors/anchor_soil.py …amples/05_Anchors/anchor_soil_layered.pyexamples/05_Anchors/anchor_soil.py renamed to examples/05_Anchors/anchor_soil_layered.py examples/05_Anchors/anchor_soil.py renamed to examples/05_Anchors/anchor_soil_layered.py

@@ -23,7 +23,7 @@
 # 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
+anchor.setSoilProfile([{'name': soil_id, 'layers': soil_profile}])  # ensures `anchor.soil_profile` is set
 
 # Step 4: Assign loads and line
 anchor.loads = {'Hm': 2e6, 'Vm': 1.5e6}

examples/05_Anchors/anchor_soil_uniform.py

@@ -0,0 +1,73 @@
+
+import sys
+sys.path.append(r'C:\Code\FAModel_anchors\famodel')
+
+from project import Project
+from anchors.anchor import Anchor
+
+# Step 1: Initialize and load soil
+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}")
+    
+# 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
+
+# Step 4: Assign loads and line
+anchor.loads = {'Hm': 1e6, 'Vm': 5e4}
+anchor.line_type = 'chain'
+anchor.d = 0.16
+anchor.w = 5000.0
+
+# Step 5: Run capacity check and optimization
+anchor.getLugForces(
+    Hm=anchor.loads['Hm'], Vm=anchor.loads['Vm'], 
+    zlug = anchor.dd['design']['zlug'], 
+    d=anchor.d, w=anchor.w,  
+    plot=True)
+
+anchor.getCapacityAnchor(
+    Hm=anchor.loads['Hm'], Vm=anchor.loads['Vm'], 
+    zlug = anchor.dd['design']['zlug'], 
+    line_type=anchor.line_type, d=anchor.d, w=anchor.w,   
+    mass_update=True,
+    plot=True)
+anchor.getCostAnchor()
+print(f'Material cost: {anchor.cost["Material cost"]:.2f} USD [2024]')
+
+results = anchor.getSizeAnchor(
+    geom = [anchor.dd['design']['L'], anchor.dd['design']['D']],
+    geomKeys = ['L', 'D'],
+    geomBounds = [(8.0, 15.0), (2.0, 4.0)],
+    loads = None,
+    lambdap_con = [3, 6],
+    zlug_fix = False,
+    safety_factor = {'SF_combined': 1},
+    plot = True)
+
+# Step 6: Report
+print('\nFinal Optimized Anchor:')
+print('Design:', anchor.dd['design'])
+print('Capacity Results:', anchor.anchorCapacity)
+anchor.getCostAnchor()
+print(f'Material cost: {anchor.cost["Material cost"]:.2f} USD [2024]')
+
+

famodel/anchors/anchor.py

@@ -205,7 +205,7 @@ def interpolateSoilProfile(self, profile_map):
         self.soilProps = dict(soilProps)
 
         print(f"[Anchor] Interpolated soil profile: {self.profile_name} with soil types {self.soil_type_list}")
-
+        
     def makeMoorPyAnchor(self, ms):
         '''
         Create a MoorPy anchor object in a MoorPy system.

famodel/project.py

@@ -1166,6 +1166,31 @@ def getSoilAtLocation(self, x, y):
             print(f"[DEBUG] Available soilProps keys: {list(self.soilProps.keys())}")
         else:
             raise ValueError("No soil grid defined")
+            
+    def convertUniformToLayered(self, default_layer=50.0):
+        '''
+        Converts self.soilProps (uniform format) into profile_map (layered format)
+        using a default thickness and assuming uniform clay profile.
+        Matches the structure of layered CPT-based soil profiles.
+        '''
+        self.profile_map = {}
+
+        for name, props in self.soilProps.items():
+            name = str(name)
+            gamma = float(props['gamma'][0])
+            Su0 = float(props['Su0'][0])
+            k = float(props['k'][0])
+
+            layer = {
+                'soil_type': 'clay',
+                'top': 0.0,
+                'bottom': default_layer,
+                'gamma_top': gamma,
+                'gamma_bot': gamma,
+                'Su_top': Su0,
+                'Su_bot': Su0 + k * default_layer}
+
+            self.profile_map[name] = [layer]  # just layers!
 
     # # ----- Anchor 
     def updateAnchor(self,anch='all',update_loc=True):