期刊文献+
共找到21篇文章
< 1 2 >
每页显示 20 50 100
DEVELOPMENT OF COMPUTER SIMULATION SOFTWAREON PHASE TRANSFORMATION DURING QUENCHING 被引量:5
1
作者 L. Chen, W. Xia and X. J. Chen Department of Materials Science and Engineering, Nanchang Institute of Aeronautical Technology, Nanchang 330034, China 《Acta Metallurgica Sinica(English Letters)》 SCIE EI CAS CSCD 2000年第2期618-622,共5页
A computer software to simulate the phase transformation during quenching is designed based on Avrami equation and Scheil additivity principle of incubation period. The isothermal transformation diagrams of supercool... A computer software to simulate the phase transformation during quenching is designed based on Avrami equation and Scheil additivity principle of incubation period. The isothermal transformation diagrams of supercooled austenite are described by cubic spline functions. This software is possess of a good graphic interface of Windows style, can simulate the whole process in austenite decomposition during continuous cooling. If the cooling rate was given, the fraction of various microstructures transformed of austenite decomposition during continuous cooling at any temperature can be calculated. The simulation results are checked with the quenching experiment of 45 steel. The results indicate that the simulation results are comparatively close to the experimental results. 展开更多
关键词 quenching phase transformation computer simulation software
在线阅读 下载PDF
COUPLED NUMERICAL SIMULATION ON COLD ROLLER'S TEMPERATURE FIFLD-PHASE TRANSFORMATION - STRESS FIELD DURING ITS QUENCHING PSOCESS 被引量:15
2
作者 J. F. Gu J. S. Pan M. J. Hu and F. F. Shen (School of Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai 200030, China) 《Acta Metallurgica Sinica(English Letters)》 SCIE EI CAS CSCD 2000年第1期254-262,共9页
Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short ti... Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short time, and these variables take a contact interactions in the whole process. In this paper, a three dimensional non - linear mathematical model for queeching process has been founded and the numerical simulation on temperature field,microstructre and stress field has been realized.In the FEM analysis, the incremental iteration method is used to deal with such complicated nonlinear as boundary nonlinear, physical property nonlinear,transformation nonlinear etc.The effect of stress on transformation kinetics has been considered in the calculation of microstructure. In the stress field anal- ysis,a thermo- elasto - plastic model has been founded, which considers such factors as transforma- tion strain,transformation plastic strain, themal strain and the effect of temperature and transforma- tion on mechanical propertier etc. The transient temperature field, microstructure distribution and stress field of the roller on any time can be displayed vividly,and the cooling curve and the changes of stress on any position can also be given. 展开更多
关键词 quenching phase transformation temperature field stress field finite element method (FEM) numerical simulation
在线阅读 下载PDF
大型锻件喷雾冷却过程数值模拟 被引量:19
3
作者 杜凤山 张芳 +1 位作者 黄华贵 李广睿 《金属热处理》 CAS CSCD 北大核心 2008年第5期35-38,共4页
为研究锻件喷雾淬火冷却过程,研制了喷雾淬火热处理装置,并对轴类件进行喷雾淬火热处理试验。根据实际情况建立包含变物性、连续过冷沸腾换热边界的温度、组织、应力三场耦合非线性有限元模型,对喷雾冷却过程进行了模拟分析计算,温度计... 为研究锻件喷雾淬火冷却过程,研制了喷雾淬火热处理装置,并对轴类件进行喷雾淬火热处理试验。根据实际情况建立包含变物性、连续过冷沸腾换热边界的温度、组织、应力三场耦合非线性有限元模型,对喷雾冷却过程进行了模拟分析计算,温度计算结果与测试结果吻合较好,并分析了大型非淬透件热处理应力的尺寸效应。 展开更多
关键词 喷雾淬火 数值模拟 温度场 相变
在线阅读 下载PDF
淬火过程的计算机模拟及其应用 被引量:35
4
作者 顾剑锋 潘键生 胡明娟 《金属热处理》 EI CAS CSCD 北大核心 2000年第5期35-37,44,共4页
籍助于温度场 -相变 -应力和应变的有限元分析方法 (FEM) ,实现了复杂形状工件淬火过程的计算机模拟。模拟结果直观地显示任一时刻工件内部的温度分布、组织分布、应力分布、应变状况及工件任意位置上的冷却曲线 ,为使用者提供了用虚拟... 籍助于温度场 -相变 -应力和应变的有限元分析方法 (FEM) ,实现了复杂形状工件淬火过程的计算机模拟。模拟结果直观地显示任一时刻工件内部的温度分布、组织分布、应力分布、应变状况及工件任意位置上的冷却曲线 ,为使用者提供了用虚拟现实方法详尽分析和研究淬火工艺的工具。本文列举了几个实际应用的例子 。 展开更多
关键词 淬火 温度场 有限元法 计算机模拟 虚拟制造
在线阅读 下载PDF
1045钢淬火时温度场的数值模拟 被引量:13
5
作者 谢建斌 程赫明 +1 位作者 何天淳 文宏光 《甘肃工业大学学报》 北大核心 2003年第4期33-37,共5页
基于实验测定的1045钢淬火冷却曲线,应用有限差分原理和非线性估计法对非线性导热方程的逆问题进行了求解,给出了一种求解钢淬火时非线性表面换热系数的方法;应用数学转换方法计算了1045钢在连续冷却时奥氏体、珠光体、贝氏体和马氏体... 基于实验测定的1045钢淬火冷却曲线,应用有限差分原理和非线性估计法对非线性导热方程的逆问题进行了求解,给出了一种求解钢淬火时非线性表面换热系数的方法;应用数学转换方法计算了1045钢在连续冷却时奥氏体、珠光体、贝氏体和马氏体的体积百分数.应用有限元分析软件ANSYS计算了1045钢淬火时具有相变和非线性表面换热系数的温度场.研究表明,钢淬火时的温度场、导热与换热是非线性的,温度场的计算结果和实验结果比较吻合. 展开更多
关键词 淬火工艺 温度场 数值模拟 换热系数 有限元分析
在线阅读 下载PDF
冷轧辊淬冷过程数值模拟的研究 被引量:18
6
作者 顾剑锋 潘健生 +1 位作者 胡明娟 沈甫法 《金属热处理学报》 CSCD 1999年第2期1-7,共7页
采用三维非线性有限元分析方法模拟小型冷轧辊淬火冷却过程,直观地显示出任一时刻轧辊上的瞬态温度场、组织分布和应力场,可以给出轧辊上任意位置的冷却曲线。在所采用的数学模型中,用增量叠代法处理边界条件非线性、物理参数非线性... 采用三维非线性有限元分析方法模拟小型冷轧辊淬火冷却过程,直观地显示出任一时刻轧辊上的瞬态温度场、组织分布和应力场,可以给出轧辊上任意位置的冷却曲线。在所采用的数学模型中,用增量叠代法处理边界条件非线性、物理参数非线性、相变潜热非线性等复杂问题。在相变量计算中引入应力状态的影响。在应力场分析中采用热弹塑性模型,考虑了相变应变、相变塑性应变、热应变、材料力学性能的温度效应和相变影响等因素。计算机模拟结果和实测结果吻合较好。由此提供了一种良好的淬火工艺虚拟生产试验手段。 展开更多
关键词 淬火 相变 温度场 应力场 计算机模拟 冷轧辊
在线阅读 下载PDF
38MnVS6非调质钢棒材轧后冷却过程中的组织转变 被引量:2
7
作者 张驰 张立文 +2 位作者 顾森东 申文飞 夏英男 《金属热处理》 CAS CSCD 北大核心 2017年第2期207-210,共4页
以MSC.Marc有限元分析软件为平台,利用其二次开发功能,结合38Mn VS6非调质钢等温转变曲线,建立了38Mn VS6非调质钢棒材轧后冷却过程的有限元模型,实现了38Mn VS6非调质钢棒材轧后冷却过程温度和相变的耦合计算。通过模拟,得到了棒材轧... 以MSC.Marc有限元分析软件为平台,利用其二次开发功能,结合38Mn VS6非调质钢等温转变曲线,建立了38Mn VS6非调质钢棒材轧后冷却过程的有限元模型,实现了38Mn VS6非调质钢棒材轧后冷却过程温度和相变的耦合计算。通过模拟,得到了棒材轧后冷却过程温度和组织的分布和演变情况。计算得到的组织转变情况与实测结果吻合较好,这对优化38Mn VS6非调质钢棒材轧后冷却工艺参数有一定指导意义。 展开更多
关键词 38MnVS6非调质钢 轧后冷却 数值模拟 相转变
原文传递
预硬塑胶模具钢端部淬火工艺有限元模拟 被引量:1
8
作者 李辉平 赵国群 +1 位作者 王广春 栾贻国 《山东大学学报(工学版)》 CAS 2004年第1期5-8,共4页
介绍了对淬火工艺组织场进行计算的数学模型 ,用Scheil叠加方法解决了TTT曲线中孕育率、虚拟时间、扩散型相变量的计算 ,利用杠杆定律及铁碳相图解决了P2 0钢的先共析铁素体的计算 .编写了淬火工艺温度场、组织场模拟软件 ,利用模拟软件... 介绍了对淬火工艺组织场进行计算的数学模型 ,用Scheil叠加方法解决了TTT曲线中孕育率、虚拟时间、扩散型相变量的计算 ,利用杠杆定律及铁碳相图解决了P2 0钢的先共析铁素体的计算 .编写了淬火工艺温度场、组织场模拟软件 ,利用模拟软件对P2 0塑胶模具钢的端部淬火试验进行了模拟 。 展开更多
关键词 淬火 相变 有限元法 模拟
在线阅读 下载PDF
高压气淬过程的数值模拟 被引量:1
9
作者 王婧 陈乃录 +1 位作者 戴华 张伟民 《金属热处理》 EI CAS CSCD 北大核心 2007年第8期72-76,共5页
应用有限体积法建立了数学模型,对高压气淬过程中气体的速度场、气体和工件的温度场以及工件内部的组织场进行了耦合计算,针对一组小棒料高压气淬的模拟结果,分析了气体流动、工件温度与相变各自的影响因素以及它们之间的相互影响,指出... 应用有限体积法建立了数学模型,对高压气淬过程中气体的速度场、气体和工件的温度场以及工件内部的组织场进行了耦合计算,针对一组小棒料高压气淬的模拟结果,分析了气体流动、工件温度与相变各自的影响因素以及它们之间的相互影响,指出几何模型准确建立的重要性。另外通过模拟结果与实测冷却曲线的比较,该模型基本反映了高压气淬过程,其中加入相变计算能更符合实际,并且还能预测出工件的最终组织。应用本文中的相变计算模型,相变的计算结果与工件金相照片中的实测结果基本吻合。 展开更多
关键词 高压气淬 数值模拟 计算流体动力学 相变
在线阅读 下载PDF
Q345E钢焊接热影响区相变的计算机模拟 被引量:3
10
作者 王薇 金成 史春元 《兵器材料科学与工程》 CAS CSCD 北大核心 2015年第4期80-84,共5页
以Kirkaldy和Koistinen-Marburger相变模型为基础,建立焊接热影响区相变模型。以Q345E钢平板对接焊为例,用Goldak双椭球功率密度分布热源模型计算焊接热循环,用ABAQUS有限元分析软件对焊接过程中的热影响区的相变进行数值模拟。以相同... 以Kirkaldy和Koistinen-Marburger相变模型为基础,建立焊接热影响区相变模型。以Q345E钢平板对接焊为例,用Goldak双椭球功率密度分布热源模型计算焊接热循环,用ABAQUS有限元分析软件对焊接过程中的热影响区的相变进行数值模拟。以相同焊接规范进行焊接试验,利用Image-pro-plus软件进行组织测定。将试验结果与模拟结果进行比较分析,二者吻合良好,证明本焊热影响区相变模型可以准确分析焊接热影响区相变情况,对优化焊接工艺参数、改善焊接接头性能、提高焊接产品质量具有重要指导作用。 展开更多
关键词 双椭球功率密度分布热源模型 相变模型 有限元软件 数值模拟
在线阅读 下载PDF
热作模具钢锻后组织模拟及软件设计 被引量:1
11
作者 杨玲玲 邓子玉 曾乐 《锻压技术》 CAS CSCD 北大核心 2007年第1期6-9,共4页
联合使用MC和CA方法,结合金属学的基本原理,模拟了5CrNiMo热作模具钢锻后微观组织演变;并采用MATLAB为开发平台,结合C语言编程完成模拟软件的开发。软件可以实现5CrNiMo钢锻后冷却过程中发生的奥氏体晶粒长大及奥氏体至铁素体相变过程... 联合使用MC和CA方法,结合金属学的基本原理,模拟了5CrNiMo热作模具钢锻后微观组织演变;并采用MATLAB为开发平台,结合C语言编程完成模拟软件的开发。软件可以实现5CrNiMo钢锻后冷却过程中发生的奥氏体晶粒长大及奥氏体至铁素体相变过程的动态演示,模拟结果直观,并可以定量计算晶粒尺寸的变化,从而实现对锻后冷却过程中不同温度下的组织形态和晶粒尺寸的预测。此外,通过实验,验证了模拟软件的正确性和可靠性。 展开更多
关键词 热作模具钢 组织模拟 软件设计 晶粒长大 相变
在线阅读 下载PDF
40Cr钢连续冷却相变的数值模拟 被引量:2
12
作者 杨书申 《热加工工艺》 CSCD 北大核心 2001年第4期51-52,共2页
淬火冷却过程中的相变对瞬态温度场和应力有较大影响。本文以 40 Cr钢为研究对象 ,用多项式拟合了等温转变图 ,采用孕期叠加法确定冷却过程相变的开始时间 ,给出了连续冷却相变动力学的数学模拟计算式 ,计算了淬火过程中的相组织的百分比。
关键词 淬火 相变 数值模拟 等温转变图 40r钢 计算机模拟
在线阅读 下载PDF
热处理模拟软件Thermal Prophet的开发与应用 被引量:16
13
作者 刘晓晖 韩利战 +1 位作者 顾剑锋 潘健生 《金属热处理》 CAS CSCD 北大核心 2012年第2期115-118,共4页
Thermal Prophet是最近开发的一款具有自主知识产权的热处理数值模拟软件。本文介绍了其基本架构、技术特点和基本功能。以热作模具钢模块淬火为例应用本软件进行了模拟,得到的畸变趋势与实际测量结果一致。通过考虑和忽略相变塑性两种... Thermal Prophet是最近开发的一款具有自主知识产权的热处理数值模拟软件。本文介绍了其基本架构、技术特点和基本功能。以热作模具钢模块淬火为例应用本软件进行了模拟,得到的畸变趋势与实际测量结果一致。通过考虑和忽略相变塑性两种模拟结果的对比,分析了相变塑性对淬火畸变所起的作用,指出了相变塑性是模拟淬火过程不可忽略的重要因素。 展开更多
关键词 THERMAL PROPHET 热处理数值模拟软件 有限元方法 相变塑性 淬火畸变
原文传递
基于扫描振镜的激光淬火温度特性模拟与实验研究 被引量:4
14
作者 陈凯烨 王晔 +3 位作者 吴国龙 王淼 张群莉 姚建华 《表面技术》 EI CAS CSCD 北大核心 2020年第5期251-258,共8页
目的与传统激光淬火(CLQ)进行对比研究,获得使用飞行光路的振镜扫描式激光淬火(GLQ)的温度变化特性和规律。方法使用两种淬火方式对45钢进行表面处理,将组织形貌和硬度分布进行比较,并结合数值模拟,研究两种激光淬火方式的温度变化特性... 目的与传统激光淬火(CLQ)进行对比研究,获得使用飞行光路的振镜扫描式激光淬火(GLQ)的温度变化特性和规律。方法使用两种淬火方式对45钢进行表面处理,将组织形貌和硬度分布进行比较,并结合数值模拟,研究两种激光淬火方式的温度变化特性和规律。结果在相同的总能量输入、扫描区域和扫描时间下,两种淬火方式硬化层横截面形貌呈现巨大差异:CLQ试样硬化层深度为753.66μm,宽度为3787.21μm,横截面形貌为"月牙形";GLQ试样硬化层深度为256.61μm,宽度为5808.77μm,形貌为"均匀性更好的近似月牙形"。截面硬度分布相似:均为中间主体区域高(GLQ试样硬度810~900HV,CLQ试样硬度790~830HV),两侧边缘区域低(均约为760HV)。模拟结果表明,CLQ试样特征点激光辐照时长共1.25 s,温度从508℃变为1063℃,奥氏体相变时间为1.17 s;GLQ试样特征点激光辐照时长共0.628 s,温度从500℃变成718℃,其波峰温度奥氏体相变时间为0.38 s,波谷温度均保持在马氏体转变温度以上。沿激光扫描方向,GLQ的温度变化率比CLQ大1个数量级;垂直激光扫描方向,GLQ温度变化率比CLQ大2个数量级。结论GLQ具有温度变化率更大、升降温多次循环、热处理相变时间更短、热积累更小等特点,有望在薄层大面积淬火上获得应用。 展开更多
关键词 激光淬火 振镜扫描 温度场 数值模拟 相变 45钢
在线阅读 下载PDF
重轨端部淬火温度场与相变的计算机模拟 被引量:3
15
作者 张敏 赵刚 +4 位作者 胡瑞海 李文凯 陈先舟 楚双学 董茂松 《特殊钢》 北大核心 2010年第5期14-17,共4页
运用有限元法对U71Mn钢(%:0.65~0.76C、0.15~0.35Si、1.10~1.40Mn、≤0.030V、≤0.025Ti)重轨端部电磁感应加热55 s,空冷8 s,喷风冷却42 s,空冷至常温的热处理过程进行了计算机模拟分析,计算机图形学理论动态显示轨端内部任意时刻的... 运用有限元法对U71Mn钢(%:0.65~0.76C、0.15~0.35Si、1.10~1.40Mn、≤0.030V、≤0.025Ti)重轨端部电磁感应加热55 s,空冷8 s,喷风冷却42 s,空冷至常温的热处理过程进行了计算机模拟分析,计算机图形学理论动态显示轨端内部任意时刻的温度分布、组织分布。模拟结果表明,在该冷却条件下淬火后端部横断面组织为细片状珠光体(索氏体+托氏体),无马氏体、贝氏体等组织,与实验结果相吻合。 展开更多
关键词 重轨 端部淬火 温度场 相变 计算机模拟
在线阅读 下载PDF
计算热力学及其在材料科学中的应用(英文) 被引量:4
16
作者 Bo Sundman Matthias Stratmann +1 位作者 张利军 杜勇 《中国材料进展》 CAS CSCD 北大核心 2015年第1期15-29,63,共16页
简要综述了材料科学领域中计算热力学与动力学的模型、数据库、软件及应用。介绍了当前最常用的软件。详细解释了用于描述不同相热力学性质的重要模型,以及获取热力学数据库中模型参数的优化过程。此外,还阐述了平衡计算中最常用的一种... 简要综述了材料科学领域中计算热力学与动力学的模型、数据库、软件及应用。介绍了当前最常用的软件。详细解释了用于描述不同相热力学性质的重要模型,以及获取热力学数据库中模型参数的优化过程。此外,还阐述了平衡计算中最常用的一种算法。最后本文给出了一些使用软件和数据库计算相平衡和相图,以及结合动力学数据用于相变及其过程模拟的示例。 展开更多
关键词 计算热力学 动力学模拟 相变 软件 数据库
在线阅读 下载PDF
超高强度钢板热冲压成形同步淬火相变不均匀性影响因素分析 被引量:9
17
作者 刘雪飞 贾卓翰 +3 位作者 张谦 张震 曾凡晶 权国政 《塑性工程学报》 CAS CSCD 北大核心 2019年第2期245-252,共8页
以超高强度钢BR1500HS为研究对象,利用材料的应力-应变数据、热物理性参数、相变模型构建了U形件热冲压成形过程的热-力-相3场动态耦合有限元数值模型。有限元模拟结果表明,法兰区域和底部区域中的马氏体体积含量显著多于侧壁区域中的... 以超高强度钢BR1500HS为研究对象,利用材料的应力-应变数据、热物理性参数、相变模型构建了U形件热冲压成形过程的热-力-相3场动态耦合有限元数值模型。有限元模拟结果表明,法兰区域和底部区域中的马氏体体积含量显著多于侧壁区域中的马氏体体积含量,揭示了钢板相变的不均匀性;同时对成形过程中板料和模具的温度场进行研究,揭示了钢板冷却速度的不均匀性。通过在整个热冲压过程中法兰区域、圆角区域、侧壁区域和底部区域选取观测点,分析了冷却速率不均匀性的影响因素,即温差、接触间隙、换热系数。最后通过实验验证了模拟结果具有较高的可靠性。 展开更多
关键词 超高强度钢板 热冲压成形 同步淬火 冷却速度 相变不均匀性 模拟仿真
在线阅读 下载PDF
一种非调质钢组织转变及性能的研究 被引量:6
18
作者 杨伟光 刘年富 +1 位作者 李祥龙 李国葵 《热处理技术与装备》 2019年第1期27-32,共6页
一种非调质钢在Gleeble 3800热模拟机上进行试验,测定了该非调质钢在不同冷却速度下膨胀曲线,同时结合金相-硬度法获得了该非调质钢的连续冷却转变时静态CCT曲线、显微组织和维氏硬度,以及等温转变时显微组织和维氏硬度。奥氏体在连续... 一种非调质钢在Gleeble 3800热模拟机上进行试验,测定了该非调质钢在不同冷却速度下膨胀曲线,同时结合金相-硬度法获得了该非调质钢的连续冷却转变时静态CCT曲线、显微组织和维氏硬度,以及等温转变时显微组织和维氏硬度。奥氏体在连续冷却速度分别为0.1、0.5、1℃/s以及在560℃等温转变前后的连续冷却速度为2℃/s,等温时间分别为10、20 min时,转变组织为铁素体+珠光体;奥氏体在560℃等温转变前后的连续冷却速度分别为3、4℃/s,等温时间分别为10、20 min时,转变组织为铁素体+珠光体+微量或少量贝氏体。本研究结果,对下游用户使用圆钢锻造成零件后,获得铁素体+珠光体组织或铁素体+珠光体+微量或少量贝氏体组织的冷却工艺提供了参考。 展开更多
关键词 非调质钢 热模拟 相转变 CCT曲线 等温转变
在线阅读 下载PDF
淬透型圆柱件热处理残余应力的数值模拟 被引量:3
19
作者 张雪姣 宋肖阳 朱琳 《一重技术》 2017年第6期39-44,共6页
对不同尺寸的小型圆柱件进行淬火过程的数值模拟,通过测量试验件的淬火残余应力验证模型的正确性。在此基础上,研究冷却速度、相变塑性、材料马氏体转变温度对淬火残余应力的影响,对淬火过程中工件内沿半径方向各点的残余应力随时间的... 对不同尺寸的小型圆柱件进行淬火过程的数值模拟,通过测量试验件的淬火残余应力验证模型的正确性。在此基础上,研究冷却速度、相变塑性、材料马氏体转变温度对淬火残余应力的影响,对淬火过程中工件内沿半径方向各点的残余应力随时间的变化规律及应力在横截面内的分布随时间的变化规律进行分析,讨论并总结形成均温淬透型圆柱件淬火残余应力的原因及规律。 展开更多
关键词 数值模拟 淬火 残余应力 相变塑性
在线阅读 下载PDF
全淬透圆柱件淬火应力的有限元模拟及实验验证 被引量:16
20
作者 刘玉 秦盛伟 +2 位作者 左训伟 陈乃录 戎咏华 《金属学报》 SCIE EI CAS CSCD 北大核心 2017年第6期733-742,共10页
利用温度场、组织场和应力场相互耦合的数学模型对直径60 mm的40CrNiMo水淬圆棒的淬火应力进行研究。结果表明,有限元模拟的淬火应力及其分布与XRD测量结果相符,论证了所运用的温度场、组织场和应力场相互耦合的数学模型(包括建立的相... 利用温度场、组织场和应力场相互耦合的数学模型对直径60 mm的40CrNiMo水淬圆棒的淬火应力进行研究。结果表明,有限元模拟的淬火应力及其分布与XRD测量结果相符,论证了所运用的温度场、组织场和应力场相互耦合的数学模型(包括建立的相变塑性函数)的正确性。通过计算机模拟分离了淬火马氏体钢中的热应力和组织应力,揭示了不同直径淬透圆棒试样的应力分布规律及其起因以及不同淬火介质对淬火应力的影响规律。 展开更多
关键词 全淬透圆柱件 淬火应力 热应力 组织应力 有限元模拟
原文传递
上一页 1 2 下一页 到第
使用帮助 返回顶部