From 7e4205fa63c0474739fc15b61a52afc344b84e1d Mon Sep 17 00:00:00 2001 From: Uwe Fechner Date: Sat, 23 Nov 2024 15:34:08 +0100 Subject: [PATCH] do not use the norm of the pitch rate --- src/KPS4.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/KPS4.jl b/src/KPS4.jl index 5b564ed9..842d5fe2 100644 --- a/src/KPS4.jl +++ b/src/KPS4.jl @@ -370,9 +370,9 @@ Updates the vector s.forces of the first parameter. s.forces[s.set.segments + 3] .+= (L2 + D2) end # if the impact of cmq is just drag - f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * norm(s.pitch_rate) * s.set.cord_length) * va_xz1 + # f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * norm(s.pitch_rate) * s.set.cord_length) * va_xz1 # otherwise - # f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * s.pitch_rate * s.set.cord_length) * va_xz1 + f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * s.pitch_rate * s.set.cord_length) * va_xz1 s.forces[s.set.segments + 2] .+= f_d s.forces[s.set.segments + 4] .+= (L3 + D3 -0.5*f_d) s.forces[s.set.segments + 5] .+= (L4 + D4 -0.5*f_d)