書名 著者 刊行 備考 コンピュータによる流体力学 ファーツィガー、ペリッチ 2002 FDMおよびFDM コンピュータによる熱移動と流れの数値解析 パタンカー 1984 FVM、SIMPLE。刊行時点でのFEMの流体計算への応用は未知 乱流の数値シミュレーション 梶島 1999 FD…

www.youtube.com

クロスパワースペクトル密度Φxyをパワースペクトル密度Φxxで割ると伝達関数が求まる。小型無人航空機制御用アクチュエータ伝達関数測定法の研究 cir.nii.ac.jp

モーメント方程式（回転運動方程式） 航法方程式（地球座標上の移動速度） 力の方程式（並進運動方程式） 運動学方程式（姿勢角速度） 関係図

以下の関数が便利 ・quaternion(q0,q1,q2,qe)：コンストラクタ ・conj()：共役クォータニオン。単位クォータニオンの場合は、逆クォータニオンと等しくなる ・parts()：クォータニオンから要素を取り出す。[q0,q1,q2,q3]=parts(quat)のように使う。[,q1,q2,q…

対流圏と中間圏では温度が高度の上昇にともなって下がる。 成層圏では逆に温度が高度の上昇にともなって上がる。 気圧は指数的に減衰している

clear; close all; g=32.174;%ヤード・ポンド法なので9.8ではない Za=-236; Zu=-0.227; Xa=14.7; Xu=-0.0215; Madot=-0.28; Ma=-3.78; Mu=0; Mq=-0.992; Zde=-12.9; Mde=-2.48; U0=293.8; Zq=-5.76; global A; global B; A=[Xu Xa 0 -g; Zu/U0,Za/U0 (U0+Zq)…

縦の運動方程式 行列表示に変形 安定軸に対する安定微係数（可変特性研究機: ロッキードP2V-7をベース） clear; close all; g=32.174;%ヤード・ポンド法なので9.8ではない Za=-236; Zu=-0.227; Xa=14.7; Xu=-0.0215; Madot=-0.28; Ma=-3.78; Mu=0; Mq=-0.992…

局所水平座標系（Local level coordinate/Local horizontal coordinate） NED（north-east-down）座標系、ENU(east-north-up)座標系などがある。 通常航空機の解析ではNED座標系を用いる。

川口、航空・宇宙における制御 p.121 NED座標系 機体の並進,回転,曲げ振動の式 宮澤、宇宙ロケット工学入門 p.42 質点の運動方程式 冨田 et al.、ロケット工学基礎講義 p.104：質点の運動方程式 p.106：エンジン推力,空気力がなく,重力のみが作用する質点の…

以下の２つを先のシミュレーションにマージ butterfly-effect.hatenablog.com butterfly-effect.hatenablog.com clear;cla global A; global U0; %有次元安定微係数 Xu = -0.01; Zu = -0.1; Mu = 0.001; Xa = 30; Za = -200; Ma = -4; Xq = 0.3; Zq = -5; Mq…

ここのコードを改変した。 butterfly-effect.hatenablog.com clear;cla global A; %有次元安定微係数 Xu = -0.01; Zu = -0.1; Mu = 0.001; Xa = 30; Za = -200; Ma = -4; Xq = 0.3; Zq = -5; Mq = -1; Yb = -45; Lb_= -2; Nb_= 1; Yp = 0.5; Lp_= -1; Np_= -…

www.mathworks.com load simdata h = Aero.Animation; h.createBody('pa24-250_orange.ac','Ac3d'); h.Bodies{1}.TimeSeriesSource = simdata; h.Camera.PositionFcn = @staticCameraPosition; h.Figure.Position(1) = h.Figure.Position(1) + 572/2; h.show…

１）慣性主軸（Principal axes） X軸を慣性主軸の１つに一致させる。このとき、Ixz=0となる。 ２）安定軸（Stability axes） X軸を左右対称な定常釣り合い飛行時の機体速度ベクトルVc0に一致させる。 運動方程式を線形化する際に、無次元安定微係数の計算が…

butterfly-effect.hatenablog.com

ロケットの動的不安定：ロケットの推力重量比が5:1以下の場合に発生する ・N-1ロケット 重量: 約640g=0.64kg 推力(N-1 motor v10): 平均推力75N~7.5kg推力重量比〜11 動的不安定性はN-1では問題なさそう

メモ ・Zoom(授業後に録画をyoutubeで共有) ・Keynote ・Ink2Go（+Wacomタブレット） ・ScreenPointer（マウスカーソルのハイライト）（左右両方のcommandキー⌘を長押しでマウスのハイライトをトグル）