期刊文献+
共找到75篇文章
< 1 2 4 >
每页显示 20 50 100
高饱和度含气土动力问题的有限体积法分步解法
1
作者 邹诗环 黄茂松 时振昊 《岩土力学》 北大核心 2025年第4期1286-1293,共8页
我国东南沿海海域广泛分布含封闭小气泡的海洋软土,小气泡的存在将对其压缩性产生很大影响。有限体积法有助于建立与流体模拟一致的海床动力响应计算方法,进而实现波浪-海床相互作用的模拟。矢量分解方法能显著提高分步解法的稳定时间步... 我国东南沿海海域广泛分布含封闭小气泡的海洋软土,小气泡的存在将对其压缩性产生很大影响。有限体积法有助于建立与流体模拟一致的海床动力响应计算方法,进而实现波浪-海床相互作用的模拟。矢量分解方法能显著提高分步解法的稳定时间步长,且该解法易于在开源有限体积法平台OpenFOAM上实现,故其具有高效与通用的特点。算例分析表明:所提出的分步求解算法能合理反映高饱和度含气土的动力响应;小气泡的存在将极大降低高饱和度含气土的瞬时孔压响应与孔压幅值,且考虑重力作用时,饱和度的下降会导致动孔压相位滞后。 展开更多
关键词 饱和度含气土 矢量分解方 分步解 有限体积 OPENFOAM
在线阅读 下载PDF
基于非结构网格有限体积法的二维高精度物质输运模拟 被引量:13
2
作者 王志力 陆永军 耿艳芬 《水科学进展》 EI CAS CSCD 北大核心 2008年第4期531-536,共6页
在二维任意多边形无结构网格有限体积水流模型的基础上,通过无插值的通量重构和引入限制因子建立物质输运的高精度无数值震荡数学模型,并从理论上分析了模型的单调和稳定条件。通过旋转流中圆形分布的对流输运算例对模型进行检验,结果... 在二维任意多边形无结构网格有限体积水流模型的基础上,通过无插值的通量重构和引入限制因子建立物质输运的高精度无数值震荡数学模型,并从理论上分析了模型的单调和稳定条件。通过旋转流中圆形分布的对流输运算例对模型进行检验,结果表明模型具有守恒、高精度无数值震荡等优点。最后将模型应用到瓯江河口的盐度输运模拟,计算的潮位、流速和盐度过程与原型吻合良好。 展开更多
关键词 非结构网格 有限体积 精度格式 物质输运
在线阅读 下载PDF
基于有限体积法碱金属高温热管冷态启动流动换热数值研究
3
作者 赵昊城 张泽秦 +3 位作者 王成龙 秋穗正 田文喜 苏光辉 《原子能科学技术》 EI CSCD 北大核心 2024年第1期50-59,共10页
为建立碱金属高温热管启动瞬态和运行稳态工作特性预测方法,本研究采用有限体积法(FVM)建立管壁导热模型、吸液芯流动传热模型和蒸气区模型。基于C语言开发并验证了碱金属高温热管的冷态启动瞬态分析程序,最大相对偏差为9.8%。仿真模拟... 为建立碱金属高温热管启动瞬态和运行稳态工作特性预测方法,本研究采用有限体积法(FVM)建立管壁导热模型、吸液芯流动传热模型和蒸气区模型。基于C语言开发并验证了碱金属高温热管的冷态启动瞬态分析程序,最大相对偏差为9.8%。仿真模拟了单根水平钠热管启动瞬态并开展敏感性分析,结果显示:对于本研究中使用的热管,在固定输入功率为1 000 W的环境条件下,启动开始后700 s热管蒸气区完全进入连续流态,到达稳态总用时为3 000 s,启动过程中工质熔化阶段吸液芯内部压力相对值逐渐增大,熔化完成后压力相对值略有降低;稳态运行下热管等温性良好,外壁面轴向温差稳定在22.5 K,吸液芯内部压降约为47 Pa;环境温度升高会延长热管到达稳态所需时间,并对稳态蒸气压力和流速分布产生一定影响;绝热段长度增加同样延长了热管启动到达稳态时间,同时对吸液芯内流体压力和速度分布存在一定影响。 展开更多
关键词 温热管 数值模拟 冷态启动 有限体积
在线阅读 下载PDF
高精度有限体积法与间断有限元法的比较 被引量:3
4
作者 范进之 李桦 《国防科技大学学报》 EI CAS CSCD 北大核心 2014年第5期33-38,共6页
通过数值算例,比较了高精度有限体积法和间断有限元法在求解不同问题时的表现。研究发现:在精度相同的条件下,间断有限元法的计算误差要明显小于有限体积法;间断有限元法的重构过程与高精度有限体积法相比较为简单,但高阶情形下解多项... 通过数值算例,比较了高精度有限体积法和间断有限元法在求解不同问题时的表现。研究发现:在精度相同的条件下,间断有限元法的计算误差要明显小于有限体积法;间断有限元法的重构过程与高精度有限体积法相比较为简单,但高阶情形下解多项式的自由度较多并且需要计算体积分,因此整个求解时间较长。降低时间积分时解多项式的自由度数目是实现高精度算法在实际问题中应用的重要手段。 展开更多
关键词 精度格式 有限体积 间断有限
在线阅读 下载PDF
双曲型守恒律的高阶、高分辨有限体积法 被引量:23
5
作者 李荫藩 宋松和 周铁 《力学进展》 EI CSCD 北大核心 2001年第2期245-263,共19页
有限体积法是一种离散积分形式守恒律的数值方法.它可以吸收有限元法和有限差分法的一些重要思想与技巧.由于它可方便地利用多种类型的网格(结构网格和非结构网格),从而非常适用于处理复杂计算区域,目前已成为一种在计算流体力学... 有限体积法是一种离散积分形式守恒律的数值方法.它可以吸收有限元法和有限差分法的一些重要思想与技巧.由于它可方便地利用多种类型的网格(结构网格和非结构网格),从而非常适用于处理复杂计算区域,目前已成为一种在计算流体力学中十分重要的方法.本文将针对二维双曲守衡律,对高精度、高分辨的有限体积法及其近年来的进展做一简要介绍. 展开更多
关键词 双曲守恒律 有限体积 结构网络 无结构网格 精度 分辨格式 计算流体力学
在线阅读 下载PDF
有限单元法在高拱坝坝肩动力稳定分析中的应用 被引量:14
6
作者 苏超 李俊宏 任青文 《河海大学学报(自然科学版)》 CAS CSCD 北大核心 2003年第2期144-147,共4页
利用有限元法及三维实体任意截面数据处理技术,建立了高拱坝坝肩稳定安全系数计算分析方法,并编制了相应的计算和处理程序.通过对某高拱坝动力稳定性的分析计算,给出了拱坝坝肩在3个典型高程分别沿两滑移面交线方向、侧滑面及底滑面滑... 利用有限元法及三维实体任意截面数据处理技术,建立了高拱坝坝肩稳定安全系数计算分析方法,并编制了相应的计算和处理程序.通过对某高拱坝动力稳定性的分析计算,给出了拱坝坝肩在3个典型高程分别沿两滑移面交线方向、侧滑面及底滑面滑动的最小安全系数. 展开更多
关键词 有限单元 拱坝 坝肩 动力稳定分析
在线阅读 下载PDF
基于有限单元法的高拱坝动力优化设计方法及其应用 被引量:5
7
作者 苏超 余天堂 姜弘道 《河海大学学报(自然科学版)》 CAS CSCD 北大核心 2002年第1期1-5,共5页
利用分布式网络并行平台 ,建立了并行环境下有约束条件的非线性优化问题的复合形算法 ,解决了基于有限单元法的高拱坝动力优化设计问题 ,并对小湾拱坝进行了以主拉应力最小为目标函数的优化设计 .
关键词 拱坝 有限单元 优化设计 地震 网络 荷载
在线阅读 下载PDF
高精度有限体积格式在三维曲线坐标系下的应用 被引量:2
8
作者 徐丹 王东方 +1 位作者 陈亚铭 邓小刚 《国防科技大学学报》 EI CAS CSCD 北大核心 2016年第2期56-60,共5页
为了构造在光滑区具有较高分辨率并且可以无振荡捕捉激波的高精度有限体积格式,同时降低格式在模板选择上遇到的困难,基于逐维重构方法,发展了结构网格下的高精度有限体积格式,并将这一格式推广到三维曲线坐标系下,从而可以适应相对复... 为了构造在光滑区具有较高分辨率并且可以无振荡捕捉激波的高精度有限体积格式,同时降低格式在模板选择上遇到的困难,基于逐维重构方法,发展了结构网格下的高精度有限体积格式,并将这一格式推广到三维曲线坐标系下,从而可以适应相对复杂外形下的计算。为充分验证格式的有效性,选取一系列典型算例进行计算:在等熵涡输运和二维喷管流动中验证了格式的精度可以达到设计精度;在双马赫反射问题中格式也表现出良好的捕捉激波的能力。数值计算表明,上述格式在曲线网格上具有较高的数值精度和鲁棒的激波捕捉能力,适用于流体力学方程的计算。 展开更多
关键词 有限体积 精度格式 曲线坐标系 逐维重构方
在线阅读 下载PDF
有限体积法中面积分离散格式的精度分析 被引量:1
9
作者 马亮 李亭鹤 《北京航空航天大学学报》 EI CAS CSCD 北大核心 2000年第5期530-534,共5页
讨论了应用顶点中心格式对欧拉方程进行有限体积离散时的精度 ,详细推导了顶点中心格式在平滑网格和扭曲网格上的截断误差 ,同时与中心格式和顶点格式进行了比较 .格式精度分析采用泰勒级数展开的方法 .分析结果表明 :顶点中心格式和顶... 讨论了应用顶点中心格式对欧拉方程进行有限体积离散时的精度 ,详细推导了顶点中心格式在平滑网格和扭曲网格上的截断误差 ,同时与中心格式和顶点格式进行了比较 .格式精度分析采用泰勒级数展开的方法 .分析结果表明 :顶点中心格式和顶点格式类似 ,总能达到至少是一阶的精度 ,并且可以通过加密网格的办法减少误差 ,而且顶点中心格式更适合于分区计算 .通过圆弧亚音速绕流和NACA0 0 1 2翼型绕流流场的数值计算 。 展开更多
关键词 精度分析 泰勒级数 有限体积 面积分离散格式
在线阅读 下载PDF
计算机精度和网格大小对大地电磁有限单元法正演的影响 被引量:6
10
作者 欧东新 《桂林工学院学报》 北大核心 2007年第3期329-332,共4页
以一维大地电磁有限单元模拟为例分析了计算机精度、网格剖分大小对高、低频计算结果的影响.发现在单精度条件下,如果要保证高频段计算精度,需采用较小网格,但此时低频段计算误差较大;采用较大的网格能使低频段精度提高,但降低了高频段... 以一维大地电磁有限单元模拟为例分析了计算机精度、网格剖分大小对高、低频计算结果的影响.发现在单精度条件下,如果要保证高频段计算精度,需采用较小网格,但此时低频段计算误差较大;采用较大的网格能使低频段精度提高,但降低了高频段计算精度.从电磁场的传播特性,以及大地电磁有限单元法正演模拟的计算过程分析了低频段误差的原因:由于低频段电磁波衰减很慢,变分方程中含有差分项,所以网格过小会引起低频段误差. 展开更多
关键词 有限单元 大地电磁 正演 计算精度 网格
在线阅读 下载PDF
利用Ringleb流动对有限体积法精度的研究
11
作者 陈海昕 李凤蔚 +1 位作者 鄂秦 李杰 《西北工业大学学报》 EI CAS CSCD 北大核心 1999年第4期653-658,共6页
利用有精确解的Ringleb 流动构造了用于检验流场计算精度的工具,并用这种工具对中心有限体积法的精度进行了研究。深入讨论了网格、边界处理、人工粘性等对计算精度的影响。结果定量地证明了改善边界处理的精度、网格的正交性... 利用有精确解的Ringleb 流动构造了用于检验流场计算精度的工具,并用这种工具对中心有限体积法的精度进行了研究。深入讨论了网格、边界处理、人工粘性等对计算精度的影响。结果定量地证明了改善边界处理的精度、网格的正交性,包括边界正交性和场内正交性,都可以提高流场计算的精度。为了减少人工粘性对无粘流计算的干扰,应该尽量降低人工粘性系数。 展开更多
关键词 精度 检验 有限体积 Ringleb流动
在线阅读 下载PDF
基于声学的高精度非接触式体积测量方法 被引量:1
12
作者 张馨龙 田光宇 黄勇 《重庆理工大学学报(自然科学)》 CAS 2016年第7期38-44,58,共8页
提出一种基于声学原理的方法以实现精确的非接触物体体积测量。基于分布参数模型,该方法被进一步优化并据此设计出一套声学体积测量仪。借助有限元软件Virtual.Lab构建声学体积计声场的有限元模型。为了得到最优几何参数,改变不同的尺寸... 提出一种基于声学原理的方法以实现精确的非接触物体体积测量。基于分布参数模型,该方法被进一步优化并据此设计出一套声学体积测量仪。借助有限元软件Virtual.Lab构建声学体积计声场的有限元模型。为了得到最优几何参数,改变不同的尺寸,计算相应的测量误差。参数涉及旁通孔的半径、高度、位置,垫片高度、放置方向以及驱动信号等因素。根据仿真得到的最佳形状参数设计出声学体积测量仪进行实验验证。实验测量5 kg,2 kg,100 g的砝码体积,最大误差不超过标称体积的0.6%,而且测量结果稳定,多次测量的标准差不超过0.1 cm3。该方法快速高效,能无损测量,造价低廉,不但适合大规模工业应用,而且能满足实验室使用的精度要求。 展开更多
关键词 声学 精度 体积测量 非接触式测量
在线阅读 下载PDF
有限单元内部场量的高精度求解
13
作者 次惠岭 白冰 +1 位作者 雷宏武 崔银祥 《岩土力学》 EI CAS CSCD 北大核心 2021年第11期3137-3146,共10页
有限单元法是目前最为广泛使用的数值分析方法,能够有效求解各类工程和科学问题。但在一些需要求单元内部场变量的情况下,传统形函数插值方法往往不能保证其结果精度。针对此问题,提出了一种高精度求解方案,基本思路是将有限单元与Taylo... 有限单元法是目前最为广泛使用的数值分析方法,能够有效求解各类工程和科学问题。但在一些需要求单元内部场变量的情况下,传统形函数插值方法往往不能保证其结果精度。针对此问题,提出了一种高精度求解方案,基本思路是将有限单元与Taylor展开相结合,建立并求解单元内部未知场量关于已知单元结点场量的线性方程组,各方程中系数取决于点的相对位置关系。研究发现:所提方案算法简单,可适用于多维度、多阶次单元,在非线性场和高阶场中性能卓越,是一种高效高精度的新型求解算法,可服务于非线性计算的中间步骤,及满足高精度后处理的需求。该方法并不局限于有限元方法,可推广用于包括无网格方法在内的各类离散数值方法。 展开更多
关键词 有限单元 场量 TAYLOR展开 精度
在线阅读 下载PDF
功能梯度材料热弹性分析的高阶格点型有限体积法 被引量:1
14
作者 龚京风 郑文利 +1 位作者 宣领宽 徐宗著 《武汉科技大学学报》 CAS 北大核心 2023年第1期33-39,共7页
为提高功能梯度材料(FGMs)热弹性分析的准确性,本文发展了一种二维高阶格点型有限体积法(CV-FVM)。该方法采用八节点四边形单元(Q8)和六节点三角形单元(T6)离散控制方程,并利用交错网格技术,将材料属性定义于单元中心,将待求解的未知变... 为提高功能梯度材料(FGMs)热弹性分析的准确性,本文发展了一种二维高阶格点型有限体积法(CV-FVM)。该方法采用八节点四边形单元(Q8)和六节点三角形单元(T6)离散控制方程,并利用交错网格技术,将材料属性定义于单元中心,将待求解的未知变量定义于单元节点。应用该方法求解FGMs的热传导、应力和热应力问题,结果与文献中的解析值吻合良好,验证了其正确性。在计算存在大温度梯度的FGMs热传导问题时,高阶CV-FVM能够有效避免数值振荡;在进行复杂FGMs热弹性分析时,高阶CV-FVM能够有效避免数值不连续问题;同时,在相同的网格数下,采用Q8单元的高阶CV-FVM比低阶CV-FVM和采用T6单元的高阶CV-FVM计算精度要高。 展开更多
关键词 格点型有限体积 单元 功能梯度材料 热弹性分析 交错网格技术 数值振荡
在线阅读 下载PDF
曲线网格下精确四阶精度有限体积紧致方法
15
作者 廖飞 叶正寅 《西北工业大学学报》 EI CAS CSCD 北大核心 2012年第6期836-840,共5页
研究了一种求解可压缩欧拉方程的精确四阶精度有限体积紧致方法。通过引入坐标变换,构造了精确四阶精度的体平均量近似和面平均量近似方法,以解决有限体积方法中的积分近似问题,并在曲线网格上辅助四阶精度Padé型紧致格式对欧拉方... 研究了一种求解可压缩欧拉方程的精确四阶精度有限体积紧致方法。通过引入坐标变换,构造了精确四阶精度的体平均量近似和面平均量近似方法,以解决有限体积方法中的积分近似问题,并在曲线网格上辅助四阶精度Padé型紧致格式对欧拉方程进行空间离散。构造了积分型高精度紧致滤波方法代替人工粘性耗散,使计算过程收敛。通过计算欧拉圆柱绕流和Ringleb流动,验证了方法的正确性和有效性。 展开更多
关键词 计算流体力学 欧拉方程 有限体积 积分近似 紧致格式 曲线网格 坐标变换 精确四阶精度
在线阅读 下载PDF
非稳态导热问题有限体积法 被引量:9
16
作者 秦跃平 孟君 +2 位作者 贾敬艳 杨小彬 刘伟 《辽宁工程技术大学学报(自然科学版)》 CAS 北大核心 2013年第5期577-581,共5页
为了研究非稳态导热问题中的常用有限体积方案(FVM1)和新提出的有限体积方案(FVM2)的精度差异问题,采用理论推导和实例验证的方法,建立基于有限体积法原理的热平衡积分方程,利用矩形网格进行离散化,建立了平面非稳态热传导问题有限体积... 为了研究非稳态导热问题中的常用有限体积方案(FVM1)和新提出的有限体积方案(FVM2)的精度差异问题,采用理论推导和实例验证的方法,建立基于有限体积法原理的热平衡积分方程,利用矩形网格进行离散化,建立了平面非稳态热传导问题有限体积法数学模型,推导出与有限单元法(FEM)方案相类似的线性方程组,对比分析所建立的各个方程组中参数的差异.研究结果表明:方案FVM2与FEM精度相同,为较优方案.最后结合一个典型实例,用VB编制了数值计算软件,对比分析这两种方案,验证了上述结果.该结果对工程计算数值模拟有一定的参考价值. 展开更多
关键词 有限体积 有限单元 控制体积 圈划方 非稳态导热 矩形网格 精度差异 数值计算
在线阅读 下载PDF
弹性力学中三种有限体积法方案的对比 被引量:6
17
作者 秦跃平 孙全 +1 位作者 刘伟 张国玉 《辽宁工程技术大学学报(自然科学版)》 CAS 北大核心 2012年第3期349-353,共5页
为了研究弹性力学中三种典型的有限体积法方案(FVM1,FVM2和FVM3)的精度差异问题,采用理论推导和算例验证的方法,建立了基于有限体积法原理的力平衡积分方程,利用三角形网格进行离散化,在计算出内部单元和边界单元对力平衡贡献的基础上,... 为了研究弹性力学中三种典型的有限体积法方案(FVM1,FVM2和FVM3)的精度差异问题,采用理论推导和算例验证的方法,建立了基于有限体积法原理的力平衡积分方程,利用三角形网格进行离散化,在计算出内部单元和边界单元对力平衡贡献的基础上,对每一种有限体积法方案确立了与有限单元法方案(FEM)刚度方程组相类似的线性方程组,对比所建立的各个方程组的待定系数.结果表明:方案FVM2与FEM精度相同,即为最优方案;在不考虑集中力的情况下,方案FVM3的精度也与FEM相同.最后以一个数值算例,对这三种有限体积法方案进行了对比分析,验证了上述分析.该结果对工程计算数值模拟具有一定的参考价值和指导意义. 展开更多
关键词 弹性力学 有限单元 有限体积 有限体积 平面应力问题 最小势能原理 三角形网格 数值计算
在线阅读 下载PDF
基于有限单元法的平纹织物复合材料强度预测:1.RVE的有限元模型 被引量:12
18
作者 左中鹅 王瑞 徐磊 《纺织学报》 EI CAS CSCD 北大核心 2009年第12期45-49,共5页
为更准确预测平纹织物复合材料的强度,基于平纹织物层合板内纱线的细观分布,建立具有双凸结构的特征体积单元RVE(representative volume element)有限元模型。模型内纱线浸渍体为弹性体,物性值随纱线走向变化,树脂为塑性体。利用有限元... 为更准确预测平纹织物复合材料的强度,基于平纹织物层合板内纱线的细观分布,建立具有双凸结构的特征体积单元RVE(representative volume element)有限元模型。模型内纱线浸渍体为弹性体,物性值随纱线走向变化,树脂为塑性体。利用有限元分析软件ANSYS,使用20节点solid单元分别对纱线、树脂进行离散,通过数值分析计算复合材料模量,并根据复合材料内纤维与基体的失效准则预测单胞强度。结果表明:通过有限元模拟得到的弹性模量值同理论计算值及实测值均一致,且预测的单胞强度同复合材料实际强度相吻合;RVE有限元模型支持参数化建模,对平纹织物复合材料的强度预测有很好的应用价值。 展开更多
关键词 有限单元 平纹织物层合板 复合材料 力学性能 强度预测 特征体积单元
在线阅读 下载PDF
基于有限杆单元法的陡坡基桩非线性分析 被引量:15
19
作者 赵明华 陈耀浩 +1 位作者 杨超炜 肖尧 《岩土力学》 EI CAS CSCD 北大核心 2018年第8期3020-3028,共9页
根据陡坡段桥梁基桩承载特性,提出了一种可考虑桩土非线性作用的基桩内力与位移分析有限杆单元方法。首先,基于单元划分结果,结合桩侧摩阻力和坡体侧向推力的分布规律,得到了相应的等效结点荷载向量表示方法;其次,在线弹性地基反力法的... 根据陡坡段桥梁基桩承载特性,提出了一种可考虑桩土非线性作用的基桩内力与位移分析有限杆单元方法。首先,基于单元划分结果,结合桩侧摩阻力和坡体侧向推力的分布规律,得到了相应的等效结点荷载向量表示方法;其次,在线弹性地基反力法的基础上,引入p-y曲线开展了桩周土抗力非线性分析,并结合桩身P-(35)效应(P为桩顶轴向荷载,(35)为桩顶水平位移)计算方法给出了单元刚度矩阵修正方法,进而提出了适用于基桩内力位移非线性分析的有限杆单元法并编制了MATLAB计算程序;最后,结合某工程实例,将计算结果与工程实测值及已有理论值进行对比。结果表明:考虑桩土非线性作用的计算方法是合理的;当桩身具有自由段时,P-(35)效应对基桩内力位移的影响较大,在实际工程设计中不可忽视。 展开更多
关键词 陡斜坡 桩基础 有限单元 P-Y曲线 P-△效应
在线阅读 下载PDF
适用于任意网格拓扑和质量的格心有限体积法 被引量:9
20
作者 蒋跃文 叶正寅 《力学学报》 EI CSCD 北大核心 2010年第5期830-837,共8页
针对格心有限体积法的离散精度易受网格类型影响的问题,基于最小二乘原理,提出了一种适用于任意网格拓扑和网格质量的有限体积方法.在解的光滑区能保证二阶精度,可光滑、陡峭地捕捉激波等强间断面.精度与网格无关,且算法统一的特性使其... 针对格心有限体积法的离散精度易受网格类型影响的问题,基于最小二乘原理,提出了一种适用于任意网格拓扑和网格质量的有限体积方法.在解的光滑区能保证二阶精度,可光滑、陡峭地捕捉激波等强间断面.精度与网格无关,且算法统一的特性使其非常适用于网格自适应和多重网格等计算应用,同时也降低了网格生成和流场求解的复杂性.跨音速算例的网格含有多种网格拓扑,计算结果表明发展的线性重构方法(linearreconstruction method,LRM)适用于不同的网格拓扑,计算得到的激波位置准确、陡峭,未产生数值振荡.运用Ringleb流动考查了该方法对低质量网格的收敛性,与传统方法相比,线性重构方法(LRM)不仅平均误差较小,而且误差随网格尺度的收敛性也更好,其精度接近二阶.三段翼型的黏性绕流计算进一步表明网格质量对其精度的影响较小. 展开更多
关键词 网格拓扑 网格质量 有限体积 线性重构方 精度
在线阅读 下载PDF
上一页 1 2 4 下一页 到第
使用帮助 返回顶部