快速检索         高级检索
2020年  第42卷  第2期      刊出日期 2020-04-28   目次   上一期   下一期  
     
封面和目次
专题:高性能计算
  • 大规模科学数据体绘制技术综述
  • 王华维,何柳,曹轶,肖丽
  • 2020, 42 (2):1-12[摘要](5522)  [PDF  8173KB](4382) [HTML](41)
    体绘制是刻画大规模科学数据中复杂物理特征的有效途径,然而,数据量极大、特征难以捕捉等问题依然是目前体绘制研究的主要挑战。为此,研究者们从三个方面对体绘制算法进行了深入研究,以提高大规模数据体绘制的效率和效果。一方面,依托硬件通过多处理器核来分担计算,降低单处理器核所要完成的计算量,是提高体绘制效率的一个有效途径。另一方面,充分发掘数据场内在特性对三维数据场进行约简,大幅减少绘制处理数据量从而降低算法开销,也是提高体绘制效率的一个有效途径。同时,在体绘制算法中融入特征分析和特征增强方法,让复杂物理特征从数据场中突显出来,以实现对科学数据的高质量绘制。本文对国内外体绘制技术相关研究进展进行了调研、综述,并分析了不同的研究方法,最后展望了未来体绘制技术研究的可能发展方向,包括应用驱动的特征体绘制、基于特征的约简体绘制、适应硬件的体绘制多级加速以及原位智能化体绘制等。
  • 面向图计算应用的处理器访存通路优化设计与实现
  • 张旭,常轶松,张科,陈明宇
  • 2020, 42 (2):13-22[摘要](6097)  [PDF  6860KB](4221) [HTML](41)
    针对图计算应用的访存特点,提出并实现一种支持高并发、乱序和异步访存的高并发访存模块(High Concurrency and high Performance Fetcher, HCPF)。通过软-硬件协同的设计方法,HCPF可同时处理192条共8种类型的内存访问请求,且访存粒度可由用户定义,满足图计算应用对海量低延迟细粒度数据访问的需求。同时,HCPF扩展了基于内存语义的跨计算节点定制互连技术,支持远程内存的细粒度直接访问,为后续实现分布式图计算框架提供技术基础。结合上述两个核心研究内容,基于流水线RISC-V处理器核,设计并实现了可支持HCPF的RISC-V片上系统(System-on-Chip,SoC)架构,搭建基于FPGA的原型验证平台,并使用自研测试程序对HCPF进行初步性能评测。实验结果表明,HCPF相比原有访存通路,最高可将基于数组和随机地址的两种随机内存访问性能分别提升至3.5倍和2.7倍。远程内存直接访问4 Byte数据的延时仅为1.63 μs。
  • 面向大规模集群的并行I/O用户层配置优化策略
  • 田鸿运,武林平,董勇,景翠萍,罗红兵,莫则尧
  • 2020, 42 (2):23-30[摘要](5405)  [PDF  5706KB](3870) [HTML](38)
    影响应用I/O性能的关键因素主要有三个层次:包括应用的I/O接口实现、体系结构和文件系统组件的性能以及应用的I/O参数配置。从应用I/O配置优化的视角,分析了大规模集群并行I/O的配置调优空间,在此基础上,给出了一套大规模集群并行I/O性能特征测试分析方法。基于该方法,在某国产超级计算集群上开展了一系列I/O测试分析来刻画系统的I/O性能特征,进而指导并行应用程序的I/O配置优化。基于优化后的配置参数,在两类典型的并行I/O场景中,针对某类生产应用程序,8192进程下的重启动数据写操作时间下降了15%,4096核的程序作业加载时间从10 min缩短到了5 s。
  • 高精度CFD程序的内外子区域划分异构并行算法
  • 王巍,徐传福,车永刚
  • 2020, 42 (2):31-40[摘要](7417)  [PDF  7495KB](3970) [HTML](37)
    对计算流体力学(Computational Fluid Dynamics, CFD)程序CNS提出一种Offload模式下对任务内外子区域划分的异构并行算法,结合结构化网格下有限差分计算和四阶龙格-库塔方法的特点,引入ghost网格点区域,设计了一种ghost区域收缩计算策略,显著降低了异构计算资源之间的数据传输开销,负载均衡时CPU端的计算与MPI通信完全和加速器端的计算重叠,提高了异构协同并行性。推导了保证计算正确性的ghost区域的参数,分析了负载均衡的条件。在“CPU(Intel Haswell Xeon E5-2670 12 cores ×2)+加速器(Xeon Phi 7120A ×2)”的服务器上测得该算法较直接将任务子块整体迁至加速器端计算的异构算法性能平均提升至5.9倍,较MPI/OpenMP两级并行算法使用24个纯CPU核的性能,该算法使用单加速器时加速至1.27倍,使用双加速器加速至1.45倍。讨论和分析了性能瓶颈与存在的问题。
  • 流水的浮点倒数近似值运算部件的设计与实现
  • 何军,王丽
  • 2020, 42 (2):41-46[摘要](5065)  [PDF  4142KB](3864) [HTML](43)
    在部分低精度浮点运算应用中,需要流水的浮点倒数近似值运算。本文基于SRT-4算法设计并实现了一种流水的浮点倒数近似值运算部件。该部件采用6级流水线结构,运算结果精度至少为8位有效尾数。为了支持对非规格化浮点数的硬件处理,还设计并实现了改进版,有利于进一步提高浮点倒数近似值运算的性能。改进版采用8级流水线结构,新增了源操作数预规格化和结果后规格化功能模块,可以实现对非规格化浮点数的硬件处理。经过逻辑综合评估,改进版的硬件开销是面积在合理范围内增加19.23%,且对时序没有明显影响,可以满足预期的1.6 GHz频率设计目标。
  • MPI并行程序中通信等待问题的诊断方法及其应用
  • 武林平,景翠萍,刘旭,田鸿运
  • 2020, 42 (2):47-54[摘要](4883)  [PDF  6214KB](4073) [HTML](37)
    随着并行规模的扩大,现有通信等待问题的诊断方法存在内存开销大、测量时间开销大等问题。通过对现有通信等待问题诊断方法的深入分析,同时考虑测量开销可控的实际需求,建立基于热点函数的通信等待问题诊断模型。基于上述模型,总结出一种更精简、更实用的通信等待问题诊断方法。将该诊断方法分别应用到二维LARED集成、LARED-S、LAP3D等大规模MPI并行程序的通信等待问题诊断过程,应用效果表明本诊断方法可精确定位导致通信等待问题的关键代码段,给出的优化方案及性能提升空间对于后续的程序改进具有参考价值,其中根据诊断结果优化后的LARED-S程序性能提升32%,通信等待时间减少44%。
  • 高铁大风预警模式挖掘
  • 滕飞,刘鉴竹,祝锦烨,勾红叶
  • 2020, 42 (2):55-63[摘要](5706)  [PDF  5757KB](3943) [HTML](37)
    高铁大风预警的传统方法基于风速预测,当瞬时值高于限速阈值时触发报警,存在大量的误报警,不必要的限速控制影响了高铁行车效率。创新地提出了基于序列模式的预警方法,旨在挖掘报警事件前序数据中的频繁模式,找出报警事件的变化规律,通过滤除与非预警序列共有的频繁模式,得到预警序列独有的序列特征,构建了预警模式库。经兰新高铁沿线的监测数据验证,该方法在提高预测准确率的基础上降低了漏报率,同时有效地减少了模式匹配所需的时间,为提前预警预留充分的时间窗口,更加符合实际应用的需求。
  • 异构混合存储的软硬件协同数据放置策略
  • 李鸿飞,杜溢墨,曾熠,王磊
  • 2020, 42 (2):64-71[摘要](6014)  [PDF  6086KB](4169) [HTML](39)
    分析比较当前大数据中心典型混合存储架构,针对其不能综合运用存储管理系统和存储设备优势的问题,提出软硬件协同的数据放置策略,同时考虑软件层混合存储管理系统和硬件层混合存储设备的特点,根据应用特性为数据选择合适的存储管理系统和设备。面向不同应用场景,提出运行前规划存储路径的静态放置模式和运行中规划存储路径的动态放置模式。基于存储管理系统和设备性能参数建模,采用仿真方法实现各数据放置策略,并运行实际应用中的三类负载进行实验测试,结果表明采用软硬件协同的数据放置策略相比只考虑存储管理系统和存储设备的数据放置策略,性能可以提高近30%。
航天工程·材料科学与工程
  • 二维可压缩流动抽吸流量模型及验证
  • 张宝虎,赵玉新,柳军
  • 2020, 42 (2):72-77[摘要](5967)  [PDF  4183KB](4350) [HTML](34)
    根据超声速与亚声速流动的特征,依据Prandtl-Meyer膨胀波理论,根据背压与来流静压确定膨胀波后的马赫数,获得了质量抽吸流率,并导出了壅塞发生的条件,建立了二维超声速抽吸模型;根据冲量定理,假定气体从主流区域到抽吸腔流态的变化是由于压差的冲量作用,建立了亚声速抽吸模型,并且进一步阐述了亚声速抽吸壅塞发生的机制。计算流体力学验证表明亚声速抽吸冲量模型的偏差较大,因此通过数值拟合模型来描述亚声速抽吸流动。该研究为三维抽吸研究提供了参考。
  • 导弹结构热防护一体化设计与优化
  • 徐世南,吴催生
  • 2020, 42 (2):78-84[摘要](6355)  [PDF  5233KB](4334) [HTML](39)
    针对导弹面临的严酷气动热力环境,基于传统波纹夹芯结构,为导弹设计了一种兼具承载和热防护功能的一体化结构。建立热-结构耦合计算模型,对一体化结构进行结构承载与热防护性能分析,结果表明:一体化结构在高温下具有较好的承载能力;增加一体化结构隔热层与内壁厚度有助于提高导弹的热防护效率。以一体化结构热防护效率最高和质量最小为目标,采用第二代非支配排序遗传算法实现了多目标优化,并认识到:一体化结构热防护效率与质量相互影响,设计时需要需综合考虑。
  • 基于位置预测的靶场图像实时判读方法
  • 钟立军,于起峰,周颉鑫,郭鹏宇,黄维
  • 2020, 42 (2):85-91[摘要](6489)  [PDF  5487KB](4077) [HTML](37)
    在靶场经纬仪对目标实时跟踪测量时,会发生相机随机抖动的情况,引起目标在图像中大幅度运动。应对大幅度运动时,基于搜索窗口的跟踪方法容易丢失目标,而基于全图搜索的跟踪方法时效性差。针对以上问题,提出一种结合核相关滤波算法(Kernelized Correlation Filter, KCF)和目标位置预测的改进的跟踪学习检测算法(Tracking-Learning-Detection, TLD)跟踪框架。利用正交多项式最优线性滤波器及相机角度信息预测目标下一帧位置,在此区域利用KCF进行快速跟踪,可以提高跟踪的成功率和时效性,跟踪失败时再进行检测。仿真实验表明,最优线性滤波器能较准确预测目标位置,给KCF提供较准确的搜索位置,算法每帧耗时仅为1.1 ms,且定位精度优于TLD和KCF,能有效应对相机抖动的问题。靶场实际试验证明该方法可提高靶场自动判读水平,减少人工干预。
  • 混合不确定条件下航空不安全事件量化分析方法
  • 崔利杰,陈浩然,任博,张贾奎
  • 2020, 42 (2):92-97[摘要](5586)  [PDF  4078KB](3938) [HTML](36)
    为解决不确定条件下航空不安全事件量化分析的难题,基于随机和模糊理论对影响飞行安全的主客观变量进行了混合不确定描述,将难以用概率表示的主观变量处理成为隶属度λ水平截集上服从均匀分布的随机变量。在此基础上,提出了混合不确定条件下航空安全指标和数值计算方法,所提指标和方法能够简洁、直观、准确地反映航空安全水平。通过航空安全领域普遍采用的Bow-tie模型对轮胎爆破事件安全性指标和隶属度关系进行仿真分析。结果表明: 将不确定变量控制在一定隶属度水平,能显著提高航空安全指标。
  • 空间双体系绳系统展开阶段末端星姿态动力学分析
  • 陈书敏,王长青,扎伯罗特诺夫·尤里,李爱军
  • 2020, 42 (2):98-106[摘要](5506)  [PDF  6225KB](3948) [HTML](36)
    为了研究空间系绳系统展开过程中末端星的姿态运动,采用第二类拉格朗日方程建立系统展开及末端星角运动的数学模型,分析展开阶段末端星的姿态动力学特性。建模时将系绳末端连接的母星与子星视为尺寸不可忽略的刚体,且母星质量远大于子星质量,系绳视为有质量的刚性杆。利用该数学模型可以分析系绳系统的展开过程,研究影响末端星姿态运动的主要因素,包括姿态角初始扰动以及末端星动/静不对称性。仿真结果表明,系绳展开时末端星姿态角初始扰动及动静不对称性可能会导致姿态角失稳,出现系绳缠绕星体的情况。仿真结果可以为展开阶段末端星姿态控制提供参考。
  • 弹(箭)体级间连接结构准静载失效实验与数值预示
  • 袁杰红,田彤辉,关振群,王青文,陈柏生
  • 2020, 42 (2):107-113, 126[摘要](5362)  [PDF  6749KB](3817) [HTML](45)
    级间螺栓法兰连接是导弹(火箭)常见的连接方式,但破坏了整体结构连续性,且承载能力薄弱,在外荷载作用下易发生失效而使整弹(箭)结构强度丧失。根据实际导弹(火箭)连接结构特点,简化设计并制作了一组原理性实验件,设计了准静载加载实验测试系统和螺栓响应信号传感器,进行了两次准静载失效实验,并利用ABAQUS软件建立了对应的有限元仿真模型。根据实验效果和实测数据,分析了连接结构在准静载荷载作用下的失效机理,并对比验证发现有限元模型数值模拟效果和精度与实验吻合较好。研究结论可为弹(箭)体级间连接结构承载能力和失效实验设计提供参考。
  • 利用预报误差分布拟合实现卫星历史轨道机动检测的方法
  • 李涛,黄昊,陈磊
  • 2020, 42 (2):114-120[摘要](6311)  [PDF  4993KB](4143) [HTML](43)
    轨道机动会引起在轨卫星的轨道异常,是空间态势感知的重点关注对象之一。为此,提出一种从卫星的两行轨道要素编目数据中检测历史机动的方法。通过分析预报误差检测轨道参数的异常编目值,进而获得相应的历史机动信息。对卫星的两行轨道要素编目数据两两预报求差得到预报误差的样本数据后,利用期望最大化算法可以拟合得到以高斯混合模型表达的预报误差的概率分布模型;以此为基础,可以确定利用预报误差检测轨道参数异常编目值的门限;利用轨道机动与异常编目值之间的对应关系,最终实现机动检测。对典型活动卫星的机动检测结果表明,所提方法能够在准确检出历史机动事件的同时保持较低的误检率。
  • 采用优化加点Kriging模型的助推火箭残骸安全区预示方法
  • 祝学军,卜奎晨,王浩,高峰,赵长见
  • 2020, 42 (2):121-126[摘要](5165)  [PDF  4210KB](3980) [HTML](35)
    建立火箭及其分离残骸弹道计算动力学模型,并采用四元数方法对姿态角解算进行处理。提出基于优化加点Kriging模型的安全区预示方法,结合Monte Carlo和Kriging代理模型的特点,给出安全区预示流程。以某型助推火箭残骸安全区计算为例,对提出的安全区预示方法进行仿真验证。仿真结果表明,提出的基于优化加点Kriging模型安全区预示方法与Monte Carlo方法相比,在不损失计算精度的前提下,具有更高的计算效率,满足快速迭代的工程需求;相比传统极限偏差叠加方法,可显著降低安全区覆盖面积,具有较强的工程应用价值。
机械工程·控制科学与工程
  • 多重优化的分布式无线覆盖探测算法
  • 王红军,向庭立,潘继飞
  • 2020, 42 (2):127-134[摘要](5859)  [PDF  5707KB](3941) [HTML](39)
    为满足无线通信网络信号覆盖有效性的实时实地可重复探测的需求,提出一种基于传感器网络的分布式无线覆盖探测算法。通过随机部署于目标区域内的无线传感器节点对无线通信网接收信号强度进行感知和预处理;利用变异函数构造新的BP神经网络目标函数,通过改进粒子群算法优化其初始权值和阈值;利用训练好的网络模型对存在探测盲区的目标区域进行插值估计,并联合传感器节点采集到的数据生成无线通信网络等信号强度线。仿真结果表明,所提算法比其他经典算法具有更高的精度,可有效探测目标区域无线通信网络的信号覆盖情况。
  • 混响背景下基于高阶统计量的动目标检测方法
  • 王晓彤,蔡志明
  • 2020, 42 (2):135-141[摘要](4996)  [PDF  4804KB](3885) [HTML](44)
    针对混响背景中的动目标检测问题,根据基阵接收数据经过波束形成与匹配滤波后的输出结果计算高阶统计量,并将其视作观测空间。基于此空间中混响和目标回波的差异,利用多ping的高阶统计量构造特征向量,计算特征向量之间的马氏距离作为混响和目标差异的量化标准,再依据最大一致条件功效检测准则选择门限检测方法。波形数据仿真与海上实录数据检验均表明该方法的检测性能优于单ping波束形成及匹配滤波方法。通过蒙特卡洛仿真获得不同信混比下的接收机工作特性曲线,与单ping检测相比,在保证虚警概率小于0.01、检测概率大于0.5的条件下,最小可检测信混比降低约6 dB。
  • 基于统一理论框架的故障诊断模型及学习进化方法
  • 张士刚,罗旭,骆彦廷,杨拥民
  • 2020, 42 (2):142-149[摘要](6463)  [PDF  5725KB](4014) [HTML](62)
    以概率图理论为基础,系统研究基于这一理论框架的故障诊断模型,对模型的构建方法以及在不同场景下的模型演化方案进行探讨,使得在统一理论框架下可实现多模式系统故障、耦合故障、动态故障、故障预测等复杂情形的诊断。为了弥补单独利用基于模型的方法和基于数据的方法的缺陷,研究了诊断模型的学习进化策略,实现了诊断效果的改进和优化。对模型后续的能力扩展和可能的研究方向进行了展望,为后续理论研究提供了参考。
  • 舰炮身管的模态分析与多目标优化
  • 刘国强,陈维义,陈华东,程晗
  • 2020, 42 (2):150-155[摘要](5795)  [PDF  4163KB](3971) [HTML](34)
    舰炮的炮口扰动是影响舰炮射击精度的重要因素。为了减少舰炮的炮口扰动,优化身管结构,建立了柔性身管的ABAQUS有限元模型,将有限元模型的模态计算值与模态试验值进行比较,发现建立的柔性身管有限元模型是合理可行的。提出了混沌量子粒子群算法与动力学联合优化的方法,进行身管和炮口的多目标优化。优化结果表明,优化后炮口中心的线速度、角速度和角位移与优化前相比显著减小,身管质量有所降低,结构更加合理,该优化方法有效可行,为下一步全炮总体优化设计提供了一定的借鉴。
  • 4电平S类数字功放的电平合成策略
  • 周强,朱蕾,陈剑斌
  • 2020, 42 (2):156-161[摘要](4959)  [PDF  4501KB](3761) [HTML](32)
    通过增加量化电平,Delta-Sigma调制器输出序列能够获得更好的调制性能。但在S类数字功放中,受限于开关功放的实现结构,多电平Delta-Sigma调制器输出脉冲序列无法直接驱动开关功放,需采用由多个开关功放单元构成的多电平开关功放。针对2比特Delta-Sigma调制器输出序列,分析了电平合成策略对4电平S类数字功放性能的影响,提出了两种优化的电平合成策略。理论分析和仿真结果表明,采用合适的电平合成策略,4电平S类数字功放可以获得更优的输出性能,尤其在输出功率回退状态下,电平合成策略对功放性能的影响更为显著。
  • ASLR机制脆弱性自动分析方法
  • 黄宁,黄曙光,潘祖烈,常超
  • 2020, 42 (2):162-170, 185[摘要](5694)  [PDF  6419KB](4173) [HTML](43)
    地址随机化是一种针对控制流劫持漏洞的防御机制。已有的漏洞自动分析与利用技术缺少对地址随机化机制影响的分析,导致生成的测试用例在实际环境中的运行效果受到极大限制。针对地址随机化的缺陷及其绕过技术的特点,提出了一种地址随机化脆弱性分析方法。该方法使用有限状态机描述程序运行路径中各关键节点的状态;针对常见的内存泄漏与控制流劫持场景建立约束条件;通过求解内存泄漏状态约束与控制流劫持状态约束的兼容性,分析地址随机化机制在特定场景下的脆弱性。实验结果表明,该方法可有效检测通过内存泄漏导致的地址随机化绕过及控制流劫持攻击,实现自动化的地址随机化脆弱性分析,提高针对软件安全性分析的效率。
  • 数据基故障诊断算法更新问题研究
  • 赵晨旭,涂遗,邱静,刘冠军
  • 2020, 42 (2):171-176[摘要](6025)  [PDF  4366KB](3807) [HTML](40)
    机内测试被广泛应用于故障诊断、装备健康管理与预测等领域。针对机内测试设备在设计和升级时遇到的分类器更新、样本数量不平衡、硬件条件限制问题,提出初步解决方案。利用基于密度的聚类和人工免疫方法处理原始数据;提出基于代表样本点的混合学习方法;利用支持向量机和仿真案例验证所提方法。结果表明所提方法能够解决上述问题,有助于基于数据的机内测试设备设计与升级。
  • 考虑非线性退化与随机失效阈值的剩余寿命预测
  • 王泽洲,陈云翔,蔡忠义,罗承昆
  • 2020, 42 (2):177-185[摘要](5126)  [PDF  6506KB](3968) [HTML](36)
    针对广泛存在的非线性退化设备,现有方法尚未考虑随机失效阈值对剩余寿命预测结果的影响。因此,通过对设备性能退化过程进行分析,提出了一种综合考虑非线性退化与随机失效阈值的剩余寿命预测方法。基于Wiener过程构建了考虑个体差异与测量误差的非线性退化模型;基于卡尔曼滤波算法建立状态空间模型以实现对退化状态的在线更新;基于极大似然法估计失效阈值分布系数估计方法,得到随机失效阈值的概率分布;基于随机失效阈值推导出剩余寿命的概率分布,实现对剩余寿命的在线预测。算例研究表明,所提方法可以有效地提升剩余寿命预测的准确性,具备一定工程应用价值。
  • 面向作战任务的最小风险传感器调度方法
  • 单甘霖,庞策,段修生
  • 2020, 42 (2):186-193[摘要](5238)  [PDF  5912KB](3797) [HTML](37)
    基于风险理论提出一种面向不同目标探测任务的传感器调度方法。将主动传感器辐射被截获风险和目标探测风险结合起来,建立一般目标探测框架下的传感器最小风险调度模型。分目标跟踪、目标识别和目标威胁等级评估三种情况将传感器最小风险调度模型具体化,给出不同情况下目标探测风险值的计算方法。针对模型的求解提出一种基于混沌思想、反向学习和双向轮盘赌的改进人工蜂群算法。通过仿真实验证明了模型的可行性和算法的有效性。
  • 采用改进遗传算法的舰载机保障调度方法
  • 刘珏,王能建,罗旭,孟祥雷
  • 2020, 42 (2):194-205[摘要](5289)  [PDF  9439KB](3992) [HTML](40)
    舰载机保障作业过程受到多种资源约束,同时可能存在突发事件的干扰。针对此特点,分析并制定干扰情况下的模型修正策略。在遗传算法中引入禁忌搜索算子改进传统遗传算法的变异操作,并通过具体实例进行仿真验证,其结果证明改进后的遗传算法在优化含干扰事件的多机保障问题时效率更高,并通过甘特图直观地反映重调度方案,为真实情况下有效处理舰载机保障过程中的干扰事件奠定理论基础。
  • 弹性轨道磁悬浮控制系统模型降阶研究
  • 刘武斌,王泉
  • 2020, 42 (2):206-210[摘要](5267)  [PDF  3158KB](3961) [HTML](37)
    在磁悬浮控制系统中,常常忽略轨道的弹性来设计控制算法。在这种控制算法作用下,当轨道刚度较小时,系统容易产生振动。为解决该问题,可以将轨道弹性加入悬浮模型,然后设计控制算法。考虑轨道弹性之后,悬浮系统的模型会比较复杂,控制算法难以在工程实现。为此,采用Hankel范数近似法对考虑轨道弹性后的模型进行降阶,并且在降阶模型的基础上设计控制算法,解决了轨道弹性引发的系统振动问题。并且,这种方法容易在工程中实现。文章最后利用仿真结果验证了降阶方法的可行性。
扫一扫关注公众号
国际标准连续出版物号
ISSN 1001-2486
国内统一连续出版物号
CN 43-1067/T
主管单位
国防科技大学
主办单位
国防科技大学教研保障中心
论文下载排行更多>>
论文点击排行更多>>
友情链接