hongshaorou/openrocket
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Closes #478

Browse Source 
Clarifies variable names

Removes previous attempt to limit damping moments by using minimum
velocity of 1.0

Clamps damping moments so they cannot exceed actual moments, avoiding
numerical instability
This commit is contained in:
JoePfeiffer 2018-12-02 14:17:33 -07:00
parent 4bd3e4ff18
commit 20fa9925cb
1 changed files with 12 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
core/src/net/sf/openrocket/aerodynamics/BarrowmanCalculator.java
View File

@ -712,17 +712,21 @@ public class BarrowmanCalculator extends AbstractAerodynamicCalculator {
// Calculate pitch and yaw damping moments // Calculate pitch and yaw damping moments
double mul = getDampingMultiplier(configuration, conditions, double mul = getDampingMultiplier(configuration, conditions,
conditions.getPitchCenter().x); conditions.getPitchCenter().x);
double pitch = conditions.getPitchRate(); double pitchRate = conditions.getPitchRate();
double yaw = conditions.getYawRate(); double yawRate = conditions.getYawRate();
double vel = conditions.getVelocity(); double velocity = conditions.getVelocity();
vel = MathUtil.max(vel, 1);
mul *= 3; // TODO: Higher damping yields much more realistic apogee turn mul *= 3; // TODO: Higher damping yields much more realistic apogee turn
total.setPitchDampingMoment(mul * MathUtil.sign(pitch) * pow2(pitch / vel)); // find magnitude of damping moments, and clamp so they can't
total.setYawDampingMoment(mul * MathUtil.sign(yaw) * pow2(yaw / vel)); // exceed magnitude of pitch and yaw moments
double pitchDampingMomentMagnitude = MathUtil.min(mul * pow2(pitchRate / velocity), total.getCm());
double yawDampingMomentMagnitude = MathUtil.min(mul * pow2(yawRate / velocity), total.getCyaw());
// multiply by sign of pitch and yaw rates
total.setPitchDampingMoment(MathUtil.sign(pitchRate) * pitchDampingMomentMagnitude);
total.setYawDampingMoment(MathUtil.sign(yawRate) * yawDampingMomentMagnitude);
} }
// TODO: MEDIUM: Are the rotation etc. being added correctly? sin/cos theta? // TODO: MEDIUM: Are the rotation etc. being added correctly? sin/cos theta?