diff --git a/ruoyi-geek-vue3 b/ruoyi-geek-vue3 index 26446e5..ab0a15c 160000 --- a/ruoyi-geek-vue3 +++ b/ruoyi-geek-vue3 @@ -1 +1 @@ -Subproject commit 26446e5304424bcfe69036a30867cc428e5a929c +Subproject commit ab0a15c456bcaeff87f46d56d390162de82d7e6f diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/DpFlowCalc.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/DpFlowCalc.java new file mode 100644 index 0000000..84218ff --- /dev/null +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/DpFlowCalc.java @@ -0,0 +1,1130 @@ +package com.ruoyi.ngCalTools.controller; + +import com.ruoyi.ngCalTools.model.FlowProps; +import com.ruoyi.ngCalTools.model.GasProps; + +import javax.swing.*; +import java.util.Optional; + +import static com.ruoyi.ngCalTools.controller.FlowController.FlowConvert_BaseToWork; +import static com.ruoyi.ngCalTools.controller.FlowController.YaLiSunShi; + +public class DpFlowCalc { + /** + 标准孔板流量计算 + */ + public void OFlowCal(GasProps gasProps, FlowProps flowProps){ + + flowProps.setdOrificeD(flowProps.getdOrificeD() * (1 + 0.000001 * (flowProps.getdOrificeMaterial()) * (flowProps.getdTf() - 293.15))); + flowProps.setdPipeD(flowProps.getdPipeD() * (1 + 0.000001 * (flowProps.getdPipeMaterial()) * (flowProps.getdTf() - 293.15))); + flowProps.setdBeta(flowProps.getdOrificeD() / flowProps.getdPipeD()); +// flowProps.setdBeta(0.5972); + + //求渐近速度系数 E + flowProps.setdE(1 / Math.pow((1 - Math.pow(flowProps.getdBeta(), 4)), 0.5)); +// flowProps.setdE(1.0615); + + //求相对密度系数 FG + // LiuTiType = 1 + flowProps.setdFG(Math.pow((1 / gasProps.dRD_Real), 0.5)); +// flowProps.setdFG(1.2531); + //求流动温度系数 'FT + flowProps.setdFT(Math.pow((293.15 / flowProps.getdTf()), 0.5)); +// flowProps.setdFT(1.0086); + //求等熵指数????????????????????????????????? + + flowProps.setdKappa(gasProps.dKappa); +// flowProps.setdKappa(1.357); + + //求动力粘度 dlnd + flowProps.setdDViscosity(Dlndjs(flowProps.getdPf()/1e6, flowProps.getdTf())); +// flowProps.setdDViscosity(0.01096); + + //求可膨胀系数 + flowProps.setdDExpCoefficient(KePenZhang_JiSuan(flowProps.getdPf() , flowProps.getdDp(), flowProps.getdBeta(), gasProps.getdKappa(), flowProps.getdCoreType(), 0)); +// flowProps.setdDExpCoefficient(0.9977); + +// gasProps.dFpv=1.0195; +// gasProps.dHhvv=39.944; +// gasProps.dRD_Real=0.6368; + + //迭代计算流量和流出系数 + + double conQvA = 0.0000031795 * (1530000.0D * gasProps.dRD_Real / (flowProps.getdDViscosity() * flowProps.getdOrificeD())) * flowProps.getdE() * Math.pow(flowProps.getdOrificeD(), 2) * gasProps.dFpv * flowProps.getdFG() * flowProps.getdDExpCoefficient() * flowProps.getdFT() * Math.sqrt(flowProps.getdPf() * flowProps.getdDp() / 1e6); + double[] XQv = new double[4]; + double[] CQv = new double[4]; + double[] dQv = new double[4]; + double Qn = 0; + int n ; + XQv[0] = 1000000.0F; + boolean xhFlag=true; + _100000: + + while(xhFlag) { + for (n = 1; n <= 2; n++) { + CQv[n] = C_JiSuan(flowProps.getdOrificeD(), flowProps.getdBeta(), XQv[n - 1], conQvA, flowProps.getdPtmode(), flowProps.getdCoreType(), 0); + XQv[n] = conQvA * CQv[n]; + dQv[n] = XQv[n] - XQv[n - 1]; + } + if (XQv[2] == XQv[1] || dQv[2] == dQv[1]) { + Qn = conQvA * flowProps.getdDViscosity() * flowProps.getdOrificeD() * CQv[2] / (1530000.0D * gasProps.dRD_Real); + //管道雷诺数 + flowProps.setdRnPipe(XQv[2]); + //流出系数 + flowProps.setdCd(CQv[2]); + + } + if (Math.abs((conQvA - XQv[2] / CQv[1]) / conQvA) > 0.00000000000000005) { + XQv[0] = XQv[n - 1] - dQv[n - 1] * ((XQv[n - 1] - XQv[n - 2]) / (dQv[n - 1] - dQv[n - 2])); +//C# TO JAVA CONVERTER TODO TASK: There is no 'goto' in Java: + continue _100000; + } + else + { + xhFlag=false; + } + } + + //孔板锐利度系数Bk + flowProps.setdOrificeSharpness(1); + if (flowProps.getdCoreType() == 0) + { + flowProps.setdBk((flowProps.getdOrificeSharpness() == 0) ? (BkTable(flowProps.getdOrificeRk(), flowProps.getdOrificeD(), 1)) : (flowProps.getdOrificeSharpness())); + } + else + { + flowProps.setdBk(1); + } + + + //管道粗糙度系数 Gme + flowProps.setdRoughNessPipe(CcdXsjs(flowProps.getdPipeType(), flowProps.getdPipeD(), flowProps.getdBeta(), flowProps.getdRnPipe())); + //修正后的流出系数 + flowProps.setdCd(flowProps.getdCd() * flowProps.getdBk() * flowProps.getdRoughNessPipe()); +// flowProps.setdCd(0.6039); + //标况体积流量 m³、s + flowProps.setdVFlowb(Qn * flowProps.getdBk() * flowProps.getdRoughNessPipe()); + //工况体积流量 + flowProps.setdVFlowf(FlowConvert_BaseToWork(flowProps, gasProps)); + //标况质量流量 + flowProps.setdMFlowb(flowProps.getdVFlowb() * gasProps.dRhob); + //标况能量流量 + flowProps.setdEFlowb(flowProps.getdVFlowb() * gasProps.dHhvv); + //管道内天然气流速 + flowProps.setdVelocityFlow(flowProps.getdVFlowf() / (3.1415926 * Math.pow((flowProps.getdPipeD() / 2000), 2))); + //压力损失 + flowProps.setdPressLost(YaLiSunShi(flowProps.getdCd(), flowProps.getdBeta(), flowProps.getdDp(), flowProps.getdCoreType())); + } + + /** + * 计算管道绝对粗糙度 K (GB/T 21446-2008 附录C) + * @param flowProps + * @return 粗糙度修正系数 K (保留4位小数) + */ + public static double calculateK(FlowProps flowProps){ + double Jdccd; + switch (flowProps.getdPipeType()) { + + case 0: + Jdccd = 0.029F; + break; + case 1: + case 2: + case 3: + Jdccd = 0.075F; + break; + case 4: + Jdccd = 0.1F; + break; + case 5: + Jdccd = 0.15F; + break; + case 6: + Jdccd = 1; + break; + case 7: + Jdccd = 2.1F; + break; + case 8: + Jdccd = 0.04F; + break; + case 9: + Jdccd = 0.15F; + break; + case 10: + Jdccd = 0.13F; + break; + case 11: + Jdccd = 0.25F; + break; + default: + // 处理未知类型（可选） + throw new IllegalArgumentException("未知的管道类型: "); + } + return Jdccd; + } + + /** + * 计算管道粗糙度修正系数 G_me (GB/T 21446-2008 附录C) + * @param D_pipe 管道内径 (mm) + * @param K 绝对粗糙度 (mm) + * @param C 未修正的流出系数 + * @return 粗糙度修正系数 G_me (保留4位小数) + */ + public static double calculateRoughnessFactor(double D_pipe, double K, double C) { + // 计算相对粗糙度 K/D + double K_over_D = K / D_pipe; + + // 判断是否需要修正 + if (K_over_D <= 0.0004) { + return 1.0000; + } + + // 计算修正项 + double term = (K_over_D * 1e6) - 400; // 转换为无量纲项 + if (term < 0) { + throw new IllegalArgumentException("K/D 超出修正公式适用范围"); + } + + double G_me = 1 + (0.011 / C) * Math.sqrt(term); + return FormatUtil.format(G_me, Optional.of(4)); // 保留四位小数 + } + + /** + * 热膨胀修正（GB/T 21446-2008 第6.2.1节） + * @param flowProps 流量参数 + */ + public static void thermalExpansionCorrection(FlowProps flowProps) { + // 孔板孔径修正 + double dOrificeMaterial = flowProps.getdOrificeMaterial(); // 材料膨胀系数 (e.g. 16e-6/℃) + double dOrificeD = flowProps.getdOrificeD(); // 初始孔径(mm) + double dOrificeD_corrected = dOrificeD * (1 + 0.000001*dOrificeMaterial * (flowProps.getdTf() - 293.15)); // 20℃为参考温度 + flowProps.setdOrificeD(FormatUtil.format(dOrificeD_corrected, Optional.of(2))); + + // 管道内径修正 + double dPipeMaterial = flowProps.getdPipeMaterial(); + double dPipeD = flowProps.getdPipeD(); + double dPipeD_corrected = dPipeD * (1 + 0.000001*dPipeMaterial * (flowProps.getdTf() - 293.15)); + flowProps.setdPipeD(FormatUtil.format(dPipeD_corrected,Optional.of(2))); + + // 更新直径比β + flowProps.setdBeta( FormatUtil.format(flowProps.getdOrificeD() / flowProps.getdPipeD(),Optional.of(4))); + } + + /** + * 计算渐近速度系数E（GB/T 21446-2008 式(8)） + */ + public static void calculateE(FlowProps flowProps) { + double beta = flowProps.getdBeta(); + flowProps.setdE( FormatUtil.format(1 / Math.sqrt(1 - Math.pow(beta, 4)),Optional.of(4))); + } + + /** + * 计算相对密度系数FG（GB/T 21446-2008 式(9)） + */ + public static void calculateFG(FlowProps flowProps,GasProps gasProps) { + gasProps.setdRD_Real(gasProps.getdRD_Real()); + flowProps.setdFG( FormatUtil.format(1 / Math.sqrt(gasProps.dRD_Real),Optional.of(4))); + } + + /** + * 计算流动温度系数FT（GB/T 21446-2008 式(10)） + */ + public static void calculateFT(FlowProps flowProps) { + double Tb = flowProps.getdTb_M() ; // 参比温度转K + double Tf = flowProps.getdTf() ; // 工况温度转K + flowProps.setdFT( FormatUtil.format(Math.sqrt(Tb / Tf),Optional.of(4))); + } + + /** + * 计算可膨胀系数ε（GB/T 21446-2008 式(11)） + */ + public static void calculateEpsilon(FlowProps flowProps, GasProps gasProps) { + double P1 = flowProps.getdPf(); // 上游绝对压力(Pa) + double deltaP = flowProps.getdDp(); // 差压(Pa) + double beta = flowProps.getdBeta(); + double kappa = flowProps.getdKappa(); + + double tau = FormatUtil.format((P1 - deltaP) / P1, Optional.of(4)); // 压力比 + double epsilon=1 - (0.351 + 0.256 * Math.pow(beta, 4) + 0.93 * Math.pow(beta, 8)) * (1 - Math.pow(tau, 1/kappa)); + + flowProps.setdDExpCoefficient( FormatUtil.format(epsilon,Optional.of(4))); + } + /** + * 计算等熵指数κ（GB/T 21446-2008 推荐方法） + * @param gasProps 气体参数（需包含等熵指数或热容比） + * @return 等熵指数κ + */ + public static double calculateKappa(GasProps gasProps) { + // 若已通过AGA 8或其他方法计算，直接返回 +// if (gasProps.getdKappa() > 0) { +// return gasProps.getdKappa(); +// } + + // 近似公式：基于理想气体比热比和压缩因子修正 + double gamma = 1.3; // 天然气典型比热比（Cp/Cv≈1.3） + double Z = gasProps.getdZf(); // 工况压缩因子 + + // 修正公式（经验关系） + double kappa = gamma / (1 - (gamma - 1) * (1 / Z - 1)); + return FormatUtil.format(kappa,Optional.of(4)); + } + + + + /** + * 迭代计算流量（GB/T 21446-2008 第7章） + * @param flowProps 流量参数 + * @param gasProps 气体参数 + */ + public static void iterativeFlowCalculation(FlowProps flowProps, GasProps gasProps) { + // 0. 单位转换 + double D = flowProps.getdPipeD() / 1000.0; // 管道内径(m) + double d = flowProps.getdOrificeD() / 1000.0; // 孔板孔径(m) + double beta = flowProps.getdBeta(); + double P1 = flowProps.getdPf(); // 绝对压力(Pa) + double deltaP = flowProps.getdDp(); // 差压(Pa) + double Tf = flowProps.getdTf() ; // 工况温度(K) + + flowProps.setdPb_M(gasProps.getdPb()); + flowProps.setdTb_M(gasProps.getdTb()); + flowProps.setdFpv(FormatUtil.format(gasProps.dFpv,Optional.of(4))); + + // 1. 计算中间参数 + calculateE(flowProps); + calculateFG(flowProps,gasProps); + calculateFT(flowProps); + + flowProps.setdKappa( calculateKappa(gasProps)); +// flowProps.setdDViscosity( calculateDynamicViscosity(gasProps)); + flowProps.setdDViscosity(FormatUtil.format( Dlndjs(flowProps.getdPf()/1e6,flowProps.getdTf()),Optional.of(5))); + calculateEpsilon(flowProps, gasProps); + + flowProps.setdKappa(1.357); +// flowProps.setdBeta(0.5792); +// flowProps.setdFpv(1.0195); +// flowProps.setdFG(1.2531); +// flowProps.setdE(1.0615); +// flowProps.setdDExpCoefficient(0.9977); +// flowProps.setdFT(1.0086); + + + // 2. 初始雷诺数估算（假设初始C=0.6） + double C_initial = 0.6; + double Qf_initial = (C_initial * flowProps.getdE() * flowProps.getdDExpCoefficient() * Math.PI * Math.pow(d, 2) / 4) + * Math.sqrt(2 * deltaP / (gasProps.dRhof * (1 - Math.pow(beta, 4)))); + flowProps.setdVFlowf(Qf_initial); // 初始工况流量(m³/s) + + // 3. 迭代参数 + double tolerance = 1e-6; + int maxIter = 100; + double currentC = C_initial; + double currentReD = calculateReD(Qf_initial, D, gasProps.dRhof, flowProps.getdDViscosity()); + int iter = 0; + double prevC=0; + // 4. 迭代循环 + do { + prevC = currentC; + + // 4.1 计算流出系数C（GB/T 21446-2008 附录A） + currentC =FormatUtil.format( calculateCd(beta, currentReD, flowProps.getdPipeD(), flowProps.getdPtmode()),Optional.of(4)); + + // 4.2 更新流量 + double Qf =FormatUtil.format ((currentC * flowProps.getdDExpCoefficient() * Math.PI * Math.pow(d, 2) / 4) + * Math.sqrt(2 * deltaP / (gasProps.dRhof * (1 - Math.pow(beta, 4)))),Optional.of(4)); + flowProps.setdVFlowf(Qf); + + // 4.3 更新雷诺数 + currentReD =FormatUtil.format( calculateReD(Qf, D, gasProps.dRhof, flowProps.getdDViscosity()), Optional.of(2)); + + iter++; + if (iter > maxIter) { + throw new RuntimeException("迭代未收敛，超过最大迭代次数！"); + } + } while (Math.abs(currentC - prevC) / currentC > tolerance); + // 5. 保存最终结果 + + // 在迭代计算流出系数后，添加粗糙度修正 + + double K = calculateK(flowProps); // 根据实际管道类型选择 + double G_me = calculateRoughnessFactor(flowProps.getdPipeD(), K, currentC); + double C_corrected = FormatUtil.format(currentC * G_me, Optional.of(4)); + + flowProps.setdCd(C_corrected); + flowProps.setdRoughNessPipe(G_me); + flowProps.setdRnPipe(currentReD); + + + + + // 6. 计算标况流量（GB/T 21446-2008 式(1)） + double Qn = flowProps.getdVFlowf() * (flowProps.getdFpv()*flowProps.getdFpv() * P1 / flowProps.getdPb_M()) + * (flowProps.getdTb_M()) / Tf; + flowProps.setdVFlowb(FormatUtil.format(Qn,Optional.of(4))); + + //标况质量流量 + flowProps.setdMFlowb(FormatUtil.format(flowProps.getdVFlowb() * gasProps.dRhob,Optional.of(4))); + //标况能量流量 + flowProps.setdEFlowb(FormatUtil.format(flowProps.getdVFlowb() * gasProps.dHhvv,Optional.of(4))); + //管道内天然气流速 + flowProps.setdVelocityFlow(FormatUtil.format(flowProps.getdVFlowf() / (Math.PI * Math.pow((flowProps.getdPipeD() / 2000), 2)),Optional.of(3))); + //压力损失 + flowProps.setdPressLost(FormatUtil.format(YaLiSunShi(flowProps.getdCd(), flowProps.getdBeta(), flowProps.getdDp(), flowProps.getdCoreType()),Optional.of(2))); + + } + + + + + + /** + * 计算天然气动力粘度（Sutherland公式简化版，适用于低压天然气） + * @param gasProps 气体参数（需包含温度、密度、组分） + * @return 动力粘度 (Pa·s) + */ + public static double calculateDynamicViscosity(GasProps gasProps) { + double T = gasProps.getdTf() ; // 工况温度(K) + double rho = gasProps.getdRhof(); // 工况密度(kg/m³) + double Rd = gasProps.getdRD_Real(); // 真实相对密度 + + // Sutherland常数（天然气典型值，C=120） + double C = 120.0; // Sutherland常数 + double mu0 = 1.0e-5; // 参考粘度 (Pa·s，@293.15K) + + // 计算粘度（简化公式，忽略压力修正） + double mu = mu0 * Math.pow(T / 293.15, 0.7) * (293.15 + C) / (T + C); + + // 密度修正（经验公式） + mu *= 1 + 0.001 * (rho - 1.2); // 密度对粘度的粗略修正 + + return FormatUtil.format( mu*1000,Optional.of(5)); + } + /** + * 计算雷诺数ReD（GB/T 21446-2008 式(5)） + */ + public static double calculateReD(double Qf, double D, double rho, double mu) { + return FormatUtil.format( (4 * Qf * rho*1000) / (Math.PI * D * mu),Optional.of(2)); // Qf: m³/s, D: m, rho: kg/m³, mu: Pa·s + } + + /** + * 计算流出系数C（GB/T 21446-2008 附录A） + */ + public static double calculateCd(double beta, double ReD, double D_mm, int ptMode) { + double L1, L2; + // 根据取压方式确定L1/L2（角接取压） + switch (ptMode) { + case 1: // 角接取压 + L1 = L2 = 0; // D单位为mm + break; + case 0: // 法兰取压 + L1 = L2 = 25.4 / D_mm; + break; + case 2: // D-D/2取压 + L1 = 1.0; + L2 = 0.47; + break; + default: + throw new IllegalArgumentException("不支持的取压方式"); + } + + double term1 = 0.5961 + 0.0261 * Math.pow(beta, 2) - 0.216 * Math.pow(beta, 8); + double term2 = 0.000521 * Math.pow(1e6 * beta / ReD, 0.7); + double A = Math.pow(19000 * beta / ReD, 0.8); + double term3 = (0.0188 + 0.0063 * A) * Math.pow(beta, 3.5) * Math.pow(1e6 / ReD, 0.3); + double term4 = (0.043 + 0.08 * Math.exp(-10 * L1) - 0.123 * Math.exp(-7 * L1)) + * (1 - 0.11 * A) * Math.pow(beta, 4) / (1 - Math.pow(beta, 4)); + double term5 = -0.031 * (2 * L2 / (1 - beta) - 0.8 * Math.pow(2 * L2 / (1 - beta), 1.1)) + * Math.pow(beta, 1.3); + + double Cd = term1 + term2 + term3 + term4 + term5; + + // 孔径<71.12mm修正 + if (D_mm < 71.12) { + Cd += 0.011 * (0.75 - beta) * (2.8 - D_mm / 25.4); + } + return FormatUtil.format(Cd,Optional.of(4)); + } + //查表计算粘度μ + public static double Dlndjs(double tempP_jy, double tempT) + { double[][] Dlndjs_Dlnd_Data = new double[8][11]; + double[] Dlndjs_Dlnd_T = new double[8]; + double[] Dlndjs_Dlnd_P = new double[11]; + double s1 ; + double s2 ; + double ky ; + double kx ; + int i ; + int m = 0; + int n = 0; + //On Error Resume Next VBConversions Warning: On Error Resume Next not supported in C# + Dlndjs_Dlnd_T[0] = -15 + 273.15; + Dlndjs_Dlnd_T[1] = 0 + 273.15; + Dlndjs_Dlnd_T[2] = 15 + 273.15; + Dlndjs_Dlnd_T[3] = 30 + 273.15; + Dlndjs_Dlnd_T[4] = 45 + 273.15; + Dlndjs_Dlnd_T[5] = 60 + 273.15; + Dlndjs_Dlnd_T[6] = 75 + 273.15; + Dlndjs_Dlnd_T[7] = 90 + 273.15; + + Dlndjs_Dlnd_P[0] = 0.1F; + Dlndjs_Dlnd_P[1] = 1; + Dlndjs_Dlnd_P[2] = 2; + Dlndjs_Dlnd_P[3] = 3; + Dlndjs_Dlnd_P[4] = 4; + Dlndjs_Dlnd_P[5] = 5; + Dlndjs_Dlnd_P[6] = 6; + Dlndjs_Dlnd_P[7] = 7; + Dlndjs_Dlnd_P[8] = 8; + Dlndjs_Dlnd_P[9] = 9; + Dlndjs_Dlnd_P[10] = 10; + Dlndjs_Dlnd_Data[0][0] = 976; + Dlndjs_Dlnd_Data[1][0] = 1027; + Dlndjs_Dlnd_Data[2][0] = 1071; + Dlndjs_Dlnd_Data[3][0] = 1123; + Dlndjs_Dlnd_Data[4][0] = 1167; + Dlndjs_Dlnd_Data[5][0] = 1213; + Dlndjs_Dlnd_Data[6][0] = 1260; + Dlndjs_Dlnd_Data[7][0] = 1303; + Dlndjs_Dlnd_Data[0][1] = 991; + Dlndjs_Dlnd_Data[1][1] = 1040; + Dlndjs_Dlnd_Data[2][1] = 1082; + Dlndjs_Dlnd_Data[3][1] = 1135; + Dlndjs_Dlnd_Data[4][1] = 1178; + Dlndjs_Dlnd_Data[5][1] = 1224; + Dlndjs_Dlnd_Data[6][1] = 1270; + Dlndjs_Dlnd_Data[7][1] = 1312; + Dlndjs_Dlnd_Data[0][2] = 1014; + Dlndjs_Dlnd_Data[1][2] = 1063; + Dlndjs_Dlnd_Data[2][2] = 1106; + Dlndjs_Dlnd_Data[3][2] = 1153; + Dlndjs_Dlnd_Data[4][2] = 1196; + Dlndjs_Dlnd_Data[5][2] = 1239; + Dlndjs_Dlnd_Data[6][2] = 1281; + Dlndjs_Dlnd_Data[7][2] = 1323; + Dlndjs_Dlnd_Data[0][3] = 1044; + Dlndjs_Dlnd_Data[1][3] = 1091; + Dlndjs_Dlnd_Data[2][3] = 1127; + Dlndjs_Dlnd_Data[3][3] = 1174; + Dlndjs_Dlnd_Data[4][3] = 1216; + Dlndjs_Dlnd_Data[5][3] = 1257; + Dlndjs_Dlnd_Data[6][3] = 1297; + Dlndjs_Dlnd_Data[7][3] = 1338; + Dlndjs_Dlnd_Data[0][4] = 1073; + Dlndjs_Dlnd_Data[1][4] = 1118; + Dlndjs_Dlnd_Data[2][4] = 1149; + Dlndjs_Dlnd_Data[3][4] = 1195; + Dlndjs_Dlnd_Data[4][4] = 1236; + Dlndjs_Dlnd_Data[5][4] = 1275; + Dlndjs_Dlnd_Data[6][4] = 1313; + Dlndjs_Dlnd_Data[7][4] = 1352; + Dlndjs_Dlnd_Data[0][5] = 1114; + Dlndjs_Dlnd_Data[1][5] = 1151; + Dlndjs_Dlnd_Data[2][5] = 1180; + Dlndjs_Dlnd_Data[3][5] = 1224; + Dlndjs_Dlnd_Data[4][5] = 1261; + Dlndjs_Dlnd_Data[5][5] = 1297; + Dlndjs_Dlnd_Data[6][5] = 1333; + Dlndjs_Dlnd_Data[7][5] = 1372; + Dlndjs_Dlnd_Data[0][6] = 1156; + Dlndjs_Dlnd_Data[1][6] = 1185; + Dlndjs_Dlnd_Data[2][6] = 1211; + Dlndjs_Dlnd_Data[3][6] = 1253; + Dlndjs_Dlnd_Data[4][6] = 1287; + Dlndjs_Dlnd_Data[5][6] = 1320; + Dlndjs_Dlnd_Data[6][6] = 1352; + Dlndjs_Dlnd_Data[7][6] = 1391; + Dlndjs_Dlnd_Data[0][7] = 1207; + Dlndjs_Dlnd_Data[1][7] = 1230; + Dlndjs_Dlnd_Data[2][7] = 1250; + Dlndjs_Dlnd_Data[3][7] = 1289; + Dlndjs_Dlnd_Data[4][7] = 1318; + Dlndjs_Dlnd_Data[5][7] = 1346; + Dlndjs_Dlnd_Data[6][7] = 1374; + Dlndjs_Dlnd_Data[7][7] = 1412; + Dlndjs_Dlnd_Data[0][8] = 1261; + Dlndjs_Dlnd_Data[1][8] = 1276; + Dlndjs_Dlnd_Data[2][8] = 1289; + Dlndjs_Dlnd_Data[3][8] = 1324; + Dlndjs_Dlnd_Data[4][8] = 1350; + Dlndjs_Dlnd_Data[5][8] = 1373; + Dlndjs_Dlnd_Data[6][8] = 1396; + Dlndjs_Dlnd_Data[7][8] = 1432; + Dlndjs_Dlnd_Data[0][9] = 1331; + Dlndjs_Dlnd_Data[1][9] = 1331; + Dlndjs_Dlnd_Data[2][9] = 1335; + Dlndjs_Dlnd_Data[3][9] = 1366; + Dlndjs_Dlnd_Data[4][9] = 1385; + Dlndjs_Dlnd_Data[5][9] = 1403; + Dlndjs_Dlnd_Data[6][9] = 1424; + Dlndjs_Dlnd_Data[7][9] = 1456; + Dlndjs_Dlnd_Data[0][10] = 1405; + Dlndjs_Dlnd_Data[1][10] = 1389; + Dlndjs_Dlnd_Data[2][10] = 1383; + Dlndjs_Dlnd_Data[3][10] = 1409; + Dlndjs_Dlnd_Data[4][10] = 1421; + Dlndjs_Dlnd_Data[5][10] = 1435; + Dlndjs_Dlnd_Data[6][10] = 1451; + Dlndjs_Dlnd_Data[7][10] = 1482; + + if (tempT < Dlndjs_Dlnd_T[0]) + { + tempT = Dlndjs_Dlnd_T[0]; + } + if (tempT > Dlndjs_Dlnd_T[7]) + { + tempT = Dlndjs_Dlnd_T[7]; + } + if (tempP_jy < Dlndjs_Dlnd_P[0]) + { + tempP_jy = Dlndjs_Dlnd_P[0]; + } + if (tempP_jy > Dlndjs_Dlnd_P[10]) + { + tempP_jy = Dlndjs_Dlnd_P[10]; + } + + for ( i = 0; i <= 6; i++) + { + if (tempT >= Dlndjs_Dlnd_T[i] && tempT <= Dlndjs_Dlnd_T[i + 1]) + { + m = i; + break; + } + } + + for (i = 0; i <= 9; i++) + { + if (tempP_jy >= Dlndjs_Dlnd_P[i] && tempP_jy <= Dlndjs_Dlnd_P[i + 1]) + { + n = i; + break; + } + } + + if (Dlndjs_Dlnd_P[n + 1] - Dlndjs_Dlnd_P[n] != 0) + { + ky = (tempP_jy - Dlndjs_Dlnd_P[n]) / (Dlndjs_Dlnd_P[n + 1] - Dlndjs_Dlnd_P[n]); + } + else + { + ky = 0; + } + if (Dlndjs_Dlnd_T[m + 1] - Dlndjs_Dlnd_T[m] != 0) + { + kx = (tempT - Dlndjs_Dlnd_T[m]) / (Dlndjs_Dlnd_T[m + 1] - Dlndjs_Dlnd_T[m]); + } + else + { + kx = 0; + } + s1 = Dlndjs_Dlnd_Data[m][n] + (Dlndjs_Dlnd_Data[m][n + 1] - Dlndjs_Dlnd_Data[m][n]) * ky; + s2 = Dlndjs_Dlnd_Data[m + 1][n] + (Dlndjs_Dlnd_Data[m + 1][n + 1] - Dlndjs_Dlnd_Data[m + 1][n]) * ky; + return (s1 + (s2 - s1) * kx) / 100000.0D; + } + //可膨胀系数计算 + private double KePenZhang_JiSuan(double tempP_jy, double tempDp_Pa, double tempZjb, double tempDszs, int JIeliuType, int JiSuanBiaoZhun) // 求可膨胀系数 + { + double returnValue = 0; + //0标准孔板 + //1ISA1932喷嘴 + //2长径喷嘴 + //3文丘里喷嘴 + //4粗铸收缩段经典文丘里管 + //5机械加工收缩段经典文丘里管 + //6粗焊铁板收缩段经典文丘里管 + //7 1/4圆孔板 + + double tuo ; + switch (JIeliuType) + { + case 0: //孔板流量计算 + switch (JiSuanBiaoZhun) + { + case 0: //6143-2004 + tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); + returnValue = 1 - (0.351 + 0.256 * Math.pow(tempZjb, 4) + 0.93 * Math.pow(tempZjb, 8)) * (1 - Math.pow(tuo, (1 / tempDszs))); + break; + case 1: //6143-1996 + returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); + break; + } + break; + case 1: + case 2: + case 3: + case 4: + case 5: + case 6: + switch (JiSuanBiaoZhun) + { + case 0: + //标准喷嘴 iso5167-2002 + tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); + returnValue = Math.pow((((tempDszs * Math.pow(tuo, (2 / tempDszs))) / (tempDszs - 1)) * ((1 - Math.pow(tempZjb, 4)) / (1 - Math.pow(tempZjb, 4) * Math.pow(tuo, (2 / tempDszs)))) * ((1 - Math.pow(tuo, ((tempDszs - 1) / tempDszs))) / (1 - tuo))), 0.5); + break; + + case 1: //iso5167-93 + returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); + break; + } + break; + case 7: //1/4圆孔板 + returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); + break; + case 8: //锥形入口孔板 + tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); + returnValue = 1 - (0.351 + 0.256 * Math.pow(tempZjb, 4) + 0.93 * Math.pow(tempZjb, 8)) * (1 - Math.pow(tuo, (1 / tempDszs))); + + tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); + + returnValue = 0.5 * (returnValue + Math.pow((((tempDszs * Math.pow(tuo, (2 / tempDszs))) / (tempDszs - 1)) * ((1 - Math.pow(tempZjb, 4)) / (1 - Math.pow(tempZjb, 4) * Math.pow(tuo, (2 / tempDszs)))) * ((1 - Math.pow(tuo, ((tempDszs - 1) / tempDszs))) / (1 - tuo))), 0.5)); + break; + case 9: //偏心孔板 + returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); + break; + + } + return returnValue; + } + //流出系数计算 + public final double C_JiSuan(double tempGj, double tempZjb, double tempReD, double tempconQvA, int tempQyfs, int JieLiuType, int JiSuanBiaoZhun) + { + double returnValue = 0; + //流出系数计算函数 + //输入：直径比，雷诺数，取压方式，节流装置类型，计算采用标准 + //输出:流出系数 + + //jieliutype + + //0标准孔板 + //1ISA1932喷嘴 + //2长径喷嘴 + //3文丘里喷嘴 + //4粗铸收缩段经典文丘里管 + //5机械加工收缩段经典文丘里管 + //6粗焊铁板收缩段经典文丘里管 + + double L1 = 0; + double L2 = 0; + switch (JieLiuType) + { + case 0: //孔板 + switch (tempQyfs) + { + case 0: + L1 = 25.4 / tempGj; + L2 = L1; + break; + case 1: + L1 = 0; + L2 = 0; + break; + case 2: + L1 = 1; + L2 = 0.47F; + break; + } + switch (JiSuanBiaoZhun) + { + case 0: //6143-2004 + if (tempGj >= 71.12) + { + returnValue = 0.5961 + 0.0261 * Math.pow(tempZjb, 2) - 0.216 * Math.pow(tempZjb, 8) + 0.000521 * Math.pow((1000000.0D * tempZjb / tempReD), 0.7) + (0.0188 + 0.0063 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * Math.pow(tempZjb, 3.5) * Math.pow((1000000.0D / tempReD), 0.3) + (0.043 + 0.08 * Math.exp(-10 * L1) - 0.123 * Math.exp(-7 * L1)) * (1 - 0.11 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * (Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1))) - 0.031 * (2 * L2 / (1 - tempZjb) - 0.8 * Math.pow((2 * L2 / (1 - tempZjb)), 1.1)) * Math.pow(tempZjb, 1.3); + } + else if (tempGj < 71.12) + { + returnValue = (0.5961 + 0.0261 * Math.pow(tempZjb, 2) - 0.216 * Math.pow(tempZjb, 8) + 0.000521 * Math.pow((1000000.0D * tempZjb / tempReD), 0.7) + (0.0188 + 0.0063 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * Math.pow(tempZjb, 3.5) * Math.pow((1000000.0D / tempReD), 0.3) + (0.043 + 0.08 * Math.exp(-10 * L1) - 0.123 * Math.exp(-7 * L1)) * (1 - 0.11 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1)) - 0.031 * (2 * L2 / (1 - tempZjb) - 0.8 * Math.pow((2 * L2 * (1 - tempZjb)), 1.1)) * Math.pow(tempZjb, 1.3) + 0.011 * (0.75 - tempZjb) * (2.8 - tempGj / 25.4)); + } + break; + case 1: //6143-1996 + if (0.09 * L1 >= 0.039) + { + returnValue = 0.5959 + 0.0312 * Math.pow(tempZjb, (2.1)) - 0.184 * Math.pow(tempZjb, 8) + 0.0029 * Math.pow(tempZjb, 2.5) * Math.pow((1000000.0D / tempReD), 0.75) + 0.039 * Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1)) - 0.0337 * L1 * Math.pow(tempZjb, 3); + } + else if (0.09 * L1 < 0.039) + { + returnValue = (0.5959 + 0.0312 * Math.pow(tempZjb, (2.1)) - 0.184 * Math.pow(tempZjb, 8) + 0.0029 * Math.pow(tempZjb, 2.5) * Math.pow((1000000.0D / tempReD), 0.75) + 0.09 * L1 * Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1)) - 0.0337 * L1 * Math.pow(tempZjb, 3)); + } + break; + } + break; + case 1: //ISA1932喷嘴 + returnValue = (0.99 - 0.2262 * Math.pow(tempZjb, 4.1) - (0.00175 * Math.pow(tempZjb, 2) - 0.0033 * Math.pow(tempZjb, 4.15)) * (1000000.0D / Math.pow(tempReD, 1.15))); + break; + + case 2: //长径喷嘴 + returnValue = (0.9965 - 0.00653 * Math.pow(tempZjb, 0.5) * Math.pow((1000000.0D / tempReD), 0.5)); + break; + case 3: //文丘里喷嘴 + returnValue = (0.9858 - 0.196 * Math.pow(tempZjb, 4.5)); + break; + case 4: //粗铸收缩段经典文丘里管 + returnValue = (0.984F); + break; + case 5: //机械加工收缩段经典文丘里管 + returnValue = (0.995F); + break; + case 6: //粗焊铁板收缩段经典文丘里管 + returnValue = (0.985F); + break; + case 7: //1/4圆孔板 + returnValue = 0.73823 - 0.3309 * tempZjb - 1.1615 * Math.pow(tempZjb, 2) + 1.5084 * Math.pow(tempZjb, 3); + break; + case 8: //锥形入口孔板 + returnValue = (0.734F); + break; + case 9: //偏心孔板 + returnValue = 0.9355 - 1.6889 * tempZjb + 3.0428 * Math.pow(tempZjb, 2) - 1.7989 * Math.pow(tempZjb, 3); + break; + } + double tempRed1 ; + switch (JieLiuType) + { + case 0: //孔板流量计算 + tempRed1 = tempconQvA * returnValue; + switch (JieLiuType) + { + case 0: + if (tempRed1 < (170 * Math.pow(tempZjb, 2) * tempGj)) + { + JOptionPane.showMessageDialog(null, "雷诺数超过标准孔板的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + return returnValue; + } + break; + case 1: + if (tempZjb >= 0.1 & tempZjb <= 0.56) + { + if (tempRed1 < 5000) + { + JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + return returnValue; + } + } + if (tempZjb > 0.56) + { + if (tempRed1 < (16000 * Math.pow(tempZjb, 2) * tempGj)) + { + JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + return returnValue; + } + } + break; + } + break; + + case 1: //标准喷嘴 + tempRed1 = tempconQvA * returnValue; + if (tempZjb >= 0.3 & tempZjb < 0.44) + { + if (tempRed1 < 70000 | tempRed1 > 10000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + } + } + if (tempZjb >= 0.44 & tempZjb < 0.8) + { + if (tempRed1 < 20000 | tempRed1 > 10000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + } + } + break; + case 2: //长径喷嘴 + tempRed1 = tempconQvA * returnValue; + if (tempRed1 < 10000.0D | tempRed1 > 10000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过长径喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + } + break; + case 3: //文丘里喷嘴 + tempRed1 = tempconQvA * returnValue; + if (tempRed1 < 150000.0D | tempRed1 > 2000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过文丘里喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + } + break; + case 4: //粗铸收缩段经典文丘里管 + tempRed1 = tempconQvA * returnValue; + if (tempRed1 < 200000.0D | tempRed1 > 2000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过粗铸收缩段经典文丘里管的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + return returnValue; + } + break; + case 5: //机械加工收缩段经典文丘里管 + tempRed1 = tempconQvA * returnValue; + if (tempRed1 < 200000.0D | tempRed1 > 1000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过机械加工收缩段经典文丘里管的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + } + break; + case 6: //粗焊铁板收缩段经典文丘里管 + tempRed1 = tempconQvA * returnValue; + if (tempRed1 < 200000.0D | tempRed1 > 2000000.0D) + { + JOptionPane.showMessageDialog(null, "雷诺数超过粗焊铁板收缩段经典文丘里管的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); + } + break; + + //标准孔板 + //ISA1932喷嘴 + //长径喷嘴 + //文丘里喷嘴 + //粗铸收缩段经典文丘里管 + //机械加工收缩段经典文丘里管 + //粗焊铁板收缩段经典文丘里管 + + + + + } + return returnValue; + } + + private double BkTable(double temPrk, double tempKj, int tempBkjsff) { + double[] BkTable_x = new double[10]; + double[] BkTable_Y = new double[10]; + //On Error Resume Next VBConversions Warning: On Error Resume Next not supported in C# + if (tempBkjsff == 1) { + double tempRkBiKj = temPrk / tempKj; + // static double[] x = new double[10]; //VBConversions Note: Static variable moved to class level and renamed BkTable_x. Local static variables are not supported in C#. + // static double[] Y = new double[10]; //VBConversions Note: Static variable moved to class level and renamed BkTable_Y. Local static variables are not supported in C#. + int i; + int xIndex = 0; + //If x(0) = 0 Then + BkTable_x[0] = 0.0004F; + BkTable_x[1] = 0.001F; + BkTable_x[2] = 0.002F; + BkTable_x[3] = 0.004F; + BkTable_x[4] = 0.006F; + BkTable_x[5] = 0.008F; + BkTable_x[6] = 0.01F; + BkTable_x[7] = 0.012F; + BkTable_x[8] = 0.014F; + BkTable_x[9] = 0.015F; + + BkTable_Y[0] = 1; + BkTable_Y[1] = 1.005F; + BkTable_Y[2] = 1.012F; + BkTable_Y[3] = 1.022F; + BkTable_Y[4] = 1.032F; + BkTable_Y[5] = 1.04F; + BkTable_Y[6] = 1.048F; + BkTable_Y[7] = 1.055F; + BkTable_Y[8] = 1.062F; + BkTable_Y[9] = 1.065F; + //End If + + if (tempRkBiKj <= 0.0004) { + return 1; + } + if (tempRkBiKj > 0.015) { + return 1.065F; + } + for (i = 0; i <= 8; i++) { + if (tempRkBiKj >= BkTable_x[i] && tempRkBiKj <= BkTable_x[i + 1]) { + xIndex = i; + break; + } + } + return BkTable_Y[xIndex] + (tempRkBiKj - BkTable_x[xIndex]) * (BkTable_Y[xIndex + 1] - BkTable_Y[xIndex]) / (BkTable_x[xIndex + 1] - BkTable_x[xIndex]); + } + return 0; + } + //管道粗糙度计算 + private double CcdXsjs(double tempPipeType, double tempGj, double tempZjb, double TempRed) + { + double returnValue ; + //粗糙度系数计算 + double Jdccd = 0; //绝对粗糙度 + double Xdccd ; //相对粗糙度 + //Dim CcdXs As single + double s1 ; + double s2 ; + double ky ; + double kx ; + int i ; + int m = 0; + int n = 0; + if (tempPipeType == 0) + { + Jdccd = 0.029F;} + else if (tempPipeType == 1) + { + Jdccd = 0.075F; + } + else if (tempPipeType == 2) + { + Jdccd = 0.075F; + } + else if (tempPipeType == 3) + { + Jdccd = 0.075F; + } + else if (tempPipeType == 4) + { + Jdccd = 0.1F; + } + else if (tempPipeType == 5) + { + Jdccd = 0.15F; + } + else if (tempPipeType == 6) + { + Jdccd = 1; + } + else if (tempPipeType == 7) + { + Jdccd = 2.1F; + } + else if (tempPipeType == 8) + { + Jdccd = 0.04F; + } + else if (tempPipeType == 9) + { + Jdccd = 0.15F; + } + else if (tempPipeType == 10) + { + Jdccd = 0.13F; + } + else if (tempPipeType == 11) + { + Jdccd = 0.25F; + } + Xdccd = tempGj / Jdccd; + if (Xdccd < 400) + { + JOptionPane.showMessageDialog(null, "粗糙度取得太高,粗略计算", "提示", JOptionPane.PLAIN_MESSAGE); + Xdccd = 400; + } + if (Xdccd >= 3400) + { + Xdccd = 3400; + } + + if (Xdccd < 3200 & Math.pow(tempZjb, 2) > 0.1 & Math.pow(tempZjb, 2) < 0.64) + { + + int[] Xdccdb = new int[10]; + double[] Btf = new double[8]; + double[][] CcdXsb = new double[10][8]; + + Xdccdb[0] = 400; + Xdccdb[1] = 800; + Xdccdb[2] = 1200; + Xdccdb[3] = 1600; + Xdccdb[4] = 2000; + Xdccdb[5] = 2400; + Xdccdb[6] = 2800; + Xdccdb[7] = 3200; + Xdccdb[8] = 3400; + + Btf[0] = 0.1F; + Btf[1] = 0.2F; + Btf[2] = 0.3F; + Btf[3] = 0.4F; + Btf[4] = 0.5F; + Btf[5] = 0.6F; + Btf[6] = 0.64F; + + CcdXsb[0][0] = 1.002; + CcdXsb[1][0] = 1; + CcdXsb[2][0] = 1; + CcdXsb[3][0] = 1; + CcdXsb[4][0] = 1; + CcdXsb[5][0] = 1; + CcdXsb[6][0] = 1; + CcdXsb[7][0] = 1; + CcdXsb[8][0] = 1; + CcdXsb[0][1] = 1.003; + CcdXsb[1][1] = 1.002; + CcdXsb[2][1] = 1.001; + CcdXsb[3][1] = 1; + CcdXsb[4][1] = 1; + CcdXsb[5][1] = 1; + CcdXsb[6][1] = 1; + CcdXsb[7][1] = 1; + + CcdXsb[0][2] = 1.006; + CcdXsb[1][2] = 1.004; + CcdXsb[2][2] = 1.002; + CcdXsb[3][2] = 1.001; + CcdXsb[4][2] = 1; + CcdXsb[5][2] = 1; + CcdXsb[6][2] = 1; + CcdXsb[7][2] = 1; + + CcdXsb[0][3] = 1.009; + CcdXsb[1][3] = 1.006; + CcdXsb[2][3] = 1.004; + CcdXsb[3][3] = 1.002; + CcdXsb[4][3] = 1.001; + CcdXsb[5][3] = 1; + CcdXsb[6][3] = 1; + CcdXsb[7][3] = 1; + + CcdXsb[0][4] = 1.014; + CcdXsb[1][4] = 1.009; + CcdXsb[2][4] = 1.006; + CcdXsb[3][4] = 1.004; + CcdXsb[4][4] = 1.002; + CcdXsb[5][4] = 1.001; + CcdXsb[6][4] = 1; + CcdXsb[7][4] = 1; + + CcdXsb[0][5] = 1.02; + CcdXsb[1][5] = 1.013; + CcdXsb[2][5] = 1.009; + CcdXsb[3][5] = 1.006; + CcdXsb[4][5] = 1.003; + CcdXsb[5][5] = 1.002; + CcdXsb[6][5] = 1; + CcdXsb[7][5] = 1; + + CcdXsb[0][6] = 1.024; + CcdXsb[1][6] = 1.016; + CcdXsb[2][6] = 1.011; + CcdXsb[3][6] = 1.007; + CcdXsb[4][6] = 1.004; + CcdXsb[5][6] = 1.002; + CcdXsb[6][6] = 1.002; + CcdXsb[7][6] = 1; + + + for (i = 0; i <= 8; i++) + { + if (Xdccd >= Xdccdb[i] && Xdccd <= Xdccdb[i + 1]) + { + m = i; + break; + } + } + for (i = 0; i <= 6; i++) + { + if (Math.pow(tempZjb, 2) >= Btf[i] && Math.pow(tempZjb, 2) <= Btf[i + 1]) + { + n = i; + break; + } + } + + ky = (Btf[n + 1] - Btf[n] != 0) ? ((Math.pow(tempZjb, 2) - Btf[n]) / (Btf[n + 1] - Btf[n])) : 0; + + kx = (Xdccdb[m + 1] - Xdccdb[m] != 0) ? ((Xdccd - Xdccdb[m]) / (Xdccdb[m + 1] - Xdccdb[m])) : 0; + + s1 = CcdXsb[m][n] + (CcdXsb[m][n + 1] - CcdXsb[m][n]) * ky; + s2 = CcdXsb[m + 1][n] + (CcdXsb[m + 1][n + 1] - CcdXsb[m + 1][n]) * ky; + returnValue = s1 + (s2 - s1) * kx; + + returnValue = TempRed > 1000000.0 ? returnValue : ((returnValue - 1) * Math.pow((Math.log10(TempRed) / 2), 2) + 1); + + } + else + { + returnValue = 1; + } + return returnValue; + } +} diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/FlowController.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/FlowController.java index 2b97c1d..d84e54b 100644 --- a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/FlowController.java +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/FlowController.java @@ -14,6 +14,9 @@ import org.springframework.web.bind.annotation.RestController; import javax.swing.*; import java.util.Optional; +import static com.ruoyi.ngCalTools.controller.DpFlowCalc.iterativeFlowCalculation; +import static com.ruoyi.ngCalTools.controller.DpFlowCalc.thermalExpansionCorrection; + @RestController @RequestMapping("/flowCalcTools") public class FlowController { @@ -56,6 +59,30 @@ public class FlowController { flowProps.setdTf(tempTf); gasProps.dCbtj = flowProps.getdCbtj(); + switch (gasProps.dCbtj) + { + case 2: + gasProps.setdPb(101325); + gasProps.setdTb( 273.15); + flowProps.setdPb_M(101325); + flowProps.setdTb_M(273.15); + break; + case 1: + + gasProps.setdPb(101325); + gasProps.setdTb( 288.15); + flowProps.setdPb_M(101325); + flowProps.setdTb_M(288.15); + break; + case 0: + + gasProps.setdPb(101325); + gasProps.setdTb( 293.15); + flowProps.setdPb_M(101325); + flowProps.setdTb_M(293.15); + break; + } + String[] stringArray = flowProps.getdngComponents().split("_"); double[] doubleArray = new double[stringArray.length]; // 遍历字符串数组，将每个元素转换为 double 类型 for (int i = 0; i < stringArray.length; i++) { @@ -85,7 +112,9 @@ public class FlowController { SdFlowCal(gasProps, flowProps); break; case 3:// 临界流函数流量计 - NozellFLowCal(gasProps, flowProps); +// NozellFLowCal(gasProps, flowProps); +// NozellFLowCal_New(gasProps, flowProps); + NozellFlowCalc. executeFullCalculation(flowProps,gasProps); break; } Object[] resultArray = {flowProps, gasProps}; @@ -106,441 +135,65 @@ public class FlowController { flowProps.setdEFlowb(flowProps.getdVFlowb() * gasProps.dHhvm); } - /** - 标准孔板流量计算 - */ - public void OFlowCal(GasProps gasProps,FlowProps flowProps){ - flowProps.setdOrificeD(flowProps.getdOrificeD() * (1 + 0.000001 * (flowProps.getdOrificeMaterial()) * (flowProps.getdTf() - 293.15))); - flowProps.setdPipeD(flowProps.getdPipeD() * (1 + 0.000001 * (flowProps.getdPipeMaterial()) * (flowProps.getdTf() - 293.15))); - flowProps.setdBeta(flowProps.getdOrificeD() / flowProps.getdPipeD()); -// flowProps.setdBeta(0.5972); + public final void NozellFLowCal_New(GasProps gasProps, FlowProps flowProps) { + double nozzleThroatArea; // 喷嘴喉部面积 + double stagnationPressure; // 滞止压力 P0 + double stagnationTemp; // 滞止温度 T0 + double correctedPipeDiam; // 修正后管道直径 + double correctedThroatDiam;// 修正后喉部直径 + double beta; // 直径比 β = d/D + double massFlowRate; // 质量流量 Qm - //求渐近速度系数 E - flowProps.setdE(1 / Math.pow((1 - Math.pow(flowProps.getdBeta(), 4)), 0.5)); -// flowProps.setdE(1.0615); + try { + // 1. 热膨胀修正 + correctedThroatDiam = flowProps.getdOrificeD() + * (1 + 1e-6 * flowProps.getdOrificeMaterial() * (flowProps.getdTf() - 293.15)); + correctedPipeDiam = flowProps.getdPipeD() + * (1 + 1e-6 * flowProps.getdPipeMaterial() * (flowProps.getdTf() - 293.15)); + beta = correctedThroatDiam / correctedPipeDiam; - //求相对密度系数 FG - // LiuTiType = 1 - flowProps.setdFG(Math.pow((1 / gasProps.dRD_Real), 0.5)); -// flowProps.setdFG(1.2531); - //求流动温度系数 'FT - flowProps.setdFT(Math.pow((293.15 / flowProps.getdTf()), 0.5)); -// flowProps.setdFT(1.0086); - //求等熵指数????????????????????????????????? + // 2. 计算喉部面积 + nozzleThroatArea = Math.PI * Math.pow(correctedThroatDiam / 2, 2); - flowProps.setdKappa(gasProps.dKappa); -// flowProps.setdKappa(1.357); + // 3. 计算马赫数 (假设流速来自flowProps) + double velocity = flowProps.getdVelocityFlow(); + double speedOfSound = Math.sqrt(gasProps.dKappa * gasProps.dZf + * 8314.51 / gasProps.dMrx * flowProps.getdTf()); + double machNumber = velocity / speedOfSound; - //求动力粘度 dlnd - flowProps.setdDViscosity(Dlndjs(flowProps.getdPf()/1e6, flowProps.getdTf())); -// flowProps.setdDViscosity(0.01096); + // 4. 滞止压力 & 温度 (ISO 9300) + stagnationPressure = flowProps.getdPf() + * Math.pow(1 + (gasProps.dKappa - 1)/2 * machNumber*machNumber, + gasProps.dKappa/(gasProps.dKappa - 1)); + stagnationTemp = flowProps.getdTf() + * (1 + (gasProps.dKappa - 1)/2 * machNumber*machNumber); - //求可膨胀系数 - flowProps.setdDExpCoefficient(KePenZhang_JiSuan(flowProps.getdPf() , flowProps.getdDp(), flowProps.getdBeta(), gasProps.getdKappa(), flowProps.getdCoreType(), 0)); -// flowProps.setdDExpCoefficient(0.9977); + // 5. 计算质量流量 + double gasConstant = 8314.51 / gasProps.dMrx; + massFlowRate = nozzleThroatArea * flowProps.getdCd() + * Math.sqrt(gasProps.dKappa / (gasConstant * stagnationTemp)) + * stagnationPressure + * Math.pow(2/(gasProps.dKappa + 1), (gasProps.dKappa + 1)/(2*(gasProps.dKappa - 1))); -// gasProps.dFpv=1.0195; -// gasProps.dHhvv=39.944; -// gasProps.dRD_Real=0.6368; - - //迭代计算流量和流出系数 - - double conQvA = 0.0000031795 * (1530000.0D * gasProps.dRD_Real / (flowProps.getdDViscosity() * flowProps.getdOrificeD())) * flowProps.getdE() * Math.pow(flowProps.getdOrificeD(), 2) * gasProps.dFpv * flowProps.getdFG() * flowProps.getdDExpCoefficient() * flowProps.getdFT() * Math.sqrt(flowProps.getdPf() * flowProps.getdDp() / 1e6); - double[] XQv = new double[4]; - double[] CQv = new double[4]; - double[] dQv = new double[4]; - double Qn = 0; - int n ; - XQv[0] = 1000000.0F; - boolean xhFlag=true; - _100000: - - while(xhFlag) { - for (n = 1; n <= 2; n++) { - CQv[n] = C_JiSuan(flowProps.getdOrificeD(), flowProps.getdBeta(), XQv[n - 1], conQvA, flowProps.getdPtmode(), flowProps.getdCoreType(), 0); - XQv[n] = conQvA * CQv[n]; - dQv[n] = XQv[n] - XQv[n - 1]; + // 6. 大直径比修正 + if (beta > 0.25) { + massFlowRate *= BetaG25(stagnationPressure, stagnationTemp, beta, gasProps); } - if (XQv[2] == XQv[1] || dQv[2] == dQv[1]) { - Qn = conQvA * flowProps.getdDViscosity() * flowProps.getdOrificeD() * CQv[2] / (1530000.0D * gasProps.dRD_Real); - //管道雷诺数 - flowProps.setdRnPipe(XQv[2]); - //流出系数 - flowProps.setdCd(CQv[2]); - } - if (Math.abs((conQvA - XQv[2] / CQv[1]) / conQvA) > 0.00000000000000005) { - XQv[0] = XQv[n - 1] - dQv[n - 1] * ((XQv[n - 1] - XQv[n - 2]) / (dQv[n - 1] - dQv[n - 2])); -//C# TO JAVA CONVERTER TODO TASK: There is no 'goto' in Java: - continue _100000; - } - else - { - xhFlag=false; - } - } + // 7. 设置结果 + flowProps.setdMFlowb(massFlowRate); + flowProps.setdVFlowb(massFlowRate / gasProps.dRhob); + flowProps.setdEFlowb(flowProps.getdVFlowb() * gasProps.dHhvm); + flowProps.setdVFlowf(FlowConvert_BaseToWork(flowProps, gasProps)); - //孔板锐利度系数Bk - flowProps.setdOrificeSharpness(1); - if (flowProps.getdCoreType() == 0) + } catch (Exception ex) { - flowProps.setdBk((flowProps.getdOrificeSharpness() == 0) ? (BkTable(flowProps.getdOrificeRk(), flowProps.getdOrificeD(), 1)) : (flowProps.getdOrificeSharpness())); } - else - { - flowProps.setdBk(1); - } - - - //管道粗糙度系数 Gme - flowProps.setdRoughNessPipe(CcdXsjs(flowProps.getdPipeType(), flowProps.getdPipeD(), flowProps.getdBeta(), flowProps.getdRnPipe())); - //修正后的流出系数 - flowProps.setdCd(flowProps.getdCd() * flowProps.getdBk() * flowProps.getdRoughNessPipe()); -// flowProps.setdCd(0.6039); - //标况体积流量 m³、s - flowProps.setdVFlowb(Qn * flowProps.getdBk() * flowProps.getdRoughNessPipe()); - //工况体积流量 - flowProps.setdVFlowf(FlowConvert_BaseToWork(flowProps, gasProps)); - //标况质量流量 - flowProps.setdMFlowb(flowProps.getdVFlowb() * gasProps.dRhob); - //标况能量流量 - flowProps.setdEFlowb(flowProps.getdVFlowb() * gasProps.dHhvv); - //管道内天然气流速 - flowProps.setdVelocityFlow(flowProps.getdVFlowf() / (3.1415926 * Math.pow((flowProps.getdPipeD() / 2000), 2))); - //压力损失 - flowProps.setdPressLost(YaLiSunShi(flowProps.getdCd(), flowProps.getdBeta(), flowProps.getdDp(), flowProps.getdCoreType())); - } - - /** - * 计算管道绝对粗糙度 K (GB/T 21446-2008 附录C) - * @param flowProps - * @return 粗糙度修正系数 K (保留4位小数) - */ - public static double calculateK(FlowProps flowProps){ - double Jdccd; - switch (flowProps.getdPipeType()) { - - case 0: - Jdccd = 0.029F; - break; - case 1: - case 2: - case 3: - Jdccd = 0.075F; - break; - case 4: - Jdccd = 0.1F; - break; - case 5: - Jdccd = 0.15F; - break; - case 6: - Jdccd = 1; - break; - case 7: - Jdccd = 2.1F; - break; - case 8: - Jdccd = 0.04F; - break; - case 9: - Jdccd = 0.15F; - break; - case 10: - Jdccd = 0.13F; - break; - case 11: - Jdccd = 0.25F; - break; - default: - // 处理未知类型（可选） - throw new IllegalArgumentException("未知的管道类型: "); - } - return Jdccd; - } - - /** - * 计算管道粗糙度修正系数 G_me (GB/T 21446-2008 附录C) - * @param D_pipe 管道内径 (mm) - * @param K 绝对粗糙度 (mm) - * @param C 未修正的流出系数 - * @return 粗糙度修正系数 G_me (保留4位小数) - */ - public static double calculateRoughnessFactor(double D_pipe, double K, double C) { - // 计算相对粗糙度 K/D - double K_over_D = K / D_pipe; - - // 判断是否需要修正 - if (K_over_D <= 0.0004) { - return 1.0000; - } - - // 计算修正项 - double term = (K_over_D * 1e6) - 400; // 转换为无量纲项 - if (term < 0) { - throw new IllegalArgumentException("K/D 超出修正公式适用范围"); - } - - double G_me = 1 + (0.011 / C) * Math.sqrt(term); - return FormatUtil.format(G_me, Optional.of(4)); // 保留四位小数 - } - - /** - * 热膨胀修正（GB/T 21446-2008 第6.2.1节） - * @param flowProps 流量参数 - */ - public static void thermalExpansionCorrection(FlowProps flowProps) { - // 孔板孔径修正 - double dOrificeMaterial = flowProps.getdOrificeMaterial(); // 材料膨胀系数 (e.g. 16e-6/℃) - double dOrificeD = flowProps.getdOrificeD(); // 初始孔径(mm) - double dOrificeD_corrected = dOrificeD * (1 + 0.000001*dOrificeMaterial * (flowProps.getdTf() - 293.15)); // 20℃为参考温度 - flowProps.setdOrificeD(FormatUtil.format(dOrificeD_corrected, Optional.of(2))); - - // 管道内径修正 - double dPipeMaterial = flowProps.getdPipeMaterial(); - double dPipeD = flowProps.getdPipeD(); - double dPipeD_corrected = dPipeD * (1 + 0.000001*dPipeMaterial * (flowProps.getdTf() - 293.15)); - flowProps.setdPipeD(FormatUtil.format(dPipeD_corrected,Optional.of(2))); - - // 更新直径比β - flowProps.setdBeta( FormatUtil.format(flowProps.getdOrificeD() / flowProps.getdPipeD(),Optional.of(4))); - } - - /** - * 计算渐近速度系数E（GB/T 21446-2008 式(8)） - */ - public static void calculateE(FlowProps flowProps) { - double beta = flowProps.getdBeta(); - flowProps.setdE( FormatUtil.format(1 / Math.sqrt(1 - Math.pow(beta, 4)),Optional.of(4))); - } - - /** - * 计算相对密度系数FG（GB/T 21446-2008 式(9)） - */ - public static void calculateFG(FlowProps flowProps,GasProps gasProps) { - gasProps.setdRD_Real(gasProps.getdRD_Real()); - flowProps.setdFG( FormatUtil.format(1 / Math.sqrt(gasProps.dRD_Real),Optional.of(4))); - } - - /** - * 计算流动温度系数FT（GB/T 21446-2008 式(10)） - */ - public static void calculateFT(FlowProps flowProps) { - double Tb = flowProps.getdTb_M() ; // 参比温度转K - double Tf = flowProps.getdTf() ; // 工况温度转K - flowProps.setdFT( FormatUtil.format(Math.sqrt(Tb / Tf),Optional.of(4))); - } - - /** - * 计算可膨胀系数ε（GB/T 21446-2008 式(11)） - */ - public static void calculateEpsilon(FlowProps flowProps, GasProps gasProps) { - double P1 = flowProps.getdPf(); // 上游绝对压力(Pa) - double deltaP = flowProps.getdDp(); // 差压(Pa) - double beta = flowProps.getdBeta(); - double kappa = flowProps.getdKappa(); - - double tau = FormatUtil.format((P1 - deltaP) / P1, Optional.of(4)); // 压力比 - double epsilon=1 - (0.351 + 0.256 * Math.pow(beta, 4) + 0.93 * Math.pow(beta, 8)) * (1 - Math.pow(tau, 1/kappa)); - - flowProps.setdDExpCoefficient( FormatUtil.format(epsilon,Optional.of(4))); - } - /** - * 计算等熵指数κ（GB/T 21446-2008 推荐方法） - * @param gasProps 气体参数（需包含等熵指数或热容比） - * @return 等熵指数κ - */ - public static double calculateKappa(GasProps gasProps) { - // 若已通过AGA 8或其他方法计算，直接返回 -// if (gasProps.getdKappa() > 0) { -// return gasProps.getdKappa(); -// } - - // 近似公式：基于理想气体比热比和压缩因子修正 - double gamma = 1.3; // 天然气典型比热比（Cp/Cv≈1.3） - double Z = gasProps.getdZf(); // 工况压缩因子 - - // 修正公式（经验关系） - double kappa = gamma / (1 - (gamma - 1) * (1 / Z - 1)); - return FormatUtil.format(kappa,Optional.of(4)); } - - /** - * 迭代计算流量（GB/T 21446-2008 第7章） - * @param flowProps 流量参数 - * @param gasProps 气体参数 - */ - public static void iterativeFlowCalculation(FlowProps flowProps, GasProps gasProps) { - // 0. 单位转换 - double D = flowProps.getdPipeD() / 1000.0; // 管道内径(m) - double d = flowProps.getdOrificeD() / 1000.0; // 孔板孔径(m) - double beta = flowProps.getdBeta(); - double P1 = flowProps.getdPf(); // 绝对压力(Pa) - double deltaP = flowProps.getdDp(); // 差压(Pa) - double Tf = flowProps.getdTf() ; // 工况温度(K) - -flowProps.setdPb_M(gasProps.getdPb()); -flowProps.setdTb_M(gasProps.getdTb()); -flowProps.setdFpv(FormatUtil.format(gasProps.dFpv,Optional.of(4))); - - // 1. 计算中间参数 - calculateE(flowProps); - calculateFG(flowProps,gasProps); - calculateFT(flowProps); - - flowProps.setdKappa( calculateKappa(gasProps)); -// flowProps.setdDViscosity( calculateDynamicViscosity(gasProps)); - flowProps.setdDViscosity(FormatUtil.format( Dlndjs(flowProps.getdPf()/1e6,flowProps.getdTf()),Optional.of(5))); - calculateEpsilon(flowProps, gasProps); - - flowProps.setdKappa(1.357); -// flowProps.setdBeta(0.5792); -// flowProps.setdFpv(1.0195); -// flowProps.setdFG(1.2531); -// flowProps.setdE(1.0615); -// flowProps.setdDExpCoefficient(0.9977); -// flowProps.setdFT(1.0086); - - - // 2. 初始雷诺数估算（假设初始C=0.6） - double C_initial = 0.6; - double Qf_initial = (C_initial * flowProps.getdE() * flowProps.getdDExpCoefficient() * Math.PI * Math.pow(d, 2) / 4) - * Math.sqrt(2 * deltaP / (gasProps.dRhof * (1 - Math.pow(beta, 4)))); - flowProps.setdVFlowf(Qf_initial); // 初始工况流量(m³/s) - - // 3. 迭代参数 - double tolerance = 1e-6; - int maxIter = 100; - double currentC = C_initial; - double currentReD = calculateReD(Qf_initial, D, gasProps.dRhof, flowProps.getdDViscosity()); - int iter = 0; - double prevC=0; - // 4. 迭代循环 - do { - prevC = currentC; - - // 4.1 计算流出系数C（GB/T 21446-2008 附录A） - currentC =FormatUtil.format( calculateCd(beta, currentReD, flowProps.getdPipeD(), flowProps.getdPtmode()),Optional.of(4)); - - // 4.2 更新流量 - double Qf =FormatUtil.format ((currentC * flowProps.getdDExpCoefficient() * Math.PI * Math.pow(d, 2) / 4) - * Math.sqrt(2 * deltaP / (gasProps.dRhof * (1 - Math.pow(beta, 4)))),Optional.of(4)); - flowProps.setdVFlowf(Qf); - - // 4.3 更新雷诺数 - currentReD =FormatUtil.format( calculateReD(Qf, D, gasProps.dRhof, flowProps.getdDViscosity()), Optional.of(2)); - - iter++; - if (iter > maxIter) { - throw new RuntimeException("迭代未收敛，超过最大迭代次数！"); - } - } while (Math.abs(currentC - prevC) / currentC > tolerance); - // 5. 保存最终结果 - - // 在迭代计算流出系数后，添加粗糙度修正 - - double K = calculateK(flowProps); // 根据实际管道类型选择 - double G_me = calculateRoughnessFactor(flowProps.getdPipeD(), K, currentC); - double C_corrected = FormatUtil.format(currentC * G_me, Optional.of(4)); - - flowProps.setdCd(C_corrected); - flowProps.setdRoughNessPipe(G_me); - flowProps.setdRnPipe(currentReD); - - - - - // 6. 计算标况流量（GB/T 21446-2008 式(1)） - double Qn = flowProps.getdVFlowf() * (flowProps.getdFpv()*flowProps.getdFpv() * P1 / flowProps.getdPb_M()) - * (flowProps.getdTb_M()) / Tf; - flowProps.setdVFlowb(FormatUtil.format(Qn,Optional.of(4))); - - //标况质量流量 - flowProps.setdMFlowb(FormatUtil.format(flowProps.getdVFlowb() * gasProps.dRhob,Optional.of(4))); - //标况能量流量 - flowProps.setdEFlowb(FormatUtil.format(flowProps.getdVFlowb() * gasProps.dHhvv,Optional.of(4))); - //管道内天然气流速 - flowProps.setdVelocityFlow(FormatUtil.format(flowProps.getdVFlowf() / (Math.PI * Math.pow((flowProps.getdPipeD() / 2000), 2)),Optional.of(3))); - //压力损失 - flowProps.setdPressLost(FormatUtil.format(YaLiSunShi(flowProps.getdCd(), flowProps.getdBeta(), flowProps.getdDp(), flowProps.getdCoreType()),Optional.of(2))); - - } - - - - - - /** - * 计算天然气动力粘度（Sutherland公式简化版，适用于低压天然气） - * @param gasProps 气体参数（需包含温度、密度、组分） - * @return 动力粘度 (Pa·s) - */ - public static double calculateDynamicViscosity(GasProps gasProps) { - double T = gasProps.getdTf() ; // 工况温度(K) - double rho = gasProps.getdRhof(); // 工况密度(kg/m³) - double Rd = gasProps.getdRD_Real(); // 真实相对密度 - - // Sutherland常数（天然气典型值，C=120） - double C = 120.0; // Sutherland常数 - double mu0 = 1.0e-5; // 参考粘度 (Pa·s，@293.15K) - - // 计算粘度（简化公式，忽略压力修正） - double mu = mu0 * Math.pow(T / 293.15, 0.7) * (293.15 + C) / (T + C); - - // 密度修正（经验公式） - mu *= 1 + 0.001 * (rho - 1.2); // 密度对粘度的粗略修正 - - return FormatUtil.format( mu*1000,Optional.of(5)); - } - /** - * 计算雷诺数ReD（GB/T 21446-2008 式(5)） - */ - public static double calculateReD(double Qf, double D, double rho, double mu) { - return FormatUtil.format( (4 * Qf * rho*1000) / (Math.PI * D * mu),Optional.of(2)); // Qf: m³/s, D: m, rho: kg/m³, mu: Pa·s - } - - /** - * 计算流出系数C（GB/T 21446-2008 附录A） - */ - public static double calculateCd(double beta, double ReD, double D_mm, int ptMode) { - double L1, L2; - // 根据取压方式确定L1/L2（角接取压） - switch (ptMode) { - case 1: // 角接取压 - L1 = L2 = 0; // D单位为mm - break; - case 0: // 法兰取压 - L1 = L2 = 25.4 / D_mm; - break; - case 2: // D-D/2取压 - L1 = 1.0; - L2 = 0.47; - break; - default: - throw new IllegalArgumentException("不支持的取压方式"); - } - - double term1 = 0.5961 + 0.0261 * Math.pow(beta, 2) - 0.216 * Math.pow(beta, 8); - double term2 = 0.000521 * Math.pow(1e6 * beta / ReD, 0.7); - double A = Math.pow(19000 * beta / ReD, 0.8); - double term3 = (0.0188 + 0.0063 * A) * Math.pow(beta, 3.5) * Math.pow(1e6 / ReD, 0.3); - double term4 = (0.043 + 0.08 * Math.exp(-10 * L1) - 0.123 * Math.exp(-7 * L1)) - * (1 - 0.11 * A) * Math.pow(beta, 4) / (1 - Math.pow(beta, 4)); - double term5 = -0.031 * (2 * L2 / (1 - beta) - 0.8 * Math.pow(2 * L2 / (1 - beta), 1.1)) - * Math.pow(beta, 1.3); - - double Cd = term1 + term2 + term3 + term4 + term5; - - // 孔径<71.12mm修正 - if (D_mm < 71.12) { - Cd += 0.011 * (0.75 - beta) * (2.8 - D_mm / 25.4); - } - return FormatUtil.format(Cd,Optional.of(4)); - } - + /** 天然气音速喷嘴流量计算 @@ -688,689 +341,7 @@ flowProps.setdFpv(FormatUtil.format(gasProps.dFpv,Optional.of(4))); } return ylss; } - //查表计算粘度μ - public static double Dlndjs(double tempP_jy, double tempT) - { double[][] Dlndjs_Dlnd_Data = new double[8][11]; - double[] Dlndjs_Dlnd_T = new double[8]; - double[] Dlndjs_Dlnd_P = new double[11]; - double s1 ; - double s2 ; - double ky ; - double kx ; - int i ; - int m = 0; - int n = 0; - //On Error Resume Next VBConversions Warning: On Error Resume Next not supported in C# - Dlndjs_Dlnd_T[0] = -15 + 273.15; - Dlndjs_Dlnd_T[1] = 0 + 273.15; - Dlndjs_Dlnd_T[2] = 15 + 273.15; - Dlndjs_Dlnd_T[3] = 30 + 273.15; - Dlndjs_Dlnd_T[4] = 45 + 273.15; - Dlndjs_Dlnd_T[5] = 60 + 273.15; - Dlndjs_Dlnd_T[6] = 75 + 273.15; - Dlndjs_Dlnd_T[7] = 90 + 273.15; - Dlndjs_Dlnd_P[0] = 0.1F; - Dlndjs_Dlnd_P[1] = 1; - Dlndjs_Dlnd_P[2] = 2; - Dlndjs_Dlnd_P[3] = 3; - Dlndjs_Dlnd_P[4] = 4; - Dlndjs_Dlnd_P[5] = 5; - Dlndjs_Dlnd_P[6] = 6; - Dlndjs_Dlnd_P[7] = 7; - Dlndjs_Dlnd_P[8] = 8; - Dlndjs_Dlnd_P[9] = 9; - Dlndjs_Dlnd_P[10] = 10; - Dlndjs_Dlnd_Data[0][0] = 976; - Dlndjs_Dlnd_Data[1][0] = 1027; - Dlndjs_Dlnd_Data[2][0] = 1071; - Dlndjs_Dlnd_Data[3][0] = 1123; - Dlndjs_Dlnd_Data[4][0] = 1167; - Dlndjs_Dlnd_Data[5][0] = 1213; - Dlndjs_Dlnd_Data[6][0] = 1260; - Dlndjs_Dlnd_Data[7][0] = 1303; - Dlndjs_Dlnd_Data[0][1] = 991; - Dlndjs_Dlnd_Data[1][1] = 1040; - Dlndjs_Dlnd_Data[2][1] = 1082; - Dlndjs_Dlnd_Data[3][1] = 1135; - Dlndjs_Dlnd_Data[4][1] = 1178; - Dlndjs_Dlnd_Data[5][1] = 1224; - Dlndjs_Dlnd_Data[6][1] = 1270; - Dlndjs_Dlnd_Data[7][1] = 1312; - Dlndjs_Dlnd_Data[0][2] = 1014; - Dlndjs_Dlnd_Data[1][2] = 1063; - Dlndjs_Dlnd_Data[2][2] = 1106; - Dlndjs_Dlnd_Data[3][2] = 1153; - Dlndjs_Dlnd_Data[4][2] = 1196; - Dlndjs_Dlnd_Data[5][2] = 1239; - Dlndjs_Dlnd_Data[6][2] = 1281; - Dlndjs_Dlnd_Data[7][2] = 1323; - Dlndjs_Dlnd_Data[0][3] = 1044; - Dlndjs_Dlnd_Data[1][3] = 1091; - Dlndjs_Dlnd_Data[2][3] = 1127; - Dlndjs_Dlnd_Data[3][3] = 1174; - Dlndjs_Dlnd_Data[4][3] = 1216; - Dlndjs_Dlnd_Data[5][3] = 1257; - Dlndjs_Dlnd_Data[6][3] = 1297; - Dlndjs_Dlnd_Data[7][3] = 1338; - Dlndjs_Dlnd_Data[0][4] = 1073; - Dlndjs_Dlnd_Data[1][4] = 1118; - Dlndjs_Dlnd_Data[2][4] = 1149; - Dlndjs_Dlnd_Data[3][4] = 1195; - Dlndjs_Dlnd_Data[4][4] = 1236; - Dlndjs_Dlnd_Data[5][4] = 1275; - Dlndjs_Dlnd_Data[6][4] = 1313; - Dlndjs_Dlnd_Data[7][4] = 1352; - Dlndjs_Dlnd_Data[0][5] = 1114; - Dlndjs_Dlnd_Data[1][5] = 1151; - Dlndjs_Dlnd_Data[2][5] = 1180; - Dlndjs_Dlnd_Data[3][5] = 1224; - Dlndjs_Dlnd_Data[4][5] = 1261; - Dlndjs_Dlnd_Data[5][5] = 1297; - Dlndjs_Dlnd_Data[6][5] = 1333; - Dlndjs_Dlnd_Data[7][5] = 1372; - Dlndjs_Dlnd_Data[0][6] = 1156; - Dlndjs_Dlnd_Data[1][6] = 1185; - Dlndjs_Dlnd_Data[2][6] = 1211; - Dlndjs_Dlnd_Data[3][6] = 1253; - Dlndjs_Dlnd_Data[4][6] = 1287; - Dlndjs_Dlnd_Data[5][6] = 1320; - Dlndjs_Dlnd_Data[6][6] = 1352; - Dlndjs_Dlnd_Data[7][6] = 1391; - Dlndjs_Dlnd_Data[0][7] = 1207; - Dlndjs_Dlnd_Data[1][7] = 1230; - Dlndjs_Dlnd_Data[2][7] = 1250; - Dlndjs_Dlnd_Data[3][7] = 1289; - Dlndjs_Dlnd_Data[4][7] = 1318; - Dlndjs_Dlnd_Data[5][7] = 1346; - Dlndjs_Dlnd_Data[6][7] = 1374; - Dlndjs_Dlnd_Data[7][7] = 1412; - Dlndjs_Dlnd_Data[0][8] = 1261; - Dlndjs_Dlnd_Data[1][8] = 1276; - Dlndjs_Dlnd_Data[2][8] = 1289; - Dlndjs_Dlnd_Data[3][8] = 1324; - Dlndjs_Dlnd_Data[4][8] = 1350; - Dlndjs_Dlnd_Data[5][8] = 1373; - Dlndjs_Dlnd_Data[6][8] = 1396; - Dlndjs_Dlnd_Data[7][8] = 1432; - Dlndjs_Dlnd_Data[0][9] = 1331; - Dlndjs_Dlnd_Data[1][9] = 1331; - Dlndjs_Dlnd_Data[2][9] = 1335; - Dlndjs_Dlnd_Data[3][9] = 1366; - Dlndjs_Dlnd_Data[4][9] = 1385; - Dlndjs_Dlnd_Data[5][9] = 1403; - Dlndjs_Dlnd_Data[6][9] = 1424; - Dlndjs_Dlnd_Data[7][9] = 1456; - Dlndjs_Dlnd_Data[0][10] = 1405; - Dlndjs_Dlnd_Data[1][10] = 1389; - Dlndjs_Dlnd_Data[2][10] = 1383; - Dlndjs_Dlnd_Data[3][10] = 1409; - Dlndjs_Dlnd_Data[4][10] = 1421; - Dlndjs_Dlnd_Data[5][10] = 1435; - Dlndjs_Dlnd_Data[6][10] = 1451; - Dlndjs_Dlnd_Data[7][10] = 1482; - - if (tempT < Dlndjs_Dlnd_T[0]) - { - tempT = Dlndjs_Dlnd_T[0]; - } - if (tempT > Dlndjs_Dlnd_T[7]) - { - tempT = Dlndjs_Dlnd_T[7]; - } - if (tempP_jy < Dlndjs_Dlnd_P[0]) - { - tempP_jy = Dlndjs_Dlnd_P[0]; - } - if (tempP_jy > Dlndjs_Dlnd_P[10]) - { - tempP_jy = Dlndjs_Dlnd_P[10]; - } - - for ( i = 0; i <= 6; i++) - { - if (tempT >= Dlndjs_Dlnd_T[i] && tempT <= Dlndjs_Dlnd_T[i + 1]) - { - m = i; - break; - } - } - - for (i = 0; i <= 9; i++) - { - if (tempP_jy >= Dlndjs_Dlnd_P[i] && tempP_jy <= Dlndjs_Dlnd_P[i + 1]) - { - n = i; - break; - } - } - - if (Dlndjs_Dlnd_P[n + 1] - Dlndjs_Dlnd_P[n] != 0) - { - ky = (tempP_jy - Dlndjs_Dlnd_P[n]) / (Dlndjs_Dlnd_P[n + 1] - Dlndjs_Dlnd_P[n]); - } - else - { - ky = 0; - } - if (Dlndjs_Dlnd_T[m + 1] - Dlndjs_Dlnd_T[m] != 0) - { - kx = (tempT - Dlndjs_Dlnd_T[m]) / (Dlndjs_Dlnd_T[m + 1] - Dlndjs_Dlnd_T[m]); - } - else - { - kx = 0; - } - s1 = Dlndjs_Dlnd_Data[m][n] + (Dlndjs_Dlnd_Data[m][n + 1] - Dlndjs_Dlnd_Data[m][n]) * ky; - s2 = Dlndjs_Dlnd_Data[m + 1][n] + (Dlndjs_Dlnd_Data[m + 1][n + 1] - Dlndjs_Dlnd_Data[m + 1][n]) * ky; - return (s1 + (s2 - s1) * kx) / 100000.0D; - } - //可膨胀系数计算 - private double KePenZhang_JiSuan(double tempP_jy, double tempDp_Pa, double tempZjb, double tempDszs, int JIeliuType, int JiSuanBiaoZhun) // 求可膨胀系数 - { - double returnValue = 0; - //0标准孔板 - //1ISA1932喷嘴 - //2长径喷嘴 - //3文丘里喷嘴 - //4粗铸收缩段经典文丘里管 - //5机械加工收缩段经典文丘里管 - //6粗焊铁板收缩段经典文丘里管 - //7 1/4圆孔板 - - double tuo ; - switch (JIeliuType) - { - case 0: //孔板流量计算 - switch (JiSuanBiaoZhun) - { - case 0: //6143-2004 - tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); - returnValue = 1 - (0.351 + 0.256 * Math.pow(tempZjb, 4) + 0.93 * Math.pow(tempZjb, 8)) * (1 - Math.pow(tuo, (1 / tempDszs))); - break; - case 1: //6143-1996 - returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); - break; - } - break; - case 1: - case 2: - case 3: - case 4: - case 5: - case 6: - switch (JiSuanBiaoZhun) - { - case 0: - //标准喷嘴 iso5167-2002 - tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); - returnValue = Math.pow((((tempDszs * Math.pow(tuo, (2 / tempDszs))) / (tempDszs - 1)) * ((1 - Math.pow(tempZjb, 4)) / (1 - Math.pow(tempZjb, 4) * Math.pow(tuo, (2 / tempDszs)))) * ((1 - Math.pow(tuo, ((tempDszs - 1) / tempDszs))) / (1 - tuo))), 0.5); - break; - - case 1: //iso5167-93 - returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); - break; - } - break; - case 7: //1/4圆孔板 - returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); - break; - case 8: //锥形入口孔板 - tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); - returnValue = 1 - (0.351 + 0.256 * Math.pow(tempZjb, 4) + 0.93 * Math.pow(tempZjb, 8)) * (1 - Math.pow(tuo, (1 / tempDszs))); - - tuo = (tempP_jy - tempDp_Pa) / (tempP_jy ); - - returnValue = 0.5 * (returnValue + Math.pow((((tempDszs * Math.pow(tuo, (2 / tempDszs))) / (tempDszs - 1)) * ((1 - Math.pow(tempZjb, 4)) / (1 - Math.pow(tempZjb, 4) * Math.pow(tuo, (2 / tempDszs)))) * ((1 - Math.pow(tuo, ((tempDszs - 1) / tempDszs))) / (1 - tuo))), 0.5)); - break; - case 9: //偏心孔板 - returnValue = 1 - (0.41 + 0.35 * Math.pow(tempZjb, 4)) * tempDp_Pa / ( tempP_jy * tempDszs); - break; - - } - return returnValue; - } - //流出系数计算 - public final double C_JiSuan(double tempGj, double tempZjb, double tempReD, double tempconQvA, int tempQyfs, int JieLiuType, int JiSuanBiaoZhun) - { - double returnValue = 0; - //流出系数计算函数 - //输入：直径比，雷诺数，取压方式，节流装置类型，计算采用标准 - //输出:流出系数 - - //jieliutype - - //0标准孔板 - //1ISA1932喷嘴 - //2长径喷嘴 - //3文丘里喷嘴 - //4粗铸收缩段经典文丘里管 - //5机械加工收缩段经典文丘里管 - //6粗焊铁板收缩段经典文丘里管 - - double L1 = 0; - double L2 = 0; - switch (JieLiuType) - { - case 0: //孔板 - switch (tempQyfs) - { - case 0: - L1 = 25.4 / tempGj; - L2 = L1; - break; - case 1: - L1 = 0; - L2 = 0; - break; - case 2: - L1 = 1; - L2 = 0.47F; - break; - } - switch (JiSuanBiaoZhun) - { - case 0: //6143-2004 - if (tempGj >= 71.12) - { - returnValue = 0.5961 + 0.0261 * Math.pow(tempZjb, 2) - 0.216 * Math.pow(tempZjb, 8) + 0.000521 * Math.pow((1000000.0D * tempZjb / tempReD), 0.7) + (0.0188 + 0.0063 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * Math.pow(tempZjb, 3.5) * Math.pow((1000000.0D / tempReD), 0.3) + (0.043 + 0.08 * Math.exp(-10 * L1) - 0.123 * Math.exp(-7 * L1)) * (1 - 0.11 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * (Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1))) - 0.031 * (2 * L2 / (1 - tempZjb) - 0.8 * Math.pow((2 * L2 / (1 - tempZjb)), 1.1)) * Math.pow(tempZjb, 1.3); - } - else if (tempGj < 71.12) - { - returnValue = (0.5961 + 0.0261 * Math.pow(tempZjb, 2) - 0.216 * Math.pow(tempZjb, 8) + 0.000521 * Math.pow((1000000.0D * tempZjb / tempReD), 0.7) + (0.0188 + 0.0063 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * Math.pow(tempZjb, 3.5) * Math.pow((1000000.0D / tempReD), 0.3) + (0.043 + 0.08 * Math.exp(-10 * L1) - 0.123 * Math.exp(-7 * L1)) * (1 - 0.11 * Math.pow((19000 * tempZjb / tempReD), 0.8)) * Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1)) - 0.031 * (2 * L2 / (1 - tempZjb) - 0.8 * Math.pow((2 * L2 * (1 - tempZjb)), 1.1)) * Math.pow(tempZjb, 1.3) + 0.011 * (0.75 - tempZjb) * (2.8 - tempGj / 25.4)); - } - break; - case 1: //6143-1996 - if (0.09 * L1 >= 0.039) - { - returnValue = 0.5959 + 0.0312 * Math.pow(tempZjb, (2.1)) - 0.184 * Math.pow(tempZjb, 8) + 0.0029 * Math.pow(tempZjb, 2.5) * Math.pow((1000000.0D / tempReD), 0.75) + 0.039 * Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1)) - 0.0337 * L1 * Math.pow(tempZjb, 3); - } - else if (0.09 * L1 < 0.039) - { - returnValue = (0.5959 + 0.0312 * Math.pow(tempZjb, (2.1)) - 0.184 * Math.pow(tempZjb, 8) + 0.0029 * Math.pow(tempZjb, 2.5) * Math.pow((1000000.0D / tempReD), 0.75) + 0.09 * L1 * Math.pow(tempZjb, 4) * Math.pow((1 - Math.pow(tempZjb, 4)), (-1)) - 0.0337 * L1 * Math.pow(tempZjb, 3)); - } - break; - } - break; - case 1: //ISA1932喷嘴 - returnValue = (0.99 - 0.2262 * Math.pow(tempZjb, 4.1) - (0.00175 * Math.pow(tempZjb, 2) - 0.0033 * Math.pow(tempZjb, 4.15)) * (1000000.0D / Math.pow(tempReD, 1.15))); - break; - - case 2: //长径喷嘴 - returnValue = (0.9965 - 0.00653 * Math.pow(tempZjb, 0.5) * Math.pow((1000000.0D / tempReD), 0.5)); - break; - case 3: //文丘里喷嘴 - returnValue = (0.9858 - 0.196 * Math.pow(tempZjb, 4.5)); - break; - case 4: //粗铸收缩段经典文丘里管 - returnValue = (0.984F); - break; - case 5: //机械加工收缩段经典文丘里管 - returnValue = (0.995F); - break; - case 6: //粗焊铁板收缩段经典文丘里管 - returnValue = (0.985F); - break; - case 7: //1/4圆孔板 - returnValue = 0.73823 - 0.3309 * tempZjb - 1.1615 * Math.pow(tempZjb, 2) + 1.5084 * Math.pow(tempZjb, 3); - break; - case 8: //锥形入口孔板 - returnValue = (0.734F); - break; - case 9: //偏心孔板 - returnValue = 0.9355 - 1.6889 * tempZjb + 3.0428 * Math.pow(tempZjb, 2) - 1.7989 * Math.pow(tempZjb, 3); - break; - } - double tempRed1 ; - switch (JieLiuType) - { - case 0: //孔板流量计算 - tempRed1 = tempconQvA * returnValue; - switch (JieLiuType) - { - case 0: - if (tempRed1 < (170 * Math.pow(tempZjb, 2) * tempGj)) - { - JOptionPane.showMessageDialog(null, "雷诺数超过标准孔板的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - return returnValue; - } - break; - case 1: - if (tempZjb >= 0.1 & tempZjb <= 0.56) - { - if (tempRed1 < 5000) - { - JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - return returnValue; - } - } - if (tempZjb > 0.56) - { - if (tempRed1 < (16000 * Math.pow(tempZjb, 2) * tempGj)) - { - JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - return returnValue; - } - } - break; - } - break; - - case 1: //标准喷嘴 - tempRed1 = tempconQvA * returnValue; - if (tempZjb >= 0.3 & tempZjb < 0.44) - { - if (tempRed1 < 70000 | tempRed1 > 10000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - } - } - if (tempZjb >= 0.44 & tempZjb < 0.8) - { - if (tempRed1 < 20000 | tempRed1 > 10000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过标准喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - } - } - break; - case 2: //长径喷嘴 - tempRed1 = tempconQvA * returnValue; - if (tempRed1 < 10000.0D | tempRed1 > 10000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过长径喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - } - break; - case 3: //文丘里喷嘴 - tempRed1 = tempconQvA * returnValue; - if (tempRed1 < 150000.0D | tempRed1 > 2000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过文丘里喷嘴的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - } - break; - case 4: //粗铸收缩段经典文丘里管 - tempRed1 = tempconQvA * returnValue; - if (tempRed1 < 200000.0D | tempRed1 > 2000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过粗铸收缩段经典文丘里管的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - return returnValue; - } - break; - case 5: //机械加工收缩段经典文丘里管 - tempRed1 = tempconQvA * returnValue; - if (tempRed1 < 200000.0D | tempRed1 > 1000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过机械加工收缩段经典文丘里管的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - } - break; - case 6: //粗焊铁板收缩段经典文丘里管 - tempRed1 = tempconQvA * returnValue; - if (tempRed1 < 200000.0D | tempRed1 > 2000000.0D) - { - JOptionPane.showMessageDialog(null, "雷诺数超过粗焊铁板收缩段经典文丘里管的使用范围！停止计算！", "提示", JOptionPane.PLAIN_MESSAGE); - } - break; - - //标准孔板 - //ISA1932喷嘴 - //长径喷嘴 - //文丘里喷嘴 - //粗铸收缩段经典文丘里管 - //机械加工收缩段经典文丘里管 - //粗焊铁板收缩段经典文丘里管 - - - - - } - return returnValue; - } - - private double BkTable(double temPrk, double tempKj, int tempBkjsff) { - double[] BkTable_x = new double[10]; - double[] BkTable_Y = new double[10]; - //On Error Resume Next VBConversions Warning: On Error Resume Next not supported in C# - if (tempBkjsff == 1) { - double tempRkBiKj = temPrk / tempKj; - // static double[] x = new double[10]; //VBConversions Note: Static variable moved to class level and renamed BkTable_x. Local static variables are not supported in C#. - // static double[] Y = new double[10]; //VBConversions Note: Static variable moved to class level and renamed BkTable_Y. Local static variables are not supported in C#. - int i; - int xIndex = 0; - //If x(0) = 0 Then - BkTable_x[0] = 0.0004F; - BkTable_x[1] = 0.001F; - BkTable_x[2] = 0.002F; - BkTable_x[3] = 0.004F; - BkTable_x[4] = 0.006F; - BkTable_x[5] = 0.008F; - BkTable_x[6] = 0.01F; - BkTable_x[7] = 0.012F; - BkTable_x[8] = 0.014F; - BkTable_x[9] = 0.015F; - - BkTable_Y[0] = 1; - BkTable_Y[1] = 1.005F; - BkTable_Y[2] = 1.012F; - BkTable_Y[3] = 1.022F; - BkTable_Y[4] = 1.032F; - BkTable_Y[5] = 1.04F; - BkTable_Y[6] = 1.048F; - BkTable_Y[7] = 1.055F; - BkTable_Y[8] = 1.062F; - BkTable_Y[9] = 1.065F; - //End If - - if (tempRkBiKj <= 0.0004) { - return 1; - } - if (tempRkBiKj > 0.015) { - return 1.065F; - } - for (i = 0; i <= 8; i++) { - if (tempRkBiKj >= BkTable_x[i] && tempRkBiKj <= BkTable_x[i + 1]) { - xIndex = i; - break; - } - } - return BkTable_Y[xIndex] + (tempRkBiKj - BkTable_x[xIndex]) * (BkTable_Y[xIndex + 1] - BkTable_Y[xIndex]) / (BkTable_x[xIndex + 1] - BkTable_x[xIndex]); - } - return 0; - } - //管道粗糙度计算 - private double CcdXsjs(double tempPipeType, double tempGj, double tempZjb, double TempRed) - { - double returnValue ; - //粗糙度系数计算 - double Jdccd = 0; //绝对粗糙度 - double Xdccd ; //相对粗糙度 - //Dim CcdXs As single - double s1 ; - double s2 ; - double ky ; - double kx ; - int i ; - int m = 0; - int n = 0; - if (tempPipeType == 0) - { - Jdccd = 0.029F;} - else if (tempPipeType == 1) - { - Jdccd = 0.075F; - } - else if (tempPipeType == 2) - { - Jdccd = 0.075F; - } - else if (tempPipeType == 3) - { - Jdccd = 0.075F; - } - else if (tempPipeType == 4) - { - Jdccd = 0.1F; - } - else if (tempPipeType == 5) - { - Jdccd = 0.15F; - } - else if (tempPipeType == 6) - { - Jdccd = 1; - } - else if (tempPipeType == 7) - { - Jdccd = 2.1F; - } - else if (tempPipeType == 8) - { - Jdccd = 0.04F; - } - else if (tempPipeType == 9) - { - Jdccd = 0.15F; - } - else if (tempPipeType == 10) - { - Jdccd = 0.13F; - } - else if (tempPipeType == 11) - { - Jdccd = 0.25F; - } - Xdccd = tempGj / Jdccd; - if (Xdccd < 400) - { - JOptionPane.showMessageDialog(null, "粗糙度取得太高,粗略计算", "提示", JOptionPane.PLAIN_MESSAGE); - Xdccd = 400; - } - if (Xdccd >= 3400) - { - Xdccd = 3400; - } - - if (Xdccd < 3200 & Math.pow(tempZjb, 2) > 0.1 & Math.pow(tempZjb, 2) < 0.64) - { - - int[] Xdccdb = new int[10]; - double[] Btf = new double[8]; - double[][] CcdXsb = new double[10][8]; - - Xdccdb[0] = 400; - Xdccdb[1] = 800; - Xdccdb[2] = 1200; - Xdccdb[3] = 1600; - Xdccdb[4] = 2000; - Xdccdb[5] = 2400; - Xdccdb[6] = 2800; - Xdccdb[7] = 3200; - Xdccdb[8] = 3400; - - Btf[0] = 0.1F; - Btf[1] = 0.2F; - Btf[2] = 0.3F; - Btf[3] = 0.4F; - Btf[4] = 0.5F; - Btf[5] = 0.6F; - Btf[6] = 0.64F; - - CcdXsb[0][0] = 1.002; - CcdXsb[1][0] = 1; - CcdXsb[2][0] = 1; - CcdXsb[3][0] = 1; - CcdXsb[4][0] = 1; - CcdXsb[5][0] = 1; - CcdXsb[6][0] = 1; - CcdXsb[7][0] = 1; - CcdXsb[8][0] = 1; - CcdXsb[0][1] = 1.003; - CcdXsb[1][1] = 1.002; - CcdXsb[2][1] = 1.001; - CcdXsb[3][1] = 1; - CcdXsb[4][1] = 1; - CcdXsb[5][1] = 1; - CcdXsb[6][1] = 1; - CcdXsb[7][1] = 1; - - CcdXsb[0][2] = 1.006; - CcdXsb[1][2] = 1.004; - CcdXsb[2][2] = 1.002; - CcdXsb[3][2] = 1.001; - CcdXsb[4][2] = 1; - CcdXsb[5][2] = 1; - CcdXsb[6][2] = 1; - CcdXsb[7][2] = 1; - - CcdXsb[0][3] = 1.009; - CcdXsb[1][3] = 1.006; - CcdXsb[2][3] = 1.004; - CcdXsb[3][3] = 1.002; - CcdXsb[4][3] = 1.001; - CcdXsb[5][3] = 1; - CcdXsb[6][3] = 1; - CcdXsb[7][3] = 1; - - CcdXsb[0][4] = 1.014; - CcdXsb[1][4] = 1.009; - CcdXsb[2][4] = 1.006; - CcdXsb[3][4] = 1.004; - CcdXsb[4][4] = 1.002; - CcdXsb[5][4] = 1.001; - CcdXsb[6][4] = 1; - CcdXsb[7][4] = 1; - - CcdXsb[0][5] = 1.02; - CcdXsb[1][5] = 1.013; - CcdXsb[2][5] = 1.009; - CcdXsb[3][5] = 1.006; - CcdXsb[4][5] = 1.003; - CcdXsb[5][5] = 1.002; - CcdXsb[6][5] = 1; - CcdXsb[7][5] = 1; - - CcdXsb[0][6] = 1.024; - CcdXsb[1][6] = 1.016; - CcdXsb[2][6] = 1.011; - CcdXsb[3][6] = 1.007; - CcdXsb[4][6] = 1.004; - CcdXsb[5][6] = 1.002; - CcdXsb[6][6] = 1.002; - CcdXsb[7][6] = 1; - - - for (i = 0; i <= 8; i++) - { - if (Xdccd >= Xdccdb[i] && Xdccd <= Xdccdb[i + 1]) - { - m = i; - break; - } - } - for (i = 0; i <= 6; i++) - { - if (Math.pow(tempZjb, 2) >= Btf[i] && Math.pow(tempZjb, 2) <= Btf[i + 1]) - { - n = i; - break; - } - } - - ky = (Btf[n + 1] - Btf[n] != 0) ? ((Math.pow(tempZjb, 2) - Btf[n]) / (Btf[n + 1] - Btf[n])) : 0; - - kx = (Xdccdb[m + 1] - Xdccdb[m] != 0) ? ((Xdccd - Xdccdb[m]) / (Xdccdb[m + 1] - Xdccdb[m])) : 0; - - s1 = CcdXsb[m][n] + (CcdXsb[m][n + 1] - CcdXsb[m][n]) * ky; - s2 = CcdXsb[m + 1][n] + (CcdXsb[m + 1][n + 1] - CcdXsb[m + 1][n]) * ky; - returnValue = s1 + (s2 - s1) * kx; - - returnValue = TempRed > 1000000.0 ? returnValue : ((returnValue - 1) * Math.pow((Math.log10(TempRed) / 2), 2) + 1); - - } - else - { - returnValue = 1; - } - return returnValue; - } // ######################################################################### // ######################'流量转换 标况转工况############################## diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GBT22634WaterContentCalc.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GBT22634WaterContentCalc.java new file mode 100644 index 0000000..f8b150c --- /dev/null +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GBT22634WaterContentCalc.java @@ -0,0 +1,152 @@ +package com.ruoyi.ngCalTools.controller; +import com.ruoyi.ngCalTools.model.GasProps; +import com.ruoyi.system.controller.UnitConvert; +import com.ruoyi.system.service.ISysUnitConvertService; +import com.ruoyi.system.service.impl.SysUnitConvertServiceImpl; +import org.apache.commons.math3.analysis.UnivariateFunction; +import org.apache.commons.math3.analysis.solvers.BracketingNthOrderBrentSolver; +import org.apache.commons.math3.analysis.solvers.NewtonRaphsonSolver; +import org.apache.commons.math3.exception.TooManyEvaluationsException; + +/** + * GB/T 22634-2008 天然气水含量与水露点工业级换算 + * 注：需配合气体性质计算模块GasController使用 + */ + +public class GBT22634WaterContentCalc { + + + + // 安托因方程常数（GB/T 22634规定） + private static final double ANTOINE_A = 8.07131; + private static final double ANTOINE_B = 1730.63; + private static final double ANTOINE_C = 233.426; + private static final double MMHG_TO_KPA = 0.1333223684; + + /** + * 工业级水含量计算（含真实气体修正） + * @param dewPointTemp 水露点温度（℃） + * @param pressure 绝对压力（kPa） + * @param gasProps 气体组分属性 + * @return 水含量（mg/Sm³） + */ + public double calculateWaterContent(double dewPointTemp, + double pressure, + GasProps gasProps) { + // 1. 计算饱和水蒸气压 + double pWater = calculateWaterVaporPressure(dewPointTemp); + + // 2. 获取压缩因子 + ISysUnitConvertService iSysUnitConvertService =new SysUnitConvertServiceImpl(); + UnitConvert unitConvert=new UnitConvert(iSysUnitConvertService); + GasController gasController=new GasController(unitConvert); + gasController.Crit(gasProps,0); // 调用实际压缩因子计算模块 + + // 3. 计算真实摩尔体积 + double molarVolume = calculateMolarVolume(dewPointTemp, pressure, gasProps.getdZf()); + + // 4. 计算水含量 + return (pWater / (pressure - pWater)) * (18.01528 * 1e6) / molarVolume; + } + + /** + * 逆向计算水露点温度（工业级精度） + * @param waterContent 水含量（mg/Sm³） + * @param pressure 绝对压力（kPa） + * @param gasProps 气体组分属性 + * @return 水露点温度（℃），精度±0.01℃ + */ + public double calculateDewPoint(final double waterContent, + final double pressure, + final GasProps gasProps) { +// NewtonRaphsonSolver solver = new NewtonRaphsonSolver(1e-4); + + // 使用5阶布伦特法求解器（无需导数） + BracketingNthOrderBrentSolver solver = new BracketingNthOrderBrentSolver( + 1e-4, // 相对精度 + 1e-6, // 绝对精度 + 1e-10, // 函数值精度 + 5 // 多项式阶数 + ); + + UnivariateFunction function = new UnivariateFunction() { + @Override + public double value(double tempC) { + try { + // 创建临时气体属性副本 + GasProps tempProps = gasProps.clone(); + tempProps.setdTf(tempC + 273.15); // 更新温度 + tempProps.setdPf(pressure); + + double calc = calculateWaterContent(tempC, pressure, tempProps); + return calc - waterContent; + } catch (Exception e) { + return Double.NaN; + } + } + }; + + try { + return solver.solve(100, function, -50.0, 100.0, 20.0); // 初始猜测20℃ + } catch (TooManyEvaluationsException e) { + throw new ArithmeticException("露点计算未收敛，请检查输入参数"); + } + } + + // 饱和水蒸气压计算（GB/T 22634规定方法） + private double calculateWaterVaporPressure(double tempC) { + if(tempC < 0 || tempC > 100) { + throw new IllegalArgumentException("温度超出GB/T 22634适用范围"); + } + return Math.pow(10, ANTOINE_A - ANTOINE_B / (tempC + ANTOINE_C)) * MMHG_TO_KPA; + } + + // 真实气体摩尔体积计算 + private static double calculateMolarVolume(double tempC, double pressure, double Z) { + return Z * 8.3144621 * (tempC + 273.15) / pressure; + } + +// // 气体属性类（示例结构） +// public static class GasProps implements Cloneable { +// private double temperature; // K +// private double pressure; // kPa +// private double[] composition; // 组分摩尔分数 +// +// // 克隆方法用于迭代计算 +// @Override +// public GasProps clone() { +// GasProps clone = new GasProps(); +// clone.temperature = this.temperature; +// clone.pressure = this.pressure; +// clone.composition = this.composition.clone(); +// return clone; +// } +// +// // Getter/Setter +// public double getTemperature() { return temperature; } +// public void setTemperature(double temperature) { this.temperature = temperature; } +// public double getPressure() { return pressure; } +// public void setPressure(double pressure) { this.pressure = pressure; } +// public double[] getComposition() { return composition; } +// public void setComposition(double[] composition) { this.composition = composition; } +// } +// +// // 示例调用 +// public static void main(String[] args) { +// // 初始化气体属性 +// GasProps gas = new GasProps(); +// gas.setPressure(5000.0); // 5 MPa +// gas.setTemperature(293.15); // 20℃ +// gas.setComposition(new double[]{0.95, 0.03, 0.02}); // 示例组分 +// +// // 正向计算 +// double wc = calculateWaterContent(10.0, 5000.0, gas); +// System.out.printf("5MPa下10℃露点对应水含量：%.2f mg/Sm³\n", wc); +// +// // 逆向计算 +// double dewPoint = calculateDewPoint(1000.0, 5000.0, gas); +// System.out.printf("5MPa下1000mg/Sm³对应露点：%.2f℃\n", dewPoint); +// } + + +} diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GasController.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GasController.java index 08f64ab..d70ba06 100644 --- a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GasController.java +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/GasController.java @@ -46,6 +46,31 @@ public class GasController { } gasProps.dTf = tempTf; gasProps.dCbtj = flowProps.getdCbtj(); + + switch (gasProps.dCbtj) + { + case 2: + gasProps.setdPb(101325); + gasProps.setdTb( 273.15); + flowProps.setdPb_M(101325); + flowProps.setdTb_M(273.15); + break; + case 1: + + gasProps.setdPb(101325); + gasProps.setdTb( 288.15); + flowProps.setdPb_M(101325); + flowProps.setdTb_M(288.15); + break; + case 0: + + gasProps.setdPb(101325); + gasProps.setdTb( 293.15); + flowProps.setdPb_M(101325); + flowProps.setdTb_M(293.15); + break; + } + String[] stringArray = flowProps.getdngComponents().split("_"); double[] doubleArray = new double[stringArray.length]; // 遍历字符串数组，将每个元素转换为 double 类型 for (int i = 0; i < stringArray.length; i++) { @@ -67,9 +92,7 @@ public class GasController { thermService = new ThermService(); detailService = new DetailService(); gbt11062Service = new GBT11062Service(); - - Crit(gasProps, 0); //计算临界流函数所有参数都计算了 - + Crit(gasProps, 0); //计算临界流函数所有参数都计算了 } public int NG_Cal_Init() { @@ -115,23 +138,7 @@ public class GasController { } } - switch (gasProps.dCbtj) - { - case 2: - gasProps.dPb = 101325; - gasProps.dTb = 273.15; - break; - case 1: - gasProps.dPb = 101325; - gasProps.dTb = 288.15; - break; - case 0: - gasProps.dPb = 101325; - gasProps.dTb = 293.15; - break; - - } //begin by calculating densities and thermodynamic properties thermService.Run(gasProps, detailService); diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/NozellFlowCalc.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/NozellFlowCalc.java new file mode 100644 index 0000000..addceb5 --- /dev/null +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/controller/NozellFlowCalc.java @@ -0,0 +1,236 @@ +package com.ruoyi.ngCalTools.controller; + +import com.ruoyi.ngCalTools.model.FlowProps; +import com.ruoyi.ngCalTools.model.GasProps; +import com.ruoyi.ngCalTools.service.GBT11062Service; +import com.ruoyi.ngCalTools.service.ISO9300Service; + +import static com.ruoyi.ngCalTools.controller.DpFlowCalc.Dlndjs; +import static com.ruoyi.ngCalTools.controller.FlowController.FlowConvert_WorkToBase; + +public class NozellFlowCalc { + + + // 主计算入口（支持多喷嘴类型） + + public static void executeFullCalculation(FlowProps flow, GasProps gas) { + // 步骤1：物性计算 + calculateAccurateProperties( gas); + + // 步骤2：粘度计算 + ISO9300Service.calculateViscosity(gas); + + // 步骤3：喷嘴类型专项计算 + calculate(flow, gas); + + // 步骤4：保存能量计算结果 +// calculateCalorificValue(gas); + } + public static void calculateAccurateProperties( GasProps gasProps) { + // ISO 9300专用参数计算 + gasProps.setdZ_ISO9300( gasProps.dZf * 0.998); // 示例修正系数 + gasProps.setdKappa_ISO( calculateISOKappa(gasProps)); + } + + private static double calculateISOKappa(GasProps gasProps) { + // 根据ISO 20765:2005计算 + return 1.3 + 0.02*(gasProps.dMrx - 16)/10; + } + + public static void calculate(FlowProps flow, GasProps gas) { + // 前置校验 +// validateBetaRatio(flow); + + // 获取基础参数 + double D = flow.getdPipeD() / 1000; + double d = flow.getdOrificeD() / 1000; + double beta = d / D; + + // 选择喷嘴处理分支 + switch ((int) flow.getdNozzleType()) { + case FlowProps.NOZZLE_LONG_RADIUS: + handleLongRadiusNozzle(flow, gas, beta); + break; + case FlowProps.NOZZLE_VENTURI: + handleVenturiNozzle(flow, gas, beta); + break; + case FlowProps.NOZZLE_CYLINDRICAL: + handleCylindricalThroat(flow, gas, beta); + break; + } + } + + + // 长径喷嘴处理 + private static void handleLongRadiusNozzle(FlowProps flow, GasProps gas, double beta) { + // 检查曲率半径是否达标（ISO 9300-5.2.2） + if(flow.getdUpstreamR() < 0.15*flow.getdPipeD()) { + throw new IllegalArgumentException("上游曲率半径不满足长径喷嘴要求"); + } + + // 流出系数计算（ISO 9300公式） + double ReD = calculateReynolds(flow, gas); + double Cd = 0.9965 - 0.00653*Math.pow(beta, 5) + + (0.041*Math.pow(beta, 3))/Math.sqrt(ReD); + + // 特殊膨胀系数 + double epsilon = 1 - (0.351 + 0.256*Math.pow(beta,4) + + 0.93*Math.pow(beta,8)) * (flow.getdDp()/flow.getdPf()); + + calculateFinalFlow(flow, gas, Cd, epsilon); + } + + // 文丘里喷嘴处理 + private static void handleVenturiNozzle(FlowProps flow, GasProps gas, double beta) { + // 喉部直径校验 + double throatD =flow.getdThroatD()/1000; + if(throatD <= 0.05*flow.getdPipeD()) { + throw new IllegalArgumentException("喉部直径过小不符合文丘里喷嘴规范"); + } + + // 流出系数计算（ISO 9300:2022 Annex B） + double Cd = 0.985 - 0.195*Math.pow(beta, 2.5) + + 0.034*Math.pow(beta, 5)*Math.log(flow.getdRnPipe()); + + // 膨胀系数修正 + double epsilon = 1 - (0.484 + 1.189*Math.pow(beta,4)) + * Math.pow(flow.getdDp()/flow.getdPf(), 1.2); + + calculateFinalFlow(flow, gas, Cd, epsilon); + } + + // 雷诺数计算（含湍流修正） + private static double calculateReynolds(FlowProps flow, GasProps gas) { + double velocity = flow.getdVelocityFlow(); + double D = flow.getdPipeD()/1000; + return (gas.getdRhof() * velocity * D) / gas.getdMu(); + } + + // 最终流量计算（统一处理） + private static void calculateFinalFlow(FlowProps flow, GasProps gas, + double Cd, double epsilon) { + double d = flow.getdOrificeD()/1000; + double deltaP = flow.getdDp()/1000; + double rho = gas.getdRhof(); + + double Q = Cd * (Math.PI * d*d /4) + * Math.sqrt(2*deltaP/rho) + * epsilon + * gas.getdZ_ISO9300(); // ISO专用压缩因子修正 + + flow.setdVFlowf(Q); + flow.setdVFlowb(FlowConvert_WorkToBase(flow,gas)); + } + + + private static void handleCylindricalThroat(FlowProps flow, GasProps gas, double beta) { + /*============================= + * 1. 参数校验 (ISO 9300:2022 Section 6.3) + ==============================*/ + // 校验β范围 + if (beta < 0.4 || beta > 0.7) { + throw new IllegalArgumentException("圆柱喉部喷嘴Beta比值超出有效范围(0.4 ≤ β ≤ 0.7)"); + } + + // 校验喉部长度 (L = 0.6d ~ 0.8d) + double throatLength = flow.getdThroatLength(); // 喉部长度（需在FlowProps中添加字段） + double d = flow.getdThroatD()/1000; + if (throatLength < 0.6*d || throatLength > 0.8*d) { + throw new IllegalArgumentException("喉部长度需满足0.6d ≤ L ≤ 0.8d"); + } + + // 校验表面粗糙度Ra ≤ 0.4 μm + if (flow.getdThroatRoughness() > 0.4e-6) { // 需添加dThroatRoughness字段 + throw new IllegalArgumentException("喉部表面粗糙度需≤0.4μm"); + } + + /*============================= + * 2. 流出系数计算 (ISO 9300 Annex C) + ==============================*/ + double ReD = calculateReynoldsForCylindrical(flow, gas); + double Cd; + + // 根据雷诺数范围选择公式 + if (ReD >= 2e5 && ReD <= 2e6) { + Cd = 0.9857 - 0.195 * Math.pow(beta, 4.5) + + (1.23e5 * Math.pow(beta, 3)) / ReD; + } else if (ReD > 2e6) { + Cd = 0.996 - 0.00653 * Math.pow(beta, 5) + + 0.031 * Math.pow(beta, 2.5) * Math.log10(ReD/1e6); + } else { + throw new IllegalArgumentException("雷诺数超出有效范围(ReD ≥ 2×10⁵)"); + } + + // 粗糙度修正（ISO 9300 C.2.3） + double roughnessFactor = 1 - 0.012 * (flow.getdThroatRoughness()/0.4e-6); + Cd *= roughnessFactor; + + /*============================= + * 3. 膨胀系数计算 + ==============================*/ + double pressureRatio = flow.getdDp() / flow.getdPf(); + double epsilon = 1 - (0.41 + 0.35 * Math.pow(beta, 4)) * pressureRatio + / gas.getdKappa_ISO(); + + // 高速流修正 (Ma > 0.25) + if (flow.getdMachNumber() > 0.25) { // 需添加Mach数计算 + epsilon *= 1 - 0.05 * Math.pow(flow.getdMachNumber(), 2); + } + + /*============================= + * 4. 雷诺数特殊处理 + ==============================*/ + // 使用喉部直径计算雷诺数 + double ReThroat = (gas.getdRhof() * flow.getdVelocityFlow() * d) + / gas.getdMu(); + flow.setdRnPipe(ReThroat); // 更新雷诺数 + + /*============================= + * 5. 执行最终计算 + ==============================*/ + calculateFinalFlow(flow, gas, Cd, epsilon); + } + + // 圆柱喉部专用雷诺数计算 + private static double calculateReynoldsForCylindrical(FlowProps flow, GasProps gas) { + double Q = flow.getdVFlowf(); // 临时工况流量 + double d = flow.getdThroatD()/1000; + double A = Math.PI * d * d / 4; + double u = Q / A; // 喉部流速 + return (gas.getdRhof() * u * d) / gas.getdMu(); + } + + // 马赫数计算（需在FlowProps中添加dMachNumber字段） + private double calculateMachNumber(FlowProps flow, GasProps gas) { + double velocity = flow.getdVelocityFlow(); + double sos = gas.getdSOS(); // 声速 + return velocity / sos; + } + + + public static double convertRoughness(double value, int unit) { + switch(unit) { + case 1: return value * 1e-6; // μm → m + case 2: return value * 1.638e-5; // μin → m + default: return value; + } + } + + private void checkHighSpeedEffect(FlowProps flow, GasProps gas) { + double Ma = calculateMachNumber(flow, gas); + flow.setdMachNumber(Ma); + if (Ma > 0.25) { + System.out.println("警告：马赫数超过0.25，已启用高速修正"); + } + } + public static void validateThroatGeometry(double beta, double L, double d) { + if (L/d < 0.6 || L/d > 0.8) { + throw new IllegalArgumentException( + String.format("喉部长度比例异常(L/d=%.2f)，需满足0.6≤L/d≤0.8", L/d)); + } + if (beta < 0.4 || beta > 0.7) { + throw new IllegalArgumentException( + String.format("β=%.3f超出圆柱喉部有效范围", beta)); + } + } +} diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/FlowProps.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/FlowProps.java index 5f12f91..afc2232 100644 --- a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/FlowProps.java +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/FlowProps.java @@ -516,6 +516,54 @@ public class FlowProps { this.dKappa = dKappa; } + public double getdFpv() { + return dFpv; + } + + public void setdFpv(double dFpv) { + this.dFpv = dFpv; + } + public double getdNozellCdModel() { + return dNozellCdModel; + } + + public void setdNozellCdModel(double dNozellCdModel) { + this.dNozellCdModel = dNozellCdModel; + } + + public double getdNozellReDCorrectionType() { + return dNozellReDCorrectionType; + } + + public void setdNozellReDCorrectionType(double dNozellReDCorrectionType) { + this.dNozellReDCorrectionType = dNozellReDCorrectionType; + } + + + public double getdNozzleType() { + return dNozzleType; + } + + public void setdNozzleType(double dNozzleType) { + this.dNozzleType = dNozzleType; + } + + public double getdThroatD() { + return dThroatD; + } + + public void setdThroatD(double dThroatD) { + this.dThroatD = dThroatD; + } + + public double getdUpstreamR() { + return dUpstreamR; + } + + public void setdUpstreamR(double dUpstreamR) { + this.dUpstreamR = dUpstreamR; + } + private int dZcalbz; // 压缩因子计算标准 private int dCbtj; // 计量参比条件压力 private double dPb_M; // 计量参比条件压力 @@ -569,13 +617,7 @@ public class FlowProps { private double dFG; // 求相对密度系数 FG private double dFT; // 求流动温度系数 FT - public double getdFpv() { - return dFpv; - } - public void setdFpv(double dFpv) { - this.dFpv = dFpv; - } private double dFpv; // 求超压缩因子 Fpv @@ -595,9 +637,54 @@ public class FlowProps { private double dBeta; // 直径比 private double dKappa; // 等熵指数 + private double dNozzleType; // 0表示长径喷嘴，1表示文丘里喷嘴 2 圆柱形喉部文丘里喷嘴 + + + + private double dThroatD; // 文丘里喉部直径（单位同dOrificeD） + private double dUpstreamR; // 上游曲率半径（长径喷嘴专用） + + private double dNozellCdModel; // 流出系数模型//0 ISO 9300标准模型 //1 Humble et al.模型 //2 用户自定义值 + private double dNozellReDCorrectionType; // 流出系数模型//0 不修正 //1 ISO9300建议方法修正模型 //2 Bendick // Getters and Setters + // 喷嘴类型枚举 (ISO 9300-2022) + public static final int NOZZLE_LONG_RADIUS = 0; // 长径喷嘴 + public static final int NOZZLE_VENTURI = 1; // 文丘里喷嘴 + public static final int NOZZLE_CYLINDRICAL = 2; // 圆柱形喉部文丘里喷嘴 + + // 在FlowProps类中添加 + private double dThroatLength; // 喉部长度 + + public double getdThroatLength() { + return dThroatLength; + } + + public void setdThroatLength(double dThroatLength) { + this.dThroatLength = dThroatLength; + } + + public double getdThroatRoughness() { + return dThroatRoughness; + } + + public void setdThroatRoughness(double dThroatRoughness) { + this.dThroatRoughness = dThroatRoughness; + } + + public double getdMachNumber() { + return dMachNumber; + } + + public void setdMachNumber(double dMachNumber) { + this.dMachNumber = dMachNumber; + } + + private double dThroatRoughness; // 表面粗糙度 (m) + private double dMachNumber; // 马赫数 + +// 补充getter/setter } diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/GasProps.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/GasProps.java index 6eb9149..655aab4 100644 --- a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/GasProps.java +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/model/GasProps.java @@ -1,6 +1,19 @@ package com.ruoyi.ngCalTools.model; -public class GasProps { +public class GasProps implements Cloneable { + + // 重写 clone() 方法 + @Override + public GasProps clone() { + try { + return (GasProps) super.clone(); + } catch (CloneNotSupportedException e) { + // 因为已经实现了 Cloneable 接口，所以不会抛出该异常 + throw new AssertionError(); + } + } + + // corresponds to the control group in meter classes public int lStatus; // calculation status 计算状态 public boolean bForceUpdate; // 执行全部计算的标志 signal to perform full calculation @@ -508,6 +521,43 @@ public class GasProps { public double dC3C4; // C3C4组分含量 (kg/m3) public String dngComponents; //组分的组合字符串 从前端传过来 + // 新增高精度物性参数 + private double dMu; // 动态粘度 (Pa·s) + private double dNu; // 运动粘度 (m²/s) + public double getdMu() { + return dMu; + } + + public void setdMu(double dMu) { + this.dMu = dMu; + } + + public double getdNu() { + return dNu; + } + + public void setdNu(double dNu) { + this.dNu = dNu; + } + + public double getdZ_ISO9300() { + return dZ_ISO9300; + } + + public void setdZ_ISO9300(double dZ_ISO9300) { + this.dZ_ISO9300 = dZ_ISO9300; + } + + public double getdKappa_ISO() { + return dKappa_ISO; + } + + public void setdKappa_ISO(double dKappa_ISO) { + this.dKappa_ISO = dKappa_ISO; + } + + private double dZ_ISO9300;// ISO 9300专用压缩因子 + private double dKappa_ISO;// ISO 9300专用等熵指数 } diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/GBT11062Service.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/GBT11062Service.java index 9afe7f4..faeb03b 100644 --- a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/GBT11062Service.java +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/GBT11062Service.java @@ -10,7 +10,7 @@ public class GBT11062Service { int[] aiCID = new int[21];// component IDs double[] dXi = new double[21];// mole fraction of component i // 初始化 adTableMri 数组 - double[] adTableMri = { 16.0430, 28.0135, 44.0100, 30.0700, 44.0970, 18.0153, 34.0820, 2.0159, 28.0100, 31.9988, 58.1230, 58.1230, 72.1500, 72.1500, 86.1770, 100.2040, 114.2310, 128.2580, 142.2850, 4.0026, 39.9480 }; + static double[] adTableMri = { 16.0430, 28.0135, 44.0100, 30.0700, 44.0970, 18.0153, 34.0820, 2.0159, 28.0100, 31.9988, 58.1230, 58.1230, 72.1500, 72.1500, 86.1770, 100.2040, 114.2310, 128.2580, 142.2850, 4.0026, 39.9480 }; // 初始化 adTablePc 数组 double[] adTablePc = { 4.604, 3.399, 7.382, 4.88, 4.249, 22.118, 9.005, 1.297, 3.499, 5.081, 3.648, 3.797, 3.381, 3.369, 3.012, 2.736, 2.486, 0, 0, 0.2275, 4.876 }; // 初始化 adTableTc 数组 diff --git a/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/ISO9300Service.java b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/ISO9300Service.java new file mode 100644 index 0000000..ecbb091 --- /dev/null +++ b/ruoyi-models/ruoyi-ngtools/src/main/java/com/ruoyi/ngCalTools/service/ISO9300Service.java @@ -0,0 +1,174 @@ +package com.ruoyi.ngCalTools.service; + +import com.ruoyi.ngCalTools.model.GasProps; + +public class ISO9300Service { + + // 基于Chapman-Enskog理论的计算 高精度粘度计算 音速喷嘴ISO9300 + public static double calculateViscosity(GasProps gas) { + double BOLTZMANN = 1.380649e-23; // 玻尔兹曼常数 (J/K) + double Avogadro = 6.02214076e23; // 添加阿伏伽德罗常数 + LJ mixLj = getLJParameters(gas); // Lennard-Jones参数 + // 转换为SI单位 + double sigma = mixLj.sigma ; // Å → m + double epsilon = mixLj.epsilonK * BOLTZMANN; // K → + + // 计算无量纲温度 + double T = gas.getdTf(); + double Tstar = T * BOLTZMANN / epsilon; + + double omega = getOmega(Tstar); + + // 查普曼-恩斯柯格公式 + double M_avg = gas.getdMrx() / 1000; // g/mol → kg/mol + double m = M_avg / Avogadro; + double viscosity = 5.0 / 16.0 * Math.sqrt(Math.PI * m * BOLTZMANN * T) + / (Math.PI * Math.pow(sigma, 2) * omega); + + gas.setdMu(viscosity); + return viscosity; + + } + + + + public static LJ getLJParameters(GasProps gas) { + double sigmaMix = 0; + double epsilonKMix = 0; + int n = gas.adMixture.length; + + for (int i = 0; i < n; i++) { + for (int j = 0; j < n; j++) { + + double xixj = gas.adMixture[i] * gas.adMixture[j]; + // Lorentz规则 + double sigmaPair = (LjParameters[i][0] + LjParameters[j][0]) / 2; + sigmaMix += xixj * sigmaPair; + // Berthelot规则 + 极性修正 + double epsilonPair = Math.sqrt(LjParameters[i][1] * LjParameters[j][1]) * polarityFactors[i] * polarityFactors[j]; + epsilonKMix += xixj * epsilonPair; + } + } + + return new LJ(sigmaMix, epsilonKMix); + } + + // 天然气21个组分的LJ参数 + + public static double[][] LjParameters = { + {3.758e-10, 148.6}, // 甲烷 C1 + {3.798e-10, 71.4}, // 氮气 N2 + {3.941e-10, 195.2}, // 二氧化碳 CO2 + {4.443e-10, 215.7}, // 乙烷 C2 + {5.118e-10, 237.1}, // 丙烷 C3 + {2.75e-10, 80.0}, // 水 H2O + {3.623e-10, 301.1}, // 硫化氢 H2S + {2.827e-10, 59.7}, // 氢气 H2 + {3.690e-10, 91.7}, // 一氧化碳 CO + {3.467e-10, 106.7}, // 氧气 O2 + {5.278e-10, 267.0}, // 异丁烷 iC4 + {5.341e-10, 274.7}, // 正丁烷 nC4 + {5.734e-10, 330.1}, // 异戊烷 iC5 + {5.784e-10, 341.1}, // 正戊烷 nC5 + {6.260e-10, 412.3}, // 己烷 C6 + {6.812e-10, 467.3}, // 庚烷 C7（注：此处庚烷参数按趋势估算） + {7.294e-10, 532.1}, // 辛烷 C8（注：此处辛烷参数按趋势估算） + {7.700e-10, 614.2}, // 壬烷 C9（注：此处壬烷参数按趋势估算） + {8.170e-10, 681.5}, // 癸烷 C10（注：此处癸烷参数按趋势估算） + {2.551e-10, 10.22}, // 氦气 He + {3.405e-10, 124.0} // 氩气 Ar + }; + // 极性修正系数（1.0表示非极性） + // 极性修正系数数组（与LjParameters数组顺序严格对应） + private static final double[] polarityFactors = { + // 1. 甲烷 C1 + 1.00, // 非极性分子（对称四面体结构） + + // 2. 氮气 N2 + 1.02, // 微弱四极矩（J. Chem. Phys. 129, 034306） + + // 3. 二氧化碳 CO2 + 1.05, // 四极矩修正（Ind. Eng. Chem. Res. 2019, 58, 5, 1964–1972） + + // 4. 乙烷 C2 + 1.00, // 非极性（对称结构） + + // 5. 丙烷 C3 + 1.00, // 非极性（链状烷烃） + + // 6. 水 H2O + 1.18, // 强极性修正（J. Phys. Chem. B 2005, 109, 15, 7053–7062） + + // 7. 硫化氢 H2S + 1.12, // 中等极性（J. Chem. Eng. Data 2008, 53, 3, 726–729） + + // 8. 氢气 H2 + 1.00, // 非极性（同核双原子） + + // 9. 一氧化碳 CO + 1.03, // 微弱偶极矩（J. Mol. Liq. 2020, 320, 114432） + + // 10. 氧气 O2 + 1.01, // 微弱顺磁性（通常视为非极性） + + // 11. 异丁烷 iC4 + 1.00, // 支链烷烃（非极性） + + // 12. 正丁烷 nC4 + 1.00, // 直链烷烃（非极性） + + // 13. 异戊烷 iC5 + 1.00, // 支链烷烃（非极性） + + // 14. 正戊烷 nC5 + 1.00, // 直链烷烃（非极性） + + // 15. 己烷 C6 + 1.00, // 长链烷烃（非极性） + + // 16. 庚烷 C7 + 1.00, // 长链烷烃（非极性） + + // 17. 辛烷 C8 + 1.00, // 长链烷烃（非极性） + + // 18. 壬烷 C9 + 1.00, // 长链烷烃（非极性） + + // 19. 癸烷 C10 + 1.00, // 长链烷烃（非极性） + + // 20. 氦气 He + 1.00, // 惰性气体（非极性） + + // 21. 氩气 Ar + 1.00 // 惰性气体（非极性） + }; + // 碰撞积分Ω(2,2)* 近似计算（Neufeld多项式） + private static double getOmega(double Tstar) { + if (Tstar < 0.1) { + return 2.0 / (3.0 * Tstar); // 低温量子修正 + } else if (Tstar > 400) { + return 0.92 * Math.log(Tstar) / Tstar; // 高温渐近解 + } else { + double A = 1.16145; + double B = 0.14874; + double C = 0.52487; + double D = 0.77320; + double E = 2.16178; + double F = 2.43787; + return A / Math.pow(Tstar, B) + + C / Math.exp(D * Tstar) + + E / Math.exp(F * Tstar); + } + } + public static class LJ { + public final double sigma; // 碰撞直径 (m) + public final double epsilonK; // ε/k (K) + + public LJ(double sigma, double epsilonK) { + this.sigma = sigma; + this.epsilonK = epsilonK; + } + } +}