硬岩掘进机（tunnel boring machine, TBM）广泛应用在城市饮水工程、地下隧洞工程，以及公路、铁路的建设中。相比传统的钻爆法，使用TBM进行施工具有自动化程度高、掘进速度快、施工质量高、环境影响小、安全性好的特点。根据设备的结构特点，TBM一般应用在硬岩的地质条件中，其刀盘要承受硬岩破碎的巨大变载荷，通过机架、液压缸和撑靴，最终传递到岩壁表面。因为外部载荷引起的TBM掘进波动较大，严重时甚至发生偏航、堵转等情况，因此通过建立推进机械系统的动力学模型，研究TBM推进系统在刀盘变载荷和变支撑刚度条件下的动力学特性有重要意义。

随着国内高铁、地铁和隧道工程的大量开工建设，TBM的推进系统的动力学研究近来已经成为热点。谢启江等研究了岩石变形和不同撑靴支撑刚度下TBM直线和曲线推进时动态特性响应，发现曲线推进时撑靴与岩石的接触刚度变化对运动特性影响较大^[1]。周小磊等应用AMESim软件对TBM推进系统进行了建模，模拟了推进液压缸快速回收和高压直线推进过程，其响应曲线和实际分析基本相符^[2]。张良贵对TBM推进系统进行了简化，计算了复合地质条件下所需的推力和扭矩^[3]。Huo等建立了整机多自由度模型，通过仿真和现场测试，得到了TBM整机在扭转、倾覆和摆动方向的响应^[4]。Sun等对TBM的刀盘系统和旋转驱动系统，建立了多自由度非线性动力学模型，分析了内外部耦合激励条件下的刀盘动力学响应^[5-6]。Xia等通过滚刀实验台的破岩实验和数值仿真，对刀盘受力进行了分析和预测^[7-8]。Gong等研究了TBM掘进贯入参数和岩石特性之间的关系，通过对实测数据库的分析，建立了掘进贯入度预测模型^[9]。相近的研究还包括对土压平衡型掘进机推进系统的研究, 其研究方法可以部分借用到硬岩掘进机的推进系统。刘杰等对土压平衡型掘进机推进系统，建立了五自由度无阻尼动力学模型，计算了该模型的响应^[10]。王洪新等通过理论结合实验分析了土压平衡盾构机推进系统的平衡控制方法^[11]。霍军周等建立了多关节面TBM刀具系统的多自由度耦合动力学模型，并实验证明接头表面对刀具系统的振动性能有很大的影响^[19]。

在之前的研究工作中，有一些因素并没有深入考虑，例如撑靴支撑液压缸和撑靴与岩壁变刚度因素，以及在不同初始条件下，动态特性的响应等。本文将综合考虑以上因素，通过对TBM推进系统多自由度动力学模型的建立和分析计算，得到一些更加具体的结论。

1 TBM推进系统动力学建模

敞开式TBM主机结构如图 1所示，由刀盘、驱动装置，主机架、推进液压缸、支撑液压缸、撑靴和扭矩液压缸组成。其中刀盘在驱动装置的扭矩驱动下，进行旋转运动，实现滚刀对岩石的滚动切割；推进液压缸用来推进主机架和刀盘做直线进给运动，实现滚刀对岩石的贯入切割，两侧的推进液压缸可以实现不同长度的伸缩，以实现盾构机向左右不同方向的偏航，用来完成隧道转弯处的施工。

盾构机的推进系统由推进液压缸，支撑液压缸、扭矩液压缸和撑靴组成。其中扭矩液压缸只用来抵消盾构机刀盘旋转的扭转力矩，在推进过程中把它和主机架看成一体。在推进过程中，推进千斤顶的推进速度很小，可以近似为静态。在正常掘进过程中，可以认为推力系统的载荷是平稳的。因此，进行以下假设：1）刀盘、主机架和扭矩液压缸可以看成一体，作为一个刚体处理；2）相对于整机和撑靴的质量，液压缸质量被忽略不计，并且近似为弹簧。

简化的盾构机模型如图 2所示。刀盘、主机架和扭矩液压缸简化为质量块$m_1$，撑靴分别简化为$m_2$和$m_3$。在$x$方向上，各质量块的位移为$x_1、x_2、x_3$，在$y$方向上，各质量块的位移为$y_1、y_2、y_3$；主机和撑靴之间的刚度为$k_{j12}、k_{j13}、k_{12x}、k_{13x}$，撑靴和隧道壁之间的刚度为$k_{2dx}、k_{2dz}、k_{3dx}、k_{3dz}$。$θ$为绕$y$方向的逆时针转动角度，$J_x$为绕$y$方向的转动惯量。

多自由度系统的受迫振动方程式为^[12]

$ \left\{ \begin{array}{l} \mathit{\boldsymbol{M\ddot q}} + \mathit{\boldsymbol{C\dot q}} + \mathit{\boldsymbol{Kq}} = \mathit{\boldsymbol{F}}\\ \mathit{\boldsymbol{J\ddot \theta }} + \mathit{\boldsymbol{C\dot \theta }} + \mathit{\boldsymbol{K\theta }} = \mathit{\boldsymbol{T}} \end{array} \right. $

(1)

式中：$\mathit{\boldsymbol M}$为系统的质量矩阵，$\mathit{\boldsymbol C}$为系统阻尼矩阵，$\mathit{\boldsymbol K}$为系统的刚度矩阵，$\mathit{\boldsymbol J}$是系统的转动惯量矩阵，$\mathit{\boldsymbol F}$是主机受到的载荷力，$\mathit{\boldsymbol T}$是主机受到的偏载力矩。

通过拉格朗日方程可以建立TBM系统的动力学模型。其中推进系统的动能为

$ \begin{array}{*{20}{c}} {T = \frac{1}{2}{m_1}\dot x_1^2 + \frac{1}{2}{m_1}\dot z_1^2 + \frac{1}{2}{m_2}\dot x_2^2 + \frac{1}{2}{m_2}\dot z_2^2 + }\\ {\frac{1}{2}{m_3}\dot x_3^2 + \frac{1}{2}{m_3}\dot z_3^2 + \frac{1}{2}{J_1}\dot \theta _1^2 + \frac{1}{2}{J_2}\dot \theta _2^2 + \frac{1}{2}{J_3}\dot \theta _3^2} \end{array} $

(2)

推进系统的势能为

$ \begin{array}{*{20}{c}} {U = \frac{1}{2}{k_{j12}}{{\left( {\frac{{{z_2} - {z_1}}}{{\sin {{30}^ \circ }}}} \right)}^2} + \frac{1}{2}{k_{j13}}{{\left( {\frac{{{z_3} - {z_1}}}{{\sin {{30}^ \circ }}}} \right)}^2} + }\\ {\frac{1}{2}{k_{12x}}\left( {{x_2} - {x_1}} \right)2 + \frac{1}{2}{k_{2dx}}x_2^2 + \frac{1}{2}{k_{2dz}}z_2^2 + }\\ {\frac{1}{2}{k_{13x}}\left( {{x_3} - {x_1}} \right)2 + \frac{1}{2}{k_{3dx}}x_3^2 + \frac{1}{2}{k_{3dz}}z_3^2 + }\\ {\frac{1}{2}{k_{12q}}{{\left( {{{\dot \theta }_2} - {{\dot \theta }_1}} \right)}^2} + \frac{1}{2}{k_{2rq}}\dot \theta _2^2 + }\\ {\frac{1}{2}{k_{13q}}{{\left( {{{\dot \theta }_3} - {{\dot \theta }_1}} \right)}^2} + \frac{1}{2}{k_{3rq}}\dot \theta _3^2} \end{array} $

(3)

推进系统的能量耗散函数为

$ \begin{array}{*{20}{c}} {D = \frac{1}{2}{C_{j12}}{{\left( {\frac{{{{\dot z}_2} - {{\dot z}_1}}}{{\sin {{30}^ \circ }}}} \right)}^2} + \frac{1}{2}{C_{j13}}{{\left( {\frac{{{{\dot z}_3} - {{\dot z}_1}}}{{\sin {{30}^ \circ }}}} \right)}^2} + }\\ {\frac{1}{2}{C_{2dx}}\dot z_2^2 + \frac{1}{2}{C_{3dx}}\dot z_3^2 + \frac{1}{2}{C_{12x}}{{\left( {{{\dot x}_2} - {{\dot x}_1}} \right)}^2} + }\\ {\frac{1}{2}{C_{2dx}}x_2^2 + \frac{1}{2}{C_{13x}}{{\left( {{{\dot x}_3} - {{\dot x}_1}} \right)}^2} + \frac{1}{2}{C_{3dx}}x_3^2 + }\\ {\frac{1}{2}{C_{12q}}{{\left( {{{\dot \theta }_2} - {{\dot \theta }_1}} \right)}^2} + \frac{1}{2}{C_{2rq}}\dot \theta _2^2 + }\\ {\frac{1}{2}{C_{13q}}{{\left( {{{\dot \theta }_3} - {{\dot \theta }_1}} \right)}^2} + \frac{1}{2}{C_{3rq}}\dot \theta _3^2} \end{array} $

(4)

拉格朗日方程：

$ \frac{{\rm{d}}}{{{\rm{d}}t}}\frac{{\partial T}}{{\partial {{\dot q}_i}}} - \frac{{\partial T}}{{\partial {q_i}}} + \frac{{\partial U}}{{\partial {q_i}}} + \frac{{\partial D}}{{\partial {{\dot q}_i}}} = {F_i}\left( t \right) $

(5)

式中:$q_i、q·_i$为系统的广义坐标和广义速度，$T、U、D$为系统的动能、势能与能量耗散函数，$F_i(t)$为广义激振力。将式(2)~(4) 代入到式（5）中，可以得到动力学方程：

$ \begin{array}{*{20}{c}} {{m_1}{{\ddot x}_1} + \left( {{k_{12x}} + {k_{13x}}} \right){x_1} - {k_{12x}}{x_2} - {k_{13x}}{x_3} + }\\ {\left( {{c_{12x}} + {c_{13x}}} \right){{\dot x}_1} - {c_{12x}}{{\dot x}_2} - {c_{13}}{{\dot x}_3} = {F_x}} \end{array} $

(6)

$ \begin{array}{*{20}{c}} {{m_2}{{\ddot x}_2} - {k_{12x}}{x_1} + \left( {{k_{2dx}} + {k_{12x}}} \right){x_2} - {c_{12x}}{{\dot x}_1} + }\\ {\left( {{c_{2dx}} + {c_{12x}}} \right){{\dot x}_2} = 0} \end{array} $

(7)

$ \begin{array}{*{20}{c}} {{m_3}{{\ddot x}_3} - {k_{13x}}{x_1} + \left( {{k_{13x}} + {k_{3dx}}} \right){x_3} - {c_{13x}}{{\dot x}_1} + }\\ {\left( {{c_{13x}} + {c_{3dx}}} \right){{\dot x}_3} = 0} \end{array} $

(8)

$ \begin{array}{*{20}{c}} {{m_1}{{\dot z}_1} + 4\left( {{k_{j12}} + {k_{j13}}} \right){z_1} - 4{k_{j12}}{z_2} - 4{k_{j13}}{z_3} + }\\ {4\left( {{c_{j12}} + {c_{j13}}} \right){{\dot z}_1} - 4{c_{j12}}{{\dot z}_2} - 4{c_{j13}}{{\dot z}_3} = {F_z}} \end{array} $

(9)

$ \begin{array}{*{20}{c}} {{m_2}{{\ddot z}_2} - 4{k_{j12}}{z_1} + \left( {{k_{j13}} + {k_{2dx}}} \right){z_2} - 4{c_{j12}}{{\dot z}_1} + }\\ {\left( {4{c_{j13}} + {c_{2dz}}} \right){{\dot z}_2} = 0} \end{array} $

(10)

$ \begin{array}{*{20}{c}} {{m_3}{{\ddot z}_3} - 4{k_{j13}}{z_1} + \left( {4{k_{j13}} + {k_{3dz}}} \right){z_3} - 4{c_{j13}}{{\dot z}_1} + }\\ {\left( {4{c_{j13}} + {c_{3dz}}} \right){{\dot z}_3} = 0} \end{array} $

(11)

$ \begin{array}{*{20}{c}} {{J_1}{{\ddot \theta }_1} + \left( {{k_{12q}} + {k_{13q}}} \right){\theta _1} - {k_{12q}}{\theta _2} - {k_{13q}}{\theta _3} + }\\ {\left( {{c_{12q}} + {c_{13q}}} \right){{\dot \theta }_1} - {c_{12q}}{{\dot \theta }_2} - {c_{13q}}{{\dot \theta }_3} = {T_y}} \end{array} $

(12)

$ \begin{array}{*{20}{c}} {{J_2}{{\ddot \theta }_2} - {k_{12q}}{\theta _1} + \left( {{k_{12q}} + {k_{2rq}}} \right){\theta _2} - {c_{12q}}{{\dot \theta }_1} + }\\ {\left( {{c_{12q}} + {c_{2rq}}} \right){{\dot \theta }_2} = 0} \end{array} $

(13)

$ \begin{array}{*{20}{c}} {{J_3}{{\ddot \theta }_3} - {k_{13q}}{\theta _1} + \left( {{k_{13q}} + {k_{3rq}}} \right){\theta _2} - {c_{13q}}{{\dot \theta }_1} + }\\ {\left( {{c_{13q}} + {c_{3rq}}} \right){{\dot \theta }_2} = 0} \end{array} $

(14)

2 动力学模型参数计算 2.1 外部激励计算

在设计TBM时，根据不同的土壤环境会采用不同的滚刀分布。常见的滚刀分布如图 3所示，其中部布置中心滚刀，周围布置正滚刀，边缘布置边滚刀。

通过滚刀的受力分析^[13]，破岩过程中的滚刀受力可以表示为：正推力$F_v$、滚动力$F_r$和侧推力$F_s$。正滚刀和边滚刀的受力如图 4所示。通过LS-DYNA的受力仿真以及科罗拉多矿业学院的滚刀受力预测模型，结合试验可以得到单个滚刀的受力^[14]。

根据单个滚刀的受力关系和滚刀的位置分布，可以列出主机受外部载荷的受力公式：

$ {F_x} = {F_{vx\Sigma }} + {F_{rx\Sigma }} + {F_{sx\Sigma }} $

(15)

$ {F_z} = \sum\limits_{i = 1}^n {{F_{vi}}} + \sum\limits_{j = 1}^m {{F_{vj}}} + \sum\limits_{k = 1}^p {\left( {{F_{vk}} + {F_{sk}}\sin {\beta _k}} \right)} $

(16)

$ {T_y} = \sum\limits_{i = 1}^n {{F_{vi}}{l_{yi}}} + \sum\limits_{j = 1}^m {{F_{vj}}{l_{yj}}} $

(17)

式中：$n$为中心滚刀的数量，$m$为正滚刀的数量，$p$为边滚刀的数量；$F_{vx∑}$为正推力在X方向的分力，$F_{rx∑}$为滚动力在$X$方向的分力，$F_{sx∑}$为侧推力在X方向的分力。

2.2 撑靴接触面刚度计算

TBM掘进过程中，撑靴和岩石表面的摩擦力提供驱动力和扭转反向力。由于刀盘阻力和推进载荷具有较大波动，所以撑靴和岩石截面间的接触刚度特性对其动力学特性非常重要。

隧道岩石的粗糙面可以认为是由很多小凸体随机分布，小凸体在受到法向载荷和切向载荷作用时，其变形和力的加载路径有关。通过使用修正各向同性三维分形表面方法来表征撑靴和岩壁的粗糙度，从而计算出撑靴和岩石之间的方向刚度和切向刚度^[15]。

$ \begin{array}{*{20}{c}} {{k_x} = \frac{4}{{3\sqrt {2{\rm{\pi }}} }}E\frac{{\left( {4 - D} \right)\left( {D - 1} \right)}}{{\left( {3 - D} \right)\left( {2 - D} \right)}}\left( {\alpha _L^{'1/2} - } \right.}\\ {\left. {\alpha _L^{'\left( {D - 1} \right)/2}\alpha _c^{'\left( {2 - D} \right)/2}} \right)} \end{array} $

(18)

$ \left\{ \begin{array}{l} {k_z} = \frac{{8Gr'}}{{\left( {2 - v} \right)}}{\left[ {\mu \frac{{{\rm{d}}P}}{{{\rm{d}}Q}} + \left( {1 - \frac{{{\rm{d}}P}}{{{\rm{d}}Q}}} \right){{\left( {1 - \frac{Q}{{\mu P}}} \right)}^{ - 1/3}}} \right]^{ - 1}},\\ \;\;\;\;\;\;\;\;\;\;0 \le \frac{{{\rm{d}}P}}{{{\rm{d}}Q}} \le \frac{1}{\eta }\\ {k_z} = \frac{{8Gr'}}{{\left( {2 - v} \right)}},\;\;\;\;\;\frac{{{\rm{d}}P}}{{{\rm{d}}Q}} \ge \frac{1}{\eta } \end{array} \right. $

(19)

式中: $D$为分形维数，$E$为岩石的杨氏模量，$a′_{c}$为临界面积，$α′_{L}$为所有接触的小凸体中发成变形后的最大接触面的体积，$P$为法向力，$r′$为名义接触半径，$Q$为切向力，$μ$为摩擦系数，$G$为剪切模量，$v$为泊松比。TMB推进过程中，撑靴接触表面的切向力、法向力及临界面积会不断变化，法向和切向刚度也会相应的改变。

2.3 外液压缸刚度计算

在重载荷慢动作过程中，液压缸的刚度成为系统的重要影响因素。液压缸的刚度计算不仅需要考虑缸内液压油的体积，还必须考虑液压管路中的液压油。通过文献[16]可以得到液压缸的刚度计算公式为

$ {k_j} = E\left( {\frac{{A_1^2}}{{{V_{L1}} + {V_1}}} + \frac{{A_2^2}}{{{V_{L2}} + {V_2}}}} \right) $

(20)

式中：$E$为液压油体积弹性模量，$V_{L1}$为输入管路中液压油的体积，$V_{L2}$为输出管路中液压油的体积，$V_1$为液压缸大端液压油的体积，$V_2$为液压缸小端液压油的体积，$A_1$为液压缸大端活塞面积，$A_2$为液压缸小端活塞面积。

2.4 阻尼计算

系统中各部件之间均为固体表面的滑移接触，所以通过以下公式计算：

$ c = 2\xi \sqrt {{m_e}{k_e}} $

(21)

式中：ξ为阻尼比，一般取值0.02~0.05，本文中取值0.02；$m_e$为等效质量；$k_e$为等效刚度。

3 动力学模型求解

本文使用的Newmark-β法是一种隐式积分方法，是无条件稳定的方法。假定在时间间隔[t, t+Δt]内加速度呈线性变化，基本假定为

$ {{\dot x}_{t + \Delta t}} = {{\dot x}_t} + \left[ {\left( {1 - \delta } \right){{\ddot x}_t} + \delta {{\ddot x}_{t + \Delta t}}} \right]\Delta t $

(22)

$ {x_{t + \Delta t}} = {x_t} + {{\dot x}_t}\Delta t + \left[ {\left( {\frac{1}{2} - \beta } \right){{\ddot x}_t} + \beta {{\ddot x}_{t + \Delta t}}} \right]\Delta {t^2} $

(23)

式中δ、β是按积分精度和稳定性要求可以调整的参数。

根据式(22)~（23），可以用$x_{t+Δt}、\dot x_t、\ddot x_t$表示$\dot x_{t+Δt}，\ddot x_{t+Δt}$，即：

$ {{\dot x}_{t + \Delta t}} = \frac{1}{{\beta \Delta {t^2}}}\left( {{x_{t + \Delta t}} - {x_t}} \right) - \frac{1}{{\beta \Delta {t^2}}}{{\dot x}_t} - \left( {\frac{1}{{2\beta }} - 1} \right)\Delta t{{\ddot x}_t} $

(24)

$ {{\ddot x}_{t + \Delta t}} = \frac{\delta }{{\beta \Delta t}}\left( {{x_{t + \Delta t}} - {x_t}} \right) + \left( {1 - \frac{\delta }{\beta }} \right){{\dot x}_t} + \left( {1 - \frac{\delta }{{2\beta }}} \right)\Delta t{{\ddot x}_t} $

(25)

每步积分应该满足$t+Δt$时刻的动力学方程：

$ M{{\ddot x}_{t + \Delta t}} + C{{\dot x}_{t + \Delta t}} + K{x_{t + \Delta t}} = {R_{t + \Delta t}} $

(26)

将式（23）~（24）代入式（25）可得关于$x_{t+Δt}$的方程为

$ \dot K{x_{t + \Delta t}} = {{\dot R}_{t + \Delta t}} $

(27)

式中：

$ \dot K = K + \frac{\delta }{{\beta \Delta t}}C + \frac{1}{{\beta \Delta {t^2}}}M $

(28)

$ \begin{array}{l} {{\dot R}_{t + \Delta t}} = {R_{t + \Delta t}} + M\left[ {\frac{1}{{\beta \Delta {t^2}}}{x_t} + \frac{\delta }{{\beta \Delta t}}{{\dot x}_t} + \left( {\frac{1}{{2\beta }} - 1} \right){{\ddot x}_t}} \right] + \\ \;\;\;\;\;\;\;\;\;\;\;C\left[ {\frac{\delta }{{\beta \Delta t}}{x_t} + \left( {\frac{\delta }{\beta } - 1} \right){{\dot x}_t} + \left( {\frac{\delta }{{2\beta }} - 1} \right)\Delta t{{\ddot x}_t}} \right] \end{array} $

(29)

4 动力学模型的实际案例适用分析 4.1 动力学响应求解

以辽宁大伙房水库施工的地质条件和采用的Robins盾构机为研究对象，其参数如表 1所示。

对滚刀和岩石的作用过程进行有限元数值仿真，并利用实测数据进行校正^[17-18]。通过式（15）~（17）的计算，可以得到如图 5的主机外部时程载荷图。

将载荷以及计算得到的刚度和阻尼代入动力学方程，应用Newmark-β法对微分方程组进行求解，可以得到系统的响应。图 6是各个方向的响应。

4.2 变初始条件下敏捷度分析

在盾构机掘进过程中，控制系统会不断矫正推进路线，使盾构机按照规划路径前进，因此盾构机的初始位移就很小。选取如下的初始位移：0≤x₁≤0.01 m，0≤x₂≤0.01 m，0≤z₁≤0.01 m，0≤z₂≤0.01 m，0≤θ₁≤0.1°，0≤θ₂≤0.1°。为了能够清楚地比较不同初始条件下，主机位移和角度响应变化。按照如下公式定义敏感系数$ν$：

$ \left\{ \begin{array}{l} {v_{11}}\left( t \right) = \frac{{{x_1}\left( t \right)}}{{{x_0}\left( t \right)}},{v_{12}}\left( t \right) = \frac{{{x_2}\left( t \right)}}{{{x_0}\left( t \right)}}\\ {v_{21}}\left( t \right) = \frac{{{z_1}\left( t \right)}}{{{z_0}\left( t \right)}},{v_{22}}\left( t \right) = \frac{{{z_2}\left( t \right)}}{{{z_0}\left( t \right)}}\\ {v_{31}}\left( t \right) = \frac{{{\theta _1}\left( t \right)}}{{{\theta _0}\left( t \right)}},{v_{32}}\left( t \right) = \frac{{{\theta _2}\left( t \right)}}{{{\theta _0}\left( t \right)}} \end{array} \right. $

(30)

图 7为敏感系数时程图。在X方向，实线是主机在初始位移0.01 m时，刀盘的敏感系数时程曲线，虚线是撑靴在初始位移为0.01 m时，刀盘敏感系数时程曲线。主机在X方向的偏航位移对刀盘的影响要远大于撑靴的偏置影响。主机在X方向偏航0.01 m，会导致刀盘在随后的最大波动幅值达到0.18 m。在Z方向，实线是主机在初始位移为0.01 m时，刀盘敏感系数时程曲线，虚线是撑靴在初始位移为0.01 m时，刀盘敏感系数时程曲线。主机在Z方向的敏感度差值平均为0.16，差值较小，该方向对偏航不会产生影响，仅影响每次进给量的误差。在绕Y轴摆动方向，实线是主机角度为1°时，刀盘的敏感系数时程曲线，虚线是撑靴在初始条件为1°时，刀盘的敏感系数时程曲线。在初始状态下，主机偏航1°，会导致刀盘在随后的最大波动幅值达到1.4°；而撑靴偏航1°，会导致刀盘的偏航幅度达到1.08°。

5 结论

1）推进系统的稳态响应具有明显的周期性，并且Z方向响应周期短于摆动方向的周期。

2）主机在X方向和绕Y方向的偏航会导致刀盘波动幅值的较大变化，从而导致掘进过程出现较大误差。

3）撑靴在初始条件下的偏航位移和偏转角度对整机的影响较小，整机的角度调节仍要以调整两侧推进液压缸为主。