期刊文献+
共找到33篇文章
< 1 2 >
每页显示 20 50 100
沥青混合料微观力学分析综述 被引量:39
1
作者 胡霞光 《长安大学学报(自然科学版)》 EI CAS CSCD 北大核心 2005年第2期6-10,共5页
针对当前沥青混合料研究方法的局限性,分析了国内外该领域的研究现状与存在问题。介绍了沥青混合料数字化方法及常用的力学分析方法。给出了运用微观力学方法研究沥青混合料的基本思路,即以二维和三维混合料数字图像为基础,数字重构沥... 针对当前沥青混合料研究方法的局限性,分析了国内外该领域的研究现状与存在问题。介绍了沥青混合料数字化方法及常用的力学分析方法。给出了运用微观力学方法研究沥青混合料的基本思路,即以二维和三维混合料数字图像为基础,数字重构沥青混合料微观结构,研究基于微观力学的沥青混合料集料接触模型,建立基于离散元的混合料微观空间结构的力学模型,应用离散元方法和微观力学理论分析该模型在劲度模量试验中的力学行为,进行沥青混合料劲度模量的虚拟试验。通过试验检验与校准,得出具有高可靠性、高精确度的劲度模量力学预估方法。 展开更多
关键词 道路工程 沥青混合料 微观力学 数字重构 离散元 综述
下载PDF
岩石微观胶结模型及离散元数值仿真方法初探 被引量:45
2
作者 蒋明镜 陈贺 刘芳 《岩石力学与工程学报》 EI CAS CSCD 北大核心 2013年第1期15-23,共9页
基于胶结铝棒接触力学特性的系列室内微观试验,提出用于模拟岩石的微观胶结接触模型,并将该模型引入二维离散元商业软件PFC2D,模拟岩石室内单双轴压缩、直接拉伸和巴西试验,并将数值试验结果与室内试验实测结果进行对比分析,验证该数值... 基于胶结铝棒接触力学特性的系列室内微观试验,提出用于模拟岩石的微观胶结接触模型,并将该模型引入二维离散元商业软件PFC2D,模拟岩石室内单双轴压缩、直接拉伸和巴西试验,并将数值试验结果与室内试验实测结果进行对比分析,验证该数值方法在反映岩石基本力学特性方面的有效性,研究结果表明,基于水泥微观胶结模型的离散单元法能够合理地描述岩石的主要力学特性;通过胶结破坏数目与轴向应力之间的关系曲线可以识别岩石裂纹扩展的规律及各个阶段,并能够确定裂纹扩展过程中的起裂应力、裂纹损伤应力和峰值应力。 展开更多
关键词 岩石力学 离散元法 胶结模型 裂纹扩展 裂纹损伤应力
下载PDF
离散单元法在土壤机械特性动态仿真中的应用进展 被引量:30
3
作者 张锐 李建桥 李因武 《农业工程学报》 EI CAS CSCD 北大核心 2003年第1期16-19,共4页
简述了离散单元法的基本原理及其国内外发展现状 ,介绍了目前地面力学研究领域中离散单元法在土壤机械特性动态仿真中的应用情况 ,分析了离散单元法在土壤动态行为仿真中应用的可行性 ,指出离散单元法应用于土壤这样的多相不连续复合体... 简述了离散单元法的基本原理及其国内外发展现状 ,介绍了目前地面力学研究领域中离散单元法在土壤机械特性动态仿真中的应用情况 ,分析了离散单元法在土壤动态行为仿真中应用的可行性 ,指出离散单元法应用于土壤这样的多相不连续复合体系中 ,以离散单元的总体行为来描述土壤动态行为具有独特的优越性 ,提出了离散单元法在土壤动态特征研究中的发展趋势和近期需要解决的关键问题。 展开更多
关键词 机械特性 应用 离散单元法 土壤动态仿真 力学模型
下载PDF
推土板表面形态对土壤动态行为影响的离散元模拟 被引量:28
4
作者 张锐 李建桥 +1 位作者 周长海 许述财 《农业工程学报》 EI CAS CSCD 北大核心 2007年第9期13-19,共7页
土壤行为的变化过程直接受触土部件结构的影响,研究不同表面推土板作用下土壤动态行为的变化规律和影响因素对于触土部件的优化设计具有重要意义。根据液桥对土壤微观力学结构的影响,在传统离散单元法理论的基础上,建立了土壤颗粒接触... 土壤行为的变化过程直接受触土部件结构的影响,研究不同表面推土板作用下土壤动态行为的变化规律和影响因素对于触土部件的优化设计具有重要意义。根据液桥对土壤微观力学结构的影响,在传统离散单元法理论的基础上,建立了土壤颗粒接触非线性力学模型,并对光滑和波纹形表面形态推土板前端土壤的动态行为进行了离散元模拟。结果表明:离散元模拟不仅准确再现了实验中表面形态对板面前端土壤动态行为影响的定性和定量结果,而且通过离散元细观分析从土壤内部土壤颗粒运动合理解释了土壤宏观形态的变化规律以及推土板表面形态对土壤动态行为的影响,即波纹形表面使板面前端土壤波动频率变大,而波动幅度变小。该研究为土壤-机械相互作用分析提供了新的研究思路和手段。 展开更多
关键词 离散单元法 土壤动态行为 非线性力学模型 宏细观模拟 推土板表面形态
下载PDF
传播介质特性对爆破振动衰减规律的影响 被引量:23
5
作者 吕涛 李海波 +3 位作者 周青春 夏祥 刘亚群 李俊如 《防灾减灾工程学报》 CSCD 2008年第3期335-341,共7页
以《工程岩体分级标准》为依据,量化作为爆破地震波载体的传播介质的特性,并采用离散单元法分析爆破地震波在传播介质中的衰减过程。数值模拟过程中详细阐述了爆破荷载、介质阻尼及相关参数的确定方法,分析了传播介质特性对爆破振动... 以《工程岩体分级标准》为依据,量化作为爆破地震波载体的传播介质的特性,并采用离散单元法分析爆破地震波在传播介质中的衰减过程。数值模拟过程中详细阐述了爆破荷载、介质阻尼及相关参数的确定方法,分析了传播介质特性对爆破振动速度衰减规律及萨道夫斯公式中K、α值的影响。研究表明,在爆破荷载作用下,近距范围内爆破振动速度迅速衰减,中远距范围内爆破振动速度衰减趋缓;传播介质为Ⅰ~Ⅲ级岩体时,随着岩体质量降低,K值有增大趋势,对于Ⅳ~Ⅴ级岩体,K值变化没有明确的规律;α值总体上具有随岩体质量降低而增大的趋势,K、α值具有与爆破振动速度近似相同的增减关系。 展开更多
关键词 传播介质 爆破振动 衰减规律 离散单元法(dem)
下载PDF
PFC滑坡模拟二、三维建模方法研究 被引量:27
6
作者 曹文 李维朝 +2 位作者 唐斌 邓刚 李俊峰 《工程地质学报》 CSCD 北大核心 2017年第2期455-462,共8页
滑坡运动过程模拟避免了试验尺度与监测手段的限制,可以详细观察破坏过程,是定量评估滑坡灾变风险的重要研究手段。开展滑坡运动过程模拟的首要工作是建立滑坡模型。颗粒流程序(PFC)虽然是滑坡运动过程模拟应用最广泛的程序之一,但在建... 滑坡运动过程模拟避免了试验尺度与监测手段的限制,可以详细观察破坏过程,是定量评估滑坡灾变风险的重要研究手段。开展滑坡运动过程模拟的首要工作是建立滑坡模型。颗粒流程序(PFC)虽然是滑坡运动过程模拟应用最广泛的程序之一,但在建立滑坡模型的前处理方面较弱,使得其在滑坡运动过程模拟中的推广应用受到限制。对此,本文指出了BallBall和Ball-Wall两种建模方法的适用性、优缺点及滑体滑床边界确定方法,并以2014年地震触发的红石岩滑坡为例,以Brick填充法为基础,从获取地形数据、确定滑体和滑床区域、建立滑体和滑床几何模型、生成颗粒模型4个步骤出发,提出了基于数字等高线地形图建立PFC二、三维复杂滑坡模型的前处理方法,弥补了PFC软件前处理的不足,从而为今后基于PFC的滑坡运动过程模拟提供有益的帮助。 展开更多
关键词 离散元 颗粒流 滑坡 数值模拟
下载PDF
推土板切土角对干土壤动态行为影响的离散元模拟 被引量:20
7
作者 张锐 李建桥 +1 位作者 许述财 李因武 《吉林大学学报(工学版)》 EI CAS CSCD 北大核心 2007年第4期822-827,共6页
根据液桥对土壤微观力学结构的影响,在传统离散单元法理论的基础上,建立了土壤颗粒接触非线性力学模型,并对不同切土角推土板前端干土壤的动态行为进行了离散元模拟。准确再现了实验室实验中切土角对板面前端干土壤动态行为影响的定性... 根据液桥对土壤微观力学结构的影响,在传统离散单元法理论的基础上,建立了土壤颗粒接触非线性力学模型,并对不同切土角推土板前端干土壤的动态行为进行了离散元模拟。准确再现了实验室实验中切土角对板面前端干土壤动态行为影响的定性和定量结果,并通过离散元细观分析从土壤内部土壤颗粒运动合理解释了干土壤宏观形态的变化规律和影响因素。 展开更多
关键词 农业工程 土壤动态行为 离散单元法 宏细观模拟 推土板切土角
下载PDF
部件复杂表面影响土壤扰动行为的离散元宏细观分析 被引量:19
8
作者 张锐 李建桥 +2 位作者 李因武 刘子斌 陈斌 《吉林大学学报(工学版)》 EI CAS CSCD 北大核心 2009年第5期1218-1223,共6页
通过分析毛管水形成液桥对土壤微观力学结构的影响,引入并行约束来表征土壤颗粒间液桥的黏性作用,建立了土壤颗粒接触非线性力学模型,并对波纹形复杂表面推土板前端土壤的动态扰动行为进行了离散元宏细观模拟。结果表明,离散元模拟不仅... 通过分析毛管水形成液桥对土壤微观力学结构的影响,引入并行约束来表征土壤颗粒间液桥的黏性作用,建立了土壤颗粒接触非线性力学模型,并对波纹形复杂表面推土板前端土壤的动态扰动行为进行了离散元宏细观模拟。结果表明,离散元模拟不仅准确再现了实验分析中复杂表面结构对土壤动态扰动行为影响的研究结果,而且从土壤颗粒运动以及土块破碎等细观角度准确分析了触土部件复杂表面对土壤动态扰动行为的影响因素和机理。 展开更多
关键词 农业工程 土壤动态扰动行为 离散单元法 非线性力学模型 宏细观模拟 波纹形表面
下载PDF
离散单元法在框架结构地震反应分析中的应用 被引量:14
9
作者 王强 吕西林 《地震工程与工程振动》 CSCD 北大核心 2004年第5期73-78,共6页
将离散单元法应用于框架结构的地震反应分析,并建立了框架结构的离散单元计算模型。在数值求解过程中,本文采用同步动态松弛法确立与求解结构单元平衡关系,避免了动态松弛法由于路径相关性所造成的误差。本文编制了框架结构地震反应时... 将离散单元法应用于框架结构的地震反应分析,并建立了框架结构的离散单元计算模型。在数值求解过程中,本文采用同步动态松弛法确立与求解结构单元平衡关系,避免了动态松弛法由于路径相关性所造成的误差。本文编制了框架结构地震反应时程分析程序,并以一实际结构为例,进行了时程分析。算例表明,本文所提出的方法简单、有效,为工程结构地震反应分析提供了新的思路与方法。 展开更多
关键词 离散单元法 框架结构 地震反应 同步动态松弛法
下载PDF
基于离散单元法的土壤力学接触模型的建立 被引量:11
10
作者 翟力欣 姬长英 《江西农业学报》 CAS 2008年第9期108-111,共4页
简述了离散单元法的发展现状,介绍了目前离散单元法的应用情况,分析了离散单元法在土壤行为中应用的可行性,指出离散单元法应用于土壤这样的多相不连续复合体系中,以离散单元的总体行为来描述土壤动态行为具有独特的优越性,提出了在离... 简述了离散单元法的发展现状,介绍了目前离散单元法的应用情况,分析了离散单元法在土壤行为中应用的可行性,指出离散单元法应用于土壤这样的多相不连续复合体系中,以离散单元的总体行为来描述土壤动态行为具有独特的优越性,提出了在离散单元法前提下建立土壤力学接触模型的方法。 展开更多
关键词 离散单元法 力学模型 土壤细观模型
下载PDF
Distinct element method analysis and field experiment of soil resistance applied on the subsoiler 被引量:10
11
作者 Li Bo Liu Fanyi +2 位作者 Mu Junying Chen Jun Han Wenting 《International Journal of Agricultural and Biological Engineering》 SCIE EI CAS 2014年第1期54-59,共6页
Since the design of the subsoiler is a complex work,the interaction between the subsoiler and soil was investigated by using Distinct Element Method(DEM)in this study.Based on the traditional discrete element theory,t... Since the design of the subsoiler is a complex work,the interaction between the subsoiler and soil was investigated by using Distinct Element Method(DEM)in this study.Based on the traditional discrete element theory,the 3D model of soil particles and the subsoiler were established after considering the liquid bridge force between soil particles.The operating resistance curves of the subsoiler were achieved after the DEM simulation at a speed of 1 m/s,and three depths of 180 mm,220 mm and 260 mm,respectively.The simulation curves agreed well with the field experimental results based on relative errors of 2.96%,14.95%and 7.15%,respectively,at three depths.All these data proved that it was feasible and favorable to analyze the performance of the subsoiler by using the DEM and it is of important significance for studying and further optimizing the structure of the subsoiler. 展开更多
关键词 distinct element method(dem) parallel bond SUBSOILER MODELING SIMULATION
原文传递
滚石离散元数值模拟的参数反演(英文) 被引量:8
12
作者 蒋景彩 能野一美 山上拓男 《岩石力学与工程学报》 EI CAS CSCD 北大核心 2008年第12期2418-2430,共13页
离散元法(DEM)是数值模拟崩塌落石运动的一个有力方法,但确定离散元输入参数是个难题。提出一个基于已知落石轨迹的离散元参数的反演计算方法。落石轨迹可以通过现场落石试验或者野外实地调查实际落石过程获得。该方法的特征为反演得到... 离散元法(DEM)是数值模拟崩塌落石运动的一个有力方法,但确定离散元输入参数是个难题。提出一个基于已知落石轨迹的离散元参数的反演计算方法。落石轨迹可以通过现场落石试验或者野外实地调查实际落石过程获得。该方法的特征为反演得到的离散元参数可以间接反映边坡与落石冲击时的变形以及植被缓冲效应。建议的几个保证解的有效性得到验证,稳定性措施效果良好。室内落石试验结果和现场落石工程实录验证证明了该算法的有效性和实用性。 展开更多
关键词 边坡工程 落石 离散元法 反演计算 dem参数 落石轨迹
下载PDF
岩石类脆性材料的离散元建模方法研究 被引量:5
13
作者 叶勇 徐西鹏 《郑州轻工业学院学报(自然科学版)》 CAS 2008年第5期62-66,共5页
通过单轴压缩、拉伸及双轴压缩模拟试验,分析了岩石脆性材料离散单元法建模过程中各个参数的关系.模拟试验表明:通过建立杨氏模量E及泊松比v与颗粒接触刚度的关系可以反向确定接触刚度值;法向接触刚度与法向连接强度的比值对岩石离散单... 通过单轴压缩、拉伸及双轴压缩模拟试验,分析了岩石脆性材料离散单元法建模过程中各个参数的关系.模拟试验表明:通过建立杨氏模量E及泊松比v与颗粒接触刚度的关系可以反向确定接触刚度值;法向接触刚度与法向连接强度的比值对岩石离散单元模型分析的影响很小;法向与切向连接强度的比值影响裂纹的破坏方式;围压的大小控制双轴试验时裂纹的类型.基于此,建立了此类了材料基于现象的宏观参数与其微观参数之间的联系规律. 展开更多
关键词 离散单元法(dem) 颗粒流 双轴压缩实验 建模方法
下载PDF
岩石力学离散介质数值方法及应用程序研究综述 被引量:4
14
作者 马钺 方祖烈 马强 《东北煤炭技术》 1999年第3期17-22,共6页
系统综述了岩石力学中离散介质的数值方法以及应用程序研究发展现状。
关键词 数值方法 离散单元法 离散介质 岩石力学
下载PDF
岩石切削加工动态行为的离散元模拟 被引量:5
15
作者 叶勇 徐西鹏 《工具技术》 北大核心 2008年第6期26-29,共4页
岩石切削加工时必须考虑其离散性。利用颗粒流程序对两种不同成份的岩石样本建立了离散元模型,通过改变岩石内部颗粒大小和形状从而获得了各个样本的最佳理论强度值,用离散元方法对两种不同岩石样本分别进行了压缩试验和四点抗弯试验数... 岩石切削加工时必须考虑其离散性。利用颗粒流程序对两种不同成份的岩石样本建立了离散元模型,通过改变岩石内部颗粒大小和形状从而获得了各个样本的最佳理论强度值,用离散元方法对两种不同岩石样本分别进行了压缩试验和四点抗弯试验数值模拟;监控整个切削加工过程中刀具受到的法向力和切向力,通过改变刀具切削角、切削速度及切深等参量研究了加工过程中岩石的动态行为。 展开更多
关键词 离散单元法(dem) 颗粒流 数值模拟 切削加工
下载PDF
Effects of intermediate stress on deep rock strainbursts under true triaxial stresses 被引量:1
16
作者 Lihua Hu Liyuan Yu +2 位作者 Minghe Ju Xiaozhao Li Chun’an Tang 《Journal of Rock Mechanics and Geotechnical Engineering》 SCIE CSCD 2023年第3期659-682,共24页
The effect of intermediate stress(in situ tunnel axial)on a strainburst is studied with a threedimensional(3D)bonded block distinct element method(DEM).A series of simulations of strainbursts under true triaxial in si... The effect of intermediate stress(in situ tunnel axial)on a strainburst is studied with a threedimensional(3D)bonded block distinct element method(DEM).A series of simulations of strainbursts under true triaxial in situ stress conditions(i.e.high tangential stress,moderate intermediate stress and low radial stress)of near-boundary rock masses are performed.Compared with the experimental results,the DEM model is able to capture the stress-strain response,failure pattern and energy balance of strainbursts.The fracturing processes of strainbursts are also numerically reproduced.Numerical results show that,as the intermediate stress increases:(1)The peak strain of strainbursts increases,the yield stress increases,the rock strength increases linearly,and the ratio of yield stress to rock strength decreases,indicating that the precursory information on strainbursts is enhanced;(2)Tensile and shear cracks increase significantly,and slabbing and bending of rock plates are more pronounced;and(3)The stored elastic strain energy and dissipated energy increase linearly,whereas the kinetic energy of the ejected rock fragments increases approximately exponentially,implying an increase in strainburst intensity.By comparing the experimental and numerical results,the effect of intermediate stress on the rock strength of strainbursts is discussed in order to address three key issues.Then,the Mogi criterion is applied to construct new strength criteria for strainbursts by converting the one-face free true triaxial stress state of a strainburst to its equivalent true triaxial stress state.In summary,the effect of intermediate stress on strainbursts is a double-edged sword that can enhance the rock strength and the precursory information of a strainburst,but also increase its intensity. 展开更多
关键词 Strainbursts Intermediate stress distinct element method(dem) Structural failure True triaxial Strength criterion
下载PDF
Strain localization analyses of idealized sands in biaxial tests by distinct element method 被引量:5
17
作者 Mingjing JIANG Hehua ZHU Xiumei LI 《Frontiers of Structural and Civil Engineering》 SCIE EI 2010年第2期208-222,共15页
This paper presents a numerical investigation on the strain localization of an idealized sand in biaxial compression tests using the distinct element method(DEM).In addition to the dilatancy and material frictional an... This paper presents a numerical investigation on the strain localization of an idealized sand in biaxial compression tests using the distinct element method(DEM).In addition to the dilatancy and material frictional angle,the principal stress field,and distributions of void ratio,particle velocity,and the averaged pure rotation rate(APR)in the DEM specimen are examined to illustrate the link between microscopic and macroscopic variables in the case of strain localization.The study shows that strain localization of the granular material in the tests proceeds with localizations of void ratio,strain and APR,and distortions of stress field and force chains.In addition,both thickness and inclination of the shear band change with the increasing of axial swain,with the former valued around 10-14 times of mean grain diameter and the later overall described by the Mohr-Coulomb theory. 展开更多
关键词 idealized sand strain localization numerical analyses distinct element method(dem)
原文传递
坝底水浮力对重力坝稳定性的影响分析(英文) 被引量:4
18
作者 UTILI Stefano 尹振宇 蒋明镜 《岩石力学与工程学报》 EI CAS CSCD 北大核心 2008年第8期1554-1568,共15页
着重研究一个典型的混凝土重力坝的坝底水浮力对大坝稳定性的影响,此大坝位于意大利的Cumbidanovu岛。大坝的基础由含有高度开裂的岩石所构成。首先,通过把大坝视为自由体的平衡分析法来评价大坝破坏前的最大水压力和有效排水系统对大... 着重研究一个典型的混凝土重力坝的坝底水浮力对大坝稳定性的影响,此大坝位于意大利的Cumbidanovu岛。大坝的基础由含有高度开裂的岩石所构成。首先,通过把大坝视为自由体的平衡分析法来评价大坝破坏前的最大水压力和有效排水系统对大坝稳定性的影响。然后,使用离散元方法来进一步评价开裂基岩中的水流状态,得到该水流产生的浮托力的分布,最终得到此水浮力对大坝稳定性的影响。对设计而言,上述分析考虑了岩基渗透,运用离散元方法进行模拟。研究结果表面,相比保守的平衡分析法,此模型可以得到更大的水浮力荷载。 展开更多
关键词 水利工程 离散元法 水浮力 大坝稳定性 重力坝
下载PDF
孔隙填充型深海能源土的离散元成样新方法及宏观力学特性 被引量:4
19
作者 贺洁 蒋明镜 《同济大学学报(自然科学版)》 EI CAS CSCD 北大核心 2016年第5期709-717,共9页
针对孔隙填充型水合物的赋存形态,提出了一种新的制备孔隙填充型能源土试样的数值成样方法.孔隙填充型水合物的砂性能源土试样是由砂粒和水合物颗粒混合而成的特殊的散粒体材料,具有明显的非连续特征.在土骨架的孔隙中,将水合物块体视... 针对孔隙填充型水合物的赋存形态,提出了一种新的制备孔隙填充型能源土试样的数值成样方法.孔隙填充型水合物的砂性能源土试样是由砂粒和水合物颗粒混合而成的特殊的散粒体材料,具有明显的非连续特征.在土骨架的孔隙中,将水合物块体视为由颗粒通过强胶结作用凝聚而成的团簇整体,随机填充生成不同水合物饱和度的沉积物试样,开展能源土宏观力学特性的离散元固结排水三轴压缩试验模拟,并从应力应变、体变、接触组构等方面进行分析.结果表明,水合物饱和度的增大对低饱和度的孔隙填充型能源土试样的初始弹性模量和强度影响较小.有效围压相同时,孔隙填充型能源土试样的体积剪缩量随着水合物饱和度的增大而减小;而水合物饱和度相同时,能源土试样强度及体积剪缩量随着有效围压的增大而增大.能源土试样的剪胀角随有效围压的增大而近似线性减少,随水合物饱和度的增大而趋于线性增大.颗粒间接触方向随着轴向应变的增大而朝竖直方向偏转. 展开更多
关键词 离散单元法 孔隙填充型能源土 数值生成方法 水合物饱和度 力学特性 接触组构
下载PDF
岩石切削过程中破坏机制的离散元分析 被引量:2
20
作者 叶勇 沈剑云 《工具技术》 北大核心 2008年第12期29-32,共4页
利用离散单元方法(DEM)对两种不同属性的花岗石样本模型进行了切割模拟,通过调整不同的切削参数计算并分析了岩石切削过程中的切削力、切削深度、切削比能及切削力倾角的变化规律,从细观上探讨了临界切深对岩石脆性破坏模式和塑性破坏... 利用离散单元方法(DEM)对两种不同属性的花岗石样本模型进行了切割模拟,通过调整不同的切削参数计算并分析了岩石切削过程中的切削力、切削深度、切削比能及切削力倾角的变化规律,从细观上探讨了临界切深对岩石脆性破坏模式和塑性破坏模式及其转变的影响,指出了当切削深度逐渐增大时,破坏模式依次从塑性破坏模式到过渡破坏模式最后转变为脆性破坏模式,同时切削力倾角并不随切削深度变化。将离散单元法分析的结果与实际实验进行了比较,表明该方法是合理的。 展开更多
关键词 离散单元法(dem) 临界切深 破坏模式 切削力 切削比能
下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部