期刊文献+
共找到232篇文章
< 1 2 12 >
每页显示 20 50 100
人工冻土分数阶导数应力-应变指数模型参数确定及验证
1
作者 姚兆明 宋梓豪 +1 位作者 陈军浩 左维亚 《煤炭学报》 EI CAS CSCD 北大核心 2024年第S01期285-294,共10页
人工冻土可看成是理想固体和理想流体以某种比例进行的勾兑,其力学特性既不符合胡克定律,也不遵守牛顿黏性定律,而是遵守介于它们之间的某种关系,分数阶导数能够很好地描述这种勾兑效应。对合肥膨胀土在不同冻结温度下进行单轴压缩试验... 人工冻土可看成是理想固体和理想流体以某种比例进行的勾兑,其力学特性既不符合胡克定律,也不遵守牛顿黏性定律,而是遵守介于它们之间的某种关系,分数阶导数能够很好地描述这种勾兑效应。对合肥膨胀土在不同冻结温度下进行单轴压缩试验,得到冻结温度对应力-应变的影响规律。将分数阶导数引入指数模型,将它改进为人工冻土单轴压缩下应力-应变分数阶指数模型。通过对建立的模型两边取自然对数,得到不同温度下应力-应变线性方程组,求解建立方程组确定出分数阶导数模型参数。为进一步验证所建立模型的适用性,引用一组南京冻结粉质黏土三轴剪切试验,在分数阶系数中考虑围压的影响,将它改进为能考虑围压影响的应力-应变分数阶指数模型。将改进的人工冻土应力-应变分数阶指数模型的计算结果与试验结果进行比较,结果表明:计算结果与试验结果吻合度较高,能准确地预测单轴压缩、三轴压缩剪切应力-应变曲线的变化趋势。改进后的分数阶导数模型参数较少且有明确的物理意义,便于工程应用。当前的模型仅仅适用于应变硬化型,为了能进一步描述应变软化型的力学特性,下一步将在模型中考虑损伤进而建立应变软化型分数阶指数模型。同时,如何在模型中反映冻土的结构性与各向异性对应力-应变的影响也将是下一步的研究内容。 展开更多
关键词 人工冻土 指数模型 应力-应变曲线 分数导数
在线阅读 下载PDF
基于再生核和有限差分法求解变系数时间分数阶对流扩散方程
2
作者 吕学琴 何松岩 王世宇 《数学物理学报(A辑)》 北大核心 2025年第1期153-164,共12页
针对变系数的时间分数阶对流-扩散方程,首先,使用有限差分法,得到了该方程的半离散格式.之后再利用再生核方法,得到了方程的精确解u(x,t_(n)),将精确解u(x,t_(n))取m项截断,可得到近似解u_(m)(x,t_(n)).通过证明,得到该方法是稳定的.最... 针对变系数的时间分数阶对流-扩散方程,首先,使用有限差分法,得到了该方程的半离散格式.之后再利用再生核方法,得到了方程的精确解u(x,t_(n)),将精确解u(x,t_(n))取m项截断,可得到近似解u_(m)(x,t_(n)).通过证明,得到该方法是稳定的.最后,通过三个数值例子,并与其他文献中的方法在同等条件下进行了比较,证明该算法有效. 展开更多
关键词 CAPUTO分数导数 再生核方法 变系数时间分数对流扩散方程 有限差分方法
在线阅读 下载PDF
一类具有推广的Caputo分数阶导数的微分方程的Ulam稳定性
3
作者 张玲玲 王森 +1 位作者 张孝锋 周先锋 《应用数学》 北大核心 2024年第3期847-855,共9页
本文研究一类具有推广的Caputo分数阶导数的微分方程,推广的Caputo分数阶导数可由Conformable导数与经典的Caputo导数结合而得或是推广的分数阶算子.我们利用拉普拉斯变换研究了线性方程的Ulam-Hyers稳定性,分别利用Banach不动点定理和G... 本文研究一类具有推广的Caputo分数阶导数的微分方程,推广的Caputo分数阶导数可由Conformable导数与经典的Caputo导数结合而得或是推广的分数阶算子.我们利用拉普拉斯变换研究了线性方程的Ulam-Hyers稳定性,分别利用Banach不动点定理和Guonwall不等式研究了非线性方程解的存在唯一性和Ulam-Hyers-Rassis稳定性,获得了几个充分条件的定理,并给出一个例子作为所得结果的应用. 展开更多
关键词 Ulam稳定性 推广的Caputo分数导数 微分方程
在线阅读 下载PDF
基于分数阶导数模型的巷道围岩变形失稳规律研究 被引量:1
4
作者 肖书文 《晋控科学技术》 2024年第1期15-18,共4页
巷道在应力长时间作用下会发生蠕变变形,为探究蠕变变形失稳规律,以分数阶导数为基础,将Newton黏壶改进为Abel黏壶,并引入负指数形式的损伤变量,建立了分数阶Nishihara本构模型。根据晋城市锦辰煤业3采区砂岩在不同应力下的蠕变实验,验... 巷道在应力长时间作用下会发生蠕变变形,为探究蠕变变形失稳规律,以分数阶导数为基础,将Newton黏壶改进为Abel黏壶,并引入负指数形式的损伤变量,建立了分数阶Nishihara本构模型。根据晋城市锦辰煤业3采区砂岩在不同应力下的蠕变实验,验证了本构模型的有效性,所提出的新模型可以更好到描述巷道围岩在长期力学作用下的变形规律。 展开更多
关键词 蠕变 巷道变形 分数导数 本构方程
在线阅读 下载PDF
分数阶时间导数方程和反常亚扩散过程——纪念茆诗松教授
5
作者 陈振庆 《应用概率统计》 CSCD 北大核心 2024年第2期323-342,共20页
本文介绍并且改进和推广了时间导数为分数阶的方程,以及与反常亚扩散过程相关联的最近的一些结果.
关键词 分数时间导数 时间为分数的方程 隶属子 逆隶属子 强解和弱解
在线阅读 下载PDF
基于分数阶导数衍生特征和小波散射网络的交通碰撞声学信号识别
6
作者 赵国荣 祁晋峰 +2 位作者 武建军 熊威 徐秀 《电子设计工程》 2024年第20期181-186,共6页
针对现阶段交通碰撞事故报警系统实时性差、识别不准等问题,提出了一种联合分数阶导数衍生特征和小波散射网络系数的交通事故声学碰撞信号识别方法。针对车辆碰撞信号具备非平稳性和长程相关性,引入分数阶导数,从中提取可以表征车辆行... 针对现阶段交通碰撞事故报警系统实时性差、识别不准等问题,提出了一种联合分数阶导数衍生特征和小波散射网络系数的交通事故声学碰撞信号识别方法。针对车辆碰撞信号具备非平稳性和长程相关性,引入分数阶导数,从中提取可以表征车辆行驶状态的时频特征作为第一个特征数据集。采用小波散射网络提取具有多尺度、稳定性和局部不变性的散射系数作为第二个特征数据集。将这两个特征数据集进行融合,基于树模型构造交通碰撞事故声学信号识别模型。实验结果表明,在已有碰撞事故样本上的检出率为100%,误报次数控制在每天每公里0.5次以内。 展开更多
关键词 分数导数 车辆碰撞 分布式声学传感 小波散射网络
在线阅读 下载PDF
基于平均法的余弦激励分数阶振子系统动力学分析
7
作者 师玮 郭蓉 +3 位作者 解加全 张彦杰 王涛 黄庆学 《工程数学学报》 CSCD 北大核心 2024年第4期623-641,共19页
提出了一种求解余弦激励下分数阶振子系统位移响应的解析和数值算法。其中解析方法是用平均法来求出系统的稳态响应解和瞬态响应解,总位移响应解即为稳态解和瞬态解的总和。提出一种数值方法,利用Grunwald-Letnikov分数阶导数的定义将... 提出了一种求解余弦激励下分数阶振子系统位移响应的解析和数值算法。其中解析方法是用平均法来求出系统的稳态响应解和瞬态响应解,总位移响应解即为稳态解和瞬态解的总和。提出一种数值方法,利用Grunwald-Letnikov分数阶导数的定义将系统中的分数阶微分项进行离散化处理,降低了原系统阶数。在一般的周期激励下,系统的近似响应解可以采用傅里叶级数展开法和线性系统叠加原理得到。最后,利用数值模拟对所提方法的有效性和可行性进行了验证,并对分数阶阶次、线性阻尼系数和分数阶导数系数对系统稳态响应振幅和总位移响应的影响进行了分析。 展开更多
关键词 平均法 分数振子系统 grunwald-letnikov分数阶导数 位移响应
在线阅读 下载PDF
基于分数阶导数的黏弹性悬架减振模型及其数值方法 被引量:13
8
作者 李占龙 孙大刚 +2 位作者 宋勇 刘付喜 赵树萍 《振动与冲击》 EI CSCD 北大核心 2016年第16期123-129,共7页
为了准确掌握黏弹性悬架的动态响应,针对传统整数阶减振模型的不足,引入分数阶导数原理,构建了黏弹性材料FKV本构模型,建立了考虑几何参数的黏弹性悬架分数阶减振模型,利用Grumwald-Letnikov定义将模型中分数阶导数离散化,并转化为状态... 为了准确掌握黏弹性悬架的动态响应,针对传统整数阶减振模型的不足,引入分数阶导数原理,构建了黏弹性材料FKV本构模型,建立了考虑几何参数的黏弹性悬架分数阶减振模型,利用Grumwald-Letnikov定义将模型中分数阶导数离散化,并转化为状态方程形式,依据矩阵函数理论推导出模型的数值解。以某型安装黏弹性悬架的履带车辆参数为例,分别建立了悬架的动态接触有限元模型和分数阶减振模型,获得了在翻越障碍工况下两种模型响应的对比解。结果表明:分数阶减振模型体现了黏弹性悬架响应具有全局相关性和记忆性,且历史作用渐近加强;黏弹性悬架有较好的缓冲减振性能;分数阶减振模型解与有限元方法计算结果有较好的一致性。旨为下一步的实车试验和实际应用提供理论参考。 展开更多
关键词 分数导数 黏弹性悬架 减振模型 数值方法
在线阅读 下载PDF
人工冻土蠕变特性粒子群分数阶导数模型 被引量:24
9
作者 陈军浩 姚兆明 +1 位作者 徐颖 王厚良 《煤炭学报》 EI CAS CSCD 北大核心 2013年第10期1763-1768,共6页
针对传统整数阶微积分本构关系所需元件多等不足,将分数阶微积分关系运用到人工冻土蠕变计算中。在元件模型基础上,引入分数阶导数,建立稳定蠕变分数阶导数定常蠕变模型。运用粒子群优化方法,通过对不同加载应力下稳态蠕变试验数据建立... 针对传统整数阶微积分本构关系所需元件多等不足,将分数阶微积分关系运用到人工冻土蠕变计算中。在元件模型基础上,引入分数阶导数,建立稳定蠕变分数阶导数定常蠕变模型。运用粒子群优化方法,通过对不同加载应力下稳态蠕变试验数据建立的蠕变模型参数进行识别,得到蠕变模型的材料参数。通过对比分数阶导数蠕变模型计算值与试验结果,发现分数阶导数蠕变模型能很好地模拟人工冻土的蠕变规律。分数阶导数模型是人工冻土领域计算的一种新方法。 展开更多
关键词 人工冻土蠕变 粒子群 分数导数 KELVIN模型
在线阅读 下载PDF
具逐项分数阶导数的积分边值问题正解的存在性 被引量:10
10
作者 李燕 刘锡平 +1 位作者 李晓晨 张莎 《上海理工大学学报》 CAS 北大核心 2016年第6期511-516,共6页
研究了一类具有逐项分数阶导数的微分方程积分边值问题正解的存在性和多解性.利用锥上不动点定理和Leggett-Williams不动点定理,分别得到了该积分边值问题至少存在1个正解和3个正解的结论.最后给出2个例子来证明结论有效.
关键词 分数微分方程 逐项分数导数 积分边值问题 锥上不动点定理 正解
在线阅读 下载PDF
分数阶导数在地震奇异性分析中的应用 被引量:7
11
作者 宋建国 刘垒 +1 位作者 李辉 刘励耘 《石油物探》 EI CSCD 北大核心 2009年第1期72-75,共4页
传统的地震解释主要是观测地震资料的振幅及相位的变化,而振幅往往并不能反映真实的地质情况。地震界面可能是岩性分界面也可能是岩性过渡带,岩性过渡带的地震反射波是入射波的分数阶导数,因此将分数阶导数引入地震属性计算中,构建一种... 传统的地震解释主要是观测地震资料的振幅及相位的变化,而振幅往往并不能反映真实的地质情况。地震界面可能是岩性分界面也可能是岩性过渡带,岩性过渡带的地震反射波是入射波的分数阶导数,因此将分数阶导数引入地震属性计算中,构建一种对波形敏感而对振幅变化不敏感的新属性——奇异性,用以刻画反射界面的横向变化。方法的基本原理是:首先计算地震子波的不同分数阶导数;然后利用匹配追踪算法将地震数据分解成地震子波的不同分数阶导数,进而获得反射波同相轴的分数阶。对胜利油田某区块实际二维地震资料进行了试处理,结果表明分数阶导数剖面能很好地描述不整合面,反映实际界面的横向变化。 展开更多
关键词 分数导数 奇异性 匹配追踪 地震属性
在线阅读 下载PDF
复杂介质中扩散和耗散行为的分数阶导数唯象建模 被引量:9
12
作者 庞国飞 陈文 +1 位作者 张晓棣 孙洪广 《应用数学和力学》 CSCD 北大核心 2015年第11期1117-1134,共18页
复杂介质一般是多相混合物.与普通固体、液体和气体相比,其力学行为具有明显的记忆、路径依赖性特征,难以用一般的经典力学模型来描述,因而显得反常.从数学力学建模上看,整数阶导数的局部极限定义不适合描述这样的非局部力学行为.分数... 复杂介质一般是多相混合物.与普通固体、液体和气体相比,其力学行为具有明显的记忆、路径依赖性特征,难以用一般的经典力学模型来描述,因而显得反常.从数学力学建模上看,整数阶导数的局部极限定义不适合描述这样的非局部力学行为.分数阶导数实质上是微分-积分算子,能精确地刻画力学行为的全局相关特征.而且分数阶模型具有明确的统计物理解释.20世纪末至今,复杂介质反常力学行为的分数阶导数模型由于具有参数少,且参数的物理意义明确等突出优点,开始引起广泛关注.该文从唯象建模的角度,综述了分数阶导数和分形导数在复杂介质的反常扩散和频率依赖能量耗散建模中的应用与发展. 展开更多
关键词 复杂介质 分数导数 分形导数 扩散 耗散 数学力学建模
在线阅读 下载PDF
采用分数阶导数描述胶凝原油的流变模型 被引量:10
13
作者 王志方 张国忠 刘刚 《中国石油大学学报(自然科学版)》 EI CAS CSCD 北大核心 2008年第2期114-118,共5页
提出了采用分数阶导数的新流变模型理论来描述胶凝原油蠕变特性和动态粘弹特性的方法。由Maxwell体和分数阶导数的类Kelvin体组合建立了类Burgers体流变模型,并推导得到了分数阶导数类Burgers体的蠕变柔量表达式。试验验证了该模型能精... 提出了采用分数阶导数的新流变模型理论来描述胶凝原油蠕变特性和动态粘弹特性的方法。由Maxwell体和分数阶导数的类Kelvin体组合建立了类Burgers体流变模型,并推导得到了分数阶导数类Burgers体的蠕变柔量表达式。试验验证了该模型能精确描述胶凝原油的蠕变过程,而经典的标准粘弹性固体模型和粘弹性流体模型都不能很好地拟合蠕变试验曲线。试验还验证了分数阶导数的类Kelvin体和类标准线性体都能精确描述胶凝原油的动态粘弹性能,而经典的标准粘弹性固体模型同样也不能很好地拟合动态粘弹试验曲线。 展开更多
关键词 分数导数 胶凝原油 粘弹性 流变模型 蠕变模型
在线阅读 下载PDF
任意荷载下分数阶导数黏弹性饱和土体一维固结 被引量:10
14
作者 汪磊 孙德安 +1 位作者 解益 李培超 《岩土工程学报》 EI CAS CSCD 北大核心 2017年第10期1823-1831,共9页
将分数阶导数理论引入Kelvin–Voigt模型,来描述黏弹性饱和土体的力学行为。对饱和土体一维固结方程和分数阶导数Kelvin–Voigt本构方程实施Laplace变换,联立求解得到变换域内有效应力和沉降的解析解。采用Crump方法实现Laplace数值反演... 将分数阶导数理论引入Kelvin–Voigt模型,来描述黏弹性饱和土体的力学行为。对饱和土体一维固结方程和分数阶导数Kelvin–Voigt本构方程实施Laplace变换,联立求解得到变换域内有效应力和沉降的解析解。采用Crump方法实现Laplace数值反演,获得了任意荷载情况下物理空间内一维固结问题的半解析解。并将指数荷载情况下分数阶导数模型退化到黏弹性情形,结果与已有文献解析解相同,验证了本研究提出任意荷载情况下分数阶导数黏弹性解的可靠性。最后,分析了相关参数对固结沉降的影响。研究表明,任意荷载情形下分数阶导数黏弹性饱和土体一维固结发展过程与黏滞系数和分数阶次有关,分数阶次越大,固结沉降发展越快;黏滞系数越大,固结沉降变化越慢;荷载变化趋势与由荷载参数变化引起的沉降变化规律是一致的,且最终沉降量一致。本研究有助于深入认识分数阶黏弹性饱和土体的固结行为。 展开更多
关键词 任意荷载 分数导数 黏弹性 饱和土体 一维固结
在线阅读 下载PDF
盐岩的分数阶导数蠕变模型 被引量:12
15
作者 吴斐 刘建锋 +1 位作者 边宇 周志威 《四川大学学报(工程科学版)》 EI CSCD 北大核心 2014年第5期22-27,共6页
盐岩蠕变特性是影响能源深部地下储存库长期运营的重要决定因素。利用基于分数阶导数定义得到的Abel粘壶以及带应变触发的Abel粘壶建立新的分数阶导数蠕变模型。通过4个月的盐岩蠕变试验,获得全过程蠕变曲线,基于盐岩蠕变试验数据,对分... 盐岩蠕变特性是影响能源深部地下储存库长期运营的重要决定因素。利用基于分数阶导数定义得到的Abel粘壶以及带应变触发的Abel粘壶建立新的分数阶导数蠕变模型。通过4个月的盐岩蠕变试验,获得全过程蠕变曲线,基于盐岩蠕变试验数据,对分数阶导数蠕变模型的参数进行拟合分析,确定模型中的参数,并对建立的分数阶导数蠕变模型与西原正夫模型进行对比分析。结果表明:分数阶导数蠕变模型与试验数据吻合得更好,不仅克服了整数阶蠕变模型在初始阶段与试验数据不能很好吻合的弊端,又能反映盐岩全过程蠕变特征。 展开更多
关键词 盐岩 分数导数蠕变模型 西原正夫模型 Abel粘壶
在线阅读 下载PDF
一类含有分数阶导数的参数激励振动问题 被引量:5
16
作者 葛志新 陈咸奖 陈松林 《振动与冲击》 EI CSCD 北大核心 2017年第4期88-92,共5页
研究了一类具有分数阶导数阻尼的参数激励振动问题。对含有由Riemann-Liouville定义的分数阶导数的Mathieu振动方程构造渐近解。利用多重尺度法,在激励参数取不同值的情况下,求得渐近解,得到分数阶指数对解的影响。
关键词 多重尺度 分数导数 参数激励 过渡曲线
在线阅读 下载PDF
宽带噪声激励下含分数阶导数的van der Pol-Duffing振子的可靠性 被引量:4
17
作者 陈林聪 李海锋 +1 位作者 梅真 朱位秋 《西南交通大学学报》 EI CSCD 北大核心 2014年第1期45-51,共7页
为了研究宽带噪声激励下含分数阶导数的van der Pol-Duffing振子的首次穿越问题,首先应用广义谐波平衡技术,将分数阶导数表示的回复力分解为等效拟线性阻尼力和拟线性回复力,获得不含分数阶导数的等效非线性随机系统;然后,应用随机平均... 为了研究宽带噪声激励下含分数阶导数的van der Pol-Duffing振子的首次穿越问题,首先应用广义谐波平衡技术,将分数阶导数表示的回复力分解为等效拟线性阻尼力和拟线性回复力,获得不含分数阶导数的等效非线性随机系统;然后,应用随机平均法将等效非线性随机系统近似为一维扩散过程,再建立和求解相应的后向Kolmogorov方程,获得系统的条件可靠性函数和平均首次穿越时间计算式;最后,通过实验结果表明,所提方法与蒙特卡罗法模拟结果吻合得非常好;系统的可靠性随分数阶数的增加而提高;分数阶导数表示的回复力不能简单地当作一类特殊的阻尼力. 展开更多
关键词 分数导数 可靠性 随机平均法 van der Pol—Duffing振子 宽带噪声
在线阅读 下载PDF
基于蚁群算法的冻结重塑黏土分数阶导数西原模型分析 被引量:4
18
作者 姚兆明 张秋瑾 牛连僧 《长江科学院院报》 CSCD 北大核心 2016年第7期81-86,共6页
掌握冻结状态下岩土体的蠕变规律,对利用冻结法施工的井筒建设安全至关重要。对人工冻结重塑黏土在-5,-10,-15℃下进行单轴抗压强度试验,得到温度对人工冻结重塑黏土单轴抗压强度的影响规律;根据重塑黏土冻结下的单轴强度,分别进行3个... 掌握冻结状态下岩土体的蠕变规律,对利用冻结法施工的井筒建设安全至关重要。对人工冻结重塑黏土在-5,-10,-15℃下进行单轴抗压强度试验,得到温度对人工冻结重塑黏土单轴抗压强度的影响规律;根据重塑黏土冻结下的单轴强度,分别进行3个加载等级的蠕变试验,得到温度、加载等级对蠕变的影响规律。将Abel黏壶引入到西原模型,建立分数阶导数西原模型;利用蚁群算法对重塑黏土冻结状态下蠕变西原模型和分数阶西原模型进行参数辨识,通过分析2种模型的模拟结果,表明分数阶西原模型更适合于计算重塑黏土冻结状态下的蠕变规律。 展开更多
关键词 冻结重塑黏土 蠕变 蚁群算法 分数导数 西原模型
在线阅读 下载PDF
分数阶导数粘弹性模型的有限元法 被引量:8
19
作者 银花 陈宁 《计算力学学报》 EI CAS CSCD 北大核心 2012年第6期966-971,共6页
在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数... 在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。 展开更多
关键词 分数导数 粘弹性 本构模型 有限元法
在线阅读 下载PDF
基于分数阶导数的沥青胶浆动态力学特性 被引量:4
20
作者 李锐铎 乐金朝 +1 位作者 李瑞霞 姚立阳 《振动与冲击》 EI CSCD 北大核心 2016年第15期46-49,92,共5页
通过对粉胶比为0.62、0.82、1.02、1.22和1.42的沥青胶浆在20℃、30℃、40℃和50℃条件下进行动态频率扫描试验,研究了不同粉胶比及试验温度条件下沥青胶浆复模量、抗车辙因子和相位角的变化规律。基于分数阶导数理论,建立了Nutting蠕... 通过对粉胶比为0.62、0.82、1.02、1.22和1.42的沥青胶浆在20℃、30℃、40℃和50℃条件下进行动态频率扫描试验,研究了不同粉胶比及试验温度条件下沥青胶浆复模量、抗车辙因子和相位角的变化规律。基于分数阶导数理论,建立了Nutting蠕变方程与经典分数阶导数Abel黏壶蠕变模型之间的关系,从而明确了Nutting蠕变方程各参数的物理意义。对分数阶Riemann-Liouville算子黏弹性蠕变本构模型的动态力学响应进行分析,提出了利用沥青胶浆动态频率扫描试验结果确定蠕变本构模型中参数A值和γ值的新方法。 展开更多
关键词 沥青胶浆 粉胶比 动态剪切流变试验 分数导数 Nutting蠕变方程 变形因子 黏弹因子
在线阅读 下载PDF
上一页 1 2 12 下一页 到第
使用帮助 返回顶部