ArcEngine下架空线之悬链线的模拟

最近在做电力方面的东西，需要对架空线进行三维建模。经过前端时间的搜索与实践，最终实现了AE下悬链线的模拟。首先说明几个问题：

（1）实地架设杆塔时，两个相邻杆塔的高程一般是不相同的，这就需要用到不等高悬点架空线的悬链线方程，在线路三维建模时要进行点位的判断

（2）在计算悬链线高度时要进行三维空间XYZ到二维空间ZU的转换，这里的u坐标轴是杆塔连线在XY平面的投影，zu所在三维空间Z=低悬点的高程值。

（3）采样精度设置为1。采样完成后将二维点再转换成三维点即可。

这里用AE写了一个类，用于悬链线的生成，公式主要参考《架空输电线路设计- 孟遂民》这本书，没有考虑复杂的温度、风力、覆冰等条件。代码如下：

/// <summary>
    /// 根据起点和终点生成悬链线
    /// 不等高悬点架空线方程
    /// 悬链线长度方程
    /// </summary>
    class PowerLineCreator
    {
        public IPoint FromPoint //起点
        { 
            get; 
            set;
        }
        public IPoint ToPoint   //终点
        { 
            get; 
            set; 
        }

        public double HorizontalStress  //弧垂最低点应力 (已知条件、水平应力)  σ0
        {
            get;
            set;
        }

        public int SampleAccuracy   //插值精度
        {
            get;
            set;
        }
        public double r    //比载：单位长度架空线上所受的荷载折算到单位载面积上的数值
        {
            get;
            set;
        }

        public double PlaneLength//悬链线档距
        {
            get
            {
                double planelength_ = Math.Sqrt(Math.Pow(FromPoint.X - ToPoint.X,2) + Math.Pow(FromPoint.Y - ToPoint.Y,2));
                return planelength_;
            }
            
        }

        public double L  //不等高悬点架空线长度
        {
            get
            {
                return Math.Sqrt(Math.Pow(L_h0,2) + Math.Pow(H,2));
            }
        }
        public double H //   悬点高程差
        {
            get 
            {
                return Math.Abs(FromPoint.Z - ToPoint.Z);    
            }      
        }

        public double L_h0  //等高悬点架空线的档内悬链线长度
        {
            get
            {
                double L = 2 * HorizontalStress / r * Math.Sinh(r * PlaneLength / (2 * HorizontalStress));
                return L;
            }
        }

        public double a   //低悬点至弧垂最低点的水平距离 
        {
            get
            {
                return PlaneLength / 2 - HorizontalStress / r * SpecialFunction.asinh(H / L_h0);
            }
        }

        public double b  //弧垂最低点至高悬点的水平距离
        {
            get
            {
                return L_h0 - a;
            }
        }

        public PowerLineCreator ()
        {
            
        }
        /// <summary>
        /// 悬链线方程
        /// </summary>
        /// <param name="u">二维空间下u坐标</param>
        /// <returns>悬链线高度</returns>
        public double Get_Z (double u) 
        {
            double Z = H / L_h0 * (2 * HorizontalStress / r * Math.Sinh(r * u / 2 / HorizontalStress) * Math.Cosh(r * (L_h0 - u) / 2 / HorizontalStress)) -
                Math.Sqrt(1 + Math.Pow(H / L_h0,2)) * 2 * HorizontalStress / r * Math.Sinh(r * u / 2 / HorizontalStress) * Math.Sinh(r * (L_h0 - u) / 2 / HorizontalStress);
            return Z;
        }
        /// <summary>
        /// 根据起始点生成悬链线插值点
        /// </summary>
        /// <returns></returns>
        public IPointCollection CreateLinePoints ()
        {
            IPointCollection ptCollection = new PolylineClass();
            ILine line_ = new LineClass();
            line_.PutCoords(FromPoint,ToPoint);
            double pAngle = line_.Angle;
            //采样悬链线点
            for(int i = 0;i < PlaneLength-1;i+=SampleAccuracy)
            {
                //二、三维空间坐标的转换
                IPoint point_ = new PointClass();
                (point_ as IZAware).ZAware = true;
                point_.X = FromPoint.X + i * Math.Cos(pAngle);
                point_.Y = FromPoint.Y + i * Math.Sin(pAngle);
                point_.Z = FromPoint.Z + Get_Z(i);
                ptCollection.AddPoint(point_);

            }
            ptCollection.AddPoint(ToPoint);
            
            return ptCollection;
        }

    }

最终实现效果图：

转载请注明出处：http://blog.csdn.net/liushuo_whu/article/details/21323807谢谢！