多阶式层级梯度蜂窝结构的共面冲击响应

李成兵 李锐 张吉涛 叶强 李仁富

阚明先, 刘利新, 南小龙, 计策, 何勇, 段书超. 磁驱动样品实验数值模拟研究[J]. 高压物理学报, 2023, 37(6): 062301. doi: 10.11858/gywlxb.20230683
引用本文: 李成兵, 李锐, 张吉涛, 叶强, 李仁富. 多阶式层级梯度蜂窝结构的共面冲击响应[J]. 高压物理学报, 2023, 37(3): 034203. doi: 10.11858/gywlxb.20230604
KAN Mingxian, LIU Lixin, NAN Xiaolong, JI Ce, HE Yong, DUAN Shuchao. Numerical Simulation of Magnetically Driven Sample Experiment[J]. Chinese Journal of High Pressure Physics, 2023, 37(6): 062301. doi: 10.11858/gywlxb.20230683
Citation: LI Chengbing, LI Rui, ZHANG Jitao, YE Qiang, LI Renfu. In-Plane Impact Response of Multi-Order Hierarchical Gradient Honeycomb Structure[J]. Chinese Journal of High Pressure Physics, 2023, 37(3): 034203. doi: 10.11858/gywlxb.20230604

多阶式层级梯度蜂窝结构的共面冲击响应

doi: 10.11858/gywlxb.20230604
基金项目: 国家自然科学基金(552174209)
详细信息
    作者简介:

    李成兵(1977-),男,博士,教授,主要从事爆炸与冲击动力学研究. E-mail:lichbing@mail.ustc.edu.cn

  • 中图分类号: O347.1

In-Plane Impact Response of Multi-Order Hierarchical Gradient Honeycomb Structure

  • 摘要: 为改善蜂窝结构共面的力学性能,基于传统六边形蜂窝结构,建立了六边形层级蜂窝结构,并利用层级蜂窝代替传统六边形蜂窝部分胞元层,复合成一种新型多阶式层级梯度蜂窝结构。利用显式动力学有限元方法研究了层级梯度蜂窝的共面在不同冲击速度作用下的冲击响应特性和能量吸收能力。研究结果表明:层级梯度蜂窝的变形模式与塑性坍塌强度和冲击速度有关;层级梯度蜂窝冲击端和固定端在不同冲击速度作用下的名义应力-应变曲线均与其变形模式有关;不同的复合方式会导致层级梯度蜂窝具有不同的平台应力和比吸能,且在高速冲击时其平台应力比传统六边形蜂窝提高45.4%~63.8%,能量吸收提升10.8%~34.1%。相对密度会影响层级梯度蜂窝的能量吸收能力。

     

  • 自1999年美国圣地亚国家实验室开展磁驱动平面实验以来,磁驱动飞片发射、磁驱动冲击压缩、磁驱动准等熵压缩、磁驱动样品研究等平面实验取得了飞速发展,促进了高压状态方程、高能量密度物理和高功率脉冲技术等研究的发展[16]。随着磁驱动平面实验的开展,磁驱动平面实验理论取得了巨大进展。2001年,Reisman等[7]首次采用磁流体力学程序对磁驱动飞片实验进行了数值模拟。2003年,Lemke等[8]提出了磁驱动平面实验中边界磁场的初始公式;2011年,他们进一步发展了该公式,认为边界磁场强度与电极宽度和阴阳电极之间的间隙两者之和成反比[9]。基于上述边界磁场的磁流体力学程序并不能完全正确模拟磁驱动双侧飞片发射实验[10],在磁驱动实验后期,模拟的飞片自由面速度与实验测量结果的偏差较大。2018年,Kan等[11]指出,电流不是从阴、阳电极电流加载端的端面通过,而是从阴、阳电极电流加载端的液/气界面通过,从而提出具有烧蚀影响的边界磁场公式,解决了磁驱动双侧飞片实验后期模拟与测量不一致的问题。

    磁驱动样品实验主要用于样品材料的状态方程参量、强度等研究。本研究将采用考虑烧蚀影响的边界磁场公式和磁流体力学程序,对大电流脉冲功率装置上的磁驱动样品实验进行模拟研究,分析磁驱动样品实验结构系数的影响因素,确定结构系数与样品材料、样品厚度、电极宽度的依赖关系,以期为磁驱动样品实验的数值设计和预测提供理论支撑。

    磁驱动样品实验的负载结构如图1所示,其中:阴极和阳极采用铝金属长方体,厚度和宽度相同,设δ为初始厚度,W为初始宽度,g0为阴、阳极板之间的初始间隙。磁驱动样品实验采用上下结构,在阴阳极上下部嵌入两对半径为R的圆柱形飞片,飞片后界面紧贴半径为rr<R)的圆柱体样品,样品后界面紧贴半径为r的圆柱体LiF窗口,设δfcδfa分别为阴、阳极上飞片的厚度,δscδsa分别为阴、阳极上样品的厚度,δwcδwa分别为阴、阳极上窗口的厚度。采用激光速度干涉仪测量样品/窗口界面速度。若样品厚度为零,即飞片与窗口之间无样品,则实验测量的速度为飞片/窗口界面速度。磁驱动样品实验的负载参数列于表1

    图  1  磁驱动样品实验负载结构示意图
    Figure  1.  Load configuration of magnetically driven sample experiments
    表  1  磁驱动样品实验负载参数
    Table  1.  Load parameters for magnetically driven sample experiments
    Exp. No. Sample material δfa/mm δsa/mm δwa/mm δfc/mm δsc/mm δwc/mm g0/mm W/mm
    Exp1-top MgAl 0.970 0 8 0.983 0.565 8 1.200 13
    Exp1-bottom MgAl 0.984 0.968 8 0.988 0.788 8 1.200 13
    Exp2-top Sn 1.000 0.386 8 1.000 0 8 1.200 11
    Exp2-bottom Sn 1.000 0.386 8 1.000 0 8 1.200 11
    Exp3-top Sn 0.993 0.417 8 0.990 0.451 8 1.200 13
    Exp3-bottom Sn 0.993 0 8 0.990 0.643 8 1.200 13
    Exp4-top Sn 0.995 0.433 8 0.998 0.650 8 1.175 15
    Exp4-bottom Sn 0.990 0.817 8 0.991 0 8 1.175 15
    下载: 导出CSV 
    | 显示表格

    二维磁驱动数值模拟程序(two-dimensional magnetically driven simulation code,MDSC2)是基于二维四边形网格的拉格朗日磁流体力学程序,包含rz平面[12]xy平面[13]两种程序模式。rz平面模式用于求解磁驱动柱对称问题,xy平面模式用于求解磁驱动飞片等平面实验问题。MDSC2程序已成功用于Z箍缩[14]、磁驱动飞片发射[11]、带窗口磁驱动准等熵压缩[15]等磁驱动实验的模拟、设计和研究。磁驱动柱面实验和磁驱动平面实验的边界磁感应强度公式是不同的。对于磁驱动平面实验,MDSC2中的电极电流加载端的边界磁感应强度B

    B=fμ0I(t)/[g(t)+gm(t)+W]
    (1)

    式中:f为结构系数,μ0为真空磁导率,I(t)为实验电流,g(t)为阴、阳电极之间的间隙,gm(t)为阴、阳电极电流加载端厚度方向上气体宽度之和[11]

    采用MDSC2程序,对大电流脉冲装置上的磁驱动样品实验进行模拟。Exp1-top和Exp1-bottom为同一电流驱动上下两组MgAl样品的磁驱动实验(Exp1),阴阳极的极板宽度为13 mm。Exp1-top为磁驱动单侧MgAl样品实验:阳极上飞片与窗口直接相连,实验测量的是飞片/窗口界面速度;阴极上飞片与窗口之间嵌入0.565 mm厚的MgAl样品,实验测量的是MgAl/窗口界面速度。Exp1-bottom为磁驱动双侧MgAl样品实验,阴、阳极上飞片与窗口之间嵌入不同厚度的MgAl样品,厚度分别为0.788和0.968 mm,实验测量的速度为样品/窗口界面速度。磁驱动样品实验Exp1的电流变化曲线见图2图3为磁驱动单侧MgAl样品实验Exp1-top和磁驱动双侧MgAl样品实验Exp1-bottom的界面速度。可见,MDSC2程序模拟的界面速度与实验测量的界面速度基本一致,说明MDSC2程序能正确模拟磁驱动MgAl样品实验Exp1。由图2图3可知,磁驱动MgAl样品实验Exp1的结构系数为0.78。

    图  2  Exp1~Exp4的实测电流
    Figure  2.  Measured current for Exp1–Exp4
    图  3  Exp1的界面速度(f=0.78)
    Figure  3.  Interface’s velocity for Exp1 (f=0.78)

    Exp2-top和Exp2-bottom为同一电流驱动上下两组Sn样品的磁驱动实验(Exp2),初始结构完全相同,阴阳极的极板宽度均为11 mm。阳极上飞片与窗口之间都嵌入0.386 mm厚的Sn样品,阴极上飞片与窗口直接相连。Exp2实验中,仅对Exp2-bottom的界面进行了速度测量。磁驱动Sn样品实验Exp2的实测电流见图2图4为Exp2-bottom磁驱动单侧Sn样品实验的界面速度。由图4可知,MDSC2程序能正确模拟磁驱动Sn样品实验Exp2。结合图2可得,Exp2的结构系数为0.84。

    图  4  Exp2-bottom的界面速度(f=0.84)
    Figure  4.  Interface’s velocity for Exp2-bottom (f=0.84)

    Exp3-top和Exp3-bottom为同一电流驱动上下两组Sn样品的磁驱动样品实验(Exp3),阴阳极的极板宽度为13 mm。其中:Exp3-top为磁驱动双侧Sn样品实验,阴、阳电极上飞片与窗口之间嵌入不同厚度的Sn样品,阴、阳极上Sn样品的厚度分别为0.451和0.417 mm,实验测量的是0.451 mm Sn/LiF界面速度;Exp3-bottom为磁驱动单侧Sn样品实验,阳极上飞片与窗口直接相连,阴极上飞片与窗口之间嵌入0.643 mm厚的Sn样品,实验测量的是0.643 mm Sn/LiF界面速度。Exp3的实测电流见图2图5给出了Exp3-top和Exp3-bottom的界面速度,可见,MDSC2程序能正确模拟Exp3。结合图2可得,Exp3的结构系数为0.85。

    图  5  Exp3的界面速度(f=0.85)
    Figure  5.  Interface’s velocity for Exp3 (f=0.85)

    Exp4-top和Exp4-bottom为同一电流驱动上下两组Sn样品的磁驱动样品实验(Exp4),阴阳电极的极板宽度为15 mm。其中:Exp4-top为磁驱动双侧Sn样品实验,阴、阳电极上飞片与窗口之间嵌入不同厚度的Sn样品(厚度分别为0.433和0.650 mm),实测速度为Sn/LiF界面速度;Exp4-bottom为磁驱动单侧Sn样品实验,阳极上飞片与窗口之间嵌入0.817 mm厚的Sn样品,实验测量的是0.817 mm Sn/LiF界面速度,阴极上飞片与窗口直接相连,实验测量的是Al/LiF界面速度。Exp4的实测电流见图2。Exp4-top和Exp4-bottom的界面速度如图6所示,可见,MDSC2程序能正确模拟Exp4。结合图2可知,Exp4的结构系数为0.88。

    图  6  Exp4的界面速度(f=0.88)
    Figure  6.  Interface’s velocity for Exp4 (f=0.88)

    综合以上结果,MDSC2程序能正确模拟大电流脉冲功率装置磁驱动样品实验,所模拟的磁驱动样品/窗口界面速度与实验测量的界面速度一致。

    尽管利用式(1)的磁流体力学程序MDSC2能正确模拟磁驱动样品实验,但是,对于不同的磁驱动样品实验,其结构系数不同。由于结构系数不同时,飞片的边界磁场强度也会不同,致使飞片所受的洛伦兹力不同,因此,同一模型的磁流体力学模拟结果显然不同。磁驱动样品实验的结构系数直接影响磁驱动样品实验模拟的正确性,也影响磁驱动样品实验设计、预测以及实验数据解读的准确性,因而确定结构系数的影响因素和变化规律非常重要。

    由第2节可知,磁驱动样品实验的结构系数是一个常数,仅由磁驱动样品实验的初始条件决定。磁驱动样品实验的初始条件包括:电极板的宽度、阴阳电极之间的间隙、飞片材料类型、飞片厚度、样品材料类型、样品厚度、窗口材料类型、窗口厚度等。对于表1中的磁驱动样品实验:飞片均为金属铝,飞片厚度约1 mm,相对偏差不大于3%,可以看成同一厚度;窗口材料均为LiF,窗口厚度均为8 mm。因此,初始条件中变化的因素包括电极板宽度、阴阳极间隙、样品材料类型和样品厚度。表2列出了磁驱动样品实验的结构系数。

    表  2  磁驱动样品实验的结构系数
    Table  2.  Structure coefficients of magnetically driven sample experiments
    Exp. No. Sample material δsa/mm δsc/mm g0/mm W/mm f
    Exp1-top MgAl 0 0.565 1.200 13 0.78
    Exp1-bottom MgAl 0.968 0.788 1.200 13 0.78
    Exp2-top Sn 0.386 0 1.200 11 0.84
    Exp2-bottom Sn 0.386 0 1.200 11 0.84
    Exp3-top Sn 0.417 0.451 1.200 13 0.85
    Exp3-bottom Sn 0 0.643 1.200 13 0.85
    Exp4-top Sn 0.433 0.650 1.175 15 0.88
    Exp4-bottom Sn 0.817 0 1.175 15 0.88
    下载: 导出CSV 
    | 显示表格

    Exp1-top和Exp1-bottom的样品材料、阴阳极间隙、电极板宽度都相同,尽管样品厚度明显不同,但结构系数相同,Exp3、Exp4也有相似的结果,说明样品厚度对结构系数的影响很小。

    在Exp1和Exp3中,极板宽度均为13 mm,阴阳极间隙均为1.200 mm,样品分别为MgAl合金和Sn,结构系数分别为0.78和0.84,考虑到样品厚度对结构系数的影响很小,因此结构系数与样品材料有关。

    在Exp2、Exp3和Exp4中,样品材料均为Sn,阴阳极间隙分别为 1.200、1.200和1.175 mm(相对偏差很小,可看作同一间隙),电极宽度分别为11、13、15 mm(相对偏差很大),结构系数分别为0.84、0.85和0.88,说明结构系数的影响因素主要是阴阳电极板宽度,极板宽度越宽,结构系数越大。

    采用考虑烧蚀影响的边界磁场公式和MDSC2程序,对磁驱动样品实验进行了模拟,并分析了磁驱动样品实验结构系数的影响因素和规律。结果表明:考虑烧蚀影响的磁流体力学程序能正确模拟磁驱动样品实验,结构系数与实验过程无关,仅由实验的负载结构确定。结构系数与样品材料相关,样品材料不同,结构系数不同;样品厚度对结构系数的影响很小;电极宽度对结构系数的影响较大,电极宽度越宽,结构系数越大。研究结果为磁驱动样品实验的设计、预测等提供了理论依据。

  • 图  蜂窝单胞结构

    Figure  1.  Unit cell of honeycomb structure

    图  中阶层级蜂窝在不同l2下的应力-应变曲线

    Figure  2.  Stress-strain curves of middle-order hierarchical honeycomb under different l2

    图  层级梯度蜂窝有限元计算模型

    Figure  3.  Schematic diagram of finite element model of hierarchical gradient honeycomb

    图  层级梯度蜂窝结构设计模型

    Figure  4.  Structure design model of hierarchical gradient honeycomb

    图  力-位移曲线及变形模式的有限元模型验证

    Figure  5.  Finite element model verification of force-displacement curves and deformation patterns

    图  层级梯度蜂窝的变形模式

    Figure  6.  Deformation modes of hierarchical gradient honeycomb

    图  层级梯度蜂窝冲击端的名义应力-应变曲线

    Figure  7.  Nominal stress-strain curves at the impact end of hierarchical gradient honeycomb

    图  层级梯度蜂窝固定端的名义应力-应变曲线

    Figure  8.  Nominal stress-strain curves at the supporting end of hierarchical gradient honeycombs

    图  层级梯度蜂窝的名义应力-应变曲线和能量吸收效率曲线

    Figure  9.  Nominal stress-strain curve and energy absorption efficiency curve of hierarchical gradient honeycomb

    图  10  不同冲击速度下层级梯度蜂窝的平台应力

    Figure  10.  Plateau stress of hierarchical gradient honeycomb at different impact velocities

    图  11  不同相对密度下层级梯度蜂窝的平台应力

    Figure  11.  Plateau stress of hierarchical gradient honeycomb at different relative densities

    图  12  不同冲击速度下层级梯度蜂窝的比吸能特性

    Figure  12.  Specific energy absorption characteristics of hierarchical gradient honeycomb at different impact velocities

    图  13  不同相对密度下层级梯度蜂窝的比吸能特性

    Figure  13.  Specific energy absorption characteristics of hierarchical gradient honeycomb at different relative densities

    表  1  单胞结构的几何参数

    Table  1.   Geometric parameters of a unit cell

    Structurel/mmt/mmh/mmα/(°)
    HL15.000.30120
    HL23.500.301.15120
    HL32.000.301.15120
    下载: 导出CSV

    表  2  材料参数[7]

    Table  2.   Material parameters[7]

    Materialρ/(kg·m–3)E/GPaμσy/MPa
    Aluminum2 700690.376
    Rigid plate7 800210
    下载: 导出CSV

    表  3  不同速度下层级梯度蜂窝的密实应变

    Table  3.   Densification strains of hierarchical gradient honeycomb at different impact velocities

    v/(m·s−1)εd
    L123L321L213L312
    100.622 20.609 40.615 80.609 4
    300.647 70.621 50.634 60.621 6
    500.673 50.641 40.673 80.654 3
    1000.755 00.776 10.747 90.762 0
    下载: 导出CSV

    表  4  相对密度参数

    Table  4.   Relative density parameters

    IdentifierStructuret/mml/mmρcom
    RD1HL30.282.000.228
    HL20.283.50
    HL10.285.00
    RD2HL30.302.000.244
    HL20.303.50
    HL10.305.00
    RD3HL30.322.000.259
    HL20.323.50
    HL10.325.00
    下载: 导出CSV
  • [1] WANG Z G, LIU J F. Numerical and theoretical analysis of honeycomb structure filled with circular aluminum tubes subjected to axial compression [J]. Composites Part B: Engineering, 2019, 165: 626–635. doi: 10.1016/j.compositesb.2019.01.070
    [2] 陈鹏, 侯秀慧, 张凯. 面内冲击荷载下半凹角蜂窝的抗冲击特性 [J]. 高压物理学报, 2019, 33(6): 064104.

    CHEN P, HOU X H, ZHANG K. Impact resistance of semi re-entrant honeycombs under in-plane dynamic crushing [J]. Chinese Journal of High Pressure Physics, 2019, 33(6): 064104.
    [3] KHAN M K, BAIG T, MIRZA S. Experimental investigation of in-plane and out-of-plane crushing of aluminum honeycomb [J]. Materials Science and Engineering: A, 2012, 539: 135–142. doi: 10.1016/j.msea.2012.01.070
    [4] 于鹏山, 刘志芳, 李世强. 新型仿生蜂窝结构的设计与耐撞性能分析 [J]. 高压物理学报, 2022, 36(1): 014204.

    YU P S, LIU Z F, LI S Q. Design and crashworthiness analysis of new bionic honeycomb structure [J]. Chinese Journal of High Pressure Physics, 2022, 36(1): 014204.
    [5] 张新春, 沈振峰, 吴鹤翔, 等. 多段填充复合蜂窝结构的动态响应特性研究 [J]. 湖南大学学报: 自然科学版, 2020, 47(4): 67–75.

    ZHANG X C, SHEN Z F, WU H X, et al. Study on dynamic response characteristics of multi-segment filled composite honeycomb structure [J]. Journal of Hunan University: Natural Science, 2020, 47(4): 67–75.
    [6] 刘颖, 何章权, 吴鹤翔, 等. 分层递变梯度蜂窝材料的面内冲击性能 [J]. 爆炸与冲击, 2011, 31(3): 225–231.

    LIU Y, HE Z Q, WU H X, et al. In-plane impact properties of hierarchical gradient gradient honeycomb materials [J]. Explosion and Shock Waves, 2011, 31(3): 225–231.
    [7] 李坚, 孟卫华, 张大海, 等. 分层密度梯度蜂窝材料面内动态压缩及吸能特性 [J]. 应用力学学报, 2021, 38(6): 2369–2375.

    LI J, MENG W H, ZHANG D H, et al. In-plane dynamic compression and energy absorption characteristics of hierarchical density gradient honeycomb materials [J]. Chinese Journal of Applied Mechanics, 2021, 38(6): 2369–2375.
    [8] NIU X Q, XU F X, ZOU Z, et al. In-plane dynamic crashing behavior and energy absorption of novel bionic honeycomb structures [J]. Composite Structures, 2022, 299: 116064. doi: 10.1016/j.compstruct.2022.116064
    [9] QIAO J X, CHEN C Q. In-plane crushing of a hierarchical honeycomb [J]. International Journal of Solids and Structures, 2016, 85/86: 57–66.
    [10] WANG Z G, DENG J J, LIU K, et al. Hybrid hierarchical square honeycomb with widely tailorable effective in-plane elastic modulus [J]. Thin-Walled Structures, 2022, 171: 108816. doi: 10.1016/j.tws.2021.108816
    [11] LIU H, ZHANG E T, NG B F. In-plane dynamic crushing of a novel honeycomb with functionally graded fractal self-similarity [J]. Composite Structures, 2021, 270: 114106. doi: 10.1016/j.compstruct.2021.114106
    [12] LI Z, JIANG Y, WANG T, et al. In-plane crushing behaviors of piecewise linear graded honeycombs [J]. Composite Structures, 2018, 207: 425437.
    [13] WU H X, LIU Y, ZHANG X C. In-plane crushing behavior and energy absorption design of composite honeycombs [J]. Acta Mechanica Sinica, 2018, 34(6): 1108–1123. doi: 10.1007/s10409-018-0798-4
    [14] GIBSON L J, ASHBY M F. Cellular solids [M]. Cambridge University Press, 1997.
    [15] GIOVANNI T, MAICON C R. On the convergence of the primal hybrid finite element method on quadrilateral meshes [J]. Applied Numerical Mathematics, 2022, 181: 552560.
    [16] TAO Y, CHEN M J, CHEN H S, et al. Strain rate effect on the out-of-plane dynamic compressive behavior of metallic honeycombs: experiment and theory [J]. Composite Structures, 2015, 132: 644651.
    [17] 张新春, 刘颖. 密度梯度蜂窝材料动力学性能研究 [J]. 工程力学, 2012, 29(8): 372–377.

    ZHANG X C, LIU Y. Research on dynamic properties of density gradient honeycomb materials [J]. Engineering Mechanics, 2012, 29(8): 372–377.
    [18] 何强, 马大为, 张震东. 分层屈服强度梯度蜂窝材料的动力学性能研究 [J]. 工程力学, 2015, 32(4): 191–196.

    HE Q, MA D W, ZHANG Z D. Study on dynamic properties of layered yield strength gradient honeycomb materials [J]. Engineering Mechanics, 2015, 32(4): 191–196.
    [19] 沈振峰, 张新春, 白江畔, 等. 负泊松比内凹环形蜂窝结构的冲击响应特性研究 [J]. 振动与冲击, 2020, 39(18): 89–95, 117.

    SHEN Z F, ZHANG X C, BAI J P, et al. Study on impact response characteristics of negative Poisson’s ratio inner concave ring honeycomb structure [J]. Journal of Vibration and Shock, 2020, 39(18): 89–95, 117.
    [20] 白江畔, 张新春, 沈振峰, 等. 冲击载荷下多胞元薄壁结构的动态压溃行为研究 [J]. 振动与冲击, 2020, 39(18): 145–152.

    BAI J P, ZHANG X C, SHEN Z F, et al. Study on dynamic crushing behavior of multi-cell thin-walled structures under shock loading [J]. Vibration and Shock, 2020, 39(18): 145–152.
    [21] 乔及森, 李明, 苗红丽. 串联梯度蜂窝结构的面内力学性能 [J]. 塑性工程学报, 2021, 28(11): 115–123.

    QIAO J S, LI M, MIAO H L. In-plane mechanical properties of series gradient honeycomb structures [J]. Chinese Journal of Plastic Engineering, 2021, 28(11): 115–123.
  • 加载中
图(13) / 表(4)
计量
  • 文章访问数:  250
  • HTML全文浏览量:  115
  • PDF下载量:  36
出版历程
  • 收稿日期:  2023-01-09
  • 修回日期:  2023-02-27
  • 录用日期:  2023-04-03
  • 网络出版日期:  2023-05-19
  • 刊出日期:  2023-06-05

目录

/

返回文章
返回