Method for integrated optimization of heat pipe constrained component layout based on mixed integer programming

Abstract

The application discloses a heat pipe constraint component layout comprehensive optimization method based on mixed integer programming, relates to the technical field of component layout optimization, and comprises the following steps: determining component information and heat pipe information in a component layout area; determining the number of heat pipes required by the components; constructing an allocation relationship optimization constraint model of the components and the heat pipes; taking the center positions of the components as optimization variables, taking the minimization of the maximum actual heat conduction power of the heat pipes as an optimization target, taking non-interference constraints between the components and the component layout area, non-interference constraints between the components, component system centroid constraints and heat coupling constraints between the components as constraint conditions, constructing an optimization target model and a constraint model; constructing an intersection constraint model of the components and the heat pipes; simultaneously solving various models to obtain and solve a component layout optimization model, and obtaining the positions of the components in the layout area. The application can realize the fast optimization solution of the component layout under the consideration of the quality characteristic requirements and the heat pipe constraints, has low solution difficulty and high efficiency.

Description

Technical Field

[0001] This invention relates to the field of component layout optimization technology, and in particular to a comprehensive optimization method for heat pipe-constrained component layout based on mixed integer programming. Background Technology

[0002] Satellite layout plays a crucial role in determining its on-orbit performance and functionality. The purpose of satellite layout design is to arrange electronic components or equipment in appropriate locations on the satellite to meet various system performance requirements, such as mass characteristics and thermal control. As an important part of the overall satellite design, satellite layout design directly determines the overall performance, development cost, design cycle, and design level of the satellite system.

[0003] Existing satellite systems typically deploy a certain number of heat pipes within a designated area to dissipate heat from components within that area, thus transferring the heat generated by the components outside the deployment area. To obtain a component layout scheme that satisfies both quality characteristics and heat pipe heat dissipation constraints, Chinese patent document CN115859571A, entitled "Component Layout Optimization Design Method Based on Integer Programming under Heat Pipe Constraints," discloses a component layout optimization design method based on integer programming under heat pipe constraints. This method decomposes the component layout optimization problem under heat pipe heat dissipation constraints into a component-heat pipe allocation optimization problem and a detailed component layout optimization problem. By solving the two optimization problems after decomposition, a component layout scheme that satisfies both quality characteristics and heat pipe heat dissipation constraints is obtained.

[0004] However, the aforementioned method employs a heuristic layout search approach based on sequence layout sampling when solving the detailed component layout optimization problem. In practical applications, this heuristic layout search method based on sequence layout sampling suffers from low search efficiency, high search difficulty, and low optimization efficiency. Furthermore, as the size of the component increases, its search difficulty and optimization efficiency increase significantly. Moreover, if other performance constraints are added when using the above method, the design optimization problem becomes even more complex and difficult to solve, making it impossible to find a layout scheme that meets the conditions.

[0005] Furthermore, existing heat pipe-constrained layout optimization methods typically assume no direct heat exchange between components, with heat dissipation occurring solely through heat pipes. However, in real-world satellite systems, the close proximity or contact between components leads to significant thermal coupling effects, resulting in localized heat accumulation, increased temperature gradients, and impacts system thermal stability and component lifespan. Summary of the Invention

[0006] To address some or all of the technical problems existing in the prior art, this invention provides a comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming.

[0007] The technical solution of the present invention is as follows: A comprehensive optimization method for heat pipe-constrained component layout based on mixed-integer programming is provided, including: Determine the component information and heat pipe information within the component layout area; Based on the structural dimensions of each component and the arrangement of the heat pipes, determine the number of heat pipes required for each component; Using constraints such as the component spanning the heat pipe, the heat pipe's thermal conductivity, and the longitudinal dimension capacity as constraints, discrete variables are used to describe the allocation relationship between the component and the heat pipe. An optimization constraint model for the allocation relationship between the component and the heat pipe is constructed using integer programming. The center position of the component is described by continuous coordinates. The center position of the component is used as the optimization variable. The optimization objective is to minimize the maximum actual heat conduction power of the heat pipe. The constraints are non-interference between components and component layout areas, non-interference between components, centroid constraint of component system, and thermal coupling constraint between components. The optimization objective model and the constraint model corresponding to each constraint condition are constructed by linear modeling. The thermal coupling constraint between components means that the interval between components is not less than a preset distance threshold, and the temperature difference between components with an interval less than the preset distance threshold does not exceed a preset temperature difference threshold. The center position of the component is described by continuous coordinates, and the position of the heat pipe is described by coordinates. Based on the distribution relationship between the component and the heat pipe, the intersection constraint model of the component and the heat pipe is constructed by linear modeling. By combining the component and heat pipe allocation relationship optimization constraint model, optimization objective model, constraint model and component and heat pipe intersection constraint model, the component layout optimization model is obtained. Solving the component layout optimization model, the position of each component in the component layout area is obtained.

[0008] In some alternative implementations, the optimization constraint model for the allocation relationship between components and heat pipes is expressed as follows: ; in, This represents an auxiliary binary variable used to describe the first... The reference heat pipe for each component, Indicates the first The heat pipe is the first one The reference heat pipe for each component is defined as the heat pipe with the smallest serial number among all the heat pipes used in the component. Indicates intermediate variables. This represents a binary variable used to describe the distribution relationship between components and heat pipes. Indicates the first The component is placed in the first On a heat pipe, Indicates the first The component was not placed in the first... On a heat pipe, Indicates the first Transformation matrix corresponding to each component The Line number Column elements, transformation matrix For one The matrix, and the matrix's first... The first line The elements up to the first One element is 1, and the rest of the elements of the matrix are 0. Indicates the first The height of each component Indicates the maximum permissible height. Indicates the first The power of each component Indicates the first The number of heat pipes required for each component This represents the maximum thermal conductivity of a given heat pipe. Indicates the number of components. Indicates the number of heat pipes; Based on the arrangement of the heat pipes, all heat pipes are sequentially numbered starting from number 1 along the arrangement direction.

[0009] In some optional implementations, the lower left corner of the component layout area is selected as the origin of the coordinate system, and the horizontal direction of the component layout area is used as the coordinate system. The vertical direction of the component layout area is along the axis. Construct a two-dimensional Cartesian coordinate system along the axes; Based on the constructed two-dimensional Cartesian coordinate system, the optimization objective model is represented as: ; in, This indicates a preset parameter. Used for equivalent substitution to optimize the objective function .

[0010] In some optional implementations, the constraint model corresponding to the non-interference constraint between components and component layout regions is expressed as: ; in, Indicates the first The center coordinates of each component Indicates the first The width of each component, Indicates the first The height of each component Indicates the width of the component layout area. Indicates the height of the component layout area.

[0011] In some alternative implementations, the constraint model corresponding to the non-interference constraint between components is represented as: ; in, Indicates an indicator variable. Indicates the first The center coordinates of each component Indicates the first The width of each component, Indicates the first The height of each component It is a positive number.

[0012] In some optional implementations, the constraint model corresponding to the centroid constraint of the component system is represented as follows: ; in, Indicates the first The quality of each component The ordinate represents the position of the expected centroid. The vertical coordinate representing the position of the centroid of the component system. This indicates the pre-set centroid deviation in the coordinate system. The component along the axial direction.

[0013] In some optional implementations, the constraint model corresponding to the thermal coupling constraint between components is represented as: ; in, Indicates the first The first component and the first Each component in Distance along the axis Indicates the first The component and the first Each component in Distance along the axis For positive integers, This indicates a preset distance threshold. , , , , It is a binary indicator variable. , Indicates the first Operating temperature of each component Indicates the first Operating temperature of each component This indicates the preset temperature difference threshold. For positive integers, Indicates the base temperature of the heat pipe. Indicates the first The power of each component Indicates the first The equivalent thermal resistance of each component to its associated heat pipe. This represents the equivalent thermal resistance of the coupling between components. Used to represent all related to the first The components whose interval is less than a preset distance threshold are the components of the first component. Additional temperature rise caused by individual components.

[0014] In some alternative implementations, the intersection constraint model between the component and the heat pipe is represented as follows: ; in, Indicates the width of the heat pipe. Indicates the first The x-coordinate of the center position of each heat pipe It is a binary variable. .

[0015] In some alternative implementations, the component layout optimization model is represented as: ; in, Indicates the component layout scheme. .

[0016] In some alternative implementations, a mathematical programming solver or an integer programming algorithm is used to solve the component layout optimization model.

[0017] The main advantages of the technical solution of this invention are as follows: The heat pipe-constrained component layout comprehensive optimization method based on mixed integer programming of the present invention decomposes the component layout optimization problem under heat pipe constraints into a component layout optimization subproblem with respect to the allocation relationship between components and heat pipes and a component layout optimization subproblem without considering the allocation relationship, and constructs the corresponding models. Then, the intersection constraint model between components and heat pipes constructed based on the allocation relationship between components and heat pipes is used to solve the models corresponding to the optimization subproblems simultaneously. It can achieve rapid optimization of component layout under the consideration of quality characteristic requirements and heat pipe constraints. The optimization solution is easy and efficient, and the optimization solution efficiency is less affected by the component size. It can still achieve rapid optimization solution when the component size is large. Attached Figure Description

[0018] The accompanying drawings, which are included to provide a further understanding of embodiments of the invention and constitute a part of this invention, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 A schematic diagram of a component layout considering heat pipe constraints provided in an embodiment of the present invention; Figure 2 The flowchart illustrates a comprehensive optimization method for heat pipe constraint component layout based on hybrid integer programming, provided in an embodiment of the present invention. Detailed Implementation

[0019] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below in conjunction with specific embodiments and corresponding drawings. Obviously, the described embodiments are only a part of the embodiments of this invention, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0020] The technical solutions provided by the embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0021] refer to Figure 1 In satellite systems, heat pipes in the component layout area are typically arranged in a horizontally uniform manner to achieve better heat transfer. The component layout optimization problem under heat pipe constraints mainly studies how to achieve better heat dissipation performance and quality characteristics of the system through component layout design, given the location of the heat pipes. The component layout optimization problem under heat pipe constraints includes the following performance design requirements: The centroid of the component system along the vertical direction of the component layout area should be located within the given centroid range and as close as possible to the desired centroid. The components must be placed on heat pipes to ensure that the heat generated by the components themselves can be dissipated through the heat pipes; The total heat dissipation of each heat pipe must not exceed its maximum allowable heat dissipation capacity; The total heat dissipation of different heat pipes should be as close as possible to ensure that the temperature distribution within the component layout area is as uniform as possible and to avoid heat concentration.

[0022] To solve the aforementioned component layout optimization problem under heat pipe constraints, improve optimization efficiency, and enhance the feasibility and quality of the obtained component layout scheme, this invention provides a comprehensive optimization method for heat pipe-constrained component layout based on mixed integer programming.

[0023] refer to Figure 2 The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming provided in this embodiment of the invention includes the following steps S1-S6: Step S1: Determine the component information and heat pipe information within the component layout area; In this embodiment of the invention, the component information includes: the number of components, the structural dimensions of each component, and the power of each component.

[0024] In this embodiment of the invention, the heat pipe information includes: the number of heat pipes, the structural dimensions of the heat pipes, the performance parameters of the heat pipes, and the arrangement of the heat pipes. The performance parameters of the heat pipes include: the maximum thermal conductivity of the heat pipes; the arrangement of the heat pipes includes: the position of each heat pipe within the component layout area.

[0025] In this embodiment of the invention, the component information and heat pipe information within the component layout area are determined based on the actual component layout.

[0026] Step S2: Determine the number of heat pipes required for each component based on the structural dimensions of each component and the arrangement of the heat pipes. In this embodiment of the invention, the number of heat pipes required for each component is determined based on the structural dimensions of each component within the component layout area and the arrangement of the heat pipes.

[0027] In this embodiment of the invention, the number of heat pipes required for each component represents the number of heat pipes that need to be contacted when the current component is placed within the component layout area. Given the structural dimensions of the component and the arrangement of the heat pipes, this number of heat pipes can be calculated based on the width of the component, the width of the heat pipes, and the spacing between adjacent heat pipes. Here, the width of the component represents the horizontal length of the component within the component layout area, and the width of the heat pipe represents the horizontal length of the heat pipe within the component layout area.

[0028] Step S3: Using the constraints of the component spanning the heat pipe, the heat pipe thermal conductivity capacity, and the longitudinal dimension capacity as constraints, discrete variables are used to describe the allocation relationship between the component and the heat pipe. An optimization constraint model for the allocation relationship between the component and the heat pipe is constructed using integer programming. In this embodiment of the invention, the component spanning a heat pipe constraint means that the component is placed on at least one heat pipe, and the positions of the multiple heat pipes spanned by the component are continuous.

[0029] In this embodiment of the invention, the heat pipe thermal conductivity capacity constraint means that the total thermal conductivity of the heat pipe does not exceed the maximum thermal conductivity allowed by the heat pipe itself. Specifically, it can also be expressed as the total power of all components placed on the heat pipe that is heated by the current heat pipe does not exceed the maximum thermal conductivity allowed by the heat pipe itself.

[0030] In this embodiment of the invention, the longitudinal dimension capacity constraint means that the total height of all components placed on each heat pipe does not exceed the maximum allowable height. The height of a component represents its length in the vertical direction of the component layout area.

[0031] In this embodiment of the invention, the maximum allowable height is predetermined according to the actual situation, and is usually the length of the component layout area in the vertical direction.

[0032] Step S4: The center position of the component is described by continuous coordinates. The center position of the component is used as the optimization variable. The optimization objective is to minimize the maximum actual heat conduction power of the heat pipe. The constraints are non-interference between components and component layout areas, non-interference between components, centroid constraint of component system, and thermal coupling constraint between components. The optimization objective model and the constraint model corresponding to each constraint condition are constructed by linear modeling. In this embodiment of the invention, based on the performance design requirements of the component layout optimization problem under heat pipe constraints, the center position of the component is taken as the optimization variable of the component layout optimization problem under heat pipe constraints, and minimizing the maximum actual heat conduction power of the heat pipe is taken as the optimization objective of the component layout optimization problem under heat pipe constraints. Based on the determined optimization variables and optimization objective, a corresponding optimization objective model is constructed through linear modeling.

[0033] In this embodiment of the invention, based on the performance design requirements of the component layout optimization problem under heat pipe constraints, the center position of the component is used as the optimization variable of the component layout optimization problem under heat pipe constraints, and the non-interference constraints between components and component layout regions, the non-interference constraints between components, the centroid constraints of the component system, and the thermal coupling constraints between components are used as the constraints of the component layout optimization problem under heat pipe constraints. Based on the determined optimization variables and constraints, the constraint models corresponding to each constraint condition are constructed through linear modeling.

[0034] In this embodiment of the invention, the non-interference constraint between components and component layout areas means that there is no interference between the boundaries of components and component layout areas; the non-interference constraint between components means that there is no interference between any two components; the centroid constraint of the component system includes the constraint on the calculation method of the centroid position of the component system and the deviation constraint between the centroid position of the component system and the desired centroid position; the thermal coupling constraint between components means that the interval between components is not less than a preset distance threshold, and the temperature difference between components with an interval less than the preset distance threshold does not exceed a preset temperature difference threshold.

[0035] Step S5: The center position of the component is described by continuous coordinates, and the position of the heat pipe is described by coordinates. Based on the distribution relationship between the component and the heat pipe, the intersection constraint model of the component and the heat pipe is constructed by linear modeling. In this embodiment of the invention, the intersection constraint between the component and the heat pipe means that there is a geometric intersection relationship between the component and all the heat pipes it crosses, but no geometric intersection relationship between the component and the other heat pipes.

[0036] Step S6: Combine the component and heat pipe allocation relationship optimization constraint model, optimization objective model, constraint model and component and heat pipe intersection constraint model to obtain the component layout optimization model. Solve the component layout optimization model to obtain the position of each component in the component layout area.

[0037] In this embodiment of the invention, since the component and heat pipe allocation relationship optimization constraint model can characterize the allocation relationship between the component and heat pipe, and the component and heat pipe intersection constraint model is constructed based on the allocation relationship between the component and heat pipe, the component and heat pipe allocation relationship optimization constraint model can be combined with the optimization target model and each constraint model through the component and heat pipe intersection constraint model, thereby obtaining the component layout optimization model.

[0038] In this embodiment of the invention, a component layout optimization model is constructed by combining the above-mentioned component and heat pipe allocation relationship optimization constraint model, optimization target model, various constraint models, and component and heat pipe intersection constraint model. The component layout optimization model is then solved directly using an existing mathematical programming solver or integer programming algorithm to obtain the corresponding component layout scheme, and thus obtain the position of each component in the component layout area.

[0039] The heat pipe-constrained component layout optimization method provided in this invention decomposes the component layout optimization problem under heat pipe constraints into a component layout optimization subproblem with respect to the allocation relationship between components and heat pipes and a component layout optimization subproblem without considering the allocation relationship. Corresponding models are constructed, and then the intersection constraint model between components and heat pipes, constructed based on the allocation relationship between components and heat pipes, is used to solve the models corresponding to the optimization subproblems simultaneously. This method enables rapid optimization of component layout considering quality characteristics and heat pipe constraints. The optimization solution is easy and efficient, and its efficiency is less affected by the component size, allowing for rapid optimization even with large component sizes.

[0040] Furthermore, in this embodiment of the invention, for ease of description and model building, all heat pipes are sequentially numbered starting from number 1 along the arrangement direction of the heat pipes according to their arrangement.

[0041] Furthermore, in this embodiment of the invention, based on the sequentially numbered heat pipes, and using the constraints of the component spanning the heat pipe, the heat pipe thermal conductivity capacity constraint, and the longitudinal dimension capacity constraint as constraints, discrete variables are used to describe the allocation relationship between the component and the heat pipe. The optimized constraint model for the allocation relationship between the component and the heat pipe, constructed through integer programming, is expressed as follows: ; in, This represents an auxiliary binary variable used to describe the first... The reference heat pipe for each component, Indicates the first The heat pipe is the first one The reference heat pipe for each component is defined as the heat pipe with the smallest serial number among all the heat pipes used in the component. Indicates intermediate variables. This represents a binary variable used to describe the distribution relationship between components and heat pipes. Indicates the first The component is placed in the first On a heat pipe, Indicates the first The component was not placed in the first... On a heat pipe, Indicates the first Transformation matrix corresponding to each component The Line number Column elements, transformation matrix For one The matrix, and the matrix's first... The first line The elements up to the first One element is 1, and the rest of the elements of the matrix are 0. Indicates the first The height of each component Indicates the maximum permissible height. Indicates the first The power of each component Indicates the first The number of heat pipes required for each component This represents the maximum thermal conductivity of a given heat pipe. Indicates the number of components. Indicates the number of heat pipes.

[0042] In the embodiments of the present invention, in the above-described optimization constraint model for the allocation relationship between components and heat pipes, and Used to define the constraints of a component spanning heat pipes. When a component spans one or more heat pipes, the heat pipe with the smallest sequence number is selected as its unique reference heat pipe. This is used to limit each component to having one and only one reference heat pipe. Descriptive variables used to define the reference heat pipe Descriptive variables relating to actual heat pipe distribution The linear transformation relationship between them.

[0043] In this embodiment of the invention, the first Transformation matrix corresponding to each component The dimensions of the components and the heat pipes are determined in advance by solving the problem.

[0044] In an embodiment of the present invention, This is used to define the longitudinal dimension capacity constraint. If this condition is met, the longitudinal dimension capacity constraint is satisfied.

[0045] In an embodiment of the present invention, This is used to define the heat pipe thermal conductivity constraint. If this condition is met, the heat pipe thermal conductivity constraint is satisfied.

[0046] In this embodiment of the invention, binary variables This is used to determine which heat pipes the components are placed on, i.e., to determine the component-heat pipe distribution relationship. Specifically, when the... Each component is placed When using a heat pipe, a binary variable satisfy and .

[0047] Furthermore, in this embodiment of the invention, in order to facilitate description and model building, improve optimization and solution efficiency, and ensure the feasibility of the obtained component layout scheme, it is set that: all components are square structures, all components are rigid bodies with uniform mass distribution, the centroid of the component coincides with the geometric center of the component, and the components are arranged in the component layout area in a horizontal direction parallel to the component layout area.

[0048] Furthermore, in this embodiment of the invention, in order to facilitate description and model construction, a point in the component layout area is selected as the origin of the coordinate system, and a Cartesian coordinate system is constructed. The center position of the component is described based on the constructed Cartesian coordinate system.

[0049] Specifically, to further facilitate description and model building, the lower left corner of the component layout area is selected as the origin of the coordinate system, and the horizontal direction of the component layout area is taken as the coordinate system. The vertical direction of the component layout area is along the axis. Construct a two-dimensional Cartesian coordinate system along the axis.

[0050] Furthermore, in this embodiment of the invention, the two-dimensional Cartesian coordinate system constructed in the above manner is further defined as follows: The center coordinates of each component are , No. The width of each component is , No. The height of each component is The width of the component layout area is The height of the component layout area is .

[0051] It should be noted that the width of a component represents its length in the horizontal direction of the component layout area, and the height of a component represents its length in the vertical direction of the component layout area. The width of the component layout area represents the length of the component layout area in the horizontal direction, and the height of the component layout area represents the length of the component layout area in the vertical direction.

[0052] It should be noted that in this embodiment of the invention, the ability of components to rotate is not considered when optimizing the component layout.

[0053] Furthermore, based on the above settings, the following optimization objective model is constructed: ; in, This indicates a preset parameter. Used for equivalent substitution to optimize the objective function .

[0054] In this embodiment of the invention, based on the above settings, preset parameters are... Used for equivalent substitution to optimize the objective function Furthermore, using preset parameters Equivalent substitution optimization objective function At the same time, add constraints. When constraints When satisfied, it means ,when If the minimum value is reached, then the condition must be met. .

[0055] In this embodiment of the invention, by using preset parameters to equivalently replace the objective function and setting corresponding constraints, the use of maximum value calculation can be avoided, thereby improving the efficiency of optimization.

[0056] Furthermore, based on the above settings, the constraint model corresponding to the non-interference constraint between components and component layout areas is constructed as follows: ; in, Indicates the first The center coordinates of each component Indicates the first The width of each component, Indicates the first The height of each component Indicates the width of the component layout area. Indicates the height of the component layout area. Indicates the number of components.

[0057] Furthermore, in this embodiment of the invention, the Phi function is used to determine the formula for calculating the interference between any two components, an indicator variable is introduced and the method of large numbers is combined to linearize the formula for calculating the interference, thereby modeling the non-interference constraint between components.

[0058] In this embodiment of the invention, the following is defined: The component and the first The interference between the components is To ensure that there is no interference between components, the amount of interference... The following conditions must be met: .

[0059] In this embodiment of the invention, the Phi function method is used to calculate the interference between any two components.

[0060] Specifically, the two components are respectively the first The component and the first Taking one component as an example, the formula for calculating the interference between two components is as follows: ; in, Indicates the first The center coordinates of each component Indicates the first The width of each component, Indicates the first The height of each component.

[0061] Expanding the absolute value operation in the above formula for calculating the interferometric amount, the formula for calculating the interferometric amount can be equivalently transformed into: .

[0062] Furthermore, a binary indicator variable of 0 / 1 is introduced. By linearizing the above interference calculation formula using the method of large numbers, we obtain the following equivalent constraint model, which is the constraint model corresponding to the non-interference constraint between components: ; in, It should be a relatively large positive number, set according to the specific situation, for example... This requires that all inequalities in the above-mentioned non-interference constraints between components be satisfied simultaneously.

[0063] In the above-mentioned non-interference constraints between components, when At that time, Can be converted If the inequality If true, it means that the non-interference constraint between components is satisfied; when At that time, Can be converted ,because Since it is a large positive constant, this inequality applies regardless of whether interference occurs between components. This holds true consistently. In this embodiment of the invention, through constraints... It can be required that at least one of the four indicator variables is equal to 1, which can guarantee that the component layout scheme obtained by solving satisfies the non-interference constraint between components.

[0064] Furthermore, in this embodiment of the invention, the following is defined: the first The mass of each component is The expected coordinates of the centroid are .

[0065] In this embodiment of the invention, since the component layout optimization problem under heat pipe constraints only considers the vertical centroid constraint, the constraint model corresponding to the centroid constraint of the component system is constructed as follows: ; in, The vertical coordinate representing the position of the centroid of the component system. This indicates the pre-set centroid deviation in the coordinate system. The component along the axial direction.

[0066] Furthermore, considering that in actual satellite operation, components in the component layout area will experience direct heat transfer due to their geometric proximity, i.e., a thermal coupling effect. This thermal coupling effect can lead to local heat accumulation, altering the temperature distribution of the components themselves and potentially affecting the heat dissipation efficiency of the heat pipes. To ensure the thermal stability of the component layout scheme, this embodiment of the invention introduces inter-component thermal coupling constraints. These constraints are used to ensure that the interval between components is not less than a preset distance threshold, and that the temperature difference between components with intervals less than the preset distance threshold does not exceed a preset temperature difference threshold.

[0067] In this embodiment of the invention, in order to integrate the thermal coupling constraints between components into the mixed integer programming framework, a discretization method based on the component projection interval is used to determine the interval between components. Simultaneously, to simplify calculations and maintain linearity, a conservative criterion is adopted: if two components are in... Axial direction and If the distances along the axis are all less than the preset distance threshold, then the interval between the two components is determined to be less than the preset distance threshold.

[0068] Specifically, with the first The component and the first Taking one component as an example, two components in Axial direction and The distances along the axes are expressed as follows: ; ; in, Indicates the first The component and the first Each component in Distance along the axis Indicates the first The component and the first Each component in Distance along the axial direction.

[0069] In this embodiment of the invention, in order to construct thermal coupling constraints between components, three sets of binary indicator variables are introduced, specifically including: a first set of binary indicator variables. and , The first group of binary indicator variables is used to eliminate absolute value operations in the distance calculation formula; the second group of binary indicator variables... and , The second set of binary indicator variables is used to represent the two components in Axial direction and Whether the distance along the axis is less than a preset distance threshold; the third group of binary indicator variables. , The third set of binary indicator variables is used to indicate whether thermal coupling exists between two components.

[0070] Based on the first set of binary indicator variables, and by linearizing the above distance calculation formula using the method of large numbers, the following equivalent form is obtained: ; ; in, This should be a relatively large positive number. The specific value should be set according to the actual situation, for example, set to the diagonal length of the layout area.

[0071] Based on the second set of binary indicator variables, the following distance determination constraints are established: ; in, This indicates a preset distance threshold, which can be set according to the actual situation.

[0072] The distance determination constraint described above is used to ensure that if and only if hour, ; and used to ensure that if and only if hour, .

[0073] Based on the established second and third sets of binary indicator variables, the following thermal coupling decision constraints are established: ; The above thermal coupling determination constraint is used to ensure that if and only if hour, That is, if and only if the two components are Axial direction and When the distances in both axial directions are less than a preset distance threshold, it is determined that there is thermal coupling between the two components.

[0074] Based on the established third set of binary indicator variables, the following conditional temperature difference constraint is established: ; in, Indicates the first Operating temperature of each component Indicates the first Operating temperature of each component This indicates a preset temperature difference threshold, which can be set according to actual conditions. This should be a relatively large normal number. The specific value should be set according to the actual situation, such as the maximum operating temperature of the component.

[0075] The aforementioned conditional temperature difference constraint is used to ensure that if and only if hour, That is, if and only if the two components are Axial direction and When the distances in both axial directions are less than a preset distance threshold, the temperature difference between the two components must not exceed a preset temperature difference threshold.

[0076] Furthermore, each component is abstracted as a thermal node. The temperature of a component is determined by its own power and the heat exchange between the component and the outside environment. Let: The equivalent thermal resistance of each component to its corresponding heat pipe is The equivalent thermal resistance of the coupling between components is The equivalent thermal resistance from the component to its heat pipe is determined by the component's packaging form, internal structure, and the characteristics of the mounting interface (such as a thermal pad), and is used as a known constant input. The coupling equivalent thermal resistance between components is used to characterize the intensity of direct thermal interaction caused by the proximity of components, and the specific value is pre-calibrated as a constant based on thermal simulation or design experience of typical engineering scenarios.

[0077] Based on the above settings, the operating temperature of the component can be estimated using the following formula: ; in, Indicates the base temperature of the heat pipe. Indicates the first The power of each component Indicates the first The power of each component Used to represent all related to the first The components whose interval is less than a preset distance threshold are the components of the first component. Additional temperature rise caused by individual components.

[0078] In this embodiment of the invention, based on the above analysis, the following constraint model corresponding to the inter-component thermal coupling constraint is constructed: .

[0079] In this embodiment of the invention, by introducing inter-component thermal coupling constraints, the thermal reliability of the component layout scheme can be significantly enhanced. Potential local overheating problems can be avoided in the initial stage of satellite overall design, thereby improving the overall performance and design quality of the entire satellite component system. Simultaneously, when constructing inter-component thermal coupling constraints, a flexible constraint of "proximity determination → temperature difference limitation" is adopted, rather than a simple geometric distance embargo. This provides the optimization algorithm with a mathematical space for intelligent trade-offs between "layout compactness" and "thermal stability," better meeting the actual needs of engineering design. Through reasonable parameter definition and thorough linearization modeling, complex thermophysical phenomena can be characterized and solved within an efficient mixed-integer programming framework.

[0080] Furthermore, in this embodiment of the invention, it is set that: each heat pipe has the same structural dimensions, and the width of the heat pipe is... , No. The horizontal coordinate of the center position of each heat pipe is .

[0081] It should be noted that the width of the heat pipe refers to its length in the horizontal direction within the component layout area.

[0082] In this embodiment of the invention, the intersection constraint between the component and the heat pipe requires that there be a geometric intersection relationship between the component and all the heat pipes it crosses, and requires that there be no geometric intersection relationship between the component and the other heat pipes.

[0083] In the above-mentioned optimization constraint model for the allocation relationship between components and heat pipes, a binary variable is defined. Used to describe the distribution relationship between components and heat pipes. Indicates the first The component is placed in the first On the heat pipe, that is, the first The component and the first The heat pipes intersect. Indicates the first The component was not placed in the first... On the heat pipe, that is, the first The component and the first The heat pipes do not intersect.

[0084] In this embodiment of the invention, based on the defined binary variable, then the first... The component and the first The intersection constraints of the heat pipes can be modeled as follows: .

[0085] Furthermore, by linearizing the above intersection constraint model, we can obtain the following model: ; in, It should be a relatively large positive number, set according to the specific situation, for example... .

[0086] Furthermore, binary variables are introduced. Expanding the absolute value operation in the linearized intersection constraint model yields the following equivalent constraint model, which is the intersection constraint model between the component and the heat pipe: .

[0087] Furthermore, in this embodiment of the invention, based on the above settings and analysis, the constructed component layout optimization model is expressed as follows: ; in, Indicates the component layout scheme. .

[0088] It should be noted that the constraint models defined in the above component layout optimization model are the same as the constraint models constructed above.

[0089] Furthermore, in this embodiment of the invention, the mathematical programming solver used to solve the component layout optimization model is the SCIP optimization solver or the CPLEX optimization solver; the integer programming algorithm used to solve the component layout optimization model is the branch and bound method.

[0090] The heat pipe-constrained component layout optimization method provided in this invention decomposes the component layout optimization problem under heat pipe constraints into a component layout optimization subproblem with respect to the allocation relationship between components and heat pipes and a component layout optimization subproblem without considering the allocation relationship. Corresponding models are constructed, and then the intersection constraint model between components and heat pipes, constructed based on the allocation relationship between components and heat pipes, is used to solve the models corresponding to the optimization subproblems simultaneously. This method enables rapid optimization of component layout considering quality characteristics and heat pipe constraints. The optimization solution is easy and efficient, and its efficiency is less affected by the component size, allowing for rapid optimization even with large component sizes.

[0091] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Additionally, the terms "front," "back," "left," "right," "upper," and "lower" in this document refer to the placement shown in the accompanying drawings.

[0092] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming, characterized in that, include: Determine the component information and heat pipe information within the component layout area; Based on the structural dimensions of each component and the arrangement of the heat pipes, determine the number of heat pipes required for each component; Using constraints such as the component spanning the heat pipe, the heat pipe's thermal conductivity, and the longitudinal dimension capacity as constraints, discrete variables are used to describe the allocation relationship between the component and the heat pipe. An optimization constraint model for the allocation relationship between the component and the heat pipe is constructed using integer programming. The center position of the component is described by continuous coordinates. The center position of the component is used as the optimization variable. The optimization objective is to minimize the maximum actual heat conduction power of the heat pipe. The constraints are non-interference between components and component layout areas, non-interference between components, centroid constraint of component system, and thermal coupling constraint between components. The optimization objective model and the constraint model corresponding to each constraint condition are constructed by linear modeling. The thermal coupling constraint between components means that the interval between components is not less than a preset distance threshold, and the temperature difference between components with an interval less than the preset distance threshold does not exceed a preset temperature difference threshold. The center position of the component is described by continuous coordinates, and the position of the heat pipe is described by coordinates. Based on the distribution relationship between the component and the heat pipe, the intersection constraint model of the component and the heat pipe is constructed by linear modeling. By combining the component and heat pipe allocation relationship optimization constraint model, optimization objective model, constraint model and component and heat pipe intersection constraint model, the component layout optimization model is obtained. Solving the component layout optimization model, the position of each component in the component layout area is obtained.

2. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 1, characterized in that, The optimal constraint model for the allocation relationship between components and heat pipes is expressed as follows: ； in, This represents an auxiliary binary variable used to describe the first... The reference heat pipe for each component, Indicates the first The heat pipe is the first one The reference heat pipe for each component is defined as the heat pipe with the smallest serial number among all the heat pipes used in the component. Indicates intermediate variables. This represents a binary variable used to describe the distribution relationship between components and heat pipes. Indicates the first The component is placed in the first On a heat pipe, Indicates the first The component was not placed in the first... On a heat pipe, Indicates the first Transformation matrices corresponding to each component The Line number Column elements, transformation matrix For one The matrix, and the matrix's first... The first line The elements up to the first One element is 1, and the rest of the elements of the matrix are 0. Indicates the first The height of each component Indicates the maximum permissible height. Indicates the first The power of each component, Indicates the first The number of heat pipes required for each component This represents the maximum thermal conductivity of a given heat pipe. Indicates the number of components. Indicates the number of heat pipes; Based on the arrangement of the heat pipes, all heat pipes are sequentially numbered starting from number 1 along the arrangement direction.

3. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 2, characterized in that, The bottom left corner of the component layout area is selected as the origin of the coordinate system, and the horizontal direction of the component layout area is used as the coordinate system. The vertical direction of the component layout area is along the axis. Construct a two-dimensional Cartesian coordinate system along the axes; Based on the constructed two-dimensional Cartesian coordinate system, the optimization objective model is represented as: ； in, This indicates a preset parameter. Used for equivalent substitution to optimize the objective function .

4. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 3, characterized in that, The constraint model corresponding to the non-interference constraint between component layout areas is represented as follows: ； in, Indicates the first The center coordinates of each component Indicates the first The width of each component, Indicates the first The height of each component Indicates the width of the component layout area. Indicates the height of the component layout area.

5. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 4, characterized in that, The constraint model corresponding to the non-interference constraint between components is represented as follows: ； in, Indicates an indicator variable. Indicates the first The center coordinates of each component Indicates the first The width of each component, Indicates the first The height of each component It is a positive number.

6. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 5, characterized in that, The constraint model corresponding to the centroid constraint of the component system is represented as follows: ； in, Indicates the first The quality of each component The ordinate represents the position of the expected centroid. The vertical coordinate representing the position of the centroid of the component system. This indicates the pre-set centroid deviation in the coordinate system. The component along the axial direction.

7. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 6, characterized in that, The constraint model corresponding to the thermal coupling constraint between components is represented as follows: ； in, Indicates the first The component and the first Each component in Distance along the axis Indicates the first The first component and the first Each component in Distance along the axis For positive integers, This indicates a preset distance threshold. , , , , It is a binary indicator variable. , Indicates the first Operating temperature of each component Indicates the first Operating temperature of each component This indicates the preset temperature difference threshold. For positive integers, Indicates the base temperature of the heat pipe. Indicates the first The power of each component, Indicates the first The equivalent thermal resistance of each component to its associated heat pipe. This represents the equivalent thermal resistance of the coupling between components. Used to represent all related to the first The components whose interval is less than a preset distance threshold are the components of the first component. Additional temperature rise caused by individual components.

8. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 7, characterized in that, The intersection constraint model between the component and the heat pipe is expressed as: ； in, Indicates the width of the heat pipe. Indicates the first The x-coordinate of the center position of each heat pipe It is a binary variable. .

9. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 8, characterized in that, The component layout optimization model is represented as follows: ； in, Indicates the component layout scheme. .

10. The comprehensive optimization method for heat pipe constrained component layout based on mixed integer programming according to claim 1, characterized in that, Solve the component layout optimization model using a mathematical programming solver or an integer programming algorithm.

Images