结构拟静力分量对中低速磁浮车-桥系统地震响应的影响
上海交通大学 船舶海洋与建筑工程学院,上海 200240
Influence of Structural Pseudo-Static Components on Seismic Responses of Low-Medium Speed Maglev Vehicle-Bridge System
School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
通讯作者: 滕念管,男,副教授,电话 (Tel.) 025-34206334; E-mail:ngteng@sjtu.edu.cn.
责任编辑: 陈晓燕
收稿日期: 2020-11-16
基金资助: |
|
Received: 2020-11-16
作者简介 About authors
黄凤华(1991-),男,江苏省南通市人,博士生,主要从事车桥耦合振动研究.
与传统桥梁抗震分析关注系统的相对运动不同,车-桥系统的抗震分析关注车辆系统的绝对运动,需考虑结构拟静力分量的影响.为探讨结构拟静力分量对中低速磁浮车-桥系统地震响应特性的影响,以某典型中低速磁浮线路为研究对象,考虑基于比例积分微分 (Proportional Integral Derivative, PID) 主动悬浮控制的磁浮车-桥耦合关系,采用相对运动法和直接求解法,分别在相对坐标系和绝对坐标系下建立中低速磁浮车-桥系统地震响应分析模型.在此基础上,重点探讨了结构拟静力分量对磁浮间隙、车辆系统以及桥梁结构动力响应的影响规律.结果表明:结构拟静力分量对车辆系统的动力响应影响较大,忽略结构拟静力分量会严重低估车辆的加速度响应,差异值最大可达447%;结构拟静力分量对桥梁结构位移响应的影响有限;建议采用考虑结构拟静力分量的绝对位移法处理磁浮车-桥系统中的地震输入.
关键词:
Unlike the traditional seismic analysis on bridge engineering which mainly concerns the relative motion, the seismic analysis of vehicle-bridge system focuses more on the absolute motion of the vehicle system, where the pseudo-static components need to be considered. In order to study the effect of structural pseudo-static components on the seismic responses of low-medium speed maglev vehicle-bridge system, the discussion relying on a typical low-medium speed maglev line was implemented. Considering the vehicle-bridge coupling relationship based on the proportional integral derivative (PID) active suspension control, two dynamic models for the low-medium speed maglev vehicle-bridge system subjected to earthquake were relatively established in the relative coordinate and the absolute coordinate, and the seismic forces were correspondingly solved by employing the relative motion method and the direct solution method. A case study was subsequently conducted, where the effect of structural pseudo-static components on dynamic responses of the vehicle-bridge system, including the maglev gap, the vehicle system and the bridge was further analyzed. The results show that the influence of pseudo-static components on the dynamic responses of the vehicle system is quite significant. Ignoring the pseudo-static components could seriously undervalue the dynamic responses of the vehicle system, of which the difference could be as high as 447%. Compared to the vehicle system, the influence of pseudo-static components on the vertical and transversal displacement responses of the bridge is fairly limited. Therefore, it is suggested that the absolute displacement method which takes the pseudo-static components into account should be adopted to deal with the seismic excitation in the maglev vehicle-bridge system.
Keywords:
本文引用格式
黄凤华, 程斌, 滕念管.
HUANG Fenghua, CHENG Bin, TENG Nianguan.
中低速磁浮交通作为一种新兴的地面交通运输方式,以其噪声低、污染小、爬坡能力强等优点受到广泛关注,目前,国内已建成多条中低速磁浮线路并成功实现商业运营.同时,我国是地震多发国,随着磁浮线路的不断扩展,磁浮列车在桥上运行时遭遇地震的可能性越来越高,地震作用不仅会引起桥梁结构的破坏,更会影响到桥上列车的行车安全.
由于磁浮交通起步较晚,现有研究多集中于磁浮车-桥系统的耦合振动,鲜有关于磁浮车-桥系统地震响应的研究.耿杰等[1,2]针对长沙中低速磁浮交通运营线路展开现场测试,揭示了列车-简支梁系统的耦合振动特性.Hu等[3,4,5,6]学者建立了基于PID主动悬浮控制的中低速磁浮列车-桥梁系统耦合振动模型,对车辆、桥梁系统的动力响应特性展开数值模拟.相比于磁浮交通,在轮轨交通中已有不少学者展开车-桥系统地震响应特性的研究.其中,车-桥-地震系统模型建立方法主要分2类,一类是建立在相对坐标系中,忽略结构拟静力分量的影响,只考虑地震力引起的动力分量,Xia等[7,8,9,10,11,12]采用相对运动法处理地震力输入,忽略了结构拟静力分量的影响,对地震作用下车-桥系统的动力响应及桥上列车的行车安全展开研究.另一类是建立在绝对坐标系中,同时考虑地震力引起的拟静力分量和动力分量,杜宪亭等[13,14,15]在绝对坐标系和相对坐标系中分别建立车-桥-地震系统模型,深入探讨了结构拟静力分量对车-桥-地震系统动力响应的影响,结果表明,忽略结构拟静力分量会造成车-桥-地震系统动力分析结果出现较大偏差.事实上,与传统结构抗震分析关注相对运动不同,车-桥系统的抗震分析同时关注车辆的绝对运动,需要考虑拟静力分量的影响,因此,研究拟静力分量的影响对建立可靠的车-桥-地震系统模型具有重要意义.
鉴于磁浮交通的车桥耦合机理与轮轨交通有较大区别,本文考虑基于PID主动悬浮控制的磁浮车-桥耦合关系,采用相对运动法和直接求解法,分别在相对坐标系和绝对坐标系下建立中低速磁浮车-桥系统地震响应分析模型.在此基础上,给出车-桥系统地震响应的求解步骤,并通过算例重点探讨结构拟静力分量对磁浮间隙、车辆系统、桥梁结构动力响应的影响规律,以期为地震作用下磁浮车-桥系统动力响应研究以及桥上列车行车安全性评估提供理论支撑.
1 磁浮车-桥-地震系统模型
1.1 车辆模型
选取典型的5模块中低速磁浮列车作为研究对象,每节车体由5个悬浮架模块支撑,每个悬浮架模块包含左、右侧两套悬浮架及防侧滚梁,电磁铁固定于悬浮架底部,车体与悬浮架之间通过二系悬挂系统相连.单节编组列车共有20个二系悬挂、40个电磁铁.车辆模型如图1所示,图中,lm1x1为第1个悬浮架上第1个电磁铁该悬浮架质心的纵向距离,lc1x1、lb1x1分别为第1个悬浮架上第1个二系悬挂到车体质心、该悬浮架质心的纵向距离.车体和悬浮架均视为刚体,不考虑纵向运动,二系悬挂系统则简化为弹簧阻尼.车体考虑沉浮(Yc)、点头(αc)、横移(Zc)、摇头(βc)、侧滚(γc) 5个自由度.考虑到防侧滚梁使得左右两侧悬浮架的运动解耦,且抑制了悬浮架的侧滚,因此,悬浮架考虑沉浮(Yb)、点头(αb)、横移(Zb)、摇头(βb) 4个自由度.单节编组列车模型总计45个自由度.
图1
车体和悬浮架的受力状态如图2所示,图中αc、βc、γc分别为车体的点头、摇头、侧滚位移;αb、βb分别为悬浮架的点头、摇头位移.由于地震不直接作用于车辆系统,而是通过磁轨关系影响车辆系统,所以地震作用下的车辆运动方程无需特别考虑地震力.基于达朗贝尔原理,建立如下车辆系统运动方程:
图2
(1) 车体垂向运动.
(2) 车体横向运动.
(3) 车体侧滚运动.
(4) 第j个悬浮架模块左侧悬浮架.
(5) 第j个悬浮架模块右侧悬浮架.
式中:
1.2 桥梁模型
基于有限元法建立桥梁系统的动力学模型,采用一致质量矩阵,得到车辆荷载和地震作用下桥梁结构在绝对坐标系中的分块运动方程:
式中:M、C、K分别为质量矩阵、阻尼矩阵和刚度矩阵,下标ss为桥梁非支撑节点自由度项,下标gg为支撑节点自由度项,下标sg和gs为支撑节点和非支撑节点之间的耦合项;u、
展开第一行,消去未知项Pgg,并忽略方程右端的阻尼项,得:
以下分别采用相对运动法和直接求解法对上述运动方程进行求解.
(1) 相对运动法.
相对运动法 (Relative Motion Method, RMM) 在相对坐标系中进行求解,是传统桥梁抗震分析中常用的求解方法,该方法将结构绝对位移u分解为拟静力位移us和相对动力位移ud,即
将式 (4) 代入式 (3),得到相对坐标系下忽略结构拟静力分量的桥梁系统运动方程:
式中:R为影响矩阵,R=-
采用振型叠加法对式 (5) 进行解耦:
式中:
将式 (6) 代入式 (5),引入瑞雷比例阻尼,基于振型正交特性,得到相对坐标系下桥梁的模态平衡方程:
式中:
联立式 (6)、(7) 和车辆运动方程即可得到相对坐标系下忽略结构拟静力分量时车-桥耦合系统的地震响应.
(2) 直接求解法.
直接求解法 (Direct Solution Method, DSM) 在绝对坐标系中进行求解,即直接对式 (3) 进行求解.同样采用振型叠加法对运动方程进行解耦:
式中:
将式 (8) 代入式 (3),得到绝对坐标系下桥梁系统的模态平衡方程:
可以看出,相比于相对运动法,直接求解法考虑了结构拟静力分量的影响,需同时输入地震加速度时程和地震位移时程.联立式 (8)、(9) 和车辆运动方程即可得到绝对坐标系下考虑结构拟静力分量时磁浮车-桥耦合系统的地震响应.
1.3 控制系统模型
中低速磁浮车辆通过U形悬浮电磁铁与F轨相互吸引产生的电磁力实现车辆的悬浮、导向功能.根据电磁铁二维受力分析[16],磁铁模块在发生横向偏移时电磁悬浮力和导向力的计算公式如下:
式中:fmy和fmz分别为电磁悬浮力和导向力;Am为磁极面积;μ0为空气磁导率;h为悬浮间隙;Wm为磁极宽度;c为导向间隙;Nconst为恒流线圈数;Iconst为恒电流,用于平衡车辆自重;Nc控制电流线圈数;ic为控制电流.
可见,磁浮间隙是计算电磁悬浮力和导向力的重点,也是准确预测车-桥系统动力响应的关键.图3给出了相对坐标系和绝对坐标系下磁浮间隙的计算示意图,其中
图3
由于中低速磁浮车的车速较低 (通常不高于 120 km/h),列车运行时所需导向力较小,所以车辆采用悬浮主动控制、导向被动控制的悬浮导向策略.侧向上,通过电磁力横向分力促使车辆自复位以实现导向功能;垂向上,引入主动控制器使悬浮间隙稳定在一定误差范围内,即在控制系统获取悬浮间隙信号后,通过控制算法产生与间隙扰动成一定关系的控制电流,调整电磁悬浮力,保证悬浮间隙在容许范围内波动.引入PID 控制器建立悬浮控制系统,基于悬浮间隙波动Δh和悬浮间隙变化率
式中:KP、KI、KD分别为比例参数、积分参数、微分参数;k为控制时间.
1.4 数值求解流程
在求解磁浮车-桥系统地震响应时,分别建立车辆、桥梁、控制器3个子系统,地震激励直接作用于桥墩支撑处,通过车-桥相互作用关系影响车辆系统,继而影响整个车-桥耦合系统.基于上述模型,编制计算程序,具体步骤如下:
(1) 输入车辆参数、桥梁结构参数、轨道不平顺、控制系统参数以及地震波信息,建立车辆、桥梁系统运动方程及悬浮控制系统模型,设置初始条件.
(2) 确定列车在梁上的位置,根据上一时刻车辆、桥梁的振动状态确定磁浮间隙及其变化率,通过控制算法式 (11) 输出控制电流,基于电磁力计算式 (10) 得到悬浮力和导向力,并作为当前时刻的电磁力作用于车辆上,然后通过四阶龙格库塔法求解车辆运动方程得到当前时刻车辆的振动状态.
(3) 计算桥梁运动方程右端的地震荷载项,叠加上一步骤的电磁悬浮力和导向力,通过龙格库塔法求解桥梁的运动方程得到当前时刻桥梁的振动状态.
(4) 保存当前时刻车辆、桥梁的振动状态,并作为下一时刻车-桥-地震系统的输入条件.
(5) 判断列车是否完全驶离桥梁,如果“是”,则结束计算;如果“否”,则重复步骤 (2) ~ (4),直至列车完全驶离桥梁.
2 模型参数及验证
2.1 计算参数
图4
表1 磁浮车辆主要计算参数
Tab.1
符号 | 名称及单位 | 数值 |
---|---|---|
Mc | 车体质量/kg | 2.0×104 |
Mb | 悬浮侧架质量/kg | 1000 |
Jαc | 车体转动惯量 (点头)/(kg·m2) | 3.85×105 |
Jβc | 车体转动惯量 (摇头)/(kg·m2) | 3.88×105 |
Jγc | 车体转动惯量 (侧滚)/(kg·m2) | 2.21×103 |
Jαb | 悬浮架转动惯量 (点头)/(kg·m2) | 1150 |
Jβb | 悬浮架转动惯量 (摇头)/(kg·m2) | 1200 |
Ksy | 二系悬挂系竖向刚度/(N·m-1) | 8.0×104 |
Csy | 二系悬挂系竖向阻尼/(N·s·m-1) | 5.0×103 |
Ksz | 二系悬挂系横向刚度/(N·m-1) | 1.5×105 |
Csz | 二系悬挂系横向阻尼/(N·s·m-1) | 2.25×103 |
图5
2.2 仿真结果验证
鉴于目前缺乏地震作用下磁浮车-桥耦合系统动力响应的实测数据,分别从车-桥耦合振动响应和桥梁地震响应两个方面出发,对所建立的磁浮车-桥-地震系统模型进行验证,以桥梁的竖向位移时程响应为例.首先,只考虑车辆荷载,将程序计算得到的桥梁位移时程与文献[19]的现场测试结果进行对比.其次,只考虑地震激励,将桥梁位移时程的程序计算值与商业软件SAP2000 计算值进行对比.对比结果如图6所示,图中:yby为轨道梁跨中竖向位移.可以看出,桥梁振动位移数值模拟波形与实测波形有所差异,这是因为桥梁动位移的波形受轨道不平顺的影响较大[20],而数值模拟中的轨道不平顺与现场测试中的轨道不平顺难以保持一致,此外,实测值易受到外界因素的干扰.总体而言,桥梁在车辆荷载下的位移响应仿真结果与实测结果较为吻合,都经历了车辆上桥、桥上运行和离桥3个阶段,仿真和实测得到的最大动位移误差在5%以内.基于RMM和DSM分别得到的桥梁地震响应仿真结果与相应的SAP2000计算结果几乎完全重合.可见,本文所建立的中低速磁浮车-桥-地震系统数值模型是可靠的.
图6
3 拟静力分量对系统响应的影响
3.1 磁浮间隙和电磁力
图7~8给出了车辆以80 km/h速度过桥时同一个电磁铁处磁浮间隙和电磁力的时程曲线,可以看出:
图7
图8
(1) 考虑拟静力分量时悬浮间隙时程曲线的波形与不考虑时较相似,但幅值有所差异.不考虑拟静力分量时悬浮间隙的最大波动幅值为1.09 mm,考虑后悬浮间隙的最大波动幅值为1.27 mm,后者比前者高16%.拟静力分量对悬浮电磁力的影响规律与悬浮间隙的较一致,考虑拟静力分量时悬浮电磁力的最大波动幅值较不考虑时大6.3%.
(2) 拟静力分量改变了导向间隙及导向电磁力时程曲线的波形,且显著增大了曲线幅值.不考虑拟静力分量时导向间隙、导向电磁力的最大波动幅值分别为7.55 mm、1.10 kN,而考虑后最大波动幅值分别为19.14 mm、2.32 kN,后者比前者分别高153%和110%.
图9
图10
图10
电磁力幅值随车速的变化
Fig.10
Variation of electromagnetic force amplitude with vehicle speed
表2 考虑和忽略拟静力分量时车桥系统响应幅值的差异值
Tab.2
车速/ (km·h-1) | 悬浮系统 | 车辆系统 | 轨道梁跨中 | |||||||
---|---|---|---|---|---|---|---|---|---|---|
悬浮间隙 波动幅值 | 导向间隙 | 悬浮电磁力 波动幅值 | 导向电磁力 | 沉浮加速度 | 偏航加速度 | 竖向位移 | 横向位移 | |||
60 | 19.2 | 111.9 | -5.7 | 143.2 | 368.3 | 61.1 | 3.7 | -1.5 | ||
70 | 23.4 | 216.5 | 7.3 | 118.1 | 394.1 | 61.2 | 0.8 | 0.2 | ||
80 | 16.3 | 153.6 | 6.3 | 110.6 | 449.1 | 50.5 | 2.1 | 1.2 | ||
90 | 7.3 | 219.1 | 4.7 | 113.8 | 322.9 | 44.3 | 2.2 | 0.2 | ||
100 | 22.5 | 192.4 | -5.1 | 82.3 | 323.0 | 46.7 | -0.3 | 2.0 | ||
110 | 26.7 | 205.8 | 9.2 | 89.0 | 309.1 | 56.2 | -0.4 | 1.0 | ||
120 | 34.8 | 171.2 | 0.97 | 101.4 | 366.9 | 47.2 | -0.2 | -1.5 |
注:差异值=(考虑拟静力时系统幅值-忽略拟静力时系统幅值)/忽略拟静力时系统幅值
(1) 考虑拟静力分量后悬浮间隙有所增大,在计算车速范围内,忽略拟静力分量会导致悬浮间隙的差异值为7.3%~34.8%,而拟静力分量对悬浮电磁力的影响较小,最大差异值为9.23%.拟静力分量显著增大了导向间隙和导向电磁力,忽略拟静力分量会使得导向间隙和导向电磁力被严重低估,差异值分别为111%~219%和82%~143%.
(2) 相比对导向间隙、导向电磁力的影响,拟静力分量对悬浮间隙、悬浮电磁力的影响较小,原因与磁浮车辆悬浮主动控制、导向被动控制的策略有关,考虑拟静力分量后,桥梁动力响应加大 (见图6(b)~6(c) ),致使车辆与轨道梁之间的间隙显著增大,导向间隙增大,然而悬浮方向采用主动控制策略,悬浮间隙被快速调节并稳定在额定范围内,波动幅度较小.
3.2 车辆系统
图11所示为车辆以80 km/h速度过桥时第2节车体质心的加速度时程曲线.可见,考虑拟静力分量后车辆沉浮、横移加速度时程曲线波形、幅值与不考虑时差别较大,不考虑拟静力分量时车辆沉浮、横移最大加速度分别为0.13、1.67 m/s2,考虑后加速度明显增大,最大加速度分别为0.72、2.52 m/s2,显然忽略拟静力分量会严重低估车辆加速度响应.
图11
图11
车体质心处加速度时程曲线
Fig.11
Acceleration time-history curves at centroid of vehicle
图12
(1) 考虑结构拟静力分量后,车辆的沉浮、横移加速度显著增大,在车速计算范围内,比不考虑拟静力分量时分别大309%~449%和44%~61%,意味着忽略拟静力分量将严重低估桥上列车的动力响应,继而可能造成地震作用下列车行车安全的误判.
(2) 比较而言,拟静力分量对车辆沉浮加速度的影响要明显大于对车辆横向加速度的影响, 原因同样与悬浮导向控制策略有关.考虑拟静力分量后桥梁动力响应增大,为保证车辆与轨道梁间的悬浮间隙在额定范围内,在控制器作用下车辆短时间内经历较大的沉浮运动,致使车辆沉浮加速度显著增大.相比于垂向的主动控制,横向靠二系悬挂系统提供阻尼,属于被动控制,调节速率相对缓慢,因而横移加速度的增大趋势比沉浮加速度的小.
3.3 桥梁结构
图13所示为车辆以80 km/h速度过桥时第2跨简支梁的跨中位移时程曲线,图中ybz为轨道梁跨中横向位移.考虑到桥梁结构关注相对运动,此处的桥梁跨中位移取相对桥墩支撑运动的动力响应.可以看出:考虑拟静力分量时桥梁跨中的位移时程变化规律与不考虑时较一致,但幅值有细微差别,不考虑拟静力分量时桥梁跨中竖向、横向最大位移为3.02、10.42 mm,考虑后最大位移分别为3.08、10.60 mm,两种情况下的位移响应差别较小.
图13
图14
图14
桥梁跨中位移随车速的变化
Fig.14
Variation of displacement at mid-span of bridge with vehicle speed
4 结论
基于磁浮车-桥动力相互作用理论,建立了不同坐标系下中低速磁浮车-桥系统地震响应的分析模型,深入研究了结构拟静力分量对车-桥系统地震响应特性的影响规律,得到如下结论:
(1) 考虑结构拟静力分量后,悬浮间隙及悬浮电磁力均有所增大,但增大趋势不明显,而导向间隙和导向电磁力则显著增大,原因与磁浮列车悬浮主动控制、导向被动控制的策略有关.
(2) 结构拟静力分量对车辆系统的动力响应影响最大,将显著增大车辆的沉浮、横移加速度,忽略结构拟静力分量会严重低估桥上列车的动力响应,继而造成地震作用下桥上列车行车安全的误判.
(3) 相比于车辆系统,结构拟静力分量对桥梁系统动力响应的影响有限,因此,当仅需计算桥梁结构在地震激励和车辆荷载同时作用下的动力响应时,可以采用传统结构抗震分析方法RMM,即忽略结构拟静力分量的影响.
(4) 为准确合理地预测地震作用下中低速磁浮车-桥系统的动力响应,建议采用考虑结构拟静力分量的绝对位移法 (DSM) 处理车-桥-地震系统的地震输入,此时需同时输入地震波位移时程和加速度时程.
参考文献
中低速磁浮列车-简支梁系统耦合振动试验研究
[J].
Experimental study on coupled vibration of low-medium speed maglev train and simply supported girder system
[J].
Dynamic ana-lysis of the interactions between a low-to-medium-speed maglev train and a bridge: Field test results of two typical bridges
[J].DOI:10.1177/0954409718758502 URL [本文引用: 1]
Coupled dynamic analysis of low and medium speed maglev vehicle-bridge interaction using SIMPACK
[J].DOI:10.1177/0954409720925676 URL [本文引用: 1]
A numerical model to predict three-dimensional interaction dynamic of low-medium-speed maglev vehicle-guideway bridge system
[J].DOI:10.1177/0954409721996553 URL [本文引用: 1]
Coupled vibration analysis of maglev vehicle-guideway while standing still or moving at low speeds
[J].DOI:10.1080/00423114.2015.1013039 URL [本文引用: 1]
Dynamic interaction analysis of maglev-guideway system based on a 3D full vehicle model
[J].DOI:10.1142/S0219455417500067 URL [本文引用: 1]
Dynamic analysis of train-bridge system subjected to non-uniform seismic excitations
[J].
地震作用下高速列车-线路-桥梁系统动力响应
[J].
Dynamic responses of high-speed train-track-bridge system under seismic excitations
[J].
近场地震作用下高速铁路简支梁桥行车安全性影响研究
[J].
Research on driving impact on high-speed railway simply supported beam bridges under near-field earthquakes
[J].
Seismic response analysis of an interacting curved bridge-train system under frequent earthquakes
[J].
Dynamic vehicle-bridge interaction under simultaneous vertical earthquake excitation
[J].DOI:10.1007/s10518-016-9954-z URL [本文引用: 1]
Nonlinear dynamic response analysis of vehicle-bridge interactive system under strong earthquakes
[J].DOI:10.1016/j.engstruct.2018.09.014 URL [本文引用: 1]
车桥耦合动力分析中地震动输入模式的研究
[J].
Study on the input patterns of seismic ground motion in the dynamic interaction analysis of train-bridge system
[J].
地震作用下结构拟静力分量对车桥系统动力响应的影响分析
[J].
Influence of structure quasi-static components on the dynamic response of train-bridge systems during earthquakes
[J].
拟静力分量对列车-轨道-桥梁系统地震响应的影响
[J].
Effects of structural quasi-static components on seismic responses of train-track-bridge system
[J].
Lift and side forces on rectangular pole pieces in two dimensions
[J].DOI:10.1063/1.1663505 URL [本文引用: 1]
低速磁浮轨道不平顺功率谱研究
[J].
Estimation of power spectrum density track irregularities of low-speed maglev railway lines
[J].
F轨对中低速磁浮列车-桥梁系统竖向耦合振动的影响研究
[J].
Study on the influence of F-rail in vertical coupling vibration of low-medium speed maglev train-bridge system
[J].
Dynamic interaction analysis of the high-speed maglev vehicle/guideway system based on a field measurement and model updating method
[J].DOI:10.1016/j.engstruct.2018.11.031 URL [本文引用: 1]
/
〈 |
|
〉 |
![](../../../richhtml/1006-2467/richHtml_jats1_1/images/1.png)