期刊文献+
共找到54篇文章
< 1 2 3 >
每页显示 20 50 100
一族求多项式零点的高阶并行迭代法
1
作者 张昕 钟孝畴 《湖南理工学院学报(自然科学版)》 2025年第1期7-9,共3页
讨论一族同时求多项式零点的并行迭代法,证明其收敛性且收敛阶为4,并通过数值例子验证该迭代法具有较高的计算效率.
关键词 多项式零点 迭代法 收敛阶
在线阅读 下载PDF
骨铣削参数对机器人颤振稳定性影响与优化
2
作者 田和强 安金昌 +1 位作者 马洪强 孟德宝 《哈尔滨工程大学学报》 北大核心 2025年第5期988-996,共9页
为解决骨科手术机器人铣削中的颤振问题,本文提出了一种优化方法,以提升手术的安全性和质量。根据再生型颤振理论,建立了动态铣削力模型,并采用零阶频域法构建颤振稳定域模型,同时绘制了叶瓣图,以分析铣削参数的影响。以骨骼去除率为优... 为解决骨科手术机器人铣削中的颤振问题,本文提出了一种优化方法,以提升手术的安全性和质量。根据再生型颤振理论,建立了动态铣削力模型,并采用零阶频域法构建颤振稳定域模型,同时绘制了叶瓣图,以分析铣削参数的影响。以骨骼去除率为优化目标,通过实验验证确定的最佳铣削参数,显著减少了颤振,并提高了手术的稳定性和效率。这项研究为骨科手术机器人在临床应用中的安全性提供了理论依据,并为铣削工艺的优化提供了实践指导。 展开更多
关键词 骨科手术机器人 再生型颤振 动态铣削力模型 零阶频域法 颤振稳定性 叶瓣图 铣削参数 优化
在线阅读 下载PDF
基于轴向准零刚度的机器人末端隔振技术研究
3
作者 王杨杰 俞哲 +3 位作者 俞建峰 钱陈豪 蒋毅 化春键 《噪声与振动控制》 北大核心 2025年第4期285-291,共7页
针对工业机器人末端作业时轴向低频隔振难的问题,设计一种由等边六边形拉簧装置和正刚度弹簧组构成的轴向准零刚度(Quasi Zero Stiffness,QZS)隔振器。首先,通过静力学分析初始高度对等边六边形拉簧装置负刚度性能和隔振器刚度特性的影... 针对工业机器人末端作业时轴向低频隔振难的问题,设计一种由等边六边形拉簧装置和正刚度弹簧组构成的轴向准零刚度(Quasi Zero Stiffness,QZS)隔振器。首先,通过静力学分析初始高度对等边六边形拉簧装置负刚度性能和隔振器刚度特性的影响;其次,建立非线性动力学方程,应用谐波平衡法推导力传递率,通过数值仿真分析振动幅值、阻尼比、初始高度和负载对力传递率的影响。最后,在不同外部激励下对轴向QZS隔振器和线性隔振器进行实验研究,并将实验结果与四阶龙格库塔法求解的加速度响应进行对比。结果表明:轴向QZS隔振器在载荷为20~80 N时能隔离3.2~20 Hz的振动;相较于其对应的线性隔振器,该隔振器的固有频率降低到1.8 Hz,起始隔振频率降低76.6%,隔振频带增宽166.7%。 展开更多
关键词 振动与波 工业机器人 低频隔振 准零刚度 力传递率 4阶龙格库塔法
在线阅读 下载PDF
零阶与一阶优化算法在悬索桥模型修正中的应用对比分析 被引量:26
4
作者 郭彤 李爱群 +1 位作者 费庆国 王浩 《振动与冲击》 EI CSCD 北大核心 2007年第4期35-38,68,共5页
为了给健康监测研究提供一个较为准确的有限元基准模型,采用零阶与一阶优化算法作为悬索桥模型修正的基本工具,同时将优化算法与大型有限元分析软件相结合,以考虑大跨悬索桥的非线性和应力刚度。对零阶与一阶优化算法在悬索桥模型修正... 为了给健康监测研究提供一个较为准确的有限元基准模型,采用零阶与一阶优化算法作为悬索桥模型修正的基本工具,同时将优化算法与大型有限元分析软件相结合,以考虑大跨悬索桥的非线性和应力刚度。对零阶与一阶优化算法在悬索桥模型修正中的具体实现方法进行了探讨。针对润扬大桥悬索桥,利用健康监测系统的实测数据进行了优化算例分析,优化过程中结合了动力与静力两方面的约束信息,使得结果更趋于合理。对两种算法的准确性和计算效率进行了对比分析,提出了应用建议。 展开更多
关键词 零阶优化算法 一阶优化算法 模型修正 悬索桥 非线性
在线阅读 下载PDF
钢管拱肋分段吊装扣索一次张拉索力改进算法 被引量:33
5
作者 周倩 周建庭 +2 位作者 马虎 黎小刚 张兰 《交通运输工程学报》 EI CSCD 北大核心 2020年第1期92-101,共10页
为改善大跨钢管拱肋分段吊装扣索索力常用算法迭代效率低、计算时耗长,且忽略了温变影响等不足,建立了可考虑温变影响和提高计算效率的改进算法;基于材料力学和几何学相关知识,推导了吊装过程中拱肋位移变化与温变的理论关系,并在计入... 为改善大跨钢管拱肋分段吊装扣索索力常用算法迭代效率低、计算时耗长,且忽略了温变影响等不足,建立了可考虑温变影响和提高计算效率的改进算法;基于材料力学和几何学相关知识,推导了吊装过程中拱肋位移变化与温变的理论关系,并在计入温变引起索长和拱肋位移改变的情况下,推导出扣索索力变化与温变的理论关系;基于扣索一次张拉法和ANSYS零阶优化法,开发了考虑温变影响且在迭代子步中对程序自动搜索实施宏观调控的扣索索力计算程序;运用改进算法对某主跨300 m钢管混凝土拱桥开展了分段吊装施工控制分析。分析结果表明:推导的理论公式和有限元分析结果的变化规律一致,拱肋位移变化的最大相对误差为11%,索力变化的最大相对误差为18%,均能满足工程精度要求;与原算法相比,采用改进算法的迭代次数由26次缩减到17次,迭代效率提高了35%,计算索力与实测索力的最大偏差由276 kN减小到100 kN;拱肋松索成拱位移理论值与实测值的最大偏差为7 mm,成拱线形正常;建立的改进算法可实现扣索一次张拉,提高迭代效率和计算精度,运用改进算法控制大跨钢管拱肋吊装施工可使拱肋松索成拱线形满足设计要求。 展开更多
关键词 桥梁工程 拱桥 分段吊装 零阶优化法 扣索一次张拉法 索力 线形 温变
原文传递
基于粒子群算法和ANSYS的层合板自由边界处层间正应力的优化研究 被引量:4
6
作者 彭文杰 陈建桥 +1 位作者 顾明凯 涂文琼 《机械科学与技术》 CSCD 北大核心 2009年第11期1496-1500,1504,共6页
由于层合板结构相邻铺层之间材料参数的不匹配,层合板受外载作用时,板的自由边界处会产生奇异的三维层间应力场。笔者以层合板自由边界处层间正应力为目标,对层合板的铺层结构进行了优化。优化方法采用粒子群算法(PSO)与有限元相结合的... 由于层合板结构相邻铺层之间材料参数的不匹配,层合板受外载作用时,板的自由边界处会产生奇异的三维层间应力场。笔者以层合板自由边界处层间正应力为目标,对层合板的铺层结构进行了优化。优化方法采用粒子群算法(PSO)与有限元相结合的方法,开发了PSO和通用有限元软件之间的接口程序。笔者还采用ANSYS自带的零阶优化方法进行了优化,并与PSO方法做了比较。优化算例分别考虑了面内拉伸载荷和面外横向载荷。优化结果表明:零阶方法的优化结果对初始值敏感,初始值不同,得到的优化解也不相同。而PSO方法由于是利用随机生成的初始值并通过进化来搜寻全局最优解,从而避免了初始值的影响,且PSO在寻优效率上具有显著的优势。 展开更多
关键词 层间应力 自由边界 粒子群算法 零阶方法 铺层设计
在线阅读 下载PDF
考虑结合处弯矩限值的混合梁斜拉桥索力优化 被引量:9
7
作者 杨德灿 章芳芳 《华中科技大学学报(自然科学版)》 EI CAS CSCD 北大核心 2010年第8期111-114,共4页
根据混合梁斜拉桥的结构特点,考虑混凝土主梁和钢主梁结合处弯矩的限值,利用有限元方法建立混合梁斜拉桥成桥索力的优化模型.以塔梁的弯曲和拉压应变能之和为目标函数,以斜拉索的初张力为设计变量,成桥索力和主梁钢混结合处的弯矩为约... 根据混合梁斜拉桥的结构特点,考虑混凝土主梁和钢主梁结合处弯矩的限值,利用有限元方法建立混合梁斜拉桥成桥索力的优化模型.以塔梁的弯曲和拉压应变能之和为目标函数,以斜拉索的初张力为设计变量,成桥索力和主梁钢混结合处的弯矩为约束条件,结合ANSYS优化模块中的零阶和一阶方法进行求解,确定合理成桥状态下的索力.在优化过程中,改变结合处弯矩的限制范围,分析成桥下的内力和位移随限制范围变化的规律.通过一座单塔斜拉桥算例计算表明:将弯矩限值控制在0.3~0.5倍的主梁最小和最大弯矩范围内,得到的成桥状态最合理.索力从塔至边墩呈均匀增长趋势,主梁弯矩图变化平顺,结合处弯矩较小. 展开更多
关键词 混合梁斜拉桥 弯矩 约束条件 零阶方法 一阶方法 合理成桥索力
原文传递
2500t环轨起重机主臂架的优化设计 被引量:7
8
作者 周奇才 张许辉 +2 位作者 赵炯 罗明川 熊肖磊 《中国工程机械学报》 北大核心 2015年第1期49-53,共5页
采用分层优化的思想,以总重量为目标函数,以臂架强度、整体刚度、局部刚度和臂端挠度为状态变量,分别用零阶方法和一阶方法对2 500t环轨起重机主臂架进行了优化设计.主臂架重量减轻了18.5%,效果显著.对比分析了两种优化方法的优缺点及... 采用分层优化的思想,以总重量为目标函数,以臂架强度、整体刚度、局部刚度和臂端挠度为状态变量,分别用零阶方法和一阶方法对2 500t环轨起重机主臂架进行了优化设计.主臂架重量减轻了18.5%,效果显著.对比分析了两种优化方法的优缺点及适用范围,为大型桁架式臂架优化方法的选择提供了参考. 展开更多
关键词 分层优化 零阶方法 一阶方法 优化设计
在线阅读 下载PDF
钢管混凝土拱桥扣索一次张拉索力优化改进算法 被引量:12
9
作者 周倩 周建庭 +1 位作者 陈静雯 冯雨实 《福州大学学报(自然科学版)》 CAS 北大核心 2019年第3期412-416,共5页
为避免施工中反复调索的繁琐工作,同时保证钢拱肋松索成拱后线形满足期望要求,基于零阶优化法,运用大型有限元分析软件ANSYS,开展钢管拱肋扣索一次张拉索力优化研究.提出基于扣索内力法确定索力迭代初值并在各迭代子步中对程序自动搜索... 为避免施工中反复调索的繁琐工作,同时保证钢拱肋松索成拱后线形满足期望要求,基于零阶优化法,运用大型有限元分析软件ANSYS,开展钢管拱肋扣索一次张拉索力优化研究.提出基于扣索内力法确定索力迭代初值并在各迭代子步中对程序自动搜索值进行宏观调整的改进算法.算例表明:改进算法可提高迭代效率和计算精度,按优化出的索力和预抬量进行施工控制,钢拱肋松索成拱线形与设计值吻合良好,克服了传统力学法调索中存在的较大马鞍形变形问题. 展开更多
关键词 钢管混凝土拱桥 零阶优化法 切线初始位移法 扣索索力 预抬量
在线阅读 下载PDF
浮式风力机平台摇荡运动混沌分析及预测 被引量:3
10
作者 余万 李春 +3 位作者 丁勤卫 郝文星 周红杰 袁全勇 《工程热物理学报》 EI CAS CSCD 北大核心 2019年第5期1050-1057,共8页
为探究漂浮式风力机平台在多种非规则环境载荷作用下摇荡运动中的非线性动力学特征,建立基于OC3 Spar平台的NREL 5 MW风力机整机模型,采用水动力软件AQWA求解平台的摇荡运动,基于相空间重构及最大李亚普诺夫指数法从定性及定量角度判断... 为探究漂浮式风力机平台在多种非规则环境载荷作用下摇荡运动中的非线性动力学特征,建立基于OC3 Spar平台的NREL 5 MW风力机整机模型,采用水动力软件AQWA求解平台的摇荡运动,基于相空间重构及最大李亚普诺夫指数法从定性及定量角度判断平台在风、浪、流载荷作用下摇荡运动时间序列的混沌特征,且进一步采用加权零阶局域法预测平台运动.结果表明:在低频时Spar平台幅值运动响应算子较大;在多重载荷作用时平台在平衡位置处做摇荡运动;平台在纵荡、垂荡及纵摇方向上的运动时间序列具有混沌特征,且三个方向短期可预测时间最大为0.4879 s;预测数据时间超过短期可预测时间时,随着时间增大预测数据和实际数据的绝对误差增大.混沌理论是非线性动力学的重要研究方向,可作为漂浮式风力机平台非线性运动分析的一种新方法. 展开更多
关键词 漂浮式风力机 SPAR平台 相空间重构 最大李亚普诺夫指数 加权零阶局域法
原文传递
基于零阶优化方法的无人机复合材料进气道结构优化与分析 被引量:5
11
作者 寇剑锋 郑锡涛 彭兴国 《机械强度》 EI CAS CSCD 北大核心 2008年第6期1018-1022,共5页
在综述优化设计理论研究的基础上,运用ANSYS分析软件,采用零阶优化方法,对某型无人机进气道复合材料层合板结构进行初始铺层的强度和刚度分析,在此基础上对其铺层厚度和铺层比例进行优化。采用两种优化方案,一种是在载荷一定的情况下求... 在综述优化设计理论研究的基础上,运用ANSYS分析软件,采用零阶优化方法,对某型无人机进气道复合材料层合板结构进行初始铺层的强度和刚度分析,在此基础上对其铺层厚度和铺层比例进行优化。采用两种优化方案,一种是在载荷一定的情况下求得尽可能轻的重量;另一种是在厚度一定的情况下调整铺层比例,使结构的承载能力更强。在得到ANSYS优化结果后再对其进行修正。优化结果显示,通过优化能进一步减轻结构重量,增强结构的承载能力。 展开更多
关键词 复合材料 零阶优化 铺层修正 无人机 进气道
在线阅读 下载PDF
零阶方法在特高压变压器线圈电场优化中的应用 被引量:2
12
作者 王立伟 吕殿利 +1 位作者 韩婷彦 汪友华 《河北工业大学学报》 CAS 北大核心 2011年第4期5-8,共4页
采用有限元分析软件ANSYS的参数化编程语言APDL(ANSYS Parametric Design Language)建立了特高压变压器底部线圈的参数化有限元模型.计算了底部线圈的电场分布,并提取出最大电场强度值.在此基础上,针对线圈结构的参数设计,采用零阶优化... 采用有限元分析软件ANSYS的参数化编程语言APDL(ANSYS Parametric Design Language)建立了特高压变压器底部线圈的参数化有限元模型.计算了底部线圈的电场分布,并提取出最大电场强度值.在此基础上,针对线圈结构的参数设计,采用零阶优化方法,对特高压变压器底部线圈的电场进行优化.优化结果表明,采用零阶优化方法优化后,最大电场强度强度值明显降低. 展开更多
关键词 特高压变压器 零阶优化方法 电场分析 三维计算 线圈
在线阅读 下载PDF
吊具夹持机构的结构优化研究 被引量:4
13
作者 代春香 孟令宽 李三雁 《中国测试》 CAS 北大核心 2019年第10期164-168,共5页
为使吊具夹持机构满足载重工况的同时又尽可能降低自身质量及体积,对其进行接触应力、结构强度及优化设计分析,并根据最优设计参数作为制作夹持机构实物的具体数据,以实现对大直径、高质量、高温度零件的吊装。首先,该文采用Ansys Workb... 为使吊具夹持机构满足载重工况的同时又尽可能降低自身质量及体积,对其进行接触应力、结构强度及优化设计分析,并根据最优设计参数作为制作夹持机构实物的具体数据,以实现对大直径、高质量、高温度零件的吊装。首先,该文采用Ansys Workbench建立夹持机构的参数化有限元模型,并进行接触分析,得出强度相对薄弱位置,进而提出改进和优化方案。其次,建立夹持机构的结构优化模型,采用零阶优化方法进行优化处理,从而得到模型的最优设计参数。从优化结果可得,夹持机构的比重较轻,夹持臂厚度变薄,而此时的最大综合位移和最大等效应力依然满足设计要求。最后,以优化分析结果作为夹持机构的实物制作数据支撑,制作调试后的吊具现场使用状况良好,效果比较理想,可验证分析结果的准确性。 展开更多
关键词 夹持机构 有限元模型 接触分析 零阶优化方法 结构优化
在线阅读 下载PDF
D24施工便梁有限元分析及优化设计 被引量:4
14
作者 王玉璞 程文明 《机械设计与制造》 北大核心 2018年第4期5-8,共4页
施工便梁是铁路线路维修改造所需的重要设备。但目前对于该设备结构的研究却相对较为薄弱。而随着我国高铁事业的逐步推进,施工便梁必将发挥愈加重要的作用。运用有限元方法,对铁路施工便梁金属结构进行静力学分析。并运用ANSYS的零阶... 施工便梁是铁路线路维修改造所需的重要设备。但目前对于该设备结构的研究却相对较为薄弱。而随着我国高铁事业的逐步推进,施工便梁必将发挥愈加重要的作用。运用有限元方法,对铁路施工便梁金属结构进行静力学分析。并运用ANSYS的零阶优化方法,采用SUMT算法对便梁进行结构的优化设计。在保证结构满足应力许用值的前提下,提高材料使用效率。经过优化后,便梁结构自重减少12.6%。并通过ANSYS分析验证了优化后的方案可行性。研究对于施工便梁及相关同类产品的结构优化有较高的借鉴意义。 展开更多
关键词 施工便梁 有限元法 零阶优化 SUMT 优化设计
在线阅读 下载PDF
基于非线性动力学的滑坡位移预测 被引量:10
15
作者 陈卫兵 李端有 《长江科学院院报》 CSCD 北大核心 2006年第2期28-30,34,共4页
利用非线性动力学的宏观研究方法对清江库岸茅坪滑坡的位移时间序列进行了分析,发现滑坡系统在演变过程中存在混沌。在此基础上,介绍了一种混沌的平均零阶局域预测法,对滑坡位移进行了预测计算,第1个月的相对误差仅为0.66%,并且前3个月... 利用非线性动力学的宏观研究方法对清江库岸茅坪滑坡的位移时间序列进行了分析,发现滑坡系统在演变过程中存在混沌。在此基础上,介绍了一种混沌的平均零阶局域预测法,对滑坡位移进行了预测计算,第1个月的相对误差仅为0.66%,并且前3个月的相对误差都控制在5%以内,初步计算结果令人满意。 展开更多
关键词 滑坡 相空间重构 最大LYAPUNOV指数 平均零阶局域预测法
在线阅读 下载PDF
湿室型泵房整体结构优化设计 被引量:1
16
作者 蔡新 崔朕铭 +2 位作者 陈卫东 黄海田 郭兴文 《河海大学学报(自然科学版)》 CAS CSCD 北大核心 2018年第1期37-42,共6页
针对现行规范方法设计的泵房存在"粗、大、笨"的问题,从降低泵房建设成本的角度,构建以泵房结构总造价最低为目标函数,泵房关键结构几何尺寸等参数为设计变量,泵房结构基底应力、地基承载力、抗滑稳定性、结构强度、沉降量及... 针对现行规范方法设计的泵房存在"粗、大、笨"的问题,从降低泵房建设成本的角度,构建以泵房结构总造价最低为目标函数,泵房关键结构几何尺寸等参数为设计变量,泵房结构基底应力、地基承载力、抗滑稳定性、结构强度、沉降量及沉降差为约束条件的优化设计数学模型。采用有限元法对泵房结构进行三维有限元数值模拟,同时借助ANSYS软件提供的优化平台对其进行优化设计研究。结果表明:与原设计方案相比,优化方案的泵房总造价减少了9.7%,结构的缝墩厚度达到限定值,既减少了工程总造价,也降低了结构对地基的要求。 展开更多
关键词 湿室型泵房 整体结构 优化设计 有限元法 参数化模型 重分析 零阶法 一阶法
在线阅读 下载PDF
石杉碱甲渗透泵控释微丸的制备及其药动学研究 被引量:6
17
作者 尹丽娜 张雅雯 +3 位作者 陈晓晓 周云琍 陈燕 郑高利 《中国药学杂志》 CAS CSCD 北大核心 2019年第9期720-725,共6页
目的制备石杉碱甲渗透泵控释微丸,并探讨其在Beagle犬内的药动学特性。方法采用挤出滚圆法制备渗透泵含药丸芯,以乙基纤维素和PEG400用量为考察因素结合星点设计效应面法优选包衣液处方,通过流化床包衣工艺制备石杉碱甲渗透泵控释微丸... 目的制备石杉碱甲渗透泵控释微丸,并探讨其在Beagle犬内的药动学特性。方法采用挤出滚圆法制备渗透泵含药丸芯,以乙基纤维素和PEG400用量为考察因素结合星点设计效应面法优选包衣液处方,通过流化床包衣工艺制备石杉碱甲渗透泵控释微丸。按最优处方制备3批微丸并以累积释放度与时间进行零级、一级、Higuchi方程拟合考察其释药特性。以市售的石杉碱甲普通片为参比制剂,考察自制石杉碱甲渗透泵控释微丸的体内药动学特点,并比较两制剂的生物等效性。结果微丸最优处方为PEG400用量10. 5%,乙基纤维素用量61. 5%,所制备的石杉碱甲渗透泵控释微丸24 h内释药具有零级释放特征,释放动力以渗透压为主。自制石杉碱甲渗透泵控释微丸的tmax与ρmax明显低于参比制剂,其t1/2比参比制剂延长约1倍,相对生物利用度为95. 8%。结论 Hup-A渗透泵控释微丸在Beagle犬中有较好的缓释效果,体内相关性良好。 展开更多
关键词 石杉碱甲 渗透泵控释微丸 星点设计效应面 零级释放 渗透压 累积释放度 药动学
原文传递
钢管混凝土拱桥拱肋吊装线形控制的分步算法 被引量:6
18
作者 赵洋 周水兴 刘静 《重庆交通大学学报(自然科学版)》 CAS 北大核心 2010年第1期1-3,37,共4页
针对钢管混凝土拱桥拱肋分节段安装由扣索张拉引起的高程差问题,以裸拱自重变形后的拱轴线形为控制目标,基于最优化理论和一次扣索张拉法,提出先用零阶优化法按整体安装计算出各拱肋节段的预抬量和扣索索力,再以整体安装计算结果为目标... 针对钢管混凝土拱桥拱肋分节段安装由扣索张拉引起的高程差问题,以裸拱自重变形后的拱轴线形为控制目标,基于最优化理论和一次扣索张拉法,提出先用零阶优化法按整体安装计算出各拱肋节段的预抬量和扣索索力,再以整体安装计算结果为目标,根据拱肋节段安装顺序,通过迭代方法计算出各个节段安装时的预抬量和扣索索力。将该算法应用到主跨240m的巫山新龙门大桥拱肋安装线形控制中,松索成拱后的线形与目标线形吻合良好。 展开更多
关键词 钢管混凝土拱桥 分步算法 线形控制 零阶优化法 整体安装
在线阅读 下载PDF
一种同时求多项式零点的加速迭代法 被引量:3
19
作者 刘兰冬 蒙杨 《工程数学学报》 CSCD 北大核心 2007年第5期935-938,共4页
本文讨论了在无重根情况下,利用改进的Newton迭代法对一种同时求多项式零点的并行迭代法进行加速,得到了一种新的加速迭代法。首先证明了该方法是收敛的,并且理论证明出收敛阶至少是5阶;其次,分析了该方法的计算效率;最后通过实际的数... 本文讨论了在无重根情况下,利用改进的Newton迭代法对一种同时求多项式零点的并行迭代法进行加速,得到了一种新的加速迭代法。首先证明了该方法是收敛的,并且理论证明出收敛阶至少是5阶;其次,分析了该方法的计算效率;最后通过实际的数值算例表明:计算收敛阶和定理结论是一致的,且本算法具有较高的计算效率。 展开更多
关键词 多项式零点 NEWTON方法 收敛阶
在线阅读 下载PDF
上一页 1 2 3 下一页 到第
使用帮助 返回顶部