ina111

% -----------
% Optimum Design of Nonplanar wings-Minimum Induced Drag 
% for A Given Lift and Wing Root Bending Moment (NAL TR-797)
% 
% Created by Takahiro Inagawa on 2013-2-28.
% Copyright (c) 2013 Takahiro Inagawa. All rights reserved.
% -----------
clear all
tic;  %Computation time measurement

ina111

##概要 JAXA（旧NAL）の1984年の航空宇宙技術研究所報告（NAL TR-797）

「非平面翼の最適設計-揚力と翼根曲げモーメントを与えた時の最小誘導抵抗-」

の数値計算部分をMatlab（Octaveで動作確認）で実装しました． 人力飛行機を念頭にした数値を入れてあります． 論文は汎用的なものなのでConstantの部分の数字を変えれば汎用的に使えます．

論文中では非平面翼の計算を行なっていますが，ここに置いているスクリプトは平面翼です．

ina111

# -*- coding: utf-8 -*-
"""
Mac OSXのmatplotlibで日本語フォント(Osaka)を使うサンプル
python2.7(anaconda-2.1.0)
cf. http://matplotlib.org/users/customizing.html
"""
# デフォルトの文字コードを変更する．
import sys
reload(sys)
sys.setdefaultencoding('utf-8')

ina111

# -*- coding: utf-8 -*-
import sys
reload(sys)
import platform
# デフォルトの文字コードを変更する．
sys.setdefaultencoding('utf-8')

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.font_manager

ina111

ina111

# coding: UTF-8
# 気象庁の高層風のページからデータ取得してプロット
# cf. http://www.data.jma.go.jp/obd/stats/etrn/upper/index.php

import os
import urllib2
from BeautifulSoup import BeautifulSoup
import datetime
import pandas as pd
import numpy as np

ina111

# -*- coding: utf-8 -*-
# ======
# 多段ロケットの最適質量配分（サイジング）問題の計算
# 必要な軌道速度に空力損失、重力損失、推力損失、制御損失を追加し、
# トータルの⊿Vを事前に算出し、その軌道速度に必要なサイジングを行う。
# 初期検討段階にのみ使用可能。
#
# 入力:
#   各段のIsp[秒]
#   各段の構造比(0.0~1.0)（各段の全備重量と推進剤以外の割合）

ina111

# -*- coding: utf-8 -*-
# 衛星高度と可視範囲のグラフ化
#
# cf. 半揚 稔雄(2014) 「ミッション解析と軌道設計の基礎」
# Copyright (c) 2016 Takahiro Inagawa
# Released under the MIT license
import sys
reload(sys)
import platform
sys.setdefaultencoding('utf-8')

ina111

# -*- coding: utf-8 -*-
# Author : Takahiro Inagawa(@ina111)

import io
import numpy as np

from jinja2 import Template

from bokeh.embed import components
from bokeh.models import Range1d

ina111

# -*- coding: utf-8 -*-
# -----------
# Optimum Design of Nonplanar wings-Minimum Induced Drag
# for A Given Lift and Wing Root Bending Moment (NAL TR-797)
#
# Created by Takahiro Inagawa on 2018-03-24.
# Copyright (c) 2018 Takahiro Inagawa. All rights reserved.
# -----------

import numpy as np