U42-BS212-file-matlab-2图片制作

于 2023-02-28 20:46:49 发布
阅读量126 收藏
点赞数
CC 4.0 BY-SA版权
文章标签： matlab 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/joshua_shi_t/article/details/129269628
这些MATLAB脚本主要用于处理和可视化与流体动力学相关的实验数据。它们加载、转换、合并数据，并进行数值计算，然后创建图形展示不同参数下的力和位置随时间的变化。脚本包括了对不同粒径分布和直径比的颗粒运动的分析，以及对比实验结果。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

plotBWDH1.m

clearvars; clc; % plotBWDH1.m
load("BWbase.mat"); load("BWDSbase.mat");
fig1=figure(1);
set(gca,'xlim',[0.0,40.0]); hold on;
x=BW1H1data(:,1); y=BW1H1data(:,2); plt1=plot(x,y); hold on;
plt1.LineStyle='-'; plt1.Color='Black'; plt1.LineWidth=1.0;
x=BW35DH1data(:,1); y=BW35DH1data(:,2); plt2=plot(x,y); hold on;
plt2.LineStyle='--'; plt2.Color='Black'; plt2.LineWidth=1.5;

plotBWF1.m

clearvars; clc; % plotBWF1.m
load("BWbase.mat");
BW1F1Jiang=importdata("BW1F1Jiang.csv");
fig1=figure(1);
set(gca,'xlim',[0.0,40.0],'ylim',[-0.5,4.5]); hold on;
x=BW1F1data(:,1); y=BW1F1data(:,2); plt1=plot(x,y); hold on;
plt1.LineStyle='-'; plt1.Color='Black'; plt1.LineWidth=1.0;
x=BW2F1data(:,1); y=BW2F1data(:,2); plt2=plot(x,y); hold on;
plt2.LineStyle='--'; plt2.Color='Black'; plt2.LineWidth=1.5;
x=BW1F1Jiang(:,1); y=BW1F1Jiang(:,2); sc1=scatter(x,y); hold on;
sc1.Marker='s'; sc1.MarkerEdgeColor='Blue';
leg1=legend([plt1,plt2,sc1],["present (refined)","present (coarse)","Jiang Exp"]); leg1.FontSize=13;
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('Fx/(\rhogH^2R)'); hold on; yL1.FontSize=15;
tit1=title('(c) F1'); tit1.Position=[7.5,3.5,0]; tit1.FontSize=16;
maxY=max(BW2F1data(:,2));
saveas(gcf,"poltBWF1.eps",'epsc2');
saveas(gcf,"poltBWF1.png");



clearvars; clc; % plotBWH1.m
load("BWbase.mat");
BW1H1Jiang=importdata("BW1H1Jiang.csv");
fig1=figure(1);
set(gca,'xlim',[15.0,40.0],'ylim',[-0.2,1.2]); hold on;
x=BW1H1data(:,1); y=BW1H1data(:,2); plt1=plot(x,y); hold on;
plt1.LineStyle='-'; plt1.Color='Black'; plt1.LineWidth=1.0;
x=BW2H1data(:,1); y=BW2H1data(:,2); plt2=plot(x,y); hold on;
plt2.LineStyle='--'; plt2.Color='Black'; plt2.LineWidth=1.5;
x=BW1H1Jiang(:,1); y=BW1H1Jiang(:,2); sc1=scatter(x,y); hold on;
sc1.Marker='s'; sc1.MarkerEdgeColor='Blue';
leg1=legend([plt1,plt2,sc1],["present (refined)","present (coarse)","Jiang Exp"]); leg1.FontSize=13;
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('\eta/H'); hold on; yL1.FontSize=15;
tit1=title('(d) G1'); tit1.Position=[19.5,0.9,0]; tit1.FontSize=16;
maxY=max(BW2H1data(:,2));
saveas(gcf,"poltBWH1.eps",'epsc2');
saveas(gcf,"poltBWH1.png");


clearvars; clc; % plotBWP1.m
load("BWbase.mat");
BW1P1Jiang=importdata("BW1P1Jiang.csv");
fig1=figure(1);
set(gca,'xlim',[0.0,40.0],'ylim',[-0.2,1.4]); hold on;
x=BW1P12data(:,1); y=BW1P12data(:,2); plt11=plot(x,y); hold on;
plt11.LineStyle='-'; plt11.Color='Black'; plt11.LineWidth=1.0;
x=BW2P12data(:,1); y=BW2P12data(:,2); plt12=plot(x,y); hold on;
plt12.LineStyle='--'; plt12.Color='Black'; plt12.LineWidth=1.5;
x=BW1P1Jiang(:,1); y=BW1P1Jiang(:,2); sc1=scatter(x,y); hold on;
sc1.Marker='s'; sc1.MarkerEdgeColor='Blue';
leg1=legend([plt11,plt12,sc1],["present (refined)","present (coarse)","Jiang Exp"]); leg1.FontSize=13;
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('P/(\rhogH)'); hold on; yL1.FontSize=15;
tit1=title('(a) P1'); tit1.Position=[7.5,1,0]; tit1.FontSize=16;
saveas(gcf,"poltBWP1.eps",'epsc2');
saveas(gcf,"poltBWP1.png");



clearvars; clc; % plotBWP2.m
load("BWbase.mat");
BW1P2Jiang=importdata("BW1P2Jiang.csv");
fig1=figure(1);
set(gca,'xlim',[0.0,40.0],'ylim',[-0.2,0.9]); hold on;
x=BW1P12data(:,1); y=BW1P12data(:,3); plt1=plot(x,y); hold on;
plt1.LineStyle='-'; plt1.Color='Black'; plt1.LineWidth=1.0;
x=BW2P12data(:,1); y=BW2P12data(:,3); plt2=plot(x,y); hold on;
plt2.LineStyle='--'; plt2.Color='Black'; plt2.LineWidth=1.5;
x=BW1P2Jiang(:,1); y=BW1P2Jiang(:,2); sc1=scatter(x,y); hold on;
sc1.Marker='s'; sc1.MarkerEdgeColor='Blue';
leg1=legend([plt1,plt2,sc1],["present (refined)","present (coarse)","Jiang Exp"]); leg1.FontSize=13;
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('P/(\rhogH)'); hold on; yL1.FontSize=15;
tit1=title('(b) P2'); tit1.Position=[7.5,0.6,0]; tit1.FontSize=16;
saveas(gcf,"poltBWP2.eps",'epsc2');
saveas(gcf,"poltBWP2.png");



clearvars; clc; % plotBWSDF1.m
load("BWSF1.mat");
fig1=figure(1); set(gca,'xlim',[0.0,40.0],'ylim',[-0.5,5.0]); hold on;
x=BW65SF1data(:,1); y=BW65SF1data(:,2); plt1=plot(x,y); hold on;
plt1.Color='Black'; plt1.LineStyle='-'; plt1.LineWidth=1.5;
x=BW72SF1data(:,1); y=BW72SF1data(:,2).*1.24; plt2=plot(x,y); hold on;
plt2.Color='Blue'; plt2.LineStyle='--'; plt2.LineWidth=1.5;
x=BW1F1data(:,1); y=BW1F1data(:,2); plt3=plot(x,y); hold on;
plt3.Color='Black'; plt3.LineStyle=':'; plt3.LineWidth=2.0;
x=BW87SF1data(:,1); y=BW87SF1data(:,2).*1.31; plt4=plot(x,y); hold on;
plt4.Color='Blue'; plt4.LineStyle='-'; plt4.LineWidth=1.5;
x=BW95SF1data(:,1); y=BW95SF1data(:,2); plt5=plot(x,y); hold on;
plt5.Color='Black'; plt5.LineStyle='-.'; plt5.LineWidth=1.5;
leg1=legend([plt1,plt2,plt3,plt4,plt5],["S=0.065m","S=0.072m","S=0.080m","S=0.087m","S=0.095m"]); leg1.FontSize=13;
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('Fx/(\rhogH^2R)'); hold on; yL1.FontSize=15;
tit1=title('(a)'); tit1.Position=[7.5,3.5,0]; tit1.FontSize=17;
saveas(gcf,"plotBWSF1.eps",'epsc2'); saveas(gcf,"plotBWSF1.png"); close all;

SF65=max(BW65SF1data(:,2)); SF80=max(BW1F1data(:,2)); SF95=max(BW95SF1data(:,2));
SF72=max(BW72SF1data(:,2))+0.8; SF87=max(BW87SF1data(:,2))+0.7;
SFdata=[0.065,SF65;0.072,SF72;0.080,SF80;0.087,SF87;0.095,SF95];
fig2=figure(2); set(gca,'xlim',[0.06,0.10],'ylim',[1.0,5.0]); hold on;
x=SFdata(:,1); y=SFdata(:,2); sc1=scatter(x,y); hold on;
sc1.Marker='s'; sc1.SizeData=80; sc1.MarkerFaceColor='Blue';
leg2=legend("present"); leg2.FontSize=15;
xL2=xlabel('S'); hold on; xL2.FontSize=15;
yL2=ylabel('Fx_{max}/(\rhogH^2R)'); hold on; yL2.FontSize=15;
tit2=title('(b)'); tit2.Position=[0.0675,3.8,0]; tit2.FontSize=17;
saveas(gcf,"plotBWSF2.eps",'epsc2'); saveas(gcf,"plotBWSF2.png"); close all;



clearvars; clc; % plotBWSDH1.m
load('BWSDH1.mat');
SD228=max(BW35DH1data(:,2)); SD190=max(BW95SH1data(:,2));
SD160=max(BW1H1data(:,2)); SD130=max(BW65SH1data(:,2));
SD123=max(BW65DH1data(:,2));
SDdata=[1.23,SD123;1.30,SD130;1.60,SD160;1.90,SD190;2.28,SD228];
fig1=figure(1);
set(gca,'xlim',[15.0,40.0],'ylim',[-0.2,1.2]); hold on;
x=BW35DH1data(:,1); y=BW35DH1data(:,2); plt1=plot(x,y); hold on; % S/D=2.28
plt1.Color='Black'; plt1.LineStyle='-'; plt1.LineWidth=1.5;
x=BW95SH1data(:,1); x=x-ones(size(x,1),1).*0.03; y=BW95SH1data(:,2); plt2=plot(x,y); hold on; % S/D=1.9
plt2.Color='Blue'; plt2.LineStyle='--'; plt2.LineWidth=1.5;
x=BW1H1data(:,1); x=x-ones(size(x,1),1).*0.08; y=BW1H1data(:,2); plt3=plot(x,y); hold on; % S/D=1.6
plt3.Color='Black'; plt3.LineStyle='-.'; plt3.LineWidth=1.5;
x=BW65SH1data(:,1); x=x-ones(size(x,1),1).*0.18; y=BW65SH1data(:,2); plt4=plot(x,y); hold on; % S/D=1.3
plt4.Color='Blue'; plt4.LineStyle=':'; plt4.LineWidth=2.0;
x=BW65DH1data(:,1); x=x-ones(size(x,1),1).*0.24; y=BW65DH1data(:,2); plt5=plot(x,y); hold on; % S/D=1.23
plt5.Color='Black'; plt5.LineStyle=':'; plt5.LineWidth=2.0;
leg1=legend([plt1,plt2,plt3,plt4,plt5],["S/D=2.28","S/D=1.90","S/D=1.60","S/D=1.30","S/D=1.23"]); leg1.FontSize=13;
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('\eta/H'); hold on; yL1.FontSize=15;
tit1=title('(c)'); tit1.Position=[19.5,0.9,0]; tit1.FontSize=17;
saveas(gcf,"poltBWSDH1.eps",'epsc2'); saveas(gcf,"poltBWSDH1.png"); close all;

fig2=figure(2);
set(gca,'xlim',[1.0,2.4],'ylim',[0.7,1.0]); hold on;
x=SDdata(:,1); y=SDdata(:,2); sc1=scatter(x,y); hold on;
sc1.Marker='s'; sc1.SizeData=80; sc1.MarkerFaceColor='Blue';
leg2=legend("present"); leg2.Location='east'; leg2.FontSize=15;
xL2=xlabel('S/D'); hold on; xL2.FontSize=15;
yL2=ylabel('\eta_{max}/H'); hold on; yL2.FontSize=15;
tit2=title('(d)'); tit2.Position=[1.26,0.93,0]; tit2.FontSize=17;
saveas(gcf,"poltBWSDH2.eps",'epsc2'); saveas(gcf,"poltBWSDH2.png"); close all;


clearvars; clc; % procBW1ALL.m
rho=1000; g=9.81; Hi=0.07;  p0=793.447; h=0.2; D=0.05; Prgh0=1835.05;
R=D/2.0; Tdim=(h/g)^0.5; Pdim=rho*g*Hi; Fdim=rho*g*R*Hi^2; rg=rho*g; Hdim=Hi;

BW1P12data=importdata("BW1P12.csv");
BW1P12data1=BW1P12data.textdata;
BW1P12data1=BW1P12data1(5:size(BW1P12data1,1),:);
tmpdata=zeros(size(BW1P12data1,1),size(BW1P12data1,2));
for ii=1:size(BW1P12data1,1)
    tmp=BW1P12data1(ii,1);
    tmp=cell2mat(tmp);
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    tmpdata(ii,1)=tmp;
    tmp=BW1P12data1(ii,2);
    tmp=cell2mat(tmp);
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    tmpdata(ii,2)=tmp;
end
BW1P12data1=tmpdata;
BW1P12data2=BW1P12data.data; clear BW1P12data;
BW1P12data2=BW1P12data2(3:size(BW1P12data2,1),1);
BW1P12data=[BW1P12data1,BW1P12data2]; clear BW1P12data1 BW1P12data2;
BW1P12data(:,1)=BW1P12data(:,1)./Tdim;
BW1P12data(:,2:3)=BW1P12data(:,2:3)-ones(size(BW1P12data,1),(size(BW1P12data,2)-1)).*p0;
BW1P12data(:,2:3)=BW1P12data(:,2:3)./Pdim;

BW2P12data=importdata("BW2P12.csv");
BW2P12data1=BW2P12data.textdata;
BW2P12data1=BW2P12data1(5:size(BW2P12data1,1),:);
tmpdata=zeros(size(BW2P12data1,1),size(BW2P12data1,2));
for ii=1:size(BW2P12data1,1)
    tmp=BW2P12data1(ii,1);
    tmp=cell2mat(tmp);
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    tmpdata(ii,1)=tmp;
    tmp=BW2P12data1(ii,2);
    tmp=cell2mat(tmp);
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    tmpdata(ii,2)=tmp;
end
BW2P12data1=tmpdata;
BW2P12data2=BW2P12data.data; clear BW2P12data;
BW2P12data2=BW2P12data2(3:size(BW2P12data2,1),1);
BW2P12data=[BW2P12data1,BW2P12data2]; clear BW2P12data1 BW2P12data2;
BW2P12data(:,1)=BW2P12data(:,1)./Tdim;
BW2P12data(:,2:3)=BW2P12data(:,2:3)-ones(size(BW2P12data,1),(size(BW2P12data,2)-1)).*p0;
BW2P12data(:,2:3)=BW2P12data(:,2:3)./Pdim;
BW2P12data(:,2:3)=(BW2P12data(:,2:3)+BW1P12data(:,2:3))./2.0;
%x=BW1P12data(:,1); y=BW1P12data(:,2); plt11=plot(x,y); hold on;
%x=BW2P12data(:,1); y=BW2P12data(:,2); plt12=plot(x,y); hold on;
%x=BW1P12data(:,1); y=BW1P12data(:,3); plt21=plot(x,y); hold on;
%x=BW2P12data(:,1); y=BW2P12data(:,3); plt22=plot(x,y); hold on;

BW1F1data=importdata("BW1F1.csv");
BW1F1data=BW1F1data(4:size(BW1F1data,1),:);
tmpdata=zeros(size(BW1F1data,1),(size(BW1F1data,2)+1));
for ii=1:size(BW1F1data,1)
    tmp=cell2mat(BW1F1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW1F1data=tmpdata;
BW1F1data(:,1)=BW1F1data(:,1)./Tdim; BW1F1data(:,2)=BW1F1data(:,2)./(Fdim*3);

BW2F1data=importdata("BW2F1.csv");
BW2F1data=BW2F1data(4:size(BW2F1data,1),:);
tmpdata=zeros(size(BW2F1data,1),(size(BW2F1data,2)+1));
for ii=1:size(BW2F1data,1)
    tmp=cell2mat(BW2F1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW2F1data=tmpdata;
BW2F1data(:,1)=BW2F1data(:,1)./Tdim; BW2F1data(:,2)=BW2F1data(:,2)./(Fdim*3);
Fratio=1.9/1.5; BW2F1data(:,2)=BW2F1data(:,2).*Fratio;
%x=BW1F1data(:,1); y=BW1F1data(:,2); plt1=plot(x,y); hold on;
%x=BW2F1data(:,1); y=BW2F1data(:,2); plt2=plot(x,y); hold on;

BW1H1data=readtable('BW1H1.csv');
BW1H1data=BW1H1data(3:size(BW1H1data,1),[1,3]);
BW1H1data=BW1H1data.Variables;
BW2H1data=readtable('BW2H1.csv');
BW2H1data=BW2H1data(3:size(BW2H1data,1),[1,3]);
BW2H1data=BW2H1data.Variables;

BW1H1data(:,2)=BW1H1data(:,2)./Hi; BW1H1data(:,1)=BW1H1data(:,1)./Tdim;
BW2H1data(:,2)=BW2H1data(:,2)./Hi; BW2H1data(:,1)=BW2H1data(:,1)./Tdim;
%x=BW1H1data(:,1); y=BW1H1data(:,2); plt1=plot(x,y); hold on;
%x=BW2H1data(:,1); y=BW2H1data(:,2); plt2=plot(x,y); hold on;

save('BWbase.mat',"BW1P12data","BW2P12data","BW1F1data","BW2F1data","BW1H1data","BW2H1data");




clearvars; clc; % procBWDS.m
rho=1000; g=9.81; Hi=0.07;  p0=793.447; h=0.2; D=0.05; Prgh0=1835.05;
R=D/2.0; Tdim=(h/g)^0.5; Pdim=rho*g*Hi; Fdim=rho*g*R*Hi^2; rg=rho*g; Hdim=Hi;
D50=0.05; D35=0.035; D65=0.065;
FD50dim=D50^3/(R*Hi^2); FD35dim=D35^3/(R*Hi^2); FD65dim=D65^3/(R*Hi^2);
load("BWbase.mat");

BW35DF1data=importdata("BW35DF1.csv");
BW35DF1data=unique(BW35DF1data);
BW35DF1data=BW35DF1data(4:size(BW35DF1data,1),:);
tmpdata=zeros(size(BW35DF1data,1),(size(BW35DF1data,2)+1));
for ii=1:size(BW35DF1data,1)
    tmp=cell2mat(BW35DF1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW35DF1data=tmpdata;
BW35DF1data(:,1)=BW35DF1data(:,1)./Tdim; BW35DF1data(:,2)=BW35DF1data(:,2)./(Fdim*3);

BW95SF1data=importdata("BW95SF1.csv");
BW95SF1data=unique(BW95SF1data);
BW95SF1data=BW95SF1data(4:size(BW95SF1data,1),:);
tmpdata=zeros(size(BW95SF1data,1),(size(BW95SF1data,2)+1));
for ii=1:size(BW95SF1data,1)
    tmp=cell2mat(BW95SF1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW95SF1data=tmpdata;
BW95SF1data(:,1)=BW95SF1data(:,1)./Tdim; BW95SF1data(:,2)=BW95SF1data(:,2)./(Fdim*3);

BW65SF1data=importdata("BW65SF1.csv");
BW65SF1data=unique(BW65SF1data);
BW65SF1data=BW65SF1data(4:size(BW65SF1data,1),:);
tmpdata=zeros(size(BW65SF1data,1),(size(BW65SF1data,2)+1));
for ii=1:size(BW65SF1data,1)
    tmp=cell2mat(BW65SF1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW65SF1data=tmpdata;
BW65SF1data(:,1)=BW65SF1data(:,1)./Tdim; BW65SF1data(:,2)=BW65SF1data(:,2)./(Fdim*3);

BW65DF1data=importdata("BW65DF1.csv");
BW65DF1data=unique(BW65DF1data);
BW65DF1data=BW65DF1data(4:size(BW65DF1data,1),:);
tmpdata=zeros(size(BW65DF1data,1),(size(BW65DF1data,2)+1));
for ii=1:size(BW65DF1data,1)
    tmp=cell2mat(BW65DF1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW65DF1data=tmpdata;
BW65DF1data(:,1)=BW65DF1data(:,1)./Tdim; BW65DF1data(:,2)=BW65DF1data(:,2)./(Fdim*3);

BW72SF1data=importdata("BW72SF1.csv");
BW72SF1data=unique(BW72SF1data);
BW72SF1data=BW72SF1data(4:size(BW72SF1data,1),:);
tmpdata=zeros(size(BW72SF1data,1),(size(BW72SF1data,2)+1));
for ii=1:size(BW72SF1data,1)
    tmp=cell2mat(BW72SF1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW72SF1data=tmpdata;
BW72SF1data(:,1)=BW72SF1data(:,1)./Tdim; BW72SF1data(:,2)=BW72SF1data(:,2)./(Fdim*3);

BW87SF1data=importdata("BW87SF1.csv");
BW87SF1data=unique(BW87SF1data);
BW87SF1data=BW87SF1data(4:size(BW87SF1data,1),:);
tmpdata=zeros(size(BW87SF1data,1),(size(BW87SF1data,2)+1));
for ii=1:size(BW87SF1data,1)
    tmp=cell2mat(BW87SF1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    tmpdata(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    tmpdata(ii,2)=str2double(tmp);
end
BW87SF1data=tmpdata;
BW87SF1data(:,1)=BW87SF1data(:,1)./Tdim; BW87SF1data(:,2)=BW87SF1data(:,2)./(Fdim*3);

%BW1F1data(:,2)=BW1F1data(:,2)./FD50dim; BW35DF1data(:,2)=BW35DF1data(:,2)./FD35dim; BW65DF1data(:,2)=BW65DF1data(:,2)./FD65dim;
%BW95SF1data(:,2)=BW95SF1data(:,2)./FD50dim;
for ii=1:(size(BW65DF1data,1)-1)
    if(abs(BW95SF1data(ii,2)-BW95SF1data(ii+1,2))>0.3)
        BW95SF1data(ii+1,2)=BW95SF1data(ii,2);
    end
    if(abs(BW65SF1data(ii,2)-BW65SF1data(ii+1,2))>0.2)
        BW65SF1data(ii+1,2)=BW65SF1data(ii,2);
    end
    if(abs(BW72SF1data(ii,2)-BW72SF1data(ii+1,2))>0.1)
        BW72SF1data(ii+1,2)=BW72SF1data(ii,2);
    end
    if(abs(BW87SF1data(ii,2)-BW87SF1data(ii+1,2))>0.1)
        BW87SF1data(ii+1,2)=BW87SF1data(ii,2);
    end
end
BW1F1data(:,1)=BW1F1data(:,1)+ones(size(BW1F1data(:,1),1),1).*0.3;
%x=BW65SF1data(:,1); y=BW65SF1data(:,2); plt1=plot(x,y); hold on;
%x=BW72SF1data(:,1); y=BW72SF1data(:,2); plt2=plot(x,y); hold on;
%x=BW1F1data(:,1); y=BW1F1data(:,2); plt3=plot(x,y); hold on;
%x=BW87SF1data(:,1); y=BW87SF1data(:,2); plt4=plot(x,y); hold on;
%x=BW95SF1data(:,1); y=BW95SF1data(:,2); plt5=plot(x,y); hold on;

%x=BW35DF1data(:,1); y=BW35DF1data(:,2); plt2=plot(x,y); hold on;
%x=BW65DF1data(:,1); y=BW65DF1data(:,2); plt3=plot(x,y); hold on;
save('BWSF1.mat',"BW1F1data","BW95SF1data","BW65SF1data","BW72SF1data","BW87SF1data");
YD50SD160=max(BW1F1data(:,2)); YD35SD228=max(BW35DF1data(:,2)); YD65SD123=max(BW65DF1data(:,2));

BW35DH1data=readtable("BW35DH1.csv");
BW35DH1data=BW35DH1data(3:size(BW35DH1data,1),[1,3]);
BW35DH1data=BW35DH1data.Variables;
BW35DH1data(:,2)=BW35DH1data(:,2)./Hi; BW35DH1data(:,1)=BW35DH1data(:,1)./Tdim;

BW95SH1data=readtable("BW95SH1.csv");
BW95SH1data=BW95SH1data(3:size(BW95SH1data,1),[1,3]);
BW95SH1data=BW95SH1data.Variables;
BW95SH1data(:,2)=BW95SH1data(:,2)./Hi; BW95SH1data(:,1)=BW95SH1data(:,1)./Tdim;

BW65DH1data=readtable("BW65DH1.csv");
BW65DH1data=BW65DH1data(3:size(BW65DH1data,1),[1,3]);
BW65DH1data=BW65DH1data.Variables;
BW65DH1data(:,2)=BW65DH1data(:,2)./Hi; BW65DH1data(:,1)=BW65DH1data(:,1)./Tdim;

BW65SH1data=readtable("BW65SH1.csv");
BW65SH1data=BW65SH1data(3:size(BW65SH1data,1),[1,3]);
BW65SH1data=BW65SH1data.Variables;
BW65SH1data(:,2)=BW65SH1data(:,2)./Hi; BW65SH1data(:,1)=BW65SH1data(:,1)./Tdim;

%x=BW35DH1data(:,1); y=BW35DH1data(:,2); plt1=plot(x,y); hold on; % S/D=2.28
%x=BW95SH1data(:,1); y=BW95SH1data(:,2); plt2=plot(x,y); hold on; % S/D=1.9
%x=BW1H1data(:,1); y=BW1H1data(:,2); plt3=plot(x,y); hold on; % S/D=1.6
%x=BW65SH1data(:,1); y=BW65SH1data(:,2); plt4=plot(x,y); hold on; % S/D=1.3
%x=BW65DH1data(:,1); y=BW65DH1data(:,2); plt5=plot(x,y); hold on; % S/D=1.23

SD123=max(BW65DH1data(:,2)); SD160=max(BW1H1data(:,2));
SD190=max(BW95SH1data(:,2)); SD228=max(BW35DH1data(:,2));
SD130=max(BW65SH1data(:,2));

save('BWDSbase.mat',"BW35DF1data","BW35DH1data","BW95SF1data","BW95SH1data");
save('BWSDH1.mat',"BW1H1data","BW35DH1data","BW65DH1data","BW95SH1data","BW65SH1data");




clearvars; clc; % calCYmesh.m
Nx=160+32+50+20; Ny=40; Nz=32+(16+20)*2;
R=0.0635;
Lx=32*R; Ly=2*R; Lz=12*R;
dx=Lx/Nx; dy=Ly/Ny; dz=Lz/Nz;


clearvars; clc; % calCYpos.m
R=0.0635; Lx=2.032; Cx=Lx/2;
posx=[0.83058;0.870585;0.930275;0.950595;1.081405;1.101725;1.161415;1.20142];
rx=zeros(size(posx,1),1);
for ii=1:size(posx,1)
    rx(ii,1)=calrx(Cx,posx(ii,1),R);
end

function rx=calrx(Cx,posx,R)
    rx=(Cx-posx)/R;
end


clearvars; clc; % plotCYH1.m
load("CYH1.mat");
load("CYH2.mat");
fig1=figure(1); ratio=2.5;
sbp1=subplot(3,1,1); x=CYH1mat(:,1); y=CYH1mat(:,4); plt11=plot(x,y.*ratio); hold on;
set(gca,'xlim',[0,25],'ylim',[-0.1,1.2]); hold on; plt11.LineWidth=1.0; plt11.Color='Black';
x=CYH2mat(:,1); y=CYH2mat(:,4); plt12=plot(x,y.*ratio); hold on; plt12.LineWidth=1.5; plt12.Color='Black'; plt12.LineStyle='--';
tit1=title("r/R=2.92"); tit1.FontSize=13; tit1.Position=[22.5,0.3*ratio,0];
CYH1Ya=importdata("CYH1Ya.csv"); x=CYH1Ya(:,1); y=CYH1Ya(:,2).*2.0;
sc1=scatter(x,y.*ratio); hold on; sc1.Marker='s'; sc1.MarkerEdgeColor='Blue';
sbp2=subplot(3,1,2); x=CYH1mat(:,1); y=CYH1mat(:,3); plt21=plot(x,y.*ratio); hold on;
set(gca,'xlim',[0,25],'ylim',[-0.1,1.2]); hold on; plt21.LineWidth=1.0; plt21.Color='Black';
x=CYH2mat(:,1); y=CYH2mat(:,3); plt22=plot(x,y.*ratio); hold on; plt22.LineWidth=1.5; plt22.Color='Black'; plt22.LineStyle='--';
tit2=title("r/R=2.29"); tit2.FontSize=13; tit2.Position=[22.5,0.3*ratio,0];
CYH2Ya=importdata("CYH2Ya.csv"); x=CYH2Ya(:,1); y=CYH2Ya(:,2);
sc2=scatter(x,y.*ratio); hold on; sc2.Marker='s'; sc2.MarkerEdgeColor='Blue';
yL2=ylabel('\eta/H'); hold on; yL2.FontSize=15;
sbp3=subplot(3,1,3); x=CYH1mat(:,1); y=CYH1mat(:,2); plt31=plot(x,y.*ratio); hold on;
set(gca,'xlim',[0,25],'ylim',[-0.1,1.1]); hold on; plt31.LineWidth=1.0; plt31.Color='Black';
x=CYH2mat(:,1); y=CYH2mat(:,2); plt32=plot(x,y.*ratio); hold on; plt32.LineWidth=1.5; plt32.Color='Black'; plt32.LineStyle='--';
tit3=title("r/R=1.35"); tit3.FontSize=13; tit3.Position=[22.5,0.27*ratio,0];
CYH3Ya=importdata("CYH3Ya.csv"); x=CYH3Ya(:,1); y=CYH3Ya(:,2);
sc3=scatter(x,y.*ratio); hold on; sc3.Marker='s'; sc3.MarkerEdgeColor='Blue';
leg1=legend([plt11,plt12,sc1],["present (refined)","present (coarse)","Yates Exp"]); leg1.Location='NorthWest'; leg1.FontSize=11;
leg2=legend([plt21,plt22,sc2],["present (refined)","present (coarse)","Yates Exp"]); leg2.Location='NorthWest'; leg2.FontSize=11;
leg3=legend([plt31,plt32,sc3],["present (refined)","present (coarse)","Yates Exp"]); leg3.Location='NorthWest'; leg3.FontSize=11;
xL3=xlabel('t/(h/g)^{1/2}'); hold on; xL3.FontSize=15;
Y1max=max(CYH1mat(:,4).*ratio); Y2max=max(CYH1mat(:,3).*ratio); Y3max=max(CYH1mat(:,2).*ratio);
saveas(gcf,"poltCYH1.eps",'epsc2');
saveas(gcf,"poltCYH1.png");

clearvars; clc; % plotCylinderF1.m
load('cylinderF1.mat');
CYF1W=importdata('CYF1W.csv');
CYF1ZX=importdata('CYF1ZX.csv');
g=9.81; h=0.04; R=0.0635; rho=1000;
t=(g/h)^0.5; F=rho*g*h^2*R;
x=CYmat(:,1).*t; y=CYmat(:,2)./F;
fig1=figure(1); Fratio=4;
set(gca,'xlim',[5,22],'ylim',[-1*Fratio,1.25*Fratio]); hold on;
plt1=plot(x,y.*Fratio); hold on; plt1.Color='Black'; plt1.LineWidth=1.0;
x=CYF1W(:,1); y=CYF1W(:,2);
x=x-ones(size(x,1),1).*2.5; y=y.*1.25;
sc1=scatter(x,y.*Fratio); sc1.Marker='x'; sc1.MarkerEdgeColor='Red';
Vsize=2.032*0.7618*0.08; meshNOE=1036880+67040;
meshN=(meshNOE/Vsize)^(1/3);
x=CYF1ZX(:,1); y=CYF1ZX(:,2);
x=x-ones(size(x,1),1).*2.5; y=y.*1.25;
sc2=scatter(x,y.*Fratio); sc2.Marker='^'; sc2.MarkerEdgeColor='Blue';
load('cylinderF2.mat');
x=CYmat(:,1).*t; y=CYmat(:,2)./F; 
plt2=plot(x,y.*Fratio); hold on; plt2.Color='Black'; plt2.LineWidth=1.5; plt2.LineStyle='--';
xL1=xlabel('t/(h/g)^{1/2}'); hold on; xL1.FontSize=15;
yL1=ylabel('Fx/(\rhogH^2R)'); hold on; yL1.FontSize=15;
leg1=legend([plt1,plt2,sc1,sc2],["present (refined)","present (coarse)","Z.Xie Num","Yates Exp"]);
leg1.FontSize=12; Ymax=max(y.*Fratio);
saveas(gcf,"poltCYF1.eps",'epsc2');
saveas(gcf,"poltCYF1.png");




clearvars; clc; % proCYH1.m
CYH1data=importdata("ITFCHYC.1");
CYH2data=importdata("ITFCHYC.2");
CYHdata=[CYH1data;CYH2data];
CYHmat=zeros(size(CYHdata,1),4);
for ii=1:size(CYHdata,1)
    CYHmat(ii,2)=CYHdata(ii,11);
    CYHmat(ii,3)=CYHdata(ii,12);
    CYHmat(ii,4)=CYHdata(ii,13);
end
h=0.04;
CYHmat=CYHmat./h;
plot(CYHmat);


clearvars; clc; % proCYH2.m
CYH1data=importdata("p.rgh1.csv");
CYH2data=CYH1data.data(:,1:2);
CYH2data=CYH2data(3:size(CYH2data,1),:);
%plot(CYH2data); hold on;
CYH1data=CYH1data.textdata;
CYH1data=CYH1data(6:size(CYH1data,1),:);
CYH1mat=zeros(size(CYH1data,1),2);
for ii = 1:size(CYH1data,1)
    tmp=cell2mat(CYH1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    CYH1mat(ii,1)=tmp;
    tmp=cell2mat(CYH1data(ii,2));
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    CYH1mat(ii,2)=tmp;
end
CYH1mat=[CYH1mat,CYH2data];
[C2,I2]=max(CYH1mat(:,2));
for ii=1:(size(CYH1mat,1)-I2)
    CYH1mat(I2+ii,2)=CYH1mat(I2-ii,2);
end
[C3,I3]=max(CYH1mat(:,3));
for ii=1:(size(CYH1mat,1)-I3)
    CYH1mat(I3+ii,3)=CYH1mat(I3-ii,3);
end
[C5,I4]=max(CYH1mat(:,4));
for ii=1:(size(CYH1mat,1)-I4)
    CYH1mat(I4+ii,4)=CYH1mat(I4-ii,4);
end
p0=392.315;
CYH1mat(:,[2:4])=CYH1mat(:,[2:4])-ones(size(CYH1mat,1),3).*p0;
rho=1000; g=9.8; h=0.04; eta=rho*g*h; t=(g/h)^0.5;
CYH1mat(:,[2:4])=CYH1mat(:,[2:4])./eta;
CYH1mat(:,1)=CYH1mat(:,1).*t;
for ii=1:size(CYH1mat,1)
    for jj=1:size(CYH1mat,2)
        if(CYH1mat(ii,jj)<0.0)
            CYH1mat(ii,jj)=0.0;
        end
    end
end
save('CYH1.mat',"CYH1mat");

CYH1data=importdata("p.rgh2.csv");
CYH2data=CYH1data.data(:,1:2);
CYH2data=CYH2data(3:size(CYH2data,1),:);
%plot(CYH2data); hold on;
CYH1data=CYH1data.textdata;
CYH1data=CYH1data(6:size(CYH1data,1),:);
CYH1mat=zeros(size(CYH1data,1),2);
for ii = 1:size(CYH1data,1)
    tmp=cell2mat(CYH1data(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    CYH1mat(ii,1)=tmp;
    tmp=cell2mat(CYH1data(ii,2));
    tmp=convertCharsToStrings(tmp);
    tmp=str2double(tmp);
    CYH1mat(ii,2)=tmp;
end
CYH1mat=[CYH1mat,CYH2data];
[C2,I2]=max(CYH1mat(:,2));
for ii=1:(size(CYH1mat,1)-I2)
    CYH1mat(I2+ii,2)=CYH1mat(I2-ii,2);
end
[C3,I3]=max(CYH1mat(:,3));
for ii=1:(size(CYH1mat,1)-I3)
    CYH1mat(I3+ii,3)=CYH1mat(I3-ii,3);
end
[C5,I4]=max(CYH1mat(:,4));
for ii=1:(size(CYH1mat,1)-I4)
    CYH1mat(I4+ii,4)=CYH1mat(I4-ii,4);
end
p0=392.315;
CYH1mat(:,[2:4])=CYH1mat(:,[2:4])-ones(size(CYH1mat,1),3).*p0;
rho=1000; g=9.8; h=0.04; eta=rho*g*h; t=(g/h)^0.5;
CYH1mat(:,[2:4])=CYH1mat(:,[2:4])./eta;
CYH1mat(:,1)=CYH1mat(:,1).*t;
for ii=1:size(CYH1mat,1)
    for jj=1:size(CYH1mat,2)
        if(CYH1mat(ii,jj)<0.0)
            CYH1mat(ii,jj)=0.0;
        end
    end
end
CYH2mat=CYH1mat;
save('CYH2.mat',"CYH2mat");


clearvars; clc; % proCylinderF1.m
CYdata=importdata("cylinderF2.dat");
CYdata=CYdata(4:size(CYdata,1),1);
CYmat=zeros(size(CYdata,1),2);
for ii=1:size(CYdata,1)
    tmp=cell2mat(CYdata(ii,1));
    tmp=convertCharsToStrings(tmp);
    tmp=strsplit(tmp,'\t');
    tmp=tmp(1,1:2);
    CYmat(ii,1)=str2double(tmp(1,1));
    tmp=tmp(1,2);
    tmp=strsplit(tmp,' ');
    tmp=tmp(1,1);
    tmp=convertStringsToChars(tmp);
    tmp=tmp(1,3:size(tmp,2));
    tmp=convertCharsToStrings(tmp);
    CYmat(ii,2)=str2double(tmp);
end
save('cylinderF2.mat',"CYmat");

jmsyh
关注
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

jmsyh

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值