2022-03-01から1ヶ月間の記事一覧
www.cnet.com
private JComponent NakujaTab() { JSplitPane horizontal = new JSplitPane(JSplitPane.HORIZONTAL_SPLIT, true); horizontal.setResizeWeight(0.5); JPanel panel = new JPanel(new MigLayout("fill, flowy", "[grow][grow 0]","[grow]")); TextArea textA…
charleslabs.fr
clear; close all; clc; tau1 = 9; sampling_time = 0.01; s = tf('s'); G = 1/tau1*(1-exp(-tau1*s))/s; H = c2d(G,sampling_time,'zoh'); [A,B1,B2,C1,C2,D11,D12,D21,D22,E,tau] = getDelayModel(H); max_simulation_time = 20; %[s] kmax = round(max_si…
clear; close all; clc; tau1 = 9; s = tf('s'); G = 1/tau1*(1-exp(-tau1*s))/s; H = c2d(G,0.01,'zoh'); H2 = getDelayModel(H,'lft'); [A,B,C,D,Ts]=ssdata(H2) step(H); 以下のステップ応答をもつ伝達関数をモデリングした。 実際はこれは移動平均フィル…
github.comgithub.com
www.youtube.com
26:42あたりで一瞬触れられている。 youtu.be
1. eclipseをインストール 2. JDK 11をインストール(Java17では動かない) 3. OpenRocketをクローン 4. eclipseでOpenRocketフォルダをインポート 5. エラーが出るので修正する。自分の場合はjson-xxx.jarとjavax.json-xxx.jarをダウンロードして読み込んだ…
datumstudio.jp以下でインストール pip install japanize-matplotlib使い方は単にインポートするだけ import japanize_matplotlib
筑波大学 熱流体制御・計測研究室 www.kz.tsukuba.ac.jp
www.youtube.com設計ツール Quicksmith quicksmith.online
https://elearning.cerfacs.fr/combustion/index.php
discussions.apple.com
公式が作成しているらしく、かなり出来が良い。 imechanica.org
clear; close all; clc; L = 1; x = linspace(0/L,L/L,100); t = [0.01 0.02 0.05 0.10 0.20]; m = 0; sol = pdepe(m,@heatpde,@heatic,@heatbc,x,t); plot(sol(1,:)); hold on; plot(sol(2,:)); plot(sol(3,:)); plot(sol(4,:)); plot(sol(5,:)); legend(st…
clear; close all; clc; L = 1; x = linspace(0/L,L/L,100); t = [0.01 0.02 0.05 0.10 0.20]; m = 0; sol = pdepe(m,@heatpde,@heatic,@heatbc,x,t); plot(sol(1,:)); hold on; plot(sol(2,:)); plot(sol(3,:)); plot(sol(4,:)); plot(sol(5,:)); legend(st…
以下を参考に。 www.mathworks.com 方程式の定義 以下の熱伝導方程式を解く 上記を以下と比較し、パラメータm、f,、c、sを決定する。 コード(mはメインの関数で用いるのでここでは使わない) function [c,f,s] = heatpde(x,t,u,dudx) c = 1; f = dudx; s = …
燃焼ガスの熱伝達係数はNakkaの値(1000)を用いた。 www.nakka-rocketry.net clear; close all; clc; L = 2e-3; %[m]. 2mm Tin = 1000; %K Tinf = 300; %K alpha_in = 1000;% 燃焼ガスの熱伝達係数 by Nakka % https://www.nakka-rocketry.net/therm.html#:…
Fluentでの操作のポイントは、Reference valueを設定すること。 以下のページに様々な形状の抗力係数が示されている。 skomo.o.oo7.jp 計算 縦に置かれた円柱、l/d=1の場合を試してみる。うまくいけばCd=0.63になるはず。 条件 円柱のd=10mmとする。奥行きL…
似ているようで異なる。 レイリー流れ ・断面積が一定の管内 ・定常流れ ・摩擦なし ・熱の授受はあり の流れをレイリー流れ(Rayleigh flow)と呼ぶ。 摩擦がないので可逆流れである。 ※垂直衝撃波はレイリー流れとファノ流れの両方の条件を満たす レイリー…
www.youtube.com www.youtube.com
Transientを選択するだけではアニメーションが保存されない。 Solution>Calculation activitiesでアニメーションの保存タイムステップを定義する必要がある。 youtu.beyoutu.beyoutu.beyoutu.be
clear; close all; clc; % N-2ロケット 固体推進剤 KNSB Pc = 1.31; %[MPa] Nakkaのグラフの値。OpenMotorでも同様の値(KNSB)。 k = 1.177; %Propep3の計算値(KNSB), O:F=65:35 pb=0.1013; %MPa 大気圧 % e = 17.4; es = [1.5 2 2.5 5 10 15 20 30 40]; pes …
Maximum licensed number of demo users already reached. と言われてAnsys Fluentが起動しない。 他のAnsys licenseのサービスが動いているのが原因。Ansys関連のプログラムをすべて終了する。さらにタスクマネージャーを開き、残っているAnsys関連のプログ…
www.nakka-rocketry.net http://www.nakka-rocketry.net/th_pix/maxthrustderiv.gif 推力係数Cfはpe=paのとき最大になる。 pa=peのときのA*/Aeは上記のequation14で求まる。この膨張比を最適膨張比という。 Nakuja N-2ロケットの海面上での最適膨張比を求め…
「圧縮性流体力学」(松尾) 「ロケット工学基礎講義」(冨田 et al.) clear; close all; clc; k = 1.177; %Propep 3での計算値(KNSB) ストイキ Pc = 1.31; %[MPa] Nakkaのグラフの値。OpenMotorでも同様の値(KNSB)。 rstar = 12; %mm %re = 50; %mm re = 2…
https://www.researchgate.net/post/Pressure_in_ANSYA_CFD-Post_is_it_Relativeどうもゲージ圧力っぽい。