喷嘴挡板伺服阀温漂电流波动研究
李磊1, 李长春1, 延皓1, 秦嘉言1, 王鹏2
1.北京交通大学 机械与电子控制工程学院, 北京 100044
2.北京精密机电控制设备研究所, 北京 100076

第一作者:李磊(1988—),男,湖北仙桃人,博士生.研究方向为流体传动与控制.email:12116315@bjtu.edu.cn.

摘要

为了找出伺服阀的温漂电流波动的原因,考虑到温度对主阀芯的影响,从仿真中阐述了阀芯楔形膨胀的变化规律,将阀芯轴向局部温度分布简化为一维稳态导热问题,推导了阀芯与阀套间隙的楔形变化规律表达式.在此基础上,利用滑阀间隙的流动规律,推导出了滑阀楔形变形下的库伦摩擦力表达式,得出了在库伦摩擦力下伺服阀的电流变化量.最后,通过理论分析和试验进行了对比,证明了伺服阀温漂电流的波动与温升引起的伺服阀内部非线性摩擦有关.

关键词: 伺服阀; 非线性摩擦; 间隙流动; 温漂
中图分类号:TH137 文献标志码:A 文章编号:1673-0291(2017)03-0096-07
Research on temperature drift current fluctuation of nozzle flapper servo valve
LI Lei1, LI Changchun1, YAN Hao1, QIN Jiayan1, WANG Peng2
1.School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044,China
2. Beijing research institute of precision electromechanical control equipment, Beijing 100076,China
Abstract

In order to find out why temperature drift can cause current fluctuation in servo valve, the changing rules of how temperature produces a wedge effect on the spool is elaborated by simulation,in which the effect of thermal expansion on main spool is taken into consideration. The local temperature distribution in spool axial dimension is simplified to a one-dimensional steady-state heat conduction. Thus an expression which show the changing rules of the wedge-shaped clearance between the spool and spool sleeve is derived. Accordingly, a coulomb friction expression under slide valve wedge deformation is deduced based on clearance flow rules of slide valve, by which the current variable quantity of servo valve under coulomb friction is obtained. Finally, by comparing theoretical analysis and tests, it is proved that the current fluctuation in servo valve caused by temperature drift is relevant to the non-linear friction caused by temperature rise inside servo valve.

Keyword: servo valve; nonlinear friction; clearance flow; temperature drift

火箭发射过程中要承受极限高温的恶劣环境, 由于高温燃气源的影响, 随着设备的运行, 伺服机构中液压油液温度由常温上升至150 ℃.作为伺服机构中重要的控制元件, 大部分伺服阀都在其零位附近工作[1], 且它的工作零位会随着油温变化发生漂移, 特别是在大批次测试中, 个别伺服阀在局部温度范围内会出温漂电流幅值急剧波动的现象, 在工况下, 如图1(a)所示, 回油温度能升高至90 ℃.当温漂电流波动较大时, 如图1(b)所示, 行业内称此现象为伺服阀“ 温漂突跳” .这种现象对液压伺服系统的可靠性、控制精度和稳定性会有很大影响[2].故找出伺服阀的温漂电流随温度波动的原因和影响因素具有重要的理论和工程意义.

图1 实际测试中的温漂现象Fig.1 Temperature drift phenomenon in test

伺服阀的制造涉及到军工方面, 具有保密性质, 因而伺服阀温漂电流急剧波动的原因在国外文献中较难查阅到相关资料.在国内, 温漂电流阈值往往是使用经验值, 甚至忽略不计.关于温度对伺服阀影响文献中, 肖其新[3, 4]指出当阀芯槽和控制棱边配合不准确时, 阀芯和阀套的热膨胀会引起零漂.李成功等[5, 6]利用能量守恒定律建立了航空柱塞泵的传热模型.文献[7, 8, 9]采用了数值有限元仿真的方法进行了滑阀和喷嘴挡板阀的热流固耦合面仿真.吕玥婷等[10, 11, 12, 13]利用有限元分析阐述了热对滑阀卡滞的影响.Therman E.Green[14]采用了红外热成像仪进行液压系统的故障诊断.徐民[15]提出了伺服阀温度零漂特性的测试方法.从国内外文献可以看出, 关于温度对伺服阀影响研究, 大部分学者也只是集中在系统层面的热仿真.

为了测试伺服阀温度漂移, 某航天院某所采用的是位置闭环测试法, 伺服阀温度筛选原理图如图2所示, 对应的伺服阀温度筛选实物图如图3所示, 在伺服机构的参与下, 利用伺服机构自带的多余度位置传感器, 伺服作动器做闭环位置运动.整个测试系统在伺服阀的过渡块上安装有测量回油温度的温度传感器, 并且在伺服放大器上加装一个100 Ω 的电阻, 用来进行伺服阀电流取样[16].通过伺服阀放大器给伺服阀线圈输入电流, 实现计算机对伺服机构的正弦位置控制.通过调整给定位置幅值实现对线圈电流幅值的调整, 使峰值电流达到额定电流的5%[17].

从上述文献可以看出, 针对伺服阀温漂中的这种特殊电流波动现象总结和分析还比较欠缺.本文作者根据以上测试原理和计算方法, 对某航天院某所的多种不同型号的伺服阀进行数百次抽样测试分析, 整个伺服机构中油温变化缓慢(< 0.1 ℃/s), 伺服阀温漂特性可以视为稳态问题.通过在此工况下考察不同伺服阀在伺服机构上温漂变化规律, 最后得出伺服阀温漂电流的波动与温升之间的关系.

图2 伺服阀温漂测试原理图Fig.2 Servo valve temperature drift test schematic

图3 被测伺服机构Fig.3 Mechanism of tested servo

1 阀芯和阀套的间隙热变形

滑阀在微小的开口情况下, 油液流速高, 粘性加热显著, 易造成滑阀局部高温的情况, 可能会引起阀芯和阀套的轴向局部变形.

文献[10]对滑阀的表面温度场进行了仿真, 仿真结果表明在滑阀的中间段及槽口附近是明显的高温区, 阀芯两端是低温区, 在阀口的节流处产生大量的热, 这一部分热量随着油液从出口流出, 一部分通过槽口表面传递给阀芯.高温下滑阀温度分布的仿真结果如图4所示, 仿真结果表明, 由于阀口之后与阀芯相接触的壁面附近流体速度关于阀芯轴线不对称, 速度较高一侧的内摩擦力较大, 产生的热量也较多, 因而温升较高, 导致阀芯的温度分布不对称.

图4 滑阀温度场分布示意图Fig.4 Spool temperature field distribution diagram

在温筛试验中, 滑阀阀芯在滑套中做开度极小且缓慢的正弦运动, 随着开度的变化, 在阀芯的槽口附近产生较大部分热量, 在槽口附近出口之间, 温度会有明显的温度梯度.

阀芯与阀套之间的配合为间隙配合(0.002~0.004 mm), 而阀体与阀套之间的配合大多数为过盈配合(一般为0.002 mm), 如图5所示, 图中 rzrarb分别为阀芯的半径、阀套内径和阀体内径, 阀芯和阀套可取同一材质Cr12MoV.如果阀芯和阀套都为自由膨胀, 则温度对缝隙变化的影响不大.但是, 考虑到阀体和阀套之间的配合为过盈配合, 限制了阀套外径方向的膨胀自由度, 则会影响阀芯和阀套的配合间隙, 甚至造成卡死的现象.

图5 阀芯和阀套的径向尺寸Fig.5 Radial dimension between spool and spool sleeve

在实际过程中, 阀芯的槽口温度要高于阀套的温度, 所以相比较可以忽略阀套的热膨胀效应, 只考虑阀芯的热膨胀.在自由膨胀的情况下, 阀芯的外表面的径向变形为

Δuzo=2αrz0rzTz(t)rdr(1)

式中: α为阀芯的材料线膨胀系数; r为阀芯方向的积分半径; Tz(t)为阀芯沿着半径方向的温度分布函数.

由图4仿真可知, 阀芯在径向的温度分布并非为恒定值, 轴的温升函数为指数分布[18], 且当一段时间后, 温度分布趋向稳定, 但是由于阀芯本身径向温差不大, 阀芯外壁较薄且导热性好, 可以简化认为阀芯在径向的温度分布为恒定值 ΔT.

则阀芯和阀套的配合间隙的变化量为

δt=2αrzΔT(2)

则阀芯和阀套的实际工作间隙为

h0=δ0-δt(3)

式中: δ0为温升前阀芯与阀套间的初始配合间隙.

从式(2)和式(3)可以看出, 温度对阀芯与阀套的间隙影响很大, 故在轴向的温度分布也将会影响轴向的阀芯与阀套配合间隙 h(x), 其可以表示为

h(x)=δ0-2αrzΔT(x)(4)

式中: ΔT(x)为阀芯x方向的温差分布.

实际上 ΔT(x)的温度函数分布模型可以简化为一维圆柱的导热问题[19], ΔT(x)可以表示为

ΔT(x)=T2-T1Lx+T1(5)

式中: T1为油入口阀芯的温度; T2为油出口阀芯的温度, T1> T2; L为阀芯突肩的宽度.

将式(5)带入式(4), 得出阀芯和阀套的配合间隙

h(x)=2αrzT1-T2Lx+(δ0-2αrzT1)(6)

从式(6)可以得出边界条件 h1h2

h1=hx=0=δ0-2αrzT1h2=hx=L=δ0-2αrzT2(7)

将式(7)带入式(6), 式(6)可以简化为

h(x)=h1+(x/L)(h2-h1)=h1+ix(8)

式中: i=(h2-h1)/L, 表示为阀芯与阀套楔形间隙 x方向的变化梯度.

式(8)即为阀芯与阀套的间隙线性化的函数表达式, 由计算结果可知, h(x)在阀芯的槽口附近因温度上升过高, 会出现摩擦力过大甚至卡滞现象( h(x)< 0时).这种楔形变化规律与文献[10]仿真中滑阀的受热变形趋势相吻合, 图6为其仿真中放大了2 000倍的阀芯受热变形情况.

图6 阀芯变形放大图Fig.6 Valve spool deformation enlarged drawing

2 滑阀的间隙流动

由上述可知, 阀芯轴向温度分布不同, 导致了阀芯与阀套之间缝隙宽度成近似楔形变化, 图7为阀芯与阀套配合示意图, 因制造安装工艺存在缺陷而导致了阀芯与阀套间存在配合偏心问题, 令其偏心量为 e.假设油液由左边 P1高压区流入右边 P2低压区( P1> P2), c1c2为阀芯和阀套的初始间隙, 由于阀芯的楔形变化, 有 c1> c2.下面讨论阀芯的油液流动规律和阀芯上 y方向的力 Fy情况.

图7 滑阀的配合示意图Fig.7 Spool valve assembly schematic

当楔形的角度很小时, 可以忽略流动 y方向的速度分量, 阀芯与阀套之间的变化 h(x)极缓慢, 所以可以将滑阀之间的层流近似于Couette流, 如图8所示, Couette流描述的是两块平板间油液的相对流动, 一块平板静止不动, 另一块是以恒定的速度 U相对运动, 其中 h为两块平板之间的距离, 油液压力 P1> P2.

图8 Couette流模型Fig.8 Couette flow model

将纳维-斯托克斯方程[20]应用到Couette流模型中, 可得到表达式如下

dpdx=μd2udy2(9)

式中: μ为油液的粘度; u为油液x方向的速度.由图8可知, 边界条件 y=0, u=0y=h, u=U可以得出两平板之间的速度流动规律

u=Uyh-h22μdpdxyh(1-yh)(10)

设两板间流量为 q, 由于 q沿着流道不变, 则

q=0hudy(11)

将式(10)代入式(11)中得

dpdx=12μqh3(Uh2q-1)(12)

式(12)中 h可以用x方向的一次函数表示, 将式(8)中 h(x)代替式(12)中 h, 由于 x=0时, p=p1, 解微分方程(12)得

p=p1+6μqi(h1+ix)2-6μqih1(13)

此外, 由于 x=L时, h=h2p=p2, 代入式(13)中, 可得

p=p1-h22(p1-p2)h22-h12(1-h12(h1+ix)2)(14)

式(14)为 x方向上流体压强的分布表达式, 将式(14)沿着阀芯突肩宽度 L方向进行积分, 可以求出单位宽度上的阀芯受力 f, 即

f=0Lpdx=[p1-h2h1+h2(p1-p2)]L(15)

滑阀侧视图如图9所示, 图中阀芯与阀套沿圆周的缝隙为

h1=c1-ecosβh2=c2-ecosβ(16)

式中: β为阀套的圆周角.

图9 滑阀侧视图Fig.9 Spool side view

因此作用在阀芯上y方向的力Fy

Fy=02πfRcosβdβ(17)

将式(16)代入式(15)中, 将式(15)带入式(17)中, 可以得出最终阀芯上 y方向的力

Fy=π(c2-c1)LR(p1-p2)2e11-[2e/(c1+c2)2]-1(18)

式(18)描述的阀芯液压力使阀芯向更加狭窄的方向作用, 阀芯在 y方向力的作用下, 外圆周与阀套内壁接触, 最后有 c1=e, 为了计算方便, 令 λ=c2/c1, 则式(18)可以简化为

Fy=πLR(p1-p2)2(1+λ)λ-1λ+3-λ+1(19)

由于阀芯和阀套之间是缓慢的正弦运动, 可以忽略粘滞摩擦力, 主要考虑轴向的库伦摩擦力, 即

fc=Fyμ(20)

3 摩擦力导致的电流变化分析

阀芯所受液压力使阀芯向更加狭窄的方向作用, 而且只有当 P1> P2时, 才会产生多余的库伦摩擦力 fc.如图10所示, 当主阀芯在零位时, 在库仑力的摩擦作用下, 为了使主阀芯维持受力平衡, 必须给伺服阀作用纠偏电流 ΔI, 使力矩马达发生偏转角度 θ, 图中s为挡板距离弹簧管旋转中心的距离, b为反馈杆末端距离挡板的距离.

图10 力矩马达衔铁示意图Fig.10 Torque motor armature schematic diagram

忽略回油压力, 因主阀芯缓慢运动且在零位时, 主阀芯受到阀芯端面压差产生的驱动力、摩擦力和反馈杆的作用力.阀芯的力平衡方程为

plvAv=fc+Kfθ(s+b)(21)

式中: plv为主阀芯端面压差; Av为主阀芯端面截面积; Kf为反馈杆的刚度; θ为力矩马达偏转角度.

前置级喷嘴挡板阀的零位线性化后的压力特性方程为

plv=psxf0xf(22)

式中: ps为供油压力; xf为挡板偏离零位的位移; xf0为挡板与喷嘴的零位间隙.

力矩马达的输出力矩 Td

Td=KtΔI+Kmθ(23)

式中: Kt为力矩马达中位电磁力矩系数; Km为力矩马达电磁弹簧力矩.

主阀芯偏移通过反馈杆作用在旋转中心的力矩为

Tf=Kfθ(s+b)2(24)

当力矩马达平衡时, 如果考虑作用在挡板的液压力, 则有如下平衡方程

Td=Tf+Kaθ+psANs(25)

式中: Ka为弹簧管的刚度; AN为喷嘴的横截面积.

综合式(21)~式(25), 可得伺服阀的电流变化量

ΔI=[kf(s+b)2+ka-km]×xf0kt[psANs-kfxf0(s+b)]fckt(26)

由式(26)可知, 伺服阀的零偏电流值与温度导致的库伦摩擦力成正比, 随着阀芯的运动, 主阀芯的温度会有短时间的分布不均, 但是随着试验的进行, 主阀芯温度最终分布均匀, 此时, 由主阀芯的温度导致的电流变化将消失.

4 试验结果分析

阀芯两端温度差并不容易测量, 只能通过fluent仿真得到模拟值, 根据文献[10]中的数据, 取阀芯两端温差为20 ℃, 配合间隙 δ0取0.007 mm, 伺服阀的阀芯 rz取7.5 mm, α为11.4× 10-6 ℃, T2=80 ℃, T1=50 ℃, 代入式(19)推导出作用在阀芯 y方向的液压力能达到16.9 kg, 最后再代入式(26), 得出摩擦力导致电流的幅值变化结果为0.2 mA, 图11为伺服阀模拟滑阀摩擦的仿真电流曲线.在仿真中第100 s处添加死区环节, 模拟由阀芯和阀套之间增加的库伦摩擦力导致的死区, 在110 s处撤销死区环节.从仿真结果中可以看到仿真曲线的电流的在第100 s处和第110 s处的波动情况.仿真表明, 添加死区环节来模拟库伦摩擦力导致的死区, 会导致伺服阀零漂电流的波动, 这与之前分析的理论相符合.

图11 温漂波动的仿真曲线图Fig.11 Servo valve temperature drift fluctuation simulation

对某型号的力反馈两级伺服阀进行温漂特性参数测试试验, 在规定试验的时间内, 伺服阀工作的油温从31.85 ℃上升到99.78 ℃, 实际的试验结果如图12所示.在825~900 s处, 伺服阀线圈零偏电流发生较为剧烈的波动, 零偏跳动幅值为0.2~0.3 mA, 变化趋势和大小与实际能够吻合. 这是由实际中滑阀的温度分布不均导致的, 所以引起的库伦摩擦力在温度影响下也不是一成不变的, 呈现出如图12所示零漂波动情况.

图12 温漂波动的实验曲线图Fig.12 Results of temperature drift fluctuation

5 结论

1)试验表明, 在线圈温漂电流幅值出现局部急剧波动现象的时候, 温漂的波动整体趋势呈现出单调变化, 在波动处正弦电流幅值变化较大.

2)伺服阀内部由于局部温升会形成楔形结构, 增大阀芯与阀套之间的库伦摩擦力, 产生非线性死区, 最后通过理论推导和试验对比, 证明了伺服阀温漂电流的波动与温升引起的伺服阀内部非线性摩擦有关.为了减小温漂的波动, 建议在生产和使用过程中, 控制油液清洁度的同时, 可以在阀芯上开平衡槽, 这样可以减小阀芯上的压强分布的不均匀性, 也可以减小楔形的变形量, 从而减小阀芯与阀套间的库伦摩擦力.

此外, 喷嘴挡板作为伺服阀前置级, 具有更高的压力灵敏度, 温度的升高对其元件的尺寸和流体的流动特性影响将更加明显, 这也是本文没有考虑到的问题, 今后将作进一步讨论.

The authors have declared that no competing interests exist.

参考文献
[1] 孙萌, 李长春, 延皓, . 阀控非对称缸位置系统的非线性建模[J]. 北京交通大学学报, 2012, 36(4): 164-168.
SUN Meng, LI Changchun, YAN Hao, et al. Nonlinear modeling for valve controlled asymmetric cylinder position system[J]. Journal of Beijing Jiaotong University, 2012, 36(4): 164-168. (in Chinese) [本文引用:1]
[2] 肖其新. 高低温对电液伺服阀性能影响[J]. 机械工程师, 2008(10): 21-23.
XIAO Qixin. Analysis of the high and low temperature impact to the electro-hydraulic servo valve[J]. Mechanical Engineers, 2008(10): 21-23. (in Chinese) [本文引用:1]
[3] 肖其新. 温度对电液伺服阀特性影响[D]. 上海: 同济大学, 2009.
XIAO Qixin. Influence of temperature on the characteristics of electro-hydraulic servo valve[D]. Shanghai: Tongji University, 2009. (in Chinese) [本文引用:1]
[4] 訚耀保, 肖其新, 闫世敏. 温度对电液伺服阀的影响分析[J]. 流体传动与控制, 2008(6): 23-26.
YIN Yaobao, XIAO Qixin, YAN Shimin. Influence of temperature on hydraulic servovalve[J]. Fluid Power Transmission & Control, 2008(6): 23-26. (in Chinese) [本文引用:1]
[5] LI Chenggong, JIAO Zongxia. Thermal-hydraulic modeling and simulation of piston pump[J]. Chinese Journal of Aeronautics, 2006, 19(4): 354-358. [本文引用:1]
[6] LI Chenggong, JIAO Zongxia. Calculation method for thermal-hydraulic system simulation[J]. Journal of Heat Transfer, 2008, 130(8): 318-323. [本文引用:1]
[7] ZHAO Jianhua, ZHOU Songlin, LU Xianghui, et al. Numerical simulation and experimental study of heat-fluid-solid coupling of double flapper-nozzle servo valve[J]. Chinese Journal of Mechanical Engineering, 2015, 28(5): 1030-1038. [本文引用:1]
[8] 周松林. 电液伺服阀喷嘴挡板组件流固耦合分析及主阀温度场分析[D]. 秦皇岛: 燕山大学, 2014.
ZHOU Songlin. Fluid-structure interaction analysis of nozzle flapper assembly and temperature field analysis of the main valve of Electro-hydraulic servo valve[D]. Qinhuangdao: Yanshan University, 2014. (in Chinese) [本文引用:1]
[9] 周大海, 刘桓龙, 秦剑, . 基于流固热耦合的滑阀温度特性研究[J]. 机床与液压, 2014 (21): 169-171.
ZHOU Dahai, LIU Huanlong, QIN Jian, et al. Temperature characteristics study of spool valve based on thermal solid-fluid coupling[J]. Machine Tool & Hydraulics, 2014(21): 169-171. (in Chinese) [本文引用:1]
[10] 吕玥婷. 液压滑阀液固热多物理场耦合分析研究[D]. 太原: 太原理工大学, 2014.
LYU Yueting. Multi-discipline coupled liquid-solid thermal analysis of hydraulic spool valve[D]. Taiyuan: Taiyuan University of Technology, 2014. (in Chinese) [本文引用:4]
[11] 刘玉珑. φ8通径伺服阀滑阀流场的CFD仿真与分析[D]. 秦皇岛: 燕山大学, 2012.
LIU Yulong. CFD simulation and analysis of flow fluid of spool valve of φ8 diameter electro-hydraulic servo valve[D]. Qinhuangdao: Yanshan University, 2014. (in Chinese) [本文引用:1]
[12] 罗绍亮. 液压滑阀卡紧现象的理论分析和解决方案[J]. 机电工程技术, 2008, 37(6): 93-95.
LUO Shaoliang. Theoretical analysis and solutions for hydraulic slide valve jamming[J]. Electrical Engine-ering Technology, 2008, 37(6): 93-95. (in Chinese) [本文引用:1]
[13] 晏静江, 柯坚, 刘桓龙, . 液压滑阀阀芯温度场的流-固-热耦合研究[J]. 中国机械工程, 2014, 25(6): 757-760.
YAN Jingjiang, KE Jian, LIU huanlong, et al. Research on fluid-solid-thermal coupling in temperature filed of hydraulic slide valves[J]. Chinese Journal of Mechanical Engineering, 2014, 25(6): 757-760. (in Chinese) [本文引用:1]
[14] GREEN T E, SNELL J R. Thermographic inspection of hydraulic systems[C]// Aerospace/defense Sensing & Controls, 1996: 25-30. [本文引用:1]
[15] 徐民. 基于C8051F020伺服温漂测控系统[J]. 单片机与嵌入式系统应用, 2009(11): 51-53.
XU Min. Temperature null shift measurement and control system for servo valve based on C8051F020[J]. Microcontroller & Embedded Systems, 2009(11): 51-53. (in Chinese) [本文引用:1]
[16] 母东杰. 双喷嘴挡板伺服阀流固耦合特性分析及振动抑制[D]. 北京: 北京交通大学, 2015.
MU Dongjie. Study on coupled vibration characteristics of electro-hydraulic servo valve and vibration control[D]. Beijing : Beijing Jiaotong University, 2015. (in Chinese) [本文引用:1]
[17] 李文祥. 伺服阀温漂特性与测试系统研究[D]. 北京: 北京交通大学, 2015.
LI Wenxiang. Study of temperature drift characteristic of servo valve and its testing system[D]. Beijing: Beijing Jiaotong University, 2015. (in Chinese) [本文引用:1]
[18] 黄其圣, 张勇, 胡鹏浩, . 温度变化对机械零件配合精度的影响[J]. 机械科学与技术, 2001, 18(3): 45-47.
HUANG Qisheng, ZHANG Yong, HU Penghao, et al. The influence of temperature variation on the fitting accuracy of machine elements[J]. Journal of Machine Design, 2001, 18(3): 45-47. (in Chinese) [本文引用:1]
[19] 曹克强, 李永林. 现代飞机液压系统热特性建模仿真与热设计[M]. 北京: 国防工业出版社, 2013.
CAO Keqiang, LI Yonglin. Modern aircraft hydraulic system simulation and modeling of the characteristics and thermal design[M]. Beijing: National Defend Industry Press, 2013. (in Chinese) [本文引用:1]
[20] 竹中利夫, 浦田映三. 液压流体力学[M]. 北京: 科学出版社, 1980.
TAKENAKA Toshio, URATA Eizo. Hydraulic fluid dynamics[M]. Beijing: Science Press, 1980. (in Chinese) [本文引用:1]