工艺与材料 EMC
A2008,35(4)
钕铁硼永磁材料热性能的分析
姚丙雷, 林 岩, 刘秀芹
1
2
1
(1.上海电器科学研究所(集团)有限公司电机分所,上海 200063;
2.LG电子(天津)电器有限公司研发中心,天津 300402)
摘 要:介绍了钕铁硼永磁材料的热稳定性。详细说明了钕铁硼永磁材料的测试方法,并对测试结果进行了分析。得出钕铁硼永磁材料在100℃及120℃时,退磁曲线的拐点位置分布情况以及永磁材料内禀矫顽力温度系数随温度变化的趋势,对永磁电机的设计具有参考价值。
关键词:钕铁硼永磁材料;高温磁性能;不可逆退磁
中图分类号:TM273∶TM201.4+2 文献标识码:A 文章编号:1673-6540(2008)04-0052-04
AnalysisofThermalPerformanceof
NdFeBPermanentMagnet
YAOBing-lei, LINYan, LIUXiu-qin
(1.ShanghaiElectricalApparatusResearchInstitute(Group)Co.,Ltd,Shanghai200063,China;
2.LG(Tianjin)ElectronicsApparatusCo.,Ltd,Tianjin300402,China)
Abstract:ThermalstabilityofNd-Fe-Bpermanentmagnetwasintroduced.ThetestmethodsofNd-Fe-Bperma-nentmagnetmaterialweredescribedindetailandtestresultswereanalyzed.ThedistributionofkneepointofNd-Fe-Bpermanentmagnetmaterialdemagnetizationcurveon100℃and120℃wasobtainedandthechangingtrendofin-herentcoerciveforcetemperaturecoefficientastemperaturechangingwasachieved.Thesereferencevaluesarehelpfultothepermanentmagnetmachinedesign.
Keywords:Nd-Fe-Bpermanentmagnetmaterial;magnetismperformanceinthermalstate;irreversibledemagnetization
1
2
1
0 引 言
随着永磁材料性能的不断提高,稀土永磁材料被广泛应用于电机领域
[1]
决定永磁电机能否可靠运行。当永磁电机运行时,有时会出现转子温升较高而使钕铁硼永磁材
料发生不可逆退磁的现象,从而引起永磁同步电动机效率和功率因数等性能指标变差。因此,有必要对钕铁硼材料的热稳定性进行研究、分析。
本文采用测量的分析方法,在大量的试验结果基础上进行数据分析,目的在于获得高温时同一永磁材料各部分性能和同一批永磁材料性能的测试数据,并对测试数据以及由此对永磁材料退磁情况的影响进行分析,为永磁电机的设计提供参考依据。
。采用高性能的稀
土永磁材料制作电机,在结构上不需要专门的电励磁系统,有利于提高电机的工作效率和功率因数,并能节省能源,减小电机的体积
[2]
。
目前国内在设计永磁电机时,一般假设永磁体各部分磁场分布均匀,且同一批次永磁体的磁性能相同,同时需要利用永磁体供应商提供的常温下的永磁体退磁曲线、内禀矫顽力、剩磁密度和最大磁能积等磁性能参数值。但是,由于工厂制造工艺和水平等原因,在高温下同一块永磁体各部分实际的性能存在差异,且同一批永磁体的性能也是不同的。而永磁电机正常运行时的性能与永磁体的磁性能密切相关,永磁体的磁性能直接
1 永磁材料的热稳定性
温度稳定性是指永磁体由所处环境温度的改变而引起磁性能变化的程度。磁稳定性是表示在施加外加磁场条件下,引起磁性能发生变化的程
2008,35(4)
[3]
工艺与材料 EMCA
度。永磁电机在运行时,由于电机温升的影响100℃时有一块切割标准样拐点位置达到24.4%,120℃时拐点位置高达49.7%。拐点位
置的差异容易发生永磁体局部失磁。
表1 N35-SH-A在100℃和120℃时的温度特性
α(B)/r(%/K)100℃
N35-SH-A-1N35-SH-A-2N35-SH-A-5
-0.09-0.1-0.08
120℃-0.09-0.1-0.11
α(H)/cJ(%/K)100℃-0.73-0.70-0.70
拐点处磁密与剩磁密度比值/(%)
以及电枢绕组的存在,永磁电机内永磁材料的性
能由温度和磁场两个因素共同影响。因此,本文的分析中将温度稳定性和磁稳定性统称为热稳定性。即对永磁材料加热至不同的温度分别施加反向外加磁场时,永磁体所表现的磁性能的变化,也就是不同温度下永磁材料的剩磁密度的温度系数和内禀矫顽力的温度系数的大小以及退磁曲线上拐点位置的变化。
试样
120℃100℃120℃-0.67-0.67-0.66
19.618.623.4
26.532.730.4
2 测试方法及结果分析
测量仪器为稀土永磁检测专用设备,其测试
原理是通过测得的J-H曲线来计算B-H曲线。首先对试样进行充磁,然后对检测设备进行初始化,最后将试样装好,并将其加热到预定的温度,起动测试程序进行测试,测得的曲线由计算机记录并打印。
测试工作主要从以下几个方面展开:(1)对一块永磁体不同位置切割标准样进行测试;
(2)对同一厂家、同一牌号的不同永磁体(主要是 10×10mm标准样)进行测试;(3)对同为SH牌号的大块样品进行检测。2.1 同一样品不同位置热稳定性的检测分析
将2块永磁体N35-SH-A和N35-SH-B试样分别切割成8块标准试样(切割方法如图1所示),从所切割的8块标准试样中选取有代表性(不同位置)的3块(1、2、5)分别进行测试,并计算其在100℃和120℃时的温度特性,分析结果见表1、2所示
。
表2 N35-SH-B在100℃和120℃时的温度特性
α(B)/r(%/K)100℃
N35-SH-B-1N35-SH-B-2N35-SH-B-5
-0.06-0.07-0.06
120℃-0.09-0.08-0.08
α(H)/cJ(%/K)100℃-0.69-0.63-0.68
拐点处磁密与剩磁密度比值/(%)
试样
120℃100℃120℃-0.67
24.4
49.724.325.2
-0.61不弯-0.64不弯
2.2 同厂家同牌号不同永磁样品热稳定性的测
试分析
分别对来自2个厂家的、不同永磁标准样品的、高温下的性能进行测试分析,结果见表3、4所示。
表3 a厂4块样品高温下测试结果
α(B)/r(%/K)100℃
a1a2a3a4
-0.11-0.12-0.11-0.11
120℃-0.11-0.11-0.11-0.11
α(H)/cJ(%/K)100℃-0.67-0.67-0.67-0.65
120℃-0.63-0.63-0.62-0.62
拐点处磁密与剩磁密度比值/(%)100℃不弯不弯不弯不弯
120℃14.516.314.512.2
试样
图1 切割方法示意图
从分析数据可看出:在高温情况下,同一磁块切割下来的不同标准试样,在标志永磁体热稳定性好坏的2个重要参数(温度系数和拐点位置)上有着很大的差别。如:N35-SH-B的3个样品,
从对钕铁硼标准样品高温下的α(B、r)α(H)和拐点位置这三个参数的检测结果可得出cJ以下结论。
(1) α(B) 的变化范围不大,基本都在r0.11±0.01%/K范围内。2个厂家的8块样品
—
工艺与材料 EMC
A
表4 b厂4块样品高温下测试结果
试样
α(B)/r(%/K)100℃
b1b2b3b4
-0.10-0.10-0.10-0.08
120℃-0.11-0.10-0.11-0.10
α(H)/cJ(%/K)100℃-0.55-0.56-0.56-0.59
120℃-0.56-0.57-0.55-0.57
拐点处磁密
与剩磁密度比值/(%)100℃不弯不弯不弯10.1
120℃2.986.004.7023.30
2008,35(4)
表5 各厂家SH牌号NdFeB永磁体在100℃
和120℃时拐点位置情况
永磁体生产厂家
100℃时的拐点(B)/%k/Br
120℃时的拐点120℃(B/B)/%kr
2块永磁体拐点出现在
8块永磁体没有出现拐点,5块永磁体
甲厂(13)出现拐点,其拐点分
别为(15.87、22.00、24.78、6.09、5.08)
10%以下,11块永磁体出现拐点,其拐点分别为(18.94、17.33、19.43、20.34、21.46、13.05、26.70、34.60、37.00、21.10、21.90)
8块永磁体没有出现拐点,4块永磁体
乙厂(12)出现拐点,其拐点分
别为(11.71、8.29、7.89、7.19)1块永磁体没有出现拐点,4块永磁体
丙厂(5)出现拐点,其拐点分
别为(17.40、9.22、7.45、4.42)
3块永磁体拐点出现在10%以下,9块永磁体出现拐点,其拐点分别(16.30、15.25、
17.95、
10.73、
10.86、24.60、21.40、11.00、24.00)
1块永磁体拐点出现在10%以下,4块永磁体出现拐点,其拐点分别(34.1、23.2、27.3、22.3
)
中, α(B 最大的也仅为0.12%/K,最小为r)0.08%/K。
(2)2个厂家的α(H检测结果各有不同。cJ)a厂和b厂样品的α(H)较为稳定,但a厂的数cJ值均普遍偏高,120℃时, α(H均高于cj) 0.6%/K,变动范围在0.61～0.63%/K之间;b厂则在0.55～0.57%/K范围以内。
(3)从检测结果来看,2个厂家的永磁材料在高温下的拐点位置普遍偏高。a厂4块样品虽然在100℃时都为直线,但在120℃时拐点位置均高于10%,且最大值为16.3%;b厂4块样品中有3块满足100℃时为直线、120℃拐点位置低于10%,但是b4号样品与前3块偏差较大,100℃时拐点位置已高于10%,120℃时更达到了23.3%,高出120℃要求的1倍多。
2.3 大块样品的热稳定性测试结果与分析
本次试验是以内置式永磁电机中的大块永磁体为试样进行测试的。表5是永磁样品在100℃和120℃时,退磁曲线上拐点处的Bk值占Br的百分比的统计情况。
从对不同厂家的钕铁硼永磁材料的试验结果可以看出:
(1)对于每块试样来说,随着试样温度的升高,其B温度系数 α(B) 呈现上升趋势,H温rrcJ度系数 α(H 呈现下降趋势,图2是两种样品cJ)的 α(B) 和 α(H) 变化趋势图;rcJ
(2)从已检测的试样来看,大部分试样都未满足电机行业的要求(如SH牌号永磁体要求100°C时退磁曲线是直线;120°C时,退磁曲线拐点低于10%B。其中100°C时为直线的仅占r)56.7%,120°C拐点低于10%Br的仅占20%。100℃和120℃时退磁曲线拐点位置的分布情况如图3所示;
图2 两种样品 α(B) 和 α(H) 随温度变化趋势rcJ
(3)所测样品的内禀矫顽力的温度系数α(H)普遍偏高。从图4可以看出, α(H
2008,35(4)工艺与材料 EMC
A
63.2%; α(H >0.7%/K的占24.5%,最大cJ)为0.769%/K。这显然对防失磁是不利的。
3 结 语
本文在对大量钕铁硼永磁材料检测的基础上,对检测结果进行了分析。测试结果表明:目前仅根据常温下的矩形度H、H值以及BK/HcJcJr值并不能判断产品的热稳定性是否合格。如SH牌
号的NdFeB永磁体在矩形度大于或等于0.9、HA(20kOe)的情况下,并不能保证120cJ≥1600k
°C时拐点小于10%。究其原因,钕铁硼永磁材料生产企业没有把内禀矫顽力的温度系数α(H)作cJ为验收标准。总之,本文的分析对永磁电机的设计具有指导意义,有较好的工程利用价值,也为正
图3 100°C和120°C
时拐点位置分布图
确选用钕铁硼永磁材料以防止永磁电机高温退磁提供了依据,具有实际意义。
【参考文献】
[1] 彭龙,张怀武,宋远强,等.永磁材料在电机中的应
用现况分析[J].新材料产业,2006(11):30-33.[2] 王秀和.永磁电机[M].北京:中国电力出版社,
2007.
[3] 唐任远.现代永磁电机理论与设计[M].北京:机
械工业出版社,1997.
收稿日期:2008-01-10
图4 所测样品的 α(H) 分布情况cJ
(上接第51页
)
图4 系统网络总结构框图
[3] 上海交通大学罗克韦尔自动化实验室.罗克韦尔自
动化网络控制平台用户手册[G].2000.
[4] 杜清珍,李洋.PLC变频调速系统应用研究[J].工
—业仪表与自动化装置,2001(6):27-30.
收稿日期:2007-08-21
工艺与材料 EMC
A2008,35(4)
钕铁硼永磁材料热性能的分析
姚丙雷, 林 岩, 刘秀芹
1
2
1
(1.上海电器科学研究所(集团)有限公司电机分所,上海 200063;
2.LG电子(天津)电器有限公司研发中心,天津 300402)
摘 要:介绍了钕铁硼永磁材料的热稳定性。详细说明了钕铁硼永磁材料的测试方法,并对测试结果进行了分析。得出钕铁硼永磁材料在100℃及120℃时,退磁曲线的拐点位置分布情况以及永磁材料内禀矫顽力温度系数随温度变化的趋势,对永磁电机的设计具有参考价值。
关键词:钕铁硼永磁材料;高温磁性能;不可逆退磁
中图分类号:TM273∶TM201.4+2 文献标识码:A 文章编号:1673-6540(2008)04-0052-04
AnalysisofThermalPerformanceof
NdFeBPermanentMagnet
YAOBing-lei, LINYan, LIUXiu-qin
(1.ShanghaiElectricalApparatusResearchInstitute(Group)Co.,Ltd,Shanghai200063,China;
2.LG(Tianjin)ElectronicsApparatusCo.,Ltd,Tianjin300402,China)
Abstract:ThermalstabilityofNd-Fe-Bpermanentmagnetwasintroduced.ThetestmethodsofNd-Fe-Bperma-nentmagnetmaterialweredescribedindetailandtestresultswereanalyzed.ThedistributionofkneepointofNd-Fe-Bpermanentmagnetmaterialdemagnetizationcurveon100℃and120℃wasobtainedandthechangingtrendofin-herentcoerciveforcetemperaturecoefficientastemperaturechangingwasachieved.Thesereferencevaluesarehelpfultothepermanentmagnetmachinedesign.
Keywords:Nd-Fe-Bpermanentmagnetmaterial;magnetismperformanceinthermalstate;irreversibledemagnetization
1
2
1
0 引 言
随着永磁材料性能的不断提高,稀土永磁材料被广泛应用于电机领域
[1]
决定永磁电机能否可靠运行。当永磁电机运行时,有时会出现转子温升较高而使钕铁硼永磁材
料发生不可逆退磁的现象,从而引起永磁同步电动机效率和功率因数等性能指标变差。因此,有必要对钕铁硼材料的热稳定性进行研究、分析。
本文采用测量的分析方法,在大量的试验结果基础上进行数据分析,目的在于获得高温时同一永磁材料各部分性能和同一批永磁材料性能的测试数据,并对测试数据以及由此对永磁材料退磁情况的影响进行分析,为永磁电机的设计提供参考依据。
。采用高性能的稀
土永磁材料制作电机,在结构上不需要专门的电励磁系统,有利于提高电机的工作效率和功率因数,并能节省能源,减小电机的体积
[2]
。
目前国内在设计永磁电机时,一般假设永磁体各部分磁场分布均匀,且同一批次永磁体的磁性能相同,同时需要利用永磁体供应商提供的常温下的永磁体退磁曲线、内禀矫顽力、剩磁密度和最大磁能积等磁性能参数值。但是,由于工厂制造工艺和水平等原因,在高温下同一块永磁体各部分实际的性能存在差异,且同一批永磁体的性能也是不同的。而永磁电机正常运行时的性能与永磁体的磁性能密切相关,永磁体的磁性能直接
1 永磁材料的热稳定性
温度稳定性是指永磁体由所处环境温度的改变而引起磁性能变化的程度。磁稳定性是表示在施加外加磁场条件下,引起磁性能发生变化的程
2008,35(4)
[3]
工艺与材料 EMCA
度。永磁电机在运行时,由于电机温升的影响100℃时有一块切割标准样拐点位置达到24.4%,120℃时拐点位置高达49.7%。拐点位
置的差异容易发生永磁体局部失磁。
表1 N35-SH-A在100℃和120℃时的温度特性
α(B)/r(%/K)100℃
N35-SH-A-1N35-SH-A-2N35-SH-A-5
-0.09-0.1-0.08
120℃-0.09-0.1-0.11
α(H)/cJ(%/K)100℃-0.73-0.70-0.70
拐点处磁密与剩磁密度比值/(%)
以及电枢绕组的存在,永磁电机内永磁材料的性
能由温度和磁场两个因素共同影响。因此,本文的分析中将温度稳定性和磁稳定性统称为热稳定性。即对永磁材料加热至不同的温度分别施加反向外加磁场时,永磁体所表现的磁性能的变化,也就是不同温度下永磁材料的剩磁密度的温度系数和内禀矫顽力的温度系数的大小以及退磁曲线上拐点位置的变化。
试样
120℃100℃120℃-0.67-0.67-0.66
19.618.623.4
26.532.730.4
2 测试方法及结果分析
测量仪器为稀土永磁检测专用设备,其测试
原理是通过测得的J-H曲线来计算B-H曲线。首先对试样进行充磁,然后对检测设备进行初始化,最后将试样装好,并将其加热到预定的温度,起动测试程序进行测试,测得的曲线由计算机记录并打印。
测试工作主要从以下几个方面展开:(1)对一块永磁体不同位置切割标准样进行测试;
(2)对同一厂家、同一牌号的不同永磁体(主要是 10×10mm标准样)进行测试;(3)对同为SH牌号的大块样品进行检测。2.1 同一样品不同位置热稳定性的检测分析
将2块永磁体N35-SH-A和N35-SH-B试样分别切割成8块标准试样(切割方法如图1所示),从所切割的8块标准试样中选取有代表性(不同位置)的3块(1、2、5)分别进行测试,并计算其在100℃和120℃时的温度特性,分析结果见表1、2所示
。
表2 N35-SH-B在100℃和120℃时的温度特性
α(B)/r(%/K)100℃
N35-SH-B-1N35-SH-B-2N35-SH-B-5
-0.06-0.07-0.06
120℃-0.09-0.08-0.08
α(H)/cJ(%/K)100℃-0.69-0.63-0.68
拐点处磁密与剩磁密度比值/(%)
试样
120℃100℃120℃-0.67
24.4
49.724.325.2
-0.61不弯-0.64不弯
2.2 同厂家同牌号不同永磁样品热稳定性的测
试分析
分别对来自2个厂家的、不同永磁标准样品的、高温下的性能进行测试分析,结果见表3、4所示。
表3 a厂4块样品高温下测试结果
α(B)/r(%/K)100℃
a1a2a3a4
-0.11-0.12-0.11-0.11
120℃-0.11-0.11-0.11-0.11
α(H)/cJ(%/K)100℃-0.67-0.67-0.67-0.65
120℃-0.63-0.63-0.62-0.62
拐点处磁密与剩磁密度比值/(%)100℃不弯不弯不弯不弯
120℃14.516.314.512.2
试样
图1 切割方法示意图
从分析数据可看出:在高温情况下,同一磁块切割下来的不同标准试样,在标志永磁体热稳定性好坏的2个重要参数(温度系数和拐点位置)上有着很大的差别。如:N35-SH-B的3个样品,
从对钕铁硼标准样品高温下的α(B、r)α(H)和拐点位置这三个参数的检测结果可得出cJ以下结论。
(1) α(B) 的变化范围不大,基本都在r0.11±0.01%/K范围内。2个厂家的8块样品
—
工艺与材料 EMC
A
表4 b厂4块样品高温下测试结果
试样
α(B)/r(%/K)100℃
b1b2b3b4
-0.10-0.10-0.10-0.08
120℃-0.11-0.10-0.11-0.10
α(H)/cJ(%/K)100℃-0.55-0.56-0.56-0.59
120℃-0.56-0.57-0.55-0.57
拐点处磁密
与剩磁密度比值/(%)100℃不弯不弯不弯10.1
120℃2.986.004.7023.30
2008,35(4)
表5 各厂家SH牌号NdFeB永磁体在100℃
和120℃时拐点位置情况
永磁体生产厂家
100℃时的拐点(B)/%k/Br
120℃时的拐点120℃(B/B)/%kr
2块永磁体拐点出现在
8块永磁体没有出现拐点,5块永磁体
甲厂(13)出现拐点,其拐点分
别为(15.87、22.00、24.78、6.09、5.08)
10%以下,11块永磁体出现拐点,其拐点分别为(18.94、17.33、19.43、20.34、21.46、13.05、26.70、34.60、37.00、21.10、21.90)
8块永磁体没有出现拐点,4块永磁体
乙厂(12)出现拐点,其拐点分
别为(11.71、8.29、7.89、7.19)1块永磁体没有出现拐点,4块永磁体
丙厂(5)出现拐点,其拐点分
别为(17.40、9.22、7.45、4.42)
3块永磁体拐点出现在10%以下,9块永磁体出现拐点,其拐点分别(16.30、15.25、
17.95、
10.73、
10.86、24.60、21.40、11.00、24.00)
1块永磁体拐点出现在10%以下,4块永磁体出现拐点,其拐点分别(34.1、23.2、27.3、22.3
)
中, α(B 最大的也仅为0.12%/K,最小为r)0.08%/K。
(2)2个厂家的α(H检测结果各有不同。cJ)a厂和b厂样品的α(H)较为稳定,但a厂的数cJ值均普遍偏高,120℃时, α(H均高于cj) 0.6%/K,变动范围在0.61～0.63%/K之间;b厂则在0.55～0.57%/K范围以内。
(3)从检测结果来看,2个厂家的永磁材料在高温下的拐点位置普遍偏高。a厂4块样品虽然在100℃时都为直线,但在120℃时拐点位置均高于10%,且最大值为16.3%;b厂4块样品中有3块满足100℃时为直线、120℃拐点位置低于10%,但是b4号样品与前3块偏差较大,100℃时拐点位置已高于10%,120℃时更达到了23.3%,高出120℃要求的1倍多。
2.3 大块样品的热稳定性测试结果与分析
本次试验是以内置式永磁电机中的大块永磁体为试样进行测试的。表5是永磁样品在100℃和120℃时,退磁曲线上拐点处的Bk值占Br的百分比的统计情况。
从对不同厂家的钕铁硼永磁材料的试验结果可以看出:
(1)对于每块试样来说,随着试样温度的升高,其B温度系数 α(B) 呈现上升趋势,H温rrcJ度系数 α(H 呈现下降趋势,图2是两种样品cJ)的 α(B) 和 α(H) 变化趋势图;rcJ
(2)从已检测的试样来看,大部分试样都未满足电机行业的要求(如SH牌号永磁体要求100°C时退磁曲线是直线;120°C时,退磁曲线拐点低于10%B。其中100°C时为直线的仅占r)56.7%,120°C拐点低于10%Br的仅占20%。100℃和120℃时退磁曲线拐点位置的分布情况如图3所示;
图2 两种样品 α(B) 和 α(H) 随温度变化趋势rcJ
(3)所测样品的内禀矫顽力的温度系数α(H)普遍偏高。从图4可以看出, α(H
2008,35(4)工艺与材料 EMC
A
63.2%; α(H >0.7%/K的占24.5%,最大cJ)为0.769%/K。这显然对防失磁是不利的。
3 结 语
本文在对大量钕铁硼永磁材料检测的基础上,对检测结果进行了分析。测试结果表明:目前仅根据常温下的矩形度H、H值以及BK/HcJcJr值并不能判断产品的热稳定性是否合格。如SH牌
号的NdFeB永磁体在矩形度大于或等于0.9、HA(20kOe)的情况下,并不能保证120cJ≥1600k
°C时拐点小于10%。究其原因,钕铁硼永磁材料生产企业没有把内禀矫顽力的温度系数α(H)作cJ为验收标准。总之,本文的分析对永磁电机的设计具有指导意义,有较好的工程利用价值,也为正
图3 100°C和120°C
时拐点位置分布图
确选用钕铁硼永磁材料以防止永磁电机高温退磁提供了依据,具有实际意义。
【参考文献】
[1] 彭龙,张怀武,宋远强,等.永磁材料在电机中的应
用现况分析[J].新材料产业,2006(11):30-33.[2] 王秀和.永磁电机[M].北京:中国电力出版社,
2007.
[3] 唐任远.现代永磁电机理论与设计[M].北京:机
械工业出版社,1997.
收稿日期:2008-01-10
图4 所测样品的 α(H) 分布情况cJ
(上接第51页
)
图4 系统网络总结构框图
[3] 上海交通大学罗克韦尔自动化实验室.罗克韦尔自
动化网络控制平台用户手册[G].2000.
[4] 杜清珍,李洋.PLC变频调速系统应用研究[J].工
—业仪表与自动化装置,2001(6):27-30.
收稿日期:2007-08-21