Simulating Floating Platforms
In this example, we make the OWENS model more complex by simulation a floating wind turbine instead of a land-based one. This involves coupling two meshes together: one representing the wind turbine (this will be similar to meshes used in previous tutorials), and one representing the floating platform.
The wind turbine mesh (or "topside") will be created similarly to the previous tutorial, though we will overwrite some of the inputs so the simulation works better for a floating system. This will help with revealing some more internals of OWENS and the potential inputs available to users, which is helpful for explaining how the floating platform mesh (or "bottom side") is defined and differs from the topside.
This example is also available as a Jupyter notebook: D_simulatingFloatingPlatforms.ipynb.
# import PyPlot
import OWENS
import OWENSFEA
import OWENSAero
runpath = path = "/home/runner/work/OWENS.jl/OWENS.jl/examples/literate"
##runpath = path = splitdir(@__FILE__)[1] # use to run locally
nothingFirst, create the inputs for the topside mesh as done previous in tutorials A and B.
Inp = OWENS.MasterInput("$runpath/sampleOWENS.yml")
analysisType = "TNB" # this differs from previous tutorials, as it has been verified to work well with floating capabilities
turbineType = Inp.turbineType
eta = Inp.eta
Nbld = Inp.Nbld
towerHeight = Inp.towerHeight
rho = Inp.rho
Vinf = Inp.Vinf
controlStrategy = Inp.controlStrategy
RPM = Inp.RPM
Nslices = Inp.Nslices
ntheta = Inp.ntheta
structuralModel = Inp.structuralModel
ntelem = Inp.ntelem
nbelem = Inp.nbelem
ncelem = Inp.ncelem
nselem = Inp.nselem
ifw = Inp.ifw
WindType = Inp.WindType
AeroModel = "DMS"#Inp.AeroModel
windINPfilename = "$(path)$(Inp.windINPfilename)"
ifw_libfile = Inp.ifw_libfile
if ifw_libfile == "nothing"
ifw_libfile = nothing
end
Blade_Height = Inp.Blade_Height
Blade_Radius = Inp.Blade_Radius
numTS = 5 # shortened since the floating simulation can take awhile
delta_t = Inp.delta_t
NuMad_geom_xlscsv_file_twr = "$(path)$(Inp.NuMad_geom_xlscsv_file_twr)"
NuMad_mat_xlscsv_file_twr = "$(path)$(Inp.NuMad_mat_xlscsv_file_twr)"
NuMad_geom_xlscsv_file_bld = "$(path)$(Inp.NuMad_geom_xlscsv_file_bld[1:end-4])DMS.csv"
NuMad_mat_xlscsv_file_bld = "$(path)$(Inp.NuMad_mat_xlscsv_file_bld)"
NuMad_geom_xlscsv_file_strut = "$(path)$(Inp.NuMad_geom_xlscsv_file_strut)"
NuMad_mat_xlscsv_file_strut = "$(path)$(Inp.NuMad_mat_xlscsv_file_strut)"
adi_lib = Inp.adi_lib
if adi_lib == "nothing"
adi_lib = nothing
end
adi_rootname = "$(path)$(Inp.adi_rootname)"
B = Nbld
R = Blade_Radius#177.2022*0.3048 #m
H = Blade_Height#1.02*R*2 #m
shapeZ = collect(LinRange(0,H,Nslices+1))
shapeX = R.*(1.0.-4.0.*(shapeZ/H.-.5).^2)#shapeX_spline(shapeZ)
nothing┌ Warning: The old combined yaml file is being depreciated in favor of the windio yaml and the modelingoptions format
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/topRunDLC.jl:647We will continue to use helper functions here to fully define the topside mesh, sectional properties, and their connection. However, note that our naming convention will be different in order to differentiate this mesh from the bottom side mesh. The outputs to the helper functions not being used elsewhere in the tutorial are returned as "_" to improve clarity.
Note that, unlike the previous examples, we will not apply a boundary condition to this mesh. The boundary condition previously was essentially fixing the bottom of the turbine tower in place, since it is installed in the ground. However, for our floating platform, we want the loads from the topside to propagate to the platform, so our restoring hydrodynamics and mooring loads will "constain" the combined meshes as we want.
topMesh, topEl, topOrt, topJoint, topSectionProps,
_, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
aeroForces, deformAero, _, _,topSystem,topAssembly, sections, _, _ = OWENS.setupOWENS(OWENSAero,path;
rho, Nslices, ntheta, RPM, Vinf, eta, B, H, R, shapeZ, shapeX,
ifw, WindType, delta_t, numTS, adi_lib, adi_rootname, windINPfilename, ifw_libfile,
NuMad_geom_xlscsv_file_twr,# = "$path/data/NuMAD_Geom_SNL_5MW_ARCUS_Cables.csv",
NuMad_mat_xlscsv_file_twr,# = "$path/data/NuMAD_Materials_SNL_5MW_D_TaperedTower.csv",
NuMad_geom_xlscsv_file_bld,# = "$path/data/NuMAD_Geom_SNL_5MW_ARCUS.csv",
NuMad_mat_xlscsv_file_bld,# = "$path/data/NuMAD_Materials_SNL_5MW_D_Carbon_LCDT_ThickFoils_ThinSkin.csv",
NuMad_geom_xlscsv_file_strut,
NuMad_mat_xlscsv_file_strut,
Htwr_base=towerHeight,
ntelem,
nbelem,
ncelem,
nselem,
joint_type = 0,
c_mount_ratio = 0.05,
AeroModel, #AD, DMS, AC
DynamicStallModel="BV",
RPI=true,
cables_connected_to_blade_base = true,
meshtype = turbineType)
nothingHub crossing at idx 1 and radially at 0.1 with AD15 hub radius of 0.1
Moving strut point from [-0.1,6.123233995736766e-18,30.5457273] to [-0.1,6.123233995736766e-18,30.5457273]
Hub crossing at idx 1 and radially at 0.1 with AD15 hub radius of 0.1
Moving strut point from [0.04999999999999999,0.08660254037844388,30.5457273] to [0.04999999999999999,0.08660254037844388,30.5457273]
Hub crossing at idx 1 and radially at 0.1 with AD15 hub radius of 0.1
Moving strut point from [0.05000000000000003,-0.08660254037844385,30.5457273] to [0.05000000000000003,-0.08660254037844385,30.5457273]
Hub crossing at idx 1 and radially at 0.1 with AD15 hub radius of 0.1
Moving strut point from [-0.1,6.123233995736767e-18,85.6371819] to [-0.1,6.123233995736767e-18,85.6371819]
Hub crossing at idx 1 and radially at 0.1 with AD15 hub radius of 0.1
Moving strut point from [0.04999999999999999,0.08660254037844388,85.6371819] to [0.04999999999999999,0.08660254037844388,85.6371819]
Hub crossing at idx 1 and radially at 0.1 with AD15 hub radius of 0.1
Moving strut point from [0.05000000000000004,-0.08660254037844385,85.6371819] to [0.05000000000000004,-0.08660254037844385,85.6371819]
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Data for SN curve control points not found in material file columns 23:28 for stress in Mpa, 29:33 for cycles in log10
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/fileio.jl:660
┌ Warning: Not using the componetized model is being depreciated, please consider updating to return_componentized=true
└ @ OWENS ~/work/OWENS.jl/OWENS.jl/src/SetupTurbine.jl:595Now we need to define the properties of the bottom side mesh, which we will do manually here instead of reading in values from a YAML file. This mesh is very simple: it consists of a single element with a nodal connection on either end. OWENS assumes the platform is a rigid body, so the element has very high stiffness to effectively function as a point mass.
nothingSince the mesh is so simple, most of these inputs are trivial and affect little. However, a few items are of note:
- The z positions in
bottomMeshstart at zero for the first node, and the second node is at the height of the tower base above the platform reference point (aka the platform location at the still water level above the platform centroid). In the case of the OC4 platform, this is 10 meters. This is also reflected in the Length argument in bottomOrt - The
bottomSectionPropsstructure has very large values for EIyy, EIzz, GJ, and EA, which are all the relevant stiffness properties, in order to effectively model the platform as a rigid body. - The
zcmparameter inbottomSectionPropsis set to the vertical distance of the platform with respect to the platform reference point. For the OC4 platform, this is 8.6588 meters below sea level.
bottomMesh = OWENSFEA.Mesh(
[1, 2], #nodeNum
1, #numElx
2, #numNodes
[0.0, 0.0], #x
[0.0, 0.0], #y
[0.0, 10.0], #z
[1], #elNum
[1 2], #conn
[0], #type
[0, 1], #meshSeg
zeros(1,1), #structuralSpanLocNorm
zeros(Int,1,1), #structuralNodeNumbers
zeros(Int,1,1) #structuralElNumbers
)
numBottomNodes = length(bottomMesh.z)
bottomOrt = OWENSFEA.Ort(
[180.0], #Psi_d
[-90.0], #Theta_d
[90.0], #Twist_d
[10.0], #Length
ones(1,1), #elNum
zeros(3,1)#Offset
)
n_ptfm_elem = 1
bottomSectionProps = Array{OWENSFEA.SectionPropsArray,1}(undef, n_ptfm_elem)
for ii = 1:n_ptfm_elem
bottomSectionProps[ii] = OWENSFEA.SectionPropsArray(
[0.0, 0.0], #ac
[0.0, 0.0], #twist_d
[0.0, 0.0], #rhoA
[1e18, 1e18], #EIyy
[1e18, 1e18], #EIzz
[1e18, 1e18], #GJ
[1e18, 1e18], #EA
[0.0, 0.0], #rhoIyy
[0.0, 0.0], #rhoIzz
[0.0, 0.0], #rhoJ
[-8.6588, 0.0], #zcm
[0.0, 0.0], #ycm,
[0.0, 0.0], #a
[0.0, 0.0], #EIyz
[0.0, 0.0], #alpha1
[0.0, 0.0], #alpha2
[0.0, 0.0], #alpha3
[0.0, 0.0], #alpha4
[0.0, 0.0], #alpha5
[0.0, 0.0], #alpha6
[0.0, 0.0], #rhoIyz
[0.0, 0.0], #b
[0.0, 0.0], #a0
[0.0, 0.0], #aeroCenterOffset
)
end
nothingWe also want to provide some additional general inputs to enable floating simulation. platformActive being the principal and obvious one, but also the input files needed by HydroDyn and MoorDyn, which are the external libraries that calculate the hydrodynamic and mooring loads to send to the bottom side mesh. These input files include:
- interpOrder: the degree of interpolation used to predict the future states within HydroDyn and MoorDyn, used in the
HD_UpdateStatesandMD_UpdateStatesfunctions - hdinputfile: the base HydroDyn input file with a
.datextension - mdinputfile: the base MoorDyn input file with a
.datextension - ssinputfile: the sea state input file with a
.datextension used to define the environmental conditions of the sea. This is used within HydroDyn. - potfilefile: the prefix of the potential flow files. At minimum, the directory must contain the .1, .3, and .hst WAMIT output files. This is also used within HydroDyn.
See the OpenFAST documentation for HydroDyn (https://openfast.readthedocs.io/en/main/source/user/hydrodyn/input_files.html) and MoorDyn (https://moordyn.readthedocs.io/en/latest/inputs.html) for more information about how to format these input files. For simplicity here, we will use predefined input files for the OC4 semisubmersible platform, which comes with OpenFAST and is copied to the data folder here.
if AeroModel=="AD"
AD15On = true
else
AD15On = false
end
platformActive = true
interpOrder = 2
hd_input_file = "$(path)/data/HydroDyn.dat"
md_input_file = "$(path)/data/MoorDyn.dat"
ss_input_file = "$(path)/data/SeaState.dat"
potflowfile = "$(path)/data/potflowdata/marin_semi"
inputs = OWENS.Inputs(;analysisType = analysisType,
tocp = [0.0,100000.1],
Omegaocp = [RPM,RPM] ./ 60,
tocp_Vinf = [0.0,100000.1],
Vinfocp = [Vinf,Vinf],
numTS,
delta_t,
AD15On,
aeroLoadsOn = 2,
platformActive,
interpOrder,
hd_input_file,
md_input_file,
ss_input_file,
potflowfile)
nothingWe also will need to specify where the HydroDyn and MoorDyn libraries exist on the computer. We will set these to nothing so the auto-installed libraries will be used by default. Change these variables as needed if the location of these libraries are different.
hd_lib = nothing
md_lib = nothing
bin = OWENS.Bin(hd_lib, md_lib)
nothingAs in previous examples, we need to define the inputs for the finite element models, though now we need a definition for both the topside and the bottom side.
The changes to the topside from example A are the lack of boundary conditions, non-default gamma and alpha terms (for better convergence with the platform mesh).
topFEAModel = OWENS.FEAModel(;
analysisType = analysisType,
dataOutputFilename = "none",
joint = topJoint,
platformTurbineConnectionNodeNumber = 1,
nlOn = true,
numNodes = topMesh.numNodes,
RayleighAlpha = 0.05,
RayleighBeta = 0.05,
iterationType = "DI",
gamma = 1.0,
alpha = 0.5)
nothingThe bottom finite element mesh will use the same inputs with two exceptions: joints is now empty (since we're considering the platform as a single rigid body), and we add an additional concentrated term at the bottom node to represent the rigid body mass matrix of the platform. The values here are pulled from the report "Definition of the Semisubmersible Floating System for Phase II of OC4" by Robertson et al. (2014) (https://doi.org/10.2172/1155123).
numDOFPerNode = 6
ptfm_mass = 3.85218e6
ptfm_roll_iner = 2.56193e9
ptfm_pitch_iner = 2.56193e9
ptfm_yaw_iner = 4.24265e9
ptfm_Ixx = ptfm_roll_iner
ptfm_Iyy = ptfm_pitch_iner
ptfm_Izz = ptfm_yaw_iner
ptfm_Ixy = 0.0
ptfm_Iyz = 0.0
ptfm_Ixz = 0.0
ptfm_mass = [
ptfm_mass 0.0 0.0 0.0 0.0 0.0
0.0 ptfm_mass 0.0 0.0 0.0 0.0
0.0 0.0 ptfm_mass 0.0 0.0 0.0
0.0 0.0 0.0 ptfm_Ixx ptfm_Ixy ptfm_Ixz
0.0 0.0 0.0 ptfm_Ixy ptfm_Iyy -ptfm_Iyz
0.0 0.0 0.0 ptfm_Ixz -ptfm_Iyz ptfm_Izz
]
bottomConcInputs = [
1 "M6" 1 1 ptfm_mass[1, 1]
1 "M6" 1 2 ptfm_mass[1, 2]
1 "M6" 1 3 ptfm_mass[1, 3]
1 "M6" 1 4 ptfm_mass[1, 4]
1 "M6" 1 5 ptfm_mass[1, 5]
1 "M6" 1 6 ptfm_mass[1, 6]
1 "M6" 2 1 ptfm_mass[2, 1]
1 "M6" 2 2 ptfm_mass[2, 2]
1 "M6" 2 3 ptfm_mass[2, 3]
1 "M6" 2 4 ptfm_mass[2, 4]
1 "M6" 2 5 ptfm_mass[2, 5]
1 "M6" 2 6 ptfm_mass[2, 6]
1 "M6" 3 1 ptfm_mass[3, 1]
1 "M6" 3 2 ptfm_mass[3, 2]
1 "M6" 3 3 ptfm_mass[3, 3]
1 "M6" 3 4 ptfm_mass[3, 4]
1 "M6" 3 5 ptfm_mass[3, 5]
1 "M6" 3 6 ptfm_mass[3, 6]
1 "M6" 4 1 ptfm_mass[4, 1]
1 "M6" 4 2 ptfm_mass[4, 2]
1 "M6" 4 3 ptfm_mass[4, 3]
1 "M6" 4 4 ptfm_mass[4, 4]
1 "M6" 4 5 ptfm_mass[4, 5]
1 "M6" 4 6 ptfm_mass[4, 6]
1 "M6" 5 1 ptfm_mass[5, 1]
1 "M6" 5 2 ptfm_mass[5, 2]
1 "M6" 5 3 ptfm_mass[5, 3]
1 "M6" 5 4 ptfm_mass[5, 4]
1 "M6" 5 5 ptfm_mass[5, 5]
1 "M6" 5 6 ptfm_mass[5, 6]
1 "M6" 6 1 ptfm_mass[6, 1]
1 "M6" 6 2 ptfm_mass[6, 2]
1 "M6" 6 3 ptfm_mass[6, 3]
1 "M6" 6 4 ptfm_mass[6, 4]
1 "M6" 6 5 ptfm_mass[6, 5]
1 "M6" 6 6 ptfm_mass[6, 6]
]
bottomConcTerms = OWENSFEA.applyConcentratedTerms(
bottomMesh.numNodes,
numDOFPerNode,
data=bottomConcInputs,
jointData=[])
bottomFEAModel = OWENS.FEAModel(;
analysisType = analysisType,
dataOutputFilename = "none",
joint = [],
platformTurbineConnectionNodeNumber = 1,
nlOn = true,
numNodes = bottomMesh.numNodes,
RayleighAlpha = 0.05,
RayleighBeta = 0.05,
iterationType = "DI",
gamma = 1.0,
alpha = 0.5,
nodalTerms=bottomConcTerms)
nothingFinally, some additional structure definitions for the bottom mesh. These don't contain anything new, but are required for OWENS to run properly.
bottomEl = OWENSFEA.El(bottomSectionProps,
bottomOrt.Length,
bottomOrt.Psi_d,
bottomOrt.Theta_d,
bottomOrt.Twist_d,
ones(bottomMesh.numEl))
nothingWe can now run our floating simulation in OWENS. Note that we are using OWENS.Unsteady now instead of OWENS.Unsteady_Land.
println("Running Unsteady")
t, aziHist, OmegaHist, OmegaDotHist, gbHist, gbDotHist, gbDotDotHist, FReactionHist, FTwrBsHist,
genTorque, genPower, torqueDriveShaft, uHist, uHist_prp,
epsilon_x_hist, epsilon_y_hist, epsilon_z_hist, kappa_x_hist, kappa_y_hist, kappa_z_hist,
FPtfmHist, FHydroHist, FMooringHist = OWENS.Unsteady(inputs,
system=topSystem,
assembly=topAssembly,
topModel=topFEAModel,
topMesh=topMesh,
topEl=topEl,
aero=aeroForces,
deformAero=deformAero,
bottomModel=bottomFEAModel,
bottomMesh=bottomMesh,
bottomEl=bottomEl,
bin=bin)
if AD15On #TODO: move this into the run functions
OWENS.OWENSOpenFASTWrappers.endTurb()
end
nothingRunning Unsteady
Reading HydroDyn data from /home/runner/work/OWENS.jl/OWENS.jl/examples/literate/data/HydroDyn.dat.
Reading SeaState data from /home/runner/work/OWENS.jl/OWENS.jl/examples/literate/data/SeaState.dat.
**************************************************************************************************
HydroDyn+SeaState library
Copyright (C) 2024 National Renewable Energy Laboratory
Copyright (C) 2024 Envision Energy USA LTD
This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
See the "LICENSE" file distributed with this software for details.
**************************************************************************************************
HydroDyn+SeaState library-v2.5.0-3490-g24c05a74
Compile Info:
- Compiler: GCC version 8.1.0
- Architecture: 64 bit
- Precision: double
- OpenMP: No
- Date: Nov 5 2024
- Time: 02:19:23
Execution Info:
- Date: 05/21/2025
- Time: 13:49:09+0000
Running SeaState.
Generating incident wave kinematics and current time history.
Running HydroDyn.
Reading in WAMIT output with root name
"/home/runner/work/OWENS.jl/OWENS.jl/examples/literate/data/potflowdata/marin_semi".
Computing radiation impulse response functions and wave diffraction forces.
Reading MoorDyn data from /home/runner/work/OWENS.jl/OWENS.jl/examples/literate/data/MoorDyn.dat.
**************************************************************************************************
MoorDyn library
Copyright (C) 2024 National Renewable Energy Laboratory
Copyright (C) 2024 Envision Energy USA LTD
This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
See the "LICENSE" file distributed with this software for details.
**************************************************************************************************
MoorDyn library-v2.5.0-3490-g24c05a74
Compile Info:
- Compiler: GCC version 8.1.0
- Architecture: 64 bit
- Precision: double
- OpenMP: No
- Date: Nov 5 2024
- Time: 02:19:23
Execution Info:
- Date: 05/21/2025
- Time: 13:49:18+0000
Running MoorDyn (v2.0.0, 2022-12-08).
This is MoorDyn v2, with significant input file changes from v1.
Copyright: (C) 2022 National Renewable Energy Laboratory, (C) 2019 Matt Hall
Parsing MoorDyn input file: notUsed
Created mooring system: 3 lines, 6 points, 0 rods, 0 bodies.
Finalizing initial conditions using dynamic relaxation.
t=1 FairTen 1: 1.18264E+06, 1, 2
t=2 FairTen 1: 1.24078E+06, 1.18264E+06, 1
t=3 FairTen 1: 1.01721E+06, 1.24078E+06, 1.18264E+06
t=4 FairTen 1: 1.11839E+06, 1.01721E+06, 1.24078E+06
t=5 FairTen 1: 1.13732E+06, 1.11839E+06, 1.01721E+06
t=6 FairTen 1: 1.05990E+06, 1.13732E+06, 1.11839E+06
t=7 FairTen 1: 1.10862E+06, 1.05990E+06, 1.13732E+06
t=8 FairTen 1: 1.13130E+06, 1.10862E+06, 1.05990E+06
t=9 FairTen 1: 1.09750E+06, 1.13130E+06, 1.10862E+06
t=10 FairTen 1: 1.07609E+06, 1.09750E+06, 1.13130E+06
t=11 FairTen 1: 1.12607E+06, 1.07609E+06, 1.09750E+06
t=12 FairTen 1: 1.10720E+06, 1.12607E+06, 1.07609E+06
t=13 FairTen 1: 1.08413E+06, 1.10720E+06, 1.12607E+06
t=14 FairTen 1: 1.12831E+06, 1.08413E+06, 1.10720E+06
t=15 FairTen 1: 1.09024E+06, 1.12831E+06, 1.08413E+06
t=16 FairTen 1: 1.08803E+06, 1.09024E+06, 1.12831E+06
t=17 FairTen 1: 1.13554E+06, 1.08803E+06, 1.09024E+06
t=18 FairTen 1: 1.08932E+06, 1.13554E+06, 1.08803E+06
t=19 FairTen 1: 1.09424E+06, 1.08932E+06, 1.13554E+06
t=20 FairTen 1: 1.11980E+06, 1.09424E+06, 1.08932E+06
t=21 FairTen 1: 1.09232E+06, 1.11980E+06, 1.09424E+06
t=22 FairTen 1: 1.10695E+06, 1.09232E+06, 1.11980E+06
t=23 FairTen 1: 1.11280E+06, 1.10695E+06, 1.09232E+06
t=24 FairTen 1: 1.09306E+06, 1.11280E+06, 1.10695E+06
t=25 FairTen 1: 1.10221E+06, 1.09306E+06, 1.11280E+06
t=26 FairTen 1: 1.11079E+06, 1.10221E+06, 1.09306E+06
t=27 FairTen 1: 1.10566E+06, 1.11079E+06, 1.10221E+06
t=28 FairTen 1: 1.09657E+06, 1.10566E+06, 1.11079E+06
t=29 FairTen 1: 1.10655E+06, 1.09657E+06, 1.10566E+06
t=30 FairTen 1: 1.10621E+06, 1.10655E+06, 1.09657E+06
t=31 FairTen 1: 1.09677E+06, 1.10621E+06, 1.10655E+06
t=32 FairTen 1: 1.11356E+06, 1.09677E+06, 1.10621E+06
t=33 FairTen 1: 1.10210E+06, 1.11356E+06, 1.09677E+06
t=34 FairTen 1: 1.09210E+06, 1.10210E+06, 1.11356E+06
t=35 FairTen 1: 1.11583E+06, 1.09210E+06, 1.10210E+06
t=36 FairTen 1: 1.10303E+06, 1.11583E+06, 1.09210E+06
t=37 FairTen 1: 1.09724E+06, 1.10303E+06, 1.11583E+06
t=38 FairTen 1: 1.11096E+06, 1.09724E+06, 1.10303E+06
t=39 FairTen 1: 1.09837E+06, 1.11096E+06, 1.09724E+06
t=40 FairTen 1: 1.10299E+06, 1.09837E+06, 1.11096E+06
t=41 FairTen 1: 1.11057E+06, 1.10299E+06, 1.09837E+06
t=42 FairTen 1: 1.10042E+06, 1.11057E+06, 1.10299E+06
t=43 FairTen 1: 1.10069E+06, 1.10042E+06, 1.11057E+06
t=44 FairTen 1: 1.10549E+06, 1.10069E+06, 1.10042E+06
t=45 FairTen 1: 1.10619E+06, 1.10549E+06, 1.10069E+06
t=46 FairTen 1: 1.10234E+06, 1.10619E+06, 1.10549E+06
t=47 FairTen 1: 1.10330E+06, 1.10234E+06, 1.10619E+06
t=48 FairTen 1: 1.10515E+06, 1.10330E+06, 1.10234E+06 Fairlead tensions converged to 1% after 48 seconds.
MoorDyn initialization completed.
Running in specified rotor speed mode
25.0%┣██████████▊ ┫ 1/4 [00:00<Inf:Inf, InfGs/it]
50.0%┣████████████████████████ ┫ 2/4 [00:04<00:08, 4s/it]
75.0%┣████████████████████████████████████ ┫ 3/4 [00:06<00:03, 3s/it]
100.0%┣███████████████████████████████████████████████┫ 4/4 [00:09<00:00, 3s/it]
Simulation Complete.Since we ran the simulation for such a short time period for sake of demonstrating completion of the simulation, the plots below are not too interesting. Try experimenting with changing the numTS variable above to a longer simulation to see how the floating VAWT responds over time. Note that the platform modeled here is not necessarily design for the Darrieus VAWT it is supporting, so your mileage may vary on the simulation results. PyPlot.plot(t, uHist_prp[:, 1]) PyPlot.title("Platform Surge Response") PyPlot.xlabel("Time [s]") PyPlot.ylabel("Force [N]") PyPlot.plot(t, FHydroHist[:, 1]) PyPlot.xlabel("Time [s]") PyPlot.ylabel("Force [N]") PyPlot.title("Hydrodynamic Surge Loading")
nothingThis page was generated using Literate.jl.