Jump to content

Ducted Fan: Difference between revisions

From mintOC
Line 40: Line 40:
{| border="1" align="center" cellpadding="5" cellspacing="0"
{| border="1" align="center" cellpadding="5" cellspacing="0"
|- bgcolor=#c7c7c7
|- bgcolor=#c7c7c7
! Symbol !! Value !! Description
! Symbol !! Value
|-
|-
| align=center | <math>m_1</math> || align=right | 100 <math>\mathrm{kg}</math> || First mass directly affected by <math>F( t )</math>
| align=center | <math>m</math> || align=right | 2.2 <math>\mathrm{kg}</math>  
|-  
|-  
| align=center | <math>m_2</math> || align=right | 2 <math>\mathrm{kg}</math>||  Second mass influenced by damping control
| align=center | <math>J</math> || align=right | 0.05 <math>\mathrm{kg \cdot m^2}</math>
|-
|-
| align=center | <math>k_1</math> || align=right | 100 <math>\mathrm{N/m}</math>|| Spring connecting first mass to reference
| align=center | <math>r</math> || align=right | 0.2 <math>\mathrm{m}</math>
|-  
|-  
| align=center | <math>k_2</math> || align=right | 3 <math>\mathrm{N/m}</math>|| Coupling spring between the two masses
| align=center | <math>\mathrm{mg}</math> || align=right | 0.4 <math>\mathrm{N}</math>
|-  
|-  
| align=center | <math>c</math> || align=right | 0.5 <math>\mathrm{Ns/m}</math>|| Damping affecting second mass
| align=center | <math>\mu</math> || align=right | 1
|-
| align=center | <math>T</math> || align=right | <math>2 \pi</math> || Duration of the motion
|-
| align=center | <math>u</math> || align=right | - || Modulates the damping of the second mass
|}
|}



Revision as of 10:15, 24 November 2025

Ducted Fan
State dimension: 1
Differential states: 4
Discrete control functions: 1


The Ducted Fan problem is a classical nonlinear benchmark in optimal control with multiple input and state constraints. This description is taken from [1].

It models the planar motion of a ducted fan aircraft, described by its horizontal and vertical positions (x1,x2) , the angle α with respect to the vertical, and their velocities (v1,v2,vα) . The inputs are the body-fixed thrust components (u1,u2) , generated by moving flaps at the end of the duct.

The objective is to steer the fan from the origin to a horizontal position of 1m at altitude 0, with zero final velocities and attitude, in a free final time tf, while minimising a trade-off between control effort and transition time.

Mathematical formulation

We summarize the states as x:=(x1,v1,x2,v2,α,vα).

minu,tf1tf0tf(2u12(t)+u22(t))dt+μtfsubject tox1˙(t)=v1(t),v1˙(t)=1m(u1cosαu2sinα),x2˙(t)=v2(t),v2˙(t)=1m(mg+u1sinα+u2cosα),α˙=vα,v˙α=rJu1,x(0)=(0,0,0,0,0,0)T,x(tf)=(1,0,0,0,0,0)T,u1(t)[5,5] t[0,tf],u2(t)[0,17] t[0,tf],α(t)[30,30] t[0,tf]

Parameters

These fixed values are used within the model:

Symbol Value
m 2.2 kg
J 0.05 kgm2
r 0.2 m
mg 0.4 N
μ 1

Reference Solutions

Here is one local solution to the above control problem.

Miscellaneous and Further Reading

This formulation and a detailed description can be found in [1].

References

[1] Caillau, J.-B., Cots, O., Gergaud, J., & Martinon, P. OptimalControlProblems.jl: a collection of optimal control problems with ODE's in Julia. https://github.com/control-toolbox/OptimalControlProblems.jl/blob/main/ext/Descriptions/double_oscillator.md