内弹道仿真器 Burnsim 【更新至1.3版】
warmonkey2011/03/06喷气推进 IP:广东
翻箱底的东西

已经经过试车检验,可认为计算较准确
RS-A-3.1-test1-sim.jpg
仿真

实测: https://www.kechuang.org/t/19505
PS:工作时间是1.1s,此帖中的估计值1s太过于粗略,具体自己看录音波形图

设计图: https://www.kechuang.org/t/19546

上一测试版本
https://www.kechuang.org/t/30813

目前能输出推力-时间图,压力-时间图,燃速-时间图,药柱质量,比冲,总冲
【此为1.0版本,1.3版本看7楼】
程序:  【直接跳到7楼下载】https://www.kechuang.org/t/32478
matlab运行环境,使用freemat代替: XXXXXXXXXXXXXXXXXXXXXt/projects/freemat/files/
文件大小17M
使用方法:
打开freemat后,再打开burnsim.m,写好参数后按F5运行
%参数名                           单位             含义
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_propellant.m中(如果你只用KNDX/KNSU/KNSB,则不需要修改)

参考文献:  https://www.kechuang.org/t/12656

(此为旧版本,请在7楼下载 1.3版 )1.0版本
attachment icon burnsim.rar 2.20KB RAR 118次下载

[修改于 8年10个月前 - 2015/07/16 02:11:13]

+2500  科创币    xj198398xujing    2011/03/07 非常不错  鼓励
来自:航空航天 / 喷气推进
17
 
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
lilith
13年3个月前 IP:未同步
284447
甚是欣喜啊,KC的火箭爱好者们的实力在与日俱增!
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
ehco
13年3个月前 IP:未同步
284531
要是搞成exe文件该多爽
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
warmonkey作者
13年3个月前 IP:未同步
284534
据说是可以把m脚本编译成exe文件的,方法正在研究中。
PS:在VC中怎样进行作图,哪个函数有这个功能?要自动的,输入x-y数据就能画出来。
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
novakon
13年3个月前 IP:未同步
284909
回 3楼(warmonkey) 的帖子
openGL美观有爱自刷新双buffer反锯齿
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
warmonkey作者
13年3个月前 修改于 8年10个月前 IP:未同步
285985
1.1版本可以把输出数据显示在图像上
110mmKNDX发动机方案初稿-RS-A-5.1beta1.jpg
(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;

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

attachment icon burnsim.rar 2.72KB RAR 100次下载
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
phpskycn
13年2个月前 IP:未同步
287975
VC中做图请使用GDI,MFC里有也可以不用MFC里的,不过效果不咋的。。。
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
warmonkey作者
12年8个月前 修改于 8年10个月前 IP:未同步
329322
更新至1.3版,修正了一些参数

attachment icon burnsim.rar 2.78KB RAR 372次下载

燃料参数在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结构描述燃速-压力曲线,也可以直接查表返回数据。
%可以根据燃料性质,自行编写这一部分。
%或者不要改动,使用现有参数
……

https://www.kechuang.org/t/30813
更详细的解释
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
warmonkey作者
12年8个月前 IP:未同步
329451
发动机参数自动优化-burnsim功能扩展示例
burnsim是一款微型固体火箭发动机内弹道仿真软件
https://www.kechuang.org/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)

attachment icon burnsim1_3branch1.rar 4.53KB RAR 125次下载

通过对参数N的扫描和反复计算,求解出【指标1,2,3,4】与N的关系,并找出最佳值
Simulate.jpg
P_Vs_N.jpg
std_P_Vs_N.jpg
I_Vs_N.jpg
Isp_Vs_N.jpg
Pmax_Vs_N.jpg

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

通过修改程序,一样可以对其他参数进行针对性的优化。欢迎大家上传自己的分支版本与仿真报告
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
jerish
12年8个月前 IP:未同步
329702
赞一个,楼主是总师
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
8-8-
12年7个月前 IP:未同步
333843
下载看看能不能用上。谢谢。
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
火星居士
12年7个月前 IP:未同步
334014
回 楼主(warmonkey) 的帖子
LZ什么时候加入BP的支持啊?
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
guohc
12年7个月前 IP:未同步
334275
求封装,用着巨麻烦
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
qqwoshishei
11年5个月前 IP:未同步
482991
想问下软件如何使用
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
冯智琰
10年10个月前 IP:未同步
548290
下载看看能不能用上。谢谢。
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
海澄
8年4个月前 IP:江苏
805988
唔,楼主是怎么处理药柱末端和喷口相分离的地方的,求指教
引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
154454496
5年3个月前 IP:河北
853682

话说这个能算RNX吗

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

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

所属专业
所属分类
上级专业
同级专业
warmonkey
学者 机友
文章
357
回复
7648
学术分
14
2008/10/11注册,16分54秒前活动

Cubesat

主体类型:个人
所属领域:无
认证方式:手机号
IP归属地:未同步
文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}