力学学报

姓名
邮箱
手机号码
标题
留言内容
验证码

冯宝平, 米建春

报道出口条件对圆形湍流射流自保持性影响的实验研究结果. 对来自渐缩和长管两种不同结构喷嘴的射流，在相同雷诺数条件下，沿轴线进行了速度测量; 研究的统计量包括平均速度、湍流强度、高阶矩、能谱和积分尺度. 实验结果表明，渐缩喷嘴射流比长管射流发展得更快、更容易达到自保持状态. 通过对比发现，在两射流的速度(温度)场中，平均速度(温度)、湍流强度、偏斜因子和平坦度因子都存在明显的异同. 同时发现两射流的积分尺度随轴向距离的增加都成线性增长，且在渐缩喷嘴射流中增长得更快. 通过对比两射流的边界层厚度、径向与轴向湍流强度的比值、湍动能能谱图并结合前人的研究结果，对两射流湍流场所表现出的不同的统计学行为给出了合理的解释.

2009, 41(5): 609-617. doi: 10.6052/0459-1879-2009-5-2008-352

PDF

两种浸入式边界方法的比较

谢胜百单鹏

介绍两类不同的浸入式边界方法及其对它的改进. 然后采用均匀矩形交错网格和压力校正投影法，对不可压流场中的二维圆柱绕流进行了数值求解并对比了两类方法的精度.计算分析表明，连续显力法具有构造简单，适用性强的优点. 但离散隐力法在物面边界精度上要优于前者. 改进后，在二阶精度的离散格式下物面边界精度较低的显示力源法的精度有一定提高，同时发现，加密网格以提高数值精度的方法对于连续显力法并不总是有效.而同样格式下，离散隐力法具有更高精度，其中预测-校正离散隐力法可以在此基础上获得更小的计算误差和更快的收敛速度. 数值解与文献已有的数值和实验结果吻合得很好，表明边界算法及其程序是可靠和有效的.

2009, 41(5): 618-627. doi: 10.6052/0459-1879-2009-5-2008-282

PDF

叶片转动对涡轮间隙内部流动影响的数值研究

牛茂升臧述升黄名海

机匣与叶片的相对转动是影响涡轮叶顶间隙流动的重要因素之一. 对LISA1.5级轴流涡轮间隙内部流动的数值计算结果表明：叶片转动对涡轮间隙流动有阻塞作用.叶片静止时，由于阻塞作用消失，导致间隙入口速度增大，间隙流量增加，并且通过间隙的流体全部卷起形成间隙涡. 同时在叶片顶部吸力面侧前缘、中部各形成一个间隙涡，使得间隙流动损失增加. 而且转速下降会加剧动叶出口截面气流过偏/偏转不足现象. 同时叶片静止时，间隙前部各个弦长截面内静压自间隙入口开始一直呈增加趋势，直到叶片尾缘附近截面，间隙截面内静压才逐渐稳定.

2009, 41(5): 628-634. doi: 10.6052/0459-1879-2009-5-2008-062

PDF

基于热格子Boltzmann方法的自然对流数值模拟

解建飞钟诚文张勇尹大川

对Eggels和Somers提出的热格子Boltzmann格式进行了改进. 在不可压缩流动的假设下，提出了一种新的温度平衡分布函数，可以克服压缩性对温度统计的影响，并且相应地修正了统计宏观温度的方法. Eggels和Somers的方法对速度和温度均采用半步长反弹格式边界条件，适合无滑移的速度边界条件.但是对温度采用该边界条件在物理本质上显得不够准确，所以在边界上对二者统一采取算法既简单又容易实现的非平衡态外推格式，同时可以与Boltzmann格式的整体二阶精度保持一致. 最后，利用改进的热格子Boltzmann方法(TLBM)模拟了Ra=10^6和Pr=0.71(空气)的方腔中的自然对流，模拟得到的流动参数与其它数值方法的结果吻合得很好，表明改进的热格子Boltzmann方法可以有效准确地模拟非等温流动.

2009, 41(5): 635-640. doi: 10.6052/0459-1879-2009-5-2008-193

PDF

基于气动力辨识的ASE模型降阶研究

张子健徐敏陈士橹

CFD/CSD耦合计算能够准确预测跨音速段飞行器弹性振动的非定常气动力, 但其带来的巨大计算量及高阶维数不利于气动弹性系统的分析与综合. 针对于此，采用系统辨识及均衡截断技术对高阶气动伺服弹性模型进行降阶处理，并利用所得到的低阶模型进行系统综合：(1) 基于Volterra级数气动力辨识技术，得到非定常气动力的时域降阶模型(ROM)，耦合结构动力学模型及控制机构动力学模型获得气动伺服弹性(ASE)状态空间方程；(2) 利用均衡截段法对时域ASE模型进行进一步降阶，得到能够较真实反映所关心频域内系统响应的低阶ASE模型；(3) 针对建模误差和降阶误差存在造成的系统不确定性问题，结合降阶模型采用混合灵敏度$H_{\infty}$控制方法设计颤振主动抑制鲁棒控制律，保证其作用于真实系统的有效性；对控制器进行均衡阶段降阶并保持其鲁棒性，得到低阶鲁棒的颤振抑制控制器. 最后利用典型的BACT模型进行气动伺服弹性的降阶及主动颤振抑制控制，仿真结果表明，基于ROM建立的低阶气动弹性模型能够较真实地反应系统的颤振特性；而基于截断后的降阶模型所设计的低阶鲁棒控制器能够有效应用于存在不确定性摄动的实际系统，并将系统颤振速度提高36%.

2009, 41(5): 641-650. doi: 10.6052/0459-1879-2009-5-2008-223

PDF

云状空化非定常脱落机理的数值与实验研究

张博王国玉黄彪余志毅

结合数值计算和实验技术研究了云状空化的非定常脱落机理. 实验采用高速录像技术观察了绕Clark-y型水翼云状空化形态随时间的变化. 数值计算采用了汽-液两相的当地均相流模型，湍流封闭采用了一种修正的RNG k-\varepsilon方程，利用商业软件的二次开发技术，引入了一种与空化区域水汽相密度相关的系数，对湍流模型进行了修正. 计算和实验结果均表明：云状空化尾部存在着准周期性的涡状空化团的脱落；局部压强的增加是引起空穴断裂进而脱落的直接原因；压强升高是由于近壁处的反向射流引起的；空穴尾部的大规模的旋涡运动引起的近壁处的逆压梯度是引起反向射流的主要原因.

2009, 41(5): 651-659. doi: 10.6052/0459-1879-2009-5-2008-152

PDF

潜堤对波浪在直立式防浪墙上爬高的影响

张景新刘桦

对直立式防浪墙前潜堤的空间布置对波浪在直墙上爬高的影响进行了研究.建立简化模型，即直立式防浪墙前海床设计为平底，潜堤设计为直立薄板. 改变潜堤高度及其与防浪墙之间距，研究不同的潜堤布置对波浪在直墙上爬高的影响. 针对线性波浪场，利用数值波浪水槽模拟了潜堤作用下直墙上波浪的爬高现象. 同时建立了理论模型，系统地分析了潜堤布置形式对直墙上波浪爬高的影响. 研究结果表明潜堤相对于水深的高度越大,对波浪在直墙上爬高的影响越大,而在一定的相对高度条件下，潜堤与直立墙之间距对波浪爬高的影响呈现出一定的规律性.

2009, 41(5): 660-665. doi: 10.6052/0459-1879-2009-5-2008-291

PDF

结合自适应网格的基于特征线方程的分离算法

陈炎曹树良梁开洪祝宝山

提出一种新的网格自适应方法：在需要加密的网格单元中心加入新结点，并对加密后的相邻三角形网格单元进行公共边变换, 构成新的网格单元. 与传统的在网格边界中点加入新节点的自适应方法相比，新方法可以更加灵活地控制网格密度，加密后的网格继承原先的网格质量不发生畸变，并且算法编程简便，容易实现. 将自适应网格生成方法和基于特征线方程的分离算法相结合，对空腔内不可压缩黏性流动进行了计算. 在特征线方向上进行时间步离散，动量方程求解过程中采用非增量型分离算法. 计算中，把求解变量梯度值作为判定准则，在变化剧烈的区域进行网格局部加密. 计算结果表明该组合算法有很好的计算精度，并有效减少了计算时间和存储量.

2009, 41(5): 666-670. doi: 10.6052/0459-1879-2009-5-2008-446

PDF

基于西原模型的非线性黏弹塑性流变模型

李良权, 徐卫亚, 王伟

在西原模型的基础上，引入能反映加速蠕变起始点的非线性黏塑性元件和瞬时塑性变形的塑性元件，建立了一个能完整反映蠕变全过程的非线性黏弹塑性流变模型. 在其一维流变模型结构的基础上建立三维的流变模型，推导三维流变方程，接着探讨了在三轴压缩试验中确定三维蠕变参数的方法，并对参数进行了敏感性分析. 用三轴压缩试验结果验证模型的有效性，结果表明模型能很好地模拟蠕变的三阶段，并且在多级偏应力水平下可以采用同一组蠕变参数.

2009, 41(5): 671-680. doi: 10.6052/0459-1879-2009-5-2008-068

PDF

双周期平面裂纹问题的特征展开-变分方法

严鹏蒋持平

根据弹性力学的变分原理，利用双周期问题位移场的双准周期性质和应力应变场的双周期性质，构造了双周期平面问题的单胞泛函变分表达式. 然后结合针对裂纹问题的复应力函数特征展开式，发展了基于单胞模型的双周期裂纹平面问题的特征展开-变分方法. 由于该方法考虑了最一般的双周期边界条件，因而能够分析一般非对称排列的双周期裂纹问题. 通过结果的收敛性分析说明了该方法具有计算效率和精度都高的优点. 最后，对于裂纹呈平行四边形排列的情况，分析了不同的裂纹排列对应力强度因子的影响.

2009, 41(5): 681-687. doi: 10.6052/0459-1879-2009-5-2008-117

PDF

考虑运动硬化的DD3镍基合金力学行为研究

聂君锋由小川庄茁李旭东

详细介绍了镍基合金的晶体塑性本构模型，在Asaro大变形晶体塑性框架下，详细介绍了镍基合金的晶体塑性本构模型，在Asaro大变形晶体塑性框架下，引入了运动硬化规律，考虑了温度和应变率对晶体塑性变形的影响，通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导，通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应，并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性，体现了镍基单晶在循环载荷作用下的拉-压不对称性.

2009, 41(5): 688-695. doi: 10.6052/0459-1879-2009-5-2008-039

PDF

热力学相容的混凝土各向异性单边损伤模型

吴建营

微裂缝演化(损伤)引起的各向异性和单边效应，对于混凝土材料和结构的变形和内力有非常重要的影响. 在描述单边效应时，已有的各向异性损伤模型均会得出与热力学基本原理相矛盾的结论即损伤卸载时能量耗散不为零. 基于不可逆热力学和内变量理论，直接以材料柔度张量的增量作为损伤内变量，建立了各向异性单边损伤模型的一般形式，给出了热力学相容的投影算子，推导了模型的率本构关系.文中详细发展了模型的Newton-Raphson数值实现算法及其算法一致性模量，建立了合理的损伤准则和演化法则并应用于混凝土材料. 数值模拟结果初步验证了建议模型的有效性. 需要说明的是，模型完全构筑于热力学基础上，无需引入应变等效或应变能等效等经验性假定.

2009, 41(5): 696-707. doi: 10.6052/0459-1879-2009-5-2008-109

PDF

PTFE/Al含能复合材料的压缩行为研究

徐松林阳世清赵鹏铎李俊玲卢芳云

金属/氟聚合物含能复合材料是一类新型的高级含能材料. 研究了室温下Al含量和应变率对PTFE/Al含能复合材料压缩性能和反应性能的影响，所加载的应变率为6\times10^{-3}s^{-1}\sim8\times 10^{3}s^{ -1}. 材料压缩性能的应变率效应明显：与静态加载相比，动态加载下材料模量和强度明显提高，但应变降低. 材料的损伤过程主要包括塑性变形、开裂和反应3部分. 随着Al含量的增加，材料准静态和动态压缩强度均呈先升后降的趋势，在 Al含量为35\%时达到最高值102.6 MPa和154 MPa; 引发反应所需加载的应变率增加，但对应的应力值差别不明显，基本在165 MPa左右，材料引发后反应完全性降低.

2009, 41(5): 708-712. doi: 10.6052/0459-1879-2009-5-2009-253

PDF

弹塑性结构安定下限分析的无网格局部

陈莘莘刘应华岑章志

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合，建立了一种安定下限分析的新方法．为了克服移动最小二乘近似难以准确施加本质边界条件的缺点，采用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念，消除了安定下限分析中由时间参数所引起的求解困难．利用减缩基技术，将安定分析问题化为一系列未知变量较少的非线性规划子问题．在每个非线性规划子问题中，自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模拟，而这些自平衡应力场基矢量可应用弹塑性增量分析中的平衡迭代结果得到．算例结果证明了提出的分析方法的有效性．

2009, 41(5): 713-721. doi: 10.6052/0459-1879-2009-5-2008-342

PDF

多尺度材料模拟的自适应有限元方法研究

徐云陈军蔚喜军

研究了材料模拟中一类新型耦合多尺度的自适应有限元方法. 采用微观分子动力学耦合宏观有限元的桥尺度方法来模拟材料破坏的前期行为，其中宏观有限元计算推广到了一般非结构三角形网格. 材料破坏形成后，停止微观尺度的计算，它的进一步发展和演化通过一个宏观模型来描述，采用自适应有限元方法来求解这一宏观模型. 其中，后验误差估计的基础是变分多尺度理论，即自适应网格加密是基于粗尺度上残差分布和细尺度上单元Green's函数. 计算中采用了破坏准则来模拟材料的断裂. 数值实验表明了方法的有效性.

2009, 41(5): 722-729. doi: 10.6052/0459-1879-2009-5-2008-024

PDF

基于球形压入法提取材料的塑性力学参数

姜鹏张泰华杨荣梁乃刚

基于球形压入的工作方式，重点研究材料塑性力学参数的表征方法. 首先，通过对球压入可测量的综合评价，选取能量比(可释放功和压入总功之比)作为主要分析参数; 其次，利用孔洞模型、量纲分析和数值模拟等工具，定义特定的代表性应变，建立起能量比和材料塑性力学参数之间的关系式, 由此，提出一种能提取材料屈服应力和硬化指数的力学表征方法，该方法避免了对接触半径的测量，确保了方法的可操作性，并通过参数重组提高其稳定性; 最后，选用45号钢和6061铝合金进行压入试验，与拉伸试验结果比对显示: 该方法识别的塑性力学参数能满足工程测试需要.

2009, 41(5): 730-738. doi: 10.6052/0459-1879-2009-5-2008-194

PDF

动能侵彻弹体的质量侵蚀模型分析

陈小伟杨世全何丽灵

通过对高速侵彻弹体的质量侵蚀实验现象进行分析，建立了质量侵蚀的工程理论模型，可对动能侵彻弹体头部侵蚀后的头形和质量侵蚀进行预期. 讨论弹体质量侵蚀的主要影响因素，认为弹体动能(或质量和速度)以及混凝土骨料硬度对弹体的质量侵蚀有显著影响，给出了更一般化的弹体质量侵蚀与撞击函数$I$的关系. 通过图形分析明确弹体的质量侵蚀主要发生于头部，并且侵蚀后的弹头仍接近尖卵形，可以通过头形变化预期弹体的质量损失.

2009, 41(5): 739-747. doi: 10.6052/0459-1879-2009-5-2008-295

PDF

基于局部模态的约束子结构模型修正法

侯吉林欧进萍

针对局部子结构为修正对象的情况提出了约束子结构修正法，实现只利用整体结构模态中对应子结构部分的模态即可以修正子结构模型.由脉冲响应结合特征系统实现法识别出子结构的低阶模态；结合识别的模态和整体结构理论模型的高阶模态构造整体结构对应子结构位置的柔度矩阵；利用柔度矩阵的物理意义，在子结构的边界上施加数值支座，把子结构从整体结构中隔离出来成为约束子结构，同时构造出约束子结构的柔度矩阵；利用灵敏度的方法根据构造出的约束子结构柔度矩阵，优化修正约束子结构，即间接等效地修正子结构模型.通过一个平面桁架结构验证了约束子结构模型修正法的可行性与有效性，即使在5%或10%的噪声影响下，仍能得到满意的修正结果.关键词模型修正，柔度矩阵，约束子结构，灵敏度，修正单元力

2009, 41(5): 748-756. doi: 10.6052/0459-1879-2009-5-2008-319

PDF

触摸法评价物体柔软性的感知觉力学原理分析

胡吉永丁辛王如彬林炜

物体柔软性的触觉认知机理是当前感官工程研究的重点. 基于触摸法评价物体柔软性时物体与指尖软组织的力学作用过程，及指尖皮下机械刺激感受体接受机械刺激的模式，建立了与这一机械作用过程等效的生物力学模型; 基于该模型，提出并证明了表征触觉系统感应、辨别物体柔软性的能力指标------机械敏感性和感觉敏感性; 以提出的两个指标，参数化讨论了人的触觉系统在判断和辨别物体柔软性的一般力学原理. 理论研究证明，以触摸方式评价物体柔软性时，皮下触觉感受体及相应高级认知系统的感知能力取决于物体与皮肤软组织的抗压力阻，且与触摸力的动态变化有关. 也正是传递机械信号的介质------软组织的可变形性致使触觉感测系统有别于实验室压缩仪的工作模式.

2009, 41(5): 757-764. doi: 10.6052/0459-1879-2009-5-2008-301

PDF

周期激励下分段线性电路的动力学行为

张晓芳陈章耀季颖毕勤胜

基于四阶自治分段线性电路的分岔特性，探讨了两种幅值周期激励下该电路系统的复杂动力学行为. 给出了弱周期激励下系统共存的两种分岔模式及其产生的原因，讨论了不同分岔模式下动力学行为的演化过程及混沌吸引子相互作用机理. 而随着激励幅值的增大，即强激励作用下，围绕两个原自治系统平衡点的周期轨道不再分裂，从而导致共存的分岔模式消失.指出无论在弱激励还是在强激励下，由于系统的固有频率与外激励频率存在量级上的差距，其相应的各种运动模式，诸如周期运动、概周期运动甚至混沌运动均表现出明显的快慢效应，进而从分岔的角度分析了不同快慢效应的产生机制.

2009, 41(5): 765-774. doi: 10.6052/0459-1879-2009-5-2008-075

PDF

双稳系统中随机共振与相同步时间的相关性

董小娟

研究了双稳系统的随机共振(SR)与输入输出之间的相同步相关性. 首先得到系统的输出信噪比(SNR)与输入输出之间的平均相同步时间的表达式，然后讨论了随机共振与输入输出之间的平均相同步时间之间的关系. 结果表明：(1)系统出现了随机共振现象，且平均相同步时间对噪声是敏感的；由于加性与乘性噪声的相互影响，随加性噪声强度的增加，输出信噪比及平均相同步时间曲线上首先出现抑制现象，然后出现峰值，并且, 减小乘性与加性噪声强度比率，可提高输出信噪比和增长平均相同步时间. (2)系统的随机共振与平均相同步时间达到最大值不同步出现,但平均相同步时间对输出信噪比是敏感的. 该结论为信号传输中利用随机共振原则改变系统工作环境提供了依据.

2009, 41(5): 775-782. doi: 10.6052/0459-1879-2009-5-2008-549

PDF

基于文献计量的《力学学报》载文被引分析解读

王惠翔陶彩军

利用中文科技期刊引文数据库中的数据，对中文核心期刊《力学学报》载文在1990-2008年6月的19年间，被国内中文期刊引用情况(包括被引论文的频次分布、年代分布、引用期刊分布和核心著者分布等)的调查研究与分析解读，揭示《力学学报》在中国学术界的整体水平及主要影响.

2009, 41(5): 783-788. doi: 10.6052/0459-1879-2009-5-2008-563

PDF

点涡动力系统的Lyapunov指数

朱浩, 程克明

平面上理想流体的三点涡系统是可积的Hamilton系统, 但其运动仍然相当复杂, 这给研究被动微粒在三点涡系统中运动带来了很大的困难. 着眼于点涡系统的被动微粒对初始小扰动的稳定性, 通过Oseledec定理定义被动微粒的Lyapunov指数, 给出了点涡系统中被动微粒稳定性的定量刻画. 同时, 由Hamilton系统的保体积性质得到的关于Lyapunov指数的简洁表达式, 避免了计算的繁琐. 利用这个定义, 点涡系的瞬时流场可以被划分成若干区域, 被动微粒的混沌运动只能在近涡的特定区域出现.

2009, 41(5): 789-793. doi: 10.6052/0459-1879-2009-5-2008-469

PDF

内皮抑素抑制肿瘤血管生成的数值模拟

赵改平陈二云沈力行郦鸣阳吴洁许世雄丁皓

数值模拟抗血管生成药物内皮抑素对肿瘤血管生成的抑制效应. 建立内皮抑素作用下肿瘤内外血管生成的二维、三维离散数学模型，模型中考虑内皮抑素的抑制作用、内皮细胞自身的增殖、促血管生成因子TAF和Fibronectin对内皮细胞产生的趋化性和趋触性以及内皮细胞自身扩散引起的随机性运动，数值模拟肿瘤内外微血管网的生成过程. 模拟结果表明，抗血管生成药物内皮抑素对肿瘤内外血管生成的速度、成熟度以及血管分支数量均有明显的抑制作用，从而有效地抑制肿瘤新生血管的形成. 该模型能够较好地模拟内皮抑素对肿瘤血管内皮细胞迁移与增殖的抑制效应，为临床抗血管生成治疗肿瘤提供有益的信息.

2009, 41(5): 794-800. doi: 10.6052/0459-1879-2009-5-2008-122

PDF

渗透各向异性可压缩地基固结的平面应变分析

艾智勇吴超

对多层地基的平面应变固结问题进行了研究，并同时考虑了土体的渗透各向异性和孔隙流体的可压缩性. 从平面应变Biot固结的控制方程出发，对时间t, 坐标z和x进行Laplace和Fourier变换，建立了地基表面(z=0)和任意深度z处的基本量在Laplace-Fourier变换域内的传递矩阵关系. 利用传递矩阵法，结合土层连续条件和边界条件，并应用Laplace-Fourier逆变换技术，推导出渗透各向异性可压缩多层地基平面应变固结的理论解. 基于该解，编制了计算程序，并进行了数值计算. 讨论了土体的渗透各向异性、孔隙流体的可压缩性以及地基的分层特性对地基固结的影响，分析结果表明：土体的渗透各向异性、孔隙流体的可压缩性，以及地基的分层特性对地基的固结行为有着重要的影响.

2009, 41(5): 801-807. doi: 10.6052/0459-1879-2009-5-2008-477

PDF

一类指数矩阵函数及其应用

富明慧林敬华

研究了一阶常微分方程组特解的精细积分方法. 针对非齐次项为多项式、指数函数以及二者的乘积的情况，在Duhamel积分形式特解的基础上，引入了一类指数矩阵函数. 通过该类函数的线性组合即可表达出非齐次方程的特解. 建立了该类指数矩阵函数的一种高效递推算法，并在此基础上实现了特解的精细积分. 由于特解的积分过程能充分利用通解精细积分过程的中间量，因此两个精细积分过程能有机地结合起来，形成了一种高效、统一的广义精细积分法. 对上述递推算法做了进一步优化，并给出了通用的计算公式.算例结果证明了该方法的有效性.

2009, 41(5): 808-814. doi: 10.6052/0459-1879-2009-5-2008-052

PDF

带有环形阀片的阻尼阀可靠性研究

陈轶杰杨占华雷强顺王亚军宁丹

针对跑车试验时油气弹簧中节流阀片受冲击失效的现象，提出了研究阻尼阀水击压强的必要性，建立了带有环形节流阀片的阻尼阀物理模型. 运用水击理论及解偏微分方程的特征线法，创造性地对给定阻尼阀结构进行了水击建模计算. 编程分析节流阀片失效的原因，系统研究了常通孔结构参数变化对阻尼阀所受瞬态冲击及其阻力值的影响规律，提出需从屈服和疲劳强度两个方面入手对阀门可靠性进行校核，并给出计算实例，完善了现有阀系设计的算法. 所得结论可以作为可靠性研究的参考，同时避免了阀片断裂的再次发生.

2009, 41(5): 815-820. doi: 10.6052/0459-1879-2009-5-2008-546

PDF

留言板

2009年第41卷第5期

通知公告MORE+

精品文章MORE+

下载中心

作者中心

留言板

2009年 第41卷 第5期

通知公告MORE+

精品文章MORE+

下载中心

作者中心

2009年第41卷第5期