期刊文献+
共找到660篇文章
< 1 2 33 >
每页显示 20 50 100
不同雷诺数下矩形钝体尾流中悬臂梁流致振动的实验研究
1
作者 王金林 霍昕芸 +1 位作者 闵为 王昕华 《液压气动与密封》 2025年第4期1-9,共9页
应用粒子图像测速技术(Particle Image Velocimetry,PIV),对雷诺数Re范围为3018~10565及悬臂梁与钝体间距s=3D条件下的矩形钝体(长宽比为1/2)尾流场特征,以及悬臂梁的振动特性进行了研究,并利用本征正交分解法(Proper Orthogonal Decomp... 应用粒子图像测速技术(Particle Image Velocimetry,PIV),对雷诺数Re范围为3018~10565及悬臂梁与钝体间距s=3D条件下的矩形钝体(长宽比为1/2)尾流场特征,以及悬臂梁的振动特性进行了研究,并利用本征正交分解法(Proper Orthogonal Decomposition,POD)对实验结果进行分析。研究结果表明,钝体尾流场前二阶POD模态含能最大并且模态频率相同,反映了大尺度的卡门涡结构,其模态频率即为涡脱频率f,该值与悬臂梁的振动主频f b一致,表明钝体尾流区悬臂梁的振动激励来自于卡门涡;矩形柱钝体尾流中悬臂梁的均方根振幅y rms与振动主频f b基本随Re的增大呈线性增大。利用前二阶POD模态对尾流场进行模态重构,揭示了矩形钝体尾流中悬臂梁对流场结构的反作用以及悬臂梁的振动机理。 展开更多
关键词 矩形钝体 悬臂梁 本征正交分解 粒子图像测速 流致振动
在线阅读 下载PDF
基于BN稳定性的非线性振动问题时间积分方法
2
作者 季奕 邢誉峰 刘波 《力学与实践》 2025年第5期1040-1045,共6页
BN稳定型方法是求解动力学问题的一类新型时间积分方法,可稳定、准确且高效地求解复杂非线性振动问题。自2021年提出后,这一类求解方法已受到学术界关注。本文对BN稳定型方法的设计思路、求解格式及其参数优化设计原则进行了讨论,并利... BN稳定型方法是求解动力学问题的一类新型时间积分方法,可稳定、准确且高效地求解复杂非线性振动问题。自2021年提出后,这一类求解方法已受到学术界关注。本文对BN稳定型方法的设计思路、求解格式及其参数优化设计原则进行了讨论,并利用高阶非线性单自由度系统和刚性卫星系统作为案例,将其与广义-α方法和ρ∞-Bathe方法进行比较。结果表明,BN稳定型方法的精度和稳定性优势。本文工作可为振动力学课程提供创新素材,有助于培养学生创新能力和对非线性振动问题的分析能力。 展开更多
关键词 BN稳定性 非线性 振动 时间积分方法
在线阅读 下载PDF
基于双线性插值的多尺度耦合一致性与收敛性研究
3
作者 李孜 杨庆成 《力学季刊》 北大核心 2025年第3期631-648,共18页
界面一致性与误差收敛性是并发多尺度计算方法中的核心问题,尤其在原子-连续介质耦合模型中尤为关键.然而,现有理论研究工作较少并且多聚焦于一维情况.本文以多分辨率分子力学方法为研究对象,通过构建原子区域与粗化区域共存的二维方形... 界面一致性与误差收敛性是并发多尺度计算方法中的核心问题,尤其在原子-连续介质耦合模型中尤为关键.然而,现有理论研究工作较少并且多聚焦于一维情况.本文以多分辨率分子力学方法为研究对象,通过构建原子区域与粗化区域共存的二维方形与三角形原子晶格模型,系统探究了双线性单元插值条件下不同能量采样方案对界面一致性和误差收敛性的影响.研究发现:界面次采样方案可显著提升界面区域一致性,其中包含所有近邻层的采样方案性能最优;误差分析结果显示,离散误差是总误差的主要来源,增加次采样点数量可有效降低采样误差,尤其在拉伸工况下表现显著.此外,两类晶格结构在误差收敛性方面展现出一致规律,验证了方法在不同晶格结构中的普适性.本研究显示了多分辨率分子力学方法中能量采样方案在处理界面一致性与收敛性方面的优势,为发展高精度跨尺度计算力学方法提供了理论支撑. 展开更多
关键词 并发多尺度计算方法 多分辨率分子力学 界面一致性 收敛性分析 双线性单元 能量采样
原文传递
流‒固‒相多场耦合数值方法及水力压裂模拟
4
作者 王飞阳 《同济大学学报(自然科学版)》 北大核心 2025年第12期1848-1857,共10页
基于统一断裂相场理论,提出了流‒固‒相多场耦合数值方法,并将其用于探究赋存地应力和天然裂缝对水力压裂及其裂缝形态的影响。数值模拟研究表明:流‒固‒相多场耦合数值方法能准确预测拉、剪和拉‒剪复合等水力压裂的裂缝形态;起裂所需的... 基于统一断裂相场理论,提出了流‒固‒相多场耦合数值方法,并将其用于探究赋存地应力和天然裂缝对水力压裂及其裂缝形态的影响。数值模拟研究表明:流‒固‒相多场耦合数值方法能准确预测拉、剪和拉‒剪复合等水力压裂的裂缝形态;起裂所需的压力水头与赋存地应力水平正相关,且水力压裂的裂缝形态受侧压力系数影响明显;天然裂缝改变了孔压的分布,进而影响水力压裂裂缝的分叉和交汇,形成不同的裂缝形态。 展开更多
关键词 水力压裂 流‒固‒相多场耦合 水平集方法 赋存地应力 天然裂缝
在线阅读 下载PDF
冻结温度下岩体THM完全耦合的理论初步分析 被引量:25
5
作者 徐光苗 刘泉声 张秀丽 《岩石力学与工程学报》 EI CAS CSCD 北大核心 2004年第21期3709-3713,共5页
在冻结温度下,岩体受内部孔隙水、结构面及温度等的影响,其物理力学性质将发生深刻的变化,因而,一般热力学理论已经难以解释该问题。从不可逆过程热力学和连续介质力学理论出发,在若干合理的假定条件下,推导了冻结温度下岩体的质量守恒... 在冻结温度下,岩体受内部孔隙水、结构面及温度等的影响,其物理力学性质将发生深刻的变化,因而,一般热力学理论已经难以解释该问题。从不可逆过程热力学和连续介质力学理论出发,在若干合理的假定条件下,推导了冻结温度下岩体的质量守恒方程、平衡方程及能量守恒方程的最终表示形式;研究了岩体在冻结温度下冰与岩石的膨胀耦合关系,并建立了相应的耦合膨胀系数的表达式。为研究寒区工程围岩温度场、渗流场及应力场的分布规律及变形机理提供了新的理论依据。 展开更多
关键词 冻结温度 岩体 结构面 物理力学性质 孔隙水 岩石 不可逆过程热力学 连续介质力学 热力学理论 平衡方程
在线阅读 下载PDF
微尺度流动研究中的几个问题 被引量:66
6
作者 陶然 权晓波 徐建中 《工程热物理学报》 EI CAS CSCD 北大核心 2001年第5期575-577,共3页
本文综述了当前微尺度流动研究的进展,对由于尺度微小化所带来的有别于宏观尺度流动的一些问题进行了讨论.分析了在微尺度流动问题中必须注意的一些因素,例如,尺度划分与连续介质假设的适用范围,表面效应,稀薄性与压缩性的作用,... 本文综述了当前微尺度流动研究的进展,对由于尺度微小化所带来的有别于宏观尺度流动的一些问题进行了讨论.分析了在微尺度流动问题中必须注意的一些因素,例如,尺度划分与连续介质假设的适用范围,表面效应,稀薄性与压缩性的作用,以及流动通道壁面的影响.最后对今后的研究进行了展望. 展开更多
关键词 微尺度流动 尺度划分 连续介质 微装置
在线阅读 下载PDF
大型地下洞室群围岩应力-损伤-渗流耦合分析 被引量:20
7
作者 张巍 肖明 范国邦 《岩土力学》 EI CSCD 北大核心 2008年第7期1813-1818,共6页
以渗透体力来考虑渗流场的力学效应,建立了应力-损伤-渗透系数关系方程来考虑应力和损伤对渗流场的影响,结合岩体结构的三维弹塑性损伤有限元分析,建立了大型地下洞室开挖围岩应力-损伤-渗流耦合的计算模型。该计算模型求解的难度主要... 以渗透体力来考虑渗流场的力学效应,建立了应力-损伤-渗透系数关系方程来考虑应力和损伤对渗流场的影响,结合岩体结构的三维弹塑性损伤有限元分析,建立了大型地下洞室开挖围岩应力-损伤-渗流耦合的计算模型。该计算模型求解的难度主要体现在岩体材料弹塑性、损伤、渗流自由面边界、渗流溢出边界、应力-损伤-渗流相互影响关系等。提出分步迭代法对以上因素进行归纳后按一定顺序分别进行迭代求解,取得了良好的计算效果。将该方法应用于某水电站大型地下洞室群分析,得出了一系列有意义的结论。 展开更多
关键词 地下洞室群 围岩 应力 损伤 渗流 耦合 分步迭代法
在线阅读 下载PDF
弹性介质中充液管道的波衰减特性 被引量:12
8
作者 刘敬喜 李天匀 +1 位作者 刘土光 李强 《华中科技大学学报(自然科学版)》 EI CAS CSCD 北大核心 2003年第10期90-92,共3页
从壳体基本运动方程出发 ,导出了弹性介质中的充液管道在轴对称运动下声振耦合系统的频散方程 .利用数值方法得到了频散方程的完全解 ,对弹性介质中充液钢质管道和PVC塑料管道波传播衰减特性进行了研究 .
关键词 充液管道 衰减 弹性介质 频散方程
在线阅读 下载PDF
微/纳尺度接触问题计算方法研究进展 被引量:12
9
作者 刘更 刘天祥 温诗铸 《力学进展》 EI CSCD 北大核心 2008年第5期521-544,共24页
接触问题广泛存在于现实生活的众多领域,近来随着微/纳米技术的不断发展,接触力学在基础理论和研究方法上面临许多新的挑战.本文在摩擦学的范畴内,对近年发展的若干求解微/纳尺度接触问题的计算方法及理论进行了综述.按发展先后及所解... 接触问题广泛存在于现实生活的众多领域,近来随着微/纳米技术的不断发展,接触力学在基础理论和研究方法上面临许多新的挑战.本文在摩擦学的范畴内,对近年发展的若干求解微/纳尺度接触问题的计算方法及理论进行了综述.按发展先后及所解决问题的尺度范围划分,主要有3类评估微/纳尺度接触性能的计算方法:(1)连续介质力学方法;(2)分子动力学模拟;(3)多尺度方法.介绍了这3类计算方法的典型理论和主要数学描述,给出了这些方法对解决若干微/纳观接触问题如黏着效应、粗糙表面描述、表面摩擦及润滑、表面热效应、生物接触等的主要应用.最后,探讨了微/纳尺度接触问题计算方法可能的发展方向及应用领域. 展开更多
关键词 微/纳尺度 接触 计算方法 连续介质 分子动力学 多尺度
在线阅读 下载PDF
流固耦合问题的描述方法及分类简化准则 被引量:48
10
作者 朱洪来 白象忠 《工程力学》 EI CSCD 北大核心 2007年第10期92-99,共8页
在流固耦合的工程实际问题中的大多数情况下,弹性薄壁构件的变形为几何非线性,再加上流体方程的非线性,将导致流体和弹性体相互作用界面上的强非线性。在界面上便可以结合拉格朗日法和欧拉法建立方程和接触条件。其方法主要有单一拉格... 在流固耦合的工程实际问题中的大多数情况下,弹性薄壁构件的变形为几何非线性,再加上流体方程的非线性,将导致流体和弹性体相互作用界面上的强非线性。在界面上便可以结合拉格朗日法和欧拉法建立方程和接触条件。其方法主要有单一拉格朗日法、单一欧拉法、相容拉格朗日-欧拉法和任意拉格朗日-欧拉法四种方法。引入描述弹性体变形特征的数值m、n、k和描述流体运动特征的数值λ、ν,可将流固耦合问题进行分类。在流体弹性力学理论的基础上,介绍了流固耦合问题界面相互作用的描述方法,并根据诺沃日洛夫BB(НовожиловВВ)在非线性弹性力学中,从几何非线性方面对弹性力学问题的分类方法出发,对流体弹性力学中的流固耦合问题进行分类,由此,可为按类别对运动学条件、动力学条件及界面上的接触条件进行相应的简化提供可靠的依据。 展开更多
关键词 流固耦合 相容拉格朗日-欧拉法 单一拉格朗日法 任意拉格朗日-欧拉法 几何非线性 分类简化准则
在线阅读 下载PDF
近场动力学与有限元的混合建模方法 被引量:11
11
作者 章青 郁杨天 顾鑫 《计算力学学报》 CAS CSCD 北大核心 2016年第4期441-448,450,共9页
综述了近场动力学与有限元混合建模方法的研究进展,阐明了各种混合建模方法的基本原理与特点,并重点介绍本课题组在近场动力学与有限元方法混合建模方面的研究工作。现有近场动力学与有限元混合建模方法包括位移协调约束、力耦合、混合... 综述了近场动力学与有限元混合建模方法的研究进展,阐明了各种混合建模方法的基本原理与特点,并重点介绍本课题组在近场动力学与有限元方法混合建模方面的研究工作。现有近场动力学与有限元混合建模方法包括位移协调约束、力耦合、混合函数方法以及子模型方法等,除子模型方法外,都可归结为并行式多尺度分析方法,其基本思想是将计算结构划分为近场动力学子域、有限元子域以及两者的交界区域(或重叠区域、或界面单元、或过渡区域)。子模型方法可归结为显-显分析方法,先采用显式有限元进行整体分析,后采用近场动力学方法对重点区域进行分析。混合建模方法需要着重提高交界区域的计算精度,并且消除虚假力和虚假应力波问题。提出了通过力耦合的近场动力学与有限元混合建模的隐式分析方法,该方法不再设置重叠区,通过杆单元连接近场动力学子域与有限元子域,其中界面上的有限元结点不仅与其所在单元的其他结点发生作用,还通过杆单元与以其为圆心、一定半径的圆域内的其他物质点相互作用。研究表明,本文提出的混合模型和求解方法既能有效解决裂纹扩展等不连续问题,又可提高计算效率,为工程结构破坏问题的计算分析提供一种有效方法。 展开更多
关键词 近场动力学 有限元 混合建模 非局部 不连续力学问题
在线阅读 下载PDF
空间杆系结构实用几何非线性分析 被引量:12
12
作者 丁泉顺 陈艾荣 项海帆 《力学季刊》 CSCD 北大核心 2001年第3期300-306,共7页
从简单实用的角度论述了空间杆系结构的几何非线性分析方法。文中分析了非线性有限元方法的求解过程,特别强调决定几何非线性收敛结果的关键问题,即由节点位移增量计算单元的内力增量。通过引入随转坐标系,论述了平面和空间梁单元小应... 从简单实用的角度论述了空间杆系结构的几何非线性分析方法。文中分析了非线性有限元方法的求解过程,特别强调决定几何非线性收敛结果的关键问题,即由节点位移增量计算单元的内力增量。通过引入随转坐标系,论述了平面和空间梁单元小应变时单元内力增量的计算问题。针对杆系结构的大应变问题,从有限应变理论出发进行分析,提出了对该问题的有效处理方法,并且用实例进行了验证。计算结果表明,该实用几何非线性分析方法是可靠和有效的。 展开更多
关键词 空间杆系结构 几何非线性 实用分析 小应变理论 大应变 有限元 单元内力增量 随转坐标系 现代工程结构 非线性连续介质力学
在线阅读 下载PDF
应变梯度理论进展 被引量:39
13
作者 陈少华 王自强 《力学进展》 EI CSCD 北大核心 2003年第2期207-216,共10页
应变梯度理论是近10年来为解释材料在微米尺度下的尺寸效应现象而发展起来的一种新理论.首先综述了应变梯度理论近年的发展及其对材料力学行为研究方面的进展.其次主要介绍了不含高阶应力的一类应变梯度理论及其应用;最后对应变梯度理... 应变梯度理论是近10年来为解释材料在微米尺度下的尺寸效应现象而发展起来的一种新理论.首先综述了应变梯度理论近年的发展及其对材料力学行为研究方面的进展.其次主要介绍了不含高阶应力的一类应变梯度理论及其应用;最后对应变梯度理论的发展做了展望. 展开更多
关键词 连续介质力学 应变梯度理论 研究进展 尺寸效应 微米尺度
在线阅读 下载PDF
考虑位移要求的大型三维连续体结构拓扑优化数值技术研究 被引量:10
14
作者 荣见华 唐国金 +1 位作者 罗银燕 杨端生 《工程力学》 EI CSCD 北大核心 2007年第3期20-27,19,共9页
大型复杂三维结构拓扑优化设计既具有理论意义,又具有重要的应用价值。基于等效转换的非奇异的结构优化模型,研究结构位移要求的最小结构重量设计问题。首先,介绍了位移约束的三维结构优化准则和公式。而后,为了提高拥有数万个单元以上... 大型复杂三维结构拓扑优化设计既具有理论意义,又具有重要的应用价值。基于等效转换的非奇异的结构优化模型,研究结构位移要求的最小结构重量设计问题。首先,介绍了位移约束的三维结构优化准则和公式。而后,为了提高拥有数万个单元以上的三维结构的计算效率,结合结构位移计算的迭代方法,在分析用于结构特性参数计算模型的基础上,建立了一套三维结构拓扑优化的求解策略和算法。最后,给出了几个典型和复杂的三维结构的拓扑优化设计算例。算例表明求解策略和算法是正确和有效的,且具有广泛的工程应用前景。 展开更多
关键词 工程力学 结构优化 结构拓扑优化 三维结构优化 数值方法
在线阅读 下载PDF
金属橡胶弹性迟滞回线的一种新的建模方法 被引量:10
15
作者 敖宏瑞 姜洪源 +1 位作者 夏宇宏 A M Ulanov 《中国矿业大学学报》 EI CAS CSCD 北大核心 2004年第4期453-456,共4页
金属橡胶作为一种干摩擦阻尼材料,其阻尼特性在很多领域得到了应用.在振动防护领域,通常用金属橡胶阻尼元件的迟滞回线来描述金属橡胶的干摩擦阻尼特性.本文提出了一种描述迟滞回线的数学方法,该方法采用泊松分布来描述金属橡胶材料内... 金属橡胶作为一种干摩擦阻尼材料,其阻尼特性在很多领域得到了应用.在振动防护领域,通常用金属橡胶阻尼元件的迟滞回线来描述金属橡胶的干摩擦阻尼特性.本文提出了一种描述迟滞回线的数学方法,该方法采用泊松分布来描述金属橡胶材料内部金属丝之间的接触点运动状态,并针对不同的振幅对金属橡胶力学本构关系中的参数进行了求解.与现有方法相比,该方法能够更加准确地描述金属橡胶的非线性迟滞特性. 展开更多
关键词 金属橡胶 迟滞回线 泊松分布 干摩擦阻尼
在线阅读 下载PDF
考虑压头曲率半径和应变梯度的微压痕分析 被引量:6
16
作者 陶彩军 王自强 陈少华 《力学学报》 EI CSCD 北大核心 2004年第6期680-687,共8页
在压头尖端曲率半径取100nm的前提下,采用Chen和Wang的应变梯度理论,对微压痕实验进行了系统的数值分析.首先通过拟合载荷-位移实验曲线的后半段来确定材料的屈服应力和幂硬化指数值,然后用有限元方法数值模拟压痕实验,并将计算得到的... 在压头尖端曲率半径取100nm的前提下,采用Chen和Wang的应变梯度理论,对微压痕实验进行了系统的数值分析.首先通过拟合载荷-位移实验曲线的后半段来确定材料的屈服应力和幂硬化指数值,然后用有限元方法数值模拟压痕实验,并将计算得到的整段载荷-位移曲线及硬度-位移曲线和实验结果进行了比较.结果表明应变梯度理论所预测的计算结果和实验结果很好地符合,包括压痕深度在亚微米和微米范围内的整段曲线. 展开更多
关键词 压痕实验 实验曲线 曲率半径 应变梯度理论 载荷-位移曲线 屈服应力 有限元方法 压头 计算结果 数值分析
在线阅读 下载PDF
振动磨非线性碰撞振动的建模与分析 被引量:14
17
作者 王树林 李巨光 +1 位作者 胡沂清 高秋芝 《青岛海洋大学学报(自然科学版)》 CSCD 1998年第3期157-164,共8页
振动磨系统的非线性主要来自散体介质的变质量碰撞,其中主碰撞和亚谐碰撞对系统动力学影响最为显著。主碰撞不但具有时变性,而且与空间坐标有关。基于研究的新进展,考虑到散体重心的动态偏移,提出了三自由度非线性振动模型,用平均... 振动磨系统的非线性主要来自散体介质的变质量碰撞,其中主碰撞和亚谐碰撞对系统动力学影响最为显著。主碰撞不但具有时变性,而且与空间坐标有关。基于研究的新进展,考虑到散体重心的动态偏移,提出了三自由度非线性振动模型,用平均法进行了求解分析,用时间序列分析的方法估计了变质量模型参数。分析求得机器振幅与实测值吻合较好,验证了所提出的理论和方法。 展开更多
关键词 非线性振动 变质量碰撞 系统建模 振动磨
原文传递
基于非线性本构关系的软管编织增强层等效刚度修正方法研究 被引量:7
18
作者 胡牧原 郑百林 +3 位作者 贺鹏飞 马朋升 张园春 郑立新 《力学季刊》 CSCD 北大核心 2015年第1期70-80,共11页
基于Hachemi提出的金属编织层本构理论以及拉伸实验的研究方法,进行了不同尺寸金属编织软管的拉伸实验.发现金属编织层在拉伸时会出现强烈的结构非线性,而该理论不能在非线性段与实验结果相吻合.提出了对理论中编织层刚度矩阵修正方法,... 基于Hachemi提出的金属编织层本构理论以及拉伸实验的研究方法,进行了不同尺寸金属编织软管的拉伸实验.发现金属编织层在拉伸时会出现强烈的结构非线性,而该理论不能在非线性段与实验结果相吻合.提出了对理论中编织层刚度矩阵修正方法,将金属纤维间交叠产生的接触等效地转化为"修正基体",独立于金属纤维自身产生的刚度.同时,结合编织角在实验中呈现的变化趋势,提出管状金属编织层在等位移拉伸荷载下编织角呈线性变化的假设.这与Hachemi的假设恰好相反.Hachemi认为编织层受拉时,编织角减小到一定程度则会发生锁定现象.而实验证明了编织角并不会发生锁定.且直到编织层破坏为止,编织角变化的趋势基本呈线性.还提出了修正编织角变化的加速系数k,通过金属纤维间的侧向接触解释了实验与数值仿真结果中,力位移曲线与编织角变化曲线不能同时吻合的力学机理. 展开更多
关键词 编织 非线性本构关系 修正基体 加速系数
原文传递
准脆性岩石水力耦合不排水多尺度本构模型 被引量:7
19
作者 朱其志 王岩岩 +2 位作者 仇晶晶 赵伦洋 邵建富 《河海大学学报(自然科学版)》 CAS CSCD 北大核心 2018年第2期165-170,共6页
以水-力耦合条件下准脆性岩石的损伤行为为研究对象,充分考虑材料的细观特征及其演化规律,结合Mori-Tanaka线性均匀化方法和不可逆热动力学理论建立了在不排水条件下水-力耦合损伤多尺度本构模型。基于小扰动假定,通过水-力耦合分析,推... 以水-力耦合条件下准脆性岩石的损伤行为为研究对象,充分考虑材料的细观特征及其演化规律,结合Mori-Tanaka线性均匀化方法和不可逆热动力学理论建立了在不排水条件下水-力耦合损伤多尺度本构模型。基于小扰动假定,通过水-力耦合分析,推导了不排水常规三轴试验条件下多尺度模型的强度的解析表达式。通过模拟不同围压下砂岩的三轴压缩不排水试验,验证了模型的有效性,并分析了岩石在损伤演化过程中应力应变曲线和孔隙水压力的演化规律。 展开更多
关键词 细观力学 脆性岩石 水力耦合模型 不排水条件
在线阅读 下载PDF
力学专业研究生计算力学能力培养及其课程体系建设 被引量:7
20
作者 杨庆生 龙连春 +1 位作者 刘赵淼 叶红玲 《力学与实践》 CSCD 北大核心 2012年第4期66-69,共4页
提出了力学专业研究生的计算力学能力要素,介绍了培养研究生计算力学能力和相关课程设置的一些新理念和做法.初步实践表明,具有高水平计算力学能力的高层次应用型创新人才受到工业技术领域的青睐.
关键词 计算力学 能力培养 课程体系 研究生
在线阅读 下载PDF
上一页 1 2 33 下一页 到第
使用帮助 返回顶部