微型固体火箭发动机内弹道仿真软件
感谢Cim的帮助更正了一个参数,目前计算已经正常

可以使用freemat代替matlab,安装包只有17M
下载地址: http://sourceforge.net/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

代码:

attachment icon burnsim.rar 2.00KB RAR 65次下载

没装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楼补充如何设置燃料
+500  科创币    虎哥   2010-12-20   还是鼓励一下。
来自:火箭技术 / 喷气推进
 
xj198398xujing
10年4个月前
1楼
不很懂程序  只有代码很无爱 能封装一下吗
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
warmonkey作者
10年4个月前
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
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
YID2014
6年6个月前
3楼
楼主威武!
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

想参与大家的讨论?现在就 登录 或者 注册

所属专业
上级专业
同级专业
warmonkey
学者 机友 笔友
文章
328
回复
6648
学术分
16
2008/10/11注册,7 小时前活动

科创航天主席

%7B%22isDisplay%22%3Atrue%7D

仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
等待中...
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
处理中..
处理失败
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传
空空如也~
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
{{forum.displayName}}
{{forum.countThreads}}
篇文章,
{{forum.countPosts}}
条回复
{{forum.description || "暂无简介"}}
ID: {{user.uid}}
学术分隐藏
{{submitted?"":"投诉或举报"}}
请选择违规类型:
{{reason.description}}
支持的图片格式:jpg, jpeg, png
插入公式
分享回复:{{shareId}}
加载中...
评论控制
加载中...
文号:{{pid}}
加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
加入关注取消关注
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
建议修改
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}
下载资料
{{fileName}}
大小:{{size}}
下载当前附件将花费 {{costMessage}}
你当前剩余 {{holdMessage}}
{{fileName}}
大小:{{size}}
当前附件免费。
你已购买过此附件,下载当前附件不需要花费积分。
加载中...
{{errorInfo}}
附件已丢失
当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}