updates for mooring adjuster function

-project.addMooring now inputs adjuster function, method, target, and i_line
-if target is not provided, the pretension or horizontal force will be calculated based on the mooring dd and base depth
-slope is calculated based on mooring span and base depth (only used for pretension option)
-mooring.reposition function updated to work with lasted adjustMooring helper function
-if slope is not provided, defaults to the 'horizontal' method and keeps same anchor spacing
-if slope is provided, the pretension method is used
-helper adjustMooring updated to use 'horizontal' as the default (requires less inputs)

Author: erickaloz

famodel/helpers.py

@@ -776,7 +776,7 @@ def getAnchors(lineAnch, arrayAnchor, proj):
 
     return(ad, mass)
 
-def adjustMooring(mooring, method = 'pretension', r=[0,0,0], project=None, target=1e6,
+def adjustMooring(mooring, method = 'horizontal', r=[0,0,0], project=None, target=1e6,
                        i_line = 0, slope = 0.58 ):
     '''Custom function to adjust a mooring, called by
     Mooring.adjust. Fairlead point should have already

famodel/mooring/mooring.py

@@ -251,7 +251,27 @@ def reposition(self, r_center=None, heading=None, project=None,
         # Run custom function to update the mooring design (and anchor position)
         # this would also szie the anchor maybe?
         if self.adjuster:
-            self.adjuster(self, r_centerB, u, project=project, **kwargs)
+            
+            #if i_line is not defined, assumed segment 0 will be adjusted
+            if not hasattr(self,'i_line'):
+                self.i_line = 0
+       
+            
+            if hasattr(self,'slope'):
+                self.adjuster(self, method = 'pretension', r=r_centerB, project=project, target = self.target, i_line = self.i_line, slope = self.slope)
+            
+            else:
+                
+                #move anchor based on set spacing then adjust line length
+                xy_loc = r_centerB[:2] + (self.span + rad_fair[1])*u
+                if project:
+                    self.dd['zAnchor'] = -project.getDepthAtLocation(xy_loc[0],xy_loc[1])
+                    self.z_anch = self.dd['zAnchor']
+                else:
+                    print('Warning: depth of mooring line, anchor, and subsystem must be updated manually.')
+                self.setEndPosition(np.hstack([r_centerB[:2] + (self.span + rad_fair[1])*u, self.z_anch]), 'a', sink=True)
+
+                self.adjuster(self, method = 'horizontal', r=r_centerB, project=project, target = self.target, i_line = self.i_line)
 
         elif self.shared == 1: # set position of end A at platform end A
             self.setEndPosition(np.hstack([r_centerA[:2] - rad_fair[0]*u, z_fair[0] + r_centerA[2]]),'a')

famodel/project.py

@@ -1484,7 +1484,8 @@ def addPlatform(self,r=[0,0,0], id=None, phi=0, entity='',
 
     def addMooring(self, id=None, endA=None, endB=None, heading=0, dd={}, 
                    section_types=[], section_lengths=[], 
-                   connectors=[], span=0, shared=0, reposition=False):
+                   connectors=[], span=0, shared=0, reposition=False, adjuster=None,
+                   method = 'horizontal', target = None, i_line = 0):
         '''
         Function to create a mooring object  and save in mooringList
         Optionally does the following:
@@ -1521,7 +1522,14 @@ def addMooring(self, id=None, endA=None, endB=None, heading=0, dd={},
             The default is 0. Used to develop a dd, unused if dd provided.
         shared : int, optional
             Describes if a mooring line is shared (1), sahred half line (2) or anchored (0). The default is 0.
-
+        adjuster : Function, optional
+            Function to adjust mooring lines for different depths. The default is none
+        method : str, optional
+            Method 'horizontal' or 'pretension' for adjuster function
+        target : float, optional
+            Target pretension or horizontal tension for adjuster function. if none, will calculate
+        i_line: int, optional
+            The index of the line segment to adjust. default is 0
         Returns
         -------
         mooring : mooring object
@@ -1556,6 +1564,20 @@ def addMooring(self, id=None, endA=None, endB=None, heading=0, dd={},
                   'z_fair':self.platformList[id_part[0]].zFair if id_part else 0}
 
         mooring = Mooring(dd=dd, id=id) # create mooring object
+        
+        #calculate target pretension or horizontal tension if none provided
+        if adjuster != None and target == None:
+            targetdd = deepcopy(dd)
+            targetdd['zAnchor'] = self.depth
+            mooring.createSubsystem(dd)
+            
+            if method == 'horizontal':
+                mooring.target = np.linalg.norm(mooring.ss.fB_L[:2])
+            elif method =='pretension':
+                mooring.target = np.linalg.norm(mooring.ss.fB_L)
+            else:
+                raise Exception('Invalid adjustment method. Must be pretension or horizontal')
+        
         # update shared prop if needed
         if len(id_part)==2 and shared<1:
             shared = 1
@@ -1582,7 +1604,25 @@ def addMooring(self, id=None, endA=None, endB=None, heading=0, dd={},
                 mooring.z_anch = -zAnew
         else:
             mooring.heading = np.degrees(heading)
+        
+        #add mooring adjuster if porivded
+        if adjuster!= None:
+            mooring.adjuster = adjuster 
+            mooring.i_line = i_line 
+            
+            # check if method is 'pretension' then save slope
+            if method == 'pretension':
 
+                if dd:
+                    
+                    #calculate mooring slope using base depth
+                    #**** this assumes that the mooring system is designed for the base depth*****
+                    mooring.slope = self.depth / dd['span']
+                    
+                else:   
+                    
+                    mooring.slope = self.depth / span
+        
         self.mooringList[id] = mooring
 
         return(mooring)