期刊文献+
共找到254篇文章
< 1 2 13 >
每页显示 20 50 100
甲烷水合物形成对蒙脱石表面区域电学特性的影响:分子动力学和有限元数值模拟
1
作者 武晓菲 邢兰昌 +3 位作者 魏伟 韩维峰 吕笑焱 钟杰 《新能源进展》 CSCD 北大核心 2024年第5期557-568,共12页
含天然气水合物泥质沉积物的电学特性是反演水合物饱和度、渗透率等储层参数的重要物理性质。甲烷水合物形成对黏土矿物颗粒表面电学性质的影响机制尚不明晰,制约了对骨架颗粒表面电导及沉积物电性参数的准确建模和计算。提出分子动力学... 含天然气水合物泥质沉积物的电学特性是反演水合物饱和度、渗透率等储层参数的重要物理性质。甲烷水合物形成对黏土矿物颗粒表面电学性质的影响机制尚不明晰,制约了对骨架颗粒表面电导及沉积物电性参数的准确建模和计算。提出分子动力学(MD)与有限元(FE)数值模拟相结合的新方法,基于Gouy-Chapman-Stern(GCS)理论建立无/含水合物-蒙脱石-Na Cl溶液体系的MD和FE数值模型,探讨甲烷水合物形成对蒙脱石表面双电层性质和表面区域电势分布的影响规律及机理。研究结果表明:(1)利用MD模型可获得含水合物条件下蒙脱石表面双电层中Stern层和扩散层的厚度以及溶液中离子的扩散系数,为FE模型提供双电层结构参数与离子运动参数,基于FE模型可分析蒙脱石表面区域电势分布特性;(2)蒙脱石表面形成稳定致密的钠离子层与界面水层,甲烷分子无法在蒙脱石表面生成水合物,因此水合物不影响Stern层而仅侵占部分扩散层区域;水合物形成过程中“排盐效应”使溶液的盐离子浓度增大,对扩散层中Na+产生排斥作用,使扩散层中的部分Na+进入Stern层,导致扩散层变薄;水合物的生长降低了溶液中离子的扩散系数;(3)水合物对扩散层的侵占使得扩散层部分区域的相对介电常数发生变化,引起该区域中电势的非均匀分布;水合物的存在不影响Stern层内的电势分布,但被水合物侵占的扩散层部分相对介电常数变小,引起扩散层内平均电势的升高;随着水合物与黏土颗粒表面距离的增加,平均电势逐渐降低。研究结果可为认识含水合物泥质沉积物中骨架颗粒的电学特性、量化颗粒表面电性参数提供理论和模型基础。 展开更多
关键词 甲烷水合物 黏土矿物 双电层 电学特性 分子动力学 有限元数值模拟
在线阅读 下载PDF
连续电性介质线源井-地电位三维有限元数值模拟 被引量:17
2
作者 刘海飞 陈德鹏 +2 位作者 戴前伟 柳建新 阮百尧 《桂林理工大学学报》 CAS 北大核心 2011年第1期28-38,共11页
为了更好的研究油藏开发过程中剩余油的分布情况以及注水前沿的推进位置,给出了三维空间任意线电流源且电导率连续变化的井-地电位三维有限元数值模拟方法。在模拟过程中,采用异常电位法消除线电流源奇异性的影响,并且将网格节点间的电... 为了更好的研究油藏开发过程中剩余油的分布情况以及注水前沿的推进位置,给出了三维空间任意线电流源且电导率连续变化的井-地电位三维有限元数值模拟方法。在模拟过程中,采用异常电位法消除线电流源奇异性的影响,并且将网格节点间的电导率设计为连续变化,使其能够模拟实际中更为复杂的地质情况;考虑到采油井或注水井可能存在垂直、倾斜及水平等任意情况,推导了任意线电流源正常电位的计算公式;采用对称超松弛预条件共轭梯度法解大型线性方程组,解决了井-地电位三维模拟速度慢的问题,并编制了实用、高效的井-地电位三维模拟软件。通过对3例地电模型进行模拟计算,分析了水泥环胶结质量、回流井到观测井的距离以及目标层的深度对井-地电位观测结果的影响。 展开更多
关键词 井-地电位 三维有限元数值模拟 异常电位 任意线电流源 预条件共轭梯度法
在线阅读 下载PDF
电子封装件受热载荷作用有限元数值模拟分析 被引量:13
3
作者 葛增杰 顾元宪 +1 位作者 王宏伟 靳永欣 《大连理工大学学报》 EI CAS CSCD 北大核心 2005年第3期320-325,共6页
针对典型塑料方形扁平封装体PQFP在工作过程中的受热分析问题,基于热弹性力学理论建立了二维和三维有限元数值模拟分析模型.研究了封装体在功率耗散情况下受均匀和非均匀热载作用时其材料的热膨胀和导热性质.有限元数值模拟取得了与实... 针对典型塑料方形扁平封装体PQFP在工作过程中的受热分析问题,基于热弹性力学理论建立了二维和三维有限元数值模拟分析模型.研究了封装体在功率耗散情况下受均匀和非均匀热载作用时其材料的热膨胀和导热性质.有限元数值模拟取得了与实验一致的结果.数值结果表明,采用较小弹性模量和热膨胀率的材料可以有效地减小热应力,基板和芯片的厚度是影响封装体变形的主要参数.数值分析结果为提高封装件的可靠性和优化设计提供了理论依据. 展开更多
关键词 数值模拟分析 载荷作用 电子封装 受热 三维有限元数值模拟 热弹性力学 分析问题 工作过程 数值结果 导热性质 功率耗散 分析模型 理论建立 热膨胀率 弹性模量 优化设计 分析结果 封装体 非均匀 热应力 可靠性 体变形 材料
在线阅读 下载PDF
低渗透气藏介质变形-渗流耦合有限元数值模拟 被引量:7
4
作者 孙峰 李行船 +1 位作者 薛世峰 仝兴华 《大庆石油地质与开发》 CAS CSCD 北大核心 2010年第6期101-104,共4页
低渗透气藏开发过程中介质变形对储层孔隙度、渗透率及产能的影响不可忽视。综合考虑储层变形、地应力、流体渗流等因素,建立了低渗透气藏介质变形一渗流耦合模型;采用有限元数值模拟技术并开发相应程序对耦合模型进行求解,定量研究... 低渗透气藏开发过程中介质变形对储层孔隙度、渗透率及产能的影响不可忽视。综合考虑储层变形、地应力、流体渗流等因素,建立了低渗透气藏介质变形一渗流耦合模型;采用有限元数值模拟技术并开发相应程序对耦合模型进行求解,定量研究开发过程中储层变形及其对储层物性和气藏产能的影响。数值模拟结果表明:非均匀水平地应力造成近井壁区域储层变形及渗透率变化各向异性,产能评价中应综合考虑介质变形的影响。 展开更多
关键词 低渗透气藏 介质变形 耦合模型 有限元数值模拟
在线阅读 下载PDF
基于有限元数值模拟的致密储层体积压裂效果影响参数分析 被引量:7
5
作者 冯福平 雷扬 +3 位作者 陈顶峰 胡超洋 王胡振 黄芮 《油气藏评价与开发》 CSCD 2019年第1期29-33,50,共6页
致密储层体积压裂作业后形成复杂裂缝网络,通过改变压裂井附近的渗流场来增大水平井产量,目前针对致密储层体积压裂影响参数对水平井投产后不同时间段压裂效果影响的权重分析较少。通过建立适用于致密油藏油水两相有限元方程对复杂裂缝... 致密储层体积压裂作业后形成复杂裂缝网络,通过改变压裂井附近的渗流场来增大水平井产量,目前针对致密储层体积压裂影响参数对水平井投产后不同时间段压裂效果影响的权重分析较少。通过建立适用于致密油藏油水两相有限元方程对复杂裂缝网络进行有限元数值模拟,采用正交实验与灰色关联法分析水平井体积压裂影响参数在油井投产10年内对产能的影响程度。研究结果表明,裂缝导流能力决定水平井压裂初期产能,裂缝长度决定水平井日产量的递减速度以及后期阶段稳产的产量水平,随着压裂井生产不断进行,裂缝长度取代裂缝导流能力成为影响产量的关键性因素。水平井投产不同时间段内各影响参数对压裂效果的影响程度研究可以为致密储层体积压裂水平井施工设计提供理论支持。 展开更多
关键词 致密油藏 影响参数 体积压裂 有限元数值模拟 复杂裂缝
在线阅读 下载PDF
护舷碰撞过程有限元数值模拟 被引量:8
6
作者 林红吉 孟宪林 +1 位作者 严建平 张海永 《船舶》 北大核心 2004年第5期47-49,共3页
本文基于ANSYS-LSDYNA,采用有限元方法模拟了钢板撞击橡胶泡沫型护舷的全过程,分析了钢板在撞击过程中速度改变量、不同时刻护舷的应变量以及护舷等效应力,结果表明有限元是一种有效的研究护舷碰撞过程的方法。
关键词 护舷 碰撞过程 有限元数值模拟 撞击 ANSYS 速度 全过程 改变量 等效应力 时刻
在线阅读 下载PDF
基于有限元数值模拟的宽厚板热轧过程宽度控制ANN模型 被引量:2
7
作者 阮金华 张立文 +2 位作者 顾森东 何文斌 陈少慧 《大连理工大学学报》 EI CAS CSCD 北大核心 2013年第5期678-684,共7页
利用商业有限元分析软件DEFORM-3D,建立宽厚板热轧过程中立辊侧压和随后平轧过程的三维热力耦合刚塑性有限元模型.模拟了不同板坯厚度、板坯宽度、立辊压下量和平辊压下率条件下320种立辊侧压和随后平轧过程,探讨了不同工艺参数对调宽... 利用商业有限元分析软件DEFORM-3D,建立宽厚板热轧过程中立辊侧压和随后平轧过程的三维热力耦合刚塑性有限元模型.模拟了不同板坯厚度、板坯宽度、立辊压下量和平辊压下率条件下320种立辊侧压和随后平轧过程,探讨了不同工艺参数对调宽效率的影响规律.模拟结果表明,调宽效率随着板坯宽度和立辊压下量的增加而升高,随着板坯厚度和平辊压下率的增加而降低.基于上述有限元数值模拟结果,借助BP人工神经网络,建立了热轧宽厚板立轧-平轧宽度控制模型,经测试,模型预测的调控效率与有限元数值模拟结果符合很好.将宽度控制模型用于宽厚板实际热轧过程的有限元数值模拟,模拟轧件宽度与实测宽度吻合很好. 展开更多
关键词 热轧宽厚板 立轧-平轧 宽度控制 有限元数值模拟 人工神经网络
在线阅读 下载PDF
混合配筋预应力混凝土密肋楼盖三维有限元数值模拟 被引量:2
8
作者 刘忠 周长顺 杨涛 《湘潭大学自然科学学报》 CAS CSCD 北大核心 2008年第4期60-64,共5页
采用有限元软件ANSYS,模拟了均布荷载作用下,周边简支混合配筋预应力混凝土密肋楼盖从开始加载到破坏的全过程,得到了极限荷载,荷载-挠度曲线和裂缝分布情况.分析表明,计算结果与试验结果吻合较好,从而验证了在对混合配筋预应力混凝土... 采用有限元软件ANSYS,模拟了均布荷载作用下,周边简支混合配筋预应力混凝土密肋楼盖从开始加载到破坏的全过程,得到了极限荷载,荷载-挠度曲线和裂缝分布情况.分析表明,计算结果与试验结果吻合较好,从而验证了在对混合配筋预应力混凝土密肋楼盖进行非线性有限元数值模拟中,应用实体单元SOLID65和等效荷载法的合理性. 展开更多
关键词 三维有限元数值模拟 等效荷载法 实体单元 极限荷载 荷载-挠度曲线
在线阅读 下载PDF
古构造应力场有限元数值模拟的应用及展望 被引量:26
9
作者 孙晓庆 《断块油气田》 CAS 2008年第3期31-33,共3页
古构造应力场有限元数值模拟技术已成为古构造应力场研究的一种有效方法,能够实现从动力学机制上定量计算整个含油气盆地古构造应力场的时空变化,对含油气盆地储集层特征描述及油气运聚分析具有重要的理论和实际意义。文中主要介绍了古... 古构造应力场有限元数值模拟技术已成为古构造应力场研究的一种有效方法,能够实现从动力学机制上定量计算整个含油气盆地古构造应力场的时空变化,对含油气盆地储集层特征描述及油气运聚分析具有重要的理论和实际意义。文中主要介绍了古构造应力场有限元数值模拟技术在含油气盆地研究中的应用现状,并对该技术的发展前景作了进一步展望。 展开更多
关键词 古构造应力场 有限元数值模拟 裂缝 油气运聚
在线阅读 下载PDF
有限元数值模拟方法在构造裂缝预测中的应用 被引量:19
10
作者 沈国华 《油气地质与采收率》 CAS CSCD 北大核心 2008年第4期24-26,29,共4页
对地下油气储层的裂缝预测是中外油气勘探和开发领域的技术难题。尽管目前有多种方法可以对地下裂缝的分布状况进行预测,但都存在不少问题。为此,提出采用有限元数值模拟方法对构造裂缝的发育状况进行预测。该方法首先是在地质规律指导... 对地下油气储层的裂缝预测是中外油气勘探和开发领域的技术难题。尽管目前有多种方法可以对地下裂缝的分布状况进行预测,但都存在不少问题。为此,提出采用有限元数值模拟方法对构造裂缝的发育状况进行预测。该方法首先是在地质规律指导下建立客观的油藏构造精细地质模型;然后将岩石力学实验与有限元数值模拟方法相结合,并将由声发射实验测算的古构造应力值,作为应力场模拟的绝对反演标准;最后模拟裂缝主要形成期的古构造应力场,依据对应力场的分析预测储层构造裂缝的发育程度。模拟结果表明,有限元数值模拟方法不仅能够恢复构造裂缝形成期储层构造应力场的分布,而且可以依据地质力学理论对储层构造裂缝的分布及发育程度进行预测,为裂缝性油藏的勘探和开发提供了技术支持。 展开更多
关键词 构造裂缝 有限元数值模拟 声发射实验 古构造应力场
在线阅读 下载PDF
金属切削过程有限元数值模拟 被引量:3
11
作者 李初晔 王焱 孟月梅 《航空制造技术》 2010年第22期28-33,共6页
本课题针对航空结构件加工中最关心的成型零件的变形状况、加工效率(金属去除率)、刀具磨损等方面,综合了近年来国内外在切削加工方面的进展,根据这些成果在现有CAE分析平台基础上开发的仿真程序,能够实现对复杂航空结构件的变形控制。
关键词 有限元数值模拟 切削过程 金属 航空结构件 加工效率 CAE分析 成型零件 刀具磨损
在线阅读 下载PDF
桶形基础采油平台循环稳定性三维有限元数值模拟 被引量:4
12
作者 李驰 王建华 +2 位作者 袁中立 秦延龙 刘振纹 《中国海洋平台》 2008年第1期26-30,共5页
对大港油田桶形基础采油平台稳定性进行研究。通过三维有限元数值模拟,计算平台基础的静承载力,并将循环荷载作用效果等效为一种与地基土初始静应力有关的地基土强度的变化,计算平台基础的循环承载力。进而,对平台基础在自身工作荷载以... 对大港油田桶形基础采油平台稳定性进行研究。通过三维有限元数值模拟,计算平台基础的静承载力,并将循环荷载作用效果等效为一种与地基土初始静应力有关的地基土强度的变化,计算平台基础的循环承载力。进而,对平台基础在自身工作荷载以及工作荷载和环境荷载共同作用下的稳定性进行评价。结果显示,该平台基础上的工作荷载仅为竖向静承载力的14%,按照平台基础所处的海洋环境,当不考虑海洋环境荷载的动力特性时,水平静承载力安全系数3.75;当考虑频率为0.1Hz规则的波浪荷载连续作用在平台基础上3个小时,水平循环承载力安全系数3.42。该平台基础偏于安全。 展开更多
关键词 海洋工程 桶形基础平台 三维有限元数值模拟 循环稳定性
在线阅读 下载PDF
有限元数值模拟预控地表沉降在大断面地铁隧道中的应用 被引量:4
13
作者 曹俊发 《铁道建筑》 北大核心 2005年第9期42-44,共3页
根据有限元数值模拟和优化分析,结合深圳地铁一期工程香—车区间人防段的地质情况和施工工艺进行了地表沉降预控演算,有效地控制了地表沉降、优化了施工方案并且验证了有限元数值模拟结果的可信性,为城市地下工程施工控制地表沉降提供... 根据有限元数值模拟和优化分析,结合深圳地铁一期工程香—车区间人防段的地质情况和施工工艺进行了地表沉降预控演算,有效地控制了地表沉降、优化了施工方案并且验证了有限元数值模拟结果的可信性,为城市地下工程施工控制地表沉降提供可借鉴的经验。 展开更多
关键词 有限元数值模拟 大断面隧道 预控 地表沉降
在线阅读 下载PDF
有限元数值模拟技术在西湖凹陷中央反转构造带形成机制研究中的应用 被引量:7
14
作者 张田 张建培 +2 位作者 张绍亮 余逸凡 唐贤君 《海洋石油》 CAS 2012年第4期11-16,共6页
在西湖凹陷中央反转构造带区域地质背景研究成果的基础上,基于有限元数值模拟技术,设计合理的有限元应力场模型和边界条件,模拟了西湖凹陷中央反转构造带的构造演化过程,探讨其成因模式。西湖凹陷中央反转构造带的构造演化过程可划分为... 在西湖凹陷中央反转构造带区域地质背景研究成果的基础上,基于有限元数值模拟技术,设计合理的有限元应力场模型和边界条件,模拟了西湖凹陷中央反转构造带的构造演化过程,探讨其成因模式。西湖凹陷中央反转构造带的构造演化过程可划分为3个阶段:(1)背斜带的初始形成阶段(始新世末,约32 Ma):在太平洋板块NWW-SEE向俯冲作用下,西湖凹陷的北部和中部地区分别形成了"A"背斜构造带和"B"背斜构造带两个不同高度的小背斜带;(2)背斜带的生长阶段(渐新世末,约23 Ma):在太平洋板块持续俯冲作用下,北部的"A"背斜构造带开始不断的向北扩展抬升,直至福江低凸起,背斜面积不断扩大。中部的"B"背斜构造带则不断的向南扩展抬升,当其扩展至西湖凹陷南部时,衍生出"C"背斜构造带和"D"背斜构造带两个小的背斜带;(3)背斜带的成型阶段(中新世末,约10 Ma):北部的"A"背斜构造带停止扩张,但垂向上仍继续抬升,抬升高度大于其它背斜构造带;中部和南部的"B"、"C"和"D"背斜构造带彼此连通,形成一个统一的大背斜构造带,由此东海陆架盆地西湖凹陷的中央反转带演化成型。 展开更多
关键词 西湖凹陷 中央反转构造带 形成机制 有限元数值模拟
在线阅读 下载PDF
有限元数值模拟在煤炭沟煤层开采中的应用 被引量:6
15
作者 姚普 肖利佳 《现代矿业》 CAS 2009年第7期59-62,共4页
对煤炭沟井田1-1′地质剖面煤层开采覆岩变形进行了有限元数值模拟分析。从岩体的本构关系入手,选取岩体参数,概化和建立工程地质模型,从而获得煤层开采中顶板岩体位移变化规律特征,分析结果表明,随着该区煤层的开采,会引发地面塌陷、... 对煤炭沟井田1-1′地质剖面煤层开采覆岩变形进行了有限元数值模拟分析。从岩体的本构关系入手,选取岩体参数,概化和建立工程地质模型,从而获得煤层开采中顶板岩体位移变化规律特征,分析结果表明,随着该区煤层的开采,会引发地面塌陷、斜坡滑坡和地裂缝等地质灾害,其结论对生产实践具有指导作用。 展开更多
关键词 有限元数值模拟 工程地质模型 地质灾害 岩体位移
在线阅读 下载PDF
基于有限元数值模拟的新下花桥静载试验分析 被引量:4
16
作者 粟进清 《公路与汽运》 2020年第5期103-105,116,共4页
对桥梁结构进行静载试验,是各类桥梁施工质量控制与评定中的重要环节。文中利用有限元软件建立新下花桥梁格有限元模型进行数值模拟,基于有限元分析结果,开展现场静载试验,通过将挠度和应变试验结果与理论结果的比较,对该桥的承载能力... 对桥梁结构进行静载试验,是各类桥梁施工质量控制与评定中的重要环节。文中利用有限元软件建立新下花桥梁格有限元模型进行数值模拟,基于有限元分析结果,开展现场静载试验,通过将挠度和应变试验结果与理论结果的比较,对该桥的承载能力和工作状况作出评价。 展开更多
关键词 桥梁 空心板桥 静载试验 有限元数值模拟
在线阅读 下载PDF
基于有限元数值模拟基坑支护结构稳定性分析 被引量:4
17
作者 范立康 《东北水利水电》 2019年第7期41-42,59,共3页
文中通过有限元数值模拟进行的固态建模分析,进而对所研究基坑工程的稳定性及应力、变形深入分析,对工程中出现的基坑边坡进行支护结构稳定性及敏感性分析计算.
关键词 有限元数值模拟 基坑 支护 稳定性
在线阅读 下载PDF
土压盾构穿越南水北调干渠的有限元数值模拟 被引量:1
18
作者 田均举 蔡松 +1 位作者 甄成 吴岱諠 《郑州铁路职业技术学院学报》 2020年第4期19-22,共4页
运用ABAQUS有限元数值模拟软件建立隧道下穿干渠三维模型,研究常水位工况和加大水位工况下土压盾构双线穿越所引起的干渠底部变形规律,并将模拟结果与现场监测数据进行对比分析。计算结果表明:在常水位工况下,地铁隧道双线通过后,南水... 运用ABAQUS有限元数值模拟软件建立隧道下穿干渠三维模型,研究常水位工况和加大水位工况下土压盾构双线穿越所引起的干渠底部变形规律,并将模拟结果与现场监测数据进行对比分析。计算结果表明:在常水位工况下,地铁隧道双线通过后,南水北调干渠结构最大沉降量为10.32 mm;在加大水位工况下,地铁隧道双线通过后,南水北调干渠结构最大沉降量为10.51 mm;最大沉降量位于隧道轴线处且以中轴线为对称轴近似呈正态分布;两种水位工况下最大沉降量相差不大,说明隧道下穿过程中干渠水位变化对干渠结构稳定性影响较小。 展开更多
关键词 有限元数值模拟 土压盾构 隧道下穿干渠 南水北调 现场监测
在线阅读 下载PDF
三维有限元数值模拟在采空区稳定性评价中的应用 被引量:2
19
作者 晁军 雷醒民 樊育豪 《工程建设与设计》 2018年第23期15-18,共4页
采空区已成为榆林某煤矿的严重安全隐患之一。针对该矿采空区数量多、体积大的特点,在现场调查和室内力学实验研究的基础上,利用计算机数值模拟技术,采用3DEC软件对采空区围岩的破坏过程进行模拟计算。结果显示,目前该矿采空区稳定性良... 采空区已成为榆林某煤矿的严重安全隐患之一。针对该矿采空区数量多、体积大的特点,在现场调查和室内力学实验研究的基础上,利用计算机数值模拟技术,采用3DEC软件对采空区围岩的破坏过程进行模拟计算。结果显示,目前该矿采空区稳定性良好;地表作用200kPa荷载情况下,由巷道引起的基岩顶部最大位移为12mm,由巷道造成的地表倾斜度不会超过2‰,由此判定巷道对地表的影响较小;三维有限元数值模拟的方法用于采空区稳定性评价,对矿山的安全生产具有重要的指导作用。 展开更多
关键词 三维有限元数值模拟 采空区 稳定性评价
在线阅读 下载PDF
钢结构螺栓连接的有限元数值模拟
20
作者 王智 《低温建筑技术》 2013年第2期61-63,共3页
螺栓连接是钢结构连接的主要方式之一,具有施工简单、拆装方便等优点。现行钢结构规范中对螺栓连接的计算都是在一定的假设条件下进行的。本文以钢结构规范中的计算公式为基础,运用有限元数值分析的方法,采用ABAQUS来分析普通螺栓群在... 螺栓连接是钢结构连接的主要方式之一,具有施工简单、拆装方便等优点。现行钢结构规范中对螺栓连接的计算都是在一定的假设条件下进行的。本文以钢结构规范中的计算公式为基础,运用有限元数值分析的方法,采用ABAQUS来分析普通螺栓群在受剪和受拉条件下的各个螺栓的受力情况。将有限元软件模拟的结果与公式计算结果进行比较,证明规范公式的正确性。文中分别模拟分析了普通螺栓群在受剪的情况下承受轴心力、扭矩,和普通螺栓群在受拉的情况下承受拉力、弯矩的受力情况。 展开更多
关键词 螺栓连接 有限元数值模拟
在线阅读 下载PDF
上一页 1 2 13 下一页 到第
使用帮助 返回顶部