2022-03-01から1ヶ月間の記事一覧

Edgeでファイルダウンロードのブロック(SmartScreen)を無効化する

www.cnet.com

OpenRocketのテストコード

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…

遅延を含むシステムのモデリング@Simulink

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); 以下のステップ応答をもつ伝達関数をモデリングした。 実際はこれは移動平均フィル…

OpenRocketプラグインのテンプレート

github.comgithub.com

OpenRocketのチュートリアル

www.youtube.com

OpenRocketのプラグイン開発

26:42あたりで一瞬触れられている。 youtu.be

OpenRocketをビルド

1. eclipseをインストール 2. JDK 11をインストール(Java17では動かない) 3. OpenRocketをクローン 4. eclipseでOpenRocketフォルダをインポート 5. エラーが出るので修正する。自分の場合はjson-xxx.jarとjavax.json-xxx.jarをダウンロードして読み込んだ…

Matplotlibの日本語化けを解消する

datumstudio.jp以下でインストール pip install japanize-matplotlib使い方は単にインポートするだけ import japanize_matplotlib

伝熱工学・数値流体力学の講義資料

筑波大学 熱流体制御・計測研究室 www.kz.tsukuba.ac.jp

スミスチャートによるマッチング回路の設計

www.youtube.com設計ツール Quicksmith quicksmith.online

境界層運動量積分方程式の導出

燃焼工学のe-leanringサイト

https://elearning.cerfacs.fr/combustion/index.php

MacのPreview.appのSignatureを他のMacに移す

discussions.apple.com

Ansys Fluentのチュートリアル資料

公式が作成しているらしく、かなり出来が良い。 imechanica.org

いろいろな非定常熱伝導問題をMatlabで解く

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…

Matlabで平板の非定常熱伝導を解く

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…

Matlabで偏微分方程式を解く(非定常熱伝導)

以下を参考に。 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#:…

Ansys Fluentで円柱周りの抗力係数を計算する

Fluentでの操作のポイントは、Reference valueを設定すること。 以下のページに様々な形状の抗力係数が示されている。 skomo.o.oo7.jp 計算 縦に置かれた円柱、l/d=1の場合を試してみる。うまくいけばCd=0.63になるはず。 条件 円柱のd=10mmとする。奥行きL…

レイリー流れとレイリー問題

似ているようで異なる。 レイリー流れ ・断面積が一定の管内 ・定常流れ ・摩擦なし ・熱の授受はあり の流れをレイリー流れ(Rayleigh flow)と呼ぶ。 摩擦がないので可逆流れである。 ※垂直衝撃波はレイリー流れとファノ流れの両方の条件を満たす レイリー…

Ansys fluentで使える粘性モデル

Ansys Fluentで揚力(Lift)と抵抗(Drag)の計算

www.youtube.com www.youtube.com

Ansys Fluentで過渡応答のシミュレーション

Transientを選択するだけではアニメーションが保存されない。 Solution>Calculation activitiesでアニメーションの保存タイムステップを定義する必要がある。 youtu.beyoutu.beyoutu.beyoutu.be

推力係数とノズル開口比eの関係

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 …

Ansys Fluentが起動しない

Maximum licensed number of demo users already reached. と言われてAnsys Fluentが起動しない。 他のAnsys licenseのサービスが動いているのが原因。Ansys関連のプログラムをすべて終了する。さらにタスクマネージャーを開き、残っているAnsys関連のプログ…

推力係数を最大するノズル開口比e(ノズル出口圧=背圧のとき→適正膨張)

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…

Ansys FluentのCFD PostのPressureはゲージ圧力なのか?

https://www.researchgate.net/post/Pressure_in_ANSYA_CFD-Post_is_it_Relativeどうもゲージ圧力っぽい。