Nejprve dynamicke vlastnosti leteckeho nosniku
/* airfoil-dynamical properties */ double airfoil_m; double airfoil_Kh; double airfoil_Ka; double airfoil_Ia; double airfoil_Sa; /* damping */ double damping_h,damping_a; /* prescribed motions */ double frequency1,amplitude1; double frequency2,amplitude2; /* static forces */ double force_lift,force_moment; |
Dale jeho aktualni poloha.
double releasetime =5.0; double cofgheight =0.0; double angle =0.0; double yvelocity =0.0; double angularvelocity =0.0; /* init location */ double restartangle=0.0; double restarty=0.0; |
Pro ODE resic ulozime historii pohybu i aerodynamickych sil.
/* ODE SOLVER HISTORY */
/* reseni systemu ODR 1.radu implicitni metodou */
double xlast[HISTORY_LIFT+1][4];
/* Koeficienty cy a cmy */
double DragForce[HISTORY_LIFT+1],
LiftForce[HISTORY_LIFT+1],
TorsionalMoment[HISTORY_LIFT+1];
|
Modifikace oblasti - oblast je modifikovana tak, ze vstup, vystup zustavaji nemenne, a deformuje se pouze oblast nosniku. Zobrazeni mezi je kombinaci identity, otaceni a posunuti.
/* MODIFICATION OF DOMAIN */ double ALE_R1=1.0; double ALE_R2=2.0; /* soucasna elasticke osy */ double xmcenter,ymcenter; /* puvodni poloha elasticke osy */ double xhc,yhc; |