期刊文献+
共找到166篇文章
< 1 2 9 >
每页显示 20 50 100
基于Boltzmann时间函数的地表任意点沉陷动态预计 被引量:3
1
作者 徐良骥 孙志豪 +2 位作者 刘潇鹏 张坤 曹宗友 《煤炭学报》 北大核心 2025年第2期715-731,共17页
煤炭资源的地下开采会造成地表沉陷,对地表生态环境及建(构)筑物的安全使用造成一定威胁,在采前获知开采过程中的地表沉陷动态预计值,是进行开采沉陷区生态环境动态修复设计的重要基础,是该领域亟待解决的问题之一。为实现煤炭地下开采... 煤炭资源的地下开采会造成地表沉陷,对地表生态环境及建(构)筑物的安全使用造成一定威胁,在采前获知开采过程中的地表沉陷动态预计值,是进行开采沉陷区生态环境动态修复设计的重要基础,是该领域亟待解决的问题之一。为实现煤炭地下开采导致地表沉陷动态过程的准确预计,根据地表沉陷动态规律总结出理想时间函数模型的形态,据此引入Boltzmann时间函数模型,从下沉量、下沉速度、下沉加速度3个方面对该模型进行分析,发现其能够满足地表沉陷动态趋势;探究该时间函数模型各参数对模型图像的影响,确定其物理意义并分别定义为最终下沉量A、最大下沉速度出现时间t0、下沉急缓程度系数B,从而构建基于Boltzmann时间函数的动态预计模型参数体系;通过对单点实测下沉量进行拟合发现该模型拟合精度较传统动态预计模型更高,拟合优度R^(2)达到0.9988;对矿区地表监测点实测下沉量进行参数反演,根据反演结果建立了沉陷盆地内任意点动态预计参数与地表最大下沉量、回采速度及覆岩岩性系数的相关关系,给出了该模型各动态预计参数在地表任意点的计算方法,并利用收集的6个工作面数据验证其精度可靠;构建了融合Boltzmann时间函数与概率积分法的地表沉陷动态预计模型,可实现对沉陷盆地内任意点任意时间的地表沉陷预计;运用该预计模型求得多时期下沉量并对其进行精度验证,结果显示开采过程中的动态预计相对误差保持在6.0%以内,相对误差最小值为2.7%。 展开更多
关键词 开采沉陷 动态预计 boltzmann时间函数 概率积分模型 下沉速度
在线阅读 下载PDF
Solution of the Poisson-Boltzmann Equation for a Cylindrical Particle with a Limited Length: Functional Theoretical Approach
2
作者 WANG Zheng-wu GU Ming-yan +1 位作者 ZHANG Ge-xin YI Xi-zhang 《Chemical Research in Chinese Universities》 SCIE CAS CSCD 2005年第5期592-596,共5页
With the help of the method of separation of variables and the Debye-Hüchel approximation, the Poisson-Boltzmann equation that describes the distribution of the potential in the electrical double layer of a cylin... With the help of the method of separation of variables and the Debye-Hüchel approximation, the Poisson-Boltzmann equation that describes the distribution of the potential in the electrical double layer of a cylindrical particle with a limited length has been firstly solved under a very low potential condition. Then with the help of the functional analysis theory this equation has been further analytically solved under general potential conditions and consequently, the corresponding surface charge densities have been obtained. Both the potential and the surface charge densities cointide with those results obtained from the Debye-Hüchel approximation when the very low potential of zeψ〈〈kT is introduced. 展开更多
关键词 Cylindrical particle Electrical double layer Poisson-boltzmann equation Surface charge density functional analysis theory
在线阅读 下载PDF
基于对流扩散系统的不可压Navier-Stokes方程的多分布正则化格子Boltzmann方法
3
作者 赵勇 葛懿萱 +2 位作者 陈鑫梦 陈震宇 汪垒 《力学学报》 北大核心 2025年第7期1597-1610,共14页
不可压Navier-Stokes方程组在环境科学、生物医学和流体力学等众多科学领域都扮演着十分重要的角色.发展求解不可压Navier-Stokes方程组的稳定高效数值模拟方法具有重要的科学价值和实际意义.为此,提出了一种求解不可压Navier-Stokes方... 不可压Navier-Stokes方程组在环境科学、生物医学和流体力学等众多科学领域都扮演着十分重要的角色.发展求解不可压Navier-Stokes方程组的稳定高效数值模拟方法具有重要的科学价值和实际意义.为此,提出了一种求解不可压Navier-Stokes方程组的多分布正则化格子Boltzmann(MDF-RLB)模型.该模型的核心思想是将不可压Navier-Stokes方程组转换为一个耦合的对流扩散系统,并针对该系统进行正则化格子Boltzmann方法的建模,即为系统中的每一个对流扩散方程(CDE)构造一个分布函数的演化方程.接着,通过Chapman-Enskog分析证明了该模型能够准确恢复基于对流扩散系统的不可压Navier-Stokes方程组.此外,本文推导了利用分布函数的零阶矩和一阶矩直接计算速度和压力,以及利用非平衡态分布函数的一阶矩来局部计算速度梯度、速度散度、应变率张量、切应力和涡度的公式.最后,通过一系列的基准解算例:二维泊肃叶流、简化二维四辊轧机问题以及二维顶盖驱动方腔流验证了本模型和这些物理量的局部计算格式的有效性和准确性.通过数值测试发现本模型和我们所提出的一些物理量的局部格式在空间上具有二阶收敛速率.同时,与多分布多松弛的格子Boltzmann(MDF-MRTLB)模型相比,MDF-RLB模型在某些情况下模拟结果更精确,且MDF-RLB模型具有更高的计算效率,计算时间减少了7%以上. 展开更多
关键词 格子boltzmann方法 不可压Navier-Stokes方程组 对流扩散系统 多重分布函数 正则化
在线阅读 下载PDF
典型PBX基于Boltzmann函数的准静态单轴拉压非线性本构模型(英文) 被引量:3
4
作者 唐维 颜熹琳 +4 位作者 温茂萍 赵龙 李明 刘彤 章定国 《含能材料》 EI CAS CSCD 北大核心 2017年第8期689-693,共5页
准静态载荷下的拉伸压缩试验中,几乎所有高聚物粘结炸药(PBX)都体现出了明显的非线性本构行为,现常用的本构模型在描述这种非线性时适应性不甚理想。针对五种典型PBX炸药(PBX-X,PBX-9502,LX-17,PBX-9501,EDC-37)应力应变曲线的非线性及... 准静态载荷下的拉伸压缩试验中,几乎所有高聚物粘结炸药(PBX)都体现出了明显的非线性本构行为,现常用的本构模型在描述这种非线性时适应性不甚理想。针对五种典型PBX炸药(PBX-X,PBX-9502,LX-17,PBX-9501,EDC-37)应力应变曲线的非线性及非对称性,基于Boltzmann函数推导了一种拉伸和压缩一起考虑的准静态单轴非线性四参数本构模型。在讨论模型待定参数物理意义的基础上,基于文中提出的待定参数初值确定方法,采用该本构模型对五种典型PBX的应力应变曲线进行了拟合。结果表明,该本构模型可较好地描述五种典型PBX炸药不同温度和应变率下的本构行为。基于Boltzmann函数的拉压非线性本构模型,有望作为一种适用于不同PBX的通用准静态本构模型广泛运用。 展开更多
关键词 高聚物粘结炸药(PBX) 应力应变曲线 本构模型 boltzmann函数
在线阅读 下载PDF
区域干旱预警标准研究Ⅰ:基于Boltzmann函数的区域旱涝指数研究 被引量:2
5
作者 侯威 杨杰 封国林 《高原气象》 CSCD 北大核心 2014年第2期432-442,共11页
基于区域内单站的逐日气象旱涝指数,通过具有明确物理意义的Boltzmann函数来求解某一区域在某一时间尺度下位于各个气象旱涝等级的逐日站数百分比累计概率分布,构建了可以反映区域气象旱涝强度的标准化区域气象旱涝强度指数(Regional me... 基于区域内单站的逐日气象旱涝指数,通过具有明确物理意义的Boltzmann函数来求解某一区域在某一时间尺度下位于各个气象旱涝等级的逐日站数百分比累计概率分布,构建了可以反映区域气象旱涝强度的标准化区域气象旱涝强度指数(Regional meteorological Drought Intensity Index,RDI)和反映区域气象旱涝空间分布状态的标准化区域气象旱涝空间分布差异指数(Regional meteorological Drought variance Index,RDVI)。以构建河南省和江苏省月尺度的RDI和RDVI指数为例说明了区域气象旱涝指数的构建方法,并验证了两个指数的有效性。对于任一区域的不同空间尺度,两个指数可以从日尺度到年尺度乃至更长的时间尺度上进行计算。通过Boltzmann函数来构建RDI和RDVI指数的方法还具有很好的拓展应用性,其他气象干旱指数也可以采用本文方法来构建区域旱涝指数。 展开更多
关键词 区域气象旱涝强度指数 区域气象旱涝空间分布差异指数 boltzmann函数
在线阅读 下载PDF
高速压制成形粉末流动过程的格子Boltzmann方法数值模拟 被引量:4
6
作者 李俏杰 郑洲顺 +1 位作者 王爽 刘建康 《中国有色金属学报》 EI CAS CSCD 北大核心 2012年第6期1754-1762,共9页
根据高速压制工艺过程的特点,底部和侧面边界采用反弹与反射的混合格式,顶部通过改变分布函数构造能模拟动态压制过程的边界格式。利用粉体的初始密度和压制最大理想密度建立一个动态的松弛函数来代替格子Boltzmann方法的固定松弛时间... 根据高速压制工艺过程的特点,底部和侧面边界采用反弹与反射的混合格式,顶部通过改变分布函数构造能模拟动态压制过程的边界格式。利用粉体的初始密度和压制最大理想密度建立一个动态的松弛函数来代替格子Boltzmann方法的固定松弛时间。对圆柱体及单级台阶的圆柱体模具的高速压制成形过程进行二维数值模拟,直观显示压制过程中粉末的流动以及应力波的形成与传播过程。模拟结果表明:模具形状对应力波的波形有较大影响;压制完成后,圆柱体压坯最大密度差为12 mg/cm3,密度分布均匀且左右对称,压坯顶部和底部密度较中间部分密度稍大,与实验结果基本相似;圆柱体压坯中存在密度较低的横向窄带区域,揭示成形坯有可能形成分层或断裂。 展开更多
关键词 高速压制 格子boltzmann方法 松弛函数 应力波
在线阅读 下载PDF
格子Boltzmann模型的改进与流体力学方程 被引量:2
7
作者 冯士德 赵 颖 +1 位作者 季仲贞 郜宪林 《力学学报》 EI CSCD 北大核心 2002年第4期492-500,共9页
流体的流动可以看成是分子以上水平的粒子基本运动组合而成,任何一个粒子系统的Hamiltonian都是由动能和势能这两部分所组成.借助于Hamiltonian建立了微观粒子和宏观流体之间的能量守恒准则,发展了一个适合于热流场数值模拟的格子Boltzm... 流体的流动可以看成是分子以上水平的粒子基本运动组合而成,任何一个粒子系统的Hamiltonian都是由动能和势能这两部分所组成.借助于Hamiltonian建立了微观粒子和宏观流体之间的能量守恒准则,发展了一个适合于热流场数值模拟的格子Boltzmann模型.从该模型可以还原出宏观的流体力学方程,所得动量方程的黏性输运项除了具有Navier-Stokes黏性力的特征外还与非定常的、非线性的动量通量和非定常的内能相关.用该模型对Benard热对流进行了数值模拟,很好地再现了Benard cell,并且克服了热格子Boltzmann模型数值稳定性差的不足. 展开更多
关键词 格子boltzmann模型 HAMILTONIAN 分布函数 流体力学方程 Benard热对流 数值模拟
在线阅读 下载PDF
一个完全守恒的格子Boltzmann模型 被引量:1
8
作者 冯士德 任荣彩 毛江玉 《计算物理》 CSCD 北大核心 2001年第5期407-411,共5页
利用六角形格子离散的方法 ,使得每个六角形Cell里含有 3种不同运动速度的粒子微团 ,宏观物理量用这些粒子微团的矩来定义 .根据微观和宏观之间的质量、动量、能量守恒准则 ,建立了一个二维的D2Q19格子Bo ltzmann模型 ,从该D2Q19模型出... 利用六角形格子离散的方法 ,使得每个六角形Cell里含有 3种不同运动速度的粒子微团 ,宏观物理量用这些粒子微团的矩来定义 .根据微观和宏观之间的质量、动量、能量守恒准则 ,建立了一个二维的D2Q19格子Bo ltzmann模型 ,从该D2Q19模型出发可推导出宏观的流体力学方程组 .用该模型对冲击波在障碍物组表面上的折射现象进行数值模拟 ,并将计算结果与实验进行了比较 . 展开更多
关键词 boltzmann模型 分布函数 冲击波 流体力学方程 六角形格子离散法 完全守恒
在线阅读 下载PDF
Euler方程的双分布函数格子Boltzmann Godunov方法 被引量:2
9
作者 阎广武 胡守信 施卫平 《计算物理》 CSCD 北大核心 1998年第2期67-72,共6页
应用双分布函数系统,通过Godunov分解,构造了一维Euler方程的格子Boltzmann算法。解决了传统格子气固有的GC问题与能量方程之间的矛盾,实现了分布函数与宏观物理量之间的一一对应。
关键词 EULER方程 格子boltzmann方法 双分布函数 Godunov分解
在线阅读 下载PDF
基于格子Boltzmann方法的方腔内自然对流与换热的数值模拟 被引量:1
10
作者 姚寿广 王公利 +1 位作者 程清芳 周长江 《江苏科技大学学报(自然科学版)》 CAS 2013年第5期478-482,共5页
基于双分布格子Boltzmann模型,建立了适合于流体流动和换热的热格子Boltzmann模型.温度分布函数中采用D2Q9离散速度模型.以热格子Boltzmann模型,模拟了方腔内自然对流的形成及其演化,通过与相关文献的计算结果对比可以发现,热格子Boltzm... 基于双分布格子Boltzmann模型,建立了适合于流体流动和换热的热格子Boltzmann模型.温度分布函数中采用D2Q9离散速度模型.以热格子Boltzmann模型,模拟了方腔内自然对流的形成及其演化,通过与相关文献的计算结果对比可以发现,热格子Boltzmann模型在处理流体流动与传热方面存在着独特的优点,文中建立的数值模拟计算方法和程序是切实有效的. 展开更多
关键词 热格子boltzmann 双分布函数模型 流动换热
在线阅读 下载PDF
参变极值问题的信息凝聚分布与Boltzmann极大熵函数 被引量:1
11
作者 刘法胜 傅白白 夏尊铨 《数学物理学报(A辑)》 CSCD 北大核心 1999年第S1期531-536,共6页
该文利用Boltzmann 熵概念给出了参变极值问题最优解的一种积分极限表达式和极值函数的极大熵函数,讨论了它们一致收敛性的要求并给出了极大熵函数一致收敛的一个充分条件,将之应用到全局最优解问题得到了全局最优解和最优... 该文利用Boltzmann 熵概念给出了参变极值问题最优解的一种积分极限表达式和极值函数的极大熵函数,讨论了它们一致收敛性的要求并给出了极大熵函数一致收敛的一个充分条件,将之应用到全局最优解问题得到了全局最优解和最优值的一种显表示,最后还探讨了极大熵函数在一类双层规划问题求解中的应用. 展开更多
关键词 极大极小问题 boltzmann 极大熵函数 全局最优解 双层规划
在线阅读 下载PDF
格子Boltzmann方法求解一维黏性Burgers方程的适用性研究 被引量:3
12
作者 张伟 李小宝 陶建华 《天津城市建设学院学报》 CAS 2006年第3期222-226,共5页
构造了一维黏性Burgers方程的一个基于3-速格子模型的格子Boltzmann格式,分析了格式的单调性和稳定性,得到了格式的单调性条件,并证明了在此条件下,格式是L_∞稳定的;分析了平衡分布函数中含有的可调系数对格式性能的影响,研究表明:... 构造了一维黏性Burgers方程的一个基于3-速格子模型的格子Boltzmann格式,分析了格式的单调性和稳定性,得到了格式的单调性条件,并证明了在此条件下,格式是L_∞稳定的;分析了平衡分布函数中含有的可调系数对格式性能的影响,研究表明:通过引入适当形式的平衡分布函数,采用简单的3-速格子,可以恢复宏观Burgers方程,通过与标准二阶精度的有限差分解的比较,证明了本文的LB格式是简单而有效的。 展开更多
关键词 格子boltzmann方法 BURGERS方程 平衡分布函数 单调性 稳定性
在线阅读 下载PDF
幂律流体扩展流动的格子Boltzmann模拟分析 被引量:2
13
作者 何录武 薛盼 李勇 《兰州大学学报(自然科学版)》 CAS CSCD 北大核心 2016年第5期708-712,共5页
应用格子Boltzmann方法和D2Q9格子模型,结合剪切速率的局部计算方法以及非平衡外推边界处理格式,数值模拟了幂律流体在3:1扩展流道(2D)中的流动行为.获得了不同雷诺数Re和幂律指数n下扩展流动的流线分布,计算出了漩涡的涡心位置和大小,... 应用格子Boltzmann方法和D2Q9格子模型,结合剪切速率的局部计算方法以及非平衡外推边界处理格式,数值模拟了幂律流体在3:1扩展流道(2D)中的流动行为.获得了不同雷诺数Re和幂律指数n下扩展流动的流线分布,计算出了漩涡的涡心位置和大小,分析了Re和n对扩展流动特点的影响.模拟结果表明本模型和处理方法具有良好的精度和稳定性. 展开更多
关键词 格子boltzmann方法 幂律流体 分布函数 扩展流 雷诺数
原文传递
偏微分方程求解的一种新颖方法——格子Boltzmann模型 被引量:2
14
作者 乐励华 高云 刘唐伟 《大学数学》 2011年第3期75-82,共8页
介绍了一种偏微分方程求解的一种新颖方法格子Boltzmann模型,详细分析了它的基本理论和基本原理。并通过不可压Navier-Stokes方程组和二维含源项扩散方程的数值模拟计算实例,说明格子Boltzmann方法的有效性,展示了广阔的应用前景,为今... 介绍了一种偏微分方程求解的一种新颖方法格子Boltzmann模型,详细分析了它的基本理论和基本原理。并通过不可压Navier-Stokes方程组和二维含源项扩散方程的数值模拟计算实例,说明格子Boltzmann方法的有效性,展示了广阔的应用前景,为今后更深入的研究和广泛应用提供参考. 展开更多
关键词 格子boltzmann方法 平衡态分布函数 D2Q9模型 NAVIER-STOKES方程 对流—扩散方程
在线阅读 下载PDF
一维粘性泥石流运动的格子Boltzmann模拟 被引量:3
15
作者 王沁 姚令侃 《灾害学》 CSCD 2007年第4期1-5,共5页
利用格子Boltzmann方法,对一维粘性泥石流的运动进行了定量化研究和数值模拟,为进一步研究泥石流的运动机制奠定了基础。在模拟一维粘性泥石流运动时,采用了含外力项的格子Boltzmann模型,通过对外力因子、平衡态分布函数的构造,使得宏... 利用格子Boltzmann方法,对一维粘性泥石流的运动进行了定量化研究和数值模拟,为进一步研究泥石流的运动机制奠定了基础。在模拟一维粘性泥石流运动时,采用了含外力项的格子Boltzmann模型,通过对外力因子、平衡态分布函数的构造,使得宏观统计特性刻画了泥石流运动的非线性特征,成功模拟了泥石流的运动过程,获得了泥石流水深、流速两大要素的时空变化规律,为研究泥石流运动过程提供了一种新手段。 展开更多
关键词 格子boltzmann方法 泥石流 平衡态分布函数 外力因子
在线阅读 下载PDF
非牛顿流体流动的格子Boltzmann方法研究进展 被引量:8
16
作者 李勇 柳文琴 《力学与实践》 北大核心 2014年第4期383-395,共13页
格子玻尔兹曼方法 (lattice Boltzmann method,LBM)能够直接计算局部剪切速率并可以达到二次精度,因此在非牛顿流动数值模拟中展现出一定优势.尽管已证实LBM对于非牛顿流动的适用性,但是LBM需要通过即时调节BGK(Bhatnagar--Gross--Krook... 格子玻尔兹曼方法 (lattice Boltzmann method,LBM)能够直接计算局部剪切速率并可以达到二次精度,因此在非牛顿流动数值模拟中展现出一定优势.尽管已证实LBM对于非牛顿流动的适用性,但是LBM需要通过即时调节BGK(Bhatnagar--Gross--Krook)碰撞项中的松弛时间来实时反映黏度改变,当松弛时间接近1/2时,迭代会出现数值不稳定现象.该文对LBM在非牛顿流体研究中的进展进行了总结,介绍了增加数值稳定性的方法并对结果的精度进行了比较,在此基础上对LBM在非牛顿研究中的进一步发展进行了展望. 展开更多
关键词 格子玻尔兹曼方法 广义牛顿流体 黏弹性流体 分布函数
在线阅读 下载PDF
劈裂注浆桩单桩承载力试验及模型预测研究
17
作者 冯杰 《建材技术与应用》 2026年第1期5-11,共7页
设计了不同桩距的劈裂注浆现场原位试验和单桩竖向静载荷试验,得到各试桩的荷载沉降数据。采用Boltzmann函数、指数函数和幂函数模型对试桩的荷载沉降数据进行曲线拟合,并进行极限承载力预测。通过回归拟合,得到最优的劈裂注浆桩单桩荷... 设计了不同桩距的劈裂注浆现场原位试验和单桩竖向静载荷试验,得到各试桩的荷载沉降数据。采用Boltzmann函数、指数函数和幂函数模型对试桩的荷载沉降数据进行曲线拟合,并进行极限承载力预测。通过回归拟合,得到最优的劈裂注浆桩单桩荷载-沉降函数模型,并得出单桩极限承载力预测公式。 展开更多
关键词 劈裂注浆桩 非线性拟合 boltzmann函数 极限承载力试验
在线阅读 下载PDF
低浓度固液两相流Boltzmann方程的同伦分析方法解 被引量:2
18
作者 张丽 王光谦 +1 位作者 傅旭东 孙其诚 《应用基础与工程科学学报》 EI CSCD 2009年第6期811-818,共8页
同伦分析方法(Homotopy analysis method,HAM)是求解强非线性问题的有力手段.针对两相流动理学理论中的非线性微分积分方程——Boltzmann方程,本文采用HAM方法选取Maxwell速度分布函数作为初始猜测解,求解得到了低浓度固液两相流的BGK模... 同伦分析方法(Homotopy analysis method,HAM)是求解强非线性问题的有力手段.针对两相流动理学理论中的非线性微分积分方程——Boltzmann方程,本文采用HAM方法选取Maxwell速度分布函数作为初始猜测解,求解得到了低浓度固液两相流的BGK模型Boltzmann方程的一阶近似解,与传统的Chapman-Enskog方法得到的一阶近似解表达式的结构一致,显示了HAM方法求解Boltzmann方程的有效性,为一般Boltzmann方程的HAM方法求解奠定了基础. 展开更多
关键词 固液两相流 boltzmann方程 同伦方法 Maxwell速度分布函数
在线阅读 下载PDF
BOLTZMANN碰撞矩阵元的计算
19
作者 向安平 朱世德 +1 位作者 任基 谢文 《计算物理》 CSCD 北大核心 1996年第3期366-372,共7页
用Burnet函数作为展开基底表示分子速度分布函数和Boltzmann方程,并由此证明,Boltzmann碰撞矩阵元的计算归结为约化积分、ClebschGordan系数和MrodyMoshinsky系数的计算。最后。
关键词 碰撞矩阵元 波尔兹曼矩阵 分子输运
在线阅读 下载PDF
交通网络β指数增长曲线的Boltzmann模型
20
作者 陈彦光 《北京大学学报(自然科学版)》 EI CAS CSCD 北大核心 2023年第2期326-330,共5页
从已被验证的城市化水平和交通网络的经验模型出发,提出交通网络β指数的数学模型。基于城市化水平的一次logistic函数,借助数学演绎法,推导出交通网络β指数增长曲线的Boltzmann方程;基于城市化水平的二次logistic函数,推导出交通网络... 从已被验证的城市化水平和交通网络的经验模型出发,提出交通网络β指数的数学模型。基于城市化水平的一次logistic函数,借助数学演绎法,推导出交通网络β指数增长曲线的Boltzmann方程;基于城市化水平的二次logistic函数,推导出交通网络β指数增长曲线的二次Boltzmann方程。前者适用于欧美国家和中国东南沿海地区,后者则适用于中国内地,特别是中国北方地区。该模型具有两个主要功能:1)解释和预测交通网络的发展;2)有助于交通网络增长过程的分段研究。 展开更多
关键词 交通网络 城市化 boltzmann方程 logistic函数 异速标度 人均收入
在线阅读 下载PDF
上一页 1 2 9 下一页 到第
使用帮助 返回顶部