基于状态模拟方法分析高压油管压力控制问题---源代码（matlab）

最新推荐文章于 2021-08-30 01:01:54 发布
最新推荐文章于 2021-08-30 01:01:54 发布
阅读量5.5k 收藏 64
点赞数 15
分类专栏： matlab 数学建模 文章标签： 数学建模
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zhulinhao/article/details/107238842
版权

应一些读者要求，这里加上了源代码，奈何本人最近忙于数学建模培训，所以排版可能不太友好，但是里面每一个都写了序号，你们可以看一下

1.m1.m
clear;clc;
load data1
y=E;
x=[ones(size(y)) p.^2 p];
[b,bint,r,rint,stats]=regress(y,x)
P=b(1)+b(2)*(p.^2)+b(3)*p;
figure;
plot(p,E,'r--',p,P,'b');
title("tan xing mo liang he ya qiang");
xlabel("ya qiang");
ylabel("tan xing mo liang");
legend('shi ji zhi','hui gui qu xian');

2.m2.m
clear;clc;
rho=0.85:0.005:0.90;
P=-53.23+227.0*tan(1.659+6.558*log(rho));
figure;
plot(rho,P);
title('???????');
xlabel('??(mg/mm^3)');
ylabel('??(MPa)');

3.m3.m
clear;clc;
load data1;
x=100:0.1:150;
y=1./(b(2).*x.^2+b(3).*x+b(1));
rho0=0.85;
rho=exp(trapz(x,y)+log(rho0))
x=100:0.1:160;
y=1./(b(2).*x.^2+b(3).*x+b(1));
rho0=0.85;
rho=exp(trapz(x,y)+log(rho0))
4.m4.m
clear;clc;
options=struct();
options.yfun=@yfun_1;
options.target=@target_1;
options.isplot=0;
options.interval=0.1;
options.duration=2000;
[fm,mx]=stepsolve(@targetfun,0.26,0.31,0.0005,0.02,5,options);
options.isplot=1;
targetfun(mx,options);
ylabel('??(MPa)');
xlabel('??(ms)');
title('???????????');
5.m5.m
clear;clc;
options=struct();
options.yfun=@yfun_1;
options.target=@target_2;
options.isplot=0;
options.interval=1;
options.duration=10000;
[fm,mx]=stepsolve(@targetfun,0.6,1,0.0005,0.05,5,options);
options.isplot=1;
targetfun(mx,options);
ylabel('??(MPa)');
xlabel('??(ms)');
title('???????????');
6.m6.m
clear;clc;
options=struct();
options.yfun=@yfun_2;
options.target=@target_1;
options.isplot=0;
options.interval=1;
options.duration=10000;
[fm,mx]=stepsolve(@targetfun,0.019,0.029,0.000001,0.001,3,options);
options.isplot=1;
targetfun(mx,options);
ylabel('??(MPa)');
xlabel('??(ms)');
title('???????????');
7.m7.m
clear;clc;
options=struct();
options.yfun=@yfun_3;
options.target=@target_1;
options.isplot=0;
options.interval=1;
options.duration=10000;
[fm,mx]=stepsolve(@targetfun,0.039,0.059,0.00001,0.002,3,options);
options.isplot=1;
targetfun(mx,options);
ylabel('高压油管压力(MPa)');
xlabel('运行时间(ms)');
title('两个喷油口工作时高压油管内部压力');
8.m8.m
clear;clc;
options=struct();
options.yfun=@yfun_4;
options.target=@target_1;
options.isplot=0;
options.interval=1;
options.duration=10000;
[fm,mx]=stepsolve(@targetfun,97,103,0.1,1,3,options);
options.isplot=1;
targetfun(mx,options);
ylabel('高压油管压力(MPa)');
xlabel('运行时间(ms)');
title('增加减压阀后高压油管内部压力');
9.Pfun.m
function P=Pfun(x)
    if (x<0.7)
        P=0;
    else
        P=(-53.23+227.0*tan(1.659+6.558*log(x)));
    end
end
10.Q_in_1.m
function Q=Q_in_1(tA,t,P_in)
%tA????
%t????
%P_in
T=tA+10;
if mod(t,T)<tA
    Q=0.85*(pi/4)*(1.4)^2*sqrt(2*(160-P_in)/0.871);
else
    Q=0;
end
end
11.Q_in_2.m
function Q=Q_in_2(w,t,P_in)
%tA????
%t????
%P_in
T=tA+10;
if mod(t,T)<tA
    Q=0.85*(pi/4)*(1.4)^2*sqrt(2*(160-P_in)/0.871);
else
    Q=0;
end
end
12.Q_out_1.m
function Q=Q_out_1(t,P)
r=mod(t,100);
if r<=0.2
    Q=100*r;
elseif r<=2.2
    Q=20;
elseif r<=2.4
    Q=-100*r+240;
else 
    Q=0;
end
end
13.Q_out_2.m
function Q=Q_out_2(t,p)
%t:时间
%p:密度
r=mod(t,100);
b1=[16.3667705292548,-2.27403710032532,0.0443842537666548];
b2=[16.4024620318978,-78.0705445458018,92.8628884740531];
if r<=0.44
    h=polyval(b1,r);
elseif r<=2
    h=2;
elseif r<=2.45
    h=polyval(b2,r);
else 
    h=0;
end
dk=1.4;
dz=2.5;
alpha=9/180*pi;
A=min([pi*dk^2/4,pi/4*(4*dz*h*tan(alpha)+2*h^2*tan(alpha)^2)]);
Q=0.85*A*(2*Pfun(p)/p)^0.5;
end
14.Q_out_3.m
function Q= Q_out_3(p,rho)
 Q=0.85*0.7^2*pi*sqrt(p/rho);
end
15.stepsolve.m
function [mf,M]=stepsolve(fun,L,R,tol,h0,d,options)
    if nargin<=6
        options=struct();
    end
    optfun=@(x)fun(x,options);
    M=inf;
    fprintf("Initiating...\n");
    mf=optfun(M);
    while R-L>tol
        fprintf("Calculating.. From %.5f to %.5f step %.5f\n",L,R,h0);
        t_A=L:h0:R;
        ff=arrayfun(optfun,t_A);
        mf=min(ff);
        M=t_A(mf==ff);
        L=M-(M-L)/d;
        R=M+(R-M)/d;
        h0=h0/d;
        fprintf("Current: Min: %.5f At %.5f\n",mf,M);
    end 
end
16.target_1.m
function [f]=target_1(TT,PP)
    f=norm(PP-100,2)/(max(TT)-min(TT));
end
17.target_2.m
function [f]=target_2(TT,PP)
    if max(PP)>150
        f=sum(diff(TT).*PP(1:end-1));
    else
        f=1e12;
    end
end
18.targetfun.m
function f=targetfun(w,options)
    if nargin==1
        options=struct();
    end
    interval=initfun(options,'interval',1);
    duration=initfun(options,'duration',20000);
    isplot=initfun(options,'isplot',1);
    target=options.target;
    yfunopt=options.yfun;
    m=0;
    global flag;
    flag=0;
    yfun=@(t,y,m)yfunopt(t,y,m,w,interval);
    y=0.850;
    steps=0:interval:duration;
    TT=zeros(size(steps));
    YY=zeros(size(steps));
    count=0;
    for t=steps
        count=count+1;
        [y,m]=yfun(t,y,m);
        TT(count)=t;
        YY(count)=y;
    end
    PP=Pfun(YY);
    f=target(TT,PP);
    if isplot
        figure;
        plot(TT,PP,'r');
        xlim([min(TT),max(TT)]);
        ylim([min(PP),max(PP)]);
    end
end
19.V_in_2.m
function V=V_in_2(w,t)
    A=2.5^2*pi;
    V0=20+2.413*A*2;
    rho=-2.413*cos(w*t)+2.413;
    V=V0-rho*A;
end
20.yfun_1.m
function [y,m]=yfun_1(t,y,m,w,dt)
    Q_in=@Q_in_1;
    Q_out=@Q_out_1;
    V_0=500*(pi/4)*(10)^2;
    y=y+((1/V_0)*(0.87*Q_in(w,t,Pfun(y))-y*Q_out(t,Pfun(y))))*dt;
end
21.yfun_2.m
function [y,m]=yfun_2(t,y,m,w,dt)
    if (mod(t,2*pi/w)<dt)
        m=V_in_2(w,0)*0.850;
    end
    V_0=500*(pi/4)*(10)^2;
    V_A=V_in_2(w,t);
    rho_A=m/V_A;
    m_out=Q_out_2(t,y)*y;
    P_A=Pfun(rho_A);
    P_N=Pfun(y);
    if (P_A>P_N)
        Q_A=0.85*0.7*0.7*3.1415*sqrt(2*(P_A-P_N)/rho_A);
        y=y+(1/V_0)*(rho_A*Q_A-m_out)*dt;
        m=m-rho_A*Q_A*dt;
        if (m<0)
            m=0;
        end
    else
        if (rho_A==0)
            pause;
        end
        y=y+(1/V_0)*(-m_out)*dt;
    end
end
22.yfun_3.m
function [y,m]=yfun_3(t,y,m,w,dt)
    if (mod(t,2*pi/w)<dt)
        m=V_in_2(w,0)*0.850;
    end
    V_0=500*(pi/4)*(10)^2;
    V_A=V_in_2(w,t);
    rho_A=m/V_A;
    m_out=Q_out_2(t,y)*y*2;
    P_A=Pfun(rho_A);
    P_N=Pfun(y);
    if (P_A>P_N)
        Q_A=0.85*0.7*0.7*3.1415*sqrt(2*(P_A-P_N)/rho_A);
        y=y+(1/V_0)*(rho_A*Q_A-m_out)*dt;
        m=m-rho_A*Q_A*dt;
        if (m<0)
            m=0;
        end
    else
        if (rho_A==0)
            pause;
        end
        y=y+(1/V_0)*(-m_out)*dt;
    end
end
23.yfun_4.m
function [y,m]=yfun_4(t,y,m,P_L,dt)
    global flag;
    w=0.05045;
    if (mod(t,2*pi/w)<dt)
        m=V_in_2(w,0)*0.850;
    end
    V_0=500*(pi/4)*(10)^2;%高压油管容积
    V_A=V_in_2(w,t);%高压油泵内体积
    rho_A=m/V_A;%高压油泵内密度
    m_out_1=0;%喷油嘴堵塞
    m_out_2=Q_out_3(Pfun(y),y)*y;%减压阀出油质量
    m_out_3=Q_out_2(t,y)*y*2;%喷油嘴修复后出油质量
    P_A=Pfun(rho_A);%高压油泵内压力
    P_N=Pfun(y);%高压油管内压力
    P_H=180;
    if(P_N>P_H)
        flag=1;
    end
    if (flag==0)%喷油嘴堵塞
        if (P_A>P_N)
            Q_A=0.85*0.7*0.7*pi*sqrt(2*(P_A-P_N)/rho_A);
            P_A=Pfun(rho_A);
            y=y+(1/V_0)*(rho_A*Q_A-m_out_1)*dt;
            P_N=Pfun(y);
            m=m-rho_A*Q_A*dt;
            if (m<0)
                m=0;
            end
        else
            y=y+(1/V_0)*(-m_out_1)*dt;
        end
    end
    %喷油嘴和减压阀同时打开
    if (flag==1&&P_N>=P_L)
        if (P_A>P_N)
            Q_A=0.85*0.7*0.7*pi*sqrt(2*(P_A-P_N)/rho_A);
            y=y+(1/V_0)*(rho_A*Q_A-m_out_2-m_out_3)*dt;
            m=m-rho_A*Q_A*dt;
            if (m<0)
                m=0;
            end
        else
            y=y+(1/V_0)*(-m_out_2-m_out_3)*dt;
        end
    end
    %减压阀关闭    
    if (P_N<P_L)
        if (P_A>P_N)
            Q_A=0.85*0.7*0.7*pi*sqrt(2*(P_A-P_N)/rho_A);
            y=y+(1/V_0)*(rho_A*Q_A-m_out_3)*dt;
            m=m-rho_A*Q_A*dt;
            if (m<0)
                m=0;
            end
        else
            y=y+(1/V_0)*(-m_out_3)*dt;
        end
    end
end
NUAA&XMU---朱林昊
关注
  • 15
    点赞
  • 踩
  • 64
    收藏
    觉得还不错? 一键收藏
  • 24
    评论
专栏目录
2019国赛A题高压油管的压力控制 Python代码_2019高压油管建模_2019国赛a_python_
10-02
2019国赛A题高压油管的压力控制 Python代码
[独有源码]java-jsp烯烃厂压力管道管理平台n6x6v从不会做毕业设计到成功完成的过程与方法
最新发布
海灵学姐
07-06 53
选题背景： 烯烃厂作为化工行业的重要组成部分，其生产过程中涉及到大量的压力管道。这些管道承载着高温、高压的流体，一旦发生泄漏或故障，可能会导致严重的事故和损失。因此，对于烯烃厂压力管道的管理至关重要。传统的管道管理方式存在着许多问题，如信息不透明、数据收集困难、维护保养不及时等，亟需建立一个科学、高效的管道管理平台。 选题意义： 建立烯烃厂压力管道管理平台具有重要的意义。首先，通过该平台可以实现对压力管道的全面监控和管理，及时发现潜在的问题并采取相应的措施，从而降低事故发生的风险。其次，该平台可以实现对管
基于状态模拟方法分析高压油管压力控制问题
zhulinhao的博客
07-06 8798
基于状态模拟方法分析高压油管压力控制问题 摘要 燃油进入和喷出是许多燃油发动机工作的基础，但燃油进入和喷出的间歇性工作过程会导致高压油管内压力的变化，使得所喷出的燃油量出现偏差，从而影响发动机的工作效率。本文主要研究的是高压油管的压力控制问题,核心是高压油管内部压力随时间的变化。由注1得油管内部压力与燃油密度的关系，进而将问题转化为高压油管内燃油密度随时间的变化。建立以油管内燃油密度与时间的微分方程并求出燃油密度与时间的关系，进一步得到油管内部压力与燃油密度的关系，从而得到油管内部压力与时间的关系。 针对问
2019高教社杯全国大学生数学建模竞赛题目A题 高压油管的压力控制
热门推荐
亲，如果文章对你有帮助的话请点个关注吧！
09-12 6万+
和往年一样，A题还是如此。。。。 A题 高压油管的压力控制 燃油进入和喷出高压油管是许多燃油发动机工作的基础，图1给出了某高压燃油系统的工作原理，燃油经过高压油泵从A处进入高压油管，再由喷口B喷出。燃油进入和喷出的间歇性工作过程会导致高压油管内压力的变化，使得所喷出的燃油量出现偏差，从而影响发动机的工作效率。 图1 高压油管示意图 问题1. 某型号高压油管的内腔长度为500m...
【数学建模】2019A题 高压油管的压力控制器（解题思路）
XiaoChen~Coder
09-19 1万+
最近收到很多小伙伴的私信，说要A题的解题思路。最近忙于考证的事情，没有太多时间来整理，以下整理思路如下： 问题分析： emmm...要去上课了，暂时先写这么多....... ...
wzzzzzw_基于matlab的高压油管压力分析_油管压力_
09-28
在本项目中，我们主要探讨的是“基于MATLAB的高压油管压力分析”。MATLAB，全称Matrix Laboratory，是一款强大的数学计算和数据可视化软件，广泛应用于工程、科学和经济领域。在高压油管的压力分析中，MATLAB可以...
（A023）基于数值模拟的高压油管压力控制模型1
08-03
总的来说，本研究通过数值模拟方法对高压油管压力控制模型进行了深入探究，不仅解决了如何在不同工况下维持油管压力稳定的问题，还为发动机性能优化提供了理论支持。这种方法论对于燃油发动机的设计和控制具有很高的...
2019年A题 数学建模 高压油管压力控制-1.pdf
08-18
该问题的核心是理解和控制高压油管内压力的变化，以确保系统的稳定性和效率。 在数学建模过程中，首先对高压油管内压强与燃油密度随时间的变化规律进行了深入研究。通过对题目中给出的压力变化量与密度变化量关系的...
论文研究 - 高压油管压力控制策略调整
05-20
基于此问题给出的数据进行相应的分析，根据不同的燃料进入和运行条件，计算高压油管中的压力变化。喷出，以确定燃油喷射系统的一些运行参数，从而提高发动机效率和经济效率。 通过在相应条件下的质量守恒公式构造...
论文研究 - 高压油管的压力控制
05-20
本文重点介绍如何控制高压管中的压力。 首先，对系统进行整体分析，以确定高压燃油管的压力平衡条件； 也就是说，流出高压燃油管的燃油量等于流入的燃油量。本文使用Excel在所有数据处理中拟合和组织数据。 在求解...
2019年全国大学生数学建模高压油管的压力控制
05-16
本文档由三人参赛时编写，荣获全国大学生数学建模省级三等奖，此次上传文档，供大家参考学习。采用的语言是MatLab，用于数据分析，绘图。
2019数学建模A题高压油管
09-21
三个问题的代码，分别把problem1，problem2，problem3文件夹复制到matlab存放代码的文件夹里面，然后分别运行main函数就可以得到每一个问题的结果了，注释里面有我调试的答案
2019A高压油管的压力控制.zip
01-03
2019A高压油管的压力控制.zip
基于智能仪表的气罐压力控制系统设计
05-18
　　　所谓压力控制系统就是利用管道或容器中的介质压力作为被控制量，从而保证输出一个恒定的气压的反馈控制系统。目前生产中应用的压力控制系统，主要以传统的PID控制算法为主。但对于复杂的大型系统，其数学模型往往难以获得，传统的PID控制方式显得无能为力。为适应复杂控制系统的控制要求，人们研究了很多智能控制方法，模糊PID控制便是其中之一。本文主要研究了模糊PID控制及其改进方法在压力系统中的应用。通过使用PID控制技术与模糊控制理论控制该压力系统，并利用MATLAB仿真软件对系统进行了仿真研究。 　
2011全国大学生数学建模竞赛A题题目及参考答案
04-22
2011全国大学生数学建模竞赛A题题目及参考答案
2019全国大学生数学建模竞赛题目A题 高压油管的压力控制
AtomQQ_2022313691的博客
09-13 1万+
A题 高压油管的压力控制 稳定模型 如何设置单向阀每次开启的时长，使得高压油管内的压力尽可能稳定在100MPa100MPa100MPa左右。首先要明确高压油管的工作原理以及过程。燃油经过高压油泵从A处进入高压油管，再由喷口B喷出。燃油进入和喷出的间歇性工作过程会导致高压油管内压力的变化。 压力与密度关系 题目中给出燃油的压力变化量与密度变化量成正比，比例系数为Eρ{\frac{E}{\rho} }ρE​ ，其中ρ{\rho}ρ为为燃油的密度。 Pk+1−Pk=E(Pk)ρk(ρk+1−ρk),k=0,1,
高压油管的压力控制
zsl1373的博客
08-30 2593
2019年数学建模A题 高压油管的压力控制
qq_42777804的博客
09-13 3万+
A题 高压油管的压力控制 燃油进入和喷出高压油管是许多燃油发动机工作的基础，图1给出了某高压燃油系统的工作原理，燃油经过高压油泵从A处进入高压油管，再由喷口B喷出。燃油进入和喷出的间歇性工作过程会导致高压油管内压力的变化，使得所喷出的燃油量出现偏差，从而影响发动机的工作效率。 图1 高压油管示意图 问题1. 某型号高压油管的内腔长度为500mm，内直径为10mm，供油入口A处小...
2019数学建模国赛Apython代码
qq_45737206的博客
08-06 1501
第二问： 凸轮角速度0.042，和优秀论文有差距，emmm，调了半天，大概率是前面拟合的问题。 想要图好看一点，可以自己悄悄改一下pb 代码有注释 写了很多没用的列表，可以去掉 import math import numpy as np import matplotlib.pyplot as plt t_rad = 0.027 ll = [] t_all = [] lt = [] l = 0 sign = 0 list_pb = [] list_pa = [] list_rb = [] l

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 24
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值