期刊文献+
共找到18篇文章
< 1 >
每页显示 20 50 100
流体与结构相互作用问题的FE-SPH耦合模拟 被引量:9
1
作者 肖毅华 韩旭 胡德安 《应用力学学报》 CAS CSCD 北大核心 2011年第1期13-18,106,共6页
应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题。流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理。为了避免隐式计算压力,通过引入人... 应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题。流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理。为了避免隐式计算压力,通过引入人工压缩率,将不可压流体近似为人工可压缩流体。采用FE-SPH耦合法对弹性板在随时间变化的水压作用下的变形以及倒塌水柱冲击弹性结构两个问题进行了模拟。模拟结果与实验结果以及其他已有数值结果符合良好,说明FE-SPH耦合法用于流体与结构相互作用问题的模拟是可行和有效的。 展开更多
关键词 流体结构相互作用 有限元-光滑粒子流体动力学耦合 接触算法 人工压缩率
在线阅读 下载PDF
土-结构-流体动力相互作用的实时耦联动力试验 被引量:10
2
作者 迟福东 王进廷 +1 位作者 金峰 汪强 《岩土力学》 EI CAS CSCD 北大核心 2010年第12期3765-3770,共6页
针对振动台试验中无限地基难以模拟和数值分析中流-固耦合作用难以计算两个难题,将最近发展的实时耦联动力试验方法引入土-结构-流体动力相互作用问题的研究。以一个渡槽结构为例,其中渡槽-水体作为物理子结构,采用振动台进行物理试验,... 针对振动台试验中无限地基难以模拟和数值分析中流-固耦合作用难以计算两个难题,将最近发展的实时耦联动力试验方法引入土-结构-流体动力相互作用问题的研究。以一个渡槽结构为例,其中渡槽-水体作为物理子结构,采用振动台进行物理试验,而无限地基作为数值子结构,采用集总参数模型进行数值模拟。两个子结构之间实时交换数据,联合评估整个耦合体系的动力响应。试验结果和有限元数值模拟结果吻合良好,表明该试验方法具有较高精度。对不同特性地基土进行的试验对比分析结果表明:对于软土地基,考虑土-结构相互作用(SSI)的结构反应幅值明显减小,周期延长;随着地基土变硬,SSI效应逐渐变弱,结构反应最终收敛至刚性地基解。 展开更多
关键词 土-结构-流体动力相互作用 实时耦联动力试验 渡槽结构 集总参数模型
在线阅读 下载PDF
肾动脉变窄和流-固结构相互作用对非Newton脉动流的影响--一个真实的腹部主动脉和肾动脉模型 被引量:2
3
作者 Z·摩塔扎维尼亚 A·扎尔 +2 位作者 A·麦迪查德 海治(译) 张禄坤(校) 《应用数学和力学》 CSCD 北大核心 2012年第2期164-176,共13页
研究肾动脉狭窄(RAS)对血液流动和血管壁的影响.根据CT扫描图像,重建腹部主动脉和肾动脉的解剖模型,通过模型的脉动流进行了仿真计算,计算中考虑了流体-固体结构的相互作用(FSI).研究RAS对血管壁剪切应力和位移的影响,RAS使得肾动脉中... 研究肾动脉狭窄(RAS)对血液流动和血管壁的影响.根据CT扫描图像,重建腹部主动脉和肾动脉的解剖模型,通过模型的脉动流进行了仿真计算,计算中考虑了流体-固体结构的相互作用(FSI).研究RAS对血管壁剪切应力和位移的影响,RAS使得肾动脉中流量减少,肾素-血管紧缩素系统可能被激活,从而导致严重的高血压. 展开更多
关键词 肾动脉狭窄 脉动 流体-固体结构相互作用 非Newton流体 高血压
在线阅读 下载PDF
瓦斯/煤尘爆炸冲击波作用下结构的动力响应探讨 被引量:4
4
作者 袁芙蓉 《煤矿机械》 2017年第6期65-67,共3页
总结了国内外瓦斯/煤尘爆炸事故特点,研究了当前形势下瓦斯/煤尘爆炸作用下的灾害作用特性,揭示了瓦斯煤尘单一相、两相共混爆炸的结构动力响应机理及爆炸冲击波与结构间的耦合作用,对比分析了当前国内外爆炸结构相应问题的研究现状,探... 总结了国内外瓦斯/煤尘爆炸事故特点,研究了当前形势下瓦斯/煤尘爆炸作用下的灾害作用特性,揭示了瓦斯煤尘单一相、两相共混爆炸的结构动力响应机理及爆炸冲击波与结构间的耦合作用,对比分析了当前国内外爆炸结构相应问题的研究现状,探讨了瓦斯煤尘爆炸防控设备及装置动态安全性的研究方法,提出了瓦斯煤尘爆炸作用下结构动力问题的研究展望。研究工作对于煤矿瓦斯/煤尘爆炸致灾防控具有一定指导意义。 展开更多
关键词 瓦斯/煤尘爆炸 流体结构相互作用 动力响应 多场耦合作用 共混爆炸 致灾机理
在线阅读 下载PDF
简单直管直接水击的液固耦合作用分析 被引量:4
5
作者 吴昊 马跃先 马驭贞 《水电能源科学》 2004年第2期71-72,85,共3页
以能量原理为基础,在考虑水体可压缩性和管壁弹性变形对水击运动影响的基础上,对简单直管内发生的水击波动与管壁结构相互作用问题进行了探讨,得到了水击波与管壁相互作用的能量公式。依据该公式,可对在FSI效应下简单直管内的直接水击... 以能量原理为基础,在考虑水体可压缩性和管壁弹性变形对水击运动影响的基础上,对简单直管内发生的水击波动与管壁结构相互作用问题进行了探讨,得到了水击波与管壁相互作用的能量公式。依据该公式,可对在FSI效应下简单直管内的直接水击压强进行计算。 展开更多
关键词 流体结构相互作用 能量原理 动能 弹性势能 水击压强
在线阅读 下载PDF
结构振动对湍流近尾迹的影响 被引量:6
6
作者 刘明侯 Zhou Y. +1 位作者 陈义良 So R.M.C. 《力学学报》 EI CSCD 北大核心 2001年第2期163-172,共10页
研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响。用激光测振仪测量圆柱在升力方向的位移;用热线和LDA(二维)测量湍流的近尾迹。通过变化自由流的速度和圆柱体直径(特征尺寸)来变化雷诺数。用两个振动特... 研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响。用激光测振仪测量圆柱在升力方向的位移;用热线和LDA(二维)测量湍流的近尾迹。通过变化自由流的速度和圆柱体直径(特征尺寸)来变化雷诺数。用两个振动特性不同的(一个相对刚性,一个相对弹性)圆柱来产生尾迹,研究固体结构振动对湍流近尾迹的平均速度场和湍流场的影响。结果表明,结构自由振动对湍流近尾这场影响明显,该影响随雷诺数的变化不明显。 展开更多
关键词 流体诱发振动 湍流 尾迹结构 流体-结构相互作用 结构振动
在线阅读 下载PDF
水下冲击载荷作用下复合材料圆柱壳的弹性响应 被引量:1
7
作者 胡刚义 夏飞 《振动与冲击》 EI CSCD 北大核心 2010年第3期67-71,88,共6页
研究了平面阶跃冲击波作用下沉浸无限长复合材料圆柱壳的瞬态响应。首先推导基于Flugge薄壳理论的复合材料圆柱壳的控制运动方程,然后采用反射尾流虚源法建立流体-结构的相互作用,最后采用有限差分法来求解平面冲击波作用下无限长复合... 研究了平面阶跃冲击波作用下沉浸无限长复合材料圆柱壳的瞬态响应。首先推导基于Flugge薄壳理论的复合材料圆柱壳的控制运动方程,然后采用反射尾流虚源法建立流体-结构的相互作用,最后采用有限差分法来求解平面冲击波作用下无限长复合材料圆柱壳的控制运动微分方程。详细考察了复合材料纤维方向和壳的半径对复合材料圆柱壳的无因次中面应变、第0阶模态径向位移和第1阶模态径向速度的影响。 展开更多
关键词 复合材料圆柱壳 水下冲击 反射尾流虚源法 流体-结构相互作用
在线阅读 下载PDF
水下悬浮隧道的技术现状和未来展望
8
作者 Marcel‘t Hart Arianna Minoretti 李英 《隧道建设(中英文)》 CSCD 北大核心 2024年第2期225-232,I0025-I0032,共16页
阐述了水下悬浮隧道(又称水下悬浮管桥)的概念;探讨了悬浮隧道结构的4种支撑形式,即岸锚式、桥墩支撑式、浮筒固定式和锚索固定式,并对4种支撑形式的适用性进行了对比;对悬浮隧道的跨海形式及主隧道结构断面形式进行了讨论,分析了多管系... 阐述了水下悬浮隧道(又称水下悬浮管桥)的概念;探讨了悬浮隧道结构的4种支撑形式,即岸锚式、桥墩支撑式、浮筒固定式和锚索固定式,并对4种支撑形式的适用性进行了对比;对悬浮隧道的跨海形式及主隧道结构断面形式进行了讨论,分析了多管系统+横通道连通形式在火灾及事故发生时的安全性。介绍了国际上悬浮隧道的相关标准:国际结构混凝土联合会(fib)在2020年发布了《水下悬浮管桥(SFTB)设计指南》(fib公告第96号),旨在为技术界提供水下悬浮管桥方案设计指南;国际隧道和地下空间协会(ITA)在2023年发布了《水下悬浮隧道(SFT)业主指南》,旨在向基础设施业主说明,水下悬浮隧道相比传统隧道的适用性;另外,挪威国家标准已率先将悬浮隧道结构纳入国家桥梁手册(2020)。梳理了悬浮隧道已有的理论研究,介绍了目前水下悬浮隧道比较热门的研究课题,比如流体-结构相互作用和冲击等。分析了水下悬浮隧道的优势,不仅能够减小主要的海洋荷载对结构的影响,而且可以降低隧道和相连结构的坡度,缩短穿越距离,保护海底地形以及减少对周围区域的噪声影响。随着人们对环境问题的日益重视,以及新材料和新技术为减少气候影响带来的机遇,悬浮隧道有望获得进一步的推广。 展开更多
关键词 水下悬浮隧道 水下悬浮管桥 指南 流体-结构相互作用 荷载 冲击 支撑形式
在线阅读 下载PDF
四角架筒型基础气浮过程中运动特性研究 被引量:1
9
作者 刘宪庆 赵明阶 +2 位作者 张浦阳 罗盛 钟韡 《太阳能学报》 EI CAS CSCD 北大核心 2022年第4期428-433,共6页
以四角架筒型基础为研究对象,通过在静水及规则波作用下的模型试验对影响结构气浮过程中运动特性的因素进行分析。结果表明,结构摇荡运动的附加质量系数随吃水的增加呈下降趋势,都大于船舶动力学的建议值1.2;随着吃水的增加,垂荡运动呈... 以四角架筒型基础为研究对象,通过在静水及规则波作用下的模型试验对影响结构气浮过程中运动特性的因素进行分析。结果表明,结构摇荡运动的附加质量系数随吃水的增加呈下降趋势,都大于船舶动力学的建议值1.2;随着吃水的增加,垂荡运动呈先减小后增大的趋势,而纵摇运动呈减小的趋势;较浅水深增强了结构的竖向运动,垂荡和纵摇运动的最大幅值随水深的变化呈相反的变化趋势;四角架筒型基础结构在4.5 m吃水和12.5 m水深下的运动响应特性最佳。 展开更多
关键词 海上风力机 传递函数 流体结构相互作用 筒型基础 附加质量
在线阅读 下载PDF
船舶碰撞动力学过程的数值仿真研究 被引量:104
10
作者 王自力 顾永宁 《爆炸与冲击》 EI CAS CSCD 北大核心 2001年第1期29-34,共6页
在分析显式非线性有限元基本理论和关键技术的基础上 ,探讨了船舶动力学过程的数值仿真方法 ,并对碰撞过程中船体的力学性能进行了分析 ,获得并讨论了碰撞力、能量吸收和损伤变形的时序结果 ,所得结论具有一般性。碰撞仿真分析中 ,被撞... 在分析显式非线性有限元基本理论和关键技术的基础上 ,探讨了船舶动力学过程的数值仿真方法 ,并对碰撞过程中船体的力学性能进行了分析 ,获得并讨论了碰撞力、能量吸收和损伤变形的时序结果 ,所得结论具有一般性。碰撞仿真分析中 ,被撞船的舷侧受撞区域和撞击船的首部区域作为可变形结构处理 ,而其余区域则视为刚体。被撞船周围的水划分成Euler有限体积网格 ,船体结构则划分成Lagrange有限元网格 ,并采用耦合技术将它们联系在一起。 展开更多
关键词 船舶碰撞 数值仿真 结构耐撞性 流体-结构相互作用 动力学过程
在线阅读 下载PDF
开阔海洋环境中的鲍鱼潜式网箱系统的工程分析
11
《渔业现代化》 北大核心 2014年第6期72-73,共2页
在还未被利用的开阔海洋环境中开发鲍鱼养殖存在极大的可能性.鲍鱼是一种摄食巨藻的贝类.本研究对具有组装式框架结构的一个商业规模的鲍鱼潜式养殖网箱养成系统的开发进行了描述.使用计算流体力学软件分析了围养鲍鱼的结构物内的流场特... 在还未被利用的开阔海洋环境中开发鲍鱼养殖存在极大的可能性.鲍鱼是一种摄食巨藻的贝类.本研究对具有组装式框架结构的一个商业规模的鲍鱼潜式养殖网箱养成系统的开发进行了描述.使用计算流体力学软件分析了围养鲍鱼的结构物内的流场特性.使用一个模拟流体结构相互作用的Morison方程类的有限元模型,调研了锚碇的围闭结构的水动力响应. 展开更多
关键词 鲍鱼养殖 海洋环境 养殖网箱 系统 流体结构相互作用 Morison方程 工程 计算流体力学
在线阅读 下载PDF
浮式跨海大桥技术的研究进展 被引量:1
12
作者 陈文创 张永良 郭鹏 《黑龙江大学工程学报》 2019年第3期1-7,14,共8页
连岛跨海大桥是实现基础设施互联互通的关键工程之一,特别在宽广深邃的峡湾或海峡上建设跨海大桥将面临着巨大的挑战,采用浮式桥梁结构使得横跨这种峡湾或海峡成为可能。回顾浮式跨海大桥技术发展,总结了两种不同结构型式(离散型和连续... 连岛跨海大桥是实现基础设施互联互通的关键工程之一,特别在宽广深邃的峡湾或海峡上建设跨海大桥将面临着巨大的挑战,采用浮式桥梁结构使得横跨这种峡湾或海峡成为可能。回顾浮式跨海大桥技术发展,总结了两种不同结构型式(离散型和连续型浮筒)的工作性能及两种支承方式(系泊于海底和端部支承)的工作特性,阐明现有浮桥模拟理论的局限性。提出需要建立黏性流体—多结构体—节间阻尼系统耦合的数学模型,探索流体—多结构体相互作用机理;提出了浮动结构的分析灵敏度矩阵中需要充分考虑黏性流体和结构的动力相互作用,探索多模态下模型的敏感性特征,建立一个充分考虑黏性流体和结构动力相互作用的浮桥结构的分析灵敏度矩阵,提高数值结果的置信度。最后提出了浮式跨海大桥技术研究中所要注意的一些问题。 展开更多
关键词 浮式跨海大桥 流体-结构相互作用 回顾 波浪
在线阅读 下载PDF
Dynamic Adaptive Finite Element Analysis of Acoustic Wave Propagation Due to Underwater Explosion for Fluid-structure Interaction Problems 被引量:4
13
作者 Seyed Shahab Emamzadeh Mohammad Taghi Ahmadi +1 位作者 Soheil Mohammadi Masoud Biglarkhani 《Journal of Marine Science and Application》 CSCD 2015年第3期302-315,共14页
In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the s... In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes:a) calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b) fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c) adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved. 展开更多
关键词 adaptive mesh fluid-structure interaction acoustic wave finite element analysis underwater explosion
在线阅读 下载PDF
A CFD Based Investigation of the Unsteady Hydrodynamic Coefficients of 3-D Fins in Viscous Flow 被引量:3
14
作者 Nazir Zulfiqar 苏玉民 王兆立 《Journal of Marine Science and Application》 2010年第3期250-255,共6页
The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly ... The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly unsteady, generating vortices and requiring detailed analysis of fluid-structure interactions.An understanding of the complexities of such flows is of interest to engineers developing vehicles capable of high dynamic performance in their propulsion and maneuvering.In the present study, a CFD based RANS simulation of a 3-D fin body moving in a viscous fluid was developed.It investigated hydrodynamic performance by evaluating the hydrodynamic coefficients (lift, drag and moment) at two different oscillating frequencies.A parametric analysis of the factors that affect the hydrodynamic performance of the fin body was done, along with a comparison of results from experiments.The results of the simulation were found in close agreement with experimental results and this validated the simulation as an effective tool for evaluation of the unsteady hydrodynamic coefficients of 3-D fins.This work can be further be used for analysis of the stability and maneuverability of fin actuated underwater vehicles. 展开更多
关键词 oscillating 3-D fin RANS hydrodynamic performance viscous flow
在线阅读 下载PDF
Unstructured Grid Immersed Boundary Method for Numerical Simulation of Fluid Structure Interaction 被引量:2
15
作者 明平剑 孙扬哲 +1 位作者 段文洋 张文平 《Journal of Marine Science and Application》 2010年第2期181-186,共6页
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ... This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder. 展开更多
关键词 fluid structure interaction immersed boundary method VOS unstructured grids finite volume method
在线阅读 下载PDF
不同质量粒子分布SPH方法及其应用 被引量:2
16
作者 艾孜海尔.哈力克 热合买提江.依明 +1 位作者 开依沙尔.热合曼 买买提明.艾尼 《振动与冲击》 EI CSCD 北大核心 2015年第22期62-67,共6页
光滑粒子流体动力学(SPH)方法是一种新发展起来的纯拉格朗日数值方法,其最大的优点是不需要对自由表面进行特殊处理,有利于模拟飞溅、破碎等高度非线性流。SPH方法的一个缺点是计算的时间复杂性高,因为需要大量粒子详细定义流体,使之计... 光滑粒子流体动力学(SPH)方法是一种新发展起来的纯拉格朗日数值方法,其最大的优点是不需要对自由表面进行特殊处理,有利于模拟飞溅、破碎等高度非线性流。SPH方法的一个缺点是计算的时间复杂性高,因为需要大量粒子详细定义流体,使之计算时间相应增加。在核半径保持不变的情况下为了提高计算效率和精度提出不同质量粒子分布的方法,即在关键计算区域分布较小质量粒子,而在其它区域分布较大质量粒子。首先通过静水数值模拟,初步验证了此方法的有效性。之后,对半浸泡振荡圆筒所产生的小振幅水波进行了数值模拟。将数值模拟结果与已有的实验结果进行了对比分析。分析结果表明此方法能够明显提高计算精度的同时可以节省大量计算时间。 展开更多
关键词 光滑流体动力学方法 波浪 振荡圆筒 不同质量粒子 流体结构相互作用
在线阅读 下载PDF
Construction of Wave-free Potentials and Multipoles in a Two-layer Fluid Having Free-surface Boundary Condition with Higher-order Derivatives 被引量:1
17
作者 Dilip Das 《Journal of Marine Science and Application》 CSCD 2015年第3期270-282,共13页
There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle th... There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace’s equation (or Helmholtz equation) arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace’s equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential. 展开更多
关键词 two-layer fluid wave-free potentials Laplace’s equation modified Helmholtz equations higher order boundary conditions MULTIPOLES
在线阅读 下载PDF
Two-dimensional Numerical Simulation of an Elastic Wedge Water Entry by a Coupled FDM-FEM Method
18
作者 Kangping Liao Changhong Hu Wenyang Duan 《Journal of Marine Science and Application》 2013年第2期163-169,共7页
Hydroelastic behavior of an elastic wedge impacting on calm water surface was investigated. A partitioned approach by coupling finite difference method (FDM) and finite element method (FEM) was developed to analyz... Hydroelastic behavior of an elastic wedge impacting on calm water surface was investigated. A partitioned approach by coupling finite difference method (FDM) and finite element method (FEM) was developed to analyze the fluid structure interaction (FSI) problem. The FDM, in which the Constraint Interpolation Profile (CIP) method was applied, was used for solving the flow field in a fixed regular Cartesian grid system. Free surface was captured by the Tangent of Hyperbola for Interface Capturing with Slope Weighting (THINC/SW) scheme. The FEM was applied for calculating the structural deformation. A volume weighted method, which was based on the immersed boundary (IB) method, was adopted for coupling the FDM and the FEM together. An elastic wedge water entry problem was calculated by the coupled FDM-FEM method. Also a comparison between the current numerical results and the published results indicate that the coupled FDM-FEM method has reasonably good accuracy in predicting the impact force. 展开更多
关键词 elastic wedge water entry coupled FDM-FEM method volume weighted method CIP method THINC/SW scheme hydroelastic behavior
在线阅读 下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部