期刊文献+
共找到87篇文章
< 1 2 5 >
每页显示 20 50 100
管式发射燃气流动计算方法及特性研究
1
作者 谢军虎 赵英豪 +2 位作者 杜文强 李浩镕 傅德彬 《兵器装备工程学报》 北大核心 2025年第7期162-167,共6页
燃气流动对火箭、导弹发射安全有着重要影响,采用有效的数值方法模拟发射时燃气流场和载荷特性具有重要价值。为探究火箭、导弹管式热发射时燃气流场特性,围绕发射过程中自由射流和发射管内受限流动状态,分别采用局部时间步长定常模型... 燃气流动对火箭、导弹发射安全有着重要影响,采用有效的数值方法模拟发射时燃气流场和载荷特性具有重要价值。为探究火箭、导弹管式热发射时燃气流场特性,围绕发射过程中自由射流和发射管内受限流动状态,分别采用局部时间步长定常模型、全局时间步长定常模型和非定常模型进行求解计算。结果表明:局部时间步长和全局时间步长能够获得相吻合的燃气自由射流结构,但不同模型获得的发射管内受限燃气流场存在显著差异;局部时间步长定常模型收敛于局部平衡状态,发射管内出现明显雍塞;全局时间步长定常模型收敛于完全发展的无雍塞流动状态;非定常模型则出现先雍塞,而后雍塞区逐渐后移并消失的时变过程。 展开更多
关键词 燃气流动 发射管 时间推进法 局部时间步长 雍塞流动
在线阅读 下载PDF
流体计算中基于流场解预估的时间步长约束方法
2
作者 陈科锦 罗灿炎 +2 位作者 高鹤 毕林 王皓 《力学学报》 北大核心 2025年第6期1351-1361,共11页
时间推进算法的稳定性是决定流体计算成败的关键.然而,即便是隐式推进算法,在启动计算的早期阶段,或者流动中存在激波、分离等大梯度的流动特征的情况,也只能采用较小的时间步长进行推进,以免计算发散.文章从守恒形式出发,通过近似线性... 时间推进算法的稳定性是决定流体计算成败的关键.然而,即便是隐式推进算法,在启动计算的早期阶段,或者流动中存在激波、分离等大梯度的流动特征的情况,也只能采用较小的时间步长进行推进,以免计算发散.文章从守恒形式出发,通过近似线性化发展了原始变量预估方法,并基于预估结果建立了一种新的时间步长约束方法,实现对每个计算单元时间步长的动态调整,无需人工手动调参,确保了启动计算初期和大梯度流动模拟的稳定性,同时能提高流动充分发展稳定后的计算效率.算例结果表明,文章方法能够将流场变量推进值控制在合理范围内,有效防止了强非线性特性等造成的计算发散情况,提升了数值求解的鲁棒性和可靠性,能够更快更稳定地收敛. 展开更多
关键词 计算流体力学 时间推进 解预估
在线阅读 下载PDF
基于时间推进的轴流涡轮一维性能分析
3
作者 唐晴 杨晨 吴虎 《推进技术》 北大核心 2025年第3期91-100,共10页
为实现轴流涡轮性能的快速分析,同时模拟叶片通道内的流动分布,发展了时间推进一维控制方程模型,并进一步开发了轴流涡轮一维性能计算程序。为模拟叶片对气流偏转作用,引入一种鲁棒的无粘叶片力模型,同时融合损失经验模型模拟粘性损失... 为实现轴流涡轮性能的快速分析,同时模拟叶片通道内的流动分布,发展了时间推进一维控制方程模型,并进一步开发了轴流涡轮一维性能计算程序。为模拟叶片对气流偏转作用,引入一种鲁棒的无粘叶片力模型,同时融合损失经验模型模拟粘性损失效应。利用该程序对KJ66发动机单级涡轮以及E^(3)发动机两级高压涡轮展开特性分析,并将预测值与相关实验数据或三维数值仿真进行了对比验证。两款涡轮单工况点的计算时长均不超过1 s。结果表明,该模型能快速准确评估涡轮特性,具备自然模拟涡轮堵塞工况的能力,在涡轮设计阶段性能评估分析中具备较好的应用潜力。 展开更多
关键词 轴流涡轮 时间推进 一维模型 无粘叶片力模型 性能分析
在线阅读 下载PDF
基于扩散飞行时间的复杂裂缝系统数值试井方法
4
作者 杜雅雯 罗万静 +2 位作者 滕柏路 李佳 韩家强 《深圳大学学报(理工版)》 北大核心 2025年第3期291-298,共8页
为精确描述水力压裂后复杂裂缝系统的压力动态特征,提出一种基于扩散飞行时间(diffusive time of flight,DTOF)的快速数值试井分析方法.将嵌入式离散裂缝模型(embedded discrete fracture model,EDFM)与波前快速推进法(fast marching me... 为精确描述水力压裂后复杂裂缝系统的压力动态特征,提出一种基于扩散飞行时间(diffusive time of flight,DTOF)的快速数值试井分析方法.将嵌入式离散裂缝模型(embedded discrete fracture model,EDFM)与波前快速推进法(fast marching method,FMM)相结合,建立以一维扩散飞行时间为空间坐标的渗流控制方程,并以此分析复杂裂缝系统试井曲线特征.研究表明,天然裂缝对非稳态渗流有较大影响,在天然裂缝和人工裂缝组合下,试井曲线表现出多种不同特征流动段的组合:仅考虑天然裂缝时,双对数曲线表现为井筒储集阶段、径向流阶段、过渡阶段和边界流阶段;考虑单一人工裂缝时,双对数曲线表现为井筒储集阶段、双线性流阶段、径向流阶段和边界流阶段;随着人工裂缝条数增加,储层的连通性得到较大改善,双对数曲线则表现为井筒储集阶段、线性流阶段和边界流阶段.本研究构建的数值模型在现场实例中得到了良好的应用,为复杂裂缝系统试井提供了新的方法. 展开更多
关键词 油气田开发 扩散飞行时间 波前快速推进法 数值试井 裂缝性储层 嵌入式离散裂缝模型 试井曲线
在线阅读 下载PDF
基于.NET Remoting分布式计算电大目标瞬态散射特性
5
作者 王文举 任猛 +1 位作者 周东明 何建国 《微波学报》 CSCD 北大核心 2009年第2期11-14,共4页
电磁场时域积分方程方法(TDIE)在解决瞬态电磁问题和宽带电磁问题等方面显示出独特的优势,然而TDIE的数值计算方法——时间步进算法(MOT)的计算规模较大,严重阻碍了TDIE方法在求解电大目标瞬态电磁特性中的应用。文章在研究MOT并行算法... 电磁场时域积分方程方法(TDIE)在解决瞬态电磁问题和宽带电磁问题等方面显示出独特的优势,然而TDIE的数值计算方法——时间步进算法(MOT)的计算规模较大,严重阻碍了TDIE方法在求解电大目标瞬态电磁特性中的应用。文章在研究MOT并行算法的基础上,开发了基于.NET Remoting的电磁场分布式数值计算方案。数值结果表明,该分布式并行方案显著提高了TDIE的计算效率,为解决电大目标瞬态电磁散射问题提供了一条有效途径。 展开更多
关键词 时域积分方程方法 时间步进算法 分布式计算 .NETRemoting
在线阅读 下载PDF
时域积分方程MOT算法的推迟位时间卷积数值积分新方法
6
作者 贾苗苗 赵延文 +1 位作者 郑宇腾 蔡强明 《电子学报》 EI CAS CSCD 北大核心 2015年第1期151-159,共9页
通过变量代换平滑三角形上推迟位(标量位函数和矢量位函数)并消除推迟矢量位旋度的奇异性,使得采用数值积分法就能够精确快速地计算任意正则时间基函数与推迟位函数及推迟矢量位旋度之间的时间卷积运算,可用于基于任意类型时间基函数的... 通过变量代换平滑三角形上推迟位(标量位函数和矢量位函数)并消除推迟矢量位旋度的奇异性,使得采用数值积分法就能够精确快速地计算任意正则时间基函数与推迟位函数及推迟矢量位旋度之间的时间卷积运算,可用于基于任意类型时间基函数的时域电场、时域磁场及其混合场积分方程时间步进(MOT)算法.与时间卷积运算的解析法对比分析表明,该时间卷积数值积分方法能够精确快速地计算基于任意类型时间基函数和不同时间步长条件下时域积分方程MOT算法的阻抗矩阵元素;而具体的计算实例也表明,阻抗矩阵的精确计算显著地提升了时域积分方程MOT算法的后时稳定性和求解精度. 展开更多
关键词 时域积分方程 时间步进算法 后时不稳定性
在线阅读 下载PDF
利用时间步进算法精确稳定求解时域积分方程 被引量:5
7
作者 赵庆广 赵延文 +1 位作者 毕海燕 聂在平 《电子学报》 EI CAS CSCD 北大核心 2008年第6期1135-1139,共5页
时域阻抗矩阵元素的计算需要分别计算场单元和源单元上的空时积分,由于时间基函数的分域性以及时间基函数(如三角型时间基函数)导数的不连续性,使得采用高斯积分方法计算源单元上空时积分的计算精度较差且误差随着时间步长的减小而增大... 时域阻抗矩阵元素的计算需要分别计算场单元和源单元上的空时积分,由于时间基函数的分域性以及时间基函数(如三角型时间基函数)导数的不连续性,使得采用高斯积分方法计算源单元上空时积分的计算精度较差且误差随着时间步长的减小而增大.本文通过将源单元上空时积分转变成为1D时间卷积分和1D空间解析积分来精确计算时域阻抗矩阵元素,并在此基础上利用时间步进算法求解了时域电场、磁场和混合场积分方程.通过计算实例表明该方法在较大的时间步长取值范围内均能确保时域积分方程时间步进算法求解的精度和后时稳定性. 展开更多
关键词 时域积分方程(TDIE) 时间步进算法(mot) 混合场积分方程(CFIE)
在线阅读 下载PDF
精确快速计算时域积分方程中奇异性积分的新方法 被引量:2
8
作者 赵延文 徐建华 +1 位作者 聂在平 武胜波 《电子与信息学报》 EI CSCD 北大核心 2005年第11期1821-1824,共4页
该文首先利用参数坐标和广义Duffy 坐标变换将时域电场积分方程(TDEFIE)的奇异性积分转换成非奇异性积分,然后根据时间基函数的特点将该积分转换成可以快速精确计算的分区域积分。数值计算实例表明,该方法可以大幅度提高求解TDEFIE 的... 该文首先利用参数坐标和广义Duffy 坐标变换将时域电场积分方程(TDEFIE)的奇异性积分转换成非奇异性积分,然后根据时间基函数的特点将该积分转换成可以快速精确计算的分区域积分。数值计算实例表明,该方法可以大幅度提高求解TDEFIE 的后时稳定性和解的精度,而不必采用任何求平均的过程。该方法适用于任意类型的时间基函数并可方便地推广到高阶曲面拟合和高阶空间基函数情形。 展开更多
关键词 时域电磁散射 时域积分方程 时间步进算法 奇异性积分
在线阅读 下载PDF
TD-EFIE分析任意几何形状导体的瞬态特性 被引量:4
9
作者 姜光兴 曹伟 朱洪波 《电子学报》 EI CAS CSCD 北大核心 2006年第4期741-745,共5页
应用时域电场积分方程方法分析了各类天线和散射体的瞬态电磁响应.在分析中,激励源的形式和导体的几何形状可以任意.首先探讨了不同时域基函数对计算结果的影响,同时为了消除后期振荡,而且保持计算结果的准确性,还提出了一种新的均衡方... 应用时域电场积分方程方法分析了各类天线和散射体的瞬态电磁响应.在分析中,激励源的形式和导体的几何形状可以任意.首先探讨了不同时域基函数对计算结果的影响,同时为了消除后期振荡,而且保持计算结果的准确性,还提出了一种新的均衡方案.计算了一系列的实例,并将计算结果与采用离散傅立叶反变换法计算得到的结果进行了比较.它们之间良好的一致性说明了本文所提方法的正确性和有效性. 展开更多
关键词 电磁散射 电磁辐射 时域基函数 时间步进法 时域电场积分方程
在线阅读 下载PDF
精确稳定求解时域电场积分方程的一种新方法 被引量:6
10
作者 赵延文 聂在平 +1 位作者 徐建华 武胜波 《电子学报》 EI CAS CSCD 北大核心 2006年第6期1104-1108,共5页
时域电场、磁场和混合场积分方程已被广泛用来分析散射体的时域散射响应.基于适当的空间积分方法和隐式的时间步进算(MOT)法在求解时域磁场和混合场积分方程时总是稳定的,然而在求解TDEFIE时则是不稳定的.在本文中,时域电场积分方程的... 时域电场、磁场和混合场积分方程已被广泛用来分析散射体的时域散射响应.基于适当的空间积分方法和隐式的时间步进算(MOT)法在求解时域磁场和混合场积分方程时总是稳定的,然而在求解TDEFIE时则是不稳定的.在本文中,时域电场积分方程的非奇异性积分采用标准的高斯求积法来计算;而利用参数坐标变换和极坐标变换将其奇异性积分转换成为可以分区域精确快速计算的非奇异性积分.通过数值实验表明,利用该方法可以非常精确稳定地求解时域电场积分方程,即使是在时间迭代后期也不必采用任何求平均的过程;另外,该方法可以用于任意时间基函数并可以推广到高阶空间基函数的情形. 展开更多
关键词 时域积分方程 时域电磁散射 时间步进算法 奇异性积分
在线阅读 下载PDF
两种自由尾迹模型在风力机气动性能预测中的应用 被引量:7
11
作者 沈昕 竺晓程 杜朝辉 《太阳能学报》 EI CAS CSCD 北大核心 2010年第7期923-927,共5页
分别采用空间迭代法以及时间步进法的自由尾迹计算模型对风力机的尾迹进行预测。结合相应的算例以及模型特点预测风力机在正对来流、偏航以及变工况下的风力机尾迹形状并与相应的实验值进行比较。结果表明,在计算稳态工况时,空间迭代法... 分别采用空间迭代法以及时间步进法的自由尾迹计算模型对风力机的尾迹进行预测。结合相应的算例以及模型特点预测风力机在正对来流、偏航以及变工况下的风力机尾迹形状并与相应的实验值进行比较。结果表明,在计算稳态工况时,空间迭代法可以快速稳定的预测出尾迹形状,而时间步进法在计算稳态工况时可以很好的收敛到初始值,而其瞬态特性保证该方法可以计算风力机在变工况下的气动特性。两种方法得到的尾迹形状基本反映叶尖涡结构的特性,可以用来作为升力线、升力面等工程模型算法的尾迹部分,用以进行水平轴风力机的气动性能预测。 展开更多
关键词 自由尾迹 空间迭代 时间步进 水平轴风力机
在线阅读 下载PDF
基于Laguerre多项式的边界积分方程时域求解 被引量:1
12
作者 柴舜连 王生水 +1 位作者 陈忠宽 毛钧杰 《电子与信息学报》 EI CSCD 北大核心 2008年第3期742-745,共4页
针对时域积分方程中存在的晚时震荡问题,介绍了基于Laguerre多项式的电场、磁场和混合场积分方程,求解了导体球和导体圆柱的时域电流分布和后向散射场以及单站RCS。结果表明,3种积分方程很好地解决了晚时震荡问题,混合场积分方程具有更... 针对时域积分方程中存在的晚时震荡问题,介绍了基于Laguerre多项式的电场、磁场和混合场积分方程,求解了导体球和导体圆柱的时域电流分布和后向散射场以及单站RCS。结果表明,3种积分方程很好地解决了晚时震荡问题,混合场积分方程具有更高的计算精度。 展开更多
关键词 单站雷达散射截面 拉盖尔多项式 混合场积分方程 阶数步进 时间步进
在线阅读 下载PDF
时间推进方法在叶轮机械内部流场计算中的进展 被引量:5
13
作者 杨策 蒋滋康 索沂生 《力学进展》 EI CSCD 北大核心 2000年第1期83-94,共12页
时间推进方法目前广泛应用于叶轮机械内部流场的计算,这种方法可以划分成显式时间推进方法和隐式时间推进方法.本文简要介绍了时间推进方法用于叶轮机械内部流场计算的起源及这种方法的优点和缺点,重点记述了国内外叶轮机界应用显式... 时间推进方法目前广泛应用于叶轮机械内部流场的计算,这种方法可以划分成显式时间推进方法和隐式时间推进方法.本文简要介绍了时间推进方法用于叶轮机械内部流场计算的起源及这种方法的优点和缺点,重点记述了国内外叶轮机界应用显式时间推进方法和隐式时间推进方法计算叶轮机内部流场的最新进展。 展开更多
关键词 叶轮机械 时间推进 数值计算 三维流场
在线阅读 下载PDF
一致性绕射理论加速的时域积分方程方法 被引量:2
14
作者 周东明 覃宇建 何建国 《微波学报》 CSCD 北大核心 2010年第3期34-37,共4页
为计算电大平台上天线间的互耦提出一种一致性绕射理论加速的时域积分方程算法。引入一致性绕射理论中射线寻迹的方法,对占时间步进法主要计算量的迟滞积分项进行加速,以牺牲部分精度为代价,提高了时域积分方程算法的计算速度。数值算... 为计算电大平台上天线间的互耦提出一种一致性绕射理论加速的时域积分方程算法。引入一致性绕射理论中射线寻迹的方法,对占时间步进法主要计算量的迟滞积分项进行加速,以牺牲部分精度为代价,提高了时域积分方程算法的计算速度。数值算例证明了有效性。 展开更多
关键词 互耦 电磁兼容 一致性绕射理论 时域积分方程 时间步进法
在线阅读 下载PDF
利用时域电场积分方程分析天线辐射问题 被引量:5
15
作者 王生水 毛钧杰 +2 位作者 姜永金 柴舜连 陈忠宽 《电波科学学报》 EI CSCD 北大核心 2007年第5期837-841,共5页
研究了基于矩量法和RWG三角基函数的隐式电场积分方程的时域算法,引入了一种激励源的时域设置方法。利用直接在时域加激励源的方法来分析天线辐射问题,所得时域数据经傅立叶变换可得到很宽频带的频域数据,与在频域逐个频点求解相比大大... 研究了基于矩量法和RWG三角基函数的隐式电场积分方程的时域算法,引入了一种激励源的时域设置方法。利用直接在时域加激励源的方法来分析天线辐射问题,所得时域数据经傅立叶变换可得到很宽频带的频域数据,与在频域逐个频点求解相比大大节省了计算时间。通过对几种典型天线的分析计算,验证了方法的正确性和算法的稳定性。 展开更多
关键词 矩量法 时域电场积分方程 隐式时间步进算法
在线阅读 下载PDF
充气结构与流场的耦合求解方法 被引量:5
16
作者 蒋跃文 叶正寅 张正科 《力学学报》 EI CSCD 北大核心 2010年第1期1-7,共7页
针对一种充气前缘增升技术,建立了其充气结构与流场耦合作用的运动方程.将方程写成状态空间形式,采用时域推进方法求解.对使用了变前缘增升技术的NACA 63-212翼型进行了充气结构静变形的数值计算,结果表明充气结构的刚度对翼型的气动特... 针对一种充气前缘增升技术,建立了其充气结构与流场耦合作用的运动方程.将方程写成状态空间形式,采用时域推进方法求解.对使用了变前缘增升技术的NACA 63-212翼型进行了充气结构静变形的数值计算,结果表明充气结构的刚度对翼型的气动特性有明显影响.与原翼型相比,在不考虑充气结构变形时,该增升技术大约能使翼型的失速迎角增加30%,最大升力系数增加22%;考虑结构变形后增升效果有所降低.刚度较低的薄膜在前缘吸力峰的作用下会隆起形成鼓包,容易引起流动分离. 展开更多
关键词 充气结构 结构运动方程 耦合求解 时域推进
在线阅读 下载PDF
求解时域电场积分方程的稳定时间步进算法 被引量:1
17
作者 李会容 赵延文 +1 位作者 张雪峰 武昭妤 《电波科学学报》 EI CSCD 北大核心 2009年第6期1129-1136,共8页
利用面积坐标变换、相对坐标变换、积分区域分解和广义Duffy坐标变换将时域电场积分方程中奇异性积分(共面、共边和共单结点的场源三角形单元上)转化成可精确计算的非奇异性积分。在不同时间基函数(导数连续和导数不连续)、不同时间步... 利用面积坐标变换、相对坐标变换、积分区域分解和广义Duffy坐标变换将时域电场积分方程中奇异性积分(共面、共边和共单结点的场源三角形单元上)转化成可精确计算的非奇异性积分。在不同时间基函数(导数连续和导数不连续)、不同时间步长情况下对比分析了该方法和现有的常用方法计算奇异性积分的精度。计算实例表明:时域阻抗矩阵的精确计算有效地改善了时间步进算法的后时稳定性。 展开更多
关键词 时域电场积分方程 奇异性积分 时间步进算法 后时不稳定性
在线阅读 下载PDF
基于时间推进自由尾迹法的水平轴风力机气动特性分析 被引量:9
18
作者 周文平 唐胜利 吕红 《中国电机工程学报》 EI CSCD 北大核心 2011年第29期124-130,共7页
建立适合于水平轴风力机气动特性分析的自由尾迹计算方法。该方法中,桨叶模型采用Weissinger-L升力面模型,转子自由尾迹的求解采用时间推进方法,数值求解方法采用4阶Adams-Moulton预估–校正算法以提高计算精度。应用所建立的方法对模... 建立适合于水平轴风力机气动特性分析的自由尾迹计算方法。该方法中,桨叶模型采用Weissinger-L升力面模型,转子自由尾迹的求解采用时间推进方法,数值求解方法采用4阶Adams-Moulton预估–校正算法以提高计算精度。应用所建立的方法对模型风力机处于轴向来流工况时的尾迹形状和气动特性进行计算,并与文献实验数据进行比较,验证了该方法的有效性。最后计算了来流风偏转和桨叶桨距角增加时尾迹形状的非定常畸变及气动载荷特性,得出一些有意义的结论。 展开更多
关键词 水平轴风力机 时间推进 自由尾迹 气动特性
在线阅读 下载PDF
一种新的界面映射推进方法及其在气动弹性力学中的应用 被引量:5
19
作者 安效民 徐敏 陈士橹 《宇航学报》 EI CAS CSCD 北大核心 2008年第5期1473-1479,共7页
非线性气动弹性研究中,涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD)耦合问题,在耦合界面上不仅要进行两场之间信息的传递,而且要有匹配的时间推进格式。针对计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法,基于... 非线性气动弹性研究中,涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD)耦合问题,在耦合界面上不仅要进行两场之间信息的传递,而且要有匹配的时间推进格式。针对计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法,基于边界能量守恒,发展了一种新的界面处理方法,包括:①界面映射方法,该方法基于耦合边界上局部的网格信息,将两场之间载荷信息和位移信息的转换放在同一个映射矩阵中来处理,克服了占用大量CPU时间和内存的需求;②界面推进方法,该法基于松耦合计算流程,通过引入半步交错推进,达到了二阶时间精度,提高了计算精度和效率。最后将该界面处理方法应用于柔性大展弦比机翼的气动弹性计算和AGARD445.6机翼的动响应分析中。结果表明该方法能够高效、高精度地处理不同物理场之间的数据交换和时间推进,并具有处理复杂几何体非线性气动弹性问题的能力。 展开更多
关键词 气动弹性 CFD/CSD耦合系统 界面技术 数据转换 时间推进
在线阅读 下载PDF
用粘性体力方法计算轴流叶轮机械内部流场 被引量:2
20
作者 杨策 蒋滋康 +1 位作者 索沂生 老大中 《应用力学学报》 CAS CSCD 北大核心 2000年第4期40-44,共5页
阐述了一种用于叶轮机械内部三维粘性流场计算程序。该程序以有限体积显式时间推进方法为基础,用Baldwin-Lomax模型模拟湍流流动,用局部时间步长和多重网格方法提高计算效率。使用上述方法计算了NASA 37#跨音速压气机转子流场,并与... 阐述了一种用于叶轮机械内部三维粘性流场计算程序。该程序以有限体积显式时间推进方法为基础,用Baldwin-Lomax模型模拟湍流流动,用局部时间步长和多重网格方法提高计算效率。使用上述方法计算了NASA 37#跨音速压气机转子流场,并与实验结果进行了比较,从而证明文中提出的方法能够快速、准确地模拟轴流叶轮机械内部复杂三维流场,该方法具有较强的工程实用意义。 展开更多
关键词 叶轮机械 湍流模型 多重网格 有限体积 时间推进方法 压气机 粘性体力方法 流场
在线阅读 下载PDF
上一页 1 2 5 下一页 到第
使用帮助 返回顶部