一種面向低頻隔振的兩自由度球面并聯(lián)機構(gòu)

張曉偉，段學(xué)超，冷國俊，?！『?/p>

(1.西安電子科技大學(xué)電子裝備結(jié)構(gòu)設(shè)計教育部重點實驗室，陜西西安　710071；2.中國電子科技集團公司第二十九研究所，四川成都　610036)

?

一種面向低頻隔振的兩自由度球面并聯(lián)機構(gòu)

張曉偉1，段學(xué)超1，冷國俊2，保宏1

(1.西安電子科技大學(xué)電子裝備結(jié)構(gòu)設(shè)計教育部重點實驗室，陜西西安710071；2.中國電子科技集團公司第二十九研究所，四川成都610036)

摘要：針對兩轉(zhuǎn)動自由度低頻隔振的需求，基于球面并聯(lián)機構(gòu)理論提出了一種兩轉(zhuǎn)動自由度球面并聯(lián)隔振裝置。從機構(gòu)學(xué)的角度對其進行了工作原理分析和自由度計算，并完成了隔振平臺的虛擬樣機設(shè)計。采用連桿參數(shù)的D-H表示法，以連桿圓弧張角為約束條件建立了機構(gòu)的約束方程，實現(xiàn)了其運動學(xué)逆解，并基于數(shù)值方法給出了運動學(xué)正解的求解過程及正逆解算例。分析和仿真結(jié)果表明：所設(shè)計的兩自由度球面并聯(lián)機構(gòu)可操控性強、運動學(xué)模型合理有效，能夠滿足隔振平臺所需的低頻隔振性能需求。

關(guān)鍵詞：機械原理與機構(gòu)學(xué)；隔振平臺；球面并聯(lián)機構(gòu)；雅可比矩陣；運動學(xué)

張曉偉，段學(xué)超，冷國俊，等.一種面向低頻隔振的兩自由度球面并聯(lián)機構(gòu)[J].河北科技大學(xué)學(xué)報,2016,37(1):13-19.

ZHANG Xiaowei, DUAN Xuechao, LENG Guojun, et al.A two-degree-of-freedom spherical parallel mechanism with low-frequency vibration isolation[J].Journal of Hebei University of Science and Technology,2016,37(1):13-19.

主動隔振技術(shù)作為振動控制領(lǐng)域的一個重要分支，具有自適應(yīng)性好、可對低頻振動進行隔離等優(yōu)點，因而成為隔振技術(shù)的研究熱點[1-2]。從隔振原理上講，主動隔振是在被控系統(tǒng)中引入次級振源，并通過一定的控制方法調(diào)節(jié)次級振源的輸出，使其產(chǎn)生的振動與主振源(干擾)的振動相抵消,從而達到隔振的目的[1]。從系統(tǒng)的功能看，隔振穩(wěn)定平臺的機械機構(gòu)是主動隔振系統(tǒng)的重要組成部分，合理的機械機構(gòu)是穩(wěn)定平臺實現(xiàn)其功能、提高穩(wěn)定精度的前提。穩(wěn)定平臺最初是串聯(lián)型機構(gòu)較多。1965年STEWART提出著名的六自由度并聯(lián)式平臺[3-4]，從此開始了基于Stewart并聯(lián)機構(gòu)的研究，目前平臺可以實現(xiàn)6個自由度的穩(wěn)定。而對于常見的少自由度的穩(wěn)定平臺，按照機構(gòu)形式可以分為2類：一類是框架機構(gòu)[5-6]，另一類為多點支撐式結(jié)構(gòu)[7]。相對于這2種結(jié)構(gòu)形式，并聯(lián)機構(gòu)具有承載能力高、剛度大、結(jié)構(gòu)簡單穩(wěn)定、精度高、易實現(xiàn)高速運動、求解運動學(xué)逆解容易等優(yōu)點[8]，從而使平臺具備更好的穩(wěn)定性。

從機構(gòu)學(xué)的角度看，機器人可分為球面機器人、平面機器人和空間機器人[9]，其中球面機器人是一特殊類型的機構(gòu)，它是聯(lián)系平面機構(gòu)與空間機構(gòu)的橋梁。當球面機構(gòu)運動時，其構(gòu)件上所有點在一個球面上運動。除具有并聯(lián)機器人的一般特性外，球面機構(gòu)還具有結(jié)構(gòu)簡單、工作空間大、不易發(fā)生干涉、運動學(xué)計算簡單、控制容易等優(yōu)點[10]。隨著機器人技術(shù)的發(fā)展，各種形式的球面并聯(lián)機構(gòu)也引起了人們的研究興趣，并在工程中得到了成功應(yīng)用[11-13]。本文提出了一種應(yīng)用于移動載體上的兩自由度主動隔振平臺機構(gòu)模型，通過隔振可減小載體由于位姿變換、干擾振動以及發(fā)動機振動等對所隔振設(shè)備的影響，使其具有一個高品質(zhì)的工作環(huán)境。

圖1　兩轉(zhuǎn)動自由度球面并聯(lián)機構(gòu)的原理圖Fig.1　Schematic of the two-degree-of-freedom spherical parallel mechanism

圖2　兩轉(zhuǎn)動自由度球面并聯(lián)機構(gòu)的虛擬樣機Fig.2　Virtual prototype of the two-degree-of-freedom spherical parallel mechanism

1機構(gòu)組成與運動原理

1.1　機構(gòu)簡介

所設(shè)計的兩轉(zhuǎn)動自由度球面并聯(lián)機構(gòu)的原理圖見圖1，機構(gòu)分為上平臺①和基座②，上平臺作為設(shè)備的承載體具有兩轉(zhuǎn)動自由度。上平臺與基座間由2個支鏈和1個支柱④連接，支鏈所有運動副都為轉(zhuǎn)動副。支柱下端與基座固連，上端通過虎克鉸③與上平臺相連。機構(gòu)符合球面機構(gòu)的一般特性[14-15]，在運動過程中支鏈連桿上的所有點均在以某一定點為球心的球面上運動，并且所有轉(zhuǎn)動副及電機軸軸線皆交匯于該定點，該定點叫做機構(gòu)的運動中心。

根據(jù)所設(shè)計機構(gòu)的特性，選擇虎克鉸中心旋轉(zhuǎn)體的型心作為機構(gòu)的轉(zhuǎn)動中心。根據(jù)機構(gòu)原理，結(jié)合機械設(shè)計相關(guān)知識，該機構(gòu)的實際設(shè)計虛擬樣機如圖2所示。

1.2　自由度計算

由機構(gòu)學(xué)知識可知，修正后的Kutzbach-Grubler自由度計算公式為[16]

(1)

式中：M為機構(gòu)的自由度數(shù)目；n為機構(gòu)中構(gòu)件的數(shù)目；g為機構(gòu)中運動副的數(shù)目；fi為第i個運動副的自由度數(shù)；d為機構(gòu)的階數(shù)；v為機構(gòu)中冗余約束的數(shù)目；ζ為機構(gòu)中局部自由度的數(shù)目。

因為平面機構(gòu)或球面機構(gòu)的公共約束為3[17]，則該機構(gòu)作為球面機構(gòu)具有3個公共約束[6]，即d=3，其他參數(shù)n=6，g=7。由于與支鏈相關(guān)的運動副都為單自由度的轉(zhuǎn)動副，所以對于這些運動副都有fi=1(i=1，2，…，6)。由于支柱與上平臺是由虎克鉸連接，故有2個自由度，所以有f7=2。從機構(gòu)的機構(gòu)形式可以看出，機構(gòu)冗余數(shù)目v=0，機構(gòu)局部自由度數(shù)目ζ=0。將上述參數(shù)代入式(1)有

M=3×(6-7-1)+6×1+2-0=2 。

所以該球面并聯(lián)機構(gòu)的自由度數(shù)目為M=2，并且通過分析可知機構(gòu)具有沿2個方向轉(zhuǎn)動的自由度。驅(qū)動位置選在支鏈與下平臺連接的轉(zhuǎn)動副處，這樣在機構(gòu)運動過程中兩驅(qū)動電機帶動支鏈運動，支鏈就可以帶動上平臺在自由度方向上完成位姿的變換。

2機構(gòu)運動學(xué)分析

2.1　坐標系的建立

按照連桿參數(shù)的D-H表示法建立坐標系：O-X0Y0Z0，O-XYZ和O-XijYijZij(ij表示第i個分支的第j個轉(zhuǎn)動副，以下i，j的取值都為1或2)，如圖3所示。由于兩自由度球面并聯(lián)機構(gòu)各個連桿系都是繞著球心轉(zhuǎn)動，因此選取球心O作為各坐標系的原點；支柱軸線作為固定坐標系的Z0軸，指向上平臺的方向為正方向；各連桿坐標系的Zij分別為第i個分支的第j個轉(zhuǎn)動副的軸線方向，方向指向球外；固定坐標系的X0軸是Z0軸和Z11軸所張成平面的法線；連桿坐標系的Xij軸分別為Zij與Zi，j+1所張成平面的法線方向；將過球心O且垂直于上平臺的垂線作為動坐標系的Z軸；X軸的方向為Z13與Z軸所張成平面的法線方向(與X13軸重合)；各坐標系的Y軸分別由右手定則確定。

圖3　各坐標系的建立Fig.3　Coordinates of the system

該機構(gòu)的結(jié)構(gòu)參數(shù)有α1，α2，β，η1i，η2i。其中η1i(i=1,2)分別為各分支的Zi1在基座的投影與第一分支的Z11投影所成的角度；η2i(i=1,2)分別為各分支的Zi3在上平臺的投影與第一分支的Z13投影所成的角度。

2.2　運動學(xué)逆解的推導(dǎo)

應(yīng)用連桿參數(shù)的D-H表示法[18]及球面解析理論[19]都可以進行機構(gòu)運動學(xué)逆解的推導(dǎo)，此處以連桿參數(shù)的D-H表示法為例進行運動學(xué)逆解推導(dǎo)。由于機構(gòu)的上平臺只有2個轉(zhuǎn)動的自由度，記上平臺位姿為橫滾角φx、俯仰角φy，則上平臺位姿可表示為

(2)

式中：c表示cos函數(shù)運算，s表示sin函數(shù)運算，下同。

2.2.1矢量ui，wi和vi的方向余弦

定義如下各單位向量：ui沿Zi1方向，wi沿Zi2方向，vi沿Zi3方向，以上i=1,2。

通過連桿參數(shù)的D-H表示法進行坐標旋轉(zhuǎn)變換可得固定坐標系O-X0Y0Z0和O-Xi2Yi2Zi2之間的轉(zhuǎn)換矩陣為

其中

(3)

因此，可得

(4)

式中：

a11=cθi2cη1icθi1+cθi2sη1icβsθi1-cα1sθi2cη1isθi1+cα1sθi2sη1icβcθi1-sη1isβsα1sθi2,

a12=-sθi2cη1icθi1-sθi2sη1icβsθi1-cα1cθi2cη1isθi1+cα1cθi2sη1icβcθi1-sη1isβsα1cθi2,

a13=-sα1cη1isθi1+sα1sη1icβcθi1+sη1isβcα1,

a21=cθi2sη1icθi1-cθi2cη1icβsθi1-cα1sθi2sη1isθi1-cα1sθi2cη1icβcθi1+cη1isβsα1sθi2,

a22=-sθi2sη1icθi1+sθi2cη1icβsθi1-cα1cθi2sη1isθi1-cα1cθi2cη1icβcθi1+cη1isβsα1cθi2,

a23=-sα1sη1isθi1-sα1cη1icβcθi1-cα1cη1isβ,

a31=sβsθi1cθi2+sβcθi1cα1sθi2+cβsα1sθi2，

a32=-sβsθi1sθi2+sβcθi1cα1cθi2+cβsα1cθi2,

a33=sβcθi1sα1-cβcα1。

ui=[sη1isβ-cη1isβ-cβ]T，

(5)

同理，根據(jù)上平臺位姿矩陣(2)，結(jié)合坐標的旋轉(zhuǎn)變換可得到vi的表達式：

(6)

2.2.2由機構(gòu)約束條件導(dǎo)出位姿逆解

各分支中間轉(zhuǎn)動副wi與相應(yīng)上平臺轉(zhuǎn)動副vi以上連桿α2相連，就有約束方程

wi·vi=cα2,i=1,2。

(7)

將式(5)和式(6)代入式(7)并整理可得

Aisθi1+Bicθi1+Ci=0,i=1,2，

式中：

Ai=-sα1cη1ivix-sα1sη1iviy,

Bi=sα1sη1icβvix-sα1cη1icβviy+sβsα1viz,

Ci=sη1isβcα1vix-cη1isβcα1viy-cβcα1viz-cα2。

解得

(8)

式(8)表明，已知動平臺的位姿，便可求輸入角θ11和θ21，并且由式(8)可以看出，每個輸出位姿對應(yīng)2組輸入解，根據(jù)機構(gòu)的運動特性易確定符合約束關(guān)系的合理解。

2.2.3逆解算例

球面并聯(lián)機構(gòu)的各設(shè)計參數(shù)分別為α1=25°,α2=30°,β=53.13°,η11=0°,η12=90°,η21=0°,η22=90°。規(guī)劃動平臺運動軌跡為與動平臺固連的動坐標系Z軸繞固定坐標系的Z0軸做錐角為30°的圓錐面勻速運動，如圖4所示。圖中圓形軌跡為動平臺單位法線(與Z軸重合)端點在運動過程中所走過的軌跡，該點也作為動平臺的標記點。

圖4　圓錐面軌跡規(guī)劃Fig.4　Motion planning of the cone trajectory

根據(jù)以上所述的機構(gòu)運動學(xué)逆解過程，對所規(guī)劃的姿態(tài)軌跡進行運動學(xué)逆解求解。圖5所示為通過逆解算法得到的2個驅(qū)動電機的角位移輸入曲線。

圖5　圓錐面軌跡規(guī)劃輸入曲線Fig.5　Input angle curves for the cone trajectory

2.3　運動學(xué)正解的推導(dǎo)

2.3.1雅克比矩陣求解

將運動學(xué)位姿反解的約束方程(7)中的輸入角和輸出角位移均看成時間t的函數(shù)，等式兩邊同時對時間t求導(dǎo)，得：

(9)

用矢量形式將式(9)表示為

Aθ′+Bφ′=0 ，

(10)

對式(10)進行整理得

φ′=Jθ′,

(11)

其中

J=-B-1A,

(12)

則J即為輸入角速度矢量θ′向輸出角轉(zhuǎn)速矢量φ′線性變換的雅克比矩陣。

2.3.2由數(shù)值形式給出機構(gòu)位姿正解

并聯(lián)機構(gòu)位姿正解是機構(gòu)學(xué)研究的一個重要課題，是研究機構(gòu)運動學(xué)、動力學(xué)和軌跡控制等方面的基礎(chǔ)。與串聯(lián)機器人相比，并聯(lián)機構(gòu)反解容易，但正解需要求解一組非線性方程組，難度較大。國內(nèi)外學(xué)者采用不同的方法對正解進行了研究[20-22]，本文用一種常用的最速下降法對機構(gòu)進行位置正解求解。

由于動平臺在其工作空間內(nèi)運動，假定它從任一已知的初始位姿P0經(jīng)過控制作用后，達到當前位姿Pd，此時各分支輸入角位置矢量為θd。在實際過程中P0靠標定來確定，θd則是由分支伺服電機的編碼器讀數(shù)值換算得出。位姿正解算法的根本任務(wù)是根據(jù)P0，θd求解出上平臺位姿Pd。

位姿正解的迭代算法如下：

1)已知上平臺初始位姿P0，分支輸入角位置矢量θd，給定算法的收斂精度ε>0，令k=0；

2)通過運動學(xué)逆解體系計算初始位姿下輸入角位置狀態(tài)：θk=InvKin(Pk)，得Δθk=θd-θk；

3)計算Sk=-J-1(θk)Δθk，若‖Sk‖≤ε，則停，Pd=Pk，否則轉(zhuǎn)4)；

4)Pk+1=Pk+Sk，k=k+1，轉(zhuǎn)2)。

由上述數(shù)值迭代過程即可獲得機構(gòu)的位置正解。

2.3.3正解算例

應(yīng)用2.2.3中的機構(gòu)參數(shù)和位姿逆解結(jié)果，將逆解所得的驅(qū)動電機輸入數(shù)據(jù)做機構(gòu)位姿的正向求解。圖6所示為所述位姿正解算法所得的動平臺橫滾角(Rollx)和俯仰角(Pitchy)的輸出曲線；圖7表示由正解輸出數(shù)據(jù)得到的動平臺標記點所走的軌跡與原規(guī)劃軌跡的對比圖。虛線表示上節(jié)中規(guī)劃的動平臺運動軌跡，星號點線代表通過正解算法獲得的動平臺運動軌跡，2組軌跡完全吻合，驗證了本文提出機構(gòu)的運動學(xué)正逆解算法的正確性和有效性。

圖6　正解輸出曲線Fig.6　Output curve of the forward kinematics

圖7　正解輸出軌跡與原規(guī)劃軌跡的對比圖Fig.7　Comparison between the solutions of forward kinematics and original planned trajectory

3結(jié)語

本文分析了兩自由度球面并聯(lián)機構(gòu)的工作原理，進行了自由度分析和虛擬樣機的設(shè)計。分別以解析形式和數(shù)值形式給出了機構(gòu)運動學(xué)正、逆解結(jié)果，并且通過仿真驗證了所提出正、逆解算法的一致性。分析表明：本機構(gòu)可操控性強，所建立的運動學(xué)模型合理有效，基于運動學(xué)方程進行合理的控制，本機構(gòu)可以在基體與設(shè)備中間實現(xiàn)振動抵消、補償?shù)淖饔?，達到主動隔振的目的。此外，根據(jù)不同的工程應(yīng)用背景，可以針對本文所提出的兩自由度球面并聯(lián)平臺進行不同的參數(shù)設(shè)定，以滿足不同場合的應(yīng)用需求。

參考文獻/References：

[1]張磊, 付永領(lǐng), 劉永光, 等. 主動隔振技術(shù)及其應(yīng)用與發(fā)展[J]. 機床與液壓, 2005 (2): 5-8.

ZHANG Lei, FU Yongling, LIU Yongguang, et al. Active vibration isolation and its application and development[J]. Machine Tool & Hydraulics, 2005(2):5-8.

[2]張春紅, 湯炳新. 主動隔振技術(shù)的回顧與展望[J]. 河海大學(xué)常州分校學(xué)報, 2002, 16(2): 1-5.

ZHANG Chunhong, TANG Bingxin. The retrospection and prospection of active vibration isolaion technique[J].Journal of Hohai University Changzhou, 2002, 16(2): 1-5.

[3]STEWART D.A platform with six degrees of freedom[J]. Proceedings of the Institution of Mechanical Engineer, 1965, 180: 371-386.

[4]HUNT K H. Kinematic Geometry of Mechanisms[M]. Oxford: Clarendon Press, 1978.

[5]王昭. 某天線穩(wěn)定平臺的結(jié)構(gòu)設(shè)計[J]. 電子機械工程, 2012, 28(4): 25-28.

WANG Zhao. Structure design of an antenna stabilizing platform[J]. Electro-Mechanical Engineering, 2012, 28(4): 25-28.

[6]張鵬. 艦載攝像穩(wěn)定平臺的結(jié)構(gòu)設(shè)計[J]. 艦船電子對抗, 2006, 29(3): 70-73.

ZHANG Peng. Structural design of shipboard shooting stabilized platform[J]. Shipboard Electronic Countermeasure, 2006, 29(3): 70-73.

[7]鄒慧君, 高峰. 現(xiàn)代機構(gòu)學(xué)進展[M]. 北京: 高等教育出版社, 2007.

[8]黃真, 孔令富. 并聯(lián)機器人機構(gòu)學(xué)理論及控制[M]. 北京: 機械工業(yè)出版社, 1997.

[9]叢爽, 尚偉偉. 并聯(lián)機器人建模、控制優(yōu)化與應(yīng)用[M]. 北京: 電子工業(yè)出版社，2010.

[10]劉艷芳, 楊隨先. 球面機構(gòu)研究動向[J]. 機械設(shè)計與研究, 2010, 27(1): 32-35.

LIU Yanfang, YANG Suixian. Advances and trends in spherical mechanisms research[J]. Machine Design and Research, 2010, 27(1): 32-35.

[11]GWINNETT J E. Amusement Devices[P]. US Patent:1,789,680,1931-01-20.

[12]方中甲. 面向頸關(guān)節(jié)的三自由度球面并聯(lián)機構(gòu)的研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2004.

FANG Zhongjia. Research on Three DOF Spherical Parallel Mechanism for The Neck Joint[D].Harbin:Harbin Institute of Technology, 2004.

[13]周玉林, 高峰. 具有偏置輸出的3自由度球面并聯(lián)仿生肩關(guān)節(jié)[P]. 中國專利：ZL200810055021.8, 2008-11-12.

[14]GOSSELIN C M,HUMEL J F.The agile eye:A high-performance three-degree-of-freedom camera-orienting device[C]//IEEE International Conference on Robotics and Automation.[S.l.]:[s.n.],1994:781-786.

[15]GOSSELIN C M, ST-PIERRE E, GAGNE M.On the development of the agile eye[J]. IEEE Robotics & Automation Magazine, 1996, 3(4): 29-37.

[16]HUANG Zhen. The kinematics and type synthesis of lower-mobility parallel manipulators[C]//Proceedings of the 11th Word Congress in Mechanism and Machine Science.Tianjin:[s.n.],2004:65-76.

[17]黃真, 劉婧芳, 李艷文. 論機構(gòu)自由度[M]. 北京: 科學(xué)出版社, 2011.

[18]戰(zhàn)麗娜. 基于球面機構(gòu)的肩關(guān)節(jié)的分析與設(shè)計[D].秦皇島:燕山大學(xué),2005.

ZHAN Lina. The Analysis and Design of a Shoulder Joint Based on Spherical Manipulator[D].Qinhuangdao:Yanshan University, 2005.

[19]黃真. 空間機構(gòu)學(xué)[M]. 北京:機械工業(yè)出版社, 1989.

[20]INNOCENTI C, PARENTI-CASTELLI V. A novel numerical approach to the closure of the 6-6 Stewart platform mechanism[C]//International Conference on Advanced Robotics, Robots in Unstructured Environments, 91 ICAR.[S.l.]:IEEE, 1991:851-855.

[21]曲義遠, 黃真. 空間六自由度多回路機構(gòu)位置的三維搜索方法[J]. 機器人, 1989, 3(5): 25-29.

QU Yiyuan, HUANG Zhen. Three-dimensional search method of a six DOF multi-loop mechanism[J]. Robot, 1989, 3(5): 25-29.

[22]李樹軍, 王陰, 王曉光. 3-RPS并聯(lián)機器人機構(gòu)位置正解的桿長逼近法[J]. 東北大學(xué)學(xué)報, 2001, 22(6): 285-287.

LI Shujun, WANG Yin, WANG Xiaoguang. Forward position analysis of 3-RPS in-parallel manipulator using self-modified successive approximation method[J]. Journal of Northeastern University, 2001, 22(6): 285-287.

A two-degree-of-freedom spherical parallel mechanism with low-frequency vibration isolation

ZHANG Xiaowei1, DUAN Xuechao1, LENG Guojun2, BAO Hong1

(1. Key Laboratory of Electronic Equipment Structure Design of Ministry of Education, Xidian University, Xi’an, Shaanxi 710071, China; 2. Southwest China Research Institute of Electronic Equipment, Chengdu, Sichuan 610036, China)

Abstract:Aiming at the requirements of the vibration isolation with two rotational degrees of freedom under specific conditions, a spherical 2-DOF parallel vibration isolation platform is proposed based on the theory of spherical parallel mechanism. From the view point of theory of mechanism, the operating principle is analyzed and the degree of freedom is calculated, and the virtual prototype design of vibration isolation platform is carried out. With D-H representation method of robot mechanism, the constraint equation is developed under the constraint conditions of the above link arc angle, and the inverse kinematics is realized. The kinematics positive solution process is given based on numerical method, and the solution examples of positive and inverse solution are given. Analysis and simulation results show that the 2-DOF spherical parallel mechanism proposed in this paper has the characteristics of good controllability, and the kinematic model is reasonable and effective, meeting the need of low-frequency vibration isolation performance.

Keywords:mechanical principle and mechanism theory; vibration isolation platform; spherical parallel mechanism; Jacobian matrix; kinematics

通訊作者：段學(xué)超博士。E-mail:xchduan@xidian.edu.cn

作者簡介：張曉偉(1988—)，男，河北承德人，碩士研究生，主要從事并聯(lián)機器人技術(shù)方面的研究。

基金項目：國家自然科學(xué)基金(51405362，51490660)

收稿日期：2015-07-10；修回日期：2015-09-30；責任編輯：馮民

中圖分類號：TD421

文獻標志碼：A

doi:10.7535/hbkd.2016yx01003

文章編號：1008-1542(2016)01-0013-07