CN1193731A - 全角跟踪振动旋转传感器 - Google Patents

全角跟踪振动旋转传感器 Download PDF

Info

Publication number
CN1193731A
CN1193731A CN98105307.6A CN98105307A CN1193731A CN 1193731 A CN1193731 A CN 1193731A CN 98105307 A CN98105307 A CN 98105307A CN 1193731 A CN1193731 A CN 1193731A
Authority
CN
China
Prior art keywords
resonator
voltage
value
angle
exciting
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
Application number
CN98105307.6A
Other languages
English (en)
Inventor
A·马修斯
G·T·凡提
李春明
D·D·凌奇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northrop Grumman Guidance and Electronics Co Inc
Original Assignee
Litton Systems Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Litton Systems Inc filed Critical Litton Systems Inc
Publication of CN1193731A publication Critical patent/CN1193731A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C19/00Gyroscopes; Turn-sensitive devices using vibrating masses; Turn-sensitive devices without moving masses; Measuring angular rate using gyroscopic effects
    • G01C19/56Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces
    • G01C19/567Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces using the phase shift of a vibration node or antinode
    • G01C19/5691Turn-sensitive devices using vibrating masses, e.g. vibratory angular rate sensors based on Coriolis forces using the phase shift of a vibration node or antinode of essentially three-dimensional vibrators, e.g. wine glass-type vibrators

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)
  • Transmission And Conversion Of Sensor Element Output (AREA)
  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

本发明提供一种振动旋转传感器,它包括:共振器,共振器为旋转对称的薄壁物体,能够以多种驻波模式振动,驻波相对共振器上参考点的方向由方位角定义,一个或多个电极附着在共振器表面,一个或多个电极与单个输出端口电学连接;具有与共振器输出端口相连的输入端口的检测电子线路,检测电子线路从共振器输出端口的共振器信号中获取驻波方位角与跟踪角之间的差值,跟踪角相对共振器的参考点而言。

Description

全角跟踪振动旋转传感器
本发明的主题与专利申请“加交流作用力和检测电子线路的振动旋转传感器”(Kumar和Foster)、“具有多路复用电子线路的振动旋转传感器”(Matthews,Darling和Varty)和“施加交流作用力电压的振动旋转传感器”(Lynch)中揭示的发明主题相同。
本发明通常涉及振动旋转传感器,特别是涉及与这类旋转传感器相连的电子线路。
图1示出了组装前的现有技术振动旋转传感器(VRS)10,它由外层部件12、半球形共振器14和内层部件16组成,所有这些部分由熔凝石英构成并通过铟结合在一起。惰性敏感元件为薄壁状直径为5.8cm的半球形共振器14,它位于外层部件12与内层部件16之间并且由杆26支承。
一个环状冲头电极20和16个分立的冲头电极22被淀积在外层部件12的内表面上。在组装好的VRS10中,环状冲头电极20和16个分立的冲头电极22紧靠在半球形共振器14的外部金属化表面32附近。在组装好的VRS中,淀积在内层部件16上的8个敏感元件电极24紧靠在半球形共振器14的内部金属化表面30附近。
电容力借助半球形共振器与环状冲头电极20之间适当的作用力电压施加在半球形共振器14上从而使其以最低阶非伸缩(或柔性)模式振动。建立起了具有沿四周以90度相隔的四个反节点的驻波,四个节点与反节点偏离45度。0度和180度反节点与90度和270度反节点形成相差90度的振荡。驻波使半球形共振器的边缘形状从圆形变化为椭圆状(半主轴通过0度/180度的反节点和90度/270度的反节点)。
VRS10围绕垂直于半球形共振器边缘34平面的轴转动,使得驻波相对VRS沿反方向转动,其转动角度正比于VRS10的转动角。这样,通过测量驻波相对VRS10的转动角度,可以确定VRS10的转动角度。
通过在半球形共振器14上施加直流偏压并在环状冲头电极20上施加交流电压激发了半球形共振器14的振动模式,其中交流电压的频率是半球形共振器14共振频率的两倍。
当半球形共振器14振动并且敏感元件电极24相对半球形共振器的电容发生变化时,通过测量流入流出敏感元件电极24的电流确定了驻波相对VRS10的模式角。x轴信号Ix源于组合I0-I90+I180-I270,这里的下标表示相对电流源电极x轴的角度定向。同样,y轴信号Iy源于组合I45-I135+I225-I315。两倍于驻波相对0度(即x)轴模式角的正切值等于Iy/Ix.。
由于半球形共振器14厚度的不均匀,第一驻波的建立将会导致第二驻波的生成,第二驻波以正交相位振荡,其反节点与第一驻波的节点重合。通过在16个分立冲头电极22上施加合适的电压可以阻止第二驻波的形成。
为了降低施加在环状冲头电极20和分立冲头电极22上的交流作用力电压的幅值并且使施加在共振器上的力相对交流驱动电压呈线性变化,直流偏压通常维持在半球形共振器14上。直流偏压的存在导致VRS电学性质变化缓慢,这归因于外层部件12和内层部件16上或内部发生的电荷迁移现象所引起的电容变化。这些缓慢变化导致时间性能令人无法接受地大幅度下降,因此需要提供特殊的装置来补偿这些效应。
本发明提供一种包含共振器和附着共振器的壳罩的振动旋转传感器,本发明还提供一种利用大小等于平均方位角的跟踪角读取这种传感器的驻波方位角的方法。共振器为旋转对称的薄壁物体,可以以多种驻波模式振动。驻波的方向由驻波某一反节点轴相对共振器上固定点的方位角限定。
在共振器表面上淀积有一个或多个电极并且与单输出端口电学连接。壳罩包括多个靠近一个或多个共振器电极的淀积电极。
读取特定驻波方位角的方法包括以下步骤:产生多个激发电压;将激发电压施加到壳罩电极上;以及通过对出现在共振器输出端口的共振器信号进行运算处理确定方位角与合成跟踪角之差。
施加在壳罩电极上的驱动电压借助壳罩电极与共振器电极之间的电容连到共振器输出口。驱动电压包括一对激发电压或作用力电压,或者包括全部二者。当激发电压到达共振器输出端口时,它对共振器的动力学基本上没有影响,但是却携带了与跟踪角有关的信息和驻波参数。作用力电压使力施加在共振器上并改变了共振器动力学性能和驻波参数。
施加在壳罩电极上的驱动电压与单个共振器信号结合在一起,这是通过壳罩-电极-共振器-电极构成的电容向共振器输出端口发送的结果。为了确定驻波方位角与跟踪角之差,激发和作用力信号被设计成可以通过对共振器信号作适当的运算而分离开来。
激发和作用力电压可以采用不同的形式构成。频分多路复用途径使得激发电压局限在分离的频带内而作用力电压的频谱局限在与激发电压有关的频带分离的频带内。
相分多路复用途径使得激发电压为频率相同而相位相差1/4周期的时间周期函数,作用力电压的频谱局限在与激发电压频率分离的频带内。
一种时分多路复用途径使得激发电压正比于取值为0和1的单一方波,而每个作用力电压包括正比于取值为0和1的方波的倍增因子,其中在任何给定时刻只有与激发和作用力电压相连的方波才取值为1。
第二种时分多路复用途径使得每个激发电压正比于具有预定频率和相位的时间周期函数与取值为0和1的独特方波之积,而每个作用力电压包括正比于取值为0和1的方波的倍增因子,其中在任何给定时刻只有与激发和作用力电压相连的一个方波才取值为1。码分多路复用途径使得激发电压正比于按照预定伪随机序列取值为1和-1的独特方波,作用力电压的频谱局限在与激发电压相关频段分离的频带内。
从共振器信号中确定驻波方位角的过程包括首先从共振器信号中提取至少两个分量并随后通过对两个分量所作的运算确定方位角与跟踪角之差。在频分多路复用情形下,两个分量占据分离的频带,通过对共振器信号所作的区分分离频带的运算提取每个分量。
在相分多路复用情形下,两个分量是频率相同而相位相差1/4周期的周期函数,通过对共振器信号所作的区分两分量相位的运算提取每个分量。
在时分多路复用情形下,两个分量存在于不同时间间隔的共振器信号内,通过对共振器信号所作的区分不同时间间隔的运算提取每个分量。
在码分多路复用情形下,两个分量是0和1的伪随机序列,伪随机序列的交叉相关度为零,通过对共振器信号所作的区分两伪随机序列的运算提取每个分量。
跟踪角不断被调整从而使方位角与跟踪角之差平均为零。通过将跟踪角加上方位角与跟踪角之差计算得到方位角。
图1示出了已有技术振动旋转传感器的结构。
图2示出了本发明的控制和读取电子线路的框图。
图3示出了本发明特定实施例的多路复用控制信号。
本发明为用多路信号实现控制和读取的振动旋转传感器。本发明的振动旋转传感器由共振器、附着共振器的壳罩和多路复用电子线路组成。共振器可以是任一具有驻波振动模式的旋转对称薄壁物体。现有技术的共振器一般为半球形。
图2示出了确定驻波参数并控制共振器动力学的简单方法。驻波相对固定于共振器的x和y轴描述。同相驻波相对共振器的方向可以用从x轴顺时针方向测量的同相反节点轴的方位角θ来定义。共振器边缘沿同相反节点轴与圆的偏离假定随cos(ωt+φ)变化,ω为振动频率,t为时间,而φ为任意相位角。四极驻波相对共振器的取向用从x轴顺时针方向测量的四极反节点轴方位角θ+π/4来定义。共振器边缘沿四极反节点轴与圆的偏离假定随sin(ωt+φ)变化。
淀积在共振器内表面上的连续圆周共振器电极42施以直流电压VB并且通过直流阻断电容器43与放大器-多路分配器44连接。附着在VRS壳罩上的8根电极46在紧靠共振器电极42附近的四周等距分开,最上面的xp电极位于x轴的中心。8根电极46的驱动电压Vxp(t)、Vxn(t)、Vyp(t)和Vyn(t)取自多路复用器48:
Vxp(t)=Vmxr(t)cos2θrcos(ωxrt+ψxr)-Vmyr(t)sin2θrcos(ωyrt+ψyr)+Vcx(t)Uxp(t)
Vxn(t)=-Vmxr(t)cos2θrcos(ωxrt+ψxr)-Vmyr(t)sin2θrcos(ωyrt+ψyr)+Vcx(t)Uxn(t)
Vyp(t)=Vmxr(t)sin2θrcos(ωxrt+ψxr)-Vmyr(t)cos2θrcos(ωyrtyr)+Vcy(t)Uyp(t)    (1)
Vyn(t)=-Vmxr(t)sin2θrcos(ωxrt+ψxr)-Vmyr(t)cos2θrcos(ωyrt+ψyr)+Vcy(t)Uyn(t)激发电压Vmxr(t)cos(ωxrt+φxr)和Vmyr(t)cos(ωyrt+φyr)是图2xr-yr跟踪角坐标系的分量(用下标r表示)。激发电压较佳实施例包括正弦cos(ωxrt+φxr)和cos(ωyrt+φyr)。有各种周期函数F(ωxrt+φxr)可以替代正弦函数,包括普通的方波。
xr轴沿x轴以跟踪角θr顺时针旋转。激发电压被设计为不会影响共振器上驻波的参量。角频率ωxr和ωyr以及相位φxr和φyr取决于所用的多路复用类型。作用力电压Vcx(t)Uxp(t)、Vcx(t)Uxn(t)、Vcy(t)Uyp(t)和Vcy(t)Uyp(t)(表示为x-y坐标系的分量)使得作用力施加在共振器上以控制共振器上一个或多个驻波的参量。函数Uxp(t)、Uxn(t)、Uyp(t)和Uyp(t)由控制单元50产生并提供给多路复用器48。电压Vcx(t)和Vcy(t)为预先确定的函数,用来隔离作用力电压与激发电压。
从共振器电极42流入放大器-多路分配器44的电流I(t)给定如下:
I(t)=Ixp(t)+Ixn(t)+Iyp(t)+Iyn(t)                     (2)这里
Ixp(t)=KI[Vmxr(t)ωxrcos2θrcos(ωxrt+ψxr)-Vmyr(t)ωyrsin2θrcos(ωyrt+ψyr)
                            +Vcx(t)ωUxpUxp(t)]Cxp
Ixn(t)=KI[-Vmxr(t)ωxrcos2θrcos(ωxrt+ψxr)+Vmyr(t)ωyrsin2θrcos(ωyrt+ψyr)
                            +Vcx(t)ωUxnUxp(t)]Cxn
Iyp(t)=KI[Vmxr(t)ωxrsin2θrcos(ωxrt+ψxr)-Vmyr(t)ωyrcos2θrcos(ωyrt+ψyr)
                            +Vcy(t)ωUypUyp(t)]Cyp    (3)
Iyn(t)=KI[-Vmxr(t)ωxrsin2θrcos(ωxrt+ψxr)+Vmyr(t)ωyrsin2θrcos(ωyrt+ψyr)
                            +Vcy(t)ωUynUyp(t)]Cyn电容Cxp、Cxn(t)、Cyp和Cyp是xp、xn、yp和yn电极46相对共振器电极42的电容。角频率ωUxp、ωUxn、ωUyp和ωUyp是与相应的Us关联的角频率并且等于或小于2ω,这里ω为共振器振动频率符号K1表示常数。驱动电压与得出电流之间的相位差没有相关性并且在上述方程中可以忽略。电容给定如下:
Cxp=Co[1+dicos2θcos(ωt+φ)-dqsin2θsin(ωt+φ)]
Cxn=Co[1-dicos2θcos(ωt+φ)+dqsin2θsin(ωt+φ)]
Cyp=Co[1+disin2θcos(ωt+φ)-dqcos2θcos(ωt+φ)]    (4)
Cyn=Co[1-disin2θcos(ωt+φ)+dqcos2θcos(ωt+φ)]这里省略了涉及di和dq的高次项。高次项的影响在后面的处理操作中考虑。量Co为共振器未激发时电极对的电容,di和dq分别为共振器未激发时共振器电极42与电极46之间缝隙除同相和四极模式的最大弯曲幅度,θ为同相驻波的反节点与x轴之间的角度,ω为共振器振动角频率,以及φ为任意的相位角。
在电流方程中代入电容表达式并求和得到的I为:
I(t)=2KICoVmxr(t)ωxrcos(ωxrt+ψxr)dicos(2θ-2θr)cos(ωt+φ)
      -2KICoVmxr(t)ωxrcos(ωxrt+ψxr)dqsin(2θ-2θr)sin(ωt+φ)
      +KIVcx(t)ωUxpUxp(t)Cxp+KIVcx(t)ωUxnUxn(t)Cxn
      +2KICoVmyr(t)ωyrcos(ωyrt+ψyr)disin(2θ-2θr)cos(ωt+φ)
      +2KICoVmyr(t)ωyrcos(ωyrt+ψyr)dqcos(2θ-2θr)sin(ωt+φ)    (5)
      +KIVcy(t)ωUypUyp(t)Cyp+KIVcy(t)ωUynUyn(t)Cyn通过放大器-多路分配器44将电流I(t)变换为电压V(t):
V(t)=KV[Vx(t)Rx(t)+Vy(t)Ry(t)]+KF[Fx(t)+Fy(t)]                     (6)这里Kv和KF为常数,
Vx(t)=Vmxr(t)ωxrcos(ωxrt+ψxr)
Vy(t)=Vmyr(t)ωyrcos(ωyrt+ψyr)
Rx(t)=dicos(2θ-2θr)cos(ωt+φ)-dqsin(2θ-2θr)sin(ωt+φ)
Ry(t)=disin(2θ-2θr)cos(ωt+φ)+dqsin(2θ-2θr)sin(ωt+φ)        (7)以及
Fx(t)=Vcx(t)[ωUxpUxp(t)CxpUxnUxn(t)Cxn]
Fy(t)=Vcy(t)[ωUypUyp(t)CypUynUyn(t)Cyn]                        (8)由于包含了驻波参量di、dq、θr-θ、ω和φ,所以信号Rx(t)和Ry(t)为所需的多路分配过程输出,该过程由施加电压V(t)的分立操作组成。
包含信号信号Rx(t)和Ry(t)的信号Sx(t)和Sy(t)由放大器-多路分配器44通过对信号Sx(t)操作Ox和对Sy(t)操作Oy提取。放大器-多路分配器44中多路分配器部分的操作原理取决于电压Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)的形式和ωxr、ωyr、φxr和φyr的数值。
对于频分多路复用来说,Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)全部为常数,ωxr、ωyr和|ωxryr|大于6ω左右,而φxr和φyr为任意的常数。包含驻波参量的信号Rx(t)和Ry(t)可以通过对两个V(t)的乘积解调得到,一个是相对cos(ωxrt+φxr),另一个是相对cos(ωyrt+φyr)。如果采用正弦以外的周期函数,则采用周期函数的副本进行解调。乘积解调包括将输入电压乘以基准正弦函数(或副本)和低通滤波乘积,低通滤波器的截止频率为3ω左右。上述处理的结果为信号SFDMx(t)和SFDMy(t):
SFDMx(t)=KFDMRx(t)
SFDMy(t)=KFDMRy(t)                                                  (9)这里KFDM为常数。由于频谱Fx(t)和Fy(t)的上限为3ω,所以通过多路分配去除了这些量。
对于相分多路复用,ωxr和ωxr具有相同的值ωo,ωo大于6ω,并且φxryr等于π/2圆弧。信号SFDMx(t)和SFDMy(t)通过相对cos(ωxrt+φxr)和cos(ωyrt+φyr)对V(t)进行乘积解调得出(或者相对所用周期函数的副本)。
SPDMx(t)=KPDMRx(t)
SPDMy(t)=KPDMRy(t)                                                 (10)这里KPDM为常数。
对于时分多路复用的一种形式,ωxr和ωxr具有相同的值ωo,ωo大于6ω,并且φxr和φyr等于任意的相位φo。电压Vmxr(t)和Vmyr(t)正比于取值为0和1的方波,只有其中一个在任意给定时刻等于1并且“1”的间隔等于2π/ω的整数倍。电压Vcx(t)和Vcy(t)等于常数。信号STDMx(t)和STDMy(t)通过相对于cos(ωot+φo)(或副本)、对V(t)执行乘积解调后与Vmxr(t)和Vmyr(t)并行相乘而得到。
STDMx(t)=KTDMVmx(t)Rx(t)
STDMy(t)=KTDMVmy(t)Ry(t)                                           (11)这里KTDM为常数。值得注意的是,只有当Vmxr(t)和Vmyr(t)非零时才可以使用Rx(t)和Ry(t)。
如果Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)正比于取值为0和1的方波,该方波只有其中一个在任意给定时刻等于1并且“1”的间隔等于2π/ω的整数倍,则可以得到同样的结果(除了常数KTDM值以外)。由于将作用力电压Vcx(t)Uxp(t)、Vcx(t)Uxn(t)、Vcy(t)Uyp(t)和Vcy(t)Uyn(t)完全相互隔离并与激发电压Vmax(t)cos(ωxrt+φxr)和Vmyr(t)cos(ωyrt+φyr)隔离开来,所以这种操作模式是期望的。
对于另一种形式的时分多路复用,ωo为零并且电压Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)正比于取值为0和1的方波,该方波只有其中一个在任意给定时刻等于1并且“1”的间隔等于2π/ω的整数倍时。将V(t)与Vmxr(t)和Vmyr(t)相乘得到与第一种时分多路复用形式相同的结果。
对于码分多路复用,ωxr、ωyr、φxr和φyr为零,Vcx(t)和Vcy(t)为常数,Vmxr(t)和Vmyr(t)正比于方波,该方波对伪随机序列取值-1/T和1/T并满足下列条件:
Figure A9810530700161
这里下标i和j表示下标mxr、myr、cx和cy中的任意一个。积分时间周期T小于2π/3ω。信号SCDMx(t)和SCDMy(t)通过将V(t)分别与Vmxr(t)和Vmyr(t)相乘并对T积分得到:
SCDMx(nT)=KCDMRx(nT)
SCDMy(nT)=KCDMRy(nT)                                               (13)这里KCDM为常数并且n为整数。值得注意的是,信号SCDMx(t)和SCDMy(t)提供了有关T间隔内驻波参量的信息。
电压Ux(t)和Uy(t)一般包括三种分量:
Uxp(t)=Uaxp(t)+Uqxp(t)+Urxp(t)
Uxn(t)=Uaxn(t)+Uqxn(t)+Urxn(t)
Uyp(t)=Uayp(t)+Uqyp(t)+Uryp(t)                                    (14)
Uyn(t)=Uayn(t)+Uqyn(t)+Uryn(t)这里下标a、q和r表示幅度、四极和速率控制电压。在所有的应用中无需将这些分量相互隔离。但是如果需要隔离,则可以在上述方程中作下列替换。
Vcax(t)Uaxp+Vcqx(t)Uqxp(t)+Vcrx(t)Urxp(t)对Vcx(t)Uxp(t)
Vcax(t)Uaxn+Vcqx(t)Uqxn(t)+Vcrx(t)Urxn(t)对Vcx(t)Uxn(t)
Vcay(t)Uayp+Vcqy(t)Uqyp(t)+Vcry(t)Uryp(t)对Vcy(t)Uyp(t)            (15)
Vcay(t)Uayn+Vcqy(t)Uqyn(t)+Vcry(t)Uryn(t)对Vcy(t)Uyn(t)经过替换,任何对Vcx(t)和Vcy(t)的限制同样可以用于Vcax(t)、Vcqx(t)、Vcrx(t)、Vcay(t)、Vcqy(t)和Vcry(t)上。例如,方程(1)变为
Vxp(t)=Vmxr(t)cos2θrcos(ωxrt+ψxr)-Vmyr(t)sin2θrcos(ωxrt+ψxr)
        +Vcax(t)Uaxp(t)+Vcqx(t)Uqxp(t)+Vcrx(t)Urxp(t)
Vxp(t)=-Vmxr(t)cos2θrcos(ωxrt+ψxr)+Vmyr(t)sin2θrcos(ωyrt+ψyr)
        +Vcax(t)Uaxn(t)+Vcqx(t)Uqxn(t)+Vcrx(t)Urxn(t)
Vyp(t)=Vmxr(t)sin2θrcos(ωxrt+ψxr)-Vmyr(t)cos2θrcos(ωyrt+ψyr)
        +Vcay(t)Uayp(t)+Vcqy(t)Uqyp(t)+Vcry(t)Uryp(t)              (16)
Vyp(t)=-Vmxr(t)sin2θrcos(ωxrt+ψxr)+Vmyr(t)sin2θrcos(ωyrt+ψyr)
        +Vcay(t)Uayn(t)+Vcqy(t)Uqyn(t)+Vcyx(t)Uryn(t)
一种可能的时分多路复用结构是与共振器柔性速率同步的间隔32π/ω的16时隙帧。多路复用控制电压如图3所示。当θ等于θr时,xr轴与反节点轴一致并且yr轴与节点轴一致。八个时隙分配给读取yr信号分量,四个时隙读取xr信号分量,并且幅度、四极和对共振器的速率作用力各分配1个时隙。对于4kHz的振动频率而言,可以分别在2kHz和1kHz速率下读取xr和yr信号分量。控制电压的施加速率为0.25kHz。
总之,从放大器-多路分配器44输出的信号Sx(t)和Sy(t)具有如下形式:
Sx(t)=KVxRx(t)
Sy(t)=KVyRxy(t)                                                   (17)这里除了时分多路复用以外,当Kvx等于KvVmx(t)而Kvy等于KvVmy(t)时,Kvx和Kvy等于Kv
为了从信号Sx(t)和Sy(t)中提取驻波参量,需要稳定和精确的共振器振动信号cos((ωt+φ)的副本。副本取自副本生成器52中的电压受控振荡器,其中电压受控振荡器对同相驻波反节点信号是锁相的。处理的第一个步骤是首先将Sx(t)和Sy(t)与副本信号cos(ωrt+φr)相乘并低通滤波结果,随后乘以相移副本sin(ωrt+φr)并低通滤波结果。该过程的结果为:
    K{dicos2θcos[(ωr-ω)t+(φr-φ)]Six(t)=
    +dqsin2θsin[(ωr-ω)t+(φr-φ)]}
    K{disin2θcos[(ωr-ω)t+(φr-φ)]Siy(t)=
    -dqcos2θsin[(ωr-ω)t+(φr-φ)]}
    K{dicos2θsin[(ωr-ω)t+(φr-φ)]Sqx(t)=                                                               (18)
    -dqsin2θcos[(ωr-ω)t+(φr-φ)]}
    K{disin2θsin[(ωr-ω)t+(φr-φ)]Sqy(t)=
    +dqcos2θcos[(ωr-ω)t+(φr-φ)]}这里K为常数。下一步骤是形成以下Six(t)、Siy(t)、Sqx(t)和Sqy(t)的乘积组合: E = S ix 2 + S qx 2 + S iy 2 + S qy 2 = K 2 ( d i 2 + d q 2 ) Q=2(SixSqy-SiySqx)=K2(2didq) R = S ix 2 + S qx 2 - S iy 2 - S qy 2 = K 2 ( d i 2 - d q 2 ) cos 4 θ - - - - ( 19 ) S = 2 ( S ix S iy + S qx ) S qy ) = K 2 ( d i 2 - d q 2 ) sin 4 θ L i = 2 ( S ix S qx + S iy S qy ) = K 2 ( d i 2 - d q 2 ) sin [ 2 ( ω r - ω ) t + 2 ( φ r - φ ) ] 这里Li(t)为误差信号,锁相环路将副本相位φr锁定为φ而ωr锁定为ω。
驻波方位角与跟踪角之间的差值θ-θr可以由下式确定: tan ( 4 θ - 4 θ r ) = S ( t ) R ( t ) - - - - ( 20 ) 并得出Six(t)与Siy(t)泊符号。量S(t)可以用作控制环路中的误差信号,该环路产生θr并使平均θ等于θr和d/dt(θ-θr)=0。数字同步的跟踪角θr用来生成sinθr和cosθr,它们被提供给多路复用器48。任意给定时刻给出的实际角度为 θ = θ r + 1 4 tan - 1 S ( t ) R ( t ) ≈ θ r + 1 4 S ( t ) R ( t ) - - - - - ( 21 )
E(t)与特定数之间的差值用作幅度控制环路中的误差信号,该环路使得正比于di 2+dq 2的复合同相和四极驻波的总能量等于特定数。
量Q(t)用作四极控制环路中的误差信号,该环路使得四极驻波幅度dq为零。当环路是闭合的时候,幅度控制环路使同相幅度di保持为特定的数值。
上述控制变量的应用经过证明是可以优化的。对于本领域内普通技术人员来说是显而易见的,即有许多控制变量的选用是次最佳的,但仍是实用的。
控制单元50的输出是函数Uxp(t)、Uxn(t)、Uyp(t)和Uyn(t)与θ的sin和cos,它们被提供给多路复用器48。
在Loper等人于1990年8月28日提交的美国专利4,951,508中包含有振动旋转传感器的其它细节,它们作为参考文献包含在这里。

Claims (46)

1.一种振动旋转传感器,其特征在于包括:
共振器,共振器为旋转对称的薄壁物体,能够至少以多种驻波模式中的一种振动,驻波相对共振器上参考点的方向由方位角定义,一个或多个电极附着在共振器表面,一个或多个电极与单个输出端口电学连接;
具有与共振器输出端口相连的输入端口的检测电子线路,检测电子线路从共振器输出端口的共振器信号中获取驻波方位角与跟踪角之间的差值,跟踪角相对共振器的参考点而言。
2.如权利要求1所述的振动旋转共振器,其特征在于信号和Vx(t)Rx(t)+Vy(t)Ry(t)包含在共振器信号之中,Vx(t)和Vy(t)是预先确定的时间t的函数而Rx(t)和Ry(t)是时间、驻波方位角与跟踪角之间的差值Δθ和其它一个或多个驻波参量的函数,检测电子线路分别对共振器信号完成Ox和Oy操作以分别获取经过标度修正的Rx(t)和Ry(t)。
3.如权利要求2所述的振动旋转传感器,其特征在于Ox包括将共振器信号与角频率为ωx而相位为φx的周期函数副本相乘并随后进行低通滤波,而Oy包括将共振器信号与角频率为ωy而相位为φy的周期函数副本相乘并随后进行低通滤波,值ωx、φx、ωy和φy都是预先确定的。
4.如权利要求2所述的振动旋转传感器,其特征在于Ox包括将共振器信号与角频率为ωo而相位为φo的周期函数副本相乘并随后进行低通滤波,而Oy包括将共振器信号与相位为(φo+π/2)的周期函数副本相乘并随后进行低通滤波,值ωo和φo都是预先确定的。
5.如权利要求2所述的振动旋转传感器,其特征在于Ox包括将共振器信号与Vmxr(t)相乘,而Oy包括将共振器信号与Vmyr(t)相乘,这里Vmxr(t)和Vmyr(t)正比于预先确定的取值为0和1的时间方波函数,方波函数在同一时刻不能同时取值为1。
6.如权利要求2所述的振动旋转传感器,其特征在于:(1)Ox包括(a)将共振器信号与周期函数副本相乘,随后进行低通滤波,以及(b)与Vmxr(t)相乘,以及(2)Oy包括(a)将共振器信号与周期函数副本相乘并随后低通滤波,以及(b)与Vmyr(t)相乘,这里周期信号副本的角频率为ωo而相位为φo,它们都是预先确定的,Vmxr(t)和Vmyr(t)正比于预先确定的取值为0和1的时间方波函数,方波函数在同一时刻不能同时取值为1。
7.如权利要求2所述的振动旋转传感器,其特征在于Ox包括将共振器信号与Vmxr(t)相乘并随后在预先确定的时间间隔T上积分,而Oy包括将共振器信号与Vmyr(t)相乘而后在时间间隔T上积分,这里Vmxr(t)和Vmyr(t)正比于在时间间隔T内取值为-1和1序列的时间方波函数。
8.如权利要求2所述的振动旋转传感器,其特征在于Rx(t)和Ry由下列方程给定:
Rx(t)=di cos(2Δθ)cos(ωt+φ)-dqsin(2Δθ)sin(ωt+φ)
Ry(t)=di sin(2Δθ)cos(ωt+φ)+dqcos(2Δθ)sin(ωt+φ)这里di和dq分别是同相和四极振动模式的柔性幅度,ω为共振器振动角频率而相位φ为振动相位。
9.如权利要求1所述的振动旋转共振器,其特征在于进一步包括:
附着共振器的壳罩,所述壳罩包括多个靠近一个或多个共振器电极的附着电极;
向第一壳罩电极提供电压Vx1(t)并向第二壳罩电极提供电压Vy2(t)的驱动电子线路,Vx1(t)包括电压分量Vmxr(t)cos(2θr)F(ωxrt+φxr)和Vmyr(t)cos(2θr)F(ωyrt+φyr)而Vy2(t)包括电压分量Vmxr(t)cos(2θr)F(ωxrt+φxr)和Vmyr(t)cos(2θr)F(ωyrt+φyr),这里Vmxr(t)和Vy2(t)是预先确定的时间t函数或者常数,θr为跟踪角,F(ωt+φ)为时间t的周期函数,频率为ω而相位为φ,ωxr、φxr、ωyr和φyr是预先确定的,Vmxr(t)F(ωxrt+φxr)和Vmyr(t)F(yryr)对共振器的驻波动力学没有大的影响。
10.如权利要求9所述的振动旋转传感器,其特征在于Vmxr和Vmyr为常数,ωxr、ωxr以及|ωxrxr|大于6ω,ω为共振器的振动角频率。
11.如权利要求9所述的振动旋转传感器,其特征在于Vmxr和Vmyr为常数,ωxr和ωxr等于预先确定的值ωo,φxryr等于π/2圆弧,ωo大于6ω,ω为共振器的振动角频率。
12.如权利要求9所述的振动旋转传感器,其特征在于ωxr和ωxr,φxryr都等于0而Vmxr(t)和Vmyr(t)正比于分别取值为0和1的第一和第二时间方波函数,第一和第二方波函数不能同时为1。
13.如权利要求9所述的振动旋转传感器,其特征在于ωxr和ωxr等于预先确定的数值ωo而φxryr等于预定值ψo,且Vmxr(t)和Vmyr(t)正比于分别取值为0和1的第一和第二方波函数,ωo大于6ω,ω为共振器的振动角频率,第一和第二方波函数不能同时为1。
14.如权利要求9所述的振动旋转传感器,其特征在于ωxr、ωxr、φxr和φyr等于0而Vmxr(t)和Vmyr(t)正比于对伪随机序列取-1和1值的第一和第二方波函数。
15.如权利要求9所述的振动旋转传感器,其特征在于Vx1(t)和Vy2(t)还分别包括电压分量Vcx(t)Ux1(t)和Vcy(t)Uy2(t),量Vcx(t)和Vcy(t)是时间t的函数或者常数,电压分量Vcx(t)Ux1(t)和Vcy(t)Uy2(t)使作用力施加到共振器上。
16.如权利要求15所述的振动旋转传感器,其特征在于Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)为常数,ωxr、ωxr和|ωxrxr|大于6ω,ω为共振器的振动角频率。
17.如权利要求15所述的振动旋转传感器,其特征在于Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)为常数,ωxr和ωyr等于预先确定的值ωo,φxryr等于π/2圆弧,ωo大于6ω,ω为共振器的振动角频率。
18.如权利要求15所述的振动旋转传感器,其特征在于ωxr、ωxr、φxr和φyr等于0而Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)正比于分别取值为0和1的第一、第二、第三和第四时间方波函数,第一、第二、第三和第四方波函数不能同时为1。
19.如权利要求15所述的振动旋转传感器,其特征在于ωxr和ωxr等于预先确定的数值ωo而φxr和φyr等于预先确定的数值φo,Vmxr(t)、Vmyr(t)、Vcx(t)和Vcy(t)正比于分别取值为0和1的第一、第二、第三和第四时间方波函数,ωo大于6ω,ω为共振器的振动角频率,第一、第二、第三和第四方波函数不能同时为1。
20.如权利要求15所述的振动旋转传感器,其特征在于Vx1(t)包括电压分量Vcax(t)Uax1(t)、Vcqx(t)Uqx1(t)、和Vcrx(t)Urx1(t),而Vy2(t)包括电压分量Vcay(t)Uay1(t)、Vcqy(t)Uqy1(t)、和Vcry(t)Ury1(t),ωxr、ωxr、φxr和φyr等于0,Vmxr(t)、Vmyr(t)、Vcax(t)、Vcqx(t)、Vcrx(t)、Vcay(t)、Vcqy(t)、和Vcry(t)正比于分别取值为0和1的第一、第二、第三、第四、第五、第六、第七和第八时间方波函数,第一、第二、第三和第四、第五、第六、第七和第八方波函数不能同时为1。
21.如权利要求15所述的振动旋转传感器,其特征在于Vx1(t)包括电压分量Vcax(t)Uax1(t)、Vcqx(t)Uqx1(t)、和Vcrx(t)Urx1(t),而Vy2(t)包括电压分量Vcay(t)Uay1(t)、Vcqy(t)Uqy1(t)、和Vcry(t)Ury1(t),ωxr和ωxr等于预先确定的数值ωo而φxr和φyr等于预先确定的数值φo,ωo大于6ω,ω为共振器的振动角频率,Vmxr(t)、Vmyr(t)、Vcax(t)、Vcqx(t)、Vcrx(t)、Vcay(t)、Vcqy(t)、和Vcry(t)正比于分别取值为0和1的第一、第二、第三、第四、第五、第六、第七和第八时间方波函数,第一、第二、第三和第四、第五、第六、第七和第八方波函数不能同时为1。
22.如权利要求15所述的振动旋转传感器,其特征在于ωxr、ωxr、φxr和φyr等于0而Vmxr(t)和Vmyr(t)正比于伪随机序列取值-1和1的的第一和第二时间方波函数,Vcx(t)和Vcy(t)为常数。
23.如权利要求1所述的振动旋转传感器,其特征在于进一步包括:
附着有共振器的壳罩,壳罩包括多个靠近一个或多个共振器电极的附着电极;
向一个或多个壳罩电极提供控制电压的驱动电子线路,驱动电子线路产生跟踪角,驱动电子线路使得跟踪角等于常数或者是多个变量中的一个或多个的函数,多个变量包括时间、驻波方向和外部电源的变量。
24.如权利要求1所述的振动旋转传感器,其特征在于进一步包括:
附着有共振器的壳罩,壳罩包括多个靠近一个或多个共振器电极的附着电极;
向一个或多个壳罩电极提供控制电压的驱动电子线路,驱动电子线路产生跟踪角,驱动电子线路通过将跟踪角加上方位角与跟踪角之差值确定驻波的方位角。
25.一种控制和读取振动旋转传感器的方法,所述传感器包括具有与单个输出端口相连的一个或多个电极的共振器和靠近共振器电极包括多个电极的壳罩,共振器能够以一种或多种驻波模式振动,每种驻波模式由多个参量定义,驻波相对共振器上参考点的方向由方位角定义,所述方法的特征在于包括以下步骤:
(a)定义相对共振器上参考点的跟踪角;
(b)生成多个作为跟踪角函数的驱动电压;
(c)将驱动电压施加到多个壳罩电极的每一个上;
(d)通过对共振器输出端口上共振器信号进行的操作确定驻波方位角与跟踪角之间的差值。
26.如权利要求25所述的方法,其特征在于每个驱动电压包括第一激发电压和第二激发电压,第一和第二激发电压的频谱局限于分立的频带。
27.如权利要求25所述的方法,其特征在于不同驱动电压的个数为2个,每个驱动电压包括第一激发电压和第二激发电压,每个激发电压都是频率和相位预先确定的时间周期函数,第一和第二激发电压的频率相同,相位相差1/4周期。
28.如权利要求25所述的方法,其特征在于每个驱动电压包括第一激发电压和第二激发电压,每个第一和第二激发电压正比于取值为0和1的独特方波,在任何给定时刻只有一个方波取值为1。
29.如权利要求25所述的方法,其特征在于每个驱动电压包括第一激发电压和第二激发电压,每个第一和第二激发电压正比于具有预先确定频率和相位的周期函数与取值为0和1的独特方波的乘积,在任何给定时刻只有一个方波取值为1。
30.如权利要求25所述的方法,其特征在于每个驱动电压包括第一激发电压和第二激发电压,每个第一和第二激发电压正比于按照预先确定的伪随机序列取值为1和-1的独特方波。
31.如权利要求25所述的方法,其特征在于每个驱动电压包括第一激发电压和第二激发电压,第一激发电压包括与两倍跟踪角的余弦值相等的倍增因子,第二激发电压包括与两倍跟踪角正弦值相等的倍增因子。
32.如权利要求25所述的方法,其特征在于至少两个驱动电压中的每一个包括第一激发电压、第二激发电压和作用力电压,第一和第二激发电压的频谱局限于分立的频带,作用力电压的频谱局限于与第一和第二激发电压有关频带相分离的频带。
33.如权利要求25所述的方法,其特征在于不同驱动电压的个数为2个,每个驱动电压包括第一激发电压、第二激发电压和作用力电压,第一和第二激发电压中的每个都是频率和相位预先确定的周期函数,第一和第二激发电压的频率相同,相位相差1/4周期,作用力电压的频谱局限于与第一和第二激发电压频率分离的频带。
34.如权利要求25所述的方法,其特征在于至少两个驱动电压中的每一个包括第一激发电压、第二激发电压和作用力电压,每个第一和第二激发电压正比于取值为0和1的独特方波,每个作用力电压包括正比于取值为0和1的方波的倍增因子,在任意时刻只有一个与激发和作用力电压相关的方波取值为1。
35.如权利要求25所述的方法,其特征在于至少两个驱动电压中的每一个包括第一激发电压、第二激发电压和作用力电压,每个第一和第二激发电压正比于具有预先确定频率和相位的周期函数与取值为0和1的独特方波的乘积,每个作用力电压包括正比于取值为0和1的方波的倍增因子,在任意时刻只有一个与激发和作用力电压相关的方波取值为1。
36.如权利要求25所述的方法,其特征在于至少两个驱动电压中的每一个包括第一激发电压、第二激发电压和作用力电压,每个第一和第二激发电压正比于按照预先确定的伪随机序列取值为1和-1的独特方波,作用力电压的频谱局限于与激发电压有关频带相分离的频带。
37.如权利要求25所述的方法,其特征在于共振器信号是两个作为驻波参量和跟踪角函数的分量之和,步骤(d)包括以下步骤:
(d1)从共振器信号中提取两个分量中的每一个;
(d2)通过对两个分量的操作确定一个驻波方位角与跟踪角之差值。
38.如权利要求37所述的方法,其特征在于两个分量占据分立的频带,每个分量通过在区分分立频带的共振器信号上进行操作而提取。
39.如权利要求37所述的方法,其特征在于两个分量为频率相同而相位相差1/4周期的时间周期函数,每个分量通过在区分两分量相位的共振器信号上进行操作而提取。
40.如权利要求37所述的方法,其特征在于两个分量存在于不同时间周期内的共振器信号中,每个分量通过在区分不同时间周期的共振器信号上进行操作而提取。
41.如权利要求37所述的方法,其特征在于两个分量是0和1的伪随机序列,两个伪随机序列的交叉相关度为零,每个分量通过在区分两伪随机序列的共振器信号上进行操作而提取。
42.如权利要求37所述的方法,其特征在于两个分量各包括两项,一项包含方位角与跟踪角差值的正弦值,另一项包含方位角与跟踪角差值的余弦值。
43.如权利要求37所述的方法,其特征在于进一步包括以下步骤:
(e)使跟踪角等于常数或者一个或多个变量的函数,多个变量包括时间、驻波方向和外部电源的变量。
44.如权利要求43所述的方法,其特征在于进一步包括以下步骤:
(f)通过将跟踪角加到方位角与跟踪角之差值上计算方位角。
45.用于实现如权利要求25所述方法的装置。
46.用于实现如权利要求37所述方法的装置。
CN98105307.6A 1997-02-18 1998-02-18 全角跟踪振动旋转传感器 Pending CN1193731A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/802,009 US5801310A (en) 1997-02-18 1997-02-18 Vibratory rotation sensor with whole-angle tracking
US08/802,009 1997-02-18

Publications (1)

Publication Number Publication Date
CN1193731A true CN1193731A (zh) 1998-09-23

Family

ID=25182606

Family Applications (1)

Application Number Title Priority Date Filing Date
CN98105307.6A Pending CN1193731A (zh) 1997-02-18 1998-02-18 全角跟踪振动旋转传感器

Country Status (6)

Country Link
US (1) US5801310A (zh)
EP (1) EP0859217B1 (zh)
JP (1) JP4100752B2 (zh)
CN (1) CN1193731A (zh)
CA (1) CA2227155A1 (zh)
DE (1) DE69841731D1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6079270A (en) * 1997-07-08 2000-06-27 Litton Systems Inc. Method and apparatus for generating driving signals for a vibratory rotation sensor
US5983719A (en) * 1997-07-31 1999-11-16 Litton Systems, Inc. Low quantization method and apparatus for vibratory rotation sensors
US6065340A (en) * 1998-07-29 2000-05-23 Litton Systems Inc. Vibratory rotation sensor with novel electrode placements
US6158282A (en) * 1999-07-07 2000-12-12 Litton Systems Inc. Extraction of double-oscillation-frequency sinusoids in vibratory rotation sensors
US6189382B1 (en) * 1999-11-05 2001-02-20 Litton Systems, Inc. Vibratory sensor with self-calibration and low noise digital conversion
US6647785B2 (en) 2001-07-27 2003-11-18 Litton Systems, Inc. Nuclear radiation hard high accuracy rotation sensor system
US7296468B2 (en) * 2005-10-26 2007-11-20 Litton Systems, Inc. Digital coriolis gyroscope
US7251900B2 (en) * 2005-10-26 2007-08-07 Guy Thomas Varty Methods and systems utilizing intermediate frequencies to control multiple coriolis gyroscopes
US8661898B2 (en) 2008-10-14 2014-03-04 Watson Industries, Inc. Vibrating structural gyroscope with quadrature control
US9494425B2 (en) 2011-05-25 2016-11-15 The Regents Of The University Of Michigan Gyroscope and method of fabricating a resonator for a gyroscope
US9742373B2 (en) 2011-10-31 2017-08-22 The Regents Of The University Of Michigan Method of manufacturing a temperature-compensated micromechanical resonator
RU2704334C1 (ru) * 2019-04-02 2019-10-28 Федеральное государственное бюджетное учреждение науки Институт проблем механики им. А.Ю. Ишлинского Российской академии наук Способ считывания и управления колебаниями волнового твердотельного гироскопа
CN110398284B (zh) * 2019-07-19 2021-05-18 天津大学 一种新型弹性波超材料结构
CN114370886B (zh) * 2021-11-23 2024-01-02 上海航天控制技术研究所 基于虚拟转动的全角模式振动陀螺测量误差自标定方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4951508A (en) * 1983-10-31 1990-08-28 General Motors Corporation Vibratory rotation sensor
CA1250458A (en) * 1983-10-31 1989-02-28 Edward J. Loper, Jr. Hemispherical resonator gyro
JP3423418B2 (ja) * 1994-06-30 2003-07-07 松下電器産業株式会社 移動体旋回角度算出装置

Also Published As

Publication number Publication date
JP4100752B2 (ja) 2008-06-11
EP0859217B1 (en) 2010-06-23
EP0859217A3 (en) 2000-01-12
CA2227155A1 (en) 1998-08-18
DE69841731D1 (de) 2010-08-05
JPH10260057A (ja) 1998-09-29
US5801310A (en) 1998-09-01
EP0859217A2 (en) 1998-08-19

Similar Documents

Publication Publication Date Title
CN1200483A (zh) 具有多路复用电路的振动旋转传感器
CN1193731A (zh) 全角跟踪振动旋转传感器
US20100011857A1 (en) Arrangement for measuring a rate of rotation using a vibration sensor
JP6180551B2 (ja) 半球共振型ジャイロ
CN101846517B (zh) 杯形波动陀螺的杯形谐振子的机械平衡方法
EP3205978A1 (en) Inertial sensors
CN1174221C (zh) 振动陀螺仪
CN1894558A (zh) 借助转速科式陀螺测量转速/加速度的方法和实现该方法的科式陀螺
CN111595313B (zh) 一种基于离散电极的轴对称谐振陀螺参数激励方法
CN1193730A (zh) 施加交流作用力电压的振动旋转传感器
JP3585959B2 (ja) 加速度と角速度との双方を検出する装置
US6065340A (en) Vibratory rotation sensor with novel electrode placements
JP3002663B2 (ja) 振動性回転センサの振動モードを制御するための方法および装置
CN1201902A (zh) 利用线性弯曲量度法的振动旋转传感器
US5902930A (en) Method for reading out a vibratory rotation sensor
US6158282A (en) Extraction of double-oscillation-frequency sinusoids in vibratory rotation sensors
JP2010160095A (ja) 角速度センサ
JP2008151728A (ja) 2軸角速度センサ
JP2010223854A (ja) 3軸検出角速度センサ
CN103017748A (zh) 一种钟形振子式角速率陀螺信号提取方法
JPH05118854A (ja) 振動ジヤイロ
JP2004233369A (ja) 加速度と角速度との双方を検出する装置
JPH0545169A (ja) 圧電振動ジヤイロ

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C01 Deemed withdrawal of patent application (patent law 1993)
WD01 Invention patent application deemed withdrawn after publication