Processing math: 100%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus
  • Indexed by Core Journals of China, Chinese S&T Journal Citation Reports
  • Chinese S&T Journal Citation Reports
  • Chinese Science Citation Database
QIN Jian, ZHANG Feikai, LI Qiying, LIU Chen. Automatic Search Method for Trestle Position of Freight Cableways[J]. Journal of Southwest Jiaotong University, 2022, 57(5): 1096-1102. doi: 10.3969/j.issn.0258-2724.20210106
Citation: DING Min, WANG Jiajia, JIANG Xiugen, CAO Qiongqiong, WANG Hongzhi. Out-of-Plane Stability Analysis Method for Circular Arch Structures[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 37-46. doi: 10.3969/j.issn.0258-2724.20191158

Out-of-Plane Stability Analysis Method for Circular Arch Structures

doi: 10.3969/j.issn.0258-2724.20191158
  • Received Date: 04 Dec 2019
  • Rev Recd Date: 20 Apr 2020
  • Available Online: 14 Dec 2020
  • Publish Date: 01 Feb 2021
  • In order to analyze the out-of-plane stability behavior of circular arch structures, according to the deformation correlation of circular arch with circular curved beam, the equilibrium equations of circular curved beam were firstly established by taking secondary moment effect into account. Combined with the geometric equation and physical equation of circular curved beam, the deflection control equation and torsion angle control equation of circular curved beam with free torsion were derived by considering large displacement. Both the general formats and corresponding simplified formats of analytical solution for circular curved beam deflection and torsion angle were obtained. Meanwhile, the expressions for deformation and internal force of circular curved beam were also derived. On this basis, the methods to analyze the out-of-plane bifurcation instability and extreme point instability of circular arch structure were presented. Critical load coefficients and their instability modes of four kinds of circular arch structure were calculated when out-of-plane bifurcation instability occurs, and the calculation results between the proposed model and the models from the literature were discussed; their load-displacement curves were calculated, and the extreme point instability of circular arch structures was analyzed. The results show that the critical load coefficient of the circular arch with two ends simply supported can be calculated by this model and has no difference with other models. In addition, this model is also useful to calculate the critical load coefficient of that with single hinge at mid-span or that with two ends inserted supported which are rarely seen in research. The out-of-plane bifurcation instability modes of all kinds of circular arch under in-plane uniformly distributed radial load are in the form of single symmetric wave. The radial load does not change the linear character of the out-of-plane load vs. displacement curve, but reduces the out-of-plane flexural rigidity. When the radial load reaches a certain value, the out-of-plane flexural rigidity becomes 0, and then the out-of-plane instability occurs.

     

  • 随着新基建进程的逐步推进和人们环保意识的逐渐提升,特高压电网建设不断对施工效率和成本提出更高的要求. 然而,由于特高压等输电线路路径走廊经常会经过交通不便甚至是人迹罕至的深山老林或崇山峻岭,输电线路建设过程中的物料运输一直是制约电网建设效率和成本的一大因素. 作为一种运输效率高、运输成本低、地形和环境适应能力强的运输方式[1],货运专用索道被越来越多地应用于输电线路施工过程中的物料运输中[2].

    为进一步降低货运索道的运输成本,提高运输效率,很多学者尝试对索道结构和设备进行研究和改进. 例如:缪谦等[3-5]研究了货运索道运输技术与设备;江明等[6-10]提出了一系列承载索和索道整体结构的方案及计算方法;白雪松等[11]开发了一种货运索道工作索计算软件,李洋等[12]设计了一种辅助索道方案优化设计平台. 而索道的架设路径也能够在很大程度上影响货运索道的物料运输效率和成本. 相对于传统的依靠人工现场勘查的索道路径规划方法,索道路径自动规划方法所考虑的备选上下料点、备选索道路径更加全面,而规划所耗费时间短,可大幅减少人工和经济成本. 在路径自动规划研究领域,王刚等[13]基于改进的A* 算法进行了输电线路路径智能选线研究;刘亮亮[14]基于蚁群算法构建了超高压输电线路路径选择的规划模型;谢景海等[15]提出了一种用于输电线路路径搜索的改进蚁群优化算法. 这些研究对货运索道的路径规划研究有较大的借鉴意义,但总体来说,现有关于货运索道的路径自动规划研究较少,且无法实现索道路径的自动规划. 例如李攀等[16]设计了一种三维GIS辅助山区输电线路货运索道选线系统,但是,由于索道支架位置自动搜索方法的缺失,该辅助选线系统仅实现了索道路径规划选线的半自动化.

    因此,本文提出了3种在二维地形剖面上进行索道支架位置自动搜索的方法:凸包点遍历法、地形自适应法和干涉点搜索法,并对这3种方法进行对比分析,提出适用于工程实际的索道支架位置搜索方法,实现输电线路货运索道路径规划的自动化.

    搜索输电线路货运索道支架位置时,需要先提供如图1所示的索道路径二维地形剖面数据. 图中坐标轴的原点O是上料点在水平面上的投影点,横坐标s经过下料点在水平面上的投影,表示水平距离,纵坐标z表示高度. 因此,索道路径二维地形剖面数据可以用坐标(st, zt)表示,其中t为地形剖面上的点. 受索道运输技术的限制,索道路径的长度一般不超过3000 m.

    图  1  索道路径二维地形剖面数据示意
    Figure  1.  Schematic of 2D terrain profile data of cableway path

    在自身重力、外部载荷和内部张力的共同作用下,相邻两个支架之间(编号为λλ + 1)的索道承载索会发生下垂,下垂后的承载索曲线g(λ,λ+1)可采用抛物线近似,如式(1)所示[17].

    z0=4fs0(ls0)l2+Cls0
    (1)

    式中:z0为承载索曲线上的点相对曲线起点的高度值;s0为承载索曲线上的点相对起点的水平距离;l为承载索曲线的终点相对起点的水平距离;C为承载索曲线的终点相对起点的高度值;f为承载索曲线的跨中垂度,其取值范围为0.050l ~ 0.080l,模拟承载索载重情况下,可采用f = 0.065l.

    在进行索道支架位置搜索时,考虑到安全性、适用性和经济性,需要对支架位置提出以下要求:

    1) 索道的支架位置必须保证承载索曲线与地面间的距离合适,以避免索道运行过程中货物触地;

    2) 相邻支架间的最大跨距不宜大于400 m;

    3) 相邻支架间的最小跨距不宜小于20 m;

    4) 中间支架数量不大于7个;

    5) 弦倾角小于45°.

    图2所示,在自动搜索索道支架位置时,需要经过以下步骤:

    图  2  索道支架位置搜索主要流程
    Figure  2.  Search process of cableway trestle position

    步骤1 索道支架位置自动搜索与冗余支架筛除. 根据索道路径二维地形剖面数据和索道承载索曲线方程,初步搜索满足索道支架位置要求1)、2)的索道支架位置,并删除多余支架以实现对支架位置的确定.

    步骤2 判断索道合理性. 根据索道支架位置要求3) ~ 5),判断索道支架位置搜索方案是否合理.

    图3所示,凸包点遍历法首先搜索出地形上所有凹度小于零的凸包点,然后在每个凸包点建立支架;再判断相邻支架间的距离,如果所有相邻支架间的距离都小于400 m,则该索道支架位置方案初步可行,后续根据其他条件判断合理性;否则,说明使用本方法在该二维地形剖面上无法建立索道.

    图  3  凸包点遍历法原理示意
    Figure  3.  Schematic of principle of convex-point traversal method

    在所有凸包点建立支架必然会造成支架数量冗余,因此,凸包遍历法还需要对冗余支架进行筛除. 对支架λλ = 1, 2, 3, …)进行冗余支架筛除的步骤如下:

    步骤1 如图4(a)所示,找出支架λ右侧且距离小于400 m的所有支架(编号为λ + 1, λ + 2, λ + 3, , γ,其中,γ为支架λ右侧所有支架的最大支架编号),并令η = λ + 2.

    图  4  冗余支架筛除示意
    Figure  4.  Schematic of deleting redundant trestles

    步骤2 如图4(b)或图4(c)所示,在支架λ和支架η之间建立承载索曲线g(λ,η),并判断曲线是否与地形干涉(当承载索曲线与地形剖面的高度差小于1 m时,认为承载索曲线与地形干涉):如果干涉(如图4(b)所示),说明支架λ和支架η之间的支架未必是冗余支架,所以直接进入下一步;如果不干涉(如图4(c)所示),说明支架λ和支架η之间的支架都是冗余支架,删除这些冗余支架(如图4(d)所示)并进入下一步.

    步骤3 令η = η + 1,重复步骤2,直至η = γ.

    地形自适应法首先在地形曲线的上方建立一条单跨索道承载索曲线(承载索曲线的起点和终点分别位于上料点和下料点正上方),然后控制索道承载索曲线起点和终点逐渐平移下落,靠近上料点和下料点. 在下落过程中,承载索曲线会与地形发生干涉. 每当干涉发生时,在干涉点的位置为索道添加一个中间支架(干涉点的位置为二维地形与承载索曲线的高度之差最大的位置),阻止干涉发生并使得承载索曲线适应地形,直至承载索曲线起点和终点分别与上料点和下料点重合.

    图5所示,该方法的原理可以通过以下步骤进行解释:

    图  5  地形自适应法原理示意
    Figure  5.  Schematic of principle of terrain adaption method

    步骤1  在上料点和下料点位置设置ad两点,建立索道承载索曲线(如图5(a)所示),判断承载索曲线是否与地形干涉. 如果不干涉,那么上料点和下料点之间只需要建立一个单跨索道g(a,d)即可;如果干涉,同时向上平移点a和点d,直到承载索曲线与地形之间仅有一个干涉点m (如图5(b)所示).

    步骤2  在该干涉点处设置中间支架m,并分别建立承载索曲线g(a,m)g(m,d). 然后,将点a和点d向下移动一定的距离,判断g(a,m)g(m,d)与地形是否干涉,并分别找到干涉位置b1c1 (如图5(c)所示).

    步骤3 分别在干涉点处设置新的中间支架b1c1,并分别建立承载索曲线g(a,b1)g(c1,d),然后,将点a和点d向下移动一定的距离,判断g(a,b1)g(c1,d)与地形是否干涉,并分别找到干涉位置b2c2 (如图5(d)所示). 重复步骤3直至点a和点d分别回到上料点和下料点.

    基于地形自适应搜索索道支架位置时,并未限制相邻支架间的距离. 为满足索道相邻支架间的最大跨距不宜大于400 m的要求,需对跨距过大的支架作进一步处理. 如图6所示,首先找出索道路径上间距大于400 m的相邻支架(如图6(a)所示),然后在这一对相邻支架的中点附近的凸包上增加新的支架(如图6(b)所示). 再找出新增支架和两侧相邻支架间的承载索曲线与地形曲线的所有干涉点,在这些干涉点位置增设新的支架(如图6(c)所示). 不断循环以上步骤优化跨距过大的相邻支架,直到所有档距都不大于400 m且所有承载索曲线与地形之间都不存在干涉点.

    图  6  支架间距优化示意
    Figure  6.  Schematic of trestle spacing optimization

    在对跨距大于400 m的相邻支架进行处理时,也会造成支架数量冗余,因此,地形自适应法也需要使用如图4所示的方法对冗余支架进行筛除.

    干涉点搜索法以已知的两个相邻支架(例如首次搜索时的上料点和下料点)为起点和终点,建立一条单跨索道承载索曲线,不断在承载索曲线与地形的最大干涉位置添加中间支架,直到承载索曲线不会与地形发生干涉;同时通过对间距大于400 m的相邻支架进行分段处理,建立满足中间支架搜索原则的货运索道.

    图7所示,该方法的原理可以通过以下步骤进行解释:

    图  7  干涉点搜索法原理示意
    Figure  7.  Schematic of principle of interference-point search method

    步骤1 在相邻支架λ和支架λ + 1之间建立承载索曲线g(λ,λ+1)(如图7(a)所示),判断曲线是否与地形干涉.

    步骤2 如果g(λ,λ+1)与地形干涉(如图7(b)所示),在干涉点处设置新的支架,编号为λ + 1 (并将原编号为λ + 1的支架及其后续支架编号顺次增大).

    步骤3 如果g(λ,λ+1)与地形不干涉,但支架λ和支架λ + 1的间距大于400 m (如图7(c)所示),那么就先找到两个支架的中点附近的凸点M,然后分别在支架λM点间以及M点和支架λ + 1间建立索道承载索曲线g(λ,M)g(M,λ+1),并分别判断g(λ,M)g(M,λ+1)曲线是否与地形干涉.

    步骤4 如果只有一条承载索曲线与地形干涉(如图7(d)所示),在相应的那个干涉点处设置新的支架,编号为λ + 1 (并将原编号为λ + 1的支架及其后续支架编号顺次增大)(如图7(e)所示);如果两条承载索曲线都与地形干涉,在相应的两个干涉点处设置两个新的支架,分别编号为λ + 1和λ + 2 (并将原编号为λ + 1的支架及其后续支架编号顺次增大);如果两条索道都不与地形干涉,则将M点设置为新的支架,编号为λ + 1 (并将原编号为λ + 1的支架及其后续支架编号顺次增大).

    步骤5 如果g(λ,λ+1)与地形不干涉,且这两点的间距小于400 m,那么支架λ和支架λ + 1之间无需设立中间支架.

    步骤6 针对所有相邻支架,重复步骤1 ~ 5,直到所有相邻支架间的承载索曲线与地形不干涉,且距离小于400 m.

    在对跨距大于400 m的相邻支架进行分段处理时,也会造成支架数量冗余,因此,干涉点搜索法也需要使用如图4所示的方法对冗余支架进行筛除.

    针对西南地区的典型地貌,共建立10万个索道路径二维地形剖面,典型特征剖面如图8中黑实线所示. 使用本文所提出的3种方法对这些二维地形剖面的索道支架位置搜索问题进行求解,求解结果样例如图8所示. 为了便于观察,在图中分别将由方法2 (地形自适应法)和方法3 (干涉点搜索法)计算得到的支架和承载索曲线在高度方向进行了平移. 在地形剖面1的计算结果中,仅有方法1 (凸包点遍历法)可以成功提供合理的索道支架位置方案,而方法2和方法3搜索的左起第3个和第4个支架之间的跨距都小于20 m,不能满足支架位置要求. 在地形剖面2 ~ 4的计算结果中,3种方法算得的大部分支架位置基本一致,只有少数支架位置有微小的差异.

    图  8  地形剖面支架位置搜索结果样例
    Figure  8.  Search results for profile examples of trestle position

    表1展示了使用本文提出的3种方法对这10万个二维地形剖面算例进行索道支架位置搜索的成功率和计算时间(计算机处理器:英特尔 Core i7-9700 @ 3.00 GHz 八核,内存为8 GB).

    表  1  本文提出的3种方法的计算结果对比
    Table  1.  Comparison of results obtained by three proposed methods
    项目方法 1 (凸包点遍历)方法 2 (地形自适应)方法 3 (干涉点搜索)
    二维地形
    剖面数/万个
    101010
    单个地形剖面平均计算时间/ms6.073.591.04
    成功搜索支架位置的索道数/条912283838255
    二维地形剖面支架位置搜索成功率/%9.128.388.26
    下载: 导出CSV 
    | 显示表格

    通过对比发现:这3种方法的支架位置搜索成功率差别不大,但单个地形剖面的平均计算时间差异较大. 方法1的支架位置搜索成功率最高,方法2次之,方法3最低;3种方法的成功率由高到低依次为9.12%、8.38%、8.26%. 方法3对单个地形剖面进行处理的平均时间最短,仅需1.04 ms;方法2次之,需3.59 ms;方法1的时间最长,为6.07 ms. 综合对比之后发现:虽然方法3的支架位置搜索成功率略低于另外两种方法,但其计算速度是方法2的3.5倍,是方法1的5.8倍. 因为3种方法的成功率并没有明显的差异,所以在实际工程应用中,推荐使用支架位置搜索速度最快的方法3.

    通过对这10万个算例的计算结果进行统计分析发现,地形剖面的水平长度、最大高度差和梯度的均方根都对支架位置搜索成功率有较大的影响,图9(a) ~ (c)分别展示了这3个因素对支架位置搜索成功率的影响. 如图9(a)所示,随着地形剖面的水平长度由200 m增加到约1600 m,3种方法的支架位置搜索成功率基本保持一致,均由50%逐渐降低到1%以下. 如图9(b)所示,随着地形剖面的最大高度差由0增加到约240 m,3种方法的支架位置搜索成功率基本保持一致,均由40%逐渐降低到1%以下. 如图9(c)所示,随着地形剖面的梯度的均方根由0.2增加到约0.8,3种方法的支架位置搜索成功率基本保持一致,均由40%逐渐降低到5%左右.

    图  9  各因素对支架位置搜索成功率的影响
    Figure  9.  Influence of various factors on success rate of trestle position search

    本文提出并对比研究了3种索道支架位置自动搜索方法:凸包点遍历法、地形自适应法和干涉点搜索法. 使用这3种方法对10万个二维地形剖面算例进行索道支架位置搜索,通过对结果进行统计分析发现:凸包点遍历法的求解速度最慢,但求解成功率最高;干涉点搜索法的求解速度最快,其求解成功率最低. 3种方法都能够在提供的二维地形剖面上自动搜索满足要求的索道支架位置,为输电线路货运索道路径自动化规划提供有效支持,进而为降低输电线路专用索道的运输成本和提高运输效率提供有力保障.

    为进一步降低输电线路施工物料的运输成本和提高运输效率,后续将以本文提出的3种索道支架位置自动搜索方法为基础,开展输电线路货运索道路径自动规划技术研究;同时结合货运索道设计计算及选型方法研究,实现货运索道的路径规划、设计计算与部件选型一体化技术研究;在此基础上,结合公路、水路等运输路径规划方法,实现输电线路施工物料运输路径的协同规划.

  • 姚玲森. 曲线梁[M]. 北京: 人民交通出版社, 1989: 1-2.
    王佳佳,丁敏,蒋秀根,等. 考虑二阶弯矩效应的自由扭转圆曲梁静力分析[J]. 中国农业大学学报,2019,24(3): 109-116. doi: 10.11841/j.issn.1007-4333.2019.03.14

    WANG Jiajia, DING Min, JIANG Xiugen, et al. Static analysis of circular curved beam with free torsion considering second-order moment effect[J]. Journal of China Agricultural University, 2019, 24(3): 109-116. doi: 10.11841/j.issn.1007-4333.2019.03.14
    铁摩辛柯 S P, 盖莱 J M. 弹性稳定理论[M]. 2版. 北京: 科学出版社, 1965: 274-281.
    项海帆, 刘光栋. 拱结构的稳定与振动[M]. 北京: 人民交通出版社, 1991: 75-87.
    窦超. 钢拱平面外稳定性能及设计方法[D]. 北京: 清华大学, 2012.
    郭彦林, 窦超. 现代拱形钢结构设计原理与应用[M]. 北京: 科学出版社, 2013: 61-86.
    窦超,郭彦林. 压弯圆弧拱平面外稳定承载力设计方法[J]. 建筑结构学报,2012,33(7): 27-36.

    DOU Chao, GUO Yanlin. Out-of-plane inelastic stability design of circular arches under combination of compression and bending[J]. Journal of Building Structures, 2012, 33(7): 27-36.
    窦超,郭彦林. 圆弧拱平面外弹性弯扭屈曲临界荷载分析[J]. 工程力学,2012,29(3): 83-89,94.

    DOU Chao, GUO Yanlin. Study on flexural-torsional buckling load of circular arches[J]. Engineering Mechanics, 2012, 29(3): 83-89,94.
    郭彦林,赵思远,窦超. 面外支撑非均匀布置的箱型与圆管截面拱平面外弹性屈曲[J]. 工程力学,2014,31(7): 29-35.

    GUO Yanlin, ZHAO Siyuan, DOU Chao. Out-of-plane elastic buckling of close-section steel arches with non-uniform distributed lateral bracings[J]. Engineering Mechanics, 2014, 31(7): 29-35.
    窦超,郭彦林. 均匀受压圆弧拱平面外弹塑性稳定设计方法[J]. 建筑结构学报,2012,33(1): 104-110.

    DOU Chao, GUO Yanlin. Out-of-plane inelastic stability and strength design of circular arches in uniform compression[J]. Journal of Building Structures, 2012, 33(1): 104-110.
    窦超,郭彦林. 受弯圆弧拱平面外稳定承载力分析[J]. 建筑结构学报,2012,33(7): 18-26.

    DOU Chao, GUO Yanlin. Out-of-plane inelastic stability of circular arches under bending moment[J]. Journal of Building Structures, 2012, 33(7): 18-26.
    邓婷. 基于大位移模型的圆拱几何非线性分析及其应用[D]. 北京: 中国农业大学, 2018.
    CHAI H Y. Flexural-torsional stability of curved beams[J]. Journal of Engineering Mechanics Division, 1982, 108(6): 1351-1369.
    CHAI H Y, PFEIFFER P A. Elastic stability of curved members[J]. Journal of Structural Engineering, 1983, 109(12): 2922-2940.
    CHAI H Y, PFEIFFER P A. Buckling of curved beams with in-plane deformation[J]. Journal of Structural Engineering, 1984, 106(2): 291-300.
    YANG Y B, KUO S R. Effect of curvature on stability of curved beams[J]. Journal of Structural Engineering, 1987, 113(6): 1185-1202. doi: 10.1061/(ASCE)0733-9445(1987)113:6(1185)
    YANG Y B, KUO S R. Curved beam elements for nonlinear analysis[J]. Journal of Engineering Mechanics, 1989, 115(4): 840-855. doi: 10.1061/(ASCE)0733-9399(1989)115:4(840)
    YANG Y B, KUO S R. Static stability of curved thin-walled beams[J]. Journal of Engineering Mechanics, 1986, 112(8): 821-841. doi: 10.1061/(ASCE)0733-9399(1986)112:8(821)
    段炼. 曲梁弯扭屈曲分析[J]. 工程力学,1989,6(3): 41-54.

    DUAN Lian. Analysis of bending and torsional buckling of curved beam[J]. Engineering Mechanics, 1989, 6(3): 41-54.
    KANG Y J, YOO C H. Thin-walled curved beams I formulation of nonlinear equations[J]. Journal of Engineering Mechanics,ASCE, 1994, 120(10): 2072-2101. doi: 10.1061/(ASCE)0733-9399(1994)120:10(2072)
    KANG Y J, YOO C H. Thin-walled curved beams II analytical solution for bucking of arches[J]. Journal of Engineering Mechanics, 1994, 120(10): 2102-2125. doi: 10.1061/(ASCE)0733-9399(1994)120:10(2102)
    许强,童根树. 任意开口薄壁截面圆弧曲梁的通用线性理论[J]. 工程力学,2002,19(6): 141-147. doi: 10.3969/j.issn.1000-4750.2002.06.028

    XU Qiang, TONG Genshu. General linear theory of curved beam with arbitrary thin-wall section[J]. Engineering Mechanics, 2002, 19(6): 141-147. doi: 10.3969/j.issn.1000-4750.2002.06.028
    杨永华. 弹性开口薄壁截面圆弧钢拱的稳定承载力研究[D]. 上海: 同济大学, 2006.
    刘磊,许克宾. 曲杆结构非线性分析中的直梁单元和曲梁单元[J]. 铁道学报,2001,23(6): 72-76. doi: 10.3321/j.issn:1001-8360.2001.06.016

    LIU Lei, XU Kebin. Curved-beam element and straight-beam element used in the nonlinear analysis of curved frame structures[J]. Journal of the China Railway Society, 2001, 23(6): 72-76. doi: 10.3321/j.issn:1001-8360.2001.06.016
    刘磊. 大跨度混凝土桥梁的双非线性分析[D]. 北京: 北京交通大学, 2000.
    程鹏. 两铰圆弧拱非线性弯曲理论和弹塑性稳定[D]. 杭州: 浙江大学, 2005.
    HEINS C P. Bending and torsional design in structural members[M]. [S.l.]: Heath and Company, 1975: 124-125.
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(5)  / Tables(4)

    Article views(692) PDF downloads(31) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return