期刊文献+
共找到41篇文章
< 1 2 3 >
每页显示 20 50 100
锥角对旋转整流罩积冰影响的模拟实验 被引量:7
1
作者 胡娅萍 吉洪湖 +4 位作者 王健 陈宁立 曹广州 童辉 庞黎刚 《航空动力学报》 EI CAS CSCD 北大核心 2014年第3期495-503,共9页
采用缩比模型对4种不同锥角的旋转整流罩进行了冰风洞积冰模拟实验.推导出旋转表面积冰相似参数,并根据重要参数的匹配确定缩比整流罩模型的积冰实验参数.在实验中采用高速摄像系统记录冰生长过程及最终冰形.结果表明:4种锥角的整流罩... 采用缩比模型对4种不同锥角的旋转整流罩进行了冰风洞积冰模拟实验.推导出旋转表面积冰相似参数,并根据重要参数的匹配确定缩比整流罩模型的积冰实验参数.在实验中采用高速摄像系统记录冰生长过程及最终冰形.结果表明:4种锥角的整流罩表面积冰均由初期生成、分布连续的明冰和后期增长迅速的白色霜冰构成;锥角小于等于74°的整流罩表面霜冰为针状或粒状,积冰厚度较小;锥角大于80°的整流罩表面霜冰为羽毛状,积冰厚度较大,并伴随冰脱落现象,锥角较大的整流罩冰脱落位置向下游移动. 展开更多
关键词 旋转整流罩 锥角 缩比模型 积冰模拟 积冰脱落
原文传递
缩比模型飞行试验相似准则研究 被引量:4
2
作者 牛文韬 高永 《兵工自动化》 2021年第8期30-34,共5页
针对全尺寸飞机直接进行飞行试验存在风险大、成本高的问题,利用缩比模型进行相关试验。在风洞试验相似准则的基础上,从飞机的六自由度运动方程入手,利用相似定理研究缩比模型飞行试验遵循的运动相似准则和飞行控制系统相似准则。以某... 针对全尺寸飞机直接进行飞行试验存在风险大、成本高的问题,利用缩比模型进行相关试验。在风洞试验相似准则的基础上,从飞机的六自由度运动方程入手,利用相似定理研究缩比模型飞行试验遵循的运动相似准则和飞行控制系统相似准则。以某型固定翼飞机为例,对缩比模型与全尺寸飞机进行开环和闭环六自由度仿真,并验证所推导的相似准则。结果表明,试验结果能准确反映原机的稳定和控制特性。 展开更多
关键词 飞行试验 缩比模型 相似准则 六自由度仿真
下载PDF
Re=3900圆柱绕流的三维大涡模拟 被引量:20
3
作者 战庆亮 周志勇 葛耀君 《哈尔滨工业大学学报》 EI CAS CSCD 北大核心 2015年第12期75-79,共5页
为研究亚临界雷诺数范围内圆柱绕流流场特性及三维大涡模拟方法的适用性,基于C++语言及有限体积法开发了三维非结构化网格的大涡模拟计算程序.采用新的高稳定性高精度二阶离散格式,及Smagorinsky亚格子模型对Re=3 900均匀来流条件下的... 为研究亚临界雷诺数范围内圆柱绕流流场特性及三维大涡模拟方法的适用性,基于C++语言及有限体积法开发了三维非结构化网格的大涡模拟计算程序.采用新的高稳定性高精度二阶离散格式,及Smagorinsky亚格子模型对Re=3 900均匀来流条件下的圆柱绕流问题进行数值模拟,并统计获得了平均流场参数及湍流流场的详细结构特性.结果表明:采用本文的网格、计算步长和高稳定性二阶离散精度大涡模拟方法计算所得的湍流场一阶统计特性和二阶统计特性与实验值吻合很好.验证了大涡模拟程序在模拟亚临界雷诺数下圆柱绕流流场平均值及脉动值的合理性. 展开更多
关键词 Smagorinsky亚格子模型 圆柱绕流 亚临界雷诺数 非结构化网格 二阶延迟修正
下载PDF
大涡模拟在水下航行体周围黏性流场计算中的初步应用 被引量:10
4
作者 黄振宇 缪国平 《水动力学研究与进展(A辑)》 CSCD 北大核心 2006年第2期190-197,共8页
在空间过滤得到的不可压缩流动的大涡模拟控制方程基础上,利用基于三阶迎风偏置格式的混合差分及二阶中心差分离散控制方程的对流项及耗散项,离散得到的线性方程组通过无矩阵计算的隐式LU-SGS(Lower-Upper SymmetricGauss Seidel)分步... 在空间过滤得到的不可压缩流动的大涡模拟控制方程基础上,利用基于三阶迎风偏置格式的混合差分及二阶中心差分离散控制方程的对流项及耗散项,离散得到的线性方程组通过无矩阵计算的隐式LU-SGS(Lower-Upper SymmetricGauss Seidel)分步迭代方法求解。利用大涡模拟(Large Eddy Simulation-LES)方法对水下航行体周围的黏性流场及水动力进行了定常数值模拟,计算得到的速度、压力分布及阻力值与试验吻合较好。 展开更多
关键词 大涡模拟(LES) 亚格子应力模型 计算流体力学(CFD) 水下航行体
下载PDF
涡轮桨搅拌槽内混合过程的大涡模拟 被引量:2
5
作者 苗一 潘家祯 +1 位作者 闵健 高正明 《华东理工大学学报(自然科学版)》 EI CAS CSCD 北大核心 2006年第5期623-628,共6页
在FLUENT 6.1软件平台和网络并行计算硬件平台上,采用大涡模拟(LES)的方法对涡轮桨搅拌槽内的混合过程进行了数值模拟。利用滤波函数对N av ier-Stokes方程进行空间滤波,对大尺寸的涡直接进行求解,而被滤掉的比网格小的旋涡通过Sam agor... 在FLUENT 6.1软件平台和网络并行计算硬件平台上,采用大涡模拟(LES)的方法对涡轮桨搅拌槽内的混合过程进行了数值模拟。利用滤波函数对N av ier-Stokes方程进行空间滤波,对大尺寸的涡直接进行求解,而被滤掉的比网格小的旋涡通过Sam agorinsky-L illy亚格子模型求解,对搅拌桨区域采用滑移网格技术。结果表明:大涡模拟对尾涡的预报优于雷诺平均(RAN S)方法,混合时间以及示踪剂响应曲线模拟结果和实验结果吻合较好,且优于RAN S方法。大涡模拟方法为准确预测搅拌槽内湍流流动的非稳态及周期性脉动特性提供了一种有效的工具。 展开更多
关键词 大涡模拟 混合时间 亚格子模型 计算流体力学(CFD) 搅拌槽
下载PDF
一种固体火箭发动机喷管摆动冷流缩比试验方法 被引量:2
6
作者 郭小帆 杨晓慧 +1 位作者 何德胜 李广武 《中国测试》 CAS 北大核心 2013年第5期39-41,共3页
针对研究喷管摆动工况下扩压器工作特性的试验中存在:喷管摆动难以实现;若以高温高压蒸汽或者燃气作为引射介质,虽引射效率高但耗费巨大的人力物力等问题在发动机研制初期很难被采纳。该文提出一种以扩压器摆动等效喷管摆动,以冷空气作... 针对研究喷管摆动工况下扩压器工作特性的试验中存在:喷管摆动难以实现;若以高温高压蒸汽或者燃气作为引射介质,虽引射效率高但耗费巨大的人力物力等问题在发动机研制初期很难被采纳。该文提出一种以扩压器摆动等效喷管摆动,以冷空气作为引射和被引射介质的简化试验方法。数值计算结果与试验结果符合较好,试验方法经济实用,操作简单,可以作为摆动喷管固体火箭发动机研制初期的试验手段。 展开更多
关键词 固体火箭发动机 冷流缩比试验 喷管摆动 试验方法
下载PDF
大型飞机的极限飞行状态研究
7
作者 颜巍 王磊 +1 位作者 李永涛 赵晶慧 《民用飞机设计与研究》 2020年第3期54-60,共7页
根据国家民用航空局的适航条款与咨询通告要求,飞机研制单位需表明飞机的失速/过失速特性以及改出特性,该项试验属于飞机的极限边界飞行和超包线飞行。原型机实施该项试验存在极大的风险,历史上出现过多起机毁人亡的重大事故。原型机缩... 根据国家民用航空局的适航条款与咨询通告要求,飞机研制单位需表明飞机的失速/过失速特性以及改出特性,该项试验属于飞机的极限边界飞行和超包线飞行。原型机实施该项试验存在极大的风险,历史上出现过多起机毁人亡的重大事故。原型机缩比模型自由飞试验可以捕捉到模型静态测力试验条件下所未发现的飞行动态特性,以及检验飞机在危险状态下的改出能力,降低原型机失速试飞的风险。某大型飞机在失速试飞前利用模型自由飞试验技术检验了飞机在极限飞行条件下的动态特性,确认了飞机在失速条件下存在偏离以及发展稳定尾旋等;同时确认了飞机具备失速改出能力,采用"三中立"、"反舵,推杆"和"实施开反尾旋伞"均能迫使飞机改出失速/过失速状态。 展开更多
关键词 缩比模型自由飞 失速试飞 过失速 反尾旋伞
下载PDF
双圆柱绕流的大涡模拟 被引量:35
8
作者 贾晓荷 刘桦 《水动力学研究与进展(A辑)》 CSCD 北大核心 2008年第6期625-632,共8页
基于FLUENT软件平台上的三维大涡模拟方法,对定常来流流场中的单圆柱及不同排列形式的双圆柱进行了水动力计算,并对圆柱周围的瞬时流场、时均流场以及圆柱受力进行了分析。分析结果表明,当Re=6×104时,串列双圆柱中上游圆柱受力明... 基于FLUENT软件平台上的三维大涡模拟方法,对定常来流流场中的单圆柱及不同排列形式的双圆柱进行了水动力计算,并对圆柱周围的瞬时流场、时均流场以及圆柱受力进行了分析。分析结果表明,当Re=6×104时,串列双圆柱中上游圆柱受力明显大于下游圆柱,下游圆柱Cd、Cl的变化频率满足fd=2fl,两圆柱的脱涡频率近似相等,涡街的相位相同;并列双圆柱的Cd基本相等,且变化趋势有一定相似性,脱涡频率均接近单个圆柱,涡街的相位相反;交错45°排列的两个圆柱Cd基本相等,脱涡频率接近,下游圆柱的涡街相位明显滞后于上游圆柱。 展开更多
关键词 大涡模拟 亚格子模式 圆柱绕流
下载PDF
大涡模拟理论进展及其在工程中的应用 被引量:15
9
作者 王汉青 王志勇 寇广孝 《流体机械》 CSCD 2004年第7期23-27,共5页
介绍了大涡模拟的理论进展和发展趋势 ,给出了大涡模拟的常用滤波函数和模型。针对国内外的研究现状 ,描述了当前大涡模拟在工程中的具体应用。
关键词 大涡模拟 湍流 数值模拟 亚格子模式
下载PDF
列车车厢内火灾烟气运动的数值模拟研究 被引量:17
10
作者 刘采峰 彭岚 刘朝 《热科学与技术》 CAS CSCD 2003年第4期352-357,共6页
使用大涡模拟和 Smagorinsky亚格子尺度模型 ,以及简单化学快速反应燃烧模型 ,对列车车厢内火灾现象进行了数值模拟。由于热烟气的弱可压缩性 ,采用了适合于低马赫数形式的三维非稳态 Navier-Stokes方程的近似形式数值研究烟气的输运过... 使用大涡模拟和 Smagorinsky亚格子尺度模型 ,以及简单化学快速反应燃烧模型 ,对列车车厢内火灾现象进行了数值模拟。由于热烟气的弱可压缩性 ,采用了适合于低马赫数形式的三维非稳态 Navier-Stokes方程的近似形式数值研究烟气的输运过程。通过对不同释热率情况下火灾过程的研究 ,初步探讨了车厢内火灾发展、烟气运动和温度分布规律。采用的研究方法和所得结论对于列车火灾的防护设计和安全管理方面有一定指导意义。 展开更多
关键词 列车车厢 火灾 烟气运动 数值模拟 大涡模拟 烟气输运 亚格子尺度模型
下载PDF
大涡模拟及其在室内空气流动换热中的应用 被引量:1
11
作者 谭羽非 《暖通空调》 北大核心 2003年第1期99-102,共4页
介绍了国内外大涡模拟的亚格子湍流模式和在室内空气流动换热中的研究现状 ,给出了大涡模拟室内空气流动的应用实例。认为大涡模拟有广阔的发展前景。
关键词 大涡模拟 亚格子湍流模式 室内 空气流动 发展前景
下载PDF
基于动态亚格子模型的方柱绕流大涡模拟 被引量:7
12
作者 张海涛 曹曙阳 《沈阳建筑大学学报(自然科学版)》 CAS 北大核心 2013年第3期434-439,共6页
目的对高雷诺数下的方柱绕流进行数值模拟,得到方柱下流场时间平均结果以及湍流流场的流场特征,用于验证亚格子模型对大涡模拟(LES)计算精确性的影响.方法基于有限体积法(FVM)的流体计算开源软件OpenFOAM,采用两种不同动态亚格子模型下... 目的对高雷诺数下的方柱绕流进行数值模拟,得到方柱下流场时间平均结果以及湍流流场的流场特征,用于验证亚格子模型对大涡模拟(LES)计算精确性的影响.方法基于有限体积法(FVM)的流体计算开源软件OpenFOAM,采用两种不同动态亚格子模型下的大涡模拟方法,模拟高雷诺数(Re=2.2×104)下的方柱绕流.结果通过数值计算得到流场完整的平均和脉动结构,并得到阻力系数、升力系数、斯托拉哈数等重要的气动力参数,经对比与试验实际情况吻合.结论高雷诺数复杂流场下亚格子应力的封闭形式对大涡模拟的计算结果有重要的影响,采用动态模拟技术一方程模型计算结果优于常规的Smagorinsky模型,在高雷诺数复杂流场下的计算是高效的,为结构抗风的数值研究提供了参考. 展开更多
关键词 大涡模拟 方柱绕流 亚格子模型 尾流结构 数值模拟
下载PDF
不同亚网格尺度应力模型在燃油喷雾大涡模拟中的应用 被引量:7
13
作者 周磊 解茂昭 +1 位作者 贾明 史俊瑞 《内燃机学报》 EI CAS CSCD 北大核心 2011年第1期29-35,共7页
亚网格尺度应力模型是影响大涡模拟结果的一个关键因素.在K-方程亚网格湍动能模型的基础上,对其它比较常用的3种代数模型(Smagoringsky模型、动态Smagoringsky模型和WALE模型)在燃油喷雾中的适用性作了探讨.在定容弹中采用3个不同的燃... 亚网格尺度应力模型是影响大涡模拟结果的一个关键因素.在K-方程亚网格湍动能模型的基础上,对其它比较常用的3种代数模型(Smagoringsky模型、动态Smagoringsky模型和WALE模型)在燃油喷雾中的适用性作了探讨.在定容弹中采用3个不同的燃油喷雾试验作为数值模拟的验证对象,分析了喷雾的演变过程.并将其应用于柴油发动机喷雾冷态研究中.结果可知,3个代数亚网格模型在蒸发和非蒸发过程中喷雾的演变过程基本与试验相符,液相贯穿距也与试验接近.K-方程模型较为准确,WALE模型在几种算例中也得到了较好的结果,Smagorinsky模型得到的液相贯穿距要稍长于其它几种模型. 展开更多
关键词 大涡模拟 亚网格模型 燃油喷雾
下载PDF
文丘里喷嘴空化流大涡模拟亚格子模型对比研究
14
作者 张楚谦 郑平 +1 位作者 陈勇刚 赵梁 《液压气动与密封》 2024年第5期32-38,共7页
针对文丘里喷嘴不同机制下可压缩空化流开展了大涡模拟(LES)研究,并对比讨论了3种常用的亚格子模型的适用性。结果显示,各模型在预测回射流机制和凝结激波机制的流速与空化云脱落周期上表现良好。在空化云演化过程中,WALE模型在两种机... 针对文丘里喷嘴不同机制下可压缩空化流开展了大涡模拟(LES)研究,并对比讨论了3种常用的亚格子模型的适用性。结果显示,各模型在预测回射流机制和凝结激波机制的流速与空化云脱落周期上表现良好。在空化云演化过程中,WALE模型在两种机制下的模拟吻合度最佳,SL模型预测溃灭过程提前,蒸汽相体积分数偏小;KET模型吻合度最差,溃灭过程以及周期内脱落存在时间延迟。初步发现回射流机制下的喉部压力功率谱密度分析服从-5/3标度律,凝结激波机制服从-7/3标度律。 展开更多
关键词 大涡模拟 亚格子模型 空化 空化流 数值模拟
下载PDF
壁面涡旋结构与湍流脉动压力的大涡模拟研究 被引量:6
15
作者 张晓龙 张楠 吴宝山 《船舶力学》 EI CSCD 北大核心 2014年第8期871-881,共11页
湍流脉动压力是重要的流噪声声源,对其进行数值计算是流声耦合领域的重要课题,开展相应的研究十分必要。文章采用大涡模拟方法(LES)结合四种亚格子涡模型与四套网格,对槽道壁面湍流脉动压力进行了数值计算,并与试验结果进行了对比分析,... 湍流脉动压力是重要的流噪声声源,对其进行数值计算是流声耦合领域的重要课题,开展相应的研究十分必要。文章采用大涡模拟方法(LES)结合四种亚格子涡模型与四套网格,对槽道壁面湍流脉动压力进行了数值计算,并与试验结果进行了对比分析,验证了数值计算方法的可靠性。首先,介绍了大涡模拟的物理内涵与基本方程,给出了常用亚格子涡模型的表达式,并给出了相应的离散求解数值方法以及边界条件的设置。其次,描述了槽道试验段的几何特征,给出了网格的剖分形式。最后,详细讨论了槽道壁面湍流脉动压力频谱计算值与试验值之间的差异,进行了定量与定性的比较分析,同时分析了涡旋结构与近壁面流速分布,研究了亚格子涡模型与网格数量对计算结果的影响,为今后复杂几何模型壁面湍流脉动压力及其频率-波数谱的计算研究工作奠定了基础。 展开更多
关键词 壁面湍流脉动压力 大涡模拟 槽道 亚格子涡模型
下载PDF
基于三种亚格子模型的空腔振荡流动计算 被引量:5
16
作者 白海涛 赖焕新 《华东理工大学学报(自然科学版)》 CAS CSCD 北大核心 2016年第1期125-131,共7页
使用三种亚格子应力模型,对长深比(L/D)为5的三维矩形开式空腔的可压缩流体进行大涡模拟计算。研究得到的空腔自激振荡频率与Rossiter公式计算结果和实验结果吻合良好,结果显示振荡能量主要集中在较低频率区域,压力幅值主要出现在前三... 使用三种亚格子应力模型,对长深比(L/D)为5的三维矩形开式空腔的可压缩流体进行大涡模拟计算。研究得到的空腔自激振荡频率与Rossiter公式计算结果和实验结果吻合良好,结果显示振荡能量主要集中在较低频率区域,压力幅值主要出现在前三阶模态。Dynamic Smagorinsky-Lilly(DSM)模型在空腔前后壁面附近区域的脉动强度分布比Smagorinsky-Lily(SM)模型更为接近实验值,Wall Adapting Local Eddy Viscosity(WALE)模型的脉动强度分布与实验值最为接近。由空腔底部监测点声压级分布及声压频谱图可以看出:WALE模型性能最佳,DSM模型结果也与实验结果相符合,SM模型的预测性能略差。 展开更多
关键词 开式空腔 自激振荡 大涡模拟 亚格子应力模型 气动噪声
下载PDF
基于大涡模拟方法的风力机气动性能和尾流干扰研究 被引量:5
17
作者 钱耀如 王同光 张震宇 《中国科学:物理学、力学、天文学》 CSCD 北大核心 2016年第12期23-30,共8页
利用致动线(Actuator Line Method,ALM)和大涡模拟(Large-Eddy Simulation,LES)结合的方法,采用4种亚格子模型,对低湍流度均匀来流中不同转速下两台串列风力机的气动性能和尾流干扰开展数值模拟研究,并探讨了亚格子模型对尾流场模拟的影... 利用致动线(Actuator Line Method,ALM)和大涡模拟(Large-Eddy Simulation,LES)结合的方法,采用4种亚格子模型,对低湍流度均匀来流中不同转速下两台串列风力机的气动性能和尾流干扰开展数值模拟研究,并探讨了亚格子模型对尾流场模拟的影响.两台风力机功率系数和推力系数,以及尾流场的轴向平均速度和雷诺应力分布的计算结果与实验值基本吻合,验证了ALM-LES方法对风力机尾流研究的可靠性.受上游风力机尾流的影响,下游风力机功率系数和推力系数大幅降低,最大功率系数仅为上游风力机最大功率系数的25%左右.与来流风况的低湍流度相比,风力机尾流场中湍流强度大幅提高.通过不同亚格子模型计算结果的对比分析,得出亚格子模型的选择对风力机气动性能和尾流场湍流特征参数的计算影响较小. 展开更多
关键词 大涡模拟 亚格子模型 风力机尾流 致动线方法
原文传递
二阶双系数动态亚格子应力模型 被引量:4
18
作者 唐学林 钱忠东 吴玉林 《水科学进展》 EI CAS CSCD 北大核心 2004年第1期50-55,共6页
基于亥姆霍兹速度分解定理和Smagorinsky模型,提出二阶双系数动态亚格子尺度应力模型,即亚格子尺度应力是可解的应变率张量和旋转率张量的函数。应用有限差分法对此模型的控制方程进行了离散,并数值模拟得到了弯管内部的流场分布和压力... 基于亥姆霍兹速度分解定理和Smagorinsky模型,提出二阶双系数动态亚格子尺度应力模型,即亚格子尺度应力是可解的应变率张量和旋转率张量的函数。应用有限差分法对此模型的控制方程进行了离散,并数值模拟得到了弯管内部的流场分布和压力分布。把雷诺数为40000的弯管流动的实验数据和此模型的模拟结果相比较,证明此模型是可行、可靠的。 展开更多
关键词 亚尺度模型 动态模型 湍流 弯管 二阶双系数 亚格子应力
下载PDF
一种改进的类DES湍流模拟方法 被引量:4
19
作者 汪洪波 孙明波 +1 位作者 吴海燕 王振国 《航空动力学报》 EI CAS CSCD 北大核心 2011年第10期2167-2173,共7页
构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格... 构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法. 展开更多
关键词 混合RANS/LES(Reynolds-averaged Navier-Stokes/large EDDY simulation) 混合函数 凹腔 超声速 亚格子模型
原文传递
不同亚格子尺度模型对气动噪声仿真精度的影响 被引量:2
20
作者 殷想 汪怡平 +2 位作者 杜敏韬 苏楚奇 孙浩 《实验流体力学》 EI CAS CSCD 北大核心 2019年第6期78-83,共6页
大涡模拟(Large Eddy Simulation,LES)湍流模型在汽车外流场及气动噪声模拟计算分析中有着广泛的应用,高精度的非定常流场计算结果是有效进行气动噪声仿真分析的前提和基础。采用现代简易汽车模型(HSM模型),基于不同亚格子模型的LES对... 大涡模拟(Large Eddy Simulation,LES)湍流模型在汽车外流场及气动噪声模拟计算分析中有着广泛的应用,高精度的非定常流场计算结果是有效进行气动噪声仿真分析的前提和基础。采用现代简易汽车模型(HSM模型),基于不同亚格子模型的LES对其进行非定常流场计算分析以及声场计算,并利用模型内部监测点的实验数据和仿真计算的结果进行对比验证。结果表明,使用WMLES模型可以得到更为准确的结果。 展开更多
关键词 气动噪声 亚格子模型 大涡模拟 HSM模型
下载PDF
上一页 1 2 3 下一页 到第
使用帮助 返回顶部