期刊文献+
共找到19篇文章
< 1 >
每页显示 20 50 100
层次分析法确定影响地下洞室围岩稳定性各因素的权值 被引量:32
1
作者 李晓静 朱维申 +1 位作者 陈卫忠 武科 《岩石力学与工程学报》 EI CAS CSCD 北大核心 2004年第z2期4731-4734,共4页
在研究大型地下洞室围岩稳定性时,综合考虑了围岩强度和刚度、洞室高度和埋深、水平向地应力及其相互间距等影响因素,采用FLAC3D软件进行了16个工况的数值模拟计算,建立了影响地下洞室稳定性的层次分析法模型,对影响稳定性的各个指标因... 在研究大型地下洞室围岩稳定性时,综合考虑了围岩强度和刚度、洞室高度和埋深、水平向地应力及其相互间距等影响因素,采用FLAC3D软件进行了16个工况的数值模拟计算,建立了影响地下洞室稳定性的层次分析法模型,对影响稳定性的各个指标因素进行判断矩阵的构造、特征向量的求解以及一致性的检验。这种分析方法可以从总体上对各因素进行排序,为将来预测地下洞室围岩稳定性提供必要的依据。 展开更多
关键词 地下工程 围岩稳定性 层次分析 数值模拟
下载PDF
库岸斜坡失稳及其涌浪灾害风险分析 被引量:4
2
作者 罗冠枝 高文伟 +2 位作者 王国卫 蒋明光 陈全明 《中国地质灾害与防治学报》 CSCD 2020年第1期8-17,共10页
库岸斜坡失稳及其引发涌浪问题是库岸斜坡灾害的重要研究部分,开展库岸斜坡失稳及其涌浪灾害风险分析研究对指导库岸斜坡防灾减灾具有重要意义。以湖南省麻阳县大水冲水库滑坡为例,探讨库岸斜坡失稳引发涌浪灾害的分析方法,采用Geostudi... 库岸斜坡失稳及其引发涌浪问题是库岸斜坡灾害的重要研究部分,开展库岸斜坡失稳及其涌浪灾害风险分析研究对指导库岸斜坡防灾减灾具有重要意义。以湖南省麻阳县大水冲水库滑坡为例,探讨库岸斜坡失稳引发涌浪灾害的分析方法,采用Geostudio软件进行10年和50年一遇5日累积降雨的滑坡稳定性模拟分析,绘制滑坡影响范围内建筑物风险变化图和室内人员风险变化图。总结了在最危险工况下,大水冲水库滑坡总经济风险为249.3万元,室内人员的风险值达到27.8,并有产生翻坝的可能,但对沟口居民和水库坝体破坏的风险低。 展开更多
关键词 库岸斜坡失稳 涌浪 风险分析 稳定性模拟分析 防灾减灾
下载PDF
湖南省麻阳县石眼潭滑坡演变机制及破坏模式分析 被引量:1
3
作者 罗冠枝 方琼 +2 位作者 徐忠 李明波 隋皓月 《东华理工大学学报(自然科学版)》 CAS 2023年第3期233-238,309,共7页
为了研究湖南省红层滑坡成灾机理和破坏模式,以湖南省麻阳县石眼潭滑坡为例,运用GeoStudio软件数值模拟方法,进行10年和20年峰值降雨的滑坡稳定性模拟分析。通过对该滑坡的演变机制和破坏模式研究发现,石眼潭滑坡属于蠕滑-拉裂破坏模式... 为了研究湖南省红层滑坡成灾机理和破坏模式,以湖南省麻阳县石眼潭滑坡为例,运用GeoStudio软件数值模拟方法,进行10年和20年峰值降雨的滑坡稳定性模拟分析。通过对该滑坡的演变机制和破坏模式研究发现,石眼潭滑坡属于蠕滑-拉裂破坏模式,提出了稳定性系数与滑面软化面积比的一元二次关系方程:在自重作用工况下,软化面积比为25.56%,斜坡稳定性较好;在20年一遇暴雨工况下,稳定性系数为0.883,软化面积比为92.05%,斜坡处于不稳定状态。当软化面积比大于61.28%时,斜坡发生破坏的概率较大,与野外调查结果基本符合,表明此关系方程适用于该类红层滑坡。研究结果可为该类红层滑坡稳定性分析和防灾减灾提供参考。 展开更多
关键词 红层滑坡 稳定性模拟分析 演变机制 破坏模式 软化面积比
下载PDF
Stability analysis of cohesive soil embankment slope based on discrete element method 被引量:4
4
作者 XU Guang-ji ZHONG Kun-zhi +2 位作者 FAN Jian-wei ZHU Ya-jing ZHANG Yu-qing 《Journal of Central South University》 SCIE EI CAS CSCD 2020年第7期1981-1991,共11页
In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of th... In order to study the safety factor and instability process of cohesive soil slope, the discrete element method(DEM) was applied. DEM software PFC2 D was used to simulate the triaxial test to study the influence of the particle micro parameters on the macroscopic characteristics of cohesive soil and calibrate the micro parameters of DEM model on this basis. Embankment slope stability analysis was carried out by strength reduction and gravity increase method, it is shown that the safety factor obtained by strength reduction method is more conservative, and the arc-shaped feature of the sliding surface under the gravity increase method is more obvious. Throughout the progressive failure process, the failure trends, maximum displacements, and velocity changes obtained by the two methods were consistent. When slope was destroyed, the upper part was cracked, the middle part was sheared, and the lower part was destroyed by extrusion. The conclusions of this paper can be applied to the safety factor calculation of cohesive soil slopes and the analysis of the instability process. 展开更多
关键词 embankment slope cohesive soil stability analysis numerical simulation PFC2D software safety factor
下载PDF
Mechanical response and stability analysis of rock mass in high geostress underground powerhouse caverns subjected to excavation 被引量:18
5
作者 LI Biao DING Quan-fu +4 位作者 XU Nu-wen LEI Yi-fan XU Yuan ZHU Zhong-ping LIU Jing-fei 《Journal of Central South University》 SCIE EI CAS CSCD 2020年第10期2971-2984,共14页
To investigate the stability of rock mass in high geostress underground powerhouse caverns subjected to excavation,a microseismic(MS)monitoring system was established and the discrete element method(DEM)-based numeric... To investigate the stability of rock mass in high geostress underground powerhouse caverns subjected to excavation,a microseismic(MS)monitoring system was established and the discrete element method(DEM)-based numerical simulation was carried out.The tempo-spatial damage characteristics of rock mass were analyzed.The evolution laws of MS source parameters during the formation of a rock collapse controlled by high geostress and geological structure were investigated.Additionally,a three-dimensional DEM model of the underground powerhouse caverns was built to reveal the deformation characteristics of rock mass.The results indicated that the MS events induced by excavation of high geostress underground powerhouse caverns occurred frequently.The large-stake crown of the main powerhouse was the main damage area.Prior to the rock collapse,the MS event count and accumulated energy release increased rapidly,while the apparent stress sharply increased and then decreased.The amount and proportion of shear and mixed MS events remarkably increased.The maximum displacement was generally located near the spandrel areas.The MS monitoring data and numerical simulation were in good agreement,which can provide significant references for damage evaluation and disaster forecasting in high geostress underground powerhouse caverns. 展开更多
关键词 high geostress underground powerhouse caverns microseismic monitoring discrete element modelling stability analysis
下载PDF
A New Car Following Model:Comprehensive Optimal Velocity Model 被引量:3
6
作者 田钧方 贾斌 李新刚 《Communications in Theoretical Physics》 SCIE CAS CSCD 2011年第6期1119-1126,共8页
In this paper, we present a new ear-following model, i.e. comprehensive optimal velocity model (COVM), whose optimal velocity function not only depends on the following distance of the preceding vehicle, but also de... In this paper, we present a new ear-following model, i.e. comprehensive optimal velocity model (COVM), whose optimal velocity function not only depends on the following distance of the preceding vehicle, but also depends on the velocity difference with preceding vehicle. Simulation results show that COVM is an improvement over the previous ones theoretically. Then, the stability condition of the model is obtained by the linear stability analysis, which has shown that the model could obtain a bigger stable region than previous models in the phase diagram. Through the nonlinear analysis, the Burgers, Korteweg-de Vries (KdV) and modified KdV (mKdV) equations are derived for the triangular shock wave, the soliton wave, and the kink-antikink soliton wave. At the same time, numerical simulations are also carried out to show that the model could simulate these density waves. 展开更多
关键词 car following traffic flow
下载PDF
Sensitivity analysis of geomechanical parameters affecting a wellbore stability 被引量:3
7
作者 Abolfazl ABDOLLAHIPOUR Hamid SOLTANIAN +2 位作者 Yaser POURMAZAHERI Ezzatollah KAZEMZADEH Mohammad FATEHI-MARJI 《Journal of Central South University》 SCIE EI CAS CSCD 2019年第3期768-778,共11页
Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive stre... Wellbore stability analysis is a growing concern in oil industries. There are many parameters affecting the stability of a wellbore including geomechanical properties (e.g., elastic modulus, uni-axial compressive strength (UCS) and cohesion) and acting forces (e.g., field stresses and mud pressure). Accurate determination of these parameters is time-consuming, expensive and sometimes even impossible. This work offers a systematic sensitivity analysis to quantify the amount of each parameter’s effect on the stability of a wellbore. Maximum wellbore wall displacement is used as a stability factor to study the stability of a wellbore. A 3D finite difference method with Mohr model is used for the numerical modeling. The numerical model is verified against an analytical solution. A dimensionless sensitivity factor is developed in order to compare the results of various parameters in the sensitivity analysis. The results show a different order of importance of parameters based on rock strength. The most sensitive properties for a weak rock are the maximum horizontal stress, internal friction angle and formation pressure, respectively, while for a strong rock, the most sensitive parameters are the maximum horizontal stress, mud pressure and pore pressure, respectively. The amount of error in wellbore stability analysis inflicted by the error in estimation of each parameter was also derived. 展开更多
关键词 wellbore stability sensitivity analysis numerical modeling
下载PDF
Numerical simulation study on the influence of the ground stress field on the stability of roadways 被引量:8
8
作者 ZHANG Hualei 《Mining Science and Technology》 EI CAS 2010年第5期707-711,共5页
We adopt the concept of generalized plane strain to model a roadway in a stress field.This can avoid limitations caused by simplifying the stress analysis as plane strain.FLAC3D was used to investigate the maximum ten... We adopt the concept of generalized plane strain to model a roadway in a stress field.This can avoid limitations caused by simplifying the stress analysis as plane strain.FLAC3D was used to investigate the maximum tensile stress and displacement of a roadway in a known stress field for angles,α,between the roadway axial direction and the maximum principal stress of 0°,30°,45°,60° and 90°.This theory was applied to the analysis of an engineering case.The results indicate that stress and displacement of the surrounding rock increase as the angle,α,increases.This provides some significant guidance for a reasonable layout of roadways in a known stress field. 展开更多
关键词 ROADWAY ground stress field generalized plane strain numerical simulation
下载PDF
Failure mechanism and stability analysis of the Zhenggang landslide in Yunnan Province of China using 3D particle flow code simulation 被引量:11
9
作者 SHI Chong LI De-jie +1 位作者 CHEN Kal-hua ZHOU Jia-wen 《Journal of Mountain Science》 SCIE CSCD 2016年第5期891-905,共15页
Based on the principle of 3D particle flow code,a numerical landslide run-out model is presented to simulate the failure process of the Zhenggang landslide(in southwestern China) under the effect of water after a rain... Based on the principle of 3D particle flow code,a numerical landslide run-out model is presented to simulate the failure process of the Zhenggang landslide(in southwestern China) under the effect of water after a rainfall.The relationship between the micro-mechanical parameters and the macro-shear strength of the grain material is determined through numerical calibrations.Then the rainfall effect is considered in numerical simulations and rain-induced sliding processes are performed,which help us to discuss the mechanism of deformation and failure of this landslide together with field observations.It shows the Zhenggang landslide would most likely be activated in Zone I and would gain momentum in Zone II.In order to prevent the potential disaster,a tailing dam is advised to be designed about 175 m downstream from the current landslide boundary of Zone II.Verified by field observations,the presented landslide model can reflect the failure mechanism after rainfall.It can also provide a method to predict the potential disaster and draft disaster prevention measures. 展开更多
关键词 Back analysis Deposit avalanche Dynamic process LANDSLIDE Particle flow code
下载PDF
Simulation and analysis of airflow stability during fire in mine belt roadway 被引量:3
10
作者 ZHANG Sheng-zhu CHENG Wei-min +2 位作者 LI Qiu-jin ZHANG Rui LUO Chuan-Iong 《Journal of Coal Science & Engineering(China)》 2010年第4期375-380,共6页
According to fluid dynamics analysis during the fire, the criteria k-ε two-equation model for solving three-dimensional turbulence was determined, the pollutants generated in the fire disaster were set by adopting Mi... According to fluid dynamics analysis during the fire, the criteria k-ε two-equation model for solving three-dimensional turbulence was determined, the pollutants generated in the fire disaster were set by adopting Mixture multiphase flow, and the SIMPLE algorithm was used for solving on the basis of comprehensive consideration on the heat radia- tion and components transmission during fire. By simulating the airflow flowing state inside the tunnel during fire disaster of downward ventilation, drift ventilation, and ascensional ventilation, respectively, with regard to the actual situation of No.l, No.3, and No.5 belt roadway in Kongzhuang Coal Mine, the velocity vector distributions of pollutants under different inlet air volumes were obtained, and the damage degree and influential factors of disaster were also clear, which is helpful to control and avoid disaster during belt roadway fire. 展开更多
关键词 belt roadway disaster airflow STABILITY numerical simulation
下载PDF
Influence of depth-thickness ratio of mining on the stability of a bedding slope with its sliding surface in concave deformation 被引量:1
11
作者 Zhang Yuning Tang Jianxin +1 位作者 Li Guodong Teng Junyang 《International Journal of Mining Science and Technology》 SCIE EI CSCD 2016年第6期1117-1123,共7页
In order to study the influence of depth-thickness ratio on bedding slope stability, whose sliding surface is flexural concave in shape under mining conditions, this paper aims to study the characteristics ofdeformari... In order to study the influence of depth-thickness ratio on bedding slope stability, whose sliding surface is flexural concave in shape under mining conditions, this paper aims to study the characteristics ofdeformarion and damage of bedding sliding with depth-thickness ratios of 200:1,150:1,120:1,100:1 and 50:1 by adopting numerical simulation analysis software combined with laboratory-made "under the influence of mining variable sliding surface slope similar simulation test bed", and to propose identification methods for slope stability under the infuence of mining. The results show that mining activities under the slope reduce slope stability. With a decrease in the mining depth ratio, the influence of mining on the slope increases gradually, and the damage to the slope gradually expands, the stability of the slope grad- ually reduces, fracture occurs on the slope toe and the central fissure gradually develops to the surface, and reaches slide threshold when the depth-thickness ratio is 50:1. 展开更多
关键词 Depth-thickness ratio Bedding slope Sliding surface Concave deflection surface Mining influence Similar simulation
下载PDF
自然崩落法超前拉底数值模拟研究 被引量:13
12
作者 梁江波 刘华武 冯兴隆 《矿业研究与开发》 CAS 北大核心 2017年第2期93-96,共4页
自然崩落法拉底与聚矿槽施工的先后顺序直接影响着底部结构的稳定性,对其进行优化研究是非常必要的。运用模拟软件FLAC^(3D)建立数值模型,假设拉底超前聚矿槽,模拟分析超前拉底条件下,底部结构的应力、塑性区变化情况,研究结果表明,在... 自然崩落法拉底与聚矿槽施工的先后顺序直接影响着底部结构的稳定性,对其进行优化研究是非常必要的。运用模拟软件FLAC^(3D)建立数值模型,假设拉底超前聚矿槽,模拟分析超前拉底条件下,底部结构的应力、塑性区变化情况,研究结果表明,在拉底推进线附近区域的底部结构处于应力集中状态,在推进线后方45m左右范围内,应力发生释放,有利于底部结构的稳定及聚矿槽的开挖。 展开更多
关键词 自然崩落法 底部结构 稳定性模拟分析 超前拉底
原文传递
Excavation-induced microseismicity: microseismic monitoring and numerical simulation 被引量:9
13
作者 Nu-wen XU Chun-an TANG +4 位作者 Hong LI Feng DAI Ke MA Jing-dong SHAO Ji-chang WU 《Journal of Zhejiang University-Science A(Applied Physics & Engineering)》 SCIE EI CAS CSCD 2012年第6期445-460,共16页
The volume of influence of excavation at the right bank slope of Dagangshan Hydropower Station, southwest China, is essentially determined from microseismic monitoring, numerical modeling and conventional measurements... The volume of influence of excavation at the right bank slope of Dagangshan Hydropower Station, southwest China, is essentially determined from microseismic monitoring, numerical modeling and conventional measurements as well as in situ observations. Microseismic monitoring is a new application technique for investigating microcrackings in rock slopes. A micro- seismic monitoring network has been systematically used to monitor rock masses unloading relaxation due to continuous exca- vation of rock slope and stress redistribution caused by dam impoundment later on, and to identify and delineate the potential slippage regions since May, 2010. An important database of seismic source locations is available. The analysis of microseismic events showed a particular tempo-spatial distribution. Seismic events predominantly occurred around the upstream slope of 1180 m elevation, especially focusing on the hanging wall of fault XL316-1. Such phenomenon was interpreted by numerical modeling using RFPA-SRM code (realistic failure process analysis-strength reduction method). By comparing microseismic activity and results of numerical simulation with in site observation and conventional measurements results, a strong correlation can he obtained between seismic source locations and excavation-induced stress distribution in the working areas. The volume of influence of the rock slope is thus determined. Engineering practices show microseismic monitoring can accurately diagnose magnitude, intensity and associated tempo-spatial characteristics of tectonic activities such as faults and unloading zones. The integrated technique combining seismic monitoring with numerical modeling, as well as in site observation and conventional surveying, leads to a better understanding of the internal effect and relationship between microseismic activity and stress field in the right bank slope from different perspectives. 展开更多
关键词 Microseismic monitoring Rock slope Numerical simulation Stability analysis Dagangshan Hydropower Station
原文传递
滇东北某高速公路多因素耦合隧道边仰坡综合设计 被引量:4
14
作者 王晓帆 《公路》 北大核心 2022年第2期395-399,共5页
滇东北地区地质构造发育,地质活动频繁,隧道洞口地形、地物、地质等建设条件复杂多变。以滇东北某高速公路特长隧道工程为依托,通过隧道口边仰坡"开挖前-开挖后-支护后"3个阶段过程模拟稳定性分析,结合施工组织等因素,进行洞... 滇东北地区地质构造发育,地质活动频繁,隧道洞口地形、地物、地质等建设条件复杂多变。以滇东北某高速公路特长隧道工程为依托,通过隧道口边仰坡"开挖前-开挖后-支护后"3个阶段过程模拟稳定性分析,结合施工组织等因素,进行洞口边仰坡"永临结合"综合设计。根据工程实践情况,洞口边仰坡设计效果良好,工程措施安全有效,可为类似建设条件下的隧道边仰坡设计提供借鉴。 展开更多
关键词 隧道边仰坡 耦合特征 过程模拟稳定性分析 永临结合 综合支护设计
原文传递
The stability of rolling motion of hypersonic vehicles with slender configuration under pitching maneuvering 被引量:4
15
作者 YE YouDa TIAN Hao ZHANG XiaFeng 《Science China(Physics,Mechanics & Astronomy)》 SCIE EI CAS CSCD 2015年第6期86-94,共9页
The configurations of near space hypersonic flying vehicles are considerably different from those of conventional aircrafts.Their configurations are relatively slender;hence their moment of inertia around the longitud... The configurations of near space hypersonic flying vehicles are considerably different from those of conventional aircrafts.Their configurations are relatively slender;hence their moment of inertia around the longitudinal axis is much smaller than those around the other two axes,resulting in strong coupling of rotations around the three axes.Thus,the stability analysis of rolling motion for such flying vehicles is more complicated than those for conventional aircrafts,and there is no available result of stability analysis which can readily be applied to such cases.This paper is mainly concerned with the stated problem.Considering the practical situation,our investigation is targeted a slightly simpler problem,namely the rolling stability of flying vehicle under known pitching motion.The stability criterion of rolling motion is obtained with and without lateral motions.We also conducted numerical simulation for the pitching-rolling coupled motions of flying vehicles by solving Navier-Stokes equations coupled with dynamic equations of flight.The results of simulation agree well with those of theoretical analysis and experiments. 展开更多
关键词 STABILITY rolling motion pitching motion hypersonic flying vehicle
原文传递
RATIO-DEPENDENT PREDATOR-PREY MODEL WITH STAGE STRUCTURE AND TIME DELAY 被引量:2
16
作者 LIJUAN ZHA JING-AN CUI XUEYONG ZHOU 《International Journal of Biomathematics》 2012年第4期15-37,共23页
Ratio-dependent predator prey models are favored by many animal ecologists recently as more suitable ones for predator-prey interactions where predation involves searching process. In this paper, a ratio-dependent pre... Ratio-dependent predator prey models are favored by many animal ecologists recently as more suitable ones for predator-prey interactions where predation involves searching process. In this paper, a ratio-dependent predator prey model with stage structure and time delay for prey is proposed and analyzed. In this model, we only consider the stage structure of immature and mature prey species and not consider the stage structure of predator species. We assume that the predator only feed on the mature prey and the time for prey from birth to maturity represented by a constant time delay. At first, we investigate the permanence and existence of the proposed model and sufficient conditions are derived. Then the global stability of the nonnegative equilibria are derived. We also get the sufficient criteria for stability switch of the positive equilibrium. Finally, some numerical simulations are carried out for supporting the analytic results. 展开更多
关键词 RATIO-DEPENDENT stage structure time delay characteristic equation globally asymptotical stability.
原文传递
Spatial patterns created by cross-diffusion for a three-species food chain model
17
作者 Canrong Tian Zhi Ling Zhigui Lin 《International Journal of Biomathematics》 2014年第2期27-49,共23页
This paper deals with the stability analysis to a three-species food chain model with crossdiffusion, the results of which show that there is no Turing instability but crossdiffusion makes the model instability possib... This paper deals with the stability analysis to a three-species food chain model with crossdiffusion, the results of which show that there is no Turing instability but crossdiffusion makes the model instability possible. We then show that the spatial patterns are spotted patterns by using numerical simulations. In order to understand why the spatial patterns happen, the existence of the nonhomogeneous steady states is investigated. Finally, using the Leray-Schauder theory, we demonstrate that cross-diffusion creates nonhomogeneous stationary patterns. 展开更多
关键词 Spatial pattern Turing instability CROSS-DIFFUSION food chain.
原文传递
Mathematical model for smoking: Effect of determination and education 被引量:1
18
作者 Anuradha Yadav Prashant K. Srivastava Anuj Kumar 《International Journal of Biomathematics》 2015年第1期1-14,共14页
We propose and analyze mathematical models to study the dynamics of smoking behavior under the influence of educational programs and also individual's determination to quit smoking. We establish the positivity and bo... We propose and analyze mathematical models to study the dynamics of smoking behavior under the influence of educational programs and also individual's determination to quit smoking. We establish the positivity and boundedness of the solutions in a biologically feasible region. A threshold value responsible for persistence of smoking is obtained and stability analysis on models is performed. We find that determination alone is not enough to eradicate smoking but it can reduce the prevalence of smoker population. Whereas the increase in education can possibly eradicate it. We performed numerical simulation for representative set of parameters to verify and discuss results obtained analytically. 展开更多
关键词 SMOKING DETERMINATION EDUCATION global stability Lyapunov function
原文传递
A fractional-order toxin producing phytoplankton and zooplankton system 被引量:1
19
作者 M. Javidi N. Nyamoradi 《International Journal of Biomathematics》 2014年第4期63-83,共21页
In this work, we investigate the dynamical behavior of a fractional-order toxin producing on a phytoplankton-zooplankton (TPPZ) system with nutrient cycling. We propose a mathematical system to model this situation.... In this work, we investigate the dynamical behavior of a fractional-order toxin producing on a phytoplankton-zooplankton (TPPZ) system with nutrient cycling. We propose a mathematical system to model this situation. All the feasible equilibria of the system are obtained and the conditions for the existence of the equilibriums are determined. Local stability analysis of the TPPZ is studied by using the fractional Routh-Hurwitz stability conditions. Numerical simulations are carried out for a hypothetical set of parameter values to substantiate our analytical findings. 展开更多
关键词 Toxic-phytoplankton ZOOPLANKTON stability numerical simulation frac-tional-order.
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部