期刊文献+
共找到27篇文章
< 1 2 >
每页显示 20 50 100
流固耦合问题中基于单元变形的有限元网格更新方法 被引量:1
1
作者 孙芳锦 张大明 殷志祥 《辽宁工程技术大学学报(自然科学版)》 CAS 北大核心 2008年第5期674-676,共3页
基于伪实体模型,根据引起单元变形的原因提出了一种网格更新方法。该方法将伪实体单元刚化矩阵中的刚化因数选为单元最大、最小主应变的指数函数,同时提出区分拉力和压力给单元变形带来的不同影响,引入修正因子修正刚化因数,最后得到刚... 基于伪实体模型,根据引起单元变形的原因提出了一种网格更新方法。该方法将伪实体单元刚化矩阵中的刚化因数选为单元最大、最小主应变的指数函数,同时提出区分拉力和压力给单元变形带来的不同影响,引入修正因子修正刚化因数,最后得到刚化因数的函数表达式。将该方法应用于典型的流固耦合问题,结果表明:本方法的正确性和优越性。提出的网格更新方法在防止网格单元极度变形,维持网格整体质量上是非常有效的。 展开更多
关键词 流固耦合 伪实体模型 单元变形 网格更新 网格质量
下载PDF
三维可变形圆化多面体离散单元法
2
作者 毛佳 余健坤 +1 位作者 邵琳玉 赵兰浩 《岩土力学》 EI CAS CSCD 北大核心 2024年第3期908-916,共9页
为真实模拟岩体的变形特性与运动形态,将圆化多面体离散元法与有限单元法结合,提出一种三维可变形圆化多面体离散单元法。该方法既能真实表征块体的不规则特征,又降低了接触判断的难度,同时能准确反映块体的变形特性。在求解切向接触力... 为真实模拟岩体的变形特性与运动形态,将圆化多面体离散元法与有限单元法结合,提出一种三维可变形圆化多面体离散单元法。该方法既能真实表征块体的不规则特征,又降低了接触判断的难度,同时能准确反映块体的变形特性。在求解切向接触力时,将接触判断对象从接触对简化为单元整体,显著提高了计算效率。为分析块体变形特性,在块体离散单元内部划分有限元网格,将最外层有限元网格作为最小接触单元。采用直接平均法,将接触力转化为等效节点力,并采用非线性有限单元法实现对单元变形特性的精确模拟,克服了圆化多面体无法反映单元变形的缺陷。通过5个算例论证了新方法在捕捉单元变形、运动形态以及其力学特征等方面的准确性和高效性。 展开更多
关键词 圆化多面体离散单元 有限元法 变形单元 等效节点力
下载PDF
一种基于单元几何变形的网格修匀方法 被引量:1
3
作者 孙树立 苑远 +1 位作者 勾志宏 张明磊 《计算力学学报》 CAS CSCD 北大核心 2016年第4期478-484,共7页
基于单元几何变形操作提出一种高效的非结构网格质量修匀方法。其基本过程是先对每个单元独立地进行拉伸-收缩操作以优化单元的形状,然后在整个网格中通过对各单元的节点位置进行加权平均来获得改善后的网格。为进一步提高修匀方法对网... 基于单元几何变形操作提出一种高效的非结构网格质量修匀方法。其基本过程是先对每个单元独立地进行拉伸-收缩操作以优化单元的形状,然后在整个网格中通过对各单元的节点位置进行加权平均来获得改善后的网格。为进一步提高修匀方法对网格质量的优化效果,并使得该方法具备一定的网格调整能力,结合动网格技术提出了对单元进行大范围和较大幅度移动的策略;在修匀过程中还通过适当算法调整单元形心位置和单元尺寸,进一步增强了修匀方法对网格局部进行疏密调节的能力。本文方法可适用于平面和三维非结构网格的质量改善及网格调整。若干算例表明了方法的有效性。 展开更多
关键词 网格修匀 网格调整 单元几何变形 拉伸-收缩操作 动网格
下载PDF
反求工程在变形吸能单元国产化中的应用研究
4
作者 吕志龙 闵阳春 刘晓波 《电力机车与城轨车辆》 2008年第3期11-12,46,共3页
阐述了反求工程的概念,介绍反求工程在HXD1型机车变形吸能单元国产化中的成功应用,提出了在我国机车车辆跨越式发展中应用反求工程的建议。
关键词 反求工程 HXD1型机车 变形吸能单元
下载PDF
可变形圆化多边形离散单元法 被引量:2
5
作者 周洋诗琦 赵兰浩 +1 位作者 邵琳玉 毛佳 《岩土力学》 EI CAS CSCD 北大核心 2022年第7期1961-1968,1977,共9页
将基于圆化多边形离散单元法与有限元方法结合,提出一种可变形圆化多边形离散单元法。此法对块体离散元进行圆化处理,可较好地表征不规则块体外形,又保留了颗粒离散元计算高效的优势。在求解接触力时,消除了角点处法向奇异等问题,同时... 将基于圆化多边形离散单元法与有限元方法结合,提出一种可变形圆化多边形离散单元法。此法对块体离散元进行圆化处理,可较好地表征不规则块体外形,又保留了颗粒离散元计算高效的优势。在求解接触力时,消除了角点处法向奇异等问题,同时增强计算的稳定性和简化接触判断。同时对切向接触力计算模型进行修正,使得接触力计算效率得到提高。此法突破了圆化多边形刚体假设的限制,可以精确计算任意形状不规则离散单元之间的相互作用,对单元的运动和变形进行模拟。通过超静定梁冲击试验、不规则块体单轴压缩试验和料斗流动“卡阻”试验3个数值模拟算例,论证此法可以有效地捕捉单元的碰撞、分离和变形等空间运动和自身特性以及其细观力学表征。 展开更多
关键词 圆化多边形离散单元 变形单元 接触力
下载PDF
新型参数化可变形单元体建筑结构研究与应用 被引量:1
6
作者 徐孝源 贾慧娟 席槿雯 《四川水泥》 2022年第12期48-50,共3页
新型参数化可变形单元体构件,基于球铰点连接进行快速组装,可通过多种搭接方式组成不同的空间结构,从而形成适用于不同环境和功能需求的建筑结构。新型参数化可变形单元体建筑具有快速拆装性能,配套了简化施工步骤节约施工时间的快速拆... 新型参数化可变形单元体构件,基于球铰点连接进行快速组装,可通过多种搭接方式组成不同的空间结构,从而形成适用于不同环境和功能需求的建筑结构。新型参数化可变形单元体建筑具有快速拆装性能,配套了简化施工步骤节约施工时间的快速拆装施工技术。该技术可通过对基本单元体角度的计算及连接排布方式参数的调整,最终达到结构、功能、环境、美感的四者统一。 展开更多
关键词 新型建筑结构 参数化 变形单元 快速施工技术
下载PDF
变形加载单元阵列的传输特性
7
作者 陈军文 黄培康 林桂森 《系统工程与电子技术》 EI CSCD 1990年第3期28-33,共6页
频率选择表面(FSS)对空间电磁波具有滤波特性,它广泛用于隐身飞行器、多频天线系统等。 本文提出了一种变形加载单元,目的在于减小一般加载单元的交叉极化电平。对使用变形加载缝隙组成的单层FSS和带有三层介质板的双层FSS进行了仿真计... 频率选择表面(FSS)对空间电磁波具有滤波特性,它广泛用于隐身飞行器、多频天线系统等。 本文提出了一种变形加载单元,目的在于减小一般加载单元的交叉极化电平。对使用变形加载缝隙组成的单层FSS和带有三层介质板的双层FSS进行了仿真计算,设计的双层FSS具有稳定的谐振频率和带宽,传输特性近似矩形,且有较窄的带宽(1GHz左右)和平坦的顶部。 单层FSS测试曲线表明了理论计算和计算机仿真的正确性。 展开更多
关键词 变形加载单元 阵列 传输特性 FSS
下载PDF
HX_D1型机车国产化变形单元非正常变形分析及解决方案
8
作者 欧阳黎健 何清和 《电力机车与城轨车辆》 2012年第3期71-73,共3页
文章针对HXD1型机车2万吨运输中出现的变形单元非正常变形的现象,从原理和实际运行等方面进行了分析,并提出了有效的结构优化解决方案。
关键词 HXD1型机车 变形单元 重载 吸能
下载PDF
97单元MEMS变形镜波前像差拟合能力分析
9
作者 徐奇 吴阳 +2 位作者 方舟 杨海波 杨慧珍 《光学技术》 CAS CSCD 北大核心 2020年第1期1-6,共6页
变形镜是自适应光学系统的关键元件,它对波前像差的拟合能力决定了自适应光学系统的校正性能。文中从变形镜对Zernike单阶像差、组合像差以及闭环校正三个方面分析97单元变形镜的拟合能力。仿真结果表明,当原始波前均方根的值为一个波长... 变形镜是自适应光学系统的关键元件,它对波前像差的拟合能力决定了自适应光学系统的校正性能。文中从变形镜对Zernike单阶像差、组合像差以及闭环校正三个方面分析97单元变形镜的拟合能力。仿真结果表明,当原始波前均方根的值为一个波长时,Zernike多项式前3~42阶残余波前像差的RMS小于0.4λ,说明变形镜对Zernike的前3~42阶像差具有较好的拟合能力。变形镜对Zernike多项式组合相差拟合以及基于随机并行梯度下降算法的闭环校正结果表明,当波前像差较小时,像差基本得到完全拟合及校正,当波前像差较大时,如D/r0=20时,残余波前像差的RMS值均小于0.14λ(初始RMS为0.63λ)。分析结果对97单元变形镜的实际应用提供了理论依据和使用参考。 展开更多
关键词 激光技术 97单元MEMS变形 校正能力 SPGD
下载PDF
杂交元假设应力模式的变形刚度分析 被引量:7
10
作者 张灿辉 黄黔 冯伟 《应用数学和力学》 CSCD 北大核心 2006年第7期757-764,共8页
通过本征变形模式提出识别杂交元零能变形模式和假设应力场中零能应力模式的新方法,同时给出了在假设应力场中增加应力模式时杂交元变形刚度的计算公式.从而从理论上阐明了在假设应力场中增加零能应力模式不仅不能抑制单元零能变形模式... 通过本征变形模式提出识别杂交元零能变形模式和假设应力场中零能应力模式的新方法,同时给出了在假设应力场中增加应力模式时杂交元变形刚度的计算公式.从而从理论上阐明了在假设应力场中增加零能应力模式不仅不能抑制单元零能变形模式而且可能增加非零能变形模式的刚度,因此不宜用来假设应力场;同时进一步指出寄生应力模式将使单元产生虚假应变能而使单元显得过刚,因此即使它能够抑制单元零能变形模式也不宜用来假设应力场,从而为假设应力场提供了合理的建议.数值算例说明了包含零能应力模式和寄生应力模式时单元的性能. 展开更多
关键词 杂交应力元 单元变形刚度 零能变形模式 零能应力模式
下载PDF
基于GPU的961单元自适应光学系统波前处理器设计及其动态性能 被引量:1
11
作者 杨乐强 王建立 +3 位作者 姚凯男 李宏壮 陈璐 邵蒙 《光学精密工程》 EI CAS CSCD 北大核心 2022年第23期3004-3012,共9页
为了满足4 m级地基大口径光学望远镜千单元自适应光学系统的波前处理规模、速度以及控制带宽要求,研究了基于GPU的千单元级自适应光学系统高速波前处理器。介绍了自适应光学系统波前处理方法,并详细论述了基于GPU的波前处理架构的实现... 为了满足4 m级地基大口径光学望远镜千单元自适应光学系统的波前处理规模、速度以及控制带宽要求,研究了基于GPU的千单元级自适应光学系统高速波前处理器。介绍了自适应光学系统波前处理方法,并详细论述了基于GPU的波前处理架构的实现与优化方法。最后,通过湍流模拟器进行桌面模拟实验,得到961单元自适应光学系统的动态性能。实验结果表明:基于GPU的961单元自适应光学系统在采样频率为1500 Hz时,系统0 dB残余误差抑制带宽能够达到100 Hz,满足4 m级地基大口径光学望远镜自适应光学系统的应用需求。 展开更多
关键词 自适应光学 波前处理 961单元变形 图形处理器
下载PDF
三维块体离散元可变形计算方法研究 被引量:14
12
作者 田振农 李世海 +1 位作者 刘晓宇 柳丙善 《岩石力学与工程学报》 EI CAS CSCD 北大核心 2008年第S1期2832-2840,共9页
在三维刚性块体离散元的基础上,研究可变形块体离散元方法,并编制相应的可变形块体离散元程序。该方法把块体单元看作弹性体,并对其划分有限元网格,对每个可变形单元的计算是假设在力边界条件下求解可变形单元的节点位移;可变形块体之... 在三维刚性块体离散元的基础上,研究可变形块体离散元方法,并编制相应的可变形块体离散元程序。该方法把块体单元看作弹性体,并对其划分有限元网格,对每个可变形单元的计算是假设在力边界条件下求解可变形单元的节点位移;可变形块体之间力和位移的传递是通过构筑的节理单元实现的,节理单元反应岩土工程中实际节理的变形性质和强度特征,结构面破坏前,体现了连续介质计算方法的特点,结构面如果发生破坏,预设节理处的节理单元直接转变化为法向和切向弹簧,体现了离散元计算方法的特点。另外,节理单元有实际的厚度,可以根据实际节理的性质确定计算参数,这一点也体现出了对刚性块体离散元方法的改进。与刚性块体离散元方法比较,该方法可以反映岩体变形性质,如泊松效应、三维岩体中的真实波速,它还可以模拟结构面的破坏演化过程;与有限元方法比较,它可以更好地描述岩体中节理的特征,允许单元有较大的位移;相对于以前的可变性离散元方法,它既克服了有限差分格式带来计算上的复杂,又建立在三维模型的基础上,为有关岩土工程的数值分析提供了一个可行的数值方法。 展开更多
关键词 数值分析 离散元法 岩体 变形单元
下载PDF
改进有限单元法动力特性方法的探讨
13
作者 江增荣 段鹏飞 +1 位作者 郭杏林 丁桦 《应用数学和力学》 CSCD 北大核心 2010年第11期1267-1275,共9页
采样定理告诉人们进行离散逼近离散尺寸的极限是问题中所涉及最高频率波长的1/2.1/2波长是一理想值,在目前的有限元分析中,通常保证精度的空间离散尺寸远小于这一理想值.文中分析了这一现象的可能原因,同时给出了一种基于有限元的在一... 采样定理告诉人们进行离散逼近离散尺寸的极限是问题中所涉及最高频率波长的1/2.1/2波长是一理想值,在目前的有限元分析中,通常保证精度的空间离散尺寸远小于这一理想值.文中分析了这一现象的可能原因,同时给出了一种基于有限元的在一定条件下能够改进空间离散效率的动力模拟方法. 展开更多
关键词 采样定理 有限元离散效率 基于变形修正的宏单元 频散效应
下载PDF
应用无网格伽辽金方法分析结构大变形问题 被引量:1
14
作者 陆新征 杨宁 江见鲸 《结构工程师》 2003年第z1期18-22,共5页
无单元法由于不需要复杂的网格划分,不存在网格畸变问题,因此在大变形分析领域有着广阔的应用前景.本文利用无网格伽辽金(EFGM)方法,对二维结构大变形问题进行了分析,得到了传统有限元方法所难以得到的结果.文中详细讨论了无网格伽辽金... 无单元法由于不需要复杂的网格划分,不存在网格畸变问题,因此在大变形分析领域有着广阔的应用前景.本文利用无网格伽辽金(EFGM)方法,对二维结构大变形问题进行了分析,得到了传统有限元方法所难以得到的结果.文中详细讨论了无网格伽辽金方法的基函数、权函数的选取及影响域的设定,并给出了各参数的具体取值.用计算实例说明了无网格伽辽金方法在解决结构大变形问题上的优势. 展开更多
关键词 无网格伽辽金方法 变形 单元
下载PDF
西西伯利亚盆地走滑断裂构造特征及其成因机制
15
作者 莫午零 计智锋 +2 位作者 姜仁 张艺琼 王素花 《天然气工业》 EI CAS CSCD 北大核心 2024年第8期1-12,共12页
西西伯利亚盆地是俄罗斯油气储量最大和产量最高含油气盆地,该盆地发育大量走滑断裂,关于其构造特征与形成机制的认识存在诸多争议。为此,基于高精度三维地震资料,在解析盆地内走滑断裂几何学、运动学特征的基础上,深入探讨了西西伯利... 西西伯利亚盆地是俄罗斯油气储量最大和产量最高含油气盆地,该盆地发育大量走滑断裂,关于其构造特征与形成机制的认识存在诸多争议。为此,基于高精度三维地震资料,在解析盆地内走滑断裂几何学、运动学特征的基础上,深入探讨了西西伯利亚盆地走滑断裂形成的动力学机制。研究结果表明:①走滑断裂主要发育北西向右行走滑断裂和北东向左行走滑断裂,2组走滑断裂之间锐夹角平均为84°,夹角平分线对应正南北方向,发育的雁列式断层与主位移带呈45°相交,基本都呈近南北向展布;②盆地存在水平运动和垂直运动,总体以南北向挤压水平运动占主导地位,其中走滑断裂性质以压剪作用为主,走滑位移幅度不显著,复杂的局部构造应力组合可能发生在断裂的末端和交叉点区域;③盆地北部发育的隆起主要为近等轴隆起(穹隆),走滑断裂以北西向为主,盆地南部发育的隆起主要为长条形(非等轴)隆起,走滑断裂发育北西向与北东向共轭断裂。结论认为:①盆地发育的走滑断裂形成机制是纯剪切模式,是同时在相互正交的水平截面中实现压缩和拉伸作用,各种性质的伴生构造和位移带为北半球沉积盆地中发育的2个稳定剪切应力场释放的结果;②在纯剪切运动模型条件下形成剪切变形单元,变形单元的一个或另一个边缘的激活,形成剪切区域的尾翼,为在不同(正交)方向的最大正应力轴上实现尾翼正交网格的不同构造激活阶段提供了条件。 展开更多
关键词 西西伯利亚盆地 构造特征 走滑断裂 变形单元 成因机制
下载PDF
型钢混凝土偏压柱粘结滑移性能的试验研究 被引量:2
16
作者 王彦宏 赵鸿铁 +1 位作者 薛建阳 杨勇 《西安建筑科技大学学报(自然科学版)》 CSCD 2003年第3期217-220,225,共5页
针对型钢混凝土粘结滑移性能,通过8根以荷载偏心距e0/h,长细比l0/h为主要参数的SRC偏心受压柱试验研究,分析了SRC柱正截面破坏的形态、特征,找出了其粘结滑移的影响因素和分布规律,并建立了合理的粘结滑移计算模型.本文可对型钢混凝土... 针对型钢混凝土粘结滑移性能,通过8根以荷载偏心距e0/h,长细比l0/h为主要参数的SRC偏心受压柱试验研究,分析了SRC柱正截面破坏的形态、特征,找出了其粘结滑移的影响因素和分布规律,并建立了合理的粘结滑移计算模型.本文可对型钢混凝土粘结滑移性能的进一步试验研究和理论分析提供一定的参考,并为了解其力学性能和型钢混凝土压弯构件的设计提供了理论依据. 展开更多
关键词 型钢混凝土结构 粘结滑移性能 偏压柱 单元变形
下载PDF
Deformation caused by the 2011 eastern Japan great earthquake monitored using the GPS single-epoch precise point positioning technique 被引量:4
17
作者 郭金运 原永东 +2 位作者 孔巧丽 李国伟 王方建 《Applied Geophysics》 SCIE CSCD 2012年第4期483-493,497,共12页
Crustal deformation can provide constraints for studying earthquake rupture and shock wave transmission for the Mw9.0 eastern Japan great earthquake. Using the single- epoch precise point positioning (PPP) method an... Crustal deformation can provide constraints for studying earthquake rupture and shock wave transmission for the Mw9.0 eastern Japan great earthquake. Using the single- epoch precise point positioning (PPP) method and the appropriate positioning flow, we process GPS data from six IGS (International GNSS Service) sites (e.g., MIZU, TSK2, USUD, MTKA, AIRA and KSMV) located in Japan and obtain the positioning results with centimeter scale precision. The displacement time series of the six sites are analyzed using the least squares spectral analysis method to estimate deformations caused by the Mw9.0 mainshock and the Mw7.9 aftershock, and the cumulative displacements after 1 day. Mainshock displacements at station MIZU, the nearest site to the mainshock in the North (N), East (E), and Up (U) directions, are -1.202 m, 2.180 m and -0.104 m, respectively, and the cumulative deformations after 1 day are -1.117 m, 2.071 m and -0.072 m, respectively. The displacements at station KSMV, the nearest site to the Mw7.9 aftershock in the N, E and U directions, are -0.032 m, 0.742 m and -0.345 m, respectively. The other sites obviously experienced eastern movements and subsidence. The deformation vectors indicate that the horizontal displacements caused by the earthquake point to the epicenter and rupture. Elastic bounds evidently took place at all sites. The results indicate that the crustal movements and earthquake were part of a megathrust caused by the Pacific Plate sinking under the North American Plate to the northeast of Japan island arc. 展开更多
关键词 Eastern Japan great earthquake GPS single-epoch precise point positioning crustal deformation
下载PDF
运动结构稳定构形切换实现景观构筑物的功能拓展 被引量:1
18
作者 王瑜 陈震 《深圳大学学报(理工版)》 EI CAS CSCD 北大核心 2017年第1期51-62,共12页
针对运动结构日益广泛用于景观构筑物的现象,从影响运动结构稳定构形状态的基本构成单元切入,根据其自身的形变,将运动结构的单元归纳为可变形单元与不可变形单元两类,不同的单元类型对应不同的稳定构形.研究表明,构成折叠结构的剪杆单... 针对运动结构日益广泛用于景观构筑物的现象,从影响运动结构稳定构形状态的基本构成单元切入,根据其自身的形变,将运动结构的单元归纳为可变形单元与不可变形单元两类,不同的单元类型对应不同的稳定构形.研究表明,构成折叠结构的剪杆单元、刚性面板单元与索杆单元,构成充气膜结构的薄膜,构成开合结构的索膜均属于可变形单元;界面开合、单元重组与整体移动则是基于不可变形单元实现的稳定构形转换.稳定构形间的灵活切换适应了不同用途下的特定要求,在较短时间内实现了构筑物功能的改变与拓展. 展开更多
关键词 建筑设计 稳定构形 运动结构 景观构筑物 变形单元 不可变形单元
下载PDF
筒型基础安装过程中筒内土塞生成的可视化模拟计算 被引量:4
19
作者 丁红岩 张浦阳 《中国海洋大学学报(自然科学版)》 CAS CSCD 北大核心 2007年第6期1034-1038,共5页
筒型基础安装过程中,筒内的负压作用会导致土塞现象的出现。以试验模型为基础建立可变形离散元计算模型,通过VC++编制的程序SPSAⅡ对筒型基础负压沉贯中土塞的生成过程进行可视化模拟计算。考虑渗流力、筒壁内外摩擦阻力和筒内负压吸力... 筒型基础安装过程中,筒内的负压作用会导致土塞现象的出现。以试验模型为基础建立可变形离散元计算模型,通过VC++编制的程序SPSAⅡ对筒型基础负压沉贯中土塞的生成过程进行可视化模拟计算。考虑渗流力、筒壁内外摩擦阻力和筒内负压吸力作为土塞产生的主要外荷载,分别模拟渤海地区3种典型土(黏土、粉土和砂土)的情况。其中,粉土情况下,程序的收敛速度最快,黏土情况下,土塞的计算高度最接近试验测量值,数值计算结果表明该计算方法作为筒型基础施工前筒内土塞高度的预测方法是可行的。 展开更多
关键词 筒型基础 土塞 变形离散单元
下载PDF
自适应光学
20
《中国光学》 EI CAS 1999年第4期99-100,共2页
O43 99042793自适应光学波前校正器=Wavefront correctors ofadaptive optics[刊,中]/凌宁(中科院光电技术研究所.四川,成都(610209))//光学技术.—1998,(3).—12—16,40报告了自适应光学波前校正器—多元分立式压电变形反射镜和两维高... O43 99042793自适应光学波前校正器=Wavefront correctors ofadaptive optics[刊,中]/凌宁(中科院光电技术研究所.四川,成都(610209))//光学技术.—1998,(3).—12—16,40报告了自适应光学波前校正器—多元分立式压电变形反射镜和两维高速压电倾斜镜的研究工作。讨论了这两类波前校正器的主要性能及其检测方法等,并给出了若干典型的检测结果。研制的19单元度形镜已成功用于ICF驱动器的波前校正,而21单元变形镜已成功用于天文观测星体成像校正。图13表1参7(方舟) 展开更多
关键词 自适应光学系统 波前校正器 高速压电倾斜镜 变形反射镜 主要性能 技术研究所 检测方法 光学技术 天文观测 单元变形
下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部