CN106556457A - 用于小型飞行器的数据记录系统 - Google Patents
用于小型飞行器的数据记录系统 Download PDFInfo
- Publication number
- CN106556457A CN106556457A CN201610997294.9A CN201610997294A CN106556457A CN 106556457 A CN106556457 A CN 106556457A CN 201610997294 A CN201610997294 A CN 201610997294A CN 106556457 A CN106556457 A CN 106556457A
- Authority
- CN
- China
- Prior art keywords
- recording system
- digital data
- data recording
- record
- aircraft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000001133 acceleration Effects 0.000 claims abstract description 8
- 238000012544 monitoring process Methods 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 2
- 238000004891 communication Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 5
- 230000005540 biological transmission Effects 0.000 abstract description 2
- 238000013500 data storage Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013523 data management Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
- G01H11/06—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
- G01H11/08—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means using piezoelectric devices
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Recording Measured Values (AREA)
Abstract
用于小型飞行器的数据记录系统,涉及数据记录领域,解决现有用于飞行器的数据记录系统存在体积大、采集信息种类单一以及通讯接口的通用差等问题,包括抗混叠滤波器、A/D采集模块、嵌入式处理器、存储器、振动传感器、过载传感器、霍尔传感器、GPS模块和电源模块;本发明通过振动传感器记录飞行器重要部位的加速度信息,记录过载传感器采集的发动机的三轴过载,通过RS232接口记录GPS模块的定位状态、位置、速度、姿态信息,数据记录系统通过CAN接口记录飞行器上通讯网络传输的控制命令与反馈信息,存储器的采集数据存储机制实现了数据记录系统的灵活多次连续工作。本发明可安装于小型飞行器上,为飞行状态监测与飞行器的设计改进提供了数据支撑。
Description
技术领域
本发明涉及数据记录领域,具体涉及一种用于小型飞行器的数据记录系统,用于监测小型飞行器(如无人机及导弹)上电压、电流、姿态、速度、控制指令及环境数据等信息的多用途数据记录系统。
背景技术
目前,用于记录飞行器在飞行试验期间的各项过程参数信息主要采用数据记录系统。但现有的数据记录系统普遍存在采集信息种类单一、体积大、重量沉、通讯接口简单及不通用等缺点。而随着科技的进步,飞行器正在向着更轻、更快、数据交换量更大的方向发展,不仅需要记录常规传感器的模拟信号、数据总线传输的指令和数据信息,还要能够记录导弹在设计优化期间所需要的飞行过程中加速度、过载、供电电压、供电电流及GPS状态等多方面信息,因此也就要求数据记录系统应具有小型化、轻量化、多接口、多功能和多用途的特点。
小型飞行器具有需求量大、应用灵活、体积小、重量轻、成本低等特点,开发适用于此类飞行器的数据记录系统应在保证抗高落地过载的同时又需要严格控制其空间尺寸和重量,因此传统的数据记录系统很难满足这种新型化的数据记录系统需要。
发明内容
本发明为解决现有用于飞行器的数据记录系统存在体积大、采集信息种类单一以及通讯接口的通用差等问题,提供一种用于小型飞行器的数据记录系统。
用于小型飞行器的数据记录系统,包括抗混叠滤波器、A/D采集模块、嵌入式处理器、存储器、振动传感器、过载传感器、霍尔传感器、GPS模块和电源模块;所述抗混叠滤波器用于对振动传感器采集的飞行器加速度变化信息进行处理,嵌入式处理器控制所述抗混叠滤波器输出的信号经A/D采集模块后存入存储器;过载传感器用于测量飞行器发射时发动机产生的三轴过载值,所述嵌入式处理器控制所述过载传感器输出的三轴过载值经A/D采集模块输出后存入存储器;所述GPS模块输出飞行器试验过程中的速度、位置及姿态信息,通过RS232接口传送至嵌入式处理器,并存入存储器;电源模块为数据记录系统供电,通过电源管理模块与内部电源连接;数据回收与解析模块,用于对存储模块输出的数据进行处理与分析。
本发明的有益效果:本发明所述的数据记录系统能够全面可靠地记录飞行器上的电源供电变化、数据通信、环境参数变化等信息,同时为达到通用性还要具备多种常见的输入和输出接口。
本发明所述的数据记录系统能够记录小型飞行器在飞行过程中的加速度、三轴过载、供电电压、供电电流、系统通讯信息及GPS模块状态等多方面的信息,为飞行器设计的优化及飞行试验过程中的问题分析提供依据;
本发明所述的数据记录系统可以实现连续多次工作,能够做到连续记录30次而记录数据不相互覆盖,这样可以在地面测试时不需要频繁地进行数据回收;
本发明所述的数据记录系统采用多层缓存与灌封工艺,可承受在导弹落地时对地面产生的不大于8000g的落地冲击,在保证抗高落地过载的同时能够满足小型飞行器对于数据记录系统的体积与重量的限制。
附图说明
图1为本发明所述的用于小型飞行器的数据记录系统的框图;
图2为本发明所述的用于小型飞行器的数据记录系统的原理框图;
图3为本发明所述的用于小型飞行器的数据记录系统中传感器安装位置示意图;
图4为采用本发明所述的用于小型飞行器的数据记录系统抗混叠滤波效果图。
具体实施方式
具体实施方式一、结合图1至图4说明本实施方式,用于小型飞行器的数据记录系统,该记录系统为具有多功能、多接口及多用途的数据记录系统,本实施方式以一种小型导弹上的数据记录系统作为具体实施例。具体包括抗混叠滤波器、A/D采集模块、嵌入式处理器、存储器、振动传感器、过载传感器、霍尔传感器、GPS模块和弹上电源模块;抗混叠滤波器用来降低振动传感器采集数据的频域混叠,A/D采集模块对采集的模拟信号进行数字量化,嵌入式处理器用于采集和存储控制,存储器用于存储各通道的采集信息,振动传感器用于采集三轴振动信息,过载传感器用于采集飞行器的三轴过载,霍尔传感器用于采集热电池电流信号,GPS模块用于获取飞行器的速度、位置及姿态信息,弹上电源模块用于与数据记录系统内部电池供电和管理。
结合图2,振动传感器1、振动传感器2和振动传感器3分别安装在导弹需要监测加速度变化的各个部位。其具体的安装位置结合图3,振动传感器1安装在导引头处,选用东华测试1A314E三向压电加速度传感器,量程为±50g,用于测量飞行试验过程中导引头处的三轴振动;振动传感器2安装在导航制导控制系统处,选用东华测试1A313E三向压电加速度传感器,量程为±125g,用于测量飞行试验过程中GPS模块和弹载计算机处的三轴振动;振动传感器3安装在舵机系统处,选用东华测试1A313E三向压电加速度传感器,量程为±125g,用于测量飞行试验过程中舵机系统的三轴振动。
振动传感器1、振动传感器2和振动传感器3采集的数据分别通过抗混叠滤波器处理以降低振动数据的频域交叉混叠,抗混叠要求的低通滤波效果,结合图3,在通频带内幅频特性误差在±3dB以内,过渡带的斜率要尽量陡峭,斜率至少为-45dB/oct,采用巴特沃斯滤波器进行低通滤波,一阶巴特沃斯滤波器可以达到-6dB/oct的衰减,使用8阶巴特沃斯滤波器进行低通滤波,衰减斜率为-48dB/oct,能够满足过渡带衰减斜率的要求。
本实施方式中所述的嵌入式处理器采用数字信号处理器DSP,型号为TMS320F28335。抗混叠滤波后的模拟输出信号在嵌入式处理器的控制下经过A/D采集后存入存储器中,所述振动传感器1、振动传感器2和振动传感器3的采样频率均为10kHz。
本实施方式中所述的过载传感器用于测量导弹发射时发动机所产生的三轴过载值,其主要以模拟量形式输出。根据飞行试验测试需求,过载传感器安装在发动机附近,结合图2,过载传感器的量程为±20g,轴向、横向和垂向的模拟过载值在嵌入式处理器的控制下经过A/D采集后存入存储器中,过载传感器输出数据的采样频率为1kHz。
电源管理模块负责数据记录系统的供电,飞行时通过机载或弹载电源供电,当其故障时切换到内部电池供电以确保数据记录系统的正常工作;在地面测试时,也可通过USB接口供电,直接通过USB接口控制数据记录系统启动采集进行地面测试。
数据回收与解析软件用于记录数据的回收、显示与擦除,可以进行简单的数据处理与分析,也可以使用数据回收与解析软件将所存储的各项信息存储成文本文件,通过MatLab对记录数据进行更深层次的数据处理分析。
导弹试验过程中的供电电源的电压信息通过数据记录系统的模拟输入通道采集并记录,在飞行试验后可以依此判断电气系统的工作状态。热电池在放电过程中电压变化缓慢,采用1kHz的采样频率即可;供电电流信息通过数据记录系统内部的霍尔传感器进行测量并记录,霍尔传感器的精度在2%左右,在飞行试验前,需要对电流传感器进行全量程标定,为了捕捉舵机动作时的脉冲电流,电流采集通道需要较高的采样频率,其采样频率设置为50kHz。
所述供电电源采用热电池供电,通过在飞行试验过程中供电电压和电流的测量及记录,试验后对电压放电曲线、平均电流、脉冲电流的峰值及持续时间的分析及处理可计算出热电池实际用电量,为热电池的后续优化设计提供重要依据。
本实施方式中所述的GPS模块输出导弹飞行试验过程中的速度、位置及姿态等动态信息,由此可推算飞行器航迹,其通常以RS232总线接口形式输出,该数据记录系统通过RS232总线接口可以完成对GPS模块输出信息的接收和记录,RS232总线通讯采用115200bps的总线波特率,8个数据位,1个停止位,无校验位,帧频为10Hz。
目前弹上通讯网络主要采用CAN总线接口形式,因此该数据记录系统配置有两路隔离CAN总线接口,可以监测并记录总线上发送的控制指令及各受控设备的反馈状态,CAN总线的波特率为1Mbps,采用标准帧数据格式,CAN标示符为11位。
本实施方式所述的数据记录系统为满足通用性要求,预留多种数字接口,其中包括SPI接口、RS232接口、RS422接口、RS485接口和CAN总线接口,可适应多种应用场合的使用。
本实施方式所述的数据记录系统可以实现连续多次工作,存储容量高达8GByte,能够做到连续记录30次而记录数据不相互覆盖,这样可以在地面测试时不需要频繁地进行数据回收;采用多层缓存与灌封工艺,可承受在导弹落地时对地面产生的不大于8000g的落地冲击。在保证抗高落地过载的同时能够满足小型飞行器对于数据记录系统的体积与重量的限制。
Claims (5)
1.用于小型飞行器的数据记录系统,包括抗混叠滤波器、A/D采集模块、嵌入式处理器、存储器、振动传感器、过载传感器、霍尔传感器、GPS模块、弹上电源模块和供电电源;其特征是;
所述抗混叠滤波器用于对振动传感器采集的飞行器加速度变化信息进行处理,嵌入式处理器控制所述抗混叠滤波器输出的信号经A/D采集模块后存入存储器;
过载传感器用于测量飞行器发射时发动机产生的三轴过载值,所述嵌入式处理器控制所述过载传感器输出的三轴过载值经A/D采集模块输出后存入存储器;
所述GPS模块输出飞行器试验过程中的速度、位置及姿态信息,通过RS232接口传送至嵌入式处理器,并存入存储器;
所述弹上电源模块为数据记录系统供电,通过电源管理模块与内部电池连接;
飞行器的供电电源为热电池,电流信息经霍尔传感器测量采集后输出至A/D采集模块记录,电压信息经A/D采集模块采集并记录;通过数据回收与解析模块,对电压放电曲线、平均电流以及脉冲电流的峰值及持续时间的分析及处理可计算热电池实际用电量;
数据回收与解析模块,用于对存储模块输出的数据进行处理与分析。
2.根据权利要求1所述的用于小型飞行器的数据记录系统,其特征在于,所述振动传感器分别安装于飞行器的导引头位置、导航制导与控制系统以及舵机系统上。
3.根据权利要求1所述的用于小型飞行器的数据记录系统,其特征在于,还包括USB接口,电源模块通过USB接口供电,控制数据记录系统启动采集并进行测试。
4.根据权利要求1所述的用于小型飞行器的数据记录系统,其特征在于,所述嵌入式处理器还设置有两路隔离CAN总线接口,用于监测并记录总线上发送的控制指令以及各受控设备的反馈状态。
5.根据权利要求1所述的用于小型飞行器的数据记录系统,其特征在于,所述嵌入式处理器预留多种数字接口,具体包括SPI接口、RS232接口、RS422接口、RS485接口和CAN总线接口,用于相应设备的连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610997294.9A CN106556457A (zh) | 2016-11-11 | 2016-11-11 | 用于小型飞行器的数据记录系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610997294.9A CN106556457A (zh) | 2016-11-11 | 2016-11-11 | 用于小型飞行器的数据记录系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106556457A true CN106556457A (zh) | 2017-04-05 |
Family
ID=58443745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610997294.9A Pending CN106556457A (zh) | 2016-11-11 | 2016-11-11 | 用于小型飞行器的数据记录系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106556457A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108613737A (zh) * | 2018-05-14 | 2018-10-02 | 南京理工大学 | 基于小波包与stft的飞行器多频振动信号的辨识方法 |
CN108845922A (zh) * | 2018-06-01 | 2018-11-20 | 深圳市道通智能航空技术有限公司 | 无人机电池信息记录方法、装置及存储介质 |
CN110430254A (zh) * | 2019-07-31 | 2019-11-08 | 西安爱生无人机技术有限公司 | 一种无人机遥控器解析信号方法、信号处理方法及系统 |
CN110758749A (zh) * | 2019-10-30 | 2020-02-07 | 中国航空救生研究所 | 一种基于三轴过载矢量和的过载启动方法 |
CN111862686A (zh) * | 2020-06-04 | 2020-10-30 | 中国人民解放军国防科技大学 | 一种飞行器运动状态测量与数据处理系统 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040220710A1 (en) * | 2003-04-30 | 2004-11-04 | General Electric Company | Vibration measurement and recording system and method |
CN1670487A (zh) * | 2004-03-15 | 2005-09-21 | 清华大学 | 用于微型飞行器的数据记录仪 |
CN102565717A (zh) * | 2012-01-18 | 2012-07-11 | 哈尔滨建成集团有限公司 | 一种热电池测试系统 |
CN102737416A (zh) * | 2012-06-08 | 2012-10-17 | 清华大学 | 无人飞行器的数据记录仪 |
CN102963533A (zh) * | 2012-12-14 | 2013-03-13 | 中国航空工业集团公司上海航空测控技术研究所 | 一种监测直升机完好性与使用状态的系统及其监测方法 |
CN103644963A (zh) * | 2013-11-28 | 2014-03-19 | 陕西千山航空电子有限责任公司 | 一种振动信号采集处理模块及处理方法 |
CN203673482U (zh) * | 2013-12-26 | 2014-06-25 | 重庆航天工业公司 | 一种弹载数据记录仪 |
CN103927798A (zh) * | 2014-04-11 | 2014-07-16 | 贵州大学 | 一种适用于中小型无人机飞行数据的记录方法及装置 |
CN104102127A (zh) * | 2014-07-17 | 2014-10-15 | 厦门大学 | 一种机载气动参数辨识系统 |
CN104816833A (zh) * | 2015-04-20 | 2015-08-05 | 中国科学院长春光学精密机械与物理研究所 | 小型飞行器电源系统及其一体化设计方法 |
CN204833392U (zh) * | 2015-08-27 | 2015-12-02 | 成都翼比特科技有限责任公司 | 一种无人机飞行数据记录仪 |
CN205192528U (zh) * | 2015-11-25 | 2016-04-27 | 中国科学院沈阳自动化研究所 | 一种面向轻型直升机的飞行数据采集系统 |
CN105785877A (zh) * | 2016-04-12 | 2016-07-20 | 清华大学 | 一种小型动态数据记录系统 |
-
2016
- 2016-11-11 CN CN201610997294.9A patent/CN106556457A/zh active Pending
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040220710A1 (en) * | 2003-04-30 | 2004-11-04 | General Electric Company | Vibration measurement and recording system and method |
CN1670487A (zh) * | 2004-03-15 | 2005-09-21 | 清华大学 | 用于微型飞行器的数据记录仪 |
CN102565717A (zh) * | 2012-01-18 | 2012-07-11 | 哈尔滨建成集团有限公司 | 一种热电池测试系统 |
CN102737416A (zh) * | 2012-06-08 | 2012-10-17 | 清华大学 | 无人飞行器的数据记录仪 |
CN102963533A (zh) * | 2012-12-14 | 2013-03-13 | 中国航空工业集团公司上海航空测控技术研究所 | 一种监测直升机完好性与使用状态的系统及其监测方法 |
CN103644963A (zh) * | 2013-11-28 | 2014-03-19 | 陕西千山航空电子有限责任公司 | 一种振动信号采集处理模块及处理方法 |
CN203673482U (zh) * | 2013-12-26 | 2014-06-25 | 重庆航天工业公司 | 一种弹载数据记录仪 |
CN103927798A (zh) * | 2014-04-11 | 2014-07-16 | 贵州大学 | 一种适用于中小型无人机飞行数据的记录方法及装置 |
CN104102127A (zh) * | 2014-07-17 | 2014-10-15 | 厦门大学 | 一种机载气动参数辨识系统 |
CN104816833A (zh) * | 2015-04-20 | 2015-08-05 | 中国科学院长春光学精密机械与物理研究所 | 小型飞行器电源系统及其一体化设计方法 |
CN204833392U (zh) * | 2015-08-27 | 2015-12-02 | 成都翼比特科技有限责任公司 | 一种无人机飞行数据记录仪 |
CN205192528U (zh) * | 2015-11-25 | 2016-04-27 | 中国科学院沈阳自动化研究所 | 一种面向轻型直升机的飞行数据采集系统 |
CN105785877A (zh) * | 2016-04-12 | 2016-07-20 | 清华大学 | 一种小型动态数据记录系统 |
Non-Patent Citations (3)
Title |
---|
纪勇等: "小型无人机CAN总线信号记录仪", 《数据采集与处理》 * |
韩勇豪等: "一种固态飞行参数记录仪的研制", 《电子设计工程》 * |
高思远等: "无人机载数据记录仪软件的实时性分析与优化", 《吉林大学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108613737A (zh) * | 2018-05-14 | 2018-10-02 | 南京理工大学 | 基于小波包与stft的飞行器多频振动信号的辨识方法 |
CN108845922A (zh) * | 2018-06-01 | 2018-11-20 | 深圳市道通智能航空技术有限公司 | 无人机电池信息记录方法、装置及存储介质 |
CN110430254A (zh) * | 2019-07-31 | 2019-11-08 | 西安爱生无人机技术有限公司 | 一种无人机遥控器解析信号方法、信号处理方法及系统 |
CN110430254B (zh) * | 2019-07-31 | 2021-11-16 | 西安爱生无人机技术有限公司 | 一种无人机遥控器解析信号方法、信号处理方法及系统 |
CN110758749A (zh) * | 2019-10-30 | 2020-02-07 | 中国航空救生研究所 | 一种基于三轴过载矢量和的过载启动方法 |
CN111862686A (zh) * | 2020-06-04 | 2020-10-30 | 中国人民解放军国防科技大学 | 一种飞行器运动状态测量与数据处理系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106556457A (zh) | 用于小型飞行器的数据记录系统 | |
CN201551845U (zh) | 固定翼航模飞行辅助系统 | |
CN205633027U (zh) | 一种多轴旋翼井下机器人系统 | |
CN105980251B (zh) | 用于飞行器的基于路径的功率产生控制 | |
CN101793591A (zh) | 飞行器气动伺服弹性地面模拟试验系统 | |
CN105807779A (zh) | 一种无人机飞行控制系统及方法 | |
CN103279126B (zh) | 小型无人直升机开发测试平台 | |
CN102288346B (zh) | 一种小型化数字式的大规模传感器阵列冲击监测系统 | |
CN102508493A (zh) | 一种小型无人飞行器飞行控制方法 | |
CN1669874A (zh) | 一种用于飞行器的自动驾驶仪 | |
CN106697292A (zh) | 一种载重量大的农林植保无人直升机 | |
CN104765377A (zh) | 基于qnx的无人直升机飞行控制平台系统 | |
CN2681997Y (zh) | 一种具有双处理器的飞行器自动驾驶仪 | |
CN102914668A (zh) | 一种研究鸟类飞行动力学的方法和装置 | |
CN107678332B (zh) | 一种基于惯性导航的快响应火箭投放系统及投放方法 | |
Nugraha et al. | Development of Rocket Telemetry in Chamber Gas Pressure Monitoring with the MPXV7002DP Gas Pressure Sensor | |
CN202483567U (zh) | 一种存储式地应力测试装置 | |
Dantsker et al. | Flight data acquisition platform development, integration, and operation on small-to medium-sized unmanned aircraft | |
CN202083397U (zh) | 一种惯性测量装置火箭橇试验数据采集系统 | |
CN105974935A (zh) | 一种四旋翼农用遥控飞行器及其控制方法 | |
CN205192528U (zh) | 一种面向轻型直升机的飞行数据采集系统 | |
KR100850896B1 (ko) | 무인 헬기의 비행 데이터 획득 시스템 및 방법 | |
CN102425978B (zh) | 一种弹载计算机 | |
CN204788412U (zh) | 一种适用于小型无人直升机的惯性传感器组件 | |
CN1282863C (zh) | 用于微型飞行器的数据记录仪 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170405 |
|
WD01 | Invention patent application deemed withdrawn after publication |