代码修改

ningshuxia

2024-03-20 785c92f5078cb7aeb05b13f6da627defb6525b36

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
using IStation.Epanet.Enums;
using System;
 
namespace IStation.Epanet.Network.Structures
{
 
    public class Pipe : Link
    {
        public Pipe(string name) : base(name) { }
 
        public override LinkType LinkType => LinkType.Pipe;
 
        public bool HasCheckValve { get; set; }
 
        public override void InitResistance(HeadLossFormulaType formflag, double hexp)
        {
            FlowResistance = Constants.CSMALL;
            double d = Diameter;
 
            switch (formflag)
            {
                case HeadLossFormulaType.HW:
                    FlowResistance = 4.727 * Lenght / Math.Pow(Kc, hexp) / Math.Pow(d, 4.871);
                    break;
                case HeadLossFormulaType.DW:
                    FlowResistance = Lenght / 2.0 / 32.2 / d / Math.Pow(Math.PI * Math.Pow(d, 2) / 4.0, 2);
                    break;
                case HeadLossFormulaType.CM:
                    FlowResistance = Math.Pow(4.0 * Kc / (1.49 * Math.PI * d * d), 2)
                                     * Math.Pow(d / 4.0, -1.333) * Lenght;
                    break;
            }
        }
 
        public override void ConvertUnits(Network nw)
        {
            FieldsMap fMap = nw.FieldsMap;
 
            if (nw.FormFlag == HeadLossFormulaType.DW)
                Kc /= 1000.0 * fMap.GetUnits(FieldType.ELEV);
 
            Diameter /= fMap.GetUnits(FieldType.DIAM);
            Lenght /= fMap.GetUnits(FieldType.LENGTH);
 
            Km = 0.02517 * Km / Math.Pow(Diameter, 2) / Math.Pow(Diameter, 2);
 
            Kb /= Constants.SECperDAY;
            Kw /= Constants.SECperDAY;
        }
 
 
    }
}