发明名称 多普勒振镜正弦调制多光束激光外差二次谐波测量电致伸缩系数的方法
摘要 多普勒振镜正弦调制多光束激光外差二次谐波测量电致伸缩系数的方法,它涉及一种测量电致伸缩系数的方法。它为解决现有采用激光外差测量技术在测量电致伸缩系数存在的采集激光差频信号质量较低,信号处理的运算速度慢的问题而提出。第二平面反射镜的反射面与薄玻璃板相互平行,振镜开始简谐振动;同时,打开H0固体激光器,调节高压电源的输出电压信号,同时信号处理系统连续采集并处理光电探测器输出的电信号获得第二平面反射镜和薄玻璃板后表面之间的距离变化量,得到待测压电陶瓷管的电磁致伸缩系数:它具有采集的激光差频信号质量高和信号处理的运算速度快的优点。它可在相关激光测风雷达等领域广泛使用。
申请公布号 CN102253002A 申请公布日期 2011.11.23
申请号 CN201110145157.X 申请日期 2011.05.31
申请人 哈尔滨工业大学 发明人 王春晖;李彦超
分类号 G01N21/45(2006.01)I 主分类号 G01N21/45(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 1.多普勒振镜正弦调制多光束激光外差二次谐波测量电致伸缩系数的方法,它是基于多普勒振镜正弦调制多光束激光外差二次谐波测量电致伸缩系数实现的,所述系统由H<sub>0</sub>固体激光器(2)、四分之一波片(12)、振镜(13)、第一平面反射镜(3)、偏振分束镜PBS(11)、会聚透镜(10)、薄玻璃板(9)、第二平面反射镜(6)、待测压电陶瓷管(7)、二维调整架(8)、高压电源、光电探测器(4)和信号处理系统(5)组成;H<sub>0</sub>固体激光器(2)发出的线偏振光经第一平面反射镜(3)反射之后入射至偏振分束镜PBS(11),经该偏振分束镜PBS(11)反射后的光束经四分之一波片(12)透射后入射至振镜(13)的光接收面,经该振镜(13)反射的光束再次经四分之一波片(12)透射后发送至偏振分束镜PBS(11),经该偏振分束镜PBS(11)透射后的光束入射至薄玻璃板(9),经该薄玻璃板(9)透射之后的光束入射至第二平面反射镜(6),该光束在相互平行的薄玻璃板(9)后表面和第二平面反射镜(6)之间反复反射和透射出薄玻璃板多次,获得多束经薄玻璃板(9)透射之后的光束和薄玻璃板前表面的反射光一起通过会聚透镜(10)汇聚至光电探测器(4)的光敏面上,所述光电探测器(4)输出电信号给信号处理系统(5);薄玻璃板(9)后表面和第二平面反射镜(6)的反射面之间的距离为d;所述第二平面反射镜(6)的背面中心与待测压电陶瓷管(7)的一端固定连接,该待测压电陶瓷管(7)的另一端固定在二维调整架(8)上,所述待测压电陶瓷管(7)的中心轴线与所述第二平面反射镜(6)的反射面相垂直;所述待测压电陶瓷管(7)的内表面(7-1)和外表面(7-2)分别通过电极(1)与高压电源的两个电压输出端连接;其特征在于所述多普勒振镜正弦调制多光束激光外差二次谐波测量电致伸缩系数的方法由如下步骤实现:首先,通过调整二维调整架(8),使与待测压电陶瓷管(7)固定连接的第二平面反射镜(6)的反射面与薄玻璃板(9)相互平行,并使第二平面反射镜(6)的反射面与薄玻璃板(9)之间的距离d为4.25mm;然后,采用高压电源为待测压电陶瓷管(7)提供驱动电压,并打开振镜(13)的驱动电源使振镜(13)开始振动;同时,打开H<sub>0</sub>固体激光器(2)。最后,调节所述高压电源的输出电压信号U,同时信号处理系统(5)连续采集光电探测器(4)输出的电信号,并对采集到的信号进行处理,进而获得第二平面反射镜(6)和薄玻璃板后表面(9)之间的距离变化量,根据该距离变化量和此时高压电源输出的电压信号获得待测压电陶瓷管(7)的电磁致伸缩系数:设用E表示待测压电陶瓷管(7)内外表面加上电压后,在待测压电陶瓷管(7)内外表面间形成的径向电场的电场强度,用ε表示待测压电陶瓷管(7)轴向的应变,α表示待测压电陶瓷管(7)在准线性区域内的电致伸缩系数,于是得到:ε=αE                                         公式1按上式有:<maths num="0001"><![CDATA[<math><mrow><mfrac><mi>&Delta;l</mi><mi>l</mi></mfrac><mo>=</mo><mi>&alpha;</mi><mfrac><mi>U</mi><msub><mi>d</mi><mn>0</mn></msub></mfrac></mrow></math>]]></maths>公式2最终得到:<maths num="0002"><![CDATA[<math><mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><msub><mi>ld</mi><mn>0</mn></msub></mrow><mi>lU</mi></mfrac></mrow></math>]]></maths>公式3公式中,Δl是待测压电陶瓷管(7)在加电前后的长度增量,即等于第二平面反射镜(6)和薄玻璃板(9)之间的距离变化量,l是待测压电陶瓷管(7)的未加电状态的原始长度;d<sub>0</sub>是待测压电陶瓷管(7)的壁厚;所述信号处理系统(5)根据连续采集光电探测器(4)输出的电信号,并对采集到的信号进行处理,进而获得第二平面反射镜(6)和薄玻璃板(9)之间的距离变化量的过程为:由于激光在薄玻璃板(9)前表面的反射光与第二平面反射镜(6)入射面反射k次和k+1次后的透射出玻璃前表面的光混频,产生两个幅度相差2~3个数量级的差频信号,所述方法的二次谐频差为检测到后表面k次反射的E<sub>k</sub>与后表面k+2次反射后的E<sub>k+2</sub>光混频所产生;不考虑薄玻璃板(9)自身厚度的情况下,当激光以入射角θ<sub>0</sub>斜入射薄玻璃板(9)前表面时的入射光场为E(t)=E<sub>l</sub>exp(iω<sub>0</sub>t),振镜(13)的简谐振动方程为x(t)=x<sub>0</sub>cos(ω<sub>c</sub>t);振镜(13)的速度方程为v(t)=-ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>t),反射光的频率变为ω=ω<sub>0</sub>(1-2ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>t)/c),上式中参数ω<sub>0</sub>为激光角频率,参数x<sub>0</sub>为振镜(13)振动的振幅,参数ω<sub>c</sub>为振镜(13)的角频率,c为光速,t为时间;则t-L/c时刻到达薄玻璃板(9)的反射光场为:E<sub>0</sub>(t)=α<sub>0</sub>E<sub>l</sub>exp{i[ω<sub>0</sub>(1-2ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>(t-L/c))/c)                公式4(t-L/c)+ω<sub>0</sub>x<sub>0</sub>cos(ω<sub>c</sub>(t-L/c))/c]}式中,参数α<sub>0</sub>=r,r为薄玻璃板(9)的反射系数;L为振镜(13)到不计厚度薄玻璃板(9)前表面的距离;E<sub>l</sub>为振幅常数;经薄玻璃板(9)透射的光在不同时刻被第二平面反射镜(6)后表面多次反射并多次透射出薄玻璃板(9)的前表面,其透射光的表达式分别写成如下形式:E<sub>1</sub>(t)=α<sub>1</sub>E<sub>l</sub>exp{i[ω<sub>0</sub>(1-2ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>(t-(L+2ndcosθ)/c))/c)(t-(L+2ndcosθ)/c)+ω<sub>0</sub>x<sub>0</sub>cos(ω<sub>c</sub>(t-(L+2ndcosθ)/c))/c]}E<sub>2</sub>(t)=α<sub>2</sub>E<sub>l</sub>exp{i[ω<sub>0</sub>(1-2ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>(t-(L+4ndcosθ)/c))/c)(t-(L+4ndcosθ)/c)+ω<sub>0</sub>x<sub>0</sub>cos(ω<sub>c</sub>(t-(L+4ndcosθ)/c))/c]}E<sub>3</sub>(t)=α<sub>3</sub>E<sub>l</sub>exp{i[ω<sub>0</sub>(1-2ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>(t-(L+6ndcosθ)/c))/c)(t-(L+6ndcosθ)/c)+ω<sub>0</sub>x<sub>0</sub>cos(ω<sub>c</sub>(t-(L+6ndcosθ)/c))/c]}              公式5E<sub>m</sub>(t)=α<sub>m</sub>E<sub>l</sub>exp{i[ω<sub>0</sub>(1-2ω<sub>c</sub>x<sub>0</sub>sin(ω<sub>c</sub>(t-(L+2mndcosθ)/c))/c)(t-(L+2mndcosθ)/c)+ω<sub>0</sub>x<sub>0</sub>cos(ω<sub>c</sub>(t-(L+2mndcosθ)/c))/c]}其中,下标m取值为0,1,2,......,n为薄玻璃板(9)和第二平面反射镜(6)之间介质的折射率,α<sub>1</sub>=β<sup>2</sup>r′,......,α<sub>m</sub>=β<sup>1</sup>r′<sup>m</sup>r<sup>m-1</sup>,β为薄玻璃板(9)的透射系数,r′为第二平面反射镜(6)的反射系数,参数d为薄玻璃板(9)到第二平面反射镜(6)的距离,θ为入射光透过薄玻璃板(9)后的折射角,光电探测器(4)接收到的总光场表示为:E(t)=E<sub>0</sub>(t)+E<sub>1</sub>(t)+…+E<sub>m</sub>(t)                                 公式6则光电探测器(4)输出的光电流表示为:<maths num="0003"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><mfrac><mi>&eta;e</mi><mi>hv</mi></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><mo>&Integral;</mo><munder><mo>&Integral;</mo><mi>S</mi></munder><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><msub><mi>E</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>E</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>]</mo><msup><mrow><mo>[</mo><msub><mi>E</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>E</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>]</mo></mrow><mo>*</mo></msup><mi>ds</mi></mrow></math>]]></maths>公式7其中,参数e为电子电量,参数Z为探测器表面介质的本征阻抗,参数η为量子效率,参数S为探测器光敏面的面积,参数h为普朗克常数,参数v为激光频率,*号表示复数共轭;整理得到激光外差二次谐波信号的中频电流为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>I</mi><mi>if</mi></msub><mo>=</mo><mfrac><mi>&eta;e</mi><mrow><mn>2</mn><mi>hv</mi></mrow></mfrac><mfrac><mn>1</mn><mi>Z</mi></mfrac><mo>&Integral;</mo><munder><mo>&Integral;</mo><mi>S</mi></munder><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>0</mn></mrow><mo>&infin;</mo></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>p</mi><mo>+</mo><mn>2</mn></mrow><mo>&infin;</mo></munderover><mrow><mo>(</mo><msub><mi>E</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><msub><mi>E</mi><mi>j</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><msub><mi>E</mi><mi>p</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>E</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>ds</mi></mrow></math>]]></maths>公式8将公式4和公式5代入公式8,结果为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>I</mi><mi>IF</mi></msub><mo>=</mo><mfrac><mi>&eta;e</mi><mi>hv</mi></mfrac><mfrac><mi>&pi;</mi><mi>Z</mi></mfrac><msubsup><mi>E</mi><mn>0</mn><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>m</mi><mo>-</mo><mi>p</mi></mrow></munderover><msub><mi>&alpha;</mi><mrow><mi>j</mi><mo>+</mo><mi>p</mi></mrow></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><mi>cos</mi><mo>[</mo><mfrac><mrow><mn>8</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi><msub><mi>&omega;</mi><mn>0</mn></msub><msubsup><mi>&omega;</mi><mi>c</mi><mn>2</mn></msubsup><msub><mi>x</mi><mn>0</mn></msub></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><mi>t</mi><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>&omega;</mi><mn>0</mn></msub><msub><mi>x</mi><mn>0</mn></msub></mrow><mi>c</mi></mfrac><mo>-</mo><mfrac><mrow><mn>4</mn><mi>nd</mi><msub><mi>&omega;</mi><mn>0</mn></msub><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac></mrow></math>]]></maths>公式9<maths num="0006"><![CDATA[<math><mrow><mo>-</mo><mfrac><mrow><mn>8</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi><msub><mi>&omega;</mi><mn>0</mn></msub><msubsup><mi>&omega;</mi><mi>c</mi><mn>2</mn></msubsup><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>L</mi><mo>+</mo><mn>2</mn><mi>pnd</mi><mi>cos</mi><mi>&theta;</mi><mo>)</mo></mrow></mrow><msup><mi>c</mi><mn>3</mn></msup></mfrac><mo>]</mo></mrow></math>]]></maths>忽略1/c<sup>3</sup>的小项之后简化为:<maths num="0007"><![CDATA[<math><mrow><msub><mi>I</mi><mi>IF</mi></msub><mo>=</mo><mfrac><mi>&eta;e</mi><mi>hv</mi></mfrac><mfrac><mi>&pi;</mi><mi>Z</mi></mfrac><msubsup><mi>E</mi><mn>0</mn><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>m</mi><mo>-</mo><mi>p</mi></mrow></munderover><msub><mi>&alpha;</mi><mrow><mi>j</mi><mo>+</mo><mi>p</mi></mrow></msub><msub><mi>&alpha;</mi><mi>j</mi></msub><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>8</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi><msub><mi>&omega;</mi><mn>0</mn></msub><msubsup><mi>&omega;</mi><mi>c</mi><mn>2</mn></msubsup><msub><mi>x</mi><mn>0</mn></msub></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><mi>t</mi><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>&omega;</mi><mn>0</mn></msub><msub><mi>x</mi><mn>0</mn></msub><mo>-</mo><mn>4</mn><mi>nd</mi><msub><mi>&omega;</mi><mn>0</mn></msub><mi>cos</mi><mi>&theta;</mi></mrow><mi>c</mi></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>公式10其中,参数p和j均为非负整数;根据公式10,把激光外差二次谐波信号的频率记为:<maths num="0008"><![CDATA[<math><mrow><mi>f</mi><mo>=</mo><mn>8</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi><msub><mi>&omega;</mi><mn>0</mn></msub><msubsup><mi>&omega;</mi><mi>c</mi><mn>2</mn></msubsup><msub><mi>x</mi><mn>0</mn></msub><mo>/</mo><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><msup><mi>c</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>=</mo><mn>4</mn><mi>nd</mi><mi>cos</mi><mi>&theta;</mi><msub><mi>&omega;</mi><mn>0</mn></msub><msubsup><mi>&omega;</mi><mi>c</mi><mn>2</mn></msubsup><msub><mi>x</mi><mn>0</mn></msub><mo>/</mo><mrow><mo>(</mo><mi>&pi;</mi><msup><mi>c</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>=</mo><mi>Kd</mi></mrow></math>]]></maths>公式11根据公式11得出,激光外差二次谐波信号的频率与薄玻璃板(9)和第二平面反射镜(6)之间的距离d成正比,比例系数为:<maths num="0009"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><mn>4</mn><mi>n</mi><mi>cos</mi><mi>&theta;</mi><msub><mi>&omega;</mi><mn>0</mn></msub><msubsup><mi>&omega;</mi><mi>c</mi><mn>2</mn></msubsup><msub><mi>x</mi><mn>0</mn></msub><mo>/</mo><mrow><mo>(</mo><mi>&pi;</mi><msup><mi>c</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></math>]]></maths>公式12通过公式11测出薄玻璃板(9)和第二平面反射镜(6)之间的距离d,当d改变时,根据公式11测出对应d的变化量Δd,得到Δd根据<img file="FDA0000065401360000044.GIF" wi="207" he="110" />获得待测压电陶瓷管(7)电致伸缩系数。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号