在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