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

Jump to bottom

Update ROM to only run rom equations if the rom model is active #19

Merged
merged 1 commit into from
Aug 9, 2024
Merged

Update ROM to only run rom equations if the rom model is active #19

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
156 changes: 79 additions & 77 deletions src/rom.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,83 +16,85 @@ This function executes a reduced order model analysis.
* `rom` object containing a reduced order feamodel
"""
function reducedOrderModel(elStorage,feamodel,mesh,el,displ)

countedNodes = []

rom0 = calculateROM(feamodel,mesh,el,displ,zeros(9),elStorage,countedNodes) #Omega = 0 #calculates system matrices for parked condition

#calculates system matrices for various acceleration, rotor speed, rotor accelration combinations
omx = 1
omy = 1
omz = 1

omxdot = 1
omydot = 1
omzdot = 1

a_x = 1
a_y = 1
a_z = 1

rom1 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 omx 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_x = 1 #OmegaDot_i = 0 accel _i = 0
rom2 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 omy 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_y = 1 #OmegaDot_i = 0 accel_i = 0
rom3 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 omz 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_z = 1 #OmegaDot_i = 0 accel_i = 0
rom4 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 omx omy 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_x,y = 1 #OmegaDot_i = 0 accel_i = 0
rom5 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 omy omz 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_y,z = 1 #OmegaDot_i = 0 accel_i = 0
rom6 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 omx 0.0 omz 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_x,z = 1 #OmegaDot_i = 0 accel_i = 0
rom7 = calculateROMGyric(feamodel,mesh,el,displ,[a_x 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_i = 0 accel_x = 1
rom8 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 a_y 0.0 0.0 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_i = 0 accel_y = 1
rom9 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 a_z 0.0 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_i = 0 accel_z = 1
rom10 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 0.0 omxdot 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_1 = 0 accel_i=0
rom11 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 0.0 0.0 omydot 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_2 = 0 accel_i=0
rom12 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 omzdot],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_3 = 0 accel_i=0

#reduced order structural stiffness, mass, and damping
Kr = rom0.Kr
Mr = rom0.Mr
Cr = rom0.Cr + rom0.CrModal
Fr = rom0.Fr

#reduced order transformation matrices
Phi = rom0.Phi
invPhi = rom0.invPhi

#reduced order spin softening coefficient matrices
SrOx2 = rom1.Kr
SrOy2 = rom2.Kr
SrOz2 = rom3.Kr
SrOxOy = rom4.Kr - rom1.Kr - rom2.Kr
SrOyOz = rom5.Kr - rom2.Kr - rom3.Kr
SrOxOz = rom6.Kr - rom1.Kr - rom3.Kr

#reduced order centrifugal load coefficient matrices
FrOx2 = rom1.Fr
FrOy2 = rom2.Fr
FrOz2 = rom3.Fr
FrOxOy = rom4.Fr - rom1.Fr - rom2.Fr
FrOyOz = rom5.Fr - rom2.Fr - rom3.Fr
FrOxOz = rom6.Fr - rom1.Fr - rom3.Fr
FrOxdot = rom10.Fr
FrOydot = rom11.Fr
FrOzdot = rom12.Fr

#need to calculate reduced order acceleration/body force coefeficients
FrAx = rom7.Fr
FrAy = rom8.Fr
FrAz = rom9.Fr

#reduced order gyric coefficient matrices
GrOx = rom1.Cr
GrOy = rom2.Cr
GrOz = rom3.Cr

#reduced order circulatory coefficient matrices
HrOx = 0.5*rom1.Cr
HrOy = 0.5*rom2.Cr
HrOz = 0.5*rom3.Cr

rom = ROM(Kr,Mr,Cr,0.0,Fr,Phi,invPhi,SrOx2,SrOy2,SrOz2,SrOxOy,SrOyOz,SrOxOz,FrOx2,FrOy2,FrOz2,FrOxOy,FrOyOz,FrOxOz,FrOxdot,FrOydot,FrOzdot,FrAx,FrAy,FrAz,GrOx,GrOy,GrOz,HrOx,HrOy,HrOz)

if feamodel.analysisType == "ROM"
countedNodes = []

rom0 = calculateROM(feamodel,mesh,el,displ,zeros(9),elStorage,countedNodes) #Omega = 0 #calculates system matrices for parked condition

#calculates system matrices for various acceleration, rotor speed, rotor accelration combinations
omx = 1
omy = 1
omz = 1

omxdot = 1
omydot = 1
omzdot = 1

a_x = 1
a_y = 1
a_z = 1

rom1 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 omx 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_x = 1 #OmegaDot_i = 0 accel _i = 0
rom2 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 omy 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_y = 1 #OmegaDot_i = 0 accel_i = 0
rom3 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 omz 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_z = 1 #OmegaDot_i = 0 accel_i = 0
rom4 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 omx omy 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_x,y = 1 #OmegaDot_i = 0 accel_i = 0
rom5 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 omy omz 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_y,z = 1 #OmegaDot_i = 0 accel_i = 0
rom6 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 omx 0.0 omz 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_x,z = 1 #OmegaDot_i = 0 accel_i = 0
rom7 = calculateROMGyric(feamodel,mesh,el,displ,[a_x 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_i = 0 accel_x = 1
rom8 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 a_y 0.0 0.0 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_i = 0 accel_y = 1
rom9 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 a_z 0.0 0.0 0.0 0.0 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_i = 0 accel_z = 1
rom10 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 0.0 omxdot 0.0 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_1 = 0 accel_i=0
rom11 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 0.0 0.0 omydot 0.0],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_2 = 0 accel_i=0
rom12 = calculateROMGyric(feamodel,mesh,el,displ,[0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 omzdot],elStorage,rom0,countedNodes) #Omega_i = 0 #OmegaDot_3 = 0 accel_i=0

#reduced order structural stiffness, mass, and damping
Kr = rom0.Kr
Mr = rom0.Mr
Cr = rom0.Cr + rom0.CrModal
Fr = rom0.Fr

#reduced order transformation matrices
Phi = rom0.Phi
invPhi = rom0.invPhi

#reduced order spin softening coefficient matrices
SrOx2 = rom1.Kr
SrOy2 = rom2.Kr
SrOz2 = rom3.Kr
SrOxOy = rom4.Kr - rom1.Kr - rom2.Kr
SrOyOz = rom5.Kr - rom2.Kr - rom3.Kr
SrOxOz = rom6.Kr - rom1.Kr - rom3.Kr

#reduced order centrifugal load coefficient matrices
FrOx2 = rom1.Fr
FrOy2 = rom2.Fr
FrOz2 = rom3.Fr
FrOxOy = rom4.Fr - rom1.Fr - rom2.Fr
FrOyOz = rom5.Fr - rom2.Fr - rom3.Fr
FrOxOz = rom6.Fr - rom1.Fr - rom3.Fr
FrOxdot = rom10.Fr
FrOydot = rom11.Fr
FrOzdot = rom12.Fr

#need to calculate reduced order acceleration/body force coefeficients
FrAx = rom7.Fr
FrAy = rom8.Fr
FrAz = rom9.Fr

#reduced order gyric coefficient matrices
GrOx = rom1.Cr
GrOy = rom2.Cr
GrOz = rom3.Cr

#reduced order circulatory coefficient matrices
HrOx = 0.5*rom1.Cr
HrOy = 0.5*rom2.Cr
HrOz = 0.5*rom3.Cr

rom = ROM(Kr,Mr,Cr,0.0,Fr,Phi,invPhi,SrOx2,SrOy2,SrOz2,SrOxOy,SrOyOz,SrOxOz,FrOx2,FrOy2,FrOz2,FrOxOy,FrOyOz,FrOxOz,FrOxdot,FrOydot,FrOzdot,FrAx,FrAy,FrAz,GrOx,GrOy,GrOz,HrOx,HrOy,HrOz)
else
rom = ROM(nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing,nothing)
end
return rom
end

Expand Down
Loading