储油罐问题.doc

问题重述 通常加油站都采用与储油罐配套的“油位计量管理系统”。采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。 使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），油箱液面不再平行于油箱底部，从而使得容积的计算公式不再准确，导致罐容表发生改变，需要重新定位。 图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。要求用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 （1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为(=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 （2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度(和横向偏转角度( ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。 问题一 无变位情形 建立坐标系 将椭圆柱体的储油箱水平放置，箱底接触地面的只有一条直线，以该条直线为x轴，油位探针所在直线为y轴，两轴所在直线相交与O点，过O点做z轴垂直于x,y轴所确定的平面。 2公式推导 在z-y平面内列写椭圆方程： ； 椭圆中阴影部分面积为； 油液容积为。 根据如上的公式，整理得到容积V与高度h之间的函数关系式为： 其中 a 0.89, b=0.6, l=2.45 (m)分别为罐体截面椭圆长半轴、短半轴和罐体长度， h 为罐内的油位高度油量/L 3659.88 3946.55 4110.15 油量/L 3659.88 3946.55 4110.15 ，因此，油平面法向变为（，油平面方程为 无变位时，油平面高度不随x坐标变化，当有变位时，油平面高度是x的函数，表达式为 计算体积的Matlab程序 function v= getv(h) a=0.89;b=0.6; alpha=4.1*3.1416/18; fun=@(x,y)a*sqrt(1-((y-b)/b).^2); hx=@(x)min(max(0,h-x*tan(alpha)),1.2); v=2*quad2d(fun,-0.4,2.05,0,hx); end 或者 function v= getv2(h) a=0.89;b=0.6;alpha=4.1*3.1416/180; fun=@(y)a*sqrt(1-((y-b)/b).^2); hx=@(x)min(max(0,h-x*tan(alpha)),1.2); f1=@(u)quadl(fun,0,hx(u)); f2=@(x)arrayfun(f1,x); v=2*quadl(f2,-0.4,2.05); end 问题二 对于实际的储油罐，坐标系的建立类似于问题一，将y轴建立在油位探针上，原点取油位探针的最低点，x轴为圆柱体水平放置时和地面的接触线，因此，储油罐左边球缺的x坐标范围为（-3，-2），右边球缺的x 坐标范围是（6，7），中间圆柱体部分的x坐标范围是（-2,6）， 圆柱体底圆半径为1.5米，容易计算出球缺的半径为1.625米 当储油罐发生纵向旋转角度，相当于坐标系关于z轴旋转了，旋转矩阵为 油平面法向由（0,1,0）变为，而当储油罐横向旋转时，相当于坐标系关于x轴旋转了，旋转矩阵为 此时方向由变为。 不管罐体如何旋转，油平面总是经过点（0,h,0） 圆柱体油量的计算 用垂直于x轴的平面截储油罐，截得是圆，圆半径为 根据油平面方程，油平面和xy面交线为 解出ｙ为，考虑到球缺上的点和 xz面的距离，记. 当储油罐横向偏转时，油体和垂直于x轴的截面为弓形，其高度为 由于我们只在储油罐所在区域考虑问题，应作修正，修正后的表达式为 由弓形的面积公式，我们得到油体和垂直于x轴的截面的截面积为 于是，油量为 变位参数的确定 数据给定了各时间段的进出油量值?及显示油位高度?，据此可以计算参数。将油位高度值代入体积表达式中，计算得到?，于是得到,通过最小二乘准则，利用数