期刊文献+
共找到21篇文章
< 1 2 >
每页显示 20 50 100
竖轴潮流能水轮机横荡和艏摇耦合水动力分析 被引量:3
1
作者 姜劲 孙科 +1 位作者 张亮 王树齐 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2020年第4期493-499,共7页
为了分析竖轴潮流能水轮机在无界均匀流下横荡与艏摇耦合运动时的水动力特性,本文采用滑移网格和动网格技术研究了轮机恒定转速下横荡和艏摇同频耦合运动及不同频耦合运动对水轮机水动力载荷的影响规律。研究结果表明:当水轮机发生横荡... 为了分析竖轴潮流能水轮机在无界均匀流下横荡与艏摇耦合运动时的水动力特性,本文采用滑移网格和动网格技术研究了轮机恒定转速下横荡和艏摇同频耦合运动及不同频耦合运动对水轮机水动力载荷的影响规律。研究结果表明:当水轮机发生横荡和艏摇同频耦合运动时,水动力载荷瞬时值明显产生比横荡(或艏摇)运动更大的波动,载荷时变频率与耦合运动频率相同,其幅值随耦合运动频率的增加而增加;当横荡频率比艏摇频率大较多时,水动力载荷的变化规律和横荡运动时的基本一致,但随着艏摇频率的逐渐增加,艏摇运动的影响就会体现出来。研究结果可为漂浮式载体的运动响应研究、水轮机更多自由度运动水动力分析及电能控制等提供相关数据。 展开更多
关键词 潮流能 竖轴轮机 耦合运动 水动力载荷 横荡 艏摇 滑移网格 动网格
在线阅读 下载PDF
穿梭油轮和浮式生产储油轮并靠作业过程中横荡方向运动和漂移力研究(英文) 被引量:1
2
作者 王洪超 王磊 邹刚 《船舶力学》 EI CSCD 北大核心 2014年第3期262-279,共18页
文章针对穿梭油轮和浮式生产储油轮在并靠作业过程中在不同浪向和不同并靠间距情况下的船体一阶耦合运动及平均漂移力等水动力干扰情况进行了研究。通过研究发现,两船之间水动力干扰情况在不同的波浪频率区域,会出现放大或遮蔽效应,这... 文章针对穿梭油轮和浮式生产储油轮在并靠作业过程中在不同浪向和不同并靠间距情况下的船体一阶耦合运动及平均漂移力等水动力干扰情况进行了研究。通过研究发现,两船之间水动力干扰情况在不同的波浪频率区域,会出现放大或遮蔽效应,这与单浮体的水动力干扰情况有着很大的不同。计算结果表明由于另一浮体的存在,横荡平均漂移力不仅在大小上和单浮体水动力情况不同,而且在方向上也发生了改变。同时,鉴于横荡方向的运动和受力在动力定位控制系统中是非常重要的,该文可以为装有动力定位系统的穿梭油轮的安全作业,船体水动力干扰情况预测以及工程实践的初步设计提供相关参考和帮助。 展开更多
关键词 水动力干扰 并靠 横荡耦合运动 平均漂移力 中场法 动力定位系统
在线阅读 下载PDF
浮式装置升沉及横荡运动下海洋立管动力响应研究 被引量:4
3
作者 娄敏 董文乙 王腾 《中国海洋平台》 2010年第4期14-18,共5页
考虑浮式装置升沉及横荡运动、立管内流体流动及立管外波浪荷载的共同作用,建立立管振动方程,采用有限元法进行求解,并分析升沉及横荡运动对海洋立管动力响应的影响规律。
关键词 海洋立管 浮式系统 升沉运动 横荡运动 波浪荷载
在线阅读 下载PDF
单驳船横荡运动对沉管管段横荡运动特性影响的试验研究 被引量:2
4
作者 左卫广 王永学 《中国海洋平台》 2015年第2期57-63,69,共8页
通过物理模型试验,探讨了在波浪条件下,系泊单驳船横荡运动对沉管管段横荡运动的影响。研究结果发现考虑驳船横荡运动后的沉管横荡运动与忽略驳船运动的沉管横荡运动相比,两者存在明显不同。比较了考虑单驳船横荡运动后的沉管横荡运动... 通过物理模型试验,探讨了在波浪条件下,系泊单驳船横荡运动对沉管管段横荡运动的影响。研究结果发现考虑驳船横荡运动后的沉管横荡运动与忽略驳船运动的沉管横荡运动相比,两者存在明显不同。比较了考虑单驳船横荡运动后的沉管横荡运动与单驳船横荡运动之间的相关性,并探讨了不同波浪周期,波高和沉管沉放深度情形下沉管和单驳船横荡运动特性。 展开更多
关键词 沉管 横荡运动 低频运动 波频运动 单驳船
在线阅读 下载PDF
船舶碎冰区横荡运动的水动力导数研究
5
作者 王超 刘正 +2 位作者 朱江波 汪春辉 郭春雨 《船舶力学》 EI CSCD 北大核心 2022年第9期1303-1314,共12页
目前敞水水动力导数的研究较为成熟,针对冰水耦合后的船舶水动力导数研究尚处于起步阶段。本文使用STAR-CCM+软件中DEM颗粒模拟碎冰粒子,求解碎冰区冰水混合介质中船舶横荡运动水动力导数。通过开启DEM模块下双向耦合模式,进行动量、能... 目前敞水水动力导数的研究较为成熟,针对冰水耦合后的船舶水动力导数研究尚处于起步阶段。本文使用STAR-CCM+软件中DEM颗粒模拟碎冰粒子,求解碎冰区冰水混合介质中船舶横荡运动水动力导数。通过开启DEM模块下双向耦合模式,进行动量、能量交换达到冰块与水耦合作用。选取低频频率f分别为0.06 Hz、0.08 Hz、0.1 Hz、0.12 Hz、0.14 Hz进行横荡运动数值模拟,忽略自由液面的影响,分别计算敞水工况以及碎冰工况下船舶各运动频率下所受的侧向力以及转艏力矩,并通过对各频率下无因次后的力以及力矩拟合来解出部分水动力导数。通过计算得知,敞水工况下船舶各部分水动力导数值与统计公式计算值相差不大,冰水耦合后大部分水动力导数在低频率下明显大于敞水工况下的水动力导数,而水动力导数Yv’值比敞水工况值小。 展开更多
关键词 横荡运动 碎冰区 冰水耦合 DEM 冰水耦合水动力导数
在线阅读 下载PDF
非线性横荡的数值计算和模型试验
6
作者 杨驰 刘应中 高木仪昌 《上海交通大学学报》 EI CAS CSCD 北大核心 1990年第4期1-7,共7页
本文用时间步进和边界元法,对在静水中由静止开始作横荡运动的三维柱体的流体动力问题进行了数值模拟,并在日本大学作了相应的模型试验.在计算中,自由面采用拉格朗日网格,物面和自由面的交界点同时满足其边界条件,外边界用不给定方向的... 本文用时间步进和边界元法,对在静水中由静止开始作横荡运动的三维柱体的流体动力问题进行了数值模拟,并在日本大学作了相应的模型试验.在计算中,自由面采用拉格朗日网格,物面和自由面的交界点同时满足其边界条件,外边界用不给定方向的单波假定.数值计算结果和模型试验结果吻合较好. 展开更多
关键词 三维柱体 非线性 横荡 流体动力
在线阅读 下载PDF
浮式风力机横荡运动下气动特性分析 被引量:3
7
作者 黄亚振 陈嘉佳 +2 位作者 沈昕 吴亚东 田杰 《太阳能学报》 EI CAS CSCD 北大核心 2021年第9期279-285,共7页
浮式平台的横荡运动会影响到风力机相对入流速度以及叶轮-尾流的相互作用,导致浮式风力机的气动特性更为复杂。基于升力面自由尾迹模型,以NREL 5 MW全尺寸浮式风力机作为研究对象,分析在横荡运动下浮式风力机整机性能以及叶片非定常气... 浮式平台的横荡运动会影响到风力机相对入流速度以及叶轮-尾流的相互作用,导致浮式风力机的气动特性更为复杂。基于升力面自由尾迹模型,以NREL 5 MW全尺寸浮式风力机作为研究对象,分析在横荡运动下浮式风力机整机性能以及叶片非定常气动特性,并探究其作用机理。研究表明,横荡运行下,远场尾迹对叶片的诱导作用减小,整机的功率及推力略有提升;横荡速度附加在旋转方向上导致功率及推力与横荡运动存在一定的相位差;叶片上切向力及法向力呈现不同的分布,频率为横荡运动的2倍,叶片变形会相应增大,叶尖与塔架碰撞风险增大。 展开更多
关键词 海上风力机 浮式平台 升力面理论 非定常气动性能 横荡
在线阅读 下载PDF
穿透自由面平板横荡和首摇水动力计算
8
作者 邹早建 《武汉交通科技大学学报》 1996年第6期641-648,共8页
介绍了一种计算穿透自由面垂直平板横荡和首摇运动水动力的三维Rankine面元法.采用在平板平均位置上及其尾流中的涡强谐变的涡分布,同时,采用在围绕平板的自由面上的源强谐变的源分布以计入自由面影响.将所介绍的方法应用于... 介绍了一种计算穿透自由面垂直平板横荡和首摇运动水动力的三维Rankine面元法.采用在平板平均位置上及其尾流中的涡强谐变的涡分布,同时,采用在围绕平板的自由面上的源强谐变的源分布以计入自由面影响.将所介绍的方法应用于小展弦比(吃水长度比)垂直平板。 展开更多
关键词 垂直平板 横荡 首摇振荡 水动力 船舶 计算
在线阅读 下载PDF
安装船横荡运动下沉管管节运动特性的数值研究 被引量:1
9
作者 刘明 张姝如 左卫广 《中国港湾建设》 2017年第6期17-21,共5页
考虑安装船做强迫横荡运动,研究沉管管节运动特性。管节运动包括仅由安装船横荡运动引起的管节运动和单独波浪影响下的管节运动两部分。应用有限水深条件下的格林函数求解控制方程,考虑安装船横荡运动,并采用四阶Runge-Kutta方法分别求... 考虑安装船做强迫横荡运动,研究沉管管节运动特性。管节运动包括仅由安装船横荡运动引起的管节运动和单独波浪影响下的管节运动两部分。应用有限水深条件下的格林函数求解控制方程,考虑安装船横荡运动,并采用四阶Runge-Kutta方法分别求解安装船瞬时运动引起的管节位移变化和管节的运动方程。通过数值算例分析,探讨了安装船做不同振幅和相位差的横荡运动和波浪联合作用下沉管管节的运动特性。 展开更多
关键词 沉管 横荡运动 低频运动 波频运动 安装船
在线阅读 下载PDF
大幅激励下液舱晃荡特性研究
10
作者 龚宸 邹昶方 傅振航 《舰船科学技术》 北大核心 2024年第20期43-49,共7页
晃荡载荷是大型液货船结构设计中的关键问题之一。本文基于数值模拟方法研究液舱在大幅激励下的晃荡特性问题,考虑的因素包括装载率、大幅横荡激励、大幅横荡-垂荡激励。采用基于有限体积法的CFD代码Fluent求解器,通过动网格技术驱动液... 晃荡载荷是大型液货船结构设计中的关键问题之一。本文基于数值模拟方法研究液舱在大幅激励下的晃荡特性问题,考虑的因素包括装载率、大幅横荡激励、大幅横荡-垂荡激励。采用基于有限体积法的CFD代码Fluent求解器,通过动网格技术驱动液舱运动并用VOF法追踪自由液面的非线性运动过程。研究结果表明,在大幅激励下,液体一阶固有频率偏离理论值;随着激励幅值的增大,晃荡压力呈线性增加到非线性增加,最终达到饱和的变化特性;对比受大幅横荡-垂荡耦合激励的液舱和受大幅横荡激励的液舱,其压力特性在低装载下相似,而中、高装载下,受大幅横荡-垂荡激励的液舱晃荡压力更高。 展开更多
关键词 冲击压力 大幅度横荡激励 大幅度横荡-垂荡耦合激励 共振频率
在线阅读 下载PDF
锚系船深水导管架上部组块吊装就位方案研究
11
作者 沈枕 王燕 +3 位作者 莫燕 连建国 李英俊 孙华峰 《石油和化工设备》 CAS 2024年第1期133-136,共4页
利用锚系船在深水海域作业时,受长周期波浪的影响,浮吊船的运动不可避免。为了保证上部组块就位安全,需要从设计方面给出参考。由于浮吊船船体呈细长型,在纵、横两个方向上运动的反应明显不同,纵荡运动衰减很快,并且对长周期波浪的反应... 利用锚系船在深水海域作业时,受长周期波浪的影响,浮吊船的运动不可避免。为了保证上部组块就位安全,需要从设计方面给出参考。由于浮吊船船体呈细长型,在纵、横两个方向上运动的反应明显不同,纵荡运动衰减很快,并且对长周期波浪的反应不敏感,运动幅值也约为横荡运动幅值的50%。在考虑了锚缆悬链线形态和锚缆破断力后,建议在进行上部组块吊装时,沿着浮吊船纵荡运动的方向靠近导管架就位。 展开更多
关键词 深水 浮吊船 纵荡 横荡
在线阅读 下载PDF
平板及立板式防波堤散射波浪力研究 被引量:5
12
作者 王科 高鑫 张犀 《船舶力学》 EI 北大核心 2011年第10期1127-1133,共7页
文章应用边界单元方法研究了波浪与平板及立板两种防波堤结构的相互作用,揭示了这两种型式防波堤的散射波浪力与板的相对板厚、相对潜深及相对板长之间的关系。研究发现:(1)有效波浪力存在的前提条件是潜深小于0.5m并且板厚为O(0.005m);... 文章应用边界单元方法研究了波浪与平板及立板两种防波堤结构的相互作用,揭示了这两种型式防波堤的散射波浪力与板的相对板厚、相对潜深及相对板长之间的关系。研究发现:(1)有效波浪力存在的前提条件是潜深小于0.5m并且板厚为O(0.005m);(2)水面和水下平板的垂荡波浪力为控制要素并有极值存在,横摇波浪力可以忽略;(3)水面或者水下立板的横荡波浪力为控制要素,并且横摇波浪力约为平板的20倍。 展开更多
关键词 边界单元法 板式防波堤 横荡 垂荡 横摇波浪力
在线阅读 下载PDF
基于实验的二维八边形液舱晃荡比尺效应研究 被引量:3
13
作者 杨志勋 徐潜岳 +3 位作者 董玉山 阮诗伦 吴尚华 岳前进 《船舶力学》 EI CSCD 北大核心 2019年第5期505-511,共7页
文章针对二维八边形液舱晃荡冲击的比尺效应开展了实验研究。实验根据弗劳德相似准则分别选择1:1原型液舱与1:2和1:3两个比尺模型液舱,在横荡单次冲击外激激励的作用下,考虑80%危险载液率分别对流体冲击侧壁、上斜板和顶板时的晃荡荷载... 文章针对二维八边形液舱晃荡冲击的比尺效应开展了实验研究。实验根据弗劳德相似准则分别选择1:1原型液舱与1:2和1:3两个比尺模型液舱,在横荡单次冲击外激激励的作用下,考虑80%危险载液率分别对流体冲击侧壁、上斜板和顶板时的晃荡荷载与冲击波形进行了比尺效应对比分析。分析结果表明,冲击过程中,如果原型液舱和缩比模型液舱中对应冲击壁面处于流体中,包裹的气泡越少,冲击压力峰值由模型还原到原型时更接近于原型实验实测值;若考虑整个冲击过程中的流场形态,应尽可能选择大比尺液舱来开展晃荡模型实验。 展开更多
关键词 二维八边形液舱 比尺效应 晃荡实验 横荡单次冲击
在线阅读 下载PDF
舰船在水平面运动中非线性模型的建立 被引量:3
14
作者 张冰 姜长生 《舰船科学技术》 北大核心 2005年第5期25-27,共3页
在选择了描述船舶运动坐标系的前提下,根据刚体动力学理论和牛顿第二运动定律,推导了舰船六自由度的非线性模型,并在此基础上建立了舰船在水平面中的非线性运动模型。同时,根据平稳随机过程理论对海浪的分析可知,利用白噪声驱动所设计... 在选择了描述船舶运动坐标系的前提下,根据刚体动力学理论和牛顿第二运动定律,推导了舰船六自由度的非线性模型,并在此基础上建立了舰船在水平面中的非线性运动模型。同时,根据平稳随机过程理论对海浪的分析可知,利用白噪声驱动所设计的成形滤波器,能够获得海浪干扰的数学模型。上述模型既能很好地反映舰船的水动力特性,又能与原始波浪力、力矩相吻合,由此证明了模型的正确性。 展开更多
关键词 船舶运动 非线性模型 横荡 横摇 首摇 成形滤波器 模糊建模
在线阅读 下载PDF
二代稳性衡准纯稳性丧失直接评估方法研究(二) 被引量:1
15
作者 鲁江 王田华 +1 位作者 储纪龙 顾民 《船舶力学》 EI CSCD 北大核心 2017年第z1期158-167,共10页
针对国际海事组织(IMO)正在制定的船舶第二代完整稳性衡准中的纯稳性丧失直接评估衡准,本文提出了纵荡-横荡-横摇-首摇4自由度运动耦合的标准数学模型,进行了纯稳性丧失直接评估.该方法首先基于MMG操纵性标准方法,构建纵荡-横荡-横摇-首... 针对国际海事组织(IMO)正在制定的船舶第二代完整稳性衡准中的纯稳性丧失直接评估衡准,本文提出了纵荡-横荡-横摇-首摇4自由度运动耦合的标准数学模型,进行了纯稳性丧失直接评估.该方法首先基于MMG操纵性标准方法,构建纵荡-横荡-横摇-首摇4自由度运动方程,同时考虑了舵控制方程,其次在纵荡、横荡、横摇、首摇方程右边考虑时域波浪力/力矩,且横摇方程右边进一步考虑了粘性横摇阻尼力矩和波浪中复原力臂变化.最后采用ONR内倾船进行了尾斜浪中纯稳性丧失直接数值计算,为纯稳性丧失直接稳性评估衡准的应用奠定了基础. 展开更多
关键词 纯稳性丧失 稳性直接评估 第二代完整稳性衡准 MMG操纵性 纵荡-横荡-横摇-首摇
在线阅读 下载PDF
深水半潜式平台的涡激运动Ⅱ.船型优化研究
16
作者 沈志平 单铁兵 +2 位作者 潘方豪 张海彬 王璞 《上海交通大学学报》 EI CAS CSCD 北大核心 2017年第12期1504-1511,共8页
基于CFD数值计算方法对深水半潜式平台的涡激运动开展优化研究,提出一种能够抑制涡激运动的半潜平台新形式,并应用模型试验进行对比验证.研究了新型平台和传统型平台的涡激运动特征,对纵向偏移、横荡响应的幅值和分布规律,锁定区间等关... 基于CFD数值计算方法对深水半潜式平台的涡激运动开展优化研究,提出一种能够抑制涡激运动的半潜平台新形式,并应用模型试验进行对比验证.研究了新型平台和传统型平台的涡激运动特征,对纵向偏移、横荡响应的幅值和分布规律,锁定区间等关键因素进行对比分析;研究了各方向耦合作用下平台的轨迹特征;详细分析了新型平台的减涡机制,从平台泄涡结构和涡激响应的均匀性、稳定性以及响应幅值等方面展开研究,阐述立柱的横剖面形状,排列形式等对平台涡激运动性能的影响.研究表明,这种"立柱的横剖面形状为梯形,呈外八字对称排列"的新型平台形式对抑制平台涡激运动具有积极影响,有助于减少或延缓强非线性涡激运动对系泊、立管系统的破坏作用. 展开更多
关键词 涡激运动 横荡响应 锁定现象 减涡机制 横剖面形状
在线阅读 下载PDF
深水半潜式平台的涡激运动Ⅰ.关键特性研究
17
作者 单铁兵 沈志平 《上海交通大学学报》 EI CAS CSCD 北大核心 2017年第12期1493-1503,共11页
采用数值计算和模型试验相结合的方法对深水半潜式平台涡激运动特性进行了系统研究.模型试验在拖曳水池内开展,采用4点水平系泊方式;数值计算选用黏性流的计算流体力学方法,建立适用于研究平台涡激运动的分析方法,开发合理的网格划分模... 采用数值计算和模型试验相结合的方法对深水半潜式平台涡激运动特性进行了系统研究.模型试验在拖曳水池内开展,采用4点水平系泊方式;数值计算选用黏性流的计算流体力学方法,建立适用于研究平台涡激运动的分析方法,开发合理的网格划分模式,引入网格变形和重构技术解决非线性的流固耦合问题;研究半潜式平台在不同来流角度下的横荡运动、纵向偏移特征;深入研究涡激锁定现象,给出不同角度下的漩涡锁定区间,获得涡激锁定的相关特性;研究平台涡激响应随折合速度所经历的演化过程;揭示了平台涡激运动随时间变化的发展阶段,分析各阶段漩涡泄放的特点. 展开更多
关键词 涡激运动 流固耦合 横荡运动 锁定现象 折合速度
在线阅读 下载PDF
复合运动下风力机翼型动态失速特性研究 被引量:3
18
作者 冯俊鑫 赵振宙 +3 位作者 陈明 江瑞芳 刘一格 王丁丁 《太阳能学报》 EI CAS CSCD 北大核心 2023年第8期430-436,共7页
采用带γ-Reθ转捩修正的SST k-ω湍流模型和动网格技术,以S809翼型为研究对象,探讨复合(俯仰+前后横荡)运动下的动态失速特性。结果表明:横荡运动扩大了翼型升阻力系数迟滞环状;上仰阶段升力系数增加,阻力系数减小,增减幅度与横荡幅度... 采用带γ-Reθ转捩修正的SST k-ω湍流模型和动网格技术,以S809翼型为研究对象,探讨复合(俯仰+前后横荡)运动下的动态失速特性。结果表明:横荡运动扩大了翼型升阻力系数迟滞环状;上仰阶段升力系数增加,阻力系数减小,增减幅度与横荡幅度成正比;降低了下俯阶段的升力系数。横荡运动通过改变翼型边界层流体动能,有效缩小翼型尾缘涡的拓扑结构,增强前缘涡的诱导效应,使得翼型上仰阶段失速程度减弱。 展开更多
关键词 海上风电 风力机 空气动力学 横荡运动 俯仰运动 动态失速 S809翼型
在线阅读 下载PDF
驳船横摇运动下波浪对沉管管段运动特性的影响 被引量:1
19
作者 刘明 张俊 左卫广 《中国港湾建设》 2017年第10期42-46,共5页
沉管管段在沉放过程中,管段运动不仅受海面风、浪、流的影响,还受到驳船运动影响,属于复杂的水动力学问题。为探讨驳船运动对沉管管段运动的影响,假设驳船作强迫横摇运动,沉管管段运动由波浪作用下的沉管管段运动和驳船横摇运动引起的... 沉管管段在沉放过程中,管段运动不仅受海面风、浪、流的影响,还受到驳船运动影响,属于复杂的水动力学问题。为探讨驳船运动对沉管管段运动的影响,假设驳船作强迫横摇运动,沉管管段运动由波浪作用下的沉管管段运动和驳船横摇运动引起的沉管管段运动两部分组成。建立波浪作用下的沉管运动时域模型,应用有限水深条件下的格林函数求解控制方程得到速度势。通过瞬时驳船横摇运动改变得到沉管管段运动新位移,在新位置计算速度势得到波浪力,带入沉管管段运动方程即可得到沉管运动。最后,探讨了不同工况下沉管管段的运动特性,并给出了相关结论。 展开更多
关键词 驳船 沉管 横摇运动 横荡运动 低频运动
在线阅读 下载PDF
乐器趣闻
20
作者 沈洪 《音乐爱好者》 1994年第5期34-34,共1页
有位意大利音乐家持一支横笛在欧洲各地巡回演奏,所到之处听众如醉如痴。他吹奏的乐曲悦耳动听,所闻的横笛更不同一般。这横笛竟是用人骨制做的!原来音乐家的一位同事逝世前立下遗嘱说,希望拿他的骨头去制做乐器,让他死后仍能发出美妙... 有位意大利音乐家持一支横笛在欧洲各地巡回演奏,所到之处听众如醉如痴。他吹奏的乐曲悦耳动听,所闻的横笛更不同一般。这横笛竟是用人骨制做的!原来音乐家的一位同事逝世前立下遗嘱说,希望拿他的骨头去制做乐器,让他死后仍能发出美妙的乐声。音乐家便用他的胫骨制成一支横笛,进行吹奏,用精彩的演出实现同事继续献身乐坛的遗愿。 荷兰的马特支城有个奇特的乐队,是由家庭主妇们组成的。她们所用的乐器也很奇特,全部是家庭厨房中的炊具和餐具,如煎锅、炒勺、菜碟、茶盘,以及刀、 展开更多
关键词 家庭厨房 人骨 悦耳动听 曲终 乐器演奏 成宪 特纳 米尼 横荡 神谷
在线阅读 下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部