期刊文献+
共找到776篇文章
< 1 2 39 >
每页显示 20 50 100
Analysis and Simulation of the Influence of Electromagnetic Fields on Living Beings near HV Power Lines Using the FDTD Method
1
作者 Anthony Bassesuka Sandoka Nzao 《Open Journal of Applied Sciences》 2023年第12期2343-2359,共17页
Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-freque... Recent decades have seen rapid advances in the field of electrical engineering, such that our environment has become a sea of electrical and magnetic signals, raising questions about the possible effects of low-frequency electromagnetic fields on the environment and which are capable of modifying and destroying our ecosystem. Particular interest was given in this article due to a massive influx of population living near high voltage lines. The analysis and simulation of the influence of low frequency electromagnetic fields on living beings in the vicinity of high voltage sources 132 kV and 220 kV in urban areas in DR Congo is the subject of our research with a view to estimating the level of exposure of humans to low frequency electromagnetic fields. To carry out our research, we used the classic method of analyzing the field produced near a high voltage line based on Maxwell’s image theory, the Maxwell-Gauss theorem and Maxwell-Ampère theorem to model and quantify low-frequency electromagnetic fields in the vicinity of high-voltage lines. The 2D FDTD numerical formulation was developed from telegraphers’ equations and allowed us to obtain models of current and voltage induced by electromagnetic fields on living beings below and near HV lines. The different simulations carried out on the proposed models illustrate the effects of the electrical and geometric parameters of the pylons on the distribution of the electromagnetic field in the vicinity of the HV lines. The results obtained were compared to the safety limits recommended by the standards. 展开更多
关键词 Low Frequency Electromagnetic Field Modeling Field and Health Maxwell’s Equations fdtd method Electromagnetic Compatibility High Voltage Line
下载PDF
Piecewise linear recursive convolution FDTD method for magnetized plasmas 被引量:4
2
作者 Liu Song Zhong Shuangying Liu Shaobin 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 2006年第2期290-295,共6页
The piecewise linear recursive convolution (PLRC) finite-different time-domain (FDTD) method greatly improves accuracy over the original recursive convolution (RC) FDTD approach but retains its speed and efficie... The piecewise linear recursive convolution (PLRC) finite-different time-domain (FDTD) method greatly improves accuracy over the original recursive convolution (RC) FDTD approach but retains its speed and efficiency advantages. A PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time is presented, enabled the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations the reflection and transmission coefficients through a magnetized plasma layer. The results show that the PLRC-FDTD method has significantly improved the accuracy over the original RC method. 展开更多
关键词 electromagnetic wave fdtd methods piecewise linear recursive convolution magnetized plasma.
下载PDF
A DIAGONAL SPLIT-CELL MODEL FOR THE OVERLAPPING YEE FDTD METHOD 被引量:3
3
作者 Jinjie Liu Moysey Brio Jerome V. Moloney 《Acta Mathematica Scientia》 SCIE CSCD 2009年第6期1670-1676,共7页
In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell's equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not req... In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell's equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not require permittivity averaging so that better accuracy can be achieved. Our numerical results on optical force computation show that the standard FDTD method converges linearly, while the proposed method achieves quadratic convergence and better accuracy. 展开更多
关键词 fdtd method Overlapping Yee method Maxwell's equations Optical force
下载PDF
A Revised Piecewise Linear Recursive Convolution FDTD Method for Magnetized Plasmas 被引量:1
4
作者 刘崧 钟双英 刘少斌 《Plasma Science and Technology》 SCIE EI CAS CSCD 2005年第6期3122-3126,共5页
The piecewise linear recursive convolution (PLRC) finite-different time-domain (FDTD) method improves accuracy over the original recursive convolution (RC) FDTD approach and current density convolution (JEC) b... The piecewise linear recursive convolution (PLRC) finite-different time-domain (FDTD) method improves accuracy over the original recursive convolution (RC) FDTD approach and current density convolution (JEC) but retains their advantages in speed and efficiency. This paper describes a revised piecewise linear recursive convolution PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time, enabling the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations of the reflection and transmission coefficients through a magnetized plasma layer. The results show that the revised PLRC-FDTD method has improved the accuracy over the original RC FDTD method and JEC FDTD method. 展开更多
关键词 Electromagnetic wave finite-different time-domain fdtd methods piecewise linear recursive convolution magnetized plasma
下载PDF
A Generalized FDTD Method with Absorbing Boundary Condition for Solving a Time-Dependent Linear Schrodinger Equation 被引量:1
5
作者 Frederick Ira Moxley III Fei Zhu Weizhong Dai 《American Journal of Computational Mathematics》 2012年第3期163-172,共10页
The Finite-Difference Time-Domain (FDTD) method is a well-known technique for the analysis of quantum devices. It solves a discretized Schrodinger equation in an iterative process. However, the method provides only a ... The Finite-Difference Time-Domain (FDTD) method is a well-known technique for the analysis of quantum devices. It solves a discretized Schrodinger equation in an iterative process. However, the method provides only a second-order accurate numerical solution and requires that the spatial grid size and time step should satisfy a very restricted condition in order to prevent the numerical solution from diverging. In this article, we present a generalized FDTD method with absorbing boundary condition for solving the one-dimensional (1D) time-dependent Schr?dinger equation and obtain a more relaxed condition for stability. The generalized FDTD scheme is tested by simulating a particle moving in free space and then hitting an energy potential. Numerical results coincide with those obtained based on the theoretical analysis. 展开更多
关键词 SCHRODINGER Equation Absorbing BOUNDARY fdtd method Stability
下载PDF
Theoretical Proof of Unconditional Stability of the 3-D ADI-FDTD Method 被引量:3
6
作者 WANG Yingjun WANG Bingzhong SHAO Wei (School of Physical Electronics,UESTC Chengdu 610054 China) 《Journal of Electronic Science and Technology of China》 2003年第1期1-5,共5页
In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is i... In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint. 展开更多
关键词 alternating-direction implicit(ADI)technique Courant-Friedrich-Levy(CFL)condition restraint finite-difference time-domain(fdtd)method stability
下载PDF
基于ADI-FDTD法的探地雷达正演数值模拟及验证
7
作者 曾良生 迟耀丹 +1 位作者 周庆才 衣淼 《吉林建筑大学学报》 CAS 2024年第1期48-55,共8页
探地雷达(GPR)是一种应用前景广泛、用于探测和定位地下物体的浅层地球物理方法。通过开展探地雷达正演模拟研究,可以获得复杂地下结构的探地雷达图像回波特征。时域有限差分(FDTD)法受到稳定性和收敛性条件的限制,导致效率和精度低。... 探地雷达(GPR)是一种应用前景广泛、用于探测和定位地下物体的浅层地球物理方法。通过开展探地雷达正演模拟研究,可以获得复杂地下结构的探地雷达图像回波特征。时域有限差分(FDTD)法受到稳定性和收敛性条件的限制,导致效率和精度低。交替方向隐式-时域有限差分(ADI-FDTD)法克服FDTD法的稳定性限制,可以选择更大的时间步长来提高计算效率。从原理出发,对ADI-FDTD法进行公式推导,在基于ADI-FDTD法的基础上进行卷积完全匹配层(CPML)的结合。通过进行ADI-FDTD法的时间稳定性仿真和不同材质、不同填充介质的双管,沙槽正演模拟实验,结果表明,ADI-FDTD法可以保证时间无条件稳定性,相同情况下可采用更大的时间步长进行正演模拟,提高正演效率。对不同管线、不同材质情况下的正演剖面曲线特征进行解译分析,最后与沙槽实测进行对比,证明得到的探测解译结果与实际状况达到较好的吻合。 展开更多
关键词 探地雷达(GPR) 正演模拟 无条件稳定 卷积完全匹配层(CPML) 交替方向隐式-时域有限差分(ADI-fdtd)法
下载PDF
基于稀疏采样的FDTD/TDPO混合优化算法
8
作者 王林茜 陈娟 牟春晖 《电波科学学报》 CSCD 北大核心 2024年第5期846-851,共6页
基于时域有限差分(finite-difference time-domain,FDTD)和时域物理光学(time-domain physical optics,TDPO)的混合算法提出用稀疏采样优化近场外推的方法,解决了大口径反射面天线电磁计算用时过长、计算效率低的问题。传统混合算法中... 基于时域有限差分(finite-difference time-domain,FDTD)和时域物理光学(time-domain physical optics,TDPO)的混合算法提出用稀疏采样优化近场外推的方法,解决了大口径反射面天线电磁计算用时过长、计算效率低的问题。传统混合算法中馈源的近场外推到反射面天线的步骤计算量较大,稀疏采样法从空间和时间两个维度进行间隔采样,减少了此步骤的计算量,从而节省了用时。给出了馈源为喇叭天线、抛物面口径为100个波长的抛物面天线算例,所得计算结果与传统方法相比最大相对误差仅有-36.03 dB,用时最快可缩短为原来的2.51%,与CST软件计算结果对比吻合度较高,用时缩减了93.18%,验证了本文优化算法可在不影响计算精度的前提下提高计算效率。 展开更多
关键词 时域有限差分(fdtd) 时域物理光学(TDPO) 混合方法 稀疏采样 抛物面天线 远区辐射场
下载PDF
FDTD Method in the Analysis of Lithium Niobate Modulator
9
作者 Jianying Zhou, Xihua Li, Xiaoping Zhou, Kejian Chen, Xiaoqin Jiang, Minghua Wang(Department of Information Science and Electronic Engineering, Zhejiang University, Hangzhou, P. R. China, Tel: 86-571-87951706, Fax: 86-571-87952867, E-mail:phd16@isee.zju.edu.cn) 《光学学报》 EI CAS CSCD 北大核心 2003年第S1期191-192,共2页
The Mach-Zehnder traveling-wave lithium niobate modulator is simulated by time domain finite difference (FDTD) method for different device geometry. The result is comparable to the reported value by finite element met... The Mach-Zehnder traveling-wave lithium niobate modulator is simulated by time domain finite difference (FDTD) method for different device geometry. The result is comparable to the reported value by finite element method. 展开更多
关键词 fdtd method in the Analysis of Lithium Niobate Modulator IS ET on of in
原文传递
Coupling Effect Analysis of Multiwaveguides by FDTD Method
10
作者 ZHOUXaojun YUZhiuan 《Semiconductor Photonics and Technology》 CAS 1998年第4期219-224,共6页
By use of finite-difference time-domain (FDTD) method,an eigenmode analysis in a multiwaveguide structure is presented.Because of difference in propagation constants of different modes,coupling effect is discussed for... By use of finite-difference time-domain (FDTD) method,an eigenmode analysis in a multiwaveguide structure is presented.Because of difference in propagation constants of different modes,coupling effect is discussed for three and five waveguide systems.The field distribution in multiwaveguides is given. 展开更多
关键词 Coupling Effect fdtd method Multiwaveguides
下载PDF
有限深度海域海面电磁散射的FDTD方法研究 被引量:2
11
作者 张童 任新成 +2 位作者 朱小敏 杨鹏举 赵晔 《微波学报》 CSCD 北大核心 2023年第2期88-94,共7页
雷达能够探测、遥感有限深度海域的海面及附近目标,针对这一需求,采用改进的文氏海谱和Monte Carlo方法模拟有限深度海域海面,选取Debye模型计算海水介电常数,运用时域有限差分法研究了有限深度海域海面电磁散射,仿真了散射系数的角分... 雷达能够探测、遥感有限深度海域的海面及附近目标,针对这一需求,采用改进的文氏海谱和Monte Carlo方法模拟有限深度海域海面,选取Debye模型计算海水介电常数,运用时域有限差分法研究了有限深度海域海面电磁散射,仿真了散射系数的角分布曲线。仿真结果表明:散射系数随散射角振荡变化,在镜反射方向处产生散射增强效应;海水深度虽然对散射系数影响较小,但是仍有规律可循,海水深度越大,散射系数越小;海面风速、入射角、风区范围对散射系数影响较大,风速越大,散射系数越小;入射角越大,散射系数越小;风区范围越大,散射系数越小。 展开更多
关键词 电磁散射 时域有限差分法 有限深度海域海面 蒙特卡罗方法 散射系数
原文传递
面向电大多尺度模型的大规模并行对称半正定亚网格FDTD算法
12
作者 陈晓洁 鲍献丰 +2 位作者 李瀚宇 张爱清 周海京 《北京理工大学学报》 EI CAS CSCD 北大核心 2023年第6期649-656,共8页
针对计算电磁学算法应用于实际工程问题时面临的计算模型电大尺寸多尺度瓶颈,基于自研的高性能并行计算框架JASMIN,通过提出融合h-自适应时间积分算法流程的对称半正定亚网格FDTD算法流程,实现了具备大规模并行能力的亚网格时域有限差分... 针对计算电磁学算法应用于实际工程问题时面临的计算模型电大尺寸多尺度瓶颈,基于自研的高性能并行计算框架JASMIN,通过提出融合h-自适应时间积分算法流程的对称半正定亚网格FDTD算法流程,实现了具备大规模并行能力的亚网格时域有限差分(FDTD)算法.该算法保留了FDTD算法高并行可扩展性优势,同时可局部使用细化网格剖分电大尺寸模型精细结构部位,支持粗网格区域和细网格区域分别采用大时间步长和小时间步长进行迭代计算.算例表明,该方法在解决电大多尺度电磁模拟问题时,具有较强的并行可扩展性,同时能够保证数值计算结果的精度和稳定性,与传统的FDTD算法相比,可显著降低内存需求,提高计算效率. 展开更多
关键词 电大多尺度模型 对称半正定亚网格fdtd算法 大规模并行
下载PDF
强电磁脉冲电离大气三维FDTD仿真
13
作者 肖科 丁亮 +2 位作者 柴舜连 伍小灿 尹文禄 《微波学报》 CSCD 北大核心 2023年第4期1-5,共5页
强电磁脉冲在大气中传播时会导致大气分子电离,且电离气体又会反作用于传播的电磁波,因此是一个传播状态动态实时变化的物理过程。要准确模拟该过程需要综合考虑电磁波传播、大气电离、带电粒子运动等物理过程。为解决该问题,文中结合... 强电磁脉冲在大气中传播时会导致大气分子电离,且电离气体又会反作用于传播的电磁波,因此是一个传播状态动态实时变化的物理过程。要准确模拟该过程需要综合考虑电磁波传播、大气电离、带电粒子运动等物理过程。为解决该问题,文中结合麦克斯韦方程、等离子体方程、磁流体动力学方程等,构建多物理场分析方程组,该方程组能够考虑强电磁脉冲传播时大气电离、扩散、附着、复合等物理效应。进而结合三维时域有限差分方法(FDTD)来构建动态仿真模型,对强电磁脉冲在大气中的传播过程进行仿真并分析其宏观大气电离效果。 展开更多
关键词 高功率微波 大气 等离子体 时域有限差分法
原文传递
FDTD Modeling of Lorentzian DNG Metamaterials by Auxiliary Differential Equation Method
14
作者 Chiranjib Goswami Saptarshi Mukherjee +2 位作者 Subrata Karmakar Manimala Pal Rowdra Ghatak 《Journal of Electromagnetic Analysis and Applications》 2014年第5期106-114,共9页
In this paper, Finite Difference Time Domain (FDTD) is utilized to simulate metamaterials of Double Negative (DNG) origin that refers to those materials having simultaneous negative permittivity and permeability. The ... In this paper, Finite Difference Time Domain (FDTD) is utilized to simulate metamaterials of Double Negative (DNG) origin that refers to those materials having simultaneous negative permittivity and permeability. The problem regarding space formulation is achieved by means of auxiliary differential equation method (ADE), which is easy, reliable and also causal process in nature thus making it proficient. It uses fair approximations to explicate the model. Mur’s boundary condition is used for 1-D problem space and convolution perfectly matched layer boundary is implemented for 2-D problem space. The properties of metamaterial conform their speculations of energy absorption, enhancement and backward propagation property with the aid of graphs engineered by Matlab simulation both in 1-D and 2-D. Also, the interaction of fields on DNG and Double Positive (DPS) layers is contrasted. The results achieved elucidate the validity and effectiveness of the ADE method and the Convolution Perfectly Match Layer (CPML) in designing DNG metamaterials. 展开更多
关键词 ABC ADE method CPML DNG fdtd method
下载PDF
Analysis of Complex Electromagnetic Structures by Hybrid FDTD/WCIP Method
15
作者 Gharbi Ramzi Zairi Hassen +1 位作者 Trabelsi Hichem Baudrand Henri 《Journal of Electromagnetic Analysis and Applications》 2012年第12期497-503,共7页
This paper proposes a hybrid full-wave analysis using Finite-Difference Time-Domain (FDTD) and Wave Concept Iterative Process (WCIP) methods, developed to analyze locally arbitrarily shaped microwave structures and Mu... This paper proposes a hybrid full-wave analysis using Finite-Difference Time-Domain (FDTD) and Wave Concept Iterative Process (WCIP) methods, developed to analyze locally arbitrarily shaped microwave structures and Multilayer Planar structure. Using the equivalence principle, the original problem can be decomposed into two sub regions and solve each sub region separately. An interpolation scheme is proposed for communicating between the FDTD fields and WCIP wave, which will not require the effort of fitting the WCIP mesh to the FDTD cells in the interface region. This method is applied to calculate the scattering parameters of arbitrary (3-D) microwave structures. Applying FDTD to 3D discontinuity and WCIP to the remaining region preserves the advantages of both WCIP flexibility and FDTD efficiency. A comparison of the results with the FDTD staircasing data verifies the accuracy of the proposed method. 展开更多
关键词 fdtd HYBRID FINITE-DIFFERENCE TIME-DOMAIN (fdtd) HYBRID Techniques WCIP method
下载PDF
等离子体散射FDTD分析的移位算子方法 被引量:44
16
作者 葛德彪 吴跃丽 朱湘琴 《电波科学学报》 EI CSCD 2003年第4期359-362,399,共5页
色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系,这就给应用FDTD方法计算色散介质中波的散射和传播带来困难。文中把色散介质的相对个电系数写成以jω为自变量的有理分式函数,用/t代替jω,过渡到时域,再引入离散时域移... 色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系,这就给应用FDTD方法计算色散介质中波的散射和传播带来困难。文中把色散介质的相对个电系数写成以jω为自变量的有理分式函数,用/t代替jω,过渡到时域,再引入离散时域移位算子代替时间微分算子来处理有理分式函数形式的介电系数,进而导出FDTD中电位移矢量D和电场强度E之间的关系。所得结果与Z变换理论的梯形近似结果一致,但与Z变换理论相比较,移位算子方法推导简单,概念简明。文中用移位算子结合FDTD方法分别计算一维等离子体平板的反射系数和二维等离子体柱的远场散射。所得结果与用Z变换理论结合FDTD方法计算的结果相符合。 展开更多
关键词 色散介质 fdtd 介电系数 移位算子 等离子体 远场散射
下载PDF
直角坐标系下非均匀FDTD网格生成系统 被引量:8
17
作者 周国祥 程萍 +1 位作者 蒋经国 杨明武 《微波学报》 CSCD 北大核心 2005年第2期56-59,共4页
基于非均匀FDTD网格的电磁场模拟是一种解决诸多实际电磁场问题的通用、有效方法。本文介绍了在直角坐标系下非均匀FDTD网格图形生成系统,可以快速、高效地产生1D、2D、3DFDTD网格,展现物体的几何图形,还可以动态地、可视化地得到检查... 基于非均匀FDTD网格的电磁场模拟是一种解决诸多实际电磁场问题的通用、有效方法。本文介绍了在直角坐标系下非均匀FDTD网格图形生成系统,可以快速、高效地产生1D、2D、3DFDTD网格,展现物体的几何图形,还可以动态地、可视化地得到检查、调节和修改。作为检验,本文以一实际结构为例,产生并显示了其网格图形。 展开更多
关键词 直角坐标系 fdtd 非均匀 网格 图形生成系统 电磁场模拟 磁场问题 有效方法 几何图形 可视化
下载PDF
ESD辐射场数值建模与FDTD分析 被引量:9
18
作者 毕增军 刘尚合 冯德玉 《高电压技术》 EI CAS CSCD 北大核心 2003年第8期26-28,共3页
采用 FDTD法研究了 ESD电流注入到细导线中所产生的辐射场 ,建立了数值模型 ,并计算和分析了产生的辐射场特征及导线半径和导线长度对辐射场的影响 ,结果表明 :在同等条件下 ,ESD电流通过细导线产生的辐射场比ESD火花产生的辐射场强的多... 采用 FDTD法研究了 ESD电流注入到细导线中所产生的辐射场 ,建立了数值模型 ,并计算和分析了产生的辐射场特征及导线半径和导线长度对辐射场的影响 ,结果表明 :在同等条件下 ,ESD电流通过细导线产生的辐射场比ESD火花产生的辐射场强的多 ,FDTD法与 Mo 展开更多
关键词 ESD 辐射场 数值建模 fdtd 分析 静电放电
下载PDF
导弹目标的FDTD建模与RCS计算 被引量:12
19
作者 郑奎松 葛德彪 魏兵 《系统工程与电子技术》 EI CSCD 北大核心 2004年第7期896-899,共4页
导弹弹头对后向散射有明显的影响。FDTD建模时 ,为了更精确的模拟弹头的外形轮廓 ,采有超椭球(Superspheriod)几何体来模拟导弹弹头的雷达罩。通过调整超椭球体方程中ν因子值的大小能够得到不同形状弹头的外形轮廓 ,其中当ν取 1 381... 导弹弹头对后向散射有明显的影响。FDTD建模时 ,为了更精确的模拟弹头的外形轮廓 ,采有超椭球(Superspheriod)几何体来模拟导弹弹头的雷达罩。通过调整超椭球体方程中ν因子值的大小能够得到不同形状弹头的外形轮廓 ,其中当ν取 1 381时可近似为传统的VonKarman雷达罩。该超椭球方程推导简单且能够模拟大多数传统雷达罩的形状。最后 ,给出了用FDTD方法计算两种导弹弹头模型的后向散射和侧向散射。结果表明 ,采用一定参数下的超椭球几何体弹头比球冠状弹头有效的减小后面RCS。 展开更多
关键词 电磁散射 雷达散射截面 fdtd方法 超椭球体 导弹
下载PDF
FDTD求解高功率微波大气传播问题的可行性研究 被引量:13
20
作者 唐涛 廖成 杨丹 《电波科学学报》 EI CSCD 北大核心 2010年第1期122-126,共5页
提出了应用时域有限差分(FDTD)方法对高功率微波(HPM)大气传播模型进行数值模拟。HPM大气传播模型是由Maxwell方程组和电子磁流体方程共同描述的。利用FDTD的特点,在适当的时间和空间位置对该方程组进行中心差分,获得其显式差分格式。... 提出了应用时域有限差分(FDTD)方法对高功率微波(HPM)大气传播模型进行数值模拟。HPM大气传播模型是由Maxwell方程组和电子磁流体方程共同描述的。利用FDTD的特点,在适当的时间和空间位置对该方程组进行中心差分,获得其显式差分格式。并结合流体力学有限元分析中的边界处理方法,建立了一套全新、高效的求解HPM大气传播的数值方法,其计算量仅与空间网格步数相当,在普通PC机上就能完成相应计算。计算结果验证了该方法的合理性和正确性。 展开更多
关键词 高功率微波(HPM) 大气传播 时域有限差分法(fdtd)
下载PDF
上一页 1 2 39 下一页 到第
使用帮助 返回顶部