水平加筋体加固岩溶路基受力机理及设计方法

张东卿; 薛元; 罗强; 刘菀茹; 郑永飞

doi:10.3969/j.issn.0258-2724.20170784

摘要: 为了解决现有设计方法未考虑水平加筋体的滑移效应，低估铁路路基面变形的问题，对水平加筋体加固铁路岩溶路基的受力机理进行了研究，提出了一种新的设计方法. 首先，明确了荷载传递机制和水平加筋体受力模型；然后，从严格控制铁路路基面变形的角度出发，在现有设计方法的基础上，考虑滑移效应并假定锚固区加筋体应变呈线性变化，通过引入加筋体松弛量建立了加筋体挠度修正公式；再次，提出了考虑滑移效应的水平加筋体加固铁路岩溶路基的设计方法，明确了验算内容和设计流程；最后，将所提方法应用于某模型试验，并与实测值及数值分析值进行了对比. 研究结果表明：考虑滑移效应能够有效修正路基面变形计算值，路基面变形值由0.18 m修正为0.25 m，接近于实测值0.23 m和数值分析值0.21 m，方法准确，可为实际设计提供依据.

关键词:

Abstract: In order to solve the problem that the existing design methods do not consider the frictional effect of the horizontal reinforcement, and the deformation of the railway surface was underestimated, the load mechanism for railway karst subgrade reinforced by horizontal geosynthetic reinforcement was studied and a new design was proposed. Firstly, the mechanism of load transfer and mechanical model were indicated. And then, from the perspective of strictly controlling the deformation of railway subgrade surface, on the basis of the existing design methods, taking into account the frictional effect, assuming that the strain in anchorage zone changes linearly, the correction formula of reinforcement deflection was established by introducing the relaxation length. Furthermore, the design method considering the frictional effect for railway karst subgrade reinforced by horizontal reinforcement was put forward; in addition, the verification content and design process were also proposed. Finally, the method proposed was applied to a model test, and the result was compared with the results of model test and numerical analysis. The research results show that: considering the frictional effect can effectively modifier the calculated deformation of subgrade surface, which can increase from 0.18 m to 0.25 m, and was close to the model test value of 0.23 m and the numerical analysis value of 0.21 m; the method proposed is reasonable and can provide reference for actual design.

Key words:

高速磁浮列车速度快、安全舒适、维保成本低，可实现大城市之间和综合经济区之间的快速通勤^[1]. 高速磁浮通过电磁力约束实现无接触运行，摆脱了传统轨道交通的弓网关系与轮轨关系，有潜力成为陆地交通系统中速度最快的运输方式^[2-3].

常导高速磁浮列车采用电励磁长定子直线同步电机（LSM）实现推力输出^[3]. 磁浮列车高速时，采用基于反电动势的无传感器控制策略，实现电流与速度闭环控制^[5]；零速与低速时，采用齿槽检测与标志板结合的定位测速系统检测磁极相角（PRW）^[1]，此时动子磁链的变化率小，且反馈电流信噪比低，通过反电动势难以估计电角度信息；当零速与低速时定位测速传感器失效，列车丢失电角度无法实现电流与速度闭环. 因此，研究磁浮列车零速与低速无传感器（以下统称零低速无传感）控制方法具有一定的实用意义.

传统的零低速无传感技术利用动子凸极性，通过增加外部激励信号跟踪动子位置，即脉冲宽度调制（pulse width modulation, PWM）激励法和高频信号注入法. 注入信号类型多样，有正弦信号、脉振信号以及随机信号（变频率、变幅值、变相位），按注入轴系主要有高频旋转电压和高频脉振电压注入^[6]. 旋转电压法注入到α、β轴，通过提取高频响应电流信号负序分量，再采用跟踪观测器、同步轴系滤波器等方法解调^[7]，得到动子位置信息. 脉振电压法分为脉振正弦和脉振方波2类信号^[8]，信号从d轴注入，通过解调高频响应信号得到基波分量和高频分量，分别实现闭环控制和角度信息估算. 脉振正弦信号分离需要多个滤波器，应用受限于系统带宽；脉振方波注入频率可接近开关频率，且信号分离可以忽略滤波器的影响. 然而，面向高速磁浮与轨道交通等大功率电传动系统，传统高频注入方法往往忽略控制与采样电路延时问题，尤其是在低开关频率下^[9]，若不对延时引起的相移补偿^[10-11]，位置估计累积误差将导致闭环控制失步或错步等故障.

针对上述基于延时的相移滞后问题，国内外学者展开了大量研究工作. 如：针对大功率电机驱动系统低采样、低开关频率导致的延时及交叉耦合增大和系统稳定性与动态性降低，采用复矢量设计电流控制器^[12]；考虑控制延时与PWM生成延时导致的比例积分（PI）控制器参数失配，通过重新设计参数区间，提升驱动系统的相位裕度和带宽^[13]；考虑开关管开关延时，采用在dq控制回路中增加前馈扰动观测器来实现补偿，提高位置估计精度^[13,15]. 为补偿位置滞后，通过计算q轴电流误差来设计位置估计算法^[16-17]，通过欧拉预测计算q轴电压误差来估计开关管工作时间^[18]；为抑制控制延时对磁浮列车电流控制性能的影响，一种预测电流控制器被提出^[19-20]，为逆变器提供弱化延时影响的开关信号；为抑制低开关频率下控制延时对轨道交通列车无传感控制精度的影响，一种基于梯度下降法的高频注入延时补偿方法被提出^[21]. 然而，针对磁浮列车零低速无传感控制，考虑延时导致的相移补偿研究较少.

本文针对磁浮列车零低速高频信号注入无传感控制忽略大功率驱动系统控制延时和采样电路延时的影响，导致高频响应电流电角度相位滞后与速度跟踪误差累积的问题，提出考虑相位滞后补偿的高频信号注入无传感器控制方法. 通过搭建短行程磁浮电机试验平台，验证所提控制方法的正确性.

1. 长定子同步直线电机高频模型

1.1 高频基波模型

常导高速磁浮列车运行推力和悬浮力都与悬浮磁场相关，该磁场大小由悬浮间隙与励磁电流决定，悬浮电磁铁既是悬浮力输出来源，也是LSM的动子部分. LSM定子部分由长1 032 mm的定子模块与三相电枢波绕组构成；动子模块由10个整磁极和2个半磁极总长3 096 mm电磁铁构成，其中头尾车厢每侧包含7.5个标准悬浮电磁铁模块，中间车厢为8个，如图1所示.

图 1 常导高速磁浮直线长定子电机结构

Figure 1. Structure of EMS high-speed linear LSM

下载: 全尺寸图片幻灯片

列车零低速采用高频信号注入无传感器控制方式，高频电流响应信号中包含有动子位置信息，由于高频段的电流信号幅值较小，频率高，则高频段的阻抗压降与速度反电势项可忽略不计，电压方程中仅考虑电流导数项，此时电机可简化为纯电感，电机高频数学模型可表示为

$\left[ {\begin{array}{*{20}{c}} {{u_{d{\text{h}}}}} \\ {{u_{q{\text{h}}}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{L_{d{\text{h}}}}}&0 \\ 0&{{L_{q{\text{h}}}}} \end{array}} \right]\frac{{\text{d}}}{{{\text{d}}t}}\left[ {\begin{array}{*{20}{c}} {{i_{d{\text{h}}}}} \\ {{i_{q{\text{h}}}}} \end{array}} \right],$

(1)

式中：u_dh和u_qh分别为d、q轴的高频电压分量，i_dh和i_qh分别为d、q轴高频电流分量，L_dh和L_qh分别为d、q轴的高频电感分量，t为运行时间.

高频电压信号注入在观测轴系，即d_e、q_e轴系，d、q为实际轴系，含延时相移. 区别于以往的高频电压注入角度估计方法，增加不含延时相移的实际轴系d_r、q_r与估计轴系d_er、q_er，用于补偿角度参考计算，如所示. 图中： $\alpha$ 、 $\beta$ 为静止轴系； ${\theta _{\text{e}}}$ 为含延时的实际电角度； ${\hat \theta _{\text{e}}}$ 为未补偿的估计电角度； ${\theta _{\text{r}}}$ 为不含延时的实际电角度； ${\hat \theta _{{\text{er}}}}$ 为补偿后的估计电角度； $\Delta \theta$ 为补偿前 ${\theta _{\text{e}}}$ 与 ${\hat \theta _{\text{e}}}$ 误差； $\Delta {\theta _{\text{r}}}$ 为 ${\theta _{\text{r}}}$ 与 ${\hat \theta _{\text{e}}}$ 误差； $\Delta {\theta _{{\text{cmp}}}}$ 为补偿后的实际角度（ ${\theta _{\text{e}}}$ + $\hat \varphi$ ）与估计角度 ${\hat \theta _{\text{e}}}$ 误差，即 $\Delta \theta$ + $\hat \varphi$ ； $\varphi$ 为实际延时电角度； $\hat \varphi$ 为估计补偿延时电角度； $\Delta \varphi$ 为 $\varphi$ 与 $\hat \varphi$ 误差，当 $\Delta \varphi$ = 0时，补偿后的实际电角度等于 ${\theta _{\text{r}}}$ .

图 2 考虑相移滞后影响的高频响应各坐标系

Figure 2. Coordinate system of high-frequency response considering effect of phase shift lag

下载: 全尺寸图片幻灯片

1.2 高频方波注入模型与响应电流模型

高频方波电压信号注入在d_eq_e轴系，注入信号的表达式为

$u_{d{\text{h}}}^{\text{e}} = \left\{ \begin{array}{*{20}{l}} {U_{\text{h}}}, & t \in ({T_0},{T_0} + {T_{\text{h}}}/2] \\ - {U_{\text{h}}}, & t \in ({T_0} + {T_{\text{h}}}/2,{T_0} + {T_{\text{h}}}] \\ \end{array} \right.$

(2)

$u_{q{\text{h}}}^{\text{e}} = 0$

(3)

式中： $u_{d{\text{h}}}^{\text{e}}$ 、 $u_{q{\text{h}}}^{\text{e}}$ 分别为观测轴系高频电压分量，U_h为注入的电压幅值，T₀为周期的初始时刻，T_h为高频方波电压周期.

基于图2所示轴系关系，d_eq_e轴系经Park变换可得到dq轴系电压方程，即式（1）左侧可写为

$\left[ {\begin{array}{*{20}{c}} {{u_{d{\text{h}}}}} \\ {{u_{q{\text{h}}}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\cos\; \Delta \theta }&{\sin\; \Delta \theta } \\ { - \sin\; \Delta \theta }&{\cos\; \Delta \theta } \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {u_{d{\text{h}}}^{\text{e}}} \\ {u_{q{\text{h}}}^{\text{e}}} \end{array}} \right].$

(4)

考虑Clarke变换不需反馈电角度，因此，在αβ轴系下对高频电流信号中的位置信息进行解调角度信息，建立起dq轴系与αβ轴系电流方程的关系，具体可表示为

$\left[ \begin{gathered} {i_{d{\text{h}}}} \\ {i_{q{\text{h}}}} \\ \end{gathered} \right] = \left[ {\begin{array}{*{20}{c}} {\cos\; {\theta _{\text{e}}}}&{\sin\; {\theta _{\text{e}}}} \\ { - \sin\; {\theta _{\text{e}}}}&{\cos\; {\theta _{\text{e}}}} \end{array}} \right]\left[ \begin{gathered} {i_{\alpha {\text{h}}}} \\ {i_{\beta {\text{h}}}} \\ \end{gathered} \right],$

(5)

式中：i_dh和i_qh分别为d、q轴的高频电流分量， ${i_{\alpha {\text{h}}}}$ 、 ${i_{\beta {\text{h}}}}$ 分别为 $\alpha$ 、 $\beta$ 轴的高频电流分量.

将式（4）、（5）代入式（1），考虑到电流导数项远大于速度项与反电动势项，因此，仅保留导数项，转换为电流状态方程，得到

$\begin{split} & \frac{{\text{d}}}{{{\text{d}}t}}\left[ \begin{gathered} {i_{\alpha {\text{h}}}} \\ {i_{\beta {\text{h}}}} \\ \end{gathered} \right] = \left[ {\begin{array}{*{20}{c}} {\cos\; {\theta _{\text{e}}}}&{ - \sin\; {\theta _{\text{e}}}} \\ {\sin\; {\theta _{\text{e}}}}&{\cos\; {\theta _{\text{e}}}} \end{array}} \right]{\left[ {\begin{array}{*{20}{c}} {{L_{d{\text{h}}}}}&0 \\ 0&{{L_{q{\text{h}}}}} \end{array}} \right]^{ - 1}}\left[ {\begin{array}{*{20}{c}} {\cos\; \Delta \theta }&{\sin \;\Delta \theta } \\ { - \sin\; \Delta \theta }&{\cos\; \Delta \theta } \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {u_{{\text{dh}}}^{\text{e}}} \\ {u_{{\text{qh}}}^{\text{e}}} \end{array}} \right] \\ &\quad = \frac{{u_{{\text{dh}}}^{\text{e}}}}{{L_{{\text{avg\_h}}}^2 - L_{{\text{dif\_h}}}^2}}\left[ \begin{gathered} {L_{{\text{avg\_h}}}}\cos\; {{\hat \theta }_{\text{e}}} - {L_{{\text{dif\_h}}}}\cos ({\theta _{\text{e}}} + \Delta \theta ) \\ {L_{{\text{avg\_h}}}}\sin\; {{\hat \theta }_{\text{e}}} - {L_{{\text{dif\_h}}}}\sin ({\theta _{\text{e}}} + \Delta \theta ) \\ \end{gathered} \right] = \frac{{u_{{\text{dh}}}^{\text{e}}}}{{{L_{{\text{dh}}}}{L_{{\text{qh}}}}}} {\boldsymbol{\varLambda}} \end{split}$

(6)

式中： ${\boldsymbol{\varLambda}} {\text{ = }}\left[ {\begin{array}{*{20}{c}} {\cos \;{\theta _e}\cos \;\Delta \theta /{L_{{\text{dh}}}} + \sin\; {\theta _e}\sin\; \Delta \theta /{L_{q{\text{h}}}}} \\ {\sin \;{\theta _e}\cos \;\Delta \theta /{L_{{\text{dh}}}} - \cos \;{\theta _e}\sin\; \Delta \theta /{L_{q{\text{h}}}}} \end{array}} \right]$ ，L_dh、L_qh分别为d、q轴的电感常数， $L_{{\text{avg\_h}}}^{}$ 与 $L_{{\text{dif\_h}}}^{}$ 分别为高频均值电感与高频偏差电感.

对式（6）左右积分，得到

$\left[ \begin{gathered} {i_{\alpha {\text{h}}}} \\ {i_{\beta {\text{h}}}} \\ \end{gathered} \right] = \varLambda \int {u_{{\text{dh}}}^e{\mathrm{d}}t}.$

(7)

根据傅里叶变换分解，注入高频方波信号为

$u_{{\text{dh}}}^{\text{e}}(t) = \frac{{4{V_{{\mathrm{in}}}}}}{{\text{π}} }\sum\limits_{k = 0}^\infty {\frac{{\sin [(2k + 1){\omega _{{\text{in}}}}t]}}{{2k + 1}}} \text{，}$

(8)

式中： ${\omega _{{\text{in}}}}$ 为注入频率.

将（8）代入（7）可得

$\left[ \begin{gathered} {i_{\alpha {\text{h}}}} \\ {i_{\beta {\text{h}}}} \\ \end{gathered} \right] = {\boldsymbol{\varLambda}} \frac{{ - 4{V_{{\text{in}}}}}}{{{\text{π}} {\omega _{{\text{in}}}}}}\sum\limits_{k = 0}^\infty {\frac{{\cos (2k + 1){\omega _{{\text{in}}}}t}}{{2k + 1}}} \text{，}$

(9)

式中： ${V_{{\text{in}}}}$ 为注入高频方波信号幅值.

从式（9）可以看出，向电流调节器输出的d轴注入高频电压，在α、β轴高频响应电流中含有动子位置信息，通过解调α、β轴响应电流的估计误差，经锁相环可实现角度位置估计. 由于Λ中含延时的电角度 ${\theta _{\text{e}}}$ 与 $\Delta \theta$ 会导致估计角度 ${\hat \theta _{\text{e}}}$ 相移，因此，需要对估计角度补偿.

2. 延时对角度估计相移偏差影响分析

2.1 延时分析与无传感控制结构

在高速磁浮交通大功率电传动控制系统中，受到开关损耗与散热的限制，高频电压注入信号的频率受到了限制，同时位置控制系统的固有延时，严重影响着位置估计精度与动态性能. 高频注入无传感控制系统中的延时问题分为控制延时与采样电路延时，控制延时会导致PWM输出电压相移与幅值误差；采样电路延时会引起注入电压偏差，导致高频响应电流劣化与相移误差，且电流控制回路动态性能下降，因此，需要对由延时引起的估计角度偏差补偿，保证无位置传感器控制系统的动态性能.

当采用单载波周期双采样实现PWM输出调制，前半周期和后半周期的脉冲分别由两个采样值计算得出. 首先，控制系统在时刻n采样相电流，通过Clarke变换转至αβ轴；其次，经信号解调与误差锁相得到电角度，结合Park将电流转换至dq轴；最后，经d、q轴电流调节器与Clarke逆变换得到静止坐标系下的α、β轴参考电压，计算出新的PWM调制信号用于时刻n + 1的调制输出，此过程存在一个采样周期T_s/2的控制延时. 考虑PWM调制占空比信号至逆变器开关管输出延时约半个采样周期T_s/4，则从总延时约3T_s/4时长.

如图3所示，根据时刻t_n的高频响应电流估计出电角度θ_e，由于控制延时，在时刻t_n + t_d将θ_e反馈闭环计算的PWM调制输出，此时动子已移动至θ_r，高频响应电流和估计的动子位置都有滞后相移，根据式（8）保留 ${\omega _{{\text{in}}}}$ ，补偿一般取基波分量，即k=0，考虑延时的高频响应电流为

图 3 考虑延时的相移滞后示意

Figure 3. Phase shift lag considering delay

下载: 全尺寸图片幻灯片

$\begin{split} & \left[ \begin{gathered} {i_{\alpha {\text{h\_d}}}} \\ {i_{\beta {\text{h\_d}}}} \\ \end{gathered} \right] = \left[ {\begin{array}{*{20}{c}} {\dfrac{{\cos \;{\theta _{\mathrm{e}}}\cos \;\Delta \theta }}{{{\omega _{{\mathrm{in}}}}{L_{{\text{dh}}}}}} + \dfrac{{\sin\; {\theta _{\mathrm{e}}}\sin \;\Delta \theta }}{{{\omega _{{\mathrm{in}}}}{L_{{\text{qh}}}}}}} \\ {\dfrac{{\sin \;{\theta _{\mathrm{e}}}\cos\; \Delta \theta }}{{{\omega _{{\mathrm{in}}}}{L_{{\text{dh}}}}}} - \dfrac{{\cos\; {\theta _{\mathrm{e}}}\sin \;\Delta \theta }}{{{\omega _{{\mathrm{in}}}}{L_{{\text{qh}}}}}}} \end{array}} \right] \\ &\quad\frac{{4{V_{{\mathrm{in}}}}}}{{\text{π}} } \cdot \sin ({\omega _{{\mathrm{in}}}}t - \varphi ), \end{split}$

(10)

式中： $\varphi = {\omega _e}{t_{\mathrm{d}}}$ ，为延时滞后角度.

考虑相位补偿的高频方波注入无传感控制框图如所示. 图中： ${v_{{\mathrm{ref}}}}$ 为给定的控制速度，即设置的目标速度； ${\hat v_{\mathrm{e}}}$ 为观测速度，由算法计算出的观测角频率 ${\hat \omega _{\mathrm{e}}}$ 乘以系统运行时间t得到的； ${\hat \omega _{\mathrm{e}}}$ 为观测角频率，由估计电角度 ${\hat \theta _{\text{e}}}$ 求导得到的； $i_{\text{q}}^*$ 为经过速度控制器ASR（automatic speed regulator）计算后得到的q轴控制电流； $i_{\mathrm{d}}^*$ 为d轴控制电流，是一个常值0； $u_{\mathrm{d}}^*$ 、 $u_{\text{q}}^*$ 分别为经过电流控制器ACR （automatic current regulator）计算后，得到的dq轴控制电压； $u_\alpha ^*$ 、 $u_\beta ^*$ 为 $u_{\mathrm{d}}^*$ 、 $u_{\text{q}}^*$ 经过Clarke反变换计算得到的αβ轴调制控制电压. Clarke反变换需要用到角度 ${\hat \theta _{{\text{er}}}}$ ； ${i_{\text{a}}}$ 、 ${i_{\text{b}}}$ 与 ${i_{\text{c}}}$ 为ABC三相电流； ${i_{\alpha \beta }}$ 为αβ轴反馈电流，包含 ${i_\alpha }$ 与 ${i_\beta }$ ； ${i_{{\text{dq}}}}$ 为dq轴反馈电流，包含 ${i_{\text{d}}}$ 与 ${i_{\text{q}}}$ ； ${i_{\alpha \beta }}$ 与 ${i_{{\text{dq}}}}$ 中都包含基频部分与高频部分； ${i_{\alpha \beta {\text{h}}}}$ 为αβ轴高频部分， ${i_{\alpha \beta {\text{h}}}}$ 分为 ${i_{\alpha {\text{h}}}}$ 与 ${i_{\beta {\text{h}}}}$ ，用于角度估计计算； ${i_{{\text{dqf}}}}$ 表示dq轴基频部分， ${i_{{\text{dqf}}}}$ 分为 $i_{{\text{df}}}^{}$ 、 $i_{{\text{qf}}}^{}$ 为dq轴基波电流，用于电流闭环控制. 电角度估算过程主要包括高频信号解调、包络处理与锁相、相移计算、相移补偿等4个阶段，其中，相移计算与相移补偿分别实现了相移角的计算与补偿实施.

图 4 考虑相位补偿的高频方波注入无传感控制框图

Figure 4. Sensorless control for HFSI considering phase shift lag compensation

下载: 全尺寸图片幻灯片

由式（10）可以看出，若能确定延时时间t_d，通过角度计算可消除延时引起的相移滞后. 为消除估计角度变量，保留滞后角度，对式（10）变换处理分析如下：1）式（10）包含 ${\theta _e}$ 、 $\varphi$ 、t与 ${i_{\alpha {\text{h}}}}$ 、 ${i_{\beta {\text{h}}}}$ 等变量；2） ${\omega _{{\mathrm{in}}}}$ 与 ${V_{{\mathrm{in}}}}$ 为高频注入信号， $\Delta \theta$ 经锁相环输出确定，即 ${\hat \theta _{\mathrm{e}}}$ 确定；3）若t已知，对 $\varphi$ 估计补偿，则应先建立只含 $\varphi$ 的成本函数，再最小化误差目标函数，最终得到延时t_d. 变换处理： ${i_{{\mathrm{td}}}}$ 为电流成本函数，通过左右两侧乘以含观测电角度 ${\hat \theta _{\mathrm{e}}}$ 的矩阵，实现 ${\theta _{\mathrm{e}}}$ 消除，得到 ${i_{{\mathrm{td}}}}$ 幅值与2 $\Delta \theta$ 相关， ${i_{{\mathrm{td}}}}$ 相移与 $\varphi$ 相关，即

$\begin{split} & {i_{{\mathrm{td}}}} = \left[ {\begin{array}{*{20}{c}} {\sin \;{{\hat \theta }_{\text{e}}}}&{ - \cos\; {{\hat \theta }_{\text{e}}}} \end{array}} \right]\left[ \begin{gathered} {i_{\alpha {\text{h}}}} \\ {i_{\beta {\text{h}}}} \\ \end{gathered} \right] \\ &\quad = \left(\frac{{ - \cos\; \Delta \theta \sin\; \Delta \theta }}{{{\omega _{{\text{in}}}}{L_{{\text{dh}}}}}} + \frac{{\sin\; \Delta \theta \cos\; \Delta \theta }}{{{\omega _{{\text{in}}}}{L_{{\text{qh}}}}}}\right) \\ &\quad \frac{{4{V_{{\text{in}}}}}}{{\text{π}} } \cdot \sin ({\omega _{{\text{in}}}}t - \varphi ). \end{split}$

(11)

${i_{{\text{td}}}}$ 的估计值

${\hat i_{{\mathrm{td}}}} = \frac{{({L_{{\text{dh}}}} - {L_{{\text{qh}}}})/2}}{{{\omega _{{\mathrm{in}}}}{L_{{\text{dh}}}}{L_{{\text{qh}}}}}}\frac{{4{V_{in}}}}{{\text{π}} }\sin 2(\Delta {\theta _{{\mathrm{cmp}}}}) \sin ({\omega _{{\mathrm{in}}}}t - \hat \varphi ).$

(12)

根据延时误差机理，不考虑延时， $\Delta {\theta _{}}{\text{ = }}{\theta _{\text{e}}} - {\hat \theta _{\text{e}}}$ ， ${\theta _{\text{e}}}$ 等于 ${\theta _{\text{r}}} - \varphi$ ；考虑延时有

$\Delta {\theta _{\text{r}}}{\text{ = }}{\theta _{\text{r}}} - {\hat \theta _{\text{e}}} = {\theta _{\text{e}}} + \varphi - {\hat \theta _{\text{e}}} = \Delta \theta + \varphi.$

(13)

为使 $\Delta {\theta _{{\text{cmp}}}}$ 逼近 $\Delta {\theta _{\text{r}}}$ ，考虑延时，补偿后角度误差表示为 $\Delta {\theta _{{\text{cmp}}}} = \Delta {\theta _{}} + \hat \varphi$ . 因此， ${\hat i_{{\text{td}}}}$ 表达式为

${\hat i_{{\text{td}}}} = \frac{{({L_{d{\text{h}}}} - {L_{q{\text{h}}}})/2}}{{{\omega _{{\text{in}}}}{L_{d{\text{h}}}}{L_{q{\text{h}}}}}} \cdot \sin [2(\Delta \theta + \hat \varphi )] \cdot \frac{{4{V_{{\text{in}}}}}}{{\text{π}}} \cdot \sin ({\omega _{{\text{in}}}}t - \hat \varphi ).$

(14)

电流误差定义为

$\begin{split} & \Delta {i_{{\text{td}}}} = {i_{{\text{td}}}} - {\hat i_{{\text{td}}}}{\text{ = }}{i_{{\text{td}}}}- \frac{{({L_{d{\text{h}}}} - {L_{q{\text{h}}}})/2}}{{{\omega _{{\text{in}}}}{L_{d{\text{h}}}}{L_{q{\text{h}}}}}} \times \\ &\quad \sin [2(\Delta \theta + \hat \varphi )] \cdot \frac{{4{V_{{\text{in}}}}}}{{\text{π}} } \cdot \sin ({\omega _{{\text{in}}}}t - \hat \varphi ) . \end{split}$

(15)

式（15）对时刻n进行离散化处理得

$\begin{split} & \Delta {i_{{\mathrm{td}}}}[n]{\text{ = }}{i_{{\text{td}}}}[n] -\frac{{({L_{d{\text{h}}}} - {L_{q{\text{h}}}})/2}}{{{\omega _{{\text{in}}}}{L_{d{\text{h}}}}{L_{q{\text{h}}}}}} \\ &\quad \sin [2(\Delta \theta + \hat \varphi )] \frac{{4{V_{{\text{in}}}}}}{{\text{π}} } \sin [n{\omega _{{\text{in}}}}{T_{\text{s}}} - {\omega _{{\text{in}}}}{t_{\text{c}}}(n)] . \end{split}$

(16)

由式（15）误差函数，设计目标函数为

$L(n,{t_{\text{c}}}(n)) = \frac{1}{2}\Delta i_{{\text{td}}}^2.$

(17)

则提出的二分法延时计算结构如图5. 图中：t_mid为二分法中间值.

图 5 考虑梯度变化的二分法延时计算结构

Figure 5. Bisection delay computational structure considering gradient change

下载: 全尺寸图片幻灯片

2.2 考虑梯度变化的二分法寻优算法

准确估计t_d等于找到使L值最小的t_c值，二分法被引入解决最小值问题，是一种高效的逐次逼近算法，这里将t_c当作索引，将 $L(n,{t_{\text{c}}}(n))$ 当作评价指标，即索引对应的数值. 以[a，b]=[T_s/2，3T_s/2]为边界，待求参数t_c的初始值设置为T_s/2以接近实际值t_d，减少迭代总次数.

考虑到 $L(n,{t_{\mathrm{c}}}(n))$ 非单调变化，采用经典二分法无法实现最优值搜索，因此，在增加二分法梯度判断条件的基础上，对搜索路径优化，计算流程如所示. 设置判断误差 $\varepsilon {\text{ = }}{10^{ - 3}}$ ，计算最优补偿指标 $L({t_{{\text{mid}}}})$ ， ${t_{{\text{mid}}}} = (a + b)/2$ . 当 $L({t_{{\text{mid}}}}) > \varepsilon$ 时：

图 6 基于二分法的延时与角度相移滞后计算流程

Figure 6. Flowchart of delay and angular phase shift lag based on bisection method

下载: 全尺寸图片幻灯片

1）若目标函数在t_mid的梯度值p_L（t_mid）<0，则梯度计算后选定的二分边界t_better=b，[a, b]=[t_mid, t_better]；

2）若p_L（t_mid）>0，则t_better=a，[a, b]=[t_better, t_mid]；

3）若p_L（t_mid）=0，p_L（a）<0且p_L（b）<0，则t_better=min{a, b}，[a, b]=[t_better, t_mid]（t_better< t_mid）或[a, b]=[t_mid, t_better]（t_better>t_mid）；

（4）若p_L（t_mid）=0，p_L（a）>0且p_L（b）>0，则完成迭代t_c=t_mid，即当 $L({t_{{\text{mid}}}}) > \varepsilon$ 时，重复上述迭代计算过程.

当 $L({t_{{\text{mid}}}}) < \varepsilon$ 时，完成迭代t_c=（a + b）/2，然后再增加一步迭代，确认最优补偿时间.

3. 试验验证与分析

3.1 试验验证

在实验室中搭建了高速磁浮长定子直线同步电机低速试验平台，如图7所示. 实际磁浮列车中，动子布置在长定子下方，动子励磁后与长定子铁芯产生吸引力克服列车重力实现悬浮. 为简化电机设计与安装难度，采用动子在上定子在下的布置方法. 其中：动子由1个整磁极与2个半磁极构成；受限于变流器容量，定子绕组采用集中式绕法；中性点采用星型连接方式；实际测速采用磁栅传感器. 驱动器试验平台，包括上位机、示波器、实时数字控制器、直流电源、三电平电压型逆变器，实时数字控制器采用RTU-Box控制器，试验逆变器开关频率2 k，注入方波频率1 k，注入电压20 V，试验平台关键参数如表1所示.

图 7 试验平台与驱动控制平台

Figure 7. Test platform and drive control platform

下载: 全尺寸图片幻灯片

表 1 磁浮电机试验平台主要参数

Table 1. Main parameters of maglev motor test platform

参数	数值
直流侧电压/V	220
定子相电阻/Ω	0.12
d轴电感/mH	1.8
q轴电感/mH	1.4
定子极距/mm	2.58
动子极距/mm	266.5
励磁电流/A	20～23
动子励磁磁链/Wb	0.324 7

下载: 导出CSV

| 显示表格

3.2 电流环试验结果

补偿前、后，当q轴控制电流设置为20、21、22 A时，实际电角度 ${\theta _{\text{e}}}$ 、估计电角度 ${\hat \theta _{\text{e}}}$ 、角度误差 $\Delta {\theta _{\text{e}}}$ 与q轴电流误差 $\Delta {i_{\text{q}}}$ 的变化分别如图8、9所示. 从图8可以看出，角度误差分别在−0.52～0.86 rad、−0.41～0.71 rad、−0.32～0.68 rad；从图9可以看出，对应的估计角度误差分别降至了−0.22～0.23 rad、−0.20～0.21 rad、−0.18～0.19 rad；根据两图电流误差显示，由于长定子定位力，反馈电流有一定的周期波动，补偿前后电流波动强度几乎保持一致.

图 8 补偿前不同给定电流控制效果

Figure 8. Control effect of different set currents before compensation

下载: 全尺寸图片幻灯片

图 9 补偿后不同给定电流控制效果

Figure 9. Control effect of different set currents after compensation

下载: 全尺寸图片幻灯片

3.3 速度环试验结果

该补偿算法在速度闭环中的控制效果、所示. 显示了补偿前，当给定速度为0.8、0.9、1.0 m/s时， ${\theta _{\text{e}}}$ 、 ${\hat \theta _{\text{e}}}$ 、 $\Delta {\theta _{\text{e}}}$ 与跟踪速度误差 $\Delta v$ 的变化，可以看出，角度误差分别在−0.23～0.53、−0.18～0.51、−0.17～0.49 rad，稳态速度最大跟踪误差为0.36 m/s左右. 图11显示了补偿后，对应的角度误差分别降至了−0.21～0.17 rad，−0.17～0.15 rad，−0.15～0.12 rad，稳态速度最大跟踪误差降至0.21 m/s左右. 可以发现，补偿后最大角度误差分别减小约67.9%、70.5%和75.5%，平均约70%左右，速度跟踪误差减小50%左右.

图 10 补偿前不同给定速度控制效果

Figure 10. Control effect of different set speeds before compensation

下载: 全尺寸图片幻灯片

综合上述，由于齿槽效应、端部效应以及长定子铁芯间气隙通断影响，当电机达到稳态时，电流环反馈电流与速度环反馈速度，在给定值上下存在一定波动，如表2、3所示. 当考虑延时对控制周期内电角度相移补偿后，电角度估计误差明显减小，且电流与速度跟踪性能也有一定的提升. 针对控制延时与采样延时导致的电角度观测误差问题，一方面可以进一步优化系统控制器性能，缩短采样周期；另一方面可以通过算法对控制系统总延时计算补偿.

表 2 不同电流下补偿前、后误差最大波动

Table 2. Maximum fluctuation of error before and after compensation under different currents

电流/A	补偿前误差		补偿后误差		波动变化
电流/A	电角/ rad	电流/ A	电角/ rad	电流/ A	电角/ %	电流/ %
20	0.86	7.5	0.23	7.5	73.3	0
21	0.71	7.5	0.21	7.5	70.4	0
22	0.68	7.5	0.19	7.5	72.1	0

下载: 导出CSV

| 显示表格

表 3 不同速度下补偿前、后误差最大波动

Table 3. Maximum fluctuation of error before and after compensation under different speeds

速度/（m·s⁻¹）	补偿前误差		补偿后误差		波动变化
速度/（m·s⁻¹）	电角/ rad	速度/ （m·s⁻¹）	电角/ rad	速度/ （m·s⁻¹）	电角/ %	速度/ %
0.8	0.53	0.36	0.17	0.21	67.9	50
0.9	0.51	0.36	0.15	0.21	70.5	50
1.0	0.49	0.36	0.12	0.21	75.5	50

下载: 导出CSV

| 显示表格

4. 结　论

针对磁浮列车零低速无传感控制高频响应电流电角度相位滞后与速度跟踪累积误差对列车控制精度的影响，提出一种考虑相移补偿的磁浮列车零低速高频信号注入无传感控制方法，得到以下结论：

1）所提出的估计角度偏差最小化寻优方法能够准确计算控制周期内的控制延时，反映相移误差.

2）补偿后的估计电角度误差减小，参与电流闭环对电流波动影响不大.

3）补偿后的角速度参与速度闭环能够在一定程度上提升速度跟踪控制性能.

本文提出的考虑相移补偿高频注入无传感方法可以准确地计及延时导致的相移滞后的影响，有利于实现列车零低速速度跟踪误差的补偿，从而提高高速磁浮无传感控制性能.

图 11 补偿后不同给定速度控制效果

Figure 11. Control effect of different set speeds after compensation

下载: 全尺寸图片幻灯片

图 1 水平加筋体受力示意

Figure 1. Sketch map of horizontal reinforcement

下载: 全尺寸图片幻灯片

图 2 水平加筋体力学计算模型

Figure 2. Calculation model of horizontal reinforcement

下载: 全尺寸图片幻灯片

图 3 水平加筋体锚固区示意

Figure 3. Sketch map of anchorage area

下载: 全尺寸图片幻灯片

图 4 变形关系示意

Figure 4. Illustration of relation between deformation

下载: 全尺寸图片幻灯片

图 5 某水平加筋体的张拉力-应变-时间特征曲线

Figure 5. Tension-strain-time curve of horizontal reinforcement

下载: 全尺寸图片幻灯片

表 1 模型试验和数值计算的参数

Table 1. Parameters of model test and numerical analysis

参数	取值
填方高度/m	1
填料重度/（kN•m^–3）	15.65
填料内摩擦角/（°）	36
溶洞直径/m	2.2
水平加筋体拉伸模量/（kN•m^–1）	2 988
附加荷载/kN	0
锚固长度/m	3

下载: 导出CSV

表 2 路基面变形对比

Table 2. Comparison of deformation of subgrade surface

模型试验	数值分析	RAFAEL方法	本文方法
0.23	0.21	0.18	0.25

下载: 导出CSV

COOPER H A, SAUNDERS M J. Road and bridge construction across gypsum karst in England[J]. Engineering Geology, 2002, 65(1): 217-223.

GIROUD J P, BONAPARTE R, BEECH J F, et al. Design of soil layer-geosynthetic systems overlying voids[J]. Geotextiles and Geomembranes, 1990, 9(1): 11-50. doi: 10.1016/0266-1144(90)90004-V

WANG M C, FENG Y X. Stability of geosynthetic-reinforced soil above a cavity[J]. Geotextiles and Geomembranes, 1996, 14(1): 95-109.

VILLARD P, BRIANCON L. Design of geosynthetic reinforcements for platforms subjected to localized sinkholes[J]. Canadian Geotechnical Journal, 2008, 45(2): 196-209. doi: 10.1139/T07-083

吕伟华,缪林昌,王非. 基于不完全土拱效应的土工格栅加固机制与设计方法[J]. 岩石力学与工程学报,2012,31(3): 632-639. doi: 10.3969/j.issn.1000-6915.2012.03.023

LÜ Weihua, MIAO Linchang, WANG Fei. Mechanism of geogrid reinforcement based on partially developed soil arch effect and design method[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(3): 632-639. doi: 10.3969/j.issn.1000-6915.2012.03.023

王非,缪林昌. 落水洞上覆路堤土工加筋设计新方法[J]. 东南大学学报(自然科学版),2012,31(3): 632-639.

WANG Fei, MIAO Linchang. New design method of geosynthetic-reinforced Embankment over sinkholes[J]. Journal of Southeast University (Natural Science Edition), 2012, 31(3): 632-639.

贺炜,付宏渊,殷苗苗. 岩溶区路堤下水平加筋体受力分析方法研究[J]. 水文地质工程地质,2010,37(6): 33-37. doi: 10.3969/j.issn.1000-3665.2010.06.007

HE Wei, FU Hongyuan, YIN Miaomiao. Study on calculating method of horizontally reinforcement beneath embankment constructed in Karst terrain[J]. Hydrogeology & Engineering Geology, 2010, 37(6): 33-37. doi: 10.3969/j.issn.1000-3665.2010.06.007

陈仁朋,汪焱卫,陈金苗. 桩网结构路基格栅加筋作用及其拉力特性研究[J]. 西南交通大学学报,2016,51(6): 1080-1086. doi: 10.3969/j.issn.0258-2724.2016.06.006

CHEN Renpeng, WANG Yanwei, CHEN Jinmiao. Reinforcing mechanism and tension behaviors of geogrid in pile-supported reinforced embankment[J]. Journal of Southwest Jiaotong University, 2016, 51(6): 1080-1086. doi: 10.3969/j.issn.0258-2724.2016.06.006

TERZAGHI K. Theoretical soil mechanics[M]. New York: John Wiley & Sons, 1943: 37 -42

郑俊杰,周燕君,曹文昭,等. 不同网孔形状格栅加筋土界面特性试验研究[J]. 西南交通大学学报,2017,52(3): 482-488. doi: 10.3969/j.issn.0258-2724.2017.03.007

ZHENG Junjie, ZHOU Yanjun, CAO Wenzhao, et al. Experimental study of interface behavior of geogrids with different aperture shapes[J]. Journal of Southwest Jiaotong University, 2017, 52(3): 482-488. doi: 10.3969/j.issn.0258-2724.2017.03.007

British Standard Institution. Code of practice for strengthened/reinforced soils and other fills: BS8006-1—2010[S]. London: British Standard Institution, 2010

LEONARD J W. Tension structures[M]. NewYork: McGraw-Hill, 1988: 35-49

朱斌,陈若曦,陈云敏. 抗沟渠型空洞水平加筋体的作用机理及设计方法[J]. 中国公路学报,2009,22(1): 11-16. doi: 10.3321/j.issn:1001-7372.2009.01.002

ZHU Bin, CHEN Ruoxi, CHEN Yunmin. Action mechanism and design method of horizontal reinforcement subjected to trench void[J]. China Journal of Highway and Transport, 2009, 22(1): 11-16. doi: 10.3321/j.issn:1001-7372.2009.01.002

HUCKERT A, BRIANON L, VILLARD P, et al. Load transfer mechanisms in geotextile-reinforced embankments overlying voids:Experimental and analytical approaches[J]. Geotextiles and Geomembranes, 2016, 44(1): 442-456.

VILLARD P, HUCKERT A, BRIANON L, et al. Load transfer mechanisms in geotextile-reinforced embankments overlying voids:Numerical approach and design[J]. Geotextiles and Geomembranes, 2016, 44(1): 381-395.

铁道第四勘察设计院. 铁路路基土工合成材料应用设计规: TB10118—2006[S]. 北京: 中国铁道出版社, 2006

中华人民共和国水利部. 土工合成材料应用技术规范: GB/T50290—2014[S]. 北京: 中国计划出版社, 2014

GEOTECHNIK. Recommendations for design and analysis of earth structures using geosynthetic reinforcements—EBGEO[M]. Berlin: Ernst & Sohn, 2011: 97-99

胡一峰, 李怒放. 高速铁路无砟轨道路基设计原理[M]. 北京: 中国铁道出版社, 2010: 194-195

期刊类型引用(15)

1.	徐秀丽，钱思博，李洪涛，吴冲，王仁贵，李琦，丁兴国，李枝军，李雪红. 新型开口肋正交异性钢桥面板疲劳性能对比试验研究. 土木工程学报. 2025(03): 50-67 . 百度学术
2.	鲁乃唯，刘静，王凯，王鸿浩. 随机轮载下钢桥面顶板与纵肋焊缝疲劳裂纹扩展特性研究. 振动与冲击. 2024(01): 193-201 . 百度学术
3.	宋旭明，罗祥威，季苏瑶，唐冕. 基于断裂力学的横隔板弧形切口CFRP加固参数研究. 中南大学学报(自然科学版). 2024(03): 1046-1055 . 百度学术
4.	邹晓琴. 正交异性钢桥面板U肋顶板全熔透焊缝制造工艺与疲劳性能研究. 公路. 2024(09): 72-78 . 百度学术
5.	彭振楠，宋神友，陈焕勇. 正交异性钢桥面顶板U肋全熔透焊接接头研发及疲劳性能研究. 广东公路交通. 2024(06): 6-11 . 百度学术
6.	袁毅，彭最. 武汉汉江湾桥主桥设计关键技术. 桥梁建设. 2023(01): 9-15 . 百度学术
7.	宋松科，张翔，杜桃明，侯宇博，权新蕊. 焊接微裂纹缺陷对顶板与纵肋构造疲劳性能影响研究. 桥梁建设. 2023(01): 24-32 . 百度学术
8.	雷丙超. 焊接缝疲劳开裂对斜拉桥钢桥面板的影响及加固方式. 现代机械. 2023(02): 98-102 . 百度学术
9.	王晓，郭亚文，梁爱霞，张莹，张雪松. 基于断裂力学的钢桥面板纵肋顶板焊缝疲劳效应分析. 公路交通科技. 2023(03): 105-114 . 百度学术
10.	赵秋，唐琨，黄冠铭，林铮哲. 对接焊缝初始缺陷对疲劳寿命影响规律研究. 合肥工业大学学报(自然科学版). 2023(04): 473-480 . 百度学术
11.	王大鹏，王景春，张先萌，张振波，徐强. 考虑列车荷载的隧道衬砌裂缝疲劳扩展可靠性分析. 中国铁道科学. 2023(06): 95-102 . 百度学术
12.	张亚海，郭宝圣，张卫国，陈方怀. 钢桥面板-纵肋双面焊缝疲劳裂纹应力强度因子. 土木与环境工程学报(中英文). 2022(03): 62-70 . 百度学术
13.	张清华，李俊，袁道云，朱金柱，崔闯. 高疲劳抗力钢桥面板的疲劳问题Ⅰ:模型试验. 中国公路学报. 2021(03): 124-135 . 百度学术
14.	李丹，张启伟，姜旭. 钢桥面板纵肋-顶板焊缝细节疲劳裂纹的应力强度因子分析. 城市道桥与防洪. 2021(07): 249-252+257+25 . 百度学术
15.	张清华，李俊，袁道云，朱金柱，崔闯. 深圳至中山跨江通道钢桥面板结构疲劳试验研究. 土木工程学报. 2020(11): 102-115 . 百度学术

其他类型引用(16)