1、油桶的变位识别与罐容器的标定 (国一论文)摘要：本文主要是通过积分方程和高度转化方程，建立在无变位，横向倾斜以及纵向倾斜在为 4.1 度的基本模型，出罐体变位后油位高度对于问题一：我们开始使用微积分，对灌容器的无变位和变位的求出它们油量的体积。首先求出无变位情况下求出它的高度和油量，然而在纵向倾斜，它有五种情况，分别对它五种积分，在五段分段函数，由于函数比较复杂，不便观察它们的变化，再求出的数据，次之使用使用 matlab 仿真出数据。再模拟出曲线。很明显能看到倾斜和无变位的变化。能较快的满意的答案。 最后得到出油桶的变化识别与灌容器的标定。 对于问题二：在于问题二中我们使用重积分的和高度转化
2、方程基础下求解截住图形的面积，借助于 matlab 工具绘制出模拟的图像，并且建立油位的高度与储油器之间的模型。所绘出在结果的图像在误差最大在一定范围上， 之间的角度在 010 度之间， 之间的角度在附件二的通过集优值求出。分别用 之间去检查者两百个油面的高度以及两百个的储油量的数据下的关系去精优最佳值。最后得出罐体变位后油位高度间隔为 10cm 的罐容表标定值。关键词：微积分方程， 高度转化方程，matlab，曲线拟合，体积标定，问题重述：通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定
3、的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为 =4
4、.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度 和横向偏转角度 ）之间的一般关系。请利用罐体变位后在进/ 出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。问题分析：油灌由于地基发生变化，显然油桶的形状发生变化。积分求解出纵向
5、倾斜，和横向倾斜以及无变位的油量，分析它们之间的影响。对于问题一：在问题一首先分了纵向倾斜为4.1度，无变位倾斜两种情况下的求解，在无变位的倾斜条件下，主要是一个椭圆的方程求解，当然需要建立一个积分方程求出它们的面积，再乘以他的宽度求出它的体积。分析出它们的变化，模拟出曲线。讨论灌容器的油量的之间在这种情况下变化。对于问题二，一般的加油站下的基本是这种类型的，意义不同的事如何让求出在两个倾斜角度求出储油量的体积。问题在加油站是很大的意义。问题二同样是属于积分问题，解决此类问题需要用多元面积分和体积积分。对已知的数据在附件二得出。进一步了解油面的高度和储油量之间的关系，首先对的预测和精优。在附件
6、二中先取值，再次去优化。在油桶的体积分为三个部分，积分求和，求出在10cm的灌容标的定值。并在此基础上进行误差分析和曲线拟合。进一步完善体积的模型。数学模型的假设与符号说明;模型假设： 1）在恒压，恒温，等不变条件下测量。2）油桶内壁无凹凸不齐，且内壁光滑。3）油的密度均匀，纯净无其它的液体。4）输出油管和输进油管在输量过程中不会衰落出油量。5）在测量时忽略误差的产生。6）忽略蒸发的液体以及油管留下的液体。符号说明：椭圆的半长轴长a椭圆的半短轴长b储油罐的总长L油位探针到油罐底部左侧l的距离储油罐的总体积TVd 为油桶的宽度储油罐的纵向倾斜角度储油罐的横向偏转角度油位高度h球冠体的半径R圆柱体
7、的底面半径0数学模型1.1 模型一;模型一是在无变位情况下产生，在油量输进和输出得到已知的数据，求出无变位的油面的高度跟油量的体积。油桶的形状是一个椭圆柱体所在求出的体积时需要V=s（x）*d，所以建立模型如下：4.1 小椭圆储油罐无变位时的模型由于此时的椭圆无变位，考虑先对二维椭圆进行积分。为方便表示油位高度，建立如图所示的坐标系，椭圆的半长轴长为 ，半短轴长为 ，则椭圆方程为abv=s(x)*d s（ x） =2(122)()=sin1（ )+12sin2sin1( )+2V=LS(X)V=sin1（ )+12sin2sin1( )+24.2 纵向倾斜的油桶模型我们将倾斜的油桶根据油面的高
8、度分成 5 个阶段进行积分。1=00.4sin1(tan )+12(tan) 22+()22+2)dydx (0abs(mv-detav)minv=abs(mv-detav);aa=a;bb=b;endendendaa bb minv%最小二乘拟合syms hb a bs=0;vh=zeros(1,10);aa=0;bb=0;mins=100000;v1=309.66 249.73 186.43 231.42 297.79 109.19 162.87 328.5 166.13 237.66 /1000;h1=11.34 9.20 6.84 8.40 10.87 3.99 5.95 11.92
9、6.03 8.59 /1000;vh1=-1/3000000000000000*(2*(22500000000- 9375-2137*cos(b)+37500*tan(a)(1/2);（表达式太长，部分省略）for a=linspace(0,0.1745,10)for b=linspace(0,0.1745,10)vh(1)=eval(vh1)vh(2)=eval(vh2);vh(3)=eval(vh3);vh(4)=eval(vh4);vh(5)=eval(vh5);vh(6)=eval(vh6);vh(7)=eval(vh7);vh(8)=eval(vh8);vh(9)=eval(vh9)
10、;vh(10)=eval(vh10);for i=1:10s=s+(v1(i)-vh(i)2endif sminsmins=s;aa=a;bb=b;endendendaa,bb%综合考虑纵向清晰和横向偏转syms x y %区域 1L=8;R=1.5;R0=1.625;H=1;l=2;a=0.0377;b=0.0785;vh1=2*int(y-R+R*(1-cos(b)+l*tan(a)/tan(a)*(R2-y2)0.5,R-R*(1-cos(b)-l*tan(a),R)r=(R02-(x-R)2)0.5vh2=int(r2*acos(R0-H)/r)-(r2-(R0-H)2)0.5)*(R
11、0-H),x,0,R*(1-cos(b)+l*tan(a)vh=vh1+vh2;eval(vh)syms x y %区域 2L=8;R=1.5;R0=1.625;H=1;l=2;a=0.0377;b=0.0785;hb2max=(L-l)*tan(a)-R*(1-cos(b)/cos(b);hb=0.2224;h=R-(R-hb)*cos(b);vh1=2*int(y-R+h+l*tan(a)/tan(a)*(R2-y2)0.5,y,R-h-l*tan(a),R);r=(R02-(x-R)2)0.5;vh2=int(r2*acos(R0-H)/r)-(r2-(R0-H)2)0.5)*(R0-H
12、),x,0,h+l*tan(a);vh=vh1+vh2;eval(vh)syms x y %区域 3L=8;R=1.5;R0=1.625;H=1;l=2;a=0.0377;b=0.0785;hb=2.93;h=R-(R-hb)*cos(b);vh1=2*L*int(R2-(y-R)2)0.5,y,0,h);vh11=2*int(-y+h+l*tan(a)/tan(a)*(R2-(y-R)2)0.5,y,h,h+l*tan(a);vh12=2*int(L-(h-y+l*tan(a)/tan(a)*(R2-(y-R)2)0.5,y,h-(L-l)*tan(a),h);r=(R02-(y-R)2)0
13、.5;vh21=int(r2*acos(R0-H)/r)-(r2-(R0-H)2)0.5)*(R0-H),y,0,h+l*tan(a);vh22=int(r2*acos(R0-H)/r)-(r2-(R0-H)2)0.5)*(R0-H),y,0,h-(L-l)*tan(a);vh=vh1+vh11-vh12+vh21+vh22;eval(vh)syms x y %区域 4L=8;R=1.5;R0=1.625;H=1;l=2;a=4*pi/180;hb=3;hb4min=(R*(1+cos(b)-l*tan(a)/cos(b)bb4max=2*Rh=R-(R-hb)*cos(b);v1=L*pi*
14、R2;vh1=2*int(y+R-h+(L-l)*tan(a)/tan(a)*(R2-y2)0.5,y,h-R-(L-l)*tan(a),R);r=(R02-(y-R)2)0.5;vh21=int(r2*acos(R0-H)/r)-(r2-(R0-H)2)0.5)*(R0-H),y,0,h-(L-l)*tan(a);vh22=pi*H2*(R0-H/3);vh=v1-vh1+vh21+vh22;eval(vh)syms x y %区域 5L=8;R=1.5;R0=1.625;H=1;l=2;a=4*pi/180;h=3;h=R-(R-hb)*cos(b);v1=L*pi*R2;vh1=2*in
15、t(y-R+(L-l)*tan(a)/tan(a)*(R2-y2)0.5,y,R-(L-l)*tan(a),R);r=(R02-(y-R)2)0.5;vh21=int(r2*acos(R0-H)/r)-(r2-(R0-H)2)0.5)*(R0-H),y,0,2*R-(L-l)*tan(a);vh22=pi*H2*(R0-H/3);vh=v1-vh1+vh21+vh22;eval(vh)%对估计的参数得到的体积变化和实际储油量的对比图syms hbh1= 2486.21 426.96 420.01 416.53 413.98 /1000;v=zeros(1,300);vs=-2431987916*(241+1312*hb)+.+(-409203/7036874664+8938596496*hb)(1/2)for i=1:300
