CN110135071A - Reliability detection method and system of single-axis fiber optic gyroscope based on multivariate performance degradation - Google Patents
Reliability detection method and system of single-axis fiber optic gyroscope based on multivariate performance degradation Download PDFInfo
- Publication number
- CN110135071A CN110135071A CN201910408484.6A CN201910408484A CN110135071A CN 110135071 A CN110135071 A CN 110135071A CN 201910408484 A CN201910408484 A CN 201910408484A CN 110135071 A CN110135071 A CN 110135071A
- Authority
- CN
- China
- Prior art keywords
- degradation
- reliability
- data
- zero bias
- fiber optic
- 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.)
- Granted
Links
- 238000006731 degradation reaction Methods 0.000 title claims abstract description 115
- 230000015556 catabolic process Effects 0.000 title claims abstract description 111
- 239000000835 fiber Substances 0.000 title claims abstract description 84
- 238000001514 detection method Methods 0.000 title claims abstract description 32
- 238000005315 distribution function Methods 0.000 claims abstract description 34
- 239000013307 optical fiber Substances 0.000 claims abstract 20
- 238000000034 method Methods 0.000 claims description 39
- 238000009792 diffusion process Methods 0.000 claims description 18
- 230000003068 static effect Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000007476 Maximum Likelihood Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 241000039077 Copula Species 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 3
- 230000007850 degeneration Effects 0.000 claims 5
- 230000003139 buffering effect Effects 0.000 claims 3
- 230000015572 biosynthetic process Effects 0.000 claims 2
- 238000003786 synthesis reaction Methods 0.000 claims 2
- 238000010276 construction Methods 0.000 claims 1
- 230000001360 synchronised effect Effects 0.000 claims 1
- 238000012360 testing method Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 230000005653 Brownian motion process Effects 0.000 description 2
- 238000005537 brownian motion Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Gyroscopes (AREA)
Abstract
本发明涉及一种基于多元性能退化的单轴光纤陀螺可靠性检测方法,包括:对光纤陀螺零偏退化量和标度因数退化量进行退化建模:构建光纤陀螺零偏性能可靠度模型;构建光纤陀螺标度因数可靠度函数模型;构造光纤陀螺的零偏和标度因数退化失效联合分布函数;计算得到光纤陀螺的可靠度综合模型,输出包括光纤陀螺可靠性信息的检测数据。本发明还公开了一种实施基于多元性能退化的单轴光纤陀螺可靠性检测方法的系统。本发明通过光纤陀螺可靠性检测系统获取所需试验数据,计算光纤陀螺可靠度综合函数,得到当前对应的被测光纤陀螺可靠性指标,直观评价当前被测光纤陀螺可靠性水平,为光纤陀螺的生产或设计,提供产品改进的技术支撑。
The invention relates to a single-axis optical fiber gyro reliability detection method based on multivariate performance degradation. The reliability function model of the FOG scale factor; construct the joint distribution function of the zero bias and the scale factor degradation failure of the FOG; calculate the reliability comprehensive model of the FOG, and output the detection data including the reliability information of the FOG. The invention also discloses a system for implementing the reliability detection method of the single-axis fiber optic gyroscope based on multivariate performance degradation. The present invention obtains the required test data through the reliability detection system of the fiber optic gyroscope, calculates the comprehensive function of the reliability of the fiber optic gyroscope, obtains the current corresponding reliability index of the fiber optic gyroscope under test, and intuitively evaluates the reliability level of the fiber optic gyroscope at present, which is the reliability level of the fiber optic gyroscope. Production or design, providing technical support for product improvement.
Description
技术领域technical field
本发明涉及可靠性分析及检测技术领域,尤其是基于多元性能退化的单轴光纤陀螺可靠性检测方法及系统。The invention relates to the technical field of reliability analysis and detection, in particular to a method and system for reliability detection of a single-axis fiber optic gyroscope based on multivariate performance degradation.
背景技术Background technique
随着我国工业生产的快速发展和装备的改进,对产品的质量和可靠性要求越来越高。光纤陀螺是惯性导航系统的关键部件,广泛应用于海陆空等军事、民用领域,具有长期贮存、一次使用、退化失效等特点。按结构分,光纤陀螺分为单轴和多轴光纤陀螺,其中单轴光纤陀螺应用最为广泛,因此对单轴光纤陀螺的性能检测、有效评估单轴光纤陀螺的可靠性水平显得尤为重要。如何通过有效的方法来检测单轴光纤陀螺的可靠性水平,不仅可以在要求进一步提高其性能水平时,指出改进设计的途径,而且在准确的了解了单轴光纤陀螺的性能水平后,可以有效的评定整个惯性系统工作的性能。With the rapid development of my country's industrial production and the improvement of equipment, the requirements for product quality and reliability are getting higher and higher. Fiber optic gyroscope is a key component of inertial navigation system, widely used in military and civil fields such as sea, land and air, and has the characteristics of long-term storage, one-time use, degradation and failure. According to the structure, fiber optic gyroscopes are divided into single-axis and multi-axis fiber-optic gyroscopes. Among them, single-axis fiber-optic gyroscopes are the most widely used. Therefore, it is particularly important to test the performance of single-axis fiber-optic gyroscopes and effectively evaluate the reliability level of single-axis fiber-optic gyroscopes. How to detect the reliability level of the single-axis fiber optic gyroscope by an effective method can not only point out the way to improve the design when it is required to further improve its performance level, but also can effectively understand the performance level of the single-axis fiber optic gyroscope. The evaluation of the performance of the entire inertial system operation.
目前针对单轴光纤陀螺的参数采集与检测系统,如基于RS232的检测系统,大多存在自动化程度低,硬件电路较复杂,测试数据的记录、处理分析等均由人工完成,耗时大,而且不能避免人为误差等;同时,在进行单轴光纤陀螺的可靠性分析时,大多针对单轴光纤陀螺的单一性能进行退化分析,而目前考虑多元性能退化进行分析的方法主要有两种:(1)假设多个性能退化过程之间相互独立或者服从多元正态分布,具有简单快速的优点,但分析不过全面,结果误差较大;(2)考虑多元性能参数之间的相关关系,如利用联合概率密度发和状态空间法,其优点是充分考虑了多元性能退化过程之间的相关性,但存在建模困难和计算量大的问题。At present, the parameter acquisition and detection systems for single-axis fiber optic gyroscopes, such as the detection systems based on RS232, mostly have a low degree of automation, complex hardware circuits, and the recording, processing and analysis of test data are all done manually, which is time-consuming and cannot be Avoid human error, etc. At the same time, when the reliability analysis of single-axis fiber optic gyroscope is carried out, most of the degradation analysis is carried out for the single performance of single-axis fiber optic gyroscope. At present, there are two main methods for analysis considering the degradation of multiple performances: (1) It is assumed that multiple performance degradation processes are independent of each other or obey the multivariate normal distribution, which has the advantage of being simple and fast, but the analysis is not comprehensive, and the result error is large; (2) Consider the correlation between multivariate performance parameters, such as using joint probability The advantage of density hair and state space method is that the correlation between multivariate performance degradation processes is fully considered, but there are problems of difficult modeling and large amount of calculation.
发明内容SUMMARY OF THE INVENTION
本发明的首要目的在于提供一种有效描述了多元性能退化量之间的相关关系,并且计算简单,外推性好的基于多元性能退化的单轴光纤陀螺可靠性检测方法。The primary purpose of the present invention is to provide a method for detecting the reliability of a single-axis fiber optic gyroscope based on the multivariate performance degradation, which effectively describes the correlation between the multivariate performance degradation quantities, is simple in calculation, and has good extrapolation.
为实现上述目的,本发明采用了以下技术方案:一种基于多元性能退化的单轴光纤陀螺可靠性检测方法,其特征在于:该方法包括以下顺序的步骤:In order to achieve the above object, the present invention adopts the following technical solutions: a method for detecting reliability of a single-axis fiber optic gyroscope based on multivariate performance degradation, characterized in that: the method comprises the steps in the following order:
(1)对光纤陀螺零偏退化量X1(t)和标度因数退化量X2(t)进行退化建模:(1) Degradation modeling is performed on the FOG bias degradation X 1 (t) and the scale factor degradation X 2 (t):
(2)构建光纤陀螺零偏性能可靠度模型;(2) Build the reliability model of the zero bias performance of the fiber optic gyroscope;
(3)构建光纤陀螺标度因数可靠度函数模型;(3) Build the reliability function model of the FOG scale factor;
(4)基于Frank Copula函数构造光纤陀螺的零偏和标度因数退化失效联合分布函数;(4) Based on the Frank Copula function, construct the joint distribution function of the zero bias and scale factor degradation failure of the fiber optic gyroscope;
(5)计算得到光纤陀螺的可靠度综合模型,由光纤陀螺可靠性综合模型生成并输出包括光纤陀螺可靠性信息的检测数据。(5) The reliability comprehensive model of the fiber optic gyroscope is obtained by calculation, and the detection data including the reliability information of the fiber optic gyroscope is generated and outputted from the reliability comprehensive model of the fiber optic gyroscope.
所述步骤(1)具体是指:The step (1) specifically refers to:
X1(t)=μ1t+σ1B(t)、X2(t)=μ2t+σ2B(t),式中,B(t)为标准布朗运动形式,μ1,σ1分别表示零偏退化率和零偏扩散系数,μ2,σ2分别表示标度因数退化率和标度因数扩散系数,由极大似然估计法得到μ1,σ1,μ2,σ2的估计值 X 1 (t)=μ 1 t+σ 1 B(t), X 2 (t)=μ 2 t+σ 2 B(t), where B(t) is the standard Brownian motion, μ 1 , σ 1 represents the zero bias degradation rate and bias diffusion coefficient, respectively, μ 2 , σ 2 represent the scale factor degradation rate and the scale factor diffusion coefficient, respectively . Estimated value of σ2
所述步骤(2)具体是指:The step (2) specifically refers to:
计算零偏退化失效分布函数其中分别为零偏退化率和零偏扩散系数的估计值,Φ(·)为标准正态分布函数,w1表示零偏退化失效阈值,构建光纤陀螺零偏性能可靠度模型 Calculate the zero-biased degradation failure distribution function in The estimated values of the zero-bias degradation rate and the zero-bias diffusion coefficient, respectively, Φ( ) is the standard normal distribution function, w 1 represents the zero-bias degradation failure threshold, and the reliability model of the zero-bias performance of the fiber optic gyroscope is constructed.
所述步骤(3)具体是指:Described step (3) specifically refers to:
计算标度因数退化失效分布函数其中分别为标度因数退化率和标度因数扩散系数的估计值,Φ(·)为标准正态分布函数,w2表示标度因数退化失效阈值,构建光纤陀螺标度因数可靠度函数模型Calculate the scale factor degradation failure distribution function in are the estimated values of the scale factor degradation rate and the scale factor diffusion coefficient, respectively, Φ( ) is the standard normal distribution function, w 2 represents the scale factor degradation failure threshold, and the FOG scale factor reliability function model is constructed
所述步骤(4)具体是指:Described step (4) specifically refers to:
其中,C(·)为具有边缘分布函数F1(t)、F2(t)的二元Copula函数,F1(t)为零偏退化量的失效分布函数,F2(t)为标度因数退化量的失效分布函数,为光纤陀螺零偏退化量和标度因数退化量之间的相关系数估计值,值由贝叶斯参数估计法得到; Among them, C(·) is a binary Copula function with marginal distribution functions F 1 (t) and F 2 (t), F 1 (t) is a failure distribution function of zero-biased degradation, and F 2 (t) is the standard The failure distribution function of the degree factor degradation amount, is the estimated value of the correlation coefficient between the offset degradation of the fiber optic gyroscope and the degradation of the scale factor, The value is obtained by the Bayesian parameter estimation method;
所述值由贝叶斯参数估计法得到具体是指:said The value obtained by the Bayesian parameter estimation method specifically refers to:
由零偏退化量的失效分布函数F1(t)、标度因数退化量的失效分布函数F2(t)计算出tj时刻函数值(F1(tj),F2(tj)),Copula函数的对数似然函数为进一步利用贝叶斯参数估计法得到相关系数估计值对应的贝叶斯公式为:P(α|(F1(tj),F2(tj)))∝π(α)·L(F1(tj),F2(tj)|α),式中P(α|(F1(tj),F2(tj)))为参数α的后验分布,π(α)为无信息先验分布。Calculate the function value (F 1 (t j ), F 2 (t j ) at time t j from the failure distribution function F 1 (t) of the zero bias degradation amount and the failure distribution function F 2 (t) of the scale factor degradation amount ), the log-likelihood function of the Copula function is Further use the Bayesian parameter estimation method to obtain the estimated value of the correlation coefficient The corresponding Bayesian formula is: P(α|(F 1 (t j ),F 2 (t j )))∝π(α)·L(F 1 (t j ),F 2 (t j )| α), where P(α|(F 1 (t j ), F 2 (t j ))) is the posterior distribution of the parameter α, and π(α) is the uninformative prior distribution.
所述步骤(5)具体是指:The step (5) specifically refers to:
式中R1(t)表示零偏性能退化量的可靠度函数,R2(t)表示标度因数性能退化量的可靠度函数,为零偏和标度因数退化失效联合分布函数;由光纤陀螺可靠性综合模型生成并输出包括光纤陀螺可靠性信息的检测数据。 where R 1 (t) represents the reliability function of the zero bias performance degradation, R 2 (t) represents the reliability function of the scale factor performance degradation, Joint distribution function of zero bias and scale factor degradation failure; generated from the comprehensive model of fiber optic gyroscope reliability and output detection data including fiber optic gyroscope reliability information.
所述由极大似然估计法得到μ1,σ1,μ2,σ2的估计值包括以下步骤:The obtaining of the estimated values of μ 1 , σ 1 , μ 2 , and σ 2 by the maximum likelihood estimation method includes the following steps:
(1a)光纤陀螺零偏退化数据形式为:(1a) The data form of FOG zero bias degradation is:
标度因数退化数据形式为:The scale factor degenerate data is of the form:
其中i=1,...,n,j=1,...,m,n为样本总数,m为测量时刻总数,初始时刻数据退化量为0,由公式Δx1,i(tj)=x1,i(tj)-x1,i(tj-1),Δx2,i(tj)=x2,i(tj)-x2,i(tj-1)分别计算得到相邻测量时刻间零偏退化增量:where i=1,...,n,j=1,...,m, n is the total number of samples, m is the total number of measurement moments, and the data degradation amount at the initial moment is 0. The formula Δx 1,i (t j ) =x 1,i (t j )-x 1,i (t j-1 ), Δx 2,i (t j )=x 2,i (t j )-x 2,i (t j-1 ), respectively Calculate the zero offset degradation increment between adjacent measurement times:
和标度因数退化增量: and scale factor degenerate increments:
(1b)根据极大似然估计法,由零偏退化增量ΔX1,i(tj)得到零偏退化率和零偏扩散系数估计值:由标度因数退化增量ΔX2,i(tj)得到标度因数退化率和标度因数扩散系数估计值:(1b) According to the maximum likelihood estimation method, the estimated values of the zero-bias degradation rate and the zero-bias diffusion coefficient are obtained from the zero-bias degradation increment ΔX 1,i (t j ): The scale factor degradation rate and scale factor diffusion coefficient estimates are obtained from the scale factor degradation increment ΔX 2,i (t j ):
本发明的另一目的在于提供一种实施基于多元性能退化的单轴光纤陀螺可靠性检测方法的系统,包括:Another object of the present invention is to provide a system for implementing a reliability detection method for a single-axis fiber optic gyroscope based on multiple performance degradation, including:
供电单元,保护光纤陀螺同时满足其工作所需+5V电压要求;The power supply unit protects the fiber optic gyroscope and meets the +5V voltage requirement for its work;
检测平台,提供光纤陀螺的检测环境,实现光纤陀螺的性能检测;The detection platform provides the detection environment of the fiber optic gyroscope and realizes the performance detection of the fiber optic gyroscope;
数据采集模块,实现光纤陀螺输出数据的采集并传送到计算机;The data acquisition module realizes the acquisition of the output data of the fiber optic gyroscope and transmits it to the computer;
计算机,对经数据采集模块传送过来的数据进行处理分析,完成光纤陀螺性能参数的标定,实现数据的显示和存储;The computer processes and analyzes the data transmitted by the data acquisition module, completes the calibration of the performance parameters of the fiber optic gyroscope, and realizes the display and storage of the data;
所述检测平台包括温控箱和单轴一体化转台,光纤陀螺放置于温控箱内;The detection platform includes a temperature control box and a single-axis integrated turntable, and the fiber optic gyroscope is placed in the temperature control box;
所述数据采集模块包括:The data acquisition module includes:
A/D转换器,将光纤陀螺输出的电压信号或电流信号转换为计算机能够识别的等效数字信号;A/D converter, which converts the voltage signal or current signal output by the fiber optic gyroscope into an equivalent digital signal that can be recognized by the computer;
数据缓冲模块,起到数据缓冲的作用,连接在A/D转换器和主机之间,A/D转换器产生的数据先存储在缓冲器中,当缓冲器中数据存满后由计算机取走数据;The data buffer module plays the role of data buffer and is connected between the A/D converter and the host. The data generated by the A/D converter is first stored in the buffer, and when the data in the buffer is full, it is taken away by the computer data;
控制电路,实现对光纤陀螺输出信号的转换以及计算机数据读取的控制;The control circuit realizes the conversion of the output signal of the fiber optic gyroscope and the control of the computer data reading;
USB接口电路,实现数据的传输,将采集数据传送到计算机。The USB interface circuit realizes data transmission and transmits the collected data to the computer.
所述控制电路包括:The control circuit includes:
USB控制器,控制数据的传输以及所有USB接口和接口上设备的正常运行;USB controller, which controls the transmission of data and the normal operation of all USB interfaces and devices on the interfaces;
逻辑控制电路,接收USB控制器发出的控制信号,产生三态缓冲电路和时钟电路的控制信号;The logic control circuit receives the control signal sent by the USB controller, and generates the control signal of the three-state buffer circuit and the clock circuit;
三态缓冲电路,为高阻态时实现光纤陀螺数据由AD9225芯片传输到静态RAM,为高、低电平态时实现静态RAM和USB控制器之间的数据传输;Tri-state buffer circuit, in the high-impedance state, the fiber optic gyroscope data is transmitted from the AD9225 chip to the static RAM, and in the high and low states, the data transmission between the static RAM and the USB controller is realized;
时钟电路,通过控制时序逻辑电路来产生AD9225芯片的时钟信号和转换开始信号、地址发生器的时钟信号以及调整A/D转换器与静态RAM之间的时钟同步,保证电路的时序正确;The clock circuit, by controlling the sequential logic circuit, generates the clock signal of the AD9225 chip and the conversion start signal, the clock signal of the address generator, and adjusts the clock synchronization between the A/D converter and the static RAM to ensure the correct timing of the circuit;
时序逻辑电路,由时钟电路控制产生AD9225芯片的时钟信号和转换开始信号、地址发生器的时钟信号;The sequential logic circuit is controlled by the clock circuit to generate the clock signal of the AD9225 chip, the conversion start signal, and the clock signal of the address generator;
地址发生器,由计数器实现,提供静态RAM的读写地址信号;The address generator, implemented by a counter, provides read and write address signals of static RAM;
AD9225芯片,提供12位精度、25Msps的高速模数转换;AD9225 chip, providing 12-bit precision, 25Msps high-speed analog-to-digital conversion;
静态RAM,为数据缓冲器,起到数据缓冲的作用。Static RAM, for data buffer, plays the role of data buffer.
由上述技术方案可知,本发明的优点在于:第一,本系统与基于RS232的检测系统相比,本系统由于采用了高速A/D芯片和USB接口,解决了传输速率慢、单位时间数据吞吐量低的问题,其数据量提高了十几个百分点;第二,同时实现了数据的的采集与存储以及基本处理功能,避免了以人工方式完成的繁琐,保证了效率;第三,本方法在进行可靠性分析时充分考虑了多元性能退化量之间的相关关系,并利用Frank Copula函数对光纤陀螺多元性能退化量进行有效耦合,与目前大多数使用的的联合概率密度法和状态空间法相比具有建模方便、计算简单的优点。As can be seen from the above technical solutions, the advantages of the present invention are: first, compared with the detection system based on RS232, the system adopts a high-speed A/D chip and a USB interface, which solves the problem of slow transmission rate and data throughput per unit time. The data volume is increased by more than ten percentage points; secondly, the collection and storage of data and basic processing functions are realized at the same time, which avoids the tediousness of manual completion and ensures the efficiency; thirdly, this method In the reliability analysis, the correlation between the multivariate performance degradation quantities is fully considered, and the Frank Copula function is used to effectively couple the multivariate performance degradation quantities of the fiber optic gyroscope. It has the advantages of convenient modeling and simple calculation.
附图说明Description of drawings
图1为本系统的组成框图;Fig. 1 is the composition block diagram of this system;
图2为控制电路的电路原理框图;Fig. 2 is the circuit principle block diagram of the control circuit;
图3为本发明的方法流程图;Fig. 3 is the method flow chart of the present invention;
图4为零偏退化数据曲线图;Figure 4 is a graph of zero bias degradation data;
图5为标度因数退化数据曲线图;Figure 5 is a graph of scale factor degradation data;
图6为零偏与标度因数退化量关系图;Figure 6 is a graph of the relationship between zero bias and scale factor degradation;
图7为贝叶斯参数估计思想流程图;Figure 7 is a flow chart of the idea of Bayesian parameter estimation;
图8为可靠度曲线图。Figure 8 is a reliability graph.
具体实施方式Detailed ways
如图3所示,一种基于多元性能退化的单轴光纤陀螺可靠性检测方法,该方法包括以下顺序的步骤:As shown in Figure 3, a method for reliability detection of single-axis fiber optic gyroscope based on multivariate performance degradation, the method includes the steps in the following sequence:
(1)对光纤陀螺零偏退化量X1(t)和标度因数退化量X2(t)进行退化建模:(1) Degradation modeling is performed on the FOG bias degradation X 1 (t) and the scale factor degradation X 2 (t):
(2)构建光纤陀螺零偏性能可靠度模型;(2) Build the reliability model of the zero bias performance of the fiber optic gyroscope;
(3)构建光纤陀螺标度因数可靠度函数模型;(3) Build the reliability function model of the FOG scale factor;
(4)基于Frank Copula函数构造光纤陀螺的零偏和标度因数退化失效联合分布函数;(4) Based on the Frank Copula function, construct the joint distribution function of the zero bias and scale factor degradation failure of the fiber optic gyroscope;
(5)计算得到光纤陀螺的可靠度综合模型,由光纤陀螺可靠性综合模型生成并输出包括光纤陀螺可靠性信息的检测数据。(5) The reliability comprehensive model of the fiber optic gyroscope is obtained by calculation, and the detection data including the reliability information of the fiber optic gyroscope is generated and outputted from the reliability comprehensive model of the fiber optic gyroscope.
所述步骤(1)具体是指:The step (1) specifically refers to:
X1(t)=μ1t+σ1B(t)、X2(t)=μ2t+σ2B(t),式中,B(t)为标准布朗运动形式,μ1,σ1分别表示零偏退化率和零偏扩散系数,μ2,σ2分别表示标度因数退化率和标度因数扩散系数,由极大似然估计法得到μ1,σ1,μ2,σ2的估计值 X 1 (t)=μ 1 t+σ 1 B(t), X 2 (t)=μ 2 t+σ 2 B(t), where B(t) is the standard Brownian motion, μ 1 , σ 1 represents the zero bias degradation rate and bias diffusion coefficient, respectively, μ 2 , σ 2 represent the scale factor degradation rate and the scale factor diffusion coefficient, respectively . Estimated value of σ2
所述步骤(2)具体是指:The step (2) specifically refers to:
计算零偏退化失效分布函数其中分别为零偏退化率和零偏扩散系数的估计值,Φ(·)为标准正态分布函数,w1表示零偏退化失效阈值,构建光纤陀螺零偏性能可靠度模型 Calculate the zero-biased degradation failure distribution function in The estimated values of the zero-bias degradation rate and the zero-bias diffusion coefficient, respectively, Φ( ) is the standard normal distribution function, w 1 represents the zero-bias degradation failure threshold, and the reliability model of the zero-bias performance of the fiber optic gyroscope is constructed.
所述步骤(3)具体是指:Described step (3) specifically refers to:
计算标度因数退化失效分布函数其中分别为标度因数退化率和标度因数扩散系数的估计值,Φ(·)为标准正态分布函数,w2表示标度因数退化失效阈值,构建光纤陀螺标度因数可靠度函数模型 Calculate the scale factor degradation failure distribution function in are the estimated values of the scale factor degradation rate and the scale factor diffusion coefficient, respectively, Φ( ) is the standard normal distribution function, w 2 represents the scale factor degradation failure threshold, and the FOG scale factor reliability function model is constructed
如图6所示,所述步骤(4)具体是指:As shown in Figure 6, the step (4) specifically refers to:
其中,C(·)为具有边缘分布函数F1(t)、F2(t)的二元Copula函数,F1(t)为零偏退化量的失效分布函数,F2(t)为标度因数退化量的失效分布函数,为光纤陀螺零偏退化量和标度因数退化量之间的相关系数估计值,值由贝叶斯参数估计法得到; Among them, C(·) is a binary Copula function with marginal distribution functions F 1 (t) and F 2 (t), F 1 (t) is a failure distribution function of zero-biased degradation, and F 2 (t) is the standard The failure distribution function of the degree factor degradation amount, is the estimated value of the correlation coefficient between the offset degradation of the fiber optic gyroscope and the degradation of the scale factor, The value is obtained by the Bayesian parameter estimation method;
如图7所示,所述值由贝叶斯参数估计法得到具体是指:As shown in Figure 7, the The value obtained by the Bayesian parameter estimation method specifically refers to:
由零偏退化量的失效分布函数F1(t)、标度因数退化量的失效分布函数F2(t)计算出tj时刻函数值(F1(tj),F2(tj)),Copula函数的对数似然函数为进一步利用贝叶斯参数估计法得到相关系数估计值对应的贝叶斯公式为:P(α|(F1(tj),F2(tj)))∝π(α)·L(F1(tj),F2(tj)|α),式中P(α|(F1(tj),F2(tj)))为参数α的后验分布,π(α)为无信息先验分布。Calculate the function value (F 1 (t j ), F 2 (t j ) at time t j from the failure distribution function F 1 (t) of the zero bias degradation amount and the failure distribution function F 2 (t) of the scale factor degradation amount ), the log-likelihood function of the Copula function is Further use the Bayesian parameter estimation method to obtain the estimated value of the correlation coefficient The corresponding Bayesian formula is: P(α|(F 1 (t j ),F 2 (t j )))∝π(α)·L(F 1 (t j ),F 2 (t j )| α), where P(α|(F 1 (t j ), F 2 (t j ))) is the posterior distribution of the parameter α, and π(α) is the uninformative prior distribution.
所述步骤(5)具体是指:The step (5) specifically refers to:
式中R1(t)表示零偏性能退化量的可靠度函数,R2(t)表示标度因数性能退化量的可靠度函数,为零偏和标度因数退化失效联合分布函数;由光纤陀螺可靠性综合模型生成并输出包括光纤陀螺可靠性信息的检测数据。 where R 1 (t) represents the reliability function of the zero bias performance degradation, R 2 (t) represents the reliability function of the scale factor performance degradation, Joint distribution function of zero bias and scale factor degradation failure; generated from the comprehensive model of fiber optic gyroscope reliability and output detection data including fiber optic gyroscope reliability information.
所述由极大似然估计法得到μ1,σ1,μ2,σ2的估计值包括以下步骤:The obtaining of the estimated values of μ 1 , σ 1 , μ 2 , and σ 2 by the maximum likelihood estimation method includes the following steps:
(1a)光纤陀螺零偏退化数据形式为:(1a) The data form of FOG zero bias degradation is:
标度因数退化数据形式为:The scale factor degenerate data is of the form:
其中i=1,...,n,j=1,...,m,n为样本总数,m为测量时刻总数,初始时刻数据退化量为0,由公式Δx1,i(tj)=x1,i(tj)-x1,i(tj-1),Δx2,i(tj)=x2,i(tj)-x2,i(tj-1)分别计算得到相邻测量时刻间零偏退化增量:和标度因数退化增量: where i=1,...,n,j=1,...,m, n is the total number of samples, m is the total number of measurement moments, and the data degradation amount at the initial moment is 0. The formula Δx 1,i (t j ) =x 1,i (t j )-x 1,i (t j-1 ), Δx 2,i (t j )=x 2,i (t j )-x 2,i (t j-1 ), respectively Calculate the zero offset degradation increment between adjacent measurement times: and scale factor degenerate increments:
(1b)根据极大似然估计法,由零偏退化增量ΔX1,i(tj)得到零偏退化率和零偏扩散系数估计值:由标度因数退化增量ΔX2,i(tj)得到标度因数退化率和标度因数扩散系数估计值:(1b) According to the maximum likelihood estimation method, the estimated values of the zero-bias degradation rate and the zero-bias diffusion coefficient are obtained from the zero-bias degradation increment ΔX 1,i (t j ): The scale factor degradation rate and scale factor diffusion coefficient estimates are obtained from the scale factor degradation increment ΔX 2,i (t j ):
如图1所示,本系统包括:As shown in Figure 1, the system includes:
供电单元,保护光纤陀螺同时满足其工作所需+5V电压要求;The power supply unit protects the fiber optic gyroscope and meets the +5V voltage requirement for its work;
检测平台,提供光纤陀螺的检测环境,实现光纤陀螺的性能检测;The detection platform provides the detection environment of the fiber optic gyroscope and realizes the performance detection of the fiber optic gyroscope;
数据采集模块,实现光纤陀螺输出数据的采集并传送到计算机;The data acquisition module realizes the acquisition of the output data of the fiber optic gyroscope and transmits it to the computer;
计算机,对经数据采集模块传送过来的数据进行处理分析,完成光纤陀螺性能参数的标定,实现数据的显示和存储;The computer processes and analyzes the data transmitted by the data acquisition module, completes the calibration of the performance parameters of the fiber optic gyroscope, and realizes the display and storage of the data;
所述检测平台包括温控箱和单轴一体化转台,光纤陀螺放置于温控箱内;The detection platform includes a temperature control box and a single-axis integrated turntable, and the fiber optic gyroscope is placed in the temperature control box;
所述数据采集模块包括:The data acquisition module includes:
A/D转换器,将光纤陀螺输出的电压信号或电流信号转换为计算机能够识别的等效数字信号;A/D converter, which converts the voltage signal or current signal output by the fiber optic gyroscope into an equivalent digital signal that can be recognized by the computer;
数据缓冲模块,起到数据缓冲的作用,连接在A/D转换器和主机之间,A/D转换器产生的数据先存储在缓冲器中,当缓冲器中数据存满后由计算机取走数据;The data buffer module plays the role of data buffer and is connected between the A/D converter and the host. The data generated by the A/D converter is first stored in the buffer, and when the data in the buffer is full, it is taken away by the computer data;
控制电路,实现对光纤陀螺输出信号的转换以及计算机数据读取的控制;The control circuit realizes the conversion of the output signal of the fiber optic gyroscope and the control of the computer data reading;
USB接口电路,实现数据的传输,将采集数据传送到计算机。The USB interface circuit realizes data transmission and transmits the collected data to the computer.
如图2所示,所述控制电路包括:As shown in Figure 2, the control circuit includes:
USB控制器,控制数据的传输以及所有USB接口和接口上设备的正常运行;USB controller, which controls the transmission of data and the normal operation of all USB interfaces and devices on the interfaces;
逻辑控制电路,接收USB控制器发出的控制信号,产生三态缓冲电路和时钟电路的控制信号;The logic control circuit receives the control signal sent by the USB controller, and generates the control signal of the three-state buffer circuit and the clock circuit;
三态缓冲电路,为高阻态时实现光纤陀螺数据由AD9225芯片传输到静态RAM,为高、低电平态时实现静态RAM和USB控制器之间的数据传输;Tri-state buffer circuit, in the high-impedance state, the fiber optic gyroscope data is transmitted from the AD9225 chip to the static RAM, and in the high and low states, the data transmission between the static RAM and the USB controller is realized;
时钟电路,通过控制时序逻辑电路来产生AD9225芯片的时钟信号和转换开始信号、地址发生器的时钟信号以及调整A/D转换器与静态RAM之间的时钟同步,保证电路的时序正确;The clock circuit, by controlling the sequential logic circuit, generates the clock signal of the AD9225 chip and the conversion start signal, the clock signal of the address generator, and adjusts the clock synchronization between the A/D converter and the static RAM to ensure the correct timing of the circuit;
时序逻辑电路,由时钟电路控制产生AD9225芯片的时钟信号和转换开始信号、地址发生器的时钟信号;The sequential logic circuit is controlled by the clock circuit to generate the clock signal of the AD9225 chip, the conversion start signal, and the clock signal of the address generator;
地址发生器,由计数器实现,提供静态RAM的读写地址信号;The address generator, implemented by a counter, provides read and write address signals of static RAM;
AD9225芯片,提供12位精度、25Msps的高速模数转换;AD9225 chip, providing 12-bit precision, 25Msps high-speed analog-to-digital conversion;
静态RAM,为数据缓冲器,起到数据缓冲的作用。Static RAM, for data buffer, plays the role of data buffer.
以下结合图1至8对本发明作进一步的说明。The present invention will be further described below with reference to FIGS. 1 to 8 .
如图1所示,光纤陀螺被置于温控箱中,通过定位夹进行固定,并保持光纤陀螺的输入基准轴IRA的方向垂直于单轴一体化转台的平台面,单轴一体化转台和温控箱为光纤陀螺提供试验环境,供电单元提供光纤陀螺工作所需的+5V的电压,光纤陀螺的输出信号经过A/D转换器进行转换,转换后的数据先存储在数据缓冲器中,当数据缓冲器的数据存满后,由计算机读取数据,控制电路控制光纤陀螺输出信号的转换以及计算机数据的读取,计算机采集到数据后进行分析处理,实现光纤陀螺零偏和标度因数的标定,并实现数据的绘图和显示。As shown in Figure 1, the fiber optic gyroscope is placed in the temperature control box, fixed by the positioning clip, and the direction of the input reference axis IRA of the fiber optic gyroscope is kept perpendicular to the platform surface of the single-axis integrated turntable. The single-axis integrated turntable and The temperature control box provides the test environment for the fiber optic gyroscope. The power supply unit provides the +5V voltage required for the operation of the fiber optic gyroscope. The output signal of the fiber optic gyroscope is converted by the A/D converter, and the converted data is first stored in the data buffer. When the data in the data buffer is full, the computer reads the data, and the control circuit controls the conversion of the output signal of the fiber optic gyro and the reading of the computer data. calibration, and realize the drawing and display of data.
选择3只同一批次的单轴光纤陀螺试验样本进行试验分析,温控箱温度控制60℃并稳定1h进行预热,然后接通光纤陀螺电源,以固定时间间隔168h测量输出值,共测量25次,光纤陀螺输出信号经计算机采集并分析处理,对其零偏和标度因数进行标定后得到其零偏和标度因数退化数据如表1、表2所示,并分别作出退化曲线如图4、图5所示,由图4、图5可看出3只光纤陀螺试验样本零偏和标度因数均有明显退化趋势。Three samples of the same batch of single-axis fiber optic gyroscopes were selected for test analysis. The temperature of the temperature control box was controlled at 60 °C and stabilized for 1 hour for preheating. Then, the fiber optic gyroscope power supply was turned on, and the output value was measured at a fixed time interval of 168 hours, with a total of 25 measurements. Second, the output signal of the fiber optic gyroscope is collected, analyzed and processed by the computer, and its zero offset and scale factor are calibrated to obtain its zero offset and scale factor degradation data as shown in Table 1 and Table 2, and the degradation curves are respectively drawn as shown in Fig. 4. As shown in Figure 5, it can be seen from Figure 4 and Figure 5 that the zero offset and scale factor of the three fiber optic gyroscope test samples have obvious degradation trends.
表1零偏退化数据(°\h)Table 1 Zero bias degradation data (°\h)
表2标度因数退化数据(10e-6)Table 2 Scale factor degradation data (10e-6)
由表1、表2数据,计算得到零偏退化增量、标度因数退化增量数据如表3、表4所示;From the data in Table 1 and Table 2, the data of zero offset degradation increment and scale factor degradation increment are calculated as shown in Table 3 and Table 4;
表3零偏退化增量Table 3 Zero offset degradation increment
表4标度因数退化增量Table 4 Scale factor degradation increments
光纤陀螺综合可靠度函数曲线如图8所示,其横坐标表示时间,单位为h,纵坐标为对应时间的可靠度,根据上述可靠度函数进行寿命分析,分析结果如表5所示:在31175h后,光纤陀螺的可靠度为0.97;在31941h后,光纤陀螺可靠度为0.95;在34693h后,光纤陀螺可靠度为0.8;在39939h后,光纤陀螺可靠度为0.3。The comprehensive reliability function curve of the fiber optic gyroscope is shown in Figure 8. The abscissa represents the time, and the unit is h, and the ordinate represents the reliability of the corresponding time. According to the above reliability function, the life analysis is carried out. After 31175h, the reliability of the FOG is 0.97; after 31941h, the reliability of the FOG is 0.95; after 34693h, the reliability of the FOG is 0.8; after 39939h, the reliability of the FOG is 0.3.
表5table 5
综上所述,本发明综合考虑零偏和标度因数性能对光纤陀螺的可靠性影响,通过光纤陀螺可靠性检测系统获取所需试验数据,计算光纤陀螺可靠度综合函数,得到当前对应的被测光纤陀螺可靠性指标,直观评价当前被测光纤陀螺可靠性水平,为光纤陀螺的生产或设计,提供产品改进的技术支撑。To sum up, the present invention comprehensively considers the influence of zero bias and scale factor performance on the reliability of the fiber optic gyroscope, obtains the required test data through the fiber optic gyroscope reliability detection system, calculates the comprehensive function of the reliability of the fiber optic gyroscope, and obtains the current corresponding Measure the reliability index of the fiber optic gyroscope, intuitively evaluate the reliability level of the current fiber optic gyroscope, and provide technical support for product improvement for the production or design of the fiber optic gyroscope.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910408484.6A CN110135071B (en) | 2019-05-16 | 2019-05-16 | Method and system for reliability detection of single-axis fiber optic gyroscope based on multivariate performance degradation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910408484.6A CN110135071B (en) | 2019-05-16 | 2019-05-16 | Method and system for reliability detection of single-axis fiber optic gyroscope based on multivariate performance degradation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110135071A true CN110135071A (en) | 2019-08-16 |
CN110135071B CN110135071B (en) | 2022-09-27 |
Family
ID=67574612
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910408484.6A Active CN110135071B (en) | 2019-05-16 | 2019-05-16 | Method and system for reliability detection of single-axis fiber optic gyroscope based on multivariate performance degradation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110135071B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111678535A (en) * | 2020-04-27 | 2020-09-18 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | Accelerated test method for reliability of fiber optic gyroscope |
CN113124899A (en) * | 2021-03-23 | 2021-07-16 | 西安航天精密机电研究所 | Method for acquiring variable-temperature scale factor of fiber optic gyroscope based on simulation technology |
CN115824264A (en) * | 2023-02-23 | 2023-03-21 | 中国船舶集团有限公司第七〇七研究所 | Method for evaluating and improving process reliability of hollow-core microstructure fiber-optic gyroscope |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014197139A1 (en) * | 2013-06-07 | 2014-12-11 | Raytheon Applied Signal Technology, Inc. | System and method for distribution free target detection in a dependent non-gaussian background |
CN109033710A (en) * | 2018-08-30 | 2018-12-18 | 电子科技大学 | A kind of momenttum wheel reliability estimation method based on more performance degradations |
-
2019
- 2019-05-16 CN CN201910408484.6A patent/CN110135071B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014197139A1 (en) * | 2013-06-07 | 2014-12-11 | Raytheon Applied Signal Technology, Inc. | System and method for distribution free target detection in a dependent non-gaussian background |
CN109033710A (en) * | 2018-08-30 | 2018-12-18 | 电子科技大学 | A kind of momenttum wheel reliability estimation method based on more performance degradations |
Non-Patent Citations (1)
Title |
---|
毛端海等: "基于多元性能退化量光纤陀螺贮存寿命综合评估", 《计算机测量与控制》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111678535A (en) * | 2020-04-27 | 2020-09-18 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | Accelerated test method for reliability of fiber optic gyroscope |
CN113124899A (en) * | 2021-03-23 | 2021-07-16 | 西安航天精密机电研究所 | Method for acquiring variable-temperature scale factor of fiber optic gyroscope based on simulation technology |
CN115824264A (en) * | 2023-02-23 | 2023-03-21 | 中国船舶集团有限公司第七〇七研究所 | Method for evaluating and improving process reliability of hollow-core microstructure fiber-optic gyroscope |
Also Published As
Publication number | Publication date |
---|---|
CN110135071B (en) | 2022-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110135071B (en) | Method and system for reliability detection of single-axis fiber optic gyroscope based on multivariate performance degradation | |
US8996928B2 (en) | Devices for indicating a physical layer error | |
CN103941579B (en) | A kind of moment for oceanographic instrumentation records and clock synchronizing method | |
US20230068927A1 (en) | Extended reality movement platform | |
CN104034351A (en) | Method for testing property of high-accuracy fiber-optic gyro on basis of dynamic condition of incremental method | |
CN116678403A (en) | Temperature compensation method, device, equipment and storage medium of inertial measurement device | |
CN102062589A (en) | Fiber-optic gyroscope based angular displacement measuring device and method | |
CN111024119B (en) | Rapid calibration method for triaxial MEMS gyroscope | |
Geist et al. | Gravity-based characterization of three-axis accelerometers in terms of intrinsic accelerometer parameters | |
CN101109637A (en) | A Method of Improving the Performance of Fiber Optic Gyroscope Bias Temperature Sensitivity | |
CN109406315B (en) | Multi-parameter cooperative measurement device for initiating explosive device impact response | |
CN103033196B (en) | Electronic star simulator for APS star sensor | |
Wedha et al. | Design and build Mini Digital Scale using internet of things | |
CN110726852A (en) | MEMS accelerometer temperature compensation method | |
CN102254569A (en) | Quad-data rate (QDR) controller and realization method thereof | |
CN201562037U (en) | A/D chip conversion error detection device | |
CN104748765B (en) | The X-ray pulsar navigation electronics read-out device of function is measured with energy level | |
US9100112B1 (en) | Latency built-in self-test | |
Untilov et al. | Results of fiber-optic gyro testing | |
CN105027037B (en) | User's input is carried out using portable electric appts | |
US9551640B2 (en) | Wire-pull test location identification on a wire of a microelectronic package | |
US12109017B2 (en) | Remotely tracking range of motion measurement | |
CN103867189A (en) | Distributed type well-temperature measuring device and measuring method | |
RU220141U1 (en) | Linear acceleration measuring device | |
US12096053B2 (en) | Content delivery platform |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |