邱紹楊, 任鴻翔, 尹金崗,2
(1.大連海事大學(xué) 航海動(dòng)態(tài)仿真和控制交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 遼寧 大連 116026;2.中國(guó)海事服務(wù)中心, 北京 100029)
船舶救生系統(tǒng)是為滿足船上人員在海難中自救的需要而設(shè)置在船上的專門設(shè)備及其附屬裝置的總稱。國(guó)際海事組織在《國(guó)際海上人命安全公約》中要求船舶必須配備救生設(shè)備的同時(shí),還在《船員培訓(xùn)值班和發(fā)證國(guó)際公約》中要求進(jìn)行船舶消防、救生艇筏操縱、海上求生和海上急救等4項(xiàng)專業(yè)訓(xùn)練。目前的船舶救生培訓(xùn)大多在真實(shí)設(shè)備上進(jìn)行,危險(xiǎn)性高。此外,受時(shí)間、設(shè)備數(shù)量和環(huán)境等因素限制,這些培訓(xùn)的效率較低。
近年來(lái),隨著虛擬現(xiàn)實(shí)技術(shù)的快速發(fā)展,虛擬現(xiàn)實(shí)在航海教學(xué)和培訓(xùn)領(lǐng)域引起廣泛關(guān)注[1-2],利用該技術(shù)開發(fā)具有逼真沉浸感和良好交互性的虛擬救生培訓(xùn)系統(tǒng)可較好地解決上述問(wèn)題。本文基于3ds Max和Unity 3D研發(fā)虛擬船舶救生培訓(xùn)系統(tǒng),提供船舶救生設(shè)備的模擬訓(xùn)練功能。
船舶救生設(shè)備主要包括共用救生設(shè)備、個(gè)人救生設(shè)備、救生視覺(jué)信號(hào)設(shè)備和救生通信與報(bào)警系統(tǒng)。救生培訓(xùn)涉及了解救生設(shè)備的主要結(jié)構(gòu)、熟悉設(shè)備的使用和操作流程及救生演習(xí)等內(nèi)容。救生演習(xí)包括棄船和人員落水,通常按照一定的流程,綜合多種救生設(shè)備的操作使用。
本文利用虛擬現(xiàn)實(shí)技術(shù)開發(fā)船舶救生培訓(xùn)系統(tǒng),構(gòu)建海況良好、船舶火災(zāi)、沉船、大風(fēng)浪天氣和人員落水等條件下的海上救生演習(xí)虛擬場(chǎng)景,在這些場(chǎng)景中實(shí)現(xiàn)救生通信和報(bào)警、救生設(shè)備操作和救生演習(xí)協(xié)同訓(xùn)練;受訓(xùn)人員在訓(xùn)練過(guò)程中能熟悉不同情況下救生設(shè)備的操作要領(lǐng)和注意事項(xiàng)[3],并提高團(tuán)隊(duì)的緊密協(xié)作能力。
系統(tǒng)由1個(gè)指揮決策中心和若干個(gè)訓(xùn)練單元組成(見圖1),采用模塊化程序設(shè)計(jì)思想,分成不同模塊進(jìn)行開發(fā),包括:指揮決策模塊、系統(tǒng)交互模塊、運(yùn)動(dòng)模型模塊和救生通信與報(bào)警系統(tǒng)模塊。各模塊的功能如下。
1.2.1指揮決策模塊
該模塊包括訓(xùn)練任務(wù)初始化、訓(xùn)練過(guò)程監(jiān)管和訓(xùn)練結(jié)果評(píng)判。
(1) 訓(xùn)練任務(wù)初始化,包括: 訓(xùn)練場(chǎng)景設(shè)定,即火災(zāi)的位置和態(tài)勢(shì)、沉船的速度、天氣狀況、人員落水的位置和船舶狀態(tài)等; 任務(wù)的分配,即船舶遇險(xiǎn)后每個(gè)受訓(xùn)人員的任務(wù); 救生設(shè)備的設(shè)定,即船舶遇險(xiǎn)后需用到哪些設(shè)備。
(2) 訓(xùn)練過(guò)程監(jiān)管,即訓(xùn)練時(shí)指揮決策模塊調(diào)整場(chǎng)景視角觀察救生現(xiàn)場(chǎng)情況并錄像,對(duì)整個(gè)過(guò)程進(jìn)行監(jiān)管并實(shí)時(shí)記錄數(shù)據(jù)。
(3) 訓(xùn)練結(jié)果評(píng)判,即訓(xùn)練結(jié)束之后對(duì)記錄的信息數(shù)據(jù)進(jìn)行讀取和回放,根據(jù)訓(xùn)練時(shí)記錄下的用戶操作順序、狀態(tài)等信息對(duì)其相應(yīng)業(yè)務(wù)處理能力進(jìn)行評(píng)判。
1.2.2系統(tǒng)交互模塊
該模塊可實(shí)現(xiàn)船舶主要救生設(shè)備的虛擬現(xiàn)實(shí)交互和相應(yīng)的應(yīng)急模擬,包括:救生艇筏、救助艇以多種方式釋放及回收;救生艇在水面操縱、救助艇在水中救人,并與大船互動(dòng);拋繩器的使用、救生圈的拋放;救生衣、浸水衣的穿戴;救生視覺(jué)信號(hào)設(shè)備的釋放和通信報(bào)警設(shè)備的使用等。
1.2.3運(yùn)動(dòng)模型模塊
該模塊主要指虛擬母船、救生艇筏、救助艇、柔性繩索和艇架吊臂的三維運(yùn)動(dòng)模型。運(yùn)動(dòng)模型包括影響運(yùn)動(dòng)的各種效應(yīng),主要有:主機(jī)、舵、纜繩、側(cè)推對(duì)船和艇的控制;風(fēng)和流對(duì)船、艇、筏及柔性繩索的影響;船與艇筏之間的效應(yīng)。
1.2.4救生通信與報(bào)警系統(tǒng)模塊
該模塊模擬船舶救生時(shí)的通信與報(bào)警系統(tǒng),主要包括全船警報(bào)、救生現(xiàn)場(chǎng)人員之間、現(xiàn)場(chǎng)人員與駕駛臺(tái)之間及救生艇筏和救助艇與大船之間的通信。報(bào)警系統(tǒng)主要包括船舶VHF/MF/HF DSC和Inmarsat設(shè)備及棄船時(shí)攜帶至救生艇筏上的VHF、SART和EPIRB。
為提高救生艇、救助艇水面運(yùn)動(dòng)仿真的逼真度,建立這2種艇的四自由度運(yùn)動(dòng)數(shù)學(xué)模型。
2.1.1坐標(biāo)系
采用大地坐標(biāo)系和附體坐標(biāo)系(見圖2),其中:O0x0y0z0為固定在地球表面的大地坐標(biāo)系,O0x0軸指向正北,O0y0軸指向正東,O0z0軸指向地心;Oxyz為原點(diǎn)位于船上某定點(diǎn)O上的附體坐標(biāo)系,Ox軸指向艏部,Oy軸指向右舷,Oz軸指向船體下方。
2.1.2四自由度運(yùn)動(dòng)方程的建立
基于分離模型(MMG)的思想計(jì)算力和力矩[4],坐標(biāo)原點(diǎn)設(shè)置在重心位置處,由于小艇的螺旋槳為導(dǎo)管槳,因此省去舵的力和力矩,推導(dǎo)出小艇四自由度動(dòng)力學(xué)方程為
(1)
式(1)中:H和P分別為作用在小艇船體和螺旋槳的力和力矩;wind和wave分別為風(fēng)、浪的力矩;m為小艇質(zhì)量;mx和my分別為縱向和橫向小艇的附加質(zhì)量;Ixx和Izz分別為繞Ox軸和Oz軸的轉(zhuǎn)動(dòng)慣量;Jxx和Jzz分別為繞Ox軸及Oz軸的附加轉(zhuǎn)動(dòng)慣量;u、v、p和r分別為縱向、橫向和繞Ox軸、Oz軸的角速度。
運(yùn)動(dòng)學(xué)方程為
(2)
式(2)中:x和y分別為縱向和橫向的位移;φ和ψ分別為縱向和垂向的歐拉角。
作用在小艇上的力和力矩的計(jì)算見文獻(xiàn)[4]~文獻(xiàn)[6],波浪和風(fēng)的干擾力計(jì)算見文獻(xiàn)[7]。利用四階龍格庫(kù)塔方法求解式(1)和式(2)。
共用救生設(shè)備的操作和救生演習(xí)需多人協(xié)作完成,該系統(tǒng)將所需的人分為不同的角色,受訓(xùn)人員先選擇角色,再進(jìn)行訓(xùn)練。為提高系統(tǒng)的靈活性,設(shè)計(jì)“真”人與虛擬人之間的協(xié)同和多人之間的協(xié)同2種工作方式。“真”人與虛擬人之間的協(xié)同是指每個(gè)受訓(xùn)人員選擇一個(gè)角色,另外的角色為虛擬輔助船員,通過(guò)系統(tǒng)提供的命令菜單來(lái)控制輔助船員操作設(shè)備,受訓(xùn)人員與輔助船員“配合”完成相關(guān)操縱;之間的協(xié)同是指系統(tǒng)中的角色全部為受訓(xùn)人員,受訓(xùn)人員通過(guò)局域網(wǎng)協(xié)同工作,并在虛擬場(chǎng)景中互見,多人協(xié)同完成相關(guān)訓(xùn)練。
在多人協(xié)同操作時(shí),同一時(shí)刻受訓(xùn)人員可能同時(shí)操作并控制1臺(tái)設(shè)備,此時(shí)若處理不當(dāng),可能會(huì)給系統(tǒng)帶來(lái)不穩(wěn)定因素,從而導(dǎo)致系統(tǒng)運(yùn)行異常,因此需有一個(gè)穩(wěn)健的協(xié)同仿真策略。[8]文獻(xiàn)[9]中的系統(tǒng)采用定時(shí)采集和發(fā)送數(shù)據(jù)的方法, 每隔一段時(shí)間采集一次系統(tǒng)中的數(shù)據(jù),進(jìn)行仿真計(jì)算。為提高系統(tǒng)的實(shí)時(shí)性,該系統(tǒng)每幀進(jìn)行一次數(shù)據(jù)采集,在每次仿真計(jì)算之前依次獲取各訓(xùn)練單元傳來(lái)的數(shù)據(jù),形成一幀完整的最新數(shù)據(jù),然后進(jìn)行仿真計(jì)算。在該情況下,若系統(tǒng)仍出現(xiàn)異常,會(huì)根據(jù)訓(xùn)練單元角色級(jí)別的高低,選擇一個(gè)訓(xùn)練單元的操作為優(yōu)先操作,優(yōu)先級(jí)順序?yàn)椋捍L(zhǎng)、大副、二副、三副和水手長(zhǎng)等,直到系統(tǒng)異常消失。
該系統(tǒng)為多訓(xùn)練單元系統(tǒng),具有多計(jì)算機(jī)、多顯示和多音頻設(shè)備的特點(diǎn),只有當(dāng)受訓(xùn)人員各自計(jì)算機(jī)上的場(chǎng)景一致時(shí),才能進(jìn)行協(xié)同操作訓(xùn)練。為實(shí)現(xiàn)場(chǎng)景的一致性和“真”人與虛擬人操作流程的一致性,系統(tǒng)構(gòu)建一個(gè)虛擬場(chǎng)景數(shù)據(jù)庫(kù),各訓(xùn)練單元直接與數(shù)據(jù)庫(kù)進(jìn)行數(shù)據(jù)交換。當(dāng)訓(xùn)練單元與系統(tǒng)交互時(shí),直接改變數(shù)據(jù)庫(kù)中的數(shù)據(jù),其他訓(xùn)練單元從數(shù)據(jù)庫(kù)中更新本機(jī)的數(shù)據(jù)即可實(shí)現(xiàn)場(chǎng)景的一致性。
在利用脫鉤手柄釋放救生艇之后,救生艇依靠重力沿滑道下滑;在脫離軌道之后,以拋物線的軌跡在空中滑行,最終入水。[10-11]為提高救生艇自由降落過(guò)程仿真的逼真度,對(duì)其建立數(shù)學(xué)模型。救生艇自由降落入水之前的過(guò)程見圖3。救生艇入水之前的運(yùn)動(dòng)分為下滑、旋轉(zhuǎn)和自由下落等3個(gè)階段。由于整個(gè)降落過(guò)程的時(shí)間較短,因此在該段時(shí)間內(nèi),救生艇滑道可視為是靜止不動(dòng)的,救生艇所受空氣阻力對(duì)救生艇運(yùn)動(dòng)的影響較小,可忽略不計(jì)。
1) 下滑:下滑過(guò)程是從救生艇被釋放到其重心靠近滑架最低點(diǎn),由位置1到位置2的過(guò)程,此時(shí)救生艇受重力mg、滑道支持力Fn及摩擦力f=μFn的作用,μ=tanΦ為摩擦系數(shù),運(yùn)動(dòng)方程為
(3)
式(3)中:I為轉(zhuǎn)動(dòng)慣性矩;θ為救生艇軸線與水平方向的夾角。
2) 旋轉(zhuǎn):旋轉(zhuǎn)過(guò)程是下滑結(jié)束到救生艇不再接觸滑架的過(guò)程,救生艇受力與下滑過(guò)程相同。此時(shí)救生艇受的重力和滑道作用于救生艇的力的作用線不在同一條直線上,使救生艇發(fā)生轉(zhuǎn)動(dòng),運(yùn)動(dòng)方程為
(4)
3) 自由下落:自由下落過(guò)程是從轉(zhuǎn)動(dòng)結(jié)束到救生艇接觸水面的過(guò)程,此時(shí)救生艇只受重力,但仍保持轉(zhuǎn)動(dòng)狀態(tài),直到接觸水面,運(yùn)動(dòng)方程為
(5)
系統(tǒng)訓(xùn)練與評(píng)估的流程為:
1) 教練員在教練員站初始化訓(xùn)練任務(wù)。
2) 系統(tǒng)將訓(xùn)練任務(wù)轉(zhuǎn)換為訓(xùn)練任務(wù)數(shù)據(jù),發(fā)送到各訓(xùn)練單元。
3) 學(xué)員在訓(xùn)練單元按照訓(xùn)練流程進(jìn)行訓(xùn)練,訓(xùn)練單元記錄操作數(shù)據(jù)。
4) 訓(xùn)練結(jié)束之后,訓(xùn)練單元將操作數(shù)據(jù)傳回教練員站,教練員站根據(jù)訓(xùn)練任務(wù)數(shù)據(jù)、操作數(shù)據(jù)及訓(xùn)練錄像進(jìn)行評(píng)判。
訓(xùn)練流程包括設(shè)備操作和指揮命令,指揮命令又包含操作命令和訓(xùn)練點(diǎn)評(píng)。訓(xùn)練點(diǎn)評(píng)的實(shí)現(xiàn)方案見圖4,系統(tǒng)以訓(xùn)練任務(wù)數(shù)據(jù)和操作記錄數(shù)據(jù)為基礎(chǔ),根據(jù)系統(tǒng)評(píng)判或主觀判斷得出訓(xùn)練過(guò)程的不足和需注意的事項(xiàng),最終給出點(diǎn)評(píng)。點(diǎn)評(píng)方式為:
1) 當(dāng)點(diǎn)評(píng)人為第一人稱時(shí),點(diǎn)評(píng)人根據(jù)主觀判斷,通過(guò)語(yǔ)音設(shè)備和鍵盤進(jìn)行語(yǔ)音及文字點(diǎn)評(píng)。
2) 當(dāng)點(diǎn)評(píng)人為第三人稱時(shí),系統(tǒng)通過(guò)評(píng)價(jià)指標(biāo)進(jìn)行計(jì)算,給出語(yǔ)音及文字點(diǎn)評(píng)。
對(duì)于訓(xùn)練評(píng)估方案(見圖5),評(píng)估方式包括機(jī)器評(píng)估、人工評(píng)估和專家評(píng)估等3種。機(jī)器評(píng)估和人工評(píng)估分別對(duì)應(yīng)設(shè)備操作及指揮命令,專家評(píng)估起到輔助評(píng)判的作用,最終修正、完善系統(tǒng)的評(píng)估。當(dāng)前系統(tǒng)的評(píng)價(jià)指標(biāo)是可量化的評(píng)價(jià)指標(biāo),評(píng)估方法為專家法與隸屬度函數(shù)相結(jié)合的方法[12],以后將采用更智能的基于云模型和機(jī)器學(xué)習(xí)的方法。
以巴拿馬型散貨船“長(zhǎng)山?!陛啚槟感痛?,利用三維模型制作軟件3ds Max建立船體、駕駛臺(tái)、救生艇筏、救助艇及其存放和登乘裝置、個(gè)人救生設(shè)備、救生視覺(jué)信號(hào)設(shè)備、救生通信與報(bào)警系統(tǒng)、救生拋繩設(shè)備的模型。將模型導(dǎo)入到Unity 3D中,構(gòu)建虛擬海上救生演習(xí)場(chǎng)景,利用音頻制作技術(shù)實(shí)現(xiàn)場(chǎng)景中船和艇的發(fā)動(dòng)機(jī)、海水、船員行走和報(bào)警等所發(fā)出聲音的模擬;利用粒子系統(tǒng)實(shí)現(xiàn)場(chǎng)景中火災(zāi)、艇入水時(shí)水花、漂浮煙霧信號(hào)、手持火焰信號(hào)、火箭降落傘信號(hào)和救生圈自發(fā)煙霧信號(hào)的模擬;利用場(chǎng)區(qū)網(wǎng)模型實(shí)現(xiàn)船舶火災(zāi)蔓延過(guò)程模擬[13];利用UGUI開發(fā)用戶界面。
本文的救生培訓(xùn)系統(tǒng)可模擬船舶火災(zāi)、沉船、大風(fēng)浪天氣和人員落水條件下進(jìn)行海上救生演習(xí)的虛擬場(chǎng)景,下面以人員落水、船舶失火和大風(fēng)浪天氣場(chǎng)景的救生演習(xí)為例進(jìn)行說(shuō)明。
3.2.1人員落水場(chǎng)景
圖6為人員落水演習(xí)。受訓(xùn)人員發(fā)現(xiàn)有人落水之后,立即報(bào)告駕駛臺(tái)并拋出救生圈(見圖6a)),值班駕駛員立即向落水者一舷操滿舵并發(fā)出人員落水警報(bào)(三長(zhǎng)聲);船長(zhǎng)到駕駛臺(tái)操縱船舶,設(shè)法接近落水者,并下達(dá)釋放救助艇救人的命令;大副收到命令之后組織人員到艇甲板集合,受訓(xùn)人員依照分配的任務(wù)釋放救助艇,系統(tǒng)能檢測(cè)船體與救助艇艇之間的碰撞(見圖6b)),受訓(xùn)人員操作失誤,導(dǎo)致救助艇撞到大船;救助艇釋放之后,受訓(xùn)人員操救助艇,使其從落水者下風(fēng)側(cè)接近落水者(見圖6c));落水人員被救起之后,回收救助艇,并向船長(zhǎng)匯報(bào);船長(zhǎng)解除警報(bào),組織講評(píng),演習(xí)結(jié)束。
3.2.2船舶失火場(chǎng)景
圖7為船舶失火演習(xí)。在駕駛臺(tái)發(fā)現(xiàn)貨艙發(fā)生大火(見圖7a)),火勢(shì)不可控,發(fā)出全船棄船警報(bào)(七短聲一長(zhǎng)聲);利用報(bào)警系統(tǒng)發(fā)出遇險(xiǎn)警報(bào),利用EPIRB和SART報(bào)警(見圖7b)和圖7c));船長(zhǎng)命令釋放救生艇和救生筏,人員集合,按照分配的任務(wù)釋放艇筏,隨后通過(guò)救生視覺(jué)信號(hào)設(shè)備求救(見圖7d)),在海面操縱救生艇(見圖7e));回收救生艇;駕駛臺(tái)解除警報(bào),現(xiàn)場(chǎng)組織講評(píng),演習(xí)結(jié)束。
3.2.3大風(fēng)浪天氣場(chǎng)景
圖8為大風(fēng)浪天氣演習(xí)。在大風(fēng)浪中,采用自由降落方式釋放救生艇(見圖8a),船舶頂流,控制航向和航速,艇長(zhǎng)選擇合適的時(shí)機(jī)釋放救生艇;艇入水之后,小心操縱,使其頂浪航行;船長(zhǎng)操縱船舶使風(fēng)舷角在30°~40°,減速,將救助艇的一舷變?yōu)橄嘛L(fēng)舷,救助艇的釋放受風(fēng)浪影響較大。在圖8b)中,由于救助艇未系艏纜和艉部止蕩繩索,導(dǎo)致救助艇搖晃劇烈,由于該救助艇處沒(méi)有羊角,因此只能在欄桿上綁系止蕩繩索來(lái)防止救助艇劇烈晃動(dòng)。在圖8c)中,艇入水前先啟動(dòng)艇機(jī),可先解開艉止蕩繩索,再解開艏纜,盡快駛離大船。
本文將船舶救生與虛擬現(xiàn)實(shí)緊密結(jié)合,使船舶救生訓(xùn)練不再局限于真實(shí)設(shè)備。在搭建船舶救生培訓(xùn)系統(tǒng)過(guò)程中,為救生艇和救助艇建立運(yùn)動(dòng)數(shù)學(xué)模型,模型滿足系統(tǒng)對(duì)艇在水面運(yùn)動(dòng)的需求;為船舶救生訓(xùn)練設(shè)計(jì)2種工作模式,提高了系統(tǒng)的靈活性;設(shè)計(jì)多人協(xié)同訓(xùn)練的仿真策略,保證了系統(tǒng)的穩(wěn)定性。經(jīng)測(cè)試,系統(tǒng)的實(shí)時(shí)性和沉浸感較好,培訓(xùn)效果較理想,能起到輔助訓(xùn)練的作用,具有一定的價(jià)值。
以后的研究將進(jìn)一步完善艇入水的數(shù)學(xué)模型、艇與大船之間的船間效應(yīng)和耦合運(yùn)動(dòng)及多角色協(xié)同訓(xùn)練的智能評(píng)估算法。
[1] 劉晶晶,任鴻翔,尹金崗,等.多平臺(tái)的船用陀螺羅經(jīng)交互仿真系統(tǒng)[J]. 大連海事大學(xué)學(xué)報(bào), 2016, 42(1):17-20.
[2] 邱紹楊,任鴻翔,尹金崗.基于虛擬現(xiàn)實(shí)技術(shù)的自由降落式救生艇培訓(xùn)系統(tǒng)[J]. 大連海事大學(xué)學(xué)報(bào),2017,43(3):14-18.
[3] 李同欽. 精通救生艇筏和救助艇[M].大連:大連海事大學(xué)出版社,2009:63-70.
[4] 賈欣樂(lè),楊鹽生. 船舶運(yùn)動(dòng)數(shù)學(xué)模型[M].大連:大連海事大學(xué)出版社,1999:49-138.
[5] 王貴彪,謝永和,許頌捷.導(dǎo)管螺旋槳四象限水動(dòng)力性能研究[J]. 船舶工程, 2015, 37(10): 26-28.
[6] LEGURSKY K. A Modified Model, Simulation, and Tests of a Full-Scale Sailing Yacht [C]//Oceans. USA: IEEE, 2012: 1-7.
[7] 張秀鳳,尹勇,金一丞.規(guī)則波中船舶運(yùn)動(dòng)六自由度數(shù)學(xué)模型[J].交通運(yùn)輸工程學(xué)報(bào),2007,7(3):40-43.
[8] 陳學(xué)文,劉玉慶,朱秀慶,等. 航天員太空協(xié)同操作虛擬訓(xùn)練仿真系統(tǒng)研究[J]. 系統(tǒng)仿真學(xué)報(bào),2013,25(10):2348-2354.
[9] 許愛軍.多用戶協(xié)同虛擬現(xiàn)實(shí)的技能訓(xùn)練系統(tǒng)[J]. 計(jì)算機(jī)系統(tǒng)應(yīng)用, 2016,25(2):65-70.
[10] 黃春平,陸軍,秦嘉岷. 自由降落救生艇試驗(yàn)過(guò)程探討[J]. 造船技術(shù),2012(3): 37-40.
[11] ARAI M, KHONDOKER M R H, INOUE Y. Water Entry Simulation of Free-Fall Lifeboat[J]. Journal of the Society Naval Architects of Japan, 1995,178: 193-201.
[12] 王德龍,任鴻翔,肖方兵. 船舶操縱模擬器單船避碰自動(dòng)評(píng)估系統(tǒng)[J].中國(guó)航海, 2015, 38(1): 44-48.
[13] 任鴻翔,金一丞,尹勇. 船舶火災(zāi)模擬訓(xùn)練系統(tǒng)研究[J]. 武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版), 2010,34(1): 19-22.