不很懂程序 只有代码很无爱 能封装一下吗
30813
%7B%22isLastPage%22%3Atrue%2C%22notes%22%3A%5B%5D%2C%22pid%22%3A%22t30813%22%2C%22tid%22%3A%2230813%22%2C%22mainForumsId%22%3A%5B%2289%22%5D%2C%22categoriesId%22%3A%5B%220%22%5D%2C%22tcId%22%3A%5B%5D%7D
%7B%22isEditMode%22%3Afalse%7D
微型固体火箭发动机内弹道仿真软件
感谢Cim的帮助更正了一个参数，目前计算已经正常
可以使用freemat代替matlab，安装包只有17M
下载地址： XXXXXXXXXXXXXXXXXXXXXt/projects/freemat/files/
自己找最新版的安装（发帖时最新版是Freemat4）
打开freemat后，再打开burnsim.m，写好参数后按F5运行，耐心等待结果，1秒的燃烧大概要5秒来仿真
暂时只能输出P-t图，燃料可以自定义无限种
药柱形状只要你自己能搞出面积公式和体积公式，就能随便改
目前是按照端面+内孔燃烧配置。
使用方法：
在burnsim.m中配置以下参数：
%参数名 单位 含义
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 大气压
burnsim_propellant.m中设置推进剂密度
然后运行burnsim.m
推进剂自己放在burnsim_propellant.m中
运行结果
![p-t.jpg]()
代码：
没装matlab的试试freemat
附画出燃速-压强图的代码：
%画出燃速-压强图
%绘图起点 点距 终点
p_start = 0.1
p_step = 0.01
p_end = 12
prop = 'KNDX' %推进剂名称
x=(p_start:p_step:p_end);
for i=1:length(x);y(i) = burnsim_propellant(x(i),prop,'burnspeed');end
plot(x,y);
xlabel P(MPa)
ylabel r(mm/s)
grid on
计算燃料质量：
V_gram(1) * de
运行完后输入
2楼补充如何设置燃料
可以使用freemat代替matlab，安装包只有17M
下载地址： XXXXXXXXXXXXXXXXXXXXXt/projects/freemat/files/
自己找最新版的安装（发帖时最新版是Freemat4）
打开freemat后，再打开burnsim.m，写好参数后按F5运行，耐心等待结果，1秒的燃烧大概要5秒来仿真
暂时只能输出P-t图，燃料可以自定义无限种
药柱形状只要你自己能搞出面积公式和体积公式，就能随便改
目前是按照端面+内孔燃烧配置。
使用方法：
在burnsim.m中配置以下参数：
%参数名 单位 含义
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 大气压
burnsim_propellant.m中设置推进剂密度
然后运行burnsim.m
推进剂自己放在burnsim_propellant.m中
运行结果
代码：
没装matlab的试试freemat
附画出燃速-压强图的代码：
%画出燃速-压强图
%绘图起点 点距 终点
p_start = 0.1
p_step = 0.01
p_end = 12
prop = 'KNDX' %推进剂名称
x=(p_start:p_step:p_end);
for i=1:length(x);y(i) = burnsim_propellant(x(i),prop,'burnspeed');end
plot(x,y);
xlabel P(MPa)
ylabel r(mm/s)
grid on
计算燃料质量：
V_gram(1) * de
运行完后输入
2楼补充如何设置燃料
在burnsim_propellant.m倒数第二个end上面添加
位置:
function result = burnsim_propellant( pressure , propname , argcname )
if( strcmp(propname,'KNSU'))
if(strcmp(argcname,'vc'))
result = 923; % m/s 特征速度
elseif(strcmp(argcname,'density'))
result = 1.81; % g/cm^3 密度
elseif( strcmp(argcname, 'burnspeed'))
a=8.263; n=0.319; % 燃速系数
result = a * pressure ^ n;
%燃速 r = a * pressure ^ n% 燃速指数
else
end
%此处添加新燃料
end
end
要添加的代码
elseif( strcmp(propname,'你的燃料名称'))
if(strcmp(argcname,'vc')) % m/s 特征速度
result = 特征速度 ; %这里填入特征速度
elseif(strcmp(argcname,'density'))
result = 密度; % g/cm^3 密度
elseif( strcmp(argcname, 'burnspeed')) %各个压力范围和燃速数据
if (pressure < 0.779) %编程语法类似C,自己填入参数
a=你的数据; n= 0.619; % 燃速系数
elseif(pressure < 2.572)
a=7.553; n=-0.009; % 燃速系数
elseif(pressure < 5.930)
a=3.841; n=0.688; % 燃速系数
elseif(pressure < 8.502)
a=17.20; n=-0.148; % 燃速系数
else
a=4.775; n=0.442; % 燃速系数
end%要更多个压力范围，就添加更多的elseif 但是必须保证最后是else。
result = a * pressure ^ n;%燃速计算方法，压力是MPa燃速是mm/s
else
end
位置:
function result = burnsim_propellant( pressure , propname , argcname )
if( strcmp(propname,'KNSU'))
if(strcmp(argcname,'vc'))
result = 923; % m/s 特征速度
elseif(strcmp(argcname,'density'))
result = 1.81; % g/cm^3 密度
elseif( strcmp(argcname, 'burnspeed'))
a=8.263; n=0.319; % 燃速系数
result = a * pressure ^ n;
%燃速 r = a * pressure ^ n% 燃速指数
else
end
%此处添加新燃料
end
end
要添加的代码
elseif( strcmp(propname,'你的燃料名称'))
if(strcmp(argcname,'vc')) % m/s 特征速度
result = 特征速度 ; %这里填入特征速度
elseif(strcmp(argcname,'density'))
result = 密度; % g/cm^3 密度
elseif( strcmp(argcname, 'burnspeed')) %各个压力范围和燃速数据
if (pressure < 0.779) %编程语法类似C,自己填入参数
a=你的数据; n= 0.619; % 燃速系数
elseif(pressure < 2.572)
a=7.553; n=-0.009; % 燃速系数
elseif(pressure < 5.930)
a=3.841; n=0.688; % 燃速系数
elseif(pressure < 8.502)
a=17.20; n=-0.148; % 燃速系数
else
a=4.775; n=0.442; % 燃速系数
end%要更多个压力范围，就添加更多的elseif 但是必须保证最后是else。
result = a * pressure ^ n;%燃速计算方法，压力是MPa燃速是mm/s
else
end
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。