diff --git a/excel_addin/Tension.cs b/excel_addin/Tension.cs index 4f7d43f..692489a 100644 --- a/excel_addin/Tension.cs +++ b/excel_addin/Tension.cs @@ -13,11 +13,11 @@ namespace excel_addin //全局用临界档距，是个单例 private static List gCriticalSpan; - [ExcelFunction(Description = "覆冰增重荷载")]//仅包含覆冰重量 public static double ice_load(double d, double iceThickness)//d导线直径，单位mm,iceThickness 覆冰厚度 单位mm { - double val = 9.80665 * 0.9 * Math.PI * iceThickness * (iceThickness + d) / 1000; + double val; + val = 9.80665 * 0.9 * Math.PI * iceThickness * (iceThickness + d) / 1000; return val; } @@ -25,33 +25,33 @@ namespace excel_addin public static double alpha_load(double basicV)//basicV 基本风速 单位m/s { double _a = 0; - if (basicV < 20) - { - _a = 1; - - } - else - { - if (basicV >= 20 && basicV < 27) + if (basicV < 20) { - _a = 0.85; + _a = 1; + } else { - if (basicV >= 27 && basicV < 31.5) + if (basicV >= 20 && basicV < 27) { - _a = 0.75; + _a = 0.85; } else { - if (basicV >= 31.5) + if (basicV >= 27 && basicV < 31.5) { - _a = 0.7; + _a = 0.75; + } + else + { + if (basicV >= 31.5) + { + _a = 0.7; + } } } - } - } + } return _a; } @@ -59,71 +59,77 @@ namespace excel_addin public static double alpha_swing(double basicV)//basicV 基本风速 单位m/s { double _a = 0; - if (basicV < 20) - { - _a = 1; - - } - else - { - if (basicV >= 20 && basicV < 27) + if (basicV < 20) { - _a = 0.75; + _a = 1; + } else { - if (basicV >= 27 && basicV < 31.5) + if (basicV >= 20 && basicV < 27) { - _a = 0.61; + _a = 0.75; } else { - if (basicV >= 31.5) + if (basicV >= 27 && basicV < 31.5) { _a = 0.61; } + else + { + if (basicV >= 31.5) + { + _a = 0.61; + } + } } - } - } + } return _a; } [ExcelFunction(Description = "覆冰后风荷载增大系数")] public static double iceB(int iceThickness)//iceThickness 覆冰厚度 单位mm { - var code = new Dictionary(); - code.Add(0, 1.0); - code.Add(5, 1.1); - code.Add(10, 1.2); - code.Add(15, 1.3); - if (iceThickness >= 20) - { - throw new Exception("缺少20mm冰以上覆冰后荷载增大系数!"); - } - return code[iceThickness]; + + var code = new Dictionary(); + code.Add(0, 1.0); + code.Add(5, 1.1); + code.Add(10, 1.2); + code.Add(15, 1.3); + if (iceThickness >= 20) + { + //throw new Exception("缺少20mm冰以上覆冰后荷载增大系数!"); + return -9999999; + } + return code[iceThickness]; + } [ExcelFunction(Description = "电线受风形体系数")] public static double mu_sc(double d,int iceThickness)//d导线直径，单位mm,iceThickness 覆冰厚度 单位mm { - if (iceThickness > 0) - { - return 1.2; - } - if (d < 17) - { - return 1.2; - } - //if d >=17 - return 1.1; + + if (iceThickness > 0) + { + return 1.2; + } + if (d < 17) + { + return 1.2; + } + //if d >=17 + return 1.1; + } [ExcelFunction(Description = "电线单位风荷载")] public static double wind_loadPN(double ave_h_v, double d, double alpha_load,double mu_sc, int iceThickness)//ave_h_v平均高处的风速 , d 导线直径 单位mm,alpha_load 风压不均匀系数，mu_sc电线体型系数,iceThickness 覆冰厚度 单位mm { - double val = 0.625 * Math.Pow(ave_h_v, 2) * (d + 2 * iceThickness) * alpha_load*mu_sc / 1000; + double val; + val = 0.625 * Math.Pow(ave_h_v, 2) * (d + 2 * iceThickness) * alpha_load * mu_sc / 1000; return val; } @@ -141,120 +147,122 @@ namespace excel_addin //临界档距 public static double[] criticalSpanForExcel(int maxSpan, double Epsilon, double Alpha, double area, double lowTemp_t, double lowTemp_load, double lowTemp_tens, double avrTemp_t, double avrTemp_load, double avrTemp_tens, double strongWind_t, double strongWind_load, double strongWind_tens, double ice_t, double ice_load, double ice_tens) { - List criticalSpanList = criticalSpan(maxSpan, Epsilon, Alpha, area, lowTemp_t, lowTemp_load, lowTemp_tens, avrTemp_t, avrTemp_load, avrTemp_tens, strongWind_t, strongWind_load, strongWind_tens, ice_t, ice_load, ice_tens); - double[] value = new double[criticalSpanList.Count * 2]; - for (int foo = 0; foo < criticalSpanList.Count; foo++) - { - value[2 * foo] = criticalSpanList[foo][0]; - value[2 * foo + 1] = criticalSpanList[foo][1]; - } - return value; + List criticalSpanList = criticalSpan(maxSpan, Epsilon, Alpha, area, lowTemp_t, lowTemp_load, lowTemp_tens, avrTemp_t, avrTemp_load, avrTemp_tens, strongWind_t, strongWind_load, strongWind_tens, ice_t, ice_load, ice_tens); + double[] value = new double[criticalSpanList.Count * 2]; + for (int foo = 0; foo < criticalSpanList.Count; foo++) + { + value[2 * foo] = criticalSpanList[foo][0]; + value[2 * foo + 1] = criticalSpanList[foo][1]; + } + return value; } //临界档距 public static List criticalSpan(int maxSpan, double Epsilon, double Alpha, double area, double lowTemp_t, double lowTemp_load, double lowTemp_tens, double avrTemp_t, double avrTemp_load, double avrTemp_tens, double strongWind_t, double strongWind_load, double strongWind_tens, double ice_t, double ice_load, double ice_tens) { - List criticalSpanList = new List(); - //计算4个工况下的Fx - //t_dangju 试探最大Fx时使用的档距 - int[] MaxFxConditionFlag = new int[maxSpan + 1];//记录不同档距下，哪种工况的Fx最大。1代表最低温，2代表年平均温，3代表大风，4代表覆冰。0位置元素无效。 - for (int span = 1; span <= maxSpan; span++) - { - double Fx_lowTemp; - double Fx_avrTemp; - double Fx_strongWind; - double Fx_Ice; - Fx_lowTemp = Fx(Epsilon, lowTemp_load, lowTemp_tens, Alpha, area, lowTemp_t, span); - Fx_avrTemp = Fx(Epsilon, avrTemp_load, avrTemp_tens, Alpha, area, avrTemp_t, span); - Fx_strongWind = Fx(Epsilon, strongWind_load, strongWind_tens, Alpha, area, strongWind_t, span); - Fx_Ice = Fx(Epsilon, ice_load, ice_tens, Alpha, area, ice_t, span); - double maxFx; - maxFx = Fx_lowTemp; - MaxFxConditionFlag[span] = 1; - if (maxFx < Fx_avrTemp) + List criticalSpanList = new List(); + //计算4个工况下的Fx + //t_dangju 试探最大Fx时使用的档距 + int[] MaxFxConditionFlag = new int[maxSpan + 1];//记录不同档距下，哪种工况的Fx最大。1代表最低温，2代表年平均温，3代表大风，4代表覆冰。0位置元素无效。 + for (int span = 1; span <= maxSpan; span++) { - maxFx = Fx_avrTemp; - MaxFxConditionFlag[span] = 2; - } - if (maxFx < Fx_strongWind) - { - maxFx = Fx_strongWind; - MaxFxConditionFlag[span] = 3; - } - if (maxFx < Fx_Ice) - { - maxFx = Fx_Ice; - MaxFxConditionFlag[span] = 4; - } - MaxFxConditionFlag[0] = MaxFxConditionFlag[1];//为了后面判断临界档距突变，不会从1米档距开始跳变。 - double currentTension = 0; - double currentLoad = 0; - double currentTemperature = 0; - switch (MaxFxConditionFlag[span]) - { - case 1: - currentLoad = lowTemp_load; - currentTension = lowTemp_tens; - currentTemperature = lowTemp_t; - break; - case 2: - currentLoad = avrTemp_load; - currentTension = avrTemp_tens; - currentTemperature = avrTemp_t; - break; - case 3: - currentLoad = strongWind_load; - currentTension = strongWind_tens; - currentTemperature = strongWind_t; - break; - case 4: - currentLoad = ice_load; - currentTension = ice_tens; - currentTemperature = ice_t; - break; - default: - break; - } - if (MaxFxConditionFlag[span] != MaxFxConditionFlag[span - 1])//有突变，找到了临界档距的范围 - { - double formerTension = 0; - double formerLoad = 0; - double formerTemperature = 0; - switch (MaxFxConditionFlag[span - 1]) + + double Fx_lowTemp; + double Fx_avrTemp; + double Fx_strongWind; + double Fx_Ice; + Fx_lowTemp = Fx(Epsilon, lowTemp_load, lowTemp_tens, Alpha, area, lowTemp_t, span); + Fx_avrTemp = Fx(Epsilon, avrTemp_load, avrTemp_tens, Alpha, area, avrTemp_t, span); + Fx_strongWind = Fx(Epsilon, strongWind_load, strongWind_tens, Alpha, area, strongWind_t, span); + Fx_Ice = Fx(Epsilon, ice_load, ice_tens, Alpha, area, ice_t, span); + double maxFx; + maxFx = Fx_lowTemp; + MaxFxConditionFlag[span] = 1; + if (maxFx < Fx_avrTemp) + { + maxFx = Fx_avrTemp; + MaxFxConditionFlag[span] = 2; + } + if (maxFx < Fx_strongWind) + { + maxFx = Fx_strongWind; + MaxFxConditionFlag[span] = 3; + } + if (maxFx < Fx_Ice) + { + maxFx = Fx_Ice; + MaxFxConditionFlag[span] = 4; + } + MaxFxConditionFlag[0] = MaxFxConditionFlag[1];//为了后面判断临界档距突变，不会从1米档距开始跳变。 + double currentTension = 0; + double currentLoad = 0; + double currentTemperature = 0; + switch (MaxFxConditionFlag[span]) { case 1: - formerLoad = lowTemp_load; - formerTension = lowTemp_tens; - formerTemperature = lowTemp_t; + currentLoad = lowTemp_load; + currentTension = lowTemp_tens; + currentTemperature = lowTemp_t; break; case 2: - formerLoad = avrTemp_load; - formerTension = avrTemp_tens; - formerTemperature = avrTemp_t; + currentLoad = avrTemp_load; + currentTension = avrTemp_tens; + currentTemperature = avrTemp_t; break; case 3: - formerLoad = strongWind_load; - formerTension = strongWind_tens; - formerTemperature = strongWind_t; + currentLoad = strongWind_load; + currentTension = strongWind_tens; + currentTemperature = strongWind_t; break; case 4: - formerLoad = ice_load; - formerTension = ice_tens; - formerTemperature = ice_t; + currentLoad = ice_load; + currentTension = ice_tens; + currentTemperature = ice_t; break; default: break; } - double t1 = 24 / Epsilon * (formerTension - currentTension) / area; - double t2 = 24 * Alpha * (formerTemperature - currentTemperature); - double t3 = Math.Pow(formerLoad / formerTension, 2) - Math.Pow(currentLoad / currentTension, 2); - double _criticalSpan = Math.Sqrt((t1 + t2) / t3); - //List中存有数组，0位置为临界档距，1位置为控制条件 - criticalSpanList.Add(new double[] { _criticalSpan, MaxFxConditionFlag[span - 1] }); + if (MaxFxConditionFlag[span] != MaxFxConditionFlag[span - 1])//有突变，找到了临界档距的范围 + { + double formerTension = 0; + double formerLoad = 0; + double formerTemperature = 0; + switch (MaxFxConditionFlag[span - 1]) + { + case 1: + formerLoad = lowTemp_load; + formerTension = lowTemp_tens; + formerTemperature = lowTemp_t; + break; + case 2: + formerLoad = avrTemp_load; + formerTension = avrTemp_tens; + formerTemperature = avrTemp_t; + break; + case 3: + formerLoad = strongWind_load; + formerTension = strongWind_tens; + formerTemperature = strongWind_t; + break; + case 4: + formerLoad = ice_load; + formerTension = ice_tens; + formerTemperature = ice_t; + break; + default: + break; + } + double t1 = 24 / Epsilon * (formerTension - currentTension) / area; + double t2 = 24 * Alpha * (formerTemperature - currentTemperature); + double t3 = Math.Pow(formerLoad / formerTension, 2) - Math.Pow(currentLoad / currentTension, 2); + double _criticalSpan = Math.Sqrt((t1 + t2) / t3); + //List中存有数组，0位置为临界档距，1位置为控制条件 + criticalSpanList.Add(new double[] { _criticalSpan, MaxFxConditionFlag[span - 1] }); + } } - } - criticalSpanList.Add(new double[] { maxSpan, MaxFxConditionFlag[maxSpan] }); - gCriticalSpan = new List(criticalSpanList); - return criticalSpanList; + criticalSpanList.Add(new double[] { maxSpan, MaxFxConditionFlag[maxSpan] }); + gCriticalSpan = new List(criticalSpanList); + return criticalSpanList; + } //从状态方程接触张力 @@ -271,7 +279,7 @@ namespace excel_addin public static double tensionFrmoStateEquation(double start_tension, double E, double alpha, double area, double control_load, double control_tension, double control_t, double load, double t, double l) { double A, B; - //A = E / 24 * Math.Pow(control_load * l / control_tension, 2) - control_tension / area + alpha * E * (t - control_t); + A = E / 24 * Math.Pow(control_load * l / control_tension, 2) - control_tension / area + alpha * E * (t - control_t); A = Fx(E, control_load, control_tension, alpha, area, control_t, l) + alpha * E * t; B = E * Math.Pow(load / area * l, 2) / 24; double tension;//牛顿法初值 @@ -281,185 +289,54 @@ namespace excel_addin { return tension; } - return -999999; + return -999999999999; + } + //拓展代表档距 + public static double[] expandSpan(int startSpan, int endSpan, int stepSpan,double [] criticalSpan) + { + List spans = new List(); + spans.Add(startSpan); + for (int foo = 1; foo < Convert.ToInt32(1.0 * (endSpan - startSpan) / stepSpan); foo++) + { + spans.Add(startSpan + stepSpan * foo); + } + spans.Add(endSpan); + for (int bar = 0; bar < criticalSpan.Length; bar++) + { + spans.Add(criticalSpan[bar]); + } + spans.Sort(); + return spans.ToArray(); + + } - //public static bool CalZhangLi(calStressLoadNew inPut, ref double ZhangLi) - //{ - // //张力=T6 - // //计算最大使用张力 - // double ZuiDaShiYongZhangLi; //最大使用张力=计算拉断力*保证破断张力取计算拉断力的/安全系数 - // ZuiDaShiYongZhangLi = inPut.EDingLaDuanLi * inPut.BaoZhengPoDuanZhangLiXiShu / 100 / inPut.AnQuan; - // //ZuiDaShiYongZhangLi = inPut.LaDuanLi / inPut.AnQuan; - // //计算温度线膨胀系数 - // double XianPengZhang = inPut.XianPengZhang; - // //计算平均运行张力 - // double PingJunYunXingZhangLi = inPut.EDingLaDuanLi * inPut.BaoZhengPoDuanZhangLiXiShu / 100 * inPut.NianPingJunXishu / 100;//年平均单位是% - // //double a, b, T1, T2, T3, T4, T5, T6; - // //计算待求工况比载 - // double DaiHeZai = 0; //待求工况比载 - // if (CalHeZai(inPut, ref DaiHeZai) == false) - // { - // return false; - // } - // //临界档距表计算 - // //一、初始计算数组,5项分别代表:1:气温，2:荷载，3:容许张力，4:P/T，5:Fo=-(T+A*E*a*t)(N) - // double[] ChuShiZuiDiQiWen = new double[5]; //最低温 - // double[] ChuShiNianPingJun = new double[5]; //年平均气温 - // double[] ChuShiDaFeng = new double[5]; //大风 - // double[] ChuShiFuBing = new double[5]; //最大覆冰 - // //1:存储各工况温度 - // ChuShiZuiDiQiWen[0] = inPut.DiWen; - // ChuShiNianPingJun[0] = inPut.PingWen; - // ChuShiDaFeng[0] = inPut.FengWen; - // ChuShiFuBing[0] = inPut.BingWen; - // //定义各个工况结构体，将每个结构体中的待求工况设置为对应工况数据，因为计算荷载时前面写的计算荷载函数中调取的是待求工况的温度等信息，所以需要对应各个工况 - // calStressLoadNew ZuiDiWen = inPut; - // calStressLoadNew NianPingJun = inPut; - // calStressLoadNew DaFeng = inPut; - // calStressLoadNew FuBing = inPut; - // ZuiDiWen.DaiBing = ZuiDiWen.DiBing; - // ZuiDiWen.DaiFeng = ZuiDiWen.DiFeng; - // ZuiDiWen.DaiZheSuanFengSu = false; - // ZuiDiWen.DaiWen = ZuiDiWen.DiWen; - // NianPingJun.DaiBing = NianPingJun.PingBing; - // NianPingJun.DaiFeng = NianPingJun.PingFeng; - // NianPingJun.DaiZheSuanFengSu = false; - // NianPingJun.DaiWen = NianPingJun.PingWen; - // DaFeng.DaiBing = DaFeng.FengBing; - // //double wind_velocity; - // //wind_velocity = Math.Round(DaFeng.FengFeng * Math.Pow(DaFeng.PingJunGaoDu / DaFeng.SheJiFengSuGaoDu, 0.16), 2); - // DaFeng.DaiFeng = DaFeng.FengFeng; - // DaFeng.DaiZheSuanFengSu = true; - // DaFeng.DaiWen = DaFeng.FengWen; - // FuBing.DaiBing = FuBing.BingBing; - // FuBing.DaiFeng = FuBing.BingFeng; - // FuBing.DaiZheSuanFengSu = false; - // FuBing.DaiWen = FuBing.BingWen; - // //2:计算并存储各个工况荷载 - // CalHeZai(ZuiDiWen, ref ChuShiZuiDiQiWen[1]); - // CalHeZai(NianPingJun, ref ChuShiNianPingJun[1]); - // CalHeZai(DaFeng, ref ChuShiDaFeng[1]); - // CalHeZai(FuBing, ref ChuShiFuBing[1]); - // //3:存储各工况容许张力 - // ChuShiZuiDiQiWen[2] = ZuiDaShiYongZhangLi; - // ChuShiNianPingJun[2] = PingJunYunXingZhangLi; - // ChuShiDaFeng[2] = ZuiDaShiYongZhangLi; - // ChuShiFuBing[2] = ZuiDaShiYongZhangLi; - // //4:计算并存储P/T 荷载/张力 - // ChuShiZuiDiQiWen[3] = ChuShiZuiDiQiWen[1] / ChuShiZuiDiQiWen[2]; - // ChuShiNianPingJun[3] = ChuShiNianPingJun[1] / ChuShiNianPingJun[2]; - // ChuShiDaFeng[3] = ChuShiDaFeng[1] / ChuShiDaFeng[2]; - // ChuShiFuBing[3] = ChuShiFuBing[1] / ChuShiFuBing[2]; + //架线表k值 + public static double tensionK(double load, double tension) //load 综合荷载N ,tension 张力 N + { + double value; + value = load / tension / 2; + return value; - // //计算4个工况下的Fx by 杜孟远2020.4.17 - // //t_dangju 试探最大Fx时使用的档距 - // int[] MaxFxConditionFlag = new int[inPut.DangJu + 1];//记录不同t_dangju下，哪种工况的Fx最大。1代表最低温，2代表年平均温，3代表大风，4代表覆冰。0位置元素无效。 - // MaxFxConditionFlag[0] = 1; - // for (int t_dangju = inPut.DangJu; t_dangju <= inPut.DangJu; t_dangju++) - // { - // double Fx_dangjuZuiDiWen; - // double Fx_dangjuNianPingJun; - // double Fx_dangjuDaFeng; - // double Fx_dangjuFuBing; - // Fx_dangjuZuiDiWen = Fx(inPut.MoLiang, ChuShiZuiDiQiWen[1], ChuShiZuiDiQiWen[2], inPut.XianPengZhang, inPut.JieMianJi, ChuShiZuiDiQiWen[0], t_dangju); - // Fx_dangjuNianPingJun = Fx(inPut.MoLiang, ChuShiNianPingJun[1], ChuShiNianPingJun[2], inPut.XianPengZhang, inPut.JieMianJi, ChuShiNianPingJun[0], t_dangju); - // Fx_dangjuDaFeng = Fx(inPut.MoLiang, ChuShiDaFeng[1], ChuShiDaFeng[2], inPut.XianPengZhang, inPut.JieMianJi, ChuShiDaFeng[0], t_dangju); - // Fx_dangjuFuBing = Fx(inPut.MoLiang, ChuShiFuBing[1], ChuShiFuBing[2], inPut.XianPengZhang, inPut.JieMianJi, ChuShiFuBing[0], t_dangju); - // double maxFx; - // maxFx = Fx_dangjuZuiDiWen; - // MaxFxConditionFlag[t_dangju] = 1; - // if (maxFx < Fx_dangjuNianPingJun) - // { - // maxFx = Fx_dangjuNianPingJun; - // MaxFxConditionFlag[t_dangju] = 2; - // } - // if (maxFx < Fx_dangjuDaFeng) - // { - // maxFx = Fx_dangjuDaFeng; - // MaxFxConditionFlag[t_dangju] = 3; - // } - // if (maxFx < Fx_dangjuFuBing) - // { - // maxFx = Fx_dangjuFuBing; - // MaxFxConditionFlag[t_dangju] = 4; - // } - // double currentTension = 0; - // double currentLoad = 0; - // double currentTemperature = 0; - // switch (MaxFxConditionFlag[t_dangju]) - // { - // case 1: - // currentLoad = ChuShiZuiDiQiWen[1]; - // currentTension = ChuShiZuiDiQiWen[2]; - // currentTemperature = ChuShiZuiDiQiWen[0]; - // break; - // case 2: - // currentLoad = ChuShiNianPingJun[1]; - // currentTension = ChuShiNianPingJun[2]; - // currentTemperature = ChuShiNianPingJun[0]; - // break; - // case 3: - // currentLoad = ChuShiDaFeng[1]; - // currentTension = ChuShiDaFeng[2]; - // currentTemperature = ChuShiDaFeng[0]; - // break; - // case 4: - // currentLoad = ChuShiFuBing[1]; - // currentTension = ChuShiFuBing[2]; - // currentTemperature = ChuShiFuBing[0]; - // break; - // default: - // break; - // } - // //TODO: 目前不进入这个判断 by 杜孟远 2020.4.18 - // MaxFxConditionFlag[t_dangju] = MaxFxConditionFlag[t_dangju - 1];//TODO: 以后删掉这句 by 杜孟远 2020.4.18 - // if (MaxFxConditionFlag[t_dangju] != MaxFxConditionFlag[t_dangju - 1])//有突变，找到了临界档距的范围 - // { - // double criticalSpan; - // double formerTension = 0; - // double formerLoad = 0; - // double formerTemperature = 0; - // switch (MaxFxConditionFlag[t_dangju - 1]) - // { - // case 1: - // formerLoad = ChuShiZuiDiQiWen[1]; - // formerTension = ChuShiZuiDiQiWen[2]; - // formerTemperature = ChuShiZuiDiQiWen[0]; - // break; - // case 2: - // formerLoad = ChuShiNianPingJun[1]; - // formerTension = ChuShiNianPingJun[2]; - // formerTemperature = ChuShiNianPingJun[0]; - // break; - // case 3: - // formerLoad = ChuShiDaFeng[1]; - // formerTension = ChuShiDaFeng[2]; - // formerTemperature = ChuShiDaFeng[0]; - // break; - // case 4: - // formerLoad = ChuShiFuBing[1]; - // formerTension = ChuShiFuBing[2]; - // formerTemperature = ChuShiFuBing[0]; - // break; - // default: - // break; - // } - // double t1 = 24 / inPut.MoLiang * (formerTension - currentTension) / inPut.JieMianJi; - // double t2 = 24 * inPut.XianPengZhang * (formerTemperature - currentTemperature); - // double t3 = Math.Pow(formerLoad / formerTension, 2) - Math.Pow(currentLoad / currentTension, 2); - // criticalSpan = Math.Sqrt((t1 + t2) / t3);//TODO:暂时不用 - // //Console.WriteLine("{0}", criticalSpan); + } - // } - // double start_tension = inPut.EDingLaDuanLi / inPut.AnQuan; - // ZhangLi = tensionFrmoStateEquation(start_tension, inPut.MoLiang, inPut.XianPengZhang, inPut.JieMianJi, currentLoad, currentTension, currentTemperature, DaiHeZai, inPut.DaiWen, inPut.DangJu); - // } + //最大弧垂 + public static double line_arc(double load, double tension, double span) + { + double value; + value = load * span * span / 8 / tension; + return value; + } - // ZhangLi = BaoLiuLiangWei(ZhangLi); - // return true; - //} + //百米架线弧垂 + public static double Meter100LineArc(double control_load, double control_tension,double control_temp ,double lower_temp,double E, double alpha, double area,double span) //control_load 综合荷载N ,control_tension 张力 N,control_temp 控制条件温度, lower_temp 需要的降温后温度,E 弹性模量,alpha 膨胀系数, area 截面积, span 档距 + { + double lower_tension = tensionFrmoStateEquation(control_tension, E, alpha, area, control_load, control_tension, control_temp, control_load, lower_temp, span);//降温后张力 + return line_arc(control_load, lower_tension, 100); + + + } } diff --git a/excel_addin/UDF.bas b/excel_addin/UDF.bas new file mode 100644 index 0000000..7d14b50 --- /dev/null +++ b/excel_addin/UDF.bas @@ -0,0 +1,23 @@ +Attribute VB_Name = "UDF" +Option Explicit +Function findControlCondition(span As Double, findStartCell As Range, count As Integer, column As Integer) +Dim foo As Integer +Dim criticalSpanN As Double +Dim criticalSpanN_1 As Double +Dim value As Double +For foo = 1 To count + criticalSpanN_1 = Cells(findStartCell.Row - 1 + foo - 1, findStartCell.column).value + criticalSpanN = Cells(findStartCell.Row - 1 + foo, findStartCell.column).value + If span <= criticalSpanN And span > criticalSpanN_1 Then + value = Cells(findStartCell.Row - 1 + foo, findStartCell.column - 1 + column).value + foo = count 'ѭ + + End If +Next foo +findControlCondition = value + +End Function + + + + diff --git a/excel_addin/力学曲线.cls b/excel_addin/力学曲线.cls new file mode 100644 index 0000000..9b1537d --- /dev/null +++ b/excel_addin/力学曲线.cls @@ -0,0 +1,119 @@ +VERSION 1.0 CLASS +BEGIN + MultiUse = -1 'True +END +Attribute VB_Name = "ѧ" +Attribute VB_GlobalNameSpace = False +Attribute VB_Creatable = False +Attribute VB_PredeclaredId = True +Attribute VB_Exposed = True +Option Explicit +Private Sub Worksheet_Change(ByVal Target As Range) + Application.EnableEvents = False + + If Target.Address = "$K$11" Or Target.Address = "$H$17" Or Target.Address = "$K$17" Or Target.Address = "$O$17" Then + + Dim criticalSpan As Variant + Dim length As Integer + Dim maxSpan As Integer + Dim Epsilon As Double + Dim Alpha As Double + Dim area As Double + Dim lowTemp_t As Double + Dim lowTemp_load As Double + Dim lowTemp_tens As Double + Dim avrTemp_t As Double + Dim avrTemp_load As Double + Dim avrTemp_tens As Double + Dim strongWind_t As Double + Dim strongWind_load As Double + Dim strongWind_tens As Double + Dim ice_t As Double + Dim ice_load As Double + Dim ice_tens As Double + maxSpan = 1000 + Epsilon = Range("J5") + Alpha = Range("J6") + area = Range("J9") + lowTemp_t = Range("I20") + lowTemp_load = Range("AD20") + lowTemp_tens = Range("AE20") + avrTemp_t = Range("I21") + avrTemp_load = Range("AD21") + avrTemp_tens = Range("AE21") + strongWind_t = Range("I22") + strongWind_load = Range("AD22") + strongWind_tens = Range("AE22") + ice_t = Range("I23") + ice_load = Range("AD23") + ice_tens = Range("AE23") + criticalSpan = Application.Run("criticalSpanForExcel", maxSpan, Epsilon, Alpha, area, lowTemp_t, lowTemp_load, lowTemp_tens, avrTemp_t, avrTemp_load, avrTemp_tens, strongWind_t, strongWind_load, strongWind_tens, ice_t, ice_load, ice_tens) + length = UBound(criticalSpan) + Dim startCell As Range + Set startCell = Range("W2") + Dim foo As Integer + Range("W2:AH2") = "" + For foo = 1 To length + Cells(startCell.Row, startCell.column - 1 + foo) = criticalSpan(foo) + Next foo + '޸Ĵ + expandSpan + + End If + + + Application.EnableEvents = True +End Sub + +Private Sub expandSpan() +Dim startSpan As Integer +Dim endSpan As Integer +Dim stepSpan As Integer +startSpan = Range("H17") +endSpan = Range("K17") +stepSpan = Range("O17") +Dim criticalSpan() As Double +Dim foo As Integer +Dim value As Double +Dim criticalSpanN As Integer +criticalSpanN = 0 +For foo = 1 To 10 + If Cells(Range("W2").Row, Range("W2").column + 2 * (foo - 1)) <> "" Then + criticalSpanN = criticalSpanN + 1 + ReDim Preserve criticalSpan(1 To criticalSpanN) + criticalSpan(foo) = Cells(Range("W2").Row, Range("W2").column + 2 * (foo - 1)) + End If +Next foo +Dim expandedSpan As Variant +expandedSpan = Application.Run("expandSpan", startSpan, endSpan, stepSpan, criticalSpan) +For foo = 1 To UBound(expandedSpan) + Cells(Range("B69").Row - 1 + foo, Range("B69").column) = expandedSpan(foo) +Next foo +' +Range("X69:Z69").Select +Selection.AutoFill Destination:=Range("X69" & ":Z" & 69 + UBound(expandedSpan) - 1), Type:=xlFillDefault +'޸ +Dim rangeA As String +Dim rangeB As String +For foo = 1 To 9 + rangeA = Chr((Asc("D") + 2 * (foo - 1))) + rangeB = Chr((Asc("D") + 2 * foo - 1)) + Range(rangeA & "69:" & rangeB & "69").Select + Selection.AutoFill Destination:=Range(rangeA & "69:" & rangeB & (69 + UBound(expandedSpan) - 1)), Type:=xlFillDefault +Next foo +'kֵ +For foo = 1 To 9 + rangeA = Chr((Asc("D") + 2 * (foo - 1))) + rangeB = Chr((Asc("D") + 2 * foo - 1)) + Range(rangeA & "161:" & rangeB & "161").Select + Selection.AutoFill Destination:=Range(rangeA & "161:" & rangeB & (161 + UBound(expandedSpan) - 1)), Type:=xlFillDefault +Next foo +'kֵ +For foo = 1 To 9 + rangeA = Chr((Asc("D") + 2 * (foo - 1))) + rangeB = Chr((Asc("D") + 2 * foo - 1)) + Range(rangeA & "209:" & rangeB & "209").Select + Selection.AutoFill Destination:=Range(rangeA & "209:" & rangeB & (209 + UBound(expandedSpan) - 1)), Type:=xlFillDefault +Next foo + +End Sub