- 1、本文档共8页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
储油罐的变位识别附件一
附件一: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))*(tan(a1))^2*l^3;
end;
文档评论(0)