期刊文献+
共找到44篇文章
< 1 2 3 >
每页显示 20 50 100
多自由度强非线性颤振分析的增量谐波平衡法 被引量:19
1
作者 蔡铭 刘济科 李军 《应用数学和力学》 EI CSCD 北大核心 2006年第7期833-838,共6页
对多个自由度上含有强非线性项系统的颤振问题,推广应用增量谐波平衡法进行分析.考虑带有强非线性立方平移和俯仰刚度项的二元机翼颤振方程,首先将方程用矩阵形式表示,然后把振动过程分解成为振动瞬态的持续增量过程,再采用振幅作为控... 对多个自由度上含有强非线性项系统的颤振问题,推广应用增量谐波平衡法进行分析.考虑带有强非线性立方平移和俯仰刚度项的二元机翼颤振方程,首先将方程用矩阵形式表示,然后把振动过程分解成为振动瞬态的持续增量过程,再采用振幅作为控制参数应用谐波平衡法,以这种推广的增量谐波平衡法求得方程解的表达式,并由此分析系统的分岔现象、极限环颤振现象和谐波项数的取值问题,最后用龙格_库塔数值方法进行验算,结果表明:分析多个自由度的强非线性颤振,增量谐波平衡法是精确有效的. 展开更多
关键词 强非线性颤振 增量谐波平衡 分岔 极限环
在线阅读 下载PDF
基于增量谐波平衡法的汽车转向轮非线性摆振的研究 被引量:8
2
作者 刘献栋 张紫广 +1 位作者 何田 单颖春 《汽车工程》 EI CSCD 北大核心 2011年第2期142-147,共6页
采用振幅作为控制参数,应用增量谐波平衡法,求得方程解的表达式,分析了系统的分岔和极限环振动现象,并用龙格库塔法进行数值计算验证,结果表明:增量谐波平衡法在汽车转向轮摆振分析中有效可靠,有足够的精度,并运用Floquet理论对周期解... 采用振幅作为控制参数,应用增量谐波平衡法,求得方程解的表达式,分析了系统的分岔和极限环振动现象,并用龙格库塔法进行数值计算验证,结果表明:增量谐波平衡法在汽车转向轮摆振分析中有效可靠,有足够的精度,并运用Floquet理论对周期解的稳定性进行分析,进而确定了极限环的稳定性,最后分析了各结构参数对摆振的影响。 展开更多
关键词 汽车转向轮 摆振 极限环 分岔 增量谐波平衡
在线阅读 下载PDF
基于增量谐波平衡法的复合行星齿轮传动系统非线性动力学 被引量:19
3
作者 刘振皓 巫世晶 +1 位作者 王晓笋 潜波 《振动与冲击》 EI CSCD 北大核心 2012年第3期117-122,共6页
建立了包含时变啮合刚度、齿侧间隙与综合啮合误差的Ravigneaux式复合行星齿轮传动系统纯扭转动力学模型。运用增量谐波平衡法对系统运动微分方程组进行求解,得到系统的基频稳态响应。研究了时变啮合刚度、外部激励、齿侧间隙等参数的... 建立了包含时变啮合刚度、齿侧间隙与综合啮合误差的Ravigneaux式复合行星齿轮传动系统纯扭转动力学模型。运用增量谐波平衡法对系统运动微分方程组进行求解,得到系统的基频稳态响应。研究了时变啮合刚度、外部激励、齿侧间隙等参数的变化对系统动力学特性的影响。研究结果表明,间隙的存在使得复合行星齿轮系统的频响曲线出现了幅值跳跃与多值解等典型非线性特征,系统参数的共同作用使得复合行星齿轮系统出现了丰富的非线性动力学行为。利用该方法可以获得系统任意精度的近似解,为控制系统的振动与噪声,实现复合行星齿轮传动系统动态设计奠定基础。 展开更多
关键词 复合行星齿轮 传动系统 非线性动力学 增量谐波平衡
在线阅读 下载PDF
基于增量谐波平衡法的流量调节器非线性频率特性 被引量:6
4
作者 刘上 刘红军 +1 位作者 孙宏明 陈宏玉 《推进技术》 EI CAS CSCD 北大核心 2012年第5期814-819,共6页
针对多个状态变量、非线性关系复杂、方程刚性强的流量调节器系统,应用增量谐波平衡法研究了调节器在入口压力扰动下的强迫振荡问题,得到了系统周期解的解析表达式。通过与龙格库塔数值仿真结果相比较,验证了增量谐波平衡法应用于该系... 针对多个状态变量、非线性关系复杂、方程刚性强的流量调节器系统,应用增量谐波平衡法研究了调节器在入口压力扰动下的强迫振荡问题,得到了系统周期解的解析表达式。通过与龙格库塔数值仿真结果相比较,验证了增量谐波平衡法应用于该系统的有效性,且计算量相对较小。进而在频域内探讨了流量调节器系统的非线性频响特性,结果表明,随着扰动振幅的增大,幅频响应在整个频率范围上都有所增大。计算结果可以用于评估线性化方法的误差范围。通过分析周期解的时均值与稳态值的差别,揭示了调节器受上游压力周期扰动时,其出口流量的时均值略低于稳态值的非线性特性。 展开更多
关键词 流量调节器 增量谐波平衡 周期解 非线性频率特性
在线阅读 下载PDF
行星齿轮系统弯扭耦合振动的增量谐波平衡法 被引量:5
5
作者 姚红良 许琦 +1 位作者 茅列前 闻邦椿 《东北大学学报(自然科学版)》 EI CAS CSCD 北大核心 2013年第10期1451-1455,共5页
针对行星齿轮传动机构弯扭耦合振动方程是半正定方程而无法直接采用增量谐波平衡法进行求解的问题,引入微小相对位移量来描述行星齿轮系统中各个部件之间的相对运动.采用运动转移矩阵将半正定方程转化为正定方程,同时考虑齿轮啮合刚度... 针对行星齿轮传动机构弯扭耦合振动方程是半正定方程而无法直接采用增量谐波平衡法进行求解的问题,引入微小相对位移量来描述行星齿轮系统中各个部件之间的相对运动.采用运动转移矩阵将半正定方程转化为正定方程,同时考虑齿轮啮合刚度时变的非线性特性,采用增量谐波平衡法分析了行星齿轮系统非线性动力学响应.通过与数值求解方法Newmark-β法相对比,计算结果一致,验证了本文方法的正确性及高效性. 展开更多
关键词 行星齿轮 弯扭耦合振动 增量谐波平衡 非线性 动力学响应
在线阅读 下载PDF
基于增量谐波平衡法的齿轮系统非线性动力学 被引量:26
6
作者 杨绍普 申永军 刘献栋 《振动与冲击》 EI CSCD 北大核心 2005年第3期40-42,95,i004,共5页
利用增量谐波平衡法研究了一类考虑时变啮合刚度和间隙的直齿轮副的非线性动力学,建立了这类模型的解的统一形式。同时研究了系统的分岔特性以及阻尼比、外激励幅值对系统幅频曲线的影响。利用本文结果可以得到任意精度的近似解,从而可... 利用增量谐波平衡法研究了一类考虑时变啮合刚度和间隙的直齿轮副的非线性动力学,建立了这类模型的解的统一形式。同时研究了系统的分岔特性以及阻尼比、外激励幅值对系统幅频曲线的影响。利用本文结果可以得到任意精度的近似解,从而可以更好地分析和控制齿轮系统的动力学。最后,指出了其它文献中未注意到的问题。 展开更多
关键词 增量谐波平衡 非线性动力学 齿轮系统 时变啮合刚度 直齿轮副 统一形式 分岔特性 幅频曲线 类模型 阻尼比 外激励 近似解
在线阅读 下载PDF
基于增量谐波平衡法的含索铰可折展桁架非线性动力学特性 被引量:3
7
作者 张静 刘荣强 +1 位作者 郭宏伟 邓宗全 《振动与冲击》 EI CSCD 北大核心 2014年第7期4-10,36,共8页
为了揭示含索铰可折展桁架的非线性动力学行为,建立了考虑铰链间隙、刚度和阻尼及索非线性特性的可折展桁架纯弯曲动力学模型。对非线性动力学方程进行一次泰勒展开和参量的多次谐波描述,实现了非线性动力学方程到代数方程的转化,通过... 为了揭示含索铰可折展桁架的非线性动力学行为,建立了考虑铰链间隙、刚度和阻尼及索非线性特性的可折展桁架纯弯曲动力学模型。对非线性动力学方程进行一次泰勒展开和参量的多次谐波描述,实现了非线性动力学方程到代数方程的转化,通过迭代进行非线性动力学系统的响应计算。并利用龙格库塔方法对非线性系统进行数值分析,与增量谐波平衡(IHB)法进行对比,验证了IHB法计算的正确性。以激振频率为变化参数,对悬臂支撑的含索铰桁架结构进行解的稳定性分析,得到铰链间隙、铰链刚度、激振力和索对结构响应稳定性的影响。基于IHB法可快速准确的进行多自由度可折展结构动力学求解,为研究大型折展桁架的动力学行为奠定了基础。 展开更多
关键词 可折展桁架 非线性特性 增量谐波平衡 铰链
在线阅读 下载PDF
优化迭代步长的两种改进增量谐波平衡法 被引量:3
8
作者 黄建亮 张兵许 陈树辉 《力学学报》 EI CAS CSCD 北大核心 2022年第5期1353-1363,共11页
增量谐波平衡法(IHB法)是一个半解析半数值的方法,其最大优点是适合于强非线性系统振动的高精度求解.然而,IHB法与其他数值方法一样,也存在如何选择初值的问题,如初值选择不当,会存在不收敛的情况.针对这一问题,本文提出了两种基于优化... 增量谐波平衡法(IHB法)是一个半解析半数值的方法,其最大优点是适合于强非线性系统振动的高精度求解.然而,IHB法与其他数值方法一样,也存在如何选择初值的问题,如初值选择不当,会存在不收敛的情况.针对这一问题,本文提出了两种基于优化算法的IHB法:一是结合回溯线搜索优化算法(BLS)的改进IHB法(GIHB1),用来调节IHB法的迭代步长,使得步长逐渐减小满足收敛条件;二是引入狗腿算法的思想并结合BLS算法的改进IHB法(GIHB2),在牛顿-拉弗森(Newton-Raphson)迭代中引入负梯度方向,并在狗腿算法中引入2个参数来调节BSL搜索方式用于调节迭代的方式,使迭代方向沿着较快的下降方向,从而减少迭代的步数,提升收敛的速度.最后,给出的两个算例表明两种改进IHB法在解决初值问题上的有效性. 展开更多
关键词 非线性振动 增量谐波平衡 初值问题 回溯线搜索 狗腿
在线阅读 下载PDF
多跨碰摩转子系统分析的降维增量谐波平衡法 被引量:1
9
作者 姚红良 丛艳 +1 位作者 许琦 闻邦椿 《东北大学学报(自然科学版)》 EI CAS CSCD 北大核心 2015年第7期991-995,共5页
多跨碰摩转子系统所建有限元动力学模型自由度很多,严重降低响应预测时的计算效率.为了提高计算效率,同时保证计算精度,提出结合固定界面模态综合法的降维增量谐波平衡法.首先利用固定界面模态综合法降低多跨碰摩转子系统的自由度数目,... 多跨碰摩转子系统所建有限元动力学模型自由度很多,严重降低响应预测时的计算效率.为了提高计算效率,同时保证计算精度,提出结合固定界面模态综合法的降维增量谐波平衡法.首先利用固定界面模态综合法降低多跨碰摩转子系统的自由度数目,再利用增量谐波平衡法求解降维后系统的稳态响应,最后还原回原系统的响应.该方法不仅能够直接获得稳态解,还可以与弧长法相结合从而求得非稳定解.对双跨转子系统进行仿真试验,结果表明该方法在主模态数较小时即可保证计算精度并且具有很高的效率. 展开更多
关键词 转子系统 碰摩 固定界面模态综合 增量谐波平衡 弧长
在线阅读 下载PDF
基于增量谐波平衡法的杜芬方程非线性区间多值自动求解算法 被引量:1
10
作者 张纯金 任勇生 纪淑娟 《山东科技大学学报(自然科学版)》 CAS 2015年第5期92-97,共6页
以典型的非线性振动系统杜芬方程为例,研究应用增量谐波平衡法自动求解非线性方程区间内多值的问题。在非线性区间内选择一个与幅频特性曲线趋向相似的谐波项系数作为主动变量,根据振幅的变化情况,自动选择不同谐波项系数,直到求出所有... 以典型的非线性振动系统杜芬方程为例,研究应用增量谐波平衡法自动求解非线性方程区间内多值的问题。在非线性区间内选择一个与幅频特性曲线趋向相似的谐波项系数作为主动变量,根据振幅的变化情况,自动选择不同谐波项系数,直到求出所有的值,所求出的值与四阶Runge-Kutta法的结果相比精确度相当,验证了算法的正确性。 展开更多
关键词 增量谐波平衡 非线性振动 杜芬方程 幅频响应
在线阅读 下载PDF
基于增量谐波平衡法的输流管道非线性振动分析(英文) 被引量:1
11
作者 孟丹 陈亮 《船舶力学》 EI 北大核心 2011年第12期1416-1428,共13页
考虑内、外流体的影响,利用Kane方法和拉格朗日应变理论建立了输流管道的二维非线性动力学模型。将动力学模型在平衡位置附近线性化,进行输流管道的线性稳定性分析。采用增量谐波平衡法求解包含三次非线性项动力学模型的稳态周期解。探... 考虑内、外流体的影响,利用Kane方法和拉格朗日应变理论建立了输流管道的二维非线性动力学模型。将动力学模型在平衡位置附近线性化,进行输流管道的线性稳定性分析。采用增量谐波平衡法求解包含三次非线性项动力学模型的稳态周期解。探讨了不同参数影响下输流管道的非线性动力响应以及内流流速的大小对管道动力特性的影响。 展开更多
关键词 非线性模型 KANE方 振动分析 增量谐波平衡
在线阅读 下载PDF
复杂矿山设备周期激励响应分析的增量谐波平衡法
12
作者 赵倩 张文兵 +1 位作者 姚红良 闻邦椿 《东北大学学报(自然科学版)》 EI CAS CSCD 北大核心 2014年第11期1588-1591,1601,共5页
针对复杂矿山设备中由于存在周期激励而引发的严重振动问题,以大型矿用给料装置为研究对象,将其钢结构机架简化为框架结构模型;考虑物料冲击是引发振动的主要原因,推导了非简谐周期激励的增量谐波平衡法;采用此方法分析周期激励下矿山... 针对复杂矿山设备中由于存在周期激励而引发的严重振动问题,以大型矿用给料装置为研究对象,将其钢结构机架简化为框架结构模型;考虑物料冲击是引发振动的主要原因,推导了非简谐周期激励的增量谐波平衡法;采用此方法分析周期激励下矿山设备的振动特性,并与Newmark-β法进行对比,验证了算法的精确性效率;研究设备与地面或基础接触时产生的非线性动响应,以及不同地面软硬程度下响应的规律.结果表明,该法用来分析复杂矿山设备在周期激励作用下的非线性稳态响应问题是有效的. 展开更多
关键词 矿山设备 框架结构 周期激励 稳态响应 增量谐波平衡
在线阅读 下载PDF
用IHB法分析双谐波激励下铰接塔-油轮系统的非线性动力学特性 被引量:1
13
作者 谢文会 唐友刚 《应用力学学报》 EI CAS CSCD 北大核心 2008年第3期498-502,共5页
研究了考虑平方阻尼情况下,铰接塔-油轮系统在双谐波激励下的非线性动力学特性。将该系统简化为单自由度分段线性恢复力,含平方阻尼的运动学分析模型,建立了铰接装载塔系统的分段非线性动力学方程。采用增量谐波平衡法获得系统周期解,使... 研究了考虑平方阻尼情况下,铰接塔-油轮系统在双谐波激励下的非线性动力学特性。将该系统简化为单自由度分段线性恢复力,含平方阻尼的运动学分析模型,建立了铰接装载塔系统的分段非线性动力学方程。采用增量谐波平衡法获得系统周期解,使用Floquet理论判断系统的运动稳定性,结合路径跟踪法跟踪系统响应曲线,获得了系统所有可能的亚谐、谐波、组合谐波共振运动。分析了不对称恢复刚度比值对系统亚谐、组合谐波共振和对系统运动倍周期分岔点的影响,比较了考虑平方阻尼和不考虑平方阻尼情况下系统非线性动力学特性,得到了系统的一些重要的非线性动力学特点。 展开更多
关键词 铰接塔 分段非线性 增量谐波平衡 平方阻尼 组合共振 亚谐共振
在线阅读 下载PDF
两端固定载流管非线性振动IHB方法研究
14
作者 王鹏 张咏鸥 +1 位作者 王晟 张涛 《振动与冲击》 EI CSCD 北大核心 2017年第20期248-253,共6页
增量平衡谐波法(IHB法)可用于求解两端简支载流管的非线性振动问题。考虑非线性约束及集中质量点的影响,利用Hamilton原理建立两端简支载流管运动微分方程,经Galerkin离散后,通过改变控制变量频率比得到系统的幅频特性曲线。讨论了系统... 增量平衡谐波法(IHB法)可用于求解两端简支载流管的非线性振动问题。考虑非线性约束及集中质量点的影响,利用Hamilton原理建立两端简支载流管运动微分方程,经Galerkin离散后,通过改变控制变量频率比得到系统的幅频特性曲线。讨论了系统运动参数例如速度、质量比、非线性约束刚度及质量点对系统幅频特性的影响。计算结果表明增量平衡谐波法是一种求解载流管非线性振动较为有效的方法。 展开更多
关键词 载流管 增量平衡谐波法 非线性约束 质量点
在线阅读 下载PDF
梁的强非线性超、次谐波共振 被引量:8
15
作者 黄建亮 黄惠仪 +1 位作者 陈恒 陈树辉 《振动与冲击》 EI CSCD 北大核心 2004年第1期1-3,共3页
本文采用改进的L -P法研究梁的强非线性振动 ,分析了超谐波共振和次谐波共振两种情况。首先建立了梁振动的Duffing方程 ,然后采用改进的L -P法分别求解超谐波响应和次谐波响应 ,最后给出两个典型算例 ,其计算结果与增量谐波平衡法和经典... 本文采用改进的L -P法研究梁的强非线性振动 ,分析了超谐波共振和次谐波共振两种情况。首先建立了梁振动的Duffing方程 ,然后采用改进的L -P法分别求解超谐波响应和次谐波响应 ,最后给出两个典型算例 ,其计算结果与增量谐波平衡法和经典的L -P法进行比较。算例表明 ,改进的L -P法得出的计算值与增量谐波平衡法所得出的结果非常吻合 ,但是经典的L -P法就失效了。 展开更多
关键词 谐波共振 谐波共振 DUFFING方程 L-P 增量谐波平衡
在线阅读 下载PDF
脉动流激励下输流管道的参数共振IHB方法研究 被引量:7
16
作者 梁峰 杨晓东 闻邦椿 《振动与冲击》 EI CSCD 北大核心 2008年第9期44-46,共3页
应用增量谐波平衡法(IHB法)研究脉动流激励下两端铰支输流管的参数共振。对无量纲的运动微分方程采用两振型Galerkin展开式离散化后,用IHB法分析了管道一阶1/2次谐波参数共振,并得到了幅-频响应曲线。将其与数值模拟解以及平均法结果比... 应用增量谐波平衡法(IHB法)研究脉动流激励下两端铰支输流管的参数共振。对无量纲的运动微分方程采用两振型Galerkin展开式离散化后,用IHB法分析了管道一阶1/2次谐波参数共振,并得到了幅-频响应曲线。将其与数值模拟解以及平均法结果比较发现,IHB法所得到的响应比平均法更加精确。研究同时表明:高阶模态对管道的一阶参数共振影响很小。 展开更多
关键词 输流管道 脉动流 增量谐波平衡(IHB) 参数共振
在线阅读 下载PDF
干摩擦LQ控制系统在简谐激励下响应的递推算法 被引量:4
17
作者 文明 邓子辰 《振动与冲击》 EI CSCD 北大核心 2004年第1期41-44,共4页
本文在用增量谐波平衡法求解系统响应和用子结构法求解Riccati方程的基础上 ,进一步研究了单自由度含有粘性阻尼双线性滞迟恢复力LQ控制系统在简谐激励下的响应计算问题 ,并给出相应算例 ,为干摩擦系统在工程中的应用提供了理论依据。
关键词 干摩擦LQ控制系统 简谐激励 递推算 增量谐波平衡 RICCATI方程 自由度 振动控制
在线阅读 下载PDF
粘弹性板的非线性动力稳定性分析方法 被引量:3
18
作者 彭凡 傅衣铭 《动力学与控制学报》 2005年第4期60-65,共6页
将微分积分型参数振动方程组转化成微分型,且基于增量谐波平衡法的一般应用途径,分析了受面内周期激励的粘弹性板的非线性动力稳定特性,揭示了主要动力不稳定区域的整体下移以及缩小和标准线性固体材料的粘性参数、板的振动频率之间的关... 将微分积分型参数振动方程组转化成微分型,且基于增量谐波平衡法的一般应用途径,分析了受面内周期激励的粘弹性板的非线性动力稳定特性,揭示了主要动力不稳定区域的整体下移以及缩小和标准线性固体材料的粘性参数、板的振动频率之间的关系.同时给出了增量谐波平衡法直接应用于非线性微分积分型参数振动方程的简化途径,并通过两种应用途径所得结果的对比,检验了这种简化途径的有效性. 展开更多
关键词 粘弹性板 非线性 动力稳定性 增量谐波平衡
在线阅读 下载PDF
IHB法在多自由度Bouc-Wen滞回非线性系统响应特性研究中的应用 被引量:4
19
作者 赵倩 刘子良 +1 位作者 姚红良 闻邦椿 《振动与冲击》 EI CSCD 北大核心 2018年第10期57-62,93,共7页
工程中常用Bouc-Wen模型来描述具有滞回特性的振动系统,此类系统是一种多值性的非解析系统,其动力学理论分析比较困难。由于Bouc-Wen滞回模型的微分形式,一般采用数值方法进行积分求解,但对于多自由度系统来说求解速度非常慢,且难以求... 工程中常用Bouc-Wen模型来描述具有滞回特性的振动系统,此类系统是一种多值性的非解析系统,其动力学理论分析比较困难。由于Bouc-Wen滞回模型的微分形式,一般采用数值方法进行积分求解,但对于多自由度系统来说求解速度非常慢,且难以求得不稳定解。故提出将滞回力引入为一个增加的自由度,重新建立振动系统的微分方程,将增量谐波平衡(IHB)法推广至求解该类含Bouc-Wen模型的多自由度滞回非线性系统,并引入弧长法解决由迟滞非线性引起的跳跃和多映射现象。利用该法分析了一些滞回系统的响应特性,通过与数值方法进行精度和效率对比,体现了该方法的优越性。 展开更多
关键词 滞回非线性 BOUC-WEN模型 悬臂梁 增量谐波平衡 弧长
在线阅读 下载PDF
准周期响应对称破缺分岔点的一种快速计算方法 被引量:1
20
作者 龚冰清 郑泽昌 +1 位作者 陈衍茂 刘济科 《力学学报》 EI CAS CSCD 北大核心 2022年第11期3181-3188,共8页
稳态响应如周期及准周期解的分岔计算,是非线性动力学研究的难点问题之一.与计算方法及分析理论相对完善的周期响应相比,准周期响应的求解只是在近些年才得到较大进展,而且其分岔分析更加棘手,仍需要更有效的理论和方法.目前,稳态响应... 稳态响应如周期及准周期解的分岔计算,是非线性动力学研究的难点问题之一.与计算方法及分析理论相对完善的周期响应相比,准周期响应的求解只是在近些年才得到较大进展,而且其分岔分析更加棘手,仍需要更有效的理论和方法.目前,稳态响应尤其是准周期响应的分岔计算,一般需采用数值方法,通过调节参数反复试算得到.为此,本文基于增量谐波平衡IHB法提出一种快速方法,可以高效地确定准周期响应的对称破缺分岔点.方法的理论基础是在准周期解的广义谐波级数表达基础上,当响应发生对称破缺分岔时,其偶次(含零次)谐波系数将逐渐由0变为小量.基于此性质,将零次谐波系数预先设定为小量,同时将分岔控制参数视为可变的迭代变量,进而通过IHB法构造迭代格式.作为算例,研究不可约频率作用下的双频激励Duffing系统以及Duffing-van der Pol耦合系统.结果表明,只要迭代格式收敛,随着预设小量减小,控制参数将逐渐接近分岔近似值;同时,通过提高谐波截断数可显著提高近似分岔值的计算精度.所提方法无需反复试算,只要迭代过程收敛、便可实现分岔点直接快速计算. 展开更多
关键词 准周期振动 对称破缺 达芬系统 分岔点 增量谐波平衡
在线阅读 下载PDF
上一页 1 2 3 下一页 到第
使用帮助 返回顶部