现在还是感觉Nakka的excel用得爽些

我以前一直用的是 那个叫做 cim 那个人做的一个 计算软件，感觉还是很好用的

不之道你的这个怎么样。

引用第3楼爱冉静的陆飞于2010-12-19 20:12发表的  :
我以前一直用的是 那个叫做 cim 那个人做的一个 计算软件，感觉还是很好用的

不之道你的这个怎么样。

你看看lz的ID。。。就知道谁是cim了

我写了个matlab版本的，支持多段药柱，正好和你的核对一下结果

LS结果如何？

%BurnSim V0.0.1
%固体火箭发动机燃烧过程仿真器

clear all;
%参数名                     单位      含义
d_nozzle_begin     = 6      % mm       喷嘴初始直径
nozzle_burn_speed  = 0      % mm/s     喷嘴烧蚀速率
d_hole(1)        =  6       % mm       药柱内孔径
D_outer        = 24             % mm       药柱外径
N             = 2              %          药柱能够燃烧的横断面数
L(1)        = 70               % mm       药柱可燃部分长度
prop         = 'KNSU'               %          推进剂名称

%初始值                                             单位
pi = 3.1415926;                                 %圆周率  
d_nozzle(1) = d_nozzle_begin;                   %   mm
r(1) = burnsim_propellant( 0.1 ,prop,'burnspeed');  %   mm/s
V_gram(1) =(D_outer^2-d_hole(1)^2)*pi/4*L(1)/1000; %   cm^3
c = burnsim_propellant(0,prop,'vc');            %   m/s
de = burnsim_propellant(0,prop,'density');      %   g/cm^3
t=0;                                            % s
dt          = 0.001                % s        分析时间片长度,数值小精度高,数值大需要内存少,修改后注意要更改图像的X坐标名称

i = 1;%循环计数
while( V_gram(i) > 0 )
%计算公式                                                       得数含义              单位
Ae(i) = d_nozzle(i)^2 * pi / 4;                                %  当前喷管直径        mm^2
S_burn(i) =( D_outer^2 - d_hole(i)^2 )*pi/4*N + d_hole(i)*pi*L(i); %  当前燃烧面积        mm^2
Kn(i) = S_burn(i) / Ae(i) ;                                    %  喷燃比
P(i) = Kn(i) * r(i) * c * de                                  %  压力                MPa
r(i+1) = burnsim_propellant(P(i),prop,'burnspeed');            %  当前燃速            mm/s
d_hole(i+1) = d_hole(i) + 2*r(i)*dt;                           %  下一时刻的药柱内孔径 mm
L(i+1) = L(i) - N * dt * r(i);                                 %  下一时刻的药柱长度  mm
V_gram(i+1) = (D_outer^2 - d_hole(i)^2)* pi /4 *L(i+1) /1000;    %  下一时刻的药柱体积  cm^3
d_nozzle(i+1) = d_nozzle(i) + 2*nozzle_burn_speed*dt;           %下一时刻的喷管内径 mm

    disp('T = ');      disp(t);
    disp('Vgram = ');  disp(V_gram(i));
    disp('[s:9] = ');      disp(P(i));
    disp('Kn = ');     disp(Kn(i));
    
i = i+1;
t = t + dt;
    
end

公式的单位有问题？我计算出来压力比自己估计的大10倍左右。

感谢Cim的帮助，原来是KNSU参数设置错误了，现在结果已经正确，发在另外一个帖子里面了，其实和LZ的差不多，不过是基于matlab的
LZ的算法比较细致，燃烧室吸热都想到了～NB。
话说LZ计算推力没有考虑排气的动量，下一版能否加上？

必须赞扬[s:251]

赞一个！！！！！！！

相当不错 赞扬
希望能加上添加其他燃料和设置燃料相关参数的选相
不懂编程来膜拜LZ了  呵呵

感谢LS捧场
第二版正在制作中主要改动如下：
1.修正KN系列燃料数据
2.添加端面燃烧情况的计算
3.添加多段药柱情况的计算
4.添加其他燃料自定义（需要燃料的参数为： 燃速系数 压强指数 特征速度 燃料密度）
5.增加动量对推力影响的讨论从而提高计算精度

预期完善计划
1.补充压强下降段计算
2.添加喷管收扩比对推力影响
3.能够输出 F－t   P－t  数据

我下载了怎么用不了那
初始化失败啊

能不能尝试加入不同的燃烧类型的计算功能，譬如多段燃料，端面燃烧等等！
很是期待第二版啊，也可以试着有一个参数输入的功能，这样对适用面就更广了！

最新版本地址
https://www.kechuang.org/t/30978

引用 usercim:
最新版本地址
https://www.kechuang.org/t/30978

非常感谢

为啥下载不了？