Temperature control optical fiber inertial navigation structure and design method thereof

Abstract

The invention relates to the technical field of optical fiber gyroscopes, in particular to a temperature control optical fiber inertial navigation structure and a design method thereof.The design method comprises the following steps of detecting heat flow direction changes, extracting offset numbers, generating frequency density, analyzing the thickness contour and deformation coincidence rate, outputting correction gradients, matching jump trends and structural connection, and generating a priority sequence; and screening offset channel numbers, integrating compression and deformation trends, and outputting a design scheme. According to the method, by capturing heat flow direction change points and extracting offset area numbers, density expression of a heat energy propagation path is achieved, the space recognition capability of thermal disturbance is enhanced, the structural adaptability of thermal response is improved, and stable conduction of heat flow between structural levels is optimized through sequence integration of a temperature jump trend and a connection relation; and extracting a stable channel number in a lateral offset state, enhancing path directivity, integrating a multi-region compression trend and structural deformation mapping, and improving thermal situation identification capability and structural thermal adaptability.

Description

Technical Field

[0001] The present invention relates to the technical field of fiber optic gyroscopes, and in particular to a temperature-controlled fiber optic inertial navigation structure and a design method thereof. Background Art

[0002] The field of fiber optic gyroscope technology includes inertial sensing technology that uses the principle of interference to measure angular velocity or angle. The core content of this technology is based on the Sagnac effect of laser interference, which interferes and compares light waves propagating in opposite directions in the optical fiber coil to determine the angular motion state of the carrier. Its overall system includes components such as light source, beam splitter, optical fiber ring, and interference detector. Fiber optic gyroscopes are widely used in a variety of high-precision attitude measurement scenarios such as aerospace, ship navigation, unmanned driving, and precision guidance. They have the characteristics of no rotating parts, fast response speed, and long life. Under long-term use or extreme environments, temperature changes will cause changes in parameters such as the refractive index and length of the optical fiber material, which in turn affects the measurement accuracy of the gyroscope. Therefore, the temperature control mechanism is crucial in the design of fiber optic gyroscopes.

[0003] Among them, the temperature-controlled fiber-optic inertial navigation structure design method refers to a method for structurally optimizing the key components in the fiber-optic inertial navigation system by means of temperature control. The technical matters involved in the subject of this patent mainly cover the layout design of the temperature control structure, the distribution method of the heat source, the construction of the thermal insulation structure and its integration with the fiber-optic ring structure. It uses thermal conductivity matching materials to design a fixed support structure, and parametrically analyzes the structural layout by simulating the heat conduction path under different working temperatures. At the same time, a temperature stabilization element is introduced to realize active regulation of the heat source. The method also optimizes the heat conduction path and reduces the temperature gradient in the fiber-optic ring area by constructing a multi-layer thermal isolation wrapping system, thereby providing a more stable thermal environment foundation for the fiber-optic gyroscope at the structural level.

[0004] In the existing technology, the thermal control structure has a weak adaptability problem when responding to multi-point spatial thermal disturbances. It lacks a dynamic tracking mechanism for the evolution of the microscopic heat flow direction, making it difficult to identify the intrinsic relationship between local high-frequency reflection behavior and the thermal direction splitting area. In terms of structural surface parameter adjustment, the thickness and deformation characteristics fail to form an effective linkage, resulting in a response hysteresis phenomenon in the heat retention area, and unable to provide a rapid control method for temperature fluctuations in a highly dynamic environment. The establishment of the heat conduction path relies on the static structural layout, lacks the coupling analysis of multi-point jump behavior and connection relationship, and is difficult to deal with the problem of heat flow coordination between complex structural layers. In the lateral offset state, no effective path number and structural feature alignment model is established, resulting in the heat channel being prone to conduction path dispersion and insufficient directionality. In the compression process of the optical fiber winding structure, the deformation area trend is not modeled as a whole, and the compression trend is disconnected from the structural offset state, affecting the matching efficiency of the temperature control response strategy and the stability of the thermal field maintenance. Summary of the Invention

[0005] The purpose of the present invention is to solve the shortcomings of the prior art and to propose a temperature-controlled fiber optic inertial navigation structure and a design method thereof.

[0006] In order to achieve the above object, the present invention adopts the following technical solution: a method for designing a temperature-controlled fiber-optic inertial navigation structure, comprising the following steps: S1: Detect the point where the heat flow direction changes at the fiber microcavity boundary, analyze the reflection direction offset trajectory and reverse expansion position, compare the reflection number distribution and the difference in the intersection area, integrate the offset area number, and generate the thermal direction splitting frequency density; S2: Extracting the thickness profile of the outer heat reflective material based on the thermal direction splitting frequency density, analyzing the consistency of the thickness fluctuation direction and the overlap rate of the retention surface deformation, identifying the mutation trend site and calculating the relationship between thickness and offset, and outputting the thermal buffer surface correction gradient; S3: calling the mutation site of the thermal buffer surface correction gradient, matching the temperature jump trend with the structural layer overlap range, screening the jump area connection and sorting the heat adsorption duration difference, and generating a heat conduction path priority sequence; S4: detecting the lateral offset morphology of the conduction direction according to the heat conduction path priority sequence, identifying the boundary area contour expansion direction and area deformation parameters, screening the continuous offset channel structure, and generating a heat channel expansion path number group; S5: Based on the thermal channel expansion path number group, the compression trend of the winding layer is compared with the degree of thermal deformation alignment, the consistency between the compression direction and the structural deformation is determined, the full structure offset trend data is integrated, and the temperature-controlled fiber optic inertial navigation design plan is output.

[0007] As a further solution of the present invention, the thermal direction splitting frequency density includes the reflection angle change sequence, the direction offset frequency distribution, and the heat flow reflection intersection area number; the thermal buffer surface correction gradient includes the material thickness difference distribution, the contour deformation offset parameter, and the abnormal mutation trend identification; the heat conduction path priority sequence includes the temperature jump area index, the structural connection relationship sequence, and the thermal adsorption continuity classification; the thermal channel expansion path number group includes the lateral offset structure number, the contour overlap deviation feature, and the area expansion channel identification; the temperature-controlled fiber optic inertial navigation design scheme includes the winding spacing compression mode, the thermal deformation alignment parameter, and the structural deformation trend distribution map.

[0008] As a further solution of the present invention, the specific steps of S1 are: S101: Obtain the heat flow direction change points at the microcavity boundary of the outer periphery of the fiber winding structure, select the boundary coordinates whose angle change rate exceeds the deflection threshold, and generate a heat flow direction mutation coordinate set; S102: Based on the heat flow direction mutation coordinate set, call the reflection direction parameter and the heat flow duration, record the reflection direction change sequence, extract the direction offset trajectory, count the rotation frequency and intersection position, and obtain the reflection trajectory intersection density distribution result; S103: According to the reflection trajectory cross density distribution result, determine the region number where the offset frequency is higher than the set threshold, extract the spatial distribution, and calculate the unit offset frequency density difference index to obtain the thermal direction splitting frequency density.

[0009] As a further solution of the present invention, the unit offset frequency density difference index calculation formula is specifically: ; in, Representative area The number is The unit offset frequency density difference index, Representative area Middle The reflection traces are numbered The frequency offset of the unit, Representative area Middle The reflection traces are numbered The effective action path distance of the unit, Indicates the The reflection traces are numbered The unit's horizontal reflection angle, Indicates the The reflection traces are numbered The vertical reflection angle of the unit, Indicates area The frequency offset threshold is set. Indicates the number The unit in the area The total number of observed reflection tracks in .

[0010] As a further solution of the present invention, the specific steps of S2 are: S201: Acquire position data of a selected area in the thermal direction splitting frequency density, extract the thickness value of the outer heat reflective material in the transverse coordinate, establish a mapping relationship between the coordinate and the thickness, determine the direction of thickness change, and obtain a thickness change trajectory value of the heat reflective layer; S202: Based on the thickness change trajectory value of the heat reflective layer, extracting the boundary coordinates of the heat retention surface contour, determining the degree of overlap between the thickness change direction and the angle between the boundary direction, counting the number of deformation point distributions, and obtaining the thickness direction consistency overlap; S203: According to the thickness direction consistency overlap, the overlapping area numbers are screened, the thickness difference and the contour offset are extracted, the thickness offset derivative is calculated and a variation function is established to obtain the thermal buffer surface correction gradient.

[0011] As a further solution of the present invention, the specific steps of S3 are: S301: Calling the mutation trend distribution position in the thermal buffer surface correction gradient, extracting the thermal retention change trajectory corresponding to the position, calculating the difference interval between the node temperature value and the average thermal retention level in the time series, establishing a mapping relationship between the coordinate sequence and the retention change, and obtaining the thermal retention trajectory offset; S302: Based on the thermal retention trajectory offset, determine the spatial region boundary corresponding to the temperature jump trend, extract the jump region and the adjacent hierarchical structure numbers, compare the connection mode of the structure number boundary with the jump coordinate intersection rate, filter the overlapping structure segments, and obtain the structure jump connection relationship group; S303: Based on the structure jump connection relationship group, the thermal adsorption duration of the structure segment is counted and the average adsorption difference is calculated. All structure combinations are sorted in ascending order of the difference, the order number of the jump node in the sorting chain is recorded, and adjacent jump structures are connected and the number sequence is combined to obtain the heat conduction path priority sequence.

[0012] As a further solution of the present invention, the specific steps of S4 are: S401: Based on the combination position in the heat conduction path priority sequence, detecting the lateral offset change of the heat conduction direction in the corresponding structure, extracting the contour extension coordinates of the path intersection area, determining the angle difference between the extension direction and the original path direction, and obtaining the lateral offset contour direction value; S402: Identifying the area deformation position of the thermal extension contour in the path connection boundary based on the lateral offset contour direction value, calculating the area change of the contour boundary in the projection plane, and comparing it with the lateral torsion rate of the heat conduction path, setting an offset overlap threshold, and filtering the numbers with overlap rates greater than the threshold to obtain an area deformation overlap sequence; S403: Call the area deformation superposition sequence, filter the path numbers in the continuous offset state, extract the continuous distribution segments of the area change and the direction offset value in the number structure, determine the stability trend of the hot channel expansion and record the corresponding structure number, and obtain the hot channel expansion path number group.

[0013] As a further solution of the present invention, the calculation formula for the area change of the contour boundary in the projection plane is specifically: ; in, Represents the area change of the contour boundary in the projection plane, The radian form represents the lateral offset profile direction value. represents the transverse torsion rate of the heat conduction path, Represents the area of the original contour boundary in the projection plane, Represents the contour boundary projection area after thermal expansion deformation, Represents the total number of segments after the contour boundary is discretized, The index number representing the discretization segment, represents the local extension strain rate of the dth boundary point along the x-axis, / Represents the ductile strain gradient at the boundary of segment d.

[0014] As a further solution of the present invention, the specific steps of S5 are: S501: Extracting the winding layer spacing value and thermal deformation offset distribution curve of the corresponding structure based on the numbered area in the hot channel expansion path number group, calculating the spacing compression ratio and offset direction angle at the numbered position, determining whether the change trends of the two are consistent, and obtaining the structural deformation alignment consistency; S502: Based on the structural deformation alignment consistency, the regions with the same compression direction under consecutive numbering are screened, the number of superimposed numbers in each region is counted, the boundary position of the structural contour after thermal action in the numbering range is extracted, and the offset trajectory vector is calculated to obtain a thermal-compression deformation superposition trend group; S503: Call the thermal compression deformation superposition trend group, integrate the superposition of regional compression direction and offset trend, establish a mapping network between numbered positions and offset vectors, identify regional sequences with stable compression directions and form structural connection combinations, and obtain a temperature-controlled fiber-optic inertial navigation design solution.

[0015] Compared with the prior art, the advantages and positive effects of the present invention are: In the present invention, by capturing the change points of heat flow direction and extracting the offset area numbers, the density expression of the heat energy propagation path is realized, the spatial recognition ability of thermal disturbance is enhanced, and the structural adaptability of thermal response is improved. By integrating the sequence of temperature jump trends and connection relationships, the stable conduction of heat flow between structural levels is optimized, the stable channel numbers under the lateral offset state are extracted, the path directionality is enhanced, and the multi-region compression trend and structural deformation mapping are integrated to improve the thermal situation recognition ability and structural thermal adaptability. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0017] Figure 1 It is a schematic diagram of the steps of the present invention. DETAILED DESCRIPTION

[0018] The technical solution of the present invention is described below in conjunction with the accompanying drawings.

[0019] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to indicate examples, illustrations, or explanations. Any embodiment or design described as an "exemplary" in the present invention should not be interpreted as being preferred or advantageous over other embodiments or designs. Rather, the use of the word "exemplary" is intended to present concepts in a concrete manner. Furthermore, in the embodiments of the present invention, "and / or" can mean both or either of the two.

[0020] In the embodiments of the present invention, the terms "image" and "picture" may sometimes be used interchangeably. It should be noted that, when the distinction is not emphasized, the meanings they convey are the same. The terms "of," "corresponding," and "corresponding" may sometimes be used interchangeably. It should be noted that, when the distinction is not emphasized, the meanings they convey are the same.

[0021] In the embodiments of the present invention, sometimes a subscript such as W1 may be written as a non-subscript such as W1. When the difference is not emphasized, the meanings to be expressed are the same.

[0022] In order to make the technical problems, technical solutions and advantages to be solved by the present invention clearer, a detailed description will be given below with reference to the accompanying drawings and specific embodiments.

[0023] See also Figure 1 A method for designing a temperature-controlled fiber-optic inertial navigation structure comprises the following steps: S1: Detect the heat flow direction change points at the microcavity boundary around the fiber-wound structure, analyze the offset trajectory of the boundary reflection direction under the action of continuous heat flow, record the reflection angle change trajectory and the reverse expansion position, compare the difference between the distribution of the number of heat flow reflections and the direction intersection area, determine the area with high direction deviation frequency and integrate the number, select the part with high intersection number as the core distribution area, and generate the thermal direction splitting frequency density; S2: Based on the position data of the selected area in the thermal direction splitting frequency density, the thickness profile of the outer heat reflective material at that location is extracted. The directional consistency of the thickness change trajectory and the degree of contour deformation of the heat retention surface are determined. The thickness fluctuation amplitude and the contour deformation coincidence rate are compared to identify the locations with abnormal mutation trends in the covering material. The relationship between the regional thickness difference and the contour change offset value is extracted to obtain the thermal buffer surface correction gradient. S3: Call the distribution position of the mutation trend in the thermal buffer surface correction gradient, extract the heat retention change trajectory of the corresponding S, compare the overlap range between the temperature jump trend and the structural layer distribution in the region, screen the structural connection relationship between the jump trend area and the adjacent S, judge the difference in thermal adsorption duration and summarize and sort it, connect the sequential combination of multiple jump structures, and generate the priority sequence of heat conduction paths; S4: Based on the combined positions in the heat conduction path priority sequence, detect the lateral offset morphology of the heat conduction direction in the structure, record the contour expansion direction of the connecting boundary area, identify the area deformation contours of the thermal extension positions of the two boundary segments, compare the degree of contour offset overlap and the lateral torsion trend of the path, screen the continuous offset states in the area expansion path, extract the structure number that forms a stable channel expansion morphology, and obtain the heat channel expansion path number group;

[0024] S5: Based on the numbered areas in the thermal channel expansion path numbering group, compare the compression trend of the winding layer spacing with the degree of alignment of the thermal deformation offset distribution, determine whether the compression amplitude direction of the winding spacing is consistent with the structural deformation trend, identify the number of superimposed continuous areas in the compression direction, summarize the actual deformation contour offset trend of the structure after heating in the same area, integrate the distribution of the compression trend within the entire structure, and obtain a temperature-controlled fiber-optic inertial navigation design plan.

[0025] The thermal direction splitting frequency density includes the reflection angle change sequence, direction offset frequency distribution, and heat flow reflection intersection area numbering; the thermal buffer surface correction gradient includes the material thickness difference distribution, contour deformation offset parameters, and abnormal mutation trend identification; the heat conduction path priority sequence includes the temperature jump area index, structural connection relationship sequence, and thermal adsorption continuity classification; the thermal channel expansion path numbering group includes the lateral offset structure number, contour superposition deviation characteristics, and area expansion channel identification; the temperature-controlled fiber optic inertial navigation design scheme includes the winding spacing compression mode, thermal deformation alignment parameters, and structural deformation trend distribution map.

[0026] The specific steps of S1 are: S101: Obtain the heat flow direction change points at the microcavity boundary of the outer periphery of the fiber winding structure, select the boundary coordinates whose angle change rate exceeds the deflection threshold, and generate a heat flow direction mutation coordinate set; When obtaining the heat flux direction change point at the microcavity boundary of the fiber-wound structure, first, heat flux sensing nodes need to be evenly distributed around the fiber, with a spacing of 2 mm between each node and a sampling period of 0.1 seconds. The heat flux value per unit time is collected one by one and the corresponding direction vector is recorded. Then, the heat flux direction difference between the sampling moments of two adjacent nodes is compared point by point. By comparing the degree of angular change in the direction of the heat flux vector, it is determined whether a mutation occurs. The threshold for the angle change to be a mutation is set to 15 degrees. This threshold is obtained by extracting the boundary heat flux angle difference from 20 groups of samples with different winding densities and counting their frequency distribution. The upper limit of the angle difference distribution within the 95% confidence interval is selected as the deflection The threshold is used, so 15 degrees is used as the judgment benchmark. During the execution process, the difference in heat flow direction recorded at adjacent time points for any two adjacent nodes is calculated. For example, if the direction of the first node is 50 degrees and the direction of the second node is 70 degrees, the angle change is 20 degrees, which exceeds the threshold and is determined to be a mutation point. The judgment operation is repeated for all nodes, and finally all boundary points with angle changes greater than or equal to 15 degrees are screened out to form a mutation coordinate set, such as (12.5, 8.4), (15.2, 9.0), (18.0, 8.9), etc. All coordinates are recorded in millimeters, forming a set of heat flow direction mutation points that can be used for subsequent directional reflection analysis.

[0027] S102: Based on the heat flow direction mutation coordinate set, call the reflection direction parameters and heat flow duration, record the reflection direction change sequence, extract the direction offset trajectory, count the rotation frequency and intersection position, and obtain the reflection trajectory intersection density distribution result; Based on the mutation coordinate set, the continuous change of the subsequent heat flow direction is recorded at each coordinate point, and the direction vector change is recorded every 0.1 second. The total heat flow reflection duration is calculated within this time range. For example, if the heat flow at a certain point continuously deviates within 0.5 seconds, the direction change sequence within this period is recorded, and the amplitude of each direction change is calculated. For example, the initial direction gradually changes from 45 degrees to 90 degrees, 135 degrees, and then to 120 degrees, etc., and the sequence 45-90-135-120 is recorded. The frequency of reverse rotation can be obtained by counting the number of changes exceeding 90 degrees. In the above sequence, 90 to 135 and 135 to 120 are both considered The direction is reversed, and the rotation frequency is counted as 2 times. In addition, in order to obtain the intersection position of the heat flow trajectory, a fixed-interval grid partitioning method is used in the heat flow path recording process. The spatial distribution area is divided into 5mm×5mm cells, and the number of times the heat flow trajectory passes through each cell is recorded. If a certain area is crossed by three trajectories in different directions more than 2 times within 0.5 seconds, it is counted as a crossing event. This process is repeated to count the number of crossings in all areas, and finally the crossing density distribution data is formed, which indicates the density of the heat flow path crossing in each cell. For example, there are 4 crossings in the (10, 10) area and 1 in the (15, 10) area, forming a spatial density comparison result.

[0028] S103: Based on the reflection trajectory cross density distribution result, determine the region number where the offset frequency is higher than the set threshold, extract the spatial distribution, and calculate the unit offset frequency density difference index to obtain the thermal direction splitting frequency density; The calculation formula of the unit offset frequency density difference index is as follows: ; in, Representative area The number is The unit offset frequency density difference index, Representative area Middle The reflection traces are numbered The frequency offset of the unit, Representative area Middle The reflection traces are numbered The effective action path distance of the unit, Indicates the The reflection traces are numbered The unit's horizontal reflection angle, Indicates the The reflection traces are numbered The vertical reflection angle of the unit, Indicates area The frequency offset threshold is set. Indicates the number The unit in the area The total number of reflection tracks observed in ; In this formula, Representative area The number is The calculation process of the unit offset frequency density difference index involves the weighted summation and normalization operations of multiple physical quantities. Each parameter involved is directly related to the physical characteristics of the reflection trajectory, such as frequency, path, and reflection angle.

[0029] The process of obtaining each parameter is as follows: :Representative area Middle The reflection traces are numbered The frequency offset of the unit.

[0030] This value is obtained by frequency domain analysis of the reflected signal, and the frequency offset is extracted by using signal spectrum analysis tools such as Fourier transform. Usually, for a reflected signal in a certain frequency band, the frequency offset is calculated by the frequency difference with the reference signal. Assume that in actual monitoring, The frequency offset of the reflection trace is (Frequency change measured based on the actual reflected signal).

[0031] :Representative area Middle The reflection traces are numbered The effective action path distance of the unit.

[0032] The path distance is calculated by measuring the propagation time of the reflected signal and combining it with the known propagation speed. The effective path distance of the reflection trajectory is .

[0033] :Indicates the The reflection traces are numbered The unit's horizontal reflection angle.

[0034] The horizontal reflection angle is usually obtained by geometric analysis of the reflection path. Specifically, it can be calculated by the angle between the radar reflection direction and the observation direction. Assuming that the angle between the reflection path and the parallel line of the ground is measured, it can be obtained. .

[0035] :Indicates the The reflection traces are numbered The vertical reflection angle of the unit.

[0036] The vertical reflection angle is also obtained through geometric analysis, similar to the measurement of the horizontal reflection angle. The reflection angle is the angle between the reflected wave and the horizontal plane. Assume that in the actual calculation, we get .

[0037] : Indicates area Set the frequency offset threshold.

[0038] This value is set by experiments or environmental standards and represents the frequency deviation threshold of the area. The setting of this parameter is usually determined by system design or environmental conditions. Assume that the setting value is .

[0039] :Indicates the number is The unit in the area The total number of observed reflection tracks in .

[0040] This value is obtained by counting all the reflection traces in the area. Assume that in a certain measurement, there are A reflection track was observed.

[0041] According to the values of the above parameters, the formula calculation process is as follows: Calculate the contribution of each reflection trajectory For each reflection trajectory , first calculate the product of its frequency offset and path distance, and then divide it by the composite value of the horizontal and vertical reflection angles (calculated using the Pythagorean theorem): ; Calculate the offset frequency density difference index Next, the difference between each reflection trace and the set frequency offset threshold is calculated, and all reflection traces are summed up. Perform averaging: ; Then divide the sum by , and get the offset frequency density index: ; The results show that the region The number is The offset frequency density index of the unit is 27.4, which reflects the average difference between the frequency offset of all reflection tracks in the unit and the set threshold.

[0042] The specific steps of S2 are: S201: Acquire the position data of the selected area in the thermal direction splitting frequency density, extract the thickness value of the outer heat reflective material in the transverse coordinate, establish a mapping relationship between the coordinate and the thickness, determine the thickness change direction, and obtain the thickness change trajectory value of the heat reflective layer; After obtaining the position data of the selected area in the thermal direction splitting frequency density, the material thickness information of the outer layer heat reflective material within the range is first extracted point by point based on the position coordinate range of the area in the heat reflective layer structure diagram. A sampling point is set every 0.5 mm to form an initial data set of coordinate-thickness value pairs. For example, the horizontal coordinate ranges from x=0mm to x=20mm, and the thickness value t(x) is recorded at 41 points x=0, 0.5, 1.0...20.0 in sequence. Its value is obtained by identifying the material boundary pixels in the cross-sectional scanning image of the reflective layer. If the thickness at x=0mm is 0.32mm, at x=5.0mm is 0.35mm, and at x=10mm is 0.42mm, it is included in the thickness value sequence. At the same time, the corresponding relationship between the horizontal coordinate and the thickness is established. The continuous thickness data is judged for directional changes. The thickness difference between two adjacent points t(i+1)-t(i) is analyzed for positive and negative signs. If the difference is greater than 0.02mm, it is considered that the thickness has increased. If it is less than -0.02mm, it is considered that the thickness decreases. If there is a point in the adjacent coordinate points where the thickness changes from increasing to decreasing, or from decreasing to increasing, it is considered that a direction change has occurred. For example, from x=6.5mm to x=7.0mm, the thickness increases from 0.38mm to 0.40mm, and then decreases to 0.36mm at x=7.5mm. This is a direction reversal point. The direction of thickness change is marked as positive or negative and each continuous interval is recorded. Finally, all the horizontal coordinates corresponding to the thickness change are connected, and the sequence trajectory of the thickness change with the horizontal coordinate is plotted to obtain the trajectory value of the heat reflective layer thickness change.

[0043] S202: Based on the thickness change trajectory value of the heat reflection layer, the boundary coordinates of the heat retention surface contour are extracted, the degree of overlap between the thickness change direction and the angle between the boundary direction is determined, the number of deformation point distributions is counted, and the thickness direction consistency overlap is obtained; Based on the obtained thickness change trajectory value, the boundary direction is calculated for each direction change interval in the selected heat reflection area. The spatial coordinates of the heat retention surface contour boundary are extracted in combination with the actual material structure diagram. A boundary curve segment is generated for every 10 sampling points to construct a complete boundary contour segment set. At the same time, the connection direction vectors of the starting and ending points of each segment are recorded. The angle between the thickness change direction vector and the boundary direction vector in each segment is compared. The angle is converted into angle data by calculating the cosine value of the angle between the two vector endpoints and judged by the set coincidence angle threshold of 15°. This threshold refers to the limit of the influence of thickness change on contour deformation in the previous experiment. The matching angle setting shows that when the angle is less than or equal to 15°, it is considered as directional overlap. For example, in segment A, the thickness direction changes horizontally to the right, and the boundary direction is also to the upper right. The angle is measured to be 12°, which is considered overlap. The point is recorded as the overlap point. If there are 4 sampling points in the segment that meet the overlap condition, they are counted as 4 deformation points. The number of overlap points in all segments is counted, and finally the ratio of the total number of overlap points to the number of all sampling points is calculated. For example, there are 60 sampling points in the area, of which 45 are overlap points. The thickness direction consistency overlap is 0.75, indicating that the thickness change direction and the boundary direction at 75% of the points in the area have directional overlap.

[0044] S203: Based on the consistency and overlap in the thickness direction, the overlapping area numbers are screened, the thickness difference and the contour offset are extracted, the thickness offset derivative is calculated, and a variation function is established to obtain the thermal buffer surface correction gradient; According to the statistically obtained thickness direction consistency overlap, each area is judged to see whether it reaches the set overlap threshold. The overlap threshold is set to 0.65. The setting basis is the minimum overlap ratio of the deformation stable area in multiple groups of thermal material structure samples. If the consistency overlap of a region is equal to or higher than 0.65, the region number is listed as the overlap region number set. For example, the overlaps of Zone_05 and Zone_07 are 0.72 and 0.68 respectively, which meet the conditions and are included in the set. Then, the thickness difference sequence corresponding to the horizontal coordinates of Zone_05 and Zone_07 is extracted respectively. The thickness difference refers to the change in the thickness value of the adjacent sampling point. For example, t(x=5.0mm)= 0.36mm, t(x=5.5mm)=0.40mm, then the thickness difference is 0.04mm, and the contour offset is the change value of the longitudinal coordinate of the actual boundary position. If the boundary longitudinal coordinate is 2.5mm when it is 5.0mm and 2.7mm at 5.5mm, then the offset is 0.2mm. The derivative of the thickness difference with respect to the offset is calculated for this set of data. With a sampling interval of 0.5mm as the step size, the difference ratio of adjacent data is calculated and a data change sequence is established. The derivative sequence operation is repeated for all areas and normalized. Finally, the derivative sequence is integrated as a function variable to construct a thermal buffer surface correction gradient, which represents the degree of influence of thickness change on boundary offset.

[0045] The specific steps of S3 are: S301: Call the mutation trend distribution position in the thermal buffer surface correction gradient, extract the thermal retention change trajectory corresponding to the position, calculate the difference interval between the node temperature value and the average thermal retention level in the time series, establish a mapping relationship between the coordinate sequence and the retention change, and obtain the thermal retention trajectory offset; After calling the mutation trend distribution position in the thermal buffer surface correction gradient, the position data of all mutation points in the thermal buffer structure coordinate system are first extracted as a sequence. For example, 5 mutation points are detected in the range of 0mm to 30mm in the x-axis direction, which are located at x=3.5mm, 7.0mm, 13.2mm, 20.5mm, and 27.8mm respectively. Then, for each position point, its temperature record value is extracted in the corresponding thermal retention time series. The time series has a sampling period of 0.2 seconds and records a total of 25 groups of node temperature data within 5 seconds. The temperature change sequence over time is constructed, and the average thermal retention level of the sequence within 5 seconds is calculated. The average thermal retention value is obtained by summing the values of all temperature sampling points and dividing it by the total number of sampling points. For example, the sequence temperatures at x=13.2mm are 45.6℃, 46.1℃ to 47.8℃, and the average value is 46.9℃ Then, the temperature value of each time node is subtracted from the average value to obtain a difference interval sequence. Combined with the time point number corresponding to the node, a mapping relationship between the coordinate sequence (t1, t2, ..., tn) and the retention change value (ΔT1, ΔT2, ..., ΔTn) is established. For example, at t = 0.6s, ΔT = -1.3°C, and at t = 1.0s, ΔT = +0.8°C. The differences of all data points are repeatedly recorded, and the trend of these mapping sequences is further judged. The length and change amplitude of the continuous positive and negative segments are analyzed. If the positive and negative values alternate frequently and the difference amplitude exceeds the set threshold of ±1.5°C, it is considered that there is a thermal retention fluctuation offset at this point. Finally, all mutation points that meet the above conditions are marked as offsets according to the maximum amplitude in their difference sequence. For example, if the maximum ΔT value at x = 20.5mm is 2.1°C, then the thermal retention trajectory offset of this point is 2.1°C.

[0046] S302: Based on the thermal retention trajectory offset, determine the spatial region boundary corresponding to the temperature jump trend, extract the jump region and the adjacent hierarchical structure numbers, compare the connection mode of the structure number boundary with the jump coordinate intersection rate, select the overlapping structure segments, and obtain the structure jump connection relationship group; Based on the offset of the thermal retention trajectory obtained above, the presence of a temperature jump trend is determined point by point. A threshold of ±2.0°C is set for the judgment condition. If the temperature change value difference exceeds this value in a continuous time period, it is marked as a jump trend point. At the mutation point x=7.0mm, the ΔT of two consecutive points is +2.3°C and -2.2°C, which meets the jump condition. Then, the corresponding thermal structure area boundary is found in the structure diagram according to the coordinates of the point, and the boundary coordinates of the spatial sub-area where the point is located are extracted. The numbers of the adjacent structure areas are found in the four directions of up, down, left and right. For example, the 7.0mm point is located in Zone_03, and the adjacent structure areas are Zone_02 and Zone_04. After recording their corresponding numbers, the relationship between Zone_03 and the adjacent Z The connection between Zone_02 and Zone_04 is analyzed. Each set of connected area boundaries is segmented, and whether the jump coordinate point falls in the middle of the connected segment is recorded. If the jump point coincides with the connected boundary more than twice, it is considered that the jump and the boundary have an intersection relationship. The intersection rate threshold is set to 0.3, and the calculation method is the number of intersection points / total number of connected boundary segments. For example, if there are 5 connecting line segments and 2 intersection points, the intersection rate is 0.4, which is greater than the threshold and is determined to be an overlapping area. The intersection relationship between all jump points and the structural boundary is determined one by one. The structural segments that meet the conditions, such as Zone_03–Zone_04 and Zone_05–Zone_06, are screened out and combined into a structural jump connection relationship group.

[0047] S303: Based on the structure jump connection relationship group, the heat adsorption duration of the structure segment is counted and the average adsorption difference is calculated. All structure combinations are sorted in ascending order of the difference, the order of the jump nodes in the sorting chain is recorded, and adjacent jump structures are connected and the number sequence is combined to obtain the heat conduction path priority sequence; According to the structure jump connection relationship group, the duration of thermal adsorption in each pair of structural segments is counted. The thermal adsorption duration is defined as the total time period when the temperature exceeds the baseline (e.g., 45°C). In the Zone_03–Zone_04 structural pair, the adsorption time is 4.2 seconds, and in Zone_05–Zone_06 it is 3.5 seconds. After recording the thermal adsorption time of each structural pair, the corresponding adsorption difference is calculated. The difference is the absolute value difference of the thermal adsorption time of the two structural segments. For example, the adsorption time of Zone_03 is 4.2 seconds, and that of Zone_04 is 3.1 seconds, then the difference is 1.1 seconds. After counting the adsorption differences of all structural jump relationship groups, they are arranged from small to large to obtain a sorted sequence and the sorting position number of each jump node in the sequence is recorded, such as Zone_03–Zone_04 is ranked third and Zone_05–Zone_06 is ranked first. Finally, the adjacent structural segments in the numbered sequence are connected according to their original spatial positions to form the final priority order of the heat conduction path, such as Zone_05–Zone_06→Zone_03–Zone_04→Zone_07–Zone_08, to obtain the heat conduction path priority sequence.

[0048] The specific steps of S4 are: S401: Based on the combined position in the heat conduction path priority sequence, detect the lateral offset change of the heat conduction direction in the corresponding structure, extract the contour extension coordinates of the path intersection area, determine the angle difference between the extension direction and the original path direction, and obtain the lateral offset contour direction value; Based on the combination position in the heat conduction path priority sequence, the lateral offset change corresponding to the heat conduction direction in the structural area of each combination number is detected one by one. First, the starting and ending positions of the path segments Zone_01–Zone_02, Zone_03–Zone_04, etc. in the two-dimensional coordinate system are extracted, and the dominant direction vector of each path segment is calculated. For example, the direction of Zone_01–Zone_02 is from (2, 3) to (8, 5), and the dominant direction is right-up. Then, at each connection point, it is expanded 2mm in the horizontal direction (along the y-axis direction), and the contour coordinate point set after the boundary expansion is extracted to form the boundary line of the extended area. The dominant direction vector of the extension line is recorded. For example, the starting point of the extended boundary is (8, 5) and the end point is (8, 7), with the direction being the positive y-axis. The angle between the main direction vector of the original path and the direction vector of the extended boundary is then calculated. The angle value is determined by the angular difference between the two vector endpoints. For example, if the main direction is 63° and the extended direction is 90°, the angle difference is 27°. The above operation is repeated for all path segments, and the angle difference threshold is set to 15°. Areas exceeding this threshold are determined as areas with significant lateral offset. For example, if the angle of the Zone_03–Zone_04 segment is 34°, its offset direction value is recorded as a significant offset state, ultimately forming a set of matching datasets of path segments and their corresponding lateral offset contour direction values.

[0049] S402: Based on the lateral offset contour direction value, the area deformation position of the thermal extension contour in the path connection boundary is identified, the area change of the contour boundary in the projection plane is calculated, and compared with the lateral torsion rate of the heat conduction path. An offset overlap threshold is set and numbers with overlap rates greater than the threshold are selected to obtain an area deformation overlap sequence; The calculation formula for the area change of the contour boundary in the projection plane is as follows: ; in, Represents the area change of the contour boundary in the projection plane, The radian form represents the lateral offset profile direction value. represents the transverse torsion rate of the heat conduction path, Represents the area of the original contour boundary in the projection plane, Represents the contour boundary projection area after thermal expansion deformation, Represents the total number of segments after the contour boundary is discretized, The index number representing the discretization segment, represents the local extension strain rate of the dth boundary point along the x-axis, / represents the ductile strain gradient at the boundary of segment d; Parameter acquisition method and specific values (Lateral offset profile direction radian): The lateral offset angle is measured by a laser profile scanner and converted into radians. The measured value =0.3rad, based on the lateral offset angle range of 0.1~0.5rad in the engineering measured data; (Transverse torsion rate of heat conduction path): The ratio of path torsion angle to path length is measured by a thermal imager. =0.12 rad / m, based on the standard range of heat conduction path torsion rate of 0.05~0.2 rad / m; (Original projection area): Use a 3D projector to scan the original contour boundary and calculate the polygon area in the projection plane =50.0m²; (Deformed Projected Area): Repeatedly scan the contour after thermal expansion to calculate the deformed area =52.5m², area difference | - |=2.5m²; (Total number of discretized segments): Based on the contour boundary length L = 60m and the discretization spacing Δx = 0.5m, the number of segments is calculated as m = L / Δx = 120 segments; / (Extension strain gradient): The strain rate in the x direction of each segment d is collected by the strain sensor =0.0015, calculate the gradient / =0.0003 / m, based on the material thermal extension strain rate range of 0.001~0.003; Formula calculation process First calculation: ; Second calculation: ; sum: ; Parameter setting basis and value association and The square root of and the denominator: Square root term Considering the nonlinear effects of lateral offset and torsion rate on area change, the denominator ( + ) Linear superposition effect of equilibrium direction and torsion; / The cumulative term: the cumulative strain gradient of each segment reflects the cumulative effect of local deformation of the overall contour; =1.957m²: This result indicates the total projected area change of the contour boundary under the coupling of thermal expansion and torsion, which is used to compare with the preset threshold to screen contour segments with excessive overlap rates, and is directly related to the generation of the area deformation overlap sequence in the step.

[0050] S403: Calling the area deformation superposition sequence, screening the path numbers in the continuous offset state, extracting the continuous distribution segments of the area change and the direction offset value in the number structure, determining the stability trend of the hot channel expansion and recording the corresponding structure number, and obtaining the hot channel expansion path number group; After calling the area deformation superposition sequence, continuous state identification is performed on each numbered path segment. The continuous state judgment condition is set as the area change in adjacent path segments does not exceed 0.5mm², and the direction offset value difference does not exceed 10°. Taking Zone_03–Zone_04 and Zone_07–Zone_08 as examples, the area changes are 2.6mm² and 2.9mm² respectively, with a difference of 0.3mm², and the direction offset angles are 27° and 32° respectively, with a difference of 5°, which all meet the continuity requirements. The continuous path number segments that meet the conditions are extracted, and their area change sequence and direction offset are recorded. The stability of the hot channel expansion trend is judged by combining the offset value sequence with the data change trend. If the area change sequence in the continuous path segment is monotonically increasing or decreasing and the direction offset angle change does not exceed 15°, then the segment is considered a stable hot channel expansion path. For example, Zone_03–Zone_04→Zone_07–Zone_08→Zone_09–Zone_10 meets the area change of 2.6→2.9→3.1mm² and the direction angle of 27°→32°→35°, and the trend is stable. Finally, these continuous and stable structure number segments are combined and recorded as a hot channel expansion path number group.

[0051] The specific steps of S5 are: S501: Extract the winding layer spacing value and thermal deformation offset distribution curve of the corresponding structure based on the numbered area in the hot channel expansion path number group, calculate the spacing compression ratio and offset direction angle at the numbered position, determine whether the change trends of the two are consistent, and obtain the structural deformation alignment consistency; According to the numbered areas in the hot channel expansion path numbering group, the winding layer data of the corresponding structures such as Zone_03, Zone_07, and Zone_09 are extracted in turn to obtain the initial winding layer spacing value of each numbered segment and the actual spacing change value after thermal action. The spacing value is achieved by lateral measurement of the microscopic section image, and sampling and recording are performed every 0.1 mm. It is calibrated before and after heat treatment. For example, the initial winding layer spacing of Zone_03 is 0.75 mm, and after heat treatment it is 0.63 mm. The compression ratio is calculated as (0.75-0.63) / 0.75=0.16, that is, the compression ratio is 16%. Then, the thermal deformation offset distribution curve of the numbered area is extracted. The curve is the offset trajectory line formed by the change of the central axis with time in the two-dimensional coordinate system during the thermal action process. The coordinate point sequence is formed by sampling and recording every 0.2 seconds. For example, the offset point is from (5.0, 4 .5) deviates to (5.6, 5.3), and then to (6.0, 5.9). A vector is constructed from it to represent the offset direction, and then the angle formed by the spacing compression direction and the offset trajectory direction is calculated. For example, if the compression direction is from top to bottom and the offset direction is from lower left to upper right, the angle between the two is 42 degrees. The angle calculation is performed for all numbered areas, and the compression ratio and angle change trend are recorded for each numbered area. The criterion for judging whether the trend is consistent is that the compression ratio and the angle increase or decrease trend are in the same direction in at least three consecutive areas. For example, from Zone_03 to Zone_07, the compression ratio increases from 16% to 18%, and the angle increases from 42° to 47°. It is judged that the trend is consistent. Finally, the structural deformation alignment consistency corresponding to each structural number segment is output, which is calculated as the ratio of the number of consistent segments to the total number of segments. For example, if 4 out of 5 segments have the same trend, the alignment consistency is 0.8.

[0052] S502: Based on the structural deformation alignment consistency, the regions with the same compression direction under consecutive numbers are screened, the number of superimposed numbers in each region is counted, the boundary position of the structural contour after thermal action in the number range is extracted, and the offset trajectory vector is calculated to obtain the thermal compression deformation superimposition trend group; Based on the consistency of structural deformation alignment, the areas with the same compression direction are screened. The compression direction is determined by taking the consistency of the sign of the change in the spacing between adjacent layers. That is, if the compression amount is positive, the direction is compression, otherwise it is expansion. For example, if Zone_03 to Zone_07 are all positive, it means that the compression direction is consistent. After determining these consecutive numbered segments, the number of numbers contained in each segment is counted. For example, there are 3 numbers in Zone_03–Zone_05 and 3 numbers in Zone_07–Zone_09 to form a structural segment number group. Then, the structural contour boundary position formed by each segment number group after thermal action is extracted. The boundary extraction is based on the coordinates of the outermost layer of the structure in the image projection. The boundary coordinates are recorded point by point to form a contour trajectory. An offset vector sequence is formed for each contour trajectory. The direction vector is constructed by the coordinate difference between two adjacent boundary points. For example, the offset vector (0.4, 0.7) is formed from the boundary points (10.2, 9.1) to (10.6, 9.8) of Zone_03. In this way, an offset vector group for the entire area is constructed. The directional median value of all offset vectors in each structural segment is taken as the thermal compression deformation offset trend value of the segment. Finally, a thermal compression deformation superposition trend group consisting of the structural segment number and its offset direction is formed.

[0053] S503: Invoke the thermal compression deformation superposition trend group, integrate the superposition of regional compression direction and offset trend, establish a mapping network between numbered positions and offset vectors, identify regional sequences with stable compression directions, form structural connection combinations, and obtain a temperature-controlled fiber-optic inertial navigation design solution. After calling the hot-pressing deformation superposition trend group, the compression direction vector and the offset trend vector of each numbered structural segment are integrated, and a two-dimensional vector coordinate mapping table is constructed for each numbered position. For example, in Zone_03, the compression direction is (0, -1) and the offset trend is (1, 1). The angle relationship of the two directions is calculated in the unit circle to compare whether there is a combination with an angle less than 20°. If the angle is less than the set threshold, it is considered to be a superposition situation. The angle overlap threshold of 20° is set based on the determination of the maximum stable segment range corresponding to the minimum angle deviation in the previous structural segment offset overlap rate experiment. The judgment is performed on all areas segment by segment to extract the numbered segments with stable directions, that is, the compression direction and offset trend are The potential direction maintains an angle change of no more than 20° in multiple consecutive numbered segments, and the angle change trend continuously rises or falls for no more than two reversal points. A stable zone number sequence is constructed. For example, Zone_03–Zone_04–Zone_05 constitutes a stable compression offset zone, and Zone_07–Zone_08–Zone_09 constitutes another stable segment. Finally, using the stable compression direction zone sequence as a unit, the regional structures are connected to form a spatially continuous structure connection combination. The path number sequence in the combination is used as the node index to generate a stable heat conduction network structure. The structure combination path is output as the core number group of the thermal conduction path layout of the temperature-controlled fiber-optic inertial navigation design scheme.

[0054] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in the present invention should be included in the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A method for designing a temperature-controlled fiber-optic inertial navigation structure, characterized in that: The following steps are involved: S1: Detect the point where the heat flow direction changes at the fiber microcavity boundary, analyze the reflection direction offset trajectory and reverse expansion position, compare the reflection number distribution and the difference in the intersection area, integrate the offset area number, and generate the thermal direction splitting frequency density; S2: Extracting the thickness profile of the outer heat reflective material based on the thermal direction splitting frequency density, analyzing the consistency of the thickness fluctuation direction and the overlap rate of the retention surface deformation, identifying the mutation trend site and calculating the relationship between thickness and offset, and outputting the thermal buffer surface correction gradient; S3: calling the mutation site of the thermal buffer surface correction gradient, matching the temperature jump trend with the structural layer overlap range, screening the jump area connection and sorting the heat adsorption duration difference, and generating a heat conduction path priority sequence; S4: detecting the lateral offset morphology of the conduction direction according to the heat conduction path priority sequence, identifying the boundary area contour expansion direction and area deformation parameters, screening the continuous offset channel structure, and generating a heat channel expansion path number group; S5: Based on the thermal channel expansion path number group, the compression trend of the winding layer is compared with the degree of thermal deformation alignment, the consistency between the compression direction and the structural deformation is determined, the full structure offset trend data is integrated, and the temperature-controlled fiber optic inertial navigation design plan is output.

2. The temperature-controlled fiber-optic inertial navigation structure design method according to claim 1, characterized in that: The thermal direction splitting frequency density includes the reflection angle change sequence, the direction offset frequency distribution, and the heat flow reflection intersection area number; the thermal buffer surface correction gradient includes the material thickness difference distribution, the contour deformation offset parameter, and the abnormal mutation trend identification; the heat conduction path priority sequence includes the temperature jump area index, the structural connection relationship sequence, and the thermal adsorption continuity classification; the thermal channel expansion path number group includes the lateral offset structure number, the contour overlap deviation feature, and the area expansion channel identification; the temperature-controlled fiber optic inertial navigation design scheme includes the winding spacing compression mode, the thermal deformation alignment parameter, and the structural deformation trend distribution map.

3. The temperature-controlled fiber-optic inertial navigation structure design method according to claim 1, characterized in that: The specific steps of S1 are: S101: Obtain the heat flow direction change points at the microcavity boundary of the outer periphery of the fiber winding structure, select the boundary coordinates whose angle change rate exceeds the deflection threshold, and generate a heat flow direction mutation coordinate set; S102: Based on the heat flow direction mutation coordinate set, call the reflection direction parameter and the heat flow duration, record the reflection direction change sequence, extract the direction offset trajectory, count the rotation frequency and intersection position, and obtain the reflection trajectory intersection density distribution result; S103: According to the reflection trajectory cross density distribution result, determine the region number where the offset frequency is higher than the set threshold, extract the spatial distribution, and calculate the unit offset frequency density difference index to obtain the thermal direction splitting frequency density.

4. The temperature-controlled fiber-optic inertial navigation structure design method according to claim 3, characterized in that: The calculation formula of the unit offset frequency density difference index is specifically as follows: ； in, Representative area The number is The unit offset frequency density difference index, Representative area Middle The reflection traces are numbered The frequency offset of the unit, Representative area Middle The reflection traces are numbered The effective action path distance of the unit, Indicates the The reflection traces are numbered The unit's horizontal reflection angle, Indicates the The reflection traces are numbered The vertical reflection angle of the unit, Indicates area The frequency offset threshold is set. Indicates the number The unit in the area The total number of observed reflection tracks in .

5. The temperature-controlled fiber-optic inertial navigation structure design method according to claim 1, characterized in that: The specific steps of S2 are: S201: Acquire position data of a selected area in the thermal direction splitting frequency density, extract the thickness value of the outer heat reflective material in the transverse coordinate, establish a mapping relationship between the coordinate and the thickness, determine the direction of thickness change, and obtain a thickness change trajectory value of the heat reflective layer; S202: Based on the thickness change trajectory value of the heat reflective layer, extracting the boundary coordinates of the heat retention surface contour, determining the degree of overlap between the thickness change direction and the angle between the boundary direction, counting the number of deformation point distributions, and obtaining the thickness direction consistency overlap; S203: According to the thickness direction consistency overlap, the overlapping area numbers are screened, the thickness difference and the contour offset are extracted, the thickness offset derivative is calculated and a variation function is established to obtain the thermal buffer surface correction gradient.

6. The temperature-controlled fiber-optic inertial navigation structure design method according to claim 1, characterized in that: The specific steps of S3 are: S301: Calling the mutation trend distribution position in the thermal buffer surface correction gradient, extracting the thermal retention change trajectory corresponding to the position, calculating the difference interval between the node temperature value and the average thermal retention level in the time series, establishing a mapping relationship between the coordinate sequence and the retention change, and obtaining the thermal retention trajectory offset; S302: Based on the thermal retention trajectory offset, determine the spatial region boundary corresponding to the temperature jump trend, extract the jump region and the adjacent hierarchical structure numbers, compare the connection mode of the structure number boundary with the jump coordinate intersection rate, filter the overlapping structure segments, and obtain the structure jump connection relationship group; S303: Based on the structure jump connection relationship group, the thermal adsorption duration of the structure segment is counted and the average adsorption difference is calculated. All structure combinations are sorted in ascending order of the difference, the order number of the jump node in the sorting chain is recorded, and adjacent jump structures are connected and the number sequence is combined to obtain the heat conduction path priority sequence.

7. The temperature-controlled fiber-optic inertial navigation structure design method according to claim 1, characterized in that: The specific steps of S4 are: S401: Based on the combination position in the heat conduction path priority sequence, detecting the lateral offset change of the heat conduction direction in the corresponding structure, extracting the contour extension coordinates of the path intersection area, determining the angle difference between the extension direction and the original path direction, and obtaining the lateral offset contour direction value; S402: Identifying the area deformation position of the thermal extension contour in the path connection boundary based on the lateral offset contour direction value, calculating the area change of the contour boundary in the projection plane, and comparing it with the lateral torsion rate of the heat conduction path, setting an offset overlap threshold, and filtering the numbers with overlap rates greater than the threshold to obtain an area deformation overlap sequence; S403: Call the area deformation superposition sequence, filter the path numbers in the continuous offset state, extract the continuous distribution segments of the area change and the direction offset value in the number structure, determine the stability trend of the hot channel expansion and record the corresponding structure number, and obtain the hot channel expansion path number group.

8. The temperature-controlled fiber-optic inertial navigation system design method according to claim 7, characterized in that: The specific calculation formula for the area change of the contour boundary in the projection plane is: ； in, Represents the area change of the contour boundary in the projection plane, The radian form represents the lateral offset profile direction value. represents the transverse torsion rate of the heat conduction path, Represents the area of the original contour boundary in the projection plane, Represents the contour boundary projection area after thermal expansion deformation, Represents the total number of segments after the contour boundary is discretized, The index number representing the discretization segment, represents the local extension strain rate of the dth boundary point along the x-axis, / Represents the ductile strain gradient at the boundary of segment d.

9. The temperature-controlled fiber-optic inertial navigation system design method according to claim 1, characterized in that: The specific steps of S5 are: S501: Extracting the winding layer spacing value and thermal deformation offset distribution curve of the corresponding structure based on the numbered area in the hot channel expansion path number group, calculating the spacing compression ratio and offset direction angle at the numbered position, determining whether the change trends of the two are consistent, and obtaining the structural deformation alignment consistency; S502: Based on the structural deformation alignment consistency, the regions with the same compression direction under consecutive numbering are screened, the number of superimposed numbers in each region is counted, the boundary position of the structural contour after thermal action in the numbering range is extracted, and the offset trajectory vector is calculated to obtain a thermal-compression deformation superposition trend group; S503: Call the thermal compression deformation superposition trend group, integrate the superposition of regional compression direction and offset trend, establish a mapping network between numbered positions and offset vectors, identify regional sequences with stable compression directions and form structural connection combinations, and obtain a temperature-controlled fiber-optic inertial navigation design solution.

10. A temperature-controlled fiber-optic inertial navigation system (FIINS) structure designed according to the temperature-controlled fiber-optic inertial navigation system (FIINS) structure design method according to any one of claims 1 to 9, the temperature-controlled fiber-optic inertial navigation system (FIINS) structure comprising a housing, a polarization controller, a fiber ring, a light source, a beam splitter, and an interferometer detector disposed within the housing, and a support structure, a feedback control circuit, and a power interface module disposed on a surface of the housing.

Images