fly/flightgear
1
0
Fork 0
Code Activity

Fade out the turbulence within two wingspans of the ground.

Browse source
This commit is contained in:
david 2003-02-12 18:34:39 +00:00
parent 4ecd146258
commit 9f5498acf4
1 changed files with 14 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

17
src/FDM/JSBSim/FGAtmosphere.cpp
View file

@ -261,7 +261,7 @@ void FGAtmosphere::Calculate(double altitude)
void FGAtmosphere::Turbulence(void) void FGAtmosphere::Turbulence(void)
{ {
switch (turbType) { switch (turbType) {
case ttBerndt: case ttBerndt: {
vDirectiondAccelDt(eX) = 1 - 2.0*(((double)(rand()))/RAND_MAX); vDirectiondAccelDt(eX) = 1 - 2.0*(((double)(rand()))/RAND_MAX);
vDirectiondAccelDt(eY) = 1 - 2.0*(((double)(rand()))/RAND_MAX); vDirectiondAccelDt(eY) = 1 - 2.0*(((double)(rand()))/RAND_MAX);
vDirectiondAccelDt(eZ) = 1 - 2.0*(((double)(rand()))/RAND_MAX); vDirectiondAccelDt(eZ) = 1 - 2.0*(((double)(rand()))/RAND_MAX);
@ -270,6 +270,16 @@ void FGAtmosphere::Turbulence(void)
MagnitudedAccelDt = 1 - 2.0*(((double)(rand()))/RAND_MAX) - Magnitude; MagnitudedAccelDt = 1 - 2.0*(((double)(rand()))/RAND_MAX) - Magnitude;
MagnitudeAccel += MagnitudedAccelDt*rate*State->Getdt(); MagnitudeAccel += MagnitudedAccelDt*rate*State->Getdt();
Magnitude += MagnitudeAccel*rate*State->Getdt(); Magnitude += MagnitudeAccel*rate*State->Getdt();
// Fade the magnitude within two wingspans
// of the ground (WAG)
double AdjustedMagnitude = Magnitude;
double AdjustedMagnitudeAccel = MagnitudeAccel;
double HOverBMAC = Position->GetHOverBMAC();
if (HOverBMAC < 2.0) {
AdjustedMagnitude *= (HOverBMAC / 2.0);
AdjustedMagnitudeAccel *= (HOverBMAC / 2.0);
}
vDirectiondAccelDt.Normalize(); vDirectiondAccelDt.Normalize();
vDirectionAccel += vDirectiondAccelDt*rate*State->Getdt(); vDirectionAccel += vDirectiondAccelDt*rate*State->Getdt();
@ -277,8 +287,8 @@ void FGAtmosphere::Turbulence(void)
vDirection += vDirectionAccel*rate*State->Getdt(); vDirection += vDirectionAccel*rate*State->Getdt();
vDirection.Normalize(); vDirection.Normalize();
vTurbulence = TurbGain*Magnitude * vDirection; vTurbulence = TurbGain*AdjustedMagnitude * vDirection;
vTurbulenceGrad = TurbGain*MagnitudeAccel * vDirection; vTurbulenceGrad = TurbGain*AdjustedMagnitudeAccel * vDirection;
vBodyTurbGrad = State->GetTl2b()*vTurbulenceGrad; vBodyTurbGrad = State->GetTl2b()*vTurbulenceGrad;
vTurbPQR(eP) = vBodyTurbGrad(eY)/Aircraft->GetWingSpan(); vTurbPQR(eP) = vBodyTurbGrad(eY)/Aircraft->GetWingSpan();
@ -293,6 +303,7 @@ void FGAtmosphere::Turbulence(void)
vTurbPQR(eR) = vBodyTurbGrad(eX)/10.0; vTurbPQR(eR) = vBodyTurbGrad(eX)/10.0;
break; break;
}
default: default:
break; break;
} }