Processing math: 100%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

基于子模型方法的地铁车体疲劳寿命评估

李新康 王苏秦 刘潮涛 张嘉 陈逸斐

李新康, 王苏秦, 刘潮涛, 张嘉, 陈逸斐. 基于子模型方法的地铁车体疲劳寿命评估[J]. 西南交通大学学报, 2022, 57(2): 295-300, 330. doi: 10.3969/j.issn.0258-2724.20200301
引用本文: 李新康, 王苏秦, 刘潮涛, 张嘉, 陈逸斐. 基于子模型方法的地铁车体疲劳寿命评估[J]. 西南交通大学学报, 2022, 57(2): 295-300, 330. doi: 10.3969/j.issn.0258-2724.20200301
LI Xinkang, WANG Suqin, LIU Chaotao, ZHANG Jia, CHEN Yifei. Fatigue Life Assessment of Metro Carbody Based on Submodel Method[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 295-300, 330. doi: 10.3969/j.issn.0258-2724.20200301
Citation: LI Xinkang, WANG Suqin, LIU Chaotao, ZHANG Jia, CHEN Yifei. Fatigue Life Assessment of Metro Carbody Based on Submodel Method[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 295-300, 330. doi: 10.3969/j.issn.0258-2724.20200301

基于子模型方法的地铁车体疲劳寿命评估

doi: 10.3969/j.issn.0258-2724.20200301
详细信息
    作者简介:

    李新康(1987—),男,高级工程师,研究方向为轨道交通车辆车体结构的仿真分析与优化,E-mail:lixinkang@cqsf.com

  • 中图分类号: U270.14

Fatigue Life Assessment of Metro Carbody Based on Submodel Method

  • 摘要:

    车体是地铁车辆的重要承载部件,针对全尺寸车体疲劳试验综合难度高的问题,基于端部底架子模型包含了牵枕缓等在车体中受力最严苛结构的特点,提出采用端部底架子模型代替全尺寸车体进行疲劳试验的方法;建立了端部底架子模型和全尺寸车体的有限元模型,并参照EN 12663标准所确定的疲劳载荷,通过设置合理的边界条件使得端部底架子模型与全尺寸车体关键位置应力分布一致;将试验测得端部底架关键位置的应力与仿真结果进行对比,验证了有限元模型的准确性,进而采用名义应力法和Eurocode 9标准规定的疲劳寿命-应力(S-N)曲线对车体和底架焊缝部位进行了疲劳损伤计算. 结果表明:端部底架3个最大损伤位置与全尺寸车体一致,并且同一位置处端部底架的损伤值均大于车体损伤值,因此采用子模型法评估全尺寸车体的疲劳寿命易于获得相对保守的结果,针对地铁车体采用子模型法进行疲劳寿命评估是可行的.

     

  • 车体作为地铁车辆关键承载零部件,在运行过程中承受来自转向架、车钩、车下设备以及乘客质量等一系列复杂载荷. 新造地铁车体由于轻量化需求,多采用铝合金型材整体焊接而成,一节车厢焊缝多达数千条. 以传统A型地铁为例,车体侧墙两边各开有5门4窗,且门窗尺寸远大于干线车体,大大增加了车身强度薄弱位置的数量.基于以上特点,地铁车体在承受复杂载荷时的疲劳问题越来越突出.

    传统A型地铁铝合金车体长约为22 m、宽度约为3 m、高度约为3 m,整备质量约为25 t,超员载荷约为50 t. 整车的疲劳试验试验周期长、经费高、风险大. 在ANSYS软件中,子模型方法通过对局部关注区域的网格进行细化,并在区域边界处施加等效的边界条件,使得计算结果更精确,计算速度更快. 故可考虑截取车体受力最恶劣的部分,对该部分的边界施加等效边界条件后进行疲劳试验. 根据车体的受力状态以及运营维修经验,车体疲劳裂纹往往产生于牵枕缓部位[1],因此可以截取车体底架端部包含牵枕缓部位的部分进行子模型疲劳试验.

    胡杰鑫等[2]利用Monte Carlo方法,对铝合金地铁车体枕梁进行了虚拟疲劳试验,并与实际疲劳试验结果进行对比,验证了该方法的有效性;王玉伟[3]基于实测载荷编制CRH2型车枕梁关键部位的应力谱,研究了运营里程与损伤之间的关系;李珊珊[4]利用子模型方法,对CRH2型车枕梁结构网格进行了细分,提出了枕梁的补强方案;臧伟锋等[5]针对飞机机身壁板设计了一种新型工装,该工装使得机身壁板的应力分布与实际全尺寸机身的应力分布相一致,并进行了内压载荷疲劳试验;四方股份联合航空工业飞机强度研究所对CRH380中间车车体全长的1/3结构进行了气密强度疲劳试验[6];Gutiérrez-Carvajal等[7]设置了多种载荷工况对铝合金车体枕梁进行了疲劳试验,利用累积损伤法评估了枕梁的疲劳寿命.

    本文将设定边界条件下的静强度试验结果与有限元计算结果进行对比分析,验证了端部底架子模型(以下简称端部底架)有限元模型的正确性. 研究了端部底架边界条件模拟方法和试验加载方法,根据标准EN12663-1:2010 + A1:2014[8](以下简称EN 12663标准)规定的载荷进行了疲劳损伤计算,得到了端部底架与全尺寸车体关键位置损伤分布的规律,验证了以端部底架子模型疲劳试验代替全尺寸车体疲劳试验的可行性. 通过端部底架疲劳试验结果对全尺寸车体的疲劳寿命进行了评估.

    在全尺寸车体模型中,车体的底架边梁与车体的侧墙连接,底架的端梁与车体的端墙连接,底架的尾部与车体地板及底架边梁连接. 截取下来的端部底架长5 m,与端墙、侧墙、地板以及边梁的连接均被截断. 所设置的边界条件需保证端部底架关键部位的应力分布与应力大小与全尺寸车体相一致.

    在有限元计算时,容易获取端部底架边界节点的力和弯矩,该力和弯矩与施加于端部底架的外部载荷使得底架处于力平衡的状态. 在疲劳试验过程中,底架车钩、空簧和中心销处的集中力可以通过作动器施加,由于底架边界在全尺寸车体中变形较小,边界处的分布力则以约束产生的反力施加.

    根据端部底架与全尺寸车体的连接方式,在端部底架的头部两侧与车体侧墙连接处、枕梁两侧与车体侧墙连接处约束端部底架的垂向自由度. 在尾部地板处约束底架的纵向和垂向自由度. 横向约束在运营过程中由车体的惯性产生,在台架试验中无法复现,因此选取与牵枕缓焊缝应力无关的边梁侧面进行约束,如图1所示,图中:X轴所示方向为纵向;Y轴所示方向为横向;Z轴所示方向为垂向.

    图  1  端部底架边界条件
    Figure  1.  End underframe boundary conditions

    根据EN 12663标准中关于车体疲劳载荷的规定:地铁车体应根据定员载荷(m1 + m3)和超员载荷(m1 + m4)分别设计疲劳载荷工况,其中:m1为车体整备质量;m3为定员乘客质量;m4为超员乘客质量,并按其在服役寿命期间所占的比例进行疲劳试验. 其他类型轨道交通车辆车体根据其载荷特点,选取恒幅定员或恒幅超员疲劳载荷进行疲劳试验.

    由于地铁车体定员载荷和超员载荷在服役期间所占比重难以准确统计,为得到较为保守的结果,认为车体在服役过程中一直承受最大载荷,车体的疲劳寿命最短,即以恒幅超员载荷(m1 + m4)进行数值计算和疲劳试验,m1 = 28.279 t,m4 = 22.440 t.

    EN 12663标准中规定整车纵向和横向疲劳载荷为 ± 0.15g,垂向载荷为g ± 0.15g,但均为单个载荷分别作用1000万次的疲劳载荷. 对于组合载荷,EN 12663标准第6.8节推荐引用VDV 152-2016[9]标准(以下简称VDV 152标准).

    VDV 152标准第5.3.2节中将X、Y、Z方向的载荷进行组合,建立了3个载荷工况,每个工况循环200万次,见表1.

    表  1  标准之间疲劳载荷工况对比
    Table  1.  Comparison of fatigue load cases between standards
    标准工况XYZ
    VDV 152 1 ± 1.20 m/s2 g ± 1.90 m/s2
    2 ± 1.60 m/s2 g ± 1.00 m/s2
    3 ± 1.50 m/s2 ± 0.60 m/s2 g ± 1.00 m/s2
    EN 12663 1 ± 0.15g
    2 ± 0.15g
    3 ± 0.15g
    下载: 导出CSV 
    | 显示表格

    VDV 152标准中载荷工况的组合与EN 12663标准中规定的3个方向载荷每个方向单独作用1000万次是等效的,即在疲劳试验中,可用VDV 152标准规定的疲劳载荷代替EN 12663标准中规定疲劳载荷[9].

    对于X方向的载荷,VDV 152标准只在一个工况下循环了200万次,因此需根据等损伤原则,将X方向200万次载荷折算成EN 12663标准中与垂向、横向共同作用的1000万次载荷.

    1.2.1   车钩载荷

    车钩载荷主要源于列车牵引或制动时的不同步,对于A型地铁车辆,一般采用4动2拖编组,拖车为首尾车,最大车钩载荷发生在拖车和动车连接的车钩位置. 按启动加速度a = 1 m/s2计算,最大车钩载荷为

    FXmax=4a(m1+m4)2=101.44kN. (1)

    根据VDV 152标准第5.3.2.1节,最大车钩力仅在牵引(制动)不同步的瞬间产生,对于动力分散列车,车钩力的疲劳载荷为最大车钩力的25%. 因此,作用200万次车钩载荷的动态幅值为${F_{X {\max} }}/{4}=25\; {\rm{KN}}$.

    根据疲劳寿命-应力曲线幂函数表达式:

    σmN=C. (2)

    式中:σ为应力;m、C为材料常数,VDV 152标准第5.3.3.3节推荐m取4;N为循环次数.

    按VDV 152标准,将0.15g作用200万次,根据式(2)折算成1000万次,则车钩载荷动态幅值$F_{\text {coupler }}= 16.96\; \mathrm{kN}.$

    1.2.2   中心销牵引(制动)载荷

    依据EN 12663标准按 ± 0.15g施加,中心销牵引(制动)载荷如式(3).

    FpivotX=0.15g(m1+m4)2=37.3kN. (3)

    按式(3)将200万次载荷折算成1000万次,载荷$F_{{\rm{pivot}}X}=24.94\; \mathrm{kN}.$

    中心销牵引(制动)载荷的加载高度根据实车确定,本文中加载高度距枕梁下表面为586.5 mm.

    1.2.3   中心销横向载荷

    列车在曲线运行时会承受来自于转向架横向止挡的横向力FY,依据EN 12663标准按± 0.15g施加,如式(4).

    FY=0.15g(m1+m4)2=37.3kN. (4)

    该载荷由中心销和空气弹簧共同承担,在试验时仅施加中心销承受的横向载荷. 中心销横向载荷加载高度距枕梁下表面为297 mm. 根据空气弹簧参数,每一个转向架单个空气弹簧承担的横向载荷${F_{{\rm{airsring}}Y}} = 3.78\;{\rm{kN}}.$

    因此中心销承担的横向载荷为

    FpivotY=FY2FairsringY=29.74kN. (5)
    1.2.4   空簧垂向载荷

    根据表1超员载荷确定空簧承受的垂向载荷. 依据EN 12663标准,按g± 0.15g施加,由空气弹簧承担.

    空簧载荷均值为

    FZm=g(m1+m4)4=124.39kN. (6)

    空簧载荷幅值为

    FZa=0.15g(m1+m4)4=18.66kN. (7)

    采用Hypermesh对端部底架进行网格划分,网格类型为Shell181. 底架所用材料为铝合金,材料密度为2.7×10−9 t/mm3,弹性模量为69 GPa,泊松比为0.3. 中心销所用材料为钢,材料密度为7.85×10−9 t/mm3,弹性模量为210 GPa,泊松比为0.3. 中心销与枕梁下表面通过梁单元(Beam188)模拟的螺栓连接. 端部底架有限元模型共划分单元213586个,节点197352个. 底架的约束方式如图2所示,FZZ方向的载荷.通过刚性元在车钩座、中心销和空簧位置施加集中力载荷.

    图  2  端部底架有限元模型
    Figure  2.  FE model of the end underframe

    全尺寸车体建模方式与端部底架相同,共划分单元1956597个,节点1779596个,如图3所示.

    车体在4个空簧位置约束垂向,在中心销横向止挡位置约束横向. 通过刚性元在车钩座、中心销空簧位置施加纵向集中力载荷,垂向载荷及横向载荷以重力场形式施加. 乘客质量通过柔性元均布在地板上,空调及车下设备等考虑其对整车刚度的加强,通过刚性元连接到设备安装座,车门及车窗通过刚性元连接到安装位置.

    图  3  全尺寸车体有限元模型
    Figure  3.  FE model of the full-size carbody

    试验过程中,利用应变片对端部底架关键部位进行应力测量,测点布置见图4S1~S9为应变片测点,R1~R4为应变花测点,应变花用相当应力进行评定.

    当载荷加到最大时,读取各个位置的应力值如图5所示.

    图  4  测点布置
    Figure  4.  Measuring points arrangement

    图5可知:测点S1、S2、S3、S4、S5、S9、R1、R2、R3、R4的实测应力与计算应力相对误差较小,测点S3误差最小为0.64%;测点S6、S7、S8误差较大,其中测点S8误差最大为122.2%.

    图  5  实测应力与计算应力对比
    Figure  5.  Comparison of measured stress and calculated stress

    测点S8位于枕内加强筋的焊趾处,该处应力梯度较大,且打磨后的几何轮廓不平整,难以准确确定焊趾位置,因此应变片位置的细微差别会导致测量结果较大的变化;测点S6与测点S3为对称点、测点S7与测点S1为对称点,其仿真应力值误差分别为8.2%和12.2%,且测点S1和测点S3的实测应力值和仿真应力值相差较小,因此测点S6和测点S7误差较大的原因可能是试验中的约束条件没有仿真中的约束条件理想,约束工装本身与车体之间存在间隙导致.

    本文中端部底架所需评估的部位主要位于枕梁内侧加强框与周边零部件的焊缝,无法直接粘贴应变片. 测点S6、S7、S8与所评估部位无直接联系,且测点部位在实际运营过程中状态良好. 除测点S6、S7、S8外,其余测点实测应力与计算应力误差较小,因此,认为端部底架有限元模型能够真实反映实际模型的应力状态.

    基于Miner线性累积损伤理论进行疲劳计算,如式(8).

    D=niNi, (8)

    式中:D为损伤值,当D<1时,认为材料没有疲劳失效的风险.;ni为第i级应力谱下循环的次数;Ni为第i级应力谱下的疲劳寿命.

    在ANSYS下计算单位静载荷下焊缝单元的名义应力,在nCode Designlife中施加恒幅疲劳载荷谱.

    端部底架同时承受来自车钩座、中心销横向、中心销纵向和空簧垂向载荷,在试验过程中其单元的主应力方向不断发生变化,属于典型的多轴疲劳. 因此采用带绝对值的最大主应力进行疲劳分析.

    在疲劳试验中,端部底架承受空簧载荷均值产生的应力,需进行平均应力修正,本文中选取FKM指导书[10]中规定的方法进行平均应力修正. 列车车体在实际服役过程中,结构承载的载荷循环次数远远超过材料疲劳极限处的循环次数,大量小幅值载荷仍能对车体造成损伤,因此将S-N曲线的第二段延伸至应力幅为0处,考虑小载荷循环的影响[11].

    本文根据Eurocode 9[12]标准建立端部底架铝合金材料不同焊缝类型的S-N曲线. $\Delta \sigma $代表循环应力范围. 根据端部底架焊缝类型选用Eurocode 9中7.2.3号焊缝、7.4.3号焊缝、9.4号焊缝和9.6号焊缝. 焊缝类型如图6所示.

    各焊缝对应的S-N曲线如图7所示. 由图7可知:各焊缝S-N曲线的拐点位于500万次循环处;S-N曲线的第一段斜率k1 = −1/3.4,第二段斜率k2 = −1/5.4.

    根据有限元模型与图6焊缝类型的对应关系,在nCode中将图7各焊缝的S-N曲线分别赋予给对应的焊缝单元,按1.2节计算的载荷进行疲劳计算. 得到端部底架和全尺寸车体的损伤云图如图8,对比如表2.

    图  6  焊缝类型示意
    Figure  6.  Weld type diagram
    图  7  焊缝S-N曲线
    Figure  7.  S-N curves of welding seams
    图  8  损伤云图
    Figure  8.  Damage nephogram
    表  2  端部底架与全尺寸车体损伤对比
    Table  2.  Damage comparison of the end underframe and the full-scale carbody
    端部底架损伤整车损伤损伤部位
    0.5028 0.1910  地板下表面与枕内加强筋底板焊缝直角拐角处
    0.1796 0.0962  地板下表面与枕内加强筋底板135°拐角附近区域
    0.0857 0.0591  外侧枕内加强筋与枕梁下盖板焊缝
    下载: 导出CSV 
    | 显示表格

    表2可以看出:端部底架损伤最大的3个部位与全尺寸车体损伤最大的3个部位位置一致,且同一位置端部底架的损伤均大于全尺寸车体损伤. 因此通过端部底架疲劳试验验证整车的疲劳强度可以获得较为保守的结果.

    端部底架包括了车体运营过程中受力最恶劣的牵枕缓结构,采用接口工装和MTS液压作动器进行疲劳加载,如图9.

    图  9  端部底架疲劳试验
    Figure  9.  End underframe fatigue testing

    端部底架的垂向和横向运动由头部及枕梁两侧边梁的“C”字形工装提供约束,纵向运动通过固定尾部地板面提供约束,如图10所示. 按1.2节计算的载荷进行1000万次循环的疲劳试验,对图8(a)所示的区域进行渗透探伤,未发现裂纹,因此由端部底架的疲劳试验结果认为整车的疲劳寿命满足EN 12663标准的要求.

    图  10  端部底架约束示意
    Figure  10.  Constraint diagram of the end underframe

    1) 全尺寸车体的疲劳试验对试验设备、试验成本均有较高要求,用小尺寸子模型代替整车进行疲劳试验能够有效地缩短研发周期,降低试验成本. 本文采用端部底架子模型代替全尺寸车体的疲劳损伤评估结果偏于保守,从数值仿真和台架试验的角度验证了方法对车体疲劳寿命进行评估是可行的.

    2) 使用子模型进行疲劳试验时,需保证子模型的边界远离可能的裂纹潜在部位,且保证子模型结构潜在危险点与全尺寸车体结构保持一致. 边界条件的设置可根据边界处刚度的方向特性,选取刚度最大的方向进行约束,释放他方向的约束. 通过有限元仿真设置的边界条件与试验进行对比,获得正确的有限元数值分析模型,为车体的寿命评估提供依据.

    3) 结合名义应力法S-N曲线和线性损伤理论,分别计算了端部底架与全尺寸车体的疲劳损伤. 发现该方法依赖于网格质量和S-N曲线的选取,针对车体和其焊缝特点,在计算中均进行了考虑,因此采用的方法适合车体的寿命评估.

    4) 采用端部底架子模型的疲劳试验与数值计算,并与全尺寸车体结构的对比分析,表明该车体能够满足EN 12663标准规定的疲劳寿命要求.

    致谢:铝合金车体疲劳寿命评估及试验方法研究(2682019CX45).

  • 图 1  端部底架边界条件

    Figure 1.  End underframe boundary conditions

    图 2  端部底架有限元模型

    Figure 2.  FE model of the end underframe

    图 3  全尺寸车体有限元模型

    Figure 3.  FE model of the full-size carbody

    图 4  测点布置

    Figure 4.  Measuring points arrangement

    图 5  实测应力与计算应力对比

    Figure 5.  Comparison of measured stress and calculated stress

    图 6  焊缝类型示意

    Figure 6.  Weld type diagram

    图 7  焊缝S-N曲线

    Figure 7.  S-N curves of welding seams

    图 8  损伤云图

    Figure 8.  Damage nephogram

    图 9  端部底架疲劳试验

    Figure 9.  End underframe fatigue testing

    图 10  端部底架约束示意

    Figure 10.  Constraint diagram of the end underframe

    表  1  标准之间疲劳载荷工况对比

    Table  1.   Comparison of fatigue load cases between standards

    标准工况XYZ
    VDV 152 1 ± 1.20 m/s2 g ± 1.90 m/s2
    2 ± 1.60 m/s2 g ± 1.00 m/s2
    3 ± 1.50 m/s2 ± 0.60 m/s2 g ± 1.00 m/s2
    EN 12663 1 ± 0.15g
    2 ± 0.15g
    3 ± 0.15g
    下载: 导出CSV

    表  2  端部底架与全尺寸车体损伤对比

    Table  2.   Damage comparison of the end underframe and the full-scale carbody

    端部底架损伤整车损伤损伤部位
    0.5028 0.1910  地板下表面与枕内加强筋底板焊缝直角拐角处
    0.1796 0.0962  地板下表面与枕内加强筋底板135°拐角附近区域
    0.0857 0.0591  外侧枕内加强筋与枕梁下盖板焊缝
    下载: 导出CSV
  • [1] 王生华. 上海轨道交通1号线列车车体裂纹原因分析及解决措施[J]. 铁道车辆,2018,56(3): 41-43. doi: 10.3969/j.issn.1002-7602.2018.03.016

    WANG Shenghua. Analysis of causes to cracking in carbodies of trains for Shanghai rail traffic No. 1 line and measures for solution[J]. Rolling Stock, 2018, 56(3): 41-43. doi: 10.3969/j.issn.1002-7602.2018.03.016
    [2] 胡杰鑫,谢里阳,喻海洋,等. 基于验证模型的枕梁疲劳寿命预测虚拟实验[J]. 西南交通大学学报,2019,54(1): 106-112.

    HU Jiexin, XIE Liyang, YU Haiyang, et al. Virtual experiments to predict bolster fatigue lifetime based on FEM model validated by static tests[J]. Journal of Southwest Jiaotong University, 2019, 54(1): 106-112.
    [3] 王玉伟. 基于实测载荷的CRH2型动车组车体枕梁疲劳强度研究[D]. 北京: 北京交通大学, 2015.
    [4] 李珊珊. CRH2046动车枕梁结构剩余寿命预测及补强方案分析[D]. 沈阳: 东北大学, 2015.
    [5] 臧伟锋,陈安,董登科. 机身壁板内压载荷试验研究[J]. 航空工程进展,2018,9(1): 69-76.

    ZANG Weifeng, CHEN An, DONG Dengke, et al. Test research on fuselage panel subjected to internal pressure load[J]. Advances in Aeronautical Science and Engineering, 2018, 9(1): 69-76.
    [6] 魏亚龙. 高速列车车体气密疲劳强度试验台气密载荷加载特性的仿真分析[D]. 北京: 北京交通大学, 2012.
    [7] GUTIÉRREZ-CARVAJAL R, BETANCUR GIRALDO G, BARBOSA J, et al. Full scale fatigue testing performed to the bolster beam of a railway vehicle[J]. International Journal on Interactive Design and Manufacturing (IJIDeM), 2016, 12(1): 253-261.
    [8] European Committee for Standardization. Railway applications—structural requirements of railway vehicle bodies, part 1: Locomotives and passenger rolling stock (and alternative method for freight wagons): EN 12663-1: 2010 + A1: 2014[S]. British: BSI Group, 2014.
    [9] Committee on Urban Rail Rolling Stock. Recommendation on the design for strength of urban rail rolling stock according to BOStrab: VDV 152-2016[S]. Germany: Verband Deutscher Verkehrsunternehmen, 2016.
    [10] Forschungskuratorium Maschinenbau. Analytical strength assessment of components : made of steel, cast iron and aluminum materials in mechanical engineering: FKM Guideline[S]. Frankfurt am Main: VDMA-Verl, 2013.
    [11] 李凡松. 服役环境下动车组车体振动与疲劳研究[D]. 成都: 西南交通大学, 2018.
    [12] European Committee for Standardization. Eurocode 9: Design of aluminium structures, part 1-3: structures susceptible to fatigue: EN 1999-1-3-2007+A1: 2011[S]. British: BSI Group, 2011.
  • 加载中
图(10) / 表(2)
计量
  • 文章访问数:  692
  • HTML全文浏览量:  339
  • PDF下载量:  43
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-05-26
  • 修回日期:  2020-09-09
  • 网络出版日期:  2020-09-30
  • 刊出日期:  2020-09-30

目录

/

返回文章
返回