ロケット工学基礎講義 p.156からの内容 計算アルゴリズム（Backward sweep）は、非線形最適制御（大塚）の数値計算の章(p.142)を参考にした。 座標の定義 状態方程式 加速度（Cは有効排気速度） 終端条件 評価関数 ハミルトニアン 随伴変数の微分方程式 最適…

初速V0＝10で鉛直上向きにボールを投げ上げる問題。 2秒間のシミュレーションを行う。 初期条件はx(0)=0, v(0)=10。 clear; close all; clc; x0 = 0; v0 = 10; %m/s t_final = 2; %s tspan = linspace(0,t_final,20); % 順方向に解く [t,xnext] = ode45(@myo…

www.youtube.com

www.nar.org

www.youtube.comferrolho.github.iogithub.com

www.youtube.com www.gpops2.com github.com www.psopt.net

clear; close all; clc; figure(); hold on; a0_g = linspace(1,3.5,50); Lambdas = [2 2.5 3 3.5 4]; %質量比の逆数 y_pks = []; a0_g_peaks = []; for idx=1:length(Lambdas) Lambda = Lambdas(idx); y_g_Isp2 = 1./a0_g.*(1+... -1/Lambda*(1+log(Lambda)…

clear; close all; clc; C = 4413; %m/s lambdas = [0.025, 0.05]; %ペイロード比 2.5%, 5% epsilon = linspace(0.85,0.95,10); %推進薬搭載率 85%~95% figure(); hold on; for idx=1:length(lambdas) lambda = lambdas(idx); deltaV = C*log(1./(1-(1-lambd…

OpenMotorで105mm x 50mmのKNSBグレインを１個〜４個に変化させたときのシミュレーションを行った。 グレインを4個以上にするとスロートでの熱流束が制限値を超えるので、４個を上限とした。 トータルインパルスはグレインに比例して線形に増加するが、３個…

ここを参考に。 https://www.nycaerospace.org/static/design1.pdf clear; close all; clc; % Rocket details z = 500; %m; g = 9.8; %m/s2 deltaV = sqrt(2*g*z); %m/s T_final= 15-6.5*z/1000; %deg C a_final = 20.0468*sqrt(T_final+273.15); %m/s M_fin…

www.uscrpl.com

www.youtube.com www.youtube.com

dronerice.jp

Survey3（MAPIR社製） www.nngists.jp

www.mapir.camera

lygte-info.dk一覧表 lygte-info.dk

ロケットの中に入ってInhibitorを装着しているのが面白かった。 www.youtube.com あと、サムネにもある表面に並べられた多数の筒状の部品が気になった。

HPR Level 1 www.nar.orgHPR Level 2 www.nar.orgHPR Level 3 www.nar.org

qd-suriken.com

「基礎から学べる機械力学」（伊藤）を参考に clear; close all; clc; e = 1; wn =1 ; t = linspace(0,10,30); x = e*wn/2.*t.*sin(wn*t); y = -e*wn/2.*t.*cos(wn*t)+e/2*sin(wn*t); plot(t,x); hold on; plot(t,y); legend('x(t)','y(t)') xlabel('time …

clear; close all; clc; E = 206e9; d = 0.1; I = pi/64*d^4; m = 30; l = 2; a = l/2; b = l/2; k = 6*E*I*l/(a*b*(l^2-a^2-b^2)); nc = 60/(2*pi)*sqrt(k/m); fprintf('危険速度: %.0f [rpm]\n',nc); 実行結果 危険速度: 4294 [rpm]

等温壁 Illingworth_Stewartson_transformation.m clear; close all; clc; eta = linspace(0, 20, 50); fprimeprime_tmps = linspace(0.45,0.48,1000); errors = []; for idx=1:length(fprimeprime_tmps) fprimeprime_tmp = fprimeprime_tmps(idx); [eta2, x…

arrayfun('UniformOutput',false)とsprintf()を使う Mes = [0, 0.8, sqrt(5), 5]; legend(arrayfun(@(x) sprintf('Me=%.2f',x),Mes, 'UniformOutput', false));

PubChemの検索結果のSynonym欄に書いてあるっぽい。 pubchem.ncbi.nlm.nih.gov 以下だとHSコードもまとめて調べられる。と思ったら有料っぽい。 www.findhs.codesHSコードはSchedule Bコードとも呼び、以下で調べられる。 uscensus.prod.3ceonline.com

以下の板に荷重を加える。荷重は端点に200N。素材は鋼（デフォルト）。 メッシュ 境界条件（端面拘束。反対側に200N） 静的解析 最大応力は51MPa. モード解析 固有周波数は以下。 Harmonic analysisでは50Hz〜910Hzまでスイープして調べる。 Harmonic analys…

clear; close all; clc; L = 50e-2; %m d = 7e-3; %m % dx = 10e-2; dy = 10e-2; nx = 20; ny = 20; dx = L/nx; dy = d/ny; mu = 1.83e-5; %Pa rho = 1.29; %kg/m3 Uinf = 20; %m/s u = zeros(nx,ny); v = zeros(nx,ny); u(1,:) = Uinf; %入り口のX速度はUin…

以下の３階常微分方程式を解く 境界条件（教科書に誤植を発見） 以下の連立一階常微分方程式を解く（Matlabではode23sを使った） clear; close all; clc; % シューティング法による境界層方程式の解法 eta_max = 7; eta = linspace(0,eta_max,100); c_initia…

clear; close all; clc; C = 1; mu = 1; rho = 1; tmax = 5; ymax = 1; ts = [0.01, 0.1]; figure(1); hold on; figure(2); hold on; for idx=1:length(ts) t = ts(idx); y = linspace(0,ymax,50); dudy = C./sqrt(mu*t/rho).*exp(-y.^2./(4*mu*t)*rho); dud…

www.mathworks.com www.youtube.comwww.youtube.comwww.youtube.com

以下で使われているGrimech3.0というソフトウェアが何なのか調べていた。 www.youtube.com UCバークレーで開発されているソフトらしい。 combustion.berkeley.edu どうもChemkinというソフトウェア上で実行されるらしいが、有料らしい。 [https://www.ansys.…