fly/flightgear
1
0
Fork 0
Code Activity

Updates to Tony's c172 model.

Browse source
This commit is contained in:
curt 1999-08-07 18:23:55 +00:00
parent b001cc55c7
commit e6aef49822
2 changed files with 12 additions and 11 deletions
src/FDM/LaRCsim
c172_aero.cc172_gear.c

19
src/FDM/LaRCsim/c172_aero.c
View file

@ -152,9 +152,11 @@ void aero( SCALAR dt, int Initialize ) {
/*Note that CLo,Cdo,Cmo will likely change with flap setting so /*Note that CLo,Cdo,Cmo will likely change with flap setting so
they may not be declared static in the future */ they may not be declared static in the future */
if (Initialize != 0) /* printf("Initialize= %d\n",Initialize); */
{ /* if (Initialize != 0)
CLadot=1.7; { */
/* printf("Initializing aero model...Initialize= %d\n", Initialize);
*/ CLadot=1.7;
CLq=3.9; CLq=3.9;
CLde=0.43; CLde=0.43;
CLo=0; CLo=0;
@ -194,7 +196,7 @@ void aero( SCALAR dt, int Initialize ) {
b=35.8; /*wing span ft */ b=35.8; /*wing span ft */
Sw=174; /*wing planform surface area ft^2*/ Sw=174; /*wing planform surface area ft^2*/
rPiARe=0.054; /*reciprocal of Pi*AR*e*/ rPiARe=0.054; /*reciprocal of Pi*AR*e*/
} /* } */
/* /*
LaRCsim uses: LaRCsim uses:
@ -249,9 +251,9 @@ void aero( SCALAR dt, int Initialize ) {
qSb=qS*b; qSb=qS*b;
/* printf("Wb: %7.4f, Ub: %7.4f, Alpha: %7.4f, elev: %7.4f, ail: %7.4f, rud: %7.4f, long_trim: %7.4f\n",W_body,U_body,Alpha*RAD_TO_DEG,elevator*RAD_TO_DEG,aileron*RAD_TO_DEG,rudder*RAD_TO_DEG,long_trim*RAD_TO_DEG); /* printf("aero: Wb: %7.4f, Ub: %7.4f, Alpha: %7.4f, elev: %7.4f, ail: %7.4f, rud: %7.4f, long_trim: %7.4f\n",W_body,U_body,Alpha*RAD_TO_DEG,elevator*RAD_TO_DEG,aileron*RAD_TO_DEG,rudder*RAD_TO_DEG,long_trim*RAD_TO_DEG);
printf("Theta: %7.4f, Gamma: %7.4f, Beta: %7.4f, Phi: %7.4f, Psi: %7.4f\n",Theta*RAD_TO_DEG,Gamma_vert_rad*RAD_TO_DEG,Beta*RAD_TO_DEG,Phi*RAD_TO_DEG,Psi*RAD_TO_DEG); */ //printf("Theta: %7.4f, Gamma: %7.4f, Beta: %7.4f, Phi: %7.4f, Psi: %7.4f\n",Theta*RAD_TO_DEG,Gamma_vert_rad*RAD_TO_DEG,Beta*RAD_TO_DEG,Phi*RAD_TO_DEG,Psi*RAD_TO_DEG);
*/
/* sum coefficients */ /* sum coefficients */
CLwbh = interp(CLtable,alpha_ind,NCL,Alpha); CLwbh = interp(CLtable,alpha_ind,NCL,Alpha);
@ -263,7 +265,7 @@ void aero( SCALAR dt, int Initialize ) {
cn = Cnbeta*Beta + (Cnp*P_body + Cnr*R_body)*b_2V + Cnda*aileron + Cndr*rudder; cn = Cnbeta*Beta + (Cnp*P_body + Cnr*R_body)*b_2V + Cnda*aileron + Cndr*rudder;
croll=Clbeta*Beta + (Clp*P_body + Clr*R_body)*b_2V + Clda*aileron + Cldr*rudder; croll=Clbeta*Beta + (Clp*P_body + Clr*R_body)*b_2V + Clda*aileron + Cldr*rudder;
/* printf("CL: %7.4f, Cd: %7.4f, Cm: %7.4f, Cy: %7.4f, Cn: %7.4f, Cl: %7.4f\n",CL,cd,cm,cy,cn,croll); /* printf("aero: CL: %7.4f, Cd: %7.4f, Cm: %7.4f, Cy: %7.4f, Cn: %7.4f, Cl: %7.4f\n",CL,cd,cm,cy,cn,croll);
*/ /*calculate wind axes forces*/ */ /*calculate wind axes forces*/
F_X_wind=-1*cd*qS; F_X_wind=-1*cd*qS;
F_Y_wind=cy*qS; F_Y_wind=cy*qS;
@ -291,4 +293,3 @@ void aero( SCALAR dt, int Initialize ) {
*//* printf("Maero: %7.4f Naero: %7.4f Raero: %7.4f\n",M_m_aero,M_n_aero,M_l_aero); *//* printf("Maero: %7.4f Naero: %7.4f Raero: %7.4f\n",M_m_aero,M_n_aero,M_l_aero);
*/ */
} }

4
src/FDM/LaRCsim/c172_gear.c
View file

@ -37,8 +37,8 @@
$Header$ $Header$
$Log$ $Log$
Revision 1.3 1999/07/31 02:57:36 curt Revision 1.4 1999/08/07 16:23:55 curt
Improvements to Tony's c172 model. Updates to Tony's c172 model.
Revision 1.1.1.1 1999/04/05 21:32:45 curt Revision 1.1.1.1 1999/04/05 21:32:45 curt
Start of 0.6.x branch. Start of 0.6.x branch.