• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

Ⅰ型裂纹中低速冲击荷载下起裂韧度测试新方法

施泽彬 朱哲明 汪小梦 王雄

刘文祥, 张庆明, 钟方平, 程帅, 张德志. 球壳塑性变形下的应变增长现象[J]. 爆炸与冲击, 2017, 37(5): 893-898. doi: 10.11883/1001-1455(2017)05-0893-06
引用本文: 施泽彬, 朱哲明, 汪小梦, 王雄. Ⅰ型裂纹中低速冲击荷载下起裂韧度测试新方法[J]. 爆炸与冲击, 2018, 38(6): 1247-1254. doi: 10.11883/bzycj-2017-0132
Liu Wenxiang, Zhang Qingming, Zhong Fangping, Cheng Shuai, Zhang Dezhi. Strain growth of spherical shell subjected to internal blast loading during plastic response[J]. Explosion And Shock Waves, 2017, 37(5): 893-898. doi: 10.11883/1001-1455(2017)05-0893-06
Citation: SHI Zebin, ZHU Zheming, WANG Xiaomeng, WANG Xiong. A new testing method for mode Ⅰ crack initiation fracture toughness under middle-low speed impacts[J]. Explosion And Shock Waves, 2018, 38(6): 1247-1254. doi: 10.11883/bzycj-2017-0132

Ⅰ型裂纹中低速冲击荷载下起裂韧度测试新方法

doi: 10.11883/bzycj-2017-0132
基金项目: 

国家自然科学基金项目 11672194

国家自然科学基金项目 11702181

四川省安全监管局安全生产科技项目 aj20170515161307

详细信息
    作者简介:

    施泽彬(1992-), 男, 硕士研究生

    通讯作者:

    朱哲明, zhemingzhu@hotmail.com

  • 中图分类号: O346;TU45

A new testing method for mode Ⅰ crack initiation fracture toughness under middle-low speed impacts

  • 摘要: 为了探寻更加合理的构型试件来研究纯Ⅰ型裂纹在冲击荷载下的起裂及扩展行为, 提出一种新构型试件, 即双倾斜底边中心裂纹试件(double inclined bottom central cracked, DIBCC)。借助于中低速落锤式冲击实验装置进行冲击实验, 通过应力波来使试件内预制裂纹起裂并扩展, 同时利用应变片测试系统监测裂纹起裂时刻, 并采用AUTODYN有限差分软件对实验过程进行数值模拟, 最后计算裂纹的动态应力强度因子, 利用实验测得的起裂时刻, 确定试件的起裂韧度。结果表明:(1)在反射拉伸波作用下, 预制裂纹两侧会产生垂直于裂纹面向外的位移, 使预制裂纹扩张, 从而使裂纹起裂。(2)数值模拟结果与实验结果在裂纹扩展路径上具有一致性, 说明本文中提出的DIBCC构型试件有效, 可以用来测试裂纹在冲击载荷下的断裂韧度。
  • 在爆炸载荷的作用下,容器壳体的最大应变没有出现在第一个应变峰上,而是出现在后期的应变峰上的现象被称为爆炸容器的应变增长现象[1]。应变增长现象中出现意料之外的应变量,会对容器安全形成威胁。自1976年被发现以来[2],应变增长现象引起了广泛关注,其形成原因可以归纳为两种:不同振动模态的叠加[1~6]以及爆炸载荷与壳体振动共振[7]。到目前为止,涉及应变增长现象的研究几乎都是在壳体弹性变形范围内。针对实验中发现的球壳塑性变形时的应变增长现象,本文中将分析该现象形成的原因以及相关规律。

    在球形容器中开展爆炸加载实验。容器内径为523 mm,壁厚为3 mm,壳体材料为Q345R,其静态屈服强度不低于345 MPa。球壳由两个半球焊接而成。容器顶部有内径125 mm的开孔,用于安装炸药。容器赤道面(即两个半球的焊缝面)上均布4个小开孔,用于安装压力传感器。实验采用球形装药,由细尼龙绳吊挂于容器中心。装药采用较严格的中心点起爆方式,由雷管(处于球壳上部的雷管腔内,与球壳内部空间隔开)引爆较细的柔爆索,柔爆索在中心起爆微型药球(直径10 mm,由微米级的PETN粉压制而成),微型药球再起爆球形主装药(由60%RDX和40%TNT浇铸而成),以形成较理想的球面冲击波。球形容器照片见图 1。容器外壁粘贴三向应变片,应变片的位置和方向见图 2

    图  1  球形容器照片
    Figure  1.  Photo of spherical vessel
    图  2  应变片位置和方向
    Figure  2.  Positions and directions of strain guages

    爆炸当量370 g TNT时球壳的应变数据见表 1。从表 1中可见,球壳应变的第一峰基本上都超过1%,远超过壳体的屈服应变。壳体材料动态屈服强度σy按620 MPa计算,参考下式:

    εy=σy(1ν)E
    (1)
    表  1  球壳应变数据
    Table  1.  Experimental strain of spherical shell
    应变片位置 方向 第一个应变峰/10-3 最大应变值/10-3 最大值应变值所处峰的位置 应变增长系数
    实测值 平均值
    S1 1 10.804 11.081±2.5% 10.804 1 1.00
    2 10.804 12.606 2 1.16
    3 11.635 13.063 2 1.12
    S2 1 10.802 10.262±6.6% 12.057 2 1.11
    2 10.397 12.057 2 1.15
    3 9.588 10.802 2 1.12
    S3 1* 11.945 11.225±6.4% 12.197 2 1.02
    (12.456) (3)
    2 10.808 11.610 2 1.07
    3 10.922 11.485 2 1.05
    S4 1 9.140 9.883±7.5% 9.140 1 1.00
    2 10.197 10.503 2 1.03
    3 10.312 10.503 2 1.01
    S5 1 11.821 13.159±14.3% 11.821 1 1.00
    2 15.036 - - -
    3 12.619 12.989 2 1.02
    S6 1 13.866 13.762±4.3% 15.211 2 1.09
    2 14.254 16.427 2 1.15
    3 13.167 14.571 2 1.10
    注:标注*的数据出现特殊情况,其第二个应变峰虽大于第一个应变峰,但应变最大值出现在第三个应变峰上。
    下载: 导出CSV 
    | 显示表格

    式中:泊松比ν=0.3,弹性模量E=200 GPa,可估算屈服应变εy=2.170×10-3图 3为球壳应变的典型曲线,壳体残余应变也超过了1%。

    图  3  典型的应变曲线
    Figure  3.  Strain-time curve of spherical shell

    图 3中还可以看出,应变最大值出现在应变曲线的第二个峰上,应变曲线出现了应变增长现象。由表 1数据可知,实验获得的有效应变数据为17组,其中14组出现了应变增长现象,13组出现了应变最大值出现在第二个应变峰上的情况,与图 3的应变曲线的特征一样。表 1中数据显示,应变增长系数最大达到1.16。

    对于振动叠加形成的应变增长现象,一般认为,弯曲波产生于扰动源(指容器上的开孔、法兰等)并传播至壳体某位置,与其他的壳体振动叠加形成应变增长现象[4-5]。由于弯曲波存在传播过程,在该应变增长现象中,不同位置球壳的应变增长现象的特征是不同的,比如应变最大值出现的时间存在差异。在本次实验出现的应变增长现象中,应变最大值均出现在第二个应变峰上,特征相同,因此其不是振动叠加形成的应变增长现象,而可能是由爆炸载荷与球壳振动共振引起的。

    图 4为实验测得的球壳内壁压力-时间曲线,曲线上周期性地出现了3个脉冲,脉冲之后容器内的压力基本趋于平静。Duffey数值模拟研究[7]表明,炸药产生的冲击波直接作用在球壳内壁形成第一个脉冲,之后冲击波返回球壳中心并反射回来再次作用在球壳内壁,形成第二个脉冲,如此反复不断产生周期性的脉冲。本次实验对应的比距离约为0.36 m/kg1/3,冲击波来回过程仅持续了3次,之后冲击波衰减至可以忽略。在当量更大、容器半径更小的情况下,冲击波更强,脉冲的数目可能大于3个。

    图  4  球壳内壁的周期性多脉冲载荷
    Figure  4.  Pressure-time curve of blast loading on inner wall of spherical shell

    借助一维球壳响应理论[8],可以由容器内壁压力-时间曲线分析出壳体响应情况。因为一维球壳除了呼吸振动外无其他振动,因此一维球壳出现的应变增长现象必定由爆炸载荷引起。

    一维球壳响应理论存在以下前提或假设:壳体材料密度变化忽略不计,即壳体材料不可压缩;球壳为薄壁,即截面上应力近似认为是均匀分布,意味着球壳半径r远大于壁厚h

    参考文献[8],假设球壳半径变化不大,壁厚不变,根据牛顿定律建立平衡方程:

    h(σθ+σφ)+ρ hrd2urdt2=rp(t)
    (2)

    式中:ρ为壳体材料密度,E为弹性模量,p(t)为容器内壁压力历程,ur为球壳径向位移,σθσφ为球壳的主应力。由于球壳为点对称模型,则

    σθ=σφ=σ
    (3)

    根据Hook定律,在壳体材料弹性变形范围内,存在关系式

    σ=E1ν ε=E1νurr
    (4)

    把式(3)、(4)代入式(2),得到

    d2urdt2+ω2ur=p(t)ρ h
    (5)

    其中

    ω2=2Eρr2(1ν)
    (6)

    设球壳材料的应力应变曲线关系为双线性弹塑性本构关系,壳体发生塑性变形时:

    σ=σy+Es(εεy)
    (7)

    式中:σyεy分别为容器材料屈服强度及屈服应变,ES为强化模量。把式(7)代入式(2),得到球壳塑性变形时的平衡方程:

    d2urdt2+ωp2ur=p(t)ρ h2(σyESεy)ρ r
    (8)

    其中

    ω2p=2ESρ r2
    (9)

    球壳弹性响应方程(5)和塑性响应方程(8)忽略了球壳半径和厚度的变化。考虑厚度的变化将更符合现实情况,变化的容器半径与厚度关系为:

    hr2=h0r20
    (10)

    借助数学软件Matlab求解以上式子。

    输入实测的球壳内壁压力-时间曲线(见图 4),计算出球壳响应曲线,见图 5。爆炸加载下壳体材料的动态屈服强度很难直接给出准确值,材料动态拉伸试验显示其屈服强度大部分在600~700 MPa之间,为了对比由压力曲线计算得到的应变曲线与真实应变曲线的特征,计算过程中壳体材料屈服强度取数值620 MPa(前文的620 MPa也来源于此),以使得计算的应变曲线的第一个应变峰与实验应变曲线的第一个应变峰相等。

    图  5  由爆炸载荷计算的应变曲线与实验应变曲线的比较
    Figure  5.  Comparison between strain-time curves obtained by experiment and calculated from pressure-time curve

    图 5比较了计算得到的球壳应变曲线和实验获得的应变曲线,计算的应变曲线与实验应变曲线的特征相似,出现了应变增长现象,且应变最大值出现在第二个应变峰上。由此可以判断,实验中应变增长现象由容器内壁的周期性多脉冲载荷引起。图 5还可以对比球壳响应曲线和球壳内壁压力-时间历程曲线的时间特征,可见仅前两个脉冲对应变增长现象有贡献。

    从响应方程(5)和(8)来看,第一个脉冲作用下壳体获得变形动力,且在内部应力和惯性力共同作用下进行振动。很显然,后续如出现第二个脉冲,且作用方向与壳体运动方向相同,将推动壳体变形在原变形基础上进一步增加,出现应变增长现象;相反,后续脉冲作用方向与壳体运动方向相反,会对壳体运动造成一定阻碍,导致壳体变形减少,不出现应变增长现象。从图 5可见,第二个脉冲持续的时间较长,要大于壳体振动的半个周期,脉冲的一部分作用方向与壳体运动方向相同,另一部分必定与壳体运动方向相反。必须注意,从压力衰减规律来看,脉冲的主要部分处于脉冲的前期。

    以球壳内壁压力-时间曲线(见图 4)为基础,人为设置第二个脉冲出现时刻不同的8种载荷,并计算得到相应的球壳响应情况,如图 6所示。以此分析第二个脉冲作用方向与壳体运动方向出现不同状况时的应变增长现象。

    图  6  后续脉冲出现时间对球壳响应的影响
    Figure  6.  Strain-time curves with second pulse appearing at different time points

    可见,第二个脉冲出现时刻在上一个应变峰的下降沿前段,脉冲的主要部分作用在该应变峰的下降沿,作用方向与壳体运动方向相反,导致下一个应变峰减少,无应变增长现象,如图 6(a)~(c);第二个脉冲出现时刻在上一个应变峰的下降沿后段,脉冲的小部分作用在该应变峰的下降沿,大部分作用在当前应变峰的上升沿,导致当前应变峰大于上一个应变峰,出现应变增长现象,如图 6(d);第二个脉冲的出现时刻在当前应变峰的上升沿,必然导致当前应变峰大于上一个应变峰,出现应变增长现象,如图 6(e)~(h),其中以第二个脉冲出现在当前应变峰的上升沿前段时的应变增长现象最严重,如图 6(e)~(f)。综上所述,后续脉冲的出现时刻决定应变增长现象的有无和强弱,本文实验中的第二个脉冲出现在应变峰的上升沿前段(见图 5),由此引起比较严重的应变增长现象。

    (1) 实验观察到球壳塑性变形量超过1%时的应变增长现象,应变增长系数最大为1.16。

    (2) 实验中应变增长现象是由球壳内壁周期性多脉冲载荷引起的,球壳内壁载荷存在3个较明显的脉冲,仅前两个脉冲对应变增长现象有贡献。

    (3) 第二个脉冲的出现时刻决定应变增长现象的有无和强弱,实验中的第二个脉冲出现在应变峰的上升沿前段,引起较严重的应变增长现象。

  • 图  1  DIBCC构型实验方法示意图

    Figure  1.  Test method of DIBCC specimen

    图  2  DIBCC构型模型示意图

    Figure  2.  Sketch of DIBCC specimen

    图  3  中低速冲击加载装置

    Figure  3.  Middle-low impact loading device

    图  4  数据采集系统

    Figure  4.  Data collection system

    图  5  入射端与透射端应变时程曲线

    Figure  5.  Curves of strain versus time recorded from incident and transmitted plates

    图  6  裂纹尖端应变片的电压信号时程曲线

    Figure  6.  Curves of voltage signal versus time recorded from strain gauges stuck at the crack tip

    图  7  3个时刻y方向应力云图

    Figure  7.  Nephograms of stress in y-direction at three moments

    图  8  裂纹尖端应力时程曲线及其右侧点的位移曲线

    Figure  8.  Curves of stress in x-direction versus time at gauge 19 and displacement in x-direction versus time at gauge 5

    图  9  数值模拟结果及实验得到的试件破坏结果

    Figure  9.  Test results of specimen fracture patterns and numerical simulation result

    图  10  裂纹张开位移示意图

    Figure  10.  Deformation of crack opening displacement

    图  11  裂纹面三个点的水平位移时程曲线

    Figure  11.  Displacement of three points on crack surface as a function of time

    图  12  动态应力强度因子时程曲线

    Figure  12.  Curve of dynamic stress intensity factor versus time

  • [1] LOEBER J F, SIH G C.Diffraction of antiplane shear waves by a finite crack[J].Acoustical Society of America, 1968, 44(1):90-98. doi: 10.1121/1.1911091
    [2] HARRIS J G.Diffraction by a crack of a cylindrical longitudinal pulse[J].Journal of Applied Mathematics and Physics, 1980, 31(3):367-383. doi: 10.1007%2FBF00236484
    [3] BERTRAM A, KALTHOFF J F.Crack propagation toughness of rock for the range of low to very high crack speeds[J].Key Engineering Materials, 2003, 251/252:423-430. doi: 10.4028/www.scientific.net/KEM.251-252
    [4] LEE D, TIPPUR H, BOGERT P.Dynamic fracture of graphite/epoxy composites stiffened by buffer strips:an experimental study[J].Composite Structures, 2012, 94(12):3538-3545. doi: 10.1016/j.compstruct.2012.05.032
    [5] 胡柳青, 李夕兵, 龚声武.冲击载荷作用下裂纹动态响应的数值模拟[J].爆炸与冲击, 2006, 26(3):214-221. doi: 10.3321/j.issn:1001-1455.2006.03.004

    HU Liuqing, LI Xibing, GONG Shengwu.Simulation on dynamic response of crack subjected to impact loading[J].Explosion and Shock Waves, 2006, 26(3):214-221. doi: 10.3321/j.issn:1001-1455.2006.03.004
    [6] 杨井瑞, 张财贵, 周妍, 等.用CSTBD试样确定砂岩的动态起裂和扩展韧度[J].爆炸与冲击, 2014, 34(3):264-271. http://www.bzycj.cn/CN/abstract/abstract8838.shtml

    YANG Jingrui, ZHANG Caigui, ZHOU Yan, et al.Determination of dynamic initiation toughness and propagation toughness of sandstone using CSTBD specimens[J].Explosion and Shock Waves, 2014, 34(3):264-271. http://www.bzycj.cn/CN/abstract/abstract8838.shtml
    [7] ZHANG Q B, ZHAO J.Effect of loading rate on fracture toughness and failure micromechanisms in marble[J].Engineering Fracture Mechanics, 2013, 102(6):288-309. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=JJ0229785505
    [8] ZHANG Z X, KOU S Q, YU J, et al.Effects of loading rate on rock fracture[J].International Journal of Rock Mechanics and Mining Sciences, 1999, 36(5):597-611. doi: 10.1016/S0148-9062(99)00031-5
    [9] 王蒙, 朱哲明, 胡荣.基于SCSCC试样的岩石复合型裂纹动态扩展特征研究[J].四川大学学报(工程科学版), 2016(2):57-65. http://d.old.wanfangdata.com.cn/Periodical/scdxxb-gckx201602009

    WANG Meng, ZHU Zheming, HU Rong.Rock experiments study of crack propagation under Ⅰmode and Ⅰ-Ⅱ mixed-mode dynamic loading using SCSCC specimens[J].Journal of Sichuan University (Engineering Science Edition), 2016(2):57-65. http://d.old.wanfangdata.com.cn/Periodical/scdxxb-gckx201602009
    [10] 王蒙, 朱哲明, 王雄.冲击荷载作用下的Ⅰ/Ⅱ复合型裂纹扩展规律研究[J].岩石力学与工程学报, 2016, 35(7):1323-1332. http://cdmd.cnki.com.cn/Article/CDMD-10674-1015635867.htm

    WANG Meng, ZHU Zheming, WANG Xiong.Study on the growth rule of mixed-modeⅠ/Ⅱ crack under impacting loads[J].Chinese Journal of Rock Mechanics and Engineering, 2016, 35(7):1323-1332. http://cdmd.cnki.com.cn/Article/CDMD-10674-1015635867.htm
    [11] 徐文涛, 朱哲明, 曾利刚.爆炸载荷下Ⅰ型裂纹动态断裂韧度测试方法初探[J].岩石力学与工程学报, 2015(增刊1):2767-2772.

    XU Wentao, ZHU Zheming, ZENG Ligang.Testing method study of mode-Ⅰ dynamic fracture toughness under blasting loads[J].Chinese Journal of Rock Mechanics And Engineering, 2015(suppl 1):2767-2772.
    [12] Zhu Z M, Xu W T, Feng R Q.A new method for measuring mode-Ⅰ dynamic fracture toughness of rock under blasting loads[J].Experimental Techniques, 2016, 40(3), 889-905. doi: 10.1007%2Fs40799-016-0093-x
    [13] Chong K P, Kuruppu M D.New specimen for fracture toughness determination of rock and other materials[J].International Journal of Fracture, 1984, 26:59-62. doi: 10.1007/BF01157555
    [14] 李清, 杨仁树.爆炸载荷裂纹扩展的应力强度因子及其断裂行为[J].煤炭学报, 2002, 27(3):290-293. doi: 10.3321/j.issn:0253-9993.2002.03.015

    LI Qing, YANG Renshu.Stress intensity factor and fracture behavior for crack propagation under blasting load[J].Journal Of China Coal Society, 2002, 27(3):290-293. doi: 10.3321/j.issn:0253-9993.2002.03.015
    [15] WANG Q Z, XING L.Determination of fracture toughness KIC by using the flattened Brazilian disc specimen for rocks[J].Engineering Fracture Mechanics, 1999, 64(2):193-201. doi: 10.1016/S0013-7944(99)00065-X
    [16] YANG S, TANG T X, ZOLLINGER D, et al.Splitting tension tests to determine concrete fracture parameters by peak-load method[J].Advanced Cement Based Materials, 1997, 5:18-28. doi: 10.1016/S1065-7355(97)90011-0
    [17] ZHU Z M, MOHANTY B, XIE H P.Numerical investigation of blasting-induced crack initiation and propagation in rocks[J].International Journal of Rock Mechanics and Mining Sciences, 2007, 44(3):412-424. doi: 10.1016/j.ijrmms.2006.09.002
    [18] ZHU Z M, XIE H P, MOHANTY B.Numerical investigation of blasting-induced damage in cylindrical rocks[J].International Journal of Rock Mechanics and Mining Sciences, 2008, 45(2):111-121. doi: 10.1016/j.ijrmms.2007.04.012
    [19] 朱哲明, 李元鑫, 周志荣, 等.爆炸荷载下缺陷岩体的动态响应[J].岩石力学与工程学报, 2011, 30(6):1157-1167. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201106009

    ZHU Zheming, LI Yuanxin, ZHOU Zhirong, et al.Dynamic response of defected rock under blasting load[J].Chinese Journal of Rock Mechanics and Engineering, 2011, 30(6):1157-1167. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201106009
    [20] ZHU Z M, WANG C, KANG J M.Study on the mechanism of zonal disintegration around an excavation[J].International Journal of Rock Mechanics and Mining Sciences, 2014, 67(4):88-95. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=JJ0232596842
    [21] 王蒙, 朱哲明, 谢军.岩石Ⅰ-Ⅱ复合型裂纹动态扩展SHPB实验及数值模拟研究[J].岩石力学与工程学报, 2015, 34(12):2474-2485.

    WANG Meng, ZHU Zheming, XIE Jun.Experimental and numerical studies of the mixed-mode Ⅰ andⅡcrack propagation under dynamic loading using SHPB[J].Chinese Journal of Rock Mechanics and Engineering, 2015, 34(12):2474-2485.
    [22] ROSSMANITH H P, DAEHNKE A, KNASSMILLNER R E K, et al.Fracture mechanics applications to drilling and blasting[J].Fatigueand Fracture Engineering Materials and Structures, 1997, 20(11):1617-1636. doi: 10.1111/ffe.1997.20.issue-11
    [23] KUTTER H K, FAIRHURST C.On the fracture process in blasting[J].International Journal of Rock Mechanics and Mining Sciences, 1971, 8:181-202. doi: 10.1016/0148-9062(71)90018-0
    [24] WANG M, ZHU Z M, DONG Y Q, et al.Study of mixed-mode Ⅰ/Ⅱ fractures using single cleavage semicircle specimens under impacting loads[J].Engineering Fracture Mechanics, 2017, 177:33-44. doi: 10.1016/j.engfracmech.2017.03.042
    [25] 宋力, 胡时胜.SHPB数据处理中的二波法与三波法[J].爆炸与冲击, 2005, 25(4):368-373. doi: 10.3321/j.issn:1001-1455.2005.04.014

    SONG Li, HU Shisheng.Two-wave and three-wave method in SHPB data processing[J].Explosion and Shock Waves, 2005, 25(4):368-373. doi: 10.3321/j.issn:1001-1455.2005.04.014
    [26] 岳中文.缺陷介质爆生裂纹扩展规律的动态焦散线试验研究[D].北京: 中国矿业大学(北京), 2009.
    [27] 郭伟国.应力波基础简明教程[M].西安:西北工业大学出版社, 2007:44-52.
  • 期刊类型引用(1)

    1. 周磊,朱哲明,董玉清,牛草原,邓帅,姜亚成. 动载荷作用下含倒U型孔洞裂纹试件的动力损伤行为研究. 岩石力学与工程学报. 2021(01): 113-125 . 百度学术

    其他类型引用(2)

  • 加载中
图(12)
计量
  • 文章访问数:  5191
  • HTML全文浏览量:  1406
  • PDF下载量:  41
  • 被引次数: 3
出版历程
  • 收稿日期:  2017-04-26
  • 修回日期:  2017-09-16
  • 刊出日期:  2018-11-25

目录

/

返回文章
返回