储油罐的变位识别附件一.doc

附件一：Matlab程序 程序1： %matlab_1.m %计算值 %计算油罐分别在偏移和没偏移时中间部分的每个高度下的实际体积并且绘出图形，同时用计算出来的偏移时的实际体积和附件中实验体积比较 close all; clear all;clc; h3=[]; %代入实测值，高度 v3=[] %代入实测值，体积 v5=[]; v5=v5+262; h5=[]; a=0.89; b=0.6; l=2.45; a1=4.1/180*pi; l1=0.825; h1=147:1:1171; %中间部分的H的范围 h1=h1/1000 %小数号 h=h1-tan(a1)*l1; % h代入for for i=1:1:1025 %有偏移，变位后，计算值 v(i)=a*b*l*( pi/2+(h(i)/b-1)*sqrt(2*h(i)/b-(h(i)/b)^2)+asin(h(i)/b-1))-(1/12*(h(i)/b-1)/sqrt(1-(h(i)/b-1)^2))*(tan(a1))^2*l^3*a/b; end; v=v'*1000; % v :1000 for i=1:1:1025 % a=0没变位，计算值 v1(i)=a*b*l*(pi/2+(h1(i)/b-1)*sqrt(2*h1(i)/b-(h1(i)/b)^2)+asin(h1(i)/b-1)) ; end; v1=v1'*1000; %1000 h1=h1*1000; % 1000 h3=h3/1000; h4=h3-tan(a1)*l1; for i=1:1:53 %把附件中实验高度代入偏移的函数中求出体积 c_v(i)=a*b*l*(pi/2+(h4(i)/b-1)*sqrt(2*h4(i)/b-(h4(i)/b)^2)+asin(h4(i)/b-1))-(1/12*(h4(i)/b-1)/sqrt(1-(h4(i)/b-1)^2))*(tan(a1))^2*l^3*a/b; end; c_v=c_v*1000 ; %理论计算值，53点 c_v=c_v'; v3; %实测值 dt_v=c_v-v3; %计算值与实测值的误差 for i=1:1:53 wc(i)=dt_v(i)/c_v(i); %计算值与实测值的误差百分比 end wc=wc'*100; ave1=sum(wc)/53; h5=h5/1000; for i=1:1:78 % a=0没变位，计算值 v6(i)=a*b*l*(pi/2+(h5(i)/b-1)*sqrt(2*h5(i)/b-(h5(i)/b)^2)+asin(h5(i)/b-1)) ; end; v6=v6'*1000; v6 h5=h5*1000; dt_v5=v6-v5; %计算值与实测值的误差 dt_v5'; for i=1:1:78 wc5(i)=dt_v5(i)/v6(i); %计算值与实测值的误差百分比 end wc5=wc5'*100 ave2=sum(wc5)/78; ave2 v_mid_low=v(1) v_mid_up=v(1025) figure plot(h5,wc5','b*'); grid on; h3=h3*1000; figure plot(h3,wc','bo'); grid on; figure plot(h1,v,'r'); %有位差，理论值 grid on; hold on; plot(h3,v3,'k'); %画实测值，53点 grid on; hold on; plot(h1,v1,'g'); %无位差值 hold on plot(h5,v5,'b'); %无位差值 hold off; 程序2： %matlab_2.m %分析变位对罐容表的影响 clear all;clc; a=0.89; b=0.6; l=2.45; a1=4.1/180*pi; l1=0.825; h1=147:1:1171; h1=h1/1000; h=h1-tan(a1)*l1; for i=1:1:1025 %变位后，理论计算值 v(i)=a*b*l*(pi/2+(h(i)/b-1)*sqrt(2*h(i)/b-(h(i)/b)^2)+asin(h(i)/b-1))-(1/12*(h(i)/b-1)/sqrt(1-(h(i)/b-1)^2))*(