力学学报 https://lxxb.www.cn100led.com/ Chinese Journal of Theoretical and Applied Mechanics day 15 2024-02-25 23:43:06 lxxb@www.cn100led.com lxxb@www.cn100led.com 2024-02-25 23:43:06 zh lxxb@www.cn100led.com 0459-1879 <![CDATA[圆柱尾流的绝对不稳定性]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1999-001 <br/><p>在水槽和低湍流度水洞中进行亚临界雷诺数圆柱尾流稳定性实验来流速度由零缓慢增长到一定值后保持不变,稳定足够长时间后,在流向某站位处给流场一个有限幅值的脉冲扰动,测量扰动前后相当长时间内下游尾流速度信号的变化情况当雷诺数处于高亚临界值时,未受扰动的尾流速度脉动很小,处于定常状态,但对近尾流进行脉冲扰动后,能够激发出不衰减的旋涡脱落发现扰动位置限制在圆柱后一定范围内才能有效,再往下游则扰动随时间衰减.说明圆柱近尾流中存在一个绝对不稳定区,在该区域内的扰动将在当地放大,经过复杂的演化,最后形成不衰减的旋涡脱落.</p> <br/>力学学报. 1999 31(1): 1-9. <br/><p>在水槽和低湍流度水洞中进行亚临界雷诺数圆柱尾流稳定性实验来流速度由零缓慢增长到一定值后保持不变,稳定足够长时间后,在流向某站位处给流场一个有限幅值的脉冲扰动,测量扰动前后相当长时间内下游尾流速度信号的变化情况当雷诺数处于高亚临界值时,未受扰动的尾流速度脉动很小,处于定常状态,但对近尾流进行脉冲扰动后,能够激发出不衰减的旋涡脱落发现扰动位置限制在圆柱后一定范围内才能有效,再往下游则扰动随时间衰减.说明圆柱近尾流中存在一个绝对不稳定区,在该区域内的扰动将在当地放大,经过复杂的演化,最后形成不衰减的旋涡脱落.</p> <br/>力学学报. 1999 31(1): 1-9. 圆柱尾流的绝对不稳定性 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 1-9. article doi:10.6052/0459-1879-1999-1-1999-001 10.6052/0459-1879-1999-1-1999-001 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1999-001 1-9 <![CDATA[非定常流函数涡量方程的一种数值解法的研究]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1999-002 <br/><p>对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解.以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数(Ra)条件下的数值研究.结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟.</p> <br/>力学学报. 1999 31(1): 10-20. <br/><p>对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解.以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数(Ra)条件下的数值研究.结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟.</p> <br/>力学学报. 1999 31(1): 10-20. 非定常流函数涡量方程的一种数值解法的研究 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 10-20. article doi:10.6052/0459-1879-1999-1-1999-002 10.6052/0459-1879-1999-1-1999-002 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1999-002 10-20 <![CDATA[亚临界雷诺数下圆柱绕流的大涡模拟]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-010 <br/><p>本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟.求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算.计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的</p> <br/>力学学报. 1999 31(1): 100-105. <br/><p>本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟.求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算.计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的</p> <br/>力学学报. 1999 31(1): 100-105. 亚临界雷诺数下圆柱绕流的大涡模拟 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 100-105. article doi:10.6052/0459-1879-1999-1-1995-010 10.6052/0459-1879-1999-1-1995-010 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-010 100-105 <![CDATA[二级加载下疲劳可靠性分析]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-011 <br/><p>基于二维概率Miner准则,建立了在二级变幅加载下进行疲劳可靠性分析的系统方法并采用八组较大子样的二级低高和高低加载下的试验数据进行了验证、结果表明;本文方法初步实现了对第一级加载后的剩余疲劳寿命和二级加载下的疲劳总寿命进行较为准确的可靠性预测,预测值与试验值吻合良好</p> <br/>力学学报. 1999 31(1): 106-112. <br/><p>基于二维概率Miner准则,建立了在二级变幅加载下进行疲劳可靠性分析的系统方法并采用八组较大子样的二级低高和高低加载下的试验数据进行了验证、结果表明;本文方法初步实现了对第一级加载后的剩余疲劳寿命和二级加载下的疲劳总寿命进行较为准确的可靠性预测,预测值与试验值吻合良好</p> <br/>力学学报. 1999 31(1): 106-112. 二级加载下疲劳可靠性分析 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 106-112. article doi:10.6052/0459-1879-1999-1-1995-011 10.6052/0459-1879-1999-1-1995-011 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-011 106-112 <![CDATA[蜂窝芯层的等效弹性参数]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-012 <br/><p>反映蜂窝材料等效弹性参数的Gibson公式结果简单,便于应用然而由于对应于蜂窝壁板伸缩变形的刚度被忽略,导致蜂窝夹层结构数值分析时,芯层材料的弹性矩阵表现出不确定性本文重新考虑了蜂窝壁板的伸缩变形对面内刚度的影响,对Gibson公式进行了修正本本文结果不但克服了Gibson公式的缺陷,同时提出了考虑蜂窝芯层面内刚度的一种简化方案,该方案可以方便地应用于蜂窝夹层结构的计算。</p> <br/>力学学报. 1999 31(1): 113-118. <br/><p>反映蜂窝材料等效弹性参数的Gibson公式结果简单,便于应用然而由于对应于蜂窝壁板伸缩变形的刚度被忽略,导致蜂窝夹层结构数值分析时,芯层材料的弹性矩阵表现出不确定性本文重新考虑了蜂窝壁板的伸缩变形对面内刚度的影响,对Gibson公式进行了修正本本文结果不但克服了Gibson公式的缺陷,同时提出了考虑蜂窝芯层面内刚度的一种简化方案,该方案可以方便地应用于蜂窝夹层结构的计算。</p> <br/>力学学报. 1999 31(1): 113-118. 蜂窝芯层的等效弹性参数 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 113-118. article doi:10.6052/0459-1879-1999-1-1995-012 10.6052/0459-1879-1999-1-1995-012 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-012 113-118 <![CDATA[多机械臂搬运同一物体的协调动态载荷分配]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-013 <br/><p>针对多机械臂共同搬运同一物体形成闭运动链的协调系统,研究了多机械臂协调动态载荷分配存在冗余情况下的实时性分配方法文中提出了“应集中惯性质量棒”的概念,进而根据各机械臂的承载能力,采用线性加权的方法发展了一种操作物体质心处零内力的载荷分配原则。最后给出了系统载荷分配的解析表达式,有效地解决了多机械臂搬运同一物体时其动态载荷分配的实时性问题。</p> <br/>力学学报. 1999 31(1): 119-125. <br/><p>针对多机械臂共同搬运同一物体形成闭运动链的协调系统,研究了多机械臂协调动态载荷分配存在冗余情况下的实时性分配方法文中提出了“应集中惯性质量棒”的概念,进而根据各机械臂的承载能力,采用线性加权的方法发展了一种操作物体质心处零内力的载荷分配原则。最后给出了系统载荷分配的解析表达式,有效地解决了多机械臂搬运同一物体时其动态载荷分配的实时性问题。</p> <br/>力学学报. 1999 31(1): 119-125. 多机械臂搬运同一物体的协调动态载荷分配 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 119-125. article doi:10.6052/0459-1879-1999-1-1995-013 10.6052/0459-1879-1999-1-1995-013 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-013 119-125 <![CDATA[温盐双扩散系统的多平衡稳定性]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-001 <br/><p>采用数值试验的方法研究了温盐双扩散系统中的多平衡稳定性在系统中引入非线性状态方程、温通量边界条件和变强迫因子在数值模拟中引入边界拟会坐标系,并按照温盐浮力比λ大于1,等于1和小于1将流场分为温度控制场、温盐均势场和盐度控制场得到了四种稳定平衡状态,特别是,大强迫因子和强非线性影响下的四胞环流对称状态.在对称的定常强迫边条件作用下的流场、温度场和盐度场均为对称分布引入温通量边条件后,对称流场在小干扰下失去对称性,发展为非对称分布,并因为干扰位置的不同而造成不同的非对称性,出现了从对称状态到非对称状态的音叉分叉,得到了a-Ras的稳定曲线及分叉图,并解释了分叉的本质</p> <br/>力学学报. 1999 31(1): 21-28. <br/><p>采用数值试验的方法研究了温盐双扩散系统中的多平衡稳定性在系统中引入非线性状态方程、温通量边界条件和变强迫因子在数值模拟中引入边界拟会坐标系,并按照温盐浮力比λ大于1,等于1和小于1将流场分为温度控制场、温盐均势场和盐度控制场得到了四种稳定平衡状态,特别是,大强迫因子和强非线性影响下的四胞环流对称状态.在对称的定常强迫边条件作用下的流场、温度场和盐度场均为对称分布引入温通量边条件后,对称流场在小干扰下失去对称性,发展为非对称分布,并因为干扰位置的不同而造成不同的非对称性,出现了从对称状态到非对称状态的音叉分叉,得到了a-Ras的稳定曲线及分叉图,并解释了分叉的本质</p> <br/>力学学报. 1999 31(1): 21-28. 温盐双扩散系统的多平衡稳定性 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 21-28. article doi:10.6052/0459-1879-1999-1-1995-001 10.6052/0459-1879-1999-1-1995-001 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-001 21-28 <![CDATA[垂直湍流液-固流中大颗粒的相对速度]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-002 <br/><p>通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式.用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数.实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。</p> <br/>力学学报. 1999 31(1): 29-37. <br/><p>通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式.用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数.实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。</p> <br/>力学学报. 1999 31(1): 29-37. 垂直湍流液-固流中大颗粒的相对速度 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 29-37. article doi:10.6052/0459-1879-1999-1-1995-002 10.6052/0459-1879-1999-1-1995-002 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-002 29-37 <![CDATA[损伤统计演化方程的性质和数值模拟]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-003 <br/><p>通过对一种含成核尺寸效应的损伤统计演化方程性质的分析和数值模拟,揭示了损伤率主要是由微损伤在二维相空间中的前沿的运动所决定的这也就是Kachanov提出的损伤率演化方程的物理基础数值结果进一步显示了含成核尺寸效应模型在损伤发展上与-维模型的区别而且,由几种形式的细观动力学算出的损伤率与损伤的关系简单,可近似拟合为宏观上封闭的形式</p> <br/>力学学报. 1999 31(1): 38-48. <br/><p>通过对一种含成核尺寸效应的损伤统计演化方程性质的分析和数值模拟,揭示了损伤率主要是由微损伤在二维相空间中的前沿的运动所决定的这也就是Kachanov提出的损伤率演化方程的物理基础数值结果进一步显示了含成核尺寸效应模型在损伤发展上与-维模型的区别而且,由几种形式的细观动力学算出的损伤率与损伤的关系简单,可近似拟合为宏观上封闭的形式</p> <br/>力学学报. 1999 31(1): 38-48. 损伤统计演化方程的性质和数值模拟 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 38-48. article doi:10.6052/0459-1879-1999-1-1995-003 10.6052/0459-1879-1999-1-1995-003 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-003 38-48 <![CDATA[弹性T项对裂尖参数和裂纹扩展路径稳定性的作用]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-004 <br/><p>研究了弹性T项在主裂纹与近裂尖微空洞干涉问题中对裂尖参数和裂纹扩展路径稳定性的影响作用.利用“伪力”法,并通过解决主裂纹与近裂尖微空洞干涉问题,对远场纯1型载荷和远场弹性T项载荷下,该干涉问题中弹性T项的影响作用进行了分析从数值结果可以看出:由于空洞的存在;释放了弹性T项所引起的应力,弹性T项对裂尖参数;应力强度因子和J积分都有直接显著的影响,因而,它对该载荷下的裂纹扩展路径的稳定性有控制作用。</p> <br/>力学学报. 1999 31(1): 49-57. <br/><p>研究了弹性T项在主裂纹与近裂尖微空洞干涉问题中对裂尖参数和裂纹扩展路径稳定性的影响作用.利用“伪力”法,并通过解决主裂纹与近裂尖微空洞干涉问题,对远场纯1型载荷和远场弹性T项载荷下,该干涉问题中弹性T项的影响作用进行了分析从数值结果可以看出:由于空洞的存在;释放了弹性T项所引起的应力,弹性T项对裂尖参数;应力强度因子和J积分都有直接显著的影响,因而,它对该载荷下的裂纹扩展路径的稳定性有控制作用。</p> <br/>力学学报. 1999 31(1): 49-57. 弹性T项对裂尖参数和裂纹扩展路径稳定性的作用 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 49-57. article doi:10.6052/0459-1879-1999-1-1995-004 10.6052/0459-1879-1999-1-1995-004 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-004 49-57 <![CDATA[特征值问题的边界形状灵敏度]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-005 <br/><p>研究连续系统振动特征值问题的边界形状灵敏度满足什么方程和边界条件,如何离散化作近似计算结果表明:如果采用相同的有限单元剖分模式,边界形状灵敏度方程和特征值问题方程具有相同的系数矩阵,但前者是非齐次方程,后者是齐次方程;前者需要施加非齐次边界条件,后者施加齐次边界条件。</p> <br/>力学学报. 1999 31(1): 58-67. <br/><p>研究连续系统振动特征值问题的边界形状灵敏度满足什么方程和边界条件,如何离散化作近似计算结果表明:如果采用相同的有限单元剖分模式,边界形状灵敏度方程和特征值问题方程具有相同的系数矩阵,但前者是非齐次方程,后者是齐次方程;前者需要施加非齐次边界条件,后者施加齐次边界条件。</p> <br/>力学学报. 1999 31(1): 58-67. 特征值问题的边界形状灵敏度 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 58-67. article doi:10.6052/0459-1879-1999-1-1995-005 10.6052/0459-1879-1999-1-1995-005 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-005 58-67 <![CDATA[有限长Timoshenko梁弹性碰撞接触瞬间的动态特性]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-006 <br/><p>给出了质点与有限长Timoshenko梁横向弹性碰撞接触问题的半解析解,分析了该碰撞问题在碰撞接触瞬间的动态响应特性:揭示了其中的波传播现象.</p> <br/>力学学报. 1999 31(1): 68-74. <br/><p>给出了质点与有限长Timoshenko梁横向弹性碰撞接触问题的半解析解,分析了该碰撞问题在碰撞接触瞬间的动态响应特性:揭示了其中的波传播现象.</p> <br/>力学学报. 1999 31(1): 68-74. 有限长Timoshenko梁弹性碰撞接触瞬间的动态特性 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 68-74. article doi:10.6052/0459-1879-1999-1-1995-006 10.6052/0459-1879-1999-1-1995-006 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-006 68-74 <![CDATA[薄圆板的后屈曲响应及圆形脱层的能量释放率]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-007 <br/><p>研究了薄圆板及圆形脱层薄膜的轴对称后屈曲问题,并首次研究了轴对称后屈曲变形之后的二次非轴对称分叉屈曲问题.采用暴级数展开的方法,结合应用傅立叶级数得出非线性及非轴对称问题方程的解井山此计算了圆形脱层薄膜与上述后屈曲变形的两个阶段对应的能量释放率</p> <br/>力学学报. 1999 31(1): 75-83. <br/><p>研究了薄圆板及圆形脱层薄膜的轴对称后屈曲问题,并首次研究了轴对称后屈曲变形之后的二次非轴对称分叉屈曲问题.采用暴级数展开的方法,结合应用傅立叶级数得出非线性及非轴对称问题方程的解井山此计算了圆形脱层薄膜与上述后屈曲变形的两个阶段对应的能量释放率</p> <br/>力学学报. 1999 31(1): 75-83. 薄圆板的后屈曲响应及圆形脱层的能量释放率 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 75-83. article doi:10.6052/0459-1879-1999-1-1995-007 10.6052/0459-1879-1999-1-1995-007 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-007 75-83 <![CDATA[安全壳中钢衬壳热屈曲问题理论与实验研究(Ⅱ)──实验部分]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-008 <br/><p>钢衬壳热屈曲问题是核工程安全壳设计中的主要问题,但实验研究方面的文章发表得不多文中以200兆瓦核电站安全壳中钢衬壳为研究对象,采用局部1:1模型,测得了钢衬壳热屈曲温度和应变载荷,给出了钢衬壳屈曲和初始后屈曲过程中挠度和温度关系、以及膜应变和温度关系,实验测得钢衬壳具有局域屈曲的现象,实验屈曲载荷与理论结果符合较好</p> <br/>力学学报. 1999 31(1): 75-83. <br/><p>钢衬壳热屈曲问题是核工程安全壳设计中的主要问题,但实验研究方面的文章发表得不多文中以200兆瓦核电站安全壳中钢衬壳为研究对象,采用局部1:1模型,测得了钢衬壳热屈曲温度和应变载荷,给出了钢衬壳屈曲和初始后屈曲过程中挠度和温度关系、以及膜应变和温度关系,实验测得钢衬壳具有局域屈曲的现象,实验屈曲载荷与理论结果符合较好</p> <br/>力学学报. 1999 31(1): 75-83. 安全壳中钢衬壳热屈曲问题理论与实验研究(Ⅱ)──实验部分 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 75-83. article doi:10.6052/0459-1879-1999-1-1995-008 10.6052/0459-1879-1999-1-1995-008 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-008 75-83 <![CDATA[Poincare型胞映射分析方法及其应用]]> https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-009 <br/><p>本文用Poincare型胞映射方法对平衡及不平.衡轴承转子非线性动力系统的全局特性进行了分析研究,同时求得了一定状态空间内系统存在的周期解及其在各不同Poincare截面上的吸引域,得到了一些新的现象和规律,并通过对平衡及不平衡轴承转子系统的全局特性异同的比较,说明了要建立既适用于平衡轴承转子系统又适用于不平衡轴承转子系统的非线性稳定性准则应注意的几个问题</p> <br/>力学学报. 1999 31(1): 91-99. <br/><p>本文用Poincare型胞映射方法对平衡及不平.衡轴承转子非线性动力系统的全局特性进行了分析研究,同时求得了一定状态空间内系统存在的周期解及其在各不同Poincare截面上的吸引域,得到了一些新的现象和规律,并通过对平衡及不平衡轴承转子系统的全局特性异同的比较,说明了要建立既适用于平衡轴承转子系统又适用于不平衡轴承转子系统的非线性稳定性准则应注意的几个问题</p> <br/>力学学报. 1999 31(1): 91-99. Poincare型胞映射分析方法及其应用 Personal use only, all commercial or other reuse prohibited . 1999 31(1): 91-99. article doi:10.6052/0459-1879-1999-1-1995-009 10.6052/0459-1879-1999-1-1995-009 力学学报 31 1 https://lxxb.www.cn100led.com/article/doi/10.6052/0459-1879-1999-1-1995-009 91-99
Baidu
map