• ISSN 1673-5722
  • CN 11-5429/P

基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析

孙柏涛 姜琪 闫培雷

王宇亮, 崔洪军, 张玉敏, 谷玉珍. 滑动长孔螺栓摩擦阻尼器力学性能试验研究[J]. 震灾防御技术, 2020, 15(1): 11-20. doi: 10.11899/zzfy20200102
引用本文: 孙柏涛, 姜琪, 闫培雷. 基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析[J]. 震灾防御技术, 2018, 13(4): 869-877. doi: 10.11899/zzfy20180414
Wang Yuliang, Cui Hongjun, Zhang Yumin, Gu Yuzhen. Experimental Research on Mechanical Property of Friction Damper with Sliding Long-hole Bolt[J]. Technology for Earthquake Disaster Prevention, 2020, 15(1): 11-20. doi: 10.11899/zzfy20200102
Citation: Sun Baitao, Jiang Qi, Yan Peilei. Optimization Design and Mechanical Property Analysis of Small Steel Frame Reaction-wall Based on Large Scale Model Structure Test[J]. Technology for Earthquake Disaster Prevention, 2018, 13(4): 869-877. doi: 10.11899/zzfy20180414

基于小比例缩尺模型结构试验的小型钢架反力墙优化设计与力学性能分析

doi: 10.11899/zzfy20180414
基金项目: 

黑龙江省应用技术研究与开发计划 国家项目省级资助GX16C009

国家自然科学青年基金 51608495

中国地震局创新团队发展计划 中国大陆地区地震灾害模拟与评估

详细信息
    作者简介:

    孙柏涛, 男, 生于1961年。研究员。主要从事地震灾害评估等方面的研究。E-mail:sunbt@iem.cn

    通讯作者:

    姜琪, 男, 生于1993年。硕士研究生。主要从事结构试验方面的研究。E-mail:jiangqi_hit@163.com

Optimization Design and Mechanical Property Analysis of Small Steel Frame Reaction-wall Based on Large Scale Model Structure Test

  • 摘要: 为了降低模型制作加工成本和节约试验空间,开展结构缩尺比例小于1:10的小模型结构抗震试验。目前,国内相关的试验设备较少,本文将参考大型反力墙的结构形式,研究小型钢架反力墙,并对其进行优化设计和力学性能分析,为后续的制作加工提供理论依据。
  • 预制装配式剪力墙结构中存在大量接缝,接缝的连接性能对剪力墙构件乃至整个结构的影响非常大。目前我国装配式剪力墙结构的设计目的可被形容为:希望把装配式剪力墙结构制作成等同于整体剪力墙结构行为。若不采用上述等同于现浇的设计理念,可考虑在竖缝中加设摩擦阻尼器进行连接,荷载较小时,阻尼器提供一定刚度,预制墙与现浇墙具有几乎相同的性质;荷载较大时,利用竖缝中的阻尼器实现耗能,提高延性,实现较好的抗震性能。

    国内外学者对摩擦阻尼器进行了大量研究工作,Gerald等(1989)研制了用于单斜拉压支撑上的长孔螺栓节点,Douglas和John等开发了用于替代钢结构梁柱节点焊接从而提高节点延性的SBCS节点和用于倒Ⅴ形支撑能产生转动滑动的SBCS节点。Martinez-Rueda(2000)针对一般交叉支撑或倒Ⅴ形支撑占据使用空间的缺点,设计出专门安装于折形支撑的转动摩擦耗能节点。Morgen等(2004)为无黏结后张拉预应力钢筋混凝土结构中梁柱节点连接专门设计了摩擦转动耗能节点。Pall等(1982)设计了安装在X形支撑交叉处的Pall摩擦耗能器。刘季等(1996)周云等(19961997)对2种不同构造Pall摩擦耗能器(十字型、交叉型)及摩擦耗能支撑进行多种工况下的低周反复荷载试验,研究其耗能性能,研究结果表明,摩擦耗能支撑具有稳定的滞回性能和极好的耗能能力,摩擦片类型与性能对摩擦耗能支撑的性能有较大影响。在试验的基础上对2种摩擦耗能器提出改进方案,设计了双摩擦耗能器。吴斌等(199719992003a2003b)通过对Pall摩擦耗能器的机构变形分析,建立摩擦耗能器芯板螺栓槽孔设计方法,并在Pall摩擦耗能器的基础上改进出构造更简单的T形芯板摩擦耗能器。Mualla(2000)设计了通过转动产生摩擦耗能的新型夹板摩擦耗能器。赵川等(2003)基于低造价、高性能、易生产、便施工的指导思想,从我国国情出发,设计了用于钢筋混凝土支撑的钢板-橡胶摩擦耗能器。张艳霞等(2015)对长孔螺栓摩擦阻尼器进行了往复荷载作用下的试验研究,研究摩擦接触面材料、长孔宽度及有无碟形垫片对阻尼器工作性能的影响。韩建强等(2018)设计了滑动长孔高强螺栓摩擦阻尼器,研究了不同摩擦板材料在不同螺栓扭力条件下阻尼器摩擦性能、耗能性能和滞回性能等,研究结果表明不同摩擦板材料阻尼器的滞回稳定性和耗能能力均较好。师骁等(2016)提出用于高层建筑的含摩擦阻尼器钢连梁,采用高强螺栓与碟形弹簧串联为摩擦片-钢板摩擦界面施加正压力,研究结果表明,选用的摩擦材料具有可靠的摩擦性能,摩擦阻尼器表现出稳定且优异的耗能能力;摩擦力相对较稳定,未出现过大的超强,有利于简化周边连接构件的设计。隋伟宁等(2018)设计了转动型摩擦阻尼器,该阻尼器通过钢板间转动摩擦消耗地震能量,并采用ABAQUS有限元软件对阻尼器力学性能进行数值模拟,通过改变阻尼器螺栓预紧力、转动摩擦系数和摩擦面个数等设计参数,确定影响阻尼器力学性能的主要参数。

    由于竖缝多位于纵横墙交接处,个别在较长剪力墙中间,因此本文提出将摩擦阻尼器设置于装配式剪力墙结构竖缝中(见图 1),在2片剪力墙中分别埋设锚栓和连接钢板,组装时通过螺栓将阻尼器与对应位置的预埋件连接,形成既能有效传递剪力,又能在竖向实现耗能的新型工艺。该工艺未采用传统预制混凝土装配式结构依靠钢筋混凝土暗柱进行连接的方法,而是以耗能连接装置取而代之,在保证水平缝可靠连接的基础上,利用竖缝阻尼器实现耗能,即满足“强水平缝、弱竖向缝”。本文借鉴已有理论和研究成果,对摩擦阻尼器进行设计,摩擦材料分别采用橡胶板、铝板、黄铜板、碳纤维板,进行拟静力循环加载试验,对阻尼器力学性能指标等进行研究。

    图 1  竖缝连接节点示意图
    Figure 1.  Schematic diagram of vertical joints

    在水平荷载作用下,竖缝两侧墙肢产生相对位移,因此考虑将摩擦片设置于阻尼器两侧预埋板和连接板之间,并用高强螺栓连接,同时高强螺栓为摩擦片提供预紧力,使用过程中通过改变螺栓预紧力及摩擦片材料调整阻尼器起滑力等参数,阻尼器连接板可根据竖缝尺寸进行调整,便于安装。结合图 1所示竖缝连接方式,根据试验条件,完成试验所用摩擦阻尼器的设计(见图 2)。试验所用阻尼器由开设长孔的滑板和夹板及高强螺栓组成,夹板的作用与图 1中预埋板的作用相同;在夹板与长孔滑板之间可嵌入不同材料的摩擦片,以获得可靠、稳定的耗能能力;中间滑板的连接位置开有水平方向的长螺孔,便于根据竖缝尺寸进行调整,以满足需要;在长孔滑板的滑动方向开设长孔,以容许可能出现的位移。摩擦片选用橡胶板、铝板、黄铜板、碳纤维板4种材料,所有材料厚度均为5mm,采用12.9级M16高强螺栓施加预紧力,并在试验前对其进行标定。为保证预紧效果,获得较稳定的预紧力,在螺栓与连接板接触位置设置碟形垫片。

    图 2  试验装置
    Figure 2.  Test equipments

    本试验根据《建筑抗震试验规程》(JGJ/T 101—2015)(中华人民共和国住房和城乡建设部,2015),采用位移控制的方式加载(见图 3),采用500kN作动器施加荷载,并以6mm为加载步幅,最大加载位移为30mm,加载速率为0.5mm/s,每级加载循环3周,若阻尼器起滑力超出作动器加载范围,则终止加载。加载方向为先向下后向上,因此定义向下为正,向上为负。沿作动器加载方向在连接板对角位置设置2个位移计(LVDT),用于量测位移,并取其平均值作为阻尼器滑动位移。将作动器荷载作为阻尼器实际承受的荷载。

    图 3  加载方案
    Figure 3.  Loading scheme

    螺栓预紧力对阻尼器工作性能参数具有重要作用,试验前首先对扭力扳手进行标定,确定扳手扭矩与螺栓预紧力之间的数量关系,得到以下拟合方程:

    $$ P = 0.7378M + 17.549 $$ (1)

    式中,P为2个12.9级M16高强螺栓预紧力;M为扭力扳手扭矩。

    开始加载后,试验装置发出“嘣嘣”声,这是由于作动器与阻尼器之间的连接螺栓在加载过程中逐渐被挤压引起的。当阻尼器开始出现滑移现象后,可观察到提供预紧力的高强螺栓发生了倾斜,铝板和黄铜板出现掉落碎屑的现象,而碳纤维板则未发现碎屑。碳纤维板进行试验时,当预紧力为70N·m,加载位移达18mm时,作动器荷载已达600kN,超出作动器加载极限,因此提前终止加载。试验后取出上述3种材料观察其表面,材料表面均有明显的划痕,说明在试验过程中钢板表面有小凸起,且由于预紧力的作用嵌入了摩擦材料,在加载过程中不仅存在摩擦力,咬合作用也起到了增大摩擦力的效果。橡胶板由于较软,试验结束后观察其表面,可知较试验前粗糙,螺栓孔部位变形明显,与连接板长孔重合部分被挤压变形,说明橡胶板强度较低,后期由于其破坏,对阻尼器的工作性能产生一定影响。加载结束后发现高强螺栓螺纹有挤压变形痕迹,螺栓与连接钢板直接接触的部位也有摩擦痕迹,也说明了高强螺栓在加载过程中发生了一定移动,摩擦片试验后的状态如图 4所示。

    图 4  试验后的摩擦片
    Figure 4.  Friction plates after test

    本次试验主要分析摩擦材料、预紧力等因素对摩擦阻尼器工作性能的影响,起滑力、起滑位移及等效初始刚度等是反映阻尼器性能的重要参数,其不仅能保证阻尼器在地震作用下为装配式剪力墙结构提供附加阻尼和足够的耗能能力,也是结构设计的重要依据。地震时结构受往复荷载作用,为考虑推拉作用下阻尼器起滑力是否稳定,引入摩擦力不均匀系数α,其值为最大正向起滑力与最大负向起滑力之比的绝对值,试验结果如表 1所示。

    表 1  试验结果
    Table 1.  Test results
    摩擦材料 预紧力
    /N·m
    起滑力
    F+/kN
    起滑力
    F-/kN
    摩擦力不均匀系数
    α
    起滑位移
    /mm
    等效初始刚度
    /N·mm-1
    黄铜板 30 245.880 -239.824 1.025 8.368 29383.365
    50 335.691 -296.601 1.132 11.113 30207.055
    70 355.744 -346.125 1.028 11.841 30043.409
    铝板 30 161.454 -159.961 1.009 6.086 26528.755
    50 215.533 -201.473 1.070 10.010 21531.768
    70 289.810 -281.060 1.031 10.630 27263.405
    碳纤维板 30 264.255 -323.789 0.816 12.009 22004.746
    50 346.364 -372.126 0.931 13.230 26180.197
    橡胶板 30 107.320 -71.060 1.510 4.002 26816.592
    50 127.390 -107.863 1.181 4.442 28678.523
    70 158.803 -124.248 1.278 5.179 30662.869
    90 172.096 -140.694 1.223 6.274 27428.280
    110 231.150 -175.435 1.318 6.678 34613.657
    下载: 导出CSV 
    | 显示表格

    摩擦阻尼器起滑力是体现其工作性能的重要指标,图 5所示为不同摩擦材料在2个加载方向起滑力变化趋势,由图 5可知起滑力近似符合线性变化趋势,且起滑力均随着预紧力的提高呈现增加的趋势,说明随着预紧力的提高,滑动摩擦力随之提高,在相同的加载位移下阻尼器可实现更大的耗能。在地震过程中,结构受往复荷载作用,若阻尼器在2个受力方向的起滑力接近,说明其工作性能良好。为此,通过摩擦力不均匀系数α反映阻尼器2个受力方向起滑力的差异。由表 1可知,黄铜板在预紧力为50N·m时α最大,为1.132;铝板在预紧力为50N·m时α最大,为1.070;橡胶板在预紧力为30N·m时α最大,为1.510,说明黄铜板、铝板在2个受力方向的起滑力较接近,而橡胶板工作性能较差,这是由于橡胶板较软、强度低,在正向加载过程中已造成材料的局部损坏,因此影响了负向起滑力。碳纤维板在预紧力为30N·m时α为0.816,正向加载时的起滑力小于负向加载,但2个受力方向起滑力之差小于20%,这是由于在前期加载过程中,造成了碳纤维板局部破坏,且高强螺栓未能完全拧紧。对4种摩擦材料进行对比分析,可知不同材料起滑力均随着预紧力的提高而增加,表现出一定规律,摩擦材料应具有一定强度,以保证其具有良好的工作性能。

    图 5  起滑力变化趋势对比
    Figure 5.  Trend contrast of sliding forces

    本文设计的摩擦阻尼器除提供刚度外,还应在较大荷载作用下尽早耗能,即阻尼器起滑位移不应太大。图 6所示为阻尼器起滑位移,由图 6表 1可知,随着预紧力的增加,不同摩擦材料起滑位移随之增加,但增幅逐渐减小。出现这种情况的主要原因是当预紧力较小时,起滑位移主要受连接板与摩擦材料产生的摩擦力影响,当预紧力增大后,与摩擦材料接触的钢板表面小凸起嵌入摩擦材料中;影响起滑位移的因素除摩擦外,还包括钢板表面凸起与摩擦材料的咬合作用,当预紧力达到一定程度后,咬合作用提供的摩擦力基本不再增加,因此,预紧力提高后起滑位移的变化逐渐减小。试验结束后观察摩擦片表面可知存在明显划痕,也说明了咬合作用的存在,对于增大摩擦力、增加阻尼器耗能而言,咬合作用有利。

    图 6  阻尼器起滑位移
    Figure 6.  Sliding displacements of the damper

    摩擦阻尼器初始刚度对于装配式剪力墙结构在小震作用下保持结构整体工作性能具有非常重要的作用,图 7所示为4种材料在不同预紧力下的刚度变化,由图 7可知,黄铜板、铝板等效初始刚度随着预紧力的增加变化不大,铝板在预紧力为50N·m时等效初始刚度有所降低,这是由于试验中螺栓未拧紧。碳纤维板和橡胶板等效初始刚度随着预紧力的增加而增加,橡胶板在预紧力为90N·m时等效初始刚度有所降低,这是由于螺栓未拧紧,使得加载过程中螺栓产生倾斜变形。受作动器加载量程的限制,碳纤维板在预紧力为70N·m已超出其最大输出荷载,因此实际获得的数据较少,需开展补充试验进行深入研究,但已有数据仍表现出一定规律性。不同材料的等效初始刚度呈现一定变化规律,表现出良好的工作性能。

    图 7  阻尼器等效初始刚度
    Figure 7.  Equivalent initial stiffness of the damper

    图 8所示为4种摩擦材料滞回曲线,由图 8可知,各材料滞回曲线均较饱满,当阻尼器达到起滑力后,滑动摩擦力先下降然后逐渐回升;铝板、黄铜板和碳纤维板在加载末期出现了滑动摩擦力上升的现象。结合试验中观察到的现象,当阻尼器开始滑移后,高强螺栓随之产生倾斜,由此造成预紧力减小;随着加载位移的增加,高强螺栓逐渐紧固,使预紧力逐渐增大,且随着预紧力的增加,摩擦片两侧钢板表面的小凸起逐渐嵌入其中,此时阻尼器除受摩擦作用影响,还受咬合作用影响,因此滑动摩擦力降低后又出现逐渐升高的趋势,且后期滑动摩擦力明显提高,有利于摩擦阻尼器耗能。而对于橡胶板,由于其较软,且弹性较大,因此在加载过程中受到的影响较小,滑动摩擦力较稳定,但随着小凸起的嵌入作用,后期滑动摩擦力略有提升。摩擦阻尼器在推拉转换过程中,在4组滞回曲线右下角及左上角存在明显阶梯状拐角,说明在加载过程中螺栓产生了偏移变形。当加载方向改变时,由于螺栓倾斜方向随之变化,因此在滞回曲线上可观察到位移变化较明显,而荷载却几乎未变化,这种情况对于阻尼器耗能不利。因此,应采取有效措施避免高强螺栓发生偏移变形,以减少预紧力的损失,保证摩擦阻尼器提供稳定的工作性能。

    图 8  阻尼器滞回曲线
    Figure 8.  Hysteretic curves of the damper

    本文进行不同摩擦材料、不同工况下滑动长孔螺栓摩擦阻尼器低周往复荷载试验研究,分析阻尼器摩擦力的产生机理、滞回性能、等效初始刚度、耗能能力及不同参数对其工作性能的影响,主要得出以下结论:

    (1)通过改变摩擦阻尼器中的摩擦材料,并对其起滑力、起滑位移及等效初始刚度等参数进行分析可知,该阻尼器工作性能良好,耗能稳定;不同摩擦材料之间存在差异,使得阻尼器表现出不同的工作性能。

    (2)黄铜板、铝板及碳纤维板摩擦力不均匀系数接近1,说明这3种摩擦材料在2个受力方向的摩擦力较接近,能承受地震作用下的往复荷载,并实现稳定耗能;而对于橡胶板,由于其自身强度较低,在正向加载时已造成了材料破坏,其摩擦力不均匀系数最大为1.510,说明其在2个受力方向上的耗能略差,因此对于摩擦材料而言,还需具备一定强度,避免受力过程中过早破坏。

    (3)高强螺栓为摩擦阻尼器提供预紧力,通过试验可知,高强螺栓偏移变形对于滑动摩擦力有很大影响,因此应尽量避免螺栓发生偏移,同时尽量采用强度较高的螺栓,避免受力过程中螺栓因破坏影响阻尼器耗能。

    (4)在工程应用中应在高强螺栓中加设碟形垫片,以保证高强螺栓始终处于紧固状态,受力过程中能提供稳定的预紧力,以实现摩擦阻尼器良好的工作性能。

  • 图  1  方案1三维模型

    Figure  1.  A three-dimensional model of design#1

    图  2  方案1侧视图

    Figure  2.  The side-map of design#1

    图  3  方案2三维模型

    Figure  3.  A three-dimensional model of design# 2

    图  4  方案2俯视图

    Figure  4.  The overlook-map of design#2

    图  5  方案3三维模型

    Figure  5.  A three-dimensional model of design#3

    图  6  方案3正视图

    Figure  6.  The positive-map of design#3

    图  7  加载制度

    Figure  7.  Rules for loading

    图  8  方案1施加荷载位置

    Figure  8.  Load location of design#1

    图  9  方案2施加荷载位置

    Figure  9.  Load location of design#2

    图  10  方案3施加荷载位置

    Figure  10.  Load location of design#3

    图  11  简化杆系受力状态

    Figure  11.  The stress state of simplified rod system

    图  12  3种设计方案的应力云图

    Figure  12.  Stress cloud diagram of three design schemes

    图  13  3种设计方案的变形云图

    Figure  13.  The deformation cloud diagram of three design schemes

    图  14  方案3的前3阶振型图

    Figure  14.  The first three-order vibration mode of design#3

    图  15  3种分析部下长肢墙的变形曲线

    Figure  15.  The deformation curve of the long limb wall under the three types of analysis

    表  1  ABAQUS统一单位

    Table  1.   Unified unit of ABAQUS

    长度 时间 荷载 质量 密度 应力
    m s N kg kg/m3 Pa
    mm s N 103kg 103kg/mm3 MPa
    下载: 导出CSV

    表  2  ABAQUS输入的钢材参数

    Table  2.   Input steel parameters in ABAQUS

    材料 密度/kN·m-3 杨氏模量/N·mm-2 泊松比
    Q235B钢 78.5 2.06×105 0.3
    下载: 导出CSV

    表  3  实际荷载值、最大挠度限值及实际最大变形值

    Table  3.   Actual load value, maximum deflection limit and actual maximum deformation value

    方案1 方案2 方案3
    根部最大剪力/kN 5 5 5
    根部最大弯矩/kN·m 7.25 6.25 7.25
    最大挠度限值/mm 1 0.83 0.94
    短肢墙最大变形/mm 0.525 0.499 0.187
    长肢墙最大变形/mm 0.734 0.540 0.441
    下载: 导出CSV

    表  4  根部最大允许剪力和最大允许弯矩

    Table  4.   Maximum allowable shear and maximum allowable bending moment of root

    方案1 方案2 方案3
    根部允许最大剪力/kN 6.8 7.7 10.7
    根部允许最大弯矩/kN·m 12.24 11.55 18.19
    下载: 导出CSV

    表  5  3种方案的自重及市场估价

    Table  5.   The weights and market valuation of the three designs

    设计方案 类型 材料 重量/kg 市场估价/万元
    方案1 板材满布形式 Q235B钢 5470 20
    方案3 型钢钢架 Q235B钢 4125 14
    方案3 型钢钢架 Q235B钢 3348 10
    下载: 导出CSV

    表  6  方案3的前3阶振型参数

    Table  6.   The first three-order vibration mode parameters of design#3

    振型 方向 特征值 固有频率/Hz 自振周期/s
    1阶振型 X轴正向平动+扭转 28072 26.666 0.0375
    2阶振型 X轴负向平动+扭转 46798 34.430 0.0290
    3阶振型 扭转 68393 41.622 0.0240
    下载: 导出CSV

    表  7  参数值

    Table  7.   Related parameters

    参数 构件和连接类别
    1 2 3 4 5 6 7 8
    C 1940×1012 861×1012 3.26×1012 2.18×1012 1.47×1012 0.96×1012 0.65×1012 0.41×1012
    β 4 4 3 3 3 3 3 3
    下载: 导出CSV

    表  8  3种分析部的最大变形值

    Table  8.   The maximum deformation value of the three types of analysis

    分析部时长 1s 10s 100s
    长肢墙最大变形值/mm 1.286 1.286 1.286
    下载: 导出CSV
  • 马永欣, 郑山锁, 2001.结构试验.北京:科学出版社.
    万宇, 2008.新疆大学结构实验室X方向反力墙有限元分析.乌鲁木齐: 新疆大学.
    王墩, 赵海琼, 吕西林, 2014.建筑结构拟静力试验方法的加载制度.四川建筑科学研究, 40(3):54-59. doi: 10.3969/j.issn.1008-1933.2014.03.014
    张敏政, 2015.地震工程的概念和应用.北京:地震出版社.
    张耀春, 2011.钢结构设计原理.北京:高等教育出版社.
    中华人民共和国建设部, 2003.钢结构设计规范GB 50017-2003.北京: 中国计划出版社.
  • 期刊类型引用(5)

    1. 付重阳,熊二刚,胡勤斌,张海荣,徐涵,樊禹江. 装配式剪力墙新型水平连接装置的力学性能数值模拟. 西安建筑科技大学学报(自然科学版). 2024(01): 82-92 . 百度学术
    2. 祁鹏飞,裴强,谯裕豪,鲁蕴华,夏超南. 基于摩擦摆支座的核安全壳结构隔震效果研究. 工程抗震与加固改造. 2024(04): 120-127 . 百度学术
    3. 隋伟宁,马勇,杨海涛,吴金国. 预制混凝土外挂墙板与钢框架摩擦耗能新型连接节点力学性能试验研究. 工业建筑. 2023(05): 109-117 . 百度学术
    4. 李成玉,马延港,陈焰周,吴东平. 柱端设置L形连接件滑移摩擦节点单层钢框架抗震性能研究. 世界地震工程. 2023(03): 56-67 . 百度学术
    5. 董标,杨进京,殷洪波,徐国林. 隧道钢拱节点让压性能分析. 隧道建设(中英文). 2021(09): 1524-1529 . 百度学术

    其他类型引用(6)

  • 加载中
图(15) / 表(8)
计量
  • 文章访问数:  96
  • HTML全文浏览量:  39
  • PDF下载量:  12
  • 被引次数: 11
出版历程
  • 收稿日期:  2018-03-15
  • 刊出日期:  2018-12-01

目录

/

返回文章
返回