期刊文献+
共找到64篇文章
< 1 2 4 >
每页显示 20 50 100
RETHINKING TO FINITE DIFFERENCE TIME-STEP INTEGRATIONS 被引量:1
1
作者 钟万勰 朱建平 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 1995年第8期705-711,共7页
The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integ... The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method. 展开更多
关键词 finite difference method time step integration numericalstability
在线阅读 下载PDF
Solving Stiff Reaction-Diffusion Equations Using Exponential Time Differences Methods
2
作者 H. A. Ashi 《American Journal of Computational Mathematics》 2018年第1期55-67,共13页
Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Th... Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Though, this paper shows that recent advance methods can be more favored. In this work, we have incorporated, throughout numerical comparison experiments, spectral methods, for the space discretization, in conjunction with second and fourth-order time integrating methods for approximating the solution of the reaction-diffusion differential equations. The results have revealed that these methods have advantages over the conventional methods, some of which to mention are: the ease of implementation, accuracy and CPU time. 展开更多
关键词 finite difference methodS EXPONENTIAL INTEGRATOR EXPONENTIAL time Differencing method REACTION-DIFFUSION System
在线阅读 下载PDF
一种适用于嵌入式导电薄层的高阶电磁波混合时域有限差分-时程精细积分法
3
作者 马亮 马西奎 +2 位作者 迟明珺 向汝 朱晓杰 《电工技术学报》 北大核心 2025年第5期1333-1343,共11页
应用单一时域数值方法,在面对嵌入式导电薄层一类多尺度问题时,都面临着建模极为困难的挑战。该文提出了一种基于时域有限差分(FDTD)法和高阶时程精细积分(PITD)法的电磁波混合数值方法。该方法对导电薄层外部进行粗网格剖分并应用FDTD... 应用单一时域数值方法,在面对嵌入式导电薄层一类多尺度问题时,都面临着建模极为困难的挑战。该文提出了一种基于时域有限差分(FDTD)法和高阶时程精细积分(PITD)法的电磁波混合数值方法。该方法对导电薄层外部进行粗网格剖分并应用FDTD法,而对薄层内部进行一维细网格剖分并应用四阶PITD法,以实现不同网格尺度的同步时间推进。为了实现粗细网格之间的信息交换,在PITD域中引入过渡区域并应用二阶PITD法,通过等效本构参数来更新交界面处的切向电场。分析了该混合算法的数值稳定性和数值反射,并通过典型数值算例验证了所提方法的有效性和准确性。 展开更多
关键词 时域有限差分法 四阶时程精细积分法 亚网格技术 矩阵指数
在线阅读 下载PDF
标量波方程时间高阶广义有限差分法及稳定性条件
4
作者 袁野 黄健良 +2 位作者 陶维祥 柳万春 吴国忱 《科学技术与工程》 北大核心 2025年第21期8796-8804,共9页
广义有限差分地震波场数值模拟方法能够适应起伏的地层界面,消除起伏界面造成的阶梯状散射现象,提高正演模拟的准确性。而使用时间二阶广义有限差分法求解波动方程时,由于时间差分精度低,时间采样间隔较大时往往会产生时间频散,影响正... 广义有限差分地震波场数值模拟方法能够适应起伏的地层界面,消除起伏界面造成的阶梯状散射现象,提高正演模拟的准确性。而使用时间二阶广义有限差分法求解波动方程时,由于时间差分精度低,时间采样间隔较大时往往会产生时间频散,影响正演模拟的精度。研究了标量波方程时间四阶广义有限差分正演模拟算法及其稳定性条件和频散特性,通过将时间四阶偏导数转嫁到空间偏导数项上实现时间四阶精度差分,时间频散得到有效的压制。此外相对于时间二阶广义有限差分法,时间四阶广义有限差分可以适应较大的时间采样间隔,一定程度上减少计算量。实验结果表明,所提出的算法能有效压制阶梯状散射和时间频散,具有更高的计算精度,将其应用于逆时偏移,可以获得高质量的成像剖面。 展开更多
关键词 广义有限差分 阶梯状散射 高阶差分格式 时间频散 稳定性条件
在线阅读 下载PDF
时域电场积分方程的稳定求解 被引量:12
5
作者 赵延文 聂在平 +1 位作者 徐建华 武胜波 《电波科学学报》 EI CSCD 2004年第2期148-152,共5页
提出了一种基于时域电场积分方程的稳定精确数值方法计算任意形状理想导体的时域散射。矢量位的时间导数采用中心差分近似 ,而标量位采用时间平均表示 ,并且标量位和矢量位在单元上随时间变化采用单元中心处的值来计算。将该方法与隐式... 提出了一种基于时域电场积分方程的稳定精确数值方法计算任意形状理想导体的时域散射。矢量位的时间导数采用中心差分近似 ,而标量位采用时间平均表示 ,并且标量位和矢量位在单元上随时间变化采用单元中心处的值来计算。将该方法与隐式方法相结合可以得到稳定精确的求解结果 ,而不必对电流密度进行时间或空间平均的处理过程。从模拟结果证明了该方法的稳定性和精确性。 展开更多
关键词 时域电场 积分方程 时域电磁散射 有限差分法 差分近似
在线阅读 下载PDF
时间有限元与保辛 被引量:30
6
作者 钟万勰 姚征 《机械强度》 EI CAS CSCD 北大核心 2005年第2期178-183,共6页
初值问题的时间积分经常采用差分近似。保守体系的特点是保辛。但通常的差分格式并不考虑保辛的性质,即使对保守体系。但有限元是自动保辛的,即使采用小参数摄动时,仍能保辛。文中先验证时间积分有限元的数值效果。再考虑齐次有阻尼Duff... 初值问题的时间积分经常采用差分近似。保守体系的特点是保辛。但通常的差分格式并不考虑保辛的性质,即使对保守体系。但有限元是自动保辛的,即使采用小参数摄动时,仍能保辛。文中先验证时间积分有限元的数值效果。再考虑齐次有阻尼Duffing方程的积分,它也可变换到时变的Hamilton系统,采用正则变换法摄动然后保辛积分可得满意结果。进一步又数值积分了二自由度齐次有阻尼Duffing方程的积分。 展开更多
关键词 时程积分 DUFFING方程 保辛 有限元 有限差分 摄动
在线阅读 下载PDF
Kirchhoff偏移法在探地雷达正演图像处理中的应用 被引量:33
7
作者 戴前伟 冯德山 何继善 《地球物理学进展》 CSCD 北大核心 2005年第3期849-853,共5页
本文首先通过对波动方程的分析,得出了声波波动方程和雷达波波动方程形式一致性,从而说明了把广泛应用在地震数据处理中的偏移技术引入到GPR资料处理中的可行性;其后说明了时域有限差分法(FDTD)的原理,并用它合成了几种常见的雷达正演剖... 本文首先通过对波动方程的分析,得出了声波波动方程和雷达波波动方程形式一致性,从而说明了把广泛应用在地震数据处理中的偏移技术引入到GPR资料处理中的可行性;其后说明了时域有限差分法(FDTD)的原理,并用它合成了几种常见的雷达正演剖面;最后利用Kirchhoff积分偏移法对正演所得的雷达剖面进行偏移处理,通过对比偏移处理前后的雷达正演剖面,可知Kirchhoff积分偏移法能使雷达正演剖面中的反射波的归位,绕射波收敛,从而大大地提高了雷达正演剖面的分辨率,更好地指导GPR剖面的地质解释和验证偏移方法的有效性. 展开更多
关键词 克希霍夫(Kirchhoff)积分法 探地雷达(GPR) 有限差分法 正演模拟 偏移
在线阅读 下载PDF
细线结构时域电场积分方程的有限差分求解 被引量:27
8
作者 谢秦川 刘顺坤 陈雨生 《电波科学学报》 EI CSCD 2000年第1期45-48,共4页
提出一种求解基于细线结构的时域电场积分方程 ( TDEFIE)的方法 -有限差分方法。与传统求解时域电场积分方程的时域矩量法相比 ,该方法易于数值实现。文中还利用该方法研究了电磁辐射、散射的三个经典问题 ,并将结果与时域矩量法的结果... 提出一种求解基于细线结构的时域电场积分方程 ( TDEFIE)的方法 -有限差分方法。与传统求解时域电场积分方程的时域矩量法相比 ,该方法易于数值实现。文中还利用该方法研究了电磁辐射、散射的三个经典问题 ,并将结果与时域矩量法的结果进行了比对 ,研究表明 ,该方法求解细线结构的 TDEFIE非常有效。 展开更多
关键词 电场积分方程 有限差分 细线结构
在线阅读 下载PDF
RTM充模过程数值模拟的隐式有限元算法 被引量:10
9
作者 谭华 晏石林 王继辉 《固体力学学报》 CAS CSCD 北大核心 2004年第3期349-354,共6页
建立了基于欧拉方法描述树脂传递模塑 (RTM)工艺充模过程的基本数学方程 ,并采用有限元隐式时间积分方法对基本方程进行了数值求解 .编制了基于隐式有限元算法及传统有限元控制体算法的程序 ,通过具体算例比较了这两种算法的优缺点 .与... 建立了基于欧拉方法描述树脂传递模塑 (RTM)工艺充模过程的基本数学方程 ,并采用有限元隐式时间积分方法对基本方程进行了数值求解 .编制了基于隐式有限元算法及传统有限元控制体算法的程序 ,通过具体算例比较了这两种算法的优缺点 .与传统的有限元控制体法相比 ,该文提出的隐式有限元算法能节省计算时间 ,特别适合于单元、节点数目多的情况 .隐式有限元算法是一种纯有限元方法 ,不需要使用控制体积技术 ,采用该算法计算出的流动前沿与时间步长无关 . 展开更多
关键词 树脂传递模塑 Darcy定律 有限元法 隐式时间积分 时间步长 材料力学
在线阅读 下载PDF
对差分法时程积分的反思 被引量:15
10
作者 钟万勰 朱建平 《应用数学和力学》 CSCD 北大核心 1995年第8期663-668,共6页
以往偏微分方程时间步的数值积分主要由有限差分法来执行,然而当时间步长较大时会引起数值不稳定性。本文给出的单点精细积分法导出的显式积分格式可证明是无条件稳定的。就扩散方程与对流─扩散方程作出了本文方法与差分法导出的格式... 以往偏微分方程时间步的数值积分主要由有限差分法来执行,然而当时间步长较大时会引起数值不稳定性。本文给出的单点精细积分法导出的显式积分格式可证明是无条件稳定的。就扩散方程与对流─扩散方程作出了本文方法与差分法导出的格式之间的对比。数值例题也表明了单点积分法的优越性。 展开更多
关键词 有限差分法 时间步积分 差分法 偏微分方程
在线阅读 下载PDF
电磁脉冲对电缆耦合问题的理论研究 被引量:16
11
作者 刘顺坤 聂鑫 陈向跃 《电波科学学报》 EI CSCD 北大核心 2010年第2期348-352,共5页
为克服传输线方程求解电磁脉冲电缆耦合问题的不足,提出了有限差分求解电场积分方程的方法,解决了架空电缆电磁脉冲耦合效应的数值模拟难题;引入了基于细线散射的时域有限差分法,解决了地面铺设电缆电磁脉冲耦合效应的数值模拟难题。利... 为克服传输线方程求解电磁脉冲电缆耦合问题的不足,提出了有限差分求解电场积分方程的方法,解决了架空电缆电磁脉冲耦合效应的数值模拟难题;引入了基于细线散射的时域有限差分法,解决了地面铺设电缆电磁脉冲耦合效应的数值模拟难题。利用时域电场积分方程方法研究了长度、线径等电缆参数变化及脉宽、极化方向等馈源参数变化对其电磁脉冲耦合效应的影响。给出了电缆耦合电流的波形特征、沿线分布规律。研究表明:在电磁脉冲作用下,自由放置的电缆,其感应皮电流为衰减振铃波形;电流最大值出现在电缆的中心处;随着电缆长度的增加,电缆上的感应皮电流非线性变大;电缆线径的增加、脉冲宽度的改变对其耦合电流峰值影响较小;电场的极化方向对电缆感应皮电流的影响较大。 展开更多
关键词 电场积分方程 时域有限差分法 电磁脉冲 电缆 耦合
在线阅读 下载PDF
时域电场积分方程在求解架空屏蔽电缆蒙皮电流中的应用 被引量:4
12
作者 刘顺坤 陈雨生 +2 位作者 孙蓓云 周辉 谢彦召 《强激光与粒子束》 EI CAS CSCD 北大核心 2005年第4期637-640,共4页
时域电场积分方程是研究电磁辐射及散射等问题的重要方程。在数值计算中,该方程方法只需将散射体进行剖分,而不必将剖分推至整个计算域内进行,计算效率较高。将时域电场积分方程方法引入到电磁脉冲作用下架高屏蔽电缆蒙皮电流的计算中,... 时域电场积分方程是研究电磁辐射及散射等问题的重要方程。在数值计算中,该方程方法只需将散射体进行剖分,而不必将剖分推至整个计算域内进行,计算效率较高。将时域电场积分方程方法引入到电磁脉冲作用下架高屏蔽电缆蒙皮电流的计算中,研究了电缆的蒙皮感应电流分布及波形特征。通过与辐射波电磁脉冲模拟器的实验结果的比对,证明了数值结果的可靠性。 展开更多
关键词 时域积分方程 有限差分 电磁脉冲 屏蔽电缆
在线阅读 下载PDF
超微光子晶体多模干涉型滤波器的设计与优化 被引量:3
13
作者 谈振兴 于天宝 +2 位作者 方利广 陈淑文 廖清华 《光子学报》 EI CAS CSCD 北大核心 2009年第3期592-596,共5页
根据对称入射时的自映像原理,设计了一种基于光子晶体多模波导的超微滤波器,对于进入该种器件的1.31μm和1.55μm两种通信波长的光信号,能够在通过一种波长的同时,对另一种波长进行阻隔,整个器件的长度不超过15μm.讨论了该滤波器的输... 根据对称入射时的自映像原理,设计了一种基于光子晶体多模波导的超微滤波器,对于进入该种器件的1.31μm和1.55μm两种通信波长的光信号,能够在通过一种波长的同时,对另一种波长进行阻隔,整个器件的长度不超过15μm.讨论了该滤波器的输出效率,并通过调整位于单模波导与多模波导联结处的介电柱位置对结构进行优化.优化后的结构具有高的输出效率和较宽的带宽.时域有限差分法的模拟结果表明该器件是可行的. 展开更多
关键词 集成光学 滤波器 光子晶体多模波导 时域有限差分法
在线阅读 下载PDF
求解三维电磁问题的自适应区域分解FDTD方法 被引量:7
14
作者 张华 洪伟 郝张成 《电波科学学报》 EI CSCD 北大核心 2006年第3期397-402,共6页
将三维电磁场问题的求解区域分解为若干独立的子区域,在各个子区域之间通过连接边界条件,进行自适应检测,当检测到波传播至检测边界上时,相应子区域被激活,并进行传统的FDTD迭代计算,否则该子区域置于休眠状态不参与场值迭代,从而减少... 将三维电磁场问题的求解区域分解为若干独立的子区域,在各个子区域之间通过连接边界条件,进行自适应检测,当检测到波传播至检测边界上时,相应子区域被激活,并进行传统的FDTD迭代计算,否则该子区域置于休眠状态不参与场值迭代,从而减少计算时间。对三维微带互连结构的信号完整性仿真结果验证了这种采用边界检测的自适应区域分解时域有限差分(ADDFDTD)方法在节省计算时间和解决复杂问题方面的有效性。 展开更多
关键词 自适应边界检测 区域分解 时域有限差分法 信号完整性
在线阅读 下载PDF
电性源瞬变电磁法油气藏动态监测模拟分析 被引量:10
15
作者 王新宇 严良俊 毛玉蓉 《石油地球物理勘探》 EI CSCD 北大核心 2022年第2期459-466,I0009,共9页
油气藏开采过程中,储层剩余油气被流体驱替前、后声阻抗变化较小,会造成时移地震监测失效,但驱替前、后储层电阻率变化较大,因而时移电磁法在油气藏动态监测中具有物性方面的优势。通过三维数值模拟,分析了电性源瞬变电磁法对剩余油气... 油气藏开采过程中,储层剩余油气被流体驱替前、后声阻抗变化较小,会造成时移地震监测失效,但驱替前、后储层电阻率变化较大,因而时移电磁法在油气藏动态监测中具有物性方面的优势。通过三维数值模拟,分析了电性源瞬变电磁法对剩余油气藏动态监测的能力。为提高数值模拟精度,基于非结构化网格矢量有限元法,采用二阶后退欧拉法(BDF2)变步长差分格式,实现了电性源瞬变电磁三维正演。对比均匀半空间电磁场解析解及三维复杂模型的数值结果,验证了该方法满足正演精度要求。基于此模拟方法,计算并分析了复杂地质背景油气藏模型的电场相对异常动态监测响应特征。基于涪陵页岩气实际地质资料进行建模,对比压裂前、后电场相对异常。研究结果表明,电性源瞬变电磁法对油气藏动态监测响应效果明显,能满足复杂三维油气藏动态监测的地质要求,应用前景广阔。 展开更多
关键词 瞬变电磁法 时移电磁法 油气藏动态监测 矢量有限元 变步长差分格式
在线阅读 下载PDF
一种多导体传输线瞬态响应时间步积分法 被引量:4
16
作者 史凌峰 成立业 +2 位作者 曹成美 孟辰 蔡成山 《电波科学学报》 EI CSCD 北大核心 2012年第5期990-996,共7页
针对时域有限差分法分析非均匀有损耗多导体传输线瞬态响应的差分振荡问题,提出一种多导体传输线瞬态响应的时间步积分方法。该方法对不等长非均匀有损多导体传输线建模时,无需对耦合传输线进行解耦,能够较方便地求解耦合状态下的不等... 针对时域有限差分法分析非均匀有损耗多导体传输线瞬态响应的差分振荡问题,提出一种多导体传输线瞬态响应的时间步积分方法。该方法对不等长非均匀有损多导体传输线建模时,无需对耦合传输线进行解耦,能够较方便地求解耦合状态下的不等长非均匀有损多导体传输线系统。所提方法首先对电波方程中的空间微分算子进行差分离散,然后采用梯形积分法对时间微分算子进行积分,最终得到分析不等长非均匀有损多导体传输线瞬态响应的时间步积分算法。该算法不仅应用于不等长非均匀有损多导体传输线模型进行数值计算和理论分析,并对传输线终端端接线性负载和非线性负载的两种情况进行仿真验证。仿真结果表明该方法是有效的。 展开更多
关键词 时间步积分 多导体传输线 瞬态分析 时域有限差分
在线阅读 下载PDF
基于慢光效应的解波分复用器的设计与研究 被引量:5
17
作者 吕淑媛 赵辉 《激光技术》 CAS CSCD 北大核心 2015年第5期617-620,共4页
为了设计一种基于正方晶格光子晶体波导的三端口解波分复用器,利用光子晶体的光子带隙特性和光在光子晶体慢光波导中的传输特性,并采用平面波展开法及时域有限差分法对解波分复用器的特性进行了理论分析,数值计算了解复用器的空间和光... 为了设计一种基于正方晶格光子晶体波导的三端口解波分复用器,利用光子晶体的光子带隙特性和光在光子晶体慢光波导中的传输特性,并采用平面波展开法及时域有限差分法对解波分复用器的特性进行了理论分析,数值计算了解复用器的空间和光谱分布特性。结果表明,调整波导中与波导紧邻的第1排柱子的位置,使三部分波导中传输的慢光频率发生改变,通过合理的结构设计,利用慢光模式传输的特点将不同频率的光从不同的输出端输出,可实现解复用的目的。 展开更多
关键词 集成光学 解波分复用器 时域有限差分法 光子晶体慢光波导
在线阅读 下载PDF
腔体内传输线耦合的电磁仿真软件与传输线方程的混合算法 被引量:1
18
作者 叶志红 廖成 +2 位作者 张敏 周海京 李瀚宇 《强激光与粒子束》 EI CAS CSCD 北大核心 2015年第10期68-71,共4页
研究了一种有限积分法软件与传输线方程相结合的混合算法,用于解决复杂电磁环境下屏蔽腔体内传输线的电磁耦合问题。利用有限积分法软件实现屏蔽腔体的建模,仿真得到腔体内部空间电磁场分布,并设置电场探针提取出传输线的激励场。利用... 研究了一种有限积分法软件与传输线方程相结合的混合算法,用于解决复杂电磁环境下屏蔽腔体内传输线的电磁耦合问题。利用有限积分法软件实现屏蔽腔体的建模,仿真得到腔体内部空间电磁场分布,并设置电场探针提取出传输线的激励场。利用传输线方程建立腔体内传输线的耦合模型,将得到的传输线激励场引入到传输线方程作为等效分布电压和电流源。利用时域有限差分(FDTD)格式对传输线方程进行离散,从而迭代求解出传输线终端负载上的电压和电流响应。通过与文献以及传统数值算法的计算结果进行对比,验证混合算法的正确性。研究表明,该混合算法在模拟电大尺寸腔体内传输线的电磁耦合方面,具有较高的精度和计算效率。 展开更多
关键词 有限积分法软件 时域有限差分法 传输线方程 激励场 电压和电流响应
在线阅读 下载PDF
FDTD法分析高速集成电路芯片内互连线 被引量:4
19
作者 袁正宇 李征帆 邹瑉珋 《电子学报》 EI CAS CSCD 北大核心 2000年第2期14-16,共3页
本文首次利用时域有限差分 (FDTD)法分析了高速集成电路芯片内半导体基片上的有耗互连传输线的电特性 .文中提出了有耗吸收边界条件 ,推导了不同媒质交界面上的边界条件通用格式 .在FDTD分析的基础上 ,得到传输线各种参数的频变特性 ,... 本文首次利用时域有限差分 (FDTD)法分析了高速集成电路芯片内半导体基片上的有耗互连传输线的电特性 .文中提出了有耗吸收边界条件 ,推导了不同媒质交界面上的边界条件通用格式 .在FDTD分析的基础上 ,得到传输线各种参数的频变特性 ,为芯片内电路模拟提供了可靠的参数 . 展开更多
关键词 高速集成电路 时域有限差分 互连线
在线阅读 下载PDF
矩形金属机壳孔缝电磁耦合特性研究 被引量:2
20
作者 吴刚 孙靖虎 +1 位作者 李兵 张晓阳 《电波科学学报》 EI CSCD 北大核心 2015年第6期1222-1227,共6页
针对矩形金属机壳上孔缝在高功率电磁环境下的电磁耦合问题,提出了一种使用模式匹配法和基于矩量法求解的混合位积分方程来分析孔缝电磁耦合特性的全波混合算法.该算法考虑了高次模、孔缝形状、孔缝厚度以及入射波极化方向等因素对电磁... 针对矩形金属机壳上孔缝在高功率电磁环境下的电磁耦合问题,提出了一种使用模式匹配法和基于矩量法求解的混合位积分方程来分析孔缝电磁耦合特性的全波混合算法.该算法考虑了高次模、孔缝形状、孔缝厚度以及入射波极化方向等因素对电磁耦合特性的影响.通过将数值仿真结果与经典孔缝电磁耦合模型的测试结果进行对比,验证了该算法具有较高的准确性,与经典的时域有限差分法相比,该算法具有很高的计算效率.研究结果表明:有孔金属机壳在外界强电磁辐射条件下具有明显的谐振效应,在谐振频率点,耦合进入机壳的电场将大大增强,且孔缝附近以及机壳中心的耦合电场峰值要高于其他位置的耦合电场峰值;随着机壳表面上孔缝厚度的增加,耦合进入其内部的电场也在减弱;当机壳上的孔缝为矩形且入射波的电场极化方向平行于矩形缝隙的短边时,对应于该极化方向的孔缝电场耦合强度是所有极化方向中最强的. 展开更多
关键词 高功率微波 电磁耦合 数值计算 金属机壳 矩形孔缝
在线阅读 下载PDF
上一页 1 2 4 下一页 到第
使用帮助 返回顶部