期刊文献+
共找到45篇文章
< 1 2 3 >
每页显示 20 50 100
Independence Polynomials and the Merrifield-Simmons Index of Mono-Layer Cylindrical Grid Graphs
1
作者 JI Lin-xing ZHANG Ke HU Wen-jun 《Chinese Quarterly Journal of Mathematics》 2024年第4期379-387,共9页
Research on the independence polynomial of graphs has been very active.However,the computational complexity of determining independence polynomials for general graphs remains NP-hard.Letα(G)be the independence number... Research on the independence polynomial of graphs has been very active.However,the computational complexity of determining independence polynomials for general graphs remains NP-hard.Letα(G)be the independence number of G and i(G;k)be the number of independent sets of order k in G,then the independence polynomial is defined as I(G;x)=∑_(k=0)^(α(G))i(G;k)x^(k),i(G;0)=1.In this paper,by utilizing the transfer matrix,we obtain an analytical expression for I(CGn;x)of mono-cylindrical grid graphs CGn and present a crucial proof of it.Moreover,we also explore the Merrifield-Simmons index and other properties of CGn. 展开更多
关键词 Independence polynomial cylindrical grid graphs Transfer matrix Merrifield-Simmons index
下载PDF
Acoustic field simulations of logging while drilling by cylindrical finite difference with variable grids 被引量:1
2
作者 Pan Yue He Xiao +1 位作者 Chen Hao Wang Xiu-Ming 《Applied Geophysics》 SCIE CSCD 2022年第1期11-28,143,共19页
This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrica... This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrical FD as the grid size gradually becomes more extensive with the increasing radius.To prevent grids from being too coarse in far fields,we compensate for the grid cell infl ation by refi ning the grid step in the azimuthal direction.The variable grid FD in the cylindrical coordinate systems has a higher effi ciency in solving acoustic logging while drilling(LWD)problems because the grid boundaries are consistent with those of the drill collar and the borehole.The proposed algorithm saves approximately 94%of the FD grids,80%of the computation time,and memory with a higher calculation accuracy than the FD on rectangular grids for the same models.We also calculate the acoustic LWD responses of the fl uid-fi lled borehole intersecting with fractures.Refl ections are generated at the fractures,which can be equivalent to an additional scattering source.The mode conversions between the collar and the Stoneley waves are revealed.The Stoneley spectra are more sensitive to the fracture.Finally,the logs in a heterogeneous formation with two refl ectors far from the borehole are modeled,and a means of estimating the azimuth of geological interfaces from refl ections is proposed. 展开更多
关键词 acoustic detection finite difference variable grid method scattering and refl ection cylindrical guided wave
下载PDF
Comparisons of some difference forms for compressible flow in cylindrical geometry on arbitrary Lagrangian and Eulerian framework 被引量:1
3
作者 Zhijun SHEN Xiao LI Jian REN 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI CSCD 2016年第11期1571-1586,共16页
The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical sym- metry properties of the numerical schemes are discus... The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical sym- metry properties of the numerical schemes are discussed, and an area weighted scheme is extended from a Lagrangian method to an arbitrary Lagrangian and Eulerian (ALE) method. Numerical results are presented to compare three discrete configurations, i.e., the control volume scheme, the area weighted scheme, and the plane scheme with the addition of a geometrical source. The fact that the singularity arises from the geometri- cal source term in the plane scheme is illustrated. A suggestion for choosing the discrete formulation is given when the strong shock wave problems are simulated. 展开更多
关键词 arbitrary Lagrangian and Eulerian (ALE) cylindrical geometry symmetrypreservation unstructured grid
下载PDF
用于高功率微波器件粒子模拟的三维柱坐标系共形网格生成技术
4
作者 张磊 翁明 +1 位作者 王玥 蒋铭 《强激光与粒子束》 CAS CSCD 北大核心 2024年第4期180-186,共7页
基于面向对象的C++语言研制三维柱坐标共形网格生成程序,对束-场互作用器件作共形网格剖分,为粒子模拟算法提供积分线元、面元。通过定义三维柱坐标网格体系如网格步长、网格索引、守护网格层、包围盒等数据,使模型的空间信息能转换成... 基于面向对象的C++语言研制三维柱坐标共形网格生成程序,对束-场互作用器件作共形网格剖分,为粒子模拟算法提供积分线元、面元。通过定义三维柱坐标网格体系如网格步长、网格索引、守护网格层、包围盒等数据,使模型的空间信息能转换成柱坐标网格信息。将轴上网格单元作特殊处理,使粒子模拟算法形式在轴上网格和在非轴上网格上保持一致。利用射线跟踪法得到属于模型子面、模型棱边的离散边界点,接着通过拓扑关系获得模型的顶点,保存上述三类离散边界点的拓扑信息和网格信息。将构建的基础网格元与边界点信息耦合,在离散网格体系中重构模型。用该共形网格剖分技术对相对论磁控管进行剖分,能够识别该磁控管的透明阴极、阳极和谐振腔等结构。 展开更多
关键词 粒子模拟 三维柱坐标 网格剖分 射线跟踪 共形技术
下载PDF
不同刚度圆柱振子涡激振动数值模拟
5
作者 宋玉泉 谭俊哲 +3 位作者 袁鹏 司先才 杜修茂 李昌明 《中国海洋大学学报(自然科学版)》 CAS CSCD 北大核心 2024年第9期133-140,共8页
为了研究不同刚度涡激振动系统对圆柱振子振动特性的影响,本文建立了单自由度涡激振动系统的仿真模型,采用计算流体力学(CFD)方法并结合重叠网格技术进行数值模拟。通过与公开发表的文献数据进行对比,验证了模型的精度和方法的可靠性。... 为了研究不同刚度涡激振动系统对圆柱振子振动特性的影响,本文建立了单自由度涡激振动系统的仿真模型,采用计算流体力学(CFD)方法并结合重叠网格技术进行数值模拟。通过与公开发表的文献数据进行对比,验证了模型的精度和方法的可靠性。本文对时变流速下圆柱振子的振幅响应、频率响应和尾涡特征的变化情况进行了详细分析。仿真结果表明:当圆柱振子的振幅随流速的改变而降低时,调整弹簧刚度的大小,可使其保持较大的无量纲振幅和振动频率,并使尾流脱落旋涡之间的间距增大,且脱落模式由“2S”模式变为“2P”模式。 展开更多
关键词 涡激振动 弹簧刚度 计算流体力学 重叠网格 圆柱振子
下载PDF
深吃水圆筒型浮式核能平台涡激运动数值模拟
6
作者 曹辰泽 何炎平 +1 位作者 王梓 刘亚东 《船舶力学》 EI CSCD 北大核心 2024年第2期239-249,共11页
深吃水圆筒型浮式核能平台是一种新型多功能高效平台,可有效解决南海岛礁开发过程中的能源供给问题。在一定来流速度下,尾流区交替泄涡进而诱导平台发生涡激运动(vortex induced motions,VIM),这将严重加速系泊和立管系统疲劳损害,同时... 深吃水圆筒型浮式核能平台是一种新型多功能高效平台,可有效解决南海岛礁开发过程中的能源供给问题。在一定来流速度下,尾流区交替泄涡进而诱导平台发生涡激运动(vortex induced motions,VIM),这将严重加速系泊和立管系统疲劳损害,同时对平台内部核反应堆运行产生不利影响。基于改进的延迟分离涡方法(improved delayed detached eddy simulation,IDDES)对平台在不同折合速度下的横荡、纵荡、艏摇运动响应进行数值模拟,并从水平面内质心运动轨迹、运动频率、三维流场特性等角度分析涡激运动关键特征。研究结果表明:当折合速度5.45<Ur<9.08时,平台横荡、艏摇振幅均逐渐增加且运动轨迹类似“香蕉”形,横荡与艏摇运动频率基本一致且横荡未出现明显的“锁定”区间变化;当7.26<Ur<9.08时,艏摇振幅近似线性递增且运动频率中主频附近出现多个峰值,运动轨迹在顺流方向上逐渐变宽;三维流场中发现尾流区三维漩涡结构相当复杂且平台底部柱靴结构对表面流动分离造成干扰,柱靴结构具有一定的减涡效果。 展开更多
关键词 深吃水圆筒型浮式核能平台 涡激运动 IDDES 重叠网格法
下载PDF
基于柱坐标系的可变密度栅格无人机终端区航路规划
7
作者 郜雅 罗喜伶 +4 位作者 周家柳 徐文翔 王宇鹏 李宇雄 刘刚 《西华大学学报(自然科学版)》 CAS 2024年第6期8-13,共6页
由于城市低空风险分布不均匀,无人机起降场终端区附近障碍物密集,多无人机在起降场终端区汇聚运行,传统基于直角坐标系的均匀密度栅格法难以解决该特定场景下的航路规划。针对以上难题,本文基于柱坐标系构建空域模型,分别提出基于A^(*)... 由于城市低空风险分布不均匀,无人机起降场终端区附近障碍物密集,多无人机在起降场终端区汇聚运行,传统基于直角坐标系的均匀密度栅格法难以解决该特定场景下的航路规划。针对以上难题,本文基于柱坐标系构建空域模型,分别提出基于A^(*)算法与遗传算法的可变密度栅格法用于无人机终端区的航路规划。仿真结果表明,基于A^(*)算法的柱坐标系可变密度栅格法比基于A^(*)传算法的柱坐标系可变密度栅格法比基于遗传算法的直角坐标系传统栅格法的航路规划效率提高了71.72%,路径总长度缩短了1.29%。本文方法解决了传统栅格法无法兼顾规划效率与环境描述精度的研究难题。 展开更多
关键词 无人机航路规划 可变密度栅格 空域建模 柱坐标系
下载PDF
含初始缺陷复合材料格栅加筋圆柱壳的鲁棒优化设计 被引量:9
8
作者 张志峰 陈浩然 白瑞祥 《固体力学学报》 CAS CSCD 北大核心 2006年第1期58-64,共7页
针对复合材料格栅加筋圆柱壳具有多约束、连续变量和非连续变量混合的特性,基于复合材料中厚圆柱壳体的非线性稳定性理论,考虑优化目标和约束鲁棒性,采用混合遗传算法,提出了一种非完善复合材料格栅加筋圆柱壳的鲁棒优化设计方法.应用... 针对复合材料格栅加筋圆柱壳具有多约束、连续变量和非连续变量混合的特性,基于复合材料中厚圆柱壳体的非线性稳定性理论,考虑优化目标和约束鲁棒性,采用混合遗传算法,提出了一种非完善复合材料格栅加筋圆柱壳的鲁棒优化设计方法.应用正交设计理论确定主要的鲁棒优化设计变量,简化了运算过程,从而提高了优化分析效率;同时又采用五点差分格式代替梯度运算,解决了在优化计算过程中不可求导的问题.通过典型算例表明,鲁棒优化设计所得到的优化分析结果与传统确定性优化方法相比差别较大,此外壳体的初始缺陷亦将对优化结果产生很大的影响. 展开更多
关键词 鲁棒优化设计 初始缺陷 复合材料格栅加筋圆柱壳
下载PDF
激光微织构滚动轴承表面润滑性能的数值分析 被引量:10
9
作者 华希俊 许洪山 +3 位作者 陈亚林 解玄 徐上 王成 《表面技术》 EI CAS CSCD 北大核心 2018年第3期36-41,共6页
目的研究圆柱滚子轴承摩擦副织构化表面润滑性能,探讨工况参数(载荷、转速)和微织构几何形貌参数(面积占有率)对油膜平均无量纲压力的影响规律。方法基于弹性流体动压润滑机理,利用有限长线接触弹流润滑理论建立理论模型,采用有限差分法... 目的研究圆柱滚子轴承摩擦副织构化表面润滑性能,探讨工况参数(载荷、转速)和微织构几何形貌参数(面积占有率)对油膜平均无量纲压力的影响规律。方法基于弹性流体动压润滑机理,利用有限长线接触弹流润滑理论建立理论模型,采用有限差分法对Reynolds方程进行离散,并运用多重网格法进行数值分析。结果在微织构面积占有率S_p=0.095、深度h_p=6μm的前提下,当载荷由100 N升为500 N时,油膜平均无量纲压力随之增大。在微织构面积占有率S_p=0.095、深度h_p=6μm的前提下,当转速由100 r/min升为500 r/min时,油膜平均无量纲压力随之增大。在载荷W=200 N、转速n=200 r/min的前提下,当微织构面积占有率由0.05升为0.15时,油膜平均无量纲压力随之增大;当微织构面积占有率由0.15升为0.25时,油膜平均无量纲压力随之减小。结论圆柱滚子轴承微织构表面的油膜平均无量纲压力随着载荷和转速的增大而增大,随着微织构面积占有率的增大而先增大后减小。 展开更多
关键词 圆柱滚子轴承 织构化表面 多重网格法 油膜平均无量纲压力
下载PDF
柱形和球形壳体内爆界面不稳定性数值分析 被引量:1
10
作者 柏劲松 王涛 +1 位作者 刘坤 钟敏 《应用力学学报》 CAS CSCD 北大核心 2012年第5期601-606,633,共6页
为了更好地研究柱形和球形构型下果冻界面不稳定性发展,避免内爆聚心反弹前后直角坐标网格计算导致的误差,提高对流场和界面位置的计算精度,通过应用考虑了MVFT程序的网格适应性,使其能够适用于柱形网格和球形网格下的界面不稳定性数值... 为了更好地研究柱形和球形构型下果冻界面不稳定性发展,避免内爆聚心反弹前后直角坐标网格计算导致的误差,提高对流场和界面位置的计算精度,通过应用考虑了MVFT程序的网格适应性,使其能够适用于柱形网格和球形网格下的界面不稳定性数值模拟,特别是能够保证内爆聚心反弹前后流场和界面计算的稳定性。应用改进的计算程序对两种构型下的界面不稳定性进行了数值模拟,并对二者界面演化规律进行了详细讨论和归纳。结果表明:对于内外半径相同的柱形和球形果冻,聚心反弹时前者半径较小,而后者反弹时刻早于前者,其向内聚心和向外运动的速度最大值大于前者,对内部气体的压缩强度强于前者。对于外边界带有正弦扰动情况,除遵循上述规律外,计算还给出了峰谷转换现象。该项研究结果为进一步深入进行复杂构型下界面不稳定性高精度数值模拟研究提供了一种分析工具。 展开更多
关键词 界面不稳定性 数值模拟 MVFT 柱形网格 球形网格
下载PDF
基于物体表面散乱三维激光扫描点的三角形格网建立 被引量:8
11
作者 郑德华 张云涛 《测绘工程》 CSCD 2004年第4期62-65,共4页
针对物体表面的散乱三维激光扫描点数据,提出了借助柱面坐标系统将三维激光扫描点转换为三角形格网构造系坐标,并根据转换后的坐标利用在平面内计算狄罗里三角形的方法建立三维物体表面格网模型。建立格网模型的效果表明,借助柱面坐标... 针对物体表面的散乱三维激光扫描点数据,提出了借助柱面坐标系统将三维激光扫描点转换为三角形格网构造系坐标,并根据转换后的坐标利用在平面内计算狄罗里三角形的方法建立三维物体表面格网模型。建立格网模型的效果表明,借助柱面坐标系统构造狄罗里三角形的方法可以作为三维物体表面三角形格网模型的一个重要的基本方法。 展开更多
关键词 物体表面 激光扫描 三维物体 效果 基本方法 数据 利用 格网 系统构造 转换
下载PDF
柱状交界面Richtmyer-Meshkov不稳定性的数值研究 被引量:2
12
作者 马东军 何兴 +1 位作者 孙德军 尹协远 《爆炸与冲击》 EI CAS CSCD 北大核心 2003年第5期398-404,共7页
采用基于有限体积法的二阶Godunov格式模拟了柱形密度交界面在同轴激波加速下的演化过程。得到了以下初步结果:在参数相同的情形下,内聚激波比中心爆炸波对界面的扰动更危险;内聚激波从轻质流体进入重质流体与从重质流体进入轻质流体相... 采用基于有限体积法的二阶Godunov格式模拟了柱形密度交界面在同轴激波加速下的演化过程。得到了以下初步结果:在参数相同的情形下,内聚激波比中心爆炸波对界面的扰动更危险;内聚激波从轻质流体进入重质流体与从重质流体进入轻质流体相比,界面有更快的增长。周向波数大小对界面增长率有很大的影响,在计算的参数下,n=8~12有最大增长率,大于和小于这个范围的周向波长,增长率均明显减弱。 展开更多
关键词 柱状交界面 RICHTMYER-MESHKOV不稳定性 二阶Godunov格式 爆炸波 激波 流体力学 爆炸力学 自适应非结构网格 数值计算
下载PDF
双层圆柱面和球面网壳的风振系数实用公式 被引量:5
13
作者 罗尧治 王彬 《空间结构》 CSCD 北大核心 2008年第1期32-39,共8页
运用频域法分析了双层圆柱面网壳和球面网壳两种典型大跨度空间结构在常见工况下的风振响应,研究了参振模态的合理截取,提出了位移和内力一致风振系数的表达式,考察了对风振系数的变化起关键作用的结构和环境的参数和因素。在两类网壳... 运用频域法分析了双层圆柱面网壳和球面网壳两种典型大跨度空间结构在常见工况下的风振响应,研究了参振模态的合理截取,提出了位移和内力一致风振系数的表达式,考察了对风振系数的变化起关键作用的结构和环境的参数和因素。在两类网壳结构的跨度、长宽比、矢跨比、地貌类别等主要参数的常见变化范围内,进行大规模算例计算,定量分析对结构风振系数的影响,并在此基础上归纳出风振系数的实用公式,以供工程设计参考. 展开更多
关键词 圆柱面网壳 球面网壳 风振响应 风振系数 频域分析 实用公式
下载PDF
分析非均匀填充圆柱介质谐振器的二阶精度FDTD公式 被引量:4
14
作者 褚庆昕 冯英 《电波科学学报》 EI CSCD 2003年第6期652-654,667,共4页
从积分形式的Maxwell方程出发 ,利用连续函数的Taylor级数展开 ,严格地给出了包含介质交界面的二阶精度时域有限差分 (FDTD)公式 ,解决了以往FDTD法处理非均匀介质填充区域问题时只有一阶精度的问题。分析表明 ,为了获得二阶精度 ,除了... 从积分形式的Maxwell方程出发 ,利用连续函数的Taylor级数展开 ,严格地给出了包含介质交界面的二阶精度时域有限差分 (FDTD)公式 ,解决了以往FDTD法处理非均匀介质填充区域问题时只有一阶精度的问题。分析表明 ,为了获得二阶精度 ,除了需要引入适当等效介电常数外 ,还必须采用适当非均匀网格。该方法被用于轴对称圆柱介质谐振器的分析。计算结果与理论值吻合良好 。 展开更多
关键词 时域有限差分法 介质交界面 二阶精度 非均匀网络 圆柱介质谐振器 FDTD
下载PDF
基于圆柱坐标非均均FDTD网格的图形可视化 被引量:1
15
作者 周国祥 侯整风 +2 位作者 郭骏 石雷 杨明武 《工程图学学报》 CSCD 2004年第4期67-71,共5页
基于非均匀 FDTD 网格的电磁场模拟是解决电磁场问题的一种有效方法。介绍一个柱坐标系下非均匀 FDTD 网格生成系统,实现了网格图形的可视化。该系统能够快速、有效地生成一维、二维、三维 FDTD 网格图形和物体几何图形。这些图形不仅... 基于非均匀 FDTD 网格的电磁场模拟是解决电磁场问题的一种有效方法。介绍一个柱坐标系下非均匀 FDTD 网格生成系统,实现了网格图形的可视化。该系统能够快速、有效地生成一维、二维、三维 FDTD 网格图形和物体几何图形。这些图形不仅直观准确,还具有动态和可调等特性,从而为 FDTD 用户提供了丰富的可视化服务。 展开更多
关键词 计算机应用 图形可视化 非均匀FDTD网格 网格图形 圆柱坐标
下载PDF
栅条圆筒筛内物料轴向运动分析 被引量:3
16
作者 连政国 张岩 姜学东 《莱阳农学院学报》 1996年第2期148-152,共5页
本文分析了栅条圆筒筛分选物料时筛筒内物料轴向运动的基本规律,提出了筛筒内物料轴向移动的计算方法。
关键词 栅条圆筒筛 筛筒 轴向位移
下载PDF
农用高压电网圆柱形殷钢电容分压器边缘效应的研究 被引量:1
17
作者 董赫 李伟凯 《黑龙江八一农垦大学学报》 2015年第2期92-96,117,共6页
通过殷钢电容分压器理论计算值与实际测量值的对比,得到边缘效应的存在导致电容分压器的实际测量值偏大的结论。以理论分析为切入点,结合微积分原理及电容串联思想,推导出计及边缘效应的分压器电容值计算公式。针对圆柱形电容器的特点,... 通过殷钢电容分压器理论计算值与实际测量值的对比,得到边缘效应的存在导致电容分压器的实际测量值偏大的结论。以理论分析为切入点,结合微积分原理及电容串联思想,推导出计及边缘效应的分压器电容值计算公式。针对圆柱形电容器的特点,提出利用增加电容分压器边缘电容环倒角半径的方法来削弱边缘效应,并通过试验验证了该方案的可行性。 展开更多
关键词 农村电网 圆柱形殷钢电容分压器 边缘效应 圆环倒角
下载PDF
单层柱面网壳的自振特性和地震响应分析 被引量:2
18
作者 黎大胜 朱召泉 《徐州建筑职业技术学院学报》 2002年第1期26-29,共4页
分析了网壳跨度、长度、矢跨比、支座刚度和网格尺寸对单斜杆型单层柱面网壳结构自振特性的影响 ,计算了这种网壳在不同方向地震作用下的动力响应 .
关键词 单层柱面网壳 自振特性 地震响应 矢跨比 梁单元 振型分解反应谱法
下载PDF
辐射热传导问题的SFVE格式与数值仿真
19
作者 聂存云 舒适 +1 位作者 杭旭登 成娟 《系统仿真学报》 CAS CSCD 北大核心 2012年第2期275-283,共9页
在柱坐标系下,首先为多介质的输运管单温辐射热传导和多介质的三温辐射热传导两类模型问题,构造了一种针对非混合网格的保对称有限体元(SFVE-I)格式;接着针对一种特殊混合网格,分别为定常扩散和三温辐射热传导模型问题,构造了一种有限体... 在柱坐标系下,首先为多介质的输运管单温辐射热传导和多介质的三温辐射热传导两类模型问题,构造了一种针对非混合网格的保对称有限体元(SFVE-I)格式;接着针对一种特殊混合网格,分别为定常扩散和三温辐射热传导模型问题,构造了一种有限体元(FVE-II)格式;进行了数值实验与数值仿真,其结果表明:对定常问题,在一致网格上有限体元解函数关于真解函数和流函数的逼近分别具有二阶、一阶精度;对上述两种模型问题,数值仿真得到与实际物理现象基本相符的热传导过程和温度分布,且能量守恒误差均小于5%。 展开更多
关键词 柱坐标 辐射热传导 保对称有限体元格式 混合网格 数值仿真
下载PDF
乏油工况下圆柱滚子线接触弹流润滑分析 被引量:1
20
作者 王靖岳 田聪 +1 位作者 武旭东 张硕 《沈阳理工大学学报》 CAS 2023年第6期75-82,共8页
建立了圆柱滚子在乏油工况下的线接触弹流润滑模型,并运用多重网格法计算得到线接触弹流润滑的压力和膜厚分布,研究乏油工况下供油膜厚、黏度等参数变化对弹流润滑特性的影响。模拟结果表明:增加供油膜厚,弹流润滑的膜厚增大且伴有颈缩... 建立了圆柱滚子在乏油工况下的线接触弹流润滑模型,并运用多重网格法计算得到线接触弹流润滑的压力和膜厚分布,研究乏油工况下供油膜厚、黏度等参数变化对弹流润滑特性的影响。模拟结果表明:增加供油膜厚,弹流润滑的膜厚增大且伴有颈缩现象,润滑压力在接触区近似为Hertz分布;黏度对压力影响较小,但对油膜厚度影响较大,随着黏度增加,油膜厚度减小;随着乏油情况的改善,接触区润滑膜厚增大,压力出现二次峰值,但供油量达到一定程度后润滑效果不再增强。 展开更多
关键词 乏油工况 圆柱滚子 弹流润滑 多重网格法
下载PDF
上一页 1 2 3 下一页 到第
使用帮助 返回顶部