力学学报 -

姓名
邮箱
手机号码
标题
留言内容
验证码

2018, 50(6): 1277-1278.

HTML PDF

吴承康

郭永怀先生是我国著名的力学家和应用数学家. 在美国留学和工作期间, 他耕耘在空气动力学前沿, 为跨声速流与应用数学做出贡献. 在回国工作的12 年间, 他积极参加规划我国力学发展蓝图, 培育力学人才, 是我国近代力学的奠基人; 他为我国原子弹、氢弹的研制贡献了全部心血乃至生命, 是我国核武器事业的先驱. 郭永怀先生热爱祖国、艰苦奋斗、谦虚质朴、无私奉献, 他学识渊博, 治学严谨, 是我国科技工作者的楷模.
郭永怀先生是《力学学报》第二任主编. 在他的主持下,《力学学报》从创刊初期就具有很高的学术起点和公正、严谨、求实的学术声誉. 这里我们向广大读者介绍郭先生的光辉一生, 并刊出他的学生们的一些论文, 以表达我们对这位主编的怀念和崇敬.

2018, 50(6): 1279-1282.

HTML PDF

高超声速飞行复现风洞理论与方法

姜宗林, 李进平, 胡宗民, 刘云峰, 俞鸿儒

针对高超声速飞行伴随的热化学反应流动,本文回顾了郭永怀先生的科研理念和学科布局,综述了他亲手成立的高温气动团队在高超声速飞行风洞实验模拟理论与方法方面的研究进展.高温气体的迅速产生与迅速应用是一种理想的风洞运行方法,而激波管就是这样一种实验装备.论文首先介绍了激波管技术的基本理论与方程,指出将其用于高超声速流动实验模拟时所具有的独特优势.然后讨论了应用激波风洞复现需要的高超声速飞行状态的可行性、基本方程和需要解决的关键问题.针对这些关键问题,进一步介绍了如何应用爆轰现象研发激波风洞驱动技术的理论,并给出了基于爆轰驱动方法的技术发展和工程应用验证.最后,论文介绍了爆轰驱动激波风洞的界面匹配条件,该条件奠定了长实验时间激波风洞运行基础,是其他驱动方法尝试解决而没能完全解决的难题.高温气动团队关于高超声速飞行复现风洞的理论与技术研究,实现了郭永怀先生的战略规划,成就了国际领先的高超声速热化学反应流动研究平台.

2018, 50(6): 1283-1291. doi:10.6052/0459-1879-18-238

HTML PDF

近空间高超声速飞行器气动特性研究的若干关键问题

叶友达, 张涵信, 蒋勤学, 张现峰

在30$\sim$70km空域机动飞行的高超声速飞行器的优点是可以耦合利用所处空域的空气产生的升力和高速飞行的离心力进行远距离机动滑翔飞行,具有重要的实用价值.尽管过去数十年在高超声速流动研究方面取得显著进展,但在设计研究近空间远程滑翔的高超声速飞行器方面仍然存在许多挑战,特别是对特定飞行条件下的流动机理了解不清楚.本文介绍了作者研究团队在开展近空间高超声速飞行器有关的关键气动问题方面的研究进展,主要包括：建立了近空间高超声速飞行的流动模型,发展了系统的相关计算空气动力学方法,针对高空高速飞行条件下稀薄气体效应和真实气体效应的耦合作用影响研究了合适的滑移边界条件,考虑了不同组分存在条件下的温度、速度和压力的滑移效应影响;提出了飞行器气动外形的动态优化方法,获得了可工程实用化的高升阻比飞行器气动外形;建立了高速飞行器动稳定性理论,在实现高超声速飞行器动态稳定飞行方面取得重大进展;最后讨论了高超声速飞行器设计中进一步需要关注的若干关键技术和科学问题、可能解决的途径及其所涉及的学科发展方向.

2018, 50(6): 1292-1310. doi:10.6052/0459-1879-18-247

HTML PDF

小行星撞击地球的超高速问题

柳森, 党雷宁, 赵君尧, 白智勇, 黄洁, 李毅, 石义雷

小行星撞击地球是人类生存面临的潜在威胁之一.在小行星进入地球大气与撞击地球表面过程中,存在烧蚀、解体、空中爆炸、火球、撞击成坑、反溅碎片云、地震以及海啸等一系列复杂的物理化学和力学现象.本文梳理和归纳了与这些现象相关的超高速空气动力学问题和超高速碰撞动力学问题.小行星进入地球大气的超高速空气动力学问题有：极高速($V = 12 ~ 20$km/s)进入条件下的气动力与轨迹,极高速进入条件下的小行星气动加热与烧蚀机理,极高速气动加热条件下的小行星结构传热与热响应,极高速进入条件下的高温气体效应,小行星进入过程的物理特征.小行星撞击地球的超高速碰撞动力学问题有：陆地撞击成坑与反溅碎片云,海洋撞击与海啸,撞击过程的地震效应.由于小行星撞击地球与超高速飞行器的再入过程在速度、材料和结构上存在较大差异,针对这些超高速问题,现有的研究手段在地面试验和数值计算两方面都存在不足.最后,从小行星进入地球大气的弹道方程、质量损失方程、解体判据和解体模型等出发,初步建立了小行星进入与撞击效应分析评估模型,并对Chelyabinsk和Tunguska两次流星事件进行了分析,重构了进入与爆炸解体过程,评估了空爆火球在地面所导致的超压和热辐射损伤.

2018, 50(6): 1311-1327. doi:10.6052/0459-1879-18-313

HTML PDF

一个改进的电磁波传输激波管实验

熊壮, 张勇勇, 王苏, 李进平, 陈宏, 范秉诚, 崔季平

在中国科学院力学研究所$\varPhi $ 800 mm高温低密度激波管上进行电磁波在等离子体中传输机理研究时,低密度和强激波条件下,由于气体解离和电离等非平衡过程,使得激波后2区宽度显著减小;同时由于边界层效应造成激波衰减和接触面加速,使得激波后2区长度进一步减小.这两个效应导致激波管2区实验观测时间减小,2区气体处于非平衡状态,增加了观察数据的不稳定性和数据分析的难度.本文提出在$\varPhi 800 $ mm高温低密度激波管中采用氩气(Ar)和空气(Air)混合气替代纯空气作为激波管实验介质气体.利用Ar不解离和难电离的特性,减小激波前后压缩比,从而增加激波后2区实验时间和气体长度. 采用Langmuir 静电探针和微波透射诊断技术测量激波后电子密度,同时利用探针测量激波后2区实验时间.结果显示,在Ar+Air混合气实验中,激波波后电子密度可达与纯Air同样的10$^{13}$cm$^{ - 3}$量级.在与纯Air相同的电子密度和碰撞频率条件下,采用95%Ar+5%Air和90%Ar+10%Air两种混合气,激波后2区实验时间和气体长度约为纯Air条件下的5$\sim $10倍,其中2区实验时间为300$\sim $800 $\mu$s,2区气体长度1$\sim $1.5 m.在$\varPhi $800 mm激波管中采用Ar+Air介质气体进行电磁波传输实验,获得了比在纯Air介质中与理论预测更一致的结果.

2018, 50(6): 1328-1366. doi:10.6052/0459-1879-18-284

HTML PDF

脉冲激光与正激波相互作用过程和减阻机理的实验研究

王殿恺, 文明, 王伟东, 卿泽旭

纳秒脉冲激光具有峰值功率密度高、易于击穿空气形成等离子体这一突出优势,在降低超声速波阻方面具有重要应用价值.以深刻揭示减阻机理为目的,针对激光与正激波相互作用这一基本物理现象开展实验研究.发展高精度纹影技术以测量复杂激波结构,时间分辨率达到 30ns,空间分辨率达到 1mm;搭建快速~PIV 实验系统以定量测量流场速度和涡量,时间分辨率达到 500ns.探明了激光等离子体引致的球面激波和高温低密度区域特性,揭示了激光等离子体在正激波冲击下的流动特性与演化规律,并结合数值模拟结果阐明了脉冲激光等离子体降低超声速波阻的根本原因.研究表明：激光等离子体引致激波的初始马赫数随着激光能量而增大,形状由水滴形逐渐发展为球面形,传播速度随着时间降低,在50$\mu$s 后接近于声速;高温低密度区域初始近似于球形,而后从激光入射方向的下游开始失稳,形成尖刺结构;在正激波冲击下,高温低密度区域演化为上下对称的双涡环结构,尺寸随着激光能量而增大.涡的卷吸和逆流可改变飞行器头部激波结构,是流场重构的重要形式,引起飞行器表面压力的大幅降低,是引起超声速飞行器波阻降低的重要机理.

2018, 50(6): 1337-1345. doi:10.6052/0459-1879-18-104

HTML PDF

内埋武器高速风洞弹射投放模型试验关键技术研究

宋威, 鲁伟, 蒋增辉, 白鹏

针对新一代战斗机超声速内埋武器弹射投放分离安全性问题,采用高速风洞投放实验技术研究内埋武器从开式武器舱弹射投放分离动态运动过程,风洞投放模型试验过程中采用除垂直加速度不足外,其余全部运动严格相似的轻模型相似设计方法,并针对轻模型法垂直加速度不足所导致的投放垂直位移偏离实物位移问题,采用一种简单易行的公式修正法进行补偿,试验给出了不同初始弹射投放分离条件下,内埋武器从载机投放分离后运动轨迹与姿态角随分离时间的变化规律,试验马赫数$Ma = 1.5$.研究结果表明：初始投放分离角速度对内埋武器投放分离后的运动轨迹及姿态角有较大的影响,当初始投放分离角速度$\omega _{z0}^s = 0^\circ/{\rm s}$时,内埋导弹出舱后先向下运动远离载机的流场干扰区,之后逐渐向载机方向抬升靠近并最终碰撞载机,高速风洞投放试验结果是不安全的,但经过公式修正后投放试验结果比较乐观,垂直方向运动仍然一直下降远离载机,这说明采用高速风洞投放试验得出的导弹不安全投放分离对真实载机来说不一定会出现,高速风洞投放试验结果比较保守. 当初始投放分离角速度$\omega _{z0}^s= 15^\circ/{\rm s}$和$\omega _{z0}^s = 30^\circ/{\rm s}$时,内埋导弹投放分离后运动趋势几乎一致, 均没出现向载机靠近的现象,内埋导弹具有一定的初始投放分离角速度有利于内埋武器的安全分离.

2018, 50(6): 1346-1355. doi:10.6052/0459-1879-18-180

HTML PDF

各向同性湍流通过正激波的演化特征研究

洪正, 叶正寅

激波与湍流相互作用(shock-turbulence interaction,STI)是空气动力学研究中的一个基础问题.基于格心有限差分法(cell-centered finite difference method,CCFDM)求解器Helios,采用五阶加权紧致非线性格式(weighted compact nonlinear scheme,WCNS)对各向同性湍流通过正激波的情形进行直接数值模拟(direct numerical simulation,DNS).对湍流相关物理量进行统计,分析结果表明,在湍流中波后的密度、温度和压力较无湍流情形下略小,而速度则略大,均在波后呈现短暂过冲然后缓慢向理论值逼近的变化趋势;波后流向雷诺应力突降随之快速增长又衰减,呈现非单调变化趋势,线性相互作用分析(linear interaction analysis,LIA)将其归结为波后能量从声模式转移为涡模式方式,与流向不同,横向雷诺应力突增后单调衰减,波后雷诺应力各向异性明显且随下游距离逐渐增强;波后湍动能突增后呈现非单调变化趋势;泰勒微尺度和Kolmogorov尺度过激波后均明显减小,说明波后湍流长度尺度变小,从而对波后网格的分辨率提出了更高的要求;密度、温度和压力过激波后脉动均方根均增加,密度和压力脉动强度减小,温度脉动强度增大.

2018, 50(6): 1356-1367. doi:10.6052/0459-1879-18-129

HTML PDF

高超声速边界层中模态转化的数值研究

高军, 李佳

在高超声速边界层中,第一模态和第二模态是与转捩有关的两个主要不稳定模态.除了不稳定模态,还存在一类稳定模态,其相速度在前缘接近快声波的相速度称为快模态.在感受性过程中,这类模态对激发边界层中不稳定模态起着很重要的作用.前缘感受性理论解释了边界层外扰动激发边界层中第一模态波的机理.针对高超声速平板边界层,利用相似性解剖面作为基本流,采用线性稳定性理论和直接数值模拟的方法研究了快模态和慢模态的稳定性行为.研究发现模态转化的位置与马赫数有关.根据线性稳定性理论的结果定义了临界频率.当扰动频率高于临界频率,第一模态与第二模态同支;而当扰动频率低于临界频率,第一模态与第二模态的共轭模态同支.借助稳定性方程的伴随方程分析了直接数值模拟的结果.直接数值模拟结果表明不论上游是快模态还是慢模态,当它们经过第二模态的不稳定区,它们都会演化成第二模态. 这可用模态在非平行流中传播的特征来解释.

2018, 50(6): 1368-1378. doi:10.6052/0459-1879-18-260

HTML PDF

聚变堆相关的液态金属膜流初步实验研究

倪明玖, 阳倦成, 任东伟, 刘佰奇, 齐天煜, 胡建生, 李建刚

在磁约束核聚变堆的面对等离子部件设计中,液态金属锂膜流因具有带走杂质、保护面对等离子固壁等优点而被认为是优选方案之一. 然而,如何克服聚变堆中强磁场环境下产生的磁流体力学效应并形成大面积均匀铺展锂膜流动是目前亟需解决的问题.本文通过搭建室温液态镓铟锡回路和高温液态锂回路,开展了两种不同特性的液态金属膜流实验, 并采用传统可视化方法获得了展向磁场存在时镓铟锡和锂在导电底板形成的液膜流动表面特征.实验结果表明: 无磁场时,两种液态金属膜流流动表面波动特性与常规流体膜流均一致, 即随着流动雷诺数的增加表面波动变得更为混乱; 而展向磁场存在时,镓铟锡膜流表面波动变得更为规则, 且沿着磁场方向平行排列,表现为拟二维波动的特征; 而锂膜流却产生了明显的磁流体力学阻力效应,表现为在流动方向局部产生锂滞留现象, 且滞留点随雷诺数增大向下游移动. 最后通过膜流受力分析,进一步阐述了锂膜流受到比镓铟锡膜流更为严重磁流体力学效应影响的原因.

2018, 50(6): 1379-1386. doi:10.6052/0459-1879-18-367

HTML PDF

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

毛洁, 王彦利, 王浩

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明：沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.

2018, 50(6): 1387-1395. doi:10.6052/0459-1879-18-276

HTML PDF

等离子体激励器控制圆柱绕流的实验研究

张鑫, 黄勇, 李华星

为了发展新型移动附面层控制技术,提升流动控制效率,采用粒子图像测速技术,开展了基于对称布局等离子体气动激励的圆柱绕流控制研究,获得了静止空气下,对称布局激励器诱导流场的演化过程,评估了来流条件下等离子体控制效果,通过等离子体诱导涡实现了虚拟移动附面层控制,分析了诱导涡随时间演化的过程,揭示了圆柱绕流等离子体控制机理.结果表明：(1)在静止空气下,对称布局激励器在刚启动瞬间,会在暴露电极两侧诱导产生一对旋转方向相反的启动涡;随着时间的推移,启动涡逐渐向远离壁面的方向运动;随后,激励器在暴露电极两侧产生了两股速度近似相等,方向相反的诱导射流,诱导射流在柯恩达效应的影响下,朝壁面方向发展.(2)当激励电压峰峰值为19.6 kV,激励频率3kHz时,施加等离子体气动激励后,圆柱脱落涡得到了较好抑制,圆柱阻力系数减小了21.8％;(3)在来流作用下,对称布局激励器在靠近来流一侧,诱导产生了较为稳定的涡结构.诱导涡通过旋转、运动,促进了壁面附近低能气流与主流之间的掺混,抑制了圆柱绕流流场分离,实现了"虚拟移动附面层控制"效果.与传统移动附面层控制技术相比,基于等离子体气动激励的新型移动附面层控制技术不需要复杂、笨重的机构,不会带来额外的阻力,具有潜在的应用前景.

2018, 50(6): 1396-1405. doi:10.6052/0459-1879-18-279

HTML PDF

Plotnikov--Toland板模型中水弹性孤立波的迎撞

巴迪M.M., 卢东强

通过奇异摄动方法研究了在薄冰层覆盖的不可压缩理想流体表面上传播的两个水弹性孤立波之间的迎面碰撞.借助特殊的 Cosserat 超弹性壳理论以及Kirchhoff--Love 板理论,冰层由 Plotnikov--Toland板模型描述.流体运动采用浅水假设和Boussinesq 近似. 应用Poincaré--Lighthill--Kuo 方法进行坐标变形,进而渐近求解控制方程及边界条件, 给出了三阶解的显式表达. 可以观察到碰撞后的孤立波不会改变它们的形状和振幅. 波浪轮廓在碰撞之前是对称的, 而在碰撞之后变成不对称的并且在波传播方向上向后倾斜. 弹性板和流体表面张力减小了波幅. 图示比较了本文与已有结果可知线性板模型可作为本文的一个特例.

2018, 50(6): 1406-1417. doi:10.6052/0459-1879-18-287

HTML PDF

奇异摄动在层流预混火焰理论研究中的应用

陈正

奇异摄动被广泛应用于求取力学问题的近似解.一个典型问题就是流体力学中的边界层.郭永怀先生曾发展了适用于平板黏性流动边界层问题的奇异摄动理论.类似于流体力学中的边界层,燃烧研究中的层流预混火焰也可以通过奇异摄动理论进行分析,在燃烧研究中通常称其为大活化能渐近分析.本文介绍了大活化能渐近分析在一维平面预混火焰和球形传播火焰中的应用及相关研究进展.首先介绍了预混火焰结构及其涉及的不同特征尺度,分析了大活化能条件下出现的特征尺度分离,并给出了关于平面预混火焰大活化能渐近分析的详细推导,讨论了热辐射对火焰传播的影响;然后介绍了大活化能渐近分析在点火与球形传播火焰方面的应用,指出了只有能够同时描述点火与球形火焰传播的理论才能准确地预测临界点火条件,并讨论了考虑链式反应的点火与火焰传播理论,分析了热辐射对球形火焰传播的影响,给了关于火焰稳定性理论研究的发展趋势.最后,基于当前研究进展对未来的研究方向进行了展望,其中涉及多步化学反应、低温冷火焰、复杂流动、辐射重吸收等.

2018, 50(6): 1418-1435. doi:10.6052/0459-1879-18-243

HTML PDF

多孔介质内黏弹性流体的热对流稳定性研究

康建宏, 谭文长

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果：(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.

2018, 50(6): 1436-1457. doi:10.6052/0459-1879-18-309

HTML PDF

三维方腔介电液体电对流的数值模拟研究

吴健, 张蒙齐, 田方宝

本文对封闭方腔内介电液体电对流进行了三维数值模拟研究.方腔的6个边界为固壁;4个侧边界为电绝缘边界;上下界面为两个电极.直流电场作用在从底部电极注入的自由电荷上,从而对液体施加库伦体积力并驱动流体流动形成电对流.为了求解这一物理问题,发展了一种二阶精度的有限体积法来求解完整的控制方程,包括Navier-Stokes方程和一组简化的Maxwell方程.考虑到电荷密度方程的强对流占优特性,采用了全逆差递减格式来求解该方程,获得了准确有界的解.通过研究发现,该流动在有限振幅区内的分叉类型为亚临界,即系统存在一个线性和非线性临界值,分别对应流动的开始和终止.由于非线性临界值比线性值小,因此两个临界值之间有一个迟滞回线.与无限大域中的自由对流相比,侧壁施加的额外约束改变了流场结构,使这两个临界值均有所增大.此外,还讨论了电荷密度和速度场的空间分布特征,发现电荷密度分布中存在电荷空白区.最后对更小空间尺寸情况计算结果表明,流动的线性分叉类型为超临界.本文的结果拓展了已有的二维有限空间内电对流的研究,并为三维电对流的线性和弱非线性理论分析提供参考.

2018, 50(6): 1458-1469. doi:10.6052/0459-1879-18-301

HTML PDF

基于非结构/混合网格模拟黏性流的高阶精度DDG/FV混合方法

邵帅, 李明, 王年华, 张来平

间断Galerkin有限元方法(discontinuous Galerkin method, DGM) 因具有计算精度高、模板紧致、易于并行等优点, 近年来已成为非结构/混合网格上广泛研究的高阶精度数值方法. 但其计算量和内存需求量巨大, 特别是对于网格规模达到百万甚至数千万的大型三维实际复杂外形问题, 其计算量和存储量对计算资源的消耗是难以承受的. 基于“混合重构”的DG/FV 格式可以有效降低DGM 的计算量和存储量. 本文将DDG 黏性项离散方法推广应用于DG/FV 混合算法, 得到新的DDG/FV混合格式, 以进一步提高DG/FV混合算法对于黏性流动模拟的计算效率. 通过Couette流动、层流平板边界层、定常圆柱绕流, 非定常圆柱绕流和NACA0012 翼型绕流等二维黏性流算例, 优化了DDG 通量公式中的参数选择, 验证了DDG/FV 混合格式对定常和非定常黏性流模拟的精度和计算效率, 并与广泛使用的BR2-DG 格式的计算结果和效率进行对比研究. 一系列数值实验结果表明, 本文构造的DDG/FV混合格式在二维非结构/混合网格的Navier-Stokes 方程求解中, 在达到相同的数值精度阶的前提下, 相比BR2-DG格式, 对于隐式时间离散的定常问题计算效率提高了2 倍以上, 对于显式时间离散的非定常问题计算效率提高1.6 倍, 并且在一些算例中, 混合格式具有更优良的计算稳定性. DDG/FV 混合格式提升了计算效率和稳定性, 具有良好的应用前景.

2018, 50(6): 1470-1482. doi:10.6052/0459-1879-18-199

HTML PDF

留言板

2018年第50卷第6期

通知公告MORE+

精品文章MORE+

下载中心

作者中心

留言板

2018年 第50卷 第6期

通知公告MORE+

精品文章MORE+

下载中心

作者中心

2018年第50卷第6期