Skip to content
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

Prepostupdate #96

Open
wants to merge 5 commits into
base: master
Choose a base branch
from
Open

Prepostupdate #96

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
f250dcf
Minor updates to pre and post processing to include more information …
kevmoor Jan 8, 2025
50f3207
Add DLCParams as an internal data part of DLC_options for internal use
kevmoor Jan 8, 2025
e681bb4
Add helper function for the single array unsteady output option
kevmoor Jan 9, 2025
630ad26
Fix fatigue calculation error
kevmoor Jan 24, 2025
1c36d0e
Also missing delta_t input for fatigue rate calculation, will need to…
kevmoor Jan 24, 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
3 changes: 2 additions & 1 deletion .gitignore
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,4 +24,5 @@ new*
*.mp4
*turbsim*/
*animation/
*.bts
*.bts
*.jld*
9 changes: 5 additions & 4 deletions examples/Optimization/windio_example.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,7 +173,8 @@ println("Percent Stress Pass: $(ipass/length(stress_U_UNIT)*100)%")
@test isapprox(topstrainout_blade_L_UNIT,topstrainout_blade_L;atol=maximum(abs.(topstrainout_blade_L_UNIT))*frac)
@test isapprox(topstrainout_tower_U_UNIT,topstrainout_tower_U;atol=maximum(abs.(topstrainout_tower_U_UNIT))*frac)
@test isapprox(topstrainout_tower_L_UNIT,topstrainout_tower_L;atol=maximum(abs.(topstrainout_tower_L_UNIT))*frac)
@test isapprox(topDamage_blade_U_UNIT,topDamage_blade_U;atol=maximum(abs.(topDamage_blade_U_UNIT))*frac)
@test isapprox(topDamage_blade_L_UNIT,topDamage_blade_L;atol=maximum(abs.(topDamage_blade_L_UNIT))*frac)
@test isapprox(topDamage_tower_U_UNIT,topDamage_tower_U;atol=maximum(abs.(topDamage_tower_U_UNIT))*frac)
@test isapprox(topDamage_tower_L_UNIT,topDamage_tower_L;atol=maximum(abs.(topDamage_tower_L_UNIT))*frac)
# TODO: write damage verification
# @test isapprox(topDamage_blade_U_UNIT,topDamage_blade_U;atol=maximum(abs.(topDamage_blade_U_UNIT))*frac)
# @test isapprox(topDamage_blade_L_UNIT,topDamage_blade_L;atol=maximum(abs.(topDamage_blade_L_UNIT))*frac)
# @test isapprox(topDamage_tower_U_UNIT,topDamage_tower_U;atol=maximum(abs.(topDamage_tower_U_UNIT))*frac)
# @test isapprox(topDamage_tower_L_UNIT,topDamage_tower_L;atol=maximum(abs.(topDamage_tower_L_UNIT))*frac)
146 changes: 1 addition & 145 deletions examples/RM2/RM2_Medium_IO.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,148 +10,4 @@ OWENS_Options = "$runpath/modeling_options_OWENS_RM2.yml"

WINDIO_filename = "$runpath/WINDIO_RM2.yaml"

OWENS.runOWENSWINDIO(OWENS_Options,WINDIO_filename,runpath)


file = "./InitialDataOutputs_windio.h5"
t_windio = HDF5.h5read(file,"t")
aziHist_windio = HDF5.h5read(file,"aziHist")
OmegaHist_windio = HDF5.h5read(file,"OmegaHist")
OmegaDotHist_windio = HDF5.h5read(file,"OmegaDotHist")
gbHist_windio = HDF5.h5read(file,"gbHist")
gbDotHist_windio = HDF5.h5read(file,"gbDotHist")
gbDotDotHist_windio = HDF5.h5read(file,"gbDotDotHist")
FReactionHist_windio = HDF5.h5read(file,"FReactionHist")
FTwrBsHist_windio = HDF5.h5read(file,"FTwrBsHist")
genTorque_windio = HDF5.h5read(file,"genTorque")
genPower_windio = HDF5.h5read(file,"genPower")
torqueDriveShaft_windio = HDF5.h5read(file,"torqueDriveShaft")
uHist_windio = HDF5.h5read(file,"uHist")
uHist_prp_windio = HDF5.h5read(file,"uHist_prp")
epsilon_x_hist_windio = HDF5.h5read(file,"epsilon_x_hist")
epsilon_y_hist_windio = HDF5.h5read(file,"epsilon_y_hist")
epsilon_z_hist_windio = HDF5.h5read(file,"epsilon_z_hist")
kappa_x_hist_windio = HDF5.h5read(file,"kappa_x_hist")
kappa_y_hist_windio = HDF5.h5read(file,"kappa_y_hist")
kappa_z_hist_windio = HDF5.h5read(file,"kappa_z_hist")
massOwens_windio = HDF5.h5read(file,"massOwens")
stress_U_windio = HDF5.h5read(file,"stress_U")
SF_ult_U_windio = HDF5.h5read(file,"SF_ult_U")
SF_buck_U_windio = HDF5.h5read(file,"SF_buck_U")
stress_L_windio = HDF5.h5read(file,"stress_L")
SF_ult_L_windio = HDF5.h5read(file,"SF_ult_L")
SF_buck_L_windio = HDF5.h5read(file,"SF_buck_L")
stress_TU_windio = HDF5.h5read(file,"stress_TU")
SF_ult_TU_windio = HDF5.h5read(file,"SF_ult_TU")
SF_buck_TU_windio = HDF5.h5read(file,"SF_buck_TU")
stress_TL_windio = HDF5.h5read(file,"stress_TL")
SF_ult_TL_windio = HDF5.h5read(file,"SF_ult_TL")
SF_buck_TL_windio = HDF5.h5read(file,"SF_buck_TL")
topstrainout_blade_U_windio = HDF5.h5read(file,"topstrainout_blade_U")
topstrainout_blade_L_windio = HDF5.h5read(file,"topstrainout_blade_L")
topstrainout_tower_U_windio = HDF5.h5read(file,"topstrainout_tower_U")
topstrainout_tower_L_windio = HDF5.h5read(file,"topstrainout_tower_L")
topDamage_blade_U_windio = HDF5.h5read(file,"topDamage_blade_U")
topDamage_blade_L_windio = HDF5.h5read(file,"topDamage_blade_L")
topDamage_tower_U_windio = HDF5.h5read(file,"topDamage_tower_U")
topDamage_tower_L_windio = HDF5.h5read(file,"topDamage_tower_L")


file = "./InitialDataOutputs_scripting.h5"
t = HDF5.h5read(file,"t")
aziHist = HDF5.h5read(file,"aziHist")
OmegaHist = HDF5.h5read(file,"OmegaHist")
OmegaDotHist = HDF5.h5read(file,"OmegaDotHist")
gbHist = HDF5.h5read(file,"gbHist")
gbDotHist = HDF5.h5read(file,"gbDotHist")
gbDotDotHist = HDF5.h5read(file,"gbDotDotHist")
FReactionHist = HDF5.h5read(file,"FReactionHist")
FTwrBsHist = HDF5.h5read(file,"FTwrBsHist")
genTorque = HDF5.h5read(file,"genTorque")
genPower = HDF5.h5read(file,"genPower")
torqueDriveShaft = HDF5.h5read(file,"torqueDriveShaft")
uHist = HDF5.h5read(file,"uHist")
uHist_prp = HDF5.h5read(file,"uHist_prp")
epsilon_x_hist = HDF5.h5read(file,"epsilon_x_hist")
epsilon_y_hist = HDF5.h5read(file,"epsilon_y_hist")
epsilon_z_hist = HDF5.h5read(file,"epsilon_z_hist")
kappa_x_hist = HDF5.h5read(file,"kappa_x_hist")
kappa_y_hist = HDF5.h5read(file,"kappa_y_hist")
kappa_z_hist = HDF5.h5read(file,"kappa_z_hist")
massOwens = HDF5.h5read(file,"massOwens")
stress_U = HDF5.h5read(file,"stress_U")
SF_ult_U = HDF5.h5read(file,"SF_ult_U")
SF_buck_U = HDF5.h5read(file,"SF_buck_U")
stress_L = HDF5.h5read(file,"stress_L")
SF_ult_L = HDF5.h5read(file,"SF_ult_L")
SF_buck_L = HDF5.h5read(file,"SF_buck_L")
stress_TU = HDF5.h5read(file,"stress_TU")
SF_ult_TU = HDF5.h5read(file,"SF_ult_TU")
SF_buck_TU = HDF5.h5read(file,"SF_buck_TU")
stress_TL = HDF5.h5read(file,"stress_TL")
SF_ult_TL = HDF5.h5read(file,"SF_ult_TL")
SF_buck_TL = HDF5.h5read(file,"SF_buck_TL")
topstrainout_blade_U = HDF5.h5read(file,"topstrainout_blade_U")
topstrainout_blade_L = HDF5.h5read(file,"topstrainout_blade_L")
topstrainout_tower_U = HDF5.h5read(file,"topstrainout_tower_U")
topstrainout_tower_L = HDF5.h5read(file,"topstrainout_tower_L")
topDamage_blade_U = HDF5.h5read(file,"topDamage_blade_U")
topDamage_blade_L = HDF5.h5read(file,"topDamage_blade_L")
topDamage_tower_U = HDF5.h5read(file,"topDamage_tower_U")
topDamage_tower_L = HDF5.h5read(file,"topDamage_tower_L")

import PyPlot
PyPlot.pygui(true)
PyPlot.figure()
PyPlot.plot(t_windio,FReactionHist_windio[:,6])
PyPlot.plot(t,FReactionHist[:,6])

frac = 1e-5
@test isapprox(t_windio,t;atol=maximum(abs.(t_windio))*frac)
@test isapprox(aziHist_windio,aziHist;atol=maximum(abs.(aziHist_windio))*frac)
@test isapprox(OmegaHist_windio,OmegaHist;atol=maximum(abs.(OmegaHist_windio))*frac)
@test isapprox(OmegaDotHist_windio,OmegaDotHist;atol=maximum(abs.(OmegaDotHist_windio))*frac)
@test isapprox(gbHist_windio,gbHist;atol=maximum(abs.(gbHist_windio))*frac)
@test isapprox(gbDotHist_windio,gbDotHist;atol=maximum(abs.(gbDotHist_windio))*frac)
@test isapprox(gbDotDotHist_windio,gbDotDotHist;atol=maximum(abs.(gbDotDotHist_windio))*frac)
@test isapprox(FReactionHist_windio,FReactionHist;atol=maximum(abs.(FReactionHist_windio))*frac)
@test isapprox(FTwrBsHist_windio,FTwrBsHist;atol=maximum(abs.(FTwrBsHist_windio))*frac)
@test isapprox(genTorque_windio,genTorque;atol=maximum(abs.(genTorque_windio))*frac)
@test isapprox(genPower_windio,genPower;atol=maximum(abs.(genPower_windio))*frac)
@test isapprox(torqueDriveShaft_windio,torqueDriveShaft;atol=maximum(abs.(torqueDriveShaft_windio))*frac)
@test isapprox(uHist_windio,uHist;atol=maximum(abs.(uHist_windio))*frac)
@test isapprox(uHist_prp_windio,uHist_prp;atol=maximum(abs.(uHist_prp_windio))*frac)
@test isapprox(epsilon_x_hist_windio,epsilon_x_hist;atol=maximum(abs.(epsilon_x_hist_windio))*frac)
@test isapprox(epsilon_y_hist_windio,epsilon_y_hist;atol=maximum(abs.(epsilon_y_hist_windio))*frac)
@test isapprox(epsilon_z_hist_windio,epsilon_z_hist;atol=maximum(abs.(epsilon_z_hist_windio))*frac)
@test isapprox(kappa_x_hist_windio,kappa_x_hist;atol=maximum(abs.(kappa_x_hist_windio))*frac)
@test isapprox(kappa_y_hist_windio,kappa_y_hist;atol=maximum(abs.(kappa_y_hist_windio))*frac)
@test isapprox(kappa_z_hist_windio,kappa_z_hist;atol=maximum(abs.(kappa_z_hist_windio))*frac)
@test isapprox(massOwens_windio,massOwens;atol=maximum(abs.(massOwens_windio))*frac)
ipass = 0
for i = 1:length(stress_U_windio)
# println("$i of $(length(stress_U_windio))")
if isapprox(stress_U_windio[i],stress_U[i];atol=maximum(abs.(stress_U_windio[i]))*frac)
global ipass += 1
end
end
println("Percent Stress Pass: $(ipass/length(stress_U_windio)*100)%")
@test isapprox(SF_ult_U_windio,SF_ult_U;atol=maximum(abs.(SF_ult_U_windio))*frac)
@test isapprox(SF_buck_U_windio,SF_buck_U;atol=maximum(abs.(SF_buck_U_windio))*frac)
@test isapprox(stress_L_windio,stress_L;atol=maximum(abs.(stress_L_windio))*frac)
@test isapprox(SF_ult_L_windio,SF_ult_L;atol=maximum(abs.(SF_ult_L_windio))*frac)
@test isapprox(SF_buck_L_windio,SF_buck_L;atol=maximum(abs.(SF_buck_L_windio))*frac)
@test isapprox(stress_TU_windio,stress_TU;atol=maximum(abs.(stress_TU_windio))*frac)
@test isapprox(SF_ult_TU_windio,SF_ult_TU;atol=maximum(abs.(SF_ult_TU_windio))*frac)
@test isapprox(SF_buck_TU_windio,SF_buck_TU;atol=maximum(abs.(SF_buck_TU_windio))*frac)
@test isapprox(stress_TL_windio,stress_TL;atol=maximum(abs.(stress_TL_windio))*frac)
@test isapprox(SF_ult_TL_windio,SF_ult_TL;atol=maximum(abs.(SF_ult_TL_windio))*frac)
@test isapprox(SF_buck_TL_windio,SF_buck_TL;atol=maximum(abs.(SF_buck_TL_windio))*frac)
@test isapprox(topstrainout_blade_U_windio,topstrainout_blade_U;atol=maximum(abs.(topstrainout_blade_U_windio))*frac)
@test isapprox(topstrainout_blade_L_windio,topstrainout_blade_L;atol=maximum(abs.(topstrainout_blade_L_windio))*frac)
@test isapprox(topstrainout_tower_U_windio,topstrainout_tower_U;atol=maximum(abs.(topstrainout_tower_U_windio))*frac)
@test isapprox(topstrainout_tower_L_windio,topstrainout_tower_L;atol=maximum(abs.(topstrainout_tower_L_windio))*frac)
@test isapprox(topDamage_blade_U_windio,topDamage_blade_U;atol=maximum(abs.(topDamage_blade_U_windio))*frac)
@test isapprox(topDamage_blade_L_windio,topDamage_blade_L;atol=maximum(abs.(topDamage_blade_L_windio))*frac)
@test isapprox(topDamage_tower_U_windio,topDamage_tower_U;atol=maximum(abs.(topDamage_tower_U_windio))*frac)
@test isapprox(topDamage_tower_L_windio,topDamage_tower_L;atol=maximum(abs.(topDamage_tower_L_windio))*frac)
OWENS.runOWENSWINDIO(OWENS_Options,WINDIO_filename,runpath)
34 changes: 34 additions & 0 deletions examples/RM2/airfoils/NACA_0021.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,34 @@
100,0.0
95,-1.412
90,-2.534
80,-4.591
70,-6.412
60,-7.986
50,-9.265
40,-10.156
30,-10.504
25,-10.397
20,-10.04
15,-9.354
10,-8.195
7.5,-7.35
5,-6.221
2.5,-4.576
1.25,-3.315
0,0
1.25,3.315
2.5,4.576
5,6.221
7.5,7.35
10,8.195
15,9.354
20,10.04
25,10.397
30,10.504
40,10.156
50,9.265
60,7.986
70,6.412
80,4.591
90,2.534
95,1.412
Loading
Loading