@@ -285,20 +285,17 @@ def compute_all_coefficients(self):
285285 self .cLf = self .compute_forcing_coefficient (
286286 self .cL_0 , self .cL_alpha , self .cL_beta
287287 )
288+ self .cLd = self .compute_damping_coefficient (self .cL_p , self .cL_q , self .cL_r )
288289
289- self .c_q_pitch_f = self .compute_forcing_coefficient (
290- self .cq_pitch_0 , self .cq_pitch_alpha , self .cq_pitch_beta
291- )
292- self .c_q_pitch_d = self .compute_damping_coefficient (
293- self .cq_pitch_p , self .cq_pitch_q , self .cq_pitch_r
290+ self .cQf = self .compute_forcing_coefficient (
291+ self .cQ_0 , self .cQ_alpha , self .cQ_beta
294292 )
293+ self .cQd = self .compute_damping_coefficient (self .cQ_p , self .cQ_q , self .cQ_r )
295294
296- self .c_d_drag_f = self .compute_forcing_coefficient (
297- self .cd_drag_0 , self .cd_drag_alpha , self .cd_drag_beta
298- )
299- self .c_d_drag_d = self .compute_damping_coefficient (
300- self .cd_drag_p , self .cd_drag_q , self .cd_drag_r
295+ self .cDf = self .compute_forcing_coefficient (
296+ self .cD_0 , self .cD_alpha , self .cD_beta
301297 )
298+ self .cDd = self .compute_damping_coefficient (self .cD_p , self .cD_q , self .cD_r )
302299
303300 self .cmf = self .compute_forcing_coefficient (
304301 self .cm_0 , self .cm_alpha , self .cm_beta
@@ -373,19 +370,19 @@ def _compute_from_coefficients(
373370 )
374371
375372 # Compute aerodynamic forces
376- lift = dyn_pressure_area * self .c_l_lift_f (
373+ lift = dyn_pressure_area * self .cLf (
377374 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
378375 ) + dyn_pressure_area_damping * self .cLd (
379376 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
380377 )
381378
382- side = dyn_pressure_area * self .c_q_pitch_f (
379+ side = dyn_pressure_area * self .cQf (
383380 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
384381 ) + dyn_pressure_area_damping * self .cQd (
385382 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
386383 )
387384
388- drag = dyn_pressure_area * self .c_d_drag_f (
385+ drag = dyn_pressure_area * self .cDf (
389386 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
390387 ) + dyn_pressure_area_damping * self .cDd (
391388 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
@@ -410,4 +407,4 @@ def _compute_from_coefficients(
410407 alpha , beta , mach , reynolds , pitch_rate , yaw_rate , roll_rate
411408 )
412409
413- return lift , side , drag , pitch , yaw , roll
410+ return lift , side , drag , pitch , yaw , roll
0 commit comments