期刊文献+
共找到16篇文章
< 1 >
每页显示 20 50 100
基于嵌入式多项式混沌展开法的随机边界下流动与传热问题不确定性量化 被引量:2
1
作者 姜昌伟 刘星 +2 位作者 石尔 李韬海 江怡 《应用数学和力学》 CSCD 北大核心 2021年第5期481-491,共11页
提出了一种嵌入式多项式混沌展开(polynomial chaos expansion,PCE)的随机边界条件下流动与传热问题不确定性量化方法及有限元程序框架.该方法利用Karhunen-Loeve展开表达随机输入边界条件,以及嵌入式多项式混沌展开法表达输出随机场;... 提出了一种嵌入式多项式混沌展开(polynomial chaos expansion,PCE)的随机边界条件下流动与传热问题不确定性量化方法及有限元程序框架.该方法利用Karhunen-Loeve展开表达随机输入边界条件,以及嵌入式多项式混沌展开法表达输出随机场;同时利用谱分解技术将控制方程转化为一组确定性控制方程,并对每个多项式混沌进行求解得到其统计特征.与Monte-Carlo法相比,该方法能够准确高效地预测随机边界条件下流动与传热问题的不确定性特征,同时可以节省大量计算资源. 展开更多
关键词 不确定性量化 嵌入式多项式混沌展开 随机边界条件 流动与传热 随机有限元
下载PDF
基于多项式混沌展开的交直流系统全纯嵌入概率潮流计算方法
2
作者 李雪 付云跃 +1 位作者 姜涛 李国庆 《电力系统自动化》 EI CSCD 北大核心 2024年第18期177-188,共12页
为快速、准确量化分析风电出力不确定性对交直流电力系统潮流分布的影响,提出一种基于多项式混沌展开(PCE)的交直流电力系统全纯嵌入概率潮流计算方法。该方法首先根据风电出力的概率分布特征选择最优正交基函数,构造近似风电出力概率... 为快速、准确量化分析风电出力不确定性对交直流电力系统潮流分布的影响,提出一种基于多项式混沌展开(PCE)的交直流电力系统全纯嵌入概率潮流计算方法。该方法首先根据风电出力的概率分布特征选择最优正交基函数,构造近似风电出力概率分布特征的PCE表达式;其次,将该PCE表达式引入交直流电力系统的全纯嵌入潮流方程中,构建基于PCE的交直流电力系统全纯嵌入概率潮流计算模型;再次,通过Galerkin投影将所构建的全纯嵌入概率潮流计算模型转化为高维确定性全纯嵌入潮流计算模型;然后,借助确定性全纯嵌入潮流模型求解方法,实现对所转化的高维确定性全纯嵌入潮流模型的求解,并根据所得PCE逼近系数计算交直流电力系统潮流的概率分布特征;最后,通过修改的PJM 5节点、IEEE 30节点和IEEE 118节点交直流测试系统算例验证所提方法的准确性和有效性。 展开更多
关键词 概率潮流 多项式混沌展开 全纯嵌入 Galerkin投影 交直流电力系统 不确定性
下载PDF
非嵌入式多项式混沌法在爆轰产物JWL参数评估中的应用 被引量:14
3
作者 王瑞利 刘全 温万治 《爆炸与冲击》 EI CAS CSCD 北大核心 2015年第1期9-15,共7页
介绍了非嵌入多项式混沌法的数学模型,给出了非嵌入式多项式混沌法进行不确定度量化的主要步骤。采用此方法研究了平面、散心爆轰问题数值模拟中,JWL模型参数R1、R2服从均匀分布的随机变量时所引起的爆轰过程计算结果的不确度性,着重分... 介绍了非嵌入多项式混沌法的数学模型,给出了非嵌入式多项式混沌法进行不确定度量化的主要步骤。采用此方法研究了平面、散心爆轰问题数值模拟中,JWL模型参数R1、R2服从均匀分布的随机变量时所引起的爆轰过程计算结果的不确度性,着重分析了爆轰传播过程中压力与密度的统计特性。研究结果表明,非嵌入式多项式混沌法可以为模型输入参数不确定性的传播对输出结果响应量的影响建立一种有效不确定度评估方法,为使用JWL模型时选取参数提供参考。 展开更多
关键词 爆炸力学 JWL状态方程 参数选取 嵌入多项式混沌 不确定度量化
下载PDF
参数随机Lamb波频散特性的非嵌入式多项式混沌方法 被引量:1
4
作者 于保华 胡小平 杨世锡 《振动工程学报》 EI CSCD 北大核心 2019年第4期712-717,共6页
Lamb波频散特性理论分析是板状结构开展导波无损检测的基础,传统的参数确定Lamb波频散特性分析方法往往忽略名义参数与实际参数之间存在的误差,而这种忽略已逐渐不能适应Lamb波高精度高效率无损检测的需求.提出一种融合谱配置方法的非... Lamb波频散特性理论分析是板状结构开展导波无损检测的基础,传统的参数确定Lamb波频散特性分析方法往往忽略名义参数与实际参数之间存在的误差,而这种忽略已逐渐不能适应Lamb波高精度高效率无损检测的需求.提出一种融合谱配置方法的非嵌入式多项式混沌分析方法,利用Galerkin映射将频散特性进行正交多项式混沌展开,得到参数随机Lamb波频散的统计特性,并与蒙特卡洛随机模拟方法进行对比,两种方法获取结果一致,而此方法具有明显的分析效率优势&研究有助于完善板状结构Lamb波传播与无损检测分析理论. 展开更多
关键词 健康监测#无损检测# Lamb波#谱配置方法#非嵌入多项式混沌
下载PDF
不确定水声场随机多项式系数解法 被引量:6
5
作者 程广利 张明敏 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2013年第1期21-25,共5页
基于非嵌入式随机多项式展开法求解了含不确定海洋环境参数的波动方程,推导了数值积分法求解多项式系数的过程.针对常规概率配点法不能准确计算随不确定输入量剧烈变化的传播损失,提出了分段概率配点法,将输入变量区间进行合理分段,基... 基于非嵌入式随机多项式展开法求解了含不确定海洋环境参数的波动方程,推导了数值积分法求解多项式系数的过程.针对常规概率配点法不能准确计算随不确定输入量剧烈变化的传播损失,提出了分段概率配点法,将输入变量区间进行合理分段,基于非嵌入式随机多项式展开法获得每段的随机多项式,继而得到整个输入变量区间对应的传播损失表达式.结果表明,数值积分法仅适合于计算随单个不确定海洋环境参数不剧烈变化的传播损失,分段概率配点法可高精度计算随不确定输入量剧烈变化的传播损失. 展开更多
关键词 不确定水声场 嵌入随机多项式展开 随机多项式系数 数值积分法 分段概率配点法
下载PDF
基于混沌多项式的RBCC飞行器上升段鲁棒轨迹快速优化
6
作者 闫循良 王培臣 +2 位作者 王舒眉 杨宇轩 王宽 《航空学报》 EI CAS CSCD 北大核心 2023年第21期299-316,共18页
针对火箭基组合循环(RBCC)高超声速飞行器上升段轨迹设计所具有的动力系统工作模态复杂、推力与飞行状态存在强耦合、模型强非线性、存在多种复杂约束限制和参数不确定性因素影响等典型特征,提出了一种基于非嵌入式混沌多项式、高斯求... 针对火箭基组合循环(RBCC)高超声速飞行器上升段轨迹设计所具有的动力系统工作模态复杂、推力与飞行状态存在强耦合、模型强非线性、存在多种复杂约束限制和参数不确定性因素影响等典型特征,提出了一种基于非嵌入式混沌多项式、高斯求积法和序列凸优化的RBCC动力上升段鲁棒轨迹优化方法,以提高轨迹的抗干扰能力和过程可靠性。首先,构建了考虑参数不确定性的RBCC高超声速上升段鲁棒轨迹优化模型,并设计了基于高斯求积与非嵌入式混沌多项式的不确定性量化传播算法,从而将其转化为维数扩展的确定性轨迹优化问题;随后,基于凸优化理论对该问题进行凸化和离散,设计了一种基于序列凸优化算法的轨迹优化求解策略,以实现对该高维确定性问题的快速求解。某空基投放上升轨迹优化结果表明,基于所构建模型和轨迹优化方法可以有效地完成RBCC高超声速飞行器上升段鲁棒轨迹优化,优化结果符合RBCC动力系统工作特点;与传统确定性轨迹优化算法相比,所提方法能够有效降低随机干扰对上升段轨迹的影响,从而提升轨迹的可靠性与鲁棒性。 展开更多
关键词 火箭基组合循环(RBCC) 高超声速飞行器 上升段轨迹 鲁棒优化 不确定性 嵌入混沌多项式 高斯求积策略 凸优化
原文传递
离心泵进口来流速度扰动不确定性对水力性能及流场的影响
7
作者 陈学炳 杨博文 +2 位作者 张人会 郭广强 杨军虎 《农业工程学报》 EI CAS CSCD 北大核心 2024年第12期77-85,共9页
离心泵运行过程中其进口来流速度随时间的波动对泵的水力性能具有不可忽略的影响,该研究采用非嵌入式多项式混沌(non-intrusive polynomial chaos, NIPC)方法对进口来流速度扰动的不确定性进行量化分析,探究其对离心泵水力性能的影响。... 离心泵运行过程中其进口来流速度随时间的波动对泵的水力性能具有不可忽略的影响,该研究采用非嵌入式多项式混沌(non-intrusive polynomial chaos, NIPC)方法对进口来流速度扰动的不确定性进行量化分析,探究其对离心泵水力性能的影响。结果表明:随机来流速度对泵的水力性能具有较大影响,且不确定度越大,其影响越大;进口来流速度扰动会引起泵叶片压力面尾缘与吸力面上压力分布及叶轮流道内流场的变化,从而造成泵的扬程及效率的波动且波动范围较大;同时,来流不确定性的影响在叶轮内流场中的传播是非对称且非均匀的。不同工况下来流速度随时间的扰动对泵性能的影响有所差异,大部分工况下来流速度扰动会造成泵性能的下降,其中不确定度为5%时,扬程最大可下降0.4 m,效率下降3%。对不同工况下泵的进口流速不确定性进行量化分析,能够对离心泵在整个运行工况下的稳健性进行综合评定,为泵稳健性设计提供一定的基础。 展开更多
关键词 离心泵 流速 水力性能 不确定性分析 嵌入多项式混沌
下载PDF
叶片安装角偏差对动叶性能影响的不确定性研究 被引量:1
8
作者 李玉 楚武利 姬田园 《西安交通大学学报》 EI CAS CSCD 北大核心 2023年第4期49-59,共11页
为了揭示叶片安装角偏差对压气机性能的影响,基于非嵌入式混沌多项式法结合数值模拟,以NASA Rotor37为研究对象,进行了叶片安装角偏差的不确定性量化分析。在峰值效率工况和近失速工况下,进行气动性能统计分析,获得了不同输出参数的均... 为了揭示叶片安装角偏差对压气机性能的影响,基于非嵌入式混沌多项式法结合数值模拟,以NASA Rotor37为研究对象,进行了叶片安装角偏差的不确定性量化分析。在峰值效率工况和近失速工况下,进行气动性能统计分析,获得了不同输出参数的均值、标准差和概率分布,探讨各性能参数与安装角偏差间的相关程度并进行流场分析,深入研究安装角偏差波动对内部流场规律的影响。结果表明:安装角偏差对气动性能的影响程度与所处工况有关,近失速工况下,绝热效率受安装角偏差影响的敏感性更大,绝热效率的相对波动幅度相比较峰值效率工况增加了22.95%;总压比与安装角偏差呈较强线性负相关,绝热效率与安装角偏差呈较强线性正相关;不同工况下流场敏感区域有所不同,峰值效率工况下,激波及波后叶尖间隙泄漏流动轨迹受误差不确定性的影响程度较大;近失速工况下,波后间隙泄漏涡破碎区域和吸力面分离流动对安装角不确定性较为敏感;流场中不同区域的损失受安装角偏差影响程度不同,其中尾缘处的损失对安装角偏差的不确定性最为敏感。 展开更多
关键词 压气机转子 安装角偏差 不确定性量化 嵌入混沌多项式
下载PDF
真实运行状态下叶顶间隙尺寸波动对跨声速转子气动性能的影响研究 被引量:7
9
作者 马驰 高丽敏 +1 位作者 李瑞宇 李杰 《推进技术》 EI CAS CSCD 北大核心 2020年第9期1958-1966,共9页
为分析跨声速转子实时波动的叶顶间隙尺寸对气动性能的影响,对跨声速压气机转子真实运行状态下一个稳定工况实时波动的叶顶间隙数据进行统计分析,获得了叶顶间隙尺寸的总体水平、波动幅值和概率分布形式。以跨声速压气机转子NASA Rotor... 为分析跨声速转子实时波动的叶顶间隙尺寸对气动性能的影响,对跨声速压气机转子真实运行状态下一个稳定工况实时波动的叶顶间隙数据进行统计分析,获得了叶顶间隙尺寸的总体水平、波动幅值和概率分布形式。以跨声速压气机转子NASA Rotor 37为研究对象,采用非嵌入式混沌多项式不确定性量化方法,对100%转速下近失速和峰值效率两个工况施加相同叶顶间隙波动对跨声速转子气动性能的影响进行了不确定性量化分析。结果表明,真实运行状态下叶顶间隙波动对气动性能的总体水平无影响,但会缩小喘振裕度3.75%;近失速工况对叶顶间隙波动更为敏感,各参数的相对波动幅值均较峰值效率工况有所增大,等熵效率受叶顶间隙波动的影响比质量流量和总压比大;近失速工况下叶顶间隙波动在叶高方向上的影响范围和强度均大于峰值效率工况,98%叶高位置处静压系数和总压损失系数最大相对波动幅值分别可达14.84%和5%。峰值效率工况下流场中的不确定性主要由叶顶泄漏流及其与激波相互作用引起;而近失速工况下流场当中的不确定性则是由激波和吸力面分离流动起主要作用。 展开更多
关键词 不确定性量化 跨声速压气机转子 叶顶间隙 气动性能 嵌入混沌多项式
下载PDF
浅海声场对不确定环境参数的灵敏度分析方法 被引量:3
10
作者 程广利 胡金华 张明敏 《兵工学报》 EI CAS CSCD 北大核心 2014年第4期567-571,共5页
为分析不确定环境参数对浅海声场的影响,用非嵌入式随机多项式展开法获得声场与不确定环境参数之间函数关系式,建立声场对单个不确定参数的局部灵敏度分析方法;当多参数不确定时,提出依据某参数变化前后声场概率密度的变化,定性表征声... 为分析不确定环境参数对浅海声场的影响,用非嵌入式随机多项式展开法获得声场与不确定环境参数之间函数关系式,建立声场对单个不确定参数的局部灵敏度分析方法;当多参数不确定时,提出依据某参数变化前后声场概率密度的变化,定性表征声场对该不确定参数的全局灵敏度分析方法,引入Sobol指数法定量计算声场对不确定环境参数的全局灵敏度。算例结果表明文中提出的分析方法可快速获得浅海声场对不确定参数的灵敏度,可用于优化不确定浅海环境参数模型。 展开更多
关键词 声学 不确定水声场 局部灵敏度分析 全局灵敏度分析 嵌入随机多项式展开 Sobol指数法
下载PDF
叶顶间隙偏差对叶片气动性能影响的不确定性研究 被引量:4
11
作者 姬田园 楚武利 +1 位作者 戴雨晨 张耀峰 《推进技术》 EI CAS CSCD 北大核心 2022年第10期129-141,共13页
在压气机叶片加工过程中,实际加工得到的叶片外形与设计叶型不可避免会存在一些偏差。为研究叶片叶顶间隙尺寸波动对性能的影响,以Rotor 37为研究对象,采用三维定常数值模拟方法,基于非嵌入式混沌多项式中二维正态分布变量产生方法,分... 在压气机叶片加工过程中,实际加工得到的叶片外形与设计叶型不可避免会存在一些偏差。为研究叶片叶顶间隙尺寸波动对性能的影响,以Rotor 37为研究对象,采用三维定常数值模拟方法,基于非嵌入式混沌多项式中二维正态分布变量产生方法,分别对叶片前缘叶顶间隙和尾缘叶顶间隙进行尺寸波动干扰,研究了叶顶间隙不确定性对气动性能和流场的不确定性影响,评估了公差带内叶顶间隙偏差与气动性能变化的相关性,并探究了叶顶间隙变化对叶片稳定裕度的影响机理。研究发现,叶顶间隙偏差对等背压条件下叶片质量流量、等熵效率和总压比的平均水平几乎不会造成影响;叶顶间隙偏差所造成的气动性能分布标准差随工况点向不稳定边界的逼近而逐渐增大。同时,叶顶间隙不确定性会导致样本稳定裕度均值有所降低,不确定性分析中稳定裕度均值相对原型叶片下降了2.42%。流场结构方面,叶顶间隙偏差主要影响80%叶高以上部分流场;叶顶间隙偏差对叶顶区域泄漏流产生的影响,导致了叶片稳定裕度的变化。 展开更多
关键词 叶顶间隙 加工偏差 嵌入混沌多项式 气动性能 不确定性量化
下载PDF
声场计算模型的环境宽容性分析
12
作者 秦凯丽 胡长青 《声学技术》 CSCD 北大核心 2019年第1期15-23,共9页
浅海水声环境容易受到各种因素的影响,导致其环境参数具有很强的不确定性,依据环境宽容性选择快速、精确的声场计算模型是保证后续研究分析正确性的重要前提。以海水声速为例,简述了基于传播损失和声场互相关系数的计算模型环境宽容性... 浅海水声环境容易受到各种因素的影响,导致其环境参数具有很强的不确定性,依据环境宽容性选择快速、精确的声场计算模型是保证后续研究分析正确性的重要前提。以海水声速为例,简述了基于传播损失和声场互相关系数的计算模型环境宽容性分析方法。为定量描述环境和声场计算模型的失配情况,在研究浅海不确定环境对于声场空间相关性影响的基础上,结合水声环境不确定性推理模型,得到声场空间相关半径和传播损失概率分布可信区间,提出利用声场空间相关半径相对值来度量声场计算模型的环境宽容性,同时利用非嵌入式随机多项式展开(NON-Polynomial Chaos Expansion, NPCE)法,结合差值评定方法对得到的环境宽容区间进行验证,结果表明,利用声场相关半径相对值可以定量分析不确定性环境下声场计算模型的宽容性。 展开更多
关键词 宽容性分析 模型失配 声场空间相关性 不确定性海洋环境 嵌入随机多项式展开
下载PDF
来流马赫数波动对扩压叶栅气动性能的影响 被引量:5
13
作者 马驰 高丽敏 +1 位作者 李瑞宇 吴瑜 《工程热物理学报》 EI CAS CSCD 北大核心 2020年第4期851-859,共9页
准确地掌握来流马赫数的波动(即不确定性)对扩压叶栅的气动性能的影响对于指导压气机设计具有重要意义。采用非嵌入式混沌多项式方法,从统计学的角度评估了随机来流马赫数对扩压叶栅气动性能的影响,并着重对比分析了不同攻角和不同随机... 准确地掌握来流马赫数的波动(即不确定性)对扩压叶栅的气动性能的影响对于指导压气机设计具有重要意义。采用非嵌入式混沌多项式方法,从统计学的角度评估了随机来流马赫数对扩压叶栅气动性能的影响,并着重对比分析了不同攻角和不同随机来流马赫数工况下来流马赫数的不确定性对扩压叶栅气动性能的整体影响和随之带来的气动性能的波动影响。研究结果表明,来流马赫数的不确定性的确对叶栅气动性能产生影响:来流攻角偏离设计攻角越多,叶栅的气动性能对来流马赫数的不确定性越敏感,且在叶栅通道中反应较敏感区域的范围及位置均随着攻角的变化而变化;随机来流马赫数越高,叶栅的气动性能对来流马赫数的不确定性越敏感。 展开更多
关键词 嵌入混沌多项式 不确定性分析 来流马赫数 气动性能 扩压叶栅
原文传递
锥导乘波体构型的气动特性不确定度分析 被引量:7
14
作者 宋赋强 阎超 +1 位作者 马宝峰 鞠胜军 《航空学报》 EI CAS CSCD 北大核心 2018年第2期97-106,共10页
为研究锥导乘波体偏离设计条件下气动特性变化情况,采用稀疏的非嵌入式混沌多项式方法,对乘波体气动特性进行了不确定性量化及全局非线性灵敏度分析。首先,采用CATIA二次开发技术对锥导乘波体进行参数化建模;其次,在来流速度、温度、密... 为研究锥导乘波体偏离设计条件下气动特性变化情况,采用稀疏的非嵌入式混沌多项式方法,对乘波体气动特性进行了不确定性量化及全局非线性灵敏度分析。首先,采用CATIA二次开发技术对锥导乘波体进行参数化建模;其次,在来流速度、温度、密度和迎角满足特定扰动的条件下,通过拉丁超立方试验设计生成样本,并采用CFD进行计算;最后,根据试验设计样本建立响应面,通过混沌多项式分析得到了乘波体气动力系数的不确定度。灵敏度分析结果表明,迎角在锥导乘波体的气动特性变化中起主导作用。对马赫数和压强的流场不确定性分析结果表明,气动特性变化主因是乘波体前缘处的压力泄漏,影响了上表面压力分布,导致了气动性能的改变。 展开更多
关键词 高超声速 气动特性 乘波体 不确定度分析 嵌入多项式混沌
原文传递
端壁倒圆半径误差对高负荷静叶栅性能影响的不确定性分析 被引量:3
15
作者 郭正涛 楚武利 +2 位作者 姬田园 杨吉博 张耀峰 《航空动力学报》 EI CAS CSCD 北大核心 2021年第10期2173-2185,共13页
为对压气机静子叶片高精度设计与制造提供有力参考,以某高负荷压气机静子叶栅为研究对象,采用基于高斯分布型随机输入的非嵌入式多项式混沌方法,量化评估了端壁倒圆半径误差对最小损失和近失速两个工况下气动性能的不确定性影响。结果表... 为对压气机静子叶片高精度设计与制造提供有力参考,以某高负荷压气机静子叶栅为研究对象,采用基于高斯分布型随机输入的非嵌入式多项式混沌方法,量化评估了端壁倒圆半径误差对最小损失和近失速两个工况下气动性能的不确定性影响。结果表明:倒圆半径误差的不确定性对气动性能的平均水平影响不大,主要反映在气动性能参数的标准差上。提高倒圆半径的加工精度可提升气动性能的鲁棒性近一倍。加工精度一定时,为对叶片进行鲁棒性优化设计,应重点关注近失速工况下气动损失的鲁棒性。根据公差带内端壁倒圆半径与近失速工况下气动性能的关联性分析,倒圆半径与气动性能呈线性关系,应避免倒圆半径与设计值相比偏小的叶片投入使用,以期获得良好的气动性能。通过近失速流场的不确定性分析,分离流动对于倒圆半径误差更为敏感,是引起气动性能不确定性变化的主要因素。 展开更多
关键词 不确定性量化 端壁倒圆 加工误差 高负荷叶栅 嵌入混沌多项式
原文传递
垂直起降运载火箭返回轨迹不确定性优化 被引量:2
16
作者 赵剑 黄悦琛 +1 位作者 李海阳 何湘粤 《航空学报》 EI CAS CSCD 北大核心 2021年第11期255-269,共15页
针对垂直起降运载火箭一子级在返回着陆的过程中存在的参数不确定性,提出了一种基于非侵入式多项式混沌展开的序列优化和可靠度评估的返回轨迹不确定性优化方法。首先,建立了返回多飞行段轨迹在确定性条件下的优化模型。然后,为同时兼... 针对垂直起降运载火箭一子级在返回着陆的过程中存在的参数不确定性,提出了一种基于非侵入式多项式混沌展开的序列优化和可靠度评估的返回轨迹不确定性优化方法。首先,建立了返回多飞行段轨迹在确定性条件下的优化模型。然后,为同时兼顾轨迹的鲁棒性和可靠性,建立了由鲁棒最优目标函数、基于可靠度的路径约束和鲁棒等式约束组成的不确定性返回轨迹优化模型。最后,基于非侵入式多项式混沌展开方法对鲁棒目标函数和等式约束进行量化处理,将原随机鲁棒优化问题转化为高维状态空间中的等价确定性优化问题;为提高路径约束的可靠度评估效率,基于非侵入式多项式混沌展开方法对最可能点法进行改进,进一步发展了序列优化和可靠度评估策略。数值仿真结果表明,所提出的不确定性优化方法具有较好的鲁棒性,可以满足工程可靠性指标要求,同时还具有较高的精度和计算效率。 展开更多
关键词 垂直起降运载火箭 返回轨迹 参数不确定性 非侵入多项式混沌展开 随机鲁棒优化 序列优化和可靠度评估
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部