{"topic":{"id":1483,"title":"平板边界层流速分布实验","slug":"topic","categoryId":13,"createdAt":"2023-06-26 05:17:39.767446","bumpedAt":"2023-06-26 08:58:48.669043","views":274,"postsCount":3,"highestPostNumber":3,"likeCount":2,"closed":false,"pinned":false,"pinnedGlobally":false,"pinnedAt":null,"pinnedUntil":null,"isBanner":false,"excerpt":"平板边界层实验在热流体实验中大二上册，适用核工能动学生 MATLAB 源程序： %原文链接 https://blog.csdn.net/dasdkhaksd/article/details/115013265?ops_request_misc=&request_id=&biz_id=102&utm_term=matla%E7%94%BB%E8%BE%B9%E7%95%8C%E5%B1%82%E9%80%9F%E5…","deletedAt":null,"isAnonymous":false,"category":{"id":13,"name":"学习资料","slug":"res","color":"FFBB00"},"author":{"id":355,"username":"Humble","name":null,"avatarUrl":"/uploads/image/f6/27/f627037a25dbe7d5b0db8824e41369f47b9ec6d12fd0032cb0fe0e7cd8cb4f64.avif"},"tags":[]},"posts":[{"id":3479,"postNumber":1,"raw":"平板边界层实验在热流体实验中大二上册，适用核工能动学生\n\n[grid mode=\"tile\"]\n![image|558x500, 50%](/uploads/image/27/14/2714ce6affda4ac5de34332c9424cd0361e8f01a9c5fac0eae91b412aa131395.avif)![image|558x500, 50%](/uploads/image/a8/9a/a89a70419f89bacb2c8d8754201ae51c754739db6ae7fe7ea70a96b2f8729265.avif)\n[/grid]\nMATLAB 源程序：\n\n```\n%原文链接 https://blog.csdn.net/dasdkhaksd/article/details/115013265?ops_request_misc=&request_id=&biz_id=102&utm_term=matla%E7%94%BB%E8%BE%B9%E7%95%8C%E5%B1%82%E9%80%9F%E5%BA%A6%E5%88%86%E5%B8%83&utm_medium=distribute.pc_search_result.none-task-blog-2~all~sobaiduweb~default-0-115013265.142^v88^control_2,239^v2^insert_chatgpt&spm=1018.2226.3001.4187\n%本例对原文有删减，依照实验数据输入即可，注释部分有示例数据\nclear,clc;\nlw1=2;lw2=2;\nmks=15;\n\nRg = 287; %气体常数（J/(kg·K)）\ng = 9.8; %重力加速度（m/s^2）\nwaterdensity = 997.125; %水的密度（kg/m^3）\n\natmosphericPressure =input('请输入大气压：(Pa)') ;%95910\ntemperature =input('请输入室温：(°C)') ;%24.8\nKVC =input('请输入运动粘性系数（查表）：(10^-7m^2/s)') ;%156\n\nK = 1;%风洞系数\ngearFactor = 0.2;%压力计倾斜因子\n\nstaticPressure0 = input('请输入机械法流场静压：(Pa)');%254.1\nstaticPressure1 = input('请输入电测法流场静压：(Pa)');%248.3\n\n%对已知实验参数的输入\n%第一截面实验数据输入\nx10 = 150;\ny10 = input('请输入第一截面边界层内距离 (记得带方括号，行向量，下同）：(mm)');\nmicromanometerP10 = input('请输入第一截面用微压计测得的读数数据：(mmC2H5OH)');\ntransmitterP10 = input('请输入第一截面用压差变送器测得的读数数据：(Pa)'); \n% y10 = [0.45 0.95 1.45 1.95 2.45 2.95 3.45 3.95 4.45];\n% micromanometerP10 = [83.9 69.1 58.5\t49 40.3 31.5 26\t20.8 17];\n% transmitterP10 = [150.4 123.3 105.6\t86 68.2\t52.9 40.6 30.5 24.8];\n\n%第二截面实验数据输入\nx20 = 250;\ny20 = input('请输入第二截面边界层内距离的数据：(mm)');\nmicromanometerP20 = input('请输入第二截面用微压计测得的读数数据：(mmC2H5OH)');\ntransmitterP20 = input('请输入第二截面用压差变送器测得的读数数据：(Pa)');\n% y20 = [0.45\t0.95 1.45 1.95 2.45 3.45 4.45 5.45 ];\n% micromanometerP20 = [85.5 74.5 66.2\t59.9 52.5 40 29 20.5];\n% transmitterP20 = [154.6\t132.6 118 104.3\t90.2 66.7 45.9 30.2];\n\n%实验过程中通过计算出来的、不变的数据\nairDensity = 1.123; %空气密度\nstaticPressure =0.5*(staticPressure0+staticPressure1) ; %流场静压取平均\nx1 = x10 * 0.001;\ny1 = y10 * 0.001;\nx2 = x20 * 0.001;\ny2 = y20 * 0.001;\nvelocityInf1 = 20.08; %V无穷\nvelocityInf2 = 20.11; %V无穷\n\nmicromanometerP1 = atmosphericPressure - micromanometerP10 * gearFactor * waterdensity * g * 0.001;\ntransmitterP1 = atmosphericPressure - transmitterP10;\nmicromanometerP2 = atmosphericPressure - micromanometerP20 * gearFactor * waterdensity * g * 0.001;\ntransmitterP2 = atmosphericPressure - transmitterP20;\n\n%对第一截面数据进行处理\nfprintf('>>>>>>>>>>**第一截面数据计算**<<<<<<<<<<\\n');\nRe1 = (velocityInf1 * x1) / (KVC * 1e-7);\nfprintf('=>雷诺数 Re1= %.3f\\n',Re1);\n\nmicromanometerVelocity1 = zeros(1,length(micromanometerP1));\nfor i=1:length(micromanometerP10)\n micromanometerVelocity1(i) =(2*g*gearFactor*(130-micromanometerP10(i))/1000*(waterdensity-airDensity)/1.123)^0.5;\nend\n% fprintf('=>用微压计测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(micromanometerVelocity1);\ntransmitterVelocity1 = zeros(1,length(transmitterP1));\nfor i=1:length(transmitterP10)\n transmitterVelocity1(i) =(2*(staticPressure1-transmitterP10(i))/airDensity)^0.5;\nend\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(transmitterVelocity1);\nmicromanometerVelocityRatio1 = micromanometerVelocity1 / velocityInf1;\n% fprintf('=>用微压计测量的数据计算得到的边界层内的速度比：\\n');\n% disp(micromanometerVelocityRatio1);\ntransmitterVelocityRatio1 = transmitterVelocity1 / velocityInf1;\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的速度比：\\n');\n% disp(transmitterVelocityRatio1);\n\n%对第二截面数据进行处理\nfprintf('>>>>>>>>>>**第二截面数据计算**<<<<<<<<<<\\n');\nRe2 = (velocityInf2 * x2) / (KVC * 1e-7);\nfprintf('=>雷诺数 Re2= %.3f\\n',Re2);\n\nmicromanometerVelocity2 = zeros(1,length(micromanometerP20));\nfor i=1:length(micromanometerP20)\n micromanometerVelocity2(i) =(2*g*gearFactor*(130-micromanometerP20(i))/1000*(waterdensity-airDensity)/1.123)^0.5;\nend\n% fprintf('=>用微压计测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(micromanometerVelocity2);\ntransmitterVelocity2 = zeros(1,length(transmitterP20));\nfor i=1:length(transmitterP20)\n transmitterVelocity2(i) = (2*(staticPressure1-transmitterP20(i))/airDensity)^0.5;\nend\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(transmitterVelocity2);\nmicromanometerVelocityRatio2 = micromanometerVelocity2 / velocityInf2;\n% fprintf('=>用微压计测量的数据计算得到的边界层内的速度比：\\n');\n% disp(micromanometerVelocityRatio2);\ntransmitterVelocityRatio2 = transmitterVelocity2 / velocityInf2;\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的速度比：\\n');\n% disp(transmitterVelocityRatio2);\n\n%数据可视化\nfigure; %第一截面\nhold on;grid on;\nplot(micromanometerVelocityRatio1,y1,'o','markersize',mks); %微压计测量所得数据点\nplot(transmitterVelocityRatio1,y1,'p','markersize',mks); %压差变送器测量所得数据点\n\nmicromanometerT1 = 0.62 : 0.01 : 0.99;\nmicromanometerU1 = spline(micromanometerVelocityRatio1,y1,micromanometerT1);\nmicromanometerBLT1 = spline(micromanometerVelocityRatio1,y1,0.99); %Boundary layer thickness\nfprintf('=>用微压计测量得出的边界层厚度为：%.3f mm\\n',micromanometerBLT1 * 1000);\nplot(micromanometerT1,micromanometerU1,'-.','linewidth',lw1); %微压计测量拟合曲线。\n\ntransmitterT1 = 0.65 : 0.01 : 0.99;\ntransmitterU1 = spline(transmitterVelocityRatio1,y1,transmitterT1);\ntransmitterBLT1 = spline(transmitterVelocityRatio1,y1,0.99);\nfprintf('=>用压差变送器测量得出的边界层厚度为：%.3f mm\\n',transmitterBLT1 * 1000);\nplot(transmitterT1,transmitterU1,'--','linewidth',lw1); %压差变送器测量拟合曲线\n\nmicromanometerY1 = 0 : 0.00001 : transmitterBLT1;\nmicromanometerX1 = ((3 * micromanometerY1 - (micromanometerY1.^3) / (transmitterBLT1^2))...\n / (2 * transmitterBLT1));\nplot(micromanometerX1,micromanometerY1,'b-','linewidth',lw2); %层流理论曲线\n\ntransmitterY1 = 0 : 0.00001 : transmitterBLT1;\ntransmitterX1 = (transmitterY1 / transmitterBLT1).^(1/7);\nplot(transmitterX1,transmitterY1,'r-','linewidth',lw2); %紊流理论曲线\n\ntitle('第一截面边界层速度分布曲线');\nxlabel('速度比','FontSize',12,'FontWeight','bold');\nylabel('离平板壁面距离 y（m）','FontSize',12,'FontWeight','bold');\naxis tight;\nlegend('微压计测量所得数据点','压差变送器测量所得数据点','微压计测量拟合曲线','压差变送器测量拟合曲线',...\n '层流理论曲线','紊流理论曲线','location','northwest');\n\nfigure; %第二截面\nhold on;grid on;\nplot(micromanometerVelocityRatio2,y2,'o','markersize',mks); %微压计测量所得数据点\nplot(transmitterVelocityRatio2,y2,'p','markersize',mks); %压差变送器测量所得数据点\n\nmicromanometerT2 = 0.62 : 0.01 : 0.99;\nmicromanometerU2 = spline(micromanometerVelocityRatio2,y2,micromanometerT2);\nmicromanometerBLT2 = spline(micromanometerVelocityRatio2,y2,0.99);\nfprintf('=>用微压计测量得出的边界层厚度为：%.3f mm\\n',micromanometerBLT2 * 1000);\nplot(micromanometerT2,micromanometerU2,'-.','linewidth',lw1); %微压计测量拟合曲线\ntransmitterT2 = 0.62 : 0.01 : 0.99;\n\ntransmitterU2 = spline(transmitterVelocityRatio2,y2,transmitterT2);\ntransmitterBLT2 = spline(transmitterVelocityRatio2,y2,0.99);\nfprintf('=>用压差变送器测量得出的边界层厚度为：%.3f mm\\n',transmitterBLT2 * 1000);\nplot(transmitterT2,transmitterU2,'--','linewidth',lw1); %压差变送器测量拟合曲线\n\nmicromanometerY2 = 0:0.00001 : transmitterBLT2;\nmicromanometerX2 = ((3 * micromanometerY2 - (micromanometerY2.^3) / (transmitterBLT2^2)))...\n / (2*transmitterBLT2);\nplot(micromanometerX2,micromanometerY2,'b-','linewidth',lw2); %层流理论曲线\n\ntransmitterY2 = 0 : 0.00001 : transmitterBLT2;\ntransmitterX2 = (transmitterY2 / transmitterBLT2).^(1/7);\nplot(transmitterX2,transmitterY2,'r-','linewidth',lw2); %紊流理论曲线\n\ntitle('第二截面边界层速度分布曲线');\nxlabel('速度比','FontSize',12,'FontWeight','bold');\nylabel('离平板壁面距离 y（m）','FontSize',12,'FontWeight','bold');\naxis tight;\nlegend('微压计测量所得数据点','压差变送器测量所得数据点','微压计测量拟合曲线','压差变送器测量拟合曲线',...\n '层流理论曲线','紊流理论曲线','location','northwest');\n\n%流态分析\nfprintf('>>>>>>>>>>**流态分析**<<<<<<<<<<\\n');\nfprintf('一、对第一截面：\\n');\n\nsigmaLin1=5 * sqrt((KVC * 1e-7 * x1) / velocityInf1);display(sigmaLin1*1000);\nsigmaTu1=0.37 * x1 * ((KVC * 1e-7) / (velocityInf1 * x1))^0.2;display(sigmaTu1*1000);\n\nfprintf('二、对第二截面：\\n');\nsigmaLin2=5 * sqrt((KVC * 1e-7 * x2) / velocityInf2);display(sigmaLin2*1000);\nsigmaTu2=0.37 * x2 * ((KVC * 1e-7) / (velocityInf2 * x2))^0.2;display(sigmaTu2*1000);\n```","cooked":"

平板边界层实验在热流体实验中大二上册，适用核工能动学生

$\"image\"$ $\"image\"$

MATLAB 源程序：

%原文链接 https://blog.csdn.net/dasdkhaksd/article/details/115013265?ops_request_misc=&request_id=&biz_id=102&utm_term=matla%E7%94%BB%E8%BE%B9%E7%95%8C%E5%B1%82%E9%80%9F%E5%BA%A6%E5%88%86%E5%B8%83&utm_medium=distribute.pc_search_result.none-task-blog-2~all~sobaiduweb~default-0-115013265.142^v88^control_2,239^v2^insert_chatgpt&spm=1018.2226.3001.4187\n%本例对原文有删减，依照实验数据输入即可，注释部分有示例数据\nclear,clc;\nlw1=2;lw2=2;\nmks=15;\n\nRg = 287; %气体常数（J/(kg·K)）\ng = 9.8; %重力加速度（m/s^2）\nwaterdensity = 997.125; %水的密度（kg/m^3）\n\natmosphericPressure =input('请输入大气压：(Pa)') ;%95910\ntemperature =input('请输入室温：(°C)')  ;%24.8\nKVC =input('请输入运动粘性系数（查表）：(10^-7m^2/s)') ;%156\n\nK = 1;%风洞系数\ngearFactor = 0.2;%压力计倾斜因子\n\nstaticPressure0 = input('请输入机械法流场静压：(Pa)');%254.1\nstaticPressure1 = input('请输入电测法流场静压：(Pa)');%248.3\n\n%对已知实验参数的输入\n%第一截面实验数据输入\nx10 = 150;\ny10 = input('请输入第一截面边界层内距离 (记得带方括号，行向量，下同）：(mm)');\nmicromanometerP10 = input('请输入第一截面用微压计测得的读数数据：(mmC2H5OH)');\ntransmitterP10 = input('请输入第一截面用压差变送器测得的读数数据：(Pa)'); \n% y10 = [0.45 0.95 1.45 1.95 2.45 2.95 3.45 3.95 4.45];\n% micromanometerP10 = [83.9 69.1 58.5\t49 40.3 31.5 26\t20.8 17];\n% transmitterP10 = [150.4 123.3 105.6\t86 68.2\t52.9 40.6 30.5 24.8];\n\n%第二截面实验数据输入\nx20 = 250;\ny20 = input('请输入第二截面边界层内距离的数据：(mm)');\nmicromanometerP20 = input('请输入第二截面用微压计测得的读数数据：(mmC2H5OH)');\ntransmitterP20 = input('请输入第二截面用压差变送器测得的读数数据：(Pa)');\n% y20 = [0.45\t0.95 1.45 1.95 2.45 3.45 4.45 5.45 ];\n% micromanometerP20 = [85.5 74.5 66.2\t59.9 52.5 40 29 20.5];\n% transmitterP20 = [154.6\t132.6 118 104.3\t90.2 66.7 45.9 30.2];\n\n%实验过程中通过计算出来的、不变的数据\nairDensity = 1.123; %空气密度\nstaticPressure =0.5*(staticPressure0+staticPressure1) ; %流场静压取平均\nx1 = x10 * 0.001;\ny1 = y10 * 0.001;\nx2 = x20 * 0.001;\ny2 = y20 * 0.001;\nvelocityInf1 = 20.08;    %V无穷\nvelocityInf2 = 20.11;    %V无穷\n\nmicromanometerP1 = atmosphericPressure - micromanometerP10 * gearFactor * waterdensity * g * 0.001;\ntransmitterP1 = atmosphericPressure - transmitterP10;\nmicromanometerP2 = atmosphericPressure - micromanometerP20 * gearFactor * waterdensity * g * 0.001;\ntransmitterP2 = atmosphericPressure - transmitterP20;\n\n%对第一截面数据进行处理\nfprintf('>>>>>>>>>>**第一截面数据计算**<<<<<<<<<<\\n');\nRe1 = (velocityInf1 * x1) / (KVC * 1e-7);\nfprintf('=>雷诺数 Re1= %.3f\\n',Re1);\n\nmicromanometerVelocity1 = zeros(1,length(micromanometerP1));\nfor i=1:length(micromanometerP10)\n    micromanometerVelocity1(i) =(2*g*gearFactor*(130-micromanometerP10(i))/1000*(waterdensity-airDensity)/1.123)^0.5;\nend\n% fprintf('=>用微压计测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(micromanometerVelocity1);\ntransmitterVelocity1 = zeros(1,length(transmitterP1));\nfor i=1:length(transmitterP10)\n    transmitterVelocity1(i) =(2*(staticPressure1-transmitterP10(i))/airDensity)^0.5;\nend\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(transmitterVelocity1);\nmicromanometerVelocityRatio1 = micromanometerVelocity1 / velocityInf1;\n% fprintf('=>用微压计测量的数据计算得到的边界层内的速度比：\\n');\n% disp(micromanometerVelocityRatio1);\ntransmitterVelocityRatio1 = transmitterVelocity1 / velocityInf1;\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的速度比：\\n');\n% disp(transmitterVelocityRatio1);\n\n%对第二截面数据进行处理\nfprintf('>>>>>>>>>>**第二截面数据计算**<<<<<<<<<<\\n');\nRe2 = (velocityInf2 * x2) / (KVC * 1e-7);\nfprintf('=>雷诺数 Re2= %.3f\\n',Re2);\n\nmicromanometerVelocity2 = zeros(1,length(micromanometerP20));\nfor i=1:length(micromanometerP20)\n    micromanometerVelocity2(i) =(2*g*gearFactor*(130-micromanometerP20(i))/1000*(waterdensity-airDensity)/1.123)^0.5;\nend\n% fprintf('=>用微压计测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(micromanometerVelocity2);\ntransmitterVelocity2 = zeros(1,length(transmitterP20));\nfor i=1:length(transmitterP20)\n    transmitterVelocity2(i) = (2*(staticPressure1-transmitterP20(i))/airDensity)^0.5;\nend\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的流速：（m/s）\\n');\n% disp(transmitterVelocity2);\nmicromanometerVelocityRatio2 = micromanometerVelocity2 / velocityInf2;\n% fprintf('=>用微压计测量的数据计算得到的边界层内的速度比：\\n');\n% disp(micromanometerVelocityRatio2);\ntransmitterVelocityRatio2 = transmitterVelocity2 / velocityInf2;\n% fprintf('=>用压差变送器测量的数据计算得到的边界层内的速度比：\\n');\n% disp(transmitterVelocityRatio2);\n\n%数据可视化\nfigure; %第一截面\nhold on;grid on;\nplot(micromanometerVelocityRatio1,y1,'o','markersize',mks); %微压计测量所得数据点\nplot(transmitterVelocityRatio1,y1,'p','markersize',mks); %压差变送器测量所得数据点\n\nmicromanometerT1 = 0.62 : 0.01 : 0.99;\nmicromanometerU1 = spline(micromanometerVelocityRatio1,y1,micromanometerT1);\nmicromanometerBLT1 = spline(micromanometerVelocityRatio1,y1,0.99); %Boundary layer thickness\nfprintf('=>用微压计测量得出的边界层厚度为：%.3f mm\\n',micromanometerBLT1 * 1000);\nplot(micromanometerT1,micromanometerU1,'-.','linewidth',lw1); %微压计测量拟合曲线。\n\ntransmitterT1 = 0.65 : 0.01 : 0.99;\ntransmitterU1 = spline(transmitterVelocityRatio1,y1,transmitterT1);\ntransmitterBLT1 = spline(transmitterVelocityRatio1,y1,0.99);\nfprintf('=>用压差变送器测量得出的边界层厚度为：%.3f mm\\n',transmitterBLT1 * 1000);\nplot(transmitterT1,transmitterU1,'--','linewidth',lw1); %压差变送器测量拟合曲线\n\nmicromanometerY1 = 0 : 0.00001 : transmitterBLT1;\nmicromanometerX1 = ((3 * micromanometerY1 - (micromanometerY1.^3) / (transmitterBLT1^2))...\n    / (2 * transmitterBLT1));\nplot(micromanometerX1,micromanometerY1,'b-','linewidth',lw2); %层流理论曲线\n\ntransmitterY1 = 0 : 0.00001 : transmitterBLT1;\ntransmitterX1 = (transmitterY1 / transmitterBLT1).^(1/7);\nplot(transmitterX1,transmitterY1,'r-','linewidth',lw2); %紊流理论曲线\n\ntitle('第一截面边界层速度分布曲线');\nxlabel('速度比','FontSize',12,'FontWeight','bold');\nylabel('离平板壁面距离 y（m）','FontSize',12,'FontWeight','bold');\naxis tight;\nlegend('微压计测量所得数据点','压差变送器测量所得数据点','微压计测量拟合曲线','压差变送器测量拟合曲线',...\n    '层流理论曲线','紊流理论曲线','location','northwest');\n\nfigure; %第二截面\nhold on;grid on;\nplot(micromanometerVelocityRatio2,y2,'o','markersize',mks); %微压计测量所得数据点\nplot(transmitterVelocityRatio2,y2,'p','markersize',mks); %压差变送器测量所得数据点\n\nmicromanometerT2 = 0.62 : 0.01 : 0.99;\nmicromanometerU2 = spline(micromanometerVelocityRatio2,y2,micromanometerT2);\nmicromanometerBLT2 = spline(micromanometerVelocityRatio2,y2,0.99);\nfprintf('=>用微压计测量得出的边界层厚度为：%.3f mm\\n',micromanometerBLT2 * 1000);\nplot(micromanometerT2,micromanometerU2,'-.','linewidth',lw1); %微压计测量拟合曲线\ntransmitterT2 = 0.62 : 0.01 : 0.99;\n\ntransmitterU2 = spline(transmitterVelocityRatio2,y2,transmitterT2);\ntransmitterBLT2 = spline(transmitterVelocityRatio2,y2,0.99);\nfprintf('=>用压差变送器测量得出的边界层厚度为：%.3f mm\\n',transmitterBLT2 * 1000);\nplot(transmitterT2,transmitterU2,'--','linewidth',lw1); %压差变送器测量拟合曲线\n\nmicromanometerY2 = 0:0.00001 : transmitterBLT2;\nmicromanometerX2 = ((3 * micromanometerY2 - (micromanometerY2.^3) / (transmitterBLT2^2)))...\n    / (2*transmitterBLT2);\nplot(micromanometerX2,micromanometerY2,'b-','linewidth',lw2); %层流理论曲线\n\ntransmitterY2 = 0 : 0.00001 : transmitterBLT2;\ntransmitterX2 = (transmitterY2 / transmitterBLT2).^(1/7);\nplot(transmitterX2,transmitterY2,'r-','linewidth',lw2); %紊流理论曲线\n\ntitle('第二截面边界层速度分布曲线');\nxlabel('速度比','FontSize',12,'FontWeight','bold');\nylabel('离平板壁面距离 y（m）','FontSize',12,'FontWeight','bold');\naxis tight;\nlegend('微压计测量所得数据点','压差变送器测量所得数据点','微压计测量拟合曲线','压差变送器测量拟合曲线',...\n    '层流理论曲线','紊流理论曲线','location','northwest');\n\n%流态分析\nfprintf('>>>>>>>>>>**流态分析**<<<<<<<<<<\\n');\nfprintf('一、对第一截面：\\n');\n\nsigmaLin1=5 * sqrt((KVC * 1e-7 * x1) / velocityInf1);display(sigmaLin1*1000);\nsigmaTu1=0.37 * x1 * ((KVC * 1e-7) / (velocityInf1 * x1))^0.2;display(sigmaTu1*1000);\n\nfprintf('二、对第二截面：\\n');\nsigmaLin2=5 * sqrt((KVC * 1e-7 * x2) / velocityInf2);display(sigmaLin2*1000);\nsigmaTu2=0.37 * x2 * ((KVC * 1e-7) / (velocityInf2 * x2))^0.2;display(sigmaTu2*1000);\n

\n","createdAt":"2023-06-26 05:17:40.295675","updatedAt":"2023-06-26 05:17:40.295675","replyToPostNumber":null,"replyCount":0,"likeCount":2,"wordCount":1503,"liked":false,"bookmarked":false,"polls":[],"reactions":[{"reaction":"❤️","count":2,"users":[{"id":309,"username":"Enriched_Uranium","name":"一小块浓缩铀"},{"id":1047,"username":"anonymous_coward_old","name":"Anonymous Coward Old"}],"hasReacted":false}],"wiki":false,"version":1,"publicVersion":1,"deletedAt":null,"isAnonymous":false,"user":{"id":355,"username":"Humble","name":null,"avatarUrl":"/uploads/image/f6/27/f627037a25dbe7d5b0db8824e41369f47b9ec6d12fd0032cb0fe0e7cd8cb4f64.avif","admin":false,"moderator":false,"signatureRaw":null,"signatureCooked":null}},{"id":3482,"postNumber":2,"raw":"非常好工作 :kissing_closed_eyes:","cooked":"

非常好工作 😚

\n","createdAt":"2023-06-26 07:35:47.812859","updatedAt":"2023-06-26 07:35:47.812859","replyToPostNumber":null,"replyCount":1,"likeCount":0,"wordCount":6,"liked":false,"bookmarked":false,"polls":[],"reactions":[],"wiki":false,"version":1,"publicVersion":1,"deletedAt":null,"isAnonymous":false,"user":{"id":309,"username":"Enriched_Uranium","name":"一小块浓缩铀","avatarUrl":"/uploads/image/af/ec/afec38c6f15759c4eea7462b7de361c2dd8fe60fa21873ace11085a7e94918eb.avif","admin":false,"moderator":false,"signatureRaw":null,"signatureCooked":null}},{"id":3485,"postNumber":3,"raw":":heart_eyes: :heart_eyes: :heart_eyes: :heart_eyes:","cooked":"

😍 😍 😍 😍

\n","createdAt":"2023-06-26 08:58:48.669043","updatedAt":"2023-06-26 08:58:48.669043","replyToPostNumber":2,"replyCount":0,"likeCount":0,"wordCount":4,"liked":false,"bookmarked":false,"polls":[],"reactions":[],"wiki":false,"version":1,"publicVersion":1,"deletedAt":null,"isAnonymous":false,"user":{"id":355,"username":"Humble","name":null,"avatarUrl":"/uploads/image/f6/27/f627037a25dbe7d5b0db8824e41369f47b9ec6d12fd0032cb0fe0e7cd8cb4f64.avif","admin":false,"moderator":false,"signatureRaw":null,"signatureCooked":null}}],"hasMore":false,"lastReadPostNumber":1,"hasUrlPostNumber":true}