期刊文献+
共找到15篇文章
< 1 >
每页显示 20 50 100
离心压气机进口径向总压畸变不确定性分析
1
作者 徐永祥 琚亚平 +3 位作者 秦瑞鸿 李震 刘艳 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2024年第8期69-79,共11页
针对离心压气机实际运行过程中的进口畸变不确定性问题,提出了一种离心压气机进口径向总压畸变不确定性参数化表征模型。以跨声速Radiver离心压气机为研究对象,采用非嵌入式多项式混沌法和Sobol敏感性分析方法,研究了压气机气动性能受... 针对离心压气机实际运行过程中的进口畸变不确定性问题,提出了一种离心压气机进口径向总压畸变不确定性参数化表征模型。以跨声速Radiver离心压气机为研究对象,采用非嵌入式多项式混沌法和Sobol敏感性分析方法,研究了压气机气动性能受叶尖畸变、轮毂畸变两类进口径向总压畸变不确定性的影响,获得了对压气机气动性能影响最为显著的不确定性因素,并分析了不确定性因素影响压气机气动性能的相关流动机理。结果表明:相比于均匀来流,叶尖畸变、轮毂畸变均使压气机级多变效率下降近1.0%;轮毂畸变、叶尖畸变分别对级总压比、扩压器总压损失系数带来的影响较为显著,级多变效率受轮毂畸变影响下降1.4%,叶尖畸变使损失系数增大11.0%;相较于轮毂侧和叶尖侧畸变峰之间的相对位置,压气机气动性能对畸变峰的数值大小更敏感;受进口径向总压畸变影响下的压气机的性能恶化与扩压器叶片进口冲角不均匀分布导致的流动分离现象密切相关。研究工作表明了进口畸变不确定性在离心压气机不确定性分析中的重要性,对未来发展离心压气机鲁棒优化设计具有一定的学术意义和工程应用价值。 展开更多
关键词 离心压气机 进口畸变 不确定性 敏感性分析 流动机理
下载PDF
多级轴流压气机末级串列静子三维优化设计与分析 被引量:1
2
作者 陈玲 琚亚平 +1 位作者 李震 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2023年第1期164-174,共11页
为充分挖掘串列叶片对压气机扩稳增效的潜能,在对多级轴流压气机末级静子串列改型的基础上,选取不同叶高处前后排叶片相对位置、负载分配、弯掠特征及端壁轮廓等25个参数为设计变量,开展某工业用6.5级轴流压气机末级串列静子叶片与端壁... 为充分挖掘串列叶片对压气机扩稳增效的潜能,在对多级轴流压气机末级静子串列改型的基础上,选取不同叶高处前后排叶片相对位置、负载分配、弯掠特征及端壁轮廓等25个参数为设计变量,开展某工业用6.5级轴流压气机末级串列静子叶片与端壁的一体化优化设计研究,并采用随机森林方法对串列叶片几何参数进行特征重要性分析,提炼出串列参数设计指导。结果表明:优化后压气机末级在设计工况、近喘工况和近堵工况下的绝热效率分别提升4.81%、3.37%和8.195%;对于末级串列静子,最优串列参数随叶高变化而变化;调整前后叶片的弦长比与弯角比以实现串列叶片负荷前加载、调整叶片后掠与反弯以降低叶片端壁处负荷、设置叶中部分更大的周向交错度使缝隙区流体充分加速,均能有效降低末级串列叶片流动损失;串列叶片缝隙区尺寸,尤其是周向交错度,较其他参数对串列叶片性能的影响更为显著。研究成果应对压气机三维串列叶片设计具有参考价值。 展开更多
关键词 轴流压气机 串列叶片 三维优化设计 随机森林 气动性能
下载PDF
多级跨声速轴流压气机几何参数化造型方法及增压特性分析 被引量:1
3
作者 关柳萍 琚亚平 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2023年第3期129-138,共10页
为发展高性能多级跨声速压气机先进设计体系,提出了压气机多圆弧叶型、多圆弧厚度分布叶型及修正NACA 6-系列叶型造型方法,利用NASA/GE研发的高效节能发动机E^(3)高压压气机前6.5级的准三维设计数据,建立了相应的压气机三维几何参数化模... 为发展高性能多级跨声速压气机先进设计体系,提出了压气机多圆弧叶型、多圆弧厚度分布叶型及修正NACA 6-系列叶型造型方法,利用NASA/GE研发的高效节能发动机E^(3)高压压气机前6.5级的准三维设计数据,建立了相应的压气机三维几何参数化模型,并通过数值模拟研究了造型压气机在多个转速下的内部流动及增压特性。结果表明:在100%设计转速近最高效率点处,造型压气机的质量流量与设计目标及与原型压气机试验结果的误差分别为0.4%和0.6%,绝热效率的误差分别为1.9%和0.9%;造型压气机模拟及原型压气机试验所得总压比均低于设计目标;前排跨声转子叶片气动负荷过大导致激波脱体严重是压气机的主要设计不足,改善叶片气动负荷分布、缓解叶片溢流是提升压气机气动性能尤其是设计点总压比的有效方法。研究成果对开展多级跨声速压气机设计优化具有参考价值。 展开更多
关键词 E^(3)高压压气机 几何参数化 增压特性 激波脱体
下载PDF
沿流向放样的轴流压气机叶道设计方法
4
作者 赵奕博 李浩 +1 位作者 琚亚平 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2023年第4期39-48,共10页
针对高负荷压气机角区流动质量差、流动损失显著的问题,以NASA Rotor37为研究对象,提出了一种轴流压气机叶道几何定义与参数化建模新方法,实现了叶片、角区和端壁结构的协同设计及叶片到端壁的高阶光滑过渡,在此基础上,开展了前缘边条... 针对高负荷压气机角区流动质量差、流动损失显著的问题,以NASA Rotor37为研究对象,提出了一种轴流压气机叶道几何定义与参数化建模新方法,实现了叶片、角区和端壁结构的协同设计及叶片到端壁的高阶光滑过渡,在此基础上,开展了前缘边条和叶身融合对压气机转子气动性能和内部流场影响的数值研究。结果表明,针对NASA Rotor37流道,采用该叶道设计方法对压气机转子角区改型后:绝热效率提高了0.3%,喘振裕度由11.85%增加至18.8%;下凹端壁可以帮助稳定压气机转子通流能力并减少改型过程的流量漂移;叶身融合结构在叶根处提供径向向上的压力梯度,将低能流体从叶片角区吹向主流区以减少角区流动损失,而近叶根处主流区流动损失略有增加。研究成果对发展轴流压气机先进设计方法应具有参考价值。 展开更多
关键词 轴流压气机 角区流动 非对称端壁 前缘边条 叶身融合
下载PDF
叶顶间隙几何不确定性对离心叶轮气动性能的影响 被引量:7
5
作者 傅珏 蒋伟 +2 位作者 柳一鸣 琚亚平 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2018年第11期51-57,共7页
为了揭示受叶顶间隙几何不确定性影响的离心叶轮气动性能和内部流场变化规律,采用非嵌入式多项式混沌(NIPC)法并结合计算流体力学方法,以Krain离心叶轮为例,研究了叶顶间隙服从高斯分布时的叶轮效率、压比性能统计变化规律,并从叶轮出... 为了揭示受叶顶间隙几何不确定性影响的离心叶轮气动性能和内部流场变化规律,采用非嵌入式多项式混沌(NIPC)法并结合计算流体力学方法,以Krain离心叶轮为例,研究了叶顶间隙服从高斯分布时的叶轮效率、压比性能统计变化规律,并从叶轮出口流动不均匀性、相对马赫数分布、载荷分布等角度探究其物理机制。结果表明:不同流量工况下,Krain叶轮效率和压比的不确定带大小相当;小流量工况下,叶轮出口流动对叶顶间隙几何不确定性的敏感度较强;大流量工况下,叶轮进口及吸力面附近流动对叶顶间隙几何不确定性的敏感程度较强。研究结果有利于深刻认识叶顶间隙几何不确定性对离心叶轮气动性能和内部流场的影响,同时为深入开展流体机械不确定性流动分析研究提供了一定的理论指导。 展开更多
关键词 离心叶轮 叶顶间隙 不确定性分析 多项式混沌法 气动性能
下载PDF
多级轴流压气机末级静子叶片串列改型方法及流动机理研究 被引量:3
6
作者 陈玲 琚亚平 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2022年第10期170-179,共10页
针对多级轴流压气机末级静子叶片转折角大、角区分离严重、流动损失显著的问题,以某工业用6.5级轴流压气机末级静子为主要研究对象,提出了一种三维串列叶片几何参数化建模方法,能够实现对串列参数沿叶高分布规律、叶片弯掠等三维特征的... 针对多级轴流压气机末级静子叶片转折角大、角区分离严重、流动损失显著的问题,以某工业用6.5级轴流压气机末级静子为主要研究对象,提出了一种三维串列叶片几何参数化建模方法,能够实现对串列参数沿叶高分布规律、叶片弯掠等三维特征的灵活控制,以支持单列叶片向串列叶片的改型设计,并结合数值计算方法,开展了末级串列静子改型设计对压气机气动性能及内部流场的影响研究。结果表明:对6.5级轴流压气机,末级静子叶片由单列改为串列后,6.5级压气机设计工况绝热效率提高了0.965%,稳定工况范围有所拓宽;串列叶片缝隙区射流的流动局部加速作用能够抑制角区分离的周向和径向发展;相比机匣附近,轮毂附近串列参数的变化对压气机气动性能的影响更为显著,且周向交错度较轴向重叠度的影响更大。研究成果对多级环境下压气机三维串列叶片设计具有参考价值。 展开更多
关键词 三维串列叶片 末级静子 轴向重叠度 周向交错度 多级轴流压气机
下载PDF
离心压缩机级通流与盘盖侧泄漏流的耦合分析 被引量:6
7
作者 王晟旻 琚亚平 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2017年第7期73-77,共5页
为揭示离心压缩机通流与轮盘及轮盖侧迷宫密封泄漏流动之间的耦合作用机理,以某天然气管线用离心压缩机首级为研究对象,数值分析了压缩机级在不考虑密封结构模型以及考虑密封结构模型下的流场分布。结果表明:在考虑密封结构的情况下,压... 为揭示离心压缩机通流与轮盘及轮盖侧迷宫密封泄漏流动之间的耦合作用机理,以某天然气管线用离心压缩机首级为研究对象,数值分析了压缩机级在不考虑密封结构模型以及考虑密封结构模型下的流场分布。结果表明:在考虑密封结构的情况下,压缩机整级的等熵效率和总压比都有所下降,且在小流量工况时下降更为明显;等厚度高低密封齿和梯形密封齿均能有效衰减泄漏气流压力能和动能,从而有效降低泄漏量;密封通道内的泄漏气流在叶轮进出口与通流区气流掺混时易造成通流分离;受轮盖密封泄漏气流的影响,在叶轮进口轮盖侧存在明显的熵增区域,使得压缩机整级性能下降。研究工作为深入了解离心压缩机密封流动结构及其对通流的影响、完善离心压缩机在实际运行状态下的气动性能预测模型奠定了坚实理论基础。 展开更多
关键词 离心压缩机 迷宫密封 级环境 通流 泄漏流动
下载PDF
叶栅叶型正反设计的伴随优化方法 被引量:3
8
作者 朱玉杰 琚亚平 张楚华 《西安交通大学学报》 EI CAS CSCD 北大核心 2019年第1期100-105,174,共7页
针对基于流场Navier-Stokes方程求解的叶栅叶型正反问题设计优化中,计算量随设计变量数目急剧增加的问题,采用伴随方法建立了集叶片几何参数化、网格生成、流场求解、伴随场求解与优化求解于一体的优化求解方法。针对反问题设计中目标... 针对基于流场Navier-Stokes方程求解的叶栅叶型正反问题设计优化中,计算量随设计变量数目急剧增加的问题,采用伴随方法建立了集叶片几何参数化、网格生成、流场求解、伴随场求解与优化求解于一体的优化求解方法。针对反问题设计中目标控制参数分布难以给定的问题,通过分析叶栅叶型正问题优化求解结果,给出了反问题优化求解所需的较优目标压力分布。利用自主程序完成了以减弱或消除流动分离、提高叶栅气动性能为目的的叶栅叶型正、反问题自动优化求解,优化后的叶型叶栅总压损失系数分别降低了5.69%、4.50%。研究表明,叶栅叶型吸力面曲率的减小、叶片前加载和中后部逆压梯度的减小,可有效抑制叶片尾缘附近的流动分离。该研究工作对发展高效宽工况叶栅叶型设计技术具有一定的参考价值。 展开更多
关键词 叶栅叶型 伴随方法 正问题 反问题 优化
下载PDF
Fluid-structure Interaction Analysis and Lifetime Estimation of a Natural Gas Pipeline Centrifugal Compressor under Near-choke and Near-surge Conditions 被引量:12
9
作者 ju yaping LIU Hui +2 位作者 YAO Ziyun XING Peng ZHANG Chuhua 《Chinese Journal of Mechanical Engineering》 SCIE EI CAS CSCD 2015年第6期1261-1268,共8页
Up to present, there have been no studies concerning the application of fluid-structure interaction(FSI) analysis to the lifetime estimation of multi-stage centrifugal compressors under dangerous unsteady aerodynami... Up to present, there have been no studies concerning the application of fluid-structure interaction(FSI) analysis to the lifetime estimation of multi-stage centrifugal compressors under dangerous unsteady aerodynamic excitations. In this paper, computational fluid dynamics(CFD) simulations of a three-stage natural gas pipeline centrifugal compressor are performed under near-choke and near-surge conditions, and the unsteady aerodynamic pressure acting on impeller blades are obtained. Then computational structural dynamics(CSD) analysis is conducted through a one-way coupling FSI model to predict alternating stresses in impeller blades. Finally, the compressor lifetime is estimated using the nominal stress approach. The FSI results show that the impellers of latter stages suffer larger fluctuation stresses but smaller mean stresses than those at preceding stages under near-choke and near-surge conditions. The most dangerous position in the compressor is found to be located near the leading edge of the last-stage impeller blade. Compressor lifetime estimation shows that the investigated compressor can run up to 102.7 h under the near-choke condition and 200.2 h under the near-surge condition. This study is expected to provide a scientific guidance for the operation safety of natural gas pipeline centrifugal compressors. 展开更多
关键词 three-stage pipeline centrifugal compressor natural gas aerodynamic performance structural durability lifetime estimation CHOKE SURGE
下载PDF
Robust Design Optimization Method for Centrifugal Impellers under Surface Roughness Uncertainties Due to Blade Fouling 被引量:9
10
作者 ju yaping ZHANG Chuhua 《Chinese Journal of Mechanical Engineering》 SCIE EI CAS CSCD 2016年第2期301-314,共14页
Blade fouling has been proved to be a great threat to compressor performance in operating stage. The current researches on fouling-induced performance degradations of centrifugal compressors are based mainly on simpli... Blade fouling has been proved to be a great threat to compressor performance in operating stage. The current researches on fouling-induced performance degradations of centrifugal compressors are based mainly on simplified roughness models without taking into account the realistic factors such as spatial non-uniformity and randomness of the fouling-induced surface roughness. Moreover, little attention has been paid to the robust design optimization of centrifugal compressor impellers with considerations of blade fouling. In this paper, a multi-objective robust design optimization method is developed for centrifugal impellers under surface roughness uncertainties due to blade fouling. A three-dimensional surface roughness map is proposed to describe the nonuniformity and randomness of realistic fouling accumulations on blades. To lower computational cost in robust design optimization, the support vector regression(SVR) metamodel is combined with the Monte Carlo simulation(MCS) method to conduct the uncertainty analysis of fouled impeller performance. The analyzed results show that the critical fouled region associated with impeller performance degradations lies at the leading edge of blade tip. The SVR metamodel has been proved to be an efficient and accurate means in the detection of impeller performance variations caused by roughness uncertainties. After design optimization, the robust optimal design is found to be more efficient and less sensitive to fouling uncertainties while maintaining good impeller performance in the clean condition. This research proposes a systematic design optimization method for centrifugal compressors with considerations of blade fouling, providing a practical guidance to the design of advanced centrifugal compressors. 展开更多
关键词 centrifugal impeller robust design optimization surface roughness uncertainty analysis
下载PDF
三级高速轴流压气机变转速导静联调特性及流动机理研究
11
作者 花雨 李聪 +1 位作者 琚亚平 张楚华 《工程热物理学报》 EI CAS CSCD 北大核心 2024年第2期418-426,共9页
变转速及变导静叶片角度是实现轴流压气机扩稳增效的两个最有效性能调节方法,对航机、燃机、工业压缩系统的高效宽工况运行具有重要影响。本文针对多级高速轴流压气机易出现低转速下裕度不够的问题,对某3.5级高速轴流压气机在100%~70%... 变转速及变导静叶片角度是实现轴流压气机扩稳增效的两个最有效性能调节方法,对航机、燃机、工业压缩系统的高效宽工况运行具有重要影响。本文针对多级高速轴流压气机易出现低转速下裕度不够的问题,对某3.5级高速轴流压气机在100%~70%设计转速范围内四种转速下的内部流动进行了全三维定常数值模拟,重点分析了低转速下的失速特性及其流动失稳机理。结果表明:100%、90%设计转速下失速裕度为16.58%、32.22%;在80%、70%设计转速下失速裕度分别为0.22%、0.34%,性能曲线表现出突变失速的不连续性。分析流动机理后发现80%转速下流动失稳现象为R1叶尖部分叶高失速,该失速无法通过增大流量解除;70%转速下流动失稳现象为S2角区分离引起S3双开型角区分离失速,该失速可通过增大流量解除。研究工作对完善多级高速轴流压气机导静叶调节扩稳设计具有指导价值。 展开更多
关键词 多级高速轴流压气机 变转速导静叶联调 突变失速 数值模拟
原文传递
三级轴流压气机复合沟槽式机匣处理扩稳研究
12
作者 张夏雯 花雨 +2 位作者 李聪 琚亚平 张楚华 《工程热物理学报》 EI CSCD 北大核心 2023年第10期2673-2684,共12页
机匣处理作为现代高负荷压气机解决裕度不足问题的重要手段,对保障燃机运行安全性及装置技术先进性具有重要影响。本文针对多级高速轴流压气机易出现低转速下裕度不足的问题,对某3.5级高速轴流压气机在新型复合沟槽式机匣处理下70%~100... 机匣处理作为现代高负荷压气机解决裕度不足问题的重要手段,对保障燃机运行安全性及装置技术先进性具有重要影响。本文针对多级高速轴流压气机易出现低转速下裕度不足的问题,对某3.5级高速轴流压气机在新型复合沟槽式机匣处理下70%~100%设计转速范围内的四种转速进行了全三维定常数值模拟,结合低转速下光壁机匣的失稳特性重点分析了复合沟槽式机匣处理的扩稳流动机理。结果表明:复合沟槽式机匣处理显著提高了80%设计转速下的失速裕度,裕度提高了13.73%。扩稳机理为:后段槽内泄漏间歇性降低叶尖两侧压差,削弱叶顶间隙泄漏射流强度;叶顶间隙泄漏流、叶尖前缘泄漏涡以及叶片表面附面层被抽吸进入后段槽阻断旋绕从而抑制泄漏涡的形成;前段槽回流在一定范围内降低了第一级动叶进气冲角。研究工作对提高压气机性能调节技术水平具有重要指导价值。 展开更多
关键词 多级高速轴流压气机 机匣处理 扩稳调节 数值模拟
原文传递
整机环境下离心压气机叶轮流固耦合不确定性量化分析
13
作者 徐永祥 宋嘉涛 +3 位作者 张本瑞 李震 琚亚平 张楚华 《西安交通大学学报》 EI CAS 2024年第10期121-132,共12页
针对流固耦合不确定性量化计算中精度和成本难以有效平衡的问题,从某涡轮增压器离心压气机双向流固耦合计算出发,提出了考虑叶片热态变形和蜗壳非对称影响的离心压气机内流场计算简化模型,并结合集成代理模型和蒙特卡罗方法,开展了整机... 针对流固耦合不确定性量化计算中精度和成本难以有效平衡的问题,从某涡轮增压器离心压气机双向流固耦合计算出发,提出了考虑叶片热态变形和蜗壳非对称影响的离心压气机内流场计算简化模型,并结合集成代理模型和蒙特卡罗方法,开展了整机环境下离心压气机叶轮流固耦合的不确定性量化分析。数值结果表明:在保证计算精度的前提下,相较于整机流固耦合模型,提出的简化模型将计算时间缩短了83%~98%;材料属性不确定性对叶轮效率的影响大于转速不确定性,密度、泊松比和弹性模量3种不确定性因素对叶轮效率的影响程度相当,转速对叶片热态变形不确定性的影响程度在转速不确定性中最为显著;叶片热态变形不确定性对叶轮效率的影响主要体现在效率均值的变化,而与热态变形无关的转速自身、叶轮出口非对称流场等不确定性因素对叶轮效率的影响主要体现在分散度的变化。研究工作突显了流固耦合分析在离心压气机不确定性量化分析中的必要性,为今后开展鲁棒优化设计提供了理论基础。 展开更多
关键词 离心压气机 叶轮 流固耦合 不确定性量化 效率
下载PDF
多块多重网格法及其跨声速转子内流并行模拟 被引量:2
14
作者 刘安 琚亚平 张楚华 《航空动力学报》 EI CAS CSCD 北大核心 2018年第7期1705-1712,共8页
为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表... 为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。 展开更多
关键词 多块结构化网格 多重网格 对接分区 内界面处理方法 并行计算流体动力学(CFD) 跨声速转子
原文传递
A Rapid and Automatic Optimal Design Method for Six-Stage Axial-Flow Industry Compressor 被引量:4
15
作者 LEI Fan ju yaping ZHANG Chuhua 《Journal of Thermal Science》 SCIE EI CAS CSCD 2021年第5期1658-1673,共16页
Existing aerodynamic design systems for multi-stage axial-flow compressor suffer from several limitations,such as experience dependent models and time costly simulations.Few attempts,however,have been devoted to the r... Existing aerodynamic design systems for multi-stage axial-flow compressor suffer from several limitations,such as experience dependent models and time costly simulations.Few attempts,however,have been devoted to the rapid and automatic optimization of aerodynamic performance at the preliminary design phase,which plays a crucial role in the final aerodynamic performance.In this work,a rapid and automatic aerodynamic optimal design method is developed for the multi-stage axial-flow compressor based on one-dimensional meanline design method,radial-equilibrium equation and genetic algorithm.The one-dimensional performance prediction model includes some popular empirical correlations to calculate the flow loss,incidence angle,deviation angle and flow blockage.The radial-equilibrium equation is solved to obtain the spanwise distribution of aerodynamic and thermodynamic parameters at the inlet and outlet of each blade row.The genetic algorithm is used for an automatic search of the global optimal compressor configuration aiming at maximizing the design efficiency.The developed method is illustrated with the aerodynamic optimal design of a 6-stage axial-flow industry compressor and verified by computational fluid dynamics simulations.The results show that the developed method is capable of improving effectively the design efficiency and predicting accurately the aerodynamic performance of the 6-stage axial-flow industry compressor in a few minutes.This work is of scientific significance to improve the axial-flow compressor design system and of engineering importance to release the designers from the heavy experience dependence especially at the preliminary design phase. 展开更多
关键词 one-dimensional meanline design two-dimensional throughflow design three-dimensional CFD simulation optimal design genetic algorithm multi-stage axial-flow compressor
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部