Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random | Template Test |  
Login
create new user
Username:
Password:
forget your password?
Main Menu
Sections

Talkback

Downloads

Information
2D with drag airship optimal control (Topic)

2D with Drag and No Wind

Disclaimer: This is a work in progress...

The Force due to drag in 2D is given by

$\displaystyle F_d = \frac{1}{2} C_d \rho A v^2$

Unlike the 1D case, which uses the absolute value, the direction of the drag force in 2D is taken into account by the magnitude of the velocity times the component velocities such that the x and y components of the drag force are

$\displaystyle F^x_d = -\frac{1}{2} C_d \rho A \dot{x} \sqrt{\dot{x}^2 + \dot{y}^2}$

$\displaystyle F^y_d = -\frac{1}{2} C_d \rho A \dot{y} \sqrt{\dot{x}^2 + \dot{y}^2}$

The equations of motion for this problem are

$\displaystyle \dot{x} = \dot{x}$

$\displaystyle \dot{y} = \dot{y}$

$\displaystyle \ddot{x} = F^x_d + \frac{T}{m} \cos \alpha$

$\displaystyle \ddot{y} = F^y_d + \frac{T}{m} \sin \alpha$

where T is the applied Torque in Newtons, $\alpha$ is the angle of applied torque counter clockwise from the negative x-axis and $m$ is the mass of the airship.

The performance measure for the minimum energy problem assuming there is no effort required to change the angle of attack

$\displaystyle J = \int_0^{t_f} \frac{T^2}{m^2} dt$ (1)

The Hamiltonian for this optimal control problem

$\displaystyle \mathcal{H} = \frac{T^2}{m^2} + p_1 \dot{x} + p_2 \dot{y} + p_3 \... ...}{m} \cos \alpha \right) + p_4 \left ( F^y_d + \frac{T}{m} \sin \alpha \right )$ (2)
The neccessary conditions for the unconstrained control inputs

$\displaystyle \frac{\partial \mathcal{H}}{\partial \alpha} = 0 $

$\displaystyle \frac{\partial \mathcal{H}}{\partial T} = 0 $

$\displaystyle \frac{\partial \mathcal{H}}{\partial \alpha} = -p_3 \frac{T}{m} \sin \alpha + p_4 \frac{T}{m} \cos \alpha = 0$

$\displaystyle \frac{p_4 T}{m}\cos \alpha = \frac{p_3 T}{m} \sin \alpha $

$\displaystyle \frac{p_4}{p_3} = \frac{\sin \alpha}{\cos \alpha} $
$\displaystyle \tan \alpha = \frac{p_4}{p_3}$ (3)

$\displaystyle \frac{\partial \mathcal{H}}{\partial T} = \frac{2T}{m^2} + \frac{p_3}{m} \cos{\alpha} + \frac{p_4}{m}{sin \alpha} = 0 $

$\displaystyle \frac{2T}{m} + p_3 \cos{\alpha} + p_4{sin \alpha} = 0 $

$\displaystyle \frac{2T}{m} = -\left (p_3 \cos{\alpha} + p_4{sin \alpha} \right ) $
$\displaystyle T = -\frac{m}{2}\left (p_3 \cos{\alpha} + p_4{sin \alpha} \right )$ (4)

The neccessary conditions for the costates are

$\displaystyle \dot{p_1} = -\frac{\partial \mathcal{H}}{\partial x} = 0$ (5)
$\displaystyle \dot{p_2} = -\frac{\partial \mathcal{H}}{\partial y} = 0$ (6)
$\displaystyle \dot{p_3} = -\frac{\partial \mathcal{H}}{\partial \dot{x}} = -p_1... ...ot{y}^2}} + \frac{p_4 C_d \rho A \dot{x} \dot{y}}{\sqrt{\dot{x}^2 + \dot{y}^2}}$ (7)
$\displaystyle \dot{p_4} = -\frac{\partial \mathcal{H}}{\partial \dot{y}} = -p_2... ...ot{y}^2}} + \frac{p_3 C_d \rho A \dot{x} \dot{y}}{\sqrt{\dot{x}^2 + \dot{y}^2}}$ (8)



"2D with drag airship optimal control" is owned by bloftin.

View style:


Cross-references: Hamiltonian, energy, mass, motion, velocity, magnitude, drag force, work

This is version 1 of 2D with drag airship optimal control, born on 2006-07-25.
Object id is 210, canonical name is 2DWithDragAirshipOptimalControl.
Accessed 1459 times total.

Classification:
Physics Classification: 02.30.Yy (Control theory)
Pending Errata and Addenda
None.
Discussion
Style: Expand: Order:

No messages.

Testing some escape charachters for html category with a generator has an injective cogenerator" now escape ” with "