甚是欣喜啊，KC的火箭爱好者们的实力在与日俱增！

要是搞成exe文件该多爽

据说是可以把m脚本编译成exe文件的，方法正在研究中。
PS:在VC中怎样进行作图，哪个函数有这个功能？要自动的，输入x-y数据就能画出来。

openGL美观有爱自刷新双buffer反锯齿

1.1版本可以把输出数据显示在图像上

(110mm KNDX发动机方案初稿)
使用方法：
打开freemat后，再打开burnsim.m，写好参数后按F5运行
%参数名                           单位             含义
solution_name = 'RS-A-5.1'  % 设计方案名称
d_nozzle_begin     = 3      % mm           喷嘴初始直径
nozzle_burn_speed  = 0      % mm/s     喷嘴烧蚀速率
d_hole(1)        =  3              % mm         药柱内孔径
D_outer        = 16                  % mm       药柱外径
N             = 0                             %          药柱能够燃烧的横断面数
L(1)        = 20                         % mm       药柱可燃部分长度（有几段可燃药柱，就填几段的长度之和）
prop         = 'KNDX'                   %          推进剂名称
ap           = 0.1                    % MPa       大气压
% 图片绘制参数，用于适应不同屏幕分辨率，请自己看情况设定，多摸索几次就知道了
line_width = 0.08;%一行文字的高度(0~0.2)
x_pos_1 = 0.03;%第一列的位置(0~1)
x_pos_2 = 0.45;%第二列的位置(0~1)
x_pos_3 = 0.55;%第三列的位置(0~1)
x_pos_4 = 0.8;%第四列的位置(0~1)
fontsize = 10;

运行后，请把图像窗口最大化，以方便观察

burnsim.rar 2.72KB RAR 123次下载

VC中做图请使用GDI，MFC里有也可以不用MFC里的，不过效果不咋的。。。

更新至1.3版，修正了一些参数

burnsim.rar 2.78KB RAR 420次下载

燃料参数在burnsim_propellant.m中修改，举例：
elseif( strcmp(propname,'KNDX'))
% /////////////这里开始是KNDX参数
    if(strcmp(argcname,'vc'))  % m/s      特征速度
           result = 897;  
%//////////////////////// 此处设定特征速度（可以理解为燃料产气产热的能力），受潮的燃料特征速度较低，897是典型值
    elseif(strcmp(argcname,'density'))
        result = 1.879*0.94;      % g/cm^3   密度
%/////////////////////理想密度1.879，实际是1.879的94%
    elseif( strcmp(argcname, 'burnspeed'))  
        if    (pressure < 0.779) %0.1MPa以下按此估算
        a=8.875; n= 0.619;      %          燃速系数
        elseif(pressure < 2.572)
%用if-elseif-else结构描述燃速-压力曲线，也可以直接查表返回数据。
%可以根据燃料性质，自行编写这一部分。
%或者不要改动，使用现有参数
……

XXXXXXXXXXXXXXXXXXXXXXXX/t/30813
更详细的解释

burnsim是一款微型固体火箭发动机内弹道仿真软件
XXXXXXXXXXXXXXXXXXXXXXXX/t/32478
该软件用matlab编写，有极强的扩展性，可以通过修改代码，增加各种功能。本帖就是一个例子。

优化问题分析：

已知参数如下
推进剂：                KNDX
喷嘴喉部初始直径 11mm
喷嘴烧蚀速率         0
药柱内孔径             8mm
药柱外径                45mm
药柱可燃部分长度 400mm
喷管效率                0.85

求解：
药柱能够燃烧的横断面数 N 的最佳值

思路：
考虑从以下【4个指标】评价参数N的合理性
1.燃烧室压力波动，即压力曲线P的标准差std(P)
2.比冲 Isp
3.发动机总冲 I
4.燃烧室峰值压力 max(P)

于是我们编写了以下代码，并命名为burnsim1.3branch1（burnsim版本1.3分支版本1）

burnsim1_3branch1.rar 4.53KB RAR 184次下载

通过对参数N的扫描和反复计算，求解出【指标1,2,3,4】与N的关系，并找出最佳值







综合图像与自动判定结果，得出结论：
N = 10是一个比较恰当的数值，N = (8~12)尚可接受，N<5则是不可接受的，N>13没有必要。

通过修改程序，一样可以对其他参数进行针对性的优化。欢迎大家上传自己的分支版本与仿真报告

赞一个，楼主是总师

下载看看能不能用上。谢谢。

LZ什么时候加入BP的支持啊？

求封装，用着巨麻烦

想问下软件如何使用

下载看看能不能用上。谢谢。

唔，楼主是怎么处理药柱末端和喷口相分离的地方的，求指教