期刊文献+
共找到114篇文章
< 1 2 6 >
每页显示 20 50 100
基于精确色散Boussinesq方程的珊瑚礁波浪传播模拟
1
作者 刘维杰 宁越 《长沙理工大学学报(自然科学版)》 CAS 2023年第4期47-54,共8页
【目的】提高激波捕捉类Boussinesq模型预测珊瑚礁地形上次重力波空间分布的准确性。【方法】基于具有精确色散性能的Boussinesq方程构建混合数值格式的数值模型,并对二维以及三维珊瑚礁地形上的波浪传播进行了数值模拟。【结果】相比... 【目的】提高激波捕捉类Boussinesq模型预测珊瑚礁地形上次重力波空间分布的准确性。【方法】基于具有精确色散性能的Boussinesq方程构建混合数值格式的数值模型,并对二维以及三维珊瑚礁地形上的波浪传播进行了数值模拟。【结果】相比目前较为流行的弱色散激波捕捉类模型,本文模型可更为准确地模拟二维珊瑚礁地形上的波能转换特征以及相应的次重力波空间分布,同时,本文模型对于三维珊瑚礁-裂口-潟湖地形上的波浪传播问题也具有较好的适应性。但由于控制方程的弱非线性,本文模型对于礁前陡坡上波浪浅化和潟湖上波高的预测仍有一定差异,且方程中的卷积项也会降低模型的计算效率。【结论】控制方程的色散性能应是限制激波捕捉类Boussinesq模型模拟珊瑚礁次重力波准确性的主要因素,如何在保持方程色散性的情况下兼顾模型的非线性性能和计算效率应是后续激波捕捉类Boussinesq模型在珊瑚礁地形上改进的主要方向。 展开更多
关键词 精确色散 BOUSSINESQ方程 珊瑚礁地形 次重力波 激波捕捉 波浪传播
下载PDF
珊瑚礁地形中波浪传播变形及增水模拟
2
作者 吕幸 季新然 《海南大学学报(自然科学版)》 CAS 2023年第3期303-311,共9页
本文基于VRANS方程建立数值模型,将珊瑚礁表层简化为具有恒定孔隙率的多孔介质,模拟分析了在带有潟湖的珊瑚礁地形中不同入射波浪参数对礁坪上及潟湖中波高分布和增水效应的影响.结果表明,规则波浪在珊瑚礁地形上发生破碎后,在礁坪上波... 本文基于VRANS方程建立数值模型,将珊瑚礁表层简化为具有恒定孔隙率的多孔介质,模拟分析了在带有潟湖的珊瑚礁地形中不同入射波浪参数对礁坪上及潟湖中波高分布和增水效应的影响.结果表明,规则波浪在珊瑚礁地形上发生破碎后,在礁坪上波高明显降低;同时,在珊瑚礁多孔介质的作用下,增水效应最明显的位置位于礁坪上靠近礁缘处,入射波高和周期越大,在礁坪及潟湖中的波高分布和增水越大.研究结果对多孔礁体上的波浪变形和增水规律作出了解释,同时也为分析礁坪上构筑物出现失效和破坏的原因提供了参考. 展开更多
关键词 多孔礁体 波浪传播 礁坪增水 潟湖
下载PDF
岸礁礁冠对波浪传播变形及增水影响的实验研究 被引量:19
3
作者 姚宇 袁万成 +1 位作者 杜睿超 蒋昌波 《热带海洋学报》 CAS CSCD 北大核心 2015年第6期19-25,共7页
通过在波浪槽中进行一系列物理模型实验来研究珊瑚岸礁礁冠的存在对波浪传播变形和礁坪上增水的影响,实验在现场观测的基础上采用了理想化的岸礁和礁冠模型,并测试了4种礁坪水深和不同的入射波工况。实验结果表明,礁冠能显著改变波浪在... 通过在波浪槽中进行一系列物理模型实验来研究珊瑚岸礁礁冠的存在对波浪传播变形和礁坪上增水的影响,实验在现场观测的基础上采用了理想化的岸礁和礁冠模型,并测试了4种礁坪水深和不同的入射波工况。实验结果表明,礁冠能显著改变波浪在礁坪边缘的传播变形过程,尤其是波浪破碎强度和高频波的产生;礁坪的增水随着入射波高和周期增加而增加,随着水深的增加而减少;礁冠的阻水效应类似于潜堤,对于其后礁坪上增水的大小起决定作用,礁冠存在时波浪增水显著增大;随后引入了一个新的无量纲参数来考虑礁冠的存在对增水的影响。 展开更多
关键词 波浪传播变形 波浪增水 岸礁礁冠 物理模型实验 无量纲参数
下载PDF
潮汐流影响下珊瑚岛礁附近波浪传播变形和增水试验 被引量:15
4
作者 姚宇 唐政江 +1 位作者 杜睿超 蒋昌波 《水科学进展》 EI CAS CSCD 北大核心 2017年第4期614-621,共8页
为揭示潮汐流影响下珊瑚岛礁附近波浪的运动特征,通过波流水槽试验对潮汐流存在下规则波的传播变形和增水规律进行了研究,测试了一系列的潮汐流流量并对比正向流、反向流以及无流的情况。研究结果表明:反向流使波浪破碎点向远海侧移动,... 为揭示潮汐流影响下珊瑚岛礁附近波浪的运动特征,通过波流水槽试验对潮汐流存在下规则波的传播变形和增水规律进行了研究,测试了一系列的潮汐流流量并对比正向流、反向流以及无流的情况。研究结果表明:反向流使波浪破碎点向远海侧移动,正向流则使其向海岸侧移动,反向流存在时破碎带内的湍流较正向流时更为剧烈;反向流时破碎点附近入射波能量由主频波向高次谐波发生转移更为显著;潮汐流的存在对波浪反射、透射和能量耗散并无明显影响;正向流造成礁坪上波浪增水减少,反向流则促进增水的增长,增水最大值与潮汐流流量间存在显著的线性关系。研究成果可望为岛礁工程的建设和维护提供一定的理论参考。 展开更多
关键词 波浪传播变形 增水 规则波 潮汐流 珊瑚岛礁
下载PDF
一般曲线坐标系下波浪传播的数值模拟 被引量:13
5
作者 张洪生 丁平兴 赵海虹 《海洋学报》 CAS CSCD 北大核心 2003年第1期110-119,共10页
在曲线坐标系下,建立了缓变水深水域波浪传播的数值模拟模型.模型适宜于复杂变化的边界形状,克服了各种代数坐标变换的局限性.在建立模型时,将原始的椭圆型缓坡方程的近似型式——依赖时间变化的抛物型方程,作为控制方程,既克服了一般... 在曲线坐标系下,建立了缓变水深水域波浪传播的数值模拟模型.模型适宜于复杂变化的边界形状,克服了各种代数坐标变换的局限性.在建立模型时,将原始的椭圆型缓坡方程的近似型式——依赖时间变化的抛物型方程,作为控制方程,既克服了一般抛物近似方法的缺点,又便利了方程的求解;从开边界条件、不同反射特性的固壁边界条件相统一的表达式出发,对边界条件进行处理;用ADI法数值求解控制方程.对模型的验证表明,数值解与物模实验值吻合良好,模型对于具有复杂边界的工程实际有较强的适应性. 展开更多
关键词 波浪传播 曲线坐标系 数值模拟 边界条件 ADI法 自适应网格 边界条件
下载PDF
海岸河口水域波浪传播数值模拟 被引量:15
6
作者 洪广文 冯卫兵 张洪生 《河海大学学报(自然科学版)》 CAS CSCD 1999年第2期1-9,共9页
提出了水深与流场缓变水域波浪传播数学模型———水流中依赖时间变量并考虑能耗的波浪“缓坡方程”及其等价的控制方程组,分析比较了无水流情况此理论模型与其相应的两种抛物型近似的差别,提供了长江口波浪变形数值模拟计算工程实例... 提出了水深与流场缓变水域波浪传播数学模型———水流中依赖时间变量并考虑能耗的波浪“缓坡方程”及其等价的控制方程组,分析比较了无水流情况此理论模型与其相应的两种抛物型近似的差别,提供了长江口波浪变形数值模拟计算工程实例.实例表明,该模型能适应河口三角洲大范围水域波浪传播数值计算. 展开更多
关键词 波浪绕射-折射 缓坡方程 波浪传播 数学模型
下载PDF
深水航道对波浪传播影响规律的研究 被引量:15
7
作者 赵智帮 刘子琪 +1 位作者 孙精石 牛恩宗 《中国港湾建设》 北大核心 2001年第6期9-15,共7页
航道对波浪的作用已经成为深水港建设的关键技术 ,本文就航道对波浪传播的影响规律进行了广泛的探讨。根据物理模型和数学模型的试验结果 ,分别阐述了航道水深、航道边坡坡度、航道宽度、海底坡度以及波浪入射角对波浪传播变形的影响。... 航道对波浪的作用已经成为深水港建设的关键技术 ,本文就航道对波浪传播的影响规律进行了广泛的探讨。根据物理模型和数学模型的试验结果 ,分别阐述了航道水深、航道边坡坡度、航道宽度、海底坡度以及波浪入射角对波浪传播变形的影响。文中还就临界入射角在工程中的意义、截短航道进行模型试验的可行性以及最大波高与水深比等进行了分析与论证。研究成果对港口及航道工程的规划设计具有指导意义。 展开更多
关键词 深水航道 波浪折射 波浪传播 波能聚集 防波堤 数学模型 物理模型 防浪设计
下载PDF
珊瑚礁地形上波浪传播变形的非静压数值模拟 被引量:4
8
作者 张其一 史宏达 +1 位作者 高伟 李金峰 《海岸工程》 2017年第4期1-9,共9页
珊瑚礁广泛分布于我国南海等海域,具有重要的生态、旅游、港口和军事价值。波浪是珊瑚礁环境最重要的动力因素之一。基于非静压浅水方程波浪模型SWASH,模拟了典型岸礁地形上的波浪传播过程,并采用Demirbilek等的水槽实验数据分别对有效... 珊瑚礁广泛分布于我国南海等海域,具有重要的生态、旅游、港口和军事价值。波浪是珊瑚礁环境最重要的动力因素之一。基于非静压浅水方程波浪模型SWASH,模拟了典型岸礁地形上的波浪传播过程,并采用Demirbilek等的水槽实验数据分别对有效波高、波浪增水与爬高进行了验证,模拟结果与实测数据符合良好,展示了非静压模型在珊瑚礁地形上良好的适用性。通过波浪频率谱和波面小波谱,分析了波浪谱在该珊瑚礁地形上的演化规律,并揭示了礁坪上波浪能量从高频向低频转移的现象。本研究结果可为珊瑚礁环境动力研究以及相关的工程建设提供有益的参考依据。 展开更多
关键词 珊瑚礁 波浪传播 非静压模型 SWASH 小波分析
下载PDF
航槽对不同入射角波浪传播的影响 被引量:14
9
作者 左其华 杨正己 《水运工程》 北大核心 1996年第1期1-8,共8页
采用数值计算方法探讨了不同入射角的波浪在典型航槽中传播时波幅的变化。计算结果表明航槽对波浪传播的影响是显著的。当入射波向与航槽轴线平行时.航槽内的波幅比随着传播距离的增长而衰减。入射角为15°~30°时.航槽... 采用数值计算方法探讨了不同入射角的波浪在典型航槽中传播时波幅的变化。计算结果表明航槽对波浪传播的影响是显著的。当入射波向与航槽轴线平行时.航槽内的波幅比随着传播距离的增长而衰减。入射角为15°~30°时.航槽迎浪一侧边坡及其附近滩地出现较大波高比区域、而另一测的波高比往往较小.航槽起到了一定的“遮掩”作用。入射角大于45°时、航槽对波浪传播的影响逐渐减弱,波浪正向穿越航槽时.波幅比可用浅水系数来计算. 展开更多
关键词 航道 航槽 入射角 波浪传播
下载PDF
基于改进缓坡方程的波浪传播数值模拟 被引量:6
10
作者 王红川 周正萍 《海洋工程》 CSCD 北大核心 2013年第3期45-53,共9页
用变分原理导出考虑底坡一阶导数平方项和二阶曲率项影响的缓坡方程,对传统缓坡方程作了改进,提高波浪在海底地形变化剧烈、水深较浅时数值模拟精度。数值计算与已有实验室试验资料比较表明,该模型可以较好地模拟有剧烈变化的海底地形... 用变分原理导出考虑底坡一阶导数平方项和二阶曲率项影响的缓坡方程,对传统缓坡方程作了改进,提高波浪在海底地形变化剧烈、水深较浅时数值模拟精度。数值计算与已有实验室试验资料比较表明,该模型可以较好地模拟有剧烈变化的海底地形的波浪传播,比传统缓坡方程模型计算结果在精度上有明显提高。 展开更多
关键词 变分原理 缓坡方程 数值计算 波浪传播
下载PDF
应用MIKE 21 BW模型分析航道对波浪传播的影响 被引量:12
11
作者 徐俊锋 方春明 +1 位作者 黄烈敏 凃洋 《中国水利水电科学研究院学报》 2011年第4期292-297,共6页
波浪在传播过程中遇到航道会发生折射、绕射以及反射现象。MIKE 21的BW模块可以很好地计算和分析港口、码头等小面积波浪场中的多方向不规则波的相互作用,浅水变形,折射,绕射和反射等各种波浪现象。本文采用该模块建立在不同航道条件下... 波浪在传播过程中遇到航道会发生折射、绕射以及反射现象。MIKE 21的BW模块可以很好地计算和分析港口、码头等小面积波浪场中的多方向不规则波的相互作用,浅水变形,折射,绕射和反射等各种波浪现象。本文采用该模块建立在不同航道条件下的波浪传播数学模型,并利用数模结果系统地分析了航道水深、航道宽度、航道边坡坡度以及波浪入射角度对波浪传播与变形的影响规律。成果对试验研究及工程实践有指导意义。 展开更多
关键词 航道 数值模拟 MIKE 21 波浪传播
下载PDF
应用Boussinesq方程对由航道开挖所造成水深变化对波浪传播所产生的异常波况的数值研究 被引量:10
12
作者 张永刚 李玉成 《海岸工程》 1997年第1期7-17,共11页
应用Boussinesq方程采用有限差分法建立了非线性数值波浪模式,并应用该模式对由航道开挖所造成水深变化对波浪传播产生的异常现象进行了数值模拟研究。从结果来看这种异常波浪局部增大现象是由于入射波与航道夹角过小,使波浪无法折射... 应用Boussinesq方程采用有限差分法建立了非线性数值波浪模式,并应用该模式对由航道开挖所造成水深变化对波浪传播产生的异常现象进行了数值模拟研究。从结果来看这种异常波浪局部增大现象是由于入射波与航道夹角过小,使波浪无法折射入航道,而在航道两侧反射叠加的结果。因此说水深的变化对波浪产生的反射现象也是不可忽视的。 展开更多
关键词 BOUSSINESQ方程 航道开挖 波浪传播 异常波况
下载PDF
曲线坐标系下非均匀水流中波浪传播的数值模拟 被引量:2
13
作者 张洪生 丁平兴 +1 位作者 吴中 缪国平 《上海交通大学学报》 EI CAS CSCD 北大核心 2003年第8期1261-1265,共5页
利用自适应网格技术 ,建立了曲线坐标系下缓变水深水域非均匀水流中波浪传播的数值模拟模型 .将椭圆型非均匀水流中波浪传播的联合折射 -绕射方程化为依赖时间变化的抛物型方程 ,并以此作为控制方程 ,克服了一般抛物近似方法的缺点 ;从... 利用自适应网格技术 ,建立了曲线坐标系下缓变水深水域非均匀水流中波浪传播的数值模拟模型 .将椭圆型非均匀水流中波浪传播的联合折射 -绕射方程化为依赖时间变化的抛物型方程 ,并以此作为控制方程 ,克服了一般抛物近似方法的缺点 ;从开边界条件、不同反射特性的固壁边界条件相统一的表达式出发 ,对边界条件进行处理 ;用 ADI法数值求解控制方程 .验证结果表明 ,模型能适应复杂的边界形状 ,可以较好地模拟复杂边界条件下波浪传播过程中的多种物理现象 ,克服了以往各种代数坐标变换的局限性 ,有效地反映了水流对波浪传播的影响 .模型的模拟结果可供实际工程参考使用 . 展开更多
关键词 数值模拟 曲线坐标系 非均匀水流 波浪传播
下载PDF
曲线坐标系下波浪传播的数学模型及其比较与验证 被引量:2
14
作者 张洪生 尤云祥 +1 位作者 朱良生 张军 《水动力学研究与进展(A辑)》 CSCD 北大核心 2005年第1期106-117,共12页
引入自适应网格技术,建立了曲线坐标系下考虑能耗影响的缓变水深水域波浪传播的数值模拟模型,通过在模型中引入能量耗散项以更好地处理全透射边界。基于原型缓坡方程,完善了等水深的轴对称双导堤水域内波浪传播的精确解表达式及其计算结... 引入自适应网格技术,建立了曲线坐标系下考虑能耗影响的缓变水深水域波浪传播的数值模拟模型,通过在模型中引入能量耗散项以更好地处理全透射边界。基于原型缓坡方程,完善了等水深的轴对称双导堤水域内波浪传播的精确解表达式及其计算结果;提供了等水深的环形河道水域内较为详细的精确解。在极坐标系下就等水深的轴对称双导堤水域内波浪的传播,推导了缓坡方程的小角度抛物近似模型的解析解表达式。就上述两水域内波浪传播的数值模拟模型的数值解、抛物近似模型的解、精确解进行了详细比较,说明数值解与精确解相吻合、数值模拟模型比小角度和大角度抛物近似模型具有更高的精度。数值模拟模型能较好地模拟水深复杂变化的地形上波浪的传播,能有效反映波浪传播过程中的多种物理现象。 展开更多
关键词 波浪传播 数学模型 曲线坐标系 能量耗散项 解析解
下载PDF
透空式防波堤周围的非线性波浪传播的数值模拟 被引量:4
15
作者 李熙 王义刚 《海洋工程》 CSCD 北大核心 2004年第3期97-101,共5页
对非线形波浪在透空式防波堤周围的波浪变形进行了数值模拟,在Boussinesq波浪方程中加入与透空建筑物有关的新的耗散项,从而界定了透空建筑物引起的部分反射和透射,波浪折射衍射的传播过程通过控制方程求解。波浪控制方程通过有限差分... 对非线形波浪在透空式防波堤周围的波浪变形进行了数值模拟,在Boussinesq波浪方程中加入与透空建筑物有关的新的耗散项,从而界定了透空建筑物引起的部分反射和透射,波浪折射衍射的传播过程通过控制方程求解。波浪控制方程通过有限差分方法求解。模型应用于模拟波浪经过具有部分反射的群桩式透空结构,结果表明透空式防波堤可以有效地衰减波浪,是重力式结构的一种替代形式。 展开更多
关键词 透空式防波 BOUSSINESQ方程 波浪传播 部分反射和透射
下载PDF
基于GPU加速的Boussinesq类波浪传播变形数值模型 被引量:2
16
作者 孙家文 朱桐 +1 位作者 房克照 刘忠波 《海洋工程》 CSCD 北大核心 2020年第2期111-119,共9页
Boussinesq波浪模型是一类相位解析模型,在时域内求解需要较高的空间和时间分辨率以保证计算精度。为提高计算效率,有必要针对该类模型开展并行算法的研究。与传统的中央处理器(CPU)相比,图形处理器(GPU)有大量的运算器,可显著提高计算... Boussinesq波浪模型是一类相位解析模型,在时域内求解需要较高的空间和时间分辨率以保证计算精度。为提高计算效率,有必要针对该类模型开展并行算法的研究。与传统的中央处理器(CPU)相比,图形处理器(GPU)有大量的运算器,可显著提高计算效率。基于统一计算设备架构CUDA C语言和图形处理器,实现了Boussinesq模型的并行运算。将本模型的计算结果同CPU数值模拟结果和解析解相比较,发现得到的结果基本一致。同时也比较了CPU端与GPU端的计算效率,结果表明,GPU数值模型的计算效率有明显提升,并且伴随数值网格的增多,提升效果更为明显。 展开更多
关键词 BOUSSINESQ方程 图形处理器 CUDA C 并行计算 计算效率 波浪传播
下载PDF
潜堤上波浪传播过程数值研究 被引量:2
17
作者 廖斌 陈善群 任超洋 《长江科学院院报》 CSCD 北大核心 2017年第8期52-58,共7页
以RANS方程结合Level-Set方法为研究基础,采用五阶WENO有限差分格式进行空间离散、三阶TVD Runge-Kutta格式进行时间推进、Level-Set方法追踪波浪与空气间自由面以及解析松弛法来实现数值水槽中的造波消波,建立了一种求解潜堤上波浪传... 以RANS方程结合Level-Set方法为研究基础,采用五阶WENO有限差分格式进行空间离散、三阶TVD Runge-Kutta格式进行时间推进、Level-Set方法追踪波浪与空气间自由面以及解析松弛法来实现数值水槽中的造波消波,建立了一种求解潜堤上波浪传播问题的数值计算模型。选取经典的潜堤上波浪传播物理试验模型对数值模型进行了验证,水位计算结果与试验值吻合较好。进一步研究了波高、潜堤顶部淹没深度以及潜堤坡度等参数对潜堤上波浪传播过程的影响。结果表明:波高越高、淹没深度越浅、潜堤向波坡度越小,波浪受浅水作用越明显;潜堤背波坡度越小,波浪受反浅水作用稍大,但并不明显。 展开更多
关键词 潜堤 Level-Set方法 波浪传播 数值研究 浅水作用
下载PDF
波流共存场中多向随机波浪传播变形数学模型 被引量:9
18
作者 郑金海 Hajime Mase 《水科学进展》 EI CAS CSCD 北大核心 2008年第1期78-83,共6页
基于波作用量守恒方程建立了波流共存场中多向随机波浪传播变形数学模型,模型中考虑了波浪绕射的影响和水流引起的波浪弥散多普勒效应,应用包含水流和地形影响的激破波模式计算波浪破碎的能量耗散,采用一阶上迎风有限差分格式离散控制... 基于波作用量守恒方程建立了波流共存场中多向随机波浪传播变形数学模型,模型中考虑了波浪绕射的影响和水流引起的波浪弥散多普勒效应,应用包含水流和地形影响的激破波模式计算波浪破碎的能量耗散,采用一阶上迎风有限差分格式离散控制方程。分别计算了有无近岸流情况下单向和多向随机波浪的波高分布,考虑水流影响的数值计算结果与物理模型实验数据吻合良好,比较分析表明,所建立的数学模型能够复演由于离岸流引起的波高增大,可用于波流共存场多向随机波浪传播变形的模拟和预报。 展开更多
关键词 波流共存场 波作用量守恒方程 多向随机波浪 数学模型 波浪传播变形
下载PDF
大范围波浪传播数学模型 被引量:2
19
作者 杨耀中 冯卫兵 《水道港口》 2010年第2期83-89,共7页
南黄海辐射沙脊群海域范围广阔、海底地形变化剧烈。在传统缓坡方程基础上建立的波浪传播数学模型难以计算该海域的波浪场。为了解决此问题,文章在洪广文推导的缓坡方程推广形式的基础上,改进了数值计算方法,建立了大范围波浪传播数学... 南黄海辐射沙脊群海域范围广阔、海底地形变化剧烈。在传统缓坡方程基础上建立的波浪传播数学模型难以计算该海域的波浪场。为了解决此问题,文章在洪广文推导的缓坡方程推广形式的基础上,改进了数值计算方法,建立了大范围波浪传播数学模型。运用波数矢无旋性方程,将缓坡方程转化为波作用守恒方程和光程函数方程,从而消除了波长对计算步长的限制。运用硬盘同步记录法,把求解时所需要的内存用量转化为硬盘用量,从而可以运用高分辨率网格计算大范围海域。在南黄海辐射沙脊群的应用表明,该模型适用于大范围、复杂水下地形的波浪场计算。 展开更多
关键词 波浪传播 数学模型 辐射沙脊群 南黄海
下载PDF
小迎浪角时长直航道对波浪传播与变形影响的数值研究 被引量:11
20
作者 吴达开 《水运工程》 北大核心 1997年第1期3-13,共11页
首先将新近研制的透浪侧界抛物线模型,结合方向谱计算,应用于航道对波浪影响的研究。其次阐明了小迎浪角时的航道散射机理;提出了航道影响线的概念。最后分析了长、直航道沿程比波高分布的趋势。得到了航道口门横断面的波能偏态盆形... 首先将新近研制的透浪侧界抛物线模型,结合方向谱计算,应用于航道对波浪影响的研究。其次阐明了小迎浪角时的航道散射机理;提出了航道影响线的概念。最后分析了长、直航道沿程比波高分布的趋势。得到了航道口门横断面的波能偏态盆形分布模式。 展开更多
关键词 小迎浪角 直航道 波浪传播 变形 航道
下载PDF
上一页 1 2 6 下一页 到第
使用帮助 返回顶部