Add anchor rerouting for cables, other small updates

Author: lsirkis

famodel/anchors/anchor.py

@@ -6,6 +6,7 @@
 from famodel.famodel_base import Node
 from famodel.mooring.mooring import Mooring
 import famodel.platform.platform 
+import shapely as sh
 
 
 class Anchor(Node):
@@ -81,6 +82,8 @@ def __init__(self, dd=None, ms=None, r=[0,0,0], aNum=None, id=None,
         # Dictionaries for additional information
         # anchor capacity
         self.anchorCapacity = {}
+        self.safety_factors = {} # calculated safety factor
+        self.safety_factors_required = {} # minimum allowable safety factor
 
         # anchor costs
         self.cost = {}
@@ -676,7 +679,10 @@ def getFS(self, loads=None, acceptance_crit=None):
         else:
             return(FS)
 
-        
+    def makeBuffer(self, buff_rad=50):
+        point = sh.Point(self.r[:2])
+        buff = point.buffer(buff_rad)
+        return buff
 
     def getCost(self,costDict='default'):
         '''find costs of anchor and store in design dictionary

famodel/cables/cable.py

@@ -4,6 +4,7 @@
 from copy import deepcopy
 from moorpy.subsystem import Subsystem
 from moorpy import helpers
+import shapely as sh
 
 from famodel.cables.dynamic_cable import DynamicCable
 from famodel.cables.static_cable import StaticCable
@@ -233,6 +234,24 @@ def getL(self):
             L += cab.L
 
         self.L = L
+        
+    def makeLine(self,buff_rad=20,include_dc=True):
+        
+        coords = []
+        for sub in self.subcomponents:
+            if isinstance(sub,Joint):
+                coords.append(sub['r'][:2])
+            elif isinstance(sub,DynamicCable) and include_dc:
+                coords.append(sub.rA[:2])
+                coords.append(sub.rB[:2])
+                pass
+            elif isinstance(sub, StaticCable):
+                if len(sub.x)>0:
+                    for i in range(len(sub.x)):
+                        coords.append([sub.x[i],sub.y[i]])
+        line = sh.LineString(coords)
+
+        return(line)
 
     def updateSpan(self,newSpan):
         '''

famodel/helpers.py

@@ -9,6 +9,7 @@
 import ruamel.yaml
 import moorpy as mp
 from moorpy.helpers import loadPointProps, getPointProps
+import shapely as sh
 
 
 def cart2pol(x, y):
@@ -211,7 +212,6 @@ def check_headings(m_headings,c_heading,rad_buff):
     # convert negative headings to positive headings
     for i,mh in enumerate(m_headings):
         if mh<0:
-            #breakpoint()
             m_headings[i] = 2*np.pi + mh
         elif mh>2*np.pi:
             m_headings[i] = mh - 2*np.pi
@@ -249,7 +249,6 @@ def head_adjust(att,heading,rad_buff=np.radians(30),endA_dir=1):
         headnew = np.pi*2 + heading
     else:
         headnew = heading
-    #breakpoint()
     attheadings = []
     flipheads = False # whether to flip headings ( for if you are looking at mooring headings of platform on the other end)
     for at in att:
@@ -849,6 +848,57 @@ def getAnchors(lineAnch, arrayAnchor, proj):
 
     return(ad, mass)
 
+def route_around_anchors(proj, anchor=True, cable=True, padding=50):
+    
+    # make anchor buffers with 50m radius
+    if anchor:
+        anchor_buffs = []
+        for anch in proj.anchorList.values():
+            anchor_buffs.append(anch.makeBuffer())
+      
+    # make cable linestrings including joint locs and static cable (no dynamic cable for simplicity)
+    if cable:
+        cable_line = {}
+        for name, cab in proj.cableList.items():
+            cable_line[name] = cab.makeLine(include_dc=False)
+
+    # Function to calculate angle of a point relative to a center
+    def angle(pt):
+        return np.arctan2(pt[1] - center[1], pt[0] - center[0])    
+
+    # Loop through each cable linestring and anchor buffer
+    for name,cab in cable_line.items():
+        for anch in anchor_buffs:
+            if cab.intersects(anch):
+                # Get the start and end of the detour (the two closest points to the buffer)
+                segments = []
+                # make additional points on the line on either side of anchor
+                dist_to_anch = cab.line_locate_point(anch.centroid)
+                if dist_to_anch > 100:
+                    segments.append(cab.interpolate(dist_to_anch - 100))
+                if cab.length - dist_to_anch > 100:
+                    segments.append(cab.interpolate(dist_to_anch + 100))
+
+                start = np.array(segments[0].coords[-1])
+
+                # Get buffer center and radius
+                center = np.array(anch.centroid.coords[0])
+                radius = anch.boundary.distance(sh.Point(center))+padding
+
+                # Calculate angle for start point relative to center
+                angle_start = angle(start)
+
+                # Generate point along the arc (detour)
+                arc_point = [center[0] + radius * np.cos(angle_start+np.pi/2), center[1] + radius * np.sin(angle_start+np.pi/2)]
+
+                # add new routing point in cable object
+                proj.cableList[name].subcomponents[2].updateRouting([list(segments[0].coords[1:]) + arc_point + list(segments[1].coords[:-1])])
+
+
+
+
+
+            
 def configureAdjuster(mooring, adjuster=None, method='horizontal',
                       i_line=0, span=None, project=None, target=None):
     '''Configures adjuster function for mooring object

famodel/project.py

@@ -34,7 +34,8 @@
 # Import select required helper functions
 from famodel.helpers import (check_headings, head_adjust, getCableDD, getDynamicCables, 
                             getMoorings, getAnchors, getFromDict, cleanDataTypes, 
-                            getStaticCables, getCableDesign,m2nm, loadYAML, configureAdjuster)
+                            getStaticCables, getCableDesign, m2nm, loadYAML, 
+                            configureAdjuster, route_around_anchors)
 
 
 class Project():
@@ -1900,7 +1901,7 @@ def addCablesConnections(self,connDict,cableType_def='dynamic_cable_66',oss=Fals
             heads = [headingA,headingB]
             # reposition cable
             cab.reposition(project=self,headings=heads,rad_fair=[5,5])
-            
+
             coords = []
             if cableConfig:
                 ref_cables = None
@@ -1930,6 +1931,7 @@ def addCablesConnections(self,connDict,cableType_def='dynamic_cable_66',oss=Fals
                 # update routing
                 cab.subcomponents[cs].updateRouting(coords) # also updates static and general cable lengths
 
+        route_around_anchors(self)
 
 
     def updatePositions(self):
@@ -2538,7 +2540,7 @@ def getMoorPyArray(self, plt=0, pristineLines=True, cables=True):
                         pass
                     elif isinstance(comp,DynamicCable):
                         # create subsystem for dynamic cable
-                        comp.createSubsystem(pristine=pristineLines)  
+                        comp.createSubsystem(pristine=pristineLines) 
 
                         if pristineLines:                           
                             ssloc = comp.ss