Skip to content

Updates to Multiphysics and Teamerlupa2 #3

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

Open
hmankle wants to merge 7 commits into
base: main
Choose a base branch
from
Open

Updates to Multiphysics and Teamerlupa2 #3

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
298da1c
Adding LUPA WEC-Sim Model with MoorDyn
hmankle Nov 27, 2024
a376bf6
Merge branch 'main' of https://github.com/PMEC-OSU/LUPA_WEC-Sim
hmankle Dec 3, 2024
7fafbd7
Revert "Adding LUPA WEC-Sim Model with MoorDyn"
hmankle Dec 3, 2024
0f800db
Merge branch 'PMEC-OSU:main' into main
hmankle Dec 11, 2024
ca4e5a7
Updating Damping loss terms and pto equilibirum position for Multiphy…
hmankle Sep 3, 2025
90108ca
Merge branch 'main' of https://github.com/hmankle/LUPA_WEC-Sim
hmankle Sep 3, 2025
180f8e4
Merge pull request #1 from hmankle/LUPA_updates
hmankle Sep 3, 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
1 change: 1 addition & 0 deletions .gitignore
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
/Additional Numerical Models/TEAMERLUPA2_inf_depth_20m/TwoBodySixDOF
92 changes: 0 additions & 92 deletions Additional Numerical Models/TEAMERLUPA2_inf_depth_20m/TwoBodySixDOF/wecSimInputFile.asv
View file
Open in desktop

This file was deleted.

12 changes: 6 additions & 6 deletions LUPATEAMER2 Fall 2023/TwoBodySixDOF/wecSimInputFile.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,7 +3,7 @@
%
% Depth is set at 3.69 m
%
% Updated by: Hannah Mankle 6/3/2025
% Updated by: Hannah Mankle 9/3/2025


%% Simulation Data
Expand Down Expand Up @@ -48,10 +48,8 @@
body(2).mass = 202.21; % [kg] Positively bouyant spar.
body(2).viz.color = [211/256 211/256 211/256];
body(2).inertia = [304.786 305.250 15.775]; % [kg-m^2] As measured from dry swing tests
% body(2).initial.axis = [1 1 0];
% body(2).initial.angle = -0.0010;
% body(2).linearDamping(4,4) = 3;
% body(2).linearDamping(5,5) = 3;
body(2).linearDamping(4,4) = 200; % Linear damping to help stablize model at beginning of simulation
body(2).linearDamping(5,5) = 200;
body(2).quadDrag.cd = [0.6 0.6 2.8 0.6 0.6 2.8]; % [-] Quadratic drag coefficient Cd as found from Beatty 2015 and Singh & Mittal 2005
body(2).quadDrag.area = [0.546 0.546 0.636 0.546 0.546 0.636]; % [m^2] Characteristic area in relevant plane
body(2).setInitDisp([0 0 0],[0 0 0 0],[0 0 -0.42]); % [m] Initial Displacement Set to engage mooring lines for pre-tension.
Expand All @@ -67,8 +65,10 @@
pto(1) = ptoClass('PTO1'); % Initialize PTO Class for PTO1
pto(1).stiffness = 0; % [N/m] PTO Stiffness
pto(1).damping = 0; % [N/(m/s)] PTO Damping Typical values fall between 0-10,000 N/(m/s)
ptoDampingLoss = -350; % [N/m/s] Linear damping Found experimentally through free decay tests. It is negative because it is a electromechanical loss.
ptoDampingLoss = 900; % [N/m/s]
pto(1).location = [0 0 0]; % [m] PTO Location
pto(1).equilibriumPosition = -0.42; % [m] initial PTO location, coresponses with spar displacement & needed if model is using stiffness terms.


%% Mooring Matrix
R = 0.325; % [m] Radius of mooring plate
Expand Down
12 changes: 4 additions & 8 deletions MULTIPHYSICSLUPA3 Spring 2024/TwoBodyHeaveOnly_D0_9m/wecSimInputFile.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@
% Heave plate diameter is 0.9 m
% Depth is set at 2.78 m
%
% Updated by: Hannah Mankle 6/6/2025
% Updated by: Hannah Mankle 9/3/2025


%% Simulation Data
Expand Down Expand Up @@ -51,16 +51,11 @@
body(2).mass = 202.21; % [kg] Positively bouyant spar.
body(2).viz.color = [211/256 211/256 211/256];
body(2).inertia = [304.786 305.250 15.775]; % [kg-m^2] As measured from dry swing tests
% body(2).initial.axis = [1 1 0];
% body(2).initial.angle = -0.0010;
% body(2).linearDamping(4,4) = 3;
% body(2).linearDamping(5,5) = 3;
body(2).linearDamping(3,3) = 700;
body(2).quadDrag.cd = [0.6 0.6 2.8 0.6 0.6 2.8]; % [-] Quadratic drag coefficient Cd as found from Beatty 2015 and Singh & Mittal 2005
body(2).quadDrag.area = [0.546 0.546 0.636 0.546 0.546 0.636]; % [m^2] Characteristic area in relevant plane
body(2).setInitDisp([0 0 0],[0 0 0 0],[0 0 -0.42]); % [m] Initial Displacement Set to engage mooring lines for pre-tension.

% add linear damping to pitch and roll

%% PTO and Constraint Parameters
% Translational Joint
constraint(1) = constraintClass('Constraint1'); % Initialize Constraint Class for Constraint1
Expand All @@ -70,8 +65,9 @@
pto(1) = ptoClass('PTO1'); % Initialize PTO Class for PTO1
pto(1).stiffness = 0; % [N/m] PTO Stiffness
pto(1).damping = 0; % [N/(m/s)] PTO Damping Typical values fall between 0-10,000 N/(m/s)
ptoDampingLoss = -350; % [N/m/s] Linear damping Found experimentally through free decay tests. It is negative because it is a electromechanical loss.
ptoDampingLoss = -900; % [N/m/s] Linear damping Found experimentally through free decay tests. It is negative because it is a electromechanical loss.
pto(1).location = [0 0 0]; % [m] PTO Location
% pto(1).location = -0.42; % [m] initial PTO location, coresponses with spar displacement & needed if model is using stiffness terms.

%% Mooring Matrix
R = 0.325; % [m] Radius of mooring plate
Expand Down
6 changes: 4 additions & 2 deletions MULTIPHYSICSLUPA3 Spring 2024/TwoBodyHeaveOnly_D1_14m/wecSimInputFile.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@
% Heave plate diameter is 1.14 m
% Depth is set at 2.78 m
%
% Updated by: Hannah Mankle 6/6/2025
% Updated by: Hannah Mankle 9/3/2025

%% Simulation Data
simu=simulationClass();
Expand Down Expand Up @@ -43,6 +43,7 @@
body(2).mass = 202.21; % [kg] Positively bouyant spar.
body(2).viz.color = [211/256 211/256 211/256];
body(2).inertia = [304.786 305.250 15.775]; % [kg-m^2] As measured from dry swing tests
body(2).linearDamping(3,3) = 700; % Based on 0.9 m heave plate, value may need to be changed
body(2).quadDrag.cd = [0.6 0.6 2.8 0.6 0.6 2.8]; % [-] Quadratic drag coefficient Cd as found from Beatty 2015 and Singh & Mittal 2005
body(2).quadDrag.area = [0.558 0.558 0.636 0.558 0.558 0.636]; % [m^2] Characteristic area in relevant plane
body(2).setInitDisp([0 0 0],[0 0 0 0],[0 0 -0.22]); % [m] Initial Displacement Set to engage mooring lines for pre-tension.
Expand All @@ -56,8 +57,9 @@
pto(1) = ptoClass('PTO1'); % Initialize PTO Class for PTO1
pto(1).stiffness = 0; % [N/m] PTO Stiffness
pto(1).damping = 0; % [N/(m/s)] PTO Damping Typical values fall between 0-10,000 N/(m/s)
ptoDampingLoss = -350; % [N/m/s] Linear damping Found experimentally through free decay tests. It is negative because it is a electromechanical loss.
ptoDampingLoss = -900; % [N/m/s] Linear damping loss. It is negative because it is a electromechanical loss. Value determined from validated 0.9 m heave plate and may need updating.
pto(1).location = [0 0 0]; % [m] PTO Location
% pto(1).location = -0.22; % [m] initial PTO location, coresponses with spar displacement & needed if model is using stiffness terms.

%% Mooring Matrix
R = 0.325; % [m] Radius of mooring plate
Expand Down