##概要 JAXA(旧NAL)の1984年の航空宇宙技術研究所報告(NAL TR-797)
「非平面翼の最適設計-揚力と翼根曲げモーメントを与えた時の最小誘導抵抗-」
の数値計算部分をMatlab(Octaveで動作確認)で実装しました. 人力飛行機を念頭にした数値を入れてあります. 論文は汎用的なものなのでConstantの部分の数字を変えれば汎用的に使えます.
論文中では非平面翼の計算を行なっていますが,ここに置いているスクリプトは平面翼です.
import os | |
import requests | |
from dotenv import load_dotenv | |
from slack_bolt import App | |
from slack_bolt.adapter.socket_mode import SocketModeHandler | |
import openai | |
load_dotenv() | |
app = App(token=os.environ["SLACK_BOT_TOKEN"]) | |
openai.api_key = os.environ["OPENAI_API_KEY"] |
import os | |
from dotenv import load_dotenv | |
from slack_bolt import App | |
from slack_bolt.adapter.socket_mode import SocketModeHandler | |
import openai | |
load_dotenv() | |
app = App(token=os.environ["SLACK_BOT_TOKEN"]) | |
openai.api_key = os.environ["OPENAI_API_KEY"] |
""" | |
軌道投入への必要ΔVの簡易計算スクリプト | |
""" | |
import sys | |
import os | |
import configparser | |
import numpy as np | |
from numpy import sqrt, cos, sin, arcsin, tan, deg2rad | |
import matplotlib.pyplot as plt |
# -*- coding: utf-8 -*- | |
# ====== | |
# 多段ロケットの最適質量配分(サイジング)問題の計算 | |
# 必要な軌道速度に空力損失、重力損失、推力損失、制御損失を追加し、 | |
# トータルの⊿Vを事前に算出し、その軌道速度に必要なサイジングを行う。 | |
# 初期検討段階にのみ使用可能。 | |
# | |
# 入力: | |
# 各段のIsp[秒] | |
# 各段の構造比(0.0~1.0)(各段の全備重量と推進剤以外の割合) |
# 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 |
{ | |
"cells": [ | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"## REFPROPのPythonラッパーの使い方\n", | |
"\n", | |
"本家NISTから公開されているラッパーの使い方のjupyter notebookはこちら\n", | |
"* https://nbviewer.jupyter.org/github/usnistgov/REFPROP-wrappers/blob/master/wrappers/python/notebooks/Tutorial.ipynb\n", |
##概要 JAXA(旧NAL)の1984年の航空宇宙技術研究所報告(NAL TR-797)
「非平面翼の最適設計-揚力と翼根曲げモーメントを与えた時の最小誘導抵抗-」
の数値計算部分をMatlab(Octaveで動作確認)で実装しました. 人力飛行機を念頭にした数値を入れてあります. 論文は汎用的なものなのでConstantの部分の数字を変えれば汎用的に使えます.
論文中では非平面翼の計算を行なっていますが,ここに置いているスクリプトは平面翼です.
# -*- 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 |
% ----------- | |
% 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 |