甘肃积石山Ms6.2级地震诱发滑坡快速评估

曾营; 张迎宾; 储峰; 柳静; 冯振海; 苏金蓉

doi:10.3969/j.issn.0258-2724.20240041

甘肃积石山Ms6.2级地震诱发滑坡快速评估

doi: 10.3969/j.issn.0258-2724.20240041

曾营^1,,
张迎宾^1, ,,
储峰²,
柳静¹,
冯振海¹,
苏金蓉³

1.
西南交通大学土木工程学院，四川成都 610031
2.
中国路桥工程有限责任公司，北京 100011
3.
四川省地震局，四川成都 610041

基金项目: 国家自然科学基金项目（52378370，52278372）

详细信息

作者简介:
曾营（1997—），男，博士研究生，研究方向为地质灾害风险性评价，E-mail：zengying@my.swjtu.edu.cn

通讯作者:
张迎宾（1983—），男，教授，博士，研究方向为岩土地震工程，E-mail：yingbinz719@swjtu.edu.cn

中图分类号: P642.2
计量
- 文章访问数: 139
- HTML全文浏览量: 60
- PDF下载量: 25
- 被引次数: 0
出版历程
- 收稿日期: 2024-01-19
- 修回日期: 2024-04-10
- 网络出版日期: 2025-01-21

Rapid Assessment of Landslides Induced by Jishishan Ms6.2 Earthquake in Gansu Province

1.
School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China
2.
China Road & Bridge Corporation, Beijing 100011, China
3.
Sichuan Earthquake Administration, Chengdu 610041, China

摘要

摘要:
震后迅速获取同震滑坡分布及灾情评估对于应急救援和重建工作至关重要，采用IDNPM （InSAR data-newmark physical fusion driver model）方法对2023年12月18日甘肃积石山地震引发的滑坡进行快速评估，以期迅速精准掌握滑坡灾害的宏观分布. 首先，通过时序SBAS-InSAR揭示该地区有着严重的冲沟发育和溯源侵蚀现象，这些地质特征为滑坡提供了有利的孕育环境；其次，运用IDNPM方法对积石山地震进行滑坡快速评估，预测出赵木川村、塔沙坡村、大河家镇等地的陡峭斜坡及沟壑两侧为地震诱发滑坡的高风险区域；最后，综合实地考察、数值模拟及卫星识别技术，验证该模型在实际应用中的可靠性. 结果表明：全区共有2.657%的高风险区，需要重点关注此类区域；对已发生崩滑的坡体紧急清理和加固，对于未发生滑移的区域，应采取监测和评估措施，以防范可能发生的震后次生滑坡事件；研究成果可为受灾区的灾后应急救援和恢复重建工作提供有力的数据支撑.
- 地震 /
- 滑坡 /
- 应急评估 /
- 预测模型 /
- InSAR技术
Abstract:
Rapidly obtaining co-seismic landslide distribution and conducting disaster assessments after earthquakes are vital for effective emergency relief and reconstruction. Therefore, in this study, the InSAR data-Newmark physical fusion driver model (IDNPM) was used to rapidly assess the landslides triggered by the earthquake in Jishishan, Gansu Province on December 18, 2023, with a view to quickly and accurately grasping the macroscopic distribution of landslide hazards. Firstly, through the time series SBAS-InSAR, it was revealed that there was serious gully development and retrogressive erosion in this area. These geological characteristics provided a favorable breeding environment for landslides. Secondly, the IDNPM was used to quickly evaluate the landslide of Jishishan earthquake, and it was predicted that the steep slopes and gully sides of Zhaomuchuan Village, Tashapo Village, and Dahejia Town were the high-risk areas for earthquake-induced landslide. Finally, based on the field investigation, numerical simulation, and satellite identification technology, the reliability of the model in practical application was proven. The results indicate that a total of 2.657% of the region is at high risk. There is a need to focus on such zones by urgently clearing and stabilizing slopes where landslides have occurred. For areas where no landslides have occurred, monitoring and assessment measures should be taken to guard against possible post-earthquake secondary landslide events. The research results can provide strong data support for emergency relief and reconstruction work after earthquakes in the affected areas.
- earthquake /
- landslide /
- emergency assessment /
- prediction model /
- InSAR technology

HTML全文

磁悬浮轴承是通过磁场产生的电磁力将转子悬浮在空中的一种无摩擦的高性能轴承^[1]，具有无接触、损耗低、运行速度高及寿命长等优点，因此被广泛应用在航空航天、飞轮储能、人工心脏等领域^[2-9]. 从工作原理的角度，磁悬浮轴承分为主动型、被动型和混合型. 主动型磁悬浮轴承通过传感器和功率放大器等单元，主动控制通电线圈的电流，以此调控电磁力，从而实现转子的主动悬浮；被动型磁悬浮轴承则利用永磁体间、永磁体和软磁材料间的斥力或电磁吸力，产生不可控的电磁力以实现转子悬浮^[10]；混合磁悬浮轴承综合两者特点，通过永磁体生成偏置磁场和电磁线圈生成控制磁场，减少功率放大器的数量，并降低其功耗^[11-12].

混合磁悬浮轴承按其结构型式可分为同极性结构和异极性结构. 同极性磁轴承的偏置磁场和控制磁场在不同的平面，工作产生的损耗较低但轴向尺寸较长^[13]，限制了其转子的最高临界转速. 异极性磁轴承磁极的极性相反，极性交替排列，偏置磁通与控制磁通在同一平面，相比于同极性轴承，其轴向长度短，对提升转子最高转速有利^[14-17]. Okada等^[18]将传统磁轴承永磁体在两径向轴承间形成磁通回路的方式改变为将永磁体置于径向磁轴承内部，使得磁路在内部闭合，减小了磁轴承漏磁，提高了轴承效率. 孙津济等^[19]提出一种外转子式的新型磁轴承结构，在定子内部引入第二气隙，使得电磁磁路与永磁磁路共面，且不经过永磁体，同时，将定子铁心磁极数由8个减为4个，大大降低了磁轴承损耗. 朱润泽等^[20]提出一种4极型磁轴承，将传统的8极结构变为4极，并引入辅助气隙以防止永磁体发生短路，使得轴承在减少损耗的同时降低径向两自由度间的电磁耦合，有效提高轴承性能. 钟志贤等^[21]提出一种新型4极径向磁轴承，将外定子在轴向上进行拓展，形成具有6条独立磁路的结构，使其电磁磁路与永磁磁路相互独立，实现径向两自由度的电磁力解耦，但由于轴向较长，在一定程度上会影响转子的临界转速.

本文提出一种新型异极径向混合磁轴承（heteropolar radial hybrid magnetic bearing，HRHMB），可以控制径向2个自由度，具有轴向长度短、空间利用率高、永磁磁路与电磁磁路耦合距离短等优点. 首先，介绍新型HRHMB的结构原理，采用等效磁路法解析得到该磁轴承的电磁力、位移刚度和电流刚度，并与仿真结果进行分析对比；为验证该磁悬浮轴承的静态性能，在同悬浮力约束条件下与传统偏置磁轴承^[22]进行对比；最后，通过有限元分析验证该新型磁轴承在径向两自由度间的电磁力耦合较小.

1. 新型HRHMB结构及原理

1.1 新型HRHMB结构

新型HRHMB和传统八极偏置磁轴承的结构及磁路如图1所示. 新型HRHMB主要由定子铁心、转子铁心、控制绕组及永磁体组成. 定子铁心上包含主动控制磁1、3、5、7和辅磁极2、4、6、8，4组控制绕组分别在4个主动控制磁极上，产生控制磁场；4个辅磁极交错分布在主动控制磁极之间. 8个永磁体采用平行充磁的方式两两分布在主动控制磁极两侧，与转子形成工作气隙，用以产生偏置磁场. 定子铁心上的主辅磁极沿圆周形成NS相间的异极极性分布，为异极型磁悬浮轴承的结构形式. 同时，在定子铁心的横向处引入辅助气隙，起到减少X方向的永磁体在定子极间的漏磁及加强偏置磁场的作用. 由于异极型磁轴承损耗较大，其定转子铁心采用硅钢片叠压而成以降低转子高速运转下的铁耗. 永磁体采用钕铁硼烧结而成，嵌装在主辅磁极之间.

图 1 两种混合型磁悬浮轴承结构

Figure 1. Two kinds of hybrid magnetic bearing structures

下载: 全尺寸图片幻灯片

1.2 新型HRHMB工作原理

由于新型HRHMB具有对称性，以Y向磁极为例，控制绕组产生的控制磁通路径如图1（a）虚线所示，主动控制磁极上的控制绕组通入电流产生Y向的磁通，经过气隙、转子铁心及磁轭形成回路. 当转子处于平衡位置时，由于定子结构对称，Y方向上永磁体产生的气隙磁密相同，此时线圈中未通入电流，没有控制磁通. 假设转子受到向下的扰动力时， + Y方向气隙长度变大，−Y方向气隙长度变小. 在差动控制方式下，线圈通入电流，控制绕组产生的控制磁通和永磁体产生的偏置磁通在 + Y方向气隙处叠加，−Y方向气隙处抵消，使 + Y方向气隙磁通大于−Y方向气隙磁通，导致转子产生一个 + Y方向的电磁力，将转子拉回平衡位置. 线圈通电产生不同方向的力以使转子受到其它方向的干扰时保持平衡位置.

新型HRHMB的优点在于磁通路径和轴向长度短，且由于偏置磁通由永磁体产生，因此，功耗较低. 另外，其永磁磁路和电励磁磁路耦合距离较短，利于降低径向二自由度之间的磁场力耦合. 与传统八极径向偏置磁悬浮轴承相比，磁极数量同为八极，但新型HRHMB磁极区域大，在飞轮储能系统和磁悬浮电机等透平机械中有潜在的应用前景.

2. HRHMB磁路分析

2.1 等效磁路模型

单独考虑永磁体产生的偏置磁场和通电线圈产生的控制磁场，对HRHMB磁场进行仿真，结果如图2所示.

图 2 新型异极径向混合磁轴承磁场仿真

Figure 2. Magnetic field simulation of novel HRHMB

下载: 全尺寸图片幻灯片

由磁场仿真结果可以看出：偏置磁通经过永磁体、主动控制磁极、气隙、转子铁心和定子辅磁极，漏磁通经过永磁体、定子主磁极、定子磁轭、定子极间空隙和定子辅磁极；Y方向偏置磁通漏磁主要分布在定子的极间，由于X方向第二气隙的引入，X方向永磁体产生的偏置磁通漏磁大为减小；控制磁通经过主动控制磁极、气隙、转子铁心和定子磁轭.

假设该新型磁轴承定转子铁心磁阻及涡流损耗为零，建立其等效磁路模型，如图3所示. 在偏置磁场等效磁路（图3（a））中，R_gb1 ～ R_gb8为永磁体下偏置磁路气隙磁阻，转子处在平衡位置时，其结构对称，此时R_gb1=R_gb2，R_gb3=R_gb4；R_k1 ～ R_k8为永磁体漏磁磁阻；F_p为永磁体磁动势；R_p为永磁体磁阻.

图 3 新型异极径向混合磁轴承等效磁路

Figure 3. Equivalent magnetic circuit of novel HRHMB

下载: 全尺寸图片幻灯片

${F_{\text{p}}} = {H_{\text{c}}}{l_{\text{m}}},$

(1)

${R_{\text{p}}} = {l_{\text{m}}}/\left( {{\mu _{\text{p}}}{S_{\text{p}}}} \right),$

(2)

式中：H_c为永磁体矫顽力，l_m为永磁体宽度，$ {\mu _{\text{p}}} $为永磁体磁导率，S_p为永磁体中性面面积.

在控制磁场等效磁路（图3（b））中，R_c1、R_c2、R_c3、R_c4分别为4个控制磁极气隙的磁阻，R_a为辅助气隙磁阻，N_x、N_y分别为水平、垂直方向的线圈匝数，i_x、i_y分别为水平、垂直方向上的电流.

$\left\{ \begin{array}{l}\left[{R}_{\text{c}1}\;{R}_{\text{c}2}\;{R}_{\text{c}3}\;{R}_{\text{c}4}\right]=\left[\dfrac{{g}_{0} + y}{{\mu }_{0}{A}_{\text{r}}}\;\dfrac{{g}_{0}-y}{{\mu }_{0}{A}_{\text{r}}}\;\dfrac{{g}_{0} + x}{{\mu }_{0}{A}_{\text{r}}}\;\dfrac{{g}_{0}-x}{{\mu }_{0}{A}_{\text{r}}}\right]\\ {R}_{\text{a}}={g}_{\text{a}}/\left({\mu }_{0}{A}_{\text{r}}\right),\end{array},\right.$

(3)

式中：μ₀为空气磁导率，g₀为转子铁心在平衡位置时的气隙长度，g_a为辅助气隙长度，x、y分别为转子在X、Y轴方向的位移量，A_r为主动控制磁极气隙的面积.

2.2 新型异极径向混合磁轴承磁通计算

在高速电机和飞轮储能系统中转子最大允许的偏心位移为气隙长度的1/10，此时可认为4个径向方向偏置磁通彼此解耦^[23-24]. 将图3（a）偏置磁路图简化为Y + 、Y−、X + 、X− 4条磁路，分别计算其气隙处磁通，由于该磁轴承结构对称，由磁网络法计算Y + 、X−处可得到其他气隙处磁通.

2.2.1 Y + 方向偏置磁路磁通

偏置磁通在Y + 方向的等效磁路如图4（a）所示，φ_pmk1、φ_pmk2为漏磁磁通，φ_pmy1、φ_pmy2为永磁体下的偏置磁通；R_gb1、R_gb2为偏置磁通气隙磁阻，R_k2为漏磁磁阻，R_k1漏磁阻可以看作是端部漏磁磁阻R_k11和极间磁阻R_k12的并联磁阻，且R_k1//R_p//R_gb1.

图 4 Y + 方向等效磁路

Figure 4. Equivalent magnetic circuit along Y + direction

下载: 全尺寸图片幻灯片

漏磁磁阻和气隙磁阻分别为

$\left\{ \begin{array}{l}{R}_{\text{k}1}={R}_{\text{k2}}={R}_{\text{k}12}{R}_{\text{k}11}/({R}_{\text{k}12} + {R}_{\text{k}11}),\\ \left[{R}_{\text{gb}1}\;{R}_{\text{gb}2}\;{R}_{\text{gb}3}\;{R}_{\text{gb}4}\right]=\left[\dfrac{{g}_{0} + y}{{\mu }_{0}{A}_{\text{ra}}}\;\dfrac{{g}_{0} + y}{{\mu }_{0}{A}_{\text{ra}}}\;\dfrac{{g}_{0}-y}{{\mu }_{0}{A}_{\text{ra}}}\;\dfrac{{g}_{0}-y}{{\mu }_{0}{A}_{\text{ra}}}\right],\end{array}\right.$

(4)

式中：$ {R_{{\text{k11}}}} = \ln ({r_{{\text{a}}2}}/{r_{{\text{a1}}}})/\left( {2{\mu _0}{\text{π}} {L_{\text{z}}}} \right) $^[25]，r_a1和r_a2为两极间半径，$ {R}_{\text{k12}}=(2/{m}_{\text{y}} + 0.5)/\left(0.64{\mu }_{0}{L}_{\text{z}}\right) $，L_z为磁轴承轴向长度，m_y为端部漏磁路半径，A_ra为辅磁极截面积.

根据磁路基尔霍夫定律可求出Y + 方向2个永磁体产生的偏置磁通和漏磁磁通，对应图4（a）的磁网络矩阵如式（5）所示. 各磁阻由图4磁阻模型得出，可得偏置磁场的磁通${\phi} _{\text{pmy1}}$、${\phi}_{\text{pmy2}} $如式（6）所示.

$\quad\quad\quad\quad\quad \left[ {\begin{array}{*{20}{c}} {{R_{{\text{k}}1}}}&0&{ - {R_{\text{p}}}}&0 \\ 0&{{R_{{\text{k}}2}} + {R_{\text{p}}}}&0&{ - {R_{\text{p}}}} \\ { - {R_{\text{p}}}}&0&{{R_{{\text{gb}}1}} + {R_{{\text{c}}1}} + {R_{\text{p}}}}&{ - {R_{{\text{c}}1}}} \\ 0&{ - {R_{\text{p}}}}&{ - {R_{{\text{c}}1}}}&{{R_{{\text{gb2}}}} + {R_{{\text{c}}1}} + {R_{\text{p}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\phi _{{\text{pmk}}1}}} \\ {{\phi _{{\text{pmk2}}}}} \\ {{\phi _{{\text{pmy}}1}}} \\ {{\phi _{{\text{pmy2}}}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{F_{\text{p}}}/\sigma } \\ { - {F_{\text{p}}}/\sigma } \\ { - {F_{\text{p}}}/\sigma } \\ {{F_{\text{p}}}/\sigma } \end{array}} \right],$

(5)

$\quad\quad\quad\quad\quad\quad\quad\quad \left\{\begin{array}{l}{\phi }_{\text{pmy}1}=\dfrac{{F}_{\text{p}}{R}_{\text{p}}({R}_{\text{c}1}-{R}_{\text{k1}})}{\sigma [{R}_{\text{p}}({R}_{\text{k1}} + {R}_{\text{gb1}} + {R}_{\text{c1}}) + {R}_{\text{k2}}({R}_{\text{gb1}} + {R}_{\text{c1}})]({R}_{\text{p}} + {R}_{\text{gb2}} + {R}_{\text{c1}})},\\ {\phi }_{\text{pmy2}}=\dfrac{{F}_{\text{p}}{R}_{\text{p}}({R}_{\text{c}1}-{R}_{\text{k2}})}{\sigma [{R}_{\text{p}}({R}_{\text{k2}} + {R}_{\text{gb2}} + {R}_{\text{c1}}) + {R}_{\text{k1}}({R}_{\text{gb2}} + {R}_{\text{c1}})]({R}_{\text{p}} + {R}_{\text{gb1}} + {R}_{\text{c1}})},\end{array}\right.$

(6)

式中：σ为永磁体漏磁系数.

2.2.2 X-方向偏置磁路磁通

偏置磁通X−方向的等效磁路如图5所示，φ_pmk5、φ_pmk6为漏磁磁通，φ_pmx1、φ_pmx2为永磁磁极下的偏置磁通；R_gb5、R_gb6为偏置磁通气隙磁阻，R_k5、R_k6为漏磁磁阻. 气隙磁阻及漏磁磁阻分别为

图 5 X−方向等效磁路

Figure 5. Equivalent magnetic circuit along X− direction

下载: 全尺寸图片幻灯片

$\left\{ \begin{array}{l}{R}_{\text{k5}}={R}_{\text{k}6}=\left(2/{m}_{x} + 0.5\right)/\left(0.64{\mu }_{0}{L}_{{\textit{z}}}\right),\\ \left[{R}_{\text{gb}5}\;{R}_{\text{gb}6}\;{R}_{\text{gb}7}\;{R}_{\text{gb}8}\right]=\left[\dfrac{{g}_{0} + x}{{{\mu }_{0}}{A}_{\text{ra}}}\;\dfrac{{g}_{0} + x}{{{\mu }_{0}}{A}_{\text{ra}}}\;\dfrac{{g}_{0}-x}{{{\mu }_{0}}{A}_{\text{ra}}}\;\dfrac{{g}_{0}-x}{{{\mu }_{0}}{A}_{\text{ra}}}\right],\end{array}\right.$

(7)

式中:$ {m_x} $为永磁体端部漏磁半径.

对应图5的磁网络矩阵如式（8）所示，式中各磁阻由图5磁阻模型得出，进而可得偏置磁通如式（9）所示.

$\quad\quad\quad\quad \left[ {\begin{array}{*{20}{c}} {{R_{{\text{gb}}5}} + {R_{\text{p}}} + {R_{{\text{c}}3}}}&{ - {R_{{\text{c}}3}}}&{ - {R_{\text{p}}}}&0 \\ { - {R_{{\text{c}}3}}}&{{R_{{\text{gb6}}}} + {R_{\text{p}}} + {R_{{\text{c}}3}}}&0&{ - {R_{\text{p}}}} \\ { - {R_{\text{p}}}}&0&{{R_{{\text{k5}}}} + {R_{\text{p}}}}&0 \\ 0&{ - {R_{\text{p}}}}&0&{{R_{{\text{gb6}}}} + {R_{\text{p}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\phi _{{\text{pmx1}}}}} \\ {{\phi _{{\text{pmx2}}}}} \\ {{\phi _{{\text{pmk5}}}}} \\ {{\phi _{{\text{pmk6}}}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{F_{\text{p}}}/\sigma } \\ { - {F_{\text{p}}}/\sigma } \\ { - {F_{\text{p}}}/\sigma } \\ {{F_{\text{p}}}/\sigma } \end{array}} \right],$

(8)

$\quad\quad\quad\quad\quad\quad\quad\quad \left\{\begin{array}{l}{\phi }_{\text{pmx}1}=\dfrac{{F}_{\text{p}}{R}_{\text{p}}({R}_{\text{c}3}-{R}_{\text{k5}})}{\sigma [{R}_{\text{p}}({R}_{\text{k5}} + {R}_{\text{gb5}} + {R}_{\text{c3}}) + {R}_{\text{k6}}({R}_{\text{gb5}} + {R}_{\text{c3}})]({R}_{\text{p}} + {R}_{\text{gb6}} + {R}_{\text{c3}})},\\ {\phi }_{\text{pmx2}}=\dfrac{{F}_{\text{p}}{R}_{\text{p}}({R}_{c3}-{R}_{\text{k5}})}{\sigma [{R}_{\text{p}}({R}_{\text{k6}} + {R}_{\text{gb6}} + {R}_{\text{c3}}) + {R}_{\text{k5}}({R}_{\text{gb6}} + {R}_{\text{c3}})]({R}_{\text{p}} + {R}_{\text{gb5}} + {R}_{\text{c3}})}.\end{array}\right.$

(9)

2.2.3 电磁磁路磁通

在图3（b）中，径向X、Y方向的控制磁路总磁阻分别为R_rx、R_ry，如式（19）所示.

$\left\{ \begin{array}{l} {R_{{\text{rx}}}} = {R_{{\text{c}}3}} + {R_{{\text{c}}4}} + 2{R_{\text{a}}} , \\ {R_{{\text{ry}}}} = {R_{{\text{c}}1}} + {R_{{\text{c}}2}}. \end{array} \right.$

(10)

径向X、Y方向的控制磁路磁通φ_cx、φ_cy分别为

$\left\{ \begin{array}{l} {\phi _{{\text{cx}}}} = 2{N_x}{i_x}/{R_{{\text{rx}}}}, \\ {\phi _{{\text{cy}}}} = 2{N_y}{i_y}/{R_{{\text{ry}}}}. \end{array} \right.$

(11)

2.3 悬浮力与刚度计算

由磁网络分析法求解该磁轴承的各磁路磁通，从而可得该磁轴承的刚度和承载力. 将笛卡尔坐标系的正方向作为力的正方向，根据麦克斯韦方程组可推导出径向两自由度的悬浮力分别为

$\left\{ \begin{array}{l} {F_y} = \dfrac{{{{({\phi _{{\text{pmy}}1}} + {\phi _{{\text{cy}}}})}^2} - {{({\phi _{{\text{pmy}}2}} - {\phi ^2_{{\text{cy}}}})}}}}{{2{\mu _0}{A_{\text{r}}}}} + \dfrac{{{\phi ^2_{{\text{pmy}}1}} - {\phi ^2_{{\text{pmy}}3}}}}{{2{\mu _0}{A_{{\text{ra}}}}}}, \\ {F_x} = \dfrac{{{{({\phi _{{\text{pmx}}1}} + {\phi _{{\text{cx}}}})}^2} - {{({\phi _{{\text{pmx}}2}} - {\phi ^2_{{\text{cx}}}})}}}}{{2{\mu _0}{A_{\text{r}}}}} + \dfrac{{{\phi ^2_{{\text{pmx}}1}} - {\phi ^2_{{\text{pmx3}}}}}}{{2{\mu _0}{A_{{\text{ra}}}}}} . \end{array} \right.$

(12)

联立式（9）～（12），求出该磁轴承在Y、X方向两自由度的电流刚度$k_{i_y} $、$k_{i_x} $和位移刚度$k_{d_y} $、$k_{d_x} $分别为

$\left\{ \begin{array}{l} k_{i_y}=\left.\dfrac{\partial F_y}{\partial i_y}\right|_{\substack{x=y=0,\\ i_x=i_y=0}}=\dfrac{2 \phi_{\mathrm{cy}} N}{\sigma \mu_0 g_0\left(R_{\mathrm{c} 2}+R_{\mathrm{c} 1} / R_{\mathrm{ry}}\right)}, \\ k_{i_x}=\left.\dfrac{\partial F_x}{\partial i_x}\right|_{\substack{x=y=0, \\ i_x=i_y=0}}=\dfrac{2 \phi_{\mathrm{cx}} N}{\sigma \mu_0 g_0\left(R_{\mathrm{rx}}+R_{\mathrm{a} 1}\right)}, \\ k_{d_y}=\left.\dfrac{\partial F_y}{\partial y}\right|_{\substack{x=y=0, \\ i_x=i_y=0}}=\dfrac{32\left(\dfrac{A_{\mathrm{ra}} g_0}{\mu_0}-\dfrac{R_{\mathrm{k} 12} R_{\mathrm{k} 11}}{R_{\mathrm{k} 12}+R_{k 11}}\right) S_{\mathrm{p}}^2 l_{\mathrm{m}}}{\sigma \mu_{\mathrm{p}}^2 A_{\mathrm{r}}^2\left(R_{\mathrm{p}} R_{\mathrm{k} 1}+6 R_{\mathrm{k} 2}+R_{\mathrm{p}}\right)^3}, \\ k_{d_x}=\left.\dfrac{\partial F_x}{\partial x}\right|_{\substack{x=y=0,\\ i_x=i_y=0}}=\dfrac{32\left(A_{\mathrm{ra}} g_0 / \mu_0-2 R_{\mathrm{c} 1}\right) S_{\mathrm{p}}^2 l_{\mathrm{m}}}{\sigma \mu_{\mathrm{p}}^2 A_{\mathrm{r}}^2\left(4 R_{\mathrm{c} 2}+2 R_{\mathrm{k} 3} R_{\mathrm{p}}+R_{\mathrm{p}}\right)^3}, \end{array}\right.$

(13)

式中：N为控制绕组的匝数.

3. 仿真分析

为验证磁路计算及悬浮力式的正确性，对新型磁轴承进行有限元仿真. 图1中2种磁轴承结构的主要模型参数如表1所示.

表 1 磁轴承主要参数

Table 1. Main parameters of magnetic bearing

参数	传统	新型
定子外直径/mm	130	139
转子外直径/mm	50	50
轴向长度/mm	45	45
主气隙长度/mm	0.5	0.5
辅助气隙长度/mm		1
永磁体厚度/mm	2	3
永磁体宽度/mm	15	8
磁极截面积/mm²	920	1560
最大控制电流/A	3.6	12/3

下载: 导出CSV

| 显示表格

3.1 磁场分布

新型磁悬浮轴承转子处于平衡位置时，给Y方向线圈通入不同电流时的磁密如图6所示，图中，序号1 ～ 8为图1（a）中对应的磁极位置. 当电流为0 A时，磁通仅由永磁体产生，磁极1、5处的磁密幅值高于磁极3、7处，这是由于引入辅助气隙导致减少了磁极1、5处的永磁体漏磁；随着控制电流的增大，永磁体和控制电流产生的磁通会在磁极1气隙处叠加，使得气隙磁密增大，在磁极5气隙处抵消，使得气隙磁密减小，其他磁极处的磁密基本不变.

图 6 气隙磁密随电流变化

Figure 6. Variation of magnetic density in air gap with current

下载: 全尺寸图片幻灯片

控制线圈未接入电流时，转子在Y方向上不同位置的气隙磁密如图7所示. 图中箭头为气隙增加趋势. 当转子向上偏转时，磁极3气隙处的磁通密度增强，而定子另一侧磁极7气隙处的磁通密度会减弱. 当转子不在平衡位置发生偏离时，此时转子和定子磁极的圆弧面不再平行，使得相同磁极不同位置上的气隙长度有所不同，出现磁极1、2、4、5、6、8处气隙磁通密度不均匀现象.

图 7 气隙磁密随位移变化

Figure 7. Variation of magnetic density in air gap with displacement

下载: 全尺寸图片幻灯片

3.2 电磁力仿真分析

悬浮力和位移刚度、电流刚度为磁悬浮轴承的基本特性，是衡量磁轴承性能的重要指标. 为验证磁路法分析的正确性，由式（13）对HRHMB的刚度特性进行计算并仿真，位移刚度和电流刚度如图8所示.

图 8 径向电流刚度与位移刚度

Figure 8. Radial current stiffness and displacement stiffness

下载: 全尺寸图片幻灯片

从图8（a）中可以看出，铁心材料处于不饱和状态时，转子所受到的悬浮力与绕组电流几乎成正比，同时也与转子在平衡位置附近的位移成正比. 当控制电流在−3 ～ 3 A时，通过磁路法计算出的最大电磁力为472.8 N，仿真得到的电磁力为447.7 N，等效磁路与仿真结果有所偏差，误差为5.6%，这是由于忽略了部分漏磁和定转子铁心中的磁阻造成的，在误差允许范围内.

在图8（b）中，转子位移在（−0.1 ～ + 0.1）mm时，等效磁路法和仿真结果的平均偏差为5.6%. 转子位置在 + 0.1 mm时，磁路法得到的电磁力为134.2 N，仿真的电磁力为123.2 N，偏差为8.9%，这是因为铁心磁通达到饱和引起的，而在高速电机等设备中转子最大允许的偏心位移为气隙长度的1/10，小于达到铁心饱和时的区间，故可忽略铁心饱和的影响. 另外，由于增加了第二气隙，该新型磁轴承X方向的位移刚度和电流刚度小于Y方向. 在应用中可以将Y轴放置于重力方向上，用以增强转子浮动的动态性能，同时，可以通过适当增加X轴方向上的控制线圈的匝数来提高其X方向上的电流刚度.

3.3 两种磁轴承对比

为比较新型和传统偏置磁轴承的悬浮性能及空间利用率，建立相同最大承载力和偏置磁通、转子尺寸等约束条件的仿真模型，传统磁轴承和新型HRHMB的永磁磁密分布如图9所示. 通过有限元仿真计算，在相同承载力下，新型磁轴承为传统偏置磁轴承体积的0.87倍，这表明新型磁轴承结构有着更高的空间利用率.

图 9 2种磁轴承的永磁磁通密度分布

Figure 9. Permanent magnetic flux density distribution of two types of magnetic bearings

下载: 全尺寸图片幻灯片

2种磁悬浮轴承的电流刚度和位移刚度如图10所示. 由图10（a）可以看出，新型HRHMB在X和Y方向上的位移刚度都相对较小，这表明当转子发生偏移时，通过控制系统回到平衡位置较为容易^[20]；图10（b）可以看出，新型HRHMB Y方向上的电流刚度大于传统结构磁轴承的电流刚度，这表明在相同的电流下，新型HRHMB能产生更大的悬浮力. 此外，由于引入辅助气隙，其X方向上的电流刚度会小于Y方向上的电流刚度，因此，新型HRHMB适用于卧式系统等应用场合.

图 10 2种磁轴承的位移刚度和电流刚度对比

Figure 10. Comparison of displacement stiffness and current stiffness of two types of magnetic bearings

下载: 全尺寸图片幻灯片

3.4 电磁力耦合分析

为分析新型异极径向混合磁轴承X、Y方向2个自由度间的电磁力耦合，对转子位移和控制电流对电磁力的影响进行仿真分析，并与图1（b）中的传统八极偏置磁轴承进行对比.

首先，定义转子处于坐标$ (x,y) $时，Y方向电磁力的相对误差为

$\delta (x,y) = \left| {\frac{{{f_Y}(x,y) - {f_Y}(0,y)}}{{{f_Y}(0,y)}}} \right| \times 100{\text{%}}$

(14)

式中：$ {f_Y}(x,y) $为转子在坐标（x，y）时Y方向上的电磁力.

在X轴、Y轴绕组不通电情况下，X和Y方向转子不同位移对Y方向电磁力的影响和电磁力相对误差如图11所示. 由图11（a）可以看出，在x位移量的影响下，新型磁轴承转子在Y方向上受到的电磁力随Y轴位移量仍成正比，变化不明显，最大相对误差为2.76%. 而从图11（b）可以看出，传统偏置磁轴承转子在Y方向上受到的电磁力最大相对误差达到了3.87%.

图 11 转子偏心时的耦合特性

Figure 11. Coupling characteristics with eccentric rotor

下载: 全尺寸图片幻灯片

定义转子位置在（0，y）时，电磁力的相对误差为

$\kappa ({I_X},y) = \left| {\frac{{{f_Y}({I_X},y) - {f_Y}(0,y)}}{{{f_Y}(0,y)}}} \right| \times 100{\text{%}}$

(15)

式中：$ {f_Y}({I_X},y) $和$ {f_Y}(0,y) $分别为转子在（0，y）处通入X方向电流I_X和不通电时Y方向的电磁力.

磁轴承x位移量和Y方向电流为0，通入X方向不同电流时，Y方向电磁力随y的关系曲线和电磁力相对误差如图12所示，由图12（a）可以看出，新型磁轴承Y方向的电磁力与Y轴位移量基本成正比关系，最大误差为3.52%，表明该新型磁轴承X方向电流对Y方向位移刚度影响较小；由图12（b）可以看出，传统偏置磁轴承X方向电流会增加转子在Y方向受的电磁力，其最大相对误差为7.6%.

图 12 转子不偏心时控制电流的耦合特性

Figure 12. Coupling characteristics of control current when rotor is not eccentric

下载: 全尺寸图片幻灯片

当转子在X方向偏心0.1 mm时，图13为磁轴承Y方向电磁力随Y轴位移量和X方向控制电流的关系和电磁力相对误差. 由图13（a）可以看出，在转子偏心条件下通入不同的X方向电流时，Y方向的电磁力仍然正比于Y轴位移量，相对误差最大值为6.5%，这表明转子在X方向发生偏移时，通入X方向电流对Y方向电磁力随Y方向位移的线性特性影响较小，有利于该新型磁悬浮轴承的控制系统设计. 由图13（b）可以看出，通入X方向电流对传统偏置磁轴承转子在Y方向上受的电磁力有较大的影响，最大相对误差为13.6%，其径向两自由度之间的运动存在一定耦合.

图 13 转子偏移0.1mm时控制电流的耦合特性

Figure 13. Coupling characteristics of control current when rotor offsets 0.1 mm

下载: 全尺寸图片幻灯片

综上可知，新型HRHMB在X方向的转子位移基本不影响Y方向的电磁力. 无论转子是否在X方向上发生偏移，X方向上的控制电流对Y方向上的电磁力基本没有影响. 其两自由度之间的电磁力耦合明显小于传统偏置磁轴承，这是由于该新型HRHMB永磁体排列组合方式使其永磁磁路与电磁磁路耦合距离短，且两磁场的相互影响效果较弱，有利于获得良好的解耦特性.

4. 结　论

对于飞轮储能系统等空间利用率要求较高的场合，本文提出了一种新型异极径向混合磁悬浮轴承，该磁轴承电磁磁路和永磁磁路耦合距离短，具有解耦性高的优点. 由等效磁路法建立了该新型轴承的磁路模型，通过解析和仿真研究得出以下结论：

1）由等效磁路法得到的新型HRHMB的承载力、位移刚度和电流刚度的理论结果，与有限元仿真结果相吻合，验证了解析结果的正确性.

2）在相同承载力、偏置磁通、转子尺寸等约束条件下，新型HRHMB结构更为紧凑，空间利用率较高，同时有着较大的电流刚度.

3）新型HRHMB在不同自由度的转子位移和控制电流变化影响下，电磁力相对误差最大值较传统偏置磁轴承小，表明该结构特殊的永磁体排列方式使得电磁磁路与永磁磁路间耦合小，径向二自由度之间电磁力解耦效果较好.

图 1 哨兵影像覆盖图

Figure 1. Sentinel image coverage map

下载: 全尺寸图片幻灯片

图 2 强震动记录台站分布

Figure 2. Distribution of stations recording strong ground motion

下载: 全尺寸图片幻灯片

图 3 地震动参数分布图

Figure 3. Distribution of ground motion parameters

下载: 全尺寸图片幻灯片

图 4 Newmark计算参数分析

Figure 4. Newmark analysis of calculation parameters

下载: 全尺寸图片幻灯片

图 5 Newmark累计位移结果及其滑坡危险性分布

Figure 5. Newmark cumulative displacement and landslide hazard zonation

下载: 全尺寸图片幻灯片

图 6 SBAS-InSAR形变速率

Figure 6. SBAS-InSAR deformation rate

下载: 全尺寸图片幻灯片

图 7 IDNPM滑坡危险性评价结果及野外考察对比（滑坡C-10、31、41、54来自于文献[11]）

Figure 7. IDNPM landslide risk assessment results and field investigation comparison （Landslides C-10, 31, 41, and 54 in the figure are from literature [11]）

下载: 全尺寸图片幻灯片

图 8 IDNPM滑坡危险性评估结果与数值模拟、卫星识别对照

Figure 8. Comparison of landslide risk assessment results of IDNPM with numerical simulations and satellite identification

下载: 全尺寸图片幻灯片

表 1 数据及来源

Table 1. Data and sources

数据名	数据格式	时间	数据源
Sentinel- 1A	5 m × 20 m	2022— 2023年	阿拉斯加卫星设备 https://search.asf.alaska.edu
精密定轨		2022— 2023年	欧洲航天局https://scihub.copernicus.eu
DEM^①	12.5 m	2009年	阿拉斯加卫星设备 https://search.asf.alaska.edu
DEM^②	30 m	2021年	日本宇宙航空研究开发机构https://global.jaxa.jp
地层岩性	1∶500000	2013年	国家地质资料库http://www.ngac.org.cn
卫星影像	10 m	2023年	欧洲航天局https://scihub.copernicus.eu

下载: 导出CSV

表 2 研究区岩性力学参数及取值

Table 2. Mechanical parameters and values of lithology in study area

序号	工程地质单元	c/kPa	$ \varphi $/（°）	$ \gamma $/（kN·m⁻³）
1	坚硬岩	38	40.0	24.5
2	较硬岩	35	35.0	23.0
3	较软岩	30	24.5	17.8
4	较弱岩	25	22.2	17.5
5	松散岩	22	21.0	15.4

下载: 导出CSV

参考文献(28)

[1]	BOMMER J J, RODRı́GUEZ C E. Earthquake-induced landslides in central America[J]. Engineering Geology, 2002, 63(3/4): 189-220.
[2]	JIBSON R W. Regression models for estimating coseismic landslide displacement[J]. Engineering Geology, 2007, 91(2/3/4): 209-218.
[3]	KORONOVSKII N V, ZAKHAROV V S, NAIMARK A A. The unpredictability of strong earthquakes: new understanding and solution of the problem[J]. Moscow University Geology Bulletin, 2021, 76(4): 366-373. doi: 10.3103/S0145875221040074
[4]	许冲,田颖颖,马思远,等. 1920年海原8.5级地震高烈度区滑坡编录与分布规律[J]. 工程地质学报,2018,26(5): 1188-1195. XU Chong, TIAN Yingying, MA Siyuan, et al. Inventory and spatial distribution of landslides in IX-XI high intensity areas of 1920 Haiyuan (China) M8.5 earthquake[J]. Journal of Engineering Geology, 2018, 26(5): 1188-1195.
[5]	XU C, XU X W, YAO X, et al. Three (nearly) complete inventories of landslides triggered by the May 12, 2008 Wenchuan Mw 7.9 earthquake of China and their spatial distribution statistical analysis[J]. Landslides, 2014, 11(3): 441-461. doi: 10.1007/s10346-013-0404-6
[6]	殷跃平. 汶川八级地震地质灾害研究[J]. 工程地质学报,2008,16(4): 433-444. doi: 10.3969/j.issn.1004-9665.2008.04.001 YIN Yueping. Researches on the geo-hazards triggered by Wenchuan earthquake, Sichuan[J]. Journal of Engineering Geology, 2008, 16(4): 433-444. doi: 10.3969/j.issn.1004-9665.2008.04.001
[7]	LAN H X, LI L P, ZHANG Y S, et al. Risk assessment of debris flow in Yushu seismic area in China: a perspective for the reconstruction[J]. Natural Hazards and Earth System Sciences, 2013, 13(11): 2957-2968. doi: 10.5194/nhess-13-2957-2013
[8]	CHEN X L, YU L, WANG M M, et al. Brief communication: landslides triggered by the Ms = 7.0 Lushan earthquake, China[J]. Natural Hazards and Earth System Sciences, 2014, 14(5): 1257-1267. doi: 10.5194/nhess-14-1257-2014
[9]	FAN X M, FANG C Y, DAI L X, et al. Near real time prediction of spatial distribution probability of earthquake-induced landslides-take the Lushan Earthquake on June 1, 2022 as an example[J]. Journal Of Engineering Geology, 2022, 30(3): 729-739.
[10]	甘肃省积石山县志编辑部. 积石山县保安族东乡族撒拉族自治县志 [M]. 北京:红旗出版社,2019:82.
[11]	陈博,宋闯,陈毅,等. 2023年甘肃积石山Ms 6.2地震同震滑坡和建筑物损毁情况应急识别与影响因素研究[J/OL]. 武汉大学学报(信息科学版),(2024-01-02)[2024-01-08], https://doi.org/10.13203/J.whugis20230497. CHEN Bo, SONG Chuang, CHEN Yi, et al. Emergency identification and influencing factor analysis of coseismic landslides and building damages induced by the 2023 Ms 6.2 Jishishan (Gansu, China) Earthquake[J/OL]. Geomatics and Information Science of Wuhan University, (2024-01-02)[2024-01-08], https://doi.org/10.13203/J.whugis20230497.
[12]	ZHUO Y Q, LIU P X, GUO Y S, et al. Cross-effects of loading rate and cumulative fault slip on pre-seismic rupture and unstable slip rate of laboratory earthquakes[J]. Tectonophysics, 2022, 826: 229266.1-229266.13.
[13]	ZENG Y, ZHANG Y B, LIU J, et al. Rapid emergency response assessment of earthquake-induced landslides driven by fusion of InSAR deformation data and newmark physical models[J]. Remote Sensing, 2023, 15(18): 4605.1-4605.21.
[14]	ZENG Y, ZHANG Y B, LIU J, et al. Assessment of earthquake-induced landslide hazard zoning using the physics-environmental coupled model[J]. Journal of Mountain Science, 2023, 20(9): 2644-2664. doi: 10.1007/s11629-023-7947-3
[15]	ZHANG Y B, LI X Q, FU J W, et al. Quick analysis model for earthquake-induced landslide movement based on energy conservation[J]. Landslides, 2024, 21(6): 1315-1331. doi: 10.1007/s10346-024-02223-x
[16]	王卫民,何建坤,丁林,等. 甘肃临夏Ms6.2级地震震源破裂过程反演初步结果[EB/OL]. (2023-12-19) [2024-01-08]. https://mp.weixin.qq.com/s/SuKbTgq8-0EV2FbNcdlYqQ. WANG Weimin,HE Jiankun,DING Lin,et al. Preliminary results of source rupture process inversion for Ms6.2 earthquake in Linxia,Gansu Province [EB/OL]. (2023-12-19) [2024-01-08]. https://mp.weixin.qq.com/s/SuKbTgq8-0EV2FbNcdlYqQ[标题&网址验证OK. 点击这里,可链出确认]
[17]	王勤彩,罗钧,陈翰林,等. 2023年12月18日甘肃积石山6.2级地震震源机制解[J]. 地震,2024,1(1): 185-188. WANG Qincai, LUO Jun, CHEN Hanlin, et al. Focal mechanism for the December 18, 2023, Jishishan Ms6.2 earthquake in Gansu Province[J]. Earthquake, 2024, 1(1): 185-188.
[18]	袁道阳,张培震,雷中生,等. 青海拉脊山断裂带新活动特征的初步研究[J]. 中国地震,2005,21(1): 93-102. YUAN Daoyang, ZHANG Peizhen, LEI Zhongsheng, et al. A preliminary study on the new activity features of the lajishan mountain fault zone in Qinghai Province[J]. Earthquake Research in China, 2005, 21(1): 93-102.
[19]	张克信,何卫红,徐亚东,等. 论从俯冲增生杂岩带重建洋板块地层主要类型与序列:以青藏特提斯二叠系为例[J]. 沉积与特提斯地质,2021,41(2):137-151. ZHANG Kexin, HE Weihong, XU Yadong, et al. , Reconstruction of main types for oceanic plate strata in the subduction accretionary complex and feature of sequence for each type: an example from the Qinghai—Tibet Tethyan Permian strata[J]. Sedimentary Geology and Tethyan Geology, 2021, 41(2):137-151.
[20]	王运生,赵波,吉锋,等. 2023年甘肃积石山Ms6.2级地震震害异常的启示[J]. 成都理工大学学报(自然科学版),2024,51(1): 1-8. WANG Yunsheng, ZHAO Bo, JI Feng, et al. Preliminary insights into the hazards triggered by the 2023 Jishishan Ms 6.2 earthquake in Gansu Province[J]. Journal of Chengdu University of Technology (Science & Technology Edition), 2024, 51(1): 1-8.
[21]	程小杰,杨为民,向灵芝,等. 基于Newmark模型的天水市北山地震黄土滑坡危险性评价[J]. 地质力学学报,2017,23(2): 296-305. CHENG Xiaojie, YANG Weimin, XIANG Lingzhi, et al. Risk assessment of seismic loess landslide based on newmark model in Beishan, Tianshui city[J]. Journal of Geomechanics, 2017, 23(2): 296-305.
[22]	李雪婧,高孟潭,徐伟进. 基于Newmark模型的概率地震滑坡危险性分析方法研究——以甘肃天水地区为例[J]. 地震学报,2019,41(6):795-807. LI Xuejing, GAO Mengtan, XU Weijin. Probabilistic seismic slope displacement hazard analysis based on Newmark displacement model: Take the area of Tianshui, Gansu Province, China as an example[J]. Acta Seismologica Sinica, 2019, 41(6):795-807.
[23]	DREYFUS D, RATHJE E M, JIBSON R W. The influence of different simplified sliding-block models and input parameters on regional predictions of seismic landslides triggered by the Northridge earthquake[J]. Engineering Geology, 2013, 163: 41-54. doi: 10.1016/j.enggeo.2013.05.015
[24]	ZHANG Y B, LIU J, CHENG Q G, et al. A new permanent displacement model considering pulse-like ground motions and its application in landslide hazard assessment[J]. Soil Dynamics and Earthquake Engineering, 2022, 163: 107556.1-107556.27.
[25]	LIU J, ZHANG Y B, WEI J T, et al. Hazard assessment of earthquake-induced landslides by using permanent displacement model considering near-fault pulse-like ground motions[J]. Bulletin of Engineering Geology and the Environment, 2021, 80(11): 8503-8518. doi: 10.1007/s10064-021-02464-3
[26]	SAMSONOV S V, TRISHCHENKO A P, TIAMPO K, et al. Removal of systematic seasonal atmospheric signal from interferometric synthetic aperture radar ground deformation time series[J]. Geophysical Research Letters, 2014, 41(17): 6123-6130. doi: 10.1002/2014GL061307
[27]	CHANG M, SUN W J, XU H Z, et al. Identification and deformation analysis of potential landslides after the Jiuzhaigou earthquake by SBAS-InSAR[J]. Environmental Science and Pollution Research, 2023, 30(13): 39093-39106. doi: 10.1007/s11356-022-25055-5
[28]	DAI C, LI W L, LU H Y, et al. Landslide hazard assessment method considering the deformation factor: a case study of Zhouqu, Gansu Province, Northwest China[J]. Remote Sensing, 2023, 15(3): 596.1-596.24.