期刊文献+
共找到39篇文章
< 1 2 >
每页显示 20 50 100
基于Infoworks ICM水力模型的圩区排涝非恒定流数值模拟研究
1
作者 王晶晶 陆雨杰 许鸣欢 《中国水能及电气化》 2025年第5期26-29,44,共5页
为了预测圩区在排涝过程中的水力行为,提出基于Infoworks ICM水力模型的圩区排涝非恒定流数值模拟研究。利用Infoworks ICM软件构建圩区排涝系统水动力模型,设定各项模型的参数,模拟不同降雨情景下的水力过程,设计降雨方案,实现圩区排... 为了预测圩区在排涝过程中的水力行为,提出基于Infoworks ICM水力模型的圩区排涝非恒定流数值模拟研究。利用Infoworks ICM软件构建圩区排涝系统水动力模型,设定各项模型的参数,模拟不同降雨情景下的水力过程,设计降雨方案,实现圩区排涝非恒定流数值模拟。模拟结果表明,该研究方法能够准确模拟圩区排涝系统的水力特性,模拟结果均与实际的水位数据基本相符,并揭示了不同降雨情景下系统性能的显著差异,具有应用价值。 展开更多
关键词 Infoworks ICM水力模型 圩区排涝 非恒定流数值模拟 水动力模型 降雨方案设计
在线阅读 下载PDF
非定常空化流动研究现状与进展 被引量:31
2
作者 黄彪 吴钦 王国玉 《排灌机械工程学报》 EI CSCD 北大核心 2018年第1期1-14,共14页
针对非定常空化流动数值计算模型,总结了不同气液相间传输模型对模拟空穴发展和空泡脱落过程的影响,发现不同空化模型对于非定常空化流场的不同区域的预测有着明显的差异.建立了一种基于混合密度分域的空化流动计算模型,合理预测了空泡... 针对非定常空化流动数值计算模型,总结了不同气液相间传输模型对模拟空穴发展和空泡脱落过程的影响,发现不同空化模型对于非定常空化流场的不同区域的预测有着明显的差异.建立了一种基于混合密度分域的空化流动计算模型,合理预测了空泡脱落时刻翼型尾部的不稳定脉动区.通过考虑空化流动的湍流多尺度与多相局部可压缩特性,提出了一种基于混合密度分域的混合湍流模型,预测结果与试验结果吻合更好.综合回顾试验与数值计算方面的研究工作,围绕空化流动特性与动力学问题,主要从非定常流场空穴形态、流动结构、空化流体动力等几个方面对非定常空化流动机理及相应研究成果进行了较为详细地评述和讨论.最后,基于非定常空化研究的发展趋势和进一步研究可能面临的突出问题,对非定常空化流动研究的发展方向进行了探讨与展望. 展开更多
关键词 非定常空化 试验测试技术 数值计算模型 空化流动机理
在线阅读 下载PDF
平原河网水动力模型及求解方法探讨 被引量:44
3
作者 卢士强 徐祖信 《水资源保护》 CAS 2003年第3期5-9,共5页
参考国内外有关资料 ,根据河网非恒定流水动力模型的控制方程组和汊点衔接条件 ,建立平原河网水动力节点 河道模型。介绍用特征线法 ,有限体积法和有限差分法中的直接解法、分级解法、汊点分组解法、矩阵标识法、非线性方法等求解的基... 参考国内外有关资料 ,根据河网非恒定流水动力模型的控制方程组和汊点衔接条件 ,建立平原河网水动力节点 河道模型。介绍用特征线法 ,有限体积法和有限差分法中的直接解法、分级解法、汊点分组解法、矩阵标识法、非线性方法等求解的基本思路。对比单元划分模型、混合模型以及人工神经网络模型等平原河网水动力模型 ,分析各个模型的优缺点 ,结果表明 ,节点 河道模型原则上可以求解任何水网的水力参数 ,单元划分模型仅适用于河道流速时空变化不大的情况 ,人工神经网络模型的验证比较困难。指出改进和设计计算方法、应用向量运算和并行算法、数值模拟可视化。 展开更多
关键词 非恒定流 水动力模型 数值模拟 人工神经网络 河网
在线阅读 下载PDF
复杂内边界长距离输水明渠的一维非恒定流数学模型 被引量:32
4
作者 张成 傅旭东 王光谦 《南水北调与水利科技》 CAS CSCD 2007年第6期16-20,共5页
大型跨流域调水工程,通常具有输水线路长、沿程过水建筑物类型和数量多的特点。节制闸、分水口及倒虹吸、渡槽等建筑物与水流的复杂作用,造成了长距离输水明渠非恒定流模拟的困难。建立了复杂内边界长距离输水明渠的一维非恒定流数学模... 大型跨流域调水工程,通常具有输水线路长、沿程过水建筑物类型和数量多的特点。节制闸、分水口及倒虹吸、渡槽等建筑物与水流的复杂作用,造成了长距离输水明渠非恒定流模拟的困难。建立了复杂内边界长距离输水明渠的一维非恒定流数学模型。模型定义和模化了节制闸、倒虹吸和多孔并联渡槽等3种类型的内边界条件,实现了对输水明渠水力过渡过程中水位、流量大幅度变化的无间断连续模拟。模型应用于南水北调中线一期工程总干渠全线和单个闸门的模拟分析,验证了模型的合理性。结果表明,模型能够模拟多内边界复杂输水渠道的水力响应过程,可作为南水北调中线一期工程等跨流域大型明渠输水工程的模拟工具。 展开更多
关键词 南水北调 非恒定流 数学模型 内边界 明渠
在线阅读 下载PDF
非恒定水流的一维、二维耦合数值模型 被引量:27
5
作者 赖锡军 汪德爟 《水利水运工程学报》 CSCD 北大核心 2002年第2期48-51,共4页
从数值波在每一时间步长内传播的距离不超过Courant数个网格点的基础出发 ,采用时间滞后条件 ,构造了虚拟的重叠计算水域 .利用重叠 -投影法 ,在连接 (点 )处建立一、二维耦合模型 ,实现了由一维计算水域向二维计算水域的准确过渡 .采... 从数值波在每一时间步长内传播的距离不超过Courant数个网格点的基础出发 ,采用时间滞后条件 ,构造了虚拟的重叠计算水域 .利用重叠 -投影法 ,在连接 (点 )处建立一、二维耦合模型 ,实现了由一维计算水域向二维计算水域的准确过渡 .采用一、二维耦合模型和全一维模型对上海市淀南片河网的验证计算表明 ,一、二维耦合模型能较好地模拟非恒定水流 。 展开更多
关键词 重叠-投影法 耦合模型 数值计算 非恒定水流
在线阅读 下载PDF
三峡水库运用初期库区水沙输移数值模拟 被引量:12
6
作者 黄仁勇 李飞 张细兵 《长江科学院院报》 CSCD 北大核心 2012年第1期7-12,共6页
考虑水沙输移过程中的非恒定性,建立了三峡水库一维非恒定水沙数学模型。此模型为一树状河网模型,可以将更多支流纳入计算范围。采用三峡水库蓄水运用后2003-2009年实测水沙资料对模型进行了验证计算,计算结果与实测资料符合较好,表明... 考虑水沙输移过程中的非恒定性,建立了三峡水库一维非恒定水沙数学模型。此模型为一树状河网模型,可以将更多支流纳入计算范围。采用三峡水库蓄水运用后2003-2009年实测水沙资料对模型进行了验证计算,计算结果与实测资料符合较好,表明模型能够用于三峡水库干支流河道水沙输移数值模拟。该模型可在今后三峡水库优化调度等相关研究工作中发挥其模拟和预测作用。 展开更多
关键词 一维 非恒定 数学模型 三峡水库 模型验证
在线阅读 下载PDF
排水管网非恒定流数值模拟新方法 被引量:3
7
作者 茅泽育 赵雪峰 +1 位作者 赵璇 罗昇 《水利水运工程学报》 CSCD 北大核心 2007年第2期42-47,共6页
提出了数值模拟城市排水管网非恒定流新的计算方法,即采用扩散波模型,将待求水位和流量分为运动波分量和扩散波分量两部分,并先后求解.在处理明满流过渡时,没有采用Pre issm ann明窄缝近似方法.用原型观测资料验证计算结果表明,计算值... 提出了数值模拟城市排水管网非恒定流新的计算方法,即采用扩散波模型,将待求水位和流量分为运动波分量和扩散波分量两部分,并先后求解.在处理明满流过渡时,没有采用Pre issm ann明窄缝近似方法.用原型观测资料验证计算结果表明,计算值与实测值吻合较好. 展开更多
关键词 扩散波模型 非恒定流 明满流 排水管网 数值计算
在线阅读 下载PDF
三峡水库干支流河道一维非恒定水沙数学模型初步研究 被引量:15
8
作者 黄仁勇 黄悦 《长江科学院院报》 CSCD 北大核心 2009年第2期9-13,共5页
考虑水沙运动过程中的非恒定性及众多支流的影响,建立了三峡水库干支流河道一维非恒定水沙数学模型,并对水流方程组采用三级解法求解。采用三峡水库蓄水运用后2003年6月至2004年12月的水库实测水沙资料对模型进行了验证,验证结果有良好... 考虑水沙运动过程中的非恒定性及众多支流的影响,建立了三峡水库干支流河道一维非恒定水沙数学模型,并对水流方程组采用三级解法求解。采用三峡水库蓄水运用后2003年6月至2004年12月的水库实测水沙资料对模型进行了验证,验证结果有良好的一致性。本模型在使用1961-1970年系列年入库水沙条件及水库正常蓄水位调度方案下进行了水库冲淤100 a预测计算,计算结果较为合理,表明本模型可用于三峡水库水沙数值模拟。 展开更多
关键词 一维 非恒定 三级解法 数学模型 三峡水库
在线阅读 下载PDF
垂直轴风机三维非定常流场与性能的数值模拟 被引量:4
9
作者 姚激 黄剑峰 +2 位作者 袁伟斌 彭明军 孙勇 《太阳能学报》 EI CAS CSCD 北大核心 2013年第7期1220-1227,共8页
应用Fluent软件求解三维非定常不可压N-S方程和RNG k-ε湍流模型,采用Simple算法,结合滑动网格技术,对垂直轴风机的三维非定常流场特性进行数值模拟研究,获得并分析不同时刻计算区域的速度场、压力场分布以及叶片上的压力分布与风轮的... 应用Fluent软件求解三维非定常不可压N-S方程和RNG k-ε湍流模型,采用Simple算法,结合滑动网格技术,对垂直轴风机的三维非定常流场特性进行数值模拟研究,获得并分析不同时刻计算区域的速度场、压力场分布以及叶片上的压力分布与风轮的转矩的变化规律。研究结果表明:模型计算区域内的速度分布、流场静压、叶片表面压力分布均呈良好周期性;垂直轴风机风轮叶片周围速度梯度明显,风轮叶片在不同迎角上的作用机理不同,每个周期内垂直轴风机的总转矩曲线上升段较下降段平缓;垂直轴风机三维非定常流场的模拟结果能有效反映风机旋转过程的流场、叶片压力和转矩变化状况。 展开更多
关键词 垂直轴风机 三维非定常流场 数值模拟 RNG k-ε湍流模型 SIMPLE算法 滑移网格
在线阅读 下载PDF
多闸坝分汊河流的洪水实时模拟和调度 被引量:5
10
作者 陈明洪 方红卫 刘军梅 《水利水电科技进展》 CSCD 北大核心 2011年第2期11-15,27,共6页
采用一维非恒定水流数学模型耦合闸坝水力学模型,构建了可用于多闸坝分汊河流中场次洪水实时模拟和调度的模型。该模型分别对河流分汊、多闸坝泄水和断流进行处理,增强了模型的实时性和稳定性。基于该模型对北运河杨洼闸以上河道进行洪... 采用一维非恒定水流数学模型耦合闸坝水力学模型,构建了可用于多闸坝分汊河流中场次洪水实时模拟和调度的模型。该模型分别对河流分汊、多闸坝泄水和断流进行处理,增强了模型的实时性和稳定性。基于该模型对北运河杨洼闸以上河道进行洪水模拟和调度,验证结果表明模拟结果令人满意。根据不同重现期场次洪水演进结果提出的调度方案为北运河重要闸坝的调度决策提供了技术支持。 展开更多
关键词 非恒定流 一维数学模型 多闸坝河流 分汊河流 洪水调度
在线阅读 下载PDF
牛顿流体圆管内非稳态Poiseuille流动特性 被引量:3
11
作者 陈雷 汤苑楠 +1 位作者 刘刚 卢兴国 《中国石油大学学报(自然科学版)》 EI CAS CSCD 北大核心 2018年第3期114-121,共8页
Poiseuille流动初始阶段存在速度发展的非稳态过程,会对测试结果造成偏差。为分析非稳态Poiseuille流动过程对测量黏度造成的偏差,以不可压缩牛顿流体为例,进行恒流量边界与非定常流量边界下的非稳态Poiseuille流动过程研究。以无量纲... Poiseuille流动初始阶段存在速度发展的非稳态过程,会对测试结果造成偏差。为分析非稳态Poiseuille流动过程对测量黏度造成的偏差,以不可压缩牛顿流体为例,进行恒流量边界与非定常流量边界下的非稳态Poiseuille流动过程研究。以无量纲黏度和无量纲时间表征非稳态过程,建立数值模型,计算给出恒平均速度边界、从0线性增加平均速度边界和恒压力边界条件下无量纲黏度数值的变化规律。结果表明:非稳态过程中无量纲黏度数值随时间逐渐减小并最终趋于1,且不同边界条件下流动达到稳定对应的无量纲时间为定值。当边界类型确定时,非稳态过程的无量纲黏度数值可视为仅与无量纲时间有关的函数;对于不同类型边界条件,从0线性增加的平均速度边界、恒压力边界、恒平均速度边界条件对应的非稳态过程逐渐缩短。 展开更多
关键词 流变学 非稳态Poiseuille流动 数值模型 边界条件 无量纲黏度 无量纲时间
在线阅读 下载PDF
低压涡轮叶型边界层相互作用的数值模拟 被引量:5
12
作者 罗华玲 乔渭阳 《推进技术》 EI CAS CSCD 北大核心 2009年第1期95-100,共6页
基于Lantry-Menter转捩模型,应用商用流体计算软件求解非定常雷诺平均N-S方程组,对进口雷诺数为2.5×104,来流湍流度为2.5%,尾迹折合频率为0.92状态下的尾迹/PAK-B叶型边界层的相互作用进行了数值模拟。数值计算揭示的尾迹在叶栅通... 基于Lantry-Menter转捩模型,应用商用流体计算软件求解非定常雷诺平均N-S方程组,对进口雷诺数为2.5×104,来流湍流度为2.5%,尾迹折合频率为0.92状态下的尾迹/PAK-B叶型边界层的相互作用进行了数值模拟。数值计算揭示的尾迹在叶栅通道中的输运特性、尾迹诱导卷起涡形成、尾迹诱导边界层转捩等物理机制大致符合相关实验定性描述。 展开更多
关键词 低压涡轮 非定常流 尾流 转捩模型+ 数值仿真
在线阅读 下载PDF
绕钝体分离流动的非定常数值模拟研究 被引量:7
13
作者 陈晓春 王元 《空气动力学学报》 CSCD 北大核心 2007年第4期513-520,共8页
利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流... 利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。 展开更多
关键词 RNG k-ε模型 钝体绕流 非定常数值模拟
在线阅读 下载PDF
基于POD方法的二维方柱低雷诺数绕流流场分析研究 被引量:14
14
作者 王掩刚 陈俊旭 先松川 《西北工业大学学报》 EI CAS CSCD 北大核心 2014年第4期612-617,共6页
以二维方柱低雷诺数绕流为研究对象,在非定常数值模拟结果的基础上,首先应用本征正交分解(proper orthogonal decomposition,POD)方法得到特征区域流场的低阶模型,获取POD基函数及时间系数分布规律;其次以POD的低阶模型为基础,实现了对... 以二维方柱低雷诺数绕流为研究对象,在非定常数值模拟结果的基础上,首先应用本征正交分解(proper orthogonal decomposition,POD)方法得到特征区域流场的低阶模型,获取POD基函数及时间系数分布规律;其次以POD的低阶模型为基础,实现了对非定常物理场的重构。研究结果表明:1应用文中的数值分析手段,捕捉到了方柱绕流非定常流动特征,且与公开文献结果吻合较好。2对于文中的研究对象,前4阶POD模态所占波动总能量99.4%,可以用于较准确描述其非定常流场特征;3各阶POD模态的时间系数之间都有明确的频率与相位关系,其中能量最大的第一对POD模态对方柱尾迹流动的波动频率和幅值起决定作用。4在所研究的雷诺数下,前4阶POD模态就可以很好地重构流场,这对以后低阶模型的建立有一定的指导意义。 展开更多
关键词 本征正交分解 低阶模型 方柱绕流 非定常
在线阅读 下载PDF
CFD在调压室涌浪水位模拟中的应用 被引量:5
15
作者 刘飞 杨建东 李进平 《中国农村水利水电》 北大核心 2011年第9期123-126,129,共5页
对某水电站的尾水隧洞及下游调压室过渡过程中的非恒定流流态进行三维数值模拟,着重观测调压室内涌浪波动及调压室内底板压差变化过程。采用3种湍流模型对下游调压室及尾水隧洞进行了模拟,将调压室涌浪波动过程与调压室底板压差变化过... 对某水电站的尾水隧洞及下游调压室过渡过程中的非恒定流流态进行三维数值模拟,着重观测调压室内涌浪波动及调压室内底板压差变化过程。采用3种湍流模型对下游调压室及尾水隧洞进行了模拟,将调压室涌浪波动过程与调压室底板压差变化过程的三维计算结果与物理模型实验、一维数值计算的相应结果予以对比分析,阐明了不同湍流模型对过渡过程模拟的影响。结合三维模拟能捕捉流场内部流态的优势,探讨了涌浪波动过程中调压室水面波动过程及水流流态。 展开更多
关键词 三维流场 湍流模型 非恒定流 调压室涌浪 调压室底板压差
在线阅读 下载PDF
旋翼桨尖气动/降噪综合优化设计研究 被引量:9
16
作者 郭旺柳 宋文萍 +1 位作者 许建华 许瑞飞 《西北工业大学学报》 EI CAS CSCD 北大核心 2012年第1期73-79,共7页
采用Kriging模型和遗传算法开展了针对桨尖形状的气动/声学综合优化设计研究。为了获得精确的气动及声压特性,得到基于声类比混合方法计算气动噪声需要的声源数据,采用非定常雷诺平均NS方程(URANS)数值模拟旋翼前飞粘性绕流。为了提高... 采用Kriging模型和遗传算法开展了针对桨尖形状的气动/声学综合优化设计研究。为了获得精确的气动及声压特性,得到基于声类比混合方法计算气动噪声需要的声源数据,采用非定常雷诺平均NS方程(URANS)数值模拟旋翼前飞粘性绕流。为了提高遗传算法的优化效率,采用一种基于EI方法的Kriging模型代替费时的数值模拟过程。以AH-1G/OLS旋翼为基准,并以气动性能为约束,噪声峰值最小为目标,进行了旋翼桨尖的降噪优化设计。优化结果表明,文中多发展的优化方法是可行的。 展开更多
关键词 旋翼 噪声 优化设计 雷诺平均NS方程 KRIGING模型
在线阅读 下载PDF
矿井涌水量预测三维数值模拟 被引量:4
17
作者 马青山 骆祖江 《西安科技大学学报》 CAS 北大核心 2015年第2期236-241,共6页
为了准确预测矿井涌水量,保障煤矿安全生产,以鄂尔多斯巴彦淖井田为例,采用地下水三维数值模拟理论和方法,通过对研究区水文地质模型的概化,建立了巴彦淖井田矿井涌水量预测的地下水三维非稳定流数值模拟模型,并结合矿井生产进度,以工... 为了准确预测矿井涌水量,保障煤矿安全生产,以鄂尔多斯巴彦淖井田为例,采用地下水三维数值模拟理论和方法,通过对研究区水文地质模型的概化,建立了巴彦淖井田矿井涌水量预测的地下水三维非稳定流数值模拟模型,并结合矿井生产进度,以工作面月回采进度为单位,模拟预测了丰水期和平水期两种情况下,各工作面不同进度期地下水位分别疏降至2煤层底时的涌水量。结果表明:在前10个进度期内,工作面最大涌水量为1 328 m3·h-1,正常涌水量为1 134m3·h-1.实践证明:该方法不但能正确刻画矿井水文地质条件,而且还能将矿井涌水量预测和矿井生产进度紧密结合起来,具有较高的精度。 展开更多
关键词 工作面 涌水量 三维非稳定流数值模型
在线阅读 下载PDF
水电站非正常运行下库区涌浪数值模拟 被引量:3
18
作者 吴时强 王煌 《水利水运科学研究》 CSCD 1999年第4期321-328,共8页
水电站非正常运行情况下,引水系统流量突变往往会引起库区及下游河道涌浪传播,对于有航运要求的水利枢纽,这种涌浪往往会影响引航道附近的水流流场,恶化通航条件。本文建立了库区涌浪计算的平面二维水流数学模型,并应用于凌津滩水... 水电站非正常运行情况下,引水系统流量突变往往会引起库区及下游河道涌浪传播,对于有航运要求的水利枢纽,这种涌浪往往会影响引航道附近的水流流场,恶化通航条件。本文建立了库区涌浪计算的平面二维水流数学模型,并应用于凌津滩水电站和琅琊山抽水蓄能电站中,表明该模型具有较强的计算模拟能力。 展开更多
关键词 非恒定流 数值模拟 水库涌浪 引航道 水电站
在线阅读 下载PDF
成层土上畦灌非恒定流方程的求解 被引量:2
19
作者 刘才良 路振广 《河海大学学报(自然科学版)》 CAS CSCD 1993年第3期59-65,共7页
农业土壤受发育过程和长期耕作的影响,具有明显的分层性.考虑到耕作层与非耕作层在土壤物理性质上的差异.用Saint-Venant方程和改进后的Green-Ampt下渗模型联立,对成层土上的畦灌非恒定流方程进行数值模拟.计算结果与实测资料吻合较好,... 农业土壤受发育过程和长期耕作的影响,具有明显的分层性.考虑到耕作层与非耕作层在土壤物理性质上的差异.用Saint-Venant方程和改进后的Green-Ampt下渗模型联立,对成层土上的畦灌非恒定流方程进行数值模拟.计算结果与实测资料吻合较好,可在此基础上确定灌水技术参数的合理组合,为灌区实行分区水管理提供依据. 展开更多
关键词 成层土 畦灌 非恒定流 数值模拟
在线阅读 下载PDF
蒙洼蓄滞洪区洪水淹没数值模拟 被引量:8
20
作者 张红萍 胡昌伟 刘舒 《中国水利水电科学研究院学报》 2005年第4期291-296,共6页
蓄滞洪区是我国江河防洪体系中重要的组成部分。蒙洼作为淮河中游一个重要的蓄洪区,在淮河洪水调度中发挥着至关重要的作用。采用二维非恒定流不规则网格模型,模拟蒙洼蓄洪区内洪水淹没过程,分析蓄洪区内洪水风险,并采用2003年淮河大洪... 蓄滞洪区是我国江河防洪体系中重要的组成部分。蒙洼作为淮河中游一个重要的蓄洪区,在淮河洪水调度中发挥着至关重要的作用。采用二维非恒定流不规则网格模型,模拟蒙洼蓄洪区内洪水淹没过程,分析蓄洪区内洪水风险,并采用2003年淮河大洪水的实测资料,对模型进行验证。结果表明,该模型能较好地模拟蒙洼蓄洪区洪水淹没过程,可用于蒙洼蓄洪区的规划、建设、管理,乃至整个淮河流域的洪水预报。 展开更多
关键词 蒙洼蓄洪区 二维非恒定流 数值模拟 淮河
在线阅读 下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部