Projeto Conceitual de Aeronaves de Transporte

Matlab code for DOC and IOC calculation
%=====================================================%
% **** PROGRAM DIRECT OPERATING COST **** %
% BASED ON ROSKAM VOL VIII (Starting at PAG 144) %
% Author: B. Mattos - June 2014 %
%=====================================================%
% Input:
% Time_Block: block time in minutes
% Cons_Block: fuel burned for the stage length (kg)
% var.Range: block length in nm
% Sal.Captain: salary Captain (US$ in a year period)
% Sal.FO: salary first officer (US$ in a year period)
% weight_empty_kg: opertional empty weight (kg)
% weight_engine_kg: weight of a single engine (kg)
% MTOW: maximum takeoff weight (kg)
% NEng: number of engines;
% T0: takeoff thrust of a single engine (lbf)
% TBO: Time between overhaul for the engine (h)
% Fuel_price= fuel price (US$/Gallon)
clc
clear all
%% Constants
kg2lb = 2.20462262; % kg to conversion factor
%%
% ********** Example of airplane (input) data
Time_Block = 120;
Cons_Block = 6500;
var.Range = 900;
Sal.Captain = 85000;
Sal.FO = 50000;
weight_empty_kg = 22000;
weight_engine_kg = 800;
MTOW = 43090;
NEng = 2;
T0 = 13700;
TBO = 2500;
Fuel_price = 2.387;
%==============================%
% MISSION DATA %
%==============================%
tbl = Time_Block/60; % [HRS] BLOCK TIME EQ 5.6
Block_Range = var.Range; % [NM] BLOCK RANGE
vbl = Block_Range / tbl; % [KTS]BLOCK SPEED EQ 5.5
Uannbl = 1e+3 * (3.4546 * tbl + 2.994 - (12.289 * (tbl^2) -
5.6626 * (tbl) + 8.964)^(1/2)); % [HRS] ANNUAL UTILIZATION
%==============================%
% DOC Page 146 %
%==============================%
%==============================%
% DOC FLT %
%==============================%
%1) CREW COST -> Ccrew