Skip to content

Adding features to MD unload function to support FFarm integration: #42

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
mattEhall merged 3 commits into from
Apr 15, 2025
Merged

Adding features to MD unload function to support FFarm integration: #42

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
c5ecf97
Adding features to MD unload function to support FFarm integration:
Yuksel-Rudy Feb 28, 2025
a32465f
changing name from phiV to just phi and writing statement to check f…
Yuksel-Rudy Apr 15, 2025
6503ad3
fixing value error
Yuksel-Rudy Apr 15, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
43 changes: 34 additions & 9 deletions moorpy/system.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1020,7 +1020,7 @@ def parseYAML(self, data):


def unload(self, fileName, MDversion=2, line_dL=0, rod_dL=0, flag='p',
outputList=[], Lm=[0], T_half=42):
outputList=[], Lm=[0], T_half=42, phi=None, MDoptionsDict={}):
'''Unloads a MoorPy system into a MoorDyn-style input file

Parameters
Expand All @@ -1040,7 +1040,10 @@ def unload(self, fileName, MDversion=2, line_dL=0, rod_dL=0, flag='p',
T_half : float, optional
For tuning response of viscoelastic model, the period when EA is
half way between the static and dynamic values [s]. Default is 42.

phi : list of floats, optional
platform's rotations variables [deg]. Defaults to None and platforms are not unrotated.
MDoptionsDict: dictionary, optional
MoorDyn Options. If not given, default options are considered.
Returns
-------
None.
Expand All @@ -1052,7 +1055,8 @@ def unload(self, fileName, MDversion=2, line_dL=0, rod_dL=0, flag='p',
#Collection of default values, each can be customized when the method is called

# Set up the dictionary that will be used to write the OPTIONS section
MDoptionsDict = dict(dtM=0.001, kb=3.0e6, cb=3.0e5, TmaxIC=60) # start by setting some key default values
if not MDoptionsDict:
MDoptionsDict = dict(dtM=0.001, kb=3.0e6, cb=3.0e5, TmaxIC=60) # start by setting some key default values
# Other available options: Echo=False, dtIC=2, CdScaleIC=10, threshIC=0.01
MDoptionsDict.update(self.MDoptions) # update the dict with any settings saved from an input file
MDoptionsDict.update(dict(g=self.g, WtrDepth=self.depth, rho=self.rho)) # lastly, apply any settings used by MoorPy
Expand Down Expand Up @@ -1125,13 +1129,20 @@ def unload(self, fileName, MDversion=2, line_dL=0, rod_dL=0, flag='p',
point_type = 'Fixed'

#Check if the point is attached to body
for body in self.bodyList:
for i, body in enumerate(self.bodyList):
for attached_Point in body.attachedP:

if attached_Point == point.number:
#point_type = "Body" + str(body.number)
point_type = "Vessel"
point_pos = body.rPointRel[body.attachedP.index(attached_Point)] # get point position in the body reference frame
point_type = "Vessel"
if phi:
c, s = np.cos(np.radians(-phi[i])), np.sin(np.radians(-phi[i]))
R = np.array([[c, -s, 0],
[s, c, 0],
[0, 0, 1]])
point_pos = R @ body.rPointRel[body.attachedP.index(attached_Point)]
else:
point_pos = body.rPointRel[body.attachedP.index(attached_Point)] # get point position in the body reference frame

elif point.type == 0: # point is coupled externally (con, free)
point_type = 'Connect'
Expand Down Expand Up @@ -1208,8 +1219,14 @@ def unload(self, fileName, MDversion=2, line_dL=0, rod_dL=0, flag='p',
#version =
#description =

# Check length of phi if given
if phi:
if len(phi) != len(self.bodyList):
raise ValueError(f"Inconsistency in the size of the given platform rotation angle phi. There are {len(self.bodyList)} bodies in MS but only {len(phi)} angles are given.")

# Set up the dictionary that will be used to write the OPTIONS section
MDoptionsDict = dict(dtM=0.001, kb=3.0e6, cb=3.0e5, TmaxIC=60) # start by setting some key default values
if not MDoptionsDict:
MDoptionsDict = dict(dtM=0.001, kb=3.0e6, cb=3.0e5, TmaxIC=60) # start by setting some key default values
MDoptionsDict.update(self.MDoptions) # update the dict with any settings saved from an input file
MDoptionsDict.update(dict(g=self.g, depth=self.depth, rho=self.rho)) # lastly, apply any settings used by MoorPy

Expand Down Expand Up @@ -1407,11 +1424,19 @@ def unload(self, fileName, MDversion=2, line_dL=0, rod_dL=0, flag='p',
point_type = 'Fixed'

#Check if the point is attached to body
for body in self.bodyList:
for i, body in enumerate(self.bodyList):
for attached_Point in body.attachedP:

if attached_Point == point.number:
point_type = "Body" + str(body.number)
point_pos = body.rPointRel[body.attachedP.index(attached_Point)] # get point position in the body reference frame
if phi:
c, s = np.cos(np.radians(-phi[i])), np.sin(np.radians(-phi[i]))
R = np.array([[c, -s, 0],
[s, c, 0],
[0, 0, 1]])
point_pos = R @ body.rPointRel[body.attachedP.index(attached_Point)]
else:
point_pos = body.rPointRel[body.attachedP.index(attached_Point)] # get point position in the body reference frame

elif point.type == 0: # point is coupled externally (con, free)
point_type = 'Free'
Expand Down
Loading