Disclosure of Invention
The invention mainly aims to provide a remote sensor thermal control index formulation method based on force-thermal coupling, which researches the coupling influence mechanism of gravity release and temperature change on the imaging quality of a system by analyzing the environmental difference between space and the ground, takes the requirement of the imaging quality index on the temperature as a constraint condition, and effectively reduces the power consumption of the thermal control system, reduces the design difficulty of the thermal control system and shortens the development period on the premise of ensuring that the design of the thermal control system meets the requirement of the imaging quality index.
The technical scheme adopted by the invention is that the method for formulating the thermal control index of the space remote sensor based on the force thermal coupling algorithm comprises the following steps:
Establishing four hypothetical conditions of force thermal coupling;
calculating the value of MTF for quantitatively evaluating imaging quality according to the force thermal coupling elastic deformation theory and the conversion relation between the optical-mechanical structure deformation and the optical system;
Determining a corresponding comprehensive imaging quality index under the coupling condition of a remote sensor thermal control subsystem and a mechanical subsystem according to the relation between the integral index of the imaging quality and the subsystem index;
And carrying out iterative computation on the coupling working conditions, and taking the corresponding temperature when the imaging quality threshold is reached as a thermal control index to finish the conversion from the optical index to the thermal control index.
Further, the assumed working conditions for establishing the four force thermal couplings include:
Before the establishment of the thermal control index, the force and heat environment of the remote sensor during on-orbit operation is described and assumed, a unified temperature range is generally defined for the whole system during the establishment of the thermal control index of the remote sensor, and four temperature working conditions are established from two aspects of uniform temperature load and uniform gradient change temperature load;
And carrying out gravity working condition assumption on the remote sensor according to the gravity direction of actual assembly and debugging, and combining the assumed four temperature working conditions with one gravity working condition to establish four coupling working conditions which are respectively (delta T) and G, X (delta T) and G, Y (delta T) and G and Z (delta T) and G.
Still further, the calculating the value of MTF for quantitatively evaluating imaging quality includes:
calculating the elastic deformation of the remote sensor under the coupling working condition by adopting a finite element analysis method of coupling elastic mechanics to ensure the formulating precision of a thermal control index, and obtaining the coordinate information of each deformed finite element node through the force-thermal coupling calculation of the remote sensor FEM, wherein the coordinate information can be used as the input condition of optical calculation;
The influence of the structural deformation of the remote sensor on the imaging quality can be analyzed through optical calculation, the optical surface nodes are rebuilt into a software usable form through a method of wavefront fitting based on Zrenike polynomials, and then all rebuilt optical surfaces are sequentially input into optical software to generate a new optical system, so that the change condition of MTF under the working condition of force-thermal coupling is analyzed;
Wherein Zrenike is mathematically described as:
Wherein K is a conic coefficient, c is a curvature, A i is a coefficient of a polynomial, Z i is a polynomial, In order to normalize the radius,Is the amplitude angle;
And comparing the calculated MTF with the MTF required by the imaging quality index, and taking the temperature corresponding to the coupling working condition when reaching the imaging quality threshold as a thermal control index.
The invention has the advantages that:
The index of the invention can reduce the heat control power consumption by about 20 percent and shortens the development period of the remote sensor to a certain extent. According to the prediction result of the MTF change condition of the remote sensor in orbit, the method can just ensure that the imaging quality in orbit meets the index, and has higher reliability compared with the index formulated by only considering the temperature.
In addition to the objects, features and advantages described above, the present invention has other objects, features and advantages. The present invention will be described in further detail with reference to the drawings.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Thermal control index formulation method
The invention provides a thermal control index formulation method taking an imaging quality index as a requirement under the condition of force thermal coupling. Fig. 1 depicts the specific process of the method, firstly, because the space environment can cause the remote sensor to generate gravity release and temperature change, four hypothetical conditions of force-heat coupling are established. And secondly, calculating the MTF value for quantitatively evaluating the imaging quality according to the force-thermal coupling elastic deformation theory and the conversion relation between the optical-mechanical structure deformation and the optical system. And then, determining the corresponding comprehensive imaging quality index under the coupling condition of the remote sensor thermal control subsystem and the mechanical subsystem according to the relation between the integral index of the imaging quality and the subsystem index. And finally, carrying out iterative computation on the coupling working conditions, and taking the corresponding temperature when the imaging quality threshold is reached as a thermal control index to finish the conversion from the optical index to the thermal control index.
Working condition hypothesis
Because the on-orbit imaging quality of the remote sensor can be influenced by the coupling of temperature change and gravity release, the invention uses an optomechanical thermal integration analysis method to calculate the structural deformation of the remote sensor caused by the force-thermal coupling environment and the imaging quality of an optical system after deformation.
The space environment (microgravity and space thermal environment) where the remote sensor is located is different from the ground environment, and the difference can cause gravity release and temperature change when the remote sensor is on track. The influence of gravity release and temperature change on the imaging quality of the remote sensor is inseparable, so that the force-heat coupling effect is comprehensively considered to formulate a heat control index. Before the establishment of the thermal control index, the force-thermal environment of the remote sensor when in orbit needs to be described and assumed:
Assume 1 temperature regime
Because the temperature distribution of each node inside the remote sensor is irregular in the space environment, a uniform temperature range is generally defined for the whole system when the thermal control index of the remote sensor is formulated. The invention establishes four temperature conditions as shown in table 1 from two aspects of uniform temperature load and uniform gradient temperature load. The initial temperature for all conditions was set to 20 ℃ taking into account the actual temperature at ground machining and conditioning. No. 1 is a uniform temperature load working condition, and the temperature value (20+/-delta T) of the working condition is changed by changing the value of the variation delta T of the temperature. 2-4 is a temperature load working condition with uniform gradient change, and the temperature gradient change has directivity due to the asymmetry of a remote sensor structure, so that the working condition is described by adopting X, Y, Z directions of a Cartesian coordinate system, and the temperature value taking 20 ℃ as a symmetry center is 20-delta T/2-20+delta T/2, wherein the positive and negative of the temperature change quantity delta T represent the change trend.
Table 1 description of the working conditions
Let 2 be the gravity condition
The remote sensor can offset the deformation caused by gravity to the structure through assembly and debugging when on the ground, but after entering a rail, new structural internal stress can be generated due to the release of the gravity so as to cause the new deformation of the remote sensor structure, so that the gravity working condition assumption is needed to be carried out on the remote sensor according to the gravity direction of actual assembly and debugging. The direction of gravity of the remote sensor in the present invention is-Y, and the gravity condition is hereinafter designated G.
The invention combines the assumed four temperature working conditions with a gravity working condition to establish four coupling working conditions which are respectively (delta T) and G, X (delta T) and G, Y (delta T) and G and Z (delta T) and G. The invention uses the coupling symbol to name the working condition in a simplified way.
Determination of thermal control indicators
Since the imaging quality is very sensitive to deformations of the opto-mechanical structure of the remote sensor, an accurate method is needed to calculate the structural deformations. The result of the force-thermal coupling method for calculating the elastic deformation is more accurate than that of a method for respectively calculating the force and the thermal deformation and then superposing, so that the method adopts a finite element analysis method of coupling elastic mechanics to calculate the elastic deformation of the remote sensor under the coupling working condition to ensure the formulating precision of the thermal control index. And (3) obtaining coordinate information of each deformed finite element node through force-thermal coupling calculation of the remote sensor FEM, wherein the coordinate information can be used as an input condition of optical calculation.
The influence of the structural deformation of the remote sensor on the imaging quality can be analyzed through optical calculation. Considering that the optical software has no data interface of finite element discrete nodes, the invention reconstructs the optical surface nodes into a software usable form by a method of wavefront fitting based on Zrenike polynomials. And then sequentially inputting all reconstructed optical surfaces into optical software to generate a new optical system, and further analyzing the variation condition of MTF under the working condition of force-thermal coupling.
Wherein Zrenike is mathematically described as:
Wherein K is a conic coefficient, c is a curvature, A i is a coefficient of a polynomial, Z i is a polynomial, In order to normalize the radius,Is the argument.
By iterating the different coupling conditions, the corresponding MTF value can be obtained. And comparing the calculated MTF with the MTF which is not required by the imaging quality, and taking the temperature corresponding to the coupling working condition when the coupling working condition is brought to the imaging quality threshold as a thermal control index.
Example calculation
Description of the optical machine structure
The invention takes a coaxial reflection type remote sensor as a research object, the imaging range is within the visible light wave band range, and the structural components of the coaxial reflection type remote sensor comprise a load body (a main mirror component, a secondary mirror component, a three-mirror component, a reflecting mirror component, a focusing mechanism, a main bearing component, a light shield component and the like), an electronic system, a thermal control system and the like. As shown in FIG. 2, the optical design of the remote sensor uses the direction of the main optical axis as +Z axis, the refractive direction of the refractive mirror as +Y axis, a Cartesian coordinate system is established by the right hand rule, and the geometric model and the finite element model of the remote sensor use the same coordinate system.
Determination of imaging quality index
The imaging quality index of the remote sensor is a precondition for the formulation of a thermal control index, wherein each subsystem of the remote sensor has a corresponding imaging quality index, and the indexes are quantitatively described through respective MTF values. From the above, the effects of temperature and gravity on imaging quality are mutually coupled, and the imaging quality index corresponding to the thermal control or mechanical subsystem cannot be independently formulated, so the invention calculates the comprehensive imaging quality index (MTF D) corresponding to the coupling condition of the thermal control and mechanical subsystem according to the relation between the MTFs of the systems.
The static MTF (MTF S) measured in the laboratory environment is a total indicator of the imaging quality of the remote sensor before emission, the present invention requires that MTF S be ≡0.12 at frequency=57.1 lp/mm. MTFs are considered to be the result of the combined actions of MTF D, process and tune transfer function (MTF m), and optoelectronic data transfer and process transfer function (MTF CCD), and have the following relationships.
In the case of a CCD device,The pixel size transfer function (MTF geometry, MTF G) and the electronic system transfer function (MTF E) are formed, and the relation is that:
MTF G can be expressed as:
Wherein f is the spatial frequency, f n is the Nyquist frequency, a is the CCD pixel size, and d is the CCD pixel spacing. When d=a, f=f n, MTF G =0.637 can be calculated.
Where MTF E =0.9 is given based on the study staff's statistics on the electronic system technology and the state of the art. MTF M was established empirically and through big data statistics, and MTF M =0.85 was obtained after the technician tested the remote sensor of the present invention.
In combination, the MTF D is more than or equal to 0.245, so that the imaging quality requirement can be met.
Determination of thermal control indicators
According to the glosser thermal integration analysis method, the step length of temperature change is 0.1 ℃, and after total iterative computation is carried out on 4 coupling working conditions for 580 times, the relation between MTF and temperature under the 4 coupling working conditions is obtained, wherein the threshold value of the MTF is 0.245. The temperature gradient sensitivity of the remote sensor to various working conditions is different, and the reasonability of the working condition type assumption is also proved. This type of assumption is reasonable because of the different sensitivity of the remote sensor to temperature gradients in different directions. The MTF of the edge view field in each working condition is most influenced by temperature change, and the threshold value of the MTF is reached first, which shows that the sensitivity of the edge view field to temperature is higher than that of other view fields, so that the temperature corresponding to the MTF=0.245 of the edge view field is used as a thermal control index.
Since the optical elements of the remote sensor are most spatially distributed along the Z direction, the maximum temperature difference between the optical elements occurs in the critical operating condition temperature distribution calculated by the total of force and thermal coupling in the G & Z (delta T) operating condition, and the maximum temperature difference between the optical elements is 4.44 ℃.
In summary, for the safety and reliability of the remote sensor, the thermal control index is formulated as follows:
a) The whole temperature range of the remote sensor is 17.3 < T < 22.6 ℃;
b) The temperature difference in the X direction is less than 5.6 ℃, the temperature difference in the Y direction is less than 6.6 ℃, and the temperature difference in the Z direction is less than 6.9 ℃;
c) The temperature difference between the optical elements is <4.4 ℃.
Simulation analysis and verification
Because the real space environment and the assumed temperature working condition are different, the temperature condition of the remote sensor during on-orbit needs to be simulated more accurately through the orbit thermal simulation analysis. The feasibility of remote sensor thermal control design and the rationality of the thermal control index formulation method can be checked by calculating and comparing the change conditions of the on-orbit imaging quality under different thermal control index formulation methods.
Rail environment analysis
The remote sensor can cause temperature change of the optical-mechanical structure due to the influences of vacuum, cold black, internal heat sources, heat radiation, heat conduction and other factors in the space. The temperature change can lead the optical element to generate displacement, distortion, thermal deformation and the like, thereby further reducing the imaging quality of the remote sensor, so that the remote sensor can perform normal imaging work only after reaching heat balance with the help of a heat control system.
The heat balance equation can be expressed as:
In the formula, Is solar radiation; is earth return; Infrared radiation of the earth; A heat source within the spacecraft; Heat radiated by the spacecraft to the space; the internal energy of the spacecraft can be changed.
The remote sensor is set to run on the sun synchronous orbit, orbit parameters are shown in table 2, and the change trend of the sun illumination angle (beta angle) can be calculated through orbit information, so that the change of the external heat flow of the remote sensor can be calculated. The change in beta angle can be described as:
Wherein I is the inclination angle between the satellite orbit plane and the earth orbit plane, I is the included angle between the earth orbit plane and the earth orbit plane, The yellow meridian of the sun, omega is the red meridian of the ascending intersection point.
Rail thermal simulation analysis
The thermal control index according to the present invention is used for designing a thermal control system. And a thermal structure analysis model of the remote sensor and the satellite platform is established. And carrying out orbit thermal simulation analysis on the whole satellite by taking an orbit environment analysis result as an input condition and 20 ℃ as a thermal control target temperature, and carrying out iterative formulation on a thermal control system by taking the simulation result as a guide until the temperature meets the index requirement.
In an actual space environment, the external heat flow of a satellite, the temperature of a platform, the attitude of the satellite and other factors can be changed, and a thermal control coating and various lacquers can be degraded with time, so that two extreme working conditions of high temperature and low temperature can be generated. According to the invention, during simulation, two extreme working conditions are analyzed, and simulation of a plurality of track periods is performed, wherein the thermal control power consumption required by the high-temperature working condition and the low-temperature working condition is respectively 98.0w and 103.5w, and the thermal control power consumption of the low-temperature working condition is larger.
FIG. 3 shows the temperature fluctuation of each optical component of the remote sensor in 10 periods, and the temperature reaches equilibrium after 3 tracks and meets the thermal control index.
The threshold value (0.245) of the MTF D is used as the index requirement of the MTF when the remote sensor is on track, the coupling relation between gravity and temperature is comprehensively considered, the temperature of each finite element node of the remote sensor obtained in the rail thermal simulation analysis is used as a temperature working condition, and the coupling working condition is built by combining with the gravity working condition. The deformation condition of the optical-mechanical structure of the remote sensor during on-orbit can be calculated by the force-thermal coupling principle. And extracting deformed optical surface node data as surface shape fitting input conditions, loading fitting results into optical software, and calculating MTF (maximum transfer function) of the remote sensor at any moment in imaging time when the remote sensor is on track.
According to the task requirement, the imaging time of each rail of the remote sensor is less than 10min, and the data in the imaging time of the 8 th rail (40176-406 s) after the remote sensor reaches the heat balance are selected for research. Wherein fig. 4 shows the MTF at a point in time when different indicators are on track in comparison with the initial MTF of the optical system. And taking 20s as a step length, carrying out iterative computation on data in the imaging time for 60 times, and obtaining the MTF change condition of different indexes in-orbit as shown in figure 4.
As can be seen from fig. 4, the present index allows the MTF to meet the index during the imaging time, to reach a minimum value (0.248) at 40616s and 40716s and just meet the threshold. This means that the index can effectively ensure the imaging quality of the remote sensor in orbit and avoid the waste of the power consumption of the thermal control system.
Thermal vacuum experiments
Remote sensors require detection of MTFs by thermal vacuum experiments prior to transmission. The thermal vacuum experiment is to simulate the orbit environment by using a space environment simulator, and calculate the energy information of the target image by using optical test equipment to reflect the imaging quality of the remote sensor. The test system comprises a camera, a collimator system, an interferometer, an infrared radiator, a satellite mounting platform, a satellite platform simulation bin and the like, and as shown in fig. 5, the systems are respectively mounted on a vacuum tank platform. The vacuum tank can provide cold black temperature below 100k and vacuum environment of 1× -4 Pa to simulate real space environment. According to the actual scheme, the remote sensor is arranged on a satellite platform, and a heating system is arranged in a satellite platform simulation bin to simulate the real temperature change of the platform. The infrared heater can be used to provide external heat flow variation on each surface of the satellite. The interferometer is arranged on the focal plane of the collimation system, the emitted light outputs parallel light beams through the collimator, and the satellite is arranged on the high-precision turntable, so that MTFs S of different fields of view of the remote sensor can be detected. The experimental results are consistent with the overall temperature indexes (17.3 ℃ and 22.6 ℃) established by the method, and further prove the rationality of the method.
In order to ensure that the space observation task is successfully completed, the temperature of the optical remote sensor needs to be maintained in a proper range through a thermal control system, so that the imaging quality is prevented from being influenced by the space thermal environment. The design and formulation of the thermal control system need reasonable thermal control indexes for guidance, however, the existing thermal control index formulation method lacks enough theoretical basis and easily causes waste of resources and time, or only considers temperature change and ignores the influence of gravity release, so that the formulated thermal control indexes are inaccurate, and accurate pre-judgment of the proper temperature level of the whole machine when the remote sensor is in orbit is difficult to form. According to the invention, the coupling influence mechanism of gravity release and temperature change on the imaging quality of the system is researched by analyzing the environmental difference between space and ground, and the requirement of the imaging quality index on the temperature is taken as a constraint condition. The track thermal simulation result shows that compared with the empirical index, the index of the invention can reduce the thermal control power consumption by about 20 percent and shortens the remote sensor development period to a certain extent. In addition, the method can just ensure that the imaging quality in the on-orbit process meets the index and has higher reliability compared with the index formulated by only considering the temperature through the prediction result of the MTF change condition of the remote sensor in the on-orbit process. Then, the thermal vacuum experimental result shows that the temperature critical value of the remote sensor under the gravity release condition accords with the thermal control index formulated by the method, and the rationality of the method is further verified. Finally, through further verification, the method can be applied to remote sensors of various types (off-axis reflection type, transmission type and the like).
The invention establishes a remote sensor thermal control index formulation method based on force-thermal coupling by researching the comprehensive influence of optics, mechanics and heat on the imaging quality of the remote sensor. The method can accurately predict the proper temperature level of the remote sensor in orbit so as to guide the design of the thermal control system to meet the requirement of imaging quality on temperature. In addition, the power consumption of the thermal control system is reduced to a certain extent, the design difficulty of the thermal control system is reduced, and the development period is shortened. The main contents are as follows:
(1) 4 kinds of force-heat coupling assumption working conditions are established for the remote sensor according to the environmental difference between the space and the ground. And analyzing the change condition of the MTF of the remote sensor under the coupling working condition according to the force-heat coupling calculation method and the surface shape fitting method.
(2) According to the relation between the integral index of imaging quality and the subsystem index, the comprehensive imaging quality index corresponding to the condition of remote sensor thermal control and mechanical subsystem coupling is determined to be MTF D >0.245, and the temperature corresponding to the condition of reaching the imaging quality index is used as a thermal control index by carrying out iterative calculation on the coupling working condition.
(3) After different indexes are applied to the remote sensor, and simulation comparison is carried out, the method provided by the invention is verified to save about 20% of power consumption of a thermal control system and shorten the development period on the premise of meeting the requirement of imaging quality. The thermal vacuum experimental result accords with the index of the invention, and further illustrates the rationality of the method of the invention.
Through further research, the method is also applicable to other types of remote sensors (off-axis reflection type, transmission type and the like). However, a remote sensor with an active optical system can resist the effect of temperature on imaging quality by adjusting the wavefront of the optical element, and thus the thermal control index formulation of the present invention is not applicable to this type of remote sensor. In conclusion, the method for formulating the thermal control index has wide applicability, can be well applied to practical engineering, and provides a certain contribution and help for the development of the aerospace industry.
The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the invention are intended to be included within the scope of the invention.