“微波多普勒测位移”用于电磁炮的弹丸位置检测
三水合番2019/08/19研究了个寂寞 IP:英国
中文摘要
想要控制多级电磁炮,一般需要知道弹丸在什么时候运动到了什么地方,也就是弹丸的“位置-时间关系”。多年以来,业余爱好者的一直都在用光电传感器来检测弹丸位置,也有一些用比如“金属探测器”的做法。不过这些做法的问题是,它们都只能检测弹丸是否运动到了某“一个”点。就是说“位置-时间关系”曲线上的每一个点都需要一个单独的传感器才能得到。想要以足够高的精度掌握整个发射过程中的“位置-时间关系”,就需要一大堆的传感器——既麻烦又贵。所以本文介绍一种“只需要一个传感器,就能在整个发射过程中一直检测弹丸位置的方案”,即“微波多普勒测位移”。
关键词
电磁炮微波多普勒

多普勒效应是波源和观察者相对运动时,观察者接收到的波的频率与波源发出的频率并不相同的现象。对于静止的波源s和运动的反射体r,在波源s处接收到的反射波的频率f满足

$$ f=\frac{v_{0}-v_{r}}{ v_{0}+v_{r}}f_{0} \tag{1}$$式中,v0为波的传播速度;vr为反射体的运动速度(当反射体远离波源时,vr为正);f0为波源发射的波的频率。对于本贴讨论的“微波”来说,它的波速为光速,即v0=c。同时,电磁炮的弹丸运动速度v远远小于光速,所以有

$$|\Delta f|=\frac{2v}{c } f_{0} \tag{2}$$其中|Δf|是波源和反射波的频率差的绝对值。通过测量这个频率差,就可以求出弹丸运动速度。

$$ v = \frac{ c |\Delta f|}{2 f_{0} }=\frac{\lambda }{2 } |\Delta f| \tag{3}$$其中,λ是微波的波长。不过我们一般需要的是弹丸的位置,而不是速度,所以还需要把这个式子再处理一下。速度是“位置对时间的导数”,而频率是“相位对时间的导数”,所以对它两端积一下分,就有

$$ \Delta x = \frac{\lambda}{4\pi} \Delta\varphi \tag{4}$$$$ 位移=波长\times\frac{相位差}{4\pi} \tag{5}$$所以只需要检测波源和反射波之间的相位差,就能知道弹丸的位移。检测相位差的方法被称为“鉴相”。对于“弹丸位置检测”这个应用,我们既不知道反射波的幅度,也不知道它的相位。对此,一般鉴相方法的是,用波源和一个“与波源带90°相差的信号”,分别与反射波相乘。测量两个相乘结果的直流分量,然后代入一个简单的公式里,就能算出波源和反射波的相位差(以及反射波的幅度)。这种方案的好处是,它输出的是一个和弹丸位移相关的模拟量,也就是说原理上可以得到任意高的分辨率。不过这种方案目前市面上没有成品,需要自行制作。现在有些比较方便的芯片可以用,比如这个

attachment icon Infineon-BGT24LTR11N16-DS-v01_03-EN.pdf 441.49KB PDF 166次下载 预览

做这种几十GHz的系统,对仪器要求非常的高,短期内似乎不会有人来做这个。所以这里介绍另一种方法——过零点计数法

这种方案可以用市面上常见的“微波多普勒模块”来做,如下图。目前这种模块在淘宝上大约是20-30元一个,工作频率为24GHz。

微波多普勒模块.png

图 1 24GHz微波多普勒模块

The Signal Path 上有关于这个模块的详细的测评。

XXXXXXXXXXXXXXXX/5vqSX40seqA

它的原理框图如下

微波多普勒雷达 函数版.png

图 2 微波多普勒模块原理框图以及各点信号波形

其中A为各点信号幅度,可以认为是未知量。发射天线发射的是一个正弦波;接收天线收到的是频率相同,但有额外相移的正弦波。整个模块输出的是相位差的函数,但是由于各点的幅度A未知,所以没法测量得到一个连续的“位置-时间关系”。

不过我们至少知道,当cos(φ)=0时,不论幅度A的值为多少,输出信号都会过一次零。所以,可以通过“数过零点的个数”来计算弹丸的位移,即前面提到的“过零点计数法”。结合图 2和式 4,可以知道弹丸每移动1/4个波长,模块的输出信号就会过一次零。不过需要注意的是,这里的波长不等于真空中的波长,而是要按“波导波长”来算,具体的我们稍后再说。


使用这种模块的时候,最简单的做法是:不加额外的结构,直接让模块的天线正对着炮管。如下图

IMG_20181129_090729.jpg

图 3 最简单的用法

这里的“炮管”是一根内径8.1mm的不锈钢管。“弹丸运动”是通过徒手拉动一根约5mm直径的不锈钢棒来模拟的。场地布置如下

IMG_20181129_090746 注释.jpg

图 4 场地布置

当不锈钢棒在炮管内运动时,模块能输出几十mV级别的信号。如下视频


VID_20181129_090630_Trim.mp4 点击下载

之后可以接一级简单的放大,送给比较器检测过零点,然后传给单片机统计脉冲个数,进而确定弹丸位移。


对于这种无线系统来说,几十mV的信号可以说是非常强了(正常情况下输出信号应当是弱到可以淹没在示波器底噪里的)。这是因为炮管起到了“波导”的作用。不锈钢制的炮管可以当作是一个“圆波导”,能以很低的损耗传递电磁波,所以模块接收到的反射波幅度也会很强。类似于对着井口喊话,听到的回音要比“对着一棵树喊话”强得多,尽管井底的面积比树小得多。

圆波导电场图.gif

图 5 24GHz的电磁波在8.1mm内径的圆波导中传播

然而“波导”也会带来两个问题。

一个是这种圆波导只允许频率足够高的电磁波在其中传播,如果电磁波频率太低(也就是波长太长),电磁波在进入波导后就会非常迅速地衰减,比如在一两个自由空间波长的距离上衰减几十dB,这个现象称为“截止”,对应的电磁波频率和自由空间波长则称为“截止频率”和“截止波长”。截止波长和波导的尺寸相关,对于圆波导,它的主模(TE11模)截止波长为1.707D,D为圆管的内直径。8.1mm的圆波导,截止波长为1.707*8.1=13.8mm,对应21.7GHz的电磁波,因此内径8.1mm的金属圆管中,只能传播频率高于21.7GHz的电磁波。

另一个问题是,波导会改变电磁波的波长,计算位移时要考虑到这个影响。电磁波在波导中的波长被称为“波导波长”,波导波长的计算方法如下

$$ \lambda_{g}=\frac{\lambda}{\sqrt{1-\left( \frac{\lambda}{\lambda_{c}} \right)^{2}}} \tag{6}$$其中,λ为自由空间(无限大真空)波长,等于“真空中光速/电磁波频率”;λc为截止波长,圆波导主模截止波长为1.707D(D为圆管的内直径)。24GHz电磁波的自由空间波长=3*10^8/(24*10^9)=12.5mm,前文已知8.1mm圆波导的主模截止波长为13.8mm。因此图 4中的波导波长为

$$ \lambda_{g}=\frac{12.5}{\sqrt{1-\left( \frac{12.5}{13.8} \right)^{2}}}=33.1\; \rm mm \tag{7}$$这个计算结果与图 5中的仿真结果(约31mm)基本相符。因此使用24GHz模块和8.1mm内径炮管,弹丸每移动8.3mm,输出信号就会过一次零。


除了波导带来的问题以外,这个方案本身也会带来一些其他问题。比如微波多普勒模块要装在炮管的尾部,弹丸也要从这个地方进入炮管,所以供弹结构会影响到模块的工作。连发炮上使用这个方案的话,可能需要特别设计的机械和微波结构来消除这个影响。

另外根据The Signal Path上的测试,图 1中的这种模块的输出频率会随很多因素改变,如图 6。

Big DopplerModule~1.gif

图 6 模块输出频率随供电电压的变化

不过好在频率变化不大,只有万分之几,影响基本可以忽略。

[修改于 4年7个月前 - 2020/05/02 20:32:39]

+1  科创币    忆昔长别    2020/02/16 涨见识了,支持一下
+1  学术分    虎哥    2019/12/06 关于小口径炮管波导测距测速的首次实验研究。
来自:物理高能技术 / 电磁炮电子信息 / 无线电严肃内容:专著/论述
10
 
6
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
zx-16533
5年4个月前 IP:广东
862633

对于速度为恒为正的情况下可以直接将模块输出的信号进行希尔伯特变换后作为Q路信号与原信号一起反正切即可完成鉴相,之后将相位角在时间上进行解缠绕可以得到位移。

安装的话大概可以弄个90°的弯波导和炮管衔接。

引用
评论
8
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
三水合番作者
5年4个月前 修改于 5年4个月前 IP:英国
862634
引用zx-16533发表于1楼的内容
对于速度为恒为正的情况下可以直接将模块输出的信号进行希尔伯特变换后作为Q路信号与原信号一起反正切即可...

这种模块输出的是已经经过混频和低通滤波的kHz级的信号😂

哦……刚理解错了……确实应该可以

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
zx-16533
5年4个月前 修改于 5年4个月前 IP:广东
862635
引用三水合番发表于2楼的内容
这种模块输出的是已经经过混频和低通滤波的kHz级的信号😂哦……刚理解错了……确实应该可以

我说的就是这个输出信号啊

这个输出应该还有个高通或者电容耦合,不然肯定会有很大的直流分量。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
三水合番作者
5年4个月前 IP:英国
862636
引用zx-16533发表于3楼的内容
我说的就是这个输出信号啊这个输出应该还有个高通或者电容耦合,不然肯定会有很大的直流分量。

额……似乎还是不行……

希尔伯特变换是“非因果”的,没法实时进行。这样就只能把波形录下来做事后分析,没法实现实时控制。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
zx-16533
5年4个月前 修改于 5年4个月前 IP:广东
862638
引用三水合番发表于4楼的内容
额……似乎还是不行……希尔伯特变换是“非因果”的,没法实时进行。这样就只能把波形录下来做事后分析,没...

数学上理想的希尔伯特变换现实上是没办法实现的,现实中一般的做法是用一个滤波器来代替,或者说类似于录下之前的一小段来做分析。这样的话就会引入一定的延迟,如果要实时控制的话似乎需要在延迟和频率响应之间权衡一下……

记得有款类似的模块是自带i/q两路正交信号输出的,那个价钱比这个贵不少。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
三水合番作者
5年4个月前 IP:欧洲
862687
引用zx-16533发表于5楼的内容
数学上理想的希尔伯特变换现实上是没办法实现的,现实中一般的做法是用一个滤波器来代替,或者说类似于录下...

刚想了一下,用一般意义上的滤波器做(线性时不变的那种),似乎也不太合适。因为这种滤波器带宽和延时是固定的,而电磁炮发射时模块的输出信号带宽很宽。想要在低频把希尔伯特变换做准的话,高频的时候延时就会大的不可接受。做“时变”的滤波器,实时改频率特性的话,似乎又有些过于麻烦了。

另外对电磁炮来说,低延迟似乎比高精度还要更有用一些。用“根据过零点求的位移”,加上简单的线性外推,或者二次抛物线外推,应该会更合适一些

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
虎哥
5年0个月前 IP:四川
867051

做一个尺寸比弹丸大得多的波导分路器,使金属部分的空间不侵占炮管径向。用塑料件来填塞多于空隙,做成介质波导,塑料上可以打孔通过弹丸。与炮管衔接处可能需要一些特殊的型面进行模式转换,加工会难一点。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
极速外卖
5年0个月前 IP:广东
867244
引用zx-16533发表于3楼的内容
我说的就是这个输出信号啊这个输出应该还有个高通或者电容耦合,不然肯定会有很大的直流分量。

是直接耦合输出的,示波器可以看到1V以上直流偏移




引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
Ball_Lightning
2年7个月前 IP:河北
903751

请问炮管必须用塑料的话怎么用这个做分析呢?

我对这种用多普勒测量弹丸运动设想很感兴趣。因为它不仅能够得到出口速度,还能得到弹丸在整个运动过程中的运动曲线,可以看出弹丸在整个过程中有没有回拉减速,各级的实际lc波长和弹丸峰值受力等参数,方便进行全面系统的分析和优化。

但是,微波多普勒很多问题,

1.在猛烈放电时不锈钢会有损耗。2.一米长的不锈钢管+打孔非常昂贵,且无法调整触发位置。3.我的触发方式是高压间隙触发,不能用不锈钢。4.微波多普勒装置里驱动线圈很近时如果触发失去时序,就会被极强的emp效应爆掉。。。。。。。。。。 

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
三水合番作者
2年7个月前 IP:四川
903752
引用Ball_Lightning发表于9楼的内容
请问炮管必须用塑料的话怎么用这个做分析呢?我对这种用多普勒测量弹丸运动设想很感兴趣。因为它不仅能够得...

炮管用塑料的话,应该是没啥简单的办法能用这个方案的,建议放弃。

1. 不锈钢的电阻率极高,是铜的四十多倍,炮管涡流的影响应该不大。可以仿真一下看看用不锈钢的时候,有炮管和没炮管的性能差别。4. 多普勒模块可以放在离线圈很远的地方,只要有金属炮管作为波导,对测试效果应该没啥影响。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

想参与大家的讨论?现在就 登录 或者 注册

文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}