Optimization and Tensile Properties of Composite Insulator at High and Low Temperature

YAN Zhitao; MA Wenjun; ZHANG Pu; QIU Zhe

doi:10.3969/j.issn.0258-2724.20181032

Volume 56 Issue 1

Jan. 2021

Turn off MathJax

Article Contents

Abstract

References

Journal of Southwest Jiaotong University > 2021 > 56(1): 56-61, 83.

YANG Yongzhi. Shape Analysis of Main Cable of Single Tower Suspension Bridge with Unilateral Spatial Cable Plane and Curved Beam[J]. Journal of Southwest Jiaotong University, 2024, 59(2): 298-306. doi: 10.3969/j.issn.0258-2724.20230197

Citation:

YAN Zhitao, MA Wenjun, ZHANG Pu, QIU Zhe. Optimization and Tensile Properties of Composite Insulator at High and Low Temperature[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 56-61, 83. doi: 10.3969/j.issn.0258-2724.20181032

Citation:

PDF( 3489 KB)

Optimization and Tensile Properties of Composite Insulator at High and Low Temperature

doi: 10.3969/j.issn.0258-2724.20181032

Received Date: 25 Dec 2018
Rev Recd Date: 16 Feb 2020

Available Online: 18 Sep 2020

Publish Date: 01 Feb 2021

Abstract

Abstract

In order to find out the influence of high and low temperature and the crimping process between the core rod and the end fittings on the bearing capacity of composite insulator, and to get the final optimized crimping mode, the finite element simulation analysis was carried out on the high and low temperature tensile test results of the core rod and the end fittings.Firstly, the stress of the fittings and mandrels was simulated at normal temperature, and the results were in good agreement with the experimental data. Then, 10 kinds of crimping conditions were designed to evaluate the optimal preloading conditions of the end fittings of composite insulator. The results show that when the preloading stress is evenly distributed along the radial direction of the fittings and 18% − 25% of the length of the fittings is reserved, the bearing capacity of the composite insulator link is the best, and the elastic limit load of the optimal model is 8.23% higher than that of the existing manufacturers. The simulation method can also be used to simulate the mechanical properties of insulators at high and low temperature, but the compression optimization cannot improve the bearing capacity significantly.
- composite insulators,
- fitting,
- preload,
- tensile properties,
- optimization

FullText(HTML)

单边空间索面曲梁悬索桥具有造型简约优雅、通透感强等优点，是景区或城市重要节点人行桥选型有竞争力的方案，国内自2013年上海迪斯尼湖畔公园景观人行桥^[1]建成后，陆续修建了多座相似桥梁^[2-3].

在单边空间索面悬索桥成桥状态主缆线形分析过程中，主缆的线形、吊杆内力及方向均为未知量，且存在耦合，这使得主缆线形求解较为复杂.

传统的主缆计算理论（如抛物线理论、分段悬链线理论等）都是围绕主缆的微分方程展开研究，通过积分或者数值方法得到主缆线形^[4]. 此过程一般假定水平力为固定值，并不适用于单边空间索面悬索桥主缆线形分析，因为斜吊杆等会引起主缆水平力改变.

关于空间索面的悬索桥主缆成桥线形分析的研究相对较少，主要是基于分段悬链线的空间分段悬链线法和数值解析法. 罗喜恒等^[5]将分段悬链线法进行拓展，用以求解空间索面的悬索桥主缆. 彭苗等^[6-9]以空间分段悬链线理论为基础，对收敛困难问题进行优化. 李传习等^[10-11]提出空间缆索悬索桥数值解析法，以索段状态方程、索段平衡方程式和分点平衡条件构成空间主缆精确找形基本方程，结合搜索合适起始端初始内力，通过迭代完成对主缆的找形计算. 张文明等^[12-13]提出基于分段悬链线理论和非线性GRG （generalized reduced gradient）法的主缆线形解析计算方法. 空间分段悬链线法和基于分段悬链线的数值解析法参数少，收敛速度快，是目前计算精度最高、应用最普遍的方法，但这类方法原理上需要引入包含索段无应力长度及变形相容条件的索段状态方程，由于缆索结构的非线性，使得初值的选择对迭代收敛性及收敛速度的影响较大.

马召宇等^[14-16]基于ANSYS有限元模型，构造坐标迭代模型，求解主缆线形. 坐标迭代法和空间分段悬链线法存在共性的问题：在迭代过程中应用变形相容条件，缆索结构强烈的非线性对初值要求较高.

孙远等^[17]提出“合理轴线”思路，基于桁架单元和结构力学的基本原理，建立单元坐标与杆端力的关系式；刘超等^[18]提出基于整体力学分析的主缆坐标体系法，分别对塔顶、地锚点及主跨跨中取矩，建立整体力学平衡方程，求解地锚支点反力，再通过迭代求解出主缆线形坐标. 这2种思路都将无应力长度及变形相容条件后置，大大降低对初值的要求.

由上述研究可知，无应力长度及变形相容条件不是求解主缆线形的必要条件，而且如果将其后置，解除缆索长度变形与空间变位的耦合关系，将大大降低对迭代初值的要求，提高算法普适性.

本文借鉴等代梁法思路^[19]，基于荷载与主缆几何线形的关联方程，用数值迭代求解主缆线形，并将计算流程拆分成粗算和精算2个阶段. 其中，粗算阶段通过合理简化，降低初值标准，既能提高算法普适性，又能得到“足够精度确保收敛”的精算阶段初值；精算阶段完成空间主缆线形的精确求解.

1. 空间主缆几何线形关联迭代法

1.1 基本假定与约定

分析主缆受力时采用如下假定：

1）主缆为理想柔性索，不考虑其抗弯刚度及抗扭刚度；

2）主缆受到的吊杆集中力只作用在节点处.

坐标系、力和力矩方向约定如图1所示，以地锚点为原点，与加劲梁的吊点连成的平面曲线的切线为x轴，竖直向上为z轴正向，按右手法则确定坐标系. 力和力矩的方向均以坐标轴同向为正. 点A为主缆在塔顶锚固点，点C为主缆跨中垂点，点B坐标（x,y,z），为主缆上任意点. 主缆跨径为L，两端点z轴方向高差为h_z，y轴方向坐标差为h_y，跨中垂度为f_z. 成桥状态悬索桥主缆承受荷载有2种：沿主缆索长分布的均布荷载函数q_t(x) （如主缆自重、主缆检修道等）及主缆受到的吊杆集中荷载函数P_t(x) （如通过吊杆传递的加劲梁恒载、索夹、吊杆及锚具自重等）. 其中：下标$ t\in \left\{x,y,{\textit{z}}\right\} $表示荷载方向， P_t(x_i)、P_t(x_j)分别为作用在x轴x_i、x_j处的主缆集中力，F_tl为地锚点O支点反力.

图 1 单跨空间主缆计算

Figure 1. Single-span spatial main cable calculation

下载: 全尺寸图片幻灯片

1.2 外荷载节点力和力矩

定义主缆上任意节点左侧外荷载在该点产生的节点力为外荷载节点力F_t(x)，产生的节点力矩为外荷载节点力矩M_t(x).

如图2所示，取任意空间主缆单元，计算单元右端的外荷载节点力为

图 2 离散单元内力计算

Figure 2. Internal force calculation of discrete element

下载: 全尺寸图片幻灯片

，

$\begin{split} &{F}_{t}({x}_{j})={F}_{t}({x}_{i}) +\frac{{q}_{t}({x}_{i}) + {q}_{t}({x}_{j})}{2}{K}_{i}\left({x}_{j}-{x}_{i}\right) + {P}_{t}({x}_{j})\text{，} \end{split}$

(1)

式中：K_i为计算单元的荷载分布长度系数，分布荷载沿主缆索长分布时，近似计算如式（2）所示，其中，y_i、y_j和z_i、z_j分别为节点i、j的y轴和z轴坐标.

${K}_{i}=\sqrt{1 + {\left(\frac{{y}_{j}-{y}_{i}}{{x}_{j}-{x}_{i}}\right)}^{2} + {\left(\frac{{{\textit{z}}}_{j}-{{\textit{z}}}_{i}}{{x}_{j}-{x}_{i}}\right)}^{2}}\text{.}$

(2)

假设作用在节点i的虚拟力V_t(x_i) （式（3））在节点j产生的力矩与节点j处的外荷载力矩相等，将此虚拟力定义为外荷载力矩等代力.

${V}_{t}({x}_{i})={F}_{t}({x}_{i}) + \frac{2{q}_{t}({x}_{i}) + {q}_{t}({x}_{j})}{6}{K}_{i}\left({x}_{j}-{x}_{i}\right).$

(3)

则外荷载节点力矩迭代式为

$\left\{ \begin{array}{c}{M}_{x}\left({x}_{j}\right)={M}_{x}\left({x}_{i}\right) + {V}_{y}\left({x}_{i}\right)\left({{\textit{z}}}_{j}-{{\textit{{\textit{z}}}}}_{i}\right)-{V}_{{\textit{z}}}\left({x}_{i}\right)\left({y}_{j}-{y}_{i}\right),\\ {M}_{y}\left({x}_{j}\right)={M}_{y}\left({x}_{i}\right) + {V}_{{\textit{z}}}\left({x}_{i}\right)\left({x}_{j}-{x}_{i}\right)-{V}_{x}\left({x}_{i}\right)\left({{\textit{z}}}_{j}-{{\textit{z}}}_{i}\right),\\ {M}_{{\textit{z}}}\left({x}_{j}\right)={M}_{{\textit{z}}}\left({x}_{i}\right) + {V}_{x}\left({x}_{i}\right)\left({y}_{j}-{y}_{i}\right)-{V}_{y}\left({x}_{i}\right)\left({x}_{j}-{x}_{i}\right). \end{array}\right.$

(4)

由式（４）可见，任意方向力矩均可表示为三部分：前节点力矩、主方向力矩和耦合方向力矩. 当计算精度要求不高时，可通过略去较小的耦合方向力矩，实现对竖向z轴和水平y轴方向的解耦.

初始端边界值为

$\left\{\begin{array}{c}{F}_{t}\left(0\right)=0,\\ {M}_{t}\left(0\right)=0. \end{array}\right.$

(5)

1.3 主缆线形方程

根据理想柔索假定，主缆上任意点力矩均为0，因此，依次对主缆上任意点（x,y,z）的x、y、z轴方向取矩，可得

${F}_{y\mathrm{l}}{\textit{z}}-{F}_{{\textit{z}}\mathrm{l}}y + {M}_{x}\left(x\right)=0\text{，}$

(6)

${F}_{{\textit{z}}\mathrm{l}}x-{F}_{x\mathrm{l}}{\textit{z}} + {M}_{y}\left(x\right)=0\text{，}$

(7)

${F}_{x\mathrm{l}}y-{F}_{y\mathrm{l}}x + {M}_{{\textit{z}}}\left(x\right)=0.$

(8)

消掉F_xl、F_yl、F_zl，得到荷载与主缆几何线形的关联方程为

${M}_{x}\left(x\right)x + {M}_{y}\left(x\right)y + {M}_{{\textit{z}}}\left(x\right){\textit{z}}=0.$

(9)

定义主缆坐标的收敛完成率为

$p=1-\frac{\left|{M}_{x}\left(x\right)x + {M}_{y}\left(x\right)y + {M}_{{\textit{z}}}\left(x\right){\textit{z}}\right|}{\left|{M}_{x}\left(x\right)x\right| + \left|{M}_{y}\left(x\right)y\right| + \left|{M}_{{\textit{z}}}\left(x\right){\textit{z}}\right|}.$

(10)

主缆坐标的收敛完成率是评价主缆坐标迭代计算完成度的指标.

如图1，主缆右端点坐标为(L,h_y,h_z)，跨中垂点坐标为(L/2, y, f_z)，将右端点和跨中垂点的x、z轴坐标分别代入式（7），可得

$\left\{\begin{array}{l}{F}_{{\textit{z}}\mathrm{l}}L-F _{x\mathrm{l}}{h}_{{\textit{z}}} + {M}_{y}\left(L\right)=0\text{，}\\ {F}_{{\textit{z}}\mathrm{l}}\dfrac{L}{2}-F _{x\mathrm{l}}\left(\dfrac{{h}_{{\textit{z}}}}{2}-{f}_{{\textit{z}}}\right) + {M}_{y}\left(\dfrac{L}{2}\right)=0. \end{array}\right.$

(11)

消元并整理后，求得地锚支点反力x轴向分力为

${F}_{x\mathrm{l}}=\frac{{M}_{y}\left(L\right)-2{M}_{y}\left(\dfrac{L}{2}\right)}{2{f}_{{\textit{z}}}}.$

(12)

把主缆右端点坐标代入式（7）、（8），分别求得

${F}_{{\textit{z}}\mathrm{l}}=\frac{{F}_{x\mathrm{l}}{h}_{{\textit{z}}}-{M}_{y}\left(L\right)}{L}\text{，}$

(13)

${F}_{y\mathrm{l}}=\frac{{F}_{x\mathrm{l}}{h}_{y} + {M}_{{\textit{z}}}\left(L\right)}{L}.$

(14)

变换式（7）、（8），得到主缆线形方程为

$\left\{\begin{array}{c}y=\dfrac{{F}_{y\mathrm{l}}x-{M}_{{\textit{z}}}\left(x\right)}{{F}_{x\mathrm{l}}}\text{，}\\ {\textit{z}}=\dfrac{{F}_{{\textit{z}}\mathrm{l}}x + {M}_{y}\left(x\right)}{{F}_{x\mathrm{l}}}. \end{array}\right.$

(15)

将式（12）～（14）代入式（15），得到主缆线形方程的另一种表达式为

$\left\{ \begin{array}{c} y=\dfrac{2{f}_{{\textit{z}}}}{{M}_{y}\left(L\right)-2{M}_{y}\left(\dfrac{L}{2}\right)}\left(\dfrac{{M}_{{\textit{z}}}\left(L\right)}{L}x-{M}_{{\textit{z}}}\left(x\right)\right) + \dfrac{{h}_{y}}{L}x\text{，}\\ {\textit{z}}=-\dfrac{2{f}_{{\textit{z}}}}{{M}_{y}\left(L\right)-2{M}_{y}\left(\dfrac{L}{2}\right)}\left(\dfrac{{M}_{y}\left(L\right)}{L}x-{M}_{y}\left(x\right)\right) + \dfrac{{h}_{{\textit{z}}}}{L}x. \end{array}\right.$

(16)

传统抛物线理论的主缆线形方程为式（16）在沿跨径方向均布荷载作用下的特例，主缆上任意点的外荷载节点力矩为

${M}_{y}\left(x\right)=\frac{1}{2}q{x}^{2}.$

(17)

将式（17）代入式（16）的z轴坐标表达式，即可得到常见的主缆线形抛物线方程为

${\textit{z}}=-\frac{4{f}_{{\textit{z}}}}{L^2}x\left(L-x\right) + \frac{{h}_{{\textit{z}}}}{L}x.$

(18)

式（15）、（16）分别适用于编程和Excel表格计算.

1.4 主缆无应力长度

1.4.1 传统方法

传统分段悬链线法及数值解析法计算主缆无应力长度方法较为复杂，引用李传习等^[20]成果如式（19）所示.

$\begin{split} & l=-\frac{H{S} _{0}}{EA}-\frac{H}{q}\Bigg[\mathrm{l}\mathrm{n}\left(V + \sqrt{{H}^{2} + {V}^{2}}\right)-\Bigg.\\ &\quad \left.\mathrm{l}\mathrm{n}\left(V-q{S} _{0} + \sqrt{{H}^{2} + {\left(V-q{S} _{0}\right)}^{2}}\right)\right]\text{，} \end{split}$

(19)

式中：l为索段左右两端点水平距离，H、V分别为索段左端张力的水平分力、竖直分力，S₀为索段的无应力长度，q为单位长度索段的重力，E为索段截面弹性模量，A为索段截面面积.

假设式（19）中除S₀外，其余均为已知量，可通过迭代求解主缆的无应力长度. 实际求解过程中，索段张力分力H、V往往是未知量，同样需要迭代求解，主缆无应力长度的求解相当复杂.

1.4.2 几何关联迭代法

几何关联迭代法不需要使用变形相容条件求解线形，因此，可以将单元的无应力长度求解过程后置于主缆线形求解之后. 这一方面可以降低对初值要求，另一方面可以在求解无应力长度时通过临时加密计算节点数，提高计算精度至满足要求.

主缆节点i (x_i, y_i, z_i)的节点张力分力为

${T}_{ti}={F}_{t\mathrm{l}} + {F}_{t}\left({x}_{i}\right).$

(20)

主缆节点i的张力为

${T}_{i}=\sqrt{{T}_{xi}^{2} + {T}_{yi}^{2} + {T}_{{\textit{z}}i}^{2}}.$

(21)

单元伸长量为

$\Delta {L}_{i}=\frac{{(T}_{i-1}{ + T}_{i}){L}_{i}}{2EA}.$

(22)

单元无应力长度为

${L}_{{\mathrm{u}}i}={L}_{i}-\Delta {L}_{i}.$

(23)

主缆无应力长度S₀计算式为

${S} _{0}=S-\sum_i \Delta {L}_{i}\text{，}$

(24)

式中：S为主缆的形状长度.

1.5 几何关联迭代法的应用

几何关联迭代法求解主缆线形需要如下参数：

1）节点坐标的初值序列，其中，x坐标应满足递增，并完整覆盖主缆跨度L，粗算阶段y和z可随意取值，精算阶段使用粗算阶段坐标结果作为初值；

2）主缆的外荷载q_t (x)和P_t (x)；

3）主缆两端z轴方向坐标差h_z、跨中垂度f_z、主缆两端y坐标差h_y （通常根据桥塔临近吊杆倾角、高度及用地情况确定）.

具体计算流程如图3所示.

图 3 几何关联迭代法求解主缆线形流程

Figure 3. Flowchart of solving main cable shape using geometric correlation iteration method

下载: 全尺寸图片幻灯片

几何关联迭代法循环流程简单，迭代收敛速度快，适用于Excel或者编程计算.

2. 单边空间索面主缆线形计算

缺少边界条件的单边空间索面主缆的合理成桥状态并非唯一，需要引入合理的约束以得到符合设计意图的索形.

吊杆的竖向分力是平衡曲线主梁结构重力最有效的途径，故可将吊杆竖向力作为找形过程中的不变量，相应地将吊杆倾角作为自变量. 吊杆张力初始值按吊杆分担桥面面积对应的一期及二期恒载确定，并根据钢箱梁结构验算适当调整. 吊杆最小倾角取为０，即垂直吊杆，最大倾角由景观设计师综合考虑用地及造型确定.

主缆线形采用几何关联迭代法，计算过程存在嵌套迭代：吊杆倾斜角度与主缆坐标匹配度的外层迭代；固定吊杆倾斜角度后，计算主缆坐标的内层迭代. 计算过程的迭代复杂度是O(n²)，对初值精度较为敏感，其中，n为迭代次数.

为降低初值需求、加速收敛，把计算流程分解成粗算阶段和精算阶段.

2.1 粗算阶段

粗算阶段通过合理简化，快速构造“足够精度确保收敛”的初值序列，主要措施如下：

1）解耦：观察主缆线形方程（式（16））及外荷载力矩（式（4））可知，外荷载与线形的耦合关系是在计算外荷载力矩过程中引入的，可以通过对式（4）进行合理简化，消除耦合关系，将迭代复杂度降低成O(n). 略去吊杆在x轴方向分力，解除x轴方向力与y、z坐标耦合后，外荷载力矩M_y(x)只与z轴方向外荷载相关，M_z(x)只与y轴方向外荷载相关，主缆线形计算拆分成竖平面线形计算和水平面线形计算2个相互独立的模块. 解耦处理后y、z坐标计算过程不再需要迭代，实质上也不需要设置初值，只需依据节点划分设定x坐标，并确保x坐标满足递增并覆盖跨径L即可.

2）集中荷载转分布荷载：粗算阶段需要通过调整吊杆倾角，达到收敛的目的，初始迭代时误差较大，如吊杆作为集中力调整，容易出现收敛计算反复横跳情况，将集中力转化为分布荷载后，大大降低荷载调幅难度，达到快速收敛的目的. 其中，吊杆竖直z向分力分散为梯形荷载，吊杆水平y向分力分散为曲边三角荷载，并均按吊杆间距分摊到主缆上.

3）沿主缆索长方向的分布荷载简化为沿主缆跨径方向的分布荷载，消除分布荷载与主缆y、z坐标的耦合性.

2.1.1 竖平面线形计算

粗算阶段简化后，计算单元的荷载分布长度系数K_i按式（25）抛物线弧长公式计算，并取代式（2）.

${K}_{i}=1 + \frac{{h}_{{\textit{z}}}^{2}}{2{L}^{2}} + \frac{8{f}_{{\textit{z}}}^{2}}{3{L}^{2}}.$

(25)

吊杆力集中荷载转化为分布荷载后，主缆的竖平面分布荷载为

${q}_{{\textit{z}}}\left(x\right)={{K}_{i}{q}_{{\mathrm{s}}} + q}_{{\mathrm{d}}}\left(x\right)\text{，}$

(26)

式中：q_s为主缆单位长度重量，按常规正桥单索面垂直吊杆模型测算主缆张力并选择主缆规格；q_d(x)为吊杆集中荷载转化的分布荷载函数.

解耦后的竖平面线形计算无需迭代，坐标初值只需满足递增，并完整覆盖主缆跨度. 按图3流程即可求解竖平面线形.

2.1.2 水平面线形计算

解耦后，固定荷载作用下的水平面线形计算同样无需迭代，原有的嵌套迭代只保留吊杆倾斜角度与主缆坐标匹配度的外层迭代. 迭代坐标初值x坐标与竖平面线形计算保持一致，y坐标随意，迭代流程如图4所示.

图 4 水平面线形计算迭代流程

Figure 4. Flowchart of planar shape calculation iteration

下载: 全尺寸图片幻灯片

水平面线形计算围绕吊杆在水平面投影长度的计算构造迭代，投影长度计算有2种方法：

1）根据y向分布荷载与竖向分布荷载比例求出吊杆倾角，再根据竖平面计算得到的吊杆高度求出吊杆在水平面投影长度L_s1；

2）已知y向分布荷载，按式（1）～（13）求出吊杆在主缆上的吊点坐标，再结合桥面钢箱梁边缘的吊点坐标求出吊杆在水平面的投影长度L_s2.

迭代的收敛条件是上述2种方法计算的吊杆在水平面的投影长度差值符合精度要求.

迭代的变量是水平面y向分布荷载. 初始迭代轮次的水平面y向分布荷载根据吊杆初始倾角及吊杆竖向分布荷载计算确定，后续迭代轮次的y向分布荷载q_y(x)均由前一迭代轮次计算采用值q_yF(x)和荷载增量 ∆q_y(x)组成，如式（27）所示.

$q_y\left(x\right)=q_{y{\mathrm{F}}}\left(x\right)+\Delta q_{y}\left(x\right).$

(27)

荷载增量 ∆q_y(x)采用前一迭代轮次计算的投影长度差值对应的荷载差值，并尽可能让当前迭代轮次的节点投影长度差低于前一迭代轮次最大误差2个数量级即可. 修改任意节点的荷载增量都会影响所有节点的y坐标，但是按上述原则保持单方向修改可以保证误差逐步收敛.

2.1.3 阶段成果

粗算阶段最终成果包括：由竖平面线形计算和水平面线形计算成果汇总成各节点坐标、竖向z轴和水平向y轴的分布荷载.

2.2 精算阶段

精算阶段应将粗算阶段成果中的吊杆分布荷载合并成集中荷载.

精算阶段需要进行嵌套迭代计算，其中，外层迭代是主缆坐标变化与吊杆倾斜度变化匹配迭代. 外层迭代流程如图5所示.

图 5 精算阶段外层迭代流程

Figure 5. Flowchart of outer iteration in precise calculation stage

下载: 全尺寸图片幻灯片

内层迭代是固定外荷载作用下的主缆线形计算，采用图3计算流程.

至此，通过粗算阶段和精算阶段的逐步迭代，最终完成主缆线形计算.

3. 算　例

如图6所示，厦门山海步道和美桥为单边索面独塔人行悬索桥，主跨216.7 m，采用V型桥塔，塔高53.56 m，桥塔菱形钢箱截面，最大尺寸2.0 m×1.6 m，根部通过球铰支座与基础连接，南北两段主缆以V塔为对称轴，主缆的低端锚固于地锚，高端锚固于V塔塔顶，V塔双柱之间设桥塔连接索，主缆采用全封闭高钒索，截面积A=129.69 cm². 南北主缆各有17根吊索，吊索间距6 m. 主缆两端z轴方向坐标差h_z=30 m、跨中垂度f_z=9 m、主缆两端y坐标差h_y=5 m.

图 6 厦门山海步道和美桥

Figure 6. Hemei bridge, Shanhai Promenade in Xiamen

下载: 全尺寸图片幻灯片

人行桥面位于半径167.65 m的圆曲线范围内，采用钢箱梁，宽度4.4 m，梁高1.2 m，钢箱吊杆吊点偏心布置于钢箱外边缘.

坐标系以北缆地锚点为坐标原点，地锚点与首根吊杆的钢箱梁吊点在水平面上投影点连线为x轴方向，竖直向上为z轴方向，坐标轴方向满足右手定则.

3.1 粗算阶段

解耦后固定外荷载作用下的竖平面线形不用迭代即可直接求出坐标.

水平面线形计算需要根据图4流程迭代，寻找满足要求的水平荷载. 吊杆倾角初值的拟定综合考虑受力、景观及用地因素，本桥设计时，希望吊杆倾角呈现有节奏的变化，最靠近塔顶的吊杆倾斜角度控制在45° 左右，最矮的吊杆为垂直吊杆，中间的倾角按线性变化，根据吊杆倾角和吊杆竖向力可以求出水平分布荷载初值及吊杆的顶端点在水平面上的投影连线.

1/2桥平面投影如图7所示，点P是主缆在塔顶锚点在水平面的投影，点E是首根吊杆，坐标轴Ox与曲线钢箱梁在点E相切，①～⑤ 是投影线. ⑤ 线是最终主缆在水平面投影线，迭代过程如下：

图 7 1/2桥平面投影

Figure 7. Projection of 1/2 bridge on plane

下载: 全尺寸图片幻灯片

第1轮迭代：根据倾角渐变原则逐一计算吊杆倾角，再通过竖平面线形计算得到主缆z坐标，即可计算出吊杆顶端点在水平面投影线 ① 线（即L_s1）；根据吊杆的竖向力和倾角求出y方向分布荷载初值，按式（1）～（16）求出吊杆在主缆上的吊点坐标，结合桥面钢箱梁边缘的吊点坐标求出吊杆顶端点在水平面的投影线 ② 线（即L_s2），这2根线最大误差为3.38 m；

第２轮迭代：按式（27）对y向荷载进行调幅，由y向分布荷载与竖向分布荷载比值求出吊杆倾角，再结合主缆z坐标可求出荷载调幅后吊杆顶端点在水平面的投影线 ③ 线（即L_s1）；再根据y向荷载计算出的主缆上吊点坐标，求出吊杆顶端点在水平面投影 ④ 线（即L_s2），2根线误差收敛到0.102 m. 后续轮次迭代模式均与第2轮次相同.

依次经过5轮迭代，最大误差依次为3380.2、101.8、15.6、1.6、0.2 mm，最终收敛到0.2 mm.

3.2 精算阶段

从粗算阶段转入精算阶段，吊杆荷载由分散荷载变成集中荷载，精算阶段首轮迭代保持荷载不变，直接进行内层迭代，修正坐标，首轮迭代坐标偏差组合值为6.2 cm.

精算阶段从第2轮迭代开始，均按图5流程先调荷载再调坐标.

精算阶段进行了5轮迭代，由于初值足够精确，所以每轮迭代都是内外层各计算一次，如表1 所示.

表 1 精算阶段迭代计算情况

Table 1. Calculation iteration in precise calculation stage

轮次	坐标偏差组合值/mm	吊杆 x 轴向分力差值/kN	吊杆 y 轴向分力差值/kN
1	62.0	0.0197	1.1095
2	1.0	0.0194	1.1095
3	0.1	0	0.0935
4	7.1×10⁻³	9.7000×10⁻⁸	0.0263
5	1.6×10⁻³	2.6000×10⁻⁸	0.0073

下载: 导出CSV

| 显示表格

主缆各索段无应力长度L_ui按式（23）计算，计算结果为S=124.2791 m，累计伸长量$\displaystyle\sum\nolimits_i^{} {} $ΔL_i=0.1943 m，S₀=124.0848 m.

3.3 与有限元软件成果对比

为验证计算方法的正确性，采用空间有限元软件Midas/Civil进行校验验算.

按精算阶段的成果输入主缆坐标，定义材料及截面，生成单元并设定索单元的初拉力，建立空间有限元模型，非线性分析收敛条件设定为10⁻⁶，进行非线性迭代分析计算，经过3轮迭代，Midas/Civil模型各节点位移收敛，最大组合位移误差小于0.1 mm.

对比几何关联迭代法结果和Midas/Civil计算结果可以看出，两者节点位置和张力差距很小，节点位置最大误差0.409 mm，节点张力最大误差0.14 kN，各节点位置和节点张力误差如图8所示.

图 8 节点的位置误差和张力误差

Figure 8. Errors of node position and node tension

下载: 全尺寸图片幻灯片

4. 结　论

本文提出基于外荷载与主缆线形的几何关联方程迭代法求解空间主缆线形，研究空间主缆线形计算新的思路，将基于分段悬链线理论的计算方法中容易引起迭代收敛失败的变形相容条件从求解必要条件中剥离，简化了求解难度，提高算法普适性.

针对于空间主缆与吊索相互耦合，提出分阶段逐步逼近的思路，粗算阶段通过解耦降低初值要求，求得高质量的精算阶段初值，精算阶段通过迭代计算得到空间主缆线形精确解.

几何关联迭代法的特点既适用于普及率较高的Excel编公式迭代，也可通过编程作为通用有限元程序的插件，为有限元程序提供成果数据.

References(16)

References

徐如恬. 瓷绝缘子强度及瓷材料种类强度水平的探讨[J]. 电瓷避雷器,1988(1): 1-8.

XU Rutian. Discussion on porcelain insulator strength and porcelain material type strength level[J]. Electric Porcelain Arrester, 1988(1): 1-8.

郭晨鋆,王科,颜冰,等. 线弹性断裂力学的支柱瓷绝缘子断裂特性及断裂机理研究[J]. 云南电力技术,2018,46(2): 12-18.

GUO Chenyun, WANG Ke, YAN Bing, et al. Fracture characteristics and fracture mechanism of pillar porcelain insulators based on linear elastic fracture mechanics[J]. Yunnan Electric Power Technology, 2018, 46(2): 12-18.

DHANESWARA D, SUHARNO B, NUGRAHA N D, et al. Effects of ceramic fibre insulation thickness on skin formation and nodule characteristics of thin wall ductile iron casting[C]//IOP Conference Series: Materials Science and Engineering, International Conference on Advanced Materials for Better Future. Surakarta: IOP Publishing Ltd., 2017, 176: 12-32

张硕,姚宁,吴继平,等. 玻璃纤维增强环氧树脂复合材料的力学性能[J]. 电工材料,2016(1): 11-14.

ZHANG Shuo, YAO Ning, WU Jiping, et al. Mechanical properties of glass fiber reinforced epoxy composites[J]. electrical materials, 2016(1): 11-14.

贾红,张鹏. 玻璃纤维增强塑料低温收缩率的测量[J]. 上海交通大学学报,2016,50(4): 502-505.

JIA Hong, ZHANG Peng. Measurement of low temperature shrinkage of glass fiber reinforced plastics[J]. Journal of Shanghai Jiaotong University, 2016, 50(4): 502-505.

PARIS L, PARGAMIN L, DUMORORA D, et al. Rating of composite suspension insulators related to the long term mechanical strength of rods[J]. IEEE Transactions on Power Delivery, 1994, 9(4): 2055-2063. doi: 10.1109/61.329538

DEMIDOV A N, KARIMBEKOV M A, MARCHENKOV A Y. Heating and cooling impact on mechanical properties of rip electric insulator for high voltage inputs[J]. Solid State Phenomena, 2017, 265: 496-500. doi: 10.4028/www.scientific.net/SSP.265.496

DE TOURREIL C, ROBERGE R, BOURDON P. Long-term mechanical properties of high voltage composite insulators[J]. Power Apparatus & Systems IEEE Transactions on, 1985, 104(10): 2918-2921.

余涛,杨新华. 压接式复合绝缘子接头处的破坏分析[J]. 水电能源科学,2010,28(10): 139-141.

YU Tao, YANG Xinhua. Failure analysis of joints of compression type composite insulators[J]. Hydropower and Energy Science, 2010, 28(10): 139-141.

袁骏. 特高压线路用复合绝缘子的机械性能[J]. 电网技术,2006(12): 29-32.

YUAN Jun. Mechanical properties of composite insulators for UHV transmission lines[J]. Power Grid Technology, 2006(12): 29-32.

钟万才. 超高压输电线路架空地线复合绝缘子断裂分析[J]. 广西电力,2014,37(5): 48-50.

ZHONG Wancai. Fracture analysis of composite insulators for overhead ground wires of EHV transmission lines[J]. Guangxi Electric Power, 2014, 37(5): 48-50.

杨祉豪. 一种通用型绝缘子卡具夹头力学及有限元分析[J]. 内燃机与配件,2019(3): 46-47.

YANG Weihao. Mechanical and finite element analysis of a universal insulator clamp[J]. Internal Combustion Engines and Accessories, 2019(3): 46-47.

陈显贻,刘志强. 棒形悬式复合绝缘子金具对芯棒不同连接结构与机械拉伸强度的关系[J]. 电瓷避雷器,2002(5): 9-13.

CHEN Xianyi, LIU Zhiqiang. The relationship between the different connecting structures of mandrel and the mechanical tensile strength of the rod suspension composite insulator fittings[J]. Electric Porcelain Arrester, 2002(5): 9-13.

谢占山. 500 kV“V”串复合绝缘子机械疲劳试验研究[D]. 北京: 华北电力大学, 2012.

MANE J V, CHANDRA S, SHARMA S, et al. Mechanical property evaluation of polyurethane foam under quasi-static and dynamic strain rates-an experimental study[J]. Procedia Engineering, 2017, 173: 726-731. doi: 10.1016/j.proeng.2016.12.160

KUMOSA M, HAN Y, KUMOSA L. Analyses of composite insulators with crimped end-fittings:part I—nonlinear finite element computations[J]. Composites Science & Technology, 2002, 62(9): 1191-1207.

Relative Articles

[1]	ZHANG Xueqin, GAO Runming, GUO Yujun, KANG Yongqiang, LI Yuansheng, WU Guangning. Hyperspectral-Based Corona Aging Evaluation for Composite Insulators[J]. Journal of Southwest Jiaotong University, 2020, 55(2): 442-449. doi: 10.3969/j.issn.0258-2724.20181062
[2]	QU Hui, PAN Hongwei, LIU Yanzhi, TAN Ding. Tensile Behaviour of Hollo-bolts Embedded in Concrete[J]. Journal of Southwest Jiaotong University, 2017, 30(5): 985-993. doi: 10.3969/j.issn.0258-2724.2017.05.020
[3]	LIU Xin, ZHONG Weizhi. Optimization and Performance Analysis for Bandwidth Spectrum Sensing in Cognitive Radio[J]. Journal of Southwest Jiaotong University, 2015, 28(1): 66-71. doi: 10.3969/j.issn.0258-2724.2015.01.010
[4]	SUN Jixing, LUO Shucai, WU Guangning, GAO Guoqiang, YANG Kunsong, HU Xueyong. Pollution Distribution Characteristic of High-Speed Train Roof Insulator and Its Optimization Strategy[J]. Journal of Southwest Jiaotong University, 2014, 27(2): 343-350. doi: 10.3969/j.issn.0258-2724.2014.02.024
[5]	WANG Gang, RONG Jian, DING Tianbao, . Optimization of Time Parameters in Static Segment in FlexRay Network[J]. Journal of Southwest Jiaotong University, 2012, 25(2): 318-324. doi: 10.3969/j.issn.0258-2724.2012.02.024
[6]	XUE Feng, WANG Ciguang, ZHANG Zhanjie, 2. Optimization Algorithm for Wagon-Flow Allocation in Marshalling Station[J]. Journal of Southwest Jiaotong University, 2010, 23(6): 932-937. doi: 10.3969/j.issn.0258-2724.2010.06.019
[7]	XIANG Hongyan, ZHANG Lin, YANG Bo. Optimization of Road Network Structure Based on Maximum-Flow Theory[J]. Journal of Southwest Jiaotong University, 2009, 22(2): 284-288.
[8]	LÜ, Xiongwei, LI Jun, LEI Ming, ZHANG Bin. Multi-objective Optimization of Stochastic Demand Inventory Routing Problem with Time Windows[J]. Journal of Southwest Jiaotong University, 2009, 22(2): 289-294.
[9]	LI Li, CUI Dabin, JIN Xuesong. State of Arts of Research on Railway Wheel Profile Optimization[J]. Journal of Southwest Jiaotong University, 2009, 22(1): 13-19.
[10]	CHEN Xiang, XU Bochu, ZHANG Weihua. Optimization of Seat Comfort of High-Speed Train[J]. Journal of Southwest Jiaotong University, 2009, 22(6): 906-911.
[11]	MOU Qifeng, WANG Ciguang. Fast Optimization of Final Approach Sequence of Aircrafts[J]. Journal of Southwest Jiaotong University, 2008, 21(3): 404-408.
[12]	LIU Xin, QIU Yuanying, SHENG Ying, LI Junhua. Optimization of Comprehensive Performance of Planar Redundant Parallel Manipulator[J]. Journal of Southwest Jiaotong University, 2008, 21(5): 626-632.
[13]	LI Zongping, JIANG Sai. Client Grouping in Distribution System Optimization[J]. Journal of Southwest Jiaotong University, 2006, 19(5): 594-598.
[14]	NIChang-jian, CUIPeng, XIANG Rui. Universal Immune Evolutionary Algorithm for Interval-Constrained Optim ization Problem s[J]. Journal of Southwest Jiaotong University, 2005, 18(4): 548-552.
[15]	LIU Ge, LI Bai-lin. New Optimization Model of Curve Fairness[J]. Journal of Southwest Jiaotong University, 2002, 15(5): 584-587.
[16]	GUANQin-chuan. Fuzzy Multi-objective Optimization Based on Neural Networks[J]. Journal of Southwest Jiaotong University, 2002, 15(3): 338-342.

Supplements(0)

Cited By

Cited by
Periodical cited type(0)
Other cited types(1)

Proportional views

Proportional views

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

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

Figures(11) / Tables(2)

Get Citation

PDF

XML

Article views(566) PDF downloads(14)

Optimization and Tensile Properties of Composite Insulator at High and Low Temperature

doi: 10.3969/j.issn.0258-2724.20181032

Abstract

1. 空间主缆几何线形关联迭代法

1.1 基本假定与约定

1.2 外荷载节点力和力矩

1.3 主缆线形方程

1.4 主缆无应力长度

1.4.1 传统方法

1.4.2 几何关联迭代法

1.5 几何关联迭代法的应用

2. 单边空间索面主缆线形计算

2.1 粗算阶段

2.1.1 竖平面线形计算

2.1.2 水平面线形计算

2.1.3 阶段成果

2.2 精算阶段

3. 算　例

3.1 粗算阶段

3.2 精算阶段

3.3 与有限元软件成果对比

4. 结　论

References

Relative Articles

Cited by

Periodical cited type(0)

Other cited types(1)

Proportional views

Catalog

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

Proportional views

Related

Optimization and Tensile Properties of Composite Insulator at High and Low Temperature

doi: 10.3969/j.issn.0258-2724.20181032

Abstract

1. 空间主缆几何线形关联迭代法

1.1 基本假定与约定

1.2 外荷载节点力和力矩

1.3 主缆线形方程

1.4 主缆无应力长度

1.4.1 传统方法

1.4.2 几何关联迭代法

1.5 几何关联迭代法的应用

2. 单边空间索面主缆线形计算

2.1 粗算阶段

2.1.1 竖平面线形计算

2.1.2 水平面线形计算

2.1.3 阶段成果

2.2 精算阶段

3. 算 例

3.1 粗算阶段

3.2 精算阶段

3.3 与有限元软件成果对比

4. 结 论

References

Relative Articles

Cited by

Periodical cited type(0)

Other cited types(1)

Proportional views

Catalog

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

Proportional views

Related

Export File

Citation

Format

Content

3. 算　例

4. 结　论