Mechanical Characteristics of Super-Long-Span Cable-Stayed Bridge with Transverse Asymmetrical Load
-
摘要:
为研究横桥向非对称恒载对大跨度斜拉桥受力行为的影响,以常泰长江大桥为工程背景,对横桥向非对称恒载桥梁结构的力学行为进行分析,并提出了合理的控制措施以实现理想的成桥状态. 提出了一种索力快速优化的方法,采用上下游索力不对称和降低铁路二期恒载的措施解决横桥向不对称恒载引起的主梁扭转、主塔偏位问题,实现了合理的成桥状态. 结果表明:通过上下游索力不对称可以有效地控制主梁扭转,斜拉索顺桥向不平衡分力和横桥向不平衡分力产生的横桥向弯曲变形大小基本一致,方向相反,有利于降低主梁横桥向弯曲变形;通过减小横桥向不对称恒载的差值可以有效地控制主塔横桥向偏位.
Abstract:In order to study the influence of transverse asymmetrical load on the mechanical characteristics of long-span cable-stayed bridges, this paper takes the Changtai Yangtze River Bridge as the engineering background and analyzes the mechanical characteristics of the bridge with a transverse asymmetrical load. In addition, reasonable control measures are put forward to realize the ideal state of the final bridge. Firstly, a fast optimization method of the cable force is proposed for the cable-stayed bridge with a transverse asymmetrical load. Secondly, asymmetrical upstream and downstream cable forces are used to control the main girder torsion caused by the transverse asymmetrical load. The influences of the unbalanced components along the cable-stayed bridge and the transverse direction of the bridge on the alignment of the main girder and the deviation of the main tower are analyzed. Finally, the method of reducing the secondary dead load of the railway is used to control the deviation of the main tower caused by the transverse asymmetrical load. The results show that the main girder torsion can be effectively controlled by the asymmetrical upstream and downstream cable forces. The transverse bending deformations generated by the unbalanced components along the cable-stayed bridge and the transverse direction of the bridge are basically consistent but opposite in direction, which is beneficial to reduce the transverse bending deformations of the main girder. The transverse deviation of the main tower can be effectively controlled by reducing the difference of the transverse asymmetrical load.
-
近年来,振动半主动控制系统在各工程领域受到广泛关注,其能提供接近主动控制装置的减振效果,且能耗低,同时还能提供如同被动装置的可靠性. 在众多半主动控制系统中,磁流变(MR)阻尼器在应用的简单性、可靠性和鲁棒性方面具有显著优势,应用磁流变阻尼器的半主动控制悬架系统受到车辆振动舒适性研究人员的重视.
对磁流变半主动悬架系统的研究内容主要包括磁流变阻尼器动力学建模[1-2]、悬架系统建模及振动传递特性研究[3]、半主动悬架控制系统仿真分析及试验验证[4-5]等. 磁流变阻尼器作为半主动悬架的执行元件,其力学特性分析及正逆模型建立一直是半主动悬架的研究热点问题. MR阻尼器的阻尼力表现出明显的非线性及滞回特性,且具有显著的场强依赖性及激励特性(振幅及频率)依赖性[6]. 阻尼力具有多值特性,不仅取决于阻尼器活塞速度,同时受活塞运动历史状态及控制电流的影响.
用于描述磁流变阻尼力的模型众多,如:包含微分环节的Bouc-Wen模型、Duhem模型、LuGre模型等;包含代数模型的sigmoid模型、双曲正切模型、魔术公式模型等;由非参数模型表示的多项式模型、神经网络模型等[7-9]. 上述模型各自具有准确性、复杂性等方面的优缺点,目前还没有一种通用的应用于半主动悬架的磁流变阻尼器力学模型.
建立准确的控制电流逆模型是使阻尼器产生期望阻尼力的重要基础,传统的逆模型建立方法包括不基于阻尼力正模型的闭环控制逆模型(如阶梯函数逆模型、Sigmoid函数逆模型、状态连续反馈逆模型等)、基于阻尼力正向模型直接逆推的逆动力学模型(如基于Bingham模型、Bouc-Wen模型、魔术公式模型模型等)[10-12]. 基于正向模型直接逆推逆动力学模型的主要难点在于容易产生难以直接计算解析解的超越方程,不可避免地需要对正向模型进行简化处理,影响建模精度. 上述建模方法无法准确描述不同工况下控制电流与阻尼器运动状态及期望阻尼力的关系,因此,基于神经网络及模糊推理系统的逆模型建模方法受到研究者们的关注.
本文在试验获取磁流变阻尼器力学性能的基础上,对磁流变阻尼器正逆动力学模型进行建模及参数识别,并将建立的正逆模型引入半主动座椅悬架控制系统模型中进行分析,对所提出正逆模型的效果进行验证.
1. 磁流变阻尼器建模及参数识别
1.1 磁流变阻尼器滞回分离模型
建立结构简单、拟合精度高且通用的磁流变阻尼器模型是研究磁流变阻尼器力学特性及半主动悬架控制系统的重点和基础. 基于魔术公式的磁流变阻尼器力学模型[13]与Bouc-Wen模型等包含微分方程的模型相比,计算复杂度明显降低,进而能够有效提高模型的参数辨识效率,如式(1)所示.
{Fd=Fk+Fc+Fm+Fb+f,Fk=k(x−x0),Fc=c˙x,Fm=m¨x,Ff=f0sign˙x,Fb=Dsin[Carctan(B˙x−E(B˙x−arctanB˙x))], (1) 式中:
$ {F}_{{\rm{d}}} $ 、$ {F}_{{\rm{k}}} $ 、$ {F}_{{\rm{c}}} $ 、$ {F}_{{\rm{m}}} $ 、$ {F}_{{\rm{f}}} $ 及$ {F}_{{\rm{b}}} $ 分别为模型输出阻尼力、弹性项、黏性项、惯性项、摩擦项及剪切项,该模型滞回特性主要由$ {F}_{{\rm{k}}} $ 及$ {F}_{{\rm{m}}} $ 表征;B、C、D、E为调整$ {F}_{{\rm{b}}} $ 曲线形状的参数,$ B={b}_{1}I + {b}_{2} $ ,$ D={d}_{1}I + {d}_{2} $ ,I为电流,b1、b2、d1、d2为与电流相关的参数;$ x $ 和$ \dot{x} $ 分别为阻尼器活塞位移及速度;$ k $ 、$ c $ 、$ {f}_{0} $ 分别为弹性项、黏性项、摩擦项系数;$ {x}_{0} $ 为弹性项初始位移;f为频率.在建立模型参数与电流的关系中,魔术公式内部参数B包含与电流相关参数,无法直接逆向推导出控制电流与期望力的关系式,对控制电流的计算精度产生影响.
本文在魔术公式模型的基础上,参考于建强等[14-15]建立的基于骨架曲线与滞回分离的建模方法,建立基于正弦、余弦型魔术公式的磁流变阻尼器滞回分离模型,如式(2). 滞回分离模型建立的原理如图1所示.
{Fd=FmaxFnorm +f0,Fnorm =Fpf+Fhyssignx+c˙x,Fpf=Dpsin[Cparctan(Bp˙x−Ep(Bp˙x−arctanBp˙x))],Fhys=Dhcos[Charctan(Bh˙x−Eh(Bh˙x−arctanBh˙x))], (2) 式中:
$ {F}_{{\rm{max}}} $ 和$ {F}_{{\rm{norm}}} $ 分别为最大阻尼力及归一化后的阻尼力;$ {F}_{{\rm{pf}}} $ 和$ {F}_{{\rm{hys}}} $ 分别对应骨架曲线与滞回分离曲线;Bp、Cp、Dp、Ep及Bh、Ch、Dh、Eh分别为影响Fpf及Fhys曲线形状的参数.1.2 滞回模型参数辨识
本文采用德国SCHENCK公司液压作动器及Lord RD-8040-1磁流变阻尼器(图2),试验工况包括振幅5~20 mm、频率0.5~9.0 Hz的正弦位移激励,阻尼器输入电流为0~1 A.
磁流变减振器力学模型参数识别问题本质上是寻找模型计算结果和减振器实际力学特性误差最小值的数学优化问题,如式(3)所示.
Fobj=min{RMS(Fcalcu−Ftest)}, (3) 式中:
$ {F}_{{\rm{obj}}} $ 、$ {F}_{{\rm{test}}} $ 和$ {F}_{{\rm{calcu}}} $ 分别为目标函数、试验值及模型计算值.为考虑磁流变阻尼力动态力不对称性,通过处理试验数据,得到不同工况下滞回环速度中心,其随电流及激励的变化规律如图3所示. 由图3可知:滞回环速度中心随施加电流增加而逐渐增加;随激励幅值及频率增加均呈现出增加趋势,当频率大于2 Hz时,该趋势更加明显.
因此,可建立滞回环速度中心
$ {\dot{x}}_{{\rm{center}}} $ 与电流及激励特性的关系,如式(4)所示.˙xcenter=(v1I+v2)(v3lnvm+v4), (4) 式中:
$ {v}_{{\rm{m}}}\mathrm{为} $ 最大速度,对于正弦激励信号,${v}_{{\rm{m}}}= 2{\text{π}} Af=\sqrt{{\dot{x}}^{2}-x\ddot{x}}$ ,$\ddot{x} $ 为阻尼器加速度,A 为激励幅值;$ {v}_{1} $ 、$ {v}_{2} $ 、$ {v}_{3} $ 及$ {v}_{4} $ 为拟合公式系数.将阻尼力非对称特性考虑进模型中,式(2)中所有速度项
$ \dot{x} $ 变为$ \dot{x}-{\dot{x}}_{{\rm{center}}} $ .图4为不同激励及电流条件下最大阻尼力变化曲线,本文试验所得阻尼力随施加电流增加而逐渐增加,并未达到完全磁饱和状态. 在相同电流条件下,最大阻尼力均随阻尼器最大活塞速度增加而增加,利用式(5)形式函数可以较准确地描述最大阻尼力和电流及激励特性的关系.
Fmax(I,vm)=(k1I2+k2I+k3)(k4vk5m), (5) 式中:
$ {k}_{1} $ ~$ {k}_{5} $ 为拟合公式系数.为对考虑阻尼力非对称性及最大阻尼力与激励特性及电流变化关系后的模型参数进行辨识,本文采用一种新型的混合优化算法——基于Sobol 序列的差分-禁忌混合算法(SS-DTHA),该算法基于差分进化算法较强的全局搜索能力,引入禁忌搜索增强算法局部搜索能力,并通过调节算法对全局搜索和局部搜索的侧重比例,进一步改进算法的探索和挖掘能力,对解决非线性系统参数辨识问题具有积极作用,如图5所示.
图6为本文提出的改进魔术公式模型与tanh模型及双Sigmod模型计算阻尼力的辨识结果对比(激励工况为15 mm-4.0 Hz-1.0 A). 由图可知:本文提出的模型在滞回环过渡区及屈服后区域均具有更好的拟合精度;3种模型与试验数据的误差均方根值分别为41.07、96.83、90.26 N.
图7为不同输入电流及激励频率条件下,模型计算阻尼力的辨识结果,并与试验测试值进行对比,图中实线和虚线分别代表试验值及模型计算值. 由图可知,在不同激励频率及电流条件下,本文提出的阻尼力模型均能准确地表征磁流变阻尼器动态阻尼力的非线性及滞回特性.
表1为不同激励幅值、频率及输入电流条件下平均百分比误差
$ {E}_{{\rm{average}}} $ 的计算结果,如式(6). 可以看出,本文提出的模型在不同工况下均能得到较低的误差,模型精度较高. 各工况平均百分比误差在3.4%附近,误差变化范围为2.7%~4.8%. 模型对试验阻尼力拟合效果较好,满足半主动控制系统建模要求.表 1 不同工况模型平均百分比误差Table 1. Comparison of average percentage error of models under different working conditions电流/A 幅值 5 mm 幅值 15 mm 3 Hz 4 Hz 3 Hz 4 Hz 0 0.0447 0.0355 0.0457 0.0452 0.2 0.0474 0.0355 0.0289 0.0396 0.4 0.0316 0.0321 0.0302 0.0319 0.6 0.0310 0.0322 0.0269 0.0284 0.8 0.0344 0.0349 0.0305 0.0315 1.0 0.0308 0.0359 0.0282 0.0273 Eaverage=RMS(Fcalcu−Ftest)/RMS(Ftest). (6) 为建立适用于阻尼器工作范围的通用模型,通过分析模型各参数的识别结果可以发现:参数Bp、Bh、c随电流及激励特性变化具有明显的规律,其他参数与电流、激励无明显变化关系,设置为常数. 参数Bp、Bh、c随电流及激励特性变化关系如图8所示. 由图可知:参数Bp、Bh随激励频率增加逐渐降低,在相同激励工况下,电流变化对参数Bp、Bh影响不大;参数c随电流变化无明显增降规律,随频率增加而逐渐降低.
上述参数与激励和电流的关系拟合为
{Bp=p1vp2m,Bh=h1vh2m,c=c1vc2m, (7) 式中:
$ {p}_{1} $ 、$ {p}_{2} $ 、$ {h}_{1} $ 、$ {h}_{2} $ 、$ {c}_{1} $ 及$ {c}_{2} $ 分别为各参数的拟合系数.基于以上参数的变化规律,对幅值为5 mm正弦激励下1.0~9.0 Hz及0~1.0 A电流共47个工况模型参数进行重新辨识,可得到不同工况通用模型的辨识结果,如表2所示.
表 2 滞回模型参数辨识结果Table 2. Parameter identification of hysteretic model变量 辨识结果 ${F}_{{\rm{max}}}$ $\left(-9.36{I}^{2} + 123.97I + 23.05\right)\left(4.29{ {v}^{0.19}_{{\rm{m}}} }\right)$ ${\dot{x} }_{{\rm{center}}}$ $\left(-3.47I + 1.91\right)(-2.35{\rm{ln} }\;{v}_{ {\rm{m} } } + 9.21)$ ${B}_{{\rm{p}}}$ $1.97{v}^{-1.09}_{{\rm{m}}}$ ${B}_{{\rm{h}}}$ $1.01{v}_{{\rm{m}}}^{-0.86}$ $ c $ $0.32{ {v}_{ {\rm{m} } } ^{-1.01}}$ Cp 1.32 Dp 0.70 Ep −7.64 Ch 1.10 Dh 0.42 Eh −11.50 f0 −80.34 图9为所建立的通用模型对10 mm正弦激励下1.0~5.0 Hz、0.2~1.0 A电流不同工况的拟合结果与试验数据对比曲线,图中实线代表试验数据,虚线代表模型计算结果. 本文所建通用模型对不同工况下的阻尼器动态力学特性均具有较好的拟合能力.
2. 基于ANFIS建立阻尼器逆模型
半主动悬架控制系统需要根据座椅运动状态及控制策略计算期望阻尼力,需建立控制电流逆模型计算获得期望阻尼力所需的电流. 本文选用自适应神经模糊系统(ANFIS)建立磁流变逆动力学模型.
为构建准确的ANFIS逆模型,建立准确全面的训练数据集是建模的重点和基础. 本文将阻尼器动态力学特性试验中所有工况案例中均匀抽取数据点整合在起,包含激励幅值5、10、15、20 mm,激励频率0.5~9.0 Hz,输入电流0~1.0 A工况,形成16400个数据点的总集合. 为避免阻尼器运动状态突变对模型的影响,选取每组数据开始及结束速度均接近0,并在每组数据后加10个全零数组,以避免先前数据集合对下一集合的输出预测的影响.
为增加ANFIS模型对变电流工况的推理能力,本文在试验获取的测试数据的基础上,根据磁流变减振器工作范围,建立激励幅值范围为2~14 mm、频率范围为0.5~10.0 Hz的扫频激励. 通过对位移信号进行微分来获得速度信号,时间步长为0.01 s,通过30 s模拟生成总共3000个数据集. 输出电流训练数据由幅值 0~1.0 A 的随机信号产生. 在激励数据不变的基础上,生成5组随机电流数组共15 000个数据点.
本文逆模型训练数据在试验恒定电流数据集及通过正向模型计算得到的随机电流数据集中各抽取7500组数据组成15000个数据点的训练集合,其余数据作为验证数据. 选取
$x\left(t\right)、\dot{x}\left(t\right)、\ddot{x}\left(t\right)、F\left(t\right)、 F\left(t-1\right)$ 作为模型输入,建立包含32 个“if-then”模糊规则的五输入单输出ANFIS模型.图10(a)为训练后逆模型的预测电流与验证目标恒定电流对比,由图可知,预测电流在目标电流值附近波动,两者之间误差均方根值为0.0869 A,电流在0附近误差较大. 图10(b)为随机电流序列模型验证结果,预测电流能够跟随目标电流的变化趋势,两者之间误差均方根值为0.1171 A.
将阻尼器逆向动力学模型与阻尼力模型结合,使该系统根据期望阻尼力产生预测电流进而产生预测阻尼力,模型验证流程如图11所示. 预测阻尼力与目标阻尼力对比结果见图12(a). 结果表明,由 ANFIS 逆模型得到的预测电流以及产生的预测阻尼力能够较好地跟随目标阻尼力. 图12(b)为一组在阻尼力活塞运动速度范围内的验证数据,可以发现,在阻尼器运动范围内,预测阻尼力与目标阻尼力变化范围基本一致. 阻尼器全工况范围内预测阻尼力误差均方根值为135.79 N,为阻尼器最大阻尼力的5.6%.
3. 半主动座椅悬架控制系统仿真
建立包含上述建立的阻尼器正逆模型的单自由度半主动座椅悬架系统,簧上质量、刚度及被动阻尼系数分为75 kg、47.37 kN/m及1200 Ns/m. 控制策略采用on-off控制、连续天棚控制及LMS控制算法进行对比分析,on-off控制根据座椅运动状态在电流为最小值(0)和最大值(1.0 A)之间切换,天棚阻尼系数设置为2000 Ns/m,LMS算法步长因子设置为0.005,除控制算法外的半主动座椅悬架模型剩余部分均设置相同. 激励信号选取幅值5 mm、频率1.0~9.0 Hz正弦扫频激励. 图13为半主动座椅悬架控制系统示意.
图14为被动和半主动座椅悬架位移传递率. 可以发现:被动座椅悬架在共振频段将座椅下地板振动激励放大;与被动座椅悬架相比,3种控制算法在低频共振频段均能够有效地衰减振动. 通过计算,在4.0 Hz处3种算法的座椅悬架位移传递率较原被动座椅分别降低了47.62%、48.72%、42.09%. 天棚控制算法及最小均方算法(least mean square, LMS)在共振频率以上仍能保持较好的振动传递特性,on-off控制对高频激励的振动衰减性能较差.
半主动控制系统的实时性对控制效果具有显著影响. 基于磁流变阻尼器的半主动座椅悬架系统各环节时间延迟如图15所示,包括传感器信号采集时延
$ \Delta {t}_{3} $ 、 控制器时延$ \Delta {t}_{2} $ (数据处理、计算期望阻尼力、输出驱动电流时延)及磁流变阻尼器响应时延$ \Delta {t}_{1} $ .本文用经典控制算法(on-off控制、天棚控制)对磁流变阻尼器响应时滞的鲁棒性进行分析,计算结果如图16所示. 可以发现:对于两种算法,增加阻尼器响应时滞均会降低振动控制的改善效果;对于采用连续控制模型的天棚控制算法,对于阻尼器响应时滞具有一定鲁棒性,时延18 ms时在仿真频域范围内仍可以实现优于被动座椅悬架的减振效果;on-off控制算法对阻尼器响应时滞较敏感,特别是超过
$ \sqrt{2} $ 倍固有频率后的高频段,当阻尼器时延18 ms时,on-off控制座椅悬架振动传递特性变差. 因此,尽可能减小半主动座椅悬架系统各个环节的时滞,对于改善半主动控制系统减振效果具有积极作用.通过上述分析可知,本文建立的阻尼力模型及控制电流逆模型能够产生满足控制策略需求的期望阻尼力,有效地衰减座椅下底板振动,对改善车辆乘坐舒适性具有积极作用.
4. 结 论
针对半主动座椅悬架控制系统,分析磁流变阻尼器非线性滞回特性,提出了一种改进魔术公式模型,并用改进差分进化算法进行模型辨识;利用自适应神经模糊系统建立磁流变阻尼器控制电流逆模型,并对包含阻尼器正逆模型的半主动座椅悬架进行验证计算. 结果表明,本文建立的正逆模性能够较好地表征磁流变阻尼器非线性滞回特性,满足座椅悬架的半主动控制需求.
1) Magic Formula 模型作为一种经典的轮胎模型,具有广泛的应用,对于非线性系统具有良好的拟合能力. 本文探索正弦型及余弦型魔术公式在磁流变减振器力学模型中的应用.
2) 在不同激励特性及电流条件下,本文建立的阻尼力模型均能较好地描述磁流变阻尼器的非线性滞回特性,对应于5 mm、3 Hz激励条件下,电流0~1.0 A变化范围内的模型平均百分比误差分别为0.0447、0.0474、0.0316、0.0310、 0.0344及0.0308.
3) 选取
$ x\left(t\right),\dot{x}\left(t\right),\ddot{x}\left(t\right),F\left(t\right),F\left(t-1\right) $ 作为模型输入,建立包含32 个“if-then”模糊规则的五输入单输出ANFIS磁流变逆动力学模型. 本文建立的逆模型对恒定电流工况及随机电流工况均可以得到较好的拟合效果,误差均方根值分别为0.0869 A及0.1171 A.4) 阻尼器逆向动力学模型与阻尼力模型串联模型对期望阻尼力的预测能力较好,阻尼器全工况范围内预测阻尼力误差均方根值与阻尼器最大阻尼力的百分比可达到5.6%.
5) 与被动座椅悬架相比,磁流变半主动座椅悬架在共振频率范围内均能够有效地衰减振动. on-off控制、连续天棚控制及LMS控制算法在4 Hz处的座椅悬架位移传递率较原被动座椅分别降低了47.62%、48.72%、42.09%. 阻尼器响应时滞增加会降低半主动控制的减振效果. 因此,应尽可能减小半主动座椅悬架系统各个环节的时滞.
-
表 1 主塔塔顶变形
Table 1. Deformations of the top of the main tower
mm 不计入收缩徐变效应 计入收缩徐变效应 顺桥向变形 横桥向变形 顺桥向变形 横桥向变形 −104.6 −88.1 −100.9 −125.2 表 2 降低铁路侧恒载措施
Table 2. Measures to reduce secondary dead load of the railway
措施 优化前 优化后 降低铺砟厚度 35 cm 厚普通
道床25 cm 厚聚氨酯固化道床 压缩道床宽度 9 m 8.2 m 调整铁路防水面板结构 6 cm 厚高性能混凝土铺装 3 mm 不锈钢与基材形成复合钢板 斜拉索 不同规格、相同
强度同种规格、不同
强度 -
[1] 秦顺全. 武汉天兴洲公铁两用长江大桥关键技术研究[J]. 工程力学,2008,25(增2): 99-105.QIN Shunquan. Key technology of Tianxingzhou Yangtze River Bridge[J]. Engineering Mechanics, 2008, 25(S2): 99-105. [2] 高宗余. 沪通长江大桥主桥技术特点[J]. 桥梁建设,2014,44(2): 1-5.GAO Zongyu. Technical characteristics of main bridge of Hutong Changjiang River Bridge[J]. Bridge Construction, 2014, 44(2): 1-5. [3] 秦顺全,徐伟,陆勤丰,等. 常泰长江大桥主航道桥总体设计与方案构思[J]. 桥梁建设,2020,50(3): 1-10. doi: 10.3969/j.issn.1003-4722.2020.03.001QIN Shunquan, XU Wei, LU Qinfeng, et al. Overall design and concept development for main navigational channel bridge of Changtai Changjiang River Bridge[J]. Bridge Construction, 2020, 50(3): 1-10. doi: 10.3969/j.issn.1003-4722.2020.03.001 [4] 秦顺全. 分阶段成形结构过程控制的无应力状态控制法[J]. 中国工程科学,2009,11(10): 72-78. doi: 10.3969/j.issn.1009-1742.2009.10.011QIN Shunquan. Unstressed state control method for process control of structure formed by stages[J]. Strategic Study of CAE, 2009, 11(10): 72-78. doi: 10.3969/j.issn.1009-1742.2009.10.011 [5] QIN S Q, WEI K, QIN J X, et al. Stress-free-state based structural analysis and construction control theory for staged construction bridges[J]. Advances in Bridge Engineering, 2020, 1(1): 1-16. doi: 10.1186/s43251-020-00001-y [6] QIN S Q, QIN J X, ZHANG D. Equilibrium equation for spatial frame structure constructed in stages and its applications[C]//IABSE Symposium Report. [S.L.]: International Association for Bridge and Structural Engineering, 2013: 1-8. [7] 但启联,秦顺全,魏凯,等. 基于平面梁单元的分阶段成形结构线形控制方程[J]. 桥梁建设,2017,47(4): 42-47. doi: 10.3969/j.issn.1003-4722.2017.04.008DAN Qilian, QIN Shunquan, WEI Kai, et al. Equation for geometric shape control of structure formed in stages based on plane beam elements[J]. Bridge Construction, 2017, 47(4): 42-47. doi: 10.3969/j.issn.1003-4722.2017.04.008 [8] 许磊平,秦顺全,马润平. 基于平面壳单元的分阶段成形结构平衡方程[J]. 西南交通大学学报,2013,48(5): 857-862. doi: 10.3969/j.issn.0258-2724.2013.05.012XU Leiping, QIN Shunquan, MA Runping. Equilibrium equation derivation of structures formed by stages based on plane shell element[J]. Journal of Southwest Jiaotong University, 2013, 48(5): 857-862. doi: 10.3969/j.issn.0258-2724.2013.05.012 [9] 许磊平,秦顺全,苑仁安. 基于Timoshenko梁的分阶段成形结构平衡方程[J]. 桥梁建设,2015,45(5): 24-29.XU Leiping, QIN Shunquan, YUAN Renan. Equilibrium equation for structure formed in stages based on Timoshenko beam[J]. Bridge Construction, 2015, 45(5): 24-29. [10] 苑仁安,秦顺全,王帆. 分阶段成形杆系结构几何非线性平衡方程[J]. 桥梁建设,2014,44(2): 50-55.YUAN Renan, QIN Shunquan, WANG Fan. Equilibrium equation for geometric nonlinearity of frame structure formed in stages[J]. Bridge Construction, 2014, 44(2): 50-55. [11] 苑仁安,秦顺全,王帆,等. 基于平面梁单元的几何非线性分阶段成形平衡方程[J]. 桥梁建设,2014,44(4): 45-49.YUAN Ren’an, QIN Shunquan, WANG Fan, et al. Equilibrium equation for geometric nonlinearity of structure formed in stages based on plane beam elements[J]. Bridge Construction, 2014, 44(4): 45-49. [12] 苑仁安,秦顺全,肖海珠. 一种斜拉桥目标状态索力快速精准确定的方法[J]. 桥梁建设,2020,50(2): 25-30. doi: 10.3969/j.issn.1003-4722.2020.02.005YUAN Renan, QIN Shunquan, XIAO Haizhu. A method to rapidly and accurately determine target cable forces for cable-stayed bridge[J]. Bridge Construction, 2020, 50(2): 25-30. doi: 10.3969/j.issn.1003-4722.2020.02.005 [13] 张清华,黄灿,卜一之,等. 大跨度钢斜拉桥制造误差的传播及其效应特性[J]. 西南交通大学学报,2015,50(5): 830-837. doi: 10.3969/j.issn.0258-2724.2015.05.010ZHANG Qinghua, HUANG Can, BU Yizhi, et al. Fabrication error propagation properties of key components of large-span cable-stayed bridges with steel box girder[J]. Journal of Southwest Jiaotong University, 2015, 50(5): 830-837. doi: 10.3969/j.issn.0258-2724.2015.05.010 [14] 谢明志,杨永清,卜一之,等. 千米级混合梁斜拉桥双目标控制施工监控体系[J]. 西南交通大学学报,2018,53(2): 244-252,321. doi: 10.3969/j.issn.0258-2724.2018.02.004XIE Mingzhi, YANG Yongqing, BU Yizhi, et al. Construction control system for thousand-meter-scale hybrid girder cable-stayed bridge based on double target control[J]. Journal of Southwest Jiaotong University, 2018, 53(2): 244-252,321. doi: 10.3969/j.issn.0258-2724.2018.02.004 -