期刊文献+
共找到141篇文章
< 1 2 8 >
每页显示 20 50 100
金属体积成形的无网格Galerkin模拟 被引量:3
1
作者 刘红生 邢忠文 杨玉英 《材料科学与工艺》 EI CAS CSCD 北大核心 2010年第2期242-244,250,共4页
为避免金属体积成形有限元法模拟中网格畸变造成网格重划和模拟精度降低,采用无网格法模拟金属体积成形.利用无网格法近似位移场,建立金属体积成形的无网格法连续性控制方程,采用罚函数法施加本征边界条件和体积不变条件,基于Markov变... 为避免金属体积成形有限元法模拟中网格畸变造成网格重划和模拟精度降低,采用无网格法模拟金属体积成形.利用无网格法近似位移场,建立金属体积成形的无网格法连续性控制方程,采用罚函数法施加本征边界条件和体积不变条件,基于Markov变分原理推导了金属体积成形的无网格Galerkin求解列式.用数值计算法求解该列式,实现金属体积成形的无网格模拟.数值结果表明,无网格法能有效处理金属体积成形中出现的大变形,避免了网格畸变和重划,具有较高的模拟精度. 展开更多
关键词 数值模拟 有限元 罚函数法 本征边界条件 无网格法
在线阅读 下载PDF
济南边界层风廓线雷达与L波段雷达大风探空测风对比 被引量:18
2
作者 王栋成 邱粲 +1 位作者 董旭光 曹洁 《气象》 CSCD 北大核心 2019年第8期1169-1180,共12页
基于相同站点的济南边界层风廓线雷达(WPR)和L波段雷达大风探空测风(以下简称LW)资料,采用相关、拟合、廓线分析等方法,系统地对比两者大风(≥10.8 m·s^-1)数据的各统计特征、相关性、时空变化规律的异同。结果表明:(1)两者u分量、... 基于相同站点的济南边界层风廓线雷达(WPR)和L波段雷达大风探空测风(以下简称LW)资料,采用相关、拟合、廓线分析等方法,系统地对比两者大风(≥10.8 m·s^-1)数据的各统计特征、相关性、时空变化规律的异同。结果表明:(1)两者u分量、v分量、风向、风速的相关系数分别为0.973、0.965、0.994、0.665,标准偏差分别为2.04 m·s^-1、2.88 m·s^-1、10.82°、2.53 m·s^-1,大风数据相关性总体较高,且相关系数风向优于风速的、降水样本优于非降水的,表明WPR观测大风数据可信、降水期间WPR水平大风数据可以使用。(2)两者大风风向、风速基本一致,但在低层差异较高层的大,尤其是340 m及以下高度层更显著,降水、非降水、全部样本在2980、1900、2740 m以下WPR风向值均略小于LW,而在1300、2740、1660 m以下WPR风速值均略大于LW。(3)两者资料互补性总体较好,大风风向、风速、v分量差值随高度变化符合对数律递增或递减规律,但u分量需分段拟合;在进行回归方程求算和资料互补时,应考虑差异随高度等的变化。 展开更多
关键词 边界层风廓线雷达 l波段雷达 大风 标准偏差法 资料互补
在线阅读 下载PDF
基于核重构的最小二乘配点法求解Helmholtz方程 被引量:4
3
作者 史宝军 袁明武 舒东伟 《力学学报》 EI CSCD 北大核心 2006年第1期125-129,共5页
基于核重构思想构造近似函数,将配点法和最小二乘原理相结合对微分方程进行离散,建立了Helmholtz 方程的最小二乘配点格式,并分别研究了Helmholtz方程的波传播问题和边界层问题.通过数值算例可以发现,给出的数值计算结果非常接近于精确... 基于核重构思想构造近似函数,将配点法和最小二乘原理相结合对微分方程进行离散,建立了Helmholtz 方程的最小二乘配点格式,并分别研究了Helmholtz方程的波传播问题和边界层问题.通过数值算例可以发现,给出的数值计算结果非常接近于精确解,计算精度明显高于SPH法的数值结果,且随着节点数目的增加, 其精确度越来越高,具有良好的收敛性. 展开更多
关键词 HElMHOlTZ方程 无网格法 重构核点法 最小二乘配点格式
在线阅读 下载PDF
无网格局部Petrov-Galerkin方法中本质边界条件的处理 被引量:3
4
作者 赵美玲 聂玉峰 袁占斌 《应用力学学报》 EI CAS CSCD 北大核心 2006年第3期493-495,517-518,共3页
鉴于无网格局部Petrov-Galerkin方法(MLPG)形函数的非插值性质,将一种新的本质边界处理方案——完全变换法与MLPG结合,通过变换矩阵修正形函数,使其满足Kronecker-δ条件,实现了本质边界的精确实施。进一步与MLPG中通常处理边界的罚方... 鉴于无网格局部Petrov-Galerkin方法(MLPG)形函数的非插值性质,将一种新的本质边界处理方案——完全变换法与MLPG结合,通过变换矩阵修正形函数,使其满足Kronecker-δ条件,实现了本质边界的精确实施。进一步与MLPG中通常处理边界的罚方法作了比较研究,数值结果表明新方法的可靠性与精确度。 展开更多
关键词 无网格方法 MlPG 本质边界 完全变换法
在线阅读 下载PDF
基于SPH-ALE方法的波浪水槽数值模拟 被引量:6
5
作者 李雪临 路宽 +4 位作者 韩林生 王静 王花梅 石建军 朱锐 《水利水电科技进展》 CSCD 北大核心 2016年第6期39-43,共5页
为了精确模拟波浪传播,基于光滑粒子流体动力学及任意拉格朗日欧拉(SPH-ALE)方法建立二维数值波浪水槽,在原SPH方法中引入近似黎曼求解器替代人工黏性项,采用排斥力边界条件防止流体粒子穿透固边界,海绵层内采用指数型衰减函数来消除水... 为了精确模拟波浪传播,基于光滑粒子流体动力学及任意拉格朗日欧拉(SPH-ALE)方法建立二维数值波浪水槽,在原SPH方法中引入近似黎曼求解器替代人工黏性项,采用排斥力边界条件防止流体粒子穿透固边界,海绵层内采用指数型衰减函数来消除水槽末端的波浪反射,并对规则波的传播进行数值模拟。结果表明:与采用人工黏性项的原SPH方法相比,SPH-ALE方法能够无衰减地模拟波浪传播,并可有效减小固边界附近的粒子压力振荡。 展开更多
关键词 波浪传播 光滑粒子 流体动力学 任意拉格朗日欧拉方法 数值模拟 波浪水槽 排斥力边界条件 海绵层
在线阅读 下载PDF
基于Fernald-PSO法反演气溶胶激光雷达比及其对斜程能见度的影响 被引量:8
6
作者 马愈昭 刘嘉琪 +3 位作者 王强强 熊兴隆 李猛 冯帅 《光子学报》 EI CAS CSCD 北大核心 2019年第3期195-204,共10页
为了提高斜程能见度的计算精度,提出了一种基于Fernald-PSO法求解气溶胶激光雷达比的方法.首先,以均匀层底部和顶部的大气消光系数相等为条件构建气溶胶激光雷达比和气溶胶消光系数边界值的方程.然后,再以激光雷达数据反演和AERONET网... 为了提高斜程能见度的计算精度,提出了一种基于Fernald-PSO法求解气溶胶激光雷达比的方法.首先,以均匀层底部和顶部的大气消光系数相等为条件构建气溶胶激光雷达比和气溶胶消光系数边界值的方程.然后,再以激光雷达数据反演和AERONET网站观测的气溶胶光学厚度相等为条件构造另一参量方程.最后,采用Fernald-PSO法解方程组,用得到的参量反演气溶胶消光系数计算斜程能见度.采用激光雷达和AERONET数据对该方法进行验证,分析气溶胶激光雷达比对斜程能见度反演的影响.结果表明,同一信号气溶胶激光雷达比假设值与计算值相差越大,气溶胶消光系数和斜程能见度的相对误差绝对值也越大,最大误差分别为18.93%和19.90%. 展开更多
关键词 斜程能见度 激光雷达 粒子群优化算法 Fernald法 激光雷达比 消光系数边界值 气溶胶光学厚度
在线阅读 下载PDF
耦合边界背景网格的无网格法泊松方程求解算法
7
作者 窦立远 孙哲 +1 位作者 谭思远 宗智 《计算力学学报》 北大核心 2025年第2期176-181,220,共7页
对于无网格方法而言,常见的边界条件处理方法是在边界外布置虚粒子或镜像粒子。这使得复杂几何边界下粒子排布十分困难,且准确度不高,此外导致基于局部支撑域的插值方法存在支撑域内离散点缺失的问题。而压力Poisson方程的求解需要Neum... 对于无网格方法而言,常见的边界条件处理方法是在边界外布置虚粒子或镜像粒子。这使得复杂几何边界下粒子排布十分困难,且准确度不高,此外导致基于局部支撑域的插值方法存在支撑域内离散点缺失的问题。而压力Poisson方程的求解需要Neumann或者Dirichlet边界条件,因此如何准确地施加边界条件至关重要。本文提出了一种针对无网格方法中边界条件处理的新方法,通过将原有Poisson方程转化为弱形式的积分方程,可以将边界条件十分方便地直接施加到方程中,避免了边界外增加虚粒子等问题。为了准确计算Poisson方程,在边界附近建立了局部规则背景网格,以规则背景插值点为新的未知变量参与运算,解决了离散点缺失的问题,并保留绝大部分流域的无网格属性。最后通过计算泰勒-格林涡旋问题,验证了所提方法的准确性和可行性。 展开更多
关键词 无网格方法 泊松方程 边界条件 虚拟粒子法 背景网格
在线阅读 下载PDF
局部Petrov-Galerkin方法在非线性边值问题中的应用 被引量:1
8
作者 胡德安 龙述尧 韩旭 《湖南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2007年第2期33-36,共4页
把一种真正的无网格局部Petrov-Galerkin方法用于求解非线性边值问题.为了克服一般局部Petrov-Galerkin方法计算工作量较大的问题,选择一个分段函数作为加权残值法的加权函数,简化了非线性问题中刚度矩阵的域积分.基于局部Petrov-Galer... 把一种真正的无网格局部Petrov-Galerkin方法用于求解非线性边值问题.为了克服一般局部Petrov-Galerkin方法计算工作量较大的问题,选择一个分段函数作为加权残值法的加权函数,简化了非线性问题中刚度矩阵的域积分.基于局部Petrov-Galerkin积分方程逐点建立的思想,推导了一种直接插值法用于施加本质边界条件.通过算例表明,这种局部Petrov-Galerkin方法是一种具有收敛快、精度高的方法. 展开更多
关键词 数值方法 局部Petrov—Galerkin方法 非线性边值问题 分段函数
在线阅读 下载PDF
小尺度封闭空间无网格Galerkin声场数值计算方法 被引量:2
9
作者 王海涛 曾向阳 陈玲 《西北工业大学学报》 EI CAS CSCD 北大核心 2012年第1期102-107,共6页
文章将无网格法引入小尺度封闭空间声场的计算。首先推导了适用于任意形状小尺度封闭空间的无网格Galerkin声场数值计算模型。在模型中,利用Galerkin型加权残量法推导了计算节点声压的系统方程,根据移动最小二乘近似法构造了无网格形函... 文章将无网格法引入小尺度封闭空间声场的计算。首先推导了适用于任意形状小尺度封闭空间的无网格Galerkin声场数值计算模型。在模型中,利用Galerkin型加权残量法推导了计算节点声压的系统方程,根据移动最小二乘近似法构造了无网格形函数,并给出了积分运算方案;然后对一个矩形封闭空间的声传递函数及混响时间进行了计算,通过与理论计算结果、SYSNOISE计算结果及其文献实验结果对比,证明了模型的正确性;最后对一个实际机舱进行了建模计算,并与测量结果进行了比较,证明了模型对较复杂结构的正确性及适用性。 展开更多
关键词 无网格法 小尺度封闭空间 声场计算 伽辽金法 实际机舱
在线阅读 下载PDF
非线性l-1模极小化问题的极大熵粒子群算法 被引量:4
10
作者 张建科 《计算机工程与应用》 CSCD 北大核心 2009年第13期62-64,共3页
针对非线性l-1模极小化问题,利用粒子群算法并结合极大熵函数法给出了此类问题的一种新混合算法。该算法首先利用极大熵函数将非线性l-1模极小化问题转化为一个光滑函数的无约束最优化问题,将此光滑函数作为粒子群算法的适应值函数;然... 针对非线性l-1模极小化问题,利用粒子群算法并结合极大熵函数法给出了此类问题的一种新混合算法。该算法首先利用极大熵函数将非线性l-1模极小化问题转化为一个光滑函数的无约束最优化问题,将此光滑函数作为粒子群算法的适应值函数;然后应用粒子群算法来优化此问题。数值结果表明,该算法收敛快、数值稳定性好,是求解非线性l-1模极小化问题的一种有效算法。 展开更多
关键词 粒子群算法 进化算法 l-1模极小化问题 极大熵函数
在线阅读 下载PDF
基于机器学习的IBIEM控制方程基本解构造模型
11
作者 刘中宪 黄珑 +1 位作者 孟思博 黄振恩 《计算力学学报》 北大核心 2025年第2期228-234,共7页
间接边界积分方程法IBIEM(indirect boundary integral equation method)求解波动问题时控制方程基本解构造依赖经验判断和试算,导致宽频散射求解不够稳定。本文通过粒子群优化-人工神经网络建立IBIEM控制方程基本解构造模型,以数据驱... 间接边界积分方程法IBIEM(indirect boundary integral equation method)求解波动问题时控制方程基本解构造依赖经验判断和试算,导致宽频散射求解不够稳定。本文通过粒子群优化-人工神经网络建立IBIEM控制方程基本解构造模型,以数据驱动代替经验判断,处理基本解构造过程中的不确定性。以二维峡谷对平面SH波散射IBIEM模拟为例验证所建模型的可靠性。结果表明,所建IBIEM控制方程基本解构造模型可对虚拟波源位置和数量的最优设置进行有效预测,兼顾计算效率和精度,大幅提高IBIEM求解波动问题时的稳定性和高效性;虚拟波源位置和数量最优设置方案受入射波频率和场地几何条件影响显著,且表现出非单调变化特征,依据经验设置基本解可靠性较差,以数据驱动的预测模型具有明显优势。本文所建方法可为IBIEM求解其他类型场地地震波动问题提供参考。 展开更多
关键词 间接边界积分方程法 虚拟波源 人工神经网络 粒子群优化 波动问题求解
在线阅读 下载PDF
Winkler地基上Kirchhoff板分析的RKPM无网格法
12
作者 曾祥勇 彭川海 邓安福 《岩土力学》 EI CAS CSCD 北大核心 2010年第S2期98-103,共6页
RKPM是基于核函数(样条函数)近似的一种典型无网格计算方法。完整推导了RKPM无网格法形函数的一阶、二阶导函数,结合Kirchhoff薄板理论,建立了Winkler地基上薄板弯曲挠度的RKPM无网格法求解控制方程,编制了相应的程序,算例分析表明,该... RKPM是基于核函数(样条函数)近似的一种典型无网格计算方法。完整推导了RKPM无网格法形函数的一阶、二阶导函数,结合Kirchhoff薄板理论,建立了Winkler地基上薄板弯曲挠度的RKPM无网格法求解控制方程,编制了相应的程序,算例分析表明,该方法有效可行。 展开更多
关键词 计算方法 RKPM 无网格法 Kirchhoff薄板 WINKlER地基
在线阅读 下载PDF
Winkler地基上中厚板计算的重构核粒子法
13
作者 曾祥勇 朱爱军 邓安福 《重庆大学学报(自然科学版)》 EI CAS CSCD 北大核心 2009年第1期44-49,共6页
重构核粒子法(reproducing kernel particle method,RKPM)是一种基于核函数近似的典型无网格方法。以RKPM法插值形函数为基础,基于Mindlin中厚板理论,建立Winkler地基上中厚板弯曲挠度的RKPM法求解控制方程,编制相应的计算程序,算例分... 重构核粒子法(reproducing kernel particle method,RKPM)是一种基于核函数近似的典型无网格方法。以RKPM法插值形函数为基础,基于Mindlin中厚板理论,建立Winkler地基上中厚板弯曲挠度的RKPM法求解控制方程,编制相应的计算程序,算例分析表明该方法有效可行。 展开更多
关键词 计算方法 重构核粒子法 无网格法 MINDlIN板 WINKlER地基
在线阅读 下载PDF
双互易杂交边界点方法求解Helmholtz方程
14
作者 杨庆年 郑俊杰 苗雨 《湖南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2010年第11期31-35,共5页
提出了一种新的边界类型无网格法——双互易杂交边界点方法,它将杂交边界点法和双互易法结合,来求解Helmholtz方程.该方法将Helmholtz方程的解分为通解和特解两部分,通解使用杂交边界点方法求解,特解则利用径向基函数近似.该方法只需要... 提出了一种新的边界类型无网格法——双互易杂交边界点方法,它将杂交边界点法和双互易法结合,来求解Helmholtz方程.该方法将Helmholtz方程的解分为通解和特解两部分,通解使用杂交边界点方法求解,特解则利用径向基函数近似.该方法只需要边界上离散的点,域内少数的点仅仅是为了径向基函数插值.通过数值算例对影响该方法性能的参数进行了研究.数值算例表明,该方法在求解Helmholtz方程时有较高的精度和数值稳定性. 展开更多
关键词 无网格法 双互易杂交边界点法 HElMHOlTZ方程 移动最小二乘近似 径向基函数
在线阅读 下载PDF
一种油冷电机计算搅油损失的方法和研究
15
作者 周航 李松松 +3 位作者 孙友情 阎阔 崔江涛 陈德鑫 《汽车技术》 北大核心 2025年第7期57-62,共6页
为解决油冷电机运行中的搅油损失问题,提出了一种油冷电机计算搅油损失的方法,通过结合流体力学和热力学理论,采用无网格粒子法三维动态计算流体力学(CFD)多物理场耦合仿真来计算搅油损失。通过与实际试验数据对比分析验证其准确性,为... 为解决油冷电机运行中的搅油损失问题,提出了一种油冷电机计算搅油损失的方法,通过结合流体力学和热力学理论,采用无网格粒子法三维动态计算流体力学(CFD)多物理场耦合仿真来计算搅油损失。通过与实际试验数据对比分析验证其准确性,为优化油冷电机设计和提升效率提供参考。 展开更多
关键词 油冷电机 搅油损失 无网格粒子法 三维动态CFD 多物理场耦合仿真
在线阅读 下载PDF
用边界轮廓法求解断裂力学L积分
16
作者 陈颂英 唐委校 +1 位作者 冷霞 郑雯 《郑州工业大学学报》 1998年第4期33-36,共4页
证明断裂力学L积分方程的被积函数的散度等于零,将面积分转化为线积分,使求解问题的维数降低两维。利用边界轮廓法的结果,使平面断裂问题L积分的求解转化为边界点的位移和面力线性迭加,避免了求解数值积分。
关键词 断裂力学 边界元法 l积分 边界轮廓法
全文增补中
负压梯度下自由面流动模拟的无网格粒子法Dirichlet边界条件 被引量:2
17
作者 韩沛东 刘启新 +2 位作者 周子棋 孙中国 席光 《西安交通大学学报》 EI CAS CSCD 北大核心 2022年第8期141-148,共8页
为改善无网格粒子法自由表面附近计算精度及稳定性,尤其在负压梯度条件下较难收敛的状况,提出了一种基于虚拟边界的狄利克雷边界条件赋值方法。以移动粒子半隐式法为例,采用最小二乘法对自由表面粒子进行多项式拟合,以拟合曲线作为虚拟... 为改善无网格粒子法自由表面附近计算精度及稳定性,尤其在负压梯度条件下较难收敛的状况,提出了一种基于虚拟边界的狄利克雷边界条件赋值方法。以移动粒子半隐式法为例,采用最小二乘法对自由表面粒子进行多项式拟合,以拟合曲线作为虚拟边界,用自由表面粒子到虚拟边界的偏差矢量对自由表面粒子赋值进行修正。通过方形液滴旋转算例验证了算法的有效性与准确性,获得了自由表面平滑的模拟结果。进一步研究了关键参数的选择对边界条件的表达和计算精度的影响,发现当拟合搜索半径为4倍初始粒子间距时,可减小拟合误差,防止边界畸变。提出的改进方法提高了自由表面核截断区域的计算精度与稳定性,可显著改善负压梯度下大变形自由曲面外缘粒子因梯度力偏差而导致的粒子飘散情况,为改善大变形复杂界面流动计算提供了新的思路。 展开更多
关键词 移动粒子半隐式法 自由面 狄利克雷边界条件 虚拟边界
在线阅读 下载PDF
无网格法中SPH法和MLS法的对比 被引量:1
18
作者 杨宇红 秦伶俐 周喆 《农机化研究》 北大核心 2004年第5期113-114,117,共3页
无网格法是近期发展起来的一种数值计算方法,其的形函数构造有多种方法。为此,着重介绍了较常用的两种方法:移动最小二乘法(MLS)和光滑粒子法(SPH),并对这两种不同的方法进行了阐述,对比了这两种插值函数的异同,最后又用数值算例对其计... 无网格法是近期发展起来的一种数值计算方法,其的形函数构造有多种方法。为此,着重介绍了较常用的两种方法:移动最小二乘法(MLS)和光滑粒子法(SPH),并对这两种不同的方法进行了阐述,对比了这两种插值函数的异同,最后又用数值算例对其计算精度进行了比较。 展开更多
关键词 无网格法 SPH法 MlS法 数值计算
在线阅读 下载PDF
弹性颗粒声散射问题直接模拟的IB-LBM研究 被引量:1
19
作者 蔡堉楠 黎胜 +1 位作者 鲁建华 纪良 《振动与冲击》 EI CSCD 北大核心 2020年第17期219-225,244,共8页
针对弹性颗粒声散射问题提出了一种基于改进IB-IBM的直接模拟方法。在该模拟方法中,流体采用格子Boltzmann方法(LBM)求解,固体采用光滑点插值法(S-PIM)求解,流固耦合采用浸入式边界方法(IBM)处理,声场通过计算求得的流场压力瞬时扰动直... 针对弹性颗粒声散射问题提出了一种基于改进IB-IBM的直接模拟方法。在该模拟方法中,流体采用格子Boltzmann方法(LBM)求解,固体采用光滑点插值法(S-PIM)求解,流固耦合采用浸入式边界方法(IBM)处理,声场通过计算求得的流场压力瞬时扰动直接获得。通过二维近似刚性颗粒声散射算例,该方法得到初步验证;与此同时,计算了二维弹性颗粒声散射问题,讨论了颗粒弹性,声传播介质类型对散射声场的影响,进一步验证了该方法的有效性。 展开更多
关键词 浸入式边界-格子Boltzmann方法(IB-lBM) 声散射 弹性颗粒 流固耦合 直接模拟
在线阅读 下载PDF
颗粒振动影响动量传递过程的格子Boltzmann方法模拟
20
作者 贾龙菲 付少童 +3 位作者 向星 张华海 张弢 王利民 《化工学报》 EI CSCD 北大核心 2023年第2期735-747,共13页
流化床内部颗粒振动对传递过程有着重要影响。格子Boltzmann方法耦合改进的浸入运动边界法模拟了不同振幅比A/D和频率比k=fe/f0下的单颗粒振动情况,并研究了不同排布和间距的双颗粒振动对传递过程中升阻力系数以及涡脱落频率的影响。结... 流化床内部颗粒振动对传递过程有着重要影响。格子Boltzmann方法耦合改进的浸入运动边界法模拟了不同振幅比A/D和频率比k=fe/f0下的单颗粒振动情况,并研究了不同排布和间距的双颗粒振动对传递过程中升阻力系数以及涡脱落频率的影响。结果表明,颗粒Reynolds数Re=100,单颗粒横向振动时,振幅增大导致锁定区间变大,颗粒锁定区间内的曳力系数大于锁定区间外,有利于传递。单颗粒流向振动,A/D=1.50时,随振动频率增大,流体流动模式为:2S模式→2P模式→2P+2S模式→混沌。相同振幅下k<1.25时,颗粒横向振动的曳力系数大于流向振动的曳力系数;k>1.25时则与之相反。因此,当k<1.25时,横向振动更有利于传递;k>1.25时,流向振动更有利于传递。串联双颗粒相互抑制涡的形成使曳力系数减小,不利于传递;相反,并联双颗粒促进传递作用,且在间距H=3D时效果最佳。以上数值结果为强化传递过程提供了一种思路。 展开更多
关键词 格子BOlTZMANN方法 浸入运动边界法 颗粒振动 曳力系数 泻涡模式 锁定现象
在线阅读 下载PDF
上一页 1 2 8 下一页 到第
使用帮助 返回顶部