期刊文献+
共找到665篇文章
< 1 2 34 >
每页显示 20 50 100
高分辨率有限差分-有限元混合方法及其在气动热计算中的应用 被引量:12
1
作者 段占元 童秉纲 姜贵庆 《空气动力学学报》 CSCD 北大核心 1997年第4期462-467,共6页
在文献[1]结合NND格式思想[2]提出的有限元格式基础上,给出一种有限差分-有限元混合方法。通过求解完全Navier-Stokes方程,得到了高超声速情况下钝头体轴对称绕流的较满意结果,同时讨论了网络Re数对驻点热流的影响。
关键词 高超声速流 热流 有限差分 有限元 气动热计算
在线阅读 下载PDF
从求学西北工业大学到参与创立有限元方法——崔俊芝院士访谈录
2
作者 崔俊芝(讲述) 王涛 魏蕾(访问整理) 《自然科学史研究》 北大核心 2025年第1期105-119,共15页
有限元方法是科学与工程计算领域最重要的算法之一,由欧美工程师与中国数学家沿着不同的学术路径独立地创立。崔俊芝院士早年求学于西北工业大学,毕业后进入中国科学院计算技术研究所工作,在冯康的指导下从事水坝计算的研究。在访谈中,... 有限元方法是科学与工程计算领域最重要的算法之一,由欧美工程师与中国数学家沿着不同的学术路径独立地创立。崔俊芝院士早年求学于西北工业大学,毕业后进入中国科学院计算技术研究所工作,在冯康的指导下从事水坝计算的研究。在访谈中,他回顾了自己的教育背景和参与创立有限元方法的详细历程;对中外有限元方法的发明与创造进行了对比分析,指出同期中国团队创立的有限元方法更具一般性;总结了有限元软件在中国的发展概况与经验教训。 展开更多
关键词 有限元方法 崔俊芝 冯康 水坝计算 变分—差分方法 直接刚度法
在线阅读 下载PDF
不同加载方法对甲板大开口分段有限元计算的误差分析 被引量:3
3
作者 石勇 侯海量 +2 位作者 白雪飞 朱锡 冯文山 《海军工程大学学报》 CAS 2002年第1期90-93,共4页
利用MSC-NASTRAN有限元程序模拟计算甲板大开口的应力.建立简单的总体和分段力学模型,采用不同的加载方法进行甲板应力计算,将局部应力计算结果与整体应力计算结果比较,得出各种加载方法有限元计算的误差,从而给出误差... 利用MSC-NASTRAN有限元程序模拟计算甲板大开口的应力.建立简单的总体和分段力学模型,采用不同的加载方法进行甲板应力计算,将局部应力计算结果与整体应力计算结果比较,得出各种加载方法有限元计算的误差,从而给出误差小、建模简单的局部分段模型加载的最佳方法. 展开更多
关键词 加载方法 甲板 大开口分段 有限元计算 差分 应力分析 局部分段加载方法 舰船设计
在线阅读 下载PDF
Navier-Stokes方程有限元并行计算方法最新进展 被引量:5
4
作者 尚月强 黄淑梅 《贵州师范大学学报(自然科学版)》 CAS 2012年第2期113-120,共8页
Navier-Stokes方程是流体力学的基本方程,其并行数值求解方法是当前计算数学和计算流体力学领域的最前沿课题之一。综述了Navier-Stokes方程有限元并行计算方法的研究现状。将已有的方法进行分类,分别介绍了其基本思想,评述了各种有限... Navier-Stokes方程是流体力学的基本方程,其并行数值求解方法是当前计算数学和计算流体力学领域的最前沿课题之一。综述了Navier-Stokes方程有限元并行计算方法的研究现状。将已有的方法进行分类,分别介绍了其基本思想,评述了各种有限元并行计算方法的优缺点,讨论了有限元并行计算方法所面临的问题,并对其发展趋势进行了展望。 展开更多
关键词 Navier—Stokes方程 有限元方法 并行计算 计算流体力学
在线阅读 下载PDF
基于有限差分-谱方法的分数阶Burgers流体的非稳态驻点流动
5
作者 白羽 王欣 +1 位作者 张艳 刘春燕 《计算力学学报》 CAS CSCD 北大核心 2024年第3期458-466,共9页
研究了分数阶Burgers流体通过拉伸平板的非稳态驻点流动问题。将分数阶导数引入Burgers流体模型可以更好地模拟流动过程,但也增加了模型的复杂性和求解难度。首次运用有限差分-谱方法求解分数阶Burgers流体模型,离散格式构造简单有效。... 研究了分数阶Burgers流体通过拉伸平板的非稳态驻点流动问题。将分数阶导数引入Burgers流体模型可以更好地模拟流动过程,但也增加了模型的复杂性和求解难度。首次运用有限差分-谱方法求解分数阶Burgers流体模型,离散格式构造简单有效。采用谱方法对控制方程中的空间项进行离散,利用有限差分方法分别结合L-1和L-2算法离散控制方程中的时间项,给出了两种离散格式,并且通过构造数值算例证明了离散格式的收敛性。结果表明,在靠近平板处,速度随着分数阶导数的增加而减小,而无穷远处的流体速度呈现出相反的趋势,体现了分数阶导数的记忆特性。此外,雷诺数越小,流体的粘度越大,导致流体速度越大。由于松弛时间参数的松弛特性,靠近平板处松弛时间参数对速度分布有抑制作用,远离平板处松弛时间促进流体流动。 展开更多
关键词 非稳态 驻点流 分数阶Burgers流体 L-2算法 有限差分-方法
在线阅读 下载PDF
复杂地质条件的间断有限元地震波数值模拟及GPU加速
6
作者 韩德超 刘卫华 +2 位作者 张春丽 袁媛 白鹏 《石油物探》 北大核心 2025年第4期639-652,共14页
间断Galerkin有限元方法(DGFEM)是一种具有较高模拟精度的有限元方法,但其算法编程难度大,其针对各类复杂介质的波动方程的算法目前未见统一的计算格式。为此,基于三角形非结构化网格以及局部Lax-Friedrichs数值流,构建了针对复杂介质... 间断Galerkin有限元方法(DGFEM)是一种具有较高模拟精度的有限元方法,但其算法编程难度大,其针对各类复杂介质的波动方程的算法目前未见统一的计算格式。为此,基于三角形非结构化网格以及局部Lax-Friedrichs数值流,构建了针对复杂介质波动方程模拟的DGFEM编程计算矩阵,并进一步得出了适用于各类复杂介质模拟的单一波场分量的通用计算格式。该通用计算格式能够有效提升DGFEM算法编程的可拓展性。基于该格式给出了DGFEM的通用CUDA核函数的构建方法,并形成CPU+GPU的二维DGFEM并行计算程序框架。通用CUDA核函数可以将DGFEM算法进一步延伸到更加复杂的介质以及三维情况。理论模型和复杂山地模型的数值实验结果表明,构建的通用计算格式和CUDA核函数可以准确模拟声波、弹性波、粘弹性波、孔隙弹性波方程描述的纵波、横波以及慢纵波等波现象。相比单核CPU模拟,二维DGFEM弹性波GPU计算加速比平均在100倍左右。同时,弹性波、粘弹性波、孔隙弹性波模拟耗时约为声波模拟的1.7,2.3,3.0倍。此结果可以指导复杂介质耦合条件模拟时多进程的负载平衡。 展开更多
关键词 间断Galerkin有限元方法 弹性波 粘弹性波 孔隙弹性波 数值模拟 GPU并行计算.
在线阅读 下载PDF
三角网格迎风差分快速扫描走时计算方法
7
作者 张大雨 黄光南 +2 位作者 曾爱平 杜贤军 李红星 《石油地球物理勘探》 北大核心 2025年第1期109-118,共10页
快速扫描法是一种基于求解程函方程的走时计算方法,通过对网格进行四个方向的扫描,可以快速计算出网格节点的走时。根据局部计算单元所使用的节点数量可以分为两种规则网格类型:五点三角网格和九点三角网格。不同网格快速扫描法计算走... 快速扫描法是一种基于求解程函方程的走时计算方法,通过对网格进行四个方向的扫描,可以快速计算出网格节点的走时。根据局部计算单元所使用的节点数量可以分为两种规则网格类型:五点三角网格和九点三角网格。不同网格快速扫描法计算走时的精度存在差异,文中对两种网格的快速扫描法走时计算进行分析对比。首先,根据不同的网格单元分别建立走时与慢度的关系式,将程函方程离散化并推导出迎风差分格式的二次方程,通过对方程进行数值求解即可得到可能的走时值;然后,讨论了两种网格快速扫描法的走时计算步骤,给出了网格边界的处理方法和九点三角网格走时计算的因果关系条件;最后,通过不同速度模型(均匀速度模型、等梯度速度模型和Sigsbee2a速度模型)进行测试,结果表明九点三角网格能有效提高快速扫描法走时计算精度。 展开更多
关键词 走时计算 因果关系条件 程函方程 有限差分 高斯-赛德尔迭代
在线阅读 下载PDF
计算射孔井水力压裂破裂压力的有限元方法 被引量:31
8
作者 胡永全 赵金洲 +2 位作者 曾庆坤 徐树汉 谢慧华 《天然气工业》 EI CAS CSCD 北大核心 2003年第2期58-59,共2页
水力压裂几乎都是在套管井中进行的 ,射孔参数必然要影响到压裂施工参数设计 ,甚至影响到压裂效果。文章将套管和岩石视为具有不同性能的两种材料 ,应用引进的CAD/CAE/CAMPro/E有限元软件在SGI工作站上计算得到了射孔孔眼周围的应力分... 水力压裂几乎都是在套管井中进行的 ,射孔参数必然要影响到压裂施工参数设计 ,甚至影响到压裂效果。文章将套管和岩石视为具有不同性能的两种材料 ,应用引进的CAD/CAE/CAMPro/E有限元软件在SGI工作站上计算得到了射孔孔眼周围的应力分布。根据主应力方位可以确定压裂裂缝方位 ,结合岩石破裂准则可以得到水力压裂施工破裂压力。 展开更多
关键词 计算 射孔井 水力压裂 油气井 破裂压力 有限元方法 应力分析
在线阅读 下载PDF
冲击-接触问题有限元仿真的并行计算 被引量:25
9
作者 亓文果 金先龙 张晓云 《振动与冲击》 EI CSCD 北大核心 2006年第4期68-72,共5页
冲击-接触问题广泛存在于汽车碰撞等的模拟计算中。简单介绍了求解该类问题的显式有限元方法,对显式有限元方法的并行性进行了讨论。根据显式有限元和冲击-接触问题的计算特点,设计并实现了接触均衡的分区算法。算例计算结果表明:该并... 冲击-接触问题广泛存在于汽车碰撞等的模拟计算中。简单介绍了求解该类问题的显式有限元方法,对显式有限元方法的并行性进行了讨论。根据显式有限元和冲击-接触问题的计算特点,设计并实现了接触均衡的分区算法。算例计算结果表明:该并行算法具有较好的加速比和并行效率。 展开更多
关键词 冲击-接触 有限元 并行计算 区域分割
在线阅读 下载PDF
GPU通用计算平台上中心差分格式显式有限元并行计算 被引量:11
10
作者 蔡勇 李光耀 王琥 《计算机研究与发展》 EI CSCD 北大核心 2013年第2期412-419,共8页
显式有限元是解决平面非线性动态问题的有效方法.由于显式有限元算法的条件稳定性,对于大规模的有限元问题的求解需要很长的计算时间.图形处理器(GPU)作为一种高度并行化的通用计算处理器,可以很好解决大规模科学计算的速度问题.统一计... 显式有限元是解决平面非线性动态问题的有效方法.由于显式有限元算法的条件稳定性,对于大规模的有限元问题的求解需要很长的计算时间.图形处理器(GPU)作为一种高度并行化的通用计算处理器,可以很好解决大规模科学计算的速度问题.统一计算架构(CUDA)为实现GPU通用计算提供了高效、简便的方法.因此,建立了基于GPU通用计算平台的中心差分格式的显式有限元并行计算方法.该方法针对GPU计算的特点,对串行算法的流程进行了优化和调整,通过采用线程与单元或节点的一一映射策略,实现了迭代过程的完全并行化.通过数值算例表明,在保证计算精度一致的前提下,采用NVIDIA GTX460显卡,该方法能够大幅度提高计算效率,是求解平面非线性动态问题的一种高效简便的数值计算方法. 展开更多
关键词 图形处理器 CUDA 中心差分 显式有限元 并行计算
在线阅读 下载PDF
基于EBE方法的三维有限元并行计算 被引量:6
11
作者 刘耀儒 周维垣 +1 位作者 陈欣 黄岩松 《岩土力学》 EI CAS CSCD 北大核心 2003年第S2期69-72,共4页
在水利工程中,施工过程的模拟、动力的时域分析、开裂计算等,都对大规模并行计算提出了迫切的需求。然而,基于高斯消去的有限元直接解法,通常会占用大量的内存,并花费大量的CPU时间。而水利工程中的问题多为大带宽问题,这些问题更为突... 在水利工程中,施工过程的模拟、动力的时域分析、开裂计算等,都对大规模并行计算提出了迫切的需求。然而,基于高斯消去的有限元直接解法,通常会占用大量的内存,并花费大量的CPU时间。而水利工程中的问题多为大带宽问题,这些问题更为突出。基于EBE-PCG方法的有限元方法,可以避免形成整体刚度矩阵,进而,显著减少内存的需求。而且,这种方法可以有效地并行实现,为大规模数值计算提供了可能。采用基于EBE策略的Jacobi预处理共轭梯度法,编制了有限元计算程序,并成功应用于溪洛渡、锦屏等工程的大规模数值分析。结果表明,对水利工程中的大带宽问题,该方法是一种很有效的并行计算方法。 展开更多
关键词 有限元 并行计算 EBE方法 水利工程 应用
在线阅读 下载PDF
基于有限元方法的大型永磁直驱同步风力发电机电磁场计算 被引量:21
12
作者 何山 王维庆 +1 位作者 张新燕 赵祥 《电网技术》 EI CSCD 北大核心 2010年第3期157-161,共5页
用有限元软件和电磁场知识对1.2MW永磁风力发电机建模,计算了发电机的静、瞬态磁场。研究了发电机静态气隙磁场的分布,通过傅里叶分解,得出了气隙磁密基波和高次谐波。采用转子斜极方法消除高次谐波,得到了转子斜极改善磁密分布的实际... 用有限元软件和电磁场知识对1.2MW永磁风力发电机建模,计算了发电机的静、瞬态磁场。研究了发电机静态气隙磁场的分布,通过傅里叶分解,得出了气隙磁密基波和高次谐波。采用转子斜极方法消除高次谐波,得到了转子斜极改善磁密分布的实际效果。对于瞬态磁场,采用电路分析方法,根据有限元方法的电路耦合功能,得出了发电机磁场和电势波形。结合实验数据,对比了静、瞬态磁场计算结果,分析了数据产生偏差的原因,验证了理论分析的正确性和有效性。 展开更多
关键词 有限元方法 永磁直驱同步风力发电机 电磁场 计算
在线阅读 下载PDF
流体-结构耦合问题的有限元并行计算研究 被引量:8
13
作者 李政 金先龙 亓文果 《计算力学学报》 EI CAS CSCD 北大核心 2007年第6期727-732,共6页
流体-结构耦合问题广泛存在于各种工程领域,本文采用ALE显式有限元法求解该类问题,并对该方法的并行性进行讨论。同时根据流体-结构耦合问题与ALE显式有限元的计算特点,在坐标递归分区方法的基础上设计并程序实现了基于流体-结构耦合均... 流体-结构耦合问题广泛存在于各种工程领域,本文采用ALE显式有限元法求解该类问题,并对该方法的并行性进行讨论。同时根据流体-结构耦合问题与ALE显式有限元的计算特点,在坐标递归分区方法的基础上设计并程序实现了基于流体-结构耦合均衡的分区算法。通过与坐标递归分区方法的计算结果相比较,对于流体-结构耦合问题的求解,耦合均衡并行分区方法具有更好的加速比和并行效率。 展开更多
关键词 流体-结构耦合 ALE有限元 显式计算 并行计算 区域分解
在线阅读 下载PDF
任意单元间断Galerkin有限元计算方法研究 被引量:16
14
作者 贺立新 张来平 张涵信 《空气动力学学报》 EI CSCD 北大核心 2007年第2期157-162,共6页
基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构... 基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构造方法,提出了非正交三棱柱单元限制器方法;利用上述方法数值模拟了球头和双椭球的高超声速粘性绕流,得到了清晰的流场结构,并得到了较好的压力和热流分布,表明该方法在复杂高超声速流动的数值模拟方面具有广阔的应用前景。 展开更多
关键词 非正交单元 间断Glerkin有限元 隐式计算方法
在线阅读 下载PDF
有限元-有限差分法二维波动逆时偏移初探 被引量:13
15
作者 董渊 杨慧珠 杜启振 《石油大学学报(自然科学版)》 CSCD 北大核心 2003年第6期25-29,共5页
提出了采用有限元有限差分实现二维波动方程的逆时偏移算法。该方法在空间上 ,联合采用有限元法和有限差分法 ;对于地表 (水平 )方向 ,使用有限元法进行离散 ,将原方程转化为一个一维 (深度和时间 )问题的方程组 ;在深度和时间方向上 ,... 提出了采用有限元有限差分实现二维波动方程的逆时偏移算法。该方法在空间上 ,联合采用有限元法和有限差分法 ;对于地表 (水平 )方向 ,使用有限元法进行离散 ,将原方程转化为一个一维 (深度和时间 )问题的方程组 ;在深度和时间方向上 ,采用有限差分法来求解。介绍了算法的基本原理 ,给出了计算实例并与使用F K(频率波数 )域相移法、频率空间域有限差分法的结果进行了比较。与采用有限元的偏移方法相比 ,本方法可以节省大量内存 ;与采用有限差分的偏移方法相比 ,可以在一定程度上提高计算精度。 展开更多
关键词 波动方程 逆时偏移 有限元 有限差分 计算精度 地震勘探 数据处理 油气勘探
在线阅读 下载PDF
基于动力有限元法的车辙预估计算方法及其应用 被引量:8
16
作者 李丽民 何兆益 +1 位作者 张国祥 雷婷 《振动与冲击》 EI CSCD 北大核心 2010年第11期143-146,共4页
基于运动车辆与路面的耦合作用,提出了更加符合实际的汽车动载模型和轮载累加作用时间、代表温度和加载修正系数的计算方法及合理的材料模型参数,采用动力有限元法,建立了沥青路面车辙预估计算方法;采用车辙预估计算方法对某环道试验柔... 基于运动车辆与路面的耦合作用,提出了更加符合实际的汽车动载模型和轮载累加作用时间、代表温度和加载修正系数的计算方法及合理的材料模型参数,采用动力有限元法,建立了沥青路面车辙预估计算方法;采用车辙预估计算方法对某环道试验柔性路面和某高速公路半刚性基层沥青路面的车辙进行了预估计算,计算值和实测值基本吻合,表明该车辙顶估计算方法是可行的,并具有较高的计算精度。该方法能在路面设计阶段就可以有针对性地解决路面使用中的车辙问题,具有重要的工程应用价值。 展开更多
关键词 动力有限元 计算方法 车辙预估 应用
在线阅读 下载PDF
钢-混凝土组合梁双层梁有限元分析方法 被引量:15
17
作者 邱文亮 张哲 黄才良 《大连理工大学学报》 EI CAS CSCD 北大核心 2003年第1期101-103,共3页
首先提出钢-混凝土组合梁的双层梁有限元分析方法,推导了考虑滑移的剪力连接件单元刚度矩阵.利用本方法计算了简支组合梁在荷载作用下的挠度、应力以及滑移分布.计算结果与实验结果一致,证明本方法分析钢-混凝土组合梁的正确性和实用性.
关键词 -混凝土组合梁 有限元分析方法 滑移 剪力连接件 双层梁模型 刚度矩阵 挠度 应力
在线阅读 下载PDF
电子器件振动特性有限元模型参数的等效计算方法 被引量:6
18
作者 杜平安 刘建涛 刘孝保 《电子学报》 EI CAS CSCD 北大核心 2010年第8期1867-1873,共7页
基于自由度缩聚技术和模态叠加原理,提出两种电子器件振动特性有限元模型参数的等效计算方法:范数极小值等效法和约束基频等效法.两种方法依据器件质量、约束基频、内部布局等特性,将具有非均匀质量和刚度分布的器件等效为均质长方形质... 基于自由度缩聚技术和模态叠加原理,提出两种电子器件振动特性有限元模型参数的等效计算方法:范数极小值等效法和约束基频等效法.两种方法依据器件质量、约束基频、内部布局等特性,将具有非均匀质量和刚度分布的器件等效为均质长方形质量块,使电子器件的网格划分和数值计算大为简化.通过印制板固有特性和响应结果的对比分析,验证了等效方法的有效性和实用性. 展开更多
关键词 电子器件 有限元模型 振动特性 等效计算方法 基频
在线阅读 下载PDF
基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用 被引量:5
19
作者 龚小权 贾洪印 +2 位作者 陈江涛 赵辉 周桂宇 《空气动力学学报》 CSCD 北大核心 2019年第1期121-132,共12页
为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进... 为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。 展开更多
关键词 隐式迭代 间断Galerkin有限元方法 GMRES 精确雅可比矩阵 计算效率
在线阅读 下载PDF
Navier-Stokes方程间断Galerkin有限元方法研究 被引量:24
20
作者 于剑 阎超 《力学学报》 EI CSCD 北大核心 2010年第5期962-970,共9页
通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤... 通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤.同时还给出了一种简单有效的计算方法来对物面边界进行高阶近似.为了能够对NS方程进行精度测试,采用对原始系统添加源项的方法构造精确解.二维Euler和NS系统的精度测试表明该方法达到了DG方法的理论精度.二维圆柱无黏绕流的计算结果表明关于物面边界的高阶近似方法能够保持DG方法原有的精度.卡门涡街数值模拟则进一步验证了该方法的正确性并且显示出DG方法较高的计算精度和分辨率. 展开更多
关键词 间断Galerkin有限元方法 NAVIER-STOKES方程 黏性项 计算流体力学
在线阅读 下载PDF
上一页 1 2 34 下一页 到第
使用帮助 返回顶部