Modifications in anchor class to include mass_update flag and minor capacity model examples updates

Author: Moreno

famodel/anchors/anchor.py

@@ -634,20 +634,18 @@ def getCapacityAnchor(self, Hm, Vm, zlug, line_type=None, d=None, w=None, mass_u
             self.anchorCapacity['Rotational displacement'] = results['Rotational displacement']
 
         # Weight calculated via dimensions
-        if mass_update == False:
+        if not mass_update:
             if 'Weight pile' in results:
                 self.anchorCapacity['Weight pile'] = results['Weight pile']
             if 'Weight plate' in results:
                 self.anchorCapacity['Weight plate'] = results['Weight plate']
         else:
             if 'Weight pile' in results:
-                if self.mass is None:
-                    self.mass = results['Weight pile']/self.g
-                self.anchorCapacity['Weight pile'] = self.mass*self.g
+                self.mass = results['Weight pile']/self.g
+                self.anchorCapacity['Weight pile'] = results['Weight pile']
             if 'Weight plate' in results:
-                if self.mass is None:
-                    self.mass = results['Weight plate']/self.g
-                self.anchorCapacity['Weight plate'] = self.mass*self.g
+                self.mass = results['Weight plate']/self.g
+                self.anchorCapacity['Weight plate'] = results['Weight plate']
 
         # print(f"[DEBUG] Stored Lateral displacement in anchorCapacity: {self.anchorCapacity['Lateral displacement']:.6f}")
 
@@ -1312,9 +1310,17 @@ def getSafetyFactor(self):
 
             return {'SF_combined': SF} 
 
-    def getCostAnchor(self, ms=None):
+    def getCostAnchor(self, ms=None, mass_update=True):
         '''
         Assign material cost using a Point object and getCost_and_MBL().
+        
+        Parameters
+        ----------
+        ms : MoorPy System, optional
+            The mooring system to which the anchor point belongs. If None, a new one is created.
+        mass_update : bool, optional
+            If True, update mpAnchor mass from self.mass.
+            If False, preserve existing mpAnchor.m if already set.
         '''
 
         # Create or use existing MoorPy system
@@ -1324,12 +1330,12 @@ def getCostAnchor(self, ms=None):
             # Create MoorPy Point using makeMoorPyAnchor
             self.makeMoorPyAnchor(ms)
 
-        # Check if mass is assigned
-        if self.mass is None:
+        # Assign self.mass if missing
+        if self.mass is None or mass_update:
             if 'Weight pile' in self.anchorCapacity:
-                self.mass = self.anchorCapacity['Weight pile'] / self.g
+                self.mass = self.anchorCapacity['Weight pile']/self.g
             elif 'Weight plate' in self.anchorCapacity:
-                self.mass = self.anchorCapacity['Weight plate'] / self.g
+                self.mass = self.anchorCapacity['Weight plate']/self.g
             else:
                 raise KeyError("Missing 'Weight pile' or 'Weight plate' in anchorCapacity. \
                 Run getCapacityAnchor() before getCostAnchor(), or define self.mass explicitly.")

famodel/anchors/anchors_famodel/capacity_dandg.py

@@ -220,7 +220,7 @@ def PileWeight(Len, Dia, tw, rho):
     Ha = 5.0e6        # Horizontal load (N)
     Va = 3.0e5        # Vertical load (N)
 
-    layers, y, z, results = getCapacityDandG(profile_map, 'CPT_rock_1', L, D, zlug, Ha, Va, plot=True)
+    layers, y, z, results = getCapacityDandG(profile_map, 'CPT_rock_1', L, D, zlug, Ha, Va, plot=True, display=0)
 
     print('\n--- Results for DandG Pile in Layered Rock ---')
     for key, val in results.items():

famodel/anchors/anchors_famodel/capacity_load.py

@@ -203,7 +203,7 @@ def getTransferLoad(profile_map, Tm, thetam, zlug, line_type, d, w=None, plot=Fa
     d = 0.25              # Chain diameter (m)
     w = 5000              # Line weight (N/m)
 
-    layers, resultsLoad = getTransferLoad(profile_map, Tm, thetam, zlug, line_type, d, w, plot=True)
+    layers, resultsLoad = getTransferLoad(profile_map, Tm, thetam, zlug, line_type, d, w, plot=True, display=1)
 
     plot_load(layers, resultsLoad['drag_values'], resultsLoad['depth_values'], 
               resultsLoad['Tm'], resultsLoad['thetam'], resultsLoad['Ta'], 

famodel/anchors/anchors_famodel/capacity_suction.py

@@ -504,7 +504,7 @@ def vertical_cross(H, M, H_target):
     zlug = (2/3)*L                    # Lug depth (m)
     theta = 5                         # Tilt misalignment angle (deg)
     psi = 7.5                         # Twist misalignment angle (deg)
-    Ha = 1e6                        # Applied horizontal load (N)
+    Ha = 1e6                          # Applied horizontal load (N)
     Va = 5.7e6                        # Applied vertical load (N)
 
     # Calculate
@@ -513,7 +513,8 @@ def vertical_cross(H, M, H_target):
         D, L, zlug,                   # Pile geometrical properties
         Ha, Va,                       # Pile loading conditions   
         thetalug=theta, psilug=psi,   # Pile misaligment tolerances 
-        plot=False
+        plot=True,
+        display=1
     )
 
     # print('\n--- Suction Pile Capacity Results ---')