期刊文献+
共找到17篇文章
< 1 >
每页显示 20 50 100
Numerical Investigation on Dynamics of the Tendon System of A TLP by Applying Absolute Nodal Coordinate Formulation 被引量:2
1
作者 Luu Quang HUNG KANG Zhuang +1 位作者 ZHANG Cheng JIE Li-shao 《China Ocean Engineering》 SCIE EI CSCD 2021年第3期384-397,共14页
In the present study,the dynamics of the tendon system of a tension-leg platform(TLP)is investigated through the absolute nodal coordinate formulation(ANCF).Based on the energy conversion principle,the stiffness,gener... In the present study,the dynamics of the tendon system of a tension-leg platform(TLP)is investigated through the absolute nodal coordinate formulation(ANCF).Based on the energy conversion principle,the stiffness,generalized elastic force,external load and mass matrices of the element are deduced to perform the element assembling by using the finite element method.Then the motion equation of the tendon/riser is established.In this study,the TLP in the International Ship Structures Committee(ISSC)model under the first and second wave forces is considered as the case study.The simulation is performed in the MATLAB environment.Moreover,the accuracy and reliability of the programs are verified for cases of beam model with theoretical solutions.It is found that the motion response of tendons is affected by the TLP movement and environmental load,simultaneously.Then,the motion response is calculated using the SESAM software and exported as the boundary of ANCF tendons.Finally,the static and dynamic characteristics of the four tendons of ISSC TLP are analyzed systematically by the ANCF method.Performed analysis proves the effectiveness and feasibility of the ANCF method.It is concluded that the proposed method is a powerful scheme for calculating the dynamics of tendon/riser in the field of ocean engineering. 展开更多
关键词 absolute nodal coordinate formulation(ancf) tendon/riser TLP platform large deformation top and middle nodes wave direction
在线阅读 下载PDF
配置节点策略的绝对节点坐标平面梁单元 被引量:1
2
作者 王佳 张钰龙 王腾飞 《哈尔滨工程大学学报》 北大核心 2025年第7期1405-1411,共7页
针对传统绝对节点坐标梁单元计算效率和收敛率较低的问题,本文提出了基于配置节点策略的绝对节点坐标有限元多节点平面梁单元的构建方法。单元位移场采用非均布的P阶勒让德多项式二阶导数的零点和单元边界点作为节点进行离散插值,单元... 针对传统绝对节点坐标梁单元计算效率和收敛率较低的问题,本文提出了基于配置节点策略的绝对节点坐标有限元多节点平面梁单元的构建方法。单元位移场采用非均布的P阶勒让德多项式二阶导数的零点和单元边界点作为节点进行离散插值,单元的弹性能和动能采用P点的龙贝格数值积分算法进行计算,这种配置使单元节点与数值积分点重合。采用P-1阶拉格朗日插值函数建立单元的位移场,使用增强的连续介质力学方法推导单元的弹性力及刚度阵,基于配置节点策略,推导了多节点绝对节点坐标平面梁单元。通过经典算例分析,验证了该单元在静力学及动力学方面的表现,讨论了单元的收敛性、精确性及计算效率。 展开更多
关键词 绝对节点坐标法 配置节点策略 配置节点绝对节点坐标梁单元 多节点梁单元 增强连续介质力学 p-细化方法 勒让德多项式 拉格朗日插值
在线阅读 下载PDF
具有随机参数的空间结构热致振动分析
3
作者 严律豪 靳艳飞 孙加亮 《力学学报》 北大核心 2025年第5期1233-1245,共13页
在轨空间结构受到周期性太阳辐射会产生热变形和热振动,考虑到设计及装配制造过程中的误差影响,研究具有随机参数的空间结构动力学模型的动态响应具有重要理论意义和应用价值.针对热环境下考虑多个随机参数的空间结构动力学模型,提出了... 在轨空间结构受到周期性太阳辐射会产生热变形和热振动,考虑到设计及装配制造过程中的误差影响,研究具有随机参数的空间结构动力学模型的动态响应具有重要理论意义和应用价值.针对热环境下考虑多个随机参数的空间结构动力学模型,提出了一种基于混沌多项式展开(PCE)的随机柔性多体系统动力学计算方法,研究了结构热致振动的动态响应,分析了不同随机参数对热致振动的影响.首先,利用绝对节点坐标法(ANCF)对涉及热效应的空间简化模型进行建模,通过傅里叶温度单元法分析结构温度响应,推导出空间热-结构耦合模型.其次,利用PCE法构建原系统的代理模型,通过单项求容积法则(MCR)选取样本点,结合随机响应面法(SRSM)计算出PCE系数,得到热致振动响应的均值和标准差,并与蒙特卡洛模拟(MCS)结果进行对比,发现本文提出的算法在满足精度要求的情况下,所需样本点仅为MCS方法的1.5%,计算效率比MCS方法提高了76倍.最后,分析了多个随机参数包括弹性模量、梁的长度、太阳辐射吸收率、表面辐射系数、太阳辐射强度和太阳辐射角度对结构热致振动响应的影响,结果表明梁的长度和弹性模量的随机离散程度对热致振动响应的影响最为明显. 展开更多
关键词 热致振动 绝对节点坐标法 混沌多项式展开 随机参数 空间热-结构耦合模型
在线阅读 下载PDF
基于多项式混沌方法的柔性多体系统不确定性分析 被引量:11
4
作者 皮霆 张云清 吴景铼 《中国机械工程》 EI CAS CSCD 北大核心 2011年第19期2341-2343,2348,共4页
研究了柔性多体系统的参数不确定性问题。采用绝对节点坐标方法来描述柔性体的变形,运用多项式混沌方法对柔性多体系统进行不确定性分析。通过理论分析与算例研究,对比了多项式混沌方法与广泛使用的Monte Carlo方法的精度和效率。结果表... 研究了柔性多体系统的参数不确定性问题。采用绝对节点坐标方法来描述柔性体的变形,运用多项式混沌方法对柔性多体系统进行不确定性分析。通过理论分析与算例研究,对比了多项式混沌方法与广泛使用的Monte Carlo方法的精度和效率。结果表明,在相同精度要求下,多项式混沌方法在复杂柔性多体系统不确定性分析中效率更高,在工程分析中更具实用价值。 展开更多
关键词 柔性多体系统 绝对节点坐标 多项式混沌 随机响应面方法
在线阅读 下载PDF
基于绝对节点坐标法的柔性多体系统灵敏度分析 被引量:2
5
作者 王铁成 陈国平 孙东阳 《振动与冲击》 EI CSCD 北大核心 2015年第24期89-92,103,共5页
针对多个设计变量情况下的柔性多体系统灵敏度分析,采用伴随变量法对基于绝对节点坐标法建立的柔性多体系统进行了研究。为了验证的该方法的计算效率,分别采用直接微分法和伴随变量法对受重力作用的柔性单摆进行了研究,结果表明:这两种... 针对多个设计变量情况下的柔性多体系统灵敏度分析,采用伴随变量法对基于绝对节点坐标法建立的柔性多体系统进行了研究。为了验证的该方法的计算效率,分别采用直接微分法和伴随变量法对受重力作用的柔性单摆进行了研究,结果表明:这两种方法计算结果的误差很小,随着设计变量数量的增加,伴随变量法有更高的计算效率。 展开更多
关键词 多体系统 绝对节点坐标法 灵敏度 伴随变量法
在线阅读 下载PDF
一种外LET柔顺半铰的动力学建模与分析 被引量:1
6
作者 李鹏飞 曹博宇 +2 位作者 汪振宇 赵晨 赵荣昌 《中国机械工程》 EI CAS CSCD 北大核心 2019年第14期1727-1733,共7页
对一种外LET柔顺半铰进行建模与动力学特性分析.采用斜率不连续的绝对节点坐标法(ANCF)建立了该柔顺半铰系统动力学方程;对外力作用下半铰中心线和外轮廓的运动变形进行仿真,与实验结果对比,证明了绝对节点坐标法对柔顺半铰动力学建模... 对一种外LET柔顺半铰进行建模与动力学特性分析.采用斜率不连续的绝对节点坐标法(ANCF)建立了该柔顺半铰系统动力学方程;对外力作用下半铰中心线和外轮廓的运动变形进行仿真,与实验结果对比,证明了绝对节点坐标法对柔顺半铰动力学建模的有效性;进一步分析了半铰外形尺寸变化对其变形规律的影响. 展开更多
关键词 柔顺机构 外LET柔顺半铰 绝对节点坐标法 斜率不连续性
在线阅读 下载PDF
基于绝对节点坐标法的空间并联机构刚柔耦合动力学建模 被引量:2
7
作者 陈修龙 崔梦强 《农业机械学报》 EI CAS CSCD 北大核心 2021年第4期392-401,共10页
为研究柔性杆件对空间并联机构动力特性的影响,以具有结构冗余的3-RRPaR空间并联机构为研究对象,提出建立冗余空间并联机构刚柔耦合动力学模型的通用方法。基于绝对节点坐标法建立了三维二节点梁单元模型,分析了机构的结构特征,利用拉... 为研究柔性杆件对空间并联机构动力特性的影响,以具有结构冗余的3-RRPaR空间并联机构为研究对象,提出建立冗余空间并联机构刚柔耦合动力学模型的通用方法。基于绝对节点坐标法建立了三维二节点梁单元模型,分析了机构的结构特征,利用拉格朗日乘子法推导了机构的刚柔耦合动力学方程,采用广义α方法在Matlab中对其动力学方程进行数值求解,得到不同弹性模量下的动力学响应曲线。结果表明,柔性杆件对机构输出特性产生重要的影响,弹性模量越小,对机构的影响越大,且弹性模量对动平台加速度的影响最为显著。本研究为空间并联机构的刚柔耦合动力学建模提供了研究思路。 展开更多
关键词 空间并联机构 绝对节点坐标法 拉格朗日乘子法 刚柔耦合动力学
在线阅读 下载PDF
含间隙多连杆机构刚柔耦合动力学建模与分析 被引量:1
8
作者 陈修龙 江守源 《计算机集成制造系统》 EI CSCD 北大核心 2023年第10期3366-3379,共14页
针对间隙和柔性构件对多连杆机构动力学响应的影响,研究了一种适用于多连杆机械式压力机的平面六杆机构。建立了转动副间隙矢量模型,通过Lankarani-Nikravesh法向接触力模型和改进的库伦摩擦力模型建立了转动副间隙处的碰撞力模型。基... 针对间隙和柔性构件对多连杆机构动力学响应的影响,研究了一种适用于多连杆机械式压力机的平面六杆机构。建立了转动副间隙矢量模型,通过Lankarani-Nikravesh法向接触力模型和改进的库伦摩擦力模型建立了转动副间隙处的碰撞力模型。基于绝对节点坐标法建立了柔性梁单元模型,通过拉格朗日乘子法建立了含间隙机构的刚柔耦合动力学方程,并通过Adams虚拟仿真软件验证了动力学模型和理论分析的正确性。研究间隙值和摩擦系数对机构动力学响应的影响,利用相图、Poincare映射图和分岔图对机构的非线性特性进行分析,结果表明,随着间隙值的增大和摩擦系数的减小,机构受到的影响越大,产生的响应峰值越高,同时转动副间隙处的混沌现象有所增强。 展开更多
关键词 多连杆机构 绝对节点坐标法 间隙 刚柔耦合动力学 非线性特性
在线阅读 下载PDF
柔性多体系统设计灵敏度的高效计算方法
9
作者 皮霆 张云清 陈立平 《华中科技大学学报(自然科学版)》 EI CAS CSCD 北大核心 2011年第7期5-8,共4页
通过使用绝对节点坐标方法和非线性应变位移关系,提出了一种组装柔性多体系统设计灵敏度方程的方法.该方法不仅可以避免推导灵敏度方程所需的大量符号计算,而且能够提高计算灵敏度的效率.给出了平面柔性曲柄滑块机构的算例,通过直接微... 通过使用绝对节点坐标方法和非线性应变位移关系,提出了一种组装柔性多体系统设计灵敏度方程的方法.该方法不仅可以避免推导灵敏度方程所需的大量符号计算,而且能够提高计算灵敏度的效率.给出了平面柔性曲柄滑块机构的算例,通过直接微分法与伴随变量法的对比,表明所提出的方法不仅能够保证精度,而且对于中等规模的设计灵敏度分析,能够使直接微分法达到与伴随变量法相当的效率,具有广泛的工程应用价值. 展开更多
关键词 动力学 灵敏度分析 柔性多体系统 绝对节点坐标 直接微分法
原文传递
随机空间柔性多体系统动力学分析 被引量:8
10
作者 郭祥 靳艳飞 田强 《力学学报》 EI CSCD 北大核心 2020年第6期1730-1742,共13页
轻质、高精度的柔性多体系统被广泛应用于实际工程领域中.由于实际设计公差、制造误差及环境温度等多种不确定因素的存在,使得柔性多体系统的结构参数(物理参数和几何参数)表现出随机性.具有随机结构参数的动力学模型能够客观地反映出... 轻质、高精度的柔性多体系统被广泛应用于实际工程领域中.由于实际设计公差、制造误差及环境温度等多种不确定因素的存在,使得柔性多体系统的结构参数(物理参数和几何参数)表现出随机性.具有随机结构参数的动力学模型能够客观地反映出真实系统的动力学行为,且结构参数的不确定性对空间柔性多体系统动力学响应的影响是不容忽视的.针对具有多个随机参数的空间柔性多体系统,提出了一种基于广义alpha算法的非侵入式随机柔性多体系统动力学计算方法.采用绝对节点坐标公式(absolute node coordinate formulation,ANCF)来描述柔性体,推导建立多体系统动力学模型.利用混沌多项式展开(polynomial chaos expansion,PCE)法构建系统随机动力学方程的代理模型,然后将随机响应面法(stochastic response surface method,SRSM)嵌入广义-alpha方法中,分别采用改进抽样的回归方法(regression method of improved sampling,RMIS)和单项求容积法则(Monte Carlo simulation,MCR)来确定样本点.将数值计算结果与蒙特卡洛模拟(Monte Carlo simulation,MCS)结果进行对比,验证了所提算法的有效性.在相同的定积分精度的条件下,根据单项求容积法则确定的样本点的计算结果稳定性更强,且其计算效率更高. 展开更多
关键词 柔性多体系统 绝对节点坐标法 混沌多项式 单项式容积法 随机响应面法
在线阅读 下载PDF
基于绝对节点坐标法的平面梁有限变形下变形重构 被引量:6
11
作者 吴懋琦 谭述君 高飞雄 《力学学报》 EI CAS CSCD 北大核心 2021年第10期2776-2789,共14页
现有的对有限变形条件下柔性结构变形重构的研究往往单纯基于曲率与应变间的几何关系,同时忽略了被测体的纵向变形及其与弯曲变形的耦合效应.为得到一种更加精确且能借助现有的力学工具进行应用方向扩展的变形重构方法,以平面梁为对象,... 现有的对有限变形条件下柔性结构变形重构的研究往往单纯基于曲率与应变间的几何关系,同时忽略了被测体的纵向变形及其与弯曲变形的耦合效应.为得到一种更加精确且能借助现有的力学工具进行应用方向扩展的变形重构方法,以平面梁为对象,借鉴变形重构逆有限元法的思想,将平面梁的变形重构问题视作一类最优化问题.首先,通过引入绝对节点坐标法(absolute nodal coordinate formulation,ANCF)对柔性结构大变形下非线性的平面梁应变-位移关系进行精确描述,构造了一种逆梯度缩减ANCF平面索梁单元.然后,对此逆ANCF单元进行改进,在简化节点自由度的同时通过引入罚函数确保单元节点处的曲率连续性,既保证了本问题的适定性,也提升了最终解的精确性.最后,基于该单元利用Newton法构造了平面梁有限变形下变形重构问题的两种求解算法,即逐单元算法和多单元整体算法,以实现不同需求下的稳定求解.数值仿真结果表明,本方法在大变形条件下的变形重构误差小于1%,而且在测点较少的情况下依然保持较高的精度,同时验证了本方法的收敛性与计算效率. 展开更多
关键词 绝对节点坐标法 变形重构 逆有限元法 有限变形 平面梁
在线阅读 下载PDF
气动齿状软体驱动器动力学建模与仿真分析 被引量:2
12
作者 苏怡仪 张力文 +1 位作者 徐齐平 刘锦阳 《振动与冲击》 EI CSCD 北大核心 2023年第8期142-151,共10页
针对发生面内弯曲、超弹性材料气动齿状软体驱动器,从非线性本构关系出发,考虑几何非线性,基于绝对节点坐标法(absolute nodal coordinate formulation,ANCF)建立二维多梁结构的动力学模型。考虑多点接触,采用罚函数法对相邻气腔之间的... 针对发生面内弯曲、超弹性材料气动齿状软体驱动器,从非线性本构关系出发,考虑几何非线性,基于绝对节点坐标法(absolute nodal coordinate formulation,ANCF)建立二维多梁结构的动力学模型。考虑多点接触,采用罚函数法对相邻气腔之间的软接触问题建模,避免相邻气腔之间的互相穿透,并能准确模拟相邻气腔从不接触到充分接触的变拓扑过程;相比于ANCF实体单元模型,二维梁模型计算规模小,更加高效。通过气动弯曲试验和ABAQUS仿真,验证理论模型的计算精度;在此基础上研究重力对软体驱动器构型的影响,分析充气速率和最终加载气压值对其动力学响应的影响。 展开更多
关键词 绝对节点坐标法(ancf) 气动齿状软体驱动器 动力学建模 软接触 材料非线性 试验验证
在线阅读 下载PDF
Dynamic analysis and regulation of the flexible pipe conveying fluid with a hard-magnetic soft segment 被引量:1
13
作者 Zilong GUO Qiao NI +2 位作者 Wei CHEN Huliang DAI Lin WANG 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI CSCD 2022年第9期1415-1430,共16页
The recently developed hard-magnetic soft(HMS)materials can play a significant role in the actuation and control of medical devices,soft robots,flexible electronics,etc.To regulate the mechanical behaviors of the cant... The recently developed hard-magnetic soft(HMS)materials can play a significant role in the actuation and control of medical devices,soft robots,flexible electronics,etc.To regulate the mechanical behaviors of the cantilevered pipe conveying fluid,the present work introduces a segment made of the HMS material located somewhere along the pipe length.Based on the absolute node coordinate formulation(ANCF),the governing equations of the pipe conveying fluid with an HMS segment are derived by the generalized Lagrange equation.By solving the derived equations with numerical methods,the static deformation,linear vibration characteristic,and nonlinear dynamic response of the pipe are analyzed.The result of the static deformation of the pipe shows that when the HMS segment is located in the middle of the pipe,the downstream portion of the pipe centerline will keep a straight shape,providing that the pipe is stable with a relatively low flow velocity.Therefore,it is possible to precisely regulate the ejection direction of the fluid flow by changing the magnetic and fluid parameters.It is also found that the intensity and direction of the external magnetic field greatly affect the stability and dynamic response of the pipe with an HMS segment.In most cases,the magnetic actuation increases the critical flow velocity for the flutter instability of the pipe system and suppresses the vibration amplitude of the pipe. 展开更多
关键词 hard-magnetic soft(HMS)material pipe conveying fluid absolute node coordinate formulation(ancf) stability dynamic response REGULATION
在线阅读 下载PDF
基于绝对节点坐标法的高次插值梁单元建模 被引量:1
14
作者 李全乐 赵春花 《上海工程技术大学学报》 CAS 2018年第4期329-333,共5页
针对传统绝对节点坐标平面梁单元存在泊松闭锁的问题,研究平面高次插值梁单元.在多项式位移模式中,基于帕斯卡三角形理论,采用与x2项共用系数的方法引入y2项,建立含有y2项的绝对节点坐标高次插值梁单元模型.实例计算结果与传统绝对节点... 针对传统绝对节点坐标平面梁单元存在泊松闭锁的问题,研究平面高次插值梁单元.在多项式位移模式中,基于帕斯卡三角形理论,采用与x2项共用系数的方法引入y2项,建立含有y2项的绝对节点坐标高次插值梁单元模型.实例计算结果与传统绝对节点坐标平面梁单元的计算结果比较说明,在泊松比不为零的情况下,绝对节点坐标高次插值梁单元的计算结果明显好于传统绝对节点坐标平面梁单元.因此,含y2项的绝对节点坐标平面梁单元可有效缓解传统绝对节点坐标平面梁单元的泊松闭锁问题. 展开更多
关键词 绝对节点坐标法 高次插值梁单元 泊松闭锁
在线阅读 下载PDF
基于绝对节点坐标的废塑料薄膜振动特性研究
15
作者 王婷 杨先海 李永波 《山东理工大学学报(自然科学版)》 CAS 2022年第4期1-6,共6页
将废塑料薄膜进行分选回收是目前高效节能的垃圾处理方式,废旧塑料薄膜的振动特性会直接影响分选的精度。基于绝对节点坐标法对废塑料薄膜进行了振动特性分析,基于薄板单元模型理论提出了一种三维柔性薄膜单元模型,引入Neo-Hookean本构... 将废塑料薄膜进行分选回收是目前高效节能的垃圾处理方式,废旧塑料薄膜的振动特性会直接影响分选的精度。基于绝对节点坐标法对废塑料薄膜进行了振动特性分析,基于薄板单元模型理论提出了一种三维柔性薄膜单元模型,引入Neo-Hookean本构模型求解单元的质量矩阵和刚度矩阵,建立废塑料薄膜动力学方程表达式。利用MATLAB软件对薄膜进行自由振动模态分析,并将数值仿真结果与ANSYS结果进行比较,验证了薄膜单元的正确性,为后续废塑料薄膜振动分选的研究与应用提供了理论基础。 展开更多
关键词 绝对节点坐标法 薄膜振动特性 模态分析 塑料分选
在线阅读 下载PDF
绝对节点坐标梁单元弹性力建模及模态仿真的常见方法对比研究
16
作者 周川 赵春花 +1 位作者 独亚平 郭嘉辉 《上海工程技术大学学报》 CAS 2022年第2期196-204,共9页
通过理论推导和数值分析,研究二维绝对节点坐标梁单元在不同弹性力建模方法下的模态参数.介绍连续介质力学、增强连续介质力学以及应变分解法,并推导应变分解法应用的局限性,从理论角度说明了3种方法的特点.基于广义特征方程得到不同弹... 通过理论推导和数值分析,研究二维绝对节点坐标梁单元在不同弹性力建模方法下的模态参数.介绍连续介质力学、增强连续介质力学以及应变分解法,并推导应变分解法应用的局限性,从理论角度说明了3种方法的特点.基于广义特征方程得到不同弹性力建模方法下的绝对节点坐标梁单元的4种固有模态,以简支梁结构为例,分析绝对节点坐标梁单元弹性力建模方法对4种固有模态的作用规律.不同弹性力建模方法下,横向低阶单元的4种固有频率相较横向高阶单元更高,表现得更“刚”.连续介质力学方法下,各单元的低阶剪切固有频率与解析解存在20%~30%误差,而应变分解法和增强连续介质力学方法能够将剪切固有频率的误差控制在4%内,提高单元收敛精度. 展开更多
关键词 绝对节点坐标法 固有频率计算 应变分解法 增强连续介质力学方法 连续介质力学方法
在线阅读 下载PDF
空间充气展开绳网捕获系统动力学建模与分析 被引量:10
17
作者 刘昊 魏承 +2 位作者 田健 谭春林 赵阳 《机械工程学报》 EI CAS CSCD 北大核心 2018年第22期145-152,共8页
空间充气展开绳网系统在空间捕获中具有自身特殊的优势:占用体积小,展开后捕获范围大,且绳网由充气杆支撑,结构更加稳固,比空间飞网系统更加具有应用潜力。基于绝对节点坐标法(Absolute node coordinate method,ANCF)建立空间充气展开... 空间充气展开绳网系统在空间捕获中具有自身特殊的优势:占用体积小,展开后捕获范围大,且绳网由充气杆支撑,结构更加稳固,比空间飞网系统更加具有应用潜力。基于绝对节点坐标法(Absolute node coordinate method,ANCF)建立空间充气展开绳网捕获系统动力学模型,应用Hertz理论建立柔性体碰撞动力学模型,对系统进行动力学仿真分析。仿真表明,静力学和振动分析能够为机构力学参数选型提供参考,相比充气杆末端激励,绳网中部受到激励时对捕获机构振动影响更大;机构收口拉力随着充气杆刚度的增加而增大,对充气杆的放气处理可以降低其刚度,从而减小收口拉力,约为原来的10%,易于收口机构的设计。通过对期望目标的捕获仿真,分析了影响收口过程的关键指标,证明了捕获策略的合理性和动力学模型的准确性。研究结果能够对空间充气展开绳网捕获系统设计提供有效的依据。 展开更多
关键词 充气展开杆 空间绳网 绝对节点坐标法 振动分析 收口捕获 碰撞计算
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部